[Azanulbizar]

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.

[apu111]

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!

[demichie]

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.

[apu111]

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!

[reza_65]

Originally Posted by demichie

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

[saturn_K]

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

[dlahaye]

Use reactingFoam with combustion de-activated?

[saturn_K]

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

[clapointe]

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

[dlahaye]

Originally Posted by saturn_K

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.

[saturn_K]

Originally Posted by clapointe

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

[clapointe]

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

[saturn_K]

Originally Posted by clapointe

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.

[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

[clapointe]

In short, you'll notice that no tutorials for rhoCentralFoam use fixedFluxPressure for walls -- for high speed compressible problems zero gradient is sufficient.

Caelan

[Marpole]

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]

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.

[Marpole]

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.