[ousoweak]

Hi,

I'm currently experimenting with the rhoSimpleFoam solver for steady-state calculations. The case is a simple 2D flat-plate flow with the following boundary conditions:
Inlet: freestream 120m/s, T = 288 K
Walls: T = 288 K
Outlet: freestream 120m/s
KOmegaSST.
Please see attached files for how I set up the boundary conditions with OpenFoam.
Calculations stop after 170+ timesteps and show error like this:
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::nutUSpaldingWallFunctionFvPatchScalarField:: calcUTau(Foam::Field<double> const&) const at ??:?
#4 Foam::nutUSpaldingWallFunctionFvPatchScalarField:: calcNut() const at ??:?
#5 Foam::nutWallFunctionFvPatchScalarField::updateCoe ffs() at ??:?
#6 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::Boundary::evaluate() at ??:?
#7 Foam::kOmegaSST<Foam::eddyViscosity<Foam::RASModel <Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Fo am::CompressibleTurbulenceModel<Foam::fluidThermo> > > > >, Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foa m::CompressibleTurbulenceModel<Foam::fluidThermo> > > >::correctNut(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#8 Foam::kOmegaSST<Foam::eddyViscosity<Foam::RASModel <Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Fo am::CompressibleTurbulenceModel<Foam::fluidThermo> > > > >, Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foa m::CompressibleTurbulenceModel<Foam::fluidThermo> > > >::correct() at ??:?
#9 ? at ??:?
#10 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11 ? at ??:?

Both the case and mesh are very simple.The case can run successfully in simpleFoam solver but failed with rhoSimpleFoam solver.I doubt whether I use the wrong boundary condition or initial condition.
Can anyone help me?Thanks.

Yiming

U code

dimensions      [0 1 -1 0 0 0 0];
internalField   uniform (120 0 0);
boundaryField
{
    inlet
    {
        type            freestream;
        freestreamValue uniform (120 0 0);
    }
    outlet
    {
        type            freestream;
        freestreamValue uniform (120 0 0);
    }
    walls
    {
        type            noSlip;
    }
    wall1
    {
        type            noSlip;
    }
    wall2
    {
        type            noSlip;
    }
    frontAndBack
    {
        type            empty;
    }
    symmdown
    {
        type            symmetryPlane;
    }
}

P code

dimensions      [1 -1 -2 0 0 0 0];
internalField   uniform 101325;
boundaryField
{
    inlet
    {
        type            zeroGradient;
        //type            freestreamPressure;
    }
    outlet
    {
        //type            zeroGradient;
        type            fixedValue;
        value           $internalField;
        //type            freestreamPressure;
    }
    walls
    {
        type            zeroGradient;
    }
    wall1
    {
        type            zeroGradient;
    }
    wall2
    {
        type            zeroGradient;
    }
    frontAndBack
    {
        type            empty;
    }
    symmdown
    {
        type            symmetryPlane;
    }
}

nut

dimensions      [0 2 -1 0 0 0 0];
internalField   uniform 1.48e-5;
boundaryField
{
    inlet
    {
        type            freestream;
        freestreamValue uniform 1.48e-5;
    }
    outlet
    {
        type            freestream;
        freestreamValue uniform 1.48e-5;
    }
    walls
    {
        type            nutUSpaldingWallFunction;
        value           uniform 0;
    }
    wall1
    {
        type            nutUSpaldingWallFunction;
        value           uniform 0;
    }
    wall2
    {
        type            nutUSpaldingWallFunction;
        value           uniform 0;
    }
    frontAndBack
    {
        type            empty;
    }
    symmdown
    {
        type            symmetryPlane;
    }
}

k

dimensions      [0 2 -2 0 0 0 0];
internalField   uniform  23.44;
boundaryField
{
    inlet
    {
        type            fixedValue;
        value           $internalField;
    }
       
    outlet
    {
        type            inletOutlet;
        inletValue      $internalField;
        value           $internalField;
    }
    walls
    {
        type            kqRWallFunction;
        value            $internalField;
    }
    wall1
    {
        type            kqRWallFunction;
        value            $internalField;
    }
    wall2
    {
        type            kqRWallFunction;
        value            $internalField;
    }
    frontAndBack
    {
        type   empty;
    }
    symmdown
    {
        type            symmetryPlane;
    }
}

omega

dimensions      [0 0 -1 0 0 0 0];
internalField   uniform 11.48;
boundaryField
{
    inlet
    {
        type            fixedValue;
        value           $internalField;
    }
    outlet
    {
        type            inletOutlet;
        inletValue      $internalField;
        value           $internalField;
    }
    walls
    {
        type            omegaWallFunction;
        value           $internalField;
    }
    wall1
    {
        type            omegaWallFunction;
        value           $internalField;
    }
    wall2
    {
        type            omegaWallFunction;
        value           $internalField;
    }
    frontAndBack
    {
        type            empty;
    }
    symmdown
    {
        type            symmetryPlane;
    }
}

attached is the complete case with out mesh file(Mesh file is too large)

[Wang Shang]

Dear Xu，

I have met the same problem, how do you sloved this problem?