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

ENH: overset: clean-up interpolation triggering. See #1041.

Browse Source
This commit is contained in:
mattijs
2018-10-17 08:54:34 +01:00
parent 4cce1d74d2
commit 03b5870c07
1 changed files with 56 additions and 46 deletions
Download Patch File Download Diff File Expand all files Collapse all files

102
src/overset/oversetPolyPatch/oversetFvPatchField.C
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2016-2017 OpenCFD Ltd.
\\ / A nd | Copyright (C) 2016-2018 OpenCFD Ltd.
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -113,60 +113,27 @@ void Foam::oversetFvPatchField<Type>::initEvaluate
<< fldName << endl;
}
}
else if
(
!fvSchemes.found("oversetInterpolation")
|| (
fvSchemes.found("oversetInterpolationRequired")
== fvSchemes.found("oversetInterpolationSuppressed")
)
)
else if (!fvSchemes.found("oversetInterpolation"))
{
IOWarningInFunction(fvSchemes)
<< "Missing required dictionary entries"
<< " 'oversetInterpolation' and 'oversetInterpolationRequired'"
<< " or 'oversetInterpolationSuppressed'"
<< "Missing required dictionary entry"
<< " 'oversetInterpolation'"
<< ". Skipping overset interpolation for field "
<< fldName << endl;
}
else if (fvSchemes.found("oversetInterpolationSuppressed"))
else if (fvSchemes.found("oversetInterpolationRequired"))
{
// Add the stencil suppression list
wordHashSet suppressed(Stencil::New(mesh).nonInterpolatedFields());
// Backwards compatibility mode: only interpolate what is
// explicitly mentioned
const dictionary* dictPtr
(
fvSchemes.findDict("oversetInterpolationSuppressed")
);
if (dictPtr)
if (fvSchemes.found("oversetInterpolationSuppressed"))
{
suppressed.insert(dictPtr->toc());
FatalIOErrorInFunction(fvSchemes)
<< "Cannot have both dictionary entry"
<< " 'oversetInterpolationSuppresed' and "
<< " 'oversetInterpolationRequired' for field "
<< fldName << exit(FatalIOError);
}
if (!suppressed.found(fldName))
{
if (debug)
{
Info<< "Interpolating non-suppressed field " << fldName
<< endl;
}
mesh.interpolate
(
const_cast<Field<Type>&>
(
this->primitiveField()
)
);
}
else if (debug)
{
Info<< "Skipping suppressed overset interpolation for field "
<< fldName << endl;
}
}
else
{
const dictionary& intDict = fvSchemes.subDict
(
"oversetInterpolationRequired"
@ -197,6 +164,49 @@ void Foam::oversetFvPatchField<Type>::initEvaluate
<< fldName << endl;
}
}
else
{
const dictionary* dictPtr
(
fvSchemes.subDictPtr("oversetInterpolationSuppressed")
);
const wordHashSet& suppress =
Stencil::New(mesh).nonInterpolatedFields();
bool skipInterpolate = suppress.found(fldName);
if (dictPtr)
{
skipInterpolate =
skipInterpolate
|| dictPtr->found(fldName);
}
if (skipInterpolate)
{
if (debug)
{
Info<< "Skipping suppressed overset interpolation"
<< " for field " << fldName << endl;
}
}
else
{
if (debug)
{
Info<< "Interpolating non-suppressed field " << fldName
<< endl;
}
mesh.interpolate
(
const_cast<Field<Type>&>
(
this->primitiveField()
)
);
}
}
}
}