exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

incompressible/icoFoam, solver + tutorials

Browse source
This commit is contained in:
Vuko Vukcevic 2016-04-05 11:23:41 +02:00
parent 5d2fb1e6f2
commit db568490b7
5 changed files with 26 additions and 31 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

37
applications/solvers/incompressible/icoFoam/icoFoam.C
View file

@ -26,6 +26,10 @@ Application
Description
Transient solver for incompressible, laminar flow of Newtonian fluids.
Consistent formulation without time-step and relaxation dependence by Jasak
Author
Hrvoje Jasak, Wikki Ltd. All rights reserved
\*---------------------------------------------------------------------------*/
@ -53,32 +57,36 @@ int main(int argc, char *argv[])
# include "readPISOControls.H"
# include "CourantNo.H"
fvVectorMatrix UEqn
// Convection-diffusion matrix
fvVectorMatrix HUEqn
(
fvm::ddt(U)
+ fvm::div(phi, U)
fvm::div(phi, U)
- fvm::laplacian(nu, U)
);
solve(UEqn == -fvc::grad(p));
// Time derivative matrix
fvVectorMatrix ddtUEqn(fvm::ddt(U));
solve(ddtUEqn + HUEqn == -fvc::grad(p));
// Prepare clean Ap without time derivative contribution
// HJ, 26/Oct/2015
volScalarField aU = HUEqn.A();
// --- PISO loop
for (int corr = 0; corr < nCorr; corr++)
{
volScalarField rUA = 1.0/UEqn.A();
U = rUA*UEqn.H();
phi = (fvc::interpolate(U) & mesh.Sf())
+ fvc::ddtPhiCorr(rUA, U, phi);
U = HUEqn.H()/aU;
phi = (fvc::interpolate(U) & mesh.Sf());
adjustPhi(phi, U, p);
for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
for (int nonOrth = 0; nonOrth <= nNonOrthCorr; nonOrth++)
{
fvScalarMatrix pEqn
(
fvm::laplacian(rUA, p) == fvc::div(phi)
fvm::laplacian(1/aU, p) == fvc::div(phi)
);
pEqn.setReference(pRefCell, pRefValue);
@ -92,7 +100,12 @@ int main(int argc, char *argv[])
# include "continuityErrs.H"
U -= rUA*fvc::grad(p);
// Note: cannot call H(U) here because the velocity is not complete
// HJ, 22/Jan/2016
U = 1.0/(aU + ddtUEqn.A())*
(
U*aU - fvc::grad(p) + ddtUEqn.H()
);
U.correctBoundaryConditions();
}

2
tutorials/incompressible/icoFoam/elbow/system/fvSchemes
View file

@ -28,7 +28,7 @@ gradSchemes
divSchemes
{
default none;
div(phi,U) Gauss limitedLinearV 1;
div(phi,U) Gauss linearUpwind leastSquares;
}
laplacianSchemes

6
tutorials/incompressible/icoFoam/mixingPlane/twoBlocksMixingPlaneMismatch_dirY_spanZ/0_orig/U
View file

@ -54,12 +54,6 @@ boundaryField
type mixingPlane;
}
top
{
type fixedValue;
value uniform (0 0 0);
}
frontAndBack
{
type empty;

6
tutorials/incompressible/icoFoam/mixingPlane/twoBlocksMixingPlane_dirY_spanZ/0_orig/U
View file

@ -54,12 +54,6 @@ boundaryField
type mixingPlane;
}
top
{
type fixedValue;
value uniform (0 0 0);
}
frontAndBack
{
type empty;

6
tutorials/incompressible/icoFoam/mixingPlane/twoBlocksMixingPlane_dirZ_spanY/0_orig/U
View file

@ -54,12 +54,6 @@ boundaryField
type mixingPlane;
}
top
{
type fixedValue;
value uniform (0 0 0);
}
frontAndBack
{
type empty;