exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Changes for passing test-loop

Browse source
This commit is contained in:
Dominik Christ 2015-08-05 12:34:57 +01:00
parent 5e05134c73
commit 5e62f6e191
56 changed files with 434 additions and 280 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

36
applications/solvers/heatTransfer/chtMultiRegionFoam/derivedFvPatchFields/solidWallHeatFluxTemperature/solidWallHeatFluxTemperatureFvPatchScalarField.C
View file

@ -39,7 +39,7 @@ solidWallHeatFluxTemperatureFvPatchScalarField
:
fixedGradientFvPatchScalarField(p, iF),
q_(p.size(), 0.0),
KName_("undefined-K")
KappaName_("undefined-Kappa")
{}
@ -54,7 +54,7 @@ solidWallHeatFluxTemperatureFvPatchScalarField
:
fixedGradientFvPatchScalarField(ptf, p, iF, mapper),
q_(ptf.q_, mapper),
KName_(ptf.KName_)
KappaName_(ptf.KappaName_)
{}
@ -68,7 +68,7 @@ solidWallHeatFluxTemperatureFvPatchScalarField
:
fixedGradientFvPatchScalarField(p, iF, dict),
q_("q", dict, p.size()),
KName_(dict.lookup("K"))
KappaName_(dict.lookup("Kappa"))
{}
@ -80,7 +80,7 @@ solidWallHeatFluxTemperatureFvPatchScalarField
:
fixedGradientFvPatchScalarField(tppsf),
q_(tppsf.q_),
KName_(tppsf.KName_)
KappaName_(tppsf.KappaName_)
{}
@ -93,7 +93,7 @@ solidWallHeatFluxTemperatureFvPatchScalarField
:
fixedGradientFvPatchScalarField(tppsf, iF),
q_(tppsf.q_),
KName_(tppsf.KName_)
KappaName_(tppsf.KappaName_)
{}
@ -125,33 +125,33 @@ void Foam::solidWallHeatFluxTemperatureFvPatchScalarField::rmap
Foam::tmp<Foam::scalarField>
Foam::solidWallHeatFluxTemperatureFvPatchScalarField::K() const
Foam::solidWallHeatFluxTemperatureFvPatchScalarField::Kappa() const
{
const fvMesh& mesh = patch().boundaryMesh().mesh();
if (mesh.objectRegistry::foundObject<volScalarField>(KName_))
if (mesh.objectRegistry::foundObject<volScalarField>(KappaName_))
{
return lookupPatchField<volScalarField, scalar>(KName_);
return lookupPatchField<volScalarField, scalar>(KappaName_);
}
else if (mesh.objectRegistry::foundObject<volSymmTensorField>(KName_))
else if (mesh.objectRegistry::foundObject<volSymmTensorField>(KappaName_))
{
const symmTensorField& KWall =
lookupPatchField<volSymmTensorField, scalar>(KName_);
const symmTensorField& KappaWall =
lookupPatchField<volSymmTensorField, scalar>(KappaName_);
vectorField n = patch().nf();
return n & KWall & n;
return n & KappaWall & n;
}
else
{
FatalErrorIn
(
"solidWallHeatFluxTemperatureFvPatchScalarField::K()"
"solidWallHeatFluxTemperatureFvPatchScalarField::Kappa()"
" const"
) << "Did not find field " << KName_
) << "Did not find field " << KappaName_
<< " on mesh " << mesh.name() << " patch " << patch().name()
<< endl
<< "Please set 'K' to a valid volScalarField"
<< "Please set 'Kappa' to a valid volScalarField"
<< " or a valid volSymmTensorField." << exit(FatalError);
return scalarField(0);
@ -166,13 +166,13 @@ void Foam::solidWallHeatFluxTemperatureFvPatchScalarField::updateCoeffs()
return;
}
gradient() = q_/K();
gradient() = q_/Kappa();
fixedGradientFvPatchScalarField::updateCoeffs();
if (debug)
{
scalar Q = gSum(K()*patch().magSf()*snGrad());
scalar Q = gSum(Kappa()*patch().magSf()*snGrad());
Info<< patch().boundaryMesh().mesh().name() << ':'
<< patch().name() << ':'
@ -194,7 +194,7 @@ void Foam::solidWallHeatFluxTemperatureFvPatchScalarField::write
{
fixedGradientFvPatchScalarField::write(os);
q_.writeEntry("q", os);
os.writeKeyword("K") << KName_ << token::END_STATEMENT << nl;
os.writeKeyword("Kappa") << KappaName_ << token::END_STATEMENT << nl;
this->writeEntry("value", os);
}

8
applications/solvers/heatTransfer/chtMultiRegionFoam/derivedFvPatchFields/solidWallHeatFluxTemperature/solidWallHeatFluxTemperatureFvPatchScalarField.H
View file

@ -31,7 +31,7 @@ Description
myWallPatch
{
type solidWallHeatFluxTemperature;
K K; // Name of K field
Kappa Kappa; // Name of Kappa field
q uniform 1000; // Heat flux / [W/m2]
value uniform 300.0; // Initial temperature / [K]
gradient uniform 0.0; // Initial gradient / [K/m]
@ -67,7 +67,7 @@ class solidWallHeatFluxTemperatureFvPatchScalarField
scalarField q_;
//- Name of thermal conductivity field
word KName_;
word KappaName_;
public:
@ -143,8 +143,8 @@ public:
// Helper
//- Get K field on this patch
tmp<scalarField> K() const;
//- Get Kappa field on this patch
tmp<scalarField> Kappa() const;
// Evaluation functions

30
applications/solvers/heatTransfer/chtMultiRegionFoam/derivedFvPatchFields/solidWallMixedTemperatureCoupled/solidWallMixedTemperatureCoupledFvPatchScalarField.C
View file

@ -41,7 +41,7 @@ solidWallMixedTemperatureCoupledFvPatchScalarField
:
mixedFvPatchScalarField(p, iF),
neighbourFieldName_("undefined-neighbourFieldName"),
KName_("undefined-K")
KappaName_("undefined-Kappa")
{
this->refValue() = 0.0;
this->refGrad() = 0.0;
@ -60,7 +60,7 @@ solidWallMixedTemperatureCoupledFvPatchScalarField
:
mixedFvPatchScalarField(ptf, p, iF, mapper),
neighbourFieldName_(ptf.neighbourFieldName_),
KName_(ptf.KName_)
KappaName_(ptf.KappaName_)
{}
@ -74,7 +74,7 @@ solidWallMixedTemperatureCoupledFvPatchScalarField
:
mixedFvPatchScalarField(p, iF),
neighbourFieldName_(dict.lookup("neighbourFieldName")),
KName_(dict.lookup("K"))
KappaName_(dict.lookup("Kappa"))
{
if (!isA<directMappedPatchBase>(this->patch().patch()))
{
@ -123,16 +123,16 @@ solidWallMixedTemperatureCoupledFvPatchScalarField
:
mixedFvPatchScalarField(wtcsf, iF),
neighbourFieldName_(wtcsf.neighbourFieldName_),
KName_(wtcsf.KName_)
KappaName_(wtcsf.KappaName_)
{}
// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
const Foam::fvPatchScalarField&
Foam::solidWallMixedTemperatureCoupledFvPatchScalarField::K() const
Foam::solidWallMixedTemperatureCoupledFvPatchScalarField::Kappa() const
{
return this->lookupPatchField<volScalarField, scalar>(KName_);
return this->lookupPatchField<volScalarField, scalar>(KappaName_);
}
@ -181,8 +181,8 @@ void Foam::solidWallMixedTemperatureCoupledFvPatchScalarField::updateCoeffs()
nbrIntFld
);
// Swap to obtain full local values of neighbour K*delta
scalarField nbrKDelta = nbrField.K()*nbrPatch.deltaCoeffs();
// Swap to obtain full local values of neighbour Kappa*delta
scalarField nbrKappaDelta = nbrField.Kappa()*nbrPatch.deltaCoeffs();
mapDistribute::distribute
(
static_cast<Pstream::commsTypes>(Pstream::defaultCommsType()),
@ -190,14 +190,14 @@ void Foam::solidWallMixedTemperatureCoupledFvPatchScalarField::updateCoeffs()
distMap.constructSize(),
distMap.subMap(), // what to send
distMap.constructMap(), // what to receive
nbrKDelta
nbrKappaDelta
);
tmp<scalarField> myKDelta = K()*patch().deltaCoeffs();
tmp<scalarField> myKappaDelta = Kappa()*patch().deltaCoeffs();
// Both sides agree on
// - temperature : (myKDelta*fld + nbrKDelta*nbrFld)/(myKDelta+nbrKDelta)
// - temperature : (myKappaDelta*fld + nbrKappaDelta*nbrFld)/(myKappaDelta+nbrKappaDelta)
// - gradient : (temperature-fld)*delta
// We've got a degree of freedom in how to implement this in a mixed bc.
// (what gradient, what fixedValue and mixing coefficient)
@ -209,21 +209,21 @@ void Foam::solidWallMixedTemperatureCoupledFvPatchScalarField::updateCoeffs()
// same on both sides. This leads to the choice of
// - refGradient = zero gradient
// - refValue = neighbour value
// - mixFraction = nbrKDelta / (nbrKDelta + myKDelta())
// - mixFraction = nbrKappaDelta / (nbrKappaDelta + myKappaDelta())
this->refValue() = nbrIntFld;
this->refGrad() = 0.0;
this->valueFraction() = nbrKDelta / (nbrKDelta + myKDelta());
this->valueFraction() = nbrKappaDelta / (nbrKappaDelta + myKappaDelta());
mixedFvPatchScalarField::updateCoeffs();
if (debug)
{
scalar Q = gSum(K()*patch().magSf()*snGrad());
scalar Q = gSum(Kappa()*patch().magSf()*snGrad());
Info<< patch().boundaryMesh().mesh().name() << ':'
<< patch().name() << ':'
@ -249,7 +249,7 @@ void Foam::solidWallMixedTemperatureCoupledFvPatchScalarField::write
mixedFvPatchScalarField::write(os);
os.writeKeyword("neighbourFieldName")<< neighbourFieldName_
<< token::END_STATEMENT << nl;
os.writeKeyword("K") << KName_ << token::END_STATEMENT << nl;
os.writeKeyword("Kappa") << KappaName_ << token::END_STATEMENT << nl;
}

8
applications/solvers/heatTransfer/chtMultiRegionFoam/derivedFvPatchFields/solidWallMixedTemperatureCoupled/solidWallMixedTemperatureCoupledFvPatchScalarField.H
View file

@ -34,7 +34,7 @@ Description
{
type solidWallMixedTemperatureCoupled;
neighbourFieldName T;
K K;
Kappa Kappa;
value uniform 300;
}
@ -72,7 +72,7 @@ class solidWallMixedTemperatureCoupledFvPatchScalarField
const word neighbourFieldName_;
//- Name of thermal conductivity field
const word KName_;
const word KappaName_;
public:
@ -142,8 +142,8 @@ public:
// Member functions
//- Get corresponding K field
const fvPatchScalarField& K() const;
//- Get corresponding Kappa field
const fvPatchScalarField& Kappa() const;
//- Update the coefficients associated with the patch field
virtual void updateCoeffs();

8
applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/createFluidFields.H
View file

@ -1,7 +1,7 @@
// Initialise fluid field pointer lists
PtrList<basicPsiThermo> thermoFluid(fluidRegions.size());
PtrList<volScalarField> rhoFluid(fluidRegions.size());
PtrList<volScalarField> KFluid(fluidRegions.size());
PtrList<volScalarField> KappaFluid(fluidRegions.size());
PtrList<volVectorField> UFluid(fluidRegions.size());
PtrList<surfaceScalarField> phiFluid(fluidRegions.size());
PtrList<uniformDimensionedVectorField> gFluid(fluidRegions.size());
@ -41,15 +41,15 @@
)
);
Info<< " Adding to KFluid\n" << endl;
KFluid.set
Info<< " Adding to KappaFluid\n" << endl;
KappaFluid.set
(
i,
new volScalarField
(
IOobject
(
"K",
"Kappa",
runTime.timeName(),
fluidRegions[i],
IOobject::NO_READ,

2
applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/pEqn.H
View file

@ -71,5 +71,5 @@
}
// Update thermal conductivity
K = thermoFluid[i].Cp()*turb.alphaEff();
Kappa = thermoFluid[i].Cp()*turb.alphaEff();
}

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

@ -2,7 +2,7 @@
basicPsiThermo& thermo = thermoFluid[i];
volScalarField& rho = rhoFluid[i];
volScalarField& K = KFluid[i];
volScalarField& Kappa = KappaFluid[i];
volVectorField& U = UFluid[i];
surfaceScalarField& phi = phiFluid[i];
const dimensionedVector& g = gFluid[i];

6
applications/solvers/heatTransfer/chtMultiRegionFoam/solid/createSolidFields.H
View file

@ -2,7 +2,7 @@
PtrList<volScalarField> rhos(solidRegions.size());
PtrList<volScalarField> cps(solidRegions.size());
PtrList<volScalarField> rhosCps(solidRegions.size());
PtrList<volScalarField> Ks(solidRegions.size());
PtrList<volScalarField> Kappas(solidRegions.size());
PtrList<volScalarField> Ts(solidRegions.size());
// Populate solid field pointer lists
@ -53,8 +53,8 @@
new volScalarField("rhosCps", rhos[i]*cps[i])
);
Info<< " Adding to Ks\n" << endl;
Ks.set
Info<< " Adding to Kappas\n" << endl;
Kappas.set
(
i,
new volScalarField

2
applications/solvers/heatTransfer/chtMultiRegionFoam/solid/setRegionSolidFields.H
View file

@ -2,5 +2,5 @@
volScalarField& rho = rhos[i];
volScalarField& cp = cps[i];
volScalarField& K = Ks[i];
volScalarField& Kappa = Kappas[i];
volScalarField& T = Ts[i];

10
applications/solvers/heatTransfer/chtMultiRegionFoam/solid/solidRegionDiffNo.C
View file

@ -31,7 +31,7 @@ Foam::scalar Foam::solidRegionDiffNo
const fvMesh& mesh,
const Time& runTime,
const volScalarField& Cprho,
const volScalarField& K
const volScalarField& Kappa
)
{
scalar DiNum = 0.0;
@ -40,14 +40,14 @@ Foam::scalar Foam::solidRegionDiffNo
//- Can have fluid domains with 0 cells so do not test.
if (mesh.nInternalFaces())
{
surfaceScalarField KrhoCpbyDelta =
surfaceScalarField KappaRhoCpbyDelta =
mesh.surfaceInterpolation::deltaCoeffs()
* fvc::interpolate(K)
* fvc::interpolate(Kappa)
/ fvc::interpolate(Cprho);
DiNum = max(KrhoCpbyDelta.internalField())*runTime.deltaT().value();
DiNum = max(KappaRhoCpbyDelta.internalField())*runTime.deltaT().value();
meanDiNum = (average(KrhoCpbyDelta)).value()*runTime.deltaT().value();
meanDiNum = (average(KappaRhoCpbyDelta)).value()*runTime.deltaT().value();
}
Info<< "Region: " << mesh.name() << " Diffusion Number mean: " << meanDiNum

2
applications/solvers/heatTransfer/chtMultiRegionFoam/solid/solidRegionDiffNo.H
View file

@ -39,7 +39,7 @@ namespace Foam
const fvMesh& mesh,
const Time& runTime,
const volScalarField& Cprho,
const volScalarField& K
const volScalarField& Kappa
);
}

2
applications/solvers/heatTransfer/chtMultiRegionFoam/solid/solidRegionDiffusionNo.H
View file

@ -9,7 +9,7 @@
solidRegions[regionI],
runTime,
rhosCps[regionI],
Ks[regionI]
Kappas[regionI]
),
DiNum
);

