Numerical analysis and accurate computation of heteroclinic orbits in the case of center manifolds

In earlier papers (see the preceding paper and the references there), Doedel and the author have developed a numerical method for the computation of branches of heteroclinic orbits for a system of autonomous ordinary differential equations in ⁿ in the case that the solution approaches the fixed points exponentially. The idea of the method is to reduce a boundary value problem on the real line to a boundary value problem on a finite interval by using linear approximation of the unstable and stable manifolds. Using the fact that the linearized operator of the problem is Fredholm in Banach spaces with exponential weights, the authors employed the general theory of approximation of nonlinear problems to show that the errors in the approximate solution decay exponentially with the length of the approximating interval. In this paper we extend the analysis in the preceding paper to the case of center manifolds which requires the refinement of the analysis in the preceding paper. The algorithm is applied to a model problem: the DC Josephson Junction. Computations are done using the software package AUTO.     

Bifurcation and coalescence of a plethora of homoclinic orbits for a Hamiltonian system

This is a further study of the set of homoclinic solutions (i.e., nonzero solutions asymptotic to 0 as ¦x¦) of the reversible Hamiltonian systemu ^iv +Pu +u–u ²=0. The present contribution is in three parts. First, rigorously for P –2, it is proved that there is a unique (up to translation) homoclinic solution of the above system, that solution is even, and on the zero-energy surface its orbit coincides with the transverse intersection of the global stable and unstable manifolds. WhenP=–2 the origin is a node on its local stable and unstable manifolds, and whenP(–2,2) it is a focus. Therefore we can infer, rigorously, from the discovery by Devaney of a Smale horseshoe in the dynamics on the zero energy set, there are infinitely many distinct infinite families of homoclinic solutions forP(–2, –2+) for some>0. Buffoni has shown globally that there are infinitely many homoclinic solutions for allP(–2,0], based on a different approach due to Champneys and Toland. Second, numerically, the development of the set of symmetric homoclinic solutions is monitored asP increases fromP=–2. It is observed that two branches extend fromP=–2 toP=+2 where their amplitudes are found to converge to 0 asP 2. All other symmetric solution branches are in the form of closed loops with a turning point betweenP=–2 andP=+2. Numerically it is observed that each such turning point is accompanied by, though not coincident with, the bifurcation of a branch of nonsymmetrical homoclinic orbits, which can, in turn, be followed back toP=–2. Finally, heuristic explanations of the numerically observed phenomena are offered in the language of geometric dynamical systems theory. One idea involves a natural ordering of homoclinic orbits on the stable and unstable manifolds, given by the Horseshoe dynamics, and goes some way to accounting for the observed order (in terms ofP-values) of the occurrence of turning points. The near-coincidence of turning and asymmetric bifurcation points is explained in terms of the nontransversality of the intersection of the stable and unstable manifolds in the zero energy set on the one hand, and the nontransversality of the intersection of the same manifolds with the symmetric section in ⁴ on the other. Some conjectures based on present understanding are recorded.     

Sequences of orbits and the boundaries of the basin of attraction for two double heteroclinic orbits

The boundaries of the basin of attraction are usually assumed to be rather elementary for Hamiltonian systems with autonomous perturbations. In the case of one saddle point, the sequences of orbits before capture are unique for each basin. However, we show that for two saddle points each with double heteroclinic orbits, there is an infinite number of different sequences of nearly homoclinic orbits before capture depending on the four heteroclinic parameters. The probabilities of capture are independent of the capture sequence.     

Piezoelectric control of Hopf bifurcations: A non-linear discrete case study

Linear and non-linear analyses of a piezoelectric controlled non-linear Ziegler column are carried out in this paper. The aim is to evaluate the effects of a linear piezoelectric element on the Hopf bifurcations of the non-linear mechanical system, triggered by the non-conservative follower force. To this end a linear stability analysis, showing the performances of the controller in shifting forward the critical load of the uncontrolled system, is carried out and the role of the electro-mechanical coupling parameter and of the mechanical damping is investigated. The beneficial effects of the controller on the amplitude of the limit cycle occurring in the post-critical field are also discussed.     

Mode jumping in the buckling of struts and plates: a comparative study

