Nonlinear wave propagation analysis in hyperelastic 1D microstructured materials constructed by homogenization

We analyze the acoustic properties of microstructured beams including a repetitive network material undergoing configuration changes leading to geometrical nonlinearities. The effective constitutive law is evaluated successively as an effective first and second order nonlinear grade 1D continuum, based on a strain driven incremental scheme written over the reference unit cell, taking into account the changes of the lattice geometry. The dynamical equations of motion are next written, leading to specific dispersion relations. The presence of second gradient order term in the nonlinear equation of motion leads to the presence of two different modes: an evanescent subsonic mode for high nonlinearity that vanishes beyond certain values of wave number, and a supersonic mode for a weak nonlinearity. This methodology is applied to analyze wave propagation within different microstructures, including the regular and reentrant hexagons, and plain weave textile pattern.     

Analysis of guided waves with a nonlinear boundary condition caused by internal resonance using the method of multiple scales

We theoretically investigated the cumulative nonlinear guided waves caused by internal resonance, using the method of multiple scales (MMS), which can construct better approximations to the solutions of perturbation problems. In this study, we consider nonlinearity only on the boundary instead of material nonlinearity or geometric nonlinearity. We showed nonlinear effects on the amplitudes of a lower mode and a higher mode depending on the propagation length. Also, we examined effects of wavenumber detuning from a phase matching condition of the two modes. If the wavenumber detuning is exactly equal to zero, the mechanical energy of the lower mode is transferred through nonlinear coupling to the energy of the higher mode, unilaterally. However, if a wavenumber detuning is not equal to zero, amplitude of the two modes change in a cyclic fashion during wave propagation. The amount of this amplitude variation and its cycle length are determined by the eigenfunctions of the two modes, the nonlinear parameter and the wavenumber detuning.     

Wave propagation along a non-principal direction in a compressible pre-stressed elastic layer

The dispersive behavior of small amplitude waves propagating along a non-principal direction in a pre-stressed, compressible elastic layer is considered. One of the principal axes of stretch is normal to the elastic layer and the direction of propagation makes an angle θ with one of the in-plane principal axes. The dispersion relations which relate wave speed and wavenumber are obtained for both symmetric and anti-symmetric motions by formulating the incremental boundary value problem for a general strain energy function. The behavior of the dispersion curves for symmetric waves is for the most part similar to that of the anti-symmetric waves at the low and high wavenumber limits. At the low wavenumber limit, depending on the pre-stress and propagation angle, it may be possible for both the fundamental mode and the next lowest mode to have finite phase speeds, while other higher modes have an infinite phase speed. At the high wavenumber limit, the phase speeds of the fundamental mode and the higher modes tend to the Rayleigh surface wave speed and the limiting wave speeds of the layer, respectively. Numerical results are presented for a Blatz–Ko material and the effect of the propagation angle is clearly illustrated.     

Weakly nonlinear wave interactions in multi-degree of freedom periodic structures

This work presents a multiple time scales perturbation analysis for analyzing weakly nonlinear wave interactions in multi-degree of freedom periodic structures. The perturbation analysis is broadly applicable to (discretized) periodic systems in any dimensional space and with a wide range of constitutive nonlinearities. Specific emphasis is placed on cubic nonlinearity, as dispersion shifts typically arise from the cubic components in nonlinear restoring forces. The procedure is first presented in general. Then, application to the diatomic chain and monoatomic two-dimensional lattice demonstrates, individually, the treatment of multiple degree of freedom systems and higher dimensional spaces. The dispersion relations are modified by weakly nonlinear wave interactions and lead to additional opportunities to control wave propagation direction, band gap size, and group velocity. Numerical simulations validate the expected dispersion shifts. An amplitude-tunable focus device demonstrates the viability of utilizing dynamically-introduced dispersion to produce beam steering that may, ultimately, lead to a phononic superprism effect as well as multiplexing/demultiplexing behavior.     

Aeroelastic analysis and active flutter control of nonlinear lattice sandwich beams

