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

rhoCentralFoam, blunt body

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 5, 2015, 15:19
Default rhoCentralFoam, blunt body
  #1
New Member
 
Join Date: Apr 2015
Posts: 3
Rep Power: 3
stack is on a distinguished road
Hello all,

I using rhoCentralFoam to run a blunt body, at Mach 5.

The simulation outputs the error below after five timesteps, with a max Courant number of 0.046. This exact setup does not crash using sonicFoam, but I donít think that solver is appropriate for this situation.

The boundary conditions on the inlet/outlets are waveTransmissive, for both p and U, and zeroGradient for T.
The BCs for the blunt are zeroGradient for p and T, and slip for U.

I have also attached pictures of the five timesteps, at the upper corner of the body, where I believe the discontinuity is causing it to crash.

If anyone could help me troubleshoot this problem, it would be greatly appreciated, thanks!

Quote:
#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::sqrt(Foam::Field<double>&, Foam::UList<double> const&) at ??:?
#4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
#5
at ??:?
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7
at ??:?
Floating point exception (core dumped)
Attached Images
File Type: jpg 2.jpg (15.6 KB, 17 views)
File Type: jpg 3.jpg (16.0 KB, 12 views)
File Type: jpg 4.jpg (16.3 KB, 9 views)
File Type: jpg 5.jpg (16.5 KB, 10 views)
File Type: jpg 6.jpg (17.1 KB, 10 views)
stack is offline   Reply With Quote

Old   April 8, 2015, 12:48
Default
  #2
New Member
 
Join Date: Apr 2015
Posts: 3
Rep Power: 3
stack is on a distinguished road
I got it working.

I changed fvSchemes, in the "interpolationSchemes" section from:

Code:
reconstruct(rho) vanLeer;
reconstruct(U) vanLeerV;
reconstruct(T) vanLeer;
to:

Code:
reconstruct(rho) upwind;
reconstruct(U) upwind;
reconstruct(T) upwind;
stack is offline   Reply With Quote

Old   May 6, 2015, 06:18
Default
  #3
New Member
 
Romania
Join Date: Dec 2014
Posts: 5
Rep Power: 3
yrganiri is on a distinguished road
Hi,

I am having the same problem with rhoCentralFoam. Can I see your files (maybe that will help me)?

The thing is: I am running a simulation on a NACA airfoil with sonicFoam (no errors there), but when I change to rhoCentralFoam I get exactly the same error as you:

Quote:
Mean and max Courant Numbers = 0.00232163 0.295612
deltaT = 3.89625e-05
Time = 0.000532113

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
ExecutionTime = 0.78 s ClockTime = 1 s

#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::sqrt(Foam::Field<double>&, Foam::UList<double> const&) at ??:?
#4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
#5
at ??:?
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7
at ??:?
Floating point exception
Actually, did you start from an OpenFoam tutorial?

Last edited by yrganiri; May 6, 2015 at 06:19. Reason: adding information
yrganiri is offline   Reply With Quote

Reply

Tags
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
Carbuncle phenomenon for blunt body flow on unstructured grid praveen Main CFD Forum 2 July 30, 2015 10:11
BC for measuring force on blunt body letthe OpenFOAM Pre-Processing 4 December 9, 2013 19:20
Vortex Shedding on Profiled Blunt Trailing Edge Body Panboy FLUENT 7 May 14, 2012 13:34
Supersonic flow over a flat plate and a blunt body. Maddy FLUENT 0 August 31, 2010 16:35
Grid singularity in blunt body simulations Stefan Main CFD Forum 8 August 4, 2001 22:33


All times are GMT -4. The time now is 02:59.