dendritic-growth simulation has succeed using prisms-pf about phasefield HOT 5 CLOSED

Hi Xin-Bo,
Great! I'll start looking through your changes to figure out how to generalize the existing framework to include this functionality.

Is there a particular paper you got your solidification model from? Having a little more familiarity with the model will help me parse the code. Thanks!

from phasefield.

Hi Steve,
The code is based on the model by Wheeler et al.
See: here
The anisotropy part refers to Karma's model, in order to handle the 3D case. And it can be deduced to the style in Wheeler's model in 2D case.

Xin-Bo

from phasefield.

Hi Xin-Bo,
After some thought, I think there are two ways to access the time derivative without a major change to the framework.

First of all, just to make sure that I'm clear on the issue, you need the old solution to calculate the time derivative of the order parameter for the heat equation, right? You should be able to calculate that from your residual for the order parameter -- to update the order parameter you use it's time derivative.

If you do absolutely need access to the old value of the order parameter, you could create it as a third variable (so you'd have temp, eta, and eta_old as variables). Then you just set eta_old = eta in residualsRHS and it'll have the value of eta from the previous time step.

Let me know if this makes sense and you think either way will work.

Best,
Steve

from phasefield.

Thanks, Steve, I think both ways are feasible.
In my pinion, the first way is not quite elegant, so I prefer the second one.
However, I don't know how to pass eta_old into residualRHS without modifying the framework. Just as the way in my code, I changed the parameter "std::vector<modelVariable>" to a Vector style "std::vector<std::vector<modelVariable>>" in order to accept eta and eta_old.

from phasefield.

In v2.0 a dendritic solidification app was added, very similar to the one here.

from phasefield.