Nonlinear aeroelastic characteristics of sandwich beams with pyramidal lattice core are investigated, and the active flutter control of the nonlinear structural system is also studied using the piezoelectric actuator/sensor pair. In the structural modeling, Reddy’s third-order shear deformation theory is applied. Aerodynamic pressure is evaluated by the supersonic piston theory. Hamilton’s principle and the assumed mode method are used to derive the equation of motion. The proportional feedback and the optimal H _∞ control methods are performed to design the controller. In the robust control, the uncertainty caused by omitting the nonlinear terms of the control equation is taken into account, and the mixed sensitivity method is used to solve the problem. The nonlinear aeroelastic property of the sandwich beam is analyzed and is compared with that of the equivalent isotropic beam with the same weight to show the superior aeroelastic characteristics of the lattice sandwich beam. Controlled vibration responses under the two different controllers are calculated and compared. Simulation results show that the robust controller is much more effective than the proportional feedback controller in the flutter suppression of the nonlinear sandwich beam.     

Localized solutions to the nonlinear wave equation for flexural motions of an elastic beam traveling in an air-filled tube

Steady-progressive-wave solutions are sought to the nonlinear wave equation derived previously [J. Fluids Struct. 16 (2002) 597] for flexural motions of an elastic beam traveling in an air-filled tube along its center axis at a subsonic speed. Fluid-structure interactions are taken into account through aerodynamic loading on the lateral surface of the beam subjected to small but finite deflection but end effects and viscous effects are neglected. Linear dispersion characteristics are first examined by exploiting the small ratio of the induced mass to the mass of the beam per unit length. Centered around the traveling speed of the beam, there exists such a narrow range of propagation velocity that the linear steady propagation is prohibited. In this range, it is revealed that some interesting nonlinear solutions exist. The periodic wavetrain is found to exist as the exact solution. Asymptotic analysis is then made by applying the method of multiple scales and the stationary nonlinear Schrödinger equation is derived for a complex amplitude. A monochromatic solution to this equation corresponds to the exact periodic solution. Imposing undisturbed boundary conditions at infinity, it is revealed that the localized solution exists as a result of balance between the linear instability and the nonlinearity. This solution is checked by solving the nonlinear equation numerically. It is further revealed that the amplitude-modulated wavetrain exists not only in the range of the velocity mentioned above but also outside of it.     

Dispersion effects of extensional waves in pre-stressed imperfectly bonded incompressible elastic layered composites

The effect of an imperfect interface, on time-harmonic extensional wave propagation in a pre-stressed symmetric layered composite is considered. The bimaterial composite consists of incompressible isotropic elastic materials. The shear spring type resistance model employed to simulate the imperfect interface can accommodate the extreme cases of perfect bonding and a fully slipping interface. The dispersion relation obtained by formulating the incremental boundary-value problem and the use of the propagator matrix technique, is analyzed at the low and high wavenumber limits. For the perfectly bonded and imperfect interface cases in the low wavenumber region, only the fundamental mode has a finite phase speed, while other higher modes have an infinite phase speed when the dimensionless wavenumber approaches zero. However, for the fully slipping interface in the low wavenumber region, both the fundamental mode and the next lowest mode have finite phase speeds. In the high wavenumber region, when the dimensionless wavenumber tends to infinity, the phase speeds of the fundamental mode and the higher modes depend on the phase speeds of the surface and interfacial waves and on the limiting phase speed of the composite. An expression to determine the cut-off frequencies is obtained from the dispersion relation. Numerical examples of dispersion curves are presented, where when the material has to be prescribed either Mooney–Rivlin material or Varga material is assumed. The effect of the imperfect interface is clearly evident in the numerical results.     

Symplectic analysis for wave propagation in one-dimensional nonlinear periodic structures

