Skip to content
GitLab
Projects
Groups
Snippets
/
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in / Register
Toggle navigation
Menu
Open sidebar
Development
openfoam
Commits
4180b685
Commit
4180b685
authored
Jun 27, 2015
by
Henry Weller
Browse files
rhoPimpleFoam: Added run-time selectable LTS support replacing rhoLTSPimpleFoam
Select LTS via the ddtScheme: ddtSchemes { default localEuler rDeltaT; }
parent
60cffb2f
Changes
37
Hide whitespace changes
Inline
Side-by-side
applications/solvers/compressible/rhoPimpleFoam/rhoLTSPimpleFoam/Make/files
deleted
100644 → 0
View file @
60cffb2f
rhoLTSPimpleFoam.C
EXE = $(FOAM_APPBIN)/rhoLTSPimpleFoam
applications/solvers/compressible/rhoPimpleFoam/rhoLTSPimpleFoam/Make/options
deleted
100644 → 0
View file @
60cffb2f
EXE_INC = \
-I.. \
-I$(LIB_SRC)/transportModels/compressible/lnInclude \
-I$(LIB_SRC)/thermophysicalModels/basic/lnInclude \
-I$(LIB_SRC)/TurbulenceModels/turbulenceModels/lnInclude \
-I$(LIB_SRC)/TurbulenceModels/compressible/lnInclude \
-I$(LIB_SRC)/finiteVolume/cfdTools \
-I$(LIB_SRC)/finiteVolume/lnInclude \
-I$(LIB_SRC)/meshTools/lnInclude \
-I$(LIB_SRC)/sampling/lnInclude \
-I$(LIB_SRC)/fvOptions/lnInclude
EXE_LIBS = \
-lcompressibleTransportModels \
-lfluidThermophysicalModels \
-lspecie \
-lturbulenceModels \
-lcompressibleTurbulenceModels \
-lfiniteVolume \
-lmeshTools \
-lsampling \
-lfvOptions
applications/solvers/compressible/rhoPimpleFoam/rhoLTSPimpleFoam/createRDeltaT.H
deleted
100644 → 0
View file @
60cffb2f
scalar
maxDeltaT
(
pimple
.
dict
().
lookupOrDefault
<
scalar
>
(
"maxDeltaT"
,
GREAT
)
);
volScalarField
rDeltaT
(
IOobject
(
"rDeltaT"
,
runTime
.
timeName
(),
mesh
,
IOobject
::
NO_READ
,
IOobject
::
AUTO_WRITE
),
mesh
,
1
/
dimensionedScalar
(
"maxDeltaT"
,
dimTime
,
maxDeltaT
),
zeroGradientFvPatchScalarField
::
typeName
);
applications/solvers/compressible/rhoPimpleFoam/rhoLTSPimpleFoam/rhoLTSPimpleFoam.C
deleted
100644 → 0
View file @
60cffb2f
/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2015 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
OpenFOAM is free software: you can redistribute it and/or modify it
under the terms of the GNU General Public License as published by
the Free Software Foundation, either version 3 of the License, or
(at your option) any later version.
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
for more details.
You should have received a copy of the GNU General Public License
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
Application
rhoLTSPorousMRFPimpleFoam
Description
Transient solver for laminar or turbulent flow of compressible fluids
with support for run-time selectable finite volume options, e.g. MRF,
explicit porosity.
Uses the flexible PIMPLE (PISO-SIMPLE) solution for time-resolved and
pseudo-transient simulations with support for local time-stepping for
efficient steady-state solution.
\*---------------------------------------------------------------------------*/
#include
"fvCFD.H"
#include
"psiThermo.H"
#include
"turbulentFluidThermoModel.H"
#include
"fvIOoptionList.H"
#include
"fvcSmooth.H"
#include
"pimpleControl.H"
#include
"bound.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
int
main
(
int
argc
,
char
*
argv
[])
{
#include
"setRootCase.H"
#include
"createTime.H"
#include
"createMesh.H"
pimpleControl
pimple
(
mesh
);
#include
"createRDeltaT.H"
#include
"createFields.H"
#include
"createMRF.H"
#include
"createFvOptions.H"
#include
"initContinuityErrs.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Info
<<
"
\n
Starting time loop
\n
"
<<
endl
;
while
(
runTime
.
run
())
{
#include
"readTimeControls.H"
#include
"compressibleCourantNo.H"
#include
"setDeltaT.H"
runTime
++
;
Info
<<
"Time = "
<<
runTime
.
timeName
()
<<
nl
<<
endl
;
#include
"setRDeltaT.H"
if
(
pimple
.
nCorrPIMPLE
()
<=
1
)
{
#include
"rhoEqn.H"
}
// --- Pressure-velocity PIMPLE corrector loop
while
(
pimple
.
loop
())
{
turbulence
->
correct
();
#include
"UEqn.H"
#include
"EEqn.H"
// --- Pressure corrector loop
while
(
pimple
.
correct
())
{
#include
"pEqn.H"
}
}
runTime
.
write
();
Info
<<
"ExecutionTime = "
<<
runTime
.
elapsedCpuTime
()
<<
" s"
<<
" ClockTime = "
<<
runTime
.
elapsedClockTime
()
<<
" s"
<<
nl
<<
endl
;
}
Info
<<
"End
\n
"
<<
endl
;
return
0
;
}
// ************************************************************************* //
applications/solvers/compressible/rhoPimpleFoam/rhoPimpleDyMFoam/rhoPimpleDyMFoam.C
View file @
4180b685
...
...
@@ -44,6 +44,8 @@ Description
#include
"pimpleControl.H"
#include
"CorrectPhi.H"
#include
"fvIOoptionList.H"
#include
"localEulerDdtScheme.H"
#include
"fvcSmooth.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
...
...
@@ -52,17 +54,23 @@ int main(int argc, char *argv[])
#include
"setRootCase.H"
#include
"createTime.H"
#include
"createDynamicFvMesh.H"
#include
"initContinuityErrs.H"
pimpleControl
pimple
(
mesh
);
#include
"readControls.H"
#include
"createRDeltaT.H"
#include
"initContinuityErrs.H"
#include
"createFields.H"
#include
"createMRF.H"
#include
"createFvOptions.H"
#include
"createRhoUf.H"
#include
"CourantNo.H"
#include
"setInitialDeltaT.H"
if
(
!
LTS
)
{
#include
"readTimeControls.H"
#include
"CourantNo.H"
#include
"setInitialDeltaT.H"
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
...
...
@@ -82,8 +90,15 @@ int main(int argc, char *argv[])
fvc
::
div
(
fvc
::
absolute
(
phi
,
rho
,
U
))
);
#include
"compressibleCourantNo.H"
#include
"setDeltaT.H"
if
(
LTS
)
{
#include
"setRDeltaT.H"
}
else
{
#include
"compressibleCourantNo.H"
#include
"setDeltaT.H"
}
runTime
++
;
...
...
applications/solvers/compressible/rhoPimpleFoam/rhoPimpleFoam.C
View file @
4180b685
...
...
@@ -39,6 +39,8 @@ Description
#include
"bound.H"
#include
"pimpleControl.H"
#include
"fvIOoptionList.H"
#include
"localEulerDdtScheme.H"
#include
"fvcSmooth.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
...
...
@@ -50,10 +52,18 @@ int main(int argc, char *argv[])
pimpleControl
pimple
(
mesh
);
#include
"createRDeltaT.H"
#include
"initContinuityErrs.H"
#include
"createFields.H"
#include
"createMRF.H"
#include
"createFvOptions.H"
#include
"initContinuityErrs.H"
if
(
!
LTS
)
{
#include
"readTimeControls.H"
#include
"CourantNo.H"
#include
"setInitialDeltaT.H"
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
...
...
@@ -62,8 +72,16 @@ int main(int argc, char *argv[])
while
(
runTime
.
run
())
{
#include
"readTimeControls.H"
#include
"compressibleCourantNo.H"
#include
"setDeltaT.H"
if
(
LTS
)
{
#include
"setRDeltaT.H"
}
else
{
#include
"compressibleCourantNo.H"
#include
"setDeltaT.H"
}
runTime
++
;
...
...
applications/solvers/compressible/rhoPimpleFoam/
rhoLTSPimpleFoam/
setRDeltaT.H
→
applications/solvers/compressible/rhoPimpleFoam/setRDeltaT.H
View file @
4180b685
{
volScalarField
&
rDeltaT
=
trDeltaT
();
const
dictionary
&
pimpleDict
=
pimple
.
dict
();
scalar
maxCo
...
...
tutorials/compressible/rhoLTSPimpleFoam/angledDuct/0/T
deleted
100644 → 0
View file @
60cffb2f
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: dev |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 1 0 0 0];
internalField uniform 293;
boundaryField
{
front
{
type zeroGradient;
}
back
{
type zeroGradient;
}
walls
{
type zeroGradient;
}
porosityWall
{
type zeroGradient;
}
inlet
{
type fixedValue;
value $internalField;
}
outlet
{
type inletOutlet;
value $internalField;
inletValue $internalField;
}
}
// ************************************************************************* //
tutorials/compressible/rhoLTSPimpleFoam/angledDuct/0/U
deleted
100644 → 0
View file @
60cffb2f
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: dev |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (47.8 0 0);
boundaryField
{
front
{
type fixedValue;
value uniform (0 0 0);
}
back
{
type fixedValue;
value uniform (0 0 0);
}
walls
{
type fixedValue;
value uniform (0 0 0);
}
porosityWall
{
type slip;
value uniform (0 0 0);
}
inlet
{
type flowRateInletVelocity;
massFlowRate constant 0.1;
value uniform (0 0 0);
}
outlet
{
type fluxCorrectedVelocity; //inletOutlet;
value uniform (0 0 0);
inletValue uniform (0 0 0);
}
}
// ************************************************************************* //
tutorials/compressible/rhoLTSPimpleFoam/angledDuct/0/alphat
deleted
100644 → 0
View file @
60cffb2f
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: dev |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object alphat;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -1 -1 0 0 0 0];
internalField uniform 0;
boundaryField
{
front
{
type compressible::alphatWallFunction;
value uniform 0;
}
back
{
type compressible::alphatWallFunction;
value uniform 0;
}
walls
{
type compressible::alphatWallFunction;
value uniform 0;
}
porosityWall
{
type compressible::alphatWallFunction;
value uniform 0;
}
inlet
{
type calculated;
value uniform 0;
}
outlet
{
type calculated;
value uniform 0;
}
}
// ************************************************************************* //
tutorials/compressible/rhoLTSPimpleFoam/angledDuct/0/epsilon
deleted
100644 → 0
View file @
60cffb2f
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: dev |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object epsilon;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -3 0 0 0 0];
internalField uniform 200;
boundaryField
{
front
{
type epsilonWallFunction;
value uniform 200;
}
back
{
type epsilonWallFunction;
value uniform 200;
}
walls
{
type epsilonWallFunction;
value uniform 200;
}
porosityWall
{
type epsilonWallFunction;
value uniform 200;
}
inlet
{
type turbulentMixingLengthDissipationRateInlet;
mixingLength 0.005;
value uniform 200;
}
outlet
{
type inletOutlet;
inletValue uniform 200;
value uniform 200;
}
}
// ************************************************************************* //
tutorials/compressible/rhoLTSPimpleFoam/angledDuct/0/k
deleted
100644 → 0
View file @
60cffb2f
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: dev |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -2 0 0 0 0];
internalField uniform 1;
boundaryField
{
front
{
type kqRWallFunction;
value uniform 1;
}
back
{
type kqRWallFunction;
value uniform 1;
}
walls
{
type kqRWallFunction;
value uniform 1;
}
porosityWall
{
type kqRWallFunction;
value uniform 1;
}
inlet
{
type turbulentIntensityKineticEnergyInlet;
intensity 0.05;
value uniform 1;
}
outlet
{
type inletOutlet;
inletValue uniform 1;
value uniform 1;
}
}
// ************************************************************************* //
tutorials/compressible/rhoLTSPimpleFoam/angledDuct/0/nut
deleted
100644 → 0
View file @
60cffb2f
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: dev |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -1 0 0 0 0];
internalField uniform 0;
boundaryField
{
front
{
type nutkWallFunction;
value uniform 0;
}
back
{
type nutkWallFunction;
value uniform 0;
}
walls
{
type nutkWallFunction;
value uniform 0;
}
porosityWall
{
type nutkWallFunction;
value uniform 0;
}
inlet
{
type calculated;
value uniform 0;
}
outlet
{