Skip to content
Snippets Groups Projects
Commit 3051d333 authored by andy's avatar andy
Browse files

STYLE: Consistency updates for LTS solver

parent 2f374d0d
Branches
Tags
No related merge requests found
......@@ -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));
// ************************************************************************* //
......@@ -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;
......
......@@ -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
......
......@@ -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;
}
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment