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.

laminarSMOKE/solvers/laminarSMOKE/properties.H

Lines 232 to 235 in 65b29ff

	for(int i=0;i<thermodynamicsMapXML->NumberOfSpecies();i++)
	massFractions[i+1] = Y[i].boundaryField()[patchi][facei];
	thermodynamicsMapXML->MoleFractions_From_MassFractions(moleFractions.GetHandle(),pMWmix[facei],massFractions.GetHandle());

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