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

compressible flow calculation error when using rhoSimpleFoam solver

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 10, 2017, 12:17
Default compressible flow calculation error when using rhoSimpleFoam solver
  #1
New Member
 
Xu Yiming
Join Date: Jan 2017
Posts: 11
Rep Power: 2
ousoweak is on a distinguished road
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
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
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
Code:
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
Code:
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
Code:
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)
Attached Files
File Type: zip 0.zip (3.9 KB, 0 views)
File Type: zip system.zip (2.7 KB, 2 views)
File Type: zip constant.zip (1.7 KB, 0 views)
ousoweak is offline   Reply With Quote

Reply

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
Fluent Adjoint Solver for compressible flow b.shuvayan FLUENT 2 January 30, 2015 08:27
Solver selection of supersonic compressible flow Leeoo Main CFD Forum 4 December 18, 2014 23:02
Star cd es-ice solver error ernarasimman STAR-CD 2 September 12, 2014 00:01
Solver for steady state compressible multispecies non-reacting flow? inf.vish OpenFOAM Running, Solving & CFD 0 October 1, 2013 01:09
Solver for compressible external flow ryan_m OpenFOAM Running, Solving & CFD 3 April 7, 2008 08:18


All times are GMT -4. The time now is 11:38.