Improvements in conjugate heat solver: non-matching interface, parallelisation

This commit is contained in:
Hrvoje Jasak 2011-09-29 02:13:00 +01:00
parent 617fcb0d32
commit 30fdce1de4
9 changed files with 152 additions and 45 deletions

View file

@ -6,10 +6,14 @@ application="conjugateHeatFoam"
runApplication blockMesh runApplication blockMesh
cp boundary constant/polyMesh cp boundary constant/polyMesh
runApplication setSet -batch setBatch
runApplication setsToZones -noFlipMap
cd ../heatedBlock cd ../heatedBlock
runApplication blockMesh runApplication blockMesh
cp boundary constant/polyMesh cp boundary constant/polyMesh
runApplication setSet -batch setBatch
runApplication setsToZones -noFlipMap
cd ../conjugateCavity cd ../conjugateCavity
runApplication $application runApplication $application

View file

@ -1,6 +1,6 @@
/*--------------------------------*- C++ -*----------------------------------*\ /*--------------------------------*- C++ -*----------------------------------*\
| ========= | | | ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open Source CFD | | \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext | | \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: www.extend-project.de | | \\ / A nd | Web: www.extend-project.de |
| \\/ M anipulation | | | \\/ M anipulation | |
@ -10,51 +10,52 @@ FoamFile
version 2.0; version 2.0;
format ascii; format ascii;
class polyBoundaryMesh; class polyBoundaryMesh;
location "constant/polyMesh";
object boundary; object boundary;
} }
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
5 5
( (
movingWall movingWall
{ {
type wall; type wall;
nFaces 10; nFaces 15;
startFace 180; startFace 420;
} }
left
{
type wall;
nFaces 15;
startFace 435;
}
right
{
type regionCouple;
nFaces 15;
startFace 450;
left zone rightZone;
{ shadowRegion solid;
type wall; shadowPatch left;
nFaces 10;
startFace 190;
}
right master yes;
{ attached on;
type regionCouple; isWall yes;
nFaces 10; bridgeOverlap no;
startFace 200; }
bottom
shadowRegion solid; {
shadowPatch left; type wall;
attached on; nFaces 15;
isWall yes; startFace 465;
} }
frontAndBack
bottom {
{ type empty;
type wall; nFaces 450;
nFaces 10; startFace 480;
startFace 210; }
}
frontAndBack
{
type empty;
nFaces 200;
startFace 220;
}
) )
// ************************************************************************* // // ************************************************************************* //

View file

@ -30,7 +30,7 @@ vertices
blocks blocks
( (
hex (0 1 2 3 4 5 6 7) (10 10 1) simpleGrading (1 1 1) hex (0 1 2 3 4 5 6 7) (15 15 1) simpleGrading (1 1 1)
); );
edges edges

View file

@ -0,0 +1,2 @@
faceSet rightZone new patchToFace right
quit

View file

@ -29,10 +29,14 @@ left
nFaces 10; nFaces 10;
startFace 200; startFace 200;
zone leftZone;
shadowRegion region0; shadowRegion region0;
shadowPatch right; shadowPatch right;
master no;
attached on; attached on;
isWall yes; // Not relevant isWall yes; // Not relevant
bridgeOverlap no;
} }
right right

View file

@ -0,0 +1,2 @@
faceSet leftZone new patchToFace left
quit

View file

@ -0,0 +1,42 @@
ln -s ../conjugateCavity/0.1/solid/ 0.1
ln -s ../conjugateCavity/0.2/solid/ 0.2
ln -s ../conjugateCavity/0.3/solid/ 0.3
ln -s ../conjugateCavity/0.4/solid/ 0.4
ln -s ../conjugateCavity/0.5/solid/ 0.5
ln -s ../conjugateCavity/0.6/solid/ 0.6
ln -s ../conjugateCavity/0.7/solid/ 0.7
ln -s ../conjugateCavity/0.8/solid/ 0.8
ln -s ../conjugateCavity/0.9/solid/ 0.9
ln -s ../conjugateCavity/0.15/solid/ 0.15
ln -s ../conjugateCavity/0.25/solid/ 0.25
ln -s ../conjugateCavity/0.35/solid/ 0.35
ln -s ../conjugateCavity/0.45/solid/ 0.45
ln -s ../conjugateCavity/0.55/solid/ 0.55
ln -s ../conjugateCavity/0.65/solid/ 0.65
ln -s ../conjugateCavity/0.75/solid/ 0.75
ln -s ../conjugateCavity/0.85/solid/ 0.85
ln -s ../conjugateCavity/0.95/solid/ 0.95
ln -s ../conjugateCavity/1.1/solid/ 1.1
ln -s ../conjugateCavity/1.2/solid/ 1.2
ln -s ../conjugateCavity/1.3/solid/ 1.3
ln -s ../conjugateCavity/1.4/solid/ 1.4
ln -s ../conjugateCavity/1.5/solid/ 1.5
ln -s ../conjugateCavity/1.6/solid/ 1.6
ln -s ../conjugateCavity/1.7/solid/ 1.7
ln -s ../conjugateCavity/1.8/solid/ 1.8
ln -s ../conjugateCavity/1.9/solid/ 1.9
ln -s ../conjugateCavity/1.15/solid/ 1.15
ln -s ../conjugateCavity/1.25/solid/ 1.25
ln -s ../conjugateCavity/1.35/solid/ 1.35
ln -s ../conjugateCavity/1.45/solid/ 1.45
ln -s ../conjugateCavity/1.55/solid/ 1.55
ln -s ../conjugateCavity/1.65/solid/ 1.65
ln -s ../conjugateCavity/1.75/solid/ 1.75
ln -s ../conjugateCavity/1.85/solid/ 1.85
ln -s ../conjugateCavity/1.95/solid/ 1.95

View file

@ -0,0 +1,52 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: www.extend-project.de |
| \\/ 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
{
lowerWall
{
type nutWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0;
}
upperWall
{
type nutWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0;
}
flow
{
type cyclic;
value uniform 0;
}
defaultFaces
{
type empty;
}
}
// ************************************************************************* //

View file

@ -22,42 +22,42 @@ solvers
preconditioner DIC; preconditioner DIC;
tolerance 1e-06; tolerance 1e-06;
relTol 0.01; relTol 0.01;
}; }
U U
{ {
solver PBiCG; solver PBiCG;
preconditioner DILU; preconditioner DILU;
tolerance 1e-05; tolerance 1e-05;
relTol 0.1; relTol 0.1;
}; }
k k
{ {
solver PBiCG; solver PBiCG;
preconditioner DILU; preconditioner DILU;
tolerance 1e-05; tolerance 1e-05;
relTol 0.1; relTol 0.1;
}; }
epsilon epsilon
{ {
solver PBiCG; solver PBiCG;
preconditioner DILU; preconditioner DILU;
tolerance 1e-05; tolerance 1e-05;
relTol 0.1; relTol 0.1;
}; }
R R
{ {
solver PBiCG; solver PBiCG;
preconditioner DILU; preconditioner DILU;
tolerance 1e-05; tolerance 1e-05;
relTol 0.1; relTol 0.1;
}; }
nuTilda nuTilda
{ {
solver PBiCG; solver PBiCG;
preconditioner DILU; preconditioner DILU;
tolerance 1e-05; tolerance 1e-05;
relTol 0.1; relTol 0.1;
}; }
} }
SIMPLE SIMPLE