The wave propagation problem in the nonlinear periodic mass-spring structure chain is analyzed using the symplectic mathematical method. The energy method is used to construct the dynamic equation, and the nonlinear dynamic equation is linearized using the small parameter perturbation method. Eigen-solutions of the symplectic matrix are used to analyze the wave propagation problem in nonlinear periodic lattices. Nonlinearity in the mass-spring chain, arising from the nonlinear spring stiffness effect, has profound effects on the overall transmission of the chain. The wave propagation characteristics are altered due to nonlinearity, and related to the incident wave intensity, which is a genuine nonlinear effect not present in the corresponding linear model. Numerical results show how the increase of nonlinearity or incident wave amplitude leads to closing of transmitting gaps. Comparison with the normal recursive approach shows effectiveness and superiority of the symplectic method for the wave propagation problem in nonlinear periodic structures.     

Time-harmonic wave propagation in a pre-stressed compressible elastic bi-material laminate

The dispersive behaviour of time-harmonic waves propagating along a principal direction in a perfectly bonded pre-stressed compressible elastic bi-material laminate is considered. The dispersion relation which relates wave speed and wavenumber is obtained by formulating the incremental boundary value problem and the use of the propagator matrix technique. At the low wavenumber limit, depending on the pre-stress, both the fundamental mode and the next lowest mode may have finite phase speeds. For the higher modes which have infinite phase speeds in the low wavenumber region, an expression to determine the cut-off frequencies is obtained. At the high wavenumber limit, the phase speeds of the fundamental mode and higher modes tend to phase speeds of the surface wave, the interfacial wave or the limiting phase speed of the composite. For numerical examples, either a two-parameter compressible neo-Hookean material or a two-parameter compressible Varga material is assumed.     

Trigonometric simplification of a class of conservative nonlinear oscillators

This paper deals with a class of conservative nonlinear oscillators of the form $\ddot x(t)+f(x(t))=0$ , where f(x) is analytic. A transformation of time from t to a new time coordinate τ is defined such that periodic solutions can be expressed in the form x(τ) = A ₀+A ₁ cos 2τ. We refer to this process of trigonometric simplification as trigonometrification. Application is given to the stability of nonlinear normal modes (NNMs) in two-degree-of-freedom systems.     

Experimental investigation of combustion mechanisms of kerosene-fueled scramjet engines with double-cavity flameholders

A scramjet combustor with double cavitybased flameholders was experimentally studied in a directconnected test bed with the inflow conditions of M = 2.64,Pt = 1.84 MPa,Tt = 1 300 K.Successful ignition and selfsustained combustion with room temperature kerosene was achieved using pilot hydrogen,and kerosene was vertically injected into the combustor through 4×φ 0.5 mm holes mounted on the wall.For different equivalence ratios and different injection schemes with both tandem cavities and parallel cavities,flow fields were obtained and compared using a high speed camera and a Schlieren system.Results revealed that the combustor inside the flow field was greatly influenced by the cavity installation scheme,cavities in tandem easily to form a single side flame distribution,and cavities in parallel are more likely to form a joint flame,forming a choked combustion mode.The supersonic combustion flame was a kind of diffusion flame and there were two kinds of combustion modes.In the unchoked combustion mode,both subsonic and supersonic combustion regions existed.While in the choked mode,the combustion region was fully subsonic with strong shock propagating upstream.Results also showed that there was a balance point between the boundary separation and shock enhanced combustion,depending on the intensity of heat release.     

Theoretical and numerical investigation of HF elastic wave propagation in two-dimensional periodic beam lattices

