# how to extend solver compressibleInterFoam to consider ideal gas

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

 February 19, 2013, 21:15 how to extend solver compressibleInterFoam to consider ideal gas #1 New Member   jianfeng zou Join Date: Feb 2013 Posts: 4 Rep Power: 6 Hi, everyone. I am a beginner to use OpenFOAM and my interest is computing the interface between two immiscible compressible gases. It seems that current compressibleInterFoam solver is designed to solve 2 compressible, isothermal immiscible fluids. So if I could extend this solver to consider the interface between two or more ideal gases or real gas? Thanks in advance. I have no idea what should I do at all, please give me some advice. immortality likes this.

 February 20, 2013, 14:09 #2 Senior Member     Matvey Kraposhin Join Date: Mar 2009 Location: Moscow, Russian Federation Posts: 330 Rep Power: 12 Hello, Your question sounds very complex (for me). You are interested in two-phase flow of immiscible ideal-gases, but interFoam family solvers considers for liquids with large density variation. Are you sure that assumptions, made in interFoam solvers family are suitable for your task? If yes, then you can do modifications easily: 1) relate compressibility of each gas to temperature psi = 1/((R/M)*T) 2) rho = psi*p 3) write temperature equation, like this: fvm::ddt(T) + fvm::div(phi,T) - fvm::Sp(divU,T) == (... some expansion sources ...) I solved this case and if you need, then we need to examine your problem step by step

 February 20, 2013, 14:44 #3 New Member   jianfeng zou Join Date: Feb 2013 Posts: 4 Rep Power: 6 Thank Kraposhin. I just intend to calculate the interaction of a shock wave with cylindrical Helium gas inhomogeneities, like the experiment conducted by J.*F. Haas and B. Sturtevant (JOURNAL OF FLUID MECHANICS, 1987). A front tracking code was used to perform this simulation in my previous work. And I hope OpenFOAM will work on the same case. And which solver did you make the mentioned modifications, interFoam or compressibleInterFoam? Could you please give me a detailed description about where and how in the solver should be modified?

February 22, 2013, 12:25
#4
Senior Member

Matvey Kraposhin
Join Date: Mar 2009
Location: Moscow, Russian Federation
Posts: 330
Rep Power: 12
Have a nice day, jianfeng
I have good news for you - in my archives i found old version of non-isothermal compressibleInterFoam (see attachment).

i made test case (based on original) and it behaves well.

But! dp/dt works very strange (for me)
Attached Files
 nonIsothermalInterFoam.tar.gz (5.2 KB, 72 views) rhoDepthCharge2D.tar.gz (3.3 KB, 41 views)

Last edited by mkraposhin; February 22, 2013 at 12:28. Reason: i forgot to attach files with source code of solver and test case

 February 25, 2013, 02:35 #5 New Member   jianfeng zou Join Date: Feb 2013 Posts: 4 Rep Power: 6 Kraposhin, really thank you for your share. I have compiled the attached code successfully and I will test it with my own case. Hope it work fine.

 February 28, 2013, 14:43 OF version #6 New Member   Concordia_CFD Join Date: Jul 2010 Location: Canada Posts: 24 Rep Power: 9 Hi, Is the solver developed for OF-1.5 dev.? Thanks

March 1, 2013, 14:52
#7
Senior Member

Matvey Kraposhin
Join Date: Mar 2009
Location: Moscow, Russian Federation
Posts: 330
Rep Power: 12
Quote:
 Originally Posted by marzbali Hi, Is the solver developed for OF-1.5 dev.? Thanks
this solver was made for OF-2.1

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Mohsin FLUENT 7 June 7, 2015 22:23 RyGuy FLUENT 4 June 28, 2012 02:07 iggor OpenFOAM 15 October 14, 2009 11:58 saurav pathak Main CFD Forum 0 February 18, 2008 06:27 Björn Mattsson FLUENT 5 September 5, 2005 04:03

All times are GMT -4. The time now is 15:28.