Stability of resonant interfacial waves in the presence of uniform magnetic field

The aim of this work is to study the instability of interacting waves between two immiscible magnetic liquids. The effects of gravitation and a uniform normal magnetic field are taken into account. The method of multiple scales is used to determine the stability criteria of the considered problem. The various stability criteria are discussed both analytically and graphically. According to the numerical examples, we have remarked that the increase of the ratio of the permeability of the liquids appears to be the destabilizing effect of the magnetic field. The short waves below the critical wavenumbers are stable whereas a number of long waves are unstable. The viscosity effect on the stability criteria is a dual-role one, depending on the strength of the applied magnetic field.     

Trigonometric shock waves in the Kaup–Boussinesq system

Nonlinear Dynamics - We consider the modulationally stable version of the Kaup–Boussinesq system which models propagation of nonlinear waves in various physical situations. It is shown that...     

Euler–Lagrange simulation of high pressure shock waves

This paper presents the numerical simulation of overdriven detonation (or O.D.D.) that occurs when a high velocity object impacts an explosive. The pressure and the velocity at this state are higher than those of the Chapman–Jouguet (C–J) state. First, before the simulation of this event, a study of PBX air blast by using multi-material Eulerian method is presented. Pressure peaks are computed for several distances from the explosive. Second, the O.D.D. phenomenon is modeled by the Euler–Lagrange penalty coupling, which permits to couple a Lagrangian mesh of the flyer plate to multi-material Eulerian mesh of explosives and air. This coupling gives us the high detonation velocities in the acceptor explosive and demonstrates that it is able to handle shock–structure interaction problems.     

Dispersive shock waves in three dimensional Benjamin–Ono equation

Dispersive shock waves (DSWs) in the three dimensional Benjamin–Ono (3DBO) equation are studied with step-like initial condition along a paraboloid front. By using a similarity reduction, the problem of studying DSWs in three space one time (3+1) dimensions reduces to finding DSW solution of a (1+1) dimensional equation. By using a special ansatz, the 3DBO equation exactly reduces to the spherical Benjamin–Ono (sBO) equation. Whitham modulation equations are derived which describes DSW evolution in the sBO equation by using a perturbation method. These equations are written in terms of appropriate Riemann type variables to obtain the sBO-Whitham system. DSW solution which is obtained from the numerical solutions of the Whitham system and the direct numerical solution of the sBO equation are compared. In this comparison, a good agreement is found between these solutions. Also, some physical qualitative results about DSWs in sBO equation are presented. It is concluded that DSW solutions in the reduced sBO equation provide some information about DSW behavior along the paraboloid fronts in the 3DBO equation.     

Aerodynamic characteristics of solid particles�� acceleration by shock waves

In this study, the interaction of a planar shock wave with a group of particles has been investigated using high-speed photography and dynamic pressure measurements. Experiments were carried out in a horizontal circular shock tube. The influence of the particle loading ratio, particle diameter, driving gas and shock wave Mach number on the acceleration was studied. It was found that the higher the particle loading ratio, the greater was the particle velocity. This is due to the higher driving pressure. Helium and nitrogen gases play quite different roles in acceleration. Pressure multiplication during shock wave interaction with particles also appears. Based on the experimental results, the discussion regarding partial quantitative velocities and accelerations of particle groups, as well as the attenuation factors when shock waves pass through the particles, is given.     

Approximate analytical solution to oscillations of a conductor in a magnetic field

An approximate analytical solution to a system exhibiting oscillations of a conductor in a magnetic field which is controlled by a discrete waveform is sought by means of multiple scales. The system involves the use of a solenoid driven by a RLC circuit, coupled with a solid state relay (SSRL), to generate large electromagnetic forces acting on a conductor, which oscillates within the solenoid. The steady state response of the metal bar, in terms of oscillations is described. This solution is expressed in terms of system and circuit parameters, valid in the weakly nonlinear region, which is identified to be small oscillatory displacement near the center of the solenoid. By analyzing different cases of resonance, period-1 and period-2 like motions are identified and validated through experimental studies. The solution provides a guideline to design an effective control strategy so as to guide the system to a desirable attractor.     

Linear Rayleigh-Taylor instability analysis of double-shell Kidder＇s self-similar implosion solution