The elastic wave propagation phenomena in two-dimensional periodic beam lattices are studied by using the Bloch wave transform. The numerical modeling is applied to the hexagonal and the rectangular beam lattices, in which, both the in-plane （with respect to the lattice plane） and out-of-plane waves are considered. The dispersion relations are obtained by calculating the Bloch eigenfrequencies and eigenmodes. The frequency bandgaps are observed and the influence of the elastic and geometric properties of the primitive cell on the bandgaps is studied. By analyzing the phase and the group velocities of the Bloch wave modes, the anisotropic behaviors and the dispersive characteristics of the hexagonal beam lattice with respect to the wave prop- agation are highlighted in high frequency domains. One im- portant result presented herein is the comparison between the first Bloch wave modes to the membrane and bend- ing/transverse shear wave modes of the classical equivalent homogenized orthotropic plate model of the hexagonal beam lattice. It is shown that, in low frequency ranges, the homog- enized plate model can correctly represent both the in-plane and out-of-plane dynamic behaviors of the beam lattice, its frequency validity domain can be precisely evaluated thanks to the Bloch modal analysis. As another important and original result, we have highlighted the existence of the retro- propagating Bloch wave modes with a negative group veloc- ity, and of the corresponding ＂retro-propagating＂ frequency bands.     

Second-harmonic generation of two-dimensional elastic wave propagation in an infinite layered structure with nonlinear spring-type interfaces

The two-dimensional elastic wave propagation in an infinite layered structure with nonlinear interlayer interfaces is analyzed theoretically to investigate the second-harmonic generation due to interfacial nonlinearity. The structure consists of identical isotropic linear elastic layers that are bonded to each other by spring-type interfaces possessing identical linear normal and shear stiffnesses but different quadratic nonlinearity parameters. Explicit analytical expressions are derived for the second-harmonic amplitudes when a single monochromatic Bloch mode propagates in the structure in arbitrary directions by applying the transfer-matrix approach and the Bloch theorem to the governing equations linearized by a perturbation method. The second-harmonic generation by a single nonlinear interface and by multiple consecutive nonlinear interfaces are shown to be profoundly influenced by the band structure of the layered structure, the fundamental Bloch wave mode, and its propagation direction. In particular, the second harmonics generated at multiple consecutive interfaces are found to grow cumulatively with the propagation distance when the phase matching occurs between the Bloch modes at the fundamental and double frequencies.     

Existence of supersonic zones in three-dimensional separation flows

Up till now the region of three-dimensional separation flows which occur with supersonic flow past obstacles has received insufficient study. Supersonic flow with a Mach number of 2.5 past a cylinder mounted on a plate was studied in [1]. A local zone with supersonic velocities was found in the reverse subsonic flow region ahead of the cylinder. Its presence is explained by the three-dimensional nature of the flow. Similar supersonic zones are not observed in the case of supersonic flow over plane and axisymmetric steps.The present paper presents the results of experimental studies whose objective was refinement of the flow pattern ahead of a cylinder on a plate and the study of the local supersonic zones.The experiments were performed in a supersonic wind tunnel with a freestream Mach number M₁=3.11. The 24-mm-diameter cylinder with pressure taps along the generating line was mounted perpendicular to the surface of a sharpened plate. The distance from the plate leading edge to the cylinder axis wasl ₀=140 mm. The plate was pressure tapped along the flow symmetry axis. The Reynolds number was Rl ₀=u₀ l ₀/v ₁, Rl ₀=1.87.10⁷, where u₁ andv ₁ are the freestream velocity and the kinematic viscosity, respectively. The pressures were measured using a Pilot probe with internal and external diameters of 0.15 and 0.9 mm, respectively.The probe was displaced in the flow symmetry plane at a distance of 1.6 mm from the plate surface and at a distance of 1.1 mm along the leading generator of the cylinder. The flow on the surface of the plate and cylinder was studied with the aid of a visualization composition and the flow past the model was photographed with a schlieren instrument. Typical patterns of the visualization composition distribution and the pressure distribution curves over the plate surface, and also photographs of the flow past the model, are shown in [1].     

Dynamics of light bullets in inhomogeneous cubic-quintic-septimal nonlinear media with $${\varvec{{\mathcal {}}}}$$-symmetric potentials

