Commit 441b8b08 authored by Andrew Heather's avatar Andrew Heather
Browse files

more tutorial updates

parent 81dfb596
......@@ -22,25 +22,25 @@ boundaryField
{
floor
{
type calculated;
type fixedFluxBuoyantPressure;
value uniform 100000;
}
ceiling
{
type calculated;
type fixedFluxBuoyantPressure;
value uniform 100000;
}
fixedWalls
{
type calculated;
type fixedFluxBuoyantPressure;
value uniform 100000;
}
box
{
type calculated;
type fixedFluxBuoyantPressure;
value uniform 100000;
}
}
......
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.5 |
| \\ / A nd | Web: http://www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object pd;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -1 -2 0 0 0 0];
internalField uniform 0;
boundaryField
{
floor
{
type fixedFluxBuoyantPressure;
value uniform 0;
}
ceiling
{
type fixedFluxBuoyantPressure;
value uniform 0;
}
fixedWalls
{
type fixedFluxBuoyantPressure;
value uniform 0;
}
box
{
type fixedFluxBuoyantPressure;
value uniform 0;
}
}
// ************************************************************************* //
......@@ -41,7 +41,7 @@ laplacianSchemes
{
default none;
laplacian(muEff,U) Gauss linear corrected;
laplacian((rho*(1|A(U))),pd) Gauss linear corrected;
laplacian((rho*(1|A(U))),p) Gauss linear corrected;
laplacian(alphaEff,h) Gauss linear corrected;
laplacian(DkEff,k) Gauss linear corrected;
laplacian(DepsilonEff,epsilon) Gauss linear corrected;
......@@ -62,7 +62,7 @@ snGradSchemes
fluxRequired
{
default no;
pd ;
p;
}
......
......@@ -17,7 +17,7 @@ FoamFile
solvers
{
pd
p
{
solver PCG;
preconditioner DIC;
......@@ -69,14 +69,14 @@ solvers
SIMPLE
{
nNonOrthogonalCorrectors 0;
pdRefCell 0;
pdRefValue 0;
pRefCell 0;
pRefValue 0;
}
relaxationFactors
{
rho 1;
pd 0.3;
p 0.3;
U 0.7;
h 0.7;
k 0.7;
......
......@@ -22,25 +22,25 @@ boundaryField
{
floor
{
type calculated;
type fixedFluxBuoyantPressure;
value uniform 100000;
}
ceiling
{
type calculated;
type fixedFluxBuoyantPressure;
value uniform 100000;
}
fixedWalls
{
type calculated;
type fixedFluxBuoyantPressure;
value uniform 100000;
}
box
{
type calculated;
type fixedFluxBuoyantPressure;
value uniform 100000;
}
}
......
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.5 |
| \\ / A nd | Web: http://www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object pd;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -1 -2 0 0 0 0];
internalField uniform 0;
boundaryField
{
floor
{
type fixedFluxBuoyantPressure;
value uniform 0;
}
ceiling
{
type fixedFluxBuoyantPressure;
value uniform 0;
}
fixedWalls
{
type fixedFluxBuoyantPressure;
value uniform 0;
}
box
{
type fixedFluxBuoyantPressure;
value uniform 0;
}
}
// ************************************************************************* //
......@@ -42,7 +42,7 @@ laplacianSchemes
{
default none;
laplacian(muEff,U) Gauss linear corrected;
laplacian((rho*(1|A(U))),pd) Gauss linear corrected;
laplacian((rho*(1|A(U))),p) Gauss linear corrected;
laplacian(alphaEff,h) Gauss linear corrected;
laplacian(DkEff,k) Gauss linear corrected;
laplacian(DepsilonEff,epsilon) Gauss linear corrected;
......@@ -63,7 +63,7 @@ snGradSchemes
fluxRequired
{
default no;
pd ;
p;
}
......
......@@ -17,7 +17,7 @@ FoamFile
solvers
{
pd
p
{
solver GAMG;
tolerance 1e-06;
......@@ -73,14 +73,14 @@ solvers
SIMPLE
{
nNonOrthogonalCorrectors 0;
pdRefCell 0;
pdRefValue 0;
pRefCell 0;
pRefValue 0;
}
relaxationFactors
{
rho 1.0;
pd 0.3;
p 0.3;
U 0.7;
h 0.7;
k 0.7;
......
......@@ -15,13 +15,13 @@ runApplication splitMeshRegions -cellZones
for i in bottomAir topAir heater leftSolid rightSolid
do
changeDictionary -region $i
changeDictionary -region $i >& log.changeDictionary.$i
done
# remove fluid fields from solid regions (important for post-processing)
for i in heater leftSolid rightSolid
do
rm -f 0*/$i/{mut,alphat,epsilon,k,p,pd,U}
rm -f 0*/$i/{mut,alphat,epsilon,k,p,p,U}
done
# remove solid fields from fluid regions (important for post-processing)
......@@ -37,7 +37,7 @@ done
# Decompose
for i in bottomAir topAir heater leftSolid rightSolid
do
decomposePar -region $i
decomposePar -region $i >& log.decomposePar.$i
done
# Run
......@@ -47,7 +47,7 @@ runParallel chtMultiRegionFoam 4 system/machines
# Reconstruct
for i in bottomAir topAir heater leftSolid rightSolid
do
reconstructPar -region $i
reconstructPar -region $i >& log.reconstructPar.$i
done
......
#!/bin/sh
set -x
m4 < constant/polyMesh/blockMeshDict.m4 > constant/polyMesh/blockMeshDict
blockMesh
cellSet
blockMesh >& log.blockMesh
cellSet >& log.cellSet
#- MRF determines its own faceZone if not supplied
#cp system/faceSetDict_rotorFaces system/faceSetDict
#faceSet
#cp system/faceSetDict_noBoundaryFaces system/faceSetDict
#faceSet
setsToZones -noFlipMap
setsToZones -noFlipMap >& log.setsToZones
......@@ -17,7 +17,7 @@ FoamFile
radiation off;
radiationModel noRadiation;
radiationModel none;
solverFreq 10;
......
......@@ -8,7 +8,7 @@
FoamFile
{
version 2.0;
format binary;
format ascii;
class dictionary;
location "system";
object cellSetDict.1;
......@@ -19,7 +19,13 @@ name c0;
action new;
topoSetSources ( boxToCell { box ( 0.004 -0.001 -1 ) ( 0.012 0.001 1 ) ; } );
topoSetSources
(
boxToCell
{
box (0.004 -0.001 -1) (0.012 0.001 1);
}
);
// ************************************************************************* //
......@@ -8,7 +8,7 @@
FoamFile
{
version 2.0;
format binary;
format ascii;
class dictionary;
location "system";
object cellSetDict.2;
......@@ -19,7 +19,13 @@ name c0;
action new;
topoSetSources ( boxToCell { box ( 0.0045 -0.00075 -1 ) ( 0.0095 0.00075 1 ) ; } );
topoSetSources
(
boxToCell
{
box (0.0045 -0.00075 -1) (0.0095 0.00075 1);
}
);
// ************************************************************************* //
......@@ -8,7 +8,7 @@
FoamFile
{
version 2.0;
format binary;
format ascii;
class dictionary;
location "system";
object cellSetDict.3;
......@@ -19,7 +19,13 @@ name c0;
action new;
topoSetSources ( boxToCell { box ( 0.00475 -0.000375 -1 ) ( 0.009 0.000375 1 ) ; } );
topoSetSources
(
boxToCell
{
box (0.00475 -0.000375 -1) (0.009 0.000375 1);
}
);
// ************************************************************************* //
......@@ -8,7 +8,7 @@
FoamFile
{
version 2.0;
format binary;
format ascii;
class dictionary;
location "system";
object cellSetDict.1;
......@@ -19,7 +19,13 @@ name c0;
action new;
topoSetSources ( boxToCell { box ( 0.004 -0.001 -1 ) ( 0.012 0.001 1 ) ; } );
topoSetSources
(
boxToCell
{
box ( 0.004 -0.001 -1 ) ( 0.012 0.001 1 );
}
);
// ************************************************************************* //
......@@ -8,7 +8,7 @@
FoamFile
{
version 2.0;
format binary;
format ascii;
class dictionary;
location "system";
object cellSetDict.2;
......@@ -19,7 +19,13 @@ name c0;
action new;
topoSetSources ( boxToCell { box ( 0.0045 -0.00075 -1 ) ( 0.0095 0.00075 1 ) ; } );
topoSetSources
(
boxToCell
{
box ( 0.0045 -0.00075 -1 ) ( 0.0095 0.00075 1 );
}
);
// ************************************************************************* //
......@@ -8,7 +8,7 @@
FoamFile
{
version 2.0;
format binary;
format ascii;
class dictionary;
location "system";
object cellSetDict.3;
......@@ -19,7 +19,13 @@ name c0;
action new;
topoSetSources ( boxToCell { box ( 0.00475 -0.000375 -1 ) ( 0.009 0.000375 1 ) ; } );
topoSetSources
(
boxToCell
{
box ( 0.00475 -0.000375 -1 ) ( 0.009 0.000375 1 );
}
);
// ************************************************************************* //
......@@ -10,7 +10,7 @@ FoamFile
version 2.0;
format ascii;
class volScalarField;
object pd;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
......
......@@ -7,6 +7,6 @@ application="compressibleInterFoam"
runApplication blockMesh
cp 0/alpha1.org 0/alpha1
cp 0/pd.org 0/pd
cp 0/p.org 0/p
runApplication setFields
runApplication $application
......@@ -30,7 +30,7 @@ divSchemes
div(rho*phi,U) Gauss upwind;
div(phi,alpha) Gauss vanLeer;
div(phirb,alpha) Gauss interfaceCompression 1;
div(phi,pd) Gauss upwind;
div(phi,p) Gauss upwind;
div(phi,k) Gauss vanLeer;
div((nuEff*dev(grad(U).T()))) Gauss linear;
}
......@@ -53,9 +53,9 @@ snGradSchemes
fluxRequired
{
default no;
pd ;
pcorr ;
gamma ;
p;
pcorr;
gamma;
}
......
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