Efficient high-order immersed interface methods for heat equations with interfaces

An efficient high-order immersed interface method （IIM） is proposed to solve two-dimensional （2D） heat problems with fixed interfaces on Cartesian grids, which has the fourth-order accuracy in the maximum norm in both time and space directions. The space variable is discretized by a high-order compact （HOC） difference scheme with correction terms added at the irregular points. The time derivative is integrated by a Crank-Nicolson and alternative direction implicit （ADI） scheme. In this case, the time accuracy is just second-order. The Richardson extrapolation method is used to improve the time accuracy to fourth-order. The numerical results confirm the convergence order and the efficiency of the method.     

STUDIES OF MELNIKOV METHOD AND TRANSVERSAL HOMOCLINIC ORBITS IN THE CIRCULAR PLANAR RESTRICTED THREE－BODY PROBLEM

ＳＴＵＤＩＥＳ　ＯＦ　ＭＥＬＮＩＫＯＶ　ＭＥＴＨＯＤ　ＡＮＤ　ＴＲＡＮＳＶＥＲＳＡＬ　ＨＯＭＯＣＬＩＮＩＣ　ＯＲＢＩＴＳ　ＩＮ　ＴＨＥ　ＣＩＲＣＵＬＡＲ　ＰＬＡＮＡＲ　ＲＥＳＴＲＩＣＴＥＤ　ＴＨＲＥＥ－ＢＯＤＹ　ＰＲＯＢＬＥＭＺｈｕＲｕｚｅｎｇ（朱如...     

Viscoelastic behaviour of partly hydrolysed polyacrylamide/chromium(III) gels

This paper reports an experimental investigation of the crosslinking process of high molecular weight partly hydrolysed polyacrylamides (HPAAm) in aqueous brine solution by trivalent chromium ions (Cr(III)). Crosslinking took place in the presence of a retardant agent (sodium citrate). First, sol-gel phase diagrams (in the polymer and concentration space) were established using the tilting tube method. Then, for a fixed composition, the gelation process was monitored systematically using dynamic viscoelastic measurements, varying the main parameters (pH, time, temperature and retardant concentration). Network formation proceeds rather slowly and an equilibrium state was not reached within 12 h. The gel was formed only at pHs between 5 and 9 and thus two gel points (i.e. at two pHs) were determined with the Winter-Chambon method. This is in agreement with the chemistry of aqueous chromium and of acrylic acid groups along the polymer backbone. Kinetics of network formation depends strongly on retardant concentration. Temperature plays an important role: network formation proceeds much faster at high temperature, in agreement with chemical kinetics.     

Geometric Implications of Linearizability

Linearizability for hyperbolic fixed points is shown to be equivalent to the existence of certain invariant foliations. The analogue for partial linearizability of nonhyperbolic fixed points is considered.     

LARGE RANGE ANALYSIS FOR NONLINEAR DYNAMIC SYSTEMS——ELEMENT MAPPING METHOD

This paper presents a new method for global analysis of nonlinear system. By means of transforming the nonlinear dynamic problems into point mapping forms which are single-valued and continuous, the state space can be regularly divided into a certain number of finitely small triangle elements on which the non-linear mapping can be approximately substituted by the linear mapping given by definition. Hence, the large range distributed problem of the mapping fixed points will be simplified as a process for solving a set of linear equations. Still further, the exact position of the fixed points can be found by the iterative technique. It is convenient to judge the stability of fixed points and the shrinkage zone in the state space by using the deformation matrix of linear mapping. In this paper, the attractive kernel for the stationary fixed points is defined, which makes great advantage for describing the attractive domains of the fixed points. The new method is more convenient and effective than the cell mapping methodl^[1]. And an example for two-dimensional mapping is given.     

LARGE RANGE ANALYSIS FOR NONLINEAR DYNAMIC SYSTEMS——ELEMENT MAPPING METHOD

This paper presents a new method for global analysis of nonlinear system. By means of transforming the nonlinear dynamic problems into point mapping forms which are single-valued and continuous, the state space can be regularly divided into a certain number of finitely small triangle elements on which the non-linear mapping can be approximately substituted by the linear mapping given by definition. Hence, the large range distributed problem of the mapping fixed points will be simplified as a process for solving a set of linear equations. Still further, the exact position of the fixed points can be found by the iterative technique. It is convenient to judge the stability of fixed points and the shrinkage zone in the state space by using the deformation matrix of linear mapping. In this paper, the attractive kernel for the stationary fixed points is defined, which makes great advantage for describing the attractive domains of the fixed points. The new method is more convenient and effective than the cell     

Non-linear in-plane analysis and buckling of pinned–fixed shallow arches subjected to a central concentrated load

This paper is concerned with an analytical study of the non-linear elastic in-plane behaviour and buckling of pinned–fixed shallow circular arches that are subjected to a central concentrated radial load. Because the boundary conditions provided by the pinned support and fixed support of a pinned–fixed arch are quite different from those of a pinned–pinned or a fixed–fixed arch, the non-linear behaviour of a pinned–fixed arch is more complicated than that of its pinned–pinned or fixed–fixed counterpart. Analytical solutions for the non-linear equilibrium path for shallow pinned–fixed circular arches are derived. The non-linear equilibrium path for a pinned–fixed arch may have one or three unstable equilibrium paths and may include two or four limit points. This is different from pinned–pinned and fixed–fixed arches that have only two limit points. The number of limit points in the non-linear equilibrium path of a pinned–fixed arch depends on the slenderness and the included angle of the arch. The switches in terms of an arch geometry parameter, which is introduced in the paper, are derived for distinguishing between arches with two limit points and those with four limit points and for distinguishing between a pinned–fixed arch and a beam curved in-elevation. It is also shown that a pinned–fixed arch under a central concentrated load can buckle in a limit point mode, but cannot buckle in a bifurcation mode. This contrasts with the buckling behaviour of pinned–pinned or fixed–fixed arches under a central concentrated load, which may buckle both in a bifurcation mode and in a limit point mode. An analytical solution for the limit point buckling load of shallow pinned–fixed circular arches is also derived. Comparisons with finite element results show that the analytical solutions can accurately predict the non-linear buckling and postbuckling behaviour of shallow pinned–fixed arches. Although the solutions are derived for shallow pinned–fixed arches, comparisons with the finite element results demonstrate that they can also provide reasonable predictions for the buckling load of deep pinned–fixed arches under a central concentrated load.     

Bifurcation and chaos in discrete-time BVP oscillator

In this paper, we investigate the discrete-time Bohöffer-Van der Pol (BVP) oscillator obtained by Euler method. We provide the sufficient conditions of existence, asymptotic stability of the fixed points, then give theoretical analysis for local bifurcations of the fixed points, and derive the conditions under which the local bifurcations such as pitchfork, saddle-node, flip and Hopf occur at the fixed points. Furthermore, we prove that the fixed point eventually evolves into a snap-back repeller which generates chaotic behavior in the sense of Marotto's chaos when certain conditions are satisfied. Finally, several numerical simulations are provided to demonstrate the theoretical results of the previous and to show the new complex dynamical behaviors of the system.     

Dynamics of two vibro-impact nonlinear energy sinks in parallel under periodic and transient excitations

A linear oscillator (LO) coupled with two vibro-impact (VI) nonlinear energy sinks (NES) in parallel is studied under periodic and transient excitations, respectively. The objective is to study response regimes and to compare their efficiency of vibration control. Through the analytical study with multiple scales method, two slow invariant manifolds (SIM) are obtained for two VI NES, and different SIM that result from different clearances analytically supports the principle of separate activation. In addition, fixed points are calculated and their positions are applied to judge response regimes. Transient responses and modulated responses can be further explained. By this way, all analysis is around the most efficient response regime. Then, numerical results demonstrate two typical responses and validate the effectiveness of analytical prediction. Finally, basic response regimes are experimentally observed and analyzed, and they can well explain the complicated variation of responses and their corresponding efficiency, not only for periodic excitations with a fixed frequency or a range of frequency, but also for transient excitation. Generally, vibration control is more effective when VI NES is activated with two impacts per cycle, whatever the types of excitation and the combinations of clearances. This observation is also well reflected by the separate activation of two VI NES with two different clearances, but at different levels of displacement amplitude of LO.     

Momentum and heat transfer of a special case of the unsteady stagnation-point flow

