[sramac3]

Hello,

I am trying to do some stefan sucking interface validation using icoReactingMultiphaseInterFoam on v1912. I am trying to set up the initial temperature field at time=0.1. But the setExprFieldsDict doesn't allow me to use the erfc() or erf() functions. I can still use sqrt(), exp(), .... and other functions in the expression, but not the error function.

I get the following error,

"Object erfc does not exist or wrong type in expression at position"

Tsup + (Tsat-Tsup)*erfc(N)/erfc(l)
^^^^ near here

I also tried using a taylor series expansion for the error function (considered five terms), but the temperature set is wrong.

I am herewith attaching the setExprFieldsDict for your reference.

Any sort of help would be appreciated.



SMR

[HPE]

I haven't tried it, but: "std::erfc(x)" would help?

[jakvb]

Were you able to solve this problem in the end? I'm also trying to use erf() in setExprFieldsDict