mod_momentum_boundary_conditions Module Reference

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. More...

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_bndc	the density (cell) boundary conditions
[in]	velocity_bndc	the velocity (face) boundary conditions
[in]	momentum_bndc	the resulting momentum (face) boundary conditions