Home
Flow Solver
Contains flow solver related guides, such as the reference and API guides.
Flow solver reference
Insert short description here
Boundary conditions
This section outlines different boundary conditions supported by the flow solver and the associated considerations.
Walls
Understand the key equations used by the flow solver to treat different types of walls used as boundary conditions.
Wall boundary conditions for turbulence models
Understand how the flow solver models wall boundary conditions for different turbulence models.
K-omega turbulence model
The flow solver applies specific boundary conditions for the k-omega turbulence model on both smooth and rough walls, using standard and hybrid wall functions to calculate turbulent kinetic energy and other parameters.

Flow Solver
Contains flow solver related guides, such as the reference and API guides.
- Flow solver guides
  Maya HTT's flow solver is included and distributed with Simcenter 3D and Simcenter Femap. This page contains links to flow solver guides, which help you understand how the flow solver works.
- Flow solver reference
  Insert short description here
  - Introduction
    The flow solver computes a solution to the non-linear, partial differential equations for the conservation of mass, momentum, energy, and general scalars in general complex 3D geometries. It uses an element-based finite volume method and iterative Krylov methods to discretize and solve the governing equations.
  - Nomenclature
    The nomenclature table offers a description and sample units or constant values of the variables used in the Flow Solver Reference Manual.
  - Governing equations
    A model can be comprised of multiple flow enclosures. Each enclosure is a separate group of three-dimensional elements that form a separate fluid domain from the remaining geometry. The physical scales of each enclosure, such as the length scale, time scale, pressure, and temperature are independent of the other enclosures. The flow solver solves governing equations separately for each flow enclosure.
  - Source term
    The source term can represent body forces or flow resistance forces to model different physical phenomena.
  - Turbulence models
    In the flow solver, the flow field can be solved as laminar or as turbulent using governing equations.
  - Boundary conditions
    This section outlines different boundary conditions supported by the flow solver and the associated considerations.
    - Walls
      Understand the key equations used by the flow solver to treat different types of walls used as boundary conditions.
      - Wall distance calculation methods
        Understand the method used by the flow solver to calculate the wall distance.
      - Wall functions
        Wall functions describe the flow velocity as a function of distance from the wall within the near-wall region in the turbulent boundary layer. By using wall functions to approximate the mean velocity field in the near-wall region, you avoid mesh requirements needed to resolve the viscous sublayer.
      - Wall functions with roughness effect
        The roughness effect is incorporated in wall functions with parameters such as the specified equivalent sand grain roughness.
      - Wall shear stress for turbulence models
        The calculation of the wall shear stress in turbulent flows depends on the applied wall treatment method. The application of wall functions and their types affects the shear stress calculation.
      - Wall boundary conditions for turbulence models
        Understand how the flow solver models wall boundary conditions for different turbulence models.
        Mixing length turbulence model
        The flow solver calculates the near wall turbulent viscosity using the eddy viscosity equation with the appropriate mixing length l and the strain rate S.
        Standard k-epsilon, RNG k-epsilon, Realizable k-epsilon turbulence models
        Understand how the flow solver calculates the near-wall turbulent viscosity for various k-epsilon turbulence models and for models with immersed boundary meshes.
        K-omega turbulence model
        The flow solver applies specific boundary conditions for the k-omega turbulence model on both smooth and rough walls, using standard and hybrid wall functions to calculate turbulent kinetic energy and other parameters.
        Shear Stress Transport model
        The flow solver applies specific boundary conditions for the Shear Stress Transport (SST) turbulence model on both smooth and rough walls, using standard wall functions, hybrid wall functions, as well as the immersed boundary method to calculate turbulent kinetic energy and other parameters.
        Spalart-Allmaras model
        The flow solver applies a specific boundary conditions for the Spalart-Allmaras turbulence model for all types of wall treatments. It also computes the near wall vorticity magnitude for smooth walls as well as for standard and hybrid wall functions.
      - Wall function for the thermal boundary layer
        The flow solver models heat flux near adiabatic and convecting walls for laminar and turbulent flows using specific boundary conditions and wall functions.
      - Thermal wall function for high speed flows and flows with viscous heating
        For flows with significant viscous heating, such as high speed flows, the heat flux depends on wall temperature and local conditions derived from thermal wall function.
      - Thermal wall function for natural convection
        For natural convection in turbulent flows, the flow solver uses a new h correlation derived from a general temperature wall function for natural convection.
      - Local heat transfer coefficient correlation for thermal wall function
        The flow solver computes the local heat transfer coefficient using the temperature scale as well as the natural convection thermal wall function.
      - Special considerations when calculating local heat transfer coefficient
        There are uncertainties regarding the validity of the correlations found in the literature when the orientation of the solid surface is such that ϕ > 150°. In addition, as the angle ϕ approaches 180° (heated surface facing downward or cooled surface facing upward), laminar flow prevails even for very high Ra numbers. In the flow solver, the methodology adopted for 150° < ϕ <180° is to evaluate laminar and turbulent heat fluxes and select the highest one
      - Slip wall treatments for gases
    - Fans
      In the flow solver, you can use Inlet Flow, Outlet Flow, Internal Fan, and Recirculation Loop fan types to model flow boundary conditions.
    - Openings
      Openings, also called vents, are pressure and temperature specified boundaries. The boundary condition imposed depends on the direction of the flow: inflow or an outflow. The direction of the flow at a vent can vary during the solve.
    - Convective outflow
      The convective outflow boundary condition models flow exiting the fluid domain without specifying the pressure. It lets the flow field exit and enter the flow domain on each element of the boundary as necessary conserving the mass.
    - Bursting membrane and flap
      The bursting membrane and flap boundary conditions model surfaces that open under specified pressure differences or time, allowing fluid flow.
    - High speed flow boundary condition
      You can use a high speed flow boundary condition for inlets or outlets, to specify the Mach number, pressure, and turbulence values at the domain boundary when the flow has velocities above Mach 0.3, which implies that the flow has density changes of more than 5% and is generally treated as compressible.
    - Screens
      The flow solver models planar resistances such as screens by calculating the pressure drop across them.
    - Symmetry boundaries
      A model is symmetric about a plane when the flow on one side of the plane is a mirror image of the flow on the opposite side of the plane. In such case, the use of a symmetry plane condition enables efficient use of computer resources by allowing the numerical solution to be obtained on a fraction of the original domain.
    - Periodic boundaries
      Periodic boundary conditions are used to simulate a flow leaving through a boundary A and entering through a boundary B under identical conditions (velocity, temperature, scalar values, etc).
    - Mixing plane
      A mixing plane boundary condition interfaces two or more fluid volumes with different flow conditions.
    - Rotating frames of reference
      The flow solver uses the frozen rotor method for rotating frames of reference (RFR). This method is useful for models where the flow distribution around the interface between the fixed and the rotating frames is non-uniform.
    - 1D network connected to 3D fluid region
      In a coupled thermal-flow analysis, you can couple a 3D fluid domain to a set of thermal streams and ducts with mass flow from a 1D fluid network that has the same fluid material. The flow solver solves the equations on the 3D fluid domain using the information from the 1D fluid network solved by the thermal solver.
  - Modeling a gas mixture
    The flow solver uses the scalar equation to model a gas mixed in any proportion with the main gas. The software supports up to five gases in the mixture. All gases are assumed to behave as ideal gases using the ideal gas law.
  - Modeling humidity
    The flow solver uses the gas mixture scalar equation to model humidity that is defined as a mixture of water vapor in air.
  - Modeling non-Newtonian fluid
  - Modeling particle tracking
    The following section describes the equations and forces that model particle trajectories, including treatment and effects of boundary conditions
  - Modeling fluid structure interaction
    Fluid structure interaction (FSI) computations account for motion or deformation of the solid structure boundaries that are interacting with the adjacent fluid flow. The flow solver handles both transient and static fluid structure interactions.
  - Discretization
    The following section describes the discretization methods used by the flow solver to discretize the governing equations.
  - Solver numerical methods
    The following section describes numerical methods used by the flow solver such as the linear equation solution, which acts as a preconditioning method to the basic iterative algorithm, and other solver methods, such as Krylov methods.
  - Flow solver files
    The following table describes the flow solver files, which are written in the run directory during the solving process.
  - Flow solver bibliography
    The following bibliography lists all sources used to explain how the flow solver works.
