No Arabic abstract
Acoustic wave propagation in a one-dimensional waveguide connected with Helmholtz resonators is studied numerically. Finite amplitude waves and viscous boundary layers are considered. The model consists of two coupled evolution equations: a nonlinear PDE describing nonlinear acoustic waves, and a linear ODE describing the oscillations in the Helmholtz resonators. The thermal and viscous losses in the tube and in the necks of the resonators are modeled by fractional derivatives. A diffusive representation is followed: the convolution kernels are replaced by a finite number of memory variables that satisfy local ordinary differential equations. A splitting method is then applied to the evolution equations: their propagative part is solved using a standard TVD scheme for hyperbolic equations, whereas their diffusive part is solved exactly. Various strategies are examined to compute the coefficients of the diffusive representation; finally, an optimization method is preferred to the usual quadrature rules. The numerical model is validated by comparisons with exact solutions. The properties of the full nonlinear solutions are investigated numerically. In particular, the existence of acoustic solitary waves is confirmed.
Heterogeneity can be accounted for by a random potential in the wave equation. For acoustic waves in a fluid with fluctuations of both density and compressibility (as well as for electromagnetic waves in a medium with fluctuation of both permittivity and permeability) the random potential entails a scalar and an operator contribution. For simplicity, the latter is usually overlooked in multiple scattering theory: whatever the type of waves, this simplification amounts to considering the Helmholtz equation with a sound speed $c$ depending on position $mathbf{r}$. In this work, a radiative transfer equation is derived from the wave equation, in order to study energy transport through a multiple scattering medium. In particular, the influence of the operator term on various transport parameters is studied, based on the diagrammatic approach of multiple scattering. Analytical results are obtained for fundamental quantities of transport theory such as the transport mean-free path $ell^*$, scattering phase function $f$ and anisotropy factor $g$. Discarding the operator term in the wave equation is shown to have a significant impact on $f$ and $g$, yet limited to the low-frequency regime i.e., when the correlation length of the disorder $ell_c$ is smaller than or comparable to the wavelength $lambda$. More surprisingly, discarding the operator part has a significant impact on the transport mean-free path $ell^*$ whatever the frequency regime. When the scalar and operator terms have identical amplitudes, the discrepancy on the transport mean-free path is around $300,%$ in the low-frequency regime, and still above $30,%$ for $ell_c/lambda=10^3$ no matter how weak fluctuations of the disorder are. Analytical results are supported by numerical simulations of the wave equation and Monte Carlo simulations.
The rising need for hybrid physical platforms has triggered a renewed interest for the development of agile radio-frequency phononic circuits with complex functionalities. The combination of travelling waves with resonant mechanical elements appears as an appealing means of harnessing elastic vibration. In this work, we demonstrate that this combination can be further enriched by the occurrence of elastic non-linearities induced travelling surface acoustic waves (SAW) interacting with a pair of otherwise linear micron-scale mechanical resonators. Reducing the resonator gap distance and increasing the SAW amplitude results in a frequency softening of the resonator pair response that lies outside the usual picture of geometrical Duffing non-linearities. The dynamics of the SAW excitation scheme allows further control of the resonator motion, notably leading to circular polarization states. These results paves the way towards versatile high-frequency phononic-MEMS/NEMS circuits fitting both classical and quantum technologies.
We investigate flow of incompressible fluid in a cylindrical tube with elastic walls. The radius of the tube may change along its length. The discussed problem is connected to the blood flow in large human arteries and especially to nonlinear wave propagation due to the pulsations of the heart. The long-wave approximation for modeling of waves in blood is applied. The obtained model Korteweg-deVries equation possessing a variable coefficient is reduced to a nonlinear dynamical system of 3 first order differential equations. The low probability of arising of a solitary wave is shown. Periodic wave solutions of the model system of equations are studied and it is shown that the waves that are consequence of the irregular heart pulsations may be modeled by a sequence of parts of such periodic wave solutions.
Assigning homogeneous boundary conditions, such as acoustic impedance, to the thermoviscous wave equations (TWE) derived by transforming the linearized Navier-Stokes equations (LNSE) to the frequency domain yields a so-called Helmholtz solver, whose output is a discrete set of complex eigenfunction and eigenvalue pairs. The proposed method -- the inverse Helmholtz solver (iHS) -- reverses such procedure by returning the value of acoustic impedance at one or more unknown impedance boundaries (IBs) of a given domain via spatial integration of the TWE for a given real-valued frequency with assigned conditions on other boundaries. The iHS procedure is applied to a second-order spatial discretization of the TWEs derived on an unstructured grid with staggered grid arrangement. The momentum equation only is extended to the center of each IB face where pressure and velocity components are co-located and treated as unknowns. One closure condition considered for the iHS is the assignment of the surface gradient of pressure phase over the IBs, corresponding to assigning the shape of the acoustic waveform at the IB. The iHS procedure is carried out independently for each frequency in order to return the complete broadband complex impedance distribution at the IBs in any desired frequency range. The iHS approach is first validated against Rotts theory for both inviscid and viscous, rectangular and circular ducts. The impedance of a geometrically complex toy cavity is then reconstructed and verified against companion full compressible unstructured Navier-Stokes simulations resolving the cavity geometry and one-dimensional impedance test tube calculations based on time-domain impedance boundary conditions (TDIBC). The iHS methodology is also shown to capture thermoacoustic effects, with reconstructed impedance values quantitatively in agreement with thermoacoustic growth rates.
We begin with the theoretical study of spectral energy cascade due to the propagation of high amplitude sound in the absence of thermal sources. To this end, a first-principles-based system of governing equations, correct up to second order in perturbation variables is derived. The exact energy corollary of such second-order system of equations is then formulated and used to elucidate the spectral energy dynamics of nonlinear acoustic waves. We then extend this analysis to thermoacoustically unstable waves -- i.e. amplified as a result of thermoacoustic instability. We drive such instability up until the generation of shock waves. We further study the nonlinear wave propagation in geometrically complex case of waves induced by the spark plasma between the electrodes. This case adds the geometrical complexity of a curved, three-dimensional shock, yielding vorticity production due to baroclinic torque. Finally, detonation waves are simulated by using a low-order approach, in a periodic setup subjected to high pressure inlet and exhaust of combustible gaseous mixture. An order adaptive fully compressible and unstructured Navier Stokes solver is currently under development to enable higher fidelity studies of both the spark plasma and detonation wave problem in the future.