Commit 939a46c8 authored by henry's avatar henry
Browse files

Added pd, changed p BCs to calculated and added pRef to thermophysicalProperties

parent 1202a245
......@@ -58,7 +58,7 @@
Info<< "Calculating field g.h\n" << endl;
volScalarField gh("gh", g & mesh.C());
dimensionedScalar pRef("pRef", p.dimensions(), 1.0e5);
dimensionedScalar pRef("pRef", p.dimensions(), thermo->lookup("pRef"));
Info<< "Creating field pd\n" << endl;
volScalarField pd
......
......@@ -54,8 +54,7 @@
Info<< "Calculating field g.h\n" << endl;
volScalarField gh("gh", g & mesh.C());
dimensionedScalar pRef("pRef", p.dimensions(), 1.0e5);
dimensionedScalar pRef("pRef", p.dimensions(), thermo->lookup("pRef"));
Info<< "Creating field pd\n" << endl;
volScalarField pd
......@@ -77,7 +76,13 @@
label pdRefCell = 0;
scalar pdRefValue = 0.0;
setRefCell(pd, mesh.solutionDict().subDict("SIMPLE"), pdRefCell, pdRefValue);
setRefCell
(
pd,
mesh.solutionDict().subDict("SIMPLE"),
pdRefCell,
pdRefValue
);
dimensionedScalar initialMass = fvc::domainIntegrate(rho);
......@@ -55,8 +55,7 @@
Info<< "Calculating field g.h\n" << endl;
volScalarField gh("gh", g & mesh.C());
dimensionedScalar pRef("pRef", p.dimensions(), 1.0e5);
dimensionedScalar pRef("pRef", p.dimensions(), thermo->lookup("pRef"));
Info<< "Creating field pd\n" << endl;
volScalarField pd
......
......@@ -11,7 +11,12 @@
List<scalar> initialMassf(fluidRegions.size());
dimensionedScalar pRef("pRef", dimensionSet(1, -1, -2, 0, 0), 1.0E5);
dimensionedScalar pRef
(
"pRef",
dimensionSet(1, -1, -2, 0, 0),
rp.lookup("pRef")
);
// Populate fluid field pointer lists
forAll(fluidRegions, i)
......
......@@ -53,7 +53,7 @@
Info<< "Calculating field g.h\n" << endl;
volScalarField gh("gh", g & mesh.C());
dimensionedScalar pRef("pRef", p.dimensions(), 0.0);
dimensionedScalar pRef("pRef", p.dimensions(), thermo->lookup("pRef"));
Info<< "Creating field pd\n" << endl;
volScalarField pd
......
......@@ -20,21 +20,21 @@ internalField uniform 100000;
boundaryField
{
floor
floor
{
type fixedFluxBuoyantPressure;
type calculated;
value uniform 100000;
}
ceiling
ceiling
{
type fixedFluxBuoyantPressure;
type calculated;
value uniform 100000;
}
fixedWalls
fixedWalls
{
type fixedFluxBuoyantPressure;
type calculated;
value uniform 100000;
}
}
......
......@@ -18,4 +18,6 @@ thermoType hThermo<pureMixture<constTransport<specieThermo<hConstThermo<per
mixture air 1 28.9 1000 0 1.8e-05 0.7;
pRef 1e5;
// ************************************************************************* //
......@@ -16,25 +16,25 @@ FoamFile
dimensions [1 -1 -2 0 0 0 0];
internalField uniform 100000;
internalField uniform 0;
boundaryField
{
floor
floor
{
type fixedFluxBuoyantPressure;
type calculated;
value uniform 100000;
}
ceiling
ceiling
{
type fixedFluxBuoyantPressure;
type calculated;
value uniform 100000;
}
fixedWalls
fixedWalls
{
type fixedFluxBuoyantPressure;
type calculated;
value uniform 100000;
}
}
......
......@@ -18,4 +18,6 @@ thermoType hThermo<pureMixture<constTransport<specieThermo<hConstThermo<per
mixture air 1 28.9 1000 0 1.8e-05 0.7;
pRef 1e5;
// ************************************************************************* //
......@@ -20,27 +20,27 @@ internalField uniform 100000;
boundaryField
{
floor
floor
{
type fixedFluxBuoyantPressure;
type calculated;
value uniform 100000;
}
ceiling
ceiling
{
type fixedFluxBuoyantPressure;
type calculated;
value uniform 100000;
}
fixedWalls
fixedWalls
{
type fixedFluxBuoyantPressure;
type calculated;
value uniform 100000;
}
box
{
type fixedFluxBuoyantPressure;
type calculated;
value uniform 100000;
}
}
......
......@@ -18,4 +18,6 @@ thermoType hThermo<pureMixture<constTransport<specieThermo<hConstThermo<per
mixture air 1 28.9 1000 0 1.8e-05 0.7;
pRef 1e5;
// ************************************************************************* //
Markdown is supported
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