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

Title: A non-intrusive parallel-in-time adjoint solver with the XBraid library

Abstract

In this paper, an adjoint solver for the multigrid-in-time software library XBraid is presented. XBraid provides a non-intrusive approach for simulating unsteady dynamics on multiple processors while parallelizing not only in space but also in the time domain (XBraid: Parallel multigrid in time, http://llnl.gov/casc/xbraid). It applies an iterative multigrid reduction in time algorithm to existing spatially parallel classical time propagators and computes the unsteady solution parallel in time. Techniques from Automatic Differentiation are used to develop a consistent discrete adjoint solver which provides sensitivity information of output quantities with respect to design parameter changes. The adjoint code runs backwards through the primal XBraid actions and accumulates gradient information parallel in time. It is highly non-intrusive as existing adjoint time propagators can easily be integrated through the adjoint interface. The adjoint code is validated on advection-dominated flow with periodic upstream boundary condition. We report it provides similar strong scaling results as the primal XBraid solver and offers great potential for speeding up the overall computational costs for sensitivity analysis using multiple processors.

Authors:
 [1];  [1];  [2]
  1. Technische Universität Kaiserslautern (Germany)
  2. Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
Publication Date:
Research Org.:
Lawrence Livermore National Laboratory (LLNL), Livermore, CA (United States)
Sponsoring Org.:
USDOE National Nuclear Security Administration (NNSA)
OSTI Identifier:
1773254
Report Number(s):
LLNL-JRNL-730159
Journal ID: ISSN 1432-9360; 881043
Grant/Contract Number:  
AC52-07NA27344
Resource Type:
Accepted Manuscript
Journal Name:
Computing and Visualization in Science
Additional Journal Information:
Journal Volume: 19; Journal Issue: 3-4; Journal ID: ISSN 1432-9360
Publisher:
Springer
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING; Parallel-in-time; Multigrid-in-time; Parareal; Optimization; Adjoint sensitivity; Unsteady adjoint; High performance computing

Citation Formats

Günther, Stefanie, Gauger, Nicolas R., and Schroder, Jacob B. A non-intrusive parallel-in-time adjoint solver with the XBraid library. United States: N. p., 2018. Web. doi:10.1007/s00791-018-0300-7.
Günther, Stefanie, Gauger, Nicolas R., & Schroder, Jacob B. A non-intrusive parallel-in-time adjoint solver with the XBraid library. United States. https://doi.org/10.1007/s00791-018-0300-7
Günther, Stefanie, Gauger, Nicolas R., and Schroder, Jacob B. Wed . "A non-intrusive parallel-in-time adjoint solver with the XBraid library". United States. https://doi.org/10.1007/s00791-018-0300-7. https://www.osti.gov/servlets/purl/1773254.
@article{osti_1773254,
title = {A non-intrusive parallel-in-time adjoint solver with the XBraid library},
author = {Günther, Stefanie and Gauger, Nicolas R. and Schroder, Jacob B.},
abstractNote = {In this paper, an adjoint solver for the multigrid-in-time software library XBraid is presented. XBraid provides a non-intrusive approach for simulating unsteady dynamics on multiple processors while parallelizing not only in space but also in the time domain (XBraid: Parallel multigrid in time, http://llnl.gov/casc/xbraid). It applies an iterative multigrid reduction in time algorithm to existing spatially parallel classical time propagators and computes the unsteady solution parallel in time. Techniques from Automatic Differentiation are used to develop a consistent discrete adjoint solver which provides sensitivity information of output quantities with respect to design parameter changes. The adjoint code runs backwards through the primal XBraid actions and accumulates gradient information parallel in time. It is highly non-intrusive as existing adjoint time propagators can easily be integrated through the adjoint interface. The adjoint code is validated on advection-dominated flow with periodic upstream boundary condition. We report it provides similar strong scaling results as the primal XBraid solver and offers great potential for speeding up the overall computational costs for sensitivity analysis using multiple processors.},
doi = {10.1007/s00791-018-0300-7},
journal = {Computing and Visualization in Science},
number = 3-4,
volume = 19,
place = {United States},
year = {Wed Jun 27 00:00:00 EDT 2018},
month = {Wed Jun 27 00:00:00 EDT 2018}
}

Works referenced in this record:

Adjoint equations in CFD - Duality, boundary conditions and solution behaviour
conference, August 1997

  • Giles, M.; Pierce, N.; Giles, M.
  • 13th Computational Fluid Dynamics Conference
  • DOI: 10.2514/6.1997-1850

Evaluating Derivatives
book, January 2008


Robust optimization – A comprehensive survey
journal, July 2007

  • Beyer, Hans-Georg; Sendhoff, Bernhard
  • Computer Methods in Applied Mechanics and Engineering, Vol. 196, Issue 33-34
  • DOI: 10.1016/j.cma.2007.03.003

Analysis of the Parareal Time‐Parallel Time‐Integration Method
journal, January 2007

  • Gander, Martin J.; Vandewalle, Stefan
  • SIAM Journal on Scientific Computing, Vol. 29, Issue 2
  • DOI: 10.1137/05064607X

Multigrid Reduction in Time for Nonlinear Parabolic Problems: A Case Study
journal, January 2017

  • Falgout, R. D.; Manteuffel, T. A.; O'Neill, B.
  • SIAM Journal on Scientific Computing, Vol. 39, Issue 5
  • DOI: 10.1137/16M1082330

Optimum Shape Design for Unsteady Flows with Time-Accurate Continuous and Discrete Adjoint Method
journal, July 2007

  • Nadarajah, Siva K.; Jameson, Antony
  • AIAA Journal, Vol. 45, Issue 7
  • DOI: 10.2514/1.24332

Parallel methods for integrating ordinary differential equations
journal, December 1964


Adjoint Recovery of Superconvergent Functionals from PDE Approximations
journal, January 2000


Explicit parallel-in-time integration of a linear acoustic-advection system
journal, April 2012


Van der Pol oscillator
journal, January 2007


Analysis of Block Parareal Preconditioners for Parabolic Optimal Control Problems
journal, January 2010

  • Mathew, Tarek P.; Sarkis, Marcus; Schaerer, Christian E.
  • SIAM Journal on Scientific Computing, Vol. 32, Issue 3
  • DOI: 10.1137/080717481

A time-domain decomposition iterative method for the solution of distributed linear quadratic optimal control problems
journal, January 2005

  • Heinkenschloss, Matthias
  • Journal of Computational and Applied Mathematics, Vol. 173, Issue 1
  • DOI: 10.1016/j.cam.2004.03.005

Unsteady Aerodynamic Design on Unstructured Meshes with Sliding Interfaces
conference, January 2013

  • Economon, Thomas; Palacios, Francisco; Alonso, Juan
  • 51st AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition
  • DOI: 10.2514/6.2013-632

Domain decomposition in time for PDE-constrained optimization
journal, December 2015


Discrete Adjoint-Based Design Optimization of Unsteady Turbulent Flows on Dynamic Unstructured Grids
journal, June 2010

  • Nielsen, Eric J.; Diskin, Boris; Yamaleev, Nail K.
  • AIAA Journal, Vol. 48, Issue 6
  • DOI: 10.2514/1.J050035

Time dependent calculations using multigrid, with applications to unsteady flows past airfoils and wings
conference, August 1991


A General Framework for the Optimal Control of Unsteady Flows with Applications
conference, June 2012

  • Rumpfkeil, Markus; Zingg, David
  • 45th AIAA Aerospace Sciences Meeting and Exhibit
  • DOI: 10.2514/6.2007-1128

Adjoint methods for PDEs: a posteriori error analysis and postprocessing by duality
journal, January 2002


Two-Level Convergence Theory for Multigrid Reduction in Time (MGRIT)
journal, January 2017

  • Dobrev, V. A.; Kolev, Tz.; Petersson, N. A.
  • SIAM Journal on Scientific Computing, Vol. 39, Issue 5
  • DOI: 10.1137/16M1074096

Résolution d'EDP par un schéma en temps «pararéel »
journal, April 2001

  • Lions, Jacques-Louis; Maday, Yvon; Turinici, Gabriel
  • Comptes Rendus de l'Académie des Sciences - Series I - Mathematics, Vol. 332, Issue 7
  • DOI: 10.1016/S0764-4442(00)01793-6

A note on MGR methods
journal, February 1983


Grid Adaptation for Functional Outputs: Application to Two-Dimensional Inviscid Flows
journal, February 2002

  • Venditti, David A.; Darmofal, David L.
  • Journal of Computational Physics, Vol. 176, Issue 1
  • DOI: 10.1006/jcph.2001.6967

Parallel Time Integration with Multigrid
journal, January 2014

  • Falgout, R. D.; Friedhoff, S.; Kolev, Tz. V.
  • SIAM Journal on Scientific Computing, Vol. 36, Issue 6
  • DOI: 10.1137/130944230

50 Years of Time Parallel Time Integration
book, January 2015


Toward an efficient parallel in time method for partial differential equations
journal, January 2012

  • Emmett, Matthew; Minion, Michael
  • Communications in Applied Mathematics and Computational Science, Vol. 7, Issue 1
  • DOI: 10.2140/camcos.2012.7.105

An Introduction to the Adjoint Approach to Design
journal, December 2000

  • Giles, Michael B.; Pierce, Niles A.
  • Flow, Turbulence and Combustion, Vol. 65, Issue 3/4, p. 393-415
  • DOI: 10.1023/A:1011430410075

Simultaneous single-step one-shot optimization with unsteady PDEs
journal, March 2016

  • Günther, Stefanie; Gauger, Nicolas R.; Wang, Qiqi
  • Journal of Computational and Applied Mathematics, Vol. 294
  • DOI: 10.1016/j.cam.2015.07.033

Solution of the Unsteady Discrete Adjoint for Three-Dimensional Problems on Dynamically Deforming Unstructured Meshes
conference, June 2012

  • Mavriplis, Dimitri
  • 46th AIAA Aerospace Sciences Meeting and Exhibit
  • DOI: 10.2514/6.2008-727

Time-decomposed parallel time-integrators: theory and feasibility studies for fluid, structure, and fluid-structure applications
journal, January 2003

  • Farhat, Charbel; Chandesris, Marion
  • International Journal for Numerical Methods in Engineering, Vol. 58, Issue 9
  • DOI: 10.1002/nme.860

Multi-level adaptive solutions to boundary-value problems
journal, May 1977


Aerodynamic design via control theory
journal, September 1988

  • Jameson, Antony
  • Journal of Scientific Computing, Vol. 3, Issue 3
  • DOI: 10.1007/BF01061285

Stabilization of Unstable Procedures: The Recursive Projection Method
journal, August 1993

  • Shroff, Gautam M.; Keller, Herbert B.
  • SIAM Journal on Numerical Analysis, Vol. 30, Issue 4
  • DOI: 10.1137/0730057

Assessment of the Recursive Projection Method for the Stabilization of Discrete Adjoint Solvers
conference, June 2017

  • Albring, Tim A.; Dick, Thomas; Gauger, Nicolas R.
  • 18th AIAA/ISSMO Multidisciplinary Analysis and Optimization Conference
  • DOI: 10.2514/6.2017-3664

A note on MGR methods
journal, February 1983


Practical and theoretical aspects of adjoint parameter estimation and identifiability in meteorology and oceanography
journal, January 1998


On optimum design in fluid mechanics
journal, June 1974


Evaluating Derivatives
book, January 2008


Parallel methods for integrating ordinary differential equations
journal, December 1964


A General Framework for the Optimal Control of Unsteady Flows with Applications
conference, June 2012

  • Rumpfkeil, Markus; Zingg, David
  • 45th AIAA Aerospace Sciences Meeting and Exhibit
  • DOI: 10.2514/6.2007-1128

Solution of the Unsteady Discrete Adjoint for Three-Dimensional Problems on Dynamically Deforming Unstructured Meshes
conference, June 2012

  • Mavriplis, Dimitri
  • 46th AIAA Aerospace Sciences Meeting and Exhibit
  • DOI: 10.2514/6.2008-727