| MtnRunBeachBum |
April 30, 2015 12:39 |
Transonic, Compressible, full flight scale Boundary Conditions?
OK, I know I've been reaching out a lot, but everyone has been very helpful, and I'm learning a lot quickly. I think I'm finally stuck on settign up boudnary conditions.
The case is rhoCentralFoam, ~0.82 Mach, in the neighborhood of 30000'. See below for blockMeshDict, U, p, and T files:
blockMeshDict
Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
convertToMeters 39.3700787; //0.0254; //(meters to inches)
vertices
/*(
(-2.54 -15.24 -2.54) // (-100 -600 -100) // full grid
(30.48 -15.24 -2.54) // (1200 -600 -100)
(30.48 15.24 -2.54) // (1200 600 -100)
(-2.54 15.24 -2.54) // (-100 600 -100)
(-2.54 -15.24 10.16) // (-100 -600 400)
(30.48 -15.24 10.16) // (1200 -600 400)
(30.48 15.24 10.16) // (1200 600 400)
(-2.54 15.24 10.16) // (-100 600 400)
);*/
(
(-2.54 0 -2.54) // (-100 -600 -100) // half grid
(30.48 0 -2.54) // (1200 -600 -100)
(30.48 15.24 -2.54) // (1200 600 -100)
(-2.54 15.24 -2.54) // (-100 600 -100)
(-2.54 0 10.16) // (-100 -600 400)
(30.48 0 10.16) // (1200 -600 400)
(30.48 15.24 10.16) // (1200 600 400)
(-2.54 15.24 10.16) // (-100 600 400)
);
blocks
(
// hex (0 1 2 3 4 5 6 7) (100 100 100) simpleGrading (1 1 1)
hex (0 3 7 4 1 2 6 5) (50 25 25) simpleGrading (1 1 1) // sets the block in order from X_min -> X_max
);
edges
(
);
boundary
(
inlet
{
type patch;
faces
(
(0 4 7 3)
);
}
outlet
{
type patch;
faces
(
(1 2 6 5)
);
}
bottom
{
type symmetryPlane;
faces
(
(0 1 5 4)
);
}
top
{
type patch;
faces
(
(2 3 7 6)
);
}
front
{
type patch;
faces
(
(0 3 2 1)
);
}
back
{
type patch;
faces
(
(4 5 6 7)
);
}
);
mergePatchPairs
(
);
// ************************************************************************* //
p-file
Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -1 -2 0 0 0 0]; //compressible
//dimensions [0 2 -2 0 0 0 0]; //incompressible
internalField uniform 32932;
boundaryField
{
inlet
{
type fixedValue;
value uniform 32932;
}
outlet
{
type zeroGradient;
value uniform 32932;
}
bottom
{
type symmetryPlane;
//type zeroGradient; //fixedValue;
//value uniform 32932;
}
top
{
//type symmetryPlane;
type zeroGradient; //fixedValue;
value uniform 32932;
}
front
{
//type symmetryPlane;
type zeroGradient; //fixedValue;
value uniform 32932;
}
back
{
//type symmetryPlane;
type zeroGradient; //fixedValue;
value uniform 32932;
}
fuselage
{
type zeroGradient;
value uniform 32932;
}
wings
{
type zeroGradient;
value uniform 32932;
}
vertTail
{
type zeroGradient;
value uniform 32932;
}
horizStab
{
type zeroGradient;
value uniform 32932;
}
engines
{
type zeroGradient;
value uniform 32932;
}
pylon
{
type zeroGradient;
value uniform 32932;
}
pylonExtension
{
type zeroGradient;
value uniform 32932;
}
pod
{
type zeroGradient;
value uniform 32932;
}
}
// ************************************************************************* //
U-file
Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (251 0 0);
boundaryField
{
inlet
{
type fixedValue;
value uniform (251 0 0);
}
outlet
{
type zeroGradient; //freestream; //inletOutlet;
//inletValue uniform (251 0 0); //
//freestreamValue uniform (251 0 0);
value uniform (251 0 0);
}
bottom
{
//type zeroGradient;
type symmetryPlane; //fixedValue;
//value uniform ( 0 0 0 ); //(0 0 0);
}
top
{
type zeroGradient;
//type symmetryPlane; //fixedValue;
value uniform ( 0 0 0 );
}
front
{
type zeroGradient;
//type symmetryPlane; //fixedValue;
value uniform ( 0 0 0 ); //(251 0 0)
}
back
{
type zeroGradient;
//type symmetryPlane; //fixedValue;
value uniform ( 0 0 0 );
}
fuselage
{
type slip;
}
wings
{
type slip;
}
vertTail
{
type slip;
}
horizStab
{
type slip;
}
engines
{
type slip;
}
pylon
{
type slip;
}
pylonExtension
{
type slip;
}
pod
{
type slip;
}
}
// ************************************************************************* //
T-file
Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 1 0 0 0];
internalField uniform 233;
boundaryField
{
inlet
{
type fixedValue;
value uniform 233;
}
outlet
{
type zeroGradient; //inletOutlet;
value uniform 233;
//inletValue uniform 233;
}
bottom
{
type symmetryPlane; //fixedValue;
//value uniform 233;
}
top
{
type fixedValue; //symmetryPlane; //fixedValue;
value uniform 233;
}
front
{
type fixedValue; //symmetryPlane; //fixedValue;
value uniform 233;
}
back
{
type fixedValue; //symmetryPlane; //fixedValue;
value uniform 233;
}
fuselage
{
type zeroGradient;
value uniform 233;
}
wings
{
type zeroGradient;
value uniform 233;
}
vertTail
{
type zeroGradient;
value uniform 233;
}
horizStab
{
type zeroGradient;
value uniform 233;
}
engines
{
type zeroGradient;
value uniform 233;
}
pylon
{
type zeroGradient;
value uniform 233;
}
pylonExtension
{
type zeroGradient;
value uniform 233;
}
pod
{
type zeroGradient;
value uniform 233;
}
}
// ************************************************************************* //
Based on the errors I get (very similar to this threat), I think most of the problem is in my boundary condition setup; however I am at wit's end...I am having trouble accounting for what looks to me to be differences in how the "patches" and so on are used for different type of solvers.
Any help or guidance would be appreciated, I've been pulling my hair out since last week, and really need it to run roday
Thanks! |
