Discussion Forum

Look "flc1hs" or "flc2hs".

Thank you very much for the reply, I tried this equation,

0.25*(flc1hs(x^2,-0.5)+flc1hs(y^2,-0.5)) and I am getting this as the initial condition (initial_comsol.jpg).

I want to achieve something like this, (initial_malab.jpg)

What happens if you write

0.4*flc1hs(0.04 - x^2 - y^2, 0.01)?

yep that worked. I am getting this white border around the initial cell. Do you know why?

I do not understand the image. How dense is the mesh around the circle? Please send a height plot.

attached are the screenshots of initial stage and after 0.6seconds

White areas probably represent function values < 0 and also shows your mesh. If you have a time dependent problem the initial steps tend to give negative values as the step function is numerically rather tricky thing to realize. Have you tried flc2hs which has also continuous second derivative? But in principle it appears to work. Defining the second parameter as 0.001 makes the transition range very sharp; values below that I have not found to have any improvement.

Thank you very much for the quick reply

I tried flc2hs, and I'm getting a similar result. I changed the second parameter to 0.001 and it is the same,

It is like this only for the initial step. at 0.1s onward white boundary is gone

You can also trying using logicals.

for instance (x^2 + y^2 > 0.4) will return 1 when true and 0 when false.

Hello Ashan,
I get the feeling that this may be a case of the situation that's been discussed in this Discussion Forum before (I cannot find the thread right now) where the user tries to impose a very rapid variation of a quantity at the nodes of a single quadratic element, not realizing that this can cause the quantity to overshoot or undershoot in between the nodes.
E.g. if you imagine a 1D quadratic element covering [-2;2] and you impose at the 3 nodes u(-2)=1, u(0)=0 and u(2)=0, then u will take negative values in [0;2]. This can be counter-intuitive but it's true. If your plot is set up so that only values of u between 0 and 1 are displayed, then you will naturally get a "hole" in your plot, like what your screenshots above show.
The solution is to use a finer mesh, one that uses more than one element to capture the gradient that you're imposing.
Jeff