基于精细积分技术的非线性动力学方程的同伦摄动法

将精细积分技术（PIM）和同伦摄动方法（HPM）相结合,给出了一种求解非线性动力学方程的新的渐近数值方法。采用精细积分法求解非线性问题时,需要将非线性项对时间参数按Taylor级数展开,在展开项少时,计算精度对时间步长敏感;随着展开项的增加,计算格式会变得越来越复杂。采用同伦摄动法,则具有相对筒单的计算格式,但计算精度较差,应用范围也限于低维非线性微分方程。将这两种方法相结合得到的新的渐近数值方法则同时具备了两者的优点,既使同伦摄动方法的应用范围推广到高维非线性动力学方程的求解,又使精细积分方法在求解非线性问题时具有较简单的计算格式。数值算例表明,该方法具有较高的数值精度和计算效率。     

Stability of nonlinear time-varying perturbed differential equations

The goal of this paper is twofold. The first part presents a converse Lyapunov theorem for the notion of uniform practical exponential stability of nonlinear differential equations in presence of small perturbation. This class of nonlinear differential equations can be viewed as parametric differential equations. The second part provides the classical perturbation method of seeking an approximate solution as a finite Taylor expansion of the exact solution. The practical asymptotic validity on the approximate is established on infinite-time interval. Finally, we give a numerical example to prove the validity of our methods.     

On the multiple scales perturbation method for difference equations

In the classical multiple scales perturbation method for ordinary difference equations (O Δ Es) as developed in 1977 by Hoppensteadt and Miranker, difference equations (describing the slow dynamics of the problem) are replaced at a certain moment in the perturbation procedure by ordinary differential equations (ODEs). Taking into account the possibly different behavior of the solutions of an O Δ E and of the solutions of a nearby ODE, one cannot always be sure that the constructed approximations by the Hoppensteadt–Miranker method indeed reflect the behavior of the exact solutions of the O Δ Es. For that reason, a version of the multiple scales perturbation method for O Δ Es will be presented and formulated in this paper completely in terms of difference equations. The goal of this paper is not only to present this method, but also to show how this method can be applied to regularly perturbed O Δ Es and to singularly perturbed, linear O Δ Es.     

Taylor expansion method for nonlinear evolution equations

Introduction Thestudyofnonlinearevolutionequationsisafascinatingproblemwhichisattheveryheart oftheunderstandingofmanyimportantproblemsinthenaturalsciences[1,2].Thenonlinear evolutionequationsandtheirnumericalapproximationareveryimportantintheareasof theoreticalmathematicsandcomputationalmathematics.Aninterestingfeatureofthe approximationtheoryofthenonlinearevolutionequationsistheapplicationsofthefunctional analyticmethodstothenumericalapproximationofthenonlinearevolutionequations. Thispaperist…     

Boundary integral equations of unique solutions in elasticity

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｂｏｕｎｄａｒｙｅｌｅｍｅｎｔｍｅｔｈｏｄ（ＢＥＭ）ｐｒｏｖｉｄｅｓａｎａｔｔｒａｃｔｉｖｅａｌｔｅｒｎａｔｉｖｅｆｏｒｔｈｅａｎａｌｙｓｉｓｏｆｅｎｇｉｎｅｅｒｉｎｇｐｒｏｂｌｅｍｓ．Ｉｔｓｍａｉｎａｄｖａｎｔａｇｅｓａｒｅｅｃｏｎｏｍｉｃａｌａｎｄｐａｒｔｉｃｕｌａｒｌｙｃｏｎｖｅｎｉｅｎｔｆｏｒｕｎｂｏｕｎｄｅｄｄｏｍａｉｎａｎｄｓｔｒｅｓｓｃｏｎｃｅｎｔｒａｔｉｏｎｐｒｏｂｌｅｍｓ．Ｔｈｅｂｏｕｎｄａｒｙｉｎｔｅｇｒａｌｅｑｕａｔｉｏｎ（ＢＩＥ）ｉｓｔｈｅ…     

Refined differential equations of deflections in axial symmetrical bending problems of spherical shell and their singular perturbation solutions

