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

ENH: better treatment of lInf in advective BC

Browse Source
This commit is contained in:
andy
2012-03-15 10:03:46 +00:00
parent 346ad0eba8
commit e887b572e9
1 changed files with 14 additions and 12 deletions
Download Patch File Download Diff File Expand all files Collapse all files

26
src/finiteVolume/fields/fvPatchFields/derived/advective/advectiveFvPatchField.C
View File

@ -45,7 +45,7 @@ Foam::advectiveFvPatchField<Type>::advectiveFvPatchField
phiName_("phi"),
rhoName_("rho"),
fieldInf_(pTraits<Type>::zero),
lInf_(0.0)
lInf_(-GREAT)
{
this->refValue() = pTraits<Type>::zero;
this->refGrad() = pTraits<Type>::zero;
@ -82,7 +82,7 @@ Foam::advectiveFvPatchField<Type>::advectiveFvPatchField
phiName_(dict.lookupOrDefault<word>("phi", "phi")),
rhoName_(dict.lookupOrDefault<word>("rho", "rho")),
fieldInf_(pTraits<Type>::zero),
lInf_(0.0)
lInf_(-GREAT)
{
if (dict.found("value"))
{
@ -110,9 +110,13 @@ Foam::advectiveFvPatchField<Type>::advectiveFvPatchField
(
"advectiveFvPatchField<Type>::"
"advectiveFvPatchField"
"(const fvPatch&, const Field<Type>&, const dictionary&)",
"("
"const fvPatch&, "
"const DimensionedField<Type, volMesh>&, "
"const dictionary&"
")",
dict
) << "unphysical lInf specified (lInf < 0)\n"
) << "unphysical lInf specified (lInf < 0)" << nl
<< " on patch " << this->patch().name()
<< " of field " << this->dimensionedInternalField().name()
<< " in file " << this->dimensionedInternalField().objectPath()
@ -217,7 +221,7 @@ void Foam::advectiveFvPatchField<Type>::updateCoeffs()
// Non-reflecting outflow boundary
// If lInf_ defined setup relaxation to the value fieldInf_.
if (lInf_ > SMALL)
if (lInf_ > 0)
{
// Calculate the field relaxation coefficient k (See notes)
const scalarField k(w*deltaT/lInf_);
@ -248,12 +252,10 @@ void Foam::advectiveFvPatchField<Type>::updateCoeffs()
}
else
{
FatalErrorIn
(
"advectiveFvPatchField<Type>::updateCoeffs()"
) << " Unsupported temporal differencing scheme : "
<< ddtScheme
<< "\n on patch " << this->patch().name()
FatalErrorIn("advectiveFvPatchField<Type>::updateCoeffs()")
<< " Unsupported temporal differencing scheme : "
<< ddtScheme << nl
<< " on patch " << this->patch().name()
<< " of field " << this->dimensionedInternalField().name()
<< " in file " << this->dimensionedInternalField().objectPath()
<< exit(FatalError);
@ -307,7 +309,7 @@ void Foam::advectiveFvPatchField<Type>::write(Ostream& os) const
this->template writeEntryIfDifferent<word>(os, "phi", "phi", phiName_);
this->template writeEntryIfDifferent<word>(os, "rho", "rho", rhoName_);
if (lInf_ > SMALL)
if (lInf_ > 0)
{
os.writeKeyword("fieldInf") << fieldInf_ << token::END_STATEMENT << nl;
os.writeKeyword("lInf") << lInf_ << token::END_STATEMENT << nl;