定常对流扩散反应方程非均匀网格上高精度紧致差分格式

本文构造了非均匀网格上求解定常对流扩散反应方程的高精度紧致差分格式.我们首先基于非均匀网格上函数的泰勒级数展开,给出了一阶导数和二阶导数的高阶近似表达式;然后将模型方程变形,借助于对流扩散方程高精度紧致格式构造的方法,结合原模型方程,得到定常对流扩散反应方程的高精度紧致差分格式;最后给出的数值算例验证了本文格式高精度和高分辨率的优点.     

一维非定常对流扩散方程非均匀网格上的高精度紧致差分格式

本文在非均匀网格上给出了求解非定常对流扩散方程的一种高精度紧致差分格式,特别适合边界层和大梯度等问题的求解.从稳态对流扩散方程入手,首先,基于非均匀网格上的泰勒级数展开对空间导数项进行离散,然后对时间项采用二阶向后欧拉差分公式,从而得到一维非定常对流扩散方程在非均匀网格上的三层全隐式紧致差分格式.新格式在时间具有二阶精度,空间具有三到四阶精度,并且是无条件稳定的.最后,通过数值实验验证了本文格式的精确性,以及在处理诸如边界层和大梯度问题上的优势.     

不可压缩可混溶驱动问题迎风区域分裂差分方法

结合区域分裂思想,本文给出了一维不可压缩可混溶驱动问题两种非重叠区域分裂迎风差分格式。由于饱和度的计算规模远大于压力方程,因此饱和度方程采用了迎风区域分裂差分法,内边界处和各子区域分别对应显隐格式。在稳定性条件下,给出了 l2 模误差估计,最后给出数值算例验证了理论结果。     

二维变系数抛物型方程高精度恒稳定的改进的Douglas格式

本文将算子方法和粘结系数方法相结合,构造出了求二维变系数抛物型方程的改进的Douglas格式,格式的截断误差阶达O(τ2+h4),并用Fourier分析方法讨论了格式的绝对稳定性.文末给出的数值例子表明了本文理论分析的正确性和所构造格式的有效性及优越性.长期以来,现有关于抛物型方程的差分格式都是对常系数方程而言的,本文较简便的给出了二维情形下变系数抛物型方程的差分解法,此方法完全可以推广到高维情形.     

非线性结构动力方程求解的显式差分格式的特性分析

本文针对软化非线性结构体系动力方程,分析了作者推导出的求解方程的显式差分格式的收敛性及稳定性,分别给出了结构体系处于非线性正刚度及屈服后的负刚度反应阶段时显式差分格式的稳定条件。稳定性分析结果表明,只要结构体系处于初始反应的粘-弹性阶段时显式差分格式满足稳定性条件,则可以保证软化非线性结构体系反应求解的整个过程中格式的计算稳定性。     

二维半线性伪抛物方程差分格式及数值模拟

本文讨论了一类具有较强应用背景的二维半线性伪抛物方程,设计了求解此类方程对应的初边值问题的隐式差分格式。基于离散泛函理论和先验估计证明了相应差分格式解的存在唯一性、收敛性和稳定性。最后针对实例给出了数值模拟结果。     

二维对流扩散方程的二阶精度特征差分格式

针对二维对流扩散方程提出了几类二阶精度特征差分格式,给出了这些格式形成的线性代数方程组可解的充分条件,分析证明了这些格式按离散L^2模是二阶收敛的。最后,具体算例表明这些格式对于对流扩散方程有良好的计算效果。     

Cahn-Hilliard方程的显式差分格式

对Ｃａｈｎ－Ｈｉｌｉａｒｄ方程的周期初值问题构造了显式差分格式，证明了该格式的收敛性和稳定性，给出了数值例子     

求解变系数对流扩散反应方程的指数型高精度紧致差分方法

本文给出了一种数值求解变系数对流扩散反应方程的指数型高精度紧致差分方法.我们首先将模型方程变形,借助常系数对流扩散方程的指数型高精度紧致差分格式,采用残量修正法得到变系数对流扩散反应方程的指数型高精度紧致差分格式;并从理论上分析了当Pelect数很大时,本文格式达到四阶计算精度时网格步长的限制条件;离散得到的代数方程组可采用追赶法直接求解.数值实验结果与理论分析完全吻合,表明了本文格式对于边界层问题或大梯度变化的物理量求解问题具有的高精度和鲁棒性的优点.     

复Ginzburg-Landau方程的数值模拟

本文对复Ginzburg-Landau方程的周期边界问题构造了三个数值格式。其中两个差分格式的精度分别为O(τ2+h2),谱方法的精度为O(τ2+hm),其中m为方程的光滑度。用线性化分析的方法给出了格式的稳定性条件,并给出了数值实验。数值实验表明,四阶紧致格式的计算效果最好,既能达到较高的计算精度又能节省大量的计算时间。     

求解带有不连续波数的二维变系数 Helmholtz方程的一种高精度紧致差分方法

很多实际物理问题都可以由带有不连续波数的变系数 Helmholtz 方程进行数值模拟。Helmholtz 方程的数值方法研究是热点问题之一,具有重要的理论和实际意义。由于波数的不连续性,使用传统的有限差分方法求解带有不连续波数的 Helmholtz 方程时通常无法达到原有差分格式的精度。结合浸入界面方法的思想,对带有不连续波数的二维变系数 Helmholtz 方程构造了一类新的四阶紧致有限差分格式,数值实验验证了新方法的可靠性和有效性。     

