99 research outputs found
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
Implications of the Choice of Quadrature Nodes for Picard Integral Deferred Corrections Methods for Ordinary Differential Equations
This paper concerns a class of deferred correction methods recently developed for initial value ordinary differential equations; such methods are based on a Picard integral form of the correction equation. These methods divide a given timestep [tn ,tn+1] into substeps, and use function values computed at these substeps to approximate the Picard integral by means of a numerical quadrature. The main purpose of this paper is to present a detailed analysis of the implications of the location of quadrature nodes on the accuracy and stability of the overall method. Comparisons between Gauss-Legendre, Gauss-Lobatto, Gauss-Radau, and uniformly spaced points are presented. Also, for a given set of quadrature nodes, quadrature rules may be formulated that include or exclude function values computed at the left-hand endpoint tn . Quadrature rules that do not depend on the left-hand endpoint (which are referred to as right-hand quadrature rules) are shown to lead to L(α)-stable implicit methods with α≈π/2. The semi-implicit analog of this property is also discussed. Numerical results suggest that the use of uniform quadrature nodes, as opposed to nodes based on Gaussian quadratures, does not significantly affect the stability or accuracy of these methods for orders less than ten. In contrast, a study of the reduction of order for stiff equations shows that when uniform quadrature nodes are used in conjunction with a right-hand quadrature rule, the form and extent of order-reduction changes considerably. Specifically, a reduction of order to is observed for uniform nodes as opposed to for non-uniform nodes, where Δt denotes the time step and ε a stiffness parameter such that ε→0 corresponds to the problem becoming increasingly stiff
A high-order spectral deferred correction strategy for low Mach number flow with complex chemistry
We present a fourth-order finite-volume algorithm in space and time for low
Mach number reacting flow with detailed kinetics and transport. Our temporal
integration scheme is based on a multi-implicit spectral deferred correction
(MISDC) strategy that iteratively couples advection, diffusion, and reactions
evolving subject to a constraint. Our new approach overcomes a stability
limitation of our previous second-order method encountered when trying to
incorporate higher-order polynomial representations of the solution in time to
increase accuracy. We have developed a new iterative scheme that naturally fits
within our MISDC framework that allows us to simultaneously conserve mass and
energy while satisfying on the equation of state. We analyse the conditions for
which the iterative schemes are guaranteed to converge to the fixed point
solution. We present numerical examples illustrating the performance of the new
method on premixed hydrogen, methane, and dimethyl ether flames.Comment: 27 pages, 5 figure
An Efficient Parallel-in-Time Method for Optimization with Parabolic PDEs
To solve optimization problems with parabolic PDE constraints, often methods
working on the reduced objective functional are used. They are computationally
expensive due to the necessity of solving both the state equation and a
backward-in-time adjoint equation to evaluate the reduced gradient in each
iteration of the optimization method. In this study, we investigate the use of
the parallel-in-time method PFASST in the setting of PDE constrained
optimization. In order to develop an efficient fully time-parallel algorithm we
discuss different options for applying PFASST to adjoint gradient computation,
including the possibility of doing PFASST iterations on both the state and
adjoint equations simultaneously. We also explore the additional gains in
efficiency from reusing information from previous optimization iterations when
solving each equation. Numerical results for both a linear and a non-linear
reaction-diffusion optimal control problem demonstrate the parallel speedup and
efficiency of different approaches
Torsion and ground state maxima: close but not the same
Could the location of the maximum point for a positive solution of a
semilinear Poisson equation on a convex domain be independent of the form of
the nonlinearity? Cima and Derrick found certain evidence for this surprising
conjecture.
We construct counterexamples on the half-disk, by working with the torsion
function and first Dirichlet eigenfunction. On an isosceles right triangle the
conjecture fails again. Yet the conjecture has merit, since the maxima of the
torsion function and eigenfunction are unexpectedly close together. It is an
open problem to quantify this closeness in terms of the domain and the
nonlinearity
A hybrid parareal spectral deferred corrections method
The parareal algorithm introduced in 2001 by Lions, Maday, and Turinici is an iterative method for the parallelization of the numerical solution of ordinary differential equations or partial differential equations discretized in the temporal direction. The temporal interval of interest is partitioned into successive domains which are assigned to separate processor units. Each iteration of the parareal algorithm consists of a high accuracy solution procedure performed in parallel on each domain using approximate initial conditions and a serial step which propagates a correction to the initial conditions through the entire time interval. The original method is designed to use classical single-step numerical methods for both of these steps. This paper investigates a variant of the parareal algorithm first outlined by Minion and Williams in 2008 that utilizes a deferred correction strategy within the parareal iterations. Here, the connections between parareal, parallel deferred corrections, and a hybrid parareal-spectral deferred correction method are further explored. The parallel speedup and efficiency of the hybrid methods are analyzed, and numerical results for ODEs and discretized PDEs are presented to demonstrate the performance of the hybrid approach
A Fourth-Order Adaptive Mesh Refinement Algorithm for the Multicomponent, Reacting Compressible Navier-Stokes Equations
In this paper we present a fourth-order in space and time block-structured
adaptive mesh refinement algorithm for the compressible multicomponent reacting
Navier-Stokes equations. The algorithm uses a finite volume approach that
incorporates a fourth-order discretization of the convective terms. The time
stepping algorithm is based on a multi-level spectral deferred corrections
method that enables explicit treatment of advection and diffusion coupled with
an implicit treatment of reactions. The temporal scheme is embedded in a
block-structured adaptive mesh refinement algorithm that includes subcycling in
time with spectral deferred correction sweeps applied on levels. Here we
present the details of the multi-level scheme paying particular attention to
the treatment of coarse-fine boundaries required to maintain fourth-order
accuracy in time. We then demonstrate the convergence properties of the
algorithm on several test cases including both nonreacting and reacting flows.
Finally we present simulations of a vitiated dimethyl ether jet in 2D and a
turbulent hydrogen jet in 3D, both with detailed kinetics and transport
A space-time parallel solver for the three-dimensional heat equation
The paper presents a combination of the time-parallel “parallel full approximation scheme in space and time” (PFASST) with a parallel multigrid method (PMG) in space, resulting in a mesh-based solver for the three-dimensional heat equation with a uniquely high degree of efficient concurrency. Parallel scaling tests are reported on the Cray XE6 machine “Monte Rosa” on up to cores and on the IBM Blue Gene/Q system “JUQUEEN” on up to cores. The efficacy of the combined spatial- and temporal parallelization is shown by demonstrating that using PFASST in addition to PMG significantly extends the strong-scaling limit. Implications of using spatial coarsening strategies in PFASST’s multi-level hierarchy in large-scale parallel simulations are discussed
- …