Interface modeling in incompressible media using level sets in Escript

Gross, Lutz, Bourgouin, Laurent, Hale, Alina and Muhlhaus, Hans (2007) Interface modeling in incompressible media using level sets in Escript. Physics of The Earth and Planetary Interiors, 163 1-4: 23-24.

Attached Files (Some files may be inaccessible until you login with your UQ eSpace credentials)
Name Description MIMEType Size Downloads
gross.pdf gross.pdf application/pdf 293.21KB 299

Author Gross, Lutz
Bourgouin, Laurent
Hale, Alina
Muhlhaus, Hans
Title Interface modeling in incompressible media using level sets in Escript
Journal name Physics of The Earth and Planetary Interiors  (ERA 2012 Listed)    (ERA 2010 Rank A)
Publication date 2007-08-15
Sub-type Article
Volume number 163
Issue number 1-4
ISSN 0031-9201; 1872-7395
Start page 23
End page 24
Total pages 2
Place of publication Amsterdam, Netherlands
Publisher Elsevier
Language eng
Subject 260299 Geophysics not elsewhere classified
Formatted abstract 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.
Keyword Partial differential equations
Numerical software
Level set method
Uzawa scheme
Rayleigh-Taylor instability
Mantle plumes
Lava dome modeling
Q-Index Code C1
Q-Index Status Provisional Code
Institutional Status Unknown
Additional Notes Issue theme: "Computational Challenges in the Earth Sciences"
 
Versions
Version Filter Type
Citation counts: Thomson Reuters Citation Count Cited 15 times in Thomson Reuters Researcher ID
Scopus Citation Count Cited 13 times in Scopus
Access Statistics: 191 Abstract Views, 299 File Downloads  -  Detailed Statistics
Created: Wed, 10 Jan 2007, 14:06:18 EST by Lutz Gross on behalf of Earth Systems Science Computational Centre  -  Detailed History