Calculation of approximate distance to nearest patch for all cells and boundary by solving the Eikonal equation in advection form with diffusion smoothing. More...
Public Member Functions | |
TypeName ("advectionDiffusion") | |
Runtime type information. More... | |
advectionDiffusion (const dictionary &dict, const fvMesh &mesh, const labelHashSet &patchIDs) | |
Construct from coefficients dictionary, mesh. More... | |
advectionDiffusion (const advectionDiffusion &)=delete | |
Disallow default bitwise copy construction. More... | |
virtual bool | correct (volScalarField &y) |
Correct the given distance-to-patch field. More... | |
virtual bool | correct (volScalarField &y, volVectorField &n) |
Correct the given distance-to-patch and normal-to-patch fields. More... | |
void | operator= (const advectionDiffusion &)=delete |
Disallow default bitwise assignment. More... | |
Public Member Functions inherited from patchDistMethod | |
TypeName ("patchDistMethod") | |
Runtime type information. More... | |
declareRunTimeSelectionTable (autoPtr, patchDistMethod, dictionary,(const dictionary &dict, const fvMesh &mesh, const labelHashSet &patchIDs),(dict, mesh, patchIDs)) | |
patchDistMethod (const fvMesh &mesh, const labelHashSet &patchIDs) | |
Construct from mesh and patch ID set. More... | |
patchDistMethod (const patchDistMethod &)=delete | |
Disallow default bitwise copy construction. More... | |
virtual | ~patchDistMethod () |
Destructor. More... | |
const labelHashSet & | patchIndices () const |
Return the patchIDs. More... | |
virtual bool | movePoints () |
Update cached geometry when the mesh moves. More... | |
virtual void | topoChange (const polyTopoChangeMap &) |
Update cached topology and geometry when the mesh changes. More... | |
virtual void | mapMesh (const polyMeshMap &) |
Update from another mesh using the given map. More... | |
virtual void | distribute (const polyDistributionMap &) |
Redistribute or update using the given distribution map. More... | |
void | operator= (const patchDistMethod &)=delete |
Disallow default bitwise assignment. More... | |
template<class Type > | |
Foam::wordList | patchTypes (const fvMesh &mesh, const labelHashSet &patchIDs) |
Additional Inherited Members | |
Static Public Member Functions inherited from patchDistMethod | |
static autoPtr< patchDistMethod > | New (const dictionary &dict, const fvMesh &mesh, const labelHashSet &patchIDs) |
template<class Type > | |
static wordList | patchTypes (const fvMesh &mesh, const labelHashSet &patchIDs) |
Return the patch types for y and n. More... | |
Protected Attributes inherited from patchDistMethod | |
const fvMesh & | mesh_ |
Reference to the mesh. More... | |
const labelHashSet | patchIndices_ |
Set of patch IDs. More... | |
Calculation of approximate distance to nearest patch for all cells and boundary by solving the Eikonal equation in advection form with diffusion smoothing.
If the diffusion coefficient is set to 0 this method is exact in principle but the numerical schemes used are not rendering the scheme approximate, but more accurate than the Poisson method. Also many models relying on the distance to the wall benefit from this field being smooth and monotonic so the addition of diffusion smoothing improves both the convergence and stability of the solution of the Eikonal equation and the behavior of the models using the distance field generated. However, it is not clear what the optimum value for the diffusion coefficient epsilon should be; a default value of 0.1 is provided but higher values may be preferable under some circumstances.
Convergence is accelerated by first generating an approximate solution using one of the simpler methods, e.g. Poisson. The method used for this prediction step is specified in the 'advectionDiffusionCoeffs' sub-dictionary, see below.
References:
P.G. Tucker, C.L. Rumsey, R.E. Bartels, R.T. Biedron, "Transport equation based wall distance computations aimed at flows with time-dependent geometry", NASA/TM-2003-212680, December 2003.
Example of the wallDist specification in fvSchemes:
laplacianSchemes { . . laplacian(yPsi) Gauss linear corrected; laplacian(yWall) Gauss linear corrected; . . } wallDist { method advectionDiffusion; // Optional entry enabling the calculation // of the normal-to-wall field nRequired false; advectionDiffusionCoeffs { method Poisson; epsilon 0.1; tolerance 1e-3; maxIter 10; } }
Also the solver specification for yWall is required in fvSolution, e.g. for simple cases:
yPsi { solver PCG; preconditioner DIC; tolerance 1e-4; relTol 0; } yWall { solver PBiCG; preconditioner DILU; tolerance 1e-4; relTol 0; } relaxationFactors { equations { . . yWall 1; . . } } or for more complex cases: yPsi { solver GAMG; smoother GaussSeidel; cacheAgglomeration true; nCellsInCoarsestLevel 10; agglomerator faceAreaPair; mergeLevels 1; tolerance 1e-4; relTol 0; } yWall { solver GAMG; smoother symGaussSeidel; cacheAgglomeration true; nCellsInCoarsestLevel 10; agglomerator faceAreaPair; mergeLevels 1; tolerance 1e-4; relTol 0; } relaxationFactors { equations { . . yWall 0.7; . . } }
Definition at line 180 of file advectionDiffusionPatchDistMethod.H.
advectionDiffusion | ( | const dictionary & | dict, |
const fvMesh & | mesh, | ||
const labelHashSet & | patchIDs | ||
) |
Construct from coefficients dictionary, mesh.
and fixed-value patch set
Definition at line 51 of file advectionDiffusionPatchDistMethod.C.
|
delete |
Disallow default bitwise copy construction.
TypeName | ( | "advectionDiffusion" | ) |
Runtime type information.
|
virtual |
Correct the given distance-to-patch field.
Implements patchDistMethod.
Definition at line 78 of file advectionDiffusionPatchDistMethod.C.
References Foam::MULES::correct(), GeometricField< Type, PatchField, GeoMesh >::null(), and y.
|
virtual |
Correct the given distance-to-patch and normal-to-patch fields.
Implements patchDistMethod.
Definition at line 84 of file advectionDiffusionPatchDistMethod.C.
References GeometricField< Type, PatchField, GeoMesh >::boundaryFieldRef(), Foam::dimless, Foam::fvc::div(), Foam::fvm::div(), forAllConstIter, Foam::fvc::grad(), Foam::fvc::interpolate(), Foam::fvm::laplacian(), Foam::mag(), n, IOobject::NO_READ, IOobject::NO_WRITE, Foam::notNull(), patches, patchi, fvMatrix< Type >::relax(), fvMatrix< Type >::solve(), Foam::fvm::Sp(), y, and Foam::Zero.
|
delete |
Disallow default bitwise assignment.