CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

rhocentralFoam error while modelling hypersonic flow

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   June 19, 2012, 11:33
Default rhocentralFoam error while modelling hypersonic flow
  #1
New Member
 
Nathan Donaldson
Join Date: Apr 2012
Location: Scotland
Posts: 6
Rep Power: 5
Nathan_Donaldson is on a distinguished road
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)
        );
    }
);

// ************************************************************************* //
Nathan_Donaldson is offline   Reply With Quote

Reply

Tags
hypersonic, openfoam 2.0.1, residuals, rhocentralfoam

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Solver for hypersonic flow in OpenFoam Santos-Dumont OpenFOAM Running, Solving & CFD 19 October 21, 2012 01:42
modelling slug flow with pressure inlet Azman FLUENT 0 May 29, 2006 06:45
modelling laminar and turbulent flow in the same Ol FLUENT 2 November 25, 2005 03:52
Modelling flow around a ship's hull using vof Manoj Kumar FLUENT 0 February 26, 2005 17:22
modelling flow in microchannel Wang FLUENT 6 October 10, 2000 09:28


All times are GMT -4. The time now is 08:13.