|
[Sponsors] |
Level Set Numerical Surface Tension Calculation |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
November 5, 2011, 00:33 |
Level Set Numerical Surface Tension Calculation
|
#1 |
New Member
Michael Franklin
Join Date: Nov 2011
Posts: 4
Rep Power: 15 |
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! |
|
Tags |
gradients, level set method, normals, phase field, surface tension |
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 |