Skip to content
Snippets Groups Projects
Commit 1161a119 authored by sergio's avatar sergio
Browse files

STY: Cleaning up tutorial settings

parent ae0b6ee7
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 75 additions and 552 deletions
applications/solvers/combustion/PDRFoam/createFieldRefs.H
+ 1
− 1
View file @ 1161a119
const volScalarField& psi = thermo.psi(); const volScalarField& psi = thermo.psi();
scalar pMin = thermo.getOrDefault<scalar>("pMin", VGREAT); scalar pMin = thermo.getOrDefault<scalar>("pMin", 0);
tutorials/combustion/PDRFoam/pipeLattice/constant/PDRProperties
+ 12
− 134
View file @ 1161a119
...@@ -15,19 +15,18 @@ FoamFile ...@@ -15,19 +15,18 @@ FoamFile
} }
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
smoothSt false;
StSmoothCoef 1000.0; StSmoothCoef 1000.0;
//No smooth if >100
schelkin schelkin
{ {
subGridSchelkin true; subGridSchelkin true;
uPrimeCoef 1.0; uPrimeCoef 1.0;
lCoef 0.336; lCoef 0.336;
maxSchFac 100.0; maxSchFac 100.0;
nrCoef 0.1; nrCoef 0.1;
nrExp 0.0; nrExp 0.0;
nrExp2 0.05; nrExp2 0.05;
} }
PDRDragModel basic; PDRDragModel basic;
...@@ -47,17 +46,8 @@ basicSchCoeffs ...@@ -47,17 +46,8 @@ basicSchCoeffs
Csk 0.0; Csk 0.0;
} }
XiModel transport; XiModel transportTwoEqs;
algebraicCoeffs
{
XiShapeCoef 1;
XpShapeCoef 1;
CpfiDot 0.0;
CpfiCross 0.0;
GEtaExp 0.0;
LOverCw 0.01;
}
transportTwoEqsCoeffs transportTwoEqsCoeffs
{ {
...@@ -69,51 +59,11 @@ transportTwoEqsCoeffs ...@@ -69,51 +59,11 @@ transportTwoEqsCoeffs
LOverCw 0.01; LOverCw 0.01;
} }
k3Coeffs
{
k3Obs 10.0;
k3Open 1.0;
}
transportThreeEqsCoeffs
{
$transportTwoEqsCoeffs;
$k3Coeffs;
}
transportFourEqsCoeffs
{
$transportTwoEqsCoeffs;
$k3Coeffs;
}
transportOneEqObsCoeffs
{
XiShapeCoef 1;
Cpfi 0.0;
GEtaExp 0.0;
}
transportXp
{
}
transportCoeffs
{
XiShapeCoef 1;
GEtaExp 0.28;
}
fixedCoeffs
{
}
/*---------------------------------------------------------------------------*\ /*---------------------------------------------------------------------------*\
XiEqModel : Model for XiEq XiEqModel : Model for XiEq
\*---------------------------------------------------------------------------*/ \*---------------------------------------------------------------------------*/
XiEqModel instability; XiEqModel instability2XiEq;
BLMcoeffs BLMcoeffs
{ {
...@@ -141,12 +91,9 @@ instability2XiEqCoeffs ...@@ -141,12 +91,9 @@ instability2XiEqCoeffs
BLMgMaXiEqCoeffs BLMgMaXiEqCoeffs
{ {
$schelkin; $schelkin;
gulderCoef 1.0; //this value is not usssed 1.0.
kaCoef 0.25; kaCoef 0.25;
lowK0 0.1; lowK0 0.1;
lowKg 0.0; lowKg 0.0;
gMaCoef 0.032; //not used
gMaCoef1 0.0; // not used
$BLMcoeffs; $BLMcoeffs;
} }
...@@ -161,45 +108,6 @@ instability2XiEqCoeffs ...@@ -161,45 +108,6 @@ instability2XiEqCoeffs
} }
} }
instabilityCoeffs
{
XiEqIn 2.5;
XiEqModel Gulder;
GulderCoeffs
{
uPrimeCoef 1.0;
subGridSchelkin true;
XiEqCoef 0.62;
}
SCOPEBlendCoeffs
{
XiEqModelL
{
XiEqModel Gulder;
GulderCoeffs
{
$schelkin;
XiEqCoef 0.62;
}
}
XiEqModelH
{
XiEqModel SCOPEXiEq;
SCOPEXiEqCoeffs
{
$schelkin;
XiEqCoef 1.6;
XiEqExp 0.33333;
}
}
}
}
/*---------------------------------------------------------------------------*\ /*---------------------------------------------------------------------------*\
XiGModel : Model for generation of Xi XiGModel : Model for generation of Xi
...@@ -207,25 +115,10 @@ instabilityCoeffs ...@@ -207,25 +115,10 @@ instabilityCoeffs
XiGModel instabilityG; XiGModel instabilityG;
instability2GCoeffs
{
lambdaIn 0.0001;
defaultCIn 10.0;
GInMult 5.0;
GInFade 4.0;
XiGModel KTS;
KTSCoeffs
{
GEtaCoef 0.28;
}
}
instabilityGCoeffs instabilityGCoeffs
{ {
lambdaIn 4.5e-3; lambdaIn lambdaIn [0 1 0 0 0 0 0] 4.5e-3;
GIn 1.917; GIn GIn [0 0 -1 0 0 0 0] 1.917;
XiGModel KTS; XiGModel KTS;
KTSCoeffs KTSCoeffs
...@@ -249,11 +142,6 @@ normBasicSubGridCoeffs ...@@ -249,11 +142,6 @@ normBasicSubGridCoeffs
Cxpe4 1.0; Cxpe4 1.0;
} }
basicSubGridCoeffs
{
}
/*---------------------------------------------------------------------------*\ /*---------------------------------------------------------------------------*\
XpGModel : Model for generation of Xp XpGModel : Model for generation of Xp
\*---------------------------------------------------------------------------*/ \*---------------------------------------------------------------------------*/
...@@ -273,14 +161,4 @@ normBasicSubGridGCoeffs ...@@ -273,14 +161,4 @@ normBasicSubGridGCoeffs
Cxpe3 400.0; Cxpe3 400.0;
Cxpe4 1.0; Cxpe4 1.0;
} }
basicSubGridGCoeffs
{
XpGModel KTS;
KTSCoeffs
{
GEtaCoef 0.28;
}
}
// ************************************************************************* // // ************************************************************************* //
tutorials/combustion/PDRFoam/pipeLattice/constant/thermophysicalProperties
+ 2
− 46
View file @ 1161a119
...@@ -15,15 +15,11 @@ FoamFile ...@@ -15,15 +15,11 @@ FoamFile
} }
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
//thermoType hhuMixtureThermo<inhomogeneousMixture<sutherlandTransport<specieThermo<janafThermo<perfectGas>>>>>;
thermoType thermoType
{ {
type heheuPsiThermo; type heheuPsiThermo;
mixture inhomogeneousMixture; mixture inhomogeneousMixture;
//mixture sprayHomogeneousMixtureDisp;
//mixture sprayInhomogeneousMixture;
//mixture inhomogeneousMixture;
transport sutherland; transport sutherland;
thermo janaf; thermo janaf;
equationOfState perfectGas; equationOfState perfectGas;
...@@ -31,26 +27,6 @@ thermoType ...@@ -31,26 +27,6 @@ thermoType
energy absoluteEnthalpy; energy absoluteEnthalpy;
} }
inertSpecie AIR;
liquids
{
//liquidComponents ( H2OPS );
H2OPS
{
defaultCoeffs yes;
}
water
{
defaultCoeffs yes;
}
}
solids
{
//solidComponents ( );
}
stoichiometricAirFuelMassRatio 17.1667; stoichiometricAirFuelMassRatio 17.1667;
fuel fuel
...@@ -119,27 +95,7 @@ burntProducts ...@@ -119,27 +95,7 @@ burntProducts
} }
} }
steam dpdt false;
{ pMin 1e3;
specie
{
nMoles 1;
molWeight 18.02;
}
thermodynamics
{
Tlow 200.00;
Thigh 5000.00;
Tcommon 1000.00;
highCpCoeffs ( 2.67215e+00 3.05629e-03 -8.73026e-07 1.20100e-10 -6.39162e-15 -29899.199219 6.862820);
lowCpCoeffs ( 3.38684e+00 3.47498e-03 -6.35470e-06 6.96858e-09 -2.50659e-12 -30208.099609 2.590230);
}
transport
{
As 1.6721e-06;
Ts 1.7067e+02;
}
}
// ************************************************************************* // // ************************************************************************* //
tutorials/combustion/PDRFoam/pipeLattice/constant/turbulenceProperties
+ 2
− 60
View file @ 1161a119
...@@ -22,18 +22,6 @@ RAS ...@@ -22,18 +22,6 @@ RAS
RASModel PDRkEpsilon; RASModel PDRkEpsilon;
PDRkEpsilonCoeffs PDRkEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 -0.33;//0; (Original : -0.33)
C4 0.1;
alphah 1;
alphak 1;
alphaEps 0.76923;
}
PDRkEpsilonFiltCoeffs
{ {
Cmu 0.09; Cmu 0.09;
C1 1.44; C1 1.44;
...@@ -49,54 +37,8 @@ RAS ...@@ -49,54 +37,8 @@ RAS
LOverMobs 0.0; LOverMobs 0.0;
maxLOverCellW 2.0; maxLOverCellW 2.0;
lCoef 0.35; lCoef 0.35;
resCoef 0.0; bMin 0.85;
maxResFac 1.0; flameFilteredG true;
filtRadMult 6.0;
bFilt 0.85;
flameFilteredG 1;
}
PDRRNGkEpsilonCoeffs
{
Cmu 0.0845;
C1 1.42;
C2 1.68;
C3 -0.33;
sigmak 0.71942;
sigmaEps 0.71942;
eta0 4.38;
beta 0.012;
LOverLobs 0.5;
LOverMobs 0.5;
maxLOverCellW 2.0;
lCoef 0.35;
C5 1.5;
C6 1.0;
alphah 1;
alphak 1;
alphaEps 0.76923;
}
PDRrealkEpsilonCoeffs
{
Cmu 0.09;
A0 4.0;
C2 1.9;
C3 -0.33;
sigmak 1;
sigmaEps 0.71942;
eta0 4.38;
beta 0.012;
LOverLobs 0.5;
LOverMobs 0.5;
maxLOverCellW 2.0;
lCoef 0.35 ;
C5 1.5;
C6 1.0;
alphah 1;
alphak 1;
alphaEps 0.76923;
} }
turbulence on; turbulence on;
......
tutorials/combustion/PDRFoam/pipeLattice/fuels/methane
+ 19
− 3
View file @ 1161a119
...@@ -4,8 +4,8 @@ name methane; ...@@ -4,8 +4,8 @@ name methane;
SCOPECoeffs SCOPECoeffs
{ {
lowerFlammabilityLimit 0.502449441017432; lowerFlamabilityLimit 0.502449441017432;
upperFlammabilityLimit 1.68468341988198; upperFlamabilityLimit 1.68468341988198;
lowerSuPolynomial lowerSuPolynomial
{ {
...@@ -72,7 +72,23 @@ SCOPECoeffs ...@@ -72,7 +72,23 @@ SCOPECoeffs
} }
Texp 2.0; Texp 2.0;
pexp -0.5; Texp0 1.5453;
Texp1 -0.4204;
Texp2 3.7503;
Texp3 1.782;
Texp4 0.0;
Texp5 0.0;
Texp6 0.0;
pexp -0.5;
pexp0 -0.314;
pexp1 0.608;
pexp2 0.0;
pexp3 0.0;
pexp4 0.0;
pexp5 0.0;
pexp6 0.0;
CIn 10.2; CIn 10.2;
} }
......
tutorials/combustion/PDRFoam/pipeLattice/system/controlDict
+ 39
− 2
View file @ 1161a119
...@@ -22,13 +22,13 @@ startTime 0; ...@@ -22,13 +22,13 @@ startTime 0;
stopAt endTime; stopAt endTime;
endTime 0.125; endTime 0.5;
deltaT 5e-06; deltaT 5e-06;
writeControl runTime; writeControl runTime;
writeInterval 0.001; writeInterval 0.01;
purgeWrite 0; purgeWrite 0;
...@@ -53,6 +53,43 @@ maxDeltaT 1; ...@@ -53,6 +53,43 @@ maxDeltaT 1;
functions functions
{ {
#include "probes" #include "probes"
normalisedGradP
{
libs (utilityFunctionObjects);
type coded;
name normalisedGradP;
writeControl none;
writeInterval -1;
executeControl timeStep;
executeInterval -1;
codeExecute
#{
const volScalarField& p = mesh().lookupObject<volScalarField>("p");
if (!mesh().foundObject<volScalarField>("normalisedGradP"))
{
volScalarField* normalisedGradPPtr =
new volScalarField("normalisedGradP", mag(fvc::grad(p)));
normalisedGradPPtr->store();
Info << "Found normalisedGradP"<< endl;
}
volScalarField& normalisedGradP =
mesh().lookupObjectRef<volScalarField>("normalisedGradP");
normalisedGradP = mag(fvc::grad(p));
Info << "Found not normalisedGradP"<< max(normalisedGradP) << endl;
#};
}
} }
// ************************************************************************* // // ************************************************************************* //
tutorials/combustion/PDRFoam/pipeLattice/system/fvSchemes
+ 0
− 17
View file @ 1161a119
...@@ -50,23 +50,6 @@ divSchemes ...@@ -50,23 +50,6 @@ divSchemes
div(phi,K) Gauss limitedLinear 1; div(phi,K) Gauss limitedLinear 1;
div(phi,ft_b_ha_hau_H2OPS_Y) Gauss multivariateSelection
{
ft limitedLinear01 1;
b limitedLinear01 1;
Xi limitedLinear 1;
Xp limitedLinear 1;
Ep limitedLinear 1;
h limitedLinear 1;
hu limitedLinear 1;
ha limitedLinear 1;
hau limitedLinear 1;
H2OPS limitedLinear 1;
water limitedLinear 1;
Y limitedLinear 1;
AIR limitedLinear 1;
};
div(phi,ft_b_h_hu) Gauss multivariateSelection div(phi,ft_b_h_hu) Gauss multivariateSelection
{ {
ft limitedLinear01 1; ft limitedLinear01 1;
......
tutorials/combustion/PDRFoam/pipeLattice/system/fvSchemesUp deleted 100644 → 0
+ 0
− 98
View file @ ae0b6ee7
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2012 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
ddtSchemes
{
default Euler;
}
gradSchemes
{
default Gauss linear;
grad(p) Gauss linear;
grad(b) fourth;
}
divSchemes
{
default Gauss upwind;
/*
div(phi,U) Gauss limitedLinearV 1;
div(phiU,p) Gauss limitedLinear 1;
div(phid,p) Gauss limitedLinear 1;
div(phi,k) Gauss limitedLinear 1;
div(phi,epsilon) Gauss limitedLinear 1;
div(phiXi,Xi) Gauss limitedLinear 1;
div(phiSt,b) Gauss limitedLinear01 1;
div(phi,K) Gauss limitedLinear 1;
*/
div(phi,ft_b_h_hu) Gauss multivariateSelection
{
ft upwind;//limitedLinear01 1;
b upwind;//limitedLinear01 1;
Xi upwind;//limitedLinear 1;
h upwind;//limitedLinear 1;
hu upwind;//limitedLinear 1;
Xp upwind;
Ep upwind;
};
div((Su*grad(b))) Gauss linear;
div((U+((Su*Xi)*grad(b)))) Gauss linear;
//div((muEff*dev2(T(grad(U))))) Gauss linear;
div(U) Gauss linear;
div((muEff*dev2(T(grad(Ubet))))) Gauss linear;
div(phi,Ubet) Gauss linear;
}
laplacianSchemes
{
default none;
laplacian(muEff,Ubet) Gauss linear limited 0.333;
laplacian(DkEff,k) Gauss linear limited 0.333;
laplacian(DepsilonEff,epsilon) Gauss linear limited 0.333;
laplacian((rho*inv((((1)*A(Ubet))+tDragDcu))),p) Gauss linear limited 0.333;
//laplacian((rho*inv((((1)*A(U))+((((0.5*rho)*CR)*mag(U))+((((Csu*(1))*betav)*muEff)*sqr(Aw)))))),p) Gauss linear limited 0.333;
laplacian(Db,b) Gauss linear limited 0.333;
laplacian(Db,ft) Gauss linear limited 0.333;
laplacian(Db,h) Gauss linear limited 0.333;
laplacian(Db,hu) Gauss linear limited 0.333;
laplacian(Db,Ep) Gauss linear limited 0.333;
laplacian(Db,GRxp) Gauss linear limited 0.333;
laplacian(Db,GRep) Gauss linear limited 0.333;
}
interpolationSchemes
{
default linear;
}
snGradSchemes
{
default limited 0.333;
}
fluxRequired
{
default no;
p;
}
// ************************************************************************* //
tutorials/combustion/PDRFoam/pipeLattice/system/fvSolutionTrans deleted 100644 → 0
+ 0
− 168
View file @ ae0b6ee7
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2012 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers
{
rho
{
solver PCG;
preconditioner DIC;
tolerance 1e-09;
relTol 0.0;
}
rhoFinal
{
$rho;
}
p
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-9;
relTol 0.01;
}
pFinal
{
$p;
relTol 0;
}
/* "(b|Xi|Xp|Ep|GRxp|GRep|ft|ha|hau|k|epsilon|H2OPS)"
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-09;
relTol 0.01;
}
"(b|Xi|Xp|Ep|GRxp|GRep|ft|ha|hau|k|epsilon|H2OPS)Final"
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-09;
relTol 0;
}
*/
"(b|Xi|Xp|Ep|GRxp|GRep|ft|h|hu|ha|hau|k|epsilon)"
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-09;
relTol 0.05;
}
"(b|Xi|Xp|Ep|GRxp|GRep|ft|h|hu|ha|hau|k|epsilon)Final"
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-09;
relTol 0;
}
}
PIMPLE
{
/* nCorrectors 3;
nOuterCorrectors 1;
nNonOrthogonalCorrectors 3;
momentumPredictor true;
transonic true;*/
nCorrectors 1;
nOuterCorrectors 10;
nNonOrthogonalCorrectors 3;
momentumPredictor true;
transonic true;
residualControl
{
b
{
tolerance 1e-5;
relTol 0;
}
p
{
tolerance 1e-5;
relTol 0;
}
Xi
{
tolerance 1e-5;
relTol 0;
}
Xp
{
tolerance 1e-5;
relTol 0;
}
Ep
{
tolerance 1e-5;
relTol 0;
}
ft
{
tolerance 1e-5;
relTol 0;
}
h
{
tolerance 1e-5;
relTol 0;
}
hu
{
tolerance 1e-5;
relTol 0;
}
k
{
tolerance 1e-5;
relTol 0;
}
epsilon
{
tolerance 1e-5;
relTol 0;
}
U
{
tolerance 1e-5;
relTol 0;
}
}
}
relaxationFactors
{
fields
{
}
equations
{
"U" 0.9;
"UFinal" 1.0;
"(b|Xi|Xp|Ep|GRxp|GRep|ft|ha|hau|k|epsilon|H2OPS)" 0.9;
"(b|Xi|Xp|Ep|GRxp|GRep|ft|ha|hau|k|epsilon|H2OPS)Final" 1.0;
}
}
// ************************************************************************* //
tutorials/combustion/PDRFoam/pipeLattice/system/topoSetDict
+ 0
− 23
View file @ 1161a119
...@@ -26,29 +26,6 @@ actions ...@@ -26,29 +26,6 @@ actions
source nearestToCell; source nearestToCell;
points ( $location ); points ( $location );
} }
// Optionally grow cellSet by attached?
/*
{
name corners;
type pointSet;
action new;
source cellToPoint;
set ignition;
option all;
}
{
name ignition;
type cellSet;
action add;
source pointToCell;
set corners;
option edge;
}
*/
); );
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment