CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

combustion, rk4 method third for time accuracy

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

LinkBack Thread Tools Search this Thread Display Modes
Old   May 11, 2021, 22:03
Default combustion, rk4 method third for time accuracy
New Member
sanjeev adhikari
Join Date: May 2018
Posts: 13
Rep Power: 8
sanjeev_adhikari is on a distinguished road
Hello everyone,

(I had posted it in main CFD forum earlier )

This is the follow up post of my earlier thread as i didn't find error after long check out. However i realized that there is no combustion(no reaction and no heat release )in the simulation so i am getting lower temperature. Combustion was working fine before upgrading time accuracy using rk4 method.

Also i checked by writing reaction rate independently it shows the concentration for other species too. Also, it concentration of initial species is computed which shows the transport of the species is working. But in overall combustion is not taking place even when the flow temperature at the stagnation region reaches more than 2700K.

i checked chemistry properties and combustion both are on and true. Also i am using "rhoCombustionModel" . can anyone suggest possible solution or error on the fomulation of species and energy equation specially involving reaction. Thank you.

I was facing the problem with time accuracy of the solution using first order Euler ddt. So, i thought of implementing rk4 third order time accuracy. Following represents the

rho = rhoOld + rkCoef[cycle] *runTime.deltaT()*(-fvc::div(Riemann.rhoFlux()))

momentum equations.
rhoU = rhoUOld + rkCoeff[cycle]*runTime.deltaT()*(-fvc::div(Riemann.rhoUFlux())+fvc::div(tauMC)+fvc:: laplacian(muEff, U));

species transport eqn:
fvScalarMatrix YiEqn
rhoYi - rhoYOldi
- rkCoeff[cycle]*runTime.deltaT()*(-fvc::div(Riemann.rhoScalarFlux(i))+
fvm::laplacian(turbulence->muEff(), Yi))
== rkCoeff[cycle]*runTime.deltaT()*reaction->R(Yi)

energy eqn

rhoE = rhoEOld + rkCoeff[cycle]*runTime.deltaT()*(
--fvc::div(Riemann.rhoEFlux() + fvc::div(sigmaDotU) + fvc::laplacian(turbulence->alphaEff(), h)


ps: courant number and timesteps are calculated as per rhocentralFoam.

Thank you all.
sanjeev_adhikari is offline   Reply With Quote


combustion, runge kuttta, species, time accuracy

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
Finite element method ANSYS CFX tomd97 CFX 2 April 4, 2018 04:10
On the alpha Eqn of VOF method when using Immersed boundary method in OpenFOAM keepfit OpenFOAM 4 January 31, 2014 14:32
multiphaseEulerFoam: method mag(...) withing method solveAlphase() maybee OpenFOAM Programming & Development 0 December 21, 2013 10:13
Hydrogen Air combustion in a combustion chamber popi CFX 7 July 11, 2007 18:40
Multiresolution + Turbulence + Combustion Raul CFX 0 March 18, 2005 07:53

All times are GMT -4. The time now is 19:28.