A (3+1)-dimensional nonlinear Schrödinger equation with variable-coefficient dispersion/diffraction and cubic-quintic-septimal nonlinearities is studied, two families of analytical light bullet solutions with two types of \({{\mathcal {PT}}}\)-symmetric potentials are obtained. The coefficient of the septimal nonlinear term strongly influences the form of light bullet. The direct numerical simulation indicates that light bullet solutions in different cubic-quintic-septimal nonlinear media exhibit different property of stability, and under different \({\mathcal {PT}}\)-symmetric potentials they also show different stability against white noise. These stabilities of evolution originate from subtle interplay among dispersion, diffraction, nonlinearity and \({\mathcal {PT}}\)-symmetric potential. Moreover, compression and expansion of light bullets in the hyperbolic dispersion/diffraction system and periodic modulation system are investigated numerically. The evolution of light bullet in periodic modulation system is more stable than that in the hyperbolic dispersion/diffraction system.     

Anti-symmetric waves in pre-stressed imperfectly bonded incompressible elastic layered composites

The effect of an imperfect interface on the dispersive behavior of in-plane time-harmonic symmetric waves in a pre-stressed incompressible symmetric layered composite, was analyzed recently by Leungvichcharoen and Wijeyewickrema (2003). In the present paper the corresponding case for time harmonic anti-symmetric waves is considered. The bi-material composite consists of incompressible isotropic elastic materials. The imperfect interface is simulated by a shear-spring type resistance model, which can also accommodate the extreme cases of perfectly bonded and fully slipping interfaces. The dispersion relation is obtained by formulating the incremental boundary-value problem and using the propagator matrix technique. The dispersion relations for anti-symmetric and symmetric waves differ from each other only through the elements of the propagator matrix associated with the inner layer. The behavior of the dispersion curves for anti-symmetric waves is for the most part similar to that of symmetric waves at the low and high wavenumber limits. At the low wavenumber limit, depending on the pre-stress for perfectly bonded and imperfect interface cases, a finite phase speed may exist only for the fundamental mode while other higher modes have an infinite phase speed. However, for a fully slipping interface in the low wavenumber region it may be possible for both the fundamental mode and the next lowest mode to have finite phase speeds. For the higher modes which have infinite phase speeds in the low wavenumber region an expression to determine the cut-off frequencies is obtained. At the high wavenumber limit, the phase speeds of the fundamental mode and the higher modes tend to the phase speeds of the surface wave or the interfacial wave or the limiting phase speed of the composite. The bifurcation equation obtained from the dispersion relation yields neutral curves that separate the stable and unstable regions associated with the fundamental mode or the next lowest mode. Numerical examples of dispersion curves are presented, where when the material has to be prescribed either Mooney–Rivlin material or Varga material is assumed. The effect of imperfect interfaces on anti-symmetric waves is clearly evident in the numerical results.     

Transonic Shocks for the Full Compressible Euler System in a General Two-Dimensional De Laval Nozzle

In this paper, we study the transonic shock problem for the full compressible Euler system in a general two-dimensional de Laval nozzle as proposed in Courant and Friedrichs (Supersonic flow and shock waves, Interscience, New York, 1948): given the appropriately large exit pressure p _e(x), if the upstream flow is still supersonic behind the throat of the nozzle, then at a certain place in the diverging part of the nozzle, a shock front intervenes and the gas is compressed and slowed down to subsonic speed so that the position and the strength of the shock front are automatically adjusted such that the end pressure at the exit becomes p _e(x). We solve this problem completely for a general class of de Laval nozzles whose divergent parts are small and arbitrary perturbations of divergent angular domains for the full steady compressible Euler system. The problem can be reduced to solve a nonlinear free boundary value problem for a mixed hyperbolic–elliptic system. One of the key ingredients in the analysis is to solve a nonlinear free boundary value problem in a weighted Hölder space with low regularities for a second order quasilinear elliptic equation with a free parameter (the position of the shock curve at one wall of the nozzle) and non-local terms involving the trace on the shock of the first order derivatives of the unknown function.     

Topology design and optimization of nonlinear periodic materials