This paper deals with the research of accuracy of differential equations of deflections.The basic idea is as follows.Firstly,considering the boundary effect the meridianmidsurface displacement u=0,thus we derive the deflection differential equations;secondly we accurately prove that by use of the deflection differential equations or theoriginal differential equations the same inner forces solutions are obtained;finally,weaccurately prove that considering the boundary effect the meridian surface displacementu=0 is an exact solution.In this paper we give the singular perturbation solution of thedeflection differential equations.Finally we check the equilibrium condition and prove theinner forces solved by perturbation method and the outer load are fully equilibrated.Itshows that perturbation solution is accurate.On the other hand,it shows again that thedeflection differential equation is an exact equation.The features of the new differential equations are as follows:1.The accuracies of the new differentia     

Least-squares finite-element method for shallow-water equations with source terms

Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of source terms in momentum equations, such as the bottom slope and friction of bed, compounds the difficulties further. In this paper, a least-squares finite-element method for the space discretization and θ-method for the time integration is developed for the 2D non-conservative SWE including the source terms. Advantages of the method include： the source terms can be approximated easily with interpolation functions, no upwind scheme is needed, as well as the resulting system equations is symmetric and positive-definite, therefore, can be solved efficiently with the conjugate gradient method. The method is applied to steady and unsteady flows, subcritical and transcritical flow over a bump, 1D and 2D circular dam-break, wave past a circular cylinder, as well as wave past a hump. Computed results show good C-property, conservation property and compare well with exact solutions and other numerical results for flows with weak and mild gradient changes, but lead to inaccurate predictions for flows with strong gradient changes and discontinuities.     

Some methods of training radial basis neural networks in solving the Navier‐Stokes equations

The purpose of this research is to analyze the application of neural networks and specific features of training radial basis functions for solving 2‐dimensional Navier‐Stokes equations. The authors developed an algorithm for solving hydrodynamic equations with representation of their solution by the method of weighted residuals upon the general neural network approximation throughout the entire computational domain. The article deals with testing of the developed algorithm through solving the 2‐dimensional Navier‐Stokes equations. Artificial neural networks are widely used for solving problems of mathematical physics; however, their use for modeling of hydrodynamic problems is very limited. At the same time, the problem of hydrodynamic modeling can be solved through neural network modeling, and our study demonstrates an example of its solution. The choice of neural networks based on radial basis functions is due to the ease of implementation and organization of the training process, the accuracy of the approximations, and smoothness of solutions. Radial basis neural networks in the solution of differential equations in partial derivatives allow obtaining a sufficiently accurate solution with a relatively small size of the neural network model. The authors propose to consider the neural network as an approximation of the unknown solution of the equation. The Gaussian distribution is used as the activation function.     

Nonlinear corrector for Reynolds-averaged Navier-Stokes equations

The scope of this paper is to present a nonlinear error estimation and correction for Navier-Stokes and Reynolds-averaged Navier-Stokes equations. This nonlinear corrector enables better solution or functional output predictions at fixed mesh complexity and can be considered in a mesh adaptation process. After solving the problem at hand, a corrected solution is obtained by solving again the problem with an added source term. This source term is deduced from the evaluation of the residual of the numerical solution interpolated on the h/2 mesh. To avoid the generation of the h/2 mesh (which is prohibitive for realistic applications), the residual at each vertex is computed by local refinement only in the neighborhood of the considered vertex. One of the main feature of this approach is that it automatically takes into account all the properties of the considered numerical method. The numerical examples point out that it successfully improves solution predictions and yields a sharp estimate of the numerical error. Moreover, we demonstrate the superiority of the nonlinear corrector with respect to linear corrector that can be found in the literature.     

Homotopy perturbation method to obtain new solitary solutions with compact support for Boussinesq‐like B(2n, 2n) equations with fully nonlinear dispersion

The aim of this paper is to obtain new solitary solutions with compact support for Boussinesq‐like B(2n, 2n) equations with fully nonlinear dispersion using the homotopy perturbation method (HPM). The special case B(2, 2) is chosen to illustrate the concrete scheme of the HPM in B(2n, 2n) equations. General formulas for the solutions of B(2n, 2n) equations are established. Copyright © 2009 John Wiley & Sons, Ltd.     

On the existence of indeterminate solutions to the equations of motion under non-associated flow

Under certain conditions, an indeterminate solution exists to the equations of motion for dynamic elastic–plastic deformation of materials using constitutive laws based on non-associated flow that suggests that an initially unbounded dynamic perturbation in the stress can develop from a quiescent state on the yield surface. The existence of this indeterminate solution has been alleged to discourage use of non-associated flow rules for both dynamic and quasi-static analysis theoretically. It is shown in this paper that the indeterminate solution that may solve the equations of motion is intrinsically dynamic, and it determinately goes to zero in the quasi-static limit regardless of other indeterminate parameters. Consequently, the existence of this unstable dynamic solution has no impact on stability and use of non-associated flow rules for analysis of the quasi-static problem. More importantly, for dynamic applications, it is also shown that the indeterminate solution solves the equations of motion only if critical restrictions are applied to the constitutive equations such that the effective modulus during loading is constant and the direction of the perturbation is unidirectional over a finite time interval. It is shown that common components of the constitutive laws used in metal forming and deformation analysis are inconsistent with these restrictions. So, these common models can be generalized to include non-associated flow for analysis of the dynamic problem without concern that the solution will become indeterminate.     

Second-order accurate difference method for the singularly perturbed problem of fourth-order ordinary differential equations

In this paper, we construct a uniform second-order difference scheme for a class of boundary value problems of fourth-order ordinary differential equations. Finally, a numerical example is given.     

Finite element solution of the equations governing the flow of electrolyte in charged microporous membranes

Electrical double-layer effects are unimportant in flows through porous media except when the Debye length k^?1 is comparable in magnitude with the pore radius a. Under these conditions the equations governing the flow of electrolyte are those of Stokes, Nernst-Planck and Poisson. These equations are non-linear and require numerical solution. The finite element method provides a useful basis for solution and various algorithms are investigated. The numerical stability and errors of each scheme are analysed together with the development of an appropriate finite element mesh. The electro-osmotic flow of a typical electrolyte (barium chloride) through a uniformly charged cylindrical membrane pore is investigated and the ion fluxes are post-computed from the numerical solutions. The ion flux is shown to be strongly dependent on both zeta potential and pore radius, ka, indicating the effects of overlapping electrical double layers.     

Evolution of a heavy rigid body rotation under the action of unsteady restoring and perturbation torques

The paper develops an approximate solution to the system of Euler’s equations with additional perturbation term for dynamically symmetric rotating rigid body. The perturbed motions of a rigid body, close to Lagrange’s case, under the action of restoring and perturbation torques that are slowly varying in time are investigated. We describe an averaging procedure for slow variables of a rigid body perturbed motion, similar to Lagrange top. Conditions for the possibility of averaging the equations of motion with respect to the nutation phase angle are presented. The averaging technique reduces the system order from 6 to 3 and does not contain fast oscillations. An example of motion of the body using linearly dissipative torques is worked out to demonstrate the use of general equations. The numerical integration of the averaged system of equations is conducted of the body motion. The graphical presentations of the solutions are represented and discussed. A new class of rotations of a dynamically symmetric rigid body about a fixed point with account for a nonstationary perturbation torque, as well as for a restoring torque that slowly varies with time, is studied. The main objective of this paper is to extend the previous results for problem of the dynamic motion of a symmetric rigid body subjected to perturbation and restoring torques. The proposed averaging method is implemented to receive the averaging system of equations of motion. The graphical representations of the solutions are presented and discussed. The attained results are a generalization of our former works where µ and M_i are independent of the slow time τ and M_i depend on the slow time only.

     

含模糊参数结构有限元方程的一种新解法

