|January 17, 2012, 11:52||
phase field interpolation
Join Date: Jun 2009
Posts: 30Rep Power: 6
I'm working in the development of my mesoscopic solidification model.
In few words I solve three coupled equations:
1 for the phase (phi) :
where phi=0 in the liquid
phi=1 in the solid
and varies as an hyperbolic arctg in between (the right hand side term imposes a fixed depth surface thickness)
Results for a growing circle are shown below (just to have an idea of what I'm trying to do):
Uploaded with ImageShack.us
2 i solve a diffusion equation for the solute :
3. i couple with an entalpy equation:
My problem is to impose the v_front velocity cause it is function of a parameter
in the interface phi ]0,1[ domain.
WHAT I NEED:
is to know the position of the isocontours 0.9 and 0.5 in the domain and to evaluate so omega.
I can try to do it by bilinear interpolation by hands but I suppose (?) that perhaps there are some "tools" already developed in OpenFOAM (and that I don't find..).
Can someone please help me or give me some ideas in order to do it in a proper way?
|Thread||Thread Starter||Forum||Replies||Last Post|
|Setting an external B0_x field in MHD module||fferroni||FLUENT||4||February 25, 2013 14:15|
|Surface interpolation schemes and parallelization||jutta||OpenFOAM Running, Solving & CFD||0||February 25, 2010 14:32|
|phase field or Cahn-Hilliard method||dusky.he||Main CFD Forum||3||February 2, 2006 16:43|
|compressible two phase flow in CFX4.4||youngan||CFX||0||July 1, 2003 23:32|
|Help!||pengsh||CFX||2||July 12, 2001 02:43|