A High-Order Difference Scheme for the Generalized Cattaneo Equation

A high-order finite difference scheme for the fractional Cattaneo equation is investigated. The $L_1$ approximation is invoked for the time fractional part, and a compact difference scheme is applied to approximate the second-order space derivative. The stability and convergence rate are discussed in the maximum norm by the energy method. Numerical examples are provided to verify the effectiveness and accuracy of the proposed difference scheme.     

A finite difference scheme for solving a three‐dimensional heat transport equation in a thin film with microscale thickness

Heat transport at the microscale is of vital importance in microtechnology applications. The heat transport equation is different from the traditional heat diffusion equation since a second‐order derivative of temperature with respect to time and a third‐order mixed derivative of temperature with respect to space and time are introduced. In this study, we develop a finite difference scheme with two levels in time for the three‐dimensional heat transport equation. It is shown by the discrete energy method that the scheme is unconditionally stable. The three‐dimensional implicit scheme is then solved by using a preconditioned Richardson iteration, so that only a tridiagonal linear system is solved each iteration. Numerical results show that the solution is accurate. Copyright © 2001 John Wiley & Sons, Ltd.     

An unconditionally stable alternating segment difference scheme of eight points for the dispersive equation

A group of new Saul'yev‐type asymmetric difference schemes to approach the dispersive equation are given here. On the basis of these schemes, an alternating difference scheme with intrinsic parallelism for solving the dispersive equation is constructed. The scheme is unconditionally stable. Numerical experiments show that the method has high accuracy. Copyright © 2006 John Wiley & Sons, Ltd.     

Finite element analysis of non-stationary one-dimensional random diffusion problems

A finite element analysis of a class of non-stationary random diffusion problems is considered. By using the one-dimensional heat equation with random initial condition and random external excitation, the statistical numerical formulation is presented. Two typical numerical examples are given for somewhat simplified problems by which the validity of the finite element scheme is discussed. The results obtained by the finite difference scheme are also shown.     

Weighted residual solutions in the time domain

The Galerkin and subdomain forms of the weighted residual method are used to generate recursive equations in time for the numerical solution of a system of ordinary differential equations. The single-step methods that result from a linear interpolation equation match currently available methods whose stability and oscillation properties are known. A three-level scheme developed by combining two linear elements is shown to be unconditionally unstable. Two of the three schemes obtained using a quadratic interpolation equation and quadratic weighting functions are also shown to be unconditionally unstable. The third scheme is unconditionally stable, but the calculated values for a numerical solution of u? + u = 0, u(0) = 1 are not as accurate as the values obtained using the single-step central difference method.     

A novel accurate and computationally efficient integration approach to viscoplastic constitutive model

A novel, accurate, and computationally efficient integration approach is developed to integrate small strain viscoplastic constitutive equations involving nonlinear coupled first-order ordinary differential equations. The developed integration scheme is achieved by a combination of the implicit backward Euler difference approximation and the implicit asymptotic integration. For the uniaxial loading case, the developed integration scheme produces accurate results irrespective of time steps. For the multiaxial loading case, the accuracy and computational efficiency of the developed integration scheme are better than those of either the implicit backward Euler difference approximation or the implicit asymptotic integration. The simplicity of the developed integration scheme is equivalent to that of the implicit backward Euler difference approximation since it also reduces the solution of integrated constitutive equations to the solution of a single nonlinear equation. The algorithm tangent constitutive matrix derived for the developed integration scheme is consistent with the integration algorithm and preserves the quadratic convergence of the Newton–Raphson method for global iterations.     

Simultaneous inversion for a diffusion coefficient and a spatially dependent source term in the SFADE

This paper deals with an inverse problem of determining a diffusion coefficient and a spatially dependent source term simultaneously in one-dimensional (1-D) space fractional advection–diffusion equation with final observations using the optimal perturbation regularization algorithm. An implicit finite difference scheme for solving the forward problem is set forth, and a fine estimation to the spectrum radius of the coefficient matrix of the difference scheme is given with which unconditional stability and convergence are proved. The simultaneous inversion problem is transformed to a minimization problem, and existence of solution to the minimum problem is proved by continuity of the input–output mapping. The optimal perturbation algorithm is introduced to solve the inverse problem, and numerical inversions are performed with the source function taking on different forms and the diffusion coefficient taking on different values, respectively. The inversion solutions give good approximations to the exact solutions demonstrating that the optimal perturbation algorithm with the Sigmoid-type regularization parameter is efficient for the simultaneous inversion problem in the space fractional diffusion equation.     

On extrapolation boundary conditions for the numerical solution of hyperbolic difference schemes

An analysis is presented for the stability of an initial and boundary value difference scheme arising from the solution of a two-dimensional hyperbolic equation by the Lax-Wendroff method. The positive and negative characteristic speeds of the linear equation are not necessarily equal in magnitude. Supplementary boundary conditions required for the finite difference solution are obtained using extrapolation for a linear combination of the dependent variables. The analysis produces stability intervals for the parameter which describes this linear combination of variables. Numerical experiments are described which support the theoretical analysis.