The Brinkman Equations Interface
The Brinkman Equations (br) interface (), found under the Porous Media and Subsurface Flow branch () when adding a physics interface, is used to compute fluid velocity and pressure fields of single-phase flow in porous media in the laminar flow regime. The physics interface extends Darcy’s law to describe the dissipation of the kinetic energy by viscous shear, similar to the Navier–Stokes equations. Fluids with varying density can be included at Mach numbers below 0.3. Also the viscosity of a fluid can vary, for example, to describe non-Newtonian fluids. To simplify the equations, select the Stokes–Brinkman flow feature to reduce the dependence on inertial effects when the Reynolds number is significantly less than 1. The physics interface can be used for stationary and time-dependent analyses.
The main node is the Porous Medium feature, which adds the Brinkman equations and provides interfaces for defining the fluid material and the porous matrix.
When this physics interface is added, the following default nodes are also added in the Model BuilderPorous Medium, Wall (the default boundary type, using No slip as the default boundary condition), and Initial Values. Then, from the Physics toolbar, add other nodes that implement, for example, boundary conditions and volume forces. You can also right-click Brinkman Equations to select physics features from the context menu.
The boundary conditions are essentially the same as for the Laminar Flow interface. Differences exist for the following boundary types: Outlet, Symmetry, Open Boundary, and Boundary Stress where the viscous part of the stress is divided by the porosity to appear as
In the COMSOL Multiphysics Reference Manual see Table 2-4 for links to common sections such as Discretization, Consistent Stabilization, and Inconsistent Stabilization, and Advanced Settings sections, all accessed by clicking the Show button () and choosing the applicable option. You can also search for information: press F1 to open the Help window or Ctrl+F1 to open the Documentation window.
Settings
The Label is the default physics interface name.
The Name is used primarily as a scope prefix for variables defined by the physics interface. Refer to such physics interface variables in expressions using the pattern <name>.<variable_name>. In order to distinguish between variables belonging to different physics interfaces, the name string must be unique. Only letters, numbers, and underscores (_) are permitted in the Name field. The first character must be a letter.
The default Name (for the first physics interface in the model) is br.
Physical Model
This node specifies the properties of the Brinkman Equations interface, which describe the overall type of fluid flow model.
Compressibility
By default the physics interface uses the Incompressible flow formulation of the Brinkman equations to model constant density flow. Alternatively, from the Compressibility list select Weakly compressible flow which should be used when the pressure dependency of the density can be neglected, or Compressible flow (Ma<0.3) to solve for the full compressible flow of which, however, the Mach number must be below 0.3.
Swirl Flow
For 2D axisymmetric models, select the Swirl flow check box to include the swirl velocity component, that is the velocity component in the azimuthal direction. While can be nonzero, there can be no gradients in the direction.
Note that this feature is only available for specific modules. See https://www.comsol.com/products/specifications/ for a detailed overview.
Neglect Inertial Term (Stokes–Brinkman Flow)
The Neglect inertial term (Stokes–Brinkman) check box is selected by default to model flow at low Reynolds numbers for which the inertial term can be neglected. This results in the linear Stokes–Brinkman equations.
Include Gravity
Gravity is not included by default. Select the Include gravity check box to activate the acceleration of gravity. This automatically adds a global Gravity feature node to the interface model tree, and the buoyancy force is included in the Equations.
Also, when the Include gravity check box is selected, the Use reduced pressure option changes the pressure formulation from using the total pressure (default) to using the reduced pressure. This option is suitable for configurations where the density changes are very small; otherwise, the default formulation can be used. For more information, see Gravity.
Porous Treatment of No Slip Condition
Choose how the No Slip condition on Wall boundaries and Interior Wall boundaries adjacent to Porous Medium should be treated. The options are Standard no slip formulation (default) and Porous slip. The latter option provides a unified treatment when the porous matrix is fully resolved as well as when it is under resolved ensuring a smooth transition between regions with different resolutions; see No Slip under Wall in the The Creeping Flow, Laminar Flow, Turbulent Flow, and Large Eddy Simulation Interfaces section.
Reference Pressure Level
Enter a Reference pressure level pref (SI unit: Pa). The default value is 1[atm].
Reference Temperature Level
Enter a Reference temperature level Tref (SI unit: K). The default value is 293.15[K].
Reference Position
If Include gravity is selected, a Reference position rref (SI unit: m) can be specified which is then used for the calculation of the hydrostatic pressure.
Turbulence
Turbulent flow can be simulated by changing the Turbulence model type to RANS (Reynolds-Averaged Navier–Stokes). If turbulent flow is activated, one of the following Turbulence models can be chosen: Algebraic yPlus, L-VEL, and k-ε, and Wall functions are preset for Wall treatment. For more information about turbulence modeling, see Theory for the Turbulent Flow Interfaces in the CFD Module User’s Guide.
Dependent Variables
The following dependent variables (fields) are defined for this physics interface — the Velocity field u (SI unit: m/s) and its components, and the Pressure p (SI unit: Pa).
Advanced Settings
To display this section, click the Show More Options button () and select Advanced Physics Options in the Show More Options dialog box. Normally these settings do not need to be changed.
The Use pseudo time stepping for stationary equation form option adds pseudo time derivatives to the equation when the Stationary equation form is used in order to speed up convergence. When selected, a CFL number expression should also be defined. For the default Automatic option, the local CFL number (from the Courant–Friedrichs–Lewy condition) is determined by a PID regulator.
Pseudo Time Stepping for Laminar Flow Models and Pseudo Time Stepping in the COMSOL Multiphysics Reference Manual