|
[Sponsors] |
compressible flow calculation error when using rhoSimpleFoam solver |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#1 |
New Member
Xu Yiming
Join Date: Jan 2017
Posts: 11
Rep Power: 9 ![]() |
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: ![]() #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; } } 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; } } 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; } } 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; } } 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; } } |
|
![]() |
![]() |
![]() |
![]() |
#2 |
New Member
Join Date: Aug 2017
Posts: 16
Rep Power: 9 ![]() |
Dear Xu,
I have met the same problem, how do you sloved this problem? |
|
![]() |
![]() |
![]() |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Fluent Adjoint Solver for compressible flow | b.shuvayan | FLUENT | 2 | January 30, 2015 07:27 |
Solver selection of supersonic compressible flow | Leeoo | Main CFD Forum | 4 | December 18, 2014 22: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 |