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

Conjugate heat transfer and reactingFoam (chtMultiRegionReactingFoam)

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

Like Tree10Likes
  • 7 Post By edaymo
  • 1 Post By melabbassi
  • 2 Post By Swagga5aur

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 29, 2016, 22:19
Default Conjugate heat transfer and reactingFoam (chtMultiRegionReactingFoam)
  #1
New Member
 
Eric Daymo
Join Date: Feb 2015
Location: Gilbert, Arizona, USA
Posts: 20
Rep Power: 7
edaymo is on a distinguished road
Hello,

If anyone is interested in a solver with conjugate heat transfer + gas phase reaction capabilities, I've posted to github a hybrid of chtMultiRegionFoam and reactingFoam called chtMultiRegionReactingFoam. Steady state and transient solvers for OpenFOAM 4 and the current OpenFOAM-dev version are available (along with tutorials) at:

https://github.com/TonkomoLLC

I hope this is useful for your work. Feel free to contact me if you have any questions or comments.

Best regards,

Eric D.
Tonkomo, LLC
Twitter: @TonkomoLLC,
email: info at tonkomo.com
edaymo is offline   Reply With Quote

Old   August 30, 2016, 12:15
Default Thanks!
  #2
Member
 
Arvind Jay
Join Date: Sep 2012
Posts: 89
Rep Power: 10
arvindpj is on a distinguished road
I was waiting for this feature. I will test and report.
-Jay
arvindpj is offline   Reply With Quote

Old   February 27, 2017, 10:57
Default
  #3
New Member
 
Mohamed el Abbassi
Join Date: Oct 2016
Location: Delft, the Netherlands
Posts: 6
Rep Power: 5
melabbassi is on a distinguished road
Hi Eric,

Thank you for your contribution. I have studied the solver (OpenFOAM 3) this month and benchmarked the implementation of reaction in the cht-version against rhoReactingFoam which is also based on p_rgh. So only 1 region and no conjugate heat transfer. As a test case I modelled the Sandia D flame with 2-step reaction mechanism of Westbrook and Dryer.

With the steady state version I initially couldn't get ignition. Even after adjusting the relaxation factors and setting the initial internal field temperature to 2000 K. The maximum temperature immediately drops to inlet values. I found out that the convention with the SIMPLE algorithm to use a timestep of 1 s was the issue. Unlike the navier stokes equations for steady state problems, the energy and species equations are largely dependent on time, which is probably taken from the 'virtual' timestep. Choosing high timesteps allows little energy release. By adjusting the time step to a much lower value (so that max Courant nr is about 0.2), reaction was finally possible.

Attached you'll find a plot of the temperature progression along the rotation axis. Maximum temperature and temperature rise are similar, but reaction starts a little bit more upstream with chtMultiRegionSimpleFoam. Still, it is a very good result.

The transient version of the solver is less robust. I needed 4 nOuterCorrectors instead of 1 with rhoReactingFoam to get the simulation going, leading to 4 times longer runtime. Can it be due to the pimple algorithm that is implemented differently?

Best,
Mohamed
Attached Images
File Type: png T_plot.png (55.2 KB, 66 views)
arvindpj likes this.
melabbassi is offline   Reply With Quote

Old   February 27, 2017, 16:57
Default
  #4
New Member
 
Eric Daymo
Join Date: Feb 2015
Location: Gilbert, Arizona, USA
Posts: 20
Rep Power: 7
edaymo is on a distinguished road
Hi, Mohamed,

Thank you very much for such a detailed review of chtMultiRegionReactingFoam! This is great feedback, especially on your solution technique for the steady state solver and the benchmark findings.

rhoReactingFoam and chtMultiRegionReactingFoam (directly based on chtMultiRegionFoam for the PIMPLE algorithm) do indeed have different formulations for PEqn, but I do know know offhand if this the root cause of the nOuterCorrectors issue you are report. I am happy to look at your test cases in more detail if you wish to send them to me. To start with, I would want to repeat with OpenFOAM 4.x or OpenFOAM-dev. The OpenFOAM 3.x version you are using is a "beta" version of the chtMultiRegionReactingFoam solver and it was not as thoroughly tested as the versions on github.

