Skip to content
Snippets Groups Projects
Commit be9f2b33 authored by Andrew Heather's avatar Andrew Heather
Browse files
parents 26e04aab 5a83c466
Branches
Tags
No related merge requests found
Showing
with 89 additions and 50 deletions
......@@ -46,3 +46,5 @@
<< " Max(alpha) = " << max(alpha).value()
<< endl;
}
rho = alpha*rhoa + beta*rhob;
......@@ -156,6 +156,16 @@
+ fvc::interpolate(beta)*phib
);
volScalarField rho
(
IOobject
(
"rho",
runTime.timeName(),
mesh
),
alpha*rhoa + beta*rhob
);
#include "createRASTurbulence.H"
......
......@@ -52,3 +52,5 @@
<< endl;
}
}
rho = alpha*rhoa + beta*rhob;
......@@ -145,6 +145,16 @@
fvc::interpolate(alpha)*phia + fvc::interpolate(beta)*phib
);
volScalarField rho
(
IOobject
(
"rho",
runTime.timeName(),
mesh
),
alpha*rhoa + beta*rhob
);
IOdictionary RASProperties
(
......
......@@ -37,12 +37,13 @@ getApplication ()
runApplication ()
{
APP_RUN=$1; shift
APP_NAME=${APP_RUN##*/}
if [ -f log.$APP_RUN ] ; then
echo "$APP_RUN already run on $PWD: remove log file to run"
if [ -f log.$APP_NAME ] ; then
echo "$APP_NAME already run on $PWD: remove log file to run"
else
echo "Running $APP_RUN on $PWD"
$APP_RUN $* > log.$APP_RUN 2>&1
$APP_RUN $* > log.$APP_NAME 2>&1
fi
}
......
......@@ -16,7 +16,6 @@ FoamFile
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -2 0 0 0 0];
value ( 0 -9.81 0 );
value (0 0 0);
// ************************************************************************* //
......@@ -16,7 +16,6 @@ FoamFile
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -2 0 0 0 0];
value ( 0 0 0 );
value (0 0 0);
// ************************************************************************* //
......@@ -22,19 +22,19 @@ boundaryField
{
walls
{
type zeroGradient;
type buoyantPressure;
value uniform 5e+06;
}
front
{
type wedge;
type wedge;
}
back
{
type wedge;
type wedge;
}
}
// ************************************************************************* //
......@@ -23,17 +23,17 @@ boundaryField
{
piston
{
type epsilonWallFunction;
type compressible::epsilonWallFunction;
value uniform 450;
}
liner
{
type epsilonWallFunction;
type compressible::epsilonWallFunction;
value uniform 450;
}
cylinderHead
{
type epsilonWallFunction;
type compressible::epsilonWallFunction;
value uniform 450;
}
}
......
......@@ -23,17 +23,17 @@ boundaryField
{
piston
{
type kQRWallFunction;
type compressible::kQRWallFunction;
value uniform 4;
}
liner
{
type kQRWallFunction;
type compressible::kQRWallFunction;
value uniform 4;
}
cylinderHead
{
type kQRWallFunction;
type compressible::kQRWallFunction;
value uniform 4;
}
}
......
../bottomAir/g
\ No newline at end of file
......@@ -20,23 +20,24 @@ internalField uniform 0;
boundaryField
{
inlet
inlet
{
type zeroGradient;
}
outlet
outlet
{
type fixedValue;
value uniform 0;
}
walls
walls
{
type zeroGradient;
type buoyantPressure;
value uniform 0;
}
defaultFaces
defaultFaces
{
type empty;
}
......
......@@ -16,7 +16,7 @@ FoamFile
dimensions [1 -1 -2 0 0 0 0];
internalField nonuniform List<scalar>
internalField nonuniform List<scalar>
1875
(
6296.13
......@@ -1899,23 +1899,24 @@ internalField nonuniform List<scalar>
boundaryField
{
inlet
inlet
{
type zeroGradient;
}
outlet
outlet
{
type fixedValue;
value uniform 0;
}
walls
walls
{
type zeroGradient;
type buoyantPressure;
value uniform 0;
}
defaultFaces
defaultFaces
{
type empty;
}
......
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.5 |
| \\ / A nd | Web: http://www.OpenFOAM.org |
| \\ / O peration | Version: dev |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
......@@ -21,25 +21,25 @@ FoamFile
{
type patch;
nFaces 51;
startFace 15151;
startFace 57362;
}
outlet
{
type patch;
nFaces 51;
startFace 15202;
startFace 57413;
}
walls
{
type wall;
nFaces 436;
startFace 15253;
nFaces 836;
startFace 57464;
}
frontBack
{
type empty;
nFaces 15420;
startFace 15689;
nFaces 57540;
startFace 58300;
}
)
......
......@@ -5,5 +5,5 @@
application=`getApplication`
./makeMesh
runApplication ./makeMesh
runApplication $application
......@@ -22,17 +22,20 @@ boundaryField
{
leftWall
{
type zeroGradient;
type buoyantPressure;
value uniform 0;
}
rightWall
{
type zeroGradient;
type buoyantPressure;
value uniform 0;
}
lowerWall
{
type zeroGradient;
type buoyantPressure;
value uniform 0;
}
atmosphere
......
......@@ -22,17 +22,20 @@ boundaryField
{
leftWall
{
type zeroGradient;
type buoyantPressure;
value uniform 0;
}
rightWall
{
type zeroGradient;
type buoyantPressure;
value uniform 0;
}
lowerWall
{
type zeroGradient;
type buoyantPressure;
value uniform 0;
}
atmosphere
......
......@@ -22,17 +22,20 @@ boundaryField
{
leftWall
{
type zeroGradient;
type buoyantPressure;
value uniform 0;
}
rightWall
{
type zeroGradient;
type buoyantPressure;
value uniform 0;
}
lowerWall
{
type zeroGradient;
type buoyantPressure;
value uniform 0;
}
atmosphere
......
......@@ -22,17 +22,20 @@ boundaryField
{
leftWall
{
type zeroGradient;
type buoyantPressure;
value uniform 0;
}
rightWall
{
type zeroGradient;
type buoyantPressure;
value uniform 0;
}
lowerWall
{
type zeroGradient;
type buoyantPressure;
value uniform 0;
}
atmosphere
......
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