2
applications/solvers/heatTransfer/chtMultiRegionFoam/solid/solveSolid.H
View file

@ -4,7 +4,7 @@
tmp<fvScalarMatrix> TEqn
(
fvm::ddt(rho*cp, T)
- fvm::laplacian(K, T)
- fvm::laplacian(Kappa, T)
);
TEqn().relax();
TEqn().solve();

18
applications/solvers/heatTransfer/chtMultiRegionSimpleFoam/derivedFvPatchFields/solidWallHeatFluxTemperature/solidWallHeatFluxTemperatureFvPatchScalarField.C
View file

@ -39,7 +39,7 @@ solidWallHeatFluxTemperatureFvPatchScalarField
:
fixedGradientFvPatchScalarField(p, iF),
q_(p.size(), 0.0),
KName_("undefined-K")
KappaName_("undefined-Kappa")
{}
@ -54,7 +54,7 @@ solidWallHeatFluxTemperatureFvPatchScalarField
:
fixedGradientFvPatchScalarField(ptf, p, iF, mapper),
q_(ptf.q_, mapper),
KName_(ptf.KName_)
KappaName_(ptf.KappaName_)
{}
@ -68,7 +68,7 @@ solidWallHeatFluxTemperatureFvPatchScalarField
:
fixedGradientFvPatchScalarField(p, iF, dict),
q_("q", dict, p.size()),
KName_(dict.lookup("K"))
KappaName_(dict.lookup("Kappa"))
{}
@ -80,7 +80,7 @@ solidWallHeatFluxTemperatureFvPatchScalarField
:
fixedGradientFvPatchScalarField(tppsf),
q_(tppsf.q_),
KName_(tppsf.KName_)
KappaName_(tppsf.KappaName_)
{}
@ -93,7 +93,7 @@ solidWallHeatFluxTemperatureFvPatchScalarField
:
fixedGradientFvPatchScalarField(tppsf, iF),
q_(tppsf.q_),
KName_(tppsf.KName_)
KappaName_(tppsf.KappaName_)
{}
@ -131,12 +131,12 @@ void Foam::solidWallHeatFluxTemperatureFvPatchScalarField::updateCoeffs()
return;
}
const scalarField& Kw = lookupPatchField<volScalarField, scalar>
const scalarField& Kappaw = lookupPatchField<volScalarField, scalar>
(
KName_
KappaName_
);
gradient() = q_/Kw;
gradient() = q_/Kappaw;
fixedGradientFvPatchScalarField::updateCoeffs();
}
@ -149,7 +149,7 @@ void Foam::solidWallHeatFluxTemperatureFvPatchScalarField::write
{
fixedGradientFvPatchScalarField::write(os);
q_.writeEntry("q", os);
os.writeKeyword("K") << KName_ << token::END_STATEMENT << nl;
os.writeKeyword("Kappa") << KappaName_ << token::END_STATEMENT << nl;
this->writeEntry("value", os);
}