This paper investigates the unsteady stagnation-point flow and heat transfer over a moving plate with mass transfer,which is also an exact solution to the unsteady Navier-Stokes(NS)equations.The boundary layer energy equation is solved with the closed form solutions for prescribed wall temperature and prescribed wall heat flux conditions.The wall temperature and heat flux have power dependence on both time and spatial distance.The solution domain,the velocity distribution,the flow field,and the temperature distribution in the fluids are studied for different controlling parameters.These parameters include the Prandtl number,the mass transfer parameter at the wall,the wall moving parameter,the time power index,and the spatial power index.It is found that two solution branches exist for certain combinations of the controlling parameters for the flow and heat transfer problems.The heat transfer solutions are given by the confluent hypergeometric function of the first kind,which can be simplified into the incomplete gamma functions for special conditions.The wall heat flux and temperature profiles show very complicated variation behaviors.The wall heat flux can have multiple poles under certain given controlling parameters,and the temperature can have significant oscillations with overshoot and negative values in the boundary layers.The relationship between the number of poles in the wall heat flux and the number of zero-crossing points is identified.The difference in the results of the prescribed wall temperature case and the prescribed wall heat flux case is analyzed.Results given in this paper provide a rare closed form analytical solution to the entire unsteady NS equations,which can be used as a benchmark problem for numerical code validation.     

一类碰撞振动系统在内伊马克沙克-音叉分岔点附近的局部两参数动力学

考虑一类具有对称性的三自由度碰撞振动系统.系统的庞加莱映射在一定条件下存在对称不动点,对应于系统的对称周期运动.根据对称性导出庞加莱映射P是另外一个隐式虚拟映射Q的二次迭代.推导了庞加莱映射对称不动点的解析表达式.根据映射不动点的稳定性及分岔理论,映射P的对称不动点发生内伊马克沙克-音叉(Neimark--Saker-pitchfork)分岔对应于映射Q发生内伊马克沙克-倍化(Neimark--Sakerflip)分岔.利用隐式虚拟映射Q,通过对范式作两参数开折分析,研究了映射P的对称不动点在内伊马克沙克-音叉分岔点附近的局部动力学行为.碰撞振动系统在这个余维二分岔点附近的局部动力学行为可能表现为投影后的庞加莱截面上的单一对称不动点、一对共轭不动点、单一对称拟周期吸引子以及一对共轭拟周期吸引子.数值模拟得到了内伊马克沙克-音叉分岔点附近的各种可能情况.内伊马克沙克-分岔和音叉分岔互相作用可能产生新的结果:对称不动点虽然首先分岔为两个共轭不动点,但是这两个共轭不动点是不稳定的,最终收敛到同一个对称拟周期吸引子.     

非线性电报方程组的3个双周期正解问题

This paper studies existence of at least three positive doubly periodic solutions of a coupled nonlinear telegraph system with doubly periodic boundary conditions. First, by using the Green function and maximum principle, existence of solutions of a nonlinear telegraph system is equivalent to existence of fixed points of an operator. By imposing growth conditions on the nonlinearities, existence of at least three fixed points in cone is obtained by using the Leggett-Williams fixed point theorem to cones in ordered Banach spaces. In other words, there exist at least three positive doubly periodic solutions of nonlinear telegraph system.     

Topological Fixed Point Theorems Do Not Hold for Random Dynamical Systems

The aim of this paper is to demonstrate that topological fixed point theorems have no canonical generalization to the case of random dynamical systems. This is done by using tools from algebraic ergodic theory. We give a criterion for the existence of invariant probability measures for group valued cocycles. With that, examples of continuous random dynamical systems on a compact interval without random invariant points, which are an appropriate generalization of fixed points, are constructed.     

The orientation dynamics of small prolate and oblate spheroids in linear shear flows

The present work deals with the stable orientation of oblate and prolate spheroids in general steady linear flows and with the mode of convergence to these stable orientations. The orientation dynamics is governed by the Jeffery equation. The stable orientations are either fixed points or limit cycles in the orientation space. The type of stable orientation depends on whether the eigenvalues of the linear part of Jeffery equation are real or complex. We define prolate and oblate spheroids to be equivalent if the aspect ratio of one is the reciprocal of the other. We show that, in a given flow, equivalent oblate and prolate spheroids possess the same number of fixed points and limit cycles of which only one is stable. If they possess only fixed points, then their corresponding stable fixed points are orthogonal. If they possess one fixed point and one limit cycle each, then the stable fixed point of one is orthogonal to the plane of the limit cycle of the other. The rate of convergence to these attractors is important to consideration of the orientations in time-space varying flow fields. We show that non-normal growth (NNG) of the distance to these attractors may delay the convergence by several characteristic shear time scales. We derive conditions for occurrence of NNG and explicit expressions for the maximal duration of the growth. We consider a specific case of which the vorticity is a stable orientation of prolate spheroids. We analyze the conditions that imply monotonic or, conversely, non-monotonic convergence to this orientation due to NNG. We thereby find the corresponding conditions for convergence of the equivalent oblate spheroids to their attractors, normal to the vorticity. We show that the convergence is monotonic if the vorticity is parallel to the strain tensor’s largest eigenvector, but that NNG occurs if the vorticity is parallel to the strain tensor’s intermediate eigenvector. The NNG duration decreases with increasing vorticity-strain ratio and with the strain intermediate eigenvalue approaching the largest eigenvalue.     

