CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

RhoCentralFoam-MultiSpecies-NonReacting

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

Like Tree2Likes
  • 1 Post By Marpole
  • 1 Post By Marpole

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 23, 2013, 19:25
Default RhoCentralFoam-MultiSpecies-NonReacting
  #1
New Member
 
Jack Lee
Join Date: Jan 2013
Location: Alberta
Posts: 5
Rep Power: 13
Azanulbizar is on a distinguished road
Hi,

I need to adapt rhocentralFoam to do mulitspecies non reacting simulation. Basically, at the inlet, I would like to inject 2-3 types of gases.

Can anyone gives me a suggestion of how to go about this matter?

The rhoCentralFoam.C looks much more complicated compared to other solvers and hence, I dont have much clue to change the codes so that I can take in multispecies/gases I need to inject in.

Thank you very much.
Azanulbizar is offline   Reply With Quote

Old   March 8, 2015, 04:53
Default
  #2
New Member
 
apurva
Join Date: Jul 2014
Posts: 2
Rep Power: 0
apu111 is on a distinguished road
Hi Azanulbizar,
Have you solved this problem?
So is it possible to introduce 2-3 non-reacting gases mixture in rhoCentralFoam?
I am working on delaval nozzle and I want to introduce them at nozzle inlet.

Thanks in advance!
apu111 is offline   Reply With Quote

Old   March 9, 2015, 13:49
Default
  #3
Member
 
Mattia de\' Michieli Vitturi
Join Date: Mar 2009
Posts: 50
Rep Power: 17
demichie is on a distinguished road
Hi,
I have published a couple of years ago a paper on a new OpenFOAM solver we have developed for a two-phase model using a central scheme approach similar to that used in rhoCetralFoam:

La Spina, G. & de' Michieli Vitturi, M. High-Resolution Finite Volume Central Schemes for a Compressible Two-Phase Model SIAM Journal on Scientific Computing, SIAM, 2012, 34, 861-880

In a more recent work we have also introduced mass transfer in the model.
If you are interested you can send an email to the first author asking for the code.

Ciao
Mattia

P.S. I have attached a figure from one of the papers comparing the results of a simulation with those of a lab experiments for the interaction of a shock wave with a cylindrical bubble.
Attached Images
File Type: jpg bubbles.jpg (46.8 KB, 90 views)
demichie is offline   Reply With Quote

Old   March 31, 2016, 15:46
Default
  #4
New Member
 
apurva
Join Date: Jul 2014
Posts: 2
Rep Power: 0
apu111 is on a distinguished road
Thanks demichie for your reply. I have gone through the paper but it's for 2 phase flow. How do I edit rhoCentralFoam.C to simlulate multi-species in the gaseous phase? Thanks in advance!
apu111 is offline   Reply With Quote

Old   December 8, 2017, 12:19
Default
  #5
Member
 
Reza khodadadi
Join Date: Apr 2011
Location: https://t.me/pump_upp
Posts: 32
Rep Power: 14
reza_65 is on a distinguished road
Send a message via ICQ to reza_65 Send a message via AIM to reza_65 Send a message via Yahoo to reza_65
Quote:
Originally Posted by demichie View Post
Hi,
I have published a couple of years ago a paper on a new OpenFOAM solver we have developed for a two-phase model using a central scheme approach similar to that used in rhoCetralFoam:

La Spina, G. & de' Michieli Vitturi, M. High-Resolution Finite Volume Central Schemes for a Compressible Two-Phase Model SIAM Journal on Scientific Computing, SIAM, 2012, 34, 861-880

In a more recent work we have also introduced mass transfer in the model.
If you are interested you can send an email to the first author asking for the code.

Ciao
Mattia

P.S. I have attached a figure from one of the papers comparing the results of a simulation with those of a lab experiments for the interaction of a shock wave with a cylindrical bubble.
It is a very nice job, well done. In 2010 I developed a two-phase flow solver, by using interFoam, for simulating shock bubble interactions as well. Also, I published an article about shock capturing which shows rhoCentralFoam has got the best shock capturing capability in openFoam solvers.

