A second-order splitting combined with orthogonal cubic spline collocation method for the Rosenau equation

A second-order splitting method is applied to a KdV-like Rosenau equation in one space variable. Then an orthogonal cubic spline collocation procedure is employed to approximate the resulting system. This semidiscrete method yields a system of differential algebraic equations (DAEs) of index 1. Error estimates in L² and L^∞ norms have been obtained for the semidiscrete approximations. For the temporal discretization, the time integrator RADAU5 is used for the resulting system. Some numerical experiments have been conducted to validate the theoretical results and to confirm the qualitative behaviors of the Rosenau equation. Finally, orthogonal cubic spline collocation method is directly applied to BBM (Benjamin–Bona–Mahony) and BBMB (Benjamin–Bona–Mahony–Burgers) equations and the well-known decay estimates are demonstrated for the computed solution. © 1998 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 14: 695–716, 1998     

开放式多体系统动力学仿真算法软件研发(Ⅱ)DAEs求解算法对比

研究求解微分-代数方程组(DAEs)的高效率、高精度和高稳定性数值积分方法一直是多体系统动力学领域的热点问题之一。本文将求解结构动力学方程的Bathe数值积分策略应用于DAEs的求解,并基于SiPESC平台开发了开放式多体系统动力学仿真算法软件,综合比较研究了Newmark法、HHT-I3法、Generalizedα方法、Bathe方法和祖冲之类Symplectic方法。通过复摆、刚-柔耦合双摆和对称陀螺三个数值算例研究了算法参数与数值阻尼的关系。数值实验表明,Newmark方法在特定参数下引入的数值阻尼通常不可控,HHT-I3方法、Generalizedα方法和Bathe方法通过选择特定步长和参数可引入可控的数值阻尼,祖冲之类Symplectic方法无数值阻尼。在求解真实高频和低频耦合问题以及高速旋转的陀螺问题时,采用祖冲之类Symplectic方法或者无耗散的Newmar方法能够对系统的高频成分进行准确模拟。     

DAEs and PDEs in elastic multibody systems

Elastic multibody systems arise in the simulation of vehicles, robots, air- and spacecrafts. They feature a mixed structure with differential-algebraic equations (DAEs) governing the gross motion and partial differential equations (PDEs) describing the elastic deformation of particular bodies. We introduce a general modelling framework for this new application field and discuss numerical simulation techniques from several points of view. Due to different time scales, singular perturbation theory and model reduction play an important role. A slider crank mechanism with a 2D FE grid for the elastic connecting rod illustrates the techniques. This revised version was published online in June 2006 with corrections to the Cover Date.     

ON A REGULARIZATION OF INDEX 2 DIFFERENTIAL-ALGEBRAIC EQUATIONS WITH PROPERLY STATED LEADING TERM

In this article, linear regular index 2 DAEs A(t)[D(t)x(t)]′+B(t)x(t)=q(t) are considered. Using a decoupling technique, initial condition and boundary condition are properly formulated. Regular index 1 DAEs are obtained by a regularization method. We study the behavior of the solution of the regularization system via asymptotic expansions. The error analysis between the solutions of the DAEs and its regularization system is given.     

Index-aware model order reduction for differential-algebraic equations

We introduce a model order reduction (MOR) procedure for differential-algebraic equations, which is based on the intrinsic differential equation contained in the starting system and on the remaining algebraic constraints. The decoupling procedure in differential and algebraic part is based on the projector and matrix chain which leads to the definition of tractability index. The differential part can be reduced by using any MOR method, we use Krylov-based projection methods to illustrate our approach. The reduction on the differential part induces a reduction on the algebraic part. In this paper, we present the method for index-1 differential-algebraic equations. We implement numerically this procedure and show numerical evidence of its validity.     

Solving differential-algebraic equations by Taylor series (II): Computing the System Jacobian

The authors have developed a Taylor series method for solving numerically an initial-value problem differential-algebraic equation (DAE) that can be of high index, high order, nonlinear, and fully implicit, BIT, 45 (2005), pp. 561–592. Numerical results have shown that this method is efficient and very accurate. Moreover, it is particularly suitable for problems that are of too high an index for present DAE solvers. This paper develops an effective method for computing a DAE’s System Jacobian, which is needed in the structural analysis of the DAE and computation of Taylor coefficients. Our method involves preprocessing of the DAE and code generation employing automatic differentiation. Theory and algorithms for preprocessing and code generation are presented. An operator-overloading approach to computing the System Jacobian is also discussed. AMS subject classification (2000)  34A09, 65L80, 65L05, 41A58     

约束多体系统动力学分析的改进的离散零空间算法

通过对已有离散零空间矩阵计算方法的改进,构造了改进的离散零空间等效变换公式,该公式可不依赖于特定的积分方法,能简洁、方便的与多种数值积分方法相结合。基于改进公式,提出了改进的离散零空间算法框架,并将该框架与一般变分积分法结合,构造了约束多体系统动力学分析的改进的离散零空间算法。通过曲柄滑块机构的数值实验,验证了改进的离散零空间等效变换公式的正确性,示例了其与数值积分算法的良好结合性,说明了改进算法的可行性和有效性。     

预测校正公式求解指标2微分代数方程

1．前言微分代数方程（EEES）是经常出现于实际问题中的一类方程．其数值求解已成为常微分方程数值求解领域十分活跃的一个方向．目前微分代数方程求解的数值方法主要是nunge－Kutta型方法及BDF方法．Runge－Kutta型方法在网,问中有详细的介绍．Hairer等人据此编制了软件RADAU,而目前使用最广泛的软件还是PetZold等编制的DASSL．DASSL使用的方法为BDF方法,它在微分代数方程中的应用最早可以追述到Gear的开创性工作问．BDF方法一个很大的优点是刚性稳定．然而对于非刚性的微分代数方程,刚性稳定已不是主要考虑的因素．因此…     

The Waveform Relaxation method for systems of differential/algebraic equations

This paper reports efforts towards establishing a parallel numerical algorithm known as Waveform Relaxation (WR) for simulating large systems of differential/algebraic equations. The WR algorithm was established as a relaxation based iterative method for the numerical integration of systems of ODEs over a finite time interval. In the WR approach, the system is broken into subsystems which are solved independently, with each subsystem using the previous iterate waveform as “guesses” about the behavior of the state variables in other subsystems. Waveforms are then exchanged between subsystems, and the subsystems are then resolved repeatedly with this improved information about the other subsystems until convergence is achieved.

In this paper, a WR algorithm is introduced for the simulation of generalized high-index DAE systems. As with ODEs, DAE systems often exhibit a multirate behavior in which the states vary as differing speeds. This can be exploited by partitioning the system into subsystems as in the WR for ODEs. One additional benefit of partitioning the DAE system into subsystems is that some of the resulting subsystems may be of lower index and, therefore, do not suffer from the numerical complications that high-index systems do. These lower index subsystems may therefore be solved by less specialized simulations. This increases the efficiency of the simulation since only a portion of the problem must be solved with specially tailored code. In addition, this paper established solvability requirements and convergence theorems for varying index DAE systems for WR simulation.     

约束多体系统独立广义坐标的数值选取

多体系统的完整约束定义了一个嵌入到欧氏空间的微分流型,多体系统独立的广义坐标选取问题等价于该流型的坐标选取问题。据此,本文提出了一种新的独立广义坐标数值选取理论,可将多体系统的微分-代数混合型运动方程转化为较易求解的纯微分方程,并且不以求解非线性方程组作为必须的手段。