A detailed investigation of the phenomenon of mode jumping in compressed struts on stiffening foundations and elastic plates of varying lengths is performed, with emphasis on the effects of altering boundary conditions. The variety of possible modal interactions is presented in a concise form using the parameter space of Arnol'd tongues, borrowed from non-linear dynamical systems theory. For the strut system, a full range of end conditions from simply supported to clamped is examined. For the plate, simply supported and clamped flexural conditions along both long (unloaded) and short (loaded) edges are considered, together with in-plane conditions ranging from free to pull in, to fully restrained. For each system, simply supported end conditions are found to provide protection against early mode jumping in a so-called “safety envelope”, but this is eroded as the end conditions are systematically altered from simply supported to clamped. For the plate system, mode jumping is induced at an earlier stage in the loading process by restricting the long (unloaded) edges against in-plane movement, but is delayed by clamping the same edges against rotation.     

Comparative analysis of chaos control methods: A mechanical system case study

Chaos may be exploited in order to design dynamical systems that may quickly react to some new situation, changing conditions and their response. In this regard, the idea that chaotic behavior may be controlled by small perturbations allows this kind of behavior to be desirable in different applications. This paper presents an overview of chaos control methods classified as follows: OGY methods - include discrete and semi-continuous approaches; multiparameter methods - also include discrete and semi-continuous approaches; and time-delayed feedback methods that are continuous approaches. These methods are employed in order to stabilize some desired UPOs establishing a comparative analysis of all methods. Essentially, a control rule is of concern and each controller needs to follow this rule. Noisy time series is treated establishing a robustness analysis of control methods. The main goal is to present a comparative analysis of the capability of each chaos control method to stabilize a desired UPO.     

A stochastic interrogation method for experimental measurements of global dynamics and basin evolution: Application to a two-well oscillator

An experimental study of local and global bifurcations in a driven two-well magneto-mechanical oscillator is presented. A detailed picture of the local bifurcation structure of the system is obtained using an automated bifurcation data acquisition system. Basins of attractions for the system are obtained using a new experimental technique: an ensemble of initial conditions is generated by switching between stochastic and deterministic excitation. Using this stochastic interrogation method, we observe the evolution of basins of attraction in the nonlinear oscillator as the forcing amplitude is increased, and find evidence for homoclinic bifurcation before the onset of chaos. Since the entire transient is collected for each initial condition, the same data can be used to obtain pictures of the flow of points in phase space. Using Liouville's Theorem, we obtain damping estimates by calculating the contraction of volumes under the action of the Poincaré map, and show that they are in good agreement with the results of more conventional damping estimation methods. Finally, the stochastic interrogation data is used to estimate transition probability matrices for finite partitions of the Poincaré section. Using these matrices, the evolution of probability densities can be studied.     

机电耦联系统失稳振荡的动态分岔研究

应用机电分析动力学的理论建立了交流电机组的机电耦联振动方程组。运用微分动力学系统理论深入分析了交流电机的机电耦联失稳振荡问题。对于该系统出现的高余维分岔问题,通过中心流形定理、多参数稳定性理论和归一化方法得到了原系统的Normal Form形式,并详细讨论了系统的分岔情况以及分岔解的稳定性,并进行了详细的数值计算分析,很好地验证了理论分析结果。取得了交流电机失稳振荡更深入一步的研究成果。     

Error analysis and iterative methods in pseudospectral semi-implicit simulations: an application to compressible convection

Compressible convection is an interesting field for numerical experiments. Rapidly varying small-scale flow structures appear as the Rayleigh number Ra increases, demanding larger spatial resolution under more and more severe Courant stability conditions. Coupling a pseudospectral approximation in space to a semi-implicit scheme in time allows one to increase the size of Δ ts, though at each time step a system of algebraic equations, whose size increases with the spatial resolution, must be solved by means of direct or iterative methods. The former allows one to minimize the consumption of CPU time but leads to unacceptable demand of memory. The efficiency and cost of the latter, on the other hand, depend heavily on the choice of the preconditioning operator and on the allowed error tolerance. In this paper we check the capabilities of iterative-like methods and we achieve the main goal of drastically reducing the memory storage with respect to direct methods, without increasing the CPU time. © 1997 John Wiley & Sons, Ltd.     

Computational Aeroacoustics: An Overview of Computational Challenges and Applications

The objective of this paper is to present an overview of recent advances in computational aeroacoustics (CAA). During the last decade, CAA has developed quite independent of computational fluid dynamics (CFD). There are computational issues that are unique to CAA and are, generally, not considered in CFD. In this paper, these issues are discussed and explained. In CAA, there is a great need to resolve high-frequency short waves with the minimum number of mesh points per wavelength. There is also a special need to minimize numerical dispersion and dissipation associated with wave propagation computation. All these have led to the development of large-stencil high-resolution schemes for CAA. A careful examination of dispersion and dissipation errors due to spatial and temporal discretization is provided. These errors are quantified and analyzed in wave number space through the use of Fourier-Laplace transforms. At this time, some of the original computational challenges to CAA have been resolved satisfactorily. A discussion of how some of these computational issues are resolved is presented. Several important CAA applications with interesting or unusual computational innovations are highlighted. Finally, a few of the most pressing outstanding computational challenges to CAA are elaborated.     

A New dynamic model for study of dislocation pattern formation

Based on the principle given in nonlinear diffusion-reaction dynamics, a new dynamic model for dislocation patterning is proposed by introducing a relaxation time to the relation between dislocation density and dislocation flux. The so-called chemical potential like quantities, which appear in the model can be derived from variation principle for free energy functional of dislocated media, where the free energy density function is expessed in terms of not only the dislocation density itself but also their spatial gradients. The linear stability analysis on the governing equations of a simple dislocation density shows that there exists an intrinsic wave number leading to bifurcation of space structure of dislocation density. At the same time, the numerical results also demonstrate the coexistence and transition between different dislocation patterns. The project supported by the National Natural Science Foundation of China, Grant No.19392300     

A study of the hybrid element methods for stress analysis

A new functional which forms the basis of an improved hybrid element formulation is proposed. The variables for the functional include stresses, strains and displacements, and the displacements and stresses are further decomposed into two parts respectively. The proposed new formulation appears to be particularly suitable for improving conforming models.The relationship between the new hybrid elements and the conventional displacement elements are also explored in this paper.     

埋地管道的横向静力分析耦合模型及其解耦条件

现有埋地管道横向静力分析中的刚性管假设可以实现横向分析与纵向分析的解耦,但往往与实际情况不符.为此,本文根据埋地管道的受力特点,建立了考虑纵向内力影响的横向静力分析耦合模型.基于此耦合模型,结合挠曲线方程与纵向内力的关系,推导出横向分析和纵向分析的解耦条件,为建立解耦模型提供了理论基础.与此同时,一方面为便于实用,另一方面为检验解耦模型的误差,文中针对圆环形管道导出了其横截面剪应力分布规律,给出了埋地管道横向静力分析的实用耦合模型,并结合弹性中心法,建议和提出了耦合模型的实现算法.最后,采用耦合模型对两端固定的埋地给水管道进行了横向内力分析及解耦模型的误差探讨.结果表明,管道两端处纵向剪力对横向内力影响显著,解耦模型误差很大;但随着与管端距离的增加纵向剪力的影响急剧减小,近跨中区域解耦误差几乎可忽略不计.     

A STUDY OF THE CATASTROPHE AND THE CAVITATION FOR A SPHERICAL CAVITY IN HOOKE’S MATERIAL WITH 1/2 POISSON’S RATIO

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

一类有限元结构分析的EBE计算方法

本文利用ＥＢＥ策略￣［１］的基本思想，给出一类有限元结构分析的ＥＢＥ计算方法，即ＥＢＥ共轭梯度法ＥＢＥ－ＣＧ和ＥＢＥ预处理共轭梯度法ＥＢＥ－ＰＣＧ，这类方法避免了传统有限元结构分析中总刚度阵的组集而可大大降低存储量要求。同时它们还特别适合在各种粒度下的多处理机系统上实现。初步数值试验结果表明：这类ＥＢＥ有限元结构分析方法对串行和并行计算都是很有效的。     

