Andersen, Peter Risby

Cutanda Henriquez, Vicente

Sánchez-Dehesa, José

Online Research Database In Technology

                       General rights Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights.  • Users may download and print one copy of any publication from the public portal for the purpose of private study or research. • You may not further distribute the material or use it for any profit-making activity or commercial gain • You may freely distribute the URL identifying the publication in the public portal   If you believe that this document breaches copyright please contact us providing details, and we will remove access to the work immediately and investigate your claim.     Downloaded from orbit.dtu.dk on: Nov 09, 2017Numerical modeling of double-negative acoustic metamaterials: Should losses beincluded?Cutanda Henriquez, Vicente; Andersen, Peter Risby; Sánchez-Dehesa, JoséPublication date:2016Document VersionPublisher's PDF, also known as Version of recordLink back to DTU OrbitCitation (APA):Cutanda Henriquez, V., Andersen, P. R., & Sánchez-Dehesa, J. (2016). Numerical modeling of double-negativeacoustic metamaterials: Should losses be included?. Abstract from EuroRegio2016, Porto, Portugal.   Numerical modeling of double-negative acoustic metamaterials: Should losses be included? Vicente Cutanda Henríquez1, Peter Risby Andersen1, José Sánchez-Dehesa2 1Centre for Acoustic-Mechanical Micro Systems Technical University of Denmark Kgs. Lyngby, Denmark {vcuhe@elektro.dtu.dk, prand@elektro.dtu.dk)  2Wave Phenomena Group Department of Electronic Engineering Universitat Politècnica de València, Camino de Vera s.n. (building 7F), ES-46022 Valencia, Spain (jsdehesa@upv.es)  Abstract Viscous and thermal losses of acoustic waves are usually neglected or accounted for as boundary impedance. It is known, however, that acoustic losses become relevant in devices with some dimension in the millimeter range or below. On the other hand, the new class of structures called acoustic metamaterials can be affected by acoustic losses, but in this case the extension of these effects is less known. Acoustic metamaterials are intricate periodic structures where, at frequencies low enough (corresponding to wavelengths much larger than the structure period), elementary units interact producing interesting unusual effects.  In this paper advanced modeling tools based on the Boundary Element Method (BEM) and the Finite Element Method (FEM) are used to study the effect of losses in an acoustic metamaterial scaled to different sizes. The conclusions are expected to give insight on the practical limitations when using acoustic metamaterials. Keywords: acoustic metamaterials, numerical methods, visco-thermal losses. PACS no. 43.20.Hq, 43.20.Rz 1 Introduction Viscous and thermal losses are present in any acoustic wave and show as attenuation over large distances or as a boundary effect. Usually the lossless wave equation is employed to describe the sound field, and loss is assigned to attenuation terms or surface impedances.[1]   However, when the dimensions of the acoustic domain are of a similar scale as the viscous and thermal boundary layers, neglecting losses can lead to unrealistic results. These dimensions are in the millimeter/micrometer range for audible frequencies. Well-known examples are acoustic couplers, acoustic transducers, or hearing aids.[2,3]  EuroRegio2016, June 13-15, Porto, Portugal     2  The periodic arrangement of specially designed structural units known as acoustic metamaterial can show unusual properties at low frequencies which can be employed to make feasible theoretical proposals like acoustic cloaking.[4,5] Acoustic metamaterials, where the acoustic wave propagates in a fluid rather than a solid, have recently joined the list of structures where acoustic fluid losses might be relevant, with the added difficulty of their necessarily intricate internal structure. The dimensions of these materials are usually designed larger than the boundary layers but much smaller than the wavelength, and therefore it is not very clear in what circumstances and to what extent viscous and thermal losses affect their very special properties.  On the other hand, new numerical modeling tools have been developed in the last 10-15 years, which can cope with viscous and thermal losses, and employ either the Finite Element Method (FEM) or the Boundary Element Method (BEM).[6,7,8] They are however much more computationally costly than their corresponding lossless implementations, and therefore are often replaced by lossless implementations when evaluating metamaterials.  The aim of this paper is evaluating how losses can affect an acoustic metamaterial at different length scales. To do this, a particular acoustic metamaterial in the literature has been chosen.[9] The metamaterial is scaled up and the effect of losses in its special properties studied.  Section 2 gives an overview of the numerical methods with losses used in the calculation. Section 3 briefly describes the acoustic metamaterial test case, taken from [9]. The scaling method and its results are shown in section 4. Section 5 is dedicated to conclusions and future work. 2 Numerical models with losses Full numerical implementations of viscous and thermal losses are based on the linearized Navier-Stokes equations:[1]  ( ) ( )( )TpcpTTCsvvptvTtsTvtopooooβγρβγγµµηρλρρρ−= −−=×∇×∇−∇∇ ++−∇=∂∂∆=∂∂=∇+∂∂21340CCCC (1a-1e)  Equations (1a) to (1e) represent respectively conservation of mass, conservation of energy, conservation of momentum and thermodynamic conditions. The symbols are: ρ, density, ρ0, static density, vC, particle velocity, T, temperature, T0, static temperature, s, entropy, λ, thermal conductivity, p, sound pressure, η, bulk viscosity, μ, coefficient of viscosity, Cp, specific heat at constant pressure, γ, ratio of specific heats and β,  rate of increase of pressure with temperature at constant volume. In addition to linear variations, no flow, homogeneous fluid and dimensions of the setup and wavelength larger than the molecular mean free path (~10-7 m) are assumed.  2.1 Finite Element Method The direct implementation of the linearized Navier-Stokes equations (1a-1e) in the Finite Element Method has been proposed and developed by several authors.[5,9] The implementation was made  EuroRegio2016, June 13-15, Porto, Portugal     3 originally in the open FEM package Elmer and it has been included later in the Comsol FEM package, which has been employed in this paper in its version 5.2.[10]  This is a computationally intensive implementation, mainly because the thin viscous and thermal boundary layers need to be properly represented. For this reason, a sparing but sufficient meshing is important, requiring an understanding of every case and how much detail is required.  FEM with losses is the most extended implementation. It was included early in a commercial package, and its implementation has been progressively refined and optimized. A growing number of examples including losses using FEM can be found in the literature and in the Comsol webpage.  2.2 Boundary Element Method The BEM implementation with losses is based on the Kirchhoff derivation of the Navier-Stokes equations (1a-1e):[1,11]   ( )( )( ) 0w000222=⋅∇=+∆=+∆=+∆vvvhhaavithvkpkpkCCC (2a-2c)     Eqs. (2a-2c) represent, the so-called acoustic (a), thermal (h) and viscous (v) modes, which can be dealt with independently in the acoustic domain. They are however linked through the temperature and velocity boundary conditions. The pressure is separated into two components p = pa + ph, while the velocity is separated into three components vha vvvvCCCC++= . The viscous velocity or rotational velocity does not have a corresponding pressure. Correspondingly, the wavenumbers ka, kh and kv are based on the lossless wavenumber k and physical properties of the fluid, such as the viscosity, bulk viscosity and thermal conductivity coefficients, air density, and specific heats. Eq. (2a) is a wave equation, while Eqs. (2b-2c) are diffusion equations, given the large imaginary part of kh and kv. Eq. (2c) is a vector equation and therefore can be split into three components, giving a total of five unknowns.  The implementation in BEM is made by discretizing equations (2a-2c) one by one and combining them into a single matrix equation using the boundary conditions. The matrix equation is solved for the acoustic pressure pa and subsequently other variables are obtained on the boundary. From the boundary values, any domain field point can be calculated. [7,8] 3 Description of the acoustic metamaterial The results shown in this paper are based on the acoustic metamaterial presented in reference [9], where details can be found. In short, the authors of [9] proposed a metamaterial that is formed by elementary units as the one shown in Fig. 1. The units are arranged in a periodic lattice, and the sound waves propagate in between two planes separated by the distance h, where the units are embedded. Such a metamaterial shows extraordinary properties within a certain  frequency band, where double negative behaviour is predicted; i.e., where both the dynamical  bulk modulus and mass density have negative values. In addition, at frequencies where the density is near zero, interesting behaviours like tunnelling of acoustic waves through narrow channels, restoring of a plane wave after scattering, and perfect transmission through bends and corners were also theoretically characterized.   EuroRegio2016, June 13-15, Porto, Portugal     4                 Figure 1 – A single unit of the example metamaterial from reference [9]. Left: schematic view with some of its defining parameters (radius Rb, height of the space between planes h, and depth of the wells, L). Right: BEM mesh version of the same unit, where only the domain boundary is meshed. The positions of the planes where the unit is embedded are shown as semi-transparent rectangles.                           Figure 2 – a) Sample of acoustic metamaterial measured in Ref. [9]; b) corresponding BEM mesh. (a) (b)  EuroRegio2016, June 13-15, Porto, Portugal     5 The results in Ref. [9] are based on an analytical model of the material and on lossless FEM simulations, both predicting the double-negative behavior. The paper also includes measurements made on a 3D printed sample, as shown in Fig. 2a. Unfortunately, the expected double-negative behavior were not supported by experimental data.    Figure 3 – Transmittance results for the metamaterial sample (sample A in Ref. [9]). Solid blue curve, lossless FEM result; dashed light blue curve, FEM with losses; dash-dotted red curve, measurements (from [9]); green-marker curve, BEM with losses. Double-negative frequency is ~2400 Hz.   In order to test the performance of the material, the transmittance (insertion loss) is simulated using the numerical methods with losses described in the previous section, and compared with the results in [9]. The transmittance is defined as the ratio of the acoustic power after an incident plane wave has traversed the material, with respect to the incident power. It should vary between 0 (no transmission) to 1 (full transmission). The transmittance results are shown in Fig. 3.   The double-negative behavior is predicted by the lossless FEM model at around 2400-2500 Hz. Lossless BEM results (not shown), also predict this behavior at a slightly higher frequency. However, measurements and models with viscous and thermal losses lack the transmittance peak at this frequency: the extraordinary behavior seems to be hindered by losses. This result was already deduced in [9] from the measurements, and it is confirmed by the numerical models. 4 Scaling of the metamaterial and results The metamaterial described in the previous section has a length (from the input to the output) of 127,5 mm. In this section models of scaled-up versions of this setup are used. The lossless description of the metamaterial should give, and it actually does for the lossless models, the same transmittance results at  EuroRegio2016, June 13-15, Porto, Portugal     6 a corresponding scaled-down frequency. Viscous and thermal losses, however, do not vary with scale in the same way.[1] The thickness of the viscous and thermal boundary layers is inversely proportional to f . Moreover, viscous losses strongly depend on the incidence angle. The scaled models are computationally demanding, but no more than the original setup. The physical dimensions are made larger, but the calculation frequencies are divided by the same factors. This means the same mesh can be reused. The transmittance is calculated with the same piston velocity excitation in all cases, which gives a growing incident power, since the section of the sample grows as the square of the scale factor.   Fig. 4 shows the transmittance at the double-negative frequency for different scale factors. The frequency at which each point is calculated is the double-negative frequency (2600 Hz in the BEM model), divided by the scale factor. The lossless calculation gives a fairly good transmission as expected, but in the models with losses the transmission improves very slowly as the scale grows. Even for the larger scale, where the metamaterial is about 2,5 m long, the transmittance is below 10-2. The FEM model shows a similar tendency, despite its evaluation of an overall lower transmittance. Curves also show that the improvement of the transmittance is less significant as the scale grows, stabilizing at a rather low value.   Figure 4 – Transmittance at the scaled double-negative frequency as a function of the scaling factor. The lossless calculation (blue '+') shows high transmission, while BEM (red 'x') and FEM (green 'o') with losses predict a quite poor transmission. Large versions of the metamaterial hardly improve. The lower figure is equivalent to the upper one, but using a logarithmic ordinate axis.  EuroRegio2016, June 13-15, Porto, Portugal     7 5 Conclusions and future work The results in section 4 seem to contradict the general belief assuming that a sufficiently large acoustic metamaterial will not be affected by viscous and thermal losses. On the contrary, results indicate that the example metamaterial cannot be made large enough so that losses become negligible. Even for large versions of the metamaterial, most of the incident power never reaches the output of the material.  This conclusion is the same whether the material is modeled using FEM or BEM with losses. These methods are based on the same physical background but have very different implementations. No measurements have been done on other than the original setup.  In the future, more work can be done in order to support the main conclusion and understand better the effect of losses on acoustic metamaterials. Besides performing verification measurements, other metamaterial setups may be tried. A detailed analysis of the simulation results can also reveal more details about what is the contribution of the different parts of the setup to the overall losses, and how the metamaterial design can minimize them. References  [1] Pierce, A.D., Acoustics. An introduction to its physical principles and applications, McGraw Hill, New York, 1981. [2] Cutanda Henriquez, V. and Juhl, P. M., Modelling measurement microphones using BEM with visco-thermal losses, Proc. Joint Baltic-Nordic Acoustics Meeting, 18-20 June 2012, Odense, Denmark. [3] Homentcovschi, D. and Miles, R. N., An analytical-numerical method for determining the mechanical response of a condenser microphone, J. Acoust. Soc. of Am., 130 (6), 3698–3705, 2011. [4] Torrent, D and Sánchez-Dehesa, J., Acoustic cloaking in two-dimension: a feasible approach, New J. Phys., 10, 063015, 2008. [5] Craster, R.V.; Guenneau, S., Eds., Acoustic Metamaterials: Negative Refraction, Lensing and Cloaking, Springer, London, 2013. [6] Malinen, M., Lyly, M., Råback, P., Kärkkäinen, A. and Kärkkäinen, L., A finite element method for the modeling of thermo-viscous effects in acoustics. Proc. 4th European Congress Computational Methods in Applied Sciences and Engineering ECCOMAS, Jyväskylä, Finland, 1-12, 2004. [7] Cutanda Henriquez, V. and Juhl, P. M., An axisymmetric boundary element formulation of sound wave propagation in fluids including viscous and thermal losses. J. Acoust. Soc. Am., 134(5), 3409-3418, 2013. [8] Cutanda Henriquez, V. and Juhl, P. M., Implementation of an Acoustic 3D BEM with Visco-Thermal Losses. Internoise 2013, 15-18 Sept, 2013, Innsbruck (Austria). [9] Graciá-Salgado, R., García-Chocano, V. M., Torrent D. and Sánchez-Dehesa, J., Negative mass density and ρ-near-zero quasi-two-dimensional metamaterials: Design and applications, Phys. Rev. B 88, 224305 (2013).  EuroRegio2016, June 13-15, Porto, Portugal     8 [10] Joly, N., Coupled equations for particle velocity and temperature variation as the fundamental formulation of linear acoustics in thermo-viscous fluids at rest, Acta Acust. Acust. 92, 202–209, 2006. [11] www.comsol.com. Consulted April 2016. [12] Bruneau, M., Herzog, Ph., Kergomard J. and Polack, J. D., General formulation of the dispersion equation in bounded visco-thermal fluid, and application to some simple geometries. Wave Motion 11,  441-451, 1989. 

Numerical modeling of double-negative acoustic metamaterials: Should losses be included?

http://orbit.dtu.dk/ws/files/124937010/EuroRegio2016_VCH_PRA_JSD.pdf

Numerical modeling of double-negative acoustic metamaterials: Should losses be included?

Abstract

Similar works

Full text

Available Versions

Online Research Database In Technology