Commit 9efaf83f authored by Andrew Heather's avatar Andrew Heather
Browse files

ENH: updated buoyantBoussinesqSimpleFoam solver to include thermal wall functions

     on (kinematic) turbulent thermal conductivity, kappat and updated tutorial
parent 8f03a325
{
volScalarField kappaEff
(
"kappaEff",
turbulence->nu()/Pr + turbulence->nut()/Prt
);
kappat = turbulence->nut()/Prt;
kappat.correctBoundaryConditions();
volScalarField kappaEff("kappaEff", turbulence->nu()/Pr + kappat);
fvScalarMatrix TEqn
(
......
......@@ -64,6 +64,21 @@
1.0 - beta*(T - TRef)
);
// kinematic turbulent thermal thermal conductivity m2/s
Info<< "Reading field kappat\n" << endl;
volScalarField kappat
(
IOobject
(
"kappat",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh
);
Info<< "Calculating field g.h\n" << endl;
volScalarField gh("gh", g & mesh.C());
surfaceScalarField ghf("ghf", g & mesh.Cf());
......
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: dev |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object kappat;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -1 0 0 0 0];
internalField uniform 0;
boundaryField
{
floor
{
type kappatJayatillekeWallFunction;
Prt 0.85;
value uniform 0;
}
ceiling
{
type kappatJayatillekeWallFunction;
Prt 0.85;
value uniform 0;
}
fixedWalls
{
type kappatJayatillekeWallFunction;
Prt 0.85;
value uniform 0;
}
}
// ************************************************************************* //
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: dev |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object kappat;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -1 0 0 0 0];
internalField uniform 0;
boundaryField
{
ground
{
type kappatJayatillekeWallFunction;
Prt 0.85;
value uniform 0;
}
igloo_region0
{
type kappatJayatillekeWallFunction;
Prt 0.85;
value uniform 0;
}
twoFridgeFreezers_seal_0
{
type kappatJayatillekeWallFunction;
Prt 0.85;
value uniform 0;
}
twoFridgeFreezers_herring_1
{
type kappatJayatillekeWallFunction;
Prt 0.85;
value uniform 0;
}
}
// ************************************************************************* //
Supports Markdown
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