I also note that since my initial August 2016 post I have placed on our github site two new solvers: "multiRegionReactingFoam" and "multiRegionReactingPimpleCentralFoam" (the latter based on https://github.com/unicfdlab/hybridCentralSolvers). I have tested multiRegionReactingFoam more extensively against reactingFoam and I would be very curious to know if this solver meets your expectations.

I look forward to further feedback. Thank you again for your insightful reply.

Best regards,

Eric
edaymo is offline   Reply With Quote

Old   March 1, 2017, 10:10
Default
  #5
New Member
 
Mohamed el Abbassi
Join Date: Oct 2016
Location: Delft, the Netherlands
Posts: 6
Rep Power: 5
melabbassi is on a distinguished road
Thanks for the clarity! That explains certain bugs I encountered in your beta version. I should have updated to OF 4 sooner, but it was good training exercise for me to be honest.

I will send you the case files for chtMultiRegionReactingFoam and rhoReactingFoam. Thank you very much for your effort! I will also run them again on OF 4.

The reason I used the p_rgh based solvers was because I wanted to make a consistent comparison and wasn't aware of the psi based multiRegionReactingFoam. Thanks for uploading it!

I'll keep you updated.
melabbassi is offline   Reply With Quote

Old   April 24, 2017, 14:13
Default Hello!
  #6
New Member
 
Jessica Smit
Join Date: Dec 2016
Posts: 1
Rep Power: 0
skibenochek is on a distinguished road
Thank you for your work! Can I use your solver with OpenFoam 2.3.0?
skibenochek is offline   Reply With Quote

Old   April 24, 2017, 15:14
Default Re: Hello!
  #7
New Member
 
Eric Daymo
Join Date: Feb 2015
Location: Gilbert, Arizona, USA
Posts: 20
Rep Power: 7
edaymo is on a distinguished road
Hi,

Thanks for the feedback!

The same programming principles apply for coding a CHT+reaction solver for OpenFOAM 2.3.0, but i think you will find that the solvers on my GitHub page will not compile as-is for OF2.3.

I've got a version of the CHT+reaction solver for OpenFOAM 3.0.x that I haven't gotten around to putting on GitHub yet, but I have not yet created one for OF 2.3/OF2.4.

Please contact me directly (e.g., private message, or info at tonkomo.com) if you would like my assistance to port this CHT+reaction solver to OF2.3. Thanks!

Best regards,

Eric
edaymo is offline   Reply With Quote

Old   October 2, 2017, 12:01
Default OFV5 multiregionsimplereactingfoam error
  #8
Member
 
Lasse Brams Vinther
Join Date: Oct 2015
Posts: 72
Rep Power: 6
Swagga5aur is on a distinguished road
Hello Eric

First of all I would like to thank you for this solver and the continuous work, as I used it for my master thesis concerning modelling a reformation process with external heating.

I noticed an error in the OF 5.0 multiRegionSimpleReactingFoam with respect to the createFluidFields.H dictionary missing the following line

Code:
    reactionFluid[i].setTurbulence(turbulence[i]);
after the following code

Code:
    Info<< "    Adding to turbulence\n" << endl;
    turbulence.set
    (
        i,
        compressible::turbulenceModel::New
        (
            rhoFluid[i],
            UFluid[i],
            phiFluid[i],
            reactionFluid[i].thermo()
        ).ptr()
    );
Prior to this addition I got a turbulence Ptr error when trying to solve a turbulence reaction model.

Just wanted to let you know

Best regards
Lasse
arvindpj and edaymo like this.
Swagga5aur is offline   Reply With Quote

Old   October 2, 2017, 16:48
Default Re: OFV5 multiregionsimplereactingfoam error
  #9
New Member
 
Eric Daymo
Join Date: Feb 2015
Location: Gilbert, Arizona, USA
Posts: 20
Rep Power: 7
edaymo is on a distinguished road
Dear Lasse,

Thank you very much for the bug report. Indeed, you are correct, and all of the multiRegionSimpleReactingFoam solvers (OF4, 5, dev, 1706, 1612+) have been updated with the reported fix to setTurbulence. I checked the other solvers (e.g., transient multiRegionReactingFoam, chtMultiRegionFoam, etc.) just to be sure, but I didn't find this bug outside of mutliRegionSimpleReactingFoam.

I really appreciate your bringing this matter to my attention. I'm also really happy to learn that this solver helped with your Master's thesis, and that you got good results with it.

Best regards,

Eric
edaymo is offline   Reply With Quote

Old   October 2, 2017, 17:15
Default
  #10
Member
 
Lasse Brams Vinther
Join Date: Oct 2015
Posts: 72
Rep Power: 6
Swagga5aur is on a distinguished road
Hello Eric,

No problem, glad to be of help.

I am experiencing some convergence issues with the multiRegionSimpleReactingFoam compared to the transient solvers multiRegionReactingFoam and reactingFoam.

Not sure if something else is missing in the solver but I'll look into it and post if I find any issues .

I believe the issue may be with the implementation of the OF5.0 EDC combustion model.

Best regards,
Lasse
Swagga5aur is offline   Reply With Quote

Old   October 2, 2017, 17:28
Default
  #11
New Member
 
Eric Daymo
Join Date: Feb 2015
Location: Gilbert, Arizona, USA
Posts: 20
Rep Power: 7
edaymo is on a distinguished road
Hello, Lasse -

I don't typically use the steady state solver, rather, I tend to use the transient solver with LTS.

One idea is to try the same problem with simpleParcelReactingFoam, which is a single region, steady state solver with chemistry, albeit with rhoReactingThermo. Assuming you can come up with reasonable BC's for a non-CHT case, this approach might narrow down whether the problem is with the multiRegion solver, or if the problem is with some other model or the setup.

Good luck and feel free to contact me if I can be of assistance.

Best regards,

Eric
edaymo is offline   Reply With Quote

Reply

Tags
chtmultiregionfoam, conjugate heat transfer, openfoam, reactingfoam, reaction

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
[swak4Foam] how can use Cp and Cv in Swak variables? immortality OpenFOAM Community Contributions 18 December 9, 2016 02:04
Multiregion heat transfer Habibfateh OpenFOAM Programming & Development 0 February 18, 2014 16:53


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