The Mach Number Limit
An important dimensionless number in fluid dynamics is the Mach number, Ma, defined by
where a is the speed of sound. A flow is formally incompressible when Ma = 0. This is theoretically achieved by letting the speed of sound tend to infinity. The Navier-Stokes equations then have the mathematical property that pressure disturbances are instantaneously propagated throughout the entire domain. This results in a parabolic equation system.
The momentum equation, Equation 4-15, is parabolic for unsteady flow and elliptic for steady flow, whereas the continuity equation, Equation 4-14, is hyperbolic for both steady and unsteady flow. The combined system of equations is thus hybrid parabolic-hyperbolic for unsteady flow and hybrid elliptic-hyperbolic for steady flow. An exception occurs when the viscous term in Equation 4-15 becomes vanishingly small, such as at an outflow boundary, in which case the momentum equation becomes locally hyperbolic. The number of boundary conditions to apply on the boundary then depends on the number of characteristics propagating into the computational domain. For the purely hyperbolic system, the number of characteristics propagating from the boundary into the domain changes as the Mach number passes through unity. Hence, the number of boundary conditions required to obtain a numerically well-posed system must also change. The compressible formulation of the laminar and turbulent interfaces uses the same boundary conditions as the incompressible formulation, which implies that the compressible interfaces are not suitable for flows with a Mach number larger than or equal to one.
The practical Mach number limit is lower than one, however. The main reason is that the numerical scheme (stabilization and boundary conditions) of the Laminar Flow interface does not recognize the direction and speed of pressure waves. The fully compressible Navier-Stokes equations do, for example, start to display very sharp gradients already at moderate Mach numbers. But the stabilization for the single-phase flow interface does not necessarily capture these gradients. It is impossible to give an exact limit where the low Mach number regime ends and the moderate Mach number regime begins, but a rule of thumb is that the Mach number effects start to appear at Ma = 0.3. For this reason, the compressible formulation is referred to as Compressible flow (Ma<0.3) in COMSOL Multiphysics. To model high Mach number flows, use one of The High Mach Number Flow Interfaces.