CFD Online URL
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

ODE solver for reactingFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 5, 2011, 12:10
Default ODE solver for reactingFoam
  #1
New Member
 
Daero Joung
Join Date: May 2009
Location: POSTECH, Republic of Korea
Posts: 3
Rep Power: 7
narak is on a distinguished road
I ran this test case with simple one step chemistry using different ODE Solvers(SIBS, KRR4, RK). I expected same result but it wasn't.
Maximum temperature is about 2252 K for SIBS, 1446 K for KRR4 and 2570 K for RK.

test case : ReactingFoam/counterFlowFlame2D

odeCoeffs
{
ODESolver RK; //SIBS, KRR4, RK
eps 0.05;
scale 1;
}

================================================== ================

But Results for simple ODE show no difference.

y" - y = 0 ; y(0) = 0, y'(0) = 1

SIBS : Numerical: y(1.0) = 2(1.175190 1.543074), hEst = 1.245250
KRR4 : Numerical: y(1.0) = 2(1.175212 1.543091), hEst = 0.186009
RK : Numerical: y(1.0) = 2(1.175180 1.543062), hEst = 0.665943

- ODEtest(test/ODETest.c)
http://www.tfd.chalmers.se/~hani/kurser/OS_CFD_2008/ZongyuanGu/reportZongyuan.pdf

- Cash-Karp Runge-kutta method in OpenFOAM
http://ftp.eq.uc.pt/books/NumericalRecipes/f16-2.pdf

Last edited by narak; April 7, 2011 at 13:02.
narak is offline   Reply With Quote

Old   April 13, 2011, 05:28
Default 15 step chemistry results
  #2
New Member
 
Daero Joung
Join Date: May 2009
Location: POSTECH, Republic of Korea
Posts: 3
Rep Power: 7
narak is on a distinguished road
I ran the same tutorial test case with 39 step chemistry(C7H16) included in OF and got almost same results.


reactingFoam SIBS KRR4 RK
T_max(K) 1892 1896 1901
dQ_max 0.128 0.130 0.133
*Ext t(s) 1818 1664 54497

* execution time for Time = 1

I don't know why 1 step chemistry showed different results yet, but it seems to run correctly with multi step chemistry including GRI 3.0.
narak is offline   Reply With Quote

Old   October 24, 2012, 16:34
Default
  #3
New Member
 
Mhsn
Join Date: Oct 2012
Posts: 24
Rep Power: 4
mhsn is on a distinguished road
Regarding your post, you are familiar with using detailed reaction mechanisms (like GriMech 3.0) in OpenFoam. What I'm going to do is to simulate premixed combustion in a small channel using OpenFoam2.1.1. I wonder what solver to use and how to implement full chemistry. I mean, how I can solve conservation equations for all species. Any help would be appreciated.
mhsn is offline   Reply With Quote

Old   November 16, 2012, 04:53
Default
  #4
New Member
 
Traib
Join Date: Sep 2012
Posts: 27
Rep Power: 4
Traib is on a distinguished road
Quote:
Originally Posted by mhsn View Post
What I'm going to do is to simulate premixed combustion in a small channel using OpenFoam2.1.1. I wonder what solver to use and how to implement full chemistry. I mean, how I can solve conservation equations for all species. Any help would be appreciated.
Hi mhsn,

You can use reactingFoam, or a density based solver rhoRactingFoam. These are solver application for which the source codes can be easily found here http://www.openfoam.org/docs/cpp/. However, if you meant to ask for a chemistry solver then you can select from the chemestry models listed in the above site.

As for solving the species transport and temperature equations, the solver applications you choose will do it for you.
Traib is offline   Reply With Quote

Old   August 18, 2013, 10:13
Default ODEsolver
  #5
Member
 
赵庆良
Join Date: Aug 2013
Posts: 56
Rep Power: 3
zqlhzx is on a distinguished road
Hi:
Foamers,Can someone tell me the way of ODE solvers what equation in reactingFoam?Is UEqn or hEqn?
zqlhzx is offline   Reply With Quote

Reply

Tags
ode, reactingfoam, solver

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
ODE solver coefficients mrangitschdowcom OpenFOAM Running, Solving & CFD 2 August 18, 2013 23:16
Working directory via command line Luiz CFX 4 March 6, 2011 21:02
Problem implementing CVODE ODE solver markusrehm OpenFOAM 20 October 13, 2010 18:02
Creating New Solver: For particle-laden compressible jets sankarv OpenFOAM 0 April 4, 2010 19:06
why the solver reject it? Anyone with experience? bearcat CFX 6 April 28, 2008 15:08


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