Wave propagation in a 2D nonlinear structural acoustic waveguide using asymptotic expansions of wavenumbers

Nonlinear acoustic wave propagation in an infinite rectangular waveguide is investigated. The upper boundary of this waveguide is a nonlinear elastic plate, whereas the lower boundary is rigid. The fluid is assumed to be inviscid with zero mean flow. The focus is restricted to non-planar modes having finite amplitudes. The approximate solution to the acoustic velocity potential of an amplitude modulated pulse is found using the method of multiple scales (MMS) involving both space and time. The calculations are presented up to the third order of the small parameter. It is found that at some frequencies the amplitude modulation is governed by the Nonlinear Schrodinger equation (NLSE). The first objective here is to study the nonlinear term in the NLSE. The sign of the nonlinear term in the NLSE plays a role in determining the stability of the amplitude modulation. Secondly, at other frequencies, the primary pulse interacts with its higher harmonics, as do two or more primary pulses with their resultant higher harmonics. This happens when the phase speeds of the waves match and the objective is to identify the frequencies of such interactions. For both the objectives, asymptotic coupled wavenumber expansions for the linear dispersion relation are required for an intermediate fluid loading. The novelty of this work lies in obtaining the asymptotic expansions and using them for predicting the sign change of the nonlinear term at various frequencies. It is found that when the coupled wavenumbers approach the uncoupled pressure-release wavenumbers, the amplitude modulation is stable. On the other hand, near the rigid-duct wavenumbers, the amplitude modulation is unstable. Also, as a further contribution, these wavenumber expansions are used to identify the frequencies of the higher harmonic interactions. And lastly, the solution for the amplitude modulation derived through the MMS is validated using these asymptotic expansions. (C) 2015 Elsevier Ltd. All rights reserved.

An alternative approach to study nonlinear inviscid flow over arbitrary bottom topography

This paper deals with a new approach to study the nonlinear inviscid flow over arbitrary bottom topography. The problem is formulated as a nonlinear boundary value problem which is reduced to a Dirichlet problem using certain transformations. The Dirichlet problem is solved by applying Plemelj-Sokhotski formulae and it is noticed that the solution of the Dirichlet problem depends on the solution of a coupled Fredholm integral equation of the second kind. These integral equations are solved numerically by using a modified method. The free-surface profile which is unknown at the outset is determined. Different kinds of bottom topographies are considered here to study the influence of bottom topography on the free-surface profile. The effects of the Froude number and the arbitrary bottom topography on the free-surface profile are demonstrated in graphical forms for the subcritical flow. Further, the nonlinear results are validated with the results available in the literature and compared with the results obtained by using linear theory. (C) 2015 Elsevier Inc. All rights reserved.

Spectral solutions to the Korteweg-de-Vries and nonlinear Schrodinger equations

In this paper, we present the solutions of 1-D and 2-D non-linear partial differential equations with initial conditions. We approach the solutions in time domain using two methods. We first solve the equations using Fourier spectral approximation in the spatial domain and secondly we compare the results with the approximation in the spatial domain using orthogonal functions such as Legendre or Chebyshev polynomials as their basis functions. The advantages and the applicability of the two different methods for different types of problems are brought out by considering 1-D and 2-D nonlinear partial differential equations namely the Korteweg-de-Vries and nonlinear Schrodinger equation with different potential function. (C) 2015 Elsevier Ltd. All rights reserved.

Weakly nonlinear acoustic wave propagation in a nonlinear orthotropic circular cylindrical waveguide

Nonlinear acoustic wave propagation is considered in an infinite orthotropic thin circular cylindrical waveguide. The modes are non-planar having small but finite amplitude. The fluid is assumed to be ideal and inviscid with no mean flow. The cylindrical waveguide is modeled using the Donnell's nonlinear theory for thin cylindrical shells. The approximate solutions for the acoustic velocity potential are found using the method of multiple scales (MMS) in space and time. The calculations are presented up to the third order of the small parameter. It is found that at some frequencies the amplitude modulation is governed by the Nonlinear Schrodinger Equation (NLSE). The first objective is to study the nonlinear term in the NLSE, as the sign of the nonlinear term determines the stability of the amplitude modulation. On the other hand, at other specific frequencies, interactions occur between the primary wave and its higher harmonics. Here, the objective is to identify the frequencies of the higher harmonic interactions. Lastly, the linear terms in the NLSE obtained using the MMS calculations are validated. All three objectives are met using an asymptotic analysis of the dispersion equation. (C) 2015 Acoustical Society of America.

Weakly corrected numerical solutions to stochastically driven nonlinear dynamical systems

