CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM

adding radiation to reactingFoam solver

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

Like Tree1Likes
  • 1 Post By megacrout

LinkBack Thread Tools Search this Thread Display Modes
Old   February 8, 2012, 20:49
Default adding radiation to reactingFoam solver
New Member
Join Date: Feb 2012
Posts: 25
Rep Power: 12
Marshak is on a distinguished road
I tried to add the radiation source term for reactingFoam in 'hs' equation and then included the header files for radiation in the reactingFoam.C and also modified the 'file' and 'options'. But when I try to run the case, everything runs fine except the radiation intensity 'G' is not calculated at each time step. My mesh is a sector mesh of a cylindrical domain and I have 'wedge' boundary condition for the periodic boundaries. Can anyone help me what can be the problem?

Has someone tried to add radiation in reactingFoam?

Last edited by Marshak; February 8, 2012 at 21:48.
Marshak is offline   Reply With Quote

Old   February 13, 2012, 06:43
Join Date: Jun 2011
Posts: 68
Rep Power: 13
megacrout is on a distinguished road
Hey Marshak,

I spent quite a lot of time playing around with reactingFoam and porousExplicitSourceReactingParcelFoam (on OF 1.7.1 though). You should have a look at the second one. It is an extension of reactingFoam and takes radiation, porosity and explicit (mass, energy and momentum) sources into account.

It can also handle multiple phases but I did not look much into it. I know however, that because of it, this solver is based on the thermo package icoPoly8 and that leads to pretty important discrepancies in the results if you donīt implement your data (very) carefully.

So, in short: I managed compiling a custom reactingFoam from pESRPF, that handles radiation and porous domains for one single (gaseous) phase using Janaf/Sutherland. However, I only ran simulations with a version in which I deactivated the radiation properties. Got reasonable results (not validated yet) but I obviously canīt say anything about radiation except that it can be done.

Your case sounds interesting and Iīd probably give your problem a shot but Iīm in the end phase of a project and just canīt do more than that at the moment.

Hope that helps (or at least relieves =D ).

babakflame likes this.
megacrout is offline   Reply With Quote


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
thobois class engineTopoChangerMesh error Peter_600 OpenFOAM 4 August 2, 2014 10:52
smoothSolver diverges - solution in using PBiCG solver? makaveli_lcf OpenFOAM Running, Solving & CFD 3 September 11, 2013 13:44
reactingFoam: ODE Solver KRR4 "Too many integration steps" pajofego OpenFOAM 3 May 18, 2013 16:15
Getting too many iterations by velocity solving (aborting). Changing U - Solver? suitup OpenFOAM Running, Solving & CFD 0 January 20, 2010 08:45
Error during Solver cfd guy CFX 4 May 8, 2001 07:04

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