ARC-length method for differential equations

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｏｒｄｉｎａｒｙａｎｄｐａｒｔｉａｌｄｉｆｅｒｅｎｔｉａｌｅｑｕａｔｉｏｎｓｏｆｃｏｎｔｉｎｕｕｍｐｒｏｂｌｅｍａｒｅｏｆｔｅｎｗｉｔｈｃｅｒｔａｉｎｔｙｐｅｓｏｆｓｉｎｇｕｌａｒｉｔｙａｓｓｔｉｆｐｒｏｐｅｒｔｙ,ｏｒ...     

求解变分型积分方程的一种新型数值方法——有限变分法

将作者提出的多虚拟裂纹扩展法(MVCE法)拓展为求解变分型积分方程问题的一种新型数值方法——有限变分法(FVM)。它的基本思想是,给定有限个(N个)局部变分模式,将所求解的未知量用适当的方法离散化,针对这N个局部变分模式列出N个方程,求解N个未知系数,从而求得未知量。单一未知变量FVM的最终方程组的系数矩阵通常是一个对称的窄带矩阵,对角元是大数,有很好的数值计算性能。用FVM求解了三维I型裂纹前缘的应力强度因子(SIF)分布。利用基于FVM的通用权函数法计算程序,可以高精度和高效率地求解表面力、体积力和温度载荷共同作用情况下三维裂纹前缘SIF的分布及其时间历程。FVM可以被推广到更广泛的领域,是一个求解变分型积分方程问题的普遍适用的新型数值方法。     

基于数值延拓方法的自振荡凝胶周期调控

自振荡凝胶是一类由自振荡反应驱动的高分子形变聚合物,其无需外界刺激就会自发地形成周期性的形变,广泛运用于仿生蠕动机器人的研究。自振荡凝胶的化学-力学耦合关系一直是研究的重点,本文基于Ru(bpy)₃催化IPAAm(N-isopropylacrylamide)高分子凝胶的力学特征及BZ(Belousov-Zhabotinsky)反应的三维Oregonator模型构造自振荡凝胶的化学-力学耦合动态模型,并构建数值延拓算法分析自振荡凝胶周期性随反应参数和力学参数变化情况。通过本文构建的微分-代数方程组数值延拓方法,可以有效求解自振荡凝胶的极限环,为自振荡凝胶周期调控提供设计基础。     

The Laplace transform method for Burgers' equation

The Laplace transform method (LTM) is introduced to solve Burgers' equation. Because of the nonlinear term in Burgers' equation, one cannot directly apply the LTM. Increment linearization technique is introduced to deal with the situation. This is a key idea in this paper. The increment linearization technique is the following: In time level t, we divide the solution u(x, t) into two parts: u(x, t^k) and w(x, t), t^k⩽t⩽t^k+1, and obtain a time‐dependent linear partial differential equation (PDE) for w(x, t). For this PDE, the LTM is applied to eliminate time dependency. The subsequent boundary value problem is solved by rational collocation method on transformed Chebyshev points. To face the well‐known computational challenge represented by the numerical inversion of the Laplace transform, Talbot's method is applied, consisting of numerically integrating the Bromwich integral on a special contour by means of trapezoidal or midpoint rules. Numerical experiments illustrate that the present method is effective and competitive. Copyright © 2009 John Wiley & Sons, Ltd.     

基于精细积分技术的非线性动力学方程的同伦摄动法

将精细积分技术（PIM）和同伦摄动方法（HPM）相结合,给出了一种求解非线性动力学方程的新的渐近数值方法。采用精细积分法求解非线性问题时,需要将非线性项对时间参数按Taylor级数展开,在展开项少时,计算精度对时间步长敏感;随着展开项的增加,计算格式会变得越来越复杂。采用同伦摄动法,则具有相对筒单的计算格式,但计算精度较差,应用范围也限于低维非线性微分方程。将这两种方法相结合得到的新的渐近数值方法则同时具备了两者的优点,既使同伦摄动方法的应用范围推广到高维非线性动力学方程的求解,又使精细积分方法在求解非线性问题时具有较简单的计算格式。数值算例表明,该方法具有较高的数值精度和计算效率。     

The fundamental solution method for incompressible Navier–Stokes equations

A complete boundary integral formulation for incompressible Navier–Stokes equations with time discretization by operator splitting is developed using the fundamental solutions of the Helmholtz operator equation with different order. The numerical results for the lift and the drag hysteresis associated with a NACA0012 aerofoil oscillating in pitch show good agreement with available experimental data. © 1998 John Wiley & Sons, Ltd.     

Sensitivity equation method for the Navier‐Stokes equations applied to uncertainty propagation

This works deals with sensitivity analysis (SA) for the Navier‐Stokes equations. The aim is to provide an estimate of the variance of the velocity field when some of the parameters are uncertain and then to use the variance to compute confidence intervals for the output of the model. First, we introduce the physical model and analyze its stability. The sensitivity equations are derived, and their stability analyzed as well. We propose a finite element‐volume numerical scheme for the state and the sensitivity, which is integrated into the open‐source industrial code TrioCFD. Finally, we present some numerical results: a steady and an unsteady test case for the channel flow problem are investigated. For the steady case, we compare the results to the Monte Carlo method and show how the SA technique succeeds in providing very accurate estimates of the variance. For the unsteady case, a new filtering procedure is proposed to deal with a sensitivity that grows in time. The filtered sensitivity is then used to compute the variance of the output and to provide confidence intervals.     

An efficient three-level explicit time-split scheme for solving two-dimensional unsteady nonlinear coupled Burgers' equations

This work deals with the numerical solutions of two-dimensional viscous coupled Burgers' equations with appropriate initial and boundary conditions using a three-level explicit time-split MacCormack approach. In this technique, the differential operators split the two-dimensional problem into two pieces so that the two-step explicit MacCormack scheme can be easily applied to each subproblem. This reduces the computational cost of the algorithm. For low Reynolds numbers, the proposed method is second-order accurate in time and fourth-order convergent in space, whereas it is second-order convergent in both time and space for high Reynolds numbers problems. This observation shows the utility and efficiency of the considered method compared with a broad range of numerical schemes widely studied in the literature for solving the two-dimensional time-dependent nonlinear coupled Burgers' equations. A large set of numerical examples that confirm the theoretical results are presented and critically discussed.     

Numerical residual elimination method of finite differential equations

In this paper,a new elimination of finite differential equations has been discussed.It applies the numerical direct iteration to obtain the residual equations,in which the number of unknowns has been reduced greatly.The solution process is simple and efficient,and the solution is exact     

精细积分法在电报方程求解中的应用

将精细积分法应用到了二维的电报方程的数值计算之中。实例计算表明，该方法具有简单、计算精度高、无条件稳定、不需要进行复杂、费时的频域一时域转换及卷积积分，直接时域分析，处理非零初始值容易等优点。与传统的ＦＦＴ法及ＮＩＬＴ法相比，其效率更高，功能更强。     

基于渐近均匀化方法的准周期梁结构等效刚度数值实现方法研究

相比周期梁结构,准周期梁结构沿轴向梯度变化,具有更大的设计自由度,能够获得更好的结构性能。由于其非均质性,一般将其均匀化为具有等效性质的均质梁结构,但现有工作很少涉及准周期梁结构等效性质的计算。本文针对由周期梁结构映射而成的准周期梁结构,通过引入雅可比矩阵,基于渐近均匀化方法推导的单胞方程及其等效性质计算列式,并建立了其单胞方程及等效刚度的有限元求解列式。该方法可以处理沿轴向变形的任意微单胞构型,数值算例验证了其正确性和有效性。     

Application of the method of lines to unsteady compressible Euler equations

In the sense of method of lines, numerical solution of the unsteady compressible Euler equations in 1D, 2D and 3D is split into three steps: First, space discretization is performed by the first‐order finite volume method using several approximate Riemann solvers. Second, smoothness and Lipschitz continuity of RHS of the arising system of ordinary dimensional equations (ODEs) is analysed and its solvability is discussed. Finally, the system of ODEs is integrated in time by means of implicit and explicit higher‐order adaptive schemes offered by ODE packages ODEPACK and DDASPK, by a backward Euler scheme based on the linearization of the RHS and by higher‐order explicit Runge–Kutta methods. Time integrators are compared from several points of view, their applicability to various types of problems is discussed, and 1D, 2D and 3D numerical examples are presented. Copyright © 2003 John Wiley & Sons, Ltd.     

用拟压缩性方法求解不可压Euler方程

