zhyang-dev/OpenFOAM-12
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

tutorials/modules/CHT/wallBoiling: corrected wall thickness

Browse Source 
Patch contributed by Juho Peltola, VTT.
This commit is contained in:
Henry Weller
2022-11-18 15:18:09 +00:00
parent b36d8fca44
commit d2378ef06e
12 changed files with 14 additions and 397 deletions
Download Patch File Download Diff File Expand all files Collapse all files

48
tutorials/modules/CHT/wallBoiling/0/fluid/dDeparture.liquid
View File

@ -1,48 +0,0 @@
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: dev
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
format ascii;
class volScalarField;
object dDeparture.liquid;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 0 0 0 0 0];
internalField uniform 0;
boundaryField
{
inlet
{
type calculated;
value uniform 0;
}
outlet
{
type calculated;
value uniform 0;
}
front
{
type wedge;
}
back
{
type wedge;
}
wall
{
type calculated;
value uniform 1e-05;
}
}
// ************************************************************************* //

48
tutorials/modules/CHT/wallBoiling/0/fluid/fDeparture.liquid
View File

@ -1,48 +0,0 @@
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: dev
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
format ascii;
class volScalarField;
object fDeparture.liquid;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 -1 0 0 0 0];
internalField uniform 0;
boundaryField
{
inlet
{
type calculated;
value uniform 0;
}
outlet
{
type calculated;
value uniform 0;
}
front
{
type wedge;
}
back
{
type wedge;
}
wall
{
type calculated;
value uniform 0;
}
}
// ************************************************************************* //

48
tutorials/modules/CHT/wallBoiling/0/fluid/nucleationSiteDensity.liquid
View File

@ -1,48 +0,0 @@
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: dev
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
format ascii;
class volScalarField;
object nucleationSiteDensity.liquid;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 -2 0 0 0 0 0];
internalField uniform 0;
boundaryField
{
inlet
{
type calculated;
value uniform 0;
}
outlet
{
type calculated;
value uniform 0;
}
front
{
type wedge;
}
back
{
type wedge;
}
wall
{
type calculated;
value uniform 0;
}
}
// ************************************************************************* //

48
tutorials/modules/CHT/wallBoiling/0/fluid/qEvaporative.liquid
View File

@ -1,48 +0,0 @@
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: dev
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
format ascii;
class volScalarField;
object qEvaporative.liquid;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 0 -3 0 0 0 0];
internalField uniform 0;
boundaryField
{
inlet
{
type calculated;
value uniform 0;
}
outlet
{
type calculated;
value uniform 0;
}
front
{
type wedge;
}
back
{
type wedge;
}
wall
{
type calculated;
value uniform 0;
}
}
// ************************************************************************* //

48
tutorials/modules/CHT/wallBoiling/0/fluid/qQuenching.liquid
View File

@ -1,48 +0,0 @@
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: dev
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
format ascii;
class volScalarField;
object qQuenching.liquid;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 0 -3 0 0 0 0];
internalField uniform 0;
boundaryField
{
inlet
{
type calculated;
value uniform 0;
}
outlet
{
type calculated;
value uniform 0;
}
front
{
type wedge;
}
back
{
type wedge;
}
wall
{
type calculated;
value uniform 0;
}
}
// ************************************************************************* //

48
tutorials/modules/CHT/wallBoiling/0/fluid/wetFraction.liquid
View File

@ -1,48 +0,0 @@
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: dev
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
format ascii;
class volScalarField;
object wetFraction.liquid;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 0 0 0 0];
internalField uniform 0;
boundaryField
{
inlet
{
type calculated;
value uniform 0;
}
outlet
{
type calculated;
value uniform 0;
}
front
{
type wedge;
}
back
{
type wedge;
}
wall
{
type calculated;
value uniform 0;
}
}
// ************************************************************************* //

2
tutorials/modules/CHT/wallBoiling/0/solid/T
View File

