zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

ENH: chtMultiRegionFoam solvers - added option to calculate using frozen flow

Browse Source
This commit is contained in:
andy
2014-10-16 10:20:19 +01:00
committed by Andrew Heather
parent e08a22afd1
commit ade2801158
6 changed files with 54 additions and 46 deletions
Download Patch File Download Diff File Expand all files Collapse all files

35
applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/fluid/createFluidFields.H
View File

@ -11,8 +11,9 @@
PtrList<radiation::radiationModel> radiation(fluidRegions.size());
List<scalar> initialMassFluid(fluidRegions.size());
List<label> pRefCellFluid(fluidRegions.size(),0);
List<scalar> pRefValueFluid(fluidRegions.size(),0.0);
List<label> pRefCellFluid(fluidRegions.size(), 0);
List<scalar> pRefValueFluid(fluidRegions.size(), 0.0);
List<bool> frozenFlowFluid(fluidRegions.size(), false);
PtrList<dimensionedScalar> rhoMax(fluidRegions.size());
PtrList<dimensionedScalar> rhoMin(fluidRegions.size());
@ -160,38 +161,22 @@
initialMassFluid[i] = fvc::domainIntegrate(rhoFluid[i]).value();
const dictionary& simpleDict =
fluidRegions[i].solutionDict().subDict("SIMPLE");
setRefCell
(
thermoFluid[i].p(),
p_rghFluid[i],
fluidRegions[i].solutionDict().subDict("SIMPLE"),
simpleDict,
pRefCellFluid[i],
pRefValueFluid[i]
);
rhoMax.set
(
i,
new dimensionedScalar
(
fluidRegions[i].solutionDict().subDict("SIMPLE").lookup
(
"rhoMax"
)
)
);
simpleDict.readIfPresent("frozenFlow", frozenFlowFluid[i]);
rhoMin.set
(
i,
new dimensionedScalar
(
fluidRegions[i].solutionDict().subDict("SIMPLE").lookup
(
"rhoMin"
)
)
);
rhoMax.set(i, new dimensionedScalar(simpleDict.lookup("rhoMax")));
rhoMin.set(i, new dimensionedScalar(simpleDict.lookup("rhoMin")));
Info<< " Adding fvOptions\n" << endl;
fluidFvOptions.set

1
applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/fluid/setRegionFluidFields.H
View File

@ -25,6 +25,7 @@
const label pRefCell = pRefCellFluid[i];
const scalar pRefValue = pRefValueFluid[i];
const bool frozenFlow = frozenFlowFluid[i];
volScalarField& p_rgh = p_rghFluid[i];
const volScalarField& gh = ghFluid[i];

22
applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/fluid/solveFluid.H
View File

@ -1,11 +1,19 @@
// Pressure-velocity SIMPLE corrector
p_rgh.storePrevIter();
rho.storePrevIter();
{
#include "UEqn.H"
#include "EEqn.H"
#include "pEqn.H"
}
if (frozenFlow)
{
#include "EEqn.H"
}
else
{
p_rgh.storePrevIter();
rho.storePrevIter();
turb.correct();
#include "UEqn.H"
#include "EEqn.H"
#include "pEqn.H"
turb.correct();
}
}

7
applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/createFluidFields.H
View File

@ -13,6 +13,7 @@
PtrList<volScalarField> dpdtFluid(fluidRegions.size());
List<scalar> initialMassFluid(fluidRegions.size());
List<bool> frozenFlowFluid(fluidRegions.size(), false);
PtrList<fv::IOoptionList> fluidFvOptions(fluidRegions.size());
@ -102,7 +103,7 @@
)
);
Info<< " Adding turbulence \n" << endl;
Info<< " Adding to turbulence\n" << endl;
turbulence.set
(
i,
@ -190,6 +191,10 @@
)
);
const dictionary& pimpleDict =
fluidRegions[i].solutionDict().subDict("PIMPLE");
pimpleDict.readIfPresent("frozenFlow", frozenFlowFluid[i]);
Info<< " Adding fvOptions\n" << endl;
fluidFvOptions.set
(

2
applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/setRegionFluidFields.H
View File

@ -28,3 +28,5 @@
dimMass,
initialMassFluid[i]
);
const bool frozenFlow = frozenFlowFluid[i];

33
applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/solveFluid.H
View File

@ -3,24 +3,31 @@ if (finalIter)
mesh.data::add("finalIteration", true);
}
if (oCorr == 0)
if (frozenFlow)
{
#include "rhoEqn.H"
#include "EEqn.H"
}
#include "UEqn.H"
#include "EEqn.H"
// --- PISO loop
for (int corr=0; corr<nCorr; corr++)
else
{
#include "pEqn.H"
if (oCorr == 0)
{
#include "rhoEqn.H"
}
#include "UEqn.H"
#include "EEqn.H"
// --- PISO loop
for (int corr=0; corr<nCorr; corr++)
{
#include "pEqn.H"
}
turb.correct();
rho = thermo.rho();
}
turb.correct();
rho = thermo.rho();
if (finalIter)
{
mesh.data::remove("finalIteration");