A Discontinuous Galerkin Method with Penalty for One-Dimensional Nonlocal Diffusion Problems

There have been many theoretical studies and numerical investigations of nonlocal diffusion(ND)problems in recent years.In this paper,we propose and analyze a new discontinuous Galerkin method for solving one-dimensional steady-state and time-dependent ND problems,based on a formulation that directly penalizes the jumps across the element interfaces in the nonlocal sense.We show that the proposed discontinuous Galerkin scheme is stable and convergent.Moreover,the local limit of such DG scheme recovers classical DG scheme for the corresponding local diff usion problem,which is a distinct feature of the new formulation and assures the asymptotic compatibility of the discretization.Numerical tests are also presented to demonstrate the eff ectiveness and the robustness of the proposed method.

High Order Semi-implicit Multistep Methods for Time-Dependent Partial Differential Equations

We consider the construction of semi-implicit linear multistep methods that can be applied to time-dependent PDEs where the separation of scales in additive form,typically used in implicit-explicit(IMEX)methods,is not possible.As shown in Boscarino et al.(J.Sci.Comput.68:975-1001,2016)for Runge-Kutta methods,these semi-implicit techniques give a great flexibility,and allow,in many cases,the construction of simple linearly implicit schemes with no need of iterative solvers.In this work,we develop a general setting for the construction of high order semi-implicit linear multistep methods and analyze their stability properties for a prototype lineal'advection-diffusion equation and in the setting of strong stability preserving(SSP)methods.Our findings are demonstrated on several examples,including nonlinear reaction-diffusion and convection-diffusion problems.

Estimation of Decision Alternatives on the Basis of Interval Pairwise Comparison Matrices

This paper deals with the calculation of a vector of reliable weights of decision alternatives on the basis of interval pairwise comparison judgments of experts. These weights are used to construct the ranking of decision alternatives and to solve selection problems, problems of ratings construction, resources allocation problems, scenarios evaluation problems, and other decision making problems. A comparative analysis of several popular models, which calculate interval weights on the basis of interval pairwise comparison matrices (IPCMs), was performed. The features of these models when they are applied to IPCMs with different inconsistency levels were identified. An algorithm is proposed which contains the stages for analyzing and increasing the IPCM inconsistency, calculating normalized interval weights, and calculating the ranking of decision alternatives on the basis of the resulting interval weights. It was found that the property of weak order preservation usually allowed identifying order-related intransitive expert pairwise comparison judgments. The correction of these elements leads to the removal of contradictions in resulting weights and increases the accuracy and reliability of results.

Strong Skew Commutativity Preserving Maps on Rings with Involution

Let R be a unital ＊-ring with the unit I. Assume that R contains a symmetric idempotent P which satisfies ARP= 0 implies A = 0 and AR（I - P） = 0 implies A = 0. In this paper, it is shown that a surjective map Ф： R →R is strong skew commutativity preserving （that is, satisfies Ф（A）Ф（B） - Ф（B）Ф（A）＊ ： AB- BA＊ for all A, B ∈R） if and only if there exist a map f ： R → ZSz（R） and an element Z ∈ ZS（R） with Z^2 =I such that Ф（A） =ZA ＋ f（A） for all A ∈ R, where ZS（R） is the symmetric center of R. As applications, the strong skew commutativity preserving maps on unital prime ＊-rings and von Neumann algebras with no central summands of type I1 are characterized.

STRONG STABILITY PRESERVING PROPERTY OF THE DEFERRED CORRECTION TIME DISCRETIZATION

In this paper, we study the strong stability preserving （SSP） property of a class of deferred correction time discretization methods, for solving the method-of-lines schemes approximating hyperbolic partial differential equations.

HERMITE WENO SCHEMES WITH STRONG STABILITY PRESERVING MULTI-STEP TEMPORAL DISCRETIZATION METHODS FOR CONSERVATION LAWS

Based on the work of Shu [SIAM J. Sci. Stat. Comput, 9 （1988）, pp.1073-1084], we construct a class of high order multi-step temporal discretization procedure for finite volume Hermite weighted essential non-oscillatory （HWENO） methods to solve hyperbolic conservation laws. The key feature of the multi-step temporal discretization procedure is to use variable time step with strong stability preserving （SSP）. The multi-step tem- poral discretization methods can make full use of computed information with HWENO spatial discretization by holding the former computational values. Extensive numerical experiments are presented to demonstrate that the finite volume HWENO schemes with multi-step diseretization can achieve high order accuracy and maintain non-oscillatory properties near discontinuous region of the solution.

A Second-Order Implicit-Explicit Scheme for the Baroclinic-Barotropic Split System of Primitive Equations

The baroclinic-barotropic mode splitting technique is commonly employed in numerical solutions of the primitive equations for ocean modeling to deal with the multiple time scales of ocean dynamics.In this paper,a second-order implicit-explicit(IMEX)scheme is proposed to advance the baroclinic-barotropic split system.Specifically,the baroclinic mode and the layer thickness of fluid are evolved explicitly via the second-order strong stability preserving Runge-Kutta scheme,while the barotropic mode is advanced implicitly using the linearized Crank-Nicolson scheme.At each time step,the baroclinic velocity is first computed using an intermediate barotropic velocity.The barotropic velocity is then corrected by re-advancing the barotropic mode with an improved barotropic forcing.Finally,the layer thickness is updated by coupling the baroclinic and barotropic velocities together.In addition,numerical inconsistencies on the discretized sea surface height caused by the mode splitting are alleviated via a reconciliation process with carefully calculated flux deficits.Temporal truncation error is also analyzed to validate the second-order accuracy of the scheme.Finally,two benchmark tests from the MPAS-Ocean platform are conducted to numerically demonstrate the performance of the proposed IMEX scheme.

ODE-Based Multistep Schemes for Backward Stochastic Differential Equations

In this paper,we explore a new approach to design and analyze numerical schemes for backward stochastic differential equations(BSDEs).By the nonlinear Feynman-Kac formula,we reformulate the BSDE into a pair of reference ordinary differential equations(ODEs),which can be directly discretized by many standard ODE solvers,yielding the corresponding numerical schemes for BSDEs.In particular,by applying strong stability preserving(SSP)time discretizations to the reference ODEs,we can propose new SSP multistep schemes for BSDEs.Theoretical analyses are rigorously performed to prove the consistency,stability and convergency of the proposed SSP multistep schemes.Numerical experiments are further carried out to verify our theoretical results and the capacity of the proposed SSP multistep schemes for solving complex associated problems.