rhoCentralFoam quitting after 1 second

Glued_To_Wonder · June 15, 2017, 17:59

So I am doing a mach tube problem and I have to have pressure, Temperature, and Velocity all initialized non-uniformly so that there are two initial quantities for those properties before and the after the shock in the middle of the tube. However the rhoCentralFoam solver quits out at an execution time of about 1~2 seconds doing maybe 20 timesteps in total.

The odd thing is too is that rhoCentralFoam solves the solution completely when U is uniform which I figured out through testing.

Has anyone run into this before or know what's up?

I am using a vanLeer discretization scheme (vanLeerV for Velocity) with a CTL number of 0.1, time step of 0.000001 and a max time of 0.002 if that helps!

ano · June 22, 2017, 06:07

Hello Cody,

what error does OF give you? A floating point exception? In which equation solver did it appear?

What are your boundary conditions especially at the inlet? I also had some problems and now avoid them by switching from a supersonic to a transsonic or sometimes even incompressible inlet boundary condition. It is unphysical and not a solution in the long run, but for my cases I obtain reasonable results.

I guess you probably already looked at the setup on this homepage:
https://curiosityfluids.com/2016/07/...hocentralfoam/

June 15, 2017, 17:59	rhoCentralFoam quitting after 1 second	#1
Glued_To_Wonder New Member Cody Knoblock Join Date: Jun 2017 Posts: 1 Rep Power: 0	So I am doing a mach tube problem and I have to have pressure, Temperature, and Velocity all initialized non-uniformly so that there are two initial quantities for those properties before and the after the shock in the middle of the tube. However the rhoCentralFoam solver quits out at an execution time of about 1~2 seconds doing maybe 20 timesteps in total. The odd thing is too is that rhoCentralFoam solves the solution completely when U is uniform which I figured out through testing. Has anyone run into this before or know what's up? I am using a vanLeer discretization scheme (vanLeerV for Velocity) with a CTL number of 0.1, time step of 0.000001 and a max time of 0.002 if that helps!

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
How to apply a turbulence model to the rhoCentralFoam solver?	arussell92	OpenFOAM Pre-Processing	18	July 25, 2022 04:26
Modify rhoCentralFoam: other equations of state	fivos	OpenFOAM Programming & Development	5	July 29, 2020 13:17
Energy equation in rhoCentralFoam (revisited)	usv001	OpenFOAM Programming & Development	3	July 5, 2018 05:52
dynamic mesh refinement and rhoCentralFoam	ChrisA	OpenFOAM Running, Solving & CFD	1	March 21, 2013 08:00
rhoCentralFoam solver with Slip BCs fails in Parallel Only	JLight	OpenFOAM Running, Solving & CFD	2	October 11, 2012 21:08

June 22, 2017, 06:07		#2
ano Member ano Join Date: Jan 2017 Location: Delft Posts: 58 Rep Power: 10	Hello Cody, what error does OF give you? A floating point exception? In which equation solver did it appear? What are your boundary conditions especially at the inlet? I also had some problems and now avoid them by switching from a supersonic to a transsonic or sometimes even incompressible inlet boundary condition. It is unphysical and not a solution in the long run, but for my cases I obtain reasonable results. I guess you probably already looked at the setup on this homepage: https://curiosityfluids.com/2016/07/...hocentralfoam/