diff --git a/tutorials/combustion/XiFoam/les/pitzDaily/0/T b/tutorials/combustion/XiFoam/les/pitzDaily/0/T index c4bae6fa13..d48433513d 100644 --- a/tutorials/combustion/XiFoam/les/pitzDaily/0/T +++ b/tutorials/combustion/XiFoam/les/pitzDaily/0/T @@ -42,7 +42,7 @@ boundaryField lowerWall { type fixedValue; - value uniform 293; + value uniform 570; } frontAndBack diff --git a/tutorials/combustion/XiFoam/les/pitzDaily/0/Tu b/tutorials/combustion/XiFoam/les/pitzDaily/0/Tu index 3e1d003deb..55c99ebb45 100644 --- a/tutorials/combustion/XiFoam/les/pitzDaily/0/Tu +++ b/tutorials/combustion/XiFoam/les/pitzDaily/0/Tu @@ -42,7 +42,7 @@ boundaryField lowerWall { type fixedValue; - value uniform 293; + value uniform 570; } frontAndBack diff --git a/tutorials/combustion/XiFoam/les/pitzDaily/0/alphaSgs b/tutorials/combustion/XiFoam/les/pitzDaily/0/alphaSgs index 24c919196b..c66c8a1534 100644 --- a/tutorials/combustion/XiFoam/les/pitzDaily/0/alphaSgs +++ b/tutorials/combustion/XiFoam/les/pitzDaily/0/alphaSgs @@ -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; } diff --git a/tutorials/combustion/XiFoam/les/pitzDaily/0/muSgs b/tutorials/combustion/XiFoam/les/pitzDaily/0/muSgs index a16e89fad9..5c85d09493 100644 --- a/tutorials/combustion/XiFoam/les/pitzDaily/0/muSgs +++ b/tutorials/combustion/XiFoam/les/pitzDaily/0/muSgs @@ -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; } diff --git a/tutorials/combustion/XiFoam/les/pitzDaily/0/p b/tutorials/combustion/XiFoam/les/pitzDaily/0/p index 3700ec72b5..b8afbd8d1a 100644 --- a/tutorials/combustion/XiFoam/les/pitzDaily/0/p +++ b/tutorials/combustion/XiFoam/les/pitzDaily/0/p @@ -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 diff --git a/tutorials/combustion/XiFoam/les/pitzDaily/constant/LESProperties b/tutorials/combustion/XiFoam/les/pitzDaily/constant/LESProperties index 56c915c924..faa025c8e2 100644 --- a/tutorials/combustion/XiFoam/les/pitzDaily/constant/LESProperties +++ b/tutorials/combustion/XiFoam/les/pitzDaily/constant/LESProperties @@ -123,6 +123,7 @@ kappa 0.4187; wallFunctionCoeffs { E 9; + Prt 0.85; } diff --git a/tutorials/combustion/XiFoam/les/pitzDaily/constant/combustionProperties b/tutorials/combustion/XiFoam/les/pitzDaily/constant/combustionProperties index 86fcfc889a..a20d6ae7ec 100644 --- a/tutorials/combustion/XiFoam/les/pitzDaily/constant/combustionProperties +++ b/tutorials/combustion/XiFoam/les/pitzDaily/constant/combustionProperties @@ -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; } ); diff --git a/tutorials/combustion/XiFoam/les/pitzDaily/constant/polyMesh/blockMeshDict b/tutorials/combustion/XiFoam/les/pitzDaily/constant/polyMesh/blockMeshDict index 1030983582..4a24aa9a9d 100644 --- a/tutorials/combustion/XiFoam/les/pitzDaily/constant/polyMesh/blockMeshDict +++ b/tutorials/combustion/XiFoam/les/pitzDaily/constant/polyMesh/blockMeshDict @@ -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) diff --git a/tutorials/combustion/XiFoam/les/pitzDaily/constant/thermophysicalProperties b/tutorials/combustion/XiFoam/les/pitzDaily/constant/thermophysicalProperties index 8e705f32e5..dc9f1dae0f 100644 --- a/tutorials/combustion/XiFoam/les/pitzDaily/constant/thermophysicalProperties +++ b/tutorials/combustion/XiFoam/les/pitzDaily/constant/thermophysicalProperties @@ -15,19 +15,44 @@ FoamFile } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // -thermoType hhuMixtureThermo>>>>; +thermoType + hhuMixtureThermo>>>>; -stoichiometricAirFuelMassRatio stoichiometricAirFuelMassRatio [ 0 0 0 0 0 0 0 ] 15.675; +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; +// phi = 0.57 +// ft = 0.0352993 +// Tad = 1650.1 -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; +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; -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; +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; -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; +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 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; +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; // ************************************************************************* // diff --git a/tutorials/combustion/XiFoam/les/pitzDaily/system/controlDict b/tutorials/combustion/XiFoam/les/pitzDaily/system/controlDict index f5fe65e622..ba28fbf560 100644 --- a/tutorials/combustion/XiFoam/les/pitzDaily/system/controlDict +++ b/tutorials/combustion/XiFoam/les/pitzDaily/system/controlDict @@ -23,7 +23,7 @@ startTime 0; stopAt endTime; -endTime 10; +endTime 0.5; deltaT 5e-06; diff --git a/tutorials/combustion/XiFoam/les/pitzDaily/system/fvSchemes b/tutorials/combustion/XiFoam/les/pitzDaily/system/fvSchemes index e74dbeb191..d4a2457583 100644 --- a/tutorials/combustion/XiFoam/les/pitzDaily/system/fvSchemes +++ b/tutorials/combustion/XiFoam/les/pitzDaily/system/fvSchemes @@ -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; diff --git a/tutorials/combustion/XiFoam/les/pitzDaily/system/fvSolution b/tutorials/combustion/XiFoam/les/pitzDaily/system/fvSolution index 6e1d1ffd3e..4dda010d66 100644 --- a/tutorials/combustion/XiFoam/les/pitzDaily/system/fvSolution +++ b/tutorials/combustion/XiFoam/les/pitzDaily/system/fvSolution @@ -124,8 +124,8 @@ solvers PISO { - nOuterCorrectors 1; - nCorrectors 2; + nOuterCorrectors 2; + nCorrectors 1; nNonOrthogonalCorrectors 0; } diff --git a/tutorials/combustion/XiFoam/les/pitzDaily3D/0/B b/tutorials/combustion/XiFoam/les/pitzDaily3D/0/B deleted file mode 100644 index 65544498af..0000000000 --- a/tutorials/combustion/XiFoam/les/pitzDaily3D/0/B +++ /dev/null @@ -1,50 +0,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; - } -} - -// ************************************************************************* // diff --git a/tutorials/combustion/XiFoam/les/pitzDaily3D/0/Su b/tutorials/combustion/XiFoam/les/pitzDaily3D/0/Su index 82eb5d80e3..024869330e 100644 --- a/tutorials/combustion/XiFoam/les/pitzDaily3D/0/Su +++ b/tutorials/combustion/XiFoam/les/pitzDaily3D/0/Su @@ -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; } diff --git a/tutorials/combustion/XiFoam/les/pitzDaily3D/0/T b/tutorials/combustion/XiFoam/les/pitzDaily3D/0/T index fdf764e06f..9c0947189b 100644 --- a/tutorials/combustion/XiFoam/les/pitzDaily3D/0/T +++ b/tutorials/combustion/XiFoam/les/pitzDaily3D/0/T @@ -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; } diff --git a/tutorials/combustion/XiFoam/les/pitzDaily3D/0/Tu b/tutorials/combustion/XiFoam/les/pitzDaily3D/0/Tu index e1299d7435..e0b9959787 100644 --- a/tutorials/combustion/XiFoam/les/pitzDaily3D/0/Tu +++ b/tutorials/combustion/XiFoam/les/pitzDaily3D/0/Tu @@ -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; } diff --git a/tutorials/combustion/XiFoam/les/pitzDaily3D/0/U b/tutorials/combustion/XiFoam/les/pitzDaily3D/0/U index ece734e9f3..c2ff1ef691 100644 --- a/tutorials/combustion/XiFoam/les/pitzDaily3D/0/U +++ b/tutorials/combustion/XiFoam/les/pitzDaily3D/0/U @@ -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; } diff --git a/tutorials/combustion/XiFoam/les/pitzDaily3D/0/Xi b/tutorials/combustion/XiFoam/les/pitzDaily3D/0/Xi index 1b3b3cc930..ac3b5633fa 100644 --- a/tutorials/combustion/XiFoam/les/pitzDaily3D/0/Xi +++ b/tutorials/combustion/XiFoam/les/pitzDaily3D/0/Xi @@ -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; } - upperWall + upperWall { type zeroGradient; } - lowerWall + lowerWall { type zeroGradient; } - frontAndBack + frontAndBack { type cyclic; } diff --git a/tutorials/combustion/XiFoam/les/pitzDaily3D/0/alphaSgs b/tutorials/combustion/XiFoam/les/pitzDaily3D/0/alphaSgs index c5331748ae..5f7deec4ec 100644 --- a/tutorials/combustion/XiFoam/les/pitzDaily3D/0/alphaSgs +++ b/tutorials/combustion/XiFoam/les/pitzDaily3D/0/alphaSgs @@ -20,27 +20,29 @@ internalField uniform 0; boundaryField { - inlet + inlet { type zeroGradient; } - outlet + outlet { type zeroGradient; } - upperWall + upperWall { - type zeroGradient; + type compressible::alphaSgsJayatillekeWallFunction; + value uniform 0; } - lowerWall + lowerWall { - type zeroGradient; + type compressible::alphaSgsJayatillekeWallFunction; + value uniform 0; } - frontAndBack + frontAndBack { type cyclic; } diff --git a/tutorials/combustion/XiFoam/les/pitzDaily3D/0/b b/tutorials/combustion/XiFoam/les/pitzDaily3D/0/b index 7be4cfc6ac..7ca894170f 100644 --- a/tutorials/combustion/XiFoam/les/pitzDaily3D/0/b +++ b/tutorials/combustion/XiFoam/les/pitzDaily3D/0/b @@ -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; } - upperWall + upperWall { type zeroGradient; } - lowerWall + lowerWall { type zeroGradient; } - frontAndBack + frontAndBack { type cyclic; } diff --git a/tutorials/combustion/XiFoam/les/pitzDaily3D/0/ft b/tutorials/combustion/XiFoam/les/pitzDaily3D/0/ft deleted file mode 100644 index 8d16f22209..0000000000 --- a/tutorials/combustion/XiFoam/les/pitzDaily3D/0/ft +++ /dev/null @@ -1,52 +0,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 volScalarField; - object ft; -} -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // - -dimensions [0 0 0 0 0 0 0]; - -internalField uniform 0; - -boundaryField -{ - inlet - { - type fixedValue; - value uniform 0; - } - - outlet - { - type inletOutlet; - inletValue uniform 0; - value uniform 0; - } - - upperWall - { - type zeroGradient; - } - - lowerWall - { - type zeroGradient; - } - - frontAndBack - { - type cyclic; - } -} - -// ************************************************************************* // diff --git a/tutorials/combustion/XiFoam/les/pitzDaily3D/0/k b/tutorials/combustion/XiFoam/les/pitzDaily3D/0/k index 69de845314..df6b081d5b 100644 --- a/tutorials/combustion/XiFoam/les/pitzDaily3D/0/k +++ b/tutorials/combustion/XiFoam/les/pitzDaily3D/0/k @@ -35,14 +35,14 @@ boundaryField upperWall { - type fixedValue; - value uniform 1e-8; + type zeroGradient; + value uniform 2e-05; } lowerWall { - type fixedValue; - value uniform 1e-8; + type zeroGradient; + value uniform 2e-05; } frontAndBack diff --git a/tutorials/combustion/XiFoam/les/pitzDaily3D/0/muSgs b/tutorials/combustion/XiFoam/les/pitzDaily3D/0/muSgs index a6befb8a9d..ff6a524d0d 100644 --- a/tutorials/combustion/XiFoam/les/pitzDaily3D/0/muSgs +++ b/tutorials/combustion/XiFoam/les/pitzDaily3D/0/muSgs @@ -20,27 +20,29 @@ internalField uniform 0; boundaryField { - inlet + inlet { type zeroGradient; } - outlet + outlet { type zeroGradient; } - upperWall + upperWall { - type zeroGradient; + type compressible::muSgsWallFunction; + value uniform 0; } - lowerWall + lowerWall { - type zeroGradient; + type compressible::muSgsWallFunction; + value uniform 0; } - frontAndBack + frontAndBack { type cyclic; } diff --git a/tutorials/combustion/XiFoam/les/pitzDaily3D/0/p b/tutorials/combustion/XiFoam/les/pitzDaily3D/0/p index 918b79c476..2a27d1e401 100644 --- a/tutorials/combustion/XiFoam/les/pitzDaily3D/0/p +++ b/tutorials/combustion/XiFoam/les/pitzDaily3D/0/p @@ -16,32 +16,39 @@ FoamFile dimensions [1 -1 -2 0 0 0 0]; -internalField uniform 100000; +internalField uniform 1e5; boundaryField { - inlet + inlet { type zeroGradient; } - outlet + 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 + upperWall { type zeroGradient; } - lowerWall + lowerWall { type zeroGradient; } - frontAndBack + frontAndBack { type cyclic; } diff --git a/tutorials/combustion/XiFoam/les/pitzDaily3D/constant/LESProperties b/tutorials/combustion/XiFoam/les/pitzDaily3D/constant/LESProperties index 95bb1c14f8..faa025c8e2 100644 --- a/tutorials/combustion/XiFoam/les/pitzDaily3D/constant/LESProperties +++ b/tutorials/combustion/XiFoam/les/pitzDaily3D/constant/LESProperties @@ -1,7 +1,7 @@ /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | -| \\ / O peration | Version: dev | +| \\ / O peration | Version: 1.5 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ @@ -123,6 +123,7 @@ kappa 0.4187; wallFunctionCoeffs { E 9; + Prt 0.85; } diff --git a/tutorials/combustion/XiFoam/les/pitzDaily3D/constant/combustionProperties b/tutorials/combustion/XiFoam/les/pitzDaily3D/constant/combustionProperties index 8d2c6962d7..f46b08c801 100644 --- a/tutorials/combustion/XiFoam/les/pitzDaily3D/constant/combustionProperties +++ b/tutorials/combustion/XiFoam/les/pitzDaily3D/constant/combustionProperties @@ -68,16 +68,16 @@ GuldersCoeffs } } -ignite yes; +ignite yes; ignitionSites ( { - location (0 0 0); + location (0.005 -0.02 0.01905); diameter 0.003; start 0; - duration 0.001; - strength 2; + duration 0.1; + strength 200; } ); diff --git a/tutorials/combustion/XiFoam/les/pitzDaily3D/constant/g b/tutorials/combustion/XiFoam/les/pitzDaily3D/constant/g index d182180d95..7b4411955d 100644 --- a/tutorials/combustion/XiFoam/les/pitzDaily3D/constant/g +++ b/tutorials/combustion/XiFoam/les/pitzDaily3D/constant/g @@ -1,7 +1,7 @@ /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | -| \\ / O peration | Version: dev | +| \\ / O peration | Version: 1.5 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ diff --git a/tutorials/combustion/XiFoam/les/pitzDaily3D/constant/thermophysicalProperties b/tutorials/combustion/XiFoam/les/pitzDaily3D/constant/thermophysicalProperties index 7c8fd76fe3..dc9f1dae0f 100644 --- a/tutorials/combustion/XiFoam/les/pitzDaily3D/constant/thermophysicalProperties +++ b/tutorials/combustion/XiFoam/les/pitzDaily3D/constant/thermophysicalProperties @@ -1,7 +1,7 @@ /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | -| \\ / O peration | Version: dev | +| \\ / O peration | Version: 1.5 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ @@ -15,19 +15,44 @@ FoamFile } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // -thermoType hhuMixtureThermo>>>>; +thermoType + hhuMixtureThermo>>>>; -stoichiometricAirFuelMassRatio stoichiometricAirFuelMassRatio [ 0 0 0 0 0 0 0 ] 15.675; +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; +// phi = 0.57 +// ft = 0.0352993 +// Tad = 1650.1 -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; +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; -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; +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; -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; +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 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; +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; // ************************************************************************* // diff --git a/tutorials/combustion/XiFoam/les/pitzDaily3D/system/controlDict b/tutorials/combustion/XiFoam/les/pitzDaily3D/system/controlDict index cd473319ff..364c823851 100644 --- a/tutorials/combustion/XiFoam/les/pitzDaily3D/system/controlDict +++ b/tutorials/combustion/XiFoam/les/pitzDaily3D/system/controlDict @@ -17,13 +17,13 @@ FoamFile application XiFoam; -startFrom startTime; +startFrom latestTime; startTime 0; stopAt endTime; -endTime 0.02; +endTime 0.5; deltaT 5e-06; diff --git a/tutorials/combustion/XiFoam/les/pitzDaily3D/system/fvSchemes b/tutorials/combustion/XiFoam/les/pitzDaily3D/system/fvSchemes index 2398134a00..d4a2457583 100644 --- a/tutorials/combustion/XiFoam/les/pitzDaily3D/system/fvSchemes +++ b/tutorials/combustion/XiFoam/les/pitzDaily3D/system/fvSchemes @@ -1,7 +1,7 @@ /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | -| \\ / O peration | Version: dev | +| \\ / O peration | Version: 1.5 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ @@ -17,28 +17,30 @@ 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,ft_b_h_hu) Gauss multivariateSelection { fu limitedLinear01 1 ; ft limitedLinear01 1 ; b limitedLinear01 1 ; h limitedLinear 1 ; hu limitedLinear 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 + { + 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; div((U+((Su*Xi)*grad(b)))) Gauss linear; @@ -52,8 +54,8 @@ laplacianSchemes laplacian(DkEff,k) Gauss linear corrected; laplacian(DBEff,B) Gauss linear corrected; laplacian((rho*(1|A(U))),p) Gauss linear corrected; - laplacian(muEff,b) Gauss linear corrected; - laplacian(muEff,ft) Gauss linear corrected; + laplacian(alphaEff,b) Gauss linear corrected; + laplacian(alphaEff,ft) Gauss linear corrected; laplacian(alphaEff,h) Gauss linear corrected; laplacian(alphaEff,hu) Gauss linear corrected; } diff --git a/tutorials/combustion/XiFoam/les/pitzDaily3D/system/fvSolution b/tutorials/combustion/XiFoam/les/pitzDaily3D/system/fvSolution index d4656531ed..4dda010d66 100644 --- a/tutorials/combustion/XiFoam/les/pitzDaily3D/system/fvSolution +++ b/tutorials/combustion/XiFoam/les/pitzDaily3D/system/fvSolution @@ -1,7 +1,7 @@ /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | -| \\ / O peration | Version: dev | +| \\ / O peration | Version: 1.5 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ @@ -124,7 +124,8 @@ solvers PISO { - nCorrectors 2; + nOuterCorrectors 2; + nCorrectors 1; nNonOrthogonalCorrectors 0; }