Problem of liquid fraction in example galliumMelting2d

[qiong]

Dear OpenLB community,

I am trying to add obstacle in the example galliumMelting2d (shown in link1). To apply the constructed obstacle for the simulation add
std::vector<T> extend1(2,T());
extend1[0] = 0.01;
extend1[1] = 0.05;
std::vector<T> origin1(2,T());
origin1[0] = 0.05;
origin1[1] = 0.02;
IndicatorCuboid2D<T> cuboid21(extend1, origin1);
superGeometry.rename(1,5,cuboid21);
From the result, I can see that the enthalpy distribution is correct (shown in link2). But it seems that the liquid fraction of the obstacle is liquid from the begining instead of solid phase(shown in link3). The enthalpy distribution is not consistent with the liquid fraction. It is because I did not set the geometry/material property correctly or the liquid fraction function is not running correctly？ Can anyone have any idea to figure it out?
link1:
https://photos.google.com/album/AF1QipPY_P4hn8sVUkZatYbr-Of5gaik8RVSqw6TdYY
link2: https://photos.google.com/album/AF1QipPY_P4hn8sVUkZatYbrOf5gaik8RVSqw6TdYY/photo/AF1QipOK8HNtGY8ILef-mnTxs218vX-DKpRbZygZ5Yk
link3:
https://photos.google.com/album/AF1QipPY_P4hn8sVUkZatYbr-Of5gaik8RVSqw6TdYY/photo/AF1QipN_aDsPXf_g4vElt91TD343nOfknfO6iaBGXAI

Thanks for your help.

Best regards,
Qiong

[Author]

Posts