research/OpenFOAM-4.x
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

advectionDiffusionPatchDistMethod: Only run the predictor method on first call

Browse source 
Improves efficiency for moving-mesh cases
This commit is contained in:
Henry 2015-01-09 22:05:32 +00:00
parent 28ec2c2a06
commit 4b4fce2c02
2 changed files with 12 additions and 3 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

12
src/finiteVolume/fvMesh/wallDist/patchDistMethods/advectionDiffusion/advectionDiffusionPatchDistMethod.C
View file

@ -68,7 +68,8 @@ Foam::patchDistMethods::advectionDiffusion::advectionDiffusion
),
epsilon_(coeffs_.lookupOrDefault<scalar>("epsilon", 0.1)),
tolerance_(coeffs_.lookupOrDefault<scalar>("tolerance", 1e-3)),
maxIter_(coeffs_.lookupOrDefault<int>("maxIter", 10))
maxIter_(coeffs_.lookupOrDefault<int>("maxIter", 10)),
predicted_(false)
{}
@ -86,7 +87,11 @@ bool Foam::patchDistMethods::advectionDiffusion::correct
volVectorField& n
)
{
pdmPredictor_->correct(y);
if (!predicted_)
{
pdmPredictor_->correct(y);
predicted_ = true;
}
volVectorField ny
(
@ -97,7 +102,7 @@ bool Foam::patchDistMethods::advectionDiffusion::correct
mesh_
),
mesh_,
dimensionedVector("nWall", dimless, vector::zero),
dimensionedVector("ny", dimless, vector::zero),
patchTypes<vector>(mesh_, patchIDs_)
);
@ -133,6 +138,7 @@ bool Foam::patchDistMethods::advectionDiffusion::correct
yEqn.relax();
initialResidual = yEqn.solve().initialResidual();
} while (initialResidual > tolerance_ && ++iter < maxIter_);
// Only calculate n if the field is defined

3
src/finiteVolume/fvMesh/wallDist/patchDistMethods/advectionDiffusion/advectionDiffusionPatchDistMethod.H
View file

@ -162,6 +162,9 @@ class advectionDiffusion
// to correct the distance-to-patch and normal-to-patch fields
int maxIter_;
//- Flag to indicate the predictor step has been executed
bool predicted_;
// Private Member Functions