A5063635634- Cryospheric studies and observations
- Computational Fluid Dynamics and Aerodynamics
- Advanced Numerical Methods in Computational Mathematics
- Meteorological Phenomena and Simulations
- Distributed and Parallel Computing Systems
- Numerical methods for differential equations
- Landslides and related hazards
- Arctic and Antarctic ice dynamics
- Parallel Computing and Optimization Techniques
- Gas Dynamics and Kinetic Theory
- Semiconductor materials and devices
- Phase Equilibria and Thermodynamics
- Climate variability and models
- Climate change and permafrost
- Embedded Systems Design Techniques
- Atmospheric and Environmental Gas Dynamics
- Simulation Techniques and Applications
- Winter Sports Injuries and Performance
- Fluid Dynamics and Turbulent Flows
- Geophysics and Gravity Measurements
- Probabilistic and Robust Engineering Design
- Advancements in Semiconductor Devices and Circuit Design
- Matrix Theory and Algorithms
- Fluid Dynamics and Thin Films
- Scientific Research and Discoveries
Sandia National Laboratories
2015-2025
Sandia National Laboratories California
2006-2024
Oak Ridge National Laboratory
1998-2022
University of Arizona
2022
United States Department of Commerce
1998-2001
System Simulation (United Kingdom)
1999
University of Minnesota
1993-1998
Parallel Consulting (United States)
1998
Sandia Research (United States)
1998
Office of Scientific and Technical Information
1996
The Trilinos Project is an effort to facilitate the design, development, integration, and ongoing support of mathematical software libraries within object-oriented framework for solution large-scale, complex multiphysics engineering scientific problems. addresses two fundamental issues developing these problems: (i) providing a streamlined process set tools development new algorithmic implementations (ii) promoting interoperability independently developed software.Trilinos uses two-level...
This work documents the first version of U.S. Department Energy (DOE) new Exascale Earth System Model (E3SMv1). We focus on standard resolution fully coupled physical model designed to address DOE mission-relevant water cycle questions. Its components include atmosphere and land (110-km grid spacing), ocean sea ice (60 km in midlatitudes 30 at equator poles), river transport (55 km) models. base configuration will also serve as a foundation for additional configurations exploring higher...
Abstract The Energy Exascale Earth System Model Atmosphere version 1, the atmospheric component of Department Energy's is described. model began as a fork well‐known Community Model, but it has evolved in new ways, and coding, performance, resolution, physical processes (primarily cloud aerosols formulations), testing development procedures now differ significantly. Vertical resolution was increased (from 30 to 72 layers), top extended 60 km (~0.1 hPa). A simple ozone photochemistry predicts...
Abstract This work documents version two of the Department Energy's Energy Exascale Earth System Model (E3SM). E3SMv2 is a significant evolution from its predecessor E3SMv1, resulting in model that nearly twice as fast and with simulated climate improved many metrics. We describe physical lower horizontal resolution configuration consisting 110 km atmosphere, 165 land, 0.5° river routing model, an ocean sea ice mesh spacing varying between 60 mid‐latitudes 30 at equator poles. The...
SUMMARY An adaptively stabilized finite element scheme is proposed for a strongly coupled hydro‐mechanical problem in fluid‐infiltrating porous solids at strain. We first present the derivation of poromechanics model via mixture theory large deformation. By exploiting assumed deformation gradient techniques, we develop numerical procedure capable simultaneously curing multiple‐locking phenomena related to shear failure, incompressibility imposed by pore fluid, and/or incompressible solid...
Abstract. We describe and evaluate version 2.1 of the Community Ice Sheet Model (CISM). CISM is a parallel, 3-D thermomechanical model, written mainly in Fortran, that solves equations for momentum balance thickness temperature evolution ice sheets. CISM's velocity solver incorporates hierarchy Stokes flow approximations, including shallow-shelf, depth-integrated higher order, order. also includes suite test cases, links to third-party libraries, parameterizations physical processes such as...
Abstract. We introduce MPAS-Albany Land Ice (MALI) v6.0, a new variable-resolution land ice model that uses unstructured Voronoi grids on plane or sphere. MALI is built using the Model for Prediction Across Scales (MPAS) framework developing Earth system components and Albany multi-physics code base solution of coupled systems partial differential equations, which itself makes use Trilinos solver libraries. includes three-dimensional first-order momentum balance (Blatter–Pattyn) by linking...
Abstract This paper describes the first implementation of Δ x = 3.25 km version Energy Exascale Earth System Model (E3SM) global atmosphere model and its behavior in a 40‐day prescribed‐sea‐surface‐temperature simulation (January 20 through February 28, 2020). was performed as part DYnamics Atmospheric general circulation Modeled On Non‐hydrostatic Domains (DYAMOND) Phase 2 intercomparison. Effective resolution is found to be horizontal dynamics grid despite using coarser for physical...
This work documents version two of the Department Energy's Energy Exascale Earth System Model (E3SM). E3SM 2 (E3SMv2) is a significant evolution from its predecessor E3SMv1, resulting in model that nearly twice as fast and with simulated climate improved many metrics. We describe physical lower horizontal resolution configuration consisting 110 km atmosphere, 165 land, 0.5° river routing model, an ocean sea ice mesh spacing varying between 60 mid-latitudes 30 at equator poles. The...
Abstract. We describe and evaluate version 2.1 of the Community Ice Sheet Model (CISM). CISM is a parallel, 3D thermomechanical model, written mainly in Fortran 90/95, that solves equations for momentum balance thickness temperature evolution ice sheets. CISM's velocity solver incorporates hierarchy Stokes-flow approximations, including shallow-shelf, depth-integrated higher-order, higher-order. also includes suite test cases, links to third-party libraries, parameterizations physical...
Abstract. This paper describes a new parallel, scalable and robust finite element based solver for the first-order Stokes momentum balance equations ice flow. The solver, known as Albany/FELIX, is constructed using component-based approach to building application codes, in which mature, modular libraries developed part of Trilinos project are combined abstract interfaces template-based generic programming, resulting final code with access dozens algorithmic advanced analysis capabilities....
This work documents version two of the Department Energy’s Energy Exascale Earth System Model (E3SM). E3SM 2 (E3SMv2) is a significant evolution from its predecessor E3SMv1, resulting in model that nearly twice as fast and with simulated climate improved many metrics. We describe physical lower horizontal resolution configuration consisting 110 km atmosphere, 165 land, 0.5° river routing model, an ocean sea ice mesh spacing varying between 60 mid-latitudes 30 at equator poles. The...
Abstract Two approaches which employ the finite element method to solve for large‐scale, coupled, incompressible flows through adjacent porous and open domains are developed evaluated in a model spontaneous ignition of coal stockpiles. Both formulations Navier‐Stokes equations do describe flow region; two different descriptions, Darcy's law Brinkman equation, employed within region. The formulation uses employs Beavers‐Joseph slip condition novel implementation interfacial conditions. other...
We present a numerical study of the structure and stability laminar isothermal flows formed by two counterflowing jets an incompressible Newtonian fluid. demonstrate that symmetric with identical mass flow rates exhibit multiple steady states and, in certain cases, time-dependent (periodic) states. Two geometric configurations were studied based on inlet jet shapes: planar axisymmetric. Stagnation both steady-state multiplicity behaviour, while axisymmetric only multiplicity. A linearized...
An approach for incorporating embedded simulation and analysis capabilities in complex codes through template-based generic programming is presented. This relies on templating operator overloading within the C++ language to transform a given calculation into one that can compute variety of additional quantities are necessary many state-of-the-art algorithms. these ideas general graph-based assembly also These have been implemented set packages Trilinos framework demonstrated simple problem...
Albany is a multiphysics code constructed by assembling set of reusable, general components. It an implicit, unstructured grid finite element that hosts advanced features are readily combined within single analysis run. uses template-based generic programming methods to provide extensibility and flexibility; it employs residual evaluation interface support the easy addition modification physics. This coupled powerful automatic differentiation utilities used implement efficient nonlinear...
We present the set of bifurcation tracking algorithms which have been developed in LOCA software library to work with large scale application codes that use fully coupled Newton's method iterative linear solvers. Turning point (fold), pitchfork, and Hopf based on implemented, particular attention scalability problem sizes parallel computers ease implementation new codes. The is accomplished by using block elimination solve Newton iterations augmented systems. applicability such for...
This paper summarizes our efforts to develop fast algorithms for density functional theory (DFT) calculations of inhomogeneous fluids. Our goal is apply DFTs a variety problems in nanotechnology and biology. To this end we have developed DFT codes treat both atomic fluid models polymeric We three-dimensional real space Fourier algorithms. The former rely on matrix-based Newton's method while the latter couple transforms with matrix-free method. Efficient computation phase diagrams...
Abstract. We present an architecture-portable and performant implementation of the atmospheric dynamical core (High-Order Methods Modeling Environment, HOMME) Energy Exascale Earth System Model (E3SM). The original Fortran is highly scalable on conventional architectures using Message Passing Interface (MPI) Open MultiProcessor (OpenMP) programming models. rewrite model in C++ use Kokkos library to express on-node parallelism a largely architecture-independent implementation. provides...
Abstract This paper presents an approach for determining the linear stability of steady states partial differential equations (PDEs) on massively parallel computers. Linearizing transient behavior around a state solution leads to eigenvalue problem. The eigenvalues with largest real part are calculated using Arnoldi's iteration driven by novel implementation Cayley transformation. transformation requires system at each Arnoldi iteration. is done iteratively so that algorithm scales problem...
A template-based generic programming approach was presented in Part I of this series papers [Sci. Program. 20 (2012), 197–219] that separates the development effort a physical model from computing additional quantities, such as derivatives, needed for embedded analysis algorithms. In paper, we describe implementation details using simulation and partial differential equations (PDEs). We detail several hurdles have encountered, some software infrastructure developed to overcome them. end with...