刚性方程组M~.(t)+C~.X(t)+KX(t)=F(t)解的探讨

讨论了刚性常微分方程组(1)的解析解和数值解,给出了解的一般形式和应用该算法的数值例子.     

A CLASS OF r-POINT (r+1)ST-ORDER A-STABLE ONE-BLOCK METHODS WITH DAMPING AT THE INFINITE POINT

1　ＩｎｔｒｏｄｕｃｔｉｏｎＦｏｒｓｏｌｖｉｎｇｓｔｉｆｆｉｎｉｔｉａｌｖａｌｕｅｐｒｏｂｌｅｍｓｆｏｒｓｙｓｔｅｍｓｏｆＯＤＥｓｙ′=ｆ(ｙ) ,ｙ(ｔ0 ) =ｙ0 ,ｔ0 <ｔ≤Ｔ ,ｙ0 ,ｙ∈Ｒｍ,ｆ :Ω Ｒｍ →Ｒｍ (1 .1 )ｍａｎｙｐａｒｔｉｃｕｌａｒｏｎｅ ｂｌｏｃｋｍｅｔｈｏｄｓｏｆｔｈｅｆｏｒｍＹｎ+1= ＡＹｎ+ｈ( Ｂ0 Ｆ(Ｙｎ) + Ｂ1Ｆ(Ｙｎ+1) ) , Ａ =Ａ Ｉｍ, Ｂｉ=Ｂｉ Ｉｍ,Ａ ,Ｂｉ∈Ｒｒ×ｒ,Ｙｎ =(ＹＴｎｒ,… ,ｙＴ(ｎ+1)ｒ- 1) Ｔ,Ｆ(Ｙｎ) =(ｆＴ(ｙｎｒ) ,… ,ｆＴ(ｙ(ｎ+1)ｒ- 1) ) Ｔ,ｙｊ≈ ｙ(ｔｊ) ,…     

New initial-value method for singularly perturbed boundary-value problems

An initial-value method is given for second-order singularly perturbed boundary-value problems with a boundary layer at one endpoint. The idea is to replace the original two-point boundary value problem by two suitable initial-value problems. The method is very easy to use and to implement. Nontrivial text problems are used to show the feasibility of the given method, its versatility, and its performance in solving linear and nonlinear singularly perturbed problems.This work was supported in part by the Consiglio Nazionale delle Ricerche, Contract No. 86.02108.01, and in part by the Ministero della Pubblica Istruzione.     

B-Theory of Runge-Kutta methods for stiff volterra functional differential equations

B-stability andB-convergence theories of Runge-Kutta methods for nonlinear stiff Volterra functional differential equations (VFDEs) are established which provide unified theoretical foundation for the study of Runge-Kutta methods when applied to nonlinear stiff initial value problems (IVPs) in ordinary differential equations (ODEs), delay differential equations (DDEs), integro-differential equations (IDEs) and VFDEs of other type which appear in practice.     

Generalizing Rosenbluth's Algorithm to Include Along‐the‐Chain Intramolecular Energies

We incorporate the Boltzmann factors for inter‐monomer bending energy into the monomer growth direction choice in Rosenbluth's algorithm to model chains of arbitrary nearest‐neighbor rigidity. This allows for the consideration of compact (bent state lower in energy), free (straight and bent state equal in energy), or extended chains (bent state higher). We validate against, and compare to, various other results, showing very good agreement with known results for short chains and demonstrate the ability to model chains up to 500 segments long, far beyond the length at which the normal Rosenbluth method becomes unstable for reasonable nonzero bending energies. This approach is easily generalizable both to other energies determinable during chain growth, for example, polymers composed of more than one type of monomer with differing monomer interaction energies, as well as to other chain production algorithms. © 2019 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2019 , 57, 1684–1691     

解强刚性块线代数方程组的两类L-收敛迭代法

对解强刚性块线代数方程组X=(A(?)J)X φ,本文提出了L-收敛的最佳单参数迭代法(L-OOPI)和L-收敛的多参数迭代直接法(L-MPID),并给出了数值例子．数例表明,对于强刚性块线代数方程组,该二迭代法是有效的．     

A time-decomposition algorithm for a stiff linear two-point boundary-value problem and its application to a nonlinear optimal control problem

An algorithm is proposed to solve a stiff linear two-point boundary-value problem (TPBVP). In a stiff problem, since some particular solutions of the system equation increase and others decrease rapidly as the independent variable changes, the integration of the system equation suffers from numerical errors. In the proposed algorithm, first, the overall interval of integration is divided into several subintervals; then, in each subinterval a sub-TPBVP with arbitrarily chosen boundary values is solved. Second, the exact boundary values which guarantee the continuity of the solution are determined algebraically. Owing to the division of the integration interval, the numerical error is effectively reduced in spite of the stiffness of the system equation. It is also shown that the algorithm is successfully imbedded into an interaction-coordination algorithm for solving a nonlinear optimal control problem.The authors would like to thank Mr. T. Sera and Mr. H. Miyake for their help with the calculations.     

Equilibrium schemes for scalar conservation laws with stiff sources

We consider a simple model case of stiff source terms in hyperbolic conservation laws, namely, the case of scalar conservation laws with a zeroth order source with low regularity. It is well known that a direct treatment of the source term by finite volume schemes gives unsatisfactory results for both the reduced CFL condition and refined meshes required because of the lack of accuracy on equilibrium states. The source term should be taken into account in the upwinding and discretized at the nodes of the grid. In order to solve numerically the problem, we introduce a so-called equilibrium schemes with the properties that (i) the maximum principle holds true; (ii) discrete entropy inequalities are satisfied; (iii) steady state solutions of the problem are maintained. One of the difficulties in studying the convergence is that there are no estimates for this problem. We therefore introduce a kinetic interpretation of upwinding taking into account the source terms. Based on the kinetic formulation we give a new convergence proof that only uses property (ii) in order to ensure desired compactness framework for a family of approximate solutions and that relies on minimal assumptions. The computational efficiency of our equilibrium schemes is demonstrated by numerical tests that show that, in comparison with an usual upwind scheme, the corresponding equilibrium version is far more accurate. Furthermore, numerical computations show that equilibrium schemes enable us to treat efficiently the sources with singularities and oscillating coefficients.

     

Lattice-Boltzmann type relaxation systems and high order relaxation schemes for the incompressible Navier-Stokes equations

A relaxation system based on a Lattice-Boltzmann type discrete velocity model is considered in the low Mach number limit. A third order relaxation scheme is developed working uniformly for all ranges of the mean free path and Mach number. In the incompressible Navier-Stokes limit the scheme reduces to an explicit high order finite difference scheme for the incompressible Navier-Stokes equations based on nonoscillatory upwind discretization. Numerical results and comparisons with other approaches are presented for several test cases in one and two space dimensions.

     

二维刚性边值问题的小波-吉尔解法

本文利用微分方程数值解的离散小波表示，讨论了此类方程在满足一定初始条件和边值条件下，在一个方向上利用小波伽辽金法，另一方向上利用吉尔方法进行求解，提出了一种解二维刚性初，边值问题的小波数值算法，计算结果表明，利用该方法所求得的数值解精度高，而且由小波特有的性质，它特别适用于求解带有奇异摄动的刚性问题。