Optimization of the parallel semi‐Lagrangian scheme based on overlapping communication with computation in the YHGSM

Author(s):  
Tao Jiang ◽  
Jianping Wu ◽  
Zhaoyang Liu ◽  
Wenpeng Zhao ◽  
Yongshun Zhang
Keyword(s):  
Author(s):  
Stefano Almi ◽  
Marco Morandotti ◽  
Francesco Solombrino

AbstractA multi-step Lagrangian scheme at discrete times is proposed for the approximation of a nonlinear continuity equation arising as a mean-field limit of spatially inhomogeneous evolutionary games, describing the evolution of a system of spatially distributed agents with strategies, or labels, whose payoff depends also on the current position of the agents. The scheme is Lagrangian, as it traces the evolution of position and labels along characteristics, and is a multi-step scheme, as it develops on the following two stages: First, the distribution of strategies or labels is updated according to a best performance criterion, and then, this is used by the agents to evolve their position. A general convergence result is provided in the space of probability measures. In the special cases of replicator-type systems and reversible Markov chains, variants of the scheme, where the explicit step in the evolution of the labels is replaced by an implicit one, are also considered and convergence results are provided.


2020 ◽  
Vol 35 (6) ◽  
pp. 355-366
Author(s):  
Vladimir V. Shashkin ◽  
Gordey S. Goyman

AbstractThis paper proposes the combination of matrix exponential method with the semi-Lagrangian approach for the time integration of shallow water equations on the sphere. The second order accuracy of the developed scheme is shown. Exponential semi-Lagrangian scheme in the combination with spatial approximation on the cubed-sphere grid is verified using the standard test problems for shallow water models. The developed scheme is as good as the conventional semi-implicit semi-Lagrangian scheme in accuracy of slowly varying flow component reproduction and significantly better in the reproduction of the fast inertia-gravity waves. The accuracy of inertia-gravity waves reproduction is close to that of the explicit time-integration scheme. The computational efficiency of the proposed exponential semi-Lagrangian scheme is somewhat lower than the efficiency of semi-implicit semi-Lagrangian scheme, but significantly higher than the efficiency of explicit, semi-implicit, and exponential Eulerian schemes.


1998 ◽  
Vol 11 (1) ◽  
pp. 376-376
Author(s):  
S.G. Moiseenko

Results of 2D numerical simulation of the magneto rotational mechanism of a supernova explosion are presented. Simulation has been done for the real equations of state and neutrino energy losses have been taken into account. Simulation has been done on the basis of an Implicit Lagrangian scheme on atriangular grid with grid reconstructuring. It is shown that, due to differential rotation of the star, a toroidal component of the magnetic field appears and grows with time. Rotational momentum transfers outwards as the toroidal component grows with time. With the evolution of the process, part of the envelope of the star is ejected. The amounts of the thrown-off mass and energy are estimated. The results of the simulation could be used as a possible explanation for the supernova explosion picture.


2005 ◽  
Vol 22 (7) ◽  
pp. 886-895 ◽  
Author(s):  
F. Ardhuin ◽  
T. H. C. Herbers

Abstract A new semi-Lagrangian advection scheme called multistep ray advection is proposed for solving the spectral energy balance equation of ocean surface gravity waves. Existing so-called piecewise ray methods advect wave energy over a single time step using “pieces” of ray trajectories, after which the spectrum is updated with source terms representing various physical processes. The generalized scheme presented here allows for an arbitrary number N of advection time steps along the same rays, thus reducing numerical diffusion, and still including source-term variations every time step. Tests are performed for alongshore uniform bottom topography, and the effects of two types of discretizations of the wave spectrum are investigated, a finite-bandwidth representation and a single frequency and direction per spectral band. In the limit of large N, both the accuracy and computation cost of the method increase, approaching a nondiffusive fully Lagrangian scheme. Even for N = 1 the semi-Lagrangian scheme test results show less numerical diffusion than predictions of the commonly used first-order upwind finite-difference scheme. Application to the refraction and shoaling of narrow swell spectra across a continental shelf illustrates the importance of controlling numerical diffusion. Numerical errors in a single-step (Δt = 600 s) scheme implemented on the North Carolina continental shelf (typical swell propagation time across the shelf is about 3 h) are shown to be comparable to the angular diffusion predicted by the wave–bottom Bragg scattering theory, in particular for narrow directional spectra, suggesting that the true directional spread of swell may not always be resolved in existing wave prediction models, because of excessive numerical diffusion. This diffusion is effectively suppressed in cases presented here with a four-step semi-Lagrangian scheme, using the same value of Δt.


2013 ◽  
Vol 6 (4) ◽  
pp. 1029-1042 ◽  
Author(s):  
B. Sørensen ◽  
E. Kaas ◽  
U. S. Korsholm

Abstract. In this paper a new advection scheme for the online coupled chemical–weather prediction model Enviro-HIRLAM is presented. The new scheme is based on the locally mass-conserving semi-Lagrangian method (LMCSL), where the original two-dimensional scheme has been extended to a fully three-dimensional version. This means that the three-dimensional semi-implicit semi-Lagrangian scheme which is currently used in Enviro-HIRLAM is largely unchanged. The HIRLAM model is a computationally efficient hydrostatic operational short-term numerical weather prediction model, which is used as the base for the online integrated Enviro-HIRLAM. The new scheme is shown to be efficient, mass conserving, and shape preserving, while only requiring minor alterations to the original code. It still retains the stability at long time steps, which the semi-Lagrangian schemes are known for, while handling the emissions of chemical species accurately. Several mass-conserving filters have been tested to assess the optimal balance of accuracy vs. efficiency.


Sign in / Sign up

Export Citation Format

Share Document