https://www.researchgate.net/publica...ds_in_OpenFOAM

DOI: 10.2298/TSCI130425048K

Best Regards,
Reza
reza_65 is offline   Reply With Quote

Old   August 27, 2020, 11:10
Default
  #6
New Member
 
Join Date: Sep 2019
Posts: 14
Rep Power: 6
saturn_K is on a distinguished road
Hi Azanulbizar

Could you finally find your answer? because I have exactly the same problem with rhoCentralFoam solver and I would be really grateful if you can help me how can I handle it.

Regards,
Keivan
saturn_K is offline   Reply With Quote

Old   August 27, 2020, 13:57
Default
  #7
Senior Member
 
Domenico Lahaye
Join Date: Dec 2013
Posts: 720
Blog Entries: 1
Rep Power: 17
dlahaye is on a distinguished road
Use reactingFoam with combustion de-activated?
dlahaye is offline   Reply With Quote

Old   August 28, 2020, 11:49
Default
  #8
New Member
 
Join Date: Sep 2019
Posts: 14
Rep Power: 6
saturn_K is on a distinguished road
Hi dlahaye

Thank you for your suggestion. But as far as I know, reactingFoam solver doesn't have the ability to capture the discontinuities like shock waves.

Regards,
Keivan
saturn_K is offline   Reply With Quote

Old   August 28, 2020, 12:03
Default
  #9
Senior Member
 
Join Date: Aug 2015
Posts: 494
Rep Power: 14
clapointe is on a distinguished road
There are various flavors of reacting rhoCentralFoam floating around (e.g. reactingPimpleCentralFoam, hy2Foam). I've put one together based on similar threads for older versions of openfoam... it can be found in this repo : https://github.com/duncanam/rhoReactingCentralFoam.

Caelan
clapointe is offline   Reply With Quote

Old   August 30, 2020, 12:38
Default
  #10
Senior Member
 
Domenico Lahaye
Join Date: Dec 2013
Posts: 720
Blog Entries: 1
Rep Power: 17
dlahaye is on a distinguished road
Quote:
Originally Posted by saturn_K View Post
Hi dlahaye

Thank you for your suggestion. But as far as I know, reactingFoam solver doesn't have the ability to capture the discontinuities like shock waves.

Regards,
Keivan
Thank you for elaborating your question.
dlahaye is offline   Reply With Quote

Old   September 2, 2020, 22:17
Default
  #11
New Member
 
Join Date: Sep 2019
Posts: 14
Rep Power: 6
saturn_K is on a distinguished road
Quote:
Originally Posted by clapointe View Post
There are various flavors of reacting rhoCentralFoam floating around (e.g. reactingPimpleCentralFoam, hy2Foam). I've put one together based on similar threads for older versions of openfoam... it can be found in this repo : https://github.com/duncanam/rhoReactingCentralFoam.

Caelan
Dear Clapointe

Thank you for your great help. Actually, I tried to use your solver and I encountered many issues and fortunately, I could solve all except one. I need to use a non-ideal equation of state (e.g. Peng-Robinson), but all the thermo types in your solver are based on perfectGas EOS. Do you know how I can add a new thermo type to your solver?

Best Regards,
Keivan
saturn_K is offline   Reply With Quote

Old   September 3, 2020, 00:41
Default
  #12
Senior Member
 
Join Date: Aug 2015
Posts: 494
Rep Power: 14
clapointe is on a distinguished road
It is simply an extension of rhoCentralFoam, so it uses existing openfoam thermo libraries. You'll need to add anything that doesn't exist natively.

Caelan
clapointe is offline   Reply With Quote

Old   September 17, 2020, 06:56
Default
  #13
New Member
 
Join Date: Sep 2019
Posts: 14
Rep Power: 6
saturn_K is on a distinguished road
Quote:
Originally Posted by clapointe View Post
It is simply an extension of rhoCentralFoam, so it uses existing openfoam thermo libraries. You'll need to add anything that doesn't exist natively.