实际工程问题中，常常会包含一些难以准确描述的系统参数。本文根据模糊集理论，将结构中的不确定参数转化为模糊参数，从而建立起含模糊参数的有限元平衡方程，应用一般模糊线性方程解的基本原理，对模糊有限元平衡方程解的概念及方法进行了讨论，并在此基础上提出了一种改进解的概念以及以摄动法为基础的求解方法。该方法假定当不确定参数相对于其清晰值的分解度不很大时，可以将其在清晰值附近作摄动展开。文中也简要讨论了改进解与现有的模糊有限元方程组各种解之间的关系。相对于现有的其它解法，本文建议的方法更易于与常规的有限元软件结合，用于处理工程实际问题。对算例结果的分析表明，本文方法可较好地符合实际要求。     

Limit cycle analysis of a class of strongly nonlinear oscillation equations

The limit cycle of a class of strongly nonlinear oscillation equations of the form % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaaiqadwhagaWaaiabgUcaRmXvP5wqonvsaeHbbjxAHXgiofMCY92D% aGqbciab-DgaNjab-HcaOiaadwhacqWFPaqkcqWF9aqpcqaH1oqzca% WGMbGaaiikaiaadwhacaGGSaGabmyDayaacaGaaiykaaaa!50B8!\[\ddot u + g(u) = \varepsilon f(u,\dot u)\] is investigated by means of a modified version of the KBM method, where is a positive small parameter. The advantage of our method is its straightforwardness and effectiveness, which is suitable for the above equation, where g(u) need not be restricted to an odd function of u, provided that the reduced equation, corresponding to =0, has a periodic solution. A specific example is presented to demonstrate the validity and accuracy of our 09 method by comparing our results with numerical ones, which are in good agreement with each other even for relatively large .     

On singular perturbation boundary-value problem of coupling type system of convection-diffusion equations

In this paper we consider the singular perturbation boundary-value problem of thefollowing coupling type system of convection-diffusion equationsWe advance two methods:the first one is the initial value solving method,by which theoriginal boundary-value problem is changed into a series of unperturbed initial-valueproblems of the first order ordinary differential equation or system so that an asymptoticexpansion is obtained;the second one is the boundary-value solving method,by which theoriginal problem is changed into a few boundary-value problems having no phenomenon ofboundary-layer so that the exact solution can be obtained and any classical numericalmethods can be used to obtain the numerical solution of consismethods can be used to obtainthe numerical solution of consistant high accuracy with respect to the perturbationparameterε     

变厚度中厚板和中厚壳的大挠度分析

采用摄动有限元法分析了变厚度中厚板和中厚壳的大挠度问题。文中借助虚功原理导出了这类板壳的一般非线性有限元方程，同时利用摄动展开求得了逐级摄动有限元的递推算式。算例表明，摄动有限元法分析变厚度中厚板壳问题同样能获得效率高精度好的结果。     

Development of a local MQ‐DQ‐based stencil adaptive method and its application to solve incompressible Navier–Stokes equations

In this paper, a local stencil adaptive method is presented, which is designed for solving computational fluid dynamics (CFD) problems with curved boundaries accurately. A local multiquadric‐differential quadrature (MQ‐DQ) method is used to discretize the governing equations, taking advantage of its meshless nature. The present method bears the properties of both local MQ‐DQ method and local stencil adaptive method and is thus named the local MQ‐DQ‐based stencil adaptive method. Two test problems with curved boundaries are solved to investigate the performance of this solution‐adaptive method. The numerical results indicate that the proposed method is effective and efficient by combining the advantages of meshless property for complex geometries and local adaptation for accuracy improvement. Copyright © 2007 John Wiley & Sons, Ltd.     

Third‐order‐accurate semi‐implicit Runge–Kutta scheme for incompressible Navier–Stokes equations

A semi‐implicit three‐step Runge–Kutta scheme for the unsteady incompressible Navier–Stokes equations with third‐order accuracy in time is presented. The higher order of accuracy as compared to the existing semi‐implicit Runge–Kutta schemes is achieved due to one additional inversion of the implicit operator I‐τγL, which requires inversion of tridiagonal matrices when using approximate factorization method. No additional solution of the pressure‐Poisson equation or evaluation of Navier–Stokes operator is needed. The scheme is supplied with a local error estimation and time‐step control algorithm. The temporal third‐order accuracy of the scheme is proved analytically and ascertained by analysing both local and global errors in a numerical example. Copyright © 2005 John Wiley & Sons, Ltd.