@ -28,7 +28,7 @@ boundaryField
wall_outer wall_outer
{ {
type externalTemperature; type externalTemperature;
q uniform 61150.34; q uniform 66919.2;
value $internalField; value $internalField;
} }
top top

2
tutorials/modules/CHT/wallBoiling/system/blockMeshDict
View File

@ -17,7 +17,7 @@ backgroundMesh
{ {
length 3.5; length 3.5;
rA 0.0096; rA 0.0096;
rB 0.0116; rB 0.0106;
lengthCells 350; lengthCells 350;
rAcells 40; rAcells 40;
rBcells 10; rBcells 10;

3
tutorials/modules/CHT/wallBoiling/system/controlDict.orig → tutorials/modules/CHT/wallBoiling/system/controlDict
View File

@ -28,7 +28,7 @@ startTime 0;
stopAt endTime; stopAt endTime;
endTime 10; endTime 8;
deltaT 0.0001; deltaT 0.0001;
@ -69,6 +69,7 @@ functions
type wallBoilingProperties; type wallBoilingProperties;
functionObjectLibs ( "libmultiphaseEulerFoamFunctionObjects.so" ); functionObjectLibs ( "libmultiphaseEulerFoamFunctionObjects.so" );
writeControl writeTime; writeControl writeTime;
executeAtStart no;
region fluid; region fluid;
phase liquid; phase liquid;
} }

30
tutorials/modules/CHT/wallBoiling/system/fvSchemes
View File

@ -9,66 +9,36 @@ FoamFile
{ {
format ascii; format ascii;
class dictionary; class dictionary;
location "system";
object fvSchemes; object fvSchemes;
} }
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
ddtSchemes ddtSchemes
{ {
default Euler;
} }
gradSchemes gradSchemes
{ {
default Gauss linear;
limited cellLimited Gauss linear 1;
} }
divSchemes divSchemes
{ {
default none;
div(phi,alpha) Gauss vanLeer;
div(phir,alpha) Gauss vanLeer;
div(alphaRhoPhi,U) Gauss linearUpwind limited;
div(phi,U) Gauss linearUpwind limited;
div(alphaRhoPhi,Yi) Gauss linearUpwind limited;
"div\(alphaRhoPhi,(h|e)\)" Gauss linearUpwind limited;
div(alphaRhoPhi,K) Gauss linearUpwind limited;
div(alphaRhoPhi,(p|rho)) Gauss linearUpwind limited;
"div\(alphaRhoPhi,(k|epsilon|omega)\)" Gauss linearUpwind limited;
"div\(phim,(k|epsilon)m\)" Gauss linearUpwind limited;
div((((alpha*rho)*nuEff)*dev2(T(grad(U))))) Gauss linear;
} }
laplacianSchemes laplacianSchemes
{ {
default Gauss linear corrected;
} }
interpolationSchemes interpolationSchemes
{ {
default linear;
} }
snGradSchemes snGradSchemes
{ {
default corrected;
} }
fluxRequired fluxRequired
{ {
default no;
}
wallDist
{
method meshWave;
nRequired yes;
} }

70
tutorials/modules/CHT/wallBoiling/system/fvSolution
View File

@ -9,79 +9,13 @@ FoamFile
{ {
format ascii; format ascii;
class dictionary; class dictionary;
location "system";
object fvSolution; object fvSolution;
} }
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers "(PIMPLE|PISO)"
{ {
"alpha.*" nOuterCorrectors 1;
{
nAlphaCorr 1;
nAlphaSubCycles 1;
}
p_rgh
{
solver GAMG;
smoother DIC;
tolerance 1e-8;
relTol 0.01;
maxIter 20;
minIter 2;
}
p_rghFinal
{
$p_rgh;
relTol 0;
}
"(e|h).*"
{
solver PBiCGStab;
preconditioner DILU;
tolerance 1e-12;
relTol 0.001;
minIter 1;
maxIter 20;
}
"(k|epsilon|omega).*"
{
solver PBiCGStab;
preconditioner DILU;
tolerance 1e-8;
relTol 0;
minIter 1;
}
} }
PIMPLE
{
nOuterCorrectors 1;
nCorrectors 1;
nNonOrthogonalCorrectors 0;
nEnergyCorrectors 1;
faceMomentum yes;
}
relaxationFactors
{
fields
{
thermalPhaseChange:dmdtf 1.0;
}
equations
{
".*" 1;
"h\..*" 1;
}
}
// ************************************************************************* // // ************************************************************************* //

