Welcome guest
Adaptive finite element methods in geodynamics: Convection dominated mid-ocean ridge and subduction zone simulations
D.R. Davies, J.H. Davies, O. Hassan, K. Morgan, P. Nithiarasu
International Journal of Numerical Methods for Heat & Fluid Flow
2008
1015 - 1035
0961-5539
10.1108/09615530810899079
Emerald Group Publishing Limited
This project was enabled with the assistance of IBM Deep Computing at Swansea University, together with the Helix facility at Cardiff University (HEFCW, SRIF). DRD would like to acknowledge support from both NERC and EPSRC as part of the Environmental Mathematics and Statistics (EMS) studentship programme (NER/S/E/2004/12725). The authors also thank Scott King for help and support with CONMAN and a number of colleagues, particularly Ben Evans, Richard Davies and Martin Wolstencroft, for helpful discussions.
Purpose – The purpose of this paper is to present an adaptive finite element procedure that improves the quality of convection dominated mid-ocean ridge (MOR) and subduction zone (SZ) simulations in geodynamics.
Design/methodology/approach – The method adapts the mesh automatically around regions of high-solution gradient, yielding enhanced resolution of the associated flow features. The approach utilizes an automatic, unstructured mesh generator and a finite element flow solver. Mesh adaptation is accomplished through mesh regeneration, employing information provided by an interpolation-based local error indicator, obtained from the computed solution on an existing mesh.
Findings – The proposed methodology works remarkably well at improving solution accuracy for both MOR and SZ simulations. Furthermore, the method is computationally highly efficient.
Originality/value – To date, successful goal-orientated/error-guided grid adaptation techniques have, to the knowledge, not been utilized within the field of geodynamics. This paper presents the first true geodynamical application of such methods.
Finite element analysis, Flow, Meshes, Oceanography, Simulation
Research paper