关于USSOR迭代法收敛的充要条件和最佳松弛因子

本文在线性方程组系数矩阵A为相容次序矩阵及A的Jacobi迭代矩阵的特征值μj均为实数的条件下,得出了USSOR迭代法收敛的充分必要性定理,并给出了USSOR迭代矩阵之谱半径ρ(ω,ω)的表达式及ρ(ω,ω)的最佳松弛因子。

Nonlinear dynamics of flexible rotor system supported on fixed-tilting pad combination journal bearing

Based on the Reynolds equation with Reynolds boundary conditions, the Castelli method was employed to solve the Reynolds equation for oil lubrication upon bearings. By doing so, a profile of nonlinear oil film force of single-pad journal bearings is established. According to the structure of combination journal bearings, nonlinear oil film force of combination journal bearing is obtained by retrieval, interpolation and assembly techniques. As for symmetrical flexible Jeffcott rotor systems supported by combination journal bearings, the nonlinear motions of the center of the rotor are calculated by the self-adaptive Runge-Kutta method and Poincar6 mapping with different rotational speeds. The numerical results show that the system performance is slightly better when the pivot ratio changes from 0.5 to 0.6, and reveals nonlinear phenomena of periodic, period-doubing, quasi-periodic motion, etc.

A High-Quality Preconditioning Technique for Multi-Length-Scale Symmetric Positive Definite Linear Systems

We study preconditioning techniques used in conjunction with the conjugate gradient method for solving multi-length-scale symmetric positive definite linear systems originating from the quantum Monte Carlo simulation of electron interaction of correlated materials. Existing preconditioning techniques are not designed to be adaptive to varying numerical properties of the multi-length-scale systems. In this paper, we propose a hybrid incomplete Cholesky （HIC） preconditioner and demonstrate its adaptivity to the multi-length-scale systems. In addition, we propose an extension of the compressed sparse column with row access （CSCR） sparse matrix storage format to efficiently accommodate the data access pattem to compute the HIC preconditioner. We show that for moderately correlated materials, the HIC preconditioner achieves the optimal linear scaling of the simulation. The development of a linear-scaling preconditioner for strongly correlated materials remains an open topic.

Positive Solutions of Boundary Value Problem for a Coupled System of Nonlinear Third-order Differential Equations

By using cone expansion-compression theorem in this paper, we study boundary value problems for a coupled system of nonlinear third-order differential equation. Some sufficient conditions are obtained which guarantee the boundary value problems for a coupled system of nonlinear third-order differential equation has at least one positive solution. Some examples are given to verify our results.

Exponential mean-square stability of the improved split-step theta methods for non-autonomous stochastic differential equations

We consider the mean-square stability of the so-called improved split-step theta method for stochastic differential equations. First, we study the mean-square stability of the method for linear test equations with real parameters. When 0 ≥ 3/2, the improved split-step theta methods can reproduce the mean-square stability of the linear test equations for any step sizes h 〉 0. Then, under a coupled condition on the drift and diffusion coefficients, we consider exponential mean-square stability of the method for nonlinear non-autonomous stochastic differential equations. Finally, the obtained results are supported by numerical experiments.

Spherical Scattered Data Quasi-interpolation by Gaussian Radial Basis Function

Since the spherical Gaussian radial function is strictly positive definite, the authors use the linear combinations of translations of the Gaussian kernel to interpolate the scattered data on spheres in this article. Seeing that target functions axe usually outside the native spaces, and that one has to solve a large scaled system of linear equations to obtain combinatorial coefficients of interpolant functions, the authors first probe into some problems about interpolation with Gaussian radial functions. Then they construct quasi- interpolation operators by Gaussian radial function, and get the degrees of approximation. Moreover, they show the error relations between quasi-interpolation and interpolation when they have the same basis functions. Finally, the authors discuss the construction and approximation of the quasi-interpolant with a local support function.

Finite element method for viscoelastic medium with damage and the application to structural analysis of solid rocket motor grain

This paper studies the damage-viscoelastic behavior of composite solid propellants of solid rocket motors(SRM).Based on viscoelastic theories and strain equivalent hypothesis in damage mechanics,a three-dimensional(3-D)nonlinear viscoelastic constitutive model incorporating with damage is developed.The resulting viscoelastic constitutive equations are numerically discretized by integration algorithm,and a stress-updating method is presented by solving nonlinear equations according to the Newton-Raphson method.A material subroutine of stress-updating is made up and embedded into commercial code of Abaqus.The material subroutine is validated through typical examples.Our results indicate that the finite element results are in good agreement with the analytical ones and have high accuracy,and the suggested method and designed subroutine are efficient and can be further applied to damage-coupling structural analysis of practical SRM grain.

Finite particle method for kinematically indeterminate bar assemblies

This study presents a structural analysis algorithm called the finite particle method （FPM） for kinematically indeterminate bar assemblies. Different from the traditional analysis method, FPM is based on the combination of the vector mechanics and numerical calculations. It models the analyzed domain composed of finite particles. Newton＇s second law is adopted to describe the motions of all particles. A convected material flame and explicit time integration for the solution procedure is also adopted in this method. By using the FPM, there is no need to solve any nonlinear equations, to calculate the stiffness matrix or equilibrium matrix, which is very helpful in the analysis of kinematically indeterminate structures. The basic formulations for the space bar are derived, following its solution procedures for bar assemblies. Three numerical examples are analyzed using the FPM. Results obtained from both the straight pretension cable and the suspension cable assembly show that the FPM can produce a more accurate analysis result. The motion simulation of the four-bar space assembly demonstrates the capability of this method in the analysis ofkinematically indeterminate structures.

The polynomial-like iterative equation for PM functions

The polynomial-like iterative equation is an important form of functional equations, in which iterates of the unknown function are linked in a linear combination. Most of known results were given for the given function to be monotone. We discuss this equation for continuous solutions in the case that the given function is a PM(piecewise monotone) function, a special class of non-monotonic functions. Using extension method, we give a general construction of solutions for the polynomial-like iterative equation.

Compatible Lie Bialgebras

A compatible Lie algebra is a pair of Lie algebras such that any linear combination of the two Lie brackets is a Lie bracket. We construct a bialgebra theory of compatible Lie Mgebras as an analogue of a piiLie bialgebra. They can also be regarded as a ＂compatible version＂ of Lie bialgebras, that is, a pair of Lie biaJgebras such that any linear combination of the two Lie bialgebras is still a Lie bialgebra. Many properties of compatible Lie bialgebras as the ＂compatible version＂ of the corresponding properties of Lie biaJgebras are presented. In particular, there is a coboundary compatible Lie bialgebra theory with a construction from the classical Yang-Baxter equation in compatible Lie algebras as a combination of two classical Yang-Baxter equations in lAe algebras. FUrthermore, a notion of compatible pre-Lie Mgebra is introduced with an interpretation of its close relation with the classical Yang-Baxter equation in compatible Lie a/gebras which leads to a construction of the solutions of the latter. As a byproduct, the compatible Lie bialgebras fit into the framework to construct non-constant solutions of the classical Yang-Baxter equation given by Golubchik and Sokolov.

Nonlinearly Coupled KdV Equations Describing the Interaction of Equatorial and Midlatitude Rossby Waves

Amplitude equations governing the nonlinear resonant interaction of equatorial baroclinic and barotropic Rossby waves were derived by Majda and Biello and used as a model for long range interactions （teleconnections） between the tropical and midlatitude troposphere. An overview of that derivation is nonlinear wave theory, but not in atmospheric presented and geared to readers versed in sciences. In the course of the derivation, two other sets of asymptotic equations are presented： the long equatorial wave equations and the weakly nonlinear, long equatorial wave equations. A linear transformation recasts the amplitude equations as nonlinear and linearly coupled KdV equations governing the amplitude of two types of modes, each of which consists of a coupled tropical/midlatitude flow. In the limit of Rossby waves with equal dispersion, the transformed amplitude equations become two KdV equations coupled only through nonlinear fluxes. Four numerical integrations are presented which show （i） the interaction of two solitons, one from either mode, （ii） and （iii） the interaction of a soliton in the presence of different mean wind shears, and （iv） the interaction of two solitons mediated by the presence of a mean wind shear.