CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

phase field interpolation

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 17, 2012, 12:52
Default phase field interpolation
  #1
Member
 
valerio
Join Date: Jun 2009
Location: Nancy
Posts: 30
Rep Power: 16
Alucard is on a distinguished road
Goodmoring
I'm working in the development of my mesoscopic solidification model.
In few words I solve three coupled equations:
1 for the phase (phi) :
d(phi)/dt+v_front*mag(grad(phi))=[lapl(phi)+phi*(1-phi)*(1-2*phi)-mag(grad(phi))*div[grad(phi)/mag(grad(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 :
dCl/dt=lapl(Cl)+S

3. i couple with an entalpy equation:

dH/dt=lapl(T)
.
My problem is to impose the v_front velocity cause it is function of a parameter
omega=Cl(phi=0.9)-Cl(phi=0.5)
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?
Thank you
Alucard is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Setting an external B0_x field in MHD module fferroni FLUENT 11 August 1, 2022 09:20
Surface interpolation schemes and parallelization jutta OpenFOAM Running, Solving & CFD 0 February 25, 2010 15:32
phase field or Cahn-Hilliard method dusky.he Main CFD Forum 3 February 2, 2006 17:43
compressible two phase flow in CFX4.4 youngan CFX 0 July 2, 2003 00:32
Help! pengsh CFX 2 July 12, 2001 03:43


All times are GMT -4. The time now is 21:05.