A method to weakly correct the solutions of stochastically driven nonlinear dynamical systems, herein numerically approximated through the Eule-Maruyama (EM) time-marching map, is proposed. An essential feature of the method is a change of measures that aims at rendering the EM-approximated solution measurable with respect to the filtration generated by an appropriately defined error process. Using Ito's formula and adopting a Monte Carlo (MC) setup, it is shown that the correction term may be additively applied to the realizations of the numerically integrated trajectories. Numerical evidence, presently gathered via applications of the proposed method to a few nonlinear mechanical oscillators and a semi-discrete form of a 1-D Burger's equation, lends credence to the remarkably improved numerical accuracy of the corrected solutions even with relatively large time step sizes. (C) 2015 Elsevier Inc. All rights reserved.

Nonlinear Dynamics of a Rotating Flexible Link

This paper deals with the study of the nonlinear dynamics of a rotating flexible link modeled as a one dimensional beam, undergoing large deformation and with geometric nonlinearities. The partial differential equation of motion is discretized using a finite element approach to yield four nonlinear, nonautonomous and coupled ordinary differential equations (ODEs). The equations are nondimensionalized using two characteristic velocities-the speed of sound in the material and a velocity associated with the transverse bending vibration of the beam. The method of multiple scales is used to perform a detailed study of the system. A set of four autonomous equations of the first-order are derived considering primary resonances of the external excitation and one-to-one internal resonances between the natural frequencies of the equations. Numerical simulations show that for certain ranges of values of these characteristic velocities, the slow flow equations can exhibit chaotic motions. The numerical simulations and the results are related to a rotating wind turbine blade and the approach can be used for the study of the nonlinear dynamics of a single link flexible manipulator.

An Ensemble Kushner-Stratonovich-Poisson Filter for Recursive Estimation in Nonlinear Dynamical Systems

We propose a Monte Carlo filter for recursive estimation of diffusive processes that modulate the instantaneous rates of Poisson measurements. A key aspect is the additive update, through a gain-like correction term, empirically approximated from the innovation integral in the time-discretized Kushner-Stratonovich equation. The additive filter-update scheme eliminates the problem of particle collapse encountered in many conventional particle filters. Through a few numerical demonstrations, the versatility of the proposed filter is brought forth.

Kramers Escape Rate in Nonlinear Diffusive Media

In this paper, we study nonlinear Kramers problem by investigating overdamped systems ruled by the one-dimensional nonlinear Fokker-Planck equation. We obtain an analytic expression for the Kramers escape rate under quasistationary conditions by employing

Instability Analysis of Nonlinear Surface Waves in a Circular Cylindrical Container Subjected to a Vertical Excitation

Singular perturbation theory of two-time scale expansions was developed both in inviscid and weak viscous fluids to investigate the motion of single surface standing wave in a liquid-filled circular cylindrical vessel, which is subject to a vertical periodical oscillation. Firstly, it is assumed that the fluid in the circular cylindrical vessel is inviscid, incompressible and the motion is irrotational, a nonlinear evolution equation of slowly varying complex amplitude, which incorporates cubic nonlinear term, external excitation and the influence of surface tension, was derived from solvability condition of high-order approximation. It shows that when forced frequency is low, the effect of surface tension on mode selection of surface wave is not important. However, when forced frequency is high, the influence of surface tension is significant, and can not be neglected. This proved that the surface tension has the function, which causes free surface returning to equilibrium location. Theoretical results much close to experimental results when the surface tension is considered. In fact, the damping will appear in actual physical system due to dissipation of viscosity of fluid. Based upon weakly viscous fluids assumption, the fluid field was divided into an outer potential flow region and an inner boundary layer region. A linear amplitude equation of slowly varying complex amplitude, which incorporates damping term and external excitation, was derived from linearized Navier-Stokes equation. The analytical expression of damping coefficient was determined and the relation between damping and other related parameters (such as viscosity, forced amplitude and depth of fluid) was presented. The nonlinear amplitude equation and a dispersion, which had been derived from the inviscid fluid approximation, were modified by adding linear damping. It was found that the modified results much reasonably close to experimental results. Moreover, the influence both of the surface tension and the weak viscosity on the mode formation was described by comparing theoretical and experimental results. The results show that when the forcing frequency is low, the viscosity of the fluid is prominent for the mode selection. However, when the forcing frequency is high, the surface tension of the fluid is prominent. Finally, instability of the surface wave is analyzed and properties of the solutions of the modified amplitude equation are determined together with phase-plane trajectories. A necessary condition of forming stable surface wave is obtained and unstable regions are illustrated. (c) 2005 Elsevier SAS. All rights reserved.

The Nonlinear Phenomena Of Thin Polydimethylsiloxane (Pdms) Films In Electrowetting

Electrowetting (EW) is an effective way to manipulate small volume liquid in micro- and nano-devices, for it can improve its wettability. Since the late 1990s, electrowetting-on-dielectric (EWOD) has been used widely in bio-MEMS, lab-on-a-chip, etc. Polydimethlsiloxane (PDMS) is extensively utilized as base materials in the fabrication of biomedical micro- and nano-devices. The properties of thin PDMS films used as dielectric layer in EW are studied in this paper. The experimental results show that the thin PDMS films exhibit good properties in EWOD. As to PDMS films with different thicknesses, a threshold voltage and a hysteresis were observed in the EIWOD experiments.