4
applications/solvers/heatTransfer/chtMultiRegionSimpleFoam/derivedFvPatchFields/solidWallHeatFluxTemperature/solidWallHeatFluxTemperatureFvPatchScalarField.H
View file

@ -31,7 +31,7 @@ Description
myWallPatch
{
type solidWallHeatFluxTemperature;
K K; // Name of K field
Kappa Kappa; // Name of Kappa field
q uniform 1000; // Heat flux / [W/m2]
value 300.0; // Initial temperature / [K]
}
@ -66,7 +66,7 @@ class solidWallHeatFluxTemperatureFvPatchScalarField
scalarField q_;
//- Name of thermal conductivity field
word KName_;
word KappaName_;
public:

8
applications/solvers/heatTransfer/chtMultiRegionSimpleFoam/fluid/createFluidFields.H
View file

@ -1,7 +1,7 @@
// Initialise fluid field pointer lists
PtrList<basicPsiThermo> thermoFluid(fluidRegions.size());
PtrList<volScalarField> rhoFluid(fluidRegions.size());
PtrList<volScalarField> KFluid(fluidRegions.size());
PtrList<volScalarField> KappaFluid(fluidRegions.size());
PtrList<volVectorField> UFluid(fluidRegions.size());
PtrList<surfaceScalarField> phiFluid(fluidRegions.size());
PtrList<uniformDimensionedVectorField> gFluid(fluidRegions.size());
@ -45,15 +45,15 @@
)
);
Info<< " Adding to KFluid\n" << endl;
KFluid.set
Info<< " Adding to KappaFluid\n" << endl;
KappaFluid.set
(
i,
new volScalarField
(
IOobject
(
"K",
"Kappa",
runTime.timeName(),
fluidRegions[i],
IOobject::NO_READ,

2
applications/solvers/heatTransfer/chtMultiRegionSimpleFoam/fluid/pEqn.H
View file

@ -70,5 +70,5 @@
<< max(rho).value() << endl;
// Update thermal conductivity
K = thermo.Cp()*turb.alphaEff();
Kappa = thermo.Cp()*turb.alphaEff();
}

2
applications/solvers/heatTransfer/chtMultiRegionSimpleFoam/fluid/setRegionFluidFields.H
View file

@ -2,7 +2,7 @@
basicPsiThermo& thermo = thermoFluid[i];
volScalarField& rho = rhoFluid[i];
volScalarField& K = KFluid[i];
volScalarField& Kappa = KappaFluid[i];
volVectorField& U = UFluid[i];
surfaceScalarField& phi = phiFluid[i];
const dimensionedVector& g = gFluid[i];

8
applications/solvers/heatTransfer/chtMultiRegionSimpleFoam/solid/createSolidFields.H
View file

@ -2,7 +2,7 @@
PtrList<volScalarField> rhos(solidRegions.size());
PtrList<volScalarField> cps(solidRegions.size());
PtrList<volScalarField> rhosCps(solidRegions.size());
PtrList<volScalarField> Ks(solidRegions.size());
PtrList<volScalarField> Kappas(solidRegions.size());
PtrList<volScalarField> Ts(solidRegions.size());
// Populate solid field pointer lists
@ -53,15 +53,15 @@
new volScalarField("rhosCps", rhos[i]*cps[i])
);
Info<< " Adding to Ks\n" << endl;
Ks.set
Info<< " Adding to Kappas\n" << endl;
Kappas.set
(
i,
new volScalarField
(
IOobject
(
"K",
"Kappa",
runTime.timeName(),
solidRegions[i],
IOobject::MUST_READ,

2
applications/solvers/heatTransfer/chtMultiRegionSimpleFoam/solid/setRegionSolidFields.H
View file

@ -2,5 +2,5 @@
volScalarField& rho = rhos[i];
volScalarField& cp = cps[i];
volScalarField& K = Ks[i];
volScalarField& Kappa = Kappas[i];
volScalarField& T = Ts[i];

2
applications/solvers/heatTransfer/chtMultiRegionSimpleFoam/solid/solveSolid.H
View file

@ -3,7 +3,7 @@
{
fvScalarMatrix tEqn
(
-fvm::laplacian(K, T)
-fvm::laplacian(Kappa, T)
);
tEqn.relax();
eqnResidual = tEqn.solve().initialResidual();

2
etc/bashrc.mingw
View file

@ -71,7 +71,7 @@ add_to_path() {
# {{{ DEFINE PROCESS FUNCTIONS
setup_foam_env() {
: ${WM_OSTYPE:=POSIX}; export WM_OSTYPE
: ${WM_OSTYPE:=MSWindows}; export WM_OSTYPE
export export FLEX_DIR=$(echo $WD\.. | sed 's/\\/\//g' | sed 's/\(.*\):/\/\1/')
FOAM_ETC_DIR=$( cd "$( dirname "${BASH_SOURCE[0]}" )" && pwd )
export FOAM_INST_DIR=$(readlink -f $FOAM_ETC_DIR/../..)

4
src/finiteVolume/finiteVolume/divSchemes/gaussDivScheme/vectorGaussDivScheme.H
View file

@ -50,13 +50,13 @@ namespace fv
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
template<> inline
template<>
tmp<BlockLduSystem<vector, scalar> > gaussDivScheme<vector>::fvmUDiv
(
const GeometricField<vector, fvPatchField, volMesh>&
) const;
template<> inline
template<>
tmp<BlockLduSystem<vector, scalar> > gaussDivScheme<vector>::fvmUDiv
(
const surfaceScalarField& flux,

46
src/turbulenceModels/compressible/RAS/derivedFvPatchFields/turbulentTemperatureCoupledBaffle/turbulentTemperatureCoupledBaffleFvPatchScalarField.C
View file

@ -112,7 +112,7 @@ turbulentTemperatureCoupledBaffleFvPatchScalarField
:
fixedValueFvPatchScalarField(p, iF),
neighbourFieldName_("undefined-neighbourFieldName"),
KName_("undefined-K")
KappaName_("undefined-Kappa")
{}
@ -127,7 +127,7 @@ turbulentTemperatureCoupledBaffleFvPatchScalarField
:
fixedValueFvPatchScalarField(ptf, p, iF, mapper),
neighbourFieldName_(ptf.neighbourFieldName_),
KName_(ptf.KName_)
KappaName_(ptf.KappaName_)
{}
@ -141,7 +141,7 @@ turbulentTemperatureCoupledBaffleFvPatchScalarField
:
fixedValueFvPatchScalarField(p, iF, dict),
neighbourFieldName_(dict.lookup("neighbourFieldName")),
KName_(dict.lookup("K"))
KappaName_(dict.lookup("Kappa"))
{
if (!isA<directMappedPatchBase>(this->patch().patch()))
{
@ -173,18 +173,18 @@ turbulentTemperatureCoupledBaffleFvPatchScalarField
:
fixedValueFvPatchScalarField(wtcsf, iF),
neighbourFieldName_(wtcsf.neighbourFieldName_),
KName_(wtcsf.KName_)
KappaName_(wtcsf.KappaName_)
{}
// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
Foam::tmp<Foam::scalarField>
Foam::turbulentTemperatureCoupledBaffleFvPatchScalarField::K() const
Foam::turbulentTemperatureCoupledBaffleFvPatchScalarField::Kappa() const
{
const fvMesh& mesh = patch().boundaryMesh().mesh();
if (KName_ == "none")
if (KappaName_ == "none")
{
const compressible::RASModel& model =
db().lookupObject<compressible::RASModel>("RASProperties");
@ -198,28 +198,28 @@ Foam::turbulentTemperatureCoupledBaffleFvPatchScalarField::K() const
talpha().boundaryField()[patch().index()]
*thermo.Cp()().boundaryField()[patch().index()];
}
else if (mesh.objectRegistry::foundObject<volScalarField>(KName_))
else if (mesh.objectRegistry::foundObject<volScalarField>(KappaName_))
{
return lookupPatchField<volScalarField, scalar>(KName_);
return lookupPatchField<volScalarField, scalar>(KappaName_);
}
else if (mesh.objectRegistry::foundObject<volSymmTensorField>(KName_))
else if (mesh.objectRegistry::foundObject<volSymmTensorField>(KappaName_))
{
const symmTensorField& KWall =
lookupPatchField<volSymmTensorField, scalar>(KName_);
const symmTensorField& KappaWall =
lookupPatchField<volSymmTensorField, scalar>(KappaName_);
vectorField n = patch().nf();
return n & KWall & n;
return n & KappaWall & n;
}
else
{
FatalErrorIn
(
"turbulentTemperatureCoupledBaffleFvPatchScalarField::K() const"
) << "Did not find field " << KName_
"turbulentTemperatureCoupledBaffleFvPatchScalarField::Kappa() const"
) << "Did not find field " << KappaName_
<< " on mesh " << mesh.name() << " patch " << patch().name()
<< endl
<< "Please set 'K' to 'none', a valid volScalarField"
<< "Please set 'Kappa' to 'none', a valid volScalarField"
<< " or a valid volSymmTensorField." << exit(FatalError);
return scalarField(0);
@ -282,8 +282,8 @@ void Foam::turbulentTemperatureCoupledBaffleFvPatchScalarField::updateCoeffs()
nbrIntFld
);
// Swap to obtain full local values of neighbour K*delta
scalarField nbrKDelta = nbrField.K()*nbrPatch.deltaCoeffs();
// Swap to obtain full local values of neighbour Kappa*delta
scalarField nbrKappaDelta = nbrField.Kappa()*nbrPatch.deltaCoeffs();
mapDistribute::distribute
(
static_cast<Pstream::commsTypes>(Pstream::defaultCommsType()),
@ -291,16 +291,16 @@ void Foam::turbulentTemperatureCoupledBaffleFvPatchScalarField::updateCoeffs()
distMap.constructSize(),
distMap.subMap(), // what to send
distMap.constructMap(), // what to receive
nbrKDelta
nbrKappaDelta
);
tmp<scalarField> myKDelta = K()*patch().deltaCoeffs();
tmp<scalarField> myKappaDelta = Kappa()*patch().deltaCoeffs();
// Calculate common wall temperature. Reuse *this to store common value.
scalarField Twall
(
(myKDelta()*intFld() + nbrKDelta*nbrIntFld)
/ (myKDelta() + nbrKDelta)
(myKappaDelta()*intFld() + nbrKappaDelta*nbrIntFld)
/ (myKappaDelta() + nbrKappaDelta)
);
// Assign to me
fvPatchScalarField::operator=(Twall);
@ -326,7 +326,7 @@ void Foam::turbulentTemperatureCoupledBaffleFvPatchScalarField::updateCoeffs()
// (*this-intFld())
// * patch().deltaCoeffs();
scalar Q = gSum(K()*patch().magSf()*snGrad());
scalar Q = gSum(Kappa()*patch().magSf()*snGrad());
Info<< patch().boundaryMesh().mesh().name() << ':'
<< patch().name() << ':'
@ -354,7 +354,7 @@ void Foam::turbulentTemperatureCoupledBaffleFvPatchScalarField::write
fixedValueFvPatchScalarField::write(os);
os.writeKeyword("neighbourFieldName")<< neighbourFieldName_
<< token::END_STATEMENT << nl;
os.writeKeyword("K") << KName_ << token::END_STATEMENT << nl;
os.writeKeyword("Kappa") << KappaName_ << token::END_STATEMENT << nl;
}

16
src/turbulenceModels/compressible/RAS/derivedFvPatchFields/turbulentTemperatureCoupledBaffle/turbulentTemperatureCoupledBaffleFvPatchScalarField.H
View file

@ -28,23 +28,23 @@ Description
Harmonic fixed value boundary condition for temperature, to be used
for heat-transfer on back-to-back baffles.
If my temperature is T1, heat conductivity K1 and neighbour is T2,K2
If my temperature is T1, heat conductivity Kappa1 and neighbour is T2,Kappa2
both sides get fixedValue (K1/dx1*T1 + K2/dx2*T2)/(K1/dx1+K2/dx2)
both sides get fixedValue (Kappa1/dx1*T1 + Kappa2/dx2*T2)/(Kappa1/dx1+Kappa2/dx2)
Example usage:
myInterfacePatchName
{
type compressible::turbulentTemperatureCoupledBaffle;
neighbourFieldName T;
K K; // or none
Kappa Kappa; // or none
value uniform 300;
}
Needs to be on underlying directMapped(Wall)FvPatch.
Note: if K is "none" looks up RASModel and basicThermo, otherwise expects
the solver to calculate a 'K' field.
Note: if Kappa is "none" looks up RASModel and basicThermo, otherwise expects
the solver to calculate a 'Kappa' field.
Note: runs in parallel with arbitrary decomposition. Uses directMapped
functionality to calculate exchange.
@ -84,7 +84,7 @@ class turbulentTemperatureCoupledBaffleFvPatchScalarField
const word neighbourFieldName_;
//- Name of thermal conductivity field
const word KName_;
const word KappaName_;
// Private Member Functions
@ -161,8 +161,8 @@ public:
// Member functions
//- Get corresponding K field
tmp<scalarField> K() const;
//- Get corresponding Kappa field
tmp<scalarField> Kappa() const;
//- Update the coefficients associated with the patch field
virtual void updateCoeffs();

48
src/turbulenceModels/compressible/RAS/derivedFvPatchFields/turbulentTemperatureCoupledBaffleMixed/turbulentTemperatureCoupledBaffleMixedFvPatchScalarField.C
View file

@ -53,7 +53,7 @@ turbulentTemperatureCoupledBaffleMixedFvPatchScalarField
:
mixedFvPatchScalarField(p, iF),
neighbourFieldName_("undefined-neighbourFieldName"),
KName_("undefined-K")
KappaName_("undefined-Kappa")
{
this->refValue() = 0.0;
this->refGrad() = 0.0;
@ -72,7 +72,7 @@ turbulentTemperatureCoupledBaffleMixedFvPatchScalarField
:
mixedFvPatchScalarField(ptf, p, iF, mapper),
neighbourFieldName_(ptf.neighbourFieldName_),
KName_(ptf.KName_)
KappaName_(ptf.KappaName_)
{}
@ -86,7 +86,7 @@ turbulentTemperatureCoupledBaffleMixedFvPatchScalarField
:
mixedFvPatchScalarField(p, iF),
neighbourFieldName_(dict.lookup("neighbourFieldName")),
KName_(dict.lookup("K"))
KappaName_(dict.lookup("Kappa"))
{
if (!isA<directMappedPatchBase>(this->patch().patch()))
{
@ -135,18 +135,18 @@ turbulentTemperatureCoupledBaffleMixedFvPatchScalarField
:
mixedFvPatchScalarField(wtcsf, iF),
neighbourFieldName_(wtcsf.neighbourFieldName_),
KName_(wtcsf.KName_)
KappaName_(wtcsf.KappaName_)
{}
// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
tmp<scalarField>
turbulentTemperatureCoupledBaffleMixedFvPatchScalarField::K() const
turbulentTemperatureCoupledBaffleMixedFvPatchScalarField::Kappa() const
{
const fvMesh& mesh = patch().boundaryMesh().mesh();
if (KName_ == "none")
if (KappaName_ == "none")
{
const compressible::RASModel& model =
db().lookupObject<compressible::RASModel>("RASProperties");
@ -160,29 +160,29 @@ turbulentTemperatureCoupledBaffleMixedFvPatchScalarField::K() const
talpha().boundaryField()[patch().index()]
*thermo.Cp()().boundaryField()[patch().index()];
}
else if (mesh.objectRegistry::foundObject<volScalarField>(KName_))
else if (mesh.objectRegistry::foundObject<volScalarField>(KappaName_))
{
return lookupPatchField<volScalarField, scalar>(KName_);
return lookupPatchField<volScalarField, scalar>(KappaName_);
}
else if (mesh.objectRegistry::foundObject<volSymmTensorField>(KName_))
else if (mesh.objectRegistry::foundObject<volSymmTensorField>(KappaName_))
{
const symmTensorField& KWall =
lookupPatchField<volSymmTensorField, scalar>(KName_);
const symmTensorField& KappaWall =
lookupPatchField<volSymmTensorField, scalar>(KappaName_);
vectorField n = patch().nf();
return n & KWall & n;
return n & KappaWall & n;
}
else
{
FatalErrorIn
(
"turbulentTemperatureCoupledBaffleMixedFvPatchScalarField::K()"
"turbulentTemperatureCoupledBaffleMixedFvPatchScalarField::Kappa()"
" const"
) << "Did not find field " << KName_
) << "Did not find field " << KappaName_
<< " on mesh " << mesh.name() << " patch " << patch().name()
<< endl
<< "Please set 'K' to 'none', a valid volScalarField"
<< "Please set 'Kappa' to 'none', a valid volScalarField"
<< " or a valid volSymmTensorField." << exit(FatalError);
return scalarField(0);
@ -241,8 +241,8 @@ void turbulentTemperatureCoupledBaffleMixedFvPatchScalarField::updateCoeffs()
nbrIntFld
);
// Swap to obtain full local values of neighbour K*delta
scalarField nbrKDelta = nbrField.K()*nbrPatch.deltaCoeffs();
// Swap to obtain full local values of neighbour Kappa*delta
scalarField nbrKappaDelta = nbrField.Kappa()*nbrPatch.deltaCoeffs();
mapDistribute::distribute
(
static_cast<Pstream::commsTypes>(Pstream::defaultCommsType()),
@ -250,14 +250,14 @@ void turbulentTemperatureCoupledBaffleMixedFvPatchScalarField::updateCoeffs()
distMap.constructSize(),
distMap.subMap(), // what to send
distMap.constructMap(), // what to receive
nbrKDelta
nbrKappaDelta
);
tmp<scalarField> myKDelta = K()*patch().deltaCoeffs();
tmp<scalarField> myKappaDelta = Kappa()*patch().deltaCoeffs();
// Both sides agree on
// - temperature : (myKDelta*fld + nbrKDelta*nbrFld)/(myKDelta+nbrKDelta)
// - temperature : (myKappaDelta*fld + nbrKappaDelta*nbrFld)/(myKappaDelta+nbrKappaDelta)
// - gradient : (temperature-fld)*delta
// We've got a degree of freedom in how to implement this in a mixed bc.
// (what gradient, what fixedValue and mixing coefficient)
@ -269,21 +269,21 @@ void turbulentTemperatureCoupledBaffleMixedFvPatchScalarField::updateCoeffs()
// same on both sides. This leads to the choice of
// - refGradient = zero gradient
// - refValue = neighbour value
// - mixFraction = nbrKDelta / (nbrKDelta + myKDelta())
// - mixFraction = nbrKappaDelta / (nbrKappaDelta + myKappaDelta())
this->refValue() = nbrIntFld;
this->refGrad() = 0.0;
this->valueFraction() = nbrKDelta / (nbrKDelta + myKDelta());
this->valueFraction() = nbrKappaDelta / (nbrKappaDelta + myKappaDelta());
mixedFvPatchScalarField::updateCoeffs();
if (debug)
{
scalar Q = gSum(K()*patch().magSf()*snGrad());
scalar Q = gSum(Kappa()*patch().magSf()*snGrad());
Info<< patch().boundaryMesh().mesh().name() << ':'
<< patch().name() << ':'
@ -309,7 +309,7 @@ void turbulentTemperatureCoupledBaffleMixedFvPatchScalarField::write
mixedFvPatchScalarField::write(os);
os.writeKeyword("neighbourFieldName")<< neighbourFieldName_
<< token::END_STATEMENT << nl;
os.writeKeyword("K") << KName_ << token::END_STATEMENT << nl;
os.writeKeyword("Kappa") << KappaName_ << token::END_STATEMENT << nl;
}

12
src/turbulenceModels/compressible/RAS/derivedFvPatchFields/turbulentTemperatureCoupledBaffleMixed/turbulentTemperatureCoupledBaffleMixedFvPatchScalarField.H
View file

@ -38,14 +38,14 @@ Description
{
type compressible::turbulentTemperatureCoupledBaffleMixed;
neighbourFieldName T;
K K; // or none
Kappa Kappa; // or none
value uniform 300;
}
Needs to be on underlying directMapped(Wall)FvPatch.
Note: if K is "none" looks up RASModel and basicThermo, otherwise expects
the solver to calculate a 'K' field.
Note: if Kappa is "none" looks up RASModel and basicThermo, otherwise expects
the solver to calculate a 'Kappa' field.
Note: runs in parallel with arbitrary decomposition. Uses directMapped
functionality to calculate exchange.
@ -87,7 +87,7 @@ class turbulentTemperatureCoupledBaffleMixedFvPatchScalarField
const word neighbourFieldName_;
//- Name of thermal conductivity field
const word KName_;
const word KappaName_;
public:
@ -162,8 +162,8 @@ public:
// Member functions
//- Get corresponding K field
tmp<scalarField> K() const;
//- Get corresponding Kappa field
tmp<scalarField> Kappa() const;
//- Update the coefficients associated with the patch field
virtual void updateCoeffs();

15
tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/Allrun
View file

@ -94,11 +94,16 @@ cp 0/* 0.001/heater/
cp 0/* 0.001/leftSolid/
cp 0/* 0.001/rightSolid/
cp 0/* 0.001/topAir/
changeDictionary -region bottomAir
changeDictionary -region topAir
changeDictionary -region heater
changeDictionary -region leftSolid
changeDictionary -region rightSolid
runApplication changeDictionary -region bottomAir
mv log.changeDictionary log.changeDictionary_bottomAir
runApplication changeDictionary -region topAir
mv log.changeDictionary log.changeDictionary_topAir
runApplication changeDictionary -region heater
mv log.changeDictionary log.changeDictionary_heater
runApplication changeDictionary -region leftSolid
mv log.changeDictionary log.changeDictionary_leftSolid
runApplication changeDictionary -region rightSolid
mv log.changeDictionary log.changeDictionary_rightSolid
runApplicationAndReportOnError chtMultiRegionFoam

6
tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/bottomAir/changeDictionaryDict
View file

@ -65,7 +65,7 @@ dictionaryReplacement
neighbourRegionName leftSolid;
neighbourPatchName leftSolid_to_bottomAir;
neighbourFieldName T;
K K;
Kappa Kappa;
value uniform 300;
}
@ -75,7 +75,7 @@ dictionaryReplacement
neighbourRegionName heater;
neighbourPatchName heater_to_bottomAir;
neighbourFieldName T;
K K;
Kappa Kappa;
value uniform 300;
}
bottomAir_to_rightSolid
@ -84,7 +84,7 @@ dictionaryReplacement
neighbourRegionName rightSolid;
neighbourPatchName rightSolid_to_bottomAir;
neighbourFieldName T;
K K;
Kappa Kappa;
value uniform 300;
}
}

10
tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/heater/changeDictionaryDict
View file

@ -43,7 +43,7 @@ dictionaryReplacement
neighbourRegionName bottomAir;
neighbourPatchName bottomAir_to_heater;
neighbourFieldName T;
K K;
Kappa Kappa;
value uniform 300;
}
@ -53,7 +53,7 @@ dictionaryReplacement
neighbourRegionName leftSolid;
neighbourPatchName leftSolid_to_heater;
neighbourFieldName T;
K K;
Kappa Kappa;
value uniform 300;
}
heater_to_rightSolid
@ -62,7 +62,7 @@ dictionaryReplacement
neighbourRegionName rightSolid;
neighbourPatchName rightSolid_to_heater;
neighbourFieldName T;
K K;
Kappa Kappa;
value uniform 300;
}
heater_to_topAir
@ -71,7 +71,7 @@ dictionaryReplacement
neighbourRegionName topAir;
neighbourPatchName topAir_to_heater;
neighbourFieldName T;
K K;
Kappa Kappa;
value uniform 300;
}
}
@ -114,7 +114,7 @@ dictionaryReplacement
}
}
K
Kappa
{
internalField uniform 80;

2
tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/heater/fvSchemes
View file

@ -32,7 +32,7 @@ divSchemes
laplacianSchemes
{
default none;
laplacian(K,T) Gauss linear limited 0.333;
laplacian(Kappa,T) Gauss linear limited 0.333;
}
interpolationSchemes

8
tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/leftSolid/changeDictionaryDict
View file

@ -40,7 +40,7 @@ dictionaryReplacement
neighbourRegionName bottomAir;
neighbourPatchName bottomAir_to_leftSolid;
neighbourFieldName T;
K K;
Kappa Kappa;
value uniform 300;
}
@ -50,7 +50,7 @@ dictionaryReplacement
neighbourRegionName heater;
neighbourPatchName heater_to_leftSolid;
neighbourFieldName T;
K K;
Kappa Kappa;
value uniform 300;
}
leftSolid_to_topAir
@ -59,7 +59,7 @@ dictionaryReplacement
neighbourRegionName topAir;
neighbourPatchName topAir_to_leftSolid;
neighbourFieldName T;
K K;
Kappa Kappa;
value uniform 300;
}
}
@ -98,7 +98,7 @@ dictionaryReplacement
}
}
K
Kappa
{
internalField uniform 80;

8
tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/rightSolid/changeDictionaryDict
View file

@ -40,7 +40,7 @@ dictionaryReplacement
neighbourRegionName heater;
neighbourPatchName heater_to_rightSolid;
neighbourFieldName T;
K K;
Kappa Kappa;
value uniform 300;
}
rightSolid_to_bottomAir
@ -49,7 +49,7 @@ dictionaryReplacement
neighbourRegionName bottomAir;
neighbourPatchName bottomAir_to_rightSolid;
neighbourFieldName T;
K K;
Kappa Kappa;
value uniform 300;
}
rightSolid_to_topAir
@ -58,7 +58,7 @@ dictionaryReplacement
neighbourRegionName topAir;
neighbourPatchName topAir_to_rightSolid;
neighbourFieldName T;
K K;
Kappa Kappa;
value uniform 300;
}
}
@ -97,7 +97,7 @@ dictionaryReplacement
}
}
K
Kappa
{
internalField uniform 80;

6
tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater/system/topAir/changeDictionaryDict
View file

@ -66,7 +66,7 @@ dictionaryReplacement
neighbourRegionName leftSolid;
neighbourPatchName leftSolid_to_topAir;
neighbourFieldName T;
K K;
Kappa Kappa;
value uniform 300;
}
@ -76,7 +76,7 @@ dictionaryReplacement
neighbourRegionName heater;
neighbourPatchName heater_to_topAir;
neighbourFieldName T;
K K;
Kappa Kappa;
value uniform 300;
}
topAir_to_rightSolid
@ -85,7 +85,7 @@ dictionaryReplacement
neighbourRegionName rightSolid;
neighbourPatchName rightSolid_to_topAir;
neighbourFieldName T;
K K;
Kappa Kappa;
value uniform 300;
}
}

240
tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/constant/polyMesh/boundary
View file

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 3.2 |
| \\ / O peration | Version: 3.1 |
| \\ / A nd | Web: http://www.foam-extend.org |
| \\/ M anipulation | |
| \\/ M anipulation | For copyright notice see file Copyright |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -15,43 +15,253 @@ FoamFile
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
6
41
(
maxY
{
type wall;
nFaces 300;
startFace 8300;
nFaces 1200;
startFace 41630;
}
minX
{
type patch;
nFaces 100;
startFace 8600;
nFaces 250;
startFace 42830;
}
maxX
{
type patch;
nFaces 100;
startFace 8700;
nFaces 250;
startFace 43080;
}
minY
{
type wall;
nFaces 300;
startFace 8800;
nFaces 1200;
startFace 43330;
}
minZ
{
type wall;
nFaces 300;
startFace 9100;
nFaces 750;
startFace 44530;
}
maxZ
{
type wall;
nFaces 300;
startFace 9400;
nFaces 750;
startFace 45280;
}
bottomAir_minX
{
type wall;
nFaces 0;
startFace 46030;
}
bottomAir_maxX
{
type wall;
nFaces 0;
startFace 46030;
}
bottomAir_minY
{
type wall;
nFaces 0;
startFace 46030;
}
bottomAir_minZ
{
type wall;
nFaces 0;
startFace 46030;
}
bottomAir_maxZ
{
type wall;
nFaces 0;
startFace 46030;
}
bottomAir_bottomAir_to_rightSolid
{
type wall;
nFaces 0;
startFace 46030;
}
bottomAir_bottomAir_to_leftSolid
{
type wall;
nFaces 0;
startFace 46030;
}
bottomAir_bottomAir_to_heater
{
type wall;
nFaces 0;
startFace 46030;
}
topAir_maxY
{
type wall;
nFaces 0;
startFace 46030;
}
topAir_minX
{
type wall;
nFaces 0;
startFace 46030;
}
topAir_maxX
{
type wall;
nFaces 0;
startFace 46030;
}
topAir_minZ
{
type wall;
nFaces 0;
startFace 46030;
}
topAir_maxZ
{
type wall;
nFaces 0;
startFace 46030;
}
topAir_topAir_to_rightSolid
{
type wall;
nFaces 0;
startFace 46030;
}
topAir_topAir_to_heater
{
type wall;
nFaces 0;
startFace 46030;
}
topAir_topAir_to_leftSolid
{
type wall;
nFaces 0;
startFace 46030;
}
leftSolid_minX
{
type wall;
nFaces 0;
startFace 46030;
}
leftSolid_minZ
{
type wall;
nFaces 0;
startFace 46030;
}
leftSolid_maxZ
{
type wall;
nFaces 0;
startFace 46030;
}
leftSolid_leftSolid_to_bottomAir
{
type wall;
nFaces 0;
startFace 46030;
}
leftSolid_leftSolid_to_heater
{
type wall;
nFaces 0;
startFace 46030;
}
leftSolid_leftSolid_to_topAir
{
type wall;
nFaces 0;
startFace 46030;
}
rightSolid_maxX
{
type wall;
nFaces 0;
startFace 46030;
}
rightSolid_minZ
{
type wall;
nFaces 0;
startFace 46030;
}
rightSolid_maxZ
{
type wall;
nFaces 0;
startFace 46030;
}
rightSolid_rightSolid_to_bottomAir
{
type wall;
nFaces 0;
startFace 46030;
}
rightSolid_rightSolid_to_topAir
{
type wall;
nFaces 0;
startFace 46030;
}
rightSolid_rightSolid_to_heater
{
type wall;
nFaces 0;
startFace 46030;
}
heater_minY
{
type wall;
nFaces 0;
startFace 46030;
}
heater_minZ
{
type wall;
nFaces 0;
startFace 46030;
}
heater_maxZ
{
type wall;
nFaces 0;
startFace 46030;
}
heater_heater_to_topAir
{
type wall;
nFaces 0;
startFace 46030;
}
heater_heater_to_rightSolid
{
type wall;
nFaces 0;
startFace 46030;
}
heater_heater_to_leftSolid
{
type wall;
nFaces 0;
startFace 46030;
}
heater_heater_to_bottomAir
{
type wall;
nFaces 0;
startFace 46030;
}
)

6
tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/system/bottomAir/changeDictionaryDict
View file

@ -102,21 +102,21 @@ dictionaryReplacement
{
type solidWallMixedTemperatureCoupled;
neighbourFieldName T;
K K;
Kappa Kappa;
value $internalField;
}
bottomAir_to_heater
{
type solidWallMixedTemperatureCoupled;
neighbourFieldName T;
K K;
Kappa Kappa;
value $internalField;
}
bottomAir_to_rightSolid
{
type solidWallMixedTemperatureCoupled;
neighbourFieldName T;
K K;
Kappa Kappa;
value $internalField;
}
}

10
tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/system/heater/changeDictionaryDict
View file

@ -43,7 +43,7 @@ dictionaryReplacement
{
type solidWallMixedTemperatureCoupled;
neighbourFieldName T;
K K;
Kappa Kappa;
value $internalField;
}
@ -51,21 +51,21 @@ dictionaryReplacement
{
type solidWallMixedTemperatureCoupled;
neighbourFieldName T;
K K;
Kappa Kappa;
value $internalField;
}
heater_to_rightSolid
{
type solidWallMixedTemperatureCoupled;
neighbourFieldName T;
K K;
Kappa Kappa;
value $internalField;
}
heater_to_topAir
{
type solidWallMixedTemperatureCoupled;
neighbourFieldName T;
K K;
Kappa Kappa;
value $internalField;
}
}
@ -115,7 +115,7 @@ dictionaryReplacement
}
}
K
Kappa
{
internalField uniform 80;

2
tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/system/heater/fvSchemes
View file

@ -32,7 +32,7 @@ divSchemes
laplacianSchemes
{
default none;
laplacian(K,T) Gauss linear limited 0.333;
laplacian(Kappa,T) Gauss linear limited 0.333;
}
interpolationSchemes

8
tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/system/leftSolid/changeDictionaryDict
View file

@ -41,7 +41,7 @@ dictionaryReplacement
{
type solidWallMixedTemperatureCoupled;
neighbourFieldName T;
K K;
Kappa Kappa;
value $internalField;
}
@ -49,14 +49,14 @@ dictionaryReplacement
{
type solidWallMixedTemperatureCoupled;
neighbourFieldName T;
K K;
Kappa Kappa;
value $internalField;
}
leftSolid_to_topAir
{
type solidWallMixedTemperatureCoupled;
neighbourFieldName T;
K K;
Kappa Kappa;
value $internalField;
}
}
@ -101,7 +101,7 @@ dictionaryReplacement
}
}
K
Kappa
{
internalField uniform 80;

2
tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/system/leftSolid/fvSchemes
View file

@ -32,7 +32,7 @@ divSchemes
laplacianSchemes
{
default none;
laplacian(K,T) Gauss linear limited 0.333;
laplacian(Kappa,T) Gauss linear limited 0.333;
}
interpolationSchemes

8
tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/system/rightSolid/changeDictionaryDict
View file

@ -41,21 +41,21 @@ dictionaryReplacement
{
type solidWallMixedTemperatureCoupled;
neighbourFieldName T;
K K;
Kappa Kappa;
value $internalField;
}
rightSolid_to_bottomAir
{
type solidWallMixedTemperatureCoupled;
neighbourFieldName T;
K K;
Kappa Kappa;
value $internalField;
}
rightSolid_to_topAir
{
type solidWallMixedTemperatureCoupled;
neighbourFieldName T;
K K;
Kappa Kappa;
value $internalField;
}
}
@ -100,7 +100,7 @@ dictionaryReplacement
}
}
K
Kappa
{
internalField uniform 80;

2
tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/system/rightSolid/fvSchemes
View file

@ -32,7 +32,7 @@ divSchemes
laplacianSchemes
{
default none;
laplacian(K,T) Gauss linear limited 0.333;
laplacian(Kappa,T) Gauss linear limited 0.333;
}
interpolationSchemes

6
tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeater/system/topAir/changeDictionaryDict
View file

@ -103,7 +103,7 @@ dictionaryReplacement
{
type solidWallMixedTemperatureCoupled;
neighbourFieldName T;
K K;
Kappa Kappa;
value $internalField;
}
@ -111,14 +111,14 @@ dictionaryReplacement
{
type solidWallMixedTemperatureCoupled;
neighbourFieldName T;
K K;
Kappa Kappa;
value $internalField;
}
topAir_to_rightSolid
{
type solidWallMixedTemperatureCoupled;
neighbourFieldName T;
K K;
Kappa Kappa;
value $internalField;
}
}

2
tutorials/heatTransfer/chtMultiRegionSimpleFoam/multiRegionHeater/Allrun
View file

@ -24,7 +24,7 @@ done
# remove solid fields from fluid regions (important for post-processing)
for i in bottomAir topAir
do
rm -f 0*/$i/{cp,K,rho}
rm -f 0*/$i/{cp,Kappa,rho}
done
for i in bottomAir topAir heater leftSolid rightSolid

2
tutorials/heatTransfer/chtMultiRegionSimpleFoam/multiRegionHeater/system/bottomAir/changeDictionaryDict
View file

@ -57,7 +57,7 @@ dictionaryReplacement
{
type compressible::turbulentTemperatureCoupledBaffle;
neighbourFieldName T;
K K;
Kappa Kappa;
value uniform 300;
}
}

4
tutorials/heatTransfer/chtMultiRegionSimpleFoam/multiRegionHeater/system/heater/changeDictionaryDict
View file

@ -47,7 +47,7 @@ dictionaryReplacement
{
type compressible::turbulentTemperatureCoupledBaffle;
neighbourFieldName T;
K K;
Kappa Kappa;
value uniform 300;
}
minY
@ -72,7 +72,7 @@ dictionaryReplacement
}
}
K
Kappa
{
internalField uniform 80;

2
tutorials/heatTransfer/chtMultiRegionSimpleFoam/multiRegionHeater/system/heater/fvSchemes
View file

@ -32,7 +32,7 @@ divSchemes
laplacianSchemes
{
default none;
laplacian(K,T) Gauss linear limited 0.333;
laplacian(Kappa,T) Gauss linear limited 0.333;
}
interpolationSchemes

4
tutorials/heatTransfer/chtMultiRegionSimpleFoam/multiRegionHeater/system/leftSolid/changeDictionaryDict
View file

@ -43,7 +43,7 @@ dictionaryReplacement
{
type compressible::turbulentTemperatureCoupledBaffle;
neighbourFieldName T;
K K;
Kappa Kappa;
value uniform 300;
}
}
@ -63,7 +63,7 @@ dictionaryReplacement
}
}
K
Kappa
{
internalField uniform 80;

4
tutorials/heatTransfer/chtMultiRegionSimpleFoam/multiRegionHeater/system/rightSolid/changeDictionaryDict
View file

@ -43,7 +43,7 @@ dictionaryReplacement
{
type compressible::turbulentTemperatureCoupledBaffle;
neighbourFieldName T;
K K;
Kappa Kappa;
value uniform 300;
}
}
@ -63,7 +63,7 @@ dictionaryReplacement
}
}
K
Kappa
{
internalField uniform 80;

2
tutorials/heatTransfer/chtMultiRegionSimpleFoam/multiRegionHeater/system/topAir/changeDictionaryDict
View file

@ -68,7 +68,7 @@ dictionaryReplacement
{
type compressible::turbulentTemperatureCoupledBaffle;
neighbourFieldName T;
K K;
Kappa Kappa;
value uniform 300;
}
}

1
tutorials/immersedBoundary/backwardStepShortTurbulentSimpleIbFoam/0_org/U
View file

@ -31,6 +31,7 @@ boundaryField
setDeadCellValue yes;
deadCellValue (0 0 0);
value uniform (0 0 0);
}
inlet

17
tutorials/immersedBoundary/backwardStepShortTurbulentSimpleIbFoam/system/controlDict
View file

@ -44,23 +44,6 @@ timePrecision 6;
runTimeModifiable yes;
libs
(
"libfoam.so"
"libincompressibleTransportModels.so"
"libcompressibleTurbulenceModel.so"
"libincompressibleRASModels.so"
// "libimmersedBoundary.so"
// "libimmersedBoundaryTurbulence.so"
"libtriSurface.so"
"libmeshTools.so"
"libsampling.so"
"libsurfMesh.so"
"libimmersedBoundary.so"
"libimmersedBoundaryTurbulence.so"
);
functions
(
forces

10
tutorials/immersedBoundary/cylinderInChannelFineIcoIbFoam/system/controlDict
View file

@ -46,16 +46,6 @@ timePrecision 6;
runTimeModifiable yes;
libs
(
"libfoam.so"
"libtriSurface.so"
"libmeshTools.so"
"libsampling.so"
"libsurfMesh.so"
"libimmersedBoundary.so"
);
functions
(
forces

10
tutorials/immersedBoundary/cylinderInChannelIcoIbFoam/system/controlDict
View file

@ -50,16 +50,6 @@ adjustTimeStep no;
maxCo 0.2;
// libs
// (
// "libfoam.so"
// "libtriSurface.so"
// "libmeshTools.so"
// "libsampling.so"
// "libsurfMesh.so"
// "libimmersedBoundary.so"
// );
functions
(
forces

10
tutorials/immersedBoundary/cylinderInChannelPotential/system/controlDict
View file

@ -48,16 +48,6 @@ adjustTimeStep no;
maxCo 0.2;
libs
(
"libfoam.so"
"libtriSurface.so"
"libmeshTools.so"
"libsampling.so"
"libsurfMesh.so"
"libimmersedBoundary.so"
);
// functions
// (
// immersedBoundaryForce

15
tutorials/immersedBoundary/pitzDailyTurbulentSimpleIbFoam/system/controlDict
View file

@ -44,21 +44,6 @@ timePrecision 6;
runTimeModifiable yes;
libs
(
"libfoam.so"
"libincompressibleTransportModels.so"
"libcompressibleTurbulenceModel.so"
"libincompressibleRASModels.so"
"libtriSurface.so"
"libmeshTools.so"
"libsampling.so"
"libsurfMesh.so"
"libimmersedBoundary.so"
"libimmersedBoundaryTurbulence.so"
);
functions
(
forces