This paper generalizes the single-shell Kidder＇s self-similar solution to the double-shell one with a discontinuity in density across the interface. An isentropic implosion model is constructed to study the Rayleigh-Taylor instability for the implosion compression. A Godunov-type method in the Lagrangian coordinates is used to compute the one-dimensional Euler equation with the initial and boundary conditions for the double-shell Kidder＇s self-similar solution in spherical geometry. Numerical results are obtained to validate the double-shell implosion model. By programming and using the linear perturbation codes, a linear stability analysis on the Rayleigh-Taylor instability for the double-shell isentropic implosion model is performed. It is found that, when the initial perturbation is concentrated much closer to the interface of the two shells, or when the spherical wave number becomes much smaller, the modal radius of the interface grows much faster, i.e., more unstable. In addition, from the spatial point of view for the compressibility effect on the perturbation evolution, the compressibility of the outer shell has a destabilization effect on the Rayleigh-Taylor instability, while the compressibility of the inner shell has a stabilization effect.     

Almost periodic solution of a non-autonomous model of phytoplankton allelopathy

In this paper, a non-autonomous two species competitive allelopathic phytoplankton model in presence of a discrete time delay is considered. We have obtained the sufficient conditions for permanence along with existence-uniqueness of an almost periodic solution. Sufficient conditions are derived for the existence of unique almost periodic solution. Analytical findings are supported through exhaustive numerical simulations. With the help of the numerical example, we have demonstrated that initial density dependent almost periodic co-existence is possible in some situations when parameter values fail to satisfy all the conditions of permanence.     

A semi-analytical solution for the confined compression of hydrated soft tissue

Confined compression is a common experimental technique aimed at gaining information on the properties of biphasic mixtures comprised of a solid saturated by a fluid, a typical example of which are soft hydrated biological tissues. When the material properties (elastic modulus, permeability) are assumed to be homogeneous, the governing equation in the axial displacement reduces to a Fourier equation which can be solved analytically. For the more realistic case of inhomogeneous material properties, the governing equation does not admit, in general, a solution in closed form. In this work, we propose a semi-analytical alternative to Finite Element analysis for the study of the confined compression of linearly elastic biphasic mixtures. The partial differential equation is discretised in the space variable and kept continuous in the time variable, by use of the Finite Difference Method, and the resulting system of ordinary differential equations is solved by means of the Laplace Transform method.     

Biomechanical effects of shock waves onEscherichia coli and λphage DNA

Escherichia coli (recombinant cells) and phage DNA in suspension liquid were exposed to pressure pulses of about 20s duration and amplitude of up to 14 MPa. These pulses were generated by a diaphragmless shock tube. The destruction of cells was monitored by the assay of phenylalanine dehydrogenase leaking from the recombinant cells and was found to increase remarkably at the peak pressure of higher than 12 MPa. A probability relation for the cell destruction expressed as a function of pressure was proposed. It is most likely that there exists a threshold pressure for the cell destruction. Fragmentation effects of shock waves on phage DNA were analyzed by electrophoresis. They were enhanced by increasing the shock wave strength and the number of shots. Probability for the DNA fragmentation as a function of pressure and molecular size was estimated with HPLC. The larger size of the DNA was more easily fragmented. A threshold pressure does not seem to exist for the DNA fragmentation.     

Basic solution of a mode-I limited-permeable crack in functionally graded piezoelectric materials

In this paper, the basic solution of a mode-I crack in functionally graded piezoelectric materials was investigated by using the generalized Almansi’s theorem. In the analysis, the electric permittivity of air inside the crack were considered. To make the analysis tractable, it was assumed that the shear modulus, piezoelectric constants and dielectric constants vary exponentially with coordinate parallel to the crack. The problem was formulated through Fourier transform into two pairs of dual integral equations, in which the unknown variables are jumps of displacements across the crack surfaces. To solve the dual integral equations, the jumps of displacements across the crack surfaces were directly expanded as a series of Jacobi polynomials. The solution of the present paper shows that the effects of the electric boundary conditions on the electric displacement fields near the crack tips can not be ignored. Simultaneously, the solution of the present paper will revert to a closed form one when the functionally graded parameter equals to zero.     

