Numerical simulations of periodic and chaotic responses in a stable duffing system

A stable Duffing system is examined by numerical simulations in order to obtain a better understanding of the behavior of periodic and chaotic responses to sinusoidal excitations. It is found that beside the multiplicity of responses, there is a duality for both periodic and chaotic responses. Period doubling does exist and this process may originate from different basic responses even with the same forcing frequency. The evolution of chaos is shown by a sequence of Poincaré maps. Finally a possible pattern for transition to chaos is suggested.     

Some phenomena for weak spring Duffing equation in subharmonic and ultrasubharmonic region

In this paper the situation in which weak spring Duffing equation gets into the chaos on account of small perturbation is discussed with Melnikov-Holmes’ method, and some phenomena in which the different subharmonic and the ultrasubharmonic coexist with the chaos are discovered.     

Numerical simulations of shock wave reflection phenomena in non-stationary flows using regularized smoothed particle hydrodynamics (RSPH)

In this paper we wish to demonstrate to what extent the numerical method regularized smoothed particle hydrodynamics (RSPH) is capable of modelling shocks and shock reflection patterns in a satisfactory manner. The use of SPH based methods to model shock wave problems has been relatively sparse, both due to historical reasons, as the method was originally developed for studies of astrophysical gas dynamics, but also due to the fact that boundary treatment in Lagrangian methods may be a difficult task. The boundary conditions have therefore been given special attention in this paper. Results presented for one quasi-stationary and three non-stationary flow tests reveal a high degree of similarity, when compared to published numerical and experimental data. The difference is found to be below 5, in the case where experimental data was found tabulated. The transition from regular reflection (RR) to Mach reflection (MR) and the opposite transition from MR to RR are studied. The results are found to be in close agreement with the results obtained from various empirical and semi-empirical formulas published in the literature. A convergence test shows a convergence rate slightly steeper than linear, comparable to what is found for other numerical methods when shocks are involved.     

Results of supercomputer simulations of meteorological mesoscale phenomena

A mesoscale model, which is suitable to simulate the distribution of meteorological variables with high resolution in space and time in a region of typically 50 km ×50 km is described. The model's applicability in complex terrain is demonstrated by several examples, namely local weather prediction, calculation of complex wind fields and dispersion, and local climate study.     

Transient phenomena in thixotropic systems

A nonlinear rheological model which accounts for the time-dependent elastic, viscous and yielding phenomena is developed in order to describe the flow behavior of thixotropic materials which exhibit yield stress. A key feature of the formulation is a smooth transition from an ‘elastically’ dominated response to a ‘viscous’ response without a discontinuity in the stress–strain curve. The model is phenomenological and is based on the kinetic processes responsible for structural changes within the thixotropic material. As such, it can predict thixotropic effects, such as stress overshoot during start-up of a steady shear flow and stress relaxation after cessation of flow. Thus this model extends a previously proposed viscoplastic model [J. Rheol. 34 (1991) 647] to include thixotropy.An analysis and comparison to experimental data involving oscillatory shear flow are provided to evaluate the accuracy of the model and to estimate the model parameters in a prototype concentrated suspension. The experiments were conducted using a series of concentrated suspensions of silicon particles and silicon carbide whiskers in polyethylene. The data obtained with this experimental system indicated much better agreement between the theory and experiments that obtained in earlier work.     

Numerical simulations of bouncing jets

Bouncing jets are fascinating phenomenon occurring under certain conditions when a jet impinges on a free surface. This effect is observed when the fluid is Newtonian and the jet falls in a bath undergoing a solid motion. It occurs also for non‐Newtonian fluids when the jets fall in a vessel at rest containing the same fluid. We investigate numerically the impact of the experimental setting and the rheological properties of the fluid on the onset of the bouncing phenomenon. Our investigations show that the occurrence of a thin lubricating layer of air separating the jet and the rest of the liquid is a key factor for the bouncing of the jet to happen. The numerical technique that is used consists of a projection method for the Navier–Stokes system coupled with a level set formulation for the representation of the interface. The space approximation is carried out with adaptive finite elements. Adaptive refinement is shown to be very important to capture the thin layer of air that is responsible for the bouncing. Copyright © 2015 John Wiley & Sons, Ltd.     

Bifurcation phenomena in non-smooth dynamical systems

