[mike.franky]

Hi,

I am trying to run a multiphase, multicomponent case using the reactingMultiphaseEulerFoam. Each phase is composed of various species. Therefore, in the boundary conditions one has to define the volume fraction of the phase and the mass fraction of the species of each phase for each cell.

In all the tutorials, the mass fraction of the internal field is uniformly set for all the cells, even if some cells have a zero volume fraction (alpha) for the phase itself. Furthermore, if I use setFieldsDict to only set the mass fraction of cells with non-zero values of alpha, the simulation crashes. So what exactly is Y? Is it simply the mass fraction of the species in the phase volume contained in the cell (i.e. zero if the volume fraction is also set to zero)? If so, why does the simulation crash if I don't define the mass fractions in all the cells?

Thanks in advance for any help

Mike