 |
 |
 |
|
|
|
[Sponsors] | ||||
November 24, 2021, 16:07
|
Solvers in OpenFOAM for LES + heat transfer
|
#1 |
|
Member
Arun subramanian
Join Date: Jun 2016
Location: Florence,Italy
Posts: 48
Rep Power: 10  |
Hello!
I am trying to run a film cooling simulation in openfoam using LES. All the tutorials I have seen so far are cold flows without changes in temperature. I would like to know if there is a 'sample test case' for forced convection problems. I have some nice results with RANS but I am struggling with LES. Is there a solver that works best for these types of problems? In the case below, I have some velocity field in the flow but no changes in the temperature field with time... I am pretty sure that I am missing something blatantly wrong. Would be nice if someone could point it out please! Code:
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 7
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
location "0";
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (10 0 0);
boundaryField
{
Body_5.Inlet
{
type surfaceNormalFixedValue;
refValue uniform -10;
}
Body_5.Inlet_coolant
{
type surfaceNormalFixedValue;
refValue uniform -10.68323;
}
Body_5.Basewall
{
type noSlip;
}
Body_5.side2
{
type noSlip;
}
Body_5.Top
{
type noSlip;
}
Body_5.Side1
{
type noSlip;
}
Body_5.Outlet
{
type inletOutlet;
inletValue uniform (0 0 0);
value uniform (0 0 0);
}
Body_5.Wall_coolant
{
type noSlip;
}
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 7
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 1 0 0 0];
internalField uniform 293;
boundaryField
{
Body_5.Inlet
{
type fixedValue;
value uniform 293;
}
Body_5.Inlet_coolant
{
type fixedValue;
value uniform 313;
}
Body_5.Basewall
{
type zeroGradient;
}
Body_5.side2
{
type zeroGradient;
}
Body_5.Top
{
type zeroGradient;
}
Body_5.Side1
{
type zeroGradient;
}
Body_5.Outlet
{
type inletOutlet;
value $internalField;
inletValue $internalField;
}
Body_5.Wall_coolant
{
type zeroGradient;
}
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 7
\\/ M anipulation |
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 7
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -1 -2 0 0 0 0];
internalField uniform 101325;
boundaryField
{
Body_5.Inlet
{
type zeroGradient;
// type mixed;
refValue uniform 101325;
refGradient uniform 0;
valueFraction uniform 0.3;
}
Body_5.Inlet_coolant
{
type zeroGradient;
// type mixed;
refValue uniform 101325;
refGradient uniform 0;
valueFraction uniform 0.3;
}
Body_5.Basewall
{
type zeroGradient;
}
Body_5.side2
{
type zeroGradient;
}
Body_5.Top
{
type zeroGradient;
}
Body_5.Side1
{
type zeroGradient;
}
Body_5.Outlet
{
type fixedValue;
value uniform 101325;
// type mixed;
// refValue uniform 110000;
// refGradient uniform 0;
// valueFraction uniform 1;
// type transonicOutletPressure;
// U U;
// phi phi;
// gamma 1.4;
// psi psi;
// pInf uniform 110000;
}
Body_5.Wall_coolant
{
type zeroGradient;
}
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 7
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -1 0 0 0 0];
internalField uniform 0;
boundaryField
{
Body_5.Inlet
{
type calculated;
value uniform 0;
}
Body_5.Inlet_coolant
{
type calculated;
value uniform 0;
}
Body_5.Basewall
{
type nutkWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0;
}
Body_5.side2
{
type nutkWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0;
}
Body_5.Top
{
type nutkWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0;
}
Body_5.Side1
{
type nutkWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0;
}
Body_5.Outlet
{
type calculated;
value uniform 0;
}
Body_5.Wall_coolant
{
type nutkWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0;
}
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 7
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object alphat;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -1 -1 0 0 0 0];
internalField uniform 0;
boundaryField
{
Body_5.Inlet
{
type calculated;
value uniform 0;
}
Body_5.Inlet_coolant
{
type calculated;
value uniform 0;
}
Body_5.Basewall
{
type compressible::alphatWallFunction;
Prt 0.85;
value uniform 0;
}
Body_5.side2
{
type compressible::alphatWallFunction;
Prt 0.85;
value uniform 0;
}
Body_5.Top
{
type compressible::alphatWallFunction;
Prt 0.85;
value uniform 0;
}
Body_5.Side1
{
type compressible::alphatWallFunction;
Prt 0.85;
value uniform 0;
}
Body_5.Outlet
{
type calculated;
value uniform 0;
}
Body_5.Wall_coolant
{
type compressible::alphatWallFunction;
Prt 0.85;
value uniform 0;
}
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 3.0.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object turbulenceProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
simulationType LES;
LES
{
LESModel Smagorinsky;
turbulence on;
printCoeffs on;
delta cubeRootVol;
cubeRootVolCoeffs
{
deltaCoeff 1;
}
PrandtlCoeffs
{
delta cubeRootVol;
cubeRootVolCoeffs
{
deltaCoeff 1;
}
smoothCoeffs
{
delta cubeRootVol;
cubeRootVolCoeffs
{
deltaCoeff 1;
}
maxDeltaRatio 1.1;
}
Cdelta 0.158;
}
vanDriestCoeffs
{
delta cubeRootVol;
cubeRootVolCoeffs
{
deltaCoeff 1;
}
smoothCoeffs
{
delta cubeRootVol;
cubeRootVolCoeffs
{
deltaCoeff 1;
}
maxDeltaRatio 1.1;
}
Aplus 26;
Cdelta 0.158;
}
smoothCoeffs
{
delta cubeRootVol;
cubeRootVolCoeffs
{
deltaCoeff 1;
}
maxDeltaRatio 1.1;
}
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 3.0.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object thermophysicalProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
thermoType
{
type hePsiThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState perfectGas;
specie specie;
energy sensibleInternalEnergy;
}
mixture
{
specie
{
nMoles 1;
molWeight 28.9;
}
thermodynamics
{
Cp 1007;
Hf 0;
}
transport
{
mu 1.8e-05; //1.2 kg/m^3 x 1.5e-05 m^2/s
Pr 1;
}
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 3.0.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers
{
"(p|rho)"
{
solver PCG;
preconditioner DIC;
tolerance 1e-6;
relTol 1e-4;
}
"(p|rho)Final"
{
$p;
relTol 0;
}
"(U|k|nuTilda)"
{
solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-6;
relTol 1e-6;
}
"(U|k|nuTilda)Final"
{
$U;
relTol 0;
}
e
{
solver smoothSolver;
smoother symGaussSeidel;
tolerance 10;
relTol 1;
}
eFinal
{
solver smoothSolver;
smoother symGaussSeidel;
tolerance 10;
relTol 1;
}
}
PIMPLE
{
momentumPredictor yes;
nOuterCorrectors 2;
nCorrectors 2;
nNonOrthogonalCorrectors 1;
rhoMin 1.2;
rhoMax 1.2;
}
relaxationFactors
{
equations
{
".*" 1;
}
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 3.0.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
ddtSchemes
{
default backward;
}
gradSchemes
{
default cellLimited Gauss linear 1;
}
divSchemes
{
default none;
div(phi,U) Gauss filteredLinear;
div(phi,e) Gauss linear;
div(phi,K) Gauss linear;
div(phiv,p) Gauss filteredLinear;
//div(phi,k) Gauss linear;
div(phi,B) Gauss linear;
div(phi,muTilda)Gauss linear;
div(B) Gauss linear;
div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}
laplacianSchemes
{
default Gauss linear corrected;
}
interpolationSchemes
{
default linear;
}
snGradSchemes
{
default corrected;
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 3.0.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
application rhoPimpleFoam;
//startFrom startTime;
startFrom latestTime;
startTime 0;
//stopAt writeNow;
stopAt endTime;
endTime 5000000;
deltaT 1e-05;
writeControl runTime;
writeInterval 4e-04;
purgeWrite 10;
writeFormat binary;
writePrecision 6;
writeCompression off;
timeFormat general;
timePrecision 6;
graphFormat raw;
runTimeModifiable true;
libs ("libLEMOS-2.4.x.so");
functions
{
/*probes
{
type probes;
functionObjectLibs ("libsampling.so");
enabled true;
outputControl timeStep;
outputInterval 1;
probeLocations
(
(0.0738 0.0015 0.0)
(0.0738 -0.0015 0.0)
);
fields
(
U
);
} */
fieldAverage
{
type fieldAverage;
functionObjectLibs ("libfieldFunctionObjects.so");
enabled true;
outputControl outputTime;
fields
(
U
{
mean on;
prime2Mean on;
base time;
}
nut
{
mean on;
prime2Mean off;
base time;
}
);
}
}
// ************************************************************************* //
|
|
|
|
|
|
November 26, 2021, 08:57
|
|
#2 |
|
New Member
JT
Join Date: Nov 2021
Posts: 1
Rep Power: 0 |
Check my former lab's validation cases: https://github.uconn.edu/CTF-UConn/chtValidationCases
Last edited by jt_turbo; November 26, 2021 at 08:58. Reason: na |
|
|
|
|
|
 |
| Tags |
| film cooling, heat transfer, les, openfoam 7, rhopimplefoam |
|
|
Similar Threads
|
||||
| Thread | Thread Starter | Forum | Replies | Last Post |
| Map of the OpenFOAM Forum - Understanding where to post your questions! | wyldckat | OpenFOAM | 10 | September 2, 2021 06:29 |
| Interphase mass transfer of a reaction | cfx_ws1992 | Main CFD Forum | 0 | May 15, 2017 22:42 |
| Question about heat transfer coefficient setting for CFX | Anna Tian | CFX | 1 | June 16, 2013 07:28 |
| Conductive Heat Transfer in OpenFoam? | FloK | OpenFOAM | 8 | March 18, 2011 14:51 |
| Convective Heat Transfer - Heat Exchanger | Mark | CFX | 6 | November 15, 2004 16:55 |
Linear Mode
