version 0.6.0
Loading...
Searching...
No Matches

Boundary conditions related routines for momentum. More...

Functions/Subroutines

subroutine generate_momentum_boundary_conditions (density_bndc, velocity_bndc, momentum_bndc)
 Generate the momentum boundary conditions based on \( \rho \) and \( \mathbf{ \mathrm{ u } } \) ones.
 

Detailed Description

Momentum in the form \( \rho \mathbf{ \mathrm{ u } } \) might need appropriate boundary conditions (ie. for the full inertial term). It is simply the product of \( \rho \) and \( \mathbf{ \mathrm{ u } } \) boundary conditions (done with some care).

Function/Subroutine Documentation

◆ generate_momentum_boundary_conditions()

subroutine mod_momentum_boundary_conditions::generate_momentum_boundary_conditions ( type(t_boundary_condition), intent(in) density_bndc,
type(t_boundary_condition_face), intent(in) velocity_bndc,
type(t_boundary_condition_face), intent(inout) momentum_bndc )

We base the treatment on the velocity bndc which guide the momentum bndc. To resume: only inlets and moving boundaries are treated differently as we have to multiply the velocity by the density at that place.

Warning
this method only holds when using classical boundary conditions. It is incompatible with non homogeneous Neumann (or round Robin) conditions. In that case, the function will quit.
Parameters
[in]density_bndcthe density (cell) boundary conditions
[in]velocity_bndcthe velocity (face) boundary conditions
[in]momentum_bndcthe resulting momentum (face) boundary conditions