- Flow solver flowcharts
  This section consists of flowcharts that illustrate the methodology of the flow solver when you run a simple steady-state or transient simulation, and does not cover all features and situations. It also highlights how the flow solver is coupled to the thermal solver for a simple steady-state or transient coupled thermal-flow simulation.

K-omega turbulence model

The flow solver applies specific boundary conditions for the k-omega turbulence model on both smooth and rough walls, using standard and hybrid wall functions to calculate turbulent kinetic energy and other parameters.

The flow solver uses the following values for the wall resolved k-omega turbulence model boundary conditions on a smooth wall:

The boundary conditions for the wall resolved k-omega turbulence model on a rough wall are given by:

where k⁺_s is the equivalent sand grain roughness height in wall units and is calculated as follow:

Standard wall function

For k-omega turbulence model that uses the standard wall function, the flow solver calculates the boundary conditions similarly as to the k-epsilon turbulence model, with exception that instead of ε value, it uses ω value:

where:

For the k equation, the convective and diffusive fluxes on the wall boundary are set to zero.

Hybrid wall function on smooth walls

The turbulent kinetic energy is zero on the wall.

The near-wall solution for ω is defined as described in [36] and can be written in the wall units as follow:

where

Hybrid wall function on rough walls

Turbulent kinetic energy is set to zero on rough walls. The boundary condition for ω is a blend of the resolved rough wall value ω⁺_rough and the log layer value ω⁺_log. Unlike the smooth resolved boundary value, ω⁺_vis, the rough wall boundary value, ω⁺_rough, is independent of y⁺.

Plot comparing blended, rough wall implementation, and log layer-based boundary values in orange, blue, and green, respectively, for an equivalent sand grain roughness height. The x-axis represents the dimensionless distance, and the y-axis represents the magnitude of the boundary value.