Nonlinear Constitutive Behavior Of Ferroelectric Materials

The ferroelectric specimen is considered as an aggregation of many randomly oriented domains. According to this mechanism, a multi-domain mechanical model is developed in this paper. Each domain is represented by one element. The applied stress and electric field are taken to be the stress and electric field in the formula of the driving force of domain switching for each element in the specimen. It means that the macroscopic switching criterion is used for calculating the volume fraction of domain switching for each element. By using the hardening relation between the driving force of domain switching and the volume fraction of domain switching calibrated, the volume fraction of domain switching for each element is calculated. Substituting the stress and electric field and the volume fraction of domain switching into the constitutive equation of ferroelectric material, one can easily get the strain and electric displacement for each element. The macroscopic behavior of the ferroelectric specimen is then directly calculated by volume averaging. Meanwhile, the nonlinear finite element analysis for the ferroelectric specimen is carried out. In the finite element simulation, the volume fraction of domain switching for each element is calculated by using the same method mentioned above. The interaction between different elements is taken into account in the finite element simulation and the local stress and electric field for each element is obtained. The macroscopic behavior of the specimen is then calculated by volume averaging. The computation results involve the electric butterfly shaped curves of axial strain versus the axial electric field and the hysteresis loops of electric displacement versus the electric field for ferroelectric specimens under the uniaxial coupled stress and electric field loading. The present theoretical prediction agrees reasonably with the experimental results.

Energy Principle And Nonlinear Electric-Mechanical Behavior Of Ferroelectric Ceramics

Many experimental observations have shown that a single domain in a ferroelectric material switches by progressive movement of domain walls, driven by a combination of electric field and stress. The mechanism of the domain switch involves the following steps: initially, the domain has a uniform spontaneous polarization; new domains with the reverse polarization direction nucleate, mainly at the surface, and grow though the crystal thickness; the new domain expands sideways as a new domain continues to form; finally, the domain switch coalesces to complete the polarization reversal. According to this mechanism, the volume fraction of the domain switching is introduced in the constitutive law of the ferroelectric material and used to study the nonlinear constitutive behavior of a ferroelectric body in this paper. The principle of stationary total potential energy is put forward in which the basic unknown quantities are the displacement u(i), electric displacement D-i and volume fraction rho(I) of the domain switching for the variant I. The mechanical field equation and a new domain switching criterion are obtained from the principle of stationary total potential energy. The domain switching criterion proposed in this paper is an expansion and development of the energy criterion established by Hwang et al. [ 1]. Based on the domain switching criterion, a set of linear algebraic equations for determining the volume fraction rho(I) of domain switching is obtained, in which the coefficients of the linear algebraic equations only contain the unknown strain and electric fields. If the volume fraction rho(I) of domain switching for each domain is prescribed, the unknown displacement and electric potential can be obtained based on the conventional finite element procedure. It is assumed that a domain switches if the reduction in potential energy exceeds a critical energy barrier. According to the experimental results, the energy barrier will strengthen when the volume fraction of the domain switching increases. The external mechanical and electric loads are increased step by step. The volume fraction rho(I) of domain switching for each element obtained from the last loading step is used as input to the constitutive equations. Then the strain and electric fields are calculated based on the conventional finite element procedure. The finite element analysis is carried out on the specimens subjected to uniaxial coupling stress and electric field. Numerical results and available experimental data are compared and discussed. The present theoretic prediction agrees reasonably with the experimental results.

Derivation Of A Nonlinear Reynolds Stress Model Using Renormalization Group Analysis And Two-Scale Expansion Technique

Adopting Yoshizawa's two-scale expansion technique, the fluctuating field is expanded around the isotropic field. The renormalization group method is applied for calculating the covariance of the fluctuating field at the lower order expansion. A nonlinear Reynolds stress model is derived and the turbulent constants inside are evaluated analytically. Compared with the two-scale direct interaction approximation analysis for turbulent shear flows proposed by Yoshizawa, the calculation is much more simple. The analytical model presented here is close to the Speziale model, which is widely applied in the numerical simulations for the complex turbulent flows.

Dynamic analysis of arrest of buckle propagation on a beam on a nonlinear elastic foundation by FEM

Based on the dynamic governing equation of propagating buckle on a beam on a nonlinear elastic foundation, this paper deals with an important problem of buckle arrest by combining the FEM with a time integration technique. A new conclusion completely different from that by the quasi-static analysis about the buckle arrestor design is drawn. This shows that the inertia of the beam cannot be ignored in the analysis under consideration, especially when the buckle propagation is suddenly stopped by the arrestors.