CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

twoPhaseEulerFoam: Exceeding iterations, residual = NaN

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

LinkBack Thread Tools Search this Thread Display Modes
Old   January 2, 2016, 13:18
Default twoPhaseEulerFoam: Exceeding iterations, residual = NaN
New Member
Join Date: Nov 2015
Posts: 6
Rep Power: 10
Flyingcircus is on a distinguished road
Hi Everyone

I've run into some trouble initializing a twoPhaseEulerFoam case.

Short description: I'm trying to eject a free water jet into open air at high velocities in order to simulate air entrainment.

The following picture, taken from a stable interFoam Simulation might help in understanding what I'm trying to achieve.

The water (red) enters the domain in a tube, then flows from left to right and passes through a slightly narrower blind, after which is traverses the bigger model area filled with air (blue) and exits the domain at the BC.

Sorry if I'm stating the obvious, I just don't want there to be any misunderstandings

Front and back BCs (in Z direction) are named "Boundary", utmost right and bottom BC are named "Drain" uppermost BC is named "Ceiling", left BC of the area where the water jet comes into contact with air is called "InletWall" (but is not actually a wall here), the utmost left BC, where water enters the domain, is named "source", the "DN150" and "Blind" BCs are for the tube and blind geometry respectively.

However when trying to run the same case in twoPhaseEulerFoam I'm running into the "NaN" issue right at the beginning.
I've been trying to juggle around with different BCs and linear solvers but had no luck so far.
I'm attaching the output file so that you can see what I'm talking about.

Furthermore I'm attaching all BC files, all "properties" files and the solution and schemes files.

Currently I'm running the case in a hexa dominant mesh, created with SHM. I've recently been forced to switch to 2.3.0 and have been working with 2.4.0 before that, which is why all of the files still bear the 2.4.0 header.

If any more info is required please let me know.

edit: I'm a bit worried whether twoPhaseEulerFoam is even the correct approach here or if I should be using a different solver.

edit1: reuploaded the archive since there was an error in the BC files which I had already fixed.
Attached Files
File Type: zip (12.0 KB, 16 views)

Last edited by Flyingcircus; January 3, 2016 at 10:39.
Flyingcircus is offline   Reply With Quote


Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
simpleFoam error - "Floating point exception" mbcx4jc2 OpenFOAM Running, Solving & CFD 12 August 4, 2015 02:20
pimpleFoam: turbulence->correct(); is not executed when using residualControl hfs OpenFOAM Running, Solving & CFD 3 October 29, 2013 08:35
calculation stops after few time steps sivakumar OpenFOAM Running, Solving & CFD 7 March 17, 2013 06:37
Orifice Plate with a fully developed flow - Problems with convergence jonmec OpenFOAM Running, Solving & CFD 3 July 28, 2011 05:24
Could anybody help me see this error and give help liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 18:07

All times are GMT -4. The time now is 06:41.