Skip to content
Snippets Groups Projects
Commit eaf50306 authored by Andrew Heather's avatar Andrew Heather
Browse files

ENH: applyBoundaryLayer - determine compressibility automatically by the...

ENH: applyBoundaryLayer - determine compressibility automatically by the presence of a thermophysicalProperties file
parent e35afa35
No related merge requests found
...@@ -3,7 +3,7 @@ ...@@ -3,7 +3,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) 2011-2015 OpenFOAM Foundation \\ / A nd | Copyright (C) 2011-2015 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation | Copyright (C) 2015 OpenCFD Ltd.
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
This file is part of OpenFOAM. This file is part of OpenFOAM.
...@@ -33,6 +33,10 @@ Description ...@@ -33,6 +33,10 @@ Description
the thickness coefficient supplied via the option -Cbl. If both options the thickness coefficient supplied via the option -Cbl. If both options
are provided -ybl is used. are provided -ybl is used.
Compressible modes is automatically selected based on the existence of the
"thermophysicalProperties" dictionary required to construct the
thermodynamics package.
\*---------------------------------------------------------------------------*/ \*---------------------------------------------------------------------------*/
#include "fvCFD.H" #include "fvCFD.H"
...@@ -331,11 +335,6 @@ int main(int argc, char *argv[]) ...@@ -331,11 +335,6 @@ int main(int argc, char *argv[])
"scalar", "scalar",
"boundary-layer thickness as Cbl * mean distance to wall" "boundary-layer thickness as Cbl * mean distance to wall"
); );
argList::addBoolOption
(
"compressible",
"apply to compressible case"
);
#include "setRootCase.H" #include "setRootCase.H"
...@@ -360,8 +359,6 @@ int main(int argc, char *argv[]) ...@@ -360,8 +359,6 @@ int main(int argc, char *argv[])
#include "createNamedMesh.H" #include "createNamedMesh.H"
#include "createFields.H" #include "createFields.H"
const bool compressible = args.optionFound("compressible");
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
// Modify velocity by applying a 1/7th power law boundary-layer // Modify velocity by applying a 1/7th power law boundary-layer
...@@ -384,7 +381,15 @@ int main(int argc, char *argv[]) ...@@ -384,7 +381,15 @@ int main(int argc, char *argv[])
U.write(); U.write();
if (compressible) if
(
IOobject
(
basicThermo::dictName,
runTime.constant(),
mesh
).headerOk()
)
{ {
calcCompressible(mesh, mask, U, y, ybl); calcCompressible(mesh, mask, U, y, ybl);
} }
......
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