守恒格式稳定性分析与耗散守恒格式

本文从守恒格式出发,建立分析稳定性和耗散性的启发性方法和Fourier分析方法,给出了耗散守恒格式的严格定义及三点耗散守恒格式的充要条件。应用本文的方法,重新分析了三点格式,得到如下结论:某些常系数耗散格式,在某些情况下,之所以会得到非物理解或发生非线性不稳定,是由于该格式在这些情况下,已经是零耗散的或是负耗散     

非正交网格上满足极值原理的扩散格式

构造了非正交网格上扩散方程新的非线性单元中心型有限体积格式,证明了该格式满足离散极值原理,且在适当条件下具有强制性、以及在离散H1范数下解的有界性和一阶收敛性.     

New approach to the Lax‐Wendroff modified differential equation for linear and nonlinear advection

The paper presents an enhanced analysis of the Lax‐Wendroff difference scheme—up to the eighth‐order with respect to time and space derivatives—of the modified‐partial differential equation (MDE) of the constant‐wind‐speed advection equation. The modified equation has been so far derived mainly as a fourth‐order equation. The Π ‐form of the first differential approximation (differential approximation or equivalent equation) derived by expressing the time derivatives in terms of the space derivatives is used for presenting the MDE. The obtained coefficients at higher order derivatives are analyzed for indications of the character of the dissipative and dispersive errors. The authors included a part of the stencil applied for determining the modified differential equation up to the eighth‐order of the analyzed modified differential equation for the second‐order Lax‐Wendroff scheme. Neither the derived coefficients at the space derivatives of order p ∈ (7 – 8) in the modified differential equation for the Lax‐Wendroff difference scheme nor the results of analyses on the basis of these coefficients of the group velocity, phase shift errors, or dispersive and dissipative features of the scheme have been published. The MDEs for 2 two‐step variants of the Lax‐Wendroff type difference schemes and the MacCormack predictor–corrector scheme (see MacCormack's study) constructed for the scalar hyperbolic conservation laws are also presented in this paper. The analysis of the inviscid Burgers equation solution with the initial condition in a form of a shock wave has been discussed on their basis. The inviscid Burgers equation with the source is also presented. The theory of MDE started to develop after the paper of C. W. Hirt was published in 1968.     

二维线性双曲型方程Neumann边值问题的紧交替方向隐格式

对二维Neumann边界条件的线性双曲型方程建立了紧交替方向的隐格式.利用方程和边界条件得到在空间上的三阶与五阶导数的边界值,进而在内点、边界内点和边界角点分别建立9点、6点和4点紧差分格式;通过引进新的范数和L₂范数估计L_∞范数;借助能量估计、Gronwall不等式和Schwarz不等式等技巧,详细分析了差分格式在无穷范数下关于时间和空间分别为二阶和四阶收敛性,并给出了稳定性结果;通过数值算例,验证了理论分析结果.     

基于概率的保持前缀地址随机化算法的安全评估

地址随机化算法通常用于在发布流量数据之前进行去隐私处理.保持前缀地址随机化算法就是其中一个常用算法.对于保持前缀地址随机化算法而言,由于引入了更多的限制,因此也面临更多的安全风险.分析了相关性攻击对保持前缀地址随机化算法的安全影响,并利用概率分析和仿真评估了不同攻击方法对其安全性能的影响.     

线性初边值问题差分格式的稳定性和收敛性(英文)

本文讨论一般的方程系数满足Lipschitz条件的变系数线性双曲型初边值问题差分格式的稳定性,并在很弱的条件下证明了几类差分格式是稳定的,本文还证明了:如果格式稳定,则在解和方程的系数足够光滑时,差分解将收敛于微分方程的解,并且g阶格式在l_2空间有g阶收敛速度。     

小参数常微分方程守恒型差分格式的一致收敛性

本文考虑自共轭常微分方程奇异摄动边值问题,构造一族带拟合因子的差分格式,给出差分格式解一致收敛于微分方程解的充分条件,由此提出几个具体格式,在条件较弱的情况下,给出较高的一致收敛阶。     

带扩散项四阶抛物方程一类新的交替分组显格式

首先给出逼近带扩散项四阶抛物方程初边值问题一类非对称差分格式,利用该组非对称格式构造了一类新的交替分组显格式算法,并给出了截断误差分析和绝对稳定性结论,最后给出数值实验.     

Unconditional convergence of DIRK schemes applied to dissipative evolution equations

In this paper we prove the convergence of algebraically stable DIRK schemes applied to dissipative evolution equations on Hilbert spaces. The convergence analysis is unconditional as we do not impose any restrictions on the initial value or assume any extra regularity of the solution. The analysis is based on the observation that the schemes are linear combinations of the Yosida approximation, which enables the usage of an abstract approximation result for dissipative maps. The analysis is also extended to the case where the dissipative vector field is perturbed by a locally Lipschitz continuous map. The efficiency and robustness of these schemes are finally illustrated by applying them to a nonlinear diffusion equation.     

一阶双曲型方程组初边值问题的差分格式及其稳定性

本文比较系统地讨论了有关数值求解两个自变量的一阶双曲型方程组初边值问题的某些问题,给出了几种能用于任何类型的初边值问题的差分格式,并在很宽的条件下证明了其中的某些变系数的初边值问题的差分格式对初值和边值是稳定的、差分格式所立出的方程组是良态的.其中的某些格式已用于解决某些复杂的实际问题(应用部分见[16]).     

自共轭常微分方程奇异摄动问题一致收敛的差分格式

本文对自共轭常微分方程奇异摄动问题,构造一族带拟合因子的差分格式,用不同于[1]的方法,通过对格式截断误差的分析,给出差分格式解一致收敛于微分方程解的充分条件;由此提出几个具体的差分格式,在较弱的条件下,给出较高的一致收敛阶,并将它们应用于例子,给出数值结果.     

On the Approximation of Linear Hamiltonian Systems

When we study the oscillation of a physical system near its equilibrium and ignore dissipative effects, we may assume it is a linear Hamiltonian system (H-system), which possesses a special symplectic structure. Thus there arises a question: how to take this structure into account in the approximation of the H-system? This question was first answered by Feng Kang for finite dimensional H-systems.We will in this paper discuss the symplectic difference schemes preserving the symplectic structure and its related properties, with emphasis on the infinite dimensional H-systems.     

Two Novel Classes of Arbitrary High-Order Structure-Preserving Algorithms for Canonical Hamiltonian Systems

In this paper, we systematically construct two classes of structure-preserving schemes with arbitrary order of accuracy for canonical Hamiltonian systems. The one class is the symplectic scheme, which contains two new families of parameterized symplectic schemes that are derived by basing on the generating function method and the symmetric composition method, respectively. Each member in these schemes is symplectic for any fixed parameter. A more general form of generating functions is introduced, which generalizes the three classical generating functions that are widely used to construct symplectic algorithms. The other class is a novel family of energy and quadratic invariants preserving schemes, which is devised by adjusting the parameter in parameterized symplectic schemes to guarantee energy conservation at each time step. The existence of the solutions of these schemes is verified. Numerical experiments demonstrate the theoretical analysis and conservation of the proposed schemes.     

非线性波动方程的弱隐式与显式差分方法

广泛出现于物理、化学、机械动力学、生物、几何学等领域的非线性波动方程已经有很多的研究工作,Sine-Gordon方程和非线性受迫振动方程就是典型的例子.周毓麟教授在[1]中研究了非线性波动方程组     

Efficient local structure‐preserving schemes for the RLW‐type equation

The multisymplectic schemes have been used in numerical simulations for the RLW‐type equation successfully. They well preserve the local geometric property, but not other local conservation laws. In this article, we propose three novel efficient local structure‐preserving schemes for the RLW‐type equation, which preserve the local energy exactly on any time‐space region and can produce richer information of the original problem. The schemes will be mass‐ and energy‐preserving as the equation is imposed on appropriate boundary conditions. Numerical experiments are presented to verify the efficiency and invariant‐preserving property of the schemes. Comparisons with the existing nonconservative schemes are made to show the behavior of the energy affects the behavior of the solution.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1678–1691, 2017     

The Conservation and Convergence of Two Finite Difference Schemes for KdV Equations with Initial and Boundary Value Conditions

Korteweg-de Vries equation is a nonlinear evolutionary partial differential equation that is of third order in space. For the approximation to this equation with the initial and boundary value conditions using the finite difference method, the difficulty is how to construct matched finite difference schemes at all the inner grid points. In this paper, two finite difference schemes are constructed for the problem. The accuracy is second-order in time and first-order in space. The first scheme is a two-level nonlinear implicit finite difference scheme and the second one is a three-level linearized finite difference scheme. The Browder fixed point theorem is used to prove the existence of the nonlinear implicit finite difference scheme. The conservation, boundedness, stability, convergence of these schemes are discussed and analyzed by the energy method together with other techniques. The two-level nonlinear finite difference scheme is proved to be unconditionally convergent and the three-level linearized one is proved to be conditionally convergent. Some numerical examples illustrate the efficiency of the proposed finite difference schemes.     

PARALLEL COMPUTING METHOD OF PURE ALTERNATIVE SEGMENT EXPLICIT-IMPLICIT DIFFERENCE SCHEME FOR NONLINEAR LELAND EQUATION

The research on the numerical solution of the nonlinear Leland equation has important theoretical significance and practical value. To solve nonlinear Leland equation, this paper offers a class of difference schemes with parallel nature which are pure alternative segment explicit-implicit(PASE-I) and implicit-explicit(PASI-E) schemes. It also gives the existence and uniqueness,the stability and the error estimate of numerical solutions for the parallel difference schemes. Theoretical analysis demonstrates that PASE-I and PASI-E schemes have obvious parallelism, unconditionally stability and second-order convergence in both space and time. The numerical experiments verify that the calculation accuracy of PASE-I and PASI-E schemes are better than that of the existing alternating segment Crank-Nicolson scheme, alternating segment explicit-implicit and implicit-explicit schemes. The speedup of PASE-I scheme is 9.89, compared to classical Crank-Nicolson scheme. Thus the schemes given by this paper are high efficient and practical for solving the nonlinear Leland equation.     

A Note on the Stability of the Integral-Differential Equation of the Hyperbolic Type in a Hilbert Space

In this paper, the initial-value problem for integral-differential equation of the hyperbolic type in a Hilbert space H is considered. The unique solvability of this problem is established. The stability estimates for the solution of this problem are obtained. The difference scheme approximately solving this problem is presented. The stability estimates for the solution of this difference scheme are obtained. In applications, the stability estimates for the solutions of the nonlocal boundary problem for one-dimensional integral-differential equation of the hyperbolic type with two dependent limits and of the local boundary problem for multidimensional integral-differential equation of the hyperbolic type with two dependent limits are obtained. The difference schemes for solving these two problems are presented. The stability estimates for the solutions of these difference schemes are obtained.     

对流扩散方程一类改进的特征线修正有限元方法

１引言在地下水污染,地下渗流驱动,核污染,半导体等问题的数值模拟中,均涉及抛物型对流扩散方程（或方程组）的数值求解问题．这些对流扩散型偏微分方程（或方程组）具有共同的特点：对流的影响远大于扩散的影响,即对流占优性,对流占优性给问题的数值求解带来许多困难,因此对流占优问题的有效数值解法一直是计算数学中重要的研究内容．用通常的差分法或有限元法进行数值求解将出现数值振荡．为了克服数值振荡,提出各种迎风方法和修正的特征方法并在这些问题上得到成功的实际应用、８０年代,Ｄｏｕｇｌａｓ和Ｒｕｓｓｅｌｌ［２］等…     

Closed Form Solution and Equivalent Equation Approximation of Linear Advection by a Non Dissipative Second Order Scheme for Step Initial Conditions

We analyse the solution of the linear advection equation on a uniform mesh by a non dissipative second order scheme for discontinuous initial condition. These schemes are known to generate parasitic oscillations in the vicinity of the discontinuity. An approximate way to predict these oscillations is provided by the equivalent equation method. More specifically, we focus on the case of advection of a step function by the leapfrog scheme. Numerical experiments show that the equivalent equation method fails to reproduce the oscillations generated by the scheme far from the discontinuity. Thus, we derive closed form exact and approximate solutions for the scheme that accurately predict these oscillations. We study the relationship between equivalent equation approximation and exact solution for the scheme, to determine its range of validity.