Commit d0c39b0e authored by sergio's avatar sergio
Browse files

ENH: Changing non-interpolation field in versetFvPatchField

for correct run time reading.
Setting up twoSimpleRotors tutorial for smooth pressure
fluctuations
parent ebdb2278
......@@ -161,6 +161,12 @@ const Foam::wordHashSet& Foam::cellCellStencil::nonInterpolatedFields() const
}
Foam::wordHashSet& Foam::cellCellStencil::nonInterpolatedFields()
{
return nonInterpolatedFields_;
}
bool Foam::cellCellStencil::localStencil(const labelUList& slots) const
{
forAll(slots, i)
......
......@@ -184,6 +184,9 @@ public:
// should not be interpolated
virtual const wordHashSet& nonInterpolatedFields() const;
//- Return non-const non-interpolating fields
virtual wordHashSet& nonInterpolatedFields();
//- Helper: is stencil fully local
bool localStencil(const labelUList&) const;
......
......@@ -25,6 +25,7 @@ License
#include "volFields.H"
#include "cellCellStencil.H"
#include "cellCellStencilObject.H"
#include "dynamicOversetFvMesh.H"
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
......@@ -130,11 +131,20 @@ void Foam::oversetFvPatchField<Type>::initEvaluate
}
else if (fvSchemes.found("oversetInterpolationSuppressed"))
{
const dictionary& intDict = fvSchemes.subDict
// Add the stencil suppression list
wordHashSet suppressed(Stencil::New(mesh).nonInterpolatedFields());
const dictionary* dictPtr
(
"oversetInterpolationSuppressed"
fvSchemes.subDictPtr("oversetInterpolationSuppressed")
);
if (!intDict.found(fldName))
if (dictPtr)
{
suppressed.insert(dictPtr->sortedToc());
}
if (!suppressed.found(fldName))
{
if (debug)
{
......
......@@ -14,7 +14,7 @@ FoamFile
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
simulationType RAS;
simulationType laminar;
RAS
{
......
......@@ -37,7 +37,7 @@ vertices
( 0.85 0.65 1)
( 0.15 0.65 1)
// movingZone2
// // movingZone2
( 0.75 0.15 0)
( 1.05 0.15 0)
( 1.05 0.85 0)
......@@ -46,6 +46,15 @@ vertices
( 1.05 0.15 1)
( 1.05 0.85 1)
( 0.75 0.85 1)
// ( 0.95 0.15 0)
// ( 1.25 0.15 0)
// ( 1.25 0.85 0)
// ( 0.95 0.85 0)
// ( 0.95 0.15 1)
// ( 1.25 0.15 1)
// ( 1.25 0.85 1)
// ( 0.95 0.85 1)
);
blocks
......
......@@ -20,23 +20,24 @@ DebugSwitches
{
overset 0;
dynamicOversetFvMesh 0;
cellVolumeWeight 0;
}
application overPimpleDyMFoam;
startFrom startTime;
startFrom latestTime;
startTime 0;
stopAt endTime;
endTime 0.08;
endTime 0.8;
deltaT 2e-5;
writeControl adjustableRunTime;
writeControl timeStep;
writeInterval 0.005;
writeInterval 30;
purgeWrite 0;
......@@ -52,10 +53,9 @@ timePrecision 6;
runTimeModifiable true;
adjustTimeStep no;
maxCo 1;
adjustTimeStep true;
maxCo 1;
functions
{
......@@ -80,7 +80,7 @@ functions
probeLocations
(
(0.015 0.005 0.005) // at inlet
(0.015 0.005 0.005)
);
}
......
......@@ -29,8 +29,8 @@ divSchemes
default none;
div(phi,U) Gauss upwind;
div(phi,epsilon) Gauss upwind;
div(phi,k) Gauss upwind;
div(phi,epsilon) Gauss limitedLinear 1;
div(phi,k) Gauss limitedLinear 1;
div((nuEff*dev2(T(grad(U))))) Gauss linear;
}
......
......@@ -30,7 +30,7 @@ solvers
{
solver PBiCGStab;
preconditioner DILU;
tolerance 1e-9;
tolerance 1e-11;
relTol 0.01;
}
......@@ -71,10 +71,10 @@ solvers
PIMPLE
{
momentumPredictor true;
momentumPredictor false;
correctPhi no;
nOuterCorrectors 1;
nCorrectors 5;
nCorrectors 3;
nNonOrthogonalCorrectors 0;
ddtCorr true;
......
......@@ -100,8 +100,9 @@ actions
{
boxes
(
(0.0025 0.0045 -100)(0.0075 0.0055 100)
(0.0085 0.0025 -100)(0.0095 0.0075 100)
(0.0025 0.0045 -100)(0.0075 0.0055 100)
(0.0085 0.0025 -100)(0.0095 0.0075 100)
// (0.0104 0.0025 -100)(0.0115 0.0075 100)
);
}
}
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment