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:
-
- Technische Universität Kaiserslautern (Germany)
- 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
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
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
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
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
Parallel methods for integrating ordinary differential equations
journal, December 1964
- Nievergelt, J.
- Communications of the ACM, Vol. 7, Issue 12
Adjoint Recovery of Superconvergent Functionals from PDE Approximations
journal, January 2000
- Pierce, Niles A.; Giles, Michael B.
- SIAM Review, Vol. 42, Issue 2
Explicit parallel-in-time integration of a linear acoustic-advection system
journal, April 2012
- Ruprecht, D.; Krause, R.
- Computers & Fluids, Vol. 59
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
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
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
Domain decomposition in time for PDE-constrained optimization
journal, December 2015
- Barker, Andrew T.; Stoll, Martin
- Computer Physics Communications, Vol. 197
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
Time dependent calculations using multigrid, with applications to unsteady flows past airfoils and wings
conference, August 1991
- Jameson, Antony
- 10th Computational Fluid Dynamics Conference
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
Adjoint methods for PDEs: a posteriori error analysis and postprocessing by duality
journal, January 2002
- Giles, Michael B.; Süli, Endre
- Acta Numerica, Vol. 11
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
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
A note on MGR methods
journal, February 1983
- Ries, Manfred; Trottenberg, Ulrich; Winter, Gerd
- Linear Algebra and its Applications, Vol. 49
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
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
50 Years of Time Parallel Time Integration
book, January 2015
- Gander, Martin J.
- Contributions in Mathematical and Computational Sciences
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
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
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
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
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
Multi-level adaptive solutions to boundary-value problems
journal, May 1977
- Brandt, Achi
- Mathematics of Computation, Vol. 31, Issue 138
Aerodynamic design via control theory
journal, September 1988
- Jameson, Antony
- Journal of Scientific Computing, Vol. 3, Issue 3
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
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
A note on MGR methods
journal, February 1983
- Ries, Manfred; Trottenberg, Ulrich; Winter, Gerd
- Linear Algebra and its Applications, Vol. 49
Practical and theoretical aspects of adjoint parameter estimation and identifiability in meteorology and oceanography
journal, January 1998
- Navon, I. M.
- Dynamics of Atmospheres and Oceans, Vol. 27, Issue 1-4
On optimum design in fluid mechanics
journal, June 1974
- Pironneau, O.
- Journal of Fluid Mechanics, Vol. 64, Issue 1
Parallel methods for integrating ordinary differential equations
journal, December 1964
- Nievergelt, J.
- Communications of the ACM, Vol. 7, Issue 12
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
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