CFD Online Discussion Forums - OF10| reactingFoam with multi-component transport

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)

- - OF10| reactingFoam with multi-component transport (https://www.cfd-online.com/Forums/openfoam-solving/245151-of10-reactingfoam-multi-component-transport.html)

OF10| reactingFoam with multi-component transport

Hi all,

I am trying to model a hydrogen diffusion flame using reactingFoam. The reference experimental flame is the DLR-H5 flame (Re_j=6600). There is noticeable preferential diffusion happening close to the jet-nozzle and I am hoping to capture this via an LES/EDC (finite-rate chem w/ TCI) simulation using OF10 solver reactingFoam.

I am fairly new to OpenFOAM and C++. I have gathered that newer versions of reactingFoam (>=OF8) has the capability to include preferential diffusion effects and thermal diffusion. I referred to the user guide section 7.1 for information on setting up my case.

According to the OF10 userguide and the Sandia_D_LTS example, I have set-up my constant/physicalProperties as:

Code:

thermoType

{

    type            hePsiThermo;

    mixture         multiComponentMixture;

    transport       sutherland;

    thermo          janaf;

    energy          sensibleEnthalpy;

    equationOfState perfectGas;

    specie          specie;

}

under the laminar examples I found a further mixture type

Code:

coefficientWilkeMultiComponentMixture

and also came across the transport type

Code:

MaxwellStefan

under /src/ThermophysicalTransportModels/laminar

I am interested in modelling the momentum transport with LES while combustion is modelled via TCI model EDC with detailed multi-component transport. (Wilke for mixture viscosity calculations and including Maxwell-Stefan diffusivities in the species transport equations). Neither of these class references are listed in the userguide under section 7.1.

My questions is if what I am trying to achieve is possible through OF10 and if so, where can I find some more information of the inputs for the coefficientWilkeMultiComponentMixture and MaxwellStefan.

Thank you in advance,

Mana

Yes, this should be possible in principle. Typically convergence is hard to achieve.

Below is a recent recent (among many other others).

With more precise indications of what you need, I can see how to provide further input.

@article{chen-19,
title = {Large Eddy Simulation of a dual swirl gas turbine combustor: Flame/flow structures and stabilisation under thermoacoustically stable and unstable conditions},
journal = {Combustion and Flame},
volume = {203},
pages = {279-300},
year = {2019},
issn = {0010-2180},
doi = {https://doi.org/10.1016/j.combustflame.2019.02.013},
url = {https://www.sciencedirect.com/science/article/pii/S001021801930077X},
author = {Zhi X. Chen and Ivan Langella and Nedunchezhian Swaminathan and Michael Stöhr and Wolfgang Meier and Hemanth Kolla},
keywords = {Large Eddy Simulation, Partially premixed combustion, Gas turbine model combustor, Dual swirl, Self-excited thermoacoustic instability}
}