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

compressibleMultiphaseInterFoam rho solver in fvSolution does not do anything!

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

LinkBack Thread Tools Search this Thread Display Modes
Old   July 9, 2019, 12:13
Default compressibleMultiphaseInterFoam rho solver in fvSolution does not do anything!
New Member
Join Date: Jun 2012
Posts: 27
Rep Power: 13
rsa is on a distinguished road

I am trying to learn a bit about the compressibleMultiphaseInterFoam solver to use it later on and build my own case. At the moment I copied the OF6 tut-case (damBreak) with a new geometry with an inlet at the bottom to inject air into a mixture of water and oil with air at the top. I managed to run the case without any problem but there were couple of issues:

1- The results looked very much like the results of incompressible solver multiphaseInterFoam!
2- I looked at the run-time output of the solver and saw the following:
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
This tells me that the solver does not calculate the density which means the density should have been constant; hence, the results looking like the incompressible case now make sense (which is veery confusing since this is suppose to be incompressible)!
3- Then to test further, I commented out the solver for rho in the fvSolution as follows:
//solver diagonal;
solver banana;

what I expected was for the solver to crash and complain BUT to my supersize, it continued withou any issue!!!
The only thing I could conclude from this is the fact that the solver acts as incompressible and does not calculate the density which according to the info in the CFD Direct webpage the density should be calculated as rho=p/(RT).

I would appreciate if someone can let me know what I am doing wrong or maybe forgetting to do.

In case anyone wants to know the BCs I am using:
U: Inlet flowRateInletVelocity 0.001 / outlet pressureInletOutletVelocity
p: initialized to 1e5 / all the boundaries calculated
p_rgh: initialized to 1e5 / inlet fixedFluxPressure 1e6 / outlet totalPressure 1e5
T: initialized to 293 / inlet fixedValue 293 / outlet inletOutlet (inletValue 293) / walls fixedValue 293 [I intentionally set the system in a isothermal state for simplicity]

rsa is offline   Reply With Quote


compressible flow, density, fvsolution, interfoam, vof

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
rho calculating and updating in buoyantPimpleFoam calf.Z OpenFOAM Programming & Development 9 March 13, 2019 10:17
SLTS+rhoPisoFoam: what is rDeltaT??? nileshjrane OpenFOAM Running, Solving & CFD 4 February 25, 2013 04:13
Quarter Burner mesh with periosic condition SamCanuck FLUENT 2 August 31, 2011 11:34
Error log OpenFOAM 6 August 7, 2009 05:44
why the solver reject it? Anyone with experience? bearcat CFX 6 April 28, 2008 14:08

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