A possible explanation of “bamboo waves” in core-annular flow of two liquids

We derive a perturbed KdV equation to model interfacial waves in viscous two-layer flows. We tentatively speculate that solutions of such an equation might be related to bamboo waves observed in recent experiments.This research was supported by the National Science Foundation under Grant DMS-9008497 and by the Office of Naval Research under Grant N00014-92-J-1664.     

Implicit and explicit integration in the solution of the absolute nodal coordinate differential/algebraic equations

This investigation is concerned with the use of an implicit integration method with adjustable numerical damping properties in the simulation of flexible multibody systems. The flexible bodies in the system are modeled using the finite element absolute nodal coordinate formulation (ANCF), which can be used in the simulation of large deformations and rotations of flexible bodies. This formulation, when used with the general continuum mechanics theory, leads to displacement modes, such as Poisson modes, that couple the cross section deformations, and bending and extension of structural elements such as beams. While these modes can be significant in the case of large deformations, and they have no significant effect on the CPU time for very flexible bodies; in the case of thin and stiff structures, the ANCF coupled deformation modes can be associated with very high frequencies that can be a source of numerical problems when explicit integration methods are used. The implicit integration method used in this investigation is the Hilber–Hughes–Taylor method applied in the context of Index 3 differential-algebraic equations (HHT-I3). The results obtained using this integration method are compared with the results obtained using an explicit Adams-predictor-corrector method, which has no adjustable numerical damping. Numerical examples that include bodies with different degrees of flexibility are solved in order to examine the performance of the HHT-I3 implicit integration method when the finite element absolute nodal coordinate formulation is used. The results obtained in this study show that for very flexible structures there is no significant difference in accuracy and CPU time between the solutions obtained using the implicit and explicit integrators. As the stiffness increases, the effect of some ANCF coupled deformation modes becomes more significant, leading to a stiff system of equations. The resulting high frequencies are filtered out when the HHT-I3 integrator is used due to its numerical damping properties. The results of this study also show that the CPU time associated with the HHT-I3 integrator does not change significantly when the stiffness of the bodies increases, while in the case of the explicit Adams method the CPU time increases exponentially. The fundamental differences between the solution procedures used with the implicit and explicit integrations are also discussed in this paper.     

Local and global exponential output synchronization of complex delayed dynamical networks

This paper proposes a new complex dynamical network model with output coupling. This model is totally different from some existing complex network models. It is well known that the node state is difficult to be observed or measured, even the node state cannot be observed or measured at all. Moreover, sometimes only part states are needed to make synchronization to come true. For these phenomena, the output synchronization is investigated in this paper. Several criteria on local and global exponential output synchronization are derived for the proposed network model. Finally, numerical simulations are given to illustrate the effectiveness of the derived results.     

Thermodynamic fluctuations in canonical shock–turbulence interaction: effect of shock strength

In this work, we use numerical simulation and linear inviscid theory to study the thermodynamic field generated by the interaction of a shock wave with homogeneous isotropic turbulence. Fluctuations in density, pressure, temperature and entropy can play an important role in shock-induced mixing, combustion and energy transfer processes. Data from shock-captured direct numerical simulations (scDNS) are used to investigate the variation of thermodynamic fluctuations for varying shock strengths, and the results are compared with linear interaction analysis (LIA). The density, pressure and temperature variances attain large values at the shock, followed by, in general, a rapid decay in the downstream flow. The rapid variation behind the shock makes it difficult to compare numerical results with theoretical predictions. A threshold method based on instantaneous shock dilatation is used to overcome this problem, and it gives excellent match between scDNS and LIA. We find cases with non-monotonic variation with Mach number as well as local peaks in density fluctuations behind the shock. These are explained in terms of the contribution of the post-shock acoustic and entropy modes in the LIA solution and their cross-correlation. Budget of the transport equations reveals interesting insight into the physics governing the thermodynamic field behind the shock wave. It is found that the variances are primarily determined by the competing effects of dilatational and dissipation mechanisms. The dominant mechanisms are identified for a range of conditions, and their implication for developing predictive models is highlighted.     

“Non-free” interaction of plane shock waves and the boundary layer in the neighborhood of the leading edge of a plate with slip

