CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

Level Set Numerical Surface Tension Calculation

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 5, 2011, 00:33
Default Level Set Numerical Surface Tension Calculation
  #1
New Member
 
Michael Franklin
Join Date: Nov 2011
Posts: 4
Rep Power: 15
mbfranklin@gmail.com is on a distinguished road
Given a phase-function 0<phi(x,y)<1 that varies smoothly from 0 to 1 in a narrow region (approx. 6 grid cells from 0.05<phi<0.95), the gradients are:

dphi/dx
dphi/dy

which are approximated via central difference. This value is located at a grid cell center (i,j) due to central difference formula. The normal vector field is then approximated as

N = gradient/mag(gradient)

This is also located at cell center (i,j). The curvature is the divergence of normal vectors, i.e.

kappa = div(normal)

also located at cell center (i,j). Now...for the surface tension, whose x and y components are located at right edges (i+1/2,j) and top edges (i,j+1/2), respectively...

When I plot gradients, normals, and curvatures, they are all symmetric, i.e., around a circle, all values point perfectly radially outward (normal) to the level set phi=0.5 interface, but when I plot my surface tension, it appears to be asymmetric! I have even taken care to stagger these two vector fields (that are located at different spatial points, remember!) when I plot the field.

Does anybody have any advice? The surface tension (presented in the paper who's work I'm following, Olsson & Kreiss, "A conservative level set method for two phase flow") is a forward difference calculation, whereas normal and gradient calculations are central differences. I'm following their method EXACTLY, I've checked everything...and yet, asymmetric surface tension...help!
mbfranklin@gmail.com is offline   Reply With Quote

Reply

Tags
gradients, level set method, normals, phase field, surface tension

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 of level set method as initial condition sooroo3 FLUENT 4 July 31, 2018 05:20
[snappyHexMesh] Snap the mesh from a level set function instead of a STL surface ginette OpenFOAM Meshing & Mesh Conversion 2 April 5, 2011 11:21
OF 1.6 | Ubuntu 9.10 (64bit) | GLIBCXX_3.4.11 not found piprus OpenFOAM Installation 22 February 25, 2010 14:43
Env variable not set gruber2 OpenFOAM Installation 5 December 30, 2005 05:27
How to set environment variables kanishka OpenFOAM Installation 1 September 4, 2005 11:15


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