zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

advectiveFvPatchField, waveTransmissiveFvPatchField: Updated docs

Browse Source 
Resolves bug-report http://bugs.openfoam.org/view.php?id=2373
This commit is contained in:
Henry Weller
2016-12-05 20:07:40 +00:00
parent 94fd633243
commit 6eeba45c21
2 changed files with 7 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

5
src/finiteVolume/fields/fvPatchFields/derived/advective/advectiveFvPatchField.H
View File

@ -29,7 +29,8 @@ Group
Description
This boundary condition provides an advective outflow condition, based on
solving DDt(psi, U) = 0 at the boundary.
solving DDt(W, field) = 0 at the boundary where \c W is the wave velocity
and \c field is the field to which this boundary condition is applied.
The standard (Euler, backward, CrankNicolson, localEuler) time schemes are
supported. Additionally an optional mechanism to relax the value at
@ -37,7 +38,7 @@ Description
switched on by specifying the relaxation length-scale \c lInf and the
far-field value \c fieldInf.
The flow/wave speed at the outlet is provided by the virtual function
The flow/wave speed \c (w) at the outlet is provided by the virtual function
advectionSpeed() the default implementation of which requires the name of
the flux field \c (phi) and optionally the density \c (rho) if the
mass-flux rather than the volumetric-flux is given.

7
src/finiteVolume/fields/fvPatchFields/derived/waveTransmissive/waveTransmissiveFvPatchField.H
View File

@ -29,18 +29,19 @@ Group
Description
This boundary condition provides a wave transmissive outflow condition,
based onsolving DDt(psi, U) = 0 at the boundary.
based on solving DDt(W, field) = 0 at the boundary \c W is the wave velocity
and \c field is the field to which this boundary condition is applied.
The wave speed is calculated using:
\f[
x_p = \frac{\phi_p}{|Sf|} + \sqrt{\frac{\gamma}{\psi_p}}
w_p = \frac{\phi_p}{|Sf|} + \sqrt{\frac{\gamma}{\psi_p}}
\f]
where
\vartable
x_p | patch values
w_p | patch wave speed
\phi_p | patch face flux
\psi_p | patch compressibility
Sf | patch face area vector