The interaction between a normally impinging shock wave and the boundary layer on a plate with slip is studied in the neighborhood of the leading edge using various experimental methods, including special laser technology, to visualize the supersonic conical gas flows. It is found that in the “non-free” interaction, when the leading edge impedes the propagation of the boundary layer separation line upstream, the structure of the disturbed flow is largely identical to that in the developed “free” interaction, but with higher parameter values and gradients in the leading part of the separation zone. The fundamental property of developed separation flows, namely, coincidence of the values of the pressure “plateau” in the separation zone and the pressure behind the oblique shock above the separation zone of the turbulent boundary layer, is conserved. Moscow. e-mail: ostap@inmech.msu.su. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 3, pp. 57–69, May–June, 2000. The work was carried out with financial support from the Russian Foundation for Basic Research (project No. 97-01-00099).     

A Theory of General Solutions of Plane Problems in Two-Dimensional Octagonal Quasicrystals

A theory of general solutions of plane problems is developed for the coupled equations in plane elasticity of two-dimensional octagonal quasicrystals. In virtue of the operator method, the general solutions of the antiplane and inplane problems are given constructively with two displacement functions. The introduced displacement functions have to satisfy higher order partial differential equations, and therefore it is difficult to obtain rigorous analytic solutions directly and is not applicable in most cases. In this case, a decomposition and superposition procedure is employed to replace the higher order displacement functions with some lower order displacement functions, and accordingly the general solutions are further simplified in terms of these functions. In consideration of different cases of characteristic roots, the general solution of the antiplane problem involves two cases, and the general solution of the inplane problem takes three cases, but all are in simple forms that are convenient to be applied. Furthermore, it is noted that the general solutions obtained here are complete in x ₃-convex domains.      

Multiple scattering from two cylindrical inclusions in a semi-infinite thin plate under flexural waves

In this study, a theoretical method is applied to investigate the multiple scattering of flexural waves from two cylindrical inclusions and dynamic stress in a semi-infinite thin plate, and the roller-supported boundary at the semi-infinite edge is considered. A computationally efficient approach utilizing the wave expansion method together with the image method is employed to formulate the problem. The addition theorem for Bessel functions is employed to accomplish the translation of wave fields between different local coordinate systems. The scattered far field amplitudes of the two inclusions are presented. Through numerical examples and analyses, it is found that the angular distribution of the dynamic stress around the inclusions shows a great difference when the positions of two inclusions vary. The effects of the elastic modulus, density, Poisson’s ratio of the inclusions on the dynamic stress are closely related to the distance between the inclusions and the semi-infinite edge and the distance between the two inclusions. The thickness and Poisson’s ratio of the inclusions show little effect on the dynamic stress. Comparisons with other existing models are also discussed.     

On delay-dependent robust exponential stability of stochastic neural networks with mixed time delays and Markovian switching

This paper deals with the global exponential stability analysis problem for a general class of uncertain stochastic neural networks with mixed time delays and Markovian switching. The mixed time delays under consideration comprise both the discrete time-varying delays and the distributed time-delays. The main purpose of this paper is to establish easily verifiable conditions under which the delayed stochastic neural network is robustly exponentially stable in the mean square in the presence of parameters uncertainties, mixed time delays, and Markovian switching. By employing new Lyapunov–Krasovskii functionals and conducting stochastic analysis, a linear matrix inequality (LMI) approach is developed to derive the criteria for the robust exponential stability, which can be readily checked by using some standard numerical packages such as the Matlab LMI Toolbox. The criteria derived are dependent on both the discrete time delay and distributed time delay, and, are therefore, less conservative. A simple example is provided to demonstrate the effectiveness and applicability of the proposed testing criteria. This work was supported in part by the Engineering and Physical Sciences Research Council (EPSRC) of the UK under Grant GR/S27658/01, the Nuffield Foundation of the UK under Grant NAL/00630/G, the National Natural Science Foundation of China under Grant 60774073, the Natural Science Foundation of Jiangsu Province of China under Grant BK2007075, the Natural Science Foundation of Jiangsu Education Committee of China under Grant 06KJD110206, the Scientific Innovation Fund of Yangzhou University of China under Grant 2006CXJ002, and the Alexander von Humboldt Foundation of Germany.