DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Large-scale optimization-based non-negative computational framework for diffusion equations: Parallel implementation and performance studies

Abstract

It is well-known that the standard Galerkin formulation, which is often the formulation of choice under the finite element method for solving self-adjoint diffusion equations, does not meet maximum principles and the non-negative constraint for anisotropic diffusion equations. Recently, optimization-based methodologies that satisfy maximum principles and the non-negative constraint for steady-state and transient diffusion-type equations have been proposed. To date, these methodologies have been tested only on small-scale academic problems. The purpose of this paper is to systematically study the performance of the non-negative methodology in the context of high performance computing (HPC). PETSc and TAO libraries are, respectively, used for the parallel environment and optimization solvers. For large-scale problems, it is important for computational scientists to understand the computational performance of current algorithms available in these scientific libraries. The numerical experiments are conducted on the state-of-the-art HPC systems, and a single-core performance model is used to better characterize the efficiency of the solvers. Furthermore, our studies indicate that the proposed non-negative computational framework for diffusion-type equations exhibits excellent strong scaling for real-world large-scale problems.

Authors:
 [1];  [2];  [1]
  1. Univ. of Houston, Houston, TX (United States)
  2. Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Publication Date:
Research Org.:
Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)
Sponsoring Org.:
USDOE
OSTI Identifier:
1296665
Report Number(s):
LA-UR-15-24900
Journal ID: ISSN 0885-7474
Grant/Contract Number:  
AC52-06NA25396
Resource Type:
Accepted Manuscript
Journal Name:
Journal of Scientific Computing
Additional Journal Information:
Journal Name: Journal of Scientific Computing; Journal ID: ISSN 0885-7474
Publisher:
Springer
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING; Computer Science; Environmental Protection

Citation Formats

Chang, Justin, Karra, Satish, and Nakshatrala, Kalyana B. Large-scale optimization-based non-negative computational framework for diffusion equations: Parallel implementation and performance studies. United States: N. p., 2016. Web. doi:10.1007/s10915-016-0250-5.
Chang, Justin, Karra, Satish, & Nakshatrala, Kalyana B. Large-scale optimization-based non-negative computational framework for diffusion equations: Parallel implementation and performance studies. United States. https://doi.org/10.1007/s10915-016-0250-5
Chang, Justin, Karra, Satish, and Nakshatrala, Kalyana B. Tue . "Large-scale optimization-based non-negative computational framework for diffusion equations: Parallel implementation and performance studies". United States. https://doi.org/10.1007/s10915-016-0250-5. https://www.osti.gov/servlets/purl/1296665.
@article{osti_1296665,
title = {Large-scale optimization-based non-negative computational framework for diffusion equations: Parallel implementation and performance studies},
author = {Chang, Justin and Karra, Satish and Nakshatrala, Kalyana B.},
abstractNote = {It is well-known that the standard Galerkin formulation, which is often the formulation of choice under the finite element method for solving self-adjoint diffusion equations, does not meet maximum principles and the non-negative constraint for anisotropic diffusion equations. Recently, optimization-based methodologies that satisfy maximum principles and the non-negative constraint for steady-state and transient diffusion-type equations have been proposed. To date, these methodologies have been tested only on small-scale academic problems. The purpose of this paper is to systematically study the performance of the non-negative methodology in the context of high performance computing (HPC). PETSc and TAO libraries are, respectively, used for the parallel environment and optimization solvers. For large-scale problems, it is important for computational scientists to understand the computational performance of current algorithms available in these scientific libraries. The numerical experiments are conducted on the state-of-the-art HPC systems, and a single-core performance model is used to better characterize the efficiency of the solvers. Furthermore, our studies indicate that the proposed non-negative computational framework for diffusion-type equations exhibits excellent strong scaling for real-world large-scale problems.},
doi = {10.1007/s10915-016-0250-5},
journal = {Journal of Scientific Computing},
number = ,
volume = ,
place = {United States},
year = {Tue Jul 26 00:00:00 EDT 2016},
month = {Tue Jul 26 00:00:00 EDT 2016}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 11 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

Flexible complementarity solvers for large-scale applications
journal, February 2006


Convex Optimization
book, January 2004


Maximum principle and uniform convergence for the finite element method
journal, February 1973


Contaminant remediation decision analysis using information gap theory
journal, February 2012

  • Harp, Dylan R.; Vesselinov, Velimir V.
  • Stochastic Environmental Research and Risk Assessment, Vol. 27, Issue 1
  • DOI: 10.1007/s00477-012-0573-1

Isotopic evidence for reduction of anthropogenic hexavalent chromium in Los Alamos National Laboratory groundwater
journal, May 2014


Three-phase numerical model for subsurface hydrology in permafrost-affected regions (PFLOTRAN-ICE v1.0)
journal, January 2014


A Fast and High Quality Multilevel Scheme for Partitioning Irregular Graphs
journal, January 1998


A simulator for modeling coupled thermo-hydro-mechanical processes in subsurface geological media
journal, September 2014


Efficient representation of computational meshes
journal, January 2009


Automated Solution of Differential Equations by the Finite Element Method
book, January 2012

  • Logg, Anders; Mardal, Kent-Andre; Wells, Garth
  • Lecture Notes in Computational Science and Engineering
  • DOI: 10.1007/978-3-642-23099-8

On mesh restrictions to satisfy comparison principles, maximum principles, and the non-negative constraint: Recent developments and new results
journal, September 2016

  • Mudunuru, Maruti; Nakshatrala, Kalyana Babu
  • Mechanics of Advanced Materials and Structures, Vol. 24, Issue 7
  • DOI: 10.1080/15502287.2016.1166160

Enforcing the non-negativity constraint and maximum principles for diffusion with decay on general computational grids
journal, August 2010

  • Nagarajan, H.; Nakshatrala, K. B.
  • International Journal for Numerical Methods in Fluids, Vol. 67, Issue 7
  • DOI: 10.1002/fld.2389

A numerical framework for diffusion-controlled bimolecular-reactive systems to enforce maximum principles and the non-negative constraint
journal, November 2013

  • Nakshatrala, K. B.; Mudunuru, M. K.; Valocchi, A. J.
  • Journal of Computational Physics, Vol. 253
  • DOI: 10.1016/j.jcp.2013.07.010

A Numerical Methodology for Enforcing Maximum Principles and the Non-Negative Constraint for Transient Diffusion Equations
journal, January 2016

  • Nakshatrala, K. B.; Nagarajan, H.; Shabouei, M.
  • Communications in Computational Physics, Vol. 19, Issue 1
  • DOI: 10.4208/cicp.180615.280815a

Non-negative mixed finite element formulations for a tensorial diffusion equation
journal, October 2009


Roofline: an insightful visual performance model for multicore architectures
journal, April 2009

  • Williams, Samuel; Waterman, Andrew; Patterson, David
  • Communications of the ACM, Vol. 52, Issue 4
  • DOI: 10.1145/1498765.1498785

Mesh Algorithms for PDE with Sieve I: Mesh Distribution
text, January 2009


Works referencing / citing this record:

Model reduction for fractured porous media: a machine learning approach for identifying main flow pathways
journal, March 2019

  • Srinivasan, Shriram; Karra, Satish; Hyman, Jeffrey
  • Computational Geosciences, Vol. 23, Issue 3
  • DOI: 10.1007/s10596-019-9811-7