Interface modeling in incompressible media using level sets in Escript

Details Interface modeling in incompressible media using level sets in Escript Interface modeling in incompressible media using level sets in Escript

Aug 2007

adsabs.harvard.edu/cgi-bin/nph-data_query?bibcode=2007pepi..163...23g&link_type=abstract

Physics of the Earth and Planetary Interiors, Volume 163, Issue 1-4, p. 23-34.

Physics

6

Scientific paper

We use a finite element (FEM) formulation of the level set method to model geological fluid flow problems involving interface propagation. Interface problems are ubiquitous in geophysics. Here we focus on a Rayleigh-Taylor instability, namely mantel plumes evolution, and the growth of lava domes. Both problems require the accurate description of the propagation of an interface between heavy and light materials (plume) or between high viscous lava and low viscous air (lava dome), respectively. The implementation of the models is based on Escript which is a Python module for the solution of partial differential equations (PDEs) using spatial discretization techniques such as FEM. It is designed to describe numerical models in the language of PDEs while using computational components implemented in C and C++ to achieve high performance for time-intensive, numerical calculations. A critical step in the solution geological flow problems is the solution of the velocity pressure problem. We describe how the Escript module can be used for a high-level implementation of an efficient variant of the well-known Uzawa scheme. We begin with a brief outline of the Escript modules and then present illustrations of its usage for the numerical solutions of the problems mentioned above.

Affiliated with

Also associated with

No associations

Rating

Interface modeling in incompressible media using level sets in Escript does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.

If you have personal experience with Interface modeling in incompressible media using level sets in Escript, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Interface modeling in incompressible media using level sets in Escript will most certainly appreciate the feedback.

Rate now

Comments { 0 }

Profile ID: LFWR-SCP-O-1797774