This paper explores optimal topologies yielding large band gap shifts in one- and two-dimensional nonlinear periodic materials. The presence of a nonlinearity in a periodic material system results in amplitude-dependent dispersion behavior, leading to novel wave-based devices such as tunable filters, resonators, and waveguides. The performance of these devices over a broad frequency range requires large, tunable band gaps, motivating the present study. Consideration of a one-dimensional bilayer system composed of alternating linear and nonlinear layers shows that optimal designs consist of thin, compliant nonlinear layers. This is at first surprising considering the source of the shift originates from only the nonlinear layer; however, thin layers lead to localized stresses that activate the nonlinear character of the system. This trend persists in two-dimensional materials where optimization studies are performed on plane-stress models discretized using bilinear Lagrange elements. A fast algorithm is introduced for computing the dispersion shifts, enabling efficient parametric analyses of two-dimensional inclusion systems. Analogous to the one-dimensional system, it is shown that thin ligaments of nonlinear material lead to large dispersion shifts and group velocity variations. Optimal topologies of the two-dimensional system are also explored using genetic algorithms aimed at producing large increases in complete band gap width and shift, or group velocity variation, without presupposing the topology. The optimal topologies that result resemble the two-dimensional inclusion systems, but with small corner features that tend to enhance the production of dispersion shift further. Finally, the study concludes with a discussion on Bloch wave modes and their important role in the production of amplitude-dependent dispersion behavior. The results of the study provide insight and guidance on selecting topologies and materials which can yield large amplitude-dependent band gap shifts and group velocity variations, thus enabling sensitive nonlinear devices.     

Global Uniqueness of Transonic Shocks in Two-Dimensional Steady Compressible Euler Flows

We prove that for the two-dimensional steady complete compressible Euler system, with given uniform upcoming supersonic flows, the following three fundamental flow patterns (special solutions) in gas dynamics involving transonic shocks are all unique in the class of piecewise C ¹ smooth functions, under appropriate conditions on the downstream subsonic flows: (i) the normal transonic shocks in a straight duct with finite or infinite length, after fixing a point the shock-front passing through; (ii) the oblique transonic shocks attached to an infinite wedge; (iii) a flat Mach configuration containing one supersonic shock, two transonic shocks, and a contact discontinuity, after fixing a point where the four discontinuities intersect. These special solutions are constructed traditionally under the assumption that they are piecewise constant, and they have played important roles in the studies of mathematical gas dynamics. Our results show that the assumption of a piecewise constant can be replaced by some weaker assumptions on the downstream subsonic flows, which are sufficient to uniquely determine these special solutions. Mathematically, these are uniqueness results on solutions of free boundary problems of a quasi-linear system of elliptic-hyperbolic composite-mixed type in bounded or unbounded planar domains, without any assumptions on smallness. The proof relies on an elliptic system of pressure p and the tangent of the flow angle w = v/u obtained by decomposition of the Euler system in Lagrangian coordinates, and a newly developed method for the L ^∞ estimate that is independent of the free boundaries, by combining the maximum principles of elliptic equations, and careful analysis of the shock polar applied on the (maybe curved) shock-fronts.     

Analysis of non-axisymmetric wave propagation in a homogeneous piezoelectric solid circular cylinder of transversely isotropic material

A study concerning the propagation of free non-axisymmetric waves in a homogeneous piezoelectric cylinder of transversely isotropic material with axial polarization is carried out on the basis of the linear theory of elasticity and linear electro-mechanical coupling. The solution of the three dimensional equations of motion and quasi-electrostatic equation is given in terms of seven mechanical and three electric potentials. The characteristic equations are obtained by the application of the mechanical and two types of electric boundary conditions at the surface of the piezoelectric cylinder. A novel method of displaying dispersion curves is described in the paper and the resulting dispersion curves are presented for propagating and evanescent waves for PZT-4 and PZT-7A piezoelectric ceramics for circumferential wave numbers m = 1, 2, and 3. It is observed that the dispersion curves are sensitive to the type of the imposed boundary conditions as well as to the measure of the electro-mechanical coupling of the material.