Adaptive time-integration for goal-oriented and coupled problems
(2021)- Abstract
- We consider efficient methods for the partitioned time-integration of multiphysics problems, which commonly exhibit a multiscale behavior, requiring independent time-grids. Examples are fluid structure interaction in e.g., the simulation of blood-flow or cooling of rocket engines, or ocean-atmosphere-vegetation interaction. The ideal method for solving these problems allows independent and adaptive time-grids, higher order time-discretizations, is fast and robust, and allows the coupling of existing subsolvers, executed in parallel.
We consider Waveform relaxation (WR) methods, which can have all of these properties. WR methods iterate on continuous-in-time interface functions, obtained via suitable interpolation. The difficulty... (More) - We consider efficient methods for the partitioned time-integration of multiphysics problems, which commonly exhibit a multiscale behavior, requiring independent time-grids. Examples are fluid structure interaction in e.g., the simulation of blood-flow or cooling of rocket engines, or ocean-atmosphere-vegetation interaction. The ideal method for solving these problems allows independent and adaptive time-grids, higher order time-discretizations, is fast and robust, and allows the coupling of existing subsolvers, executed in parallel.
We consider Waveform relaxation (WR) methods, which can have all of these properties. WR methods iterate on continuous-in-time interface functions, obtained via suitable interpolation. The difficulty is to find suitable convergence acceleration, which is required for the iteration converge quickly.
We develop a fast and highly robust, second order in time, adaptive WR method for unsteady thermal fluid structure interaction (FSI), modelled by heterogeneous coupled linear heat equations. We use a Dirichlet-Neumann coupling at the interface and an analytical optimal relaxation parameter derived for the fully-discrete scheme. While this method is sequential, it is notably faster and more robust than similar parallel methods.
We further develop a novel, parallel WR method, using asynchronous communication techniques during time-integration to accelerate convergence. Instead of exchanging interpolated time-dependent functions at the end of each time-window or iteration, we exchange time-point data immediately after each timestep. The analytical description and convergence results of this method generalize existing WR theory.
Since WR methods allow coupling of problems in a relative black-box manner, we developed adapters to PDE-subsolvers implemented using DUNE and FEniCS. We demonstrate this coupling in a thermal FSI test case.
Lastly, we consider adaptive time-integration for goal-oriented problems, where one is interested in a quantity of interest (QoI), which is a functional of the solution. The state-of-the-art method is the dual-weighted residual (DWR) method, which is extremely costly in both computation and implementation. We develop a goal oriented adaptive method based on local error estimates, which is considerably cheaper in computation. We prove convergence of the error in the QoI for tolerance to zero under a controllability assumption. By analyzing global error propagation with respect to the QoI, we can identify possible issues and make performance predictions. Numerical results verify these results and show our method to be more efficient than the DWR method. (Less)
Please use this url to cite or link to this publication:
https://lup.lub.lu.se/record/549d68ec-4f0a-47c4-8677-63fb4da3f1e9
- author
- Meisrimel, Peter LU
- supervisor
-
- Philipp Birken LU
- Claus Führer LU
- opponent
-
- Associate professor (Reader) Dedner, Andreas, University of Warwick, UK
- organization
- publishing date
- 2021-06-03
- type
- Thesis
- publication status
- published
- subject
- keywords
- time-adaptivity, goal-oriented problems, error estimation, coupled problems, waveform relaxation, asynchronous methods, conjugate heat transfer
- pages
- 176 pages
- publisher
- Lund University
- defense location
- MH:H, Lund. Join via zoom: https://lu-se.zoom.us/j/65642304248
- defense date
- 2021-09-02 14:00:00
- ISBN
- 978-91-7895-930-3
- 978-91-7895-929-7
- project
- eSSENCE@LU 4:4 - 3M: Multiphysics, multicore, multirate solution of coupled dynamic problems
- language
- English
- LU publication?
- yes
- id
- 549d68ec-4f0a-47c4-8677-63fb4da3f1e9
- date added to LUP
- 2021-06-03 13:11:47
- date last changed
- 2021-11-26 08:27:40
@phdthesis{549d68ec-4f0a-47c4-8677-63fb4da3f1e9,
abstract = {{We consider efficient methods for the partitioned time-integration of multiphysics problems, which commonly exhibit a multiscale behavior, requiring independent time-grids. Examples are fluid structure interaction in e.g., the simulation of blood-flow or cooling of rocket engines, or ocean-atmosphere-vegetation interaction. The ideal method for solving these problems allows independent and adaptive time-grids, higher order time-discretizations, is fast and robust, and allows the coupling of existing subsolvers, executed in parallel. <br/><br/>We consider Waveform relaxation (WR) methods, which can have all of these properties. WR methods iterate on continuous-in-time interface functions, obtained via suitable interpolation. The difficulty is to find suitable convergence acceleration, which is required for the iteration converge quickly. <br/><br/>We develop a fast and highly robust, second order in time, adaptive WR method for unsteady thermal fluid structure interaction (FSI), modelled by heterogeneous coupled linear heat equations. We use a Dirichlet-Neumann coupling at the interface and an analytical optimal relaxation parameter derived for the fully-discrete scheme. While this method is sequential, it is notably faster and more robust than similar parallel methods.<br/><br/>We further develop a novel, parallel WR method, using asynchronous communication techniques during time-integration to accelerate convergence. Instead of exchanging interpolated time-dependent functions at the end of each time-window or iteration, we exchange time-point data immediately after each timestep. The analytical description and convergence results of this method generalize existing WR theory.<br/><br/>Since WR methods allow coupling of problems in a relative black-box manner, we developed adapters to PDE-subsolvers implemented using DUNE and FEniCS. We demonstrate this coupling in a thermal FSI test case.<br/><br/>Lastly, we consider adaptive time-integration for goal-oriented problems, where one is interested in a quantity of interest (QoI), which is a functional of the solution. The state-of-the-art method is the dual-weighted residual (DWR) method, which is extremely costly in both computation and implementation. We develop a goal oriented adaptive method based on local error estimates, which is considerably cheaper in computation. We prove convergence of the error in the QoI for tolerance to zero under a controllability assumption. By analyzing global error propagation with respect to the QoI, we can identify possible issues and make performance predictions. Numerical results verify these results and show our method to be more efficient than the DWR method.}},
author = {{Meisrimel, Peter}},
isbn = {{978-91-7895-930-3}},
keywords = {{time-adaptivity; goal-oriented problems; error estimation; coupled problems; waveform relaxation; asynchronous methods; conjugate heat transfer}},
language = {{eng}},
month = {{06}},
publisher = {{Lund University}},
school = {{Lund University}},
title = {{Adaptive time-integration for goal-oriented and coupled problems}},
url = {{https://lup.lub.lu.se/search/files/99039864/Peter_Meisrimel_Dissertation_Summary.pdf}},
year = {{2021}},
}