updates to enable requested field activated injection

This commit is contained in:
andy
2009-05-18 19:04:04 +01:00
parent 98fc11e827
commit a723b85abd
2 changed files with 12 additions and 23 deletions

View File

@ -2,7 +2,7 @@
========= | ========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | \\ / O peration |
\\ / A nd | Copyright (C) 1991-2007 OpenCFD Ltd. \\ / A nd | Copyright (C) 1991-2009 OpenCFD Ltd.
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
@ -76,6 +76,9 @@ int main(int argc, char *argv[])
Info<< "Time = " << runTime.timeName() << nl << endl; Info<< "Time = " << runTime.timeName() << nl << endl;
rhoEffLagrangian = coalParcels.rhoEff() + limestoneParcels.rhoEff();
pDyn = 0.5*rho*magSqr(U);
Info << "Evolving coal cloud" << endl; Info << "Evolving coal cloud" << endl;
coalParcels.evolve(); coalParcels.evolve();

View File

@ -28,50 +28,36 @@
thermo->rho() thermo->rho()
); );
// lagrangian coal density field // lagrangian effective density field - used externally (optional)
/* volScalarField rholc volScalarField rhoEffLagrangian
( (
IOobject IOobject
( (
"rholc", "rhoEffLagrangian",
runTime.timeName(), runTime.timeName(),
mesh, mesh,
IOobject::NO_READ, IOobject::NO_READ,
IOobject::AUTO_WRITE IOobject::AUTO_WRITE
), ),
mesh, mesh,
dimensionedScalar("zero", dimensionSet(1, -3, 0, 0, 0, 0, 0), 0.0) dimensionedScalar("zero", dimDensity, 0.0)
); );
// lagrangian limestone density field // dynamic pressure field - used externally (optional)
volScalarField rhols volScalarField pDyn
( (
IOobject IOobject
( (
"rhols", "injectionIndicator",
runTime.timeName(), runTime.timeName(),
mesh, mesh,
IOobject::NO_READ, IOobject::NO_READ,
IOobject::AUTO_WRITE IOobject::AUTO_WRITE
), ),
mesh, mesh,
dimensionedScalar("zero", dimensionSet(1, -3, 0, 0, 0, 0, 0), 0.0) dimensionedScalar("zero", dimPressure, 0.0)
); );
// lagrangian total density field
volScalarField rhol
(
IOobject
(
"rhol",
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
mesh,
dimensionedScalar("zero", dimensionSet(1, -3, 0, 0, 0, 0, 0), 0.0)
);*/
Info<< "\nReading field U\n" << endl; Info<< "\nReading field U\n" << endl;
volVectorField U volVectorField U