To enable this setting for some physics interfaces and features, click the Show button () and select Discretization. There are two categories of discretization — a section on the physics interface node’s Settings window (described here) and adding a Discretization (Node) for global equation-based modeling.
The software adapts the order of the numerical integration to the element orders for the physics in the model. Some physics interfaces use special element types or a reduced element order for some of the field variables. Select the Shape function type and the Element order as, in most cases, Linear, Quadratic, Cubic, Quartic, or Quintic (for order 1–5, respectively).Table 16-2 is an overview of the available shape function types and the element orders supported.
Table 16-2: Shape Function Types
Some additional information is included Shape Functions and Element Types in the COMSOL Multiphysics Programming Reference Manual.
• P1+P1 means linear elements for both the velocity components and the pressure field. Linear elements are computationally cheaper than higher-order elements and are also less prone to introducing spurious oscillations, thereby improving the numerical robustness. In other words, this can be computationally efficient but requires streamline stabilization of the Navier-Stokes equations. This is the default element order for the Laminar Flow and Turbulent Flow single-phase flow interfaces and the discretization of fluids in the multiphase flow interfaces.
• P2+P1 means second-order elements for the velocity components and linear elements for the pressure field. Second-order elements work well for low flow velocities. This is the default element order for the Creeping Flow interface.
• P3+P2 means third-order elements for the velocity components and second-order elements for the pressure field. This can add additional accuracy but it also adds additional degrees of freedom compared to P2+P1 elements.
The abbreviation PmPn is often used to indicate the polynomial order of, in this case, the shape functions (elements) for the velocity components (m) and the pressure (n) when using tetrahedral or triangular elements. Here a corresponding nomenclature is used for all element shapes.The theory about this is in P.M. Gresho and R.L. Sani, Incompressible Flow and the Finite Element Method, Volume 2: Isothermal Laminar Flow, John Wiley & Sons, 2000.
The discretization of the temperature field follows that of the fluid components, so the temperature order is 1 for P1+P1, 2 for P2+P1, and 3 for P3+P2.
Some physics can create and compute variables that accurately represent the flux across all boundaries. To enable these variables, select the Compute boundary fluxes check box. Optionally, the smoothing can provide a more well-behaved flux value close to singularities. You add smoothing by selecting the Apply smoothing to boundary fluxes check box.
From the Value type when using splitting of complex variables list, you can specify the value type (Real or Complex) of dependent variables when the Split complex variables in real and imaginary parts setting is activated in the Compile Equations node of any solver sequence used. The default is the complex value type, but you can specify that the value of a dependent variable is real to make sure that it does not get affected by small imaginary contributions, which can occur, for example, when combining a Time Dependent or Stationary study with a frequency-domain study. If the split complex variables setting is not active, the value type is ignored.
By default, the general number formatting algorithm for complex-valued numbers in COMSOL Multiphysics ignores small real or imaginary parts. If you want to see also such small real or imaginary parts, open the Preferences dialog box, and under Precision on the General page, clear the Suppress small real or imaginary part check box.