Caelan
Dear Clapointe

Still, I am trying to use your solver and so far I have had good progression and I should thank you for that.
But recently, I got a new error which is related to the "fixedFluxPressure" boundary condition.

Quote:
[3] --> FOAM FATAL ERROR:
[3] updateCoeffs(const scalarField& snGradp) MUST be called before updateCoeffs() or evaluate() to set the boundary gradient.
[3]
[3] From function virtual void Foam::fixedFluxPressureFvPatchScalarField::updateC oeffs()
[3] in file fields/fvPatchFields/derived/fixedFluxPressure/fixedFluxPressureFvPatchScalarField.C at line 157.
[3]
This happens after the first iteration. Have you ever encounter this error? Do you know how I can fix it?

Regards,
Keivan
saturn_K is offline   Reply With Quote

Old   September 17, 2020, 11:48
Default
  #14
Senior Member
 
Join Date: Aug 2015
Posts: 494
Rep Power: 14
clapointe is on a distinguished road
In short, you'll notice that no tutorials for rhoCentralFoam use fixedFluxPressure for walls -- for high speed compressible problems zero gradient is sufficient.

Caelan
clapointe is offline   Reply With Quote

Old   June 16, 2022, 02:08
Default
  #15
Senior Member
 
Charles
Join Date: Aug 2016
Location: Vancouver, Canada
Posts: 148
Rep Power: 9
Marpole is on a distinguished road
Hello Caelan,

I'm very interested in rhoReactingCentralFoam you developed 2 years ago. Is there any further development such as upgrading to a new version for OpenFoam v9?

Regards,
ET3 likes this.
__________________
Charles L.
Marpole is offline   Reply With Quote

Old   July 12, 2022, 15:49
Default rhoReactingCentralFoam
  #16
ET3
New Member
 
Erik T
Join Date: Jun 2022
Posts: 5
Rep Power: 3
ET3 is on a distinguished road
Charles,



Sadly, Caelen died earlier this year.



I too am interested in updating this solver. Have you tried doing so?



I am trying it right now and I am stuck. The updates to versions 8 changed things quite a bit. I'd be interesting in learning what you have tried and if you have success. I am new to this but so far it looks like there is more than a handful of changes and based on the threads you have been involved in, it seems like reactingFoam would be the replacement to rhoCentralFoam, but reactingFoam may not be the right solution for these problems.
ET3 is offline   Reply With Quote

Old   July 14, 2022, 00:37
Default
  #17
Senior Member
 
Charles
Join Date: Aug 2016
Location: Vancouver, Canada
Posts: 148
Rep Power: 9
Marpole is on a distinguished road
Hi Erik,
Thanks for letting me know. I visited Caelan's website. What a loss!

I haven't started updating the program. Instead, I installed openfoam 7 and can run the program to address my immediate problem. It solved the problem though slow.

I tried, on openfoam 9, reactingFoam. It has everything I need but it cannot solve shock wave and supersonic speed flow. So not much useful.

Thanks again for letting me know.
ET3 likes this.
__________________
Charles L.
Marpole is offline   Reply With Quote

Reply

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
rhoCentralFoam transport equation JoaoDMiranda OpenFOAM Programming & Development 28 January 15, 2020 09:12
dynamic mesh refinement and rhoCentralFoam ChrisA OpenFOAM Running, Solving & CFD 1 March 21, 2013 09:00
calculate a nozzle flow using rhoCentralFoam hg2lf OpenFOAM 0 October 25, 2012 22:26
rhoCentralFoam solver with Slip BCs fails in Parallel Only JLight OpenFOAM Running, Solving & CFD 2 October 11, 2012 22:08
laminar compressible solver for rhoCentralFoam (OF-2.0.x) turbulentwakes OpenFOAM Running, Solving & CFD 2 March 21, 2012 04:02


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