The aim of the paper is to give an overview of bifurcation phenomena which are typical for non-smooth dynamical systems. A small number of well-chosen examples of various kinds of non-smooth systems will be presented, followed by a discussion of the bifurcation phenomena in hand and a brief introduction to the mathematical tools which have been developed to study these phenomena. The bifurcations of equilibria in two planar non-smooth continuous systems are analysed by using a generalised Jacobian matrix. A mechanical example of a non-autonomous Filippov system, belonging to the class of differential inclusions, is studied and shows a number of remarkable discontinuous bifurcations of periodic solutions. A generalisation of the Floquet theory is introduced which explains bifurcation phenomena in differential inclusions. Lastly, the dynamics of the Woodpecker Toy is analysed with a one-dimensional Poincaré map method. The dynamics is greatly influenced by simultaneous impacts which cause discontinuous bifurcations.     

Unbounded sequences of stable limit cycles in the delayed Duffing equation: an exact analysis

Nonlinear Dynamics - The delayed Duffing equation $${\ddot{x}}(t)+x(t-T)+x^3(t)=0$$ is shown to possess an infinite and unbounded sequence of rapidly oscillating, asymptotically stable periodic...     

Gas-dynamic colliders: Numerical simulations

A collision of supersonic flows of gas mixtures with disparate molecular weights, which are limited in their cross-sectional size, in vacuum leads to formation of a cloud with an elevated concentration and elevated temperature of the heavy gas. Under certain conditions, the governing factor is the collision of molecules of the heavy gas being compressed at the center of the collision of the flows. The generator of such a flow can be called a collider. Results of studying the flows in jet-type, cylindrical, and mixed two-stage colliders are described. The main attention is paid to separation of gases in terms of energy and composition. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 3, pp. 142–151, May–June, 2007.     

Critical phenomena and waves in gas-liquid systems

The results of the hydraulic studies of gas-liquid media, wave processes in two-phase media and critical phenomena are described. Some methodological foundations to describe these media and methods to obtain the basic similarity criteria for the hydraulics and gas-dynamics of bubble suspensions are discussed. A detailed consideration is given for the phase transition processes on interfaces and the interface stability. A relation has been revealed between the wave and critical phenomena in two-phase systems.Nomenclature a thermal diffusivity - Ar Archimedes number - B gas constant - C heat capacity - C _p heat capacity at constant pressure - C _v heat capacity at constant volume - c ₀ acoustic velocity in the mixture - c _l acoustic velocity in the liquid - C _f flow resistance coefficient - G mass rate of flow - g gravitational acceleration - L latent heat of evaporation - l initial perturbation width - M Mach number - Nu Nusselt number - P pressure - Pr Prandtl number - R bubble radius - (3P ₀/R ₀ ² _f )^–1 bubble resonance frequency square - T temperature - U medium motion velocity - W heavy phase velocity - W light phase velocity - We Weber number - heat release coefficient - dispersion coefficient - void fraction - adiabatic index - film thickness - dimensionless film thickness - kinematic viscosity coefficient - dynamical viscosity coefficient - dissipation coefficient in the mixture - dispersion parameter - _f liquid phase density - light phase density - heat conductivity - surface tension - frequency, ₀ ² =3P ₀/ _f R ₀ ²     

Interior layer phenomena of semilinear systems