Simulation of oscillatory flow in an aortic bifurcation using FVM and FEM: A comparative study of implementation strategies

Cardiovascular diseases are one of the major causes of long‐term morbidity and mortality in human beings. The nearly epidemic increase in prevalence of such diseases poses a serious threat to public health and calls for efficient methods of diagnosis and treatment. Non‐invasive diagnostic procedures such as MRI are often used in this context; however, these are limited in terms of spatial and temporal resolution and do not provide information on time‐dependent pressures and wall shear stresses—key quantities considered to be partially responsible for the formation and development of related pathologies. The present study is concerned with the numerical simulation of oscillatory flow through the abdominal aortic bifurcation. Computational fluid dynamics simulation of oscillatory flow in a branched geometry at high Reynolds numbers poses considerable challenges. The present study reports a detailed comparison of simulations performed with a finite volume and a finite element method, two approaches with significant differences in their discretization strategy, treatment of boundary conditions and other numerical aspects. Both solvers were parallelized, using loop parallelization of the BiCGStab linear solver for the finite volume and domain decomposition based on the Schur complement method for the finite element technique. The experience gained with these two approaches for the solution of flow in a bifurcation forms the focus of this study. Although similar results were obtained for both methods, the computation time required for convergence was found to be significantly smaller for the finite element approach. Copyright © 2010 John Wiley & Sons, Ltd.     

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.     

A spectral projection method for the analysis of autocorrelation functions and projection errors in discrete particle simulation

Discrete particle simulation is a well‐established tool for the simulation of particles and droplets suspended in turbulent flows of academic and industrial applications. The study of some properties such as the preferential concentration of inertial particles in regions of high shear and low vorticity requires the computation of autocorrelation functions. This can be a tedious task as the discrete point particles need to be projected in some manner to obtain the continuous autocorrelation functions. Projection of particle properties on to a computational grid, for instance, the grid of the carrier phase, is furthermore an issue when quantities such as particle concentrations are to be computed or source terms between the carrier phase and the particles are exchanged. The errors committed by commonly used projection methods are often unknown and are difficult to analyse. Grid and sampling size limit the possibilities in terms of precision per computational cost. Here, we present a spectral projection method that is not affected by sampling issues and addresses all of the above issues. The technique is only limited by computational resources and is easy to parallelize. The only visible drawback is the limitation to simple geometries and therefore limited to academic applications. The spectral projection method consists of a discrete Fourier‐transform of the particle locations. The Fourier‐transformed particle number density and momentum fields can then be used to compute the autocorrelation functions and the continuous physical space fields for the evaluation of the projection methods error. The number of Fourier components used to discretize the projector kernel can be chosen such that the corresponding characteristic length scale is as small as needed. This allows to study the phenomena of particle motion, for example, in a region of preferential concentration that may be smaller than the cell size of the carrier phase grid. The precision of the spectral projection method depends, therefore, only on the number of Fourier modes considered. Copyright © 2008 John Wiley & Sons, Ltd.     

SPR误差评估在动静力数值分析中的验证与应用

有限元方法是一种便捷、强大的分析方法,常用于解决工程设计和研究中的各种复杂问题,但作为一种逼近的数值分析方法,其计算结果存在一定的误差,需要利用合理的误差评估方法评价有限元解并为进一步的改进提供依据.因此,将SPR误差评估嵌入弹塑性有限元静力数值分析中,验证在该计算中SPR误差评估的可靠性,并应用于相应的流弹塑性的动力分析中,为进一步优化网格、改进有限元计算结果提供依据.     

A two scaled numerical method for global analysis of high dimensional nonlinear systems

This paper first analyzes the features of two classes of numerical methods for global analysis of nonlinear dynamical systems, which regard state space respectively as continuous and discrete ones. On basis of this understanding it then points out that the previously proposed method of point mapping under cell reference (PMUCR), has laid a frame work for the development of a two scaled numerical method suitable for the global analysis of high dimensional nonlinear systems, which may take the advantages of both classes of single scaled methods but will release the difficulties induced by the disadvantages of them. The basic ideas and main steps of implementation of the two scaled method, namely extended PMUCR, are elaborated. Finally, two examples are presented to demonstrate the capabilities of the proposed method.