I am trying to solve a case involving two phase boundaries using SteadyPhaseChange. My axisymmetric setup consists of three blocks stacked on top of each other with top and bottom solid and liquid in the center. When using only one SteadyPhaseChangesolver it works well, with the case attached to this post I ran into the following error:
Code: Select all
SteadyPhaseChange: --------------------------------------------
SteadyPhaseChange: Using steady algorithm to find the isotherm
SteadyPhaseChange: --------------------------------------------
SteadyPhaseChange: Steady state isoterm formulation
SteadyPhaseChange: Melting point found: 5.0500E+02
Note: The following floating-point exceptions are signalling: IEEE_UNDERFLOW_FLAG
STOP 1
ERROR:: SteadyPhaseChange: Isotherm is empty thus cannot map phase change surface
I played a bit around and tried different configurations but nothing is working. It can be all found here, including a more extensive description and a couple of screenshots: https://github.com/nemocrys/elmer-double-phasechange
Thanks for your help and best regards
Arved