radGradTotalPressure

Total pressure boundary condition with user-specified radial gradient for OpenFOAM. Adapted from Foam::rotatingTotalPressureFvPatchScalarField. (https://github.com/OpenFOAM/OpenFOAM-7/tree/master/src/finiteVolume/fields/fvPatchFields/derived/rotatingTotalPressure)

Installation

1. Clone the repository into $WM_PROJECT_USER_DIR/src/finiteVolume/fields/fvPatchFields/derived.
2. Copy the make directory to $WM_PROJECT_USER_DIR/src/finiteVolumne.
3. Run wmake libso in $WM_PROJECT_USER_DIR/src/finiteVolume.

TODO

• Rename and refactor .C and .h files to use radGrad instead of radEq.

• Set up test cases and testing with github actions.

• Add support for general rotation axes (currently limited to z-axis rotation only).

Usage

Tell OpenFOAM to dynamically link the library by including

libs
(
    "libradGradTotalPressureFvPatchScalarField.so"
);

in <MY_CASE>/system/controlDict.

In the file <MY_CASE>/0/p or <MY_CASE>/0/p_rgh, the boundary condition can be set using e.g.:

<PATCH_NAME>
{
    type            radGradTotalPressure;
    p0              uniform 100000;
    rRef            0.1;
    dp0dr           1000;
    rho             thermo:rho;
    psi             thermo:psi;
    gamma           1.4;
    omega           (0 0 100);
}

Property	Description	Required	Default value
U	velocity field name	no	U
phi	flux field name	no	phi
rho	density field name	no	none
psi	compressibility field name	no	none
gamma	ratio of specific heats (Cp/Cv)	yes	none
p0	static pressure reference	yes	none
omega	angular velocity of the frame [rad/s]	yes	none
dp0dr	radial total pressure gradient	yes	none
rRef	reference radial location for p0	yes	none