Dimensionless Numbers Important for Solver Stability
When solving a microfluidics problem numerically there are two critical dimensionless numbers to consider in terms of solver stability. These numbers are:
Each of these numbers is defined in Table 2-2. Both the Reynolds number and the Peclet number are associated with the relative importance of convective terms in the corresponding partial differential equation. The Peclet number describes the importance of convection in relation to diffusion (for either heat or mass transfer), and the Reynolds number describes the importance of the “convective” inertia term in relation to viscosity in the Navier–Stokes equations themselves.
Both the Reynolds number and the Peclet number can be defined on the “cell” or element level. As they are defined in COMSOL Multiphysics numerical instabilities can arise when the cell Reynolds or Peclet number is greater than one. These instabilities are usually manifested as spurious oscillations in the solution. Taking the Peclet number as an example, oscillations can occur when the cell Peclet number is greater than one in the following circumstances:
In theory the grid can be refined to bring the cell Reynolds or Peclet number below one, although this is often impractical for many problems. Several stabilization techniques are included, which enable problems with larger cell Reynolds or Peclet numbers to be solved. At the crudest level additional numerical diffusion can be added to the problem to improve its stability. This is achieved by selecting Isotropic Diffusion under Inconsistent Stabilization for any physics interface.
The stabilization options are visible when the Show More Options button () is clicked and Stabilization is selected in the Show More Options dialog box.
This method is termed “inconsistent” as a solution to the problem without numerical diffusion is not necessarily a solution to the problem with diffusion. COMSOL Multiphysics also has consistent stabilization options. A consistent stabilization technique reduces the numerical diffusion added to the problem as the solution approaches the exact solution. Both streamline diffusion and crosswind diffusion are available. Streamline diffusion adds numerical diffusion along the direction of the flow velocity (that is, the diffusion is parallel to the streamlines). Crosswind diffusion adds diffusion in the direction orthogonal to the velocity.
Generally it is best to use consistent stabilization where possible. If convergence problems are still encountered, inconsistent stabilization can be used with a parametric or time-dependent solver that slowly eliminates this term.
Selection Information and Numerical Stabilization in the COMSOL Multiphysics Reference Manual
Re
Pe
Pe
Symbol definitions: v is the characteristic velocity or cell velocity, r is the fluid density, L is the characteristic length scale for the problem, μ is the fluid viscosity, h is the element size, cp is the heat capacity at constant pressure, κ is the thermal conductivity, and D is the diffusion constant. For anisotropic diffusion or conductivity an appropriate average is computed. The variables names assume the default name spf.