Commit 4363dac7 authored by henry's avatar henry
Browse files

Further updates of the XiFoam tutorials.

parent 589a342e
......@@ -42,7 +42,7 @@ boundaryField
lowerWall
{
type fixedValue;
value uniform 293;
value uniform 570;
}
frontAndBack
......
......@@ -42,7 +42,7 @@ boundaryField
lowerWall
{
type fixedValue;
value uniform 293;
value uniform 570;
}
frontAndBack
......
......@@ -32,15 +32,13 @@ boundaryField
upperWall
{
type zeroGradient;
//type compressible::alphaSgsWallFunction;
type compressible::alphaSgsJayatillekeWallFunction;
value uniform 0;
}
lowerWall
{
type zeroGradient;
//type compressible::alphaSgsWallFunction;
type compressible::alphaSgsJayatillekeWallFunction;
value uniform 0;
}
......
......@@ -32,15 +32,13 @@ boundaryField
upperWall
{
type zeroGradient;
//type compressible::muSgsWallFunction;
type compressible::muSgsWallFunction;
value uniform 0;
}
lowerWall
{
type zeroGradient;
//type compressible::muSgsWallFunction;
type compressible::muSgsWallFunction;
value uniform 0;
}
......
......@@ -16,7 +16,7 @@ FoamFile
dimensions [1 -1 -2 0 0 0 0];
internalField uniform 100000;
internalField uniform 1e5;
boundaryField
{
......@@ -27,8 +27,15 @@ boundaryField
outlet
{
type fixedValue;
value uniform 100000;
type waveTransmissive;
field p;
phi phi;
rho rho;
psi psi;
gamma 1.3;
fieldInf 1e5;
lInf 0.3;
value uniform 1e5;
}
upperWall
......
......@@ -123,6 +123,7 @@ kappa 0.4187;
wallFunctionCoeffs
{
E 9;
Prt 0.85;
}
......
......@@ -21,7 +21,7 @@ fuel Propane;
Su Su [0 1 -1 0 0 0 0] 0.135;
SuModel unstrained;
SuModel transport;
equivalenceRatio equivalenceRatio [0 0 0 0 0 0 0] 0.6;
......@@ -68,7 +68,7 @@ GuldersCoeffs
}
}
ignite yes;
ignite yes;
ignitionSites
(
......@@ -76,8 +76,8 @@ ignitionSites
location (0.005 -0.02 0);
diameter 0.003;
start 0;
duration 0.1;
strength 40;
duration 0.05;
strength 20;
}
);
......
......@@ -16,55 +16,55 @@ FoamFile
convertToMeters 0.001;
vertices
vertices
(
(-20.6 0 -0.5)
(-20.6 3 -0.5)
(-20.6 12.7 -0.5)
(-20.6 25.4 -0.5)
(0 -25.4 -0.5)
(0 -5 -0.5)
(0 0 -0.5)
(0 3 -0.5)
(0 12.7 -0.5)
(0 25.4 -0.5)
(206 -25.4 -0.5)
(206 -8.5 -0.5)
(206 0 -0.5)
(206 6.5 -0.5)
(206 17 -0.5)
(206 25.4 -0.5)
(290 -16.6 -0.5)
(290 -6.3 -0.5)
(290 0 -0.5)
(290 4.5 -0.5)
(290 11 -0.5)
(290 16.6 -0.5)
(-20.6 0 0.5)
(-20.6 3 0.5)
(-20.6 12.7 0.5)
(-20.6 25.4 0.5)
(0 -25.4 0.5)
(0 -5 0.5)
(0 0 0.5)
(0 3 0.5)
(0 12.7 0.5)
(0 25.4 0.5)
(206 -25.4 0.5)
(206 -8.5 0.5)
(206 0 0.5)
(206 6.5 0.5)
(206 17 0.5)
(206 25.4 0.5)
(290 -16.6 0.5)
(290 -6.3 0.5)
(290 0 0.5)
(290 4.5 0.5)
(290 11 0.5)
(290 16.6 0.5)
(-20.6 0 0)
(-20.6 3 0)
(-20.6 12.7 0)
(-20.6 25.4 0)
(0 -25.4 0)
(0 -5 0)
(0 0 0)
(0 3 0)
(0 12.7 0)
(0 25.4 0)
(206 -25.4 0)
(206 -8.5 0)
(206 0 0)
(206 6.5 0)
(206 17 0)
(206 25.4 0)
(290 -16.6 0)
(290 -6.3 0)
(290 0 0)
(290 4.5 0)
(290 11 0)
(290 16.6 0)
(-20.6 0 38.1)
(-20.6 3 38.1)
(-20.6 12.7 38.1)
(-20.6 25.4 38.1)
(0 -25.4 38.1)
(0 -5 38.1)
(0 0 38.1)
(0 3 38.1)
(0 12.7 38.1)
(0 25.4 38.1)
(206 -25.4 38.1)
(206 -8.5 38.1)
(206 0 38.1)
(206 6.5 38.1)
(206 17 38.1)
(206 25.4 38.1)
(290 -16.6 38.1)
(290 -6.3 38.1)
(290 0 38.1)
(290 4.5 38.1)
(290 11 38.1)
(290 16.6 38.1)
);
blocks
blocks
(
hex (0 6 7 1 22 28 29 23) (18 7 1) simpleGrading (0.5 1.8 1)
hex (1 7 8 2 23 29 30 24) (18 10 1) simpleGrading (0.5 4 1)
......@@ -81,19 +81,19 @@ blocks
hex (14 20 21 15 36 42 43 37) (25 13 1) simpleGrading (2.5 0.25 1)
);
edges
edges
(
);
patches
patches
(
patch inlet
patch inlet
(
(0 22 23 1)
(1 23 24 2)
(2 24 25 3)
)
patch outlet
patch outlet
(
(16 17 39 38)
(17 18 40 39)
......@@ -101,13 +101,13 @@ patches
(19 20 42 41)
(20 21 43 42)
)
wall upperWall
wall upperWall
(
(3 25 31 9)
(9 31 37 15)
(15 37 43 21)
)
wall lowerWall
wall lowerWall
(
(0 6 28 22)
(6 5 27 28)
......@@ -115,7 +115,7 @@ patches
(4 10 32 26)
(10 16 38 32)
)
empty frontAndBack
empty frontAndBack
(
(22 28 29 23)
(23 29 30 24)
......
......@@ -15,19 +15,44 @@ FoamFile
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
thermoType hhuMixtureThermo<homogeneousMixture<sutherlandTransport<specieThermo<janafThermo<perfectGas>>>>>;
stoichiometricAirFuelMassRatio stoichiometricAirFuelMassRatio [ 0 0 0 0 0 0 0 ] 15.675;
fuel fuel 1 44.0962 200 5000 1000 7.53414 0.0188722 -6.27185e-06 9.14756e-10 -4.78381e-14 -16467.5 -17.8923 0.933554 0.0264246 6.10597e-06 -2.19775e-08 9.51493e-12 -13958.5 19.2017 1.67212e-06 170.672;
oxidant oxidant 1 28.8504 200 6000 1000 3.10205 0.00123963 -4.17512e-07 6.60292e-11 -3.87448e-15 -985.517 5.35187 3.58378 -0.0007269 1.66985e-06 -1.08452e-10 -4.31951e-13 -1050.53 3.11223 1.67212e-06 170.672;
reactants reactants 24.8095 29.4649 200 5000 1000 3.28069 0.00195035 -6.53483e-07 1.00239e-10 -5.64653e-15 -1609.55 4.41496 3.47696 0.000367499 1.84866e-06 -9.8993e-10 -3.10214e-14 -1570.81 3.76075 1.67212e-06 170.672;
products products 1 28.3233 200 5000 1000 3.106 0.00179682 -5.94382e-07 9.04998e-11 -5.08033e-15 -11003.7 5.11872 3.49612 0.000650364 -2.08029e-07 1.2291e-09 -7.73697e-13 -11080.3 3.18978 1.67212e-06 170.672;
burntProducts burntProducts 25.8095 28.3233 200 6000 1000 3.106 0.00179682 -5.94382e-07 9.04998e-11 -5.08033e-15 -11003.7 5.11872 3.49612 0.000650364 -2.08029e-07 1.2291e-09 -7.73697e-13 -11080.3 3.18978 1.67212e-06 170.672;
thermoType
hhuMixtureThermo<homogeneousMixture<sutherlandTransport<specieThermo<janafThermo<perfectGas>>>>>;
stoichiometricAirFuelMassRatio
stoichiometricAirFuelMassRatio [0 0 0 0 0 0 0] 15.675;
// phi = 0.57
// ft = 0.0352993
// Tad = 1650.1
fuel fuel 1 44.0962
100 5000 1000
7.53414 0.0188722 -6.27185e-06 9.14756e-10 -4.78381e-14 -16467.5 -17.8923
0.933554 0.0264246 6.10597e-06 -2.19775e-08 9.51493e-12 -13958.5 19.2017
1.67212e-06 170.672;
oxidant oxidant 1 28.8504
100 6000 1000
3.10131 0.00124137 -4.18816e-07 6.64158e-11 -3.91274e-15 -985.266 5.35597
3.58378 -0.000727005 1.67057e-06 -1.09203e-10 -4.31765e-13 -1050.53 3.11239
1.67212e-06 170.672;
reactants reactants 1 29.2068
100 5000 1000
3.20495 0.00165359 -5.55661e-07 8.62503e-11 -4.93973e-15 -1347.25 4.81241
3.52181 -9.21936e-05 1.77427e-06 -6.2049e-10 -1.99209e-13 -1352.32 3.48856
1.67212e-06 170.672;
burntProducts burntProducts 1 28.3233
100 6000 1000
3.10558 0.00179747 -5.94696e-07 9.05605e-11 -5.08443e-15 -11003.6 5.12104
3.49796 0.000638555 -1.83881e-07 1.20989e-09 -7.68691e-13 -11080.5 3.18188
1.67212e-06 170.672;
products products 1 28.5396
100 5000 1000
3.10383 0.00156927 -5.22523e-07 8.06527e-11 -4.60363e-15 -6892.54 5.21744
3.53318 7.81943e-05 5.77097e-07 6.68595e-10 -6.30433e-13 -6964.71 3.15336
1.67212e-06 170.672;
// ************************************************************************* //
......@@ -23,7 +23,7 @@ startTime 0;
stopAt endTime;
endTime 10;
endTime 0.5;
deltaT 5e-06;
......
......@@ -17,34 +17,29 @@ FoamFile
ddtSchemes
{
default Euler;
default backward;
}
gradSchemes
{
default Gauss linear;
grad(p) Gauss linear;
}
divSchemes
{
default none;
div(phi,U) Gauss limitedLinearV 1;
div(phid,p) Gauss linear;
div(phi,U) Gauss linear;
div(phiU,p) Gauss linear;
div(phi,k) Gauss limitedLinear 1;
div(phi,B) Gauss limitedLinear 1;
div(B) Gauss linear;
div(phiXi,Xi) Gauss limitedLinear 1;
div(phiXi,Su) Gauss limitedLinear 1;
div(phiSt,b) Gauss limitedLinear01 1;
div(phi,k) Gauss limitedLinear 0.1;
div(phiXi,Xi) Gauss limitedLinear01 0.1;
div(phiXi,Su) Gauss limitedLinear01 0.1;
div(phiSt,b) Gauss limitedLinear01 0.1;
div(phi,ft_b_h_hu) Gauss multivariateSelection
{
fu limitedLinear01 1;
ft limitedLinear01 1;
b limitedLinear01 1;
h limitedLinear 1;
hu limitedLinear 1;
ft limitedLinear01 0.1;
b limitedLinear01 0.1;
h limitedLinear 0.1;
hu limitedLinear 0.1;
};
div(U) Gauss linear;
div((Su*grad(b))) Gauss linear;
......
......@@ -124,8 +124,8 @@ solvers
PISO
{
nOuterCorrectors 1;
nCorrectors 2;
nOuterCorrectors 2;
nCorrectors 1;
nNonOrthogonalCorrectors 0;
}
......
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.5 |
| \\ / A nd | Web: http://www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volTensorField;
object B;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -1 -2 0 0 0 0];
internalField uniform (0 0 0 0 0 0 0 0 0);
boundaryField
{
inlet
{
type fixedValue;
value uniform (0 0 0 0 0 0 0 0 0);
}
outlet
{
type zeroGradient;
}
upperWall
{
type zeroGradient;
}
lowerWall
{
type zeroGradient;
}
frontAndBack
{
type cyclic;
}
}
// ************************************************************************* //
......@@ -20,30 +20,30 @@ internalField uniform 0.135;
boundaryField
{
inlet
inlet
{
type fixedValue;
value uniform 0.135;
}
outlet
outlet
{
type inletOutlet;
inletValue uniform 0.135;
value uniform 0.135;
}
upperWall
upperWall
{
type zeroGradient;
}
lowerWall
lowerWall
{
type zeroGradient;
}
frontAndBack
frontAndBack
{
type cyclic;
}
......
......@@ -20,32 +20,32 @@ internalField uniform 293;
boundaryField
{
inlet
inlet
{
type fixedValue;
value uniform 293;
}
outlet
outlet
{
type inletOutlet;
inletValue uniform 293;
value uniform 293;
}
upperWall
upperWall
{
type fixedValue;
value uniform 293;
}
lowerWall
lowerWall
{
type fixedValue;
value uniform 293;
value uniform 570;
}
frontAndBack
frontAndBack
{
type cyclic;
}
......
......@@ -20,32 +20,32 @@ internalField uniform 293;
boundaryField
{
inlet
inlet
{
type fixedValue;
value uniform 293;
}
outlet
outlet
{
type inletOutlet;
inletValue uniform 293;
value uniform 293;
}
upperWall
upperWall
{
type fixedValue;
value uniform 293;
}
lowerWall
lowerWall
{
type fixedValue;
value uniform 293;
value uniform 570;
}
frontAndBack
frontAndBack
{
type cyclic;
}
......
......@@ -20,33 +20,33 @@ internalField uniform (0 0 0);
boundaryField
{
inlet
inlet
{
type turbulentInlet;
referenceField uniform (13.3 0 0);
fluctuationScale (0.04 0.02 0.02);
}
outlet
outlet
{
type inletOutlet;
inletValue uniform (0 0 0);
value uniform (0 0 0);
}
upperWall
upperWall
{
type fixedValue;
value uniform (0 0 0);
}
lowerWall
lowerWall
{
type fixedValue;
value uniform (0 0 0);
}
frontAndBack
frontAndBack
{
type cyclic;
}
......
......@@ -20,30 +20,30 @@ internalField uniform 1;
boundaryField
{
inlet
inlet
{
type fixedValue;
value uniform 1;
}
outlet
outlet
{
type inletOutlet;
inletValue uniform 1;
value uniform 1;
}