In this paper,we study the interior layer phenomena of singular perturbation boundaryvalue problems for semilinear systems:εy” = f(t,y.ε)(a0 is a small parameter,y.f. A and B are n-dimensional vector functions.Thisvector boundary value peoblem does not appear to have been studied,although the scalarboundary problem has been treated extensively.Under appopriate assumptions we obtainexistence of solution as in the scalar problem and the estimate of this solution in terms ofappropriate inequalities as well.     

Impact phenomena of rotor-casing dynamical systems

Rubbing and impacting between a rotor and adjacent motion-constraining structures is a serious malfunction in rotating machinery. A shaver rotor-casing system with clearance and mass imbalance is modelled with two second-order ordinary differential equations and inelastic impact conditions. The dynamics is investigated analytically, as well as by numerical simulation. A Lyapunov exponent technique is developed to characterize the topologically different behavior as the parameters are varied. The dry friction coefficient and the eccentricity of the rotor imbalance were chosen to be the two variable parameters, the effect of which on the system dynamics is illustrated through phase plots, bifurcation diagrams, as well as Poincaré maps. The results demonstrate the existence of both rubbing and impacting behavior. Depending on values of the parameters, rubbing motion in both the clockwise and counter-clockwise directions may occur. Within the impact regime, the impact behavior could be periodic, quasi-periodic or chaotic, as confirmed by the calculation of Lyapunov exponents.     

Numerical simulations of turbulent shear flows

We simulate three-dimensional mixing layers, isotropic turbulence, and rotating turbulence. In the mixing-layer case, we show that high compressibility inhibits helical pairings obtained in the incompressible case, yielding a staggered array of large -shaped vortices. For isotropic turbulence, one shows the existence of large coherent low-pressure vortices, characterized by exponential tails of various p.d.f.'s. We develop also a new subgrid-scale model applied with success to the backward-facing step, and to the compressible boundary layer. Finally, one looks at the influence of solid-body rotation upon free-shear flows or homogeneous turbulence. At moderate Rossby numbers, cyclonic vortices are stabilized, while anticyclonic vortices are disrupted into intense Görtler-type alternate longitudinal vortices.     

Numerical simulations of a transport-aircraft configuration

An implicit low-cost Navier–Stokes solver combined with a multigrid algorithm and wall functions has been developed for efficient numerical simulations on a realistic wing-body aircraft configuration. A study of the behaviour of different transport-equation turbulence models is given. Comparisons are made with experimental data. The structure of the three-dimensional flow separation predicted by computations is described and its topological coherence is checked.     

Numerical simulations of nozzle starting process

The starting process of two-dimensional nozzle flow is investigated both experimentally and numerically. Discussions are made on the comparison between experimental and numerical results. Performances of two numerical methods which are used in the present study of unsteady flow problem are also discussed and indications for future development of numerical tools to study nozzle problems are obtained. Received 16 June 1998 / Accepted 17 August 1998     

Numerical simulations of large falling drops

The numerical simulation in a two‐phase medium of falling drops allows their velocities and shapes during the fall to be calculated. The terminal velocities and shapes for bromoform and chlorobenzene drops falling into water have been obtained. Although the method used calculates the flow inside and around the drop, it has not been possible to give results independent of the spatial discretization and the boundary effects. However, taking these influences into account, the numerical results agree with the experimental data given and the study consists of a good validation of the SURFER code used. Copyright © 2004 John Wiley & Sons, Ltd.     

Numerical simulations of energy transfer in counter-streaming plasmas

Collisionless shock formation is investigated with large scale fully electromagnetic two-dimensional Particle-in-Cell numerical simulations. Two plasmas are colliding in the center of mass reference frame at sub-relativistic velocities. Their interaction leads to collisionless stochastic electron heating, ion slowing down and formation of a shock front. We focus here on the initial stage of evolution where electron heating is due to the Weibel-like micro-instability driven by the high-speed ion flow. A two stage process is described in the detailed analysis of our simulation results. Filament generation, followed by turbulent mixing, constitute the dominant mechanism for energy repartition. The global properties are illustrated by examination of single filament evolution in terms of energy/particle density and fields.     

Stochastic jump and bifurcation of Duffing oscillator with fractional derivative damping under combined harmonic and white noise excitations

The stochastic jump and bifurcation of Duffing oscillator with fractional derivative damping of order α (0<α<1) under combined harmonic and white noise excitations are studied. First, the system state is approximately represented by two-dimensional time-homogeneous diffusive Markov process of amplitude and phase difference using the stochastic averaging method. Then, the method of reduced Fokker–Plank–Kolmogorov (FPK) equation is used to predict the stationary response of the original system. The phenomenon of stochastic jump and bifurcation as the fractional orders' change is examined.     

单模态Richtmyer-Meshkov不稳定性数值模拟

采用VOF(Volume of Fluid)方法和PPM(Piecewise Parabolic Method)方法,发展了可用于可压缩多介质粘性流体动力学问题的数值模拟方法MVPPM(Multi-Viscous-Fluid Piecewise Parabolic Method)。利用MVPPM对多个具有不同初始扰动振幅的二维和三维单模态RM(Richtmyer-Meshkov)不稳定性模型进行了数值计算,并与理论模型的计算结果进行了比较。结果表明,无论二维还是三维情况,当初始扰动振幅相对于波长较小的时候,计算的扰动振幅和增长率与理论模型的计算结果一致。当初始扰动波长不变而振幅逐渐增大时,界面振幅和增长率也逐渐增大。对于具有相同初始扰动的情况,三维计算结果在线性段与二维计算结果相同,但是在非线性段比二维结果大,说明非线性和三维效应在RM不稳定性发展过程中起着重要作用。