[kasper1301]

Hi,

I am trying to implement an equation for chemical components in twoPhaseEulerFoam and I am wondering how to discretize each term in OpenFOAM.

The transport equation:


I was thinking that I could create some new fields, such as rhoAlpha and rhoAlphaPhi (alphaPhi already exists), and implement the equation as following

fvm::ddt(rhoAlpha1, Yi)
 + fvm:div(rhoAlphaPhi1, Yi)
 - fvm::laplacian(turbulence->muEff(), Yi)
  ==
  0

The alternative would be to apply the product rule of differentation and use the fvc schemes.

Does anyone know of any pros or cons with these approaches, maybe there exists other approaches that are better?

Kasper