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

Tutorial update

Browse source
This commit is contained in:
Hrvoje Jasak 2016-05-27 10:33:19 +01:00
parent 3522b3e62b
commit d589290b43
8 changed files with 55 additions and 57 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

4
tutorials/compressible/steadyUniversalMRFFoam/axialTurbineGgiJump/0/T
View file

@ -86,8 +86,8 @@ boundaryField
DTINLET
{
type ggiEnthalpyJump;
patchType ggi;
rotating false;
patchType ggi;
rotating false;
value uniform 305.66;
}
DTOUTLET

68
tutorials/compressible/steadyUniversalMRFFoam/axialTurbineGgiJump/0/i
View file

@ -23,15 +23,15 @@ boundaryField
{
GVINLET
{
type fixedEnthalpy;
value uniform 307494;
type fixedEnthalpy;
value uniform 307494;
}
GVOUTLET
{
type ggiEnthalpyJump;
patchType ggi;
rotating false;
value uniform 307494;
type ggiEnthalpyJump;
patchType ggi;
rotating false;
value uniform 307494;
}
GVCYCLIC
{
@ -48,91 +48,91 @@ boundaryField
{
type gradientEnthalpy;
gradient uniform 0;
value uniform 307494;
value uniform 307494;
}
GVSHROUD
{
type gradientEnthalpy;
gradient uniform 0;
value uniform 307494;
value uniform 307494;
}
RUINLET
{
type ggiEnthalpyJump;
patchType ggi;
rotating true;
value uniform 307494;
{
type ggiEnthalpyJump;
patchType ggi;
rotating true;
value uniform 307494;
}
RUOUTLET
{
type ggiEnthalpyJump;
patchType ggi;
rotating true;
value uniform 307494;
type ggiEnthalpyJump;
patchType ggi;
rotating true;
value uniform 307494;
}
RUCYCLIC1
{
type cyclicGgi;
value uniform 307494;
value uniform 307494;
}
RUCYCLIC2
{
type cyclicGgi;
value uniform 307494;
value uniform 307494;
}
RUBLADE
{
type gradientEnthalpy;
gradient uniform 307494;
value uniform 307494;
}
gradient uniform 307494;
value uniform 307494;
}
RUHUB
{
type gradientEnthalpy;
gradient uniform 307494;
type gradientEnthalpy;
gradient uniform 307494;
value uniform 307494;
}
RUSHROUD
{
type gradientEnthalpy;
gradient uniform 307494;
gradient uniform 307494;
value uniform 307494;
}
DTINLET
{
type ggiEnthalpyJump;
patchType ggi;
rotating false;
value uniform 305.66;
type ggiEnthalpyJump;
patchType ggi;
rotating false;
value uniform 305.66;
}
DTOUTLET
{
type gradientEnthalpy;
gradient uniform 0;
value uniform 307494;
value uniform 307494;
}
DTCYCLIC1
{
type cyclicGgi;
value uniform 307494;
value uniform 307494;
}
DTCYCLIC2
{
type cyclicGgi;
value uniform 307494;
value uniform 307494;
}
DTHUB
{
type gradientEnthalpy;
gradient uniform 0;
value uniform 307494;
value uniform 307494;
}
DTSHROUD
{
type gradientEnthalpy;
gradient uniform 0;
value uniform 307494;
value uniform 307494;
}
}

12
tutorials/compressible/steadyUniversalMRFFoam/axialTurbineGgiJump/0/p
View file

@ -27,8 +27,7 @@ boundaryField
GVOUTLET
{
type ggi;
patchType ggi;
value uniform 100000;
value uniform 100000;
}
GVCYCLIC
{
@ -49,14 +48,12 @@ boundaryField
RUINLET
{
type ggi;
patchType ggi;
value uniform 100000;
value uniform 100000;
}
RUOUTLET
{
type ggi;
patchType ggi;
value uniform 100000;
value uniform 100000;
}
RUCYCLIC1
{
@ -81,8 +78,7 @@ boundaryField
DTINLET
{
type ggi;
patchType ggi;
value uniform 100000;
value uniform 100000;
}
DTOUTLET
{

2
tutorials/incompressible/simpleFoam/motorBike/system/fvSchemes
View file

@ -30,7 +30,7 @@ divSchemes
div(phi,U) Gauss linearUpwindV Gauss linear;
div(phi,k) Gauss upwind;
div(phi,omega) Gauss upwind;
div((nuEff*dev(grad(U).T()))) Gauss linear;
div((nuEff*dev(T(grad(U))))) Gauss linear;
}
laplacianSchemes

2
tutorials/lagrangian/rhoPisoTwinParcelFoam/rhoPisoTwinParcelFoam/UEqn.H
View file

@ -11,7 +11,7 @@
UEqn.relax();
if (momentumPredictor)
if (pimple.momentumPredictor())
{
solve(UEqn == -fvc::grad(p));
}

10
tutorials/lagrangian/rhoPisoTwinParcelFoam/rhoPisoTwinParcelFoam/pEqn.H
View file

@ -3,7 +3,7 @@ rho = thermo.rho();
volScalarField rUA = 1.0/UEqn.A();
U = rUA*UEqn.H();
if (transonic)
if (pimple.transonic())
{
surfaceScalarField phid
(
@ -15,7 +15,7 @@ if (transonic)
)
);
for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
while (pimple.correctNonOrthogonal())
{
fvScalarMatrix pEqn
(
@ -26,7 +26,7 @@ if (transonic)
pEqn.solve();
if (nonOrth == nNonOrthCorr)
if (pimple.finalNonOrthogonalIter())
{
phi == pEqn.flux();
}
@ -41,7 +41,7 @@ else
+ fvc::ddtPhiCorr(rUA, rho, U, phi)
);
for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
while (pimple.correctNonOrthogonal())
{
fvScalarMatrix pEqn
(
@ -52,7 +52,7 @@ else
pEqn.solve();
if (nonOrth == nNonOrthCorr)
if (pimple.finalNonOrthogonalIter())
{
phi += pEqn.flux();
}

12
tutorials/lagrangian/rhoPisoTwinParcelFoam/rhoPisoTwinParcelFoam/rhoPisoTwinParcelFoam.C
View file

@ -32,6 +32,7 @@ Description
#include "fvCFD.H"
#include "basicPsiThermo.H"
#include "turbulenceModel.H"
#include "pimpleControl.H"
#include "basicThermoCloud.H"
#include "basicKinematicCloud.H"
@ -44,12 +45,14 @@ int main(int argc, char *argv[])
#include "createTime.H"
#include "createMesh.H"
pimpleControl pimple(mesh);
#include "readGravitationalAcceleration.H"
#include "createFields.H"
#include "createClouds.H"
#include "readPISOControls.H"
#include "initContinuityErrs.H"
#include "readTimeControls.H"
#include "createTimeControls.H"
#include "compressibleCourantNo.H"
#include "setInitialDeltaT.H"
@ -60,7 +63,6 @@ int main(int argc, char *argv[])
while (runTime.run())
{
#include "readTimeControls.H"
#include "readPISOControls.H"
#include "compressibleCourantNo.H"
#include "setDeltaT.H"
@ -78,12 +80,12 @@ int main(int argc, char *argv[])
#include "rhoEqn.H"
// --- PIMPLE loop
for (int ocorr=1; ocorr<=nOuterCorr; ocorr++)
while (pimple.loop())
{
#include "UEqn.H"
// --- PISO loop
for (int corr=1; corr<=nCorr; corr++)
while (pimple.correct())
{
#include "hsEqn.H"
#include "pEqn.H"

2
tutorials/lagrangian/rhoPisoTwinParcelFoam/simplifiedSiwek/system/fvSolution
View file

@ -74,7 +74,7 @@ solvers
}
}
PISO
PIMPLE
{
transonic yes;
nCorrectors 2;