Skip to content
Snippets Groups Projects
Commit 216f0a0c authored by Henry's avatar Henry
Browse files

pimpleDyMFoam: Corrected caching of rUA/rAU

parent bbecb950
No related branches found
No related tags found
No related merge requests found
...@@ -9,7 +9,7 @@ tmp<fvVectorMatrix> UEqn ...@@ -9,7 +9,7 @@ tmp<fvVectorMatrix> UEqn
UEqn().relax(); UEqn().relax();
volScalarField rUA = 1.0/UEqn().A(); volScalarField rAU = 1.0/UEqn().A();
if (momentumPredictor) if (momentumPredictor)
{ {
...@@ -17,6 +17,6 @@ if (momentumPredictor) ...@@ -17,6 +17,6 @@ if (momentumPredictor)
} }
else else
{ {
U = rUA*(UEqn().H() - fvc::grad(p)); U = rAU*(UEqn().H() - fvc::grad(p));
U.correctBoundaryConditions(); U.correctBoundaryConditions();
} }
U = rUA*UEqn().H(); U = rAU*UEqn().H();
if (nCorr <= 1) if (nCorr <= 1)
{ {
...@@ -6,7 +6,7 @@ if (nCorr <= 1) ...@@ -6,7 +6,7 @@ if (nCorr <= 1)
} }
phi = (fvc::interpolate(U) & mesh.Sf()) phi = (fvc::interpolate(U) & mesh.Sf())
+ fvc::ddtPhiCorr(rUA, U, phi); + fvc::ddtPhiCorr(rAU, U, phi);
adjustPhi(phi, U, p); adjustPhi(phi, U, p);
...@@ -16,7 +16,7 @@ for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++) ...@@ -16,7 +16,7 @@ for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
// Pressure corrector // Pressure corrector
fvScalarMatrix pEqn fvScalarMatrix pEqn
( (
fvm::laplacian(rUA, p) == fvc::div(phi) fvm::laplacian(rAU, p) == fvc::div(phi)
); );
pEqn.setReference(pRefCell, pRefValue); pEqn.setReference(pRefCell, pRefValue);
...@@ -47,5 +47,5 @@ for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++) ...@@ -47,5 +47,5 @@ for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
// Explicitly relax pressure for momentum corrector // Explicitly relax pressure for momentum corrector
p.relax(); p.relax();
U -= rUA*fvc::grad(p); U -= rAU*fvc::grad(p);
U.correctBoundaryConditions(); U.correctBoundaryConditions();
...@@ -27,7 +27,7 @@ ...@@ -27,7 +27,7 @@
mesh mesh
); );
# include "createPhi.H" #include "createPhi.H"
label pRefCell = 0; label pRefCell = 0;
......
...@@ -46,6 +46,7 @@ for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++) ...@@ -46,6 +46,7 @@ for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
#include "continuityErrs.H" #include "continuityErrs.H"
// Explicitly relax pressure for momentum corrector
p.relax(); p.relax();
// Make the fluxes relative to the mesh motion // Make the fluxes relative to the mesh motion
......
...@@ -20,47 +20,47 @@ internalField uniform (0 0 0); ...@@ -20,47 +20,47 @@ internalField uniform (0 0 0);
boundaryField boundaryField
{ {
movingWall movingWall
{ {
type fixedValue; type fixedValue;
value uniform (1 0 0); value uniform (1 0 0);
} }
farFieldMoving farFieldMoving
{ {
type fixedValue; type fixedValue;
value uniform (0 0 0); value uniform (0 0 0);
} }
fixedWall fixedWall
{ {
type fixedValue; type fixedValue;
value uniform (0 0 0); value uniform (0 0 0);
} }
axis axis
{ {
type symmetryPlane; type symmetryPlane;
} }
left left
{ {
type pressureInletOutletVelocity; type pressureInletOutletVelocity;
value uniform (0 0 0); value uniform (0 0 0);
} }
farField farField
{ {
type fixedValue; type fixedValue;
value uniform (0 0 0); value uniform (0 0 0);
} }
back back
{ {
type wedge; type wedge;
} }
front front
{ {
type wedge; type wedge;
} }
......
...@@ -20,27 +20,27 @@ internalField uniform 0; ...@@ -20,27 +20,27 @@ internalField uniform 0;
boundaryField boundaryField
{ {
movingWall movingWall
{ {
type zeroGradient; type zeroGradient;
} }
farFieldMoving farFieldMoving
{ {
type zeroGradient; type zeroGradient;
} }
fixedWall fixedWall
{ {
type zeroGradient; type zeroGradient;
} }
axis axis
{ {
type symmetryPlane; type symmetryPlane;
} }
left left
{ {
type totalPressure; type totalPressure;
p0 uniform 0; p0 uniform 0;
...@@ -52,17 +52,17 @@ boundaryField ...@@ -52,17 +52,17 @@ boundaryField
value uniform 0; value uniform 0;
} }
farField farField
{ {
type zeroGradient; type zeroGradient;
} }
back back
{ {
type wedge; type wedge;
} }
front front
{ {
type wedge; type wedge;
} }
......
...@@ -21,7 +21,7 @@ solvers ...@@ -21,7 +21,7 @@ solvers
{ {
solver PCG; solver PCG;
preconditioner DIC; preconditioner DIC;
tolerance 0.01; tolerance 0.02;
relTol 0; relTol 0;
} }
...@@ -30,7 +30,7 @@ solvers ...@@ -30,7 +30,7 @@ solvers
solver PCG; solver PCG;
preconditioner DIC; preconditioner DIC;
tolerance 1e-06; tolerance 1e-06;
relTol 0.05; relTol 0.01;
} }
pFinal pFinal
...@@ -66,10 +66,16 @@ solvers ...@@ -66,10 +66,16 @@ solvers
PIMPLE PIMPLE
{ {
correctPhi yes;
nOuterCorrectors 2; nOuterCorrectors 2;
nCorrectors 1; nCorrectors 1;
nNonOrthogonalCorrectors 0; nNonOrthogonalCorrectors 0;
} }
relaxationFactors
{
U 1;
UFinal 1;
}
// ************************************************************************* // // ************************************************************************* //
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment