THE MULTI-SYMPLECTIC ALGORITHM FOR “GOOD” BOUSSINESQ EQUATION

The multi-symplectic formulations of the “Good” Boussinesq equation were considered. For the multi-symplectic formulation, a new fifteen-point difference scheme which is equivalent to the multi-symplectic Preissman integrator was derived. The numerical experiments show that the multi-symplectic scheme have excellent long-time numerical behavior. Foundation items: the Foundation for Key Laboratory of Scientific/Engineering Computing Institute of Computational Mathematics and Scientific/Engineering Computing, Chinese Academy of Sciences; the Natural Science Foundation of Huaqiao University. Biography: ZENG Wen-ping (1940-), Professor (E-mail: qmz@1sec.cc.ac.cn)     

Boundary element method for buckling eigenvalue problem and its convergence analysis

The conditions for determining solution of buckling eigenvalue problem are discussed. The corresponding system of integral equations with constraint conditions and boundary variational equations with Lagrange multiplier are established. The theorems on the existence and uniqueness of the solution for these problems are given. The corresponding boundary element method is constructed and the error estimation for the approximation solution is obtained. Finally the numerical example is given. Foundation item: the National Natural Science Foundation Pre-research Project (T4107015) Biography: Ding Rui (1969-)     

MIXED FINITE ELEMENT METHODS FOR THE SHALLOWWATER EQUATIONS INCLUDING CURRENT AND SILTSEDIMENTATION （Ⅱ）——THE DISCRETE-TIMECASE ALONG CHARACTERISTICS

The mixed finite element(MFE) methods for a shallow water equation systemconsisting of water dynamics equations, silt transport equation, and the equation of bottomtopography change were derived. A fully discrete MFE scheme for the discrete-time alongcharacteristics is presented and error estimates are established. The existence andconvergence of MFE solution of the discrete current velocity, elevation of the bottomtopography, thickness of fluid column, and mass rate of sediment is demonstrated.     

Investigation on gradient-dependent nonlocal constitutive models for elasto-plasticity coupled with damage

IntroductionIntheframeworkofconventionalplasticity ,materialinstabilityisoneoftheprincipalfactorsthatresultinginthestrainlocalizationphenomenon .Byusingtheterminology‘homogenized’ ,itisreferredtothefactthatinitialflawsandboundaryconditionsnecessarilyinduceanon_homogeneousstressstateinaspecimenduringtesting .Inparticularintheprocessofprogressivefailure ,theflawsandlocalstressconcentrationwillcausestronglyinhomogeneousdeformationofthespecimen[1,2 ].Asthedeformation_inducedfracture/damagephenom…     

IMPULSIVE CONTROL FOR THE STABILIZATION AND SYNCHRONIZATION OF LUR＇ E SYSTEMS

An impulsive control scheme of the Lur‘e system and several theorems on stability of impulsive control systems was presented, these theorems were then used to find the conditions under which the Lur‘e system can be stabilized by using impulsive control with varying impulsive intervals. The parameters of Lur‘e system and impulsive control law are given, a theory of impulsive synchronization of two Lur‘e system is also presented. A numerical example is used to verify the theoretical result.     

MIXED FINITE ELEMENT METHODS FOR THE SHALLOW WATER EQUATIONS INCLUDING CURRENT AND SILT SEDIMENTATION （Ⅰ）——THE CONTINUOUS-TIME CASE

An initial-boundary value problem for shallow equation system consisting of water dynamics equations, silt transport equation, the equation of bottom topography change, and of some boundary and initial conditions is studied, the existence of its generalized solution and semidiscrete mixed finite element (MFE) solution was discussed, and the error estimates of the semidiscrete MFE solution was derived. The error estimates are optimal.     

Difference scheme for two-phase flow

A numerical method for two-phase flow with hydrodynamics behavior was considered. The nonconservative hyperbolic governing equations proposed by Saurel and Gallout were adopted. Dissipative effects were neglected but they could be included in the model without major difficulties. Based on the opinion proposed by Abgrall that “a two phase system, uniform in velocity and pressure at t = 0 will be uniform on the same variable during its temporal evolution“, a simple accurate and fully Eulerian numerical method was presented for the simulation of multiphase compressible flows in hydrodynamic regime. The numerical method relies on Godunov-typescheme, with HLLC and Lax-Friedrichs type approximate Riemann solvers for the resolution of conservation equations, and nonconservative equation. Speed relaxation and pressure relaxation processes were introduced to account for the interaction between the phases. Test problem was presented in one space dimension which illustrated that our scheme is accurate, stable and oscillation free.     

The high accuracy explicit difference scheme for solving parabolic equations 3-dimension

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

An automatic constraint violation stabilization method for differential/ algebraic equations of motion in multibody system dynamics

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

Computation of super-convergent nodal stresses of timoshenko beam elements by EEP method

The newly proposed element energy projection (EEP) method has been applied to the computation of super-convergent nodal stresses of Timoshenko beam elements. Generalformul as based on element projection theorem were derived and illustrative numerical examples using two typical elements were given. Both the analysis and examples show that EEP method also works very well for the problems with vector function solutions. The EEP method gives super-convergent nodal stresses, which are well comparable to the nodal displacements in terms of both convergence rate and error magnitude. And in addition, it can overcome the “ shear locking“ difficulty for stresses even when the displacements are badly affected. This research paves the way for application of the EEP method to general onedimensional systems of ordinary differential equations.     

MODAL SYNTHESIS METHOD FOR NORM COMPUTATION OF H∞ DECENTRALIZED CONTROL SYSTEMS （II）

When using H∞ techniques to design decentralized controllers for large systems, the whole system is divided into subsystems, which are analysed using H∞ control theorybefore being recombined. An analogy was established with substructural analysis instructural mechanics, in which H∞ decentralized control theory corresponds to substructuralmodal synthesis theory so that the optimal H∞ norm of the whole system corresponds to thefundamental vibration frequency of the whole structure. Hence, modal synthesismethodology and the extended Wittrick-Williams algorithm were transplanted from structuralmechanics to compute the optimal H∞ norm of the control system. The orthogonality and theexpansion theorem of eigenfunctions of the subsystems H~ control are presented in part (I)of the paper. The modal synthesis method for computation of the optimal H∞ norm ofdecentralized control systems and numerical examples are presented in part (Ⅱ).     

Dynamical Behavior of Viscoelastic Cylindrical Shells Under Axial Pressures

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

Ishikawa iterative process for constructing solutions ofm-accretive operator equations

The convergence of theIshikawa iteration sequences with errors for constructing solutions of m-accretive operator equations is characterized. Moreover, the error estimates of approximate solutions for locally Lipschitzian and m-accretive operator equations are established. Foundation items: the National Natural Science Foundation of China (19801023); the Teaching and Research Award Fund for Outstanding Young Teachers in Higher Education Institutions of MOE, China Biography: ZENG Liu-chuan (1965-), Professor     

ON THE CONVERGENCE PROBLEMS OF ISHIKAWA AND MANN ITERATIVE PROCESSES WITH ERROR FOR Φ-PSEUDO CONTRACTIVE TYPE MAPPINGS

The purpose of this paper is to introduce the concept of Φ-pseudo contractive type mapping and to study the convergence problem of Ishikawa and Mann iterative processes with error for this kind of mappings. The results presented in this paper improve and extend many authors’recent results.     

NONLINEAR COMPLEX DYNAMIC PHENOMENA OF THE PERTURBED METALLIC BAR CONSIDERING DISSIPATING EFFECT

IntroductionTotheweaklydamped ,periodicallyforcedsine_Gordonequation ,A .R .Bishop[1～ 3]analyzeditssolutionunderperiodicboundaryconditionandconcludedthatitssolutionwouldshowdifferentspatialstructuresandlong_timeasymptoticstatesalongwiththevariationofpara…     

A CLASS OF ALTERNATING GROUP METHODOF BURGERS＇ EQUATION

Some new Saul‘ yev type asymmetric difference schemes for Burgers’ equation isgiven, by the use of the schemes, a kind of alternating group four points method for solvingnonlinear Burgers‘ equation is constructed here. The basic idea of the method is that thegrid points on the same time level is divided into a number of groups, the differenceequations of each group can be solved independently, hence the method with intrinsicparallelism can be used directly on parallel computer. The method is unconditionally stableby analysis of linearization procedure. The numerical experiments show that the method has good stability and accuracy.     

MODAL SYNTHESIS METHOD FOR NORM COMPUTATION OF H∞ DECENTRALIZED CONTROL SYSTEMS （ I ）

When using H∞ techniques to design decentralized controllers for large systems,the whole system is divided into subsystems, which are analysed using H∞ control theorybefore being recombined. An analogy was established with substructural analysis instructural mechanics, in which H∞ decentralized control theory corresponds to substructuralmodal synthesis theory so that the optimal H∞ norm of the whole system corresponds to thefundamental vibration frequency of the whole structure. Hence, modal synthesismethodology and the extended Wittrick-Williams algorithm were transplanted from structuralmechanics to compute the optimal H∞ norm of the control system. The orthogonality and theexpansion theorem of eigenfunctions of the subsystems H∞ control are presented in part(I) of the paper. The modal synthesis method for computation of the optimal H∞ norm ofdecentralized control systems and numerical examples are presented in part (Ⅱ).     

Lower bound limit analysis of three-dimensional elastoplastic structures by boundary element method

IntroductionThelimitanalysisofstructuresisoneofthemostpracticalandusefulbranchesinplasticity .Ithasimportantapplicationbackgroundforproblemssuchasthedeterminationofloadcarryingcapacityandplasticformingofmetal.Thepurposeofthelimitanalysisofstructuresistoprovidereliabletheoreticalbasesforengineeringdesignandsafetyassessment.Asasimplifiedmethodforelastoplasticproblems,limitanalysisneednotrequirethehistoryofloadandcancomputethelimitloadsdirectlyinsteadofelastoplasticincrementalcomputationwhichisus…     

Construction of third-order WNND scheme and its application in complex flow

IntroductionWiththedevelopmentofaeronauticsandaerospacetechnology ,moreandmorerequirementsarearisingforCFD (computationalfluiddynamics) .Oneoftheproblemsistodevelophigherorderaccuracyschemes.Forexample ,whenapplyingLES (largeeddysimulation)orDNS(directnumericalsimulation)methodtosimulatingturbulenceproblem ,theschemesneedthirdorderaccuracyormoreinspace .Anotherquestionistheinfluenceofgrid’sscaletotopologicalstructureofflowfield .Inordertosimulatecomplicatedflowswithseparationorturbulencec…     

On the formal solution of initial value problem of Navier-Stokes equation

Ｉｎｒｅｆｅｒｅｎｃｅ [1 ] ,ｗｅｄｉｓｃｕｓｓｅｄｔｈｅｉｎｓｔａｂｉｌｉｔｙｏｆＮａｖｉｅｒ＿ＳｔｏｋｅｓｅｑｕａｔｉｏｎｏｎＲ3×Ｒ ,ａｎｄ ,ｉｎｔｈｅｉｌｌｕｓｔｒａｔｉｏｎｏｆｉｔ,ｗｅｐｒｅｓｅｎｔｅｄｓｏｍｅｅｘａｍｐｌｅｓｉｎｔｈｅｆｏｒｍｏｆｅｘａｃｔｓｏｌｕｔｉｏｎ .Ｔｈｅｂａｓｉｓｏｆｔｈｅｃｏｎｓｔｒｕｃｔｉｏｎｏｆｔｈｅｍｉｓｔｈｅ“ｅｑｕａｔｉｏｎｓｅｃｏｎｄａｉｒｅ”ｏｆｔｈｅｅｑｕａｔｉｏｎ .Ｉｎｆａｃｔ,ｔｈｅｓｅｅｘａｃｔｓｏｌｕｔｉｏｎｓａｒｅｓｏｍｅｌｉｍｉｔ…