CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

combustion third order time accuracy using runge kutta method

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 10, 2021, 03:41
Default combustion third order time accuracy using runge kutta method
  #1
New Member
 
sanjeev adhikari
Join Date: May 2018
Posts: 13
Rep Power: 7
sanjeev_adhikari is on a distinguished road
Hello everyone,


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.



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
equations:

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)
+fvc::laplacian(k,T)+reaction->Sh());

OF2.3.1


courant number and timesteps are calculated as per rhocentralFoam.






Thank you all.
sanjeev_adhikari is offline   Reply With Quote

Reply

Tags
combustion, openfoam, reaction flow, runge-kutta, specie

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
How to export time series of variables for one point? mary mor OpenFOAM Post-Processing 8 July 19, 2017 11:54
pressure in incompressible solvers e.g. simpleFoam chrizzl OpenFOAM Running, Solving & CFD 13 March 28, 2017 06:49
Moving mesh Niklas Wikstrom (Wikstrom) OpenFOAM Running, Solving & CFD 122 June 15, 2014 07:20
mixerVesselAMI2D's mass is not balancing sharonyue OpenFOAM Running, Solving & CFD 6 June 10, 2013 10:34
same geometry,structured and unstructured mesh,different behaviour. sharonyue OpenFOAM Running, Solving & CFD 13 January 2, 2013 23:40


All times are GMT -4. The time now is 07:01.