U.S. Department of Energy Office of Biological and Environmental Research

BER Research Highlights

Performance Analysis of Fully Explicit and Fully Implicit Solvers Within A Spectral Element Shallow-Water Atmosphere Model
Published: November 14, 2017
Posted: January 23, 2018

Performance analysis of fully implicit Spectral Element Community Atmosphere Model (CAM-SE).

The Science
A Trilinos library based implicit solver that uses the Graphics Processing Unit (GPU)  within the residual has been implemented within the spectral element CAM. The implicit solver provides accurate solutions for a range of problem types and scales very well on high performance computers, to more than 86,400 cores. The Newton-based implicit algorithms have been evaluated at this scale of complexity for atmosphere fluid flow applications for a wide variety of configurations.

The Impact
The implicit solver is able to use large time step sizes such that no subcycling of tracers and physics is needed. Also, it shows equal performance to explicit time-stepping for strongly regionally refined configurations. There is potential for increased efficiency using the GPU, however more development of the interface between the solver library and the code is needed.

Several methods utilizing a Newton-Krylov nonlinear solver are evaluated for a range of configurations of the shallow-water dynamical core of the spectral element community atmosphere model to evaluate their computational performance. These configurations are designed to explore the attributes of each method under different but relevant model usage scenarios, including varied spectral order within an element, static regional refinement, and scaling to the largest problem sizes. The limitations and benefits to using explicit Runge-Kutta versus implicit multistep methods, with different parameters and settings, are discussed in light of the trade-offs with Message Passing Interface (MPI) communication and memory and their inherent efficiency bottlenecks. The recommendation for future work using the implicit solvers is conditional based on scale separation and the stiffness of the problem. For the regionally refined configurations, the implicit method has about the same efficiency as the explicit method, without considering efficiency gains from a preconditioner. Initial simulations with OpenACC directives to utilize a GPU when performing function evaluations show improvements locally, and that overall gains are possible with adjustments to data exchanges.

Contacts (BER PMs)
Dorothy Koch
Earth System Modeling

PI Contact
Katherine Evans
Oak Ridge National Laboratory

Support from the U.S. Department of Energy Office of Science, Biological and Environmental Research, Earth System Modeling Program and the Advanced Scientific Computing Research Office, under the SciDAC program.

Evans, K. J., R. K. Archibald, D. J. Gardner, M. R. Norman, M. A. Taylor, C. S. Woodward, and P. H. Worley. "Performance Analysis of Fully Explicit and Fully Implicit Solvers within A Spectral Element Shallow-Water Atmosphere Model." The International Journal of High Performance Computing Applications (2017). [DOI: 10.1177/1094342017736373].

Related Links

Topic Areas:

  • Research Area: Earth and Environmental Systems Modeling
  • Cross-Cutting: Scientific Computing and SciDAC

Division: SC-23.1 Climate and Environmental Sciences Division, BER


BER supports basic research and scientific user facilities to advance DOE missions in energy and environment. More about BER

Recent Highlights

Aug 24, 2019
New Approach for Studying How Microbes Influence Their Environment
A diverse group of scientists suggests a common framework and targeting of known microbial processes [more...]

Aug 08, 2019
Nutrient-Hungry Peatland Microbes Reduce Carbon Loss Under Warmer Conditions
Enzyme production in peatlands reduces carbon lost to respiration under future high temperatures. [more...]

Aug 05, 2019
Amazon Forest Response to CO2 Fertilization Dependent on Plant Phosphorus Acquisition
AmazonFACE Model Intercomparison. The Science Plant growth is dependent on the availabi [more...]

Jul 29, 2019
A Slippery Slope: Soil Carbon Destabilization
Carbon gain or loss depends on the balance between competing biological, chemical, and physical reac [more...]

Jul 15, 2019
Field Evaluation of Gas Analyzers for Measuring Ecosystem Fluxes
How gas analyzer type and correction method impact measured fluxes. The Science A side- [more...]

List all highlights (possible long download time)