Commit 9cd66598 authored by Henry's avatar Henry
Browse files

tutorials/multiphase/interPhaseChangeFoam/cavitatingBullet: Added a potentialFoam initializer step

parent b44acdb0
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: dev |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object Phi;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -1 0 0];
internalField uniform 0;
boundaryField
{
inlet
{
type zeroGradient;
}
outlet
{
type fixedValue;
value $internalField;
}
walls
{
type symmetry;
}
bullet
{
type zeroGradient;
}
}
// ************************************************************************* //
......@@ -4,9 +4,12 @@ cd ${0%/*} || exit 1 # Run from this directory
# Source tutorial clean functions
. $WM_PROJECT_DIR/bin/tools/CleanFunctions
# remove surface
# Remove surface
\rm -f constant/triSurface/bullet.stl.gz
# Remove copy of 0.org
rm -rf 0
cleanCase
# ----------------------------------------------------------------- end-of-file
......@@ -4,7 +4,7 @@ cd ${0%/*} || exit 1 # Run from this directory
# Source tutorial run functions
. $WM_PROJECT_DIR/bin/tools/RunFunctions
# copy bullet surface from resources directory
# Copy bullet surface from resources directory
cp $FOAM_TUTORIALS/resources/geometry/bullet.stl.gz constant/triSurface/
# Generate the base block mesh
......@@ -13,6 +13,11 @@ runApplication blockMesh
# Generate the snappy mesh
runApplication snappyHexMesh -overwrite
cp -r 0.org 0
# Initialise with potentialFoam solution
runApplication potentialFoam -pName p_rgh
# Run the solver
runApplication `getApplication`
......
......@@ -39,6 +39,7 @@ FoamFile
bullet
{
type wall;
inGroups 1(wall);
nFaces 37743;
startFace 1133431;
}
......
......@@ -16,43 +16,44 @@ FoamFile
ddtSchemes
{
default Euler;
default Euler;
}
interpolationSchemes
{
default linear;
default linear;
}
divSchemes
{
default none;
div(rhoPhi,U) Gauss linearUpwind grad(U);
div(phi,omega) Gauss linearUpwind grad(omega);
div(phi,k) Gauss linearUpwind grad(k);
div(phi,alpha) Gauss vanLeer;
div(phirb,alpha) Gauss linear;
default none;
div(rhoPhi,U) Gauss linearUpwind grad(U);
div(phi,omega) Gauss linearUpwind grad(omega);
div(phi,k) Gauss linearUpwind grad(k);
div(phi,alpha) Gauss vanLeer;
div(phirb,alpha) Gauss linear;
div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}
gradSchemes
{
default Gauss linear;
default Gauss linear;
}
laplacianSchemes
{
default Gauss linear limited corrected 0.5;
default Gauss linear limited corrected 0.5;
}
snGradSchemes
{
default limited corrected 0.5;
default limited corrected 0.5;
}
fluxRequired
{
default none;
default none;
Phi;
p_rgh;
pcorr;
alpha.water;
......
......@@ -59,12 +59,6 @@ solvers
maxIter 50;
};
pcorr
{
$p_rgh;
relTol 0;
};
p_rghFinal
{
solver PCG;
......@@ -91,6 +85,23 @@ solvers
relTol 0;
maxIter 50;
};
pcorr
{
$p_rgh;
relTol 0;
};
Phi
{
$p_rgh;
relTol 0;
};
}
potentialFlow
{
nNonOrthogonalCorrectors 3;
}
PIMPLE
......@@ -103,12 +114,9 @@ PIMPLE
relaxationFactors
{
fields
{
}
equations
{
"U.*" 1;
"U.*" 1;
}
}
......
Markdown is supported
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