58,559 research outputs found
A rational deferred correction approach to parabolic optimal control problems
The accurate and efficient solution of time-dependent PDE-constrained optimization problems is a challenging task, in large part due to the very high dimension of the matrix systems that need to be solved. We devise a new deferred correction method for coupled systems of time-dependent PDEs, allowing one to iteratively improve the accuracy of low-order time stepping schemes. We consider two variants of our method, a splitting and a coupling version, and analyze their convergence properties. We then test our approach on a number of PDE-constrained optimization problems. We obtain solution accuracies far superior to that achieved when solving a single discretized problem, in particular in cases where the accuracy is limited by the time discretization. Our approach allows for the direct reuse of existing solvers for the resulting matrix systems, as well as state-of-the-art preconditioning strategies
Parallel-in-Time Multi-Level Integration of the Shallow-Water Equations on the Rotating Sphere
The modeling of atmospheric processes in the context of weather and climate
simulations is an important and computationally expensive challenge. The
temporal integration of the underlying PDEs requires a very large number of
time steps, even when the terms accounting for the propagation of fast
atmospheric waves are treated implicitly. Therefore, the use of
parallel-in-time integration schemes to reduce the time-to-solution is of
increasing interest, particularly in the numerical weather forecasting field.
We present a multi-level parallel-in-time integration method combining the
Parallel Full Approximation Scheme in Space and Time (PFASST) with a spatial
discretization based on Spherical Harmonics (SH). The iterative algorithm
computes multiple time steps concurrently by interweaving parallel high-order
fine corrections and serial corrections performed on a coarsened problem. To do
that, we design a methodology relying on the spectral basis of the SH to
coarsen and interpolate the problem in space. The methods are evaluated on the
shallow-water equations on the sphere using a set of tests commonly used in the
atmospheric flow community. We assess the convergence of PFASST-SH upon
refinement in time. We also investigate the impact of the coarsening strategy
on the accuracy of the scheme, and specifically on its ability to capture the
high-frequency modes accumulating in the solution. Finally, we study the
computational cost of PFASST-SH to demonstrate that our scheme resolves the
main features of the solution multiple times faster than the serial schemes
A multigrid perspective on the parallel full approximation scheme in space and time
For the numerical solution of time-dependent partial differential equations,
time-parallel methods have recently shown to provide a promising way to extend
prevailing strong-scaling limits of numerical codes. One of the most complex
methods in this field is the "Parallel Full Approximation Scheme in Space and
Time" (PFASST). PFASST already shows promising results for many use cases and
many more is work in progress. However, a solid and reliable mathematical
foundation is still missing. We show that under certain assumptions the PFASST
algorithm can be conveniently and rigorously described as a multigrid-in-time
method. Following this equivalence, first steps towards a comprehensive
analysis of PFASST using block-wise local Fourier analysis are taken. The
theoretical results are applied to examples of diffusive and advective type
- …