mirror of
https://develop.openfoam.com/Development/openfoam.git
synced 2025-12-28 03:37:59 +00:00
345a2a42f1
- use succincter method names that more closely resemble dictionary
and HashTable method names. This improves method name consistency
between classes and also requires less typing effort:
args.found(optName) vs. args.optionFound(optName)
args.readIfPresent(..) vs. args.optionReadIfPresent(..)
...
args.opt<scalar>(optName) vs. args.optionRead<scalar>(optName)
args.read<scalar>(index) vs. args.argRead<scalar>(index)
- the older method names forms have been retained for code compatibility,
but are now deprecated
168 lines
3.3 KiB
C
168 lines
3.3 KiB
C
#include "readGravitationalAcceleration.H"
|
|
|
|
word continuousPhaseName
|
|
(
|
|
IOdictionary
|
|
(
|
|
IOobject
|
|
(
|
|
"transportProperties",
|
|
runTime.constant(),
|
|
mesh,
|
|
IOobject::MUST_READ
|
|
)
|
|
).lookup("continuousPhase")
|
|
);
|
|
|
|
Info<< "Reading field U\n" << endl;
|
|
volVectorField Uc
|
|
(
|
|
IOobject
|
|
(
|
|
IOobject::groupName("U", continuousPhaseName),
|
|
runTime.timeName(),
|
|
mesh,
|
|
IOobject::MUST_READ,
|
|
IOobject::AUTO_WRITE
|
|
),
|
|
mesh
|
|
);
|
|
|
|
Info<< "Reading field p\n" << endl;
|
|
volScalarField p
|
|
(
|
|
IOobject
|
|
(
|
|
"p",
|
|
runTime.timeName(),
|
|
mesh,
|
|
IOobject::MUST_READ,
|
|
IOobject::AUTO_WRITE
|
|
),
|
|
mesh
|
|
);
|
|
|
|
|
|
Info<< "Reading/calculating continuous-phase face flux field phic\n"
|
|
<< endl;
|
|
|
|
surfaceScalarField phic
|
|
(
|
|
IOobject
|
|
(
|
|
IOobject::groupName("phi", continuousPhaseName),
|
|
runTime.timeName(),
|
|
mesh,
|
|
IOobject::READ_IF_PRESENT,
|
|
IOobject::AUTO_WRITE
|
|
),
|
|
linearInterpolate(Uc) & mesh.Sf()
|
|
);
|
|
|
|
label pRefCell = 0;
|
|
scalar pRefValue = 0.0;
|
|
setRefCell(p, pimple.dict(), pRefCell, pRefValue);
|
|
mesh.setFluxRequired(p.name());
|
|
|
|
Info<< "Creating turbulence model\n" << endl;
|
|
|
|
singlePhaseTransportModel continuousPhaseTransport(Uc, phic);
|
|
|
|
dimensionedScalar rhocValue
|
|
(
|
|
IOobject::groupName("rho", continuousPhaseName),
|
|
dimDensity,
|
|
continuousPhaseTransport.lookup
|
|
(
|
|
IOobject::groupName("rho", continuousPhaseName)
|
|
)
|
|
);
|
|
|
|
volScalarField rhoc
|
|
(
|
|
IOobject
|
|
(
|
|
rhocValue.name(),
|
|
runTime.timeName(),
|
|
mesh,
|
|
IOobject::NO_READ,
|
|
IOobject::AUTO_WRITE
|
|
),
|
|
mesh,
|
|
rhocValue
|
|
);
|
|
|
|
volScalarField muc
|
|
(
|
|
IOobject
|
|
(
|
|
IOobject::groupName("mu", continuousPhaseName),
|
|
runTime.timeName(),
|
|
mesh,
|
|
IOobject::NO_READ,
|
|
IOobject::AUTO_WRITE
|
|
),
|
|
rhoc*continuousPhaseTransport.nu()
|
|
);
|
|
|
|
Info << "Creating field alphac\n" << endl;
|
|
// alphac must be constructed before the cloud
|
|
// so that the drag-models can find it
|
|
volScalarField alphac
|
|
(
|
|
IOobject
|
|
(
|
|
IOobject::groupName("alpha", continuousPhaseName),
|
|
runTime.timeName(),
|
|
mesh,
|
|
IOobject::READ_IF_PRESENT,
|
|
IOobject::AUTO_WRITE
|
|
),
|
|
mesh,
|
|
dimensionedScalar("0", dimless, 0)
|
|
);
|
|
|
|
word kinematicCloudName("kinematicCloud");
|
|
args.readIfPresent("cloud", kinematicCloudName);
|
|
|
|
Info<< "Constructing kinematicCloud " << kinematicCloudName << endl;
|
|
basicKinematicTypeCloud kinematicCloud
|
|
(
|
|
kinematicCloudName,
|
|
rhoc,
|
|
Uc,
|
|
muc,
|
|
g
|
|
);
|
|
|
|
// Particle fraction upper limit
|
|
scalar alphacMin
|
|
(
|
|
1.0
|
|
- readScalar
|
|
(
|
|
kinematicCloud.particleProperties().subDict("constantProperties")
|
|
.lookup("alphaMax")
|
|
)
|
|
);
|
|
|
|
// Update alphac from the particle locations
|
|
alphac = max(1.0 - kinematicCloud.theta(), alphacMin);
|
|
alphac.correctBoundaryConditions();
|
|
|
|
surfaceScalarField alphacf("alphacf", fvc::interpolate(alphac));
|
|
surfaceScalarField alphaPhic("alphaPhic", alphacf*phic);
|
|
|
|
autoPtr<PhaseIncompressibleTurbulenceModel<singlePhaseTransportModel>>
|
|
continuousPhaseTurbulence
|
|
(
|
|
PhaseIncompressibleTurbulenceModel<singlePhaseTransportModel>::New
|
|
(
|
|
alphac,
|
|
Uc,
|
|
alphaPhic,
|
|
phic,
|
|
continuousPhaseTransport
|
|
)
|
|
);
|