Existence and stability of common fixed points for systems of mappings

In this paper,we establish some common fixed point theorems and stability theorems of the sets of common fixed points for the systems of set-valued and single-valued nonlinear contractive type mappings in a finite Cartesian product of metric spaces.     

Viscosity method for hierarchical fixed point and variational inequalities with applications

A viscosity method for a hierarchical fixed point solving variational inequality problems is presented. The method is used to solve variational inequalities, where the involved mappings are non-expansive. Solutions are sought in the set of the fixed points of another non-expansive mapping. As applications, we use the results to study problems of the monotone variational inequality, the convex programming, the hierarchical minimization, and the quadratic minimization over fixed point sets.     

Viscosity method for hierarchical fixed point and variational inequalities with applications

A viscosity method for a hierarchical fixed point solving variational inequality problems is presented. The method is used to solve variational inequalities, where the involved mappings are non-expansive. Solutions are sought in the set of the fixed points of another non-expansive mapping. As applications, we use the results to study problems of the monotone variational inequality, the convex programming, the hierarchical minimization, and the quadratic minimization over fixed point sets.     

On the hyperchaotic complex Lü system

The aim of this paper is to introduce the new hyperchaotic complex Lü system. This system has complex nonlinear behavior which is studied and investigated in this work. Numerically the range of parameter values of the system at which hyperchaotic attractors exist is calculated. This new system has a whole circle of equilibria and three isolated fixed points, while the real counterpart has only three isolated ones. The stability analysis of the trivial fixed point is studied. Its dynamics is more rich in the sense that our system exhibits both chaotic and hyperchaotic attractors, as well as periodic and quasi-periodic solutions and solutions that approach fixed points. The nonlinear control method based on Lyapunov function is used to synchronize the hyperchaotic attractors. The control of these attractors is studied. Different forms of hyperchaotic complex Lü systems are constructed using the state feedback controller and complex periodic forcing.     

Interative approximation of fixed points for almost asymptotically nonexpansive type mappings in Banach spaces

IntroductionThroughoutthispaperweassumethatEisarealBanachspace ,E isthedualspaceofE ,DisanonemptysubsetofEandJ:E →2 E isthenormalizeddualitymappingdefinedbyJ(x) =f∈E :〈x ,f〉=‖x‖·‖f‖,‖f‖=‖x‖,　　 x∈E .　　Definition 1　LetT :D →Dbeamapping .1 )Tissaidtobeasymptoticallynonexpansive[1],ifthereexistsasequence kn [1 ,∞)withlimn→∞kn =1suchthat　　　　　‖Tnx-Tny‖≤kn‖x-y‖forall　　x ,y∈D ,n≥0 ;(1 )　　2 )Tissaidtobeofasymptoticallynonexpansivetype[2 ],if　　　　　lims…     

Bifurcation analysis of an electro-statically actuated micro-beam in the presence of centrifugal forces

In this paper, the influence of centrifugal forces on the stability of an electro-statically actuated clamped–clamped micro-beam has been investigated. The non-dimensional governing static and dynamic equations have been linearized using the step by step linearization method (SSLM), then, a Galerkin-based reduced order model has been used to solve the linearized equations. For constant value of a bias DC voltage and different values of angular velocity the equilibrium points of the corresponding autonomous system including stable center points, unstable saddle points and singular points have been obtained using the equivalent mass-spring model. Subsequently the bifurcation diagram has been depicted using the obtained fixed point. The static pull-in voltage value for different values of angular velocity and the static pull-in angular velocity for different values of bias voltage have been calculated. The obtained results are validated using results of previous studies and a good agreement has been observed. The effect of the centrifugal force on the fixed points has been studied using the phase portraits of the system for different initial conditions. Moreover, the effects of centrifugal forces on the dynamic pull-in behavior have been investigated using time histories and phase portraits for different angular velocities.