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

Energy equation terms for reacting flows

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 17, 2019, 03:59
Default Energy equation terms for reacting flows
  #1
New Member
 
Pavan Inguva
Join Date: Dec 2019
Location: Singapore
Posts: 7
Rep Power: 6
pavaninguva is on a distinguished road
Hi!

I am trying to implement a model for a single phase flow reaction with an energy equation to model non-isothermal effects. But the formulation of the energy equation in the system is something I am trying to figure out.

The model thus far:
Navier-Stokes:

\nabla \cdot \bold{u} = 0

\frac{\partial}{\partial t}(\bold{u}) + \nabla \cdot (\bold{u} \otimes \bold{u}) - \nabla \cdot (\nu \nabla\bold{u}) = -\nabla p

Species transport:
\frac{\partial C_{i}}{\partial t} + \bold{u}\cdot \nabla C_{i} = D_{ij}\nabla^{2}C_{i} + R_{i}

Energy Equation:

This is where I am a bit unsure. I have seen a formulation energy equation that accounts a term for "energy flux by molecular diffusion" (Froment & Bischoff) which is the last term in the equation below:

\sum_{j}MW_{j} \ C_{j} \ c_{p,j}\bigg( \frac{\partial T}{\partial t} + \bold{u}\cdot \nabla T \bigg) = \nabla \cdot (\lambda \nabla T) + \sum_{i}(-\Delta H_{i})r_{i} - \sum_{j}\bold{J}_{j}\nabla H_{j}

But I am more familiar with the following formulation which has 1) transient term, 2) advective term, 3) diffusive term and 4) source term which would be the heat of reaction:

\frac{\partial T}{\partial t} + \bold{u}\cdot \nabla T = \alpha \nabla^{2}T + S

So would the additional term (energy flux by molecular diffusion) be necessary for a good model for such a reactive flow problem.
pavaninguva is offline   Reply With Quote

Old   December 17, 2019, 07:20
Default
  #2
Senior Member
 
Eifoehn4's Avatar
 
-
Join Date: Jul 2012
Location: Germany
Posts: 184
Rep Power: 14
Eifoehn4 is on a distinguished road
The answer to your question depends on the problem you want to solve.

Your heat flux due to species diffusion in your energy equation is a very simple approximation of Onsager reciprocal relations, namely the Dufour effect, based on non-equilibrium thermodynamics.

What problem do you want to solve?

Regards
__________________
Check out my side project:

A multiphysics discontinuous Galerkin framework: Youtube, Gitlab.

Last edited by Eifoehn4; December 17, 2019 at 14:44.
Eifoehn4 is offline   Reply With Quote

Old   December 17, 2019, 13:48
Default
  #3
Senior Member
 
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,747
Rep Power: 66
LuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura about
The diffusion of energy is more generally written as: \nabla \cdot \bold{q}
where q is a generalized diffusive energy flux. For homogeneous media, q is only heat diffusion from Fourier's law: \bold{q}=-k \nabla T
You can of course use, \alpha = \frac{k}{\rho c_p}, to write it in terms of thermal diffusivity instead of thermal conductivity. For multi-component problems, there is diffusion of heat via concentration gradients, which is the last term in your equation.

If you derive your energy equation assuming diffusion of heat is only via heat conduction/diffusion then you end up with: \frac{\partial T}{\partial t} + \bold{u}\cdot \nabla T = \alpha \nabla^{2}T + S

The source is a fudge factor because you didn't account for any possible chemical reactions or alternative forms of energy transport. Note that when the equation is written in this way, S is a temperature source. It's not even an energy source!

If you properly derive an energy transport equation starting with conservation of total energy then you should get the more complicated looking thing you showed:
\sum_{j}MW_{j} \ C_{j} \ c_{p,j}\bigg( \frac{\partial T}{\partial t} + \bold{u}\cdot \nabla T \bigg) = \nabla \cdot (\lambda \nabla T) + \sum_{i}(-\Delta H_{i})r_{i} - \sum_{j}\bold{J}_{j}\nabla H_{j}

The heat of reaction is contained in the term: \sum_{i}(-\Delta H_{i})r_{i}

The contribution of heat diffusion via temperature gradients vs species gradients depends on the problem. If the reaction is nearly isothermal, then species diffusion will dominate. If the reaction is highly exothermic, then thermal diffusion dominates.

It may also be important to mention that these are the time-accurate equations (which means you are doing DNS, whether laminar or turbulent). If you do time/Reynolds-averaging, the contribution of turbulent transport also needs to be taken into account and this gets even more wonky.

Last edited by LuckyTran; December 17, 2019 at 21:17.
LuckyTran is offline   Reply With Quote

Old   December 27, 2019, 03:35
Default
  #4
New Member
 
Pavan Inguva
Join Date: Dec 2019
Location: Singapore
Posts: 7
Rep Power: 6
pavaninguva is on a distinguished road
Thanks for your really valuable insights regarding the Dufour effects. Fortunately it seems that neglecting the Dufour effect is not a bad assumption for many cases of reacting flows.

For others' reference, I found that this book was great for outlining the derivation of the model equations and also discusses the Dufour effects with some detail: https://www.springer.com/gp/book/9783540686224.
pavaninguva is offline   Reply With Quote

Reply

Tags
energy balance, reacting flows, species tranport


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
Energy equation in rhoCentralFoam (revisited) usv001 OpenFOAM Programming & Development 4 September 6, 2024 08:51
Porous Modeling of Energy equation in OpenFOAM mohammad_kordo OpenFOAM Running, Solving & CFD 9 November 22, 2020 08:18
Fluent do not use my velocity field(by UDF) to solve energy equation tangleiplus Fluent UDF and Scheme Programming 6 January 21, 2019 22:28
rhoCentralFoam: making sense of the diffusive corrector equation for energy tatu OpenFOAM Running, Solving & CFD 5 February 19, 2013 09:19
how to disable energy equation for for non-newton yueroo FLUENT 2 April 3, 2001 23:11


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