[Yassin-K]

Hey Foamers,
I am simulating ice formation on wind turbine airfoils so I am injecting water droplets at temperatures < 273.16 deg K in reactingParcelFoam (OF6).

When I do that....I received the following error:

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  ? in "/lib/x86_64-linux-gnu/libm.so.6"
#4  pow in "/lib/x86_64-linux-gnu/libm.so.6"
#5  Foam::H2O::rho(double, double) const at ??:?
#6  Foam::regionModels::surfaceFilmModels::liquidFilmThermo::rho() const at ??:?
#7  Foam::regionModels::surfaceFilmModels::thermoSingleLayer::correctThermoFields() at ??:?
#8  Foam::regionModels::surfaceFilmModels::kinematicSingleLayer::preEvolveRegion() at ??:?
#9  Foam::regionModels::surfaceFilmModels::thermoSingleLayer::preEvolveRegion() at ??:?
#10  Foam::regionModels::regionModel::evolve() at ??:?
#11  ? at ??:?
#12  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#13  ? at ??:?
Floating point exception (core dumped)

does anybody know the problem??

Thank you in advance.

Regards,
Khaled

[Yassin-K]

Hi all,
I have found a work-around that can temporarily solve the problem but not completely.
In surfaceFilmProperties, I put filmThermomodel = constant and I put constantCoeffs with the following values (values correspond to water at T = 273.16 deg K):

    constantCoeffs
    {
        specie      H2O;//all properties at 273.16 deg K
        rho0        999.8;
        mu0         1.792e-3;
	Cp0             4217;
        sigma0      0.07;
	kappa0	0.561;
	W0	18.0153;
	Tb0	100;
	pv0	611.3;
	hl0	2501e3;
	D0	2.09e-5;
    }

and it worked.
if somebody has a better idea, please share it.