A Simple Semi-Implicit Scheme for Partial Differential Equations with Obstacle Constraints

We develop a simple and efficient numerical scheme to solve a class of obstacle problems encountered in various applications. Mathematically, obstacle problems are usually formulated using nonlinear partial differential equations (PDE). To construct a computationally efficient scheme, we introduce a time derivative term and convert the PDE into a time-dependent problem. But due to its nonlinearity, the time step is in general chosen to satisfy a very restrictive stability condition. To relax such a time step constraint when solving a time dependent evolution equation, we decompose the nonlinear obstacle constraint in the PDE into a linear part and a nonlinear part and apply the semi-implicit technique. We take the linear part implicitly while treating the nonlinear part explicitly. Our method can be easily applied to solve the fractional obstacle problem and min curvature flow problem. The article will analyze the convergence of our proposed algorithm. Numerical experiments are given to demonstrate the efficiency of our algorithm.     

二维抛物型方程的高精度多重网格解法

提出了数值求解二维抛物型方程的一种新的高精度加权平均紧隐格式，利用Fourier分析方法证明了该格式是无条件稳定的，为了克服传统迭代法在求解隐格式是收敛速度慢的缺陷，利用了多重网格加速技术，大大加快了迭代收敛速度，提高了求解效率，数值实验结果验证了方法的精确性和可靠性。     

Periodic solutions for the Landau-Lifshitz-Gilbert equation

Ferromagnetic materials tend to develop very complex magnetization patterns whose time evolution is modeled by the so-called Landau-Lifshitz-Gilbert equation (LLG). In this paper, we construct time-periodic solutions for LLG in the regime of soft and small ferromagnetic particles which satisfy a certain shape condition. Roughly speaking, it is assumed that the length of the particle is greater than its hight and its width. The approach is based on a perturbation argument and the spectral analysis of the corresponding linearized problem as well as the theory of sectorial operators.     

关于一个二阶非线性中立时滞微分方程的非振荡解的存在性

A new second-order nonlinear neutral delay differential equation r（t） x（t） ＋ P（t）x（t-τ） ＋ cr（t） x（t）-x（t-τ） ＋ F t,x（t-σ1）,x（t-σ2）,...,x（t-σn） = G（t）,t ≥ t0,where τ 0,σ1,σ2,...,σn ≥ 0,P,r ∈ C（[t0,＋∞）,R）,F ∈ C（[t0,＋∞）×Rn,R）,G ∈ C（[t0,＋∞）,R） and c is a constant,is studied in this paper,and some sufficient conditions for existence of nonoscillatory solutions for this equation are established and expatiated through five theorems according to the range of value of function P（t）.Two examples are presented to illustrate that our works are proper generalizations of the other corresponding results.Furthermore,our results omit the restriction of Q1（t） dominating Q2（t）（See condition C in the text）.     

Stability of Semi-implicit Finite Volume Scheme for Level Set Like Equation

We study numerical methods for level set like equations arising in image processing and curve evolution problems. Semi-implicit finite volume-element type schemes are constructed for the general level set like equation(image selective smoothing model) given by Alvarez et al.(Alvarez L, Lions P L, Morel J M. Image selective smoothing and edge detection by nonlinear diffusion II. SIAM J. Numer.Anal., 1992, 29: 845–866). Through the reasonable semi-implicit discretization in time and co-volume method for space approximation, we give finite volume schemes,unconditionally stable in L∞ and W1,2(W1,1) sense in isotropic(anisotropic) diffusion domain.     

二阶中立型微分方程非振动解的存在性

考虑了一类变系数的具有强迫项的二阶中立型微分方程(x(t)+R(t)x(h(t)))″+P(t)x(g_1(t))-Q(t)x(g_2(t))=f(t)非振动解的存在性问题.通过Banach压缩映像原理,分别得到了方程存在满足■|x(t)|>0的非振动解x(t)的充分条件与必要条件,推广了一阶变系数方程的相应结果.     

On the Cauchy problem for the noncompact Landau-Lifshitz-Gilbert equation

Local well-posedness of the Cauchy problem for the noncompact Landau-Lifshitz-Gilbert equation is investigated via the pseudo-stereographic projection. Existence of global solutions is established for small initial data. In the case of one space dimension global existence theorems are proved for large initial data.     

非线性二阶时滞微分方程边值问题的正解

本文讨论了一类二阶时滞微分方程边值问题的正解存在性,在不要求非线性项取值恒为非负的情形下,利用锥上不动点指数的计算得到了该问题的正解.     

一种求解Burgers方程的高精度紧致差分格式

针对Burgers方程,采用余项修正法和欧拉公式,推导了一种新的四层高精度紧致差分隐格式,其截断误差为O(τ~2+τh~2+h~4),即当τ=O(h~2)时,格式空间具有四阶精度;然后通过数值实验验证了格式的精确性和可靠性.     

Unconditionally Optimal Error Analysis of the Second-Order BDF Finite Element Method for the Kuramoto-Tsuzuki Equation

This paper aims to study a second-order semi-implicit BDF finite element scheme for the Kuramoto-Tsuzuki equations in two dimensional and three dimensional spaces. The proposed scheme is stable and the nonlinear term is linearized by the extrapolation technique. Moreover, we prove that the error estimate in $L^2$-norm is unconditionally optimal which means that there has not any restriction on the time step and the mesh size. Finally, numerical results are displayed to illustrate our theoretical analysis.     

三维泊松方程的高精度多重网格解法

利用对称网格点泰勒展开式中各阶导数项明显的对称性,得到了数值求解三维泊松方程的四阶和六阶精度的紧致差分格式,其推导过程简便直接.为了克服传统迭代法在求解高维问题时计算量大、收敛速度慢的缺陷,采用了多重网格加速技术,设计了相应的多重网格算法,求解了三维泊松方程的Dirichlet边值问题.数值实验结果表明,本文所提出的高精度紧致格式达到了期望的精度并且多重网格方法的加速效果是非常显著的.     

解二维抛物型方程的高精度差分格式

本文构造了一个解二维抛物型方程的高精度三层显式差分格式,其稳定性条件为r=△t/△x2=△t/△y2≤1/4,截断误差为O(△t2+△x4).     

时间分数阶扩散方程的数值解法

分数阶微分方程在许多应用科学上比整数阶微分方程更能准确地模拟自然现象.考虑时间分数阶扩散方程,将一阶的时间导数用分数阶导数α(0<α<1)替换,给出了一种计算有效的隐式差分格式,并证明了这个隐式差分格式是无条件稳定和无条件收敛的,最后用数值例子说明差分格式是有效的.     

First-Order and Second-Order,Chaos-Free,Finite Difference Schemes for Fisher Equation

A new class of finite difference schemes is constructed for Fisher partial differential equation i.e. the reaction-diffusion equation with stiff source term: $au(1-u)$. These schemes have the properties that they reduce to high fidelity algorithms in the diffusion-free case namely in which the numerical solutions preserve the properties in the exact solutions for arbitrary time step-size and reaction coefficient α＞0 and all nonphysical spurious solutions including bifurcations and chaos that normally appear in the standard discrete models of Fisher partial differential equation will not occur. The implicit schemes so developed obtain the numerical solutions by solving a single linear algebraic system at each step. The boundness and asymptotic behaviour of numerical solutions obtained by all these schemes are given. The approach constructing the above schemes can be extended to reaction-diffusion equations with other stiff source terms.     

Error Estimates for Semi-Discrete and Fully Discrete Galerkin Finite Element Approximations of the General Linear Second-Order Hyperbolic Equation

In this article, we derive error estimates for the semi-discrete and fully discrete Galerkin approximations of a general linear second-order hyperbolic partial differential equation with general damping (which includes boundary damping). The results can be applied to a variety of cases (e.g. vibrating systems of linked elastic bodies). The results generalize pioneering work of Dupont and complement a recent article by Basson and Van Rensburg.     

高维热传导方程的高精度分支显格式

A high-order accuracy explicit difference scheme for solving 4-dimensional heatconduction equation is constructed.The stability condition is r = △t/△x2 = △t/△y2 = △t/△]z2 = △t/△w2＜3/8,and the truncation error is O(△t2 △x4).     

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

利用有界延拓法,研究了非线性波动方程周期初边值问题的显式差分解的收敛性与稳定性,避免了较难的先验估计,并放宽了非线性项的条件。     

抛物型方程的一个高精度隐格式

利用待定参数法,对一维抛物型方程构造出了一个截断误差为Ｏ（△ｘ＾４＋△ｘ＾４）的隐式差分格式,格式的稳定性条件为ｒ＝ａ△ｔ／△ｘ＾２≤１／√２,可用追赶法求解。     

解四阶抛物型方程高精度紧致差分格式

针对四阶抛物型方程周期初值问题,提出了一个两层隐式差分格式和一个三层隐式差分格式.它们的局部截断误差分别为O((Δt)2+(Δx)4)和O((Δt)2+(Δt)(Δx)2+(Δx)4),其中Δt,Δx分别为时间步长和空间步长.误差分析和数值实验均表明,本文构造的差分格式比经典的Crank-Nicolson格式和Saul’ev构造的差分格式精度更高.从精度及稳定性方面考虑,本文构造的格式也比文[5]的显式格式要好.