CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   rhocentralFoam error while modelling hypersonic flow (http://www.cfd-online.com/Forums/openfoam/103500-rhocentralfoam-error-while-modelling-hypersonic-flow.html)

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)
        );
    }
);

// ************************************************************************* //



All times are GMT -4. The time now is 04:30.