Runge–Kutta Methods for Itô Stochastic Differential Equations with Scalar Noise

A general class of stochastic Runge–Kutta methods for Itô stochastic differential equation systems w.r.t. a one-dimensional Wiener process is introduced. The colored rooted tree analysis is applied to derive conditions for the coefficients of the stochastic Runge–Kutta method assuring convergence in the weak sense with a prescribed order. Some coefficients for new stochastic Runge–Kutta schemes of order two are calculated explicitly and a simulation study reveals their good performance.     

Case of a Boltzmann gas leading to the Smoluchowski coagulation equation

A special model of a rarefied hard-sphere gas is considered. The hard-sphere particles undergo absolutely elastic collisions. It is assumed that particles can collide only if their nonzero velocities are orthogonal to each other. The model makes it possible to proceed from the Boltzmann equation to the Smoluchowski coagulation equation, where coagulation means that the kinetic energies of the colliding particles are added. A Monte Carlo scheme for simulation of the phenomenon is described, and the convergence of the simulation algorithm is proved. The convergence of numerical results to exact solutions of the Smoluchowski equation and to finite-difference solutions is tested.     

四阶方程两点边值问题Hermite有限元解的渐近展式与外推

1引言有限元解的渐近展式是提高微分方程数值解精度的重要工具,比如亏量校正和外推就是建立在有限元解的渐近展式的基础之上．许多作者对此进行了大量的研究(见[1]-[4]),特别是文[1],提出了在研究有限元解的渐近展式中十分有用的能量嵌入技巧．本文利用能量嵌入定理得到了四阶方程两点边值问题Hermite有限元解及其二阶平均导数的渐近展式,进一步我们还讨论了它们的Richardson外推公式．考虑四阶方程两点边值问题     

一类非定常对流占优扩散问题差分-流线扩散法的后处理

讨论了一类非定常对流占优扩散方程的差分-流线扩散格式(FDSD),利用插值后处理技术,提高了特殊网格下该FDSD格式在双线性元空间的精度,从而按L∞(L2(Ω) 模达到最优．     

On solving energy-dependent partitioned eigenvalue problem by genetic algorithm: The case of real symmetric Hamiltonian matrices

An energy-dependent partitioning scheme is explored for extracting a small number of eigenvalues of a real symmetric matrix with the help of genetic algorithm. The proposed method is tested with matrices of different sizes (30 × 30 to 1000 × 1000). Comparison is made with Löwdin’s strategy for solving the problem. The relative advantages and disadvantages of the GA-based method are analyzed     

Numerical solution of differential equations using Sinc method based on the interpolation of the highest derivatives

In general, we will use the numerical differentiation when dealing with the differential equations. Thus the differential equations can be transformed into algebraic equations and then we can get the numerical solutions. But as we all have known, the numerical differentiation process is very sensitive to even a small level of errors. In contrast it is expected that on average the numerical integration process is much less sensitive to errors. In this paper, based on the Sinc method we provide a new method using Sinc method incorporated with the double exponential transformation based on the interpolation of the highest derivatives (SIHD) for the differential equations. The error in the approximation of the solution is shown to converge at an exponential rate. The numerical results show that compared with the exiting results, our method is of high accuracy, of good convergence with little computational efforts. It is easy to treat nonhomogeneous mixed boundary condition for our method, which is unlike the traditional Sinc method.     

Multisplitting Iteration Schemes for Solving a Class of Nonlinear Complementarity Problems

We consider several synchronous and asynchronous multisplitting iteration schemes for solving aclass of nonlinear complementarity problems with the system matrix being an H-matrix.We establish theconvergence theorems for the schemes.The numerical experiments show that the schemes are efficient forsolving the class of nonlinear complementarity problems.     

Bending-torsion vibration of a partially submerged cylinder with an arbitrary cross-section

An analytical method is presented to investigate the bending-torsion vibration characteristics of a cylinder with an arbitrary cross-section and partially submerged in water. The compressibility and the free surface waves of the water are considered simultaneously in the analysis. The exact solution of structure–water interaction is obtained mathematically. Firstly, the analytical expression of the velocity potential of the water is derived by using the method of separation of variables. The unknown coefficients in the velocity potential are determined by the longitudinal and circumferential Fourier expansions along the outer surface of the cylinder and are expressed in the form of integral equations including the unknown dynamic bending deflection and torsional angle of the cylinder. Secondly, the force and torque acting on the cylinder per unit length, provided by the water, are obtained by integrating the water dynamic pressure along the circumference of the cylinder. The general solution of bending-torsion vibration of the cylinder under the water dynamic pressure is derived analytically. The integral equations included in the velocity potential of the water can be solved exactly. Finally, the eigenfrequency equation of cylinder–water interaction is obtained by means of the boundary conditions of the cylinder. Some numerical examples for elliptical columns partially submerged in water are provided to show the application of the present method.     

On the convergence rate and optimization of a numerical method with splitting of boundary conditions for the stokes system in a spherical layer in the axisymmetric case: Modification for thick layers

The convergence rate of a fast-converging second-order accurate iterative method with splitting of boundary conditions constructed by the authors for solving an axisymmetric Dirichlet boundary value problem for the Stokes system in a spherical gap is studied numerically. For R/r exceeding about 30, where r and R are the radii of the inner and outer boundary spheres, it is established that the convergence rate of the method is lower (and considerably lower for large R/r) than the convergence rate of its differential version. For this reason, a really simpler, more slowly converging modification of the original method is constructed on the differential level and a finite-element implementation of this modification is built. Numerical experiments have revealed that this modification has the same convergence rate as its differential counterpart for R/r of up to 5 × 10³. When the multigrid method is used to solve the split and auxiliary boundary value problems arising at iterations, the modification is more efficient than the original method starting from R/r ～ 30 and is considerably more efficient for large values of R/r. It is also established that the convergence rates of both methods depend little on the stretching coefficient η of circularly rectangular mesh cells in a range of η that is well sufficient for effective use of the multigrid method for arbitrary values of R/r smaller than ～ 5 × 10³.     

Error estimate for the approximation of nonlinear conservation laws on bounded domains by the finite volume method

In this paper we derive a priori and a posteriori error estimates for cell centered finite volume approximations of nonlinear conservation laws on polygonal bounded domains. Numerical experiments show the applicability of the a posteriori result for the derivation of local adaptive solution strategies.