1D case error? : “MoleFractions_From_MassFractions” fails when mass fractions are all zero.

[hofbeer]

Dear developers and contributors,

I try to simulate a 1D premixed flame to get the laminar flame speed. So I made 1D shape as below.

Boundary Conditions Inlet: fixedValue for mass fractions Outlet: InletOutlet for mass fractions Sides: empty (to achieve 1D simulation)

Here is a problem. When I run this case, this error appears.

I checked ~laminarSMOKE/solvers/laminarSMOKE/properties.H, and it seems to be caused by some boundary faces that have zero mass fractions initially. “MoleFractions_From_MassFractions” does not work if all input mass fractions are zero. https://github.com/acuoci/laminarSMOKE/blob/65b29ff7e957ba7522e454e5f766ccecfee98c5a/solvers/laminarSMOKE/properties.H#L232-L235

My trial solutions)

1. Run reactingFoam by one timestep and restart with a new initial solution, so that zero mass fractions does not appear at boundary faces.
2. Or edit properties.H; (ex, If all mass fractions are zero, omit above lines.)

Two ways seem to work.
If there are better ways, please let me know. It would be grateful if properties.H could work for this kind of situation.

Thank you

Best regards, Daeyoung Jun