diff --git a/src/finiteVolume/cfdTools/general/SRF/derivedFvPatchFields/SRFFreestreamVelocityFvPatchVectorField/SRFFreestreamVelocityFvPatchVectorField.H b/src/finiteVolume/cfdTools/general/SRF/derivedFvPatchFields/SRFFreestreamVelocityFvPatchVectorField/SRFFreestreamVelocityFvPatchVectorField.H index fe50e35cac8bff1060efab71259fb8f43e226be0..58c7ab9ce081d103576f51f41df1c165427eb262 100644 --- a/src/finiteVolume/cfdTools/general/SRF/derivedFvPatchFields/SRFFreestreamVelocityFvPatchVectorField/SRFFreestreamVelocityFvPatchVectorField.H +++ b/src/finiteVolume/cfdTools/general/SRF/derivedFvPatchFields/SRFFreestreamVelocityFvPatchVectorField/SRFFreestreamVelocityFvPatchVectorField.H @@ -48,12 +48,17 @@ Description \heading Patch usage + \table + Property | Description | Required | Default value + UInf | free stream velocity | yes | + \endtable + Example of the boundary condition specification: \verbatim myPatch { type SRFFreestreamVelocity; - UInf (0 0 0); // free stream velocity + UInf (0 0 0); value uniform (0 0 0); // optional value entry } \endverbatim diff --git a/src/finiteVolume/cfdTools/general/SRF/derivedFvPatchFields/SRFVelocityFvPatchVectorField/SRFVelocityFvPatchVectorField.H b/src/finiteVolume/cfdTools/general/SRF/derivedFvPatchFields/SRFVelocityFvPatchVectorField/SRFVelocityFvPatchVectorField.H index 71bb989bcd04945ec30276a7b11abb4474436340..32aaa246a76693adea2654d6b0533e9ca06cfe14 100644 --- a/src/finiteVolume/cfdTools/general/SRF/derivedFvPatchFields/SRFVelocityFvPatchVectorField/SRFVelocityFvPatchVectorField.H +++ b/src/finiteVolume/cfdTools/general/SRF/derivedFvPatchFields/SRFVelocityFvPatchVectorField/SRFVelocityFvPatchVectorField.H @@ -41,24 +41,38 @@ Description - relative = yes: inlet velocity applied 'as is': \f[ - Up = U_in + U_p = U_{in} \f] - relative = no : SRF velocity is subtracted from the inlet velocity: \f[ - Up = U_in - U_srf + U_p = U_{in} - U_{p,srf} \f] + where + \vartable + U_p = patch velocity [m/s] + U_{in} = user-specified inlet velocity + U_{p,srf} = SRF velocity + \endvartable + + \heading Patch usage + \table + Property | Description | Required | Default value + relative | relative motion to the SRF? | yes | + inletValue | inlet velocity | yes | + \endtable + Example of the boundary condition specification: \verbatim myPatch { type SRFVelocity; - relative yes; // relative motion to the SRF? - inletValue uniform (0 0 0); // inlet velocity + relative yes; + inletValue uniform (0 0 0); value uniform (0 0 0); // initial value } \endverbatim diff --git a/src/finiteVolume/fields/fvPatchFields/basic/mixed/mixedFvPatchField.H b/src/finiteVolume/fields/fvPatchFields/basic/mixed/mixedFvPatchField.H index 3b1bc57c27af9d860b534a1f71e43beec7993ca3..7c475c1111e0dfeaa645c21c391ddd5f887df11f 100644 --- a/src/finiteVolume/fields/fvPatchFields/basic/mixed/mixedFvPatchField.H +++ b/src/finiteVolume/fields/fvPatchFields/basic/mixed/mixedFvPatchField.H @@ -35,19 +35,28 @@ Description The respective contributions from each is determined by a weight field: \f[ - x_p = w*A + (1-w)*(x_c + B/\Delta) + x_p = w x_p + (1-w) \left(x_c + \frac{\nabla_\perp x}{\Delta}\right) \f] where \vartable x_p | patch values x_c | patch internal cell values - w | weight field, \c valueFraction_ - A | fixed value, \c refValue_ - B | patch normal gradient, \c refGrad_ + w | weight field \Delta| inverse distance from face centre to internal cell centre + w | weighting (0-1) \endvartable + + \heading Patch usage + + \table + Property | Description | Required | Default value + valueFraction | weight field | yes | + refValue | fixed value | yes | + refGrad | patch normal gradient | yes | + \endtable + \note This condition is not usually applied directly; instead, use a derived mixed condition such as \c inletOutlet diff --git a/src/finiteVolume/fields/fvPatchFields/basic/zeroGradient/zeroGradientFvPatchField.H b/src/finiteVolume/fields/fvPatchFields/basic/zeroGradient/zeroGradientFvPatchField.H index 99038a12402e0d341dc6b5b092a68028a15c9ab5..cb238253f1d5829ca9b281415cd0eee6c498dbda 100644 --- a/src/finiteVolume/fields/fvPatchFields/basic/zeroGradient/zeroGradientFvPatchField.H +++ b/src/finiteVolume/fields/fvPatchFields/basic/zeroGradient/zeroGradientFvPatchField.H @@ -28,7 +28,7 @@ Group grpGenericBoundaryConditions Description - This boundary condition appies a zero-gradient condition from the patch + This boundary condition applies a zero-gradient condition from the patch internal field onto the patch faces. \heading Patch usage @@ -57,7 +57,7 @@ namespace Foam { /*---------------------------------------------------------------------------*\ - Class zeroGradientFvPatch Declaration + Class zeroGradientFvPatchField Declaration \*---------------------------------------------------------------------------*/ template<class Type> diff --git a/src/finiteVolume/fields/fvPatchFields/derived/fan/fanFvPatchField.H b/src/finiteVolume/fields/fvPatchFields/derived/fan/fanFvPatchField.H index be28d244488e522d51593c6f73921f1d39630dc2..9306796d799c61f073f56c93a0fcf49bb6da7a4e 100644 --- a/src/finiteVolume/fields/fvPatchFields/derived/fan/fanFvPatchField.H +++ b/src/finiteVolume/fields/fvPatchFields/derived/fan/fanFvPatchField.H @@ -28,7 +28,7 @@ Group grpCoupledBoundaryConditions Description - This boundary condition supplies a jump condition, using the \c cyclic + This boundary condition provides a jump condition, using the \c cyclic condition as a base. The jump is specified as a \c DataEntry type, to enable the use of, e.g. diff --git a/src/finiteVolume/fields/fvPatchFields/derived/fanPressure/fanPressureFvPatchScalarField.H b/src/finiteVolume/fields/fvPatchFields/derived/fanPressure/fanPressureFvPatchScalarField.H index a8124bb38b5bfe30e34b41d311ada6d31b0d06eb..022fb48ccb74e2ef3480e486414991afd69c04b8 100644 --- a/src/finiteVolume/fields/fvPatchFields/derived/fanPressure/fanPressureFvPatchScalarField.H +++ b/src/finiteVolume/fields/fvPatchFields/derived/fanPressure/fanPressureFvPatchScalarField.H @@ -31,31 +31,36 @@ Description This boundary condition can be applied to assign either a pressure inlet or outlet total pressure condition for a fan. - User specifies: - \li pressure drop vs volumetric flow rate table (fan curve) file name; - \li direction of normal flow through the fan, in or out; - \li total pressure of the environment. + \heading Patch usage + + \table + Property | Description | Required | Default value + fileName | fan curve file name | yes | + outOfBounds | out of bounds handling | yes | + direction | direction of flow through fan [in/out] | yes | + p0 | environmental total pressure | yes | + \endtable Example of the boundary condition specification: \verbatim inlet { type fanPressure; - fileName "fanCurve"; // Fan curve file name - outOfBounds clamp; // (error|warn|clamp|repeat) - direction in; // Direction of flow through fan - p0 uniform 0; // Environmental total pressure - value uniform 0; // Initial pressure + fileName "fanCurve"; + outOfBounds clamp; + direction in; + p0 uniform 0; + value uniform 0; } outlet { type fanPressure; - fileName "fanCurve"; // Fan curve file name - outOfBounds clamp; // (error|warn|clamp|repeat) - direction out; // Direction of flow through fan - p0 uniform 0; // Environmental total pressure - value uniform 0; // Initial pressure + fileName "fanCurve"; + outOfBounds clamp; + direction out; + p0 uniform 0; + value uniform 0; } \endverbatim diff --git a/src/finiteVolume/fields/fvPatchFields/derived/mappedFixedPushedInternalValue/mappedFixedPushedInternalValueFvPatchField.H b/src/finiteVolume/fields/fvPatchFields/derived/mappedFixedPushedInternalValue/mappedFixedPushedInternalValueFvPatchField.H index 18884da980f8d795cc6d5fe67c7980db61969b2f..04ec0c4ef7e7731f60ff22d6583e81ba23fd4785 100644 --- a/src/finiteVolume/fields/fvPatchFields/derived/mappedFixedPushedInternalValue/mappedFixedPushedInternalValueFvPatchField.H +++ b/src/finiteVolume/fields/fvPatchFields/derived/mappedFixedPushedInternalValue/mappedFixedPushedInternalValueFvPatchField.H @@ -44,7 +44,7 @@ Description \verbatim myPatch { - type mappedFixedInternalValue; + type mappedFixedPushedInternalValue; fieldName T; setAverage no; average 0; @@ -78,7 +78,7 @@ namespace Foam { /*---------------------------------------------------------------------------*\ - Class mappedFixedPushedInternalValueFvPatchField Declaration + Class mappedFixedPushedInternalValueFvPatchField Declaration \*---------------------------------------------------------------------------*/ template<class Type> diff --git a/src/finiteVolume/fields/fvPatchFields/derived/mappedVelocityFluxFixedValue/mappedVelocityFluxFixedValueFvPatchField.H b/src/finiteVolume/fields/fvPatchFields/derived/mappedVelocityFluxFixedValue/mappedVelocityFluxFixedValueFvPatchField.H index a160e2c292ba6fbd090aa4e50fc9da855db22182..8dac3f6f56d35b55786211418cd0c4279e5e6c01 100644 --- a/src/finiteVolume/fields/fvPatchFields/derived/mappedVelocityFluxFixedValue/mappedVelocityFluxFixedValueFvPatchField.H +++ b/src/finiteVolume/fields/fvPatchFields/derived/mappedVelocityFluxFixedValue/mappedVelocityFluxFixedValueFvPatchField.H @@ -42,7 +42,7 @@ Description \verbatim myPatch { - type mappedFixedInternalValue; + type mappedVelocityFlux; phi phi; value uniform 0; // place holder } diff --git a/src/finiteVolume/fields/fvPatchFields/derived/multiphaseFixedFluxPressure/multiphaseFixedFluxPressureFvPatchScalarField.H b/src/finiteVolume/fields/fvPatchFields/derived/multiphaseFixedFluxPressure/multiphaseFixedFluxPressureFvPatchScalarField.H index f71feaa1e34af2b08c54787355be477c1fee3c8c..c7f13d7506a3306eb355be381fa137a1b0974ca0 100644 --- a/src/finiteVolume/fields/fvPatchFields/derived/multiphaseFixedFluxPressure/multiphaseFixedFluxPressureFvPatchScalarField.H +++ b/src/finiteVolume/fields/fvPatchFields/derived/multiphaseFixedFluxPressure/multiphaseFixedFluxPressureFvPatchScalarField.H @@ -61,7 +61,7 @@ Description \verbatim myPatch { - type externalWallHeatFluxTemperature; + type multiphaseFixedFluxPressure; phiHbyA phiHbyA; phi phi; rho rho; diff --git a/src/finiteVolume/fields/fvPatchFields/derived/surfaceNormalFixedValue/surfaceNormalFixedValueFvPatchVectorField.H b/src/finiteVolume/fields/fvPatchFields/derived/surfaceNormalFixedValue/surfaceNormalFixedValueFvPatchVectorField.H index 5906f790a65e86a8414871f83ad0c8cbaf68f808..2b6242a42d384411f3f893230c4ecf17e15ab9c6 100644 --- a/src/finiteVolume/fields/fvPatchFields/derived/surfaceNormalFixedValue/surfaceNormalFixedValueFvPatchVectorField.H +++ b/src/finiteVolume/fields/fvPatchFields/derived/surfaceNormalFixedValue/surfaceNormalFixedValueFvPatchVectorField.H @@ -42,7 +42,7 @@ Description \verbatim myPatch { - type inletOutlet; + type surfaceNormalFixedValue; refValue -10; // 10 INTO the domain } \endverbatim diff --git a/src/finiteVolume/fields/fvPatchFields/derived/translatingWallVelocity/translatingWallVelocityFvPatchVectorField.H b/src/finiteVolume/fields/fvPatchFields/derived/translatingWallVelocity/translatingWallVelocityFvPatchVectorField.H index d32a1ba46084684e538a1e7507dae293191ea1d9..c6e9c70b88b0bd2a79706387096264e55b8d1fd5 100644 --- a/src/finiteVolume/fields/fvPatchFields/derived/translatingWallVelocity/translatingWallVelocityFvPatchVectorField.H +++ b/src/finiteVolume/fields/fvPatchFields/derived/translatingWallVelocity/translatingWallVelocityFvPatchVectorField.H @@ -42,7 +42,7 @@ Description \verbatim myPatch { - type rotatingWallVelocity; + type translatingWallVelocity; U (100 0 0); } \endverbatim diff --git a/src/thermophysicalModels/basic/derivedFvPatchFields/gradientEnergy/gradientEnergyFvPatchScalarField.C b/src/thermophysicalModels/basic/derivedFvPatchFields/gradientEnergy/gradientEnergyFvPatchScalarField.C index 42ed096ea60e297e05d85600c12843d4ebf7eeca..5f0e009d13c0731faf04bb3068ee972c80a825e9 100644 --- a/src/thermophysicalModels/basic/derivedFvPatchFields/gradientEnergy/gradientEnergyFvPatchScalarField.C +++ b/src/thermophysicalModels/basic/derivedFvPatchFields/gradientEnergy/gradientEnergyFvPatchScalarField.C @@ -97,10 +97,8 @@ void Foam::gradientEnergyFvPatchScalarField::updateCoeffs() return; } - const basicThermo& thermo = db().lookupObject<basicThermo> - ( - "thermophysicalProperties" - ); + const basicThermo& thermo = + db().lookupObject<basicThermo>("thermophysicalProperties"); const label patchi = patch().index(); diff --git a/src/thermophysicalModels/basic/derivedFvPatchFields/gradientEnergy/gradientEnergyFvPatchScalarField.H b/src/thermophysicalModels/basic/derivedFvPatchFields/gradientEnergy/gradientEnergyFvPatchScalarField.H index 2ee94c08cc6b60425d9ea35938df69750d5b21d8..a510d3ac82db8f6ad2b902e9e113282227665946 100644 --- a/src/thermophysicalModels/basic/derivedFvPatchFields/gradientEnergy/gradientEnergyFvPatchScalarField.H +++ b/src/thermophysicalModels/basic/derivedFvPatchFields/gradientEnergy/gradientEnergyFvPatchScalarField.H @@ -28,9 +28,25 @@ Group grpThermoBoundaryConditions Description - This boundary condition provides a gradient condition for internal energy + This boundary condition provides a gradient condition for internal energy, + where the gradient is calculated using: + + \f[ + \nabla(e_p) = \nabla_\perp C_p(p, T) + \frac{e_p - e_c}{\Delta} + \f] + + where + \vartable + e_p | energy at patch faces [J] + e_c | energy at patch internal cells [J] + p | pressure [bar] + T | temperature [K] + C_p | specific heat [J/kg/K] + \Delta | distance between patch face and internal cell centres [m] + \endvartable \heading Patch usage + Example of the boundary condition specification: \verbatim myPatch diff --git a/src/transportModels/twoPhaseInterfaceProperties/alphaContactAngle/alphaContactAngle/alphaContactAngleFvPatchScalarField.H b/src/transportModels/twoPhaseInterfaceProperties/alphaContactAngle/alphaContactAngle/alphaContactAngleFvPatchScalarField.H index de6086c906fd6f4c67ea4850d88bcc5aecab029d..6d1f65fe73f0c42feabea6a470be8093302116b8 100644 --- a/src/transportModels/twoPhaseInterfaceProperties/alphaContactAngle/alphaContactAngle/alphaContactAngleFvPatchScalarField.H +++ b/src/transportModels/twoPhaseInterfaceProperties/alphaContactAngle/alphaContactAngle/alphaContactAngleFvPatchScalarField.H @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2011 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -24,6 +24,9 @@ License Class Foam::alphaContactAngleFvPatchScalarField +Group + grpWallBoundaryConditions grpGenericBoundaryConditions + Description Abstract base class for alphaContactAngle boundary conditions. @@ -47,11 +50,12 @@ Description on p_rgh must set to guarantee that the flux is corrected to be zero at the wall e.g. - walls - { - type fixedFluxPressure; - adjoint no; - } + \verbatim + myPatch + { + type alphaContactAngle; + } + \endverbatim If "limit zeroGradient;" is used the pressure BCs can be left as before. @@ -72,7 +76,7 @@ namespace Foam { /*---------------------------------------------------------------------------*\ - Class alphaContactAngleFvPatch Declaration + Class alphaContactAngleFvPatchScalarField Declaration \*---------------------------------------------------------------------------*/ class alphaContactAngleFvPatchScalarField diff --git a/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/convectiveHeatTransfer/convectiveHeatTransferFvPatchScalarField.H b/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/convectiveHeatTransfer/convectiveHeatTransferFvPatchScalarField.H index aaaee1c6504498e6031092b2818f25b85f4f66c7..2bc06bb0e47751bf72d8887ca9ba011fb647ed34 100644 --- a/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/convectiveHeatTransfer/convectiveHeatTransferFvPatchScalarField.H +++ b/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/convectiveHeatTransfer/convectiveHeatTransferFvPatchScalarField.H @@ -61,7 +61,7 @@ Description \verbatim myPatch { - type turbulentInlet; + type convectiveHeatTransfer; L 0.1; } \endverbatim diff --git a/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/alphatWallFunctions/alphatJayatillekeWallFunction/alphatJayatillekeWallFunctionFvPatchScalarField.H b/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/alphatWallFunctions/alphatJayatillekeWallFunction/alphatJayatillekeWallFunctionFvPatchScalarField.H index 099819304627c9888248d00eb61427233995d7a3..60271aae48a4a17fc420737a61f4696468d2efd3 100644 --- a/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/alphatWallFunctions/alphatJayatillekeWallFunction/alphatJayatillekeWallFunctionFvPatchScalarField.H +++ b/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/alphatWallFunctions/alphatJayatillekeWallFunction/alphatJayatillekeWallFunctionFvPatchScalarField.H @@ -45,7 +45,7 @@ Description \verbatim myPatch { - type alphaSgsJayatillekeWallFunction; + type alphatJayatillekeWallFunction; Prt 0.85; kappa 0.41; E 9.8; diff --git a/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/mutWallFunctions/mutURoughWallFunction/mutURoughWallFunctionFvPatchScalarField.H b/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/mutWallFunctions/mutURoughWallFunction/mutURoughWallFunctionFvPatchScalarField.H index ee89a45f6cb57e188684ef4706511a7baeac575f..479cb18b4009dc3b940cd965a6ddcae6417726cd 100644 --- a/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/mutWallFunctions/mutURoughWallFunction/mutURoughWallFunctionFvPatchScalarField.H +++ b/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/mutWallFunctions/mutURoughWallFunction/mutURoughWallFunctionFvPatchScalarField.H @@ -37,14 +37,14 @@ Description Property | Description | Required | Default value roughnessHeight | roughness height | yes | roughnessConstant | roughness constanr | yes | - roughnessFactor | scaling factor | + roughnessFactor | scaling factor | yes | \endtable Example of the boundary condition specification: \verbatim myPatch { - type nutURoughWallFunction; + type mutURoughWallFunction; roughnessHeight 1e-5; roughnessConstant 0.5; roughnessFactor 1; diff --git a/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/mutWallFunctions/mutUSpaldingWallFunction/mutUSpaldingWallFunctionFvPatchScalarField.H b/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/mutWallFunctions/mutUSpaldingWallFunction/mutUSpaldingWallFunctionFvPatchScalarField.H index 4eeefb74ccf7f75f348492778aa132c9dd894d25..2501fa5b719427f210694eb2b417ab8b18ae851f 100644 --- a/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/mutWallFunctions/mutUSpaldingWallFunction/mutUSpaldingWallFunctionFvPatchScalarField.H +++ b/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/mutWallFunctions/mutUSpaldingWallFunction/mutUSpaldingWallFunctionFvPatchScalarField.H @@ -29,9 +29,22 @@ Group Description This boundary condition provides a turbulent viscosity condition when - using wall functions for rough walls, based on velocity, using Spaldings + using wall functions for rough walls, based on velocity, using Spalding's law to give a continuous nut profile to the wall (y+ = 0) + \f[ + y^+ = u^+ + \frac{1}{E} \left[exp(\kappa u^+) - 1 - \kappa u^+\, + - 0.5 (\kappa u^+)^2 - \frac{1}{6} (\kappa u^+)^3\right] + \f] + + where + \vartable + y^+ | non-dimensional position + u^+ | non-dimensional velocity + \kappa | Von Karman constant + \endvartable + + \heading Patch usage Example of the boundary condition specification: diff --git a/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/mutWallFunctions/mutkRoughWallFunction/mutkRoughWallFunctionFvPatchScalarField.H b/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/mutWallFunctions/mutkRoughWallFunction/mutkRoughWallFunctionFvPatchScalarField.H index 8daa978a221c1b9c21a0ee3d5beca04373d1b7d9..fdf5a0a7d6364919c2cea2b7604c9951953f2473 100644 --- a/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/mutWallFunctions/mutkRoughWallFunction/mutkRoughWallFunctionFvPatchScalarField.H +++ b/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/mutWallFunctions/mutkRoughWallFunction/mutkRoughWallFunctionFvPatchScalarField.H @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2011 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -49,7 +49,7 @@ Description \verbatim myPatch { - type nutkAtmRoughWallFunction; + type mutkRoughWallFunction; Ks uniform 0; Cs uniform 0.5; } diff --git a/src/turbulenceModels/derivedFvPatchFields/fixedShearStress/fixedShearStressFvPatchVectorField.H b/src/turbulenceModels/derivedFvPatchFields/fixedShearStress/fixedShearStressFvPatchVectorField.H index e8d9cc0376814ccdb6c5005442087cdf39993689..d1ad27977f37dfbad72012823b1158ee1d0d8522 100644 --- a/src/turbulenceModels/derivedFvPatchFields/fixedShearStress/fixedShearStressFvPatchVectorField.H +++ b/src/turbulenceModels/derivedFvPatchFields/fixedShearStress/fixedShearStressFvPatchVectorField.H @@ -25,7 +25,8 @@ License Class Foam::fixedShearStressFvPatchVectorField -Group grpWallBoundaryConditions +Group + grpWallBoundaryConditions Description Set a constant shear stress as tau0 = -nuEff dU/dn. diff --git a/src/turbulenceModels/derivedFvPatchFields/porousBafflePressure/porousBafflePressureFvPatchField.H b/src/turbulenceModels/derivedFvPatchFields/porousBafflePressure/porousBafflePressureFvPatchField.H index 8e0dccd0cc824874a7f81eec921840696519ad57..b45b7a330c2f01bbad0ca1c6eb34454fe6e7f582 100644 --- a/src/turbulenceModels/derivedFvPatchFields/porousBafflePressure/porousBafflePressureFvPatchField.H +++ b/src/turbulenceModels/derivedFvPatchFields/porousBafflePressure/porousBafflePressureFvPatchField.H @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2011 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -24,18 +24,57 @@ License Class Foam::porousBafflePressureFvPatchField -Description - Foam::porousBafflePressureFvPatchField - the porous baffle operates on a cyclic patch and introduce a jump on the p - field as follow: - - deltaP = -(I*mu*U + 0.5*D*rho*magSqr(U)*L) +Group + grpCoupledBoundaryConditions - where: - - I is the inertial coefficient - D is the darcy coeafficient - L is the porous media lenght in the flow direction +Description + This boundary condition provides a jump condition, using the \cyclic + condition as a base. + + The porous baffle introduces a pressure jump defined by: + + \f[ + \Delta p = -(I \mu U + 0.5 D \rho |U|^2 L) + \f] + + where + + \vartable + p | pressure [Pa] + \rho | density [kg/m3] + \mu | viscosity [Pa s] + I | inertial coefficient + D | Darcy coefficient + L | porous media length in the flow direction + \endvartable + + + \heading Patch usage + + \table + Property | Description | Required | Default value + patchType | underlying patch type should be \c cyclic| yes | + D | Darcy coefficient | yes | + I | inertial coefficient | yes | + L | porous media length in the flow direction | yes | + \endtable + + Example of the boundary condition specification: + \verbatim + myPatch + { + type porousBafflePressure; + patchType cyclic; + jump uniform 0; + D 1000000; + I 0.001; + L 0.1; + value uniform 0; + } + \endverbatim + + \note + The underlying \c patchType should be set to \c cyclic SourceFiles porousBafflePressureFvPatchField.C @@ -63,7 +102,7 @@ class porousBafflePressureFvPatchField { // Private data - //- Darcy pressure lost coefficient + //- Darcy pressure loss coefficient scalar D_; //- Inertia pressure lost coefficient diff --git a/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/atmBoundaryLayerInletEpsilon/atmBoundaryLayerInletEpsilonFvPatchScalarField.H b/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/atmBoundaryLayerInletEpsilon/atmBoundaryLayerInletEpsilonFvPatchScalarField.H index 10e6691933fa1a6fa50c728b801efa8aeb2bc0f9..f52dbac1736b3e12fd3f7aaad257fa43902c513c 100644 --- a/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/atmBoundaryLayerInletEpsilon/atmBoundaryLayerInletEpsilonFvPatchScalarField.H +++ b/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/atmBoundaryLayerInletEpsilon/atmBoundaryLayerInletEpsilonFvPatchScalarField.H @@ -34,7 +34,7 @@ Description \c ABLInletVelocity inlet velocity boundary condition. \f[ - \epsilon = \frac{(U^*)^3}{(K(z - z_g + z_0))} + \epsilon = \frac{(U^*)^3}{K(z - z_g + z_0)} \f] where @@ -49,7 +49,7 @@ Description and: \f[ - U^* = K \frac{U_{ref}}{ln(\frac{Z_{ref} + z_0}{z_0})} + U^* = K \frac{U_{ref}}{ln\left(\frac{Z_{ref} + z_0}{z_0}\right)} \f] where: @@ -74,7 +74,7 @@ Description \verbatim myPatch { - type atmBoundaryLayerInletVelocity; + type atmBoundaryLayerInletEpsilon; z 1.0; kappa 0.41; Uref 1.0; diff --git a/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/turbulentHeatFluxTemperature/turbulentHeatFluxTemperatureFvPatchScalarField.H b/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/turbulentHeatFluxTemperature/turbulentHeatFluxTemperatureFvPatchScalarField.H index 6e686dfbdd5414ca03b2000c5df913bd68b0c3e7..60d150deb0e32d40d42d4a43bd8f90292b4cd4a2 100644 --- a/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/turbulentHeatFluxTemperature/turbulentHeatFluxTemperatureFvPatchScalarField.H +++ b/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/turbulentHeatFluxTemperature/turbulentHeatFluxTemperatureFvPatchScalarField.H @@ -37,7 +37,7 @@ Description \table Property | Description | Required | Default value heatSource | heat source type: \c flux [W/m2] or \c power [W]| yes | - q | heat source value | yea | + q | heat source value | yes | alphaEff | turbulent thermal diffusivity field name | yes | \endtable diff --git a/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutURoughWallFunction/nutURoughWallFunctionFvPatchScalarField.H b/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutURoughWallFunction/nutURoughWallFunctionFvPatchScalarField.H index c0f8066e974195b58d03d17d013339bd245ec41d..881752d9f54e21fb0565311a0756957634d5cb68 100644 --- a/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutURoughWallFunction/nutURoughWallFunctionFvPatchScalarField.H +++ b/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutURoughWallFunction/nutURoughWallFunctionFvPatchScalarField.H @@ -37,7 +37,7 @@ Description Property | Description | Required | Default value roughnessHeight | roughness height | yes | roughnessConstant | roughness constanr | yes | - roughnessFactor | scaling factor | + roughnessFactor | scaling factor | yes | \endtable Example of the boundary condition specification: diff --git a/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutUSpaldingWallFunction/nutUSpaldingWallFunctionFvPatchScalarField.H b/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutUSpaldingWallFunction/nutUSpaldingWallFunctionFvPatchScalarField.H index ba1dde4514f0ad13b345dbe937c58634b2403086..b9f75ed1b17395581b2b981b9532bf4ba7bd7257 100644 --- a/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutUSpaldingWallFunction/nutUSpaldingWallFunctionFvPatchScalarField.H +++ b/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutUSpaldingWallFunction/nutUSpaldingWallFunctionFvPatchScalarField.H @@ -30,7 +30,20 @@ Group Description This boundary condition provides a turbulent kinematic viscosity condition when using wall functions for rough walls, based on velocity, using - Spaldings law to give a continuous nut profile to the wall (y+ = 0) + Spalding's law to give a continuous nut profile to the wall (y+ = 0) + + \f[ + y^+ = u^+ + \frac{1}{E} \left[exp(\kappa u^+) - 1 - \kappa u^+\, + - 0.5 (\kappa u^+)^2 - \frac{1}{6} (\kappa u^+)^3\right] + \f] + + where + \vartable + y^+ | non-dimensional position + u^+ | non-dimensional velocity + \kappa | Von Karman constant + \endvartable + \heading Patch usage diff --git a/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutkRoughWallFunction/nutkRoughWallFunctionFvPatchScalarField.H b/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutkRoughWallFunction/nutkRoughWallFunctionFvPatchScalarField.H index 8bbe8b16b8c6cabec9caa5e3a461f29251109d03..980f55f3b12e0c993ec541a00442aa6d10bac369 100644 --- a/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutkRoughWallFunction/nutkRoughWallFunctionFvPatchScalarField.H +++ b/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutkRoughWallFunction/nutkRoughWallFunctionFvPatchScalarField.H @@ -48,7 +48,7 @@ Description \verbatim myPatch { - type nutkAtmRoughWallFunction; + type nutkRoughWallFunction; Ks uniform 0; Cs uniform 0.5; }