|November 5, 2011, 00:33||
Level Set Numerical Surface Tension Calculation
Join Date: Nov 2011
Posts: 3Rep Power: 5
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:
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!
|gradients, level set method, normals, phase field, surface tension|
|Thread||Thread Starter||Forum||Replies||Last Post|
|setting of level set method as initial condition||sooroo3||FLUENT||3||January 7, 2015 20:15|
|Snap the mesh from a level set function instead of a STL surface||ginette||OpenFOAM Native Meshers: snappyHexMesh and Others||2||April 5, 2011 10: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 10:15|