diff --git a/applications/solvers/lagrangian/LTSReactingParcelFoam/readTimeControls.H b/applications/solvers/lagrangian/LTSReactingParcelFoam/readTimeControls.H index 5513848c647124e13a8e6ebbcf4214ce2eda0e05..b73efd593e1bcc40490583f034237dbec2b728ea 100644 --- a/applications/solvers/lagrangian/LTSReactingParcelFoam/readTimeControls.H +++ b/applications/solvers/lagrangian/LTSReactingParcelFoam/readTimeControls.H @@ -27,13 +27,16 @@ License scalar maxCo(readScalar(pimple.dict().lookup("maxCo"))); // Maximum time scale -scalar maxDeltaT = readScalar(pimple.dict().lookup("maxDeltaT")); +scalar maxDeltaT(pimple.dict().lookupOrDefault<scalar>("maxDeltaT", GREAT)); // Smoothing parameter (0-1) when smoothing iterations > 0 -scalar alphaTauSmooth(pimple.dict().lookupOrDefault("alphaTauSmooth", 0.1)); +scalar rDeltaTSmoothingCoeff +( + pimple.dict().lookupOrDefault<scalar>("rDeltaTSmoothingCoeff", 0.1) +); // Maximum change in cell temperature per iteration (relative to previous value) -scalar alphaTauTemp(pimple.dict().lookupOrDefault("alphaTauTemp", 0.05)); +scalar alphaTemp(pimple.dict().lookupOrDefault("alphaTemp", 0.05)); // ************************************************************************* // diff --git a/applications/solvers/lagrangian/LTSReactingParcelFoam/timeScales.H b/applications/solvers/lagrangian/LTSReactingParcelFoam/timeScales.H index d93b5a40ed0fc18389fd9bad86528657907eed7f..e6d0fa85c4919cb322b30620c2850d408fb0fb8f 100644 --- a/applications/solvers/lagrangian/LTSReactingParcelFoam/timeScales.H +++ b/applications/solvers/lagrangian/LTSReactingParcelFoam/timeScales.H @@ -82,7 +82,7 @@ Info<< "Time scales min/max:" << endl; /rho ); - tau = alphaTauTemp*thermo.Cp()*T/(tau + ROOTVSMALL); + tau = alphaTemp*thermo.Cp()*T/(tau + ROOTVSMALL); Info<< " Temperature = " << min(maxDeltaT, gMin(tau)) << ", " << min(maxDeltaT, gMax(tau)) << endl; @@ -105,7 +105,7 @@ Info<< "Time scales min/max:" << endl; // Spatially smooth the time scale field // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ - fvc::smooth(rDeltaT, alphaTauSmooth); + fvc::smooth(rDeltaT, rDeltaTSmoothingCoeff); Info<< " Overall = " << min(1/rDeltaT).value() << ", " << max(1/rDeltaT).value() << nl << endl; diff --git a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSolution b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSolution index 69d711902521ece7bad44a33b7afcb1419f64240..aaff31c4c7e8a6d5301e142bcb10eb6273d82e1f 100644 --- a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSolution +++ b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSolution @@ -85,8 +85,9 @@ PIMPLE rhoMax rhoMax [1 -3 0 0 0] 1.5; maxCo 0.5; + rDeltaTSmoothingCoeff 0.1; maxDeltaT 1; - alphaTauTemp 0.005; + alphaTemp 0.005; } additional diff --git a/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/system/fvSolution b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/system/fvSolution index b750e0d5e5735db6b82353f9ecd091857b694bef..4dfe0a0beb61e2d751d2038467a8bedf9c5eb3a4 100644 --- a/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/system/fvSolution +++ b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/system/fvSolution @@ -72,8 +72,8 @@ PIMPLE rhoMax rhoMax [1 -3 0 0 0] 1.5; maxCo 5; - alphaTauSmooth 1; - alphaTauTemp 1; + rDeltaTSmoothingCoeff 1; + alphaTemp 1; maxDeltaT 1; }