From edcc472535b1a7314fa7ec6c1946ee31a333bbe7 Mon Sep 17 00:00:00 2001 From: Henry Weller Date: Fri, 26 Jun 2015 22:05:39 +0100 Subject: [PATCH] reactingTwoPhaseEulerFoam: Removed superfluous trSubDeltaT field --- .../reactingTwoPhaseEulerFoam/createRDeltaT.H | 47 ------------------- .../reactingTwoPhaseEulerFoam/setRDeltaT.H | 6 --- 2 files changed, 53 deletions(-) delete mode 100644 applications/solvers/multiphase/reactingTwoPhaseEulerFoam/createRDeltaT.H diff --git a/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/createRDeltaT.H b/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/createRDeltaT.H deleted file mode 100644 index b45dff137e..0000000000 --- a/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/createRDeltaT.H +++ /dev/null @@ -1,47 +0,0 @@ -bool LTS = - word(mesh.ddtScheme("default")) - == fv::localEulerDdtScheme::typeName; - -tmp trDeltaT; -tmp trSubDeltaT; - -if (LTS) -{ - scalar maxDeltaT - ( - pimple.dict().lookupOrDefault("maxDeltaT", GREAT) - ); - - trDeltaT = tmp - ( - new volScalarField - ( - IOobject - ( - "rDeltaT", - runTime.timeName(), - mesh, - IOobject::NO_READ, - IOobject::AUTO_WRITE - ), - mesh, - 1/dimensionedScalar("maxDeltaT", dimTime, maxDeltaT), - zeroGradientFvPatchScalarField::typeName - ) - ); - - trSubDeltaT = tmp - ( - new volScalarField - ( - IOobject - ( - "rSubDeltaT", - runTime.timeName(), - mesh - ), - mesh, - 1/dimensionedScalar("maxDeltaT", dimTime, maxDeltaT) - ) - ); -} diff --git a/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/setRDeltaT.H b/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/setRDeltaT.H index 99e08c3f56..f9f736f788 100644 --- a/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/setRDeltaT.H +++ b/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/setRDeltaT.H @@ -1,6 +1,5 @@ { volScalarField& rDeltaT = trDeltaT(); - volScalarField& rSubDeltaT = trSubDeltaT(); scalar rDeltaTSmoothingCoeff ( @@ -27,9 +26,4 @@ Info<< "Flow time scale min/max = " << gMin(1/rDeltaT.internalField()) << ", " << gMax(1/rDeltaT.internalField()) << endl; - - const dictionary& alphaControls = mesh.solverDict(alpha1.name()); - label nAlphaSubCycles(readLabel(alphaControls.lookup("nAlphaSubCycles"))); - - rSubDeltaT = rDeltaT*nAlphaSubCycles; }