23rd International Meshing Roundtable (IMR23) Mesh Infrastructure for Coupled Multiprocess Geophysical Simulations

We have developed a sophisticated mesh infrastructure capability to support large scale multiphysics simulations such as subsurface flow and reactive contaminant transport at storage sites as well as the analysis of the e ects of a warming climate on the terrestrial arctic. These simulations involve a wide range of coupled processes including overland flow, subsurface flow, freezing and thawing of ice rich soil, accumulation, redistribution and melting of snow, biogeochemical processes involving plant matter and finally, microtopography evolution due to melting and degradation of ice wedges below the surface. In addition to supporting the usual topological and geometric queries about the mesh, the mesh infrastructure adds capabilities such as identifying columnar structures in the mesh, enabling deforming of the mesh subject to constraints and enabling the simultaneous use of meshes of di erent dimensionality for subsurface and surface processes. The generic mesh interface is capable of using three di erent open source mesh frameworks (MSTK, MOAB and STKmesh) under the hood allowing the developers to directly compare them and choose one that is best suited for the application’s needs. We demonstrate the results of some simulations using these capabilities as well as present a comparison of the performance of the di erent mesh frameworks. c 2014 The Authors. Published by Elsevier Ltd. Peer-review under responsibility of organizing committee of the 23rd International Meshing Roundtable (IMR23).

[1]  Timothy J. Tautges,et al.  An Interoperable, Data-Structure-Neutral Component for Mesh Query and Manipulation , 2010, ACM Trans. Math. Softw..

[2]  Nathan R. Tallent,et al.  HPCTOOLKIT: tools for performance analysis of optimized parallel programs , 2010, Concurr. Comput. Pract. Exp..

[3]  Mark S. Shephard,et al.  Efficient distributed mesh data structure for parallel automated adaptive analysis , 2006, Engineering with Computers.

[4]  M. Shashkov,et al.  Adjoint operators for the natural discretizations of the divergence gradient and curl on logically rectangular grids , 1997 .

[5]  Stefan Finsterle,et al.  ADVANCED SIMULATION CAPABILITY FOR ENVIRONMENTAL MANAGEMENT (ASCEM): AN OVERVIEW OF INITIAL RESULTS , 2011 .

[6]  Nicholas Nethercote,et al.  "Building Workload Characterization Tools with Valgrind" , 2006, 2006 IEEE International Symposium on Workload Characterization.

[7]  L. A. Richards Capillary conduction of liquids through porous mediums , 1931 .

[8]  Benjamin S. Kirk,et al.  Library for Parallel Adaptive Mesh Refinement / Coarsening Simulations , 2006 .

[9]  Yves Bertrand,et al.  GMDS: A Generic Mesh Data Structure , 2008 .

[10]  Tamara G. Kolda,et al.  An overview of the Trilinos project , 2005, TOMS.

[11]  John Calcote Autotools: A Practioner's Guide to GNU Autoconf, Automake, and Libtool , 2010 .

[12]  J. Canadell,et al.  Soil organic carbon pools in the northern circumpolar permafrost region , 2009 .

[13]  Ian Gorton,et al.  Amanzi and Akuna: Two New Community Codes for Subsurface Contaminant Flow and Transport , 2011 .

[14]  Rao V. Garimella,et al.  Mesh data structure selection for mesh generation and FEA applications , 2002 .

[15]  Alan B. Williams,et al.  STK-mesh tutorial minisymposium. , 2011 .

[16]  Ümit V. Çatalyürek,et al.  The Zoltan and Isorropia parallel toolkits for combinatorial scientific computing: Partitioning, ordering and coloring , 2012 .

[17]  L. A. Schoof,et al.  EXODUS II: A finite element data model , 1994 .

[18]  Rao V. Garimella,et al.  MSTK - A Flexible Infrastructure Library for Developing Mesh Based Applications , 2004, IMR.

[19]  M. Torre Jorgenson,et al.  Abrupt increase in permafrost degradation in Arctic Alaska , 2006 .