16
tutorials/modules/CHT/wallBoiling/validation/createGraphs
View File

@ -40,31 +40,31 @@ gnuplot<<EOF
set key at graph 0.99,0.95 set key at graph 0.99,0.95
set xlabel 'Time (s)' set xlabel 'Time (s)'
set xrange [0.5:10] set xrange [0.5:8]
set ylabel 'Mass in- and outflow (kg/s)' set ylabel 'Mass in- and outflow (kg/s)'
set yrange [0.5:0.8] set yrange [0.5:0.8]
plot \ plot \
'./postProcessing/fluid/inflow/$surfaceFile' \ './postProcessing/fluid/inflow/$surfaceFile' \
every ::4::20000 u 1:(-1*360*\$3) w l lt 1 t 'liquid inflow', \ every ::10::20000 u 1:(-1*360*\$3) w l lt 1 t 'liquid inflow', \
'./postProcessing/fluid/outflow/$surfaceFile' \ './postProcessing/fluid/outflow/$surfaceFile' \
every ::4::20000 u 1:(360*\$3) w l lt 3 t 'liquid outflow', \ every ::10::20000 u 1:(360*\$3) w l lt 3 t 'liquid outflow', \
'./postProcessing/fluid/outflow/$surfaceFile' \ './postProcessing/fluid/outflow/$surfaceFile' \
every ::4::20000 u 1:(360*\$2+360*\$3) w l lt 4 t 'liquid + vapor outflow' every ::10::20000 u 1:(360*\$2+360*\$3) w l lt 4 t 'liquid + vapor outflow'
set ylabel 'Enthalpy in- and outflow (J/s)' set ylabel 'Enthalpy in- and outflow (J/s)'
set yrange [1.5e5:2.0e5] set yrange [1.5e5:2.0e5]
plot \ plot \
'./postProcessing/fluid/inletLiquid/$surfaceFile' \ './postProcessing/fluid/inletLiquid/$surfaceFile' \
every ::4::20000 u 1:(-1*360*\$2) w l lt 1 t 'liquid inflow', \ every ::10::20000 u 1:(-1*360*\$2) w l lt 1 t 'liquid inflow', \
'./postProcessing/fluid/inletLiquid/$surfaceFile' \ './postProcessing/fluid/inletLiquid/$surfaceFile' \
every ::4::20000 u 1:(-1*360*\$2+360*A*q) w l lt 8 t 'liquid inflow + wall heat flux', \ every ::10::20000 u 1:(-1*360*\$2+360*A*q) w l lt 8 t 'liquid inflow + wall heat flux', \
'./postProcessing/fluid/outletLiquid/$surfaceFile' \ './postProcessing/fluid/outletLiquid/$surfaceFile' \
every ::4::20000 u 1:(360*\$2) w l lt 3 t 'liquid outflow', \ every ::10::20000 u 1:(360*\$2) w l lt 3 t 'liquid outflow', \
"<paste \ "<paste \
./postProcessing/fluid/outletLiquid/$surfaceFile \ ./postProcessing/fluid/outletLiquid/$surfaceFile \
./postProcessing/fluid/outletGas/$surfaceFile" \ ./postProcessing/fluid/outletGas/$surfaceFile" \
every ::4::20000 u 1:(360*\$2+360*\$4) w l lt 4 t 'liquid + vapor outflow' every ::10::20000 u 1:(360*\$2+360*\$4) w l lt 4 t 'liquid + vapor outflow'
unset multiplot unset multiplot