prisms-center/phaseField

Missing cubic term in model used by alloySolidification app

Closed this issue · 1 comments

The term phi^3 in governing equation 132 from Echebarria et al. (https://doi.org/10.1103/PhysRevE.70.061604). Is missing from the model. Adding this term results in unphysical simulation results with the current settings of the app.

The term is not missing. It appears by simply rearranging the expression below in the code:

scalarvalueType eq_xi =
((phi-constV(lamda)(U+tep+constV(U_off))(constV(1.0)-phiphi))(constV(1.0)-phi*phi));