The CFD Module Physics Interfaces
The Fluid Flow interfaces in this module are based on the laws for conservation of momentum, mass, and energy in fluids. The different flow models contain different combinations and formulations of the conservation laws that apply to the physics of the flow field. These laws of physics are translated into partial differential equations and are solved together with the specified initial and boundary conditions.
A physics interface defines a number of features. These features are used to specify the fluid properties, initial conditions, boundary conditions, and possible constraints. Each feature represents an operation describing a term or condition in the conservation equations. Such a term or condition can be defined on a geometric entity of the component, such as a domain, boundary, edge (for 2D components), or point.
Figure 5 shows the Model Builder, including a Laminar Flow interface, and the Settings window for the selected Fluid Properties 1 feature node. The Fluid Properties 1 node adds the marked terms to the component equations in a selected geometry domain. Furthermore, the Fluid Properties 1 feature may link to the Materials feature node to obtain physical properties such as density and dynamic viscosity, in this case the fluid properties of water. The fluid properties, defined by the Water, liquid material, can be functions of the modeled physical quantities, such as pressure and temperature. In the same way, the Wall 1 node adds the boundary conditions at the walls of the fluid domain.
Figure 5: The Model Builder including a Laminar Flow interface (left), and the Settings window for Fluid Properties for the selected feature node (right). The Equation section in the Settings window shows the component equations and the terms added by the Fluid Properties 1 node. The added terms are underlined with a dotted line. The arrows also explain the link between the Materials node and the values for the fluid properties.
The CFD Module includes a large number of Fluid Flow interfaces for different types of flow. It also includes Chemical Species Transport interfaces for reacting flows in multicomponent solutions, and physics interfaces for heat transfer in solids, fluids, and porous media found under the Heat Transfer branch.
Figure 6 shows the Fluid Flow interfaces as they are displayed when you add a physics interface in 3D (see also Physics Interface Guide by Space Dimension and Study Type for further information). A short description of the physics interfaces follows.
Figure 6: The physics interfaces for the CFD Module as shown in the Model Wizard.
Single-Phase Flow
The Creeping Flow interface () approximates the Navier–Stokes equations for very low Reynolds numbers. This is often referred to as Stokes flow and is applicable when viscous effects are dominant, such as in very small channels or microfluidics devices.
The Laminar Flow interface () is primarily applied to flows at low to intermediate Reynolds numbers. This physics interface solves the Navier–Stokes equations for incompressible, weakly compressible, and compressible flows (up to Mach 0.3). The Laminar Flow interface also allows for simulation of non-Newtonian fluids.
The physics interfaces under the Turbulent Flow branch () model flows at high Reynolds numbers. These physics interfaces solve the Reynolds-averaged Navier–Stokes (RANS) equations for the averaged velocity and pressure fields. The turbulent flow interfaces provide different options for modeling the turbulent viscosity. There are several turbulence models available—two algebraic turbulence models, the Algebraic yPlus and L-VEL models, and seven transport-equation models, including a standard k-ε model, the Realizable k-ε model, a k-ω model, an SST (Shear Stress Transport) model, a low Reynolds number k-ε model, the Spalart–Allmaras model, and the v2-f model. Similarly to the Laminar Flow interface, incompressible flow is selected by default.
The Algebraic yPlus and L-VEL turbulence models are so-called enhanced viscosity models. A turbulent viscosity is computed from the local distance to the nearest wall. For this reason, the algebraic turbulence models are best suited for internal flows, such as in electronic cooling applications. Algebraic turbulence models are computationally economical, and more robust but, in general, less accurate than transport-equation models.
Among the transport-equation turbulence models, the standard k-ε model is the most widely used since it often is a good compromise between accuracy and computational cost (memory and CPU time). The Realizable k-ε model is similar to the standard k-ε model but has built-in realizability constraints, resulting in improved performance for certain flows, such as turbulent jets. The k-ω model is an alternative to the standard k-ε model and often gives more accurate results, especially in recirculation regions and close to solid walls. However, the k-ω model is also less robust than the standard k-ε model.
The SST model combines the robustness of the k-ε model with the accuracy of the k-ω model, making it applicable to a wide variety of turbulent flows. The Low Reynolds number k-ε model is more accurate than the standard k-ε model, especially close to walls.
The Spalart-Allmaras model is specifically designed for aerodynamic applications, such as flow around wing profiles, but is also widely used in other applications due to its high robustness and decent accuracy.
In the v2-f model, the turbulent viscosity is based on the wall-normal velocity fluctuations, whereby wall blockage effects and low Reynolds number effects are captured separately. The v2-f model also includes non-local effects of the fluctuating pressure on the turbulent fields. Higher resolution is needed in the near-wall region for the SST model, the Low Reynolds number k-ε model, the Spalart-Allmaras model, and the v2-f model. Thus, the better accuracy provided by these models comes at a higher computational cost.
The physics interfaces under the Large Eddy Simulation branch () also model flows at high Reynolds numbers. For these interfaces, the larger unsteady, three-dimensional eddies are fully resolved whereas the smaller eddies are modeled. The dependent variables, velocity and pressure, are divided into resolved and unresolved scales, and are substituted into the Navier–Stokes and continuity equations. The additional stress terms involving the unresolved scales are modeled in various ways in the three different interfaces under this branch. For the Residual Based Variational Multiscale (RBVM) interface and the Residual Based Variational Multiscale with Viscosity (RBVMWV) interface, the unresolved scales are modeled in term of the equation residuals, which means that, in the limit when all scales are resolved, the simulation reverts to a direct numerical simulation (DNS) of the Navier–Stokes and the continuity equations. The third interface uses a Smagorinsky term to model the effect of the unresolved Reynolds stress on the smallest resolved scales. This model reverts to DNS in the limit when the large resolved scales comprise the whole spectrum of turbulent scales. The Large Eddy Simulation interfaces are applicable to incompressible, three-dimensional, time-dependent flows.
The Rotating Machinery interfaces () combine the Single-Phase Flow interfaces and a Rotating Domain, and are applicable to fluid-flow problems where one or more of the boundaries rotate, for example in mixers and around propellers. The physics interfaces support incompressible, weakly compressible and compressible (Mach < 0.3) flows, laminar Newtonian and non-Newtonian flows, and turbulent flow using the standard k-ε model or either of the two algebraic turbulence models (Algebraic yPlus or L-VEL).
The Viscoelastic Flow interface () is used to simulate incompressible isothermal flow of viscoelastic fluids. It solves the continuity equation, the momentum equation and a constitutive equation that defines the elastic stresses. There are three predefined models for the elastic stresses: Oldroyd-B, FENE-P and Giesekus.
Multiphase Flow
The physics interfaces under the Bubbly Flow branch () model two-phase flow where the fluids form a gas-liquid mixture, and the content of the gas is less than 10%. There is support for both laminar flow and turbulent flow using extended versions of the turbulence models that accounts for bubble-induced turbulence. For laminar flow, the physics interface supports non-Newtonian liquids. The Bubbly Flow interfaces also allow for mass transfer between the two phases.
The physics interfaces under the Mixture Model branch () are similar to the Bubbly Flow interfaces but assume that the dispersed phase consists of solid particles or liquid droplets. The continuous phase has to be a liquid. There is support for both laminar flow and turbulent flow. The Mixture Model interfaces also allow for mass transfer between the two phases.
The Euler-Euler Model interface () for two-phase flow is able to handle the same cases as the Bubbly Flow and Mixture Model interfaces, but is not limited to low concentrations of the dispersed phase. In addition, the Euler-Euler Model interface can handle large differences in density between the phases, such as the case of solid particles in air. This makes the model suitable for simulations of, for example, fluidized beds. There is support for both laminar flow and turbulent flow using either a mixture or phase-specific k-ε turbulence model.
The predefined multiphysics couplings, Phase Transport Mixture Model (), couple a Laminar Flow or Turbulent Flow (RANS) interface to the Phase Transport interface, providing mass-averaged formulations of the mixture model interface that can handle an arbitrary number of dispersed phases.
The Nonisothermal Mixture Model () interfaces couple a Laminar Flow or Turbulent Flow (RANS) interface, the Phase Transport interface and a Heat Transfer in Fluids interface using the three multiphysics coupling nodes: Nonisothermal Flow (), Mixture Model () and Nonisothermal Mixture Model (). They can be used for studies of boiling and cavitation.
The Two-Phase Flow, Level Set interface (), the Two-Phase Flow, Phase Field interface (), and the Two Phase Flow, Moving Mesh interface () are used to model two fluids separated by a fluid-fluid interface. The moving interface is tracked in detail using either the level set method, the phase field method, or by a moving mesh, respectively. The level set and phase field methods use a fixed mesh and solve additional equations to track the interface location. The moving mesh method solves the Navier Stokes equations on a moving mesh with boundary conditions to represent the interface. In this case equations must be solved for the mesh deformation. Since a surface in the geometry is used to represent the interface between the two fluids in the Moving Mesh interface, the interface itself cannot break up into multiple disconnected surfaces. This means that the Moving Mesh interface cannot be applied to problems such as droplet formation in inkjet devices (in these applications the level set or phase field interfaces are appropriate). These physics interfaces support incompressible flows, where one or both fluids can be non-Newtonian. The physics interfaces also support turbulent flow.
The Two-Phase Thin-Film Flow, Phase Field interfaces () are multiphysics interfaces which couple a Thin Film Flow interface to a Phase Field interface to model two-phase flow in thin layers. Both of the coupled physics interfaces are boundary physics interfaces, which means that the boundary level is the highest level; they do not have a domain level.
The Laminar Three-Phase Flow, Phase Field interface () models laminar flow of three incompressible phases which may be either Newtonian or non-Newtonian. The moving fluid-fluid interfaces between the three phases are tracked in detail using the phase-field method.
The Phase Transport interface ()can be used to model transport of an arbitrary number of phases. The transport velocity must be obtained from one of the momentum transport interfaces (that is, from the single-phase or multiphase flow interfaces) or be specified manually.
The Phase Transport in Porous Media interface () can be used to model transport of an arbitrary number of phases in a porous medium. The transport velocity is derived from an absolute pressure, which must be supplied.
Porous Media and Subsurface Flow
The Brinkman Equations interface () models flow through porous media where the influence of shear stresses is significant. This physics interface supports the Stokes-Brinkman formulation, suitable for very low flow velocities, as well as the full Brinkman equations including convective terms and Forchheimer drag, which is used to account for effects at high interstitial velocities. The flow can be either incompressible or compressible, provided the Mach number is less than 0.3.
The Darcy’s Law interface () models relatively slow flows through porous media for cases where the effects of shear stresses perpendicular to the flow are small.
The Multiphase Flow in Porous Media interface () can be used to model transport of an arbitrary number of phases in a porous medium where the velocity is obtained from Darcy’s law. This is a multiphysics interface for which the couplings between fluid flow and phase transport are set up automatically.
The Two-Phase Darcy’s Law interface () sets up two Darcy-law equations, one for each fluid phase in the porous medium. It couples the two, for example using a capillary expression. It is tailored to model effects such as moisture transport in porous media.
The Free and Porous Media Flow interface ()models porous media containing open channels connected to the porous media, such as in fixed-bed reactors and catalytic converters.
Nonisothermal Flow
The Nonisothermal Flow, Laminar Flow interface () is primarily applied to model flow at low to intermediate Reynolds numbers in situations where the temperature and flow fields have to be coupled. A typical example is natural convection, where thermal buoyancy forces drive the flow. This is a multiphysics interface for which the nonlocal couplings between fluid flow and heat transfer are set up automatically.
The Nonisothermal Flow, Turbulent Flow interfaces () solve the Reynolds-Averaged Navier–Stokes (RANS) equations coupled to heat transfer in fluids and in solids. There is support for all the fluid-flow turbulence models – the Algebraic yPlus model, the L-VEL model, the standard k-ε model, the Realizable k-ε model, a k-ω model, an SST model, a low Reynolds number k-ε model, the Spalart-Allmaras model, and the v2-f model.
The Nonisothermal Flow, LES interfaces () couple one of the three LES interfaces to a Heat Transfer interface. These interfaces are applicable to incompressible flow, and buoyant convection can be modeled using the Boussinesq approximation.
The Brinkman Equations interface () couples the Brinkman Equations interface to a Heat Transfer in Porous Media interface and automatically adds a Porous Material node under Materials in the Model Tree.
The Conjugate Heat Transfer interfaces () are also included with the CFD Module and are almost identical to the Nonisothermal Flow interfaces. They only differ in the default domain feature selected - Heat transfer in Solids instead of Fluid.
High Mach Number Flow
The High Mach Number Flow, Laminar Flow interface () solves the continuity, momentum, and energy equations for fully compressible laminar flow. This physics interface is typically used to model low-pressure systems, for which the Mach number can be large but the flow remains laminar.
The High Mach Number Flow, Turbulent Flow interfaces () solve the continuity, momentum and energy equations for the averaged flow variables in fully compressible turbulent flow, coupled to transport equations for the turbulence quantities. There are two versions: one that uses to the k-ε turbulence model and one that uses to the Spalart-Allmaras turbulence model.
The Compressible Euler Equations interface () solves the equations for transient inviscid ideal gas flow where the velocity magnitude is comparable to the speed of sound, that is, compressible flows in the transonic and supersonic range.
The Rotating Machinery, High Mach Number Flow interfaces () combine the High Mach Number Flow interfaces and a Rotating Domain, and are applicable to fluid-flow problems where one or more of the boundaries rotate, for example in turbines and around propellers. The physics interfaces support compressible, laminar, and turbulent flow using the k-ε model or the Spalart-Allmaras model.
Fluid-structure interaction
The Fluid-Structure Interaction interfaces () couple a Single-Phase Flow or Two-Phase Flow, Phase Field interface to a Solid Mechanics interface for studies of deformation induced by fluid forces.
Thin-Film Flow
The Thin Film Flow interfaces () model the flow of liquids or gases confined in a thin layer on a surface. Applying equations defined on a surface, these physics interfaces compute the average velocity and pressure across the layer in narrow planar structures. The physics interfaces are thus boundary physics interfaces, which means that the boundary level is the highest level; they do not have a domain level. The simulation of the flow of a lubrication oil between two rotating cylinders is an example of a possible application for this physics interface.
Shallow Water Equations
The Shallow Water Equations interface () models free-surface flows in 1D and 2D domains using a depth-averaged formulation. The interface is based on the discontinuous Galerkin method (dG-FEM) and uses a time explicit solver.
Reacting Flow
The Laminar Flow interface () and the Laminar Flow, Diluted Species interface () under the Reacting Flow branch combine the functionality of the Single-Phase Flow and Transport of Concentrated Species interfaces or Transport of Diluted Species respectively. These physics interfaces are primarily applied to model flow at low to intermediate Reynolds numbers in situations where the mass transport and flow fields have to be coupled.
The Turbulent Flow interfaces () and the Turbulent Flow, Diluted Species interfaces () under the Reacting Flow branch apply the Reynolds-Averaged Navier–Stokes (RANS) equations together with the functionality in the Transport of Concentrated Species interface or Transport of Diluted Species respectively. They model mass and momentum transport in turbulent reacting fluid flow. The supported RANS models include the standard k-ε model, a k-ω model, the SST model, and a low Reynolds number k-ε model.
Reacting Flow in Porous Media
The Reacting Flow in Porous Media, Transport of Diluted Species interface () and the Reacting Flow in Porous Media, Transport of Concentrated Species interface () model diluted and concentrated reacting mixtures, respectively, transported by free and porous media flow. Effective diffusion coefficients in a porous matrix can be calculated from the porosity.
Physics Interface Guide by Space Dimension and Study Type
 
