Propagation and attenuation of sound through a layered phononic crystal with viscous constituents is theoretically studied. The Navier–Stokes equation with appropriate boundary conditions is solved and the dispersion relation for sound is obtained for a periodic layered heterogeneous structure where at least one of the constituents is a viscous fluid. Simplified dispersion equations are obtained when the other component of the unit is either elastic solid, viscous fluid, or ideal fluid. The limit of low frequencies when periodic structure homogenizes and the frequencies close to the band edge when propagating Bloch wave becomes a standing wave are considered and enhanced viscous dissipation is calculated. Angular dependence of the attenuation coefficient is analyzed. It is shown that transition from dissipation in the bulk to dissipation in a narrow boundary layer occurs in the region of angles close to normal incidence. Enormously high dissipation is predicted for solid–fluid structure in the region of angles where transmission practically vanishes due to appearance of so-called “transmission zeros,” according to El Hassouani, El Boudouti, Djafari-Rouhani, and Aynaou [Phys. Rev. B 78, 174306 (2008)]. For the case when the unit cell contains a narrow layer of high viscosity fluid, the anomaly related to acoustic manifestation of Borrmann effect is explained.
more »
« less
Finite boundary effects on the spherical Rayleigh–Taylor instability between viscous fluids
For the Rayleigh–Taylor unstable arrangement of a viscous fluid sphere embedded in a finite viscous fluid spherical shell with a rigid boundary and a radially directed acceleration, a dispersion relation is developed from a linear stability analysis using the method of normal modes. aR1 is the radially directed acceleration at the interface. ρi denotes the density, μi is the viscosity, and Ri is the radius, where i = 1 is the inner sphere and i = 2 is the outer sphere. The dispersion relation is a function of the following dimensionless variables: viscosity ratio s=μ1μ2, density ratio d=ρ1ρ2, spherical harmonic mode n, B=R1aR1ρ22μ221/3, H=R2R1, and the dimensionless growth rate α=σμ2aR12ρ21/3, where σ is the exponential growth rate. We show that the boundedness provided by the outer spherical shell has a strong influence on the instability behavior, which is reflected not only in the modulation of the growth rate but also in the selection of the most unstable modes that are physically possible. This outer boundary effect is quantified by the relative magnitude of the radius ratio H. We find that when H is close to unity, lower order harmonics are excluded from becoming the most unstable within a vast region of the parameter space. In other words, the effect of H has precedence over the other controlling parameters d, B, and a wide range of s in establishing what the lowest most unstable mode can be. When H ∼ 1, low order harmonics can become the most unstable only for s ≫ 1. However, in the limit when s → ∞, we show that the most unstable mode is n = 1 and derive the dispersion relation in this limit. The exclusion of most unstable low order harmonics caused by a finite outer boundary is not realized when the outer boundary extends beyond a certain threshold length-scale in which case all modes are equally possible depending on the value of B.
more »
« less
- Award ID(s):
- 1748883
- PAR ID:
- 10597334
- Publisher / Repository:
- American Institute of Physics
- Date Published:
- Journal Name:
- AIP Advances
- Volume:
- 12
- Issue:
- 4
- ISSN:
- 2158-3226
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Filamentary structures in neutral hydrogen (Hi) emission are well aligned with the interstellar magnetic field, so Hiemission morphology can be used to construct templates that strongly correlate with measurements of polarized thermal dust emission. We explore how the quantification of filament morphology affects this correlation. We introduce a new implementation of the Rolling Hough Transform (RHT) using spherical harmonic convolutions, which enables efficient quantification of filamentary structure on the sphere. We use this Spherical RHT algorithm along with a Hessian-based method to construct Hi-based polarization templates. We discuss improvements to each algorithm relative to similar implementations in the literature and compare their outputs. By exploring the parameter space of filament morphologies with the Spherical RHT, we find that the most informative Histructures for modeling the magnetic field structure are the thinnest resolved filaments. For this reason, we find a ∼10% enhancement in theB-mode correlation with polarized dust emission with higher-resolution Hiobservations. We demonstrate that certain interstellar morphologies can produce parity-violating signatures, i.e., nonzeroTBandEB, even under the assumption that filaments are locally aligned with the magnetic field. Finally, we demonstrate thatBmodes from interstellar dust filaments are mostly affected by the topology of the filaments with respect to one another and their relative polarized intensities, whereasEmodes are mostly sensitive to the shapes of individual filaments.more » « less
-
We analyse the results of direct numerical simulations of rotating convection in spherical shell geometries with stress-free boundary conditions, which develop strong zonal flows. Both the Ekman number and the Rayleigh number are varied. We find that the asymptotic theory for rapidly rotating convection can be used to predict the Ekman number dependence of each term in the governing equations, along with the convective flow speeds and the dominant length scales. Using a balance between the Reynolds stress and the viscous stress, together with the asymptotic scaling for the convective velocity, we derive an asymptotic prediction for the scaling behaviour of the zonal flow with respect to the Ekman number, which is supported by the numerical simulations. We do not find evidence of distinct asymptotic scalings for the buoyancy and viscous forces and, in agreement with previous results from asymptotic plane layer models, we find that the ratio of the viscous force to the buoyancy force increases with Rayleigh number. Thus, viscosity remains non-negligible and we do not observe a trend towards a diffusion-free scaling behaviour within the rapidly rotating regime.more » « less
-
We study the periodic buckling patterns that emerge when elastic shells are subjected to geometric confinement. Residual swelling provides access to range of shapes (saddles, rolled sheets, cylinders, and spherical sections) which vary in their extrinsic and intrinsic curvatures. Our experimental and numerical data show that when these moderately thick structures are radially confined, a single geometric parameter – the ratio of the total shell radius to the amount of unconstrained material – predicts the number of lobes formed. We present a model that interprets this scaling as the competition between radial and circumferential bending. Next, we show that reducing the transverse confinement of saddles causes the lobe number to decrease with a similar scaling analysis. Hence, one geometric parameter captures the wave number through a wide range of radial and transverse confinement, connecting the shell shape to the shape of the boundary that confines it. We expect these results to be relevant for an expanse of shell shapes, and thus applicable to the design of shape-shifting materials and the swelling and growth of soft structures.more » « less
-
For plane channel flow, thermal stratification resulting from a wall-normal temperature gradient together with an opposing gravitational field can lead to buoyancy-driven instability of three-dimensional waves. Moreover, viscosity-driven instability can lead to the amplification of two-dimensional Tollmien-Schlichting waves. Temporal stability simulations considering different combinations of Reynolds number and Rayleigh number were performed to investigate both the buoyancy and viscosity-driven instability of Rayleigh-Benard-Poiseuille flow. The investigated cases are either (1) stable, (2) unstable with respect to three-dimensional waves (buoyancy-driven instability), or (3) unstable with respect to two-dimensional waves (viscosity-driven instability). Two new and highly accurate computational fluid dynamics codes have been developed for solving the full and linearized unsteady compressible Navier-Stokes equations in Cartesian coordinates. The codes employ fifth-order-accurate upwind-biased compact finite differences for the convective terms and fourth-order-accurate compact finite differences for the viscous terms. For the case with buoyancy-driven instability, strong linear growth is observed for a broad range of spanwise wavenumbers and the wavelength of the spanwise mode with the strongest non-linear growth is gradually decreasing in time. For the case with viscosity-driven instability, the linear growth rates are lower and the first mode to experience non-linear growth is a higher harmonic with half the wavelength of the primary wave. The present results are consistent with the neutral curves from the linear stability theory analysis by Gage and Reid.more » « less
An official website of the United States government
