Computational Science Stack Exchange is a question and answer site for scientists using computers to solve scientific problems. It's 100% free, no registration required.

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

I have a 2D surface in 3D that I want to advect under a velocity field. More precisely, I have a surface $S$ and a velocity field $v$ and I want to advect $S$ under $v$ using the flow map of $v$, i.e. I want to compute $F^t(S_0)$, where $F^t$ is the flow of $v$. The velocity field is chaotic, so a direct computation leaves a bunch of scattered points. $S_0$ is very smooth.

I have been told that there are the so-called level set methods to do what I'd like. The idea sounds very nice, but this is not my background and it would be great if there are any codes available that do this.

share|improve this question

Level set methods are really only one particular case of advection PDEs which you can solve with any of the modern PDE toolboxes such as deal.II (disclaimer: that's my own library), fenics, libmesh, ... What makes it a level set method is that you the evaluate the solution function to find that surface where, for example, the solution equals zero. I'd try to look into one of these PDE toolboxes and see what you can use to solve your problem. They certainly all can solve advection equations.

share|improve this answer

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.