211 research outputs found
Hybridizable compatible finite element discretizations for numerical weather prediction: implementation and analysis
There is a current explosion of interest in new numerical methods for atmospheric modeling. A driving force behind this is the need to be able to simulate, with high efficiency, large-scale geophysical flows on increasingly more parallel computer systems. Many current operational models, including that of the UK Met Office, depend on orthogonal meshes, such as the latitude-longitude grid. This facilitates the development of finite difference discretizations with favorable numerical properties. However, such methods suffer from the ``pole problem," which prohibits the model to make efficient use of a large number of computing processors due to excessive concentration of grid-points at the poles.
Recently developed finite element discretizations, known as ``compatible" finite elements, avoid this issue while maintaining the key numerical properties essential for accurate geophysical simulations. Moreover, these properties can be obtained on arbitrary, non-orthogonal meshes. However, the efficient solution of the resulting discrete systems depend on transforming the mixed velocity-pressure (or velocity-pressure-buoyancy) system into an elliptic problem for the pressure. This is not so straightforward within the compatible finite element framework due to inter-element coupling.
This thesis supports the proposition that systems arising from compatible finite element discretizations can be solved efficiently using a technique known as ``hybridization." Hybridization removes inter-element coupling while maintaining the desired numerical properties. This permits the construction of sparse, elliptic problems, for which fast solver algorithms are known, using localized algebra. We first introduce the technique for compatible finite element discretizations of simplified atmospheric models. We then develop a general software abstraction for the rapid implementation and composition of hybridization methods, with an emphasis on preconditioning.
Finally, we extend the technique for a new compatible method for the full, compressible atmospheric equations used in operational models.Open Acces
Multilevel convergence analysis of multigrid-reduction-in-time
This paper presents a multilevel convergence framework for
multigrid-reduction-in-time (MGRIT) as a generalization of previous two-grid
estimates. The framework provides a priori upper bounds on the convergence of
MGRIT V- and F-cycles, with different relaxation schemes, by deriving the
respective residual and error propagation operators. The residual and error
operators are functions of the time stepping operator, analyzed directly and
bounded in norm, both numerically and analytically. We present various upper
bounds of different computational cost and varying sharpness. These upper
bounds are complemented by proposing analytic formulae for the approximate
convergence factor of V-cycle algorithms that take the number of fine grid time
points, the temporal coarsening factors, and the eigenvalues of the time
stepping operator as parameters.
The paper concludes with supporting numerical investigations of parabolic
(anisotropic diffusion) and hyperbolic (wave equation) model problems. We
assess the sharpness of the bounds and the quality of the approximate
convergence factors. Observations from these numerical investigations
demonstrate the value of the proposed multilevel convergence framework for
estimating MGRIT convergence a priori and for the design of a convergent
algorithm. We further highlight that observations in the literature are
captured by the theory, including that two-level Parareal and multilevel MGRIT
with F-relaxation do not yield scalable algorithms and the benefit of a
stronger relaxation scheme. An important observation is that with increasing
numbers of levels MGRIT convergence deteriorates for the hyperbolic model
problem, while constant convergence factors can be achieved for the diffusion
equation. The theory also indicates that L-stable Runge-Kutta schemes are more
amendable to multilevel parallel-in-time integration with MGRIT than A-stable
Runge-Kutta schemes.Comment: 26 pages; 17 pages Supplementary Material
- …