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

OpenFoam Joule Thomson Simulation

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

Like Tree1Likes
  • 1 Post By user_of_cfx

LinkBack Thread Tools Search this Thread Display Modes
Old   January 4, 2016, 06:53
Default OpenFoam Joule Thomson Simulation
New Member
Join Date: Oct 2015
Posts: 2
Rep Power: 0
chisoph is on a distinguished road
Hello my friends,

I'm new in OpenFoam and I'm now working 3 months on this problem and I can't find a solution for my problem in the web. So I thought I ask the elite of OpenFoam.

I would like to simulate the Joule-Thomson effect in a pipe with a sudden expansion because of a diameter change of the pipe. The simulation starts and runs till it reaches the outlet. Then suddenly the simulation crashes. The picture shows you the last step of the simulation before it crashes. I really tried to solve the problem by my own, but I have no further idea....

Attached you can find my case setup.

I use
Solver: rhoPimpleFoam
Mesh: 2D-3D makeAxialMesh
sutherland transport
janaf thermo
PengRobinsonGas equation of State

If you need more information please don't hesitate to ask.

Thanks for your help

Attached Images
File Type: png Screenshot from 2016-01-04 12_22_35.png (148.0 KB, 94 views)
Attached Files
File Type: zip Simulation (13.8 KB, 26 views)
chisoph is offline   Reply With Quote

Old   January 12, 2016, 10:47
Join Date: Apr 2011
Posts: 53
Rep Power: 15
user_of_cfx is on a distinguished road
I think one piece of additional information which would help would be your simulation log, so people can see the error you are getting and what led up to it.

I am no expert by any means, but I noticed you are using a very small number of iterations in fvSolution: nOuterCorrectors 1; It may be you have convergence issues because of this? I would increase that number by a lot, maybe above 250 or even higher. Then you can monitor the cumulative error in the log to check that it does not increase too much as then the simulation will crash.

Also, what is your velocity doing? Does it reach Mach = 1 at any point? if so maybe sonicFoam is a better solver for this? Again, no expert, this is just an idea.

Btw is this J-T cooling or heating you're modelling? because the T fluctuation all the way from 200 K up to 800-odd K is strange for gas expansion starting at 295K. Not sure if this is physically possible. So I am inclined to think that either a. your simulation does not converge or b. your thermo properties may not be correct to capture J-T.

One final thing, you have not fixed the pressure anywhere in your system?


chisoph likes this.
user_of_cfx is offline   Reply With Quote

Old   January 12, 2016, 12:01
New Member
Join Date: Oct 2015
Posts: 2
Rep Power: 0
chisoph is on a distinguished road
Hi Christa,

Thank you very much for helping me.

Attached you can find the last entries of the Simulation Log. The Warnings based on JANAF also appear at the beginning of the Simulation. This Warnings don't lead to the crash I think.

I will try to increase the nOuterCorrectos next time.

Yes I do reach Mach > 1. As I said, I'm new in OpenFoam and I thought rhoPimpleFoam sounds like it could work. But I will try it using the SonicFoam Solver.

The gas expansion should lead to a Joule-Thomson cooling. I have no idea where those high values in the temperature come from.

What do you mean with fixing the pressure? The inletflow is controlled over the velocity. Relating to this velocity, the pressure is calculated. I want to simulate the effect of Joule-Thomson in a laminar, a laminar-turbulent-transition and a turbulent flow. Thats why I want to control the simulation over the velocity.


greets chisoph
Attached Files
File Type: txt Simulation Log.txt (25.5 KB, 9 views)
chisoph is offline   Reply With Quote

Old   January 15, 2018, 07:03
New Member
Ahmed Shoeb
Join Date: Sep 2017
Location: Germany
Posts: 3
Rep Power: 8
Ahmedshoeb is on a distinguished road
Hello chisoph,

Have you found a solution to your problem ?
i am also working on simulating the Joule Thompson effect and i am having some troubles setting up the case.
Ahmedshoeb is offline   Reply With Quote

Old   February 26, 2021, 12:25
Join Date: Mar 2019
Posts: 81
Rep Power: 7
mm66 is on a distinguished road
Hi Foamers,

Any recommendations for the solver? I am having a hard time finding a suitable solver for the refrigerant expansion with the JT effect... Any help is very much appreciated.

mm66 is offline   Reply With Quote


joule thomson effect, pengrobinsongas, rhopimplefoam

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
Map of the OpenFOAM Forum - Understanding where to post your questions! wyldckat OpenFOAM 10 September 2, 2021 05:29
OpenFoam Airfoil simulation - replacing Fluent with OpenFOAM shereez234 OpenFOAM Running, Solving & CFD 1 November 3, 2015 03:54
OpenFOAM Training, London, Chicago, Munich, Sep-Oct 2015 OpenFOAM Announcements from Other Sources 2 August 31, 2015 13:36
running multiple parallel cases in openfoam slows the simulation kkpal OpenFOAM Running, Solving & CFD 2 August 21, 2015 11:08
diesel Engine simulation in OpenFOAM karam OpenFOAM Running, Solving & CFD 1 March 1, 2011 09:46

All times are GMT -4. The time now is 23:34.