Chemical Species Transport
Reacting Flow
Turbulent Flow
Turbulent Flow, Diluted Species
Reacting Flow in Porous Media
Fluid Flow
Single-Phase Flow
Turbulent Flow
Large Eddy Simulation
Rotating Machinery, Fluid Flow
Multiphase Flow
Bubbly Flow
Bubbly Flow, Turbulent Flow
Mixture Model
Mixture Model, Turbulent Flow
Euler-Euler Model
Phase Transport Mixture Model
Nonisothermal Mixture Model
Two-Phase Flow, Moving Mesh
Two-Phase Flow, Level Set
Two-Phase Flow, Phase Field
Two-Phase Thin-Film Flow Phase Flow, Phase Field
Three-Phase Flow, Phase Field
Phase Transport
Porous Media and Subsurface Flow
Nonisothermal Flow
Turbulent Flow
Large Eddy Simulation
High Mach Number Flow
Rotating Machinery, High Mach Number Flow
Fluid-Structure Interaction
Thin-Film Flow
Shallow Water Equations
Heat Transfer
Conjugate Heat Transfer
Turbulent Flow
Mathematics
Moving Interface
1 This physics interface is included with the core COMSOL package but has added functionality for this module.
2 This physics interface is a predefined multiphysics coupling that automatically adds all the physics interfaces and coupling features required.