用拟压缩性方法和Ｊａｍｅｓｏｎ的有限体积算法求解了二维和三维定常可可压Ｅｕｌｅｒ方程。分别采用显、隐式时间离散推进求解;分析了人工粘性的阶数对定常解收敛性的影响,应用该方法计算了单个翼型和翼身组合体的低速绕流,结果与实验吻合较好。     

A mixture differential quadrature method for solving two-dimensional incompressible Navier-Stokes equations

Ｔｈｅｄｉｆｆｅｒｅｎｔｉａｌｑｕａｄｒａｔｕｒｅｍｅｔｈｏｄ（ＤＱＭ）ｐｒｏｐｏｓｅｄｂｙＲ．Ｂｅｌｌｍａｎ［１,２］ｈａｓｂｅｅｎｓｕｃｃｅｓｓｆｕｌｌｙｅｍｐｌｏｙｅｄｉｎｎｕｍｅｒｉｃａｌｃｏｍｐｕｔａｔｉｏｎｓｏｆｐｒｏｂｌｅｍｓｉｎｅｎｇｉｎｅｅｒｉｎｇａｎｄｐｈｙｓｉｃａｌｓｃｉｅｎｃｅ．ＢｅｃａｕｓｅｔｈｅｉｎｆｏｒｍａｔｉｏｎｏｎａｌｌｇｒｉｄｐｏｉｎｔｓｉｓｕｓｅｄｔｏｆｉｔｔｈｅｄｅｒｉｖａｔｉｖｅｓａｔｇｒｉｄｐｏｉｎｔｓｉｎｔｈｅＤＱＭ,ｉｔｉｓｅｎｏｕｇｈｔｏｏｂｔａ…     

弹塑性接触问题的非光滑非线性方程组方法

将求解三维弹性摩擦接触问题的非光滑非线性方程组方法推广到弹塑性(Mises材料)情形，提出了两种应用方法：一种是将非光滑非线性方程组方法和求解弹塑性问题常用的Newton—Raphson迭代方法结合起来；另一种是将问题写成统一的非光滑非线性方程组，直接求解。数值算例验证了两种方法的有效性，并进行了结果比较。     

On the discontinuous Galerkin method for the numerical solution of the Navier–Stokes equations

The paper deals with the use of the discontinuous Galerkin finite element method (DGFEM) for the numerical solution of viscous compressible flows. We start with a scalar convection–diffusion equation and present a discretization with the aid of the non‐symmetric variant of DGFEM with interior and boundary penalty terms. We also mention some theoretical results. Then we extend the scheme to the system of the Navier–Stokes equations and discuss the treatment of stabilization terms. Several numerical examples are presented. Copyright © 2004 John Wiley & Sons, Ltd.     

Solution of shallow-water equations using least-squares finite-element method

A least-squares finite-element method （LSFEM） for the non-conservative shallow-water equations is presented. The model is capable of handling complex topography, steady and unsteady flows, subcritical and supercritical flows, and flows with smooth and sharp gradient changes. Advantages of the model include： （1） sources terms, such as the bottom slope, surface stresses and bed frictions, can be treated easily without any special treatment; （2） upwind scheme is no needed; （3） a single approximating space can be used for all variables, and its choice of approximating space is not subject to the Ladyzhenskaya-Babuska-Brezzi （LBB） condition; and （4） the resulting system of equations is symmetric and positive-definite （SPD） which can be solved efficiently with the preconditioned conjugate gradient method. The model is verified with flow over a bump, tide induced flow, and dam-break. Computed results are compared with analytic solutions or other numerical results, and show the model is conservative and accurate. The model is then used to simulate flow past a circular cylinder. Important flow charac-teristics, such as variation of water surface around the cylinder and vortex shedding behind the cylinder are investigated. Computed results compare well with experiment data and other numerical results.     

扩散方程单内点精细积分法与差分法比较研究

一维扩散方程初值问题可以用全域或子域精细积分求解。子域积分可以采用不同数量的内点，单内点是其最简单的情况。当单内点精细积分中的传递函数即指数函数用其泰勒展开式的一阶近似来替代时，精细积分转化为差分方程。本文研究了这一对应关系。各种常见差分格式均找到了对应的单点精细积分格式，并在单点精细积分一般公式中得到了统一表达形式