From f406c36e4a847fc51d9ab1175450d97c885b6d23 Mon Sep 17 00:00:00 2001 From: space Date: Mon, 14 Aug 2017 14:21:43 +0900 Subject: [PATCH 1/2] Schmidt Number, Outlet Name information edit --- solvers_post/SLFMFoam/Mixturefraction.H | 2 +- solvers_post/SLFMFoam/MixturefractionVar.H | 4 ++-- solvers_post/SLFMFoam/SLFMFoam.C | 14 +++++++++----- solvers_post/SLFMFoam/createFields.H | 2 ++ solvers_post/SLFMFoam/readSLFMlib.H | 4 ++-- solvers_post/SLFMFoam/updateT_RHO.H | 2 +- 6 files changed, 17 insertions(+), 11 deletions(-) diff --git a/solvers_post/SLFMFoam/Mixturefraction.H b/solvers_post/SLFMFoam/Mixturefraction.H index 79d73d1..f4c7a20 100644 --- a/solvers_post/SLFMFoam/Mixturefraction.H +++ b/solvers_post/SLFMFoam/Mixturefraction.H @@ -3,7 +3,7 @@ fvScalarMatrix mfEqn ( fvm::div(phi, mf) - - fvm::laplacian(1.47*turbulence->mut(), mf) //let 1/Sc = 1.47 + - fvm::laplacian((1/Sc)*turbulence->mut(), mf) //let 1/Sc = 1.47 ); diff --git a/solvers_post/SLFMFoam/MixturefractionVar.H b/solvers_post/SLFMFoam/MixturefractionVar.H index 2db84b5..b05b582 100644 --- a/solvers_post/SLFMFoam/MixturefractionVar.H +++ b/solvers_post/SLFMFoam/MixturefractionVar.H @@ -7,9 +7,9 @@ fvScalarMatrix mfVarEqn ( fvm::div(phi, mfVar) - - fvm::laplacian(1.47*turbulence->mut(), mfVar) //let 1/Sc = 1.47 + - fvm::laplacian((1/Sc)*turbulence->mut(), mfVar) == - 2 * ( 1.47 * turbulence->mut() ) * (Gradmf & Gradmf) + 2 * ( (1/Sc) * turbulence->mut() ) * (Gradmf & Gradmf) - 2 * rho * SDR ); diff --git a/solvers_post/SLFMFoam/SLFMFoam.C b/solvers_post/SLFMFoam/SLFMFoam.C index 70e2244..d434ede 100644 --- a/solvers_post/SLFMFoam/SLFMFoam.C +++ b/solvers_post/SLFMFoam/SLFMFoam.C @@ -25,10 +25,14 @@ Application SLFMFoam Description - Steady state solver for turbulent combustion. + Steady Laminar Flamelet Model(SLFM) solver for turbulent combustion. SLFM + assumes a turbulent flame composed of thin stretched laminar flamelets. + Each flamelet library is generated by external program and imported by + the solver. References - N. Peters, Turbulent Combustion, Cambridge University Press; 2000. + A.Y. Klimenko, R.W. Bilger, Progress in Energy and Combustion Science 25 (1999) 595-687 + N. Peters, Turbulent Combustion, Cambridge University Press (2000) Contact POSTECH combustion lab. @@ -55,7 +59,7 @@ int main(int argc, char *argv[]) #include "createControl.H" #include "createFields.H" #include "createFieldRefs.H" - #include "BetaPDF.H" //make detailed integration space and evaluate beta-pdf integration + #include "BetaPDF.H" //make detailed integration space and perform beta-pdf integration #include "createFvOptions.H" #include "initContinuityErrs.H" @@ -81,8 +85,8 @@ int main(int argc, char *argv[]) // --- Pressure-velocity SIMPLE corrector loop { #include "UEqn.H" - #include "Mixturefraction.H" //mixture fraction - #include "MixturefractionVar.H" //mixture fraction variance + #include "Mixturefraction.H" //mean mixture fraction + #include "MixturefractionVar.H" //mean mixture fraction variance #include "setSDR.H" //calculate scalar dissipation rate #include "updateT_RHO.H" //update temperature and density #include "pEqn.H" diff --git a/solvers_post/SLFMFoam/createFields.H b/solvers_post/SLFMFoam/createFields.H index e5b53c9..0b611ba 100644 --- a/solvers_post/SLFMFoam/createFields.H +++ b/solvers_post/SLFMFoam/createFields.H @@ -308,7 +308,9 @@ scalar stoiMF(readScalar(SLFMdict.lookup("stoiMF"))); scalar stretchFac(readScalar(SLFMdict.lookup("stretchFactor"))); scalar lowerN(readScalar(SLFMdict.lookup("lowerN"))); scalar upperN(readScalar(SLFMdict.lookup("upperN"))); +scalar Sc(readScalar(SLFMdict.lookup("Sc"))); scalar etamax(lowerN+upperN); +word outletName(SLFMdict.lookup("outletName")); scalarField etaValue(etamax+1,0); scalarField etaIndex(etamax+1,0); diff --git a/solvers_post/SLFMFoam/readSLFMlib.H b/solvers_post/SLFMFoam/readSLFMlib.H index dc26a22..f58e724 100644 --- a/solvers_post/SLFMFoam/readSLFMlib.H +++ b/solvers_post/SLFMFoam/readSLFMlib.H @@ -30,12 +30,12 @@ scalarFieldArray2d W_SLFM //reaction rate [1/sec] scalarField(etamax+1,0) ) ); -scalarFieldArray1d rho_SLFM //density [g/cm3 ?] +scalarFieldArray1d rho_SLFM //density [g/cm3] ( NstList.size(), scalarField(etamax+1, 0) ); -scalarFieldArray1d h_SLFM //enthalpy [?] +scalarFieldArray1d h_SLFM //enthalpy [erg/gm] ( NstList.size(), scalarField(etamax+1, 0) diff --git a/solvers_post/SLFMFoam/updateT_RHO.H b/solvers_post/SLFMFoam/updateT_RHO.H index cf6e310..3099a32 100644 --- a/solvers_post/SLFMFoam/updateT_RHO.H +++ b/solvers_post/SLFMFoam/updateT_RHO.H @@ -63,7 +63,7 @@ const fvPatchList& patches = mesh.boundary(); forAll(patches, patchI) { const fvPatch& curPatch = patches[patchI]; - if(curPatch.name() == "outlet") + if(curPatch.name() == outletName) { forAll(curPatch, faceI) { From 08b013fa35ca417a479de838c54d94dfb8aece24 Mon Sep 17 00:00:00 2001 From: space Date: Wed, 16 Aug 2017 09:30:09 +0900 Subject: [PATCH 2/2] examples updated --- examples/FlameD_Equilibrium/constant/SLFMdict | 3 +++ examples/FlameD_SLFM/constant/SLFMdict | 3 +++ 2 files changed, 6 insertions(+) diff --git a/examples/FlameD_Equilibrium/constant/SLFMdict b/examples/FlameD_Equilibrium/constant/SLFMdict index 85fc436..d85dd74 100644 --- a/examples/FlameD_Equilibrium/constant/SLFMdict +++ b/examples/FlameD_Equilibrium/constant/SLFMdict @@ -23,6 +23,9 @@ upperN 30; //number of eta-space over stoichiometric value Equilibrium true; NVar 40; //20; +Sc 0.68; //Schmidt number +outletName outlet; // outletName + detailedEta ( 1e-6 1e-5 1e-4 1e-3 1e-2 diff --git a/examples/FlameD_SLFM/constant/SLFMdict b/examples/FlameD_SLFM/constant/SLFMdict index ef97d9f..de0f9fc 100644 --- a/examples/FlameD_SLFM/constant/SLFMdict +++ b/examples/FlameD_SLFM/constant/SLFMdict @@ -24,6 +24,9 @@ upperN 30; //number of eta-space over stoichiometric value Equilibrium false; NVar 40; //20; +Sc 0.68; //Schmidt number +outletName outlet; // outletName + detailedEta ( 1e-6 1e-5 1e-4 1e-3 1e-2