Nathan_Donaldson |
June 19, 2012 11:33 |
rhocentralFoam error while modelling hypersonic flow
Hello all,
I'm a newcomer to OpenFOAM and am trying to model a hypersonic calibration model using rhoCentralFoam. After completing the Wedge15Ma5 and motorBike tutorials, I successfully created a mesh of my geometry using snappyHexMesh.
The problem I'm having is that whenever I run the case and post-process it, the results (velocity, temperature and pressure) are all completely uniform. checkMesh returns no errors, and I've used a very small time step (0.001). Also, the mean and maximum Courant numbers remain constant between iterations, and the residuals prefixed with 'rho' are perpetually 0.
The freestream boundary conditions I've tried to code are:
Mach number = 17.8
Pressure = 29.5 Pa
Temperature = 287K
I've pasted in the Ma, p, U, T and blockMeshDict files below. Are there any errors that are immediately obvious? I'm sure it'll be something basic that I'm missing.
Many thanks,
Nathan
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object Ma;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 0 0 0 0];
internalField uniform 17.8;
boundaryField
{
frontAndBack
{
type symmetryPlane;
}
outlet
{
type calculated;
value uniform 17.8;
}
inlet
{
type calculated;
value uniform 17.8;
}
lowerAndUpperWall
{
type symmetryPlane;
}
"HB2.*"
{
type calculated;
}
defaultFaces
{
type empty;
}
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -1 -2 0 0 0 0];
internalField uniform 29.5;
boundaryField
{
frontAndBack
{
type symmetryPlane;
}
inlet
{
type zeroGradient;
value uniform 29.5;
}
outlet
{
type zeroGradient;
value uniform 29.5;
}
lowerAndUpperWall
{
type symmetryPlane;
}
"HB2.*"
{
type zeroGradient;
}
defaultFaces
{
type empty;
}
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 1 0 0 0];
internalField uniform 287;
boundaryField
{
frontAndBack
{
type symmetryPlane;
}
inlet
{
type fixedValue;
value uniform 287;
}
outlet
{
type inletOutlet;
inletValue uniform 287;
value uniform 287;
}
lowerAndUpperWall
{
type symmetryPlane;
}
"HB2.*"
{
type zeroGradient;
}
defaultFaces
{
type empty;
}
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (17.8 0 0);
boundaryField
{
frontAndBack
{
type symmetryPlane;
}
inlet
{
type fixedValue;
value uniform (17.8 0 0);
}
outlet
{
type inletOutlet;
inletValue uniform (17.8 0 0);
value uniform (17.8 0 0);
}
lowerAndUpperWall
{
type symmetryPlane;
}
"HB2.*"
{
type zeroGradient;
}
defaultFaces
{
type empty;
}
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
convertToMeters 1;
vertices
(
(-100 -50 -50)
(-100 50 -50)
(-100 50 50)
(-100 -50 50)
(400 -50 -50)
(400 50 -50)
(400 50 50)
(400 -50 50)
);
blocks
(
hex (0 1 2 3 4 5 6 7) (10 10 50) simpleGrading (1 1 1)
);
edges
(
);
boundary
(
frontAndBack
{
type symmetryPlane;
faces
(
(2 3 7 6)
(0 1 5 4)
);
}
inlet
{
type patch;
faces
(
(0 3 2 1)
);
}
outlet
{
type patch;
faces
(
(4 5 6 7)
);
}
lowerAndUpperWall
{
type symmetryPlane;
faces
(
(0 4 7 3)
(1 2 6 5)
);
}
);
// ************************************************************************* //
|