Issue 
Acta Acust.
Volume 7, 2023



Article Number  30  
Number of page(s)  11  
Section  Ultrasonics  
DOI  https://doi.org/10.1051/aacus/2023024  
Published online  26 June 2023 
Scientific Article
Numerical study of beam shape adaptation by anisotropic disk covering transducer or metamaterial
Georgia Institute of Technology, George W. Woodruff School of Mechanical Engineering, Georgia Tech Europe, Laboratory for Ultrasonic Nondestructive Evaluation, IRL 2958 Georgia Tech – CNRS, 2 rue Marconi, 57070 Metz, France
^{*} Corresponding author: declercq@gatech.edu
Received:
20
January
2023
Accepted:
22
May
2023
Metamaterials are intensely explored for their capabilities to modify sound beams. In addition to frequency filtering, acoustic lenses offer intriguing possibilities for shaping sound beams. For the time being, the versatility of metamaterials remains limitless. In beamshape adaptation, however, their complexity suggests that manufacturers of transducers could benefit from combining metamaterials with more conventional materials. This paper investigates the transmission of a circumscribed beam through a stratum of anisotropic material to examine the change in beam shape after transmission. The incident sound is presumed to originate from a conventional transducer, possibly coated with a metamaterial to modify the sound field, before being transmitted through the anisotropic layer. Different incident beam shapes, such as conicallike, Gaussian, and pillar beams, are investigated. While the results are not exhaustive, they demonstrate the beam shape’s adaptability.
Key words: Beam shape adaptation / Gaussian beams / Pillar beams / Conicallike beams / Fiberreinforced composites / Metamaterials
© The Author(s), Published by EDP Sciences, 2023
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Recently, it has been demonstrated that metasurfaces are efficient and compact structures for designing arbitrary wavefronts. Metasurfaces are planar metamaterials with a subwavelength thickness that enables wavefront sculpting by introducing gradients in the spatial wave response of these flat structures. Impedancebased holographic acoustic lenses have been developed [1] to transform the sound output from a transducer into the desired field for medical applications [2]. Bessel beams, known for their diffraction resilience, can be formed equally by a metamaterial placed before a transducer [3]. The focus can be adjusted by shifting one material to another, covering the transducer surface [4]. Acoustic gradientindex (GRIN) metasurfaces, engineered from soft gradedporous silicone rubber, also permit beam steering and focusing [5]. Acoustic metamaterials are also developed for making a Mikaelian lens [6] or coding acoustic waves by sending them through the metamaterial before transmission [7]. Certain metamaterials achieve anisotropic behavior but not in the homogenized sense of actual crystals or fiberreinforced composites with constituents way below the involved acoustic wavelengths [8]. Embedded inplate applications also exist in which metamaterials steer or focus beams [9].
The versatility of metamaterials remains limitless for the time being. In beamshape modifications, however, the degree of complexity implies that transducer manufacturers might gain from merging metamaterials with traditional substances, if possible. Those substances can be crystals or fiberreinforced composites, anisotropic but homogeneous for the involved sound waves. Emitting sound through such conventional materials can further adjust beams without sophistication, and, therefore, find applications in the abovementioned topics of beamshape, steering and focus adjustments in the medical field.
It is known that sound transmits through crystals in a directionaldependent manner. However, the effect on beam shape, merely due to anisotropy and not due to impedance inhomogeneity, usually is not considered. Here, we investigate the transmission of acoustic beams perpendicularly incident on an anisotropic plate. The 3D beams are modeled by a plane wave expansion based on a Fourier series approach whereby each constituting plane wave interacts with the plate under consideration.
That sound is sensitive to the anisotropic linear stiffness of composites and has been widely exploited in the past. Indeed, ultrasonic investigation of fiberreinforced composites’ stiffnesses is based on straightforward principles involving Snell’s law, Newton’s laws, Hooke’s law, and the continuity of mechanical strains and stresses at the involved interfaces between composite layers or the surrounding water. Linear stiffness determination is practically achieved by exploiting direction dependence in experiments. It is, therefore, somewhat different from damage detection using Cscans or guided waves [10–13] or nonlinear acoustics [14]. Breakthroughs came from the team of Dale Chimenti [15–20] and others [21–23]. A direct method to probe the stiffness of composites is based on multipleangle incidence, done in the form of a transducer rotation along a sphere [24] or through focused sound with the help of phased array technology [25, 26] or otherwise [19, 27–29]. A significant challenge is to work on an aircoupled system [30, 31], which is not addressed in the current paper. In the last decades, the extensive use of fiberreinforced composites forced researchers to develop novel and practical techniques to realize nondestructive testing (NDT) on those materials. Such evolution is continuing today [32–35].
Understanding sound transmission through an anisotropic layer and how an acoustic beam deforms by it can be done in the same framework as studies aimed at nondestructive testing. It has been shown before by Rehman et al. [36] that the reflected beam pattern may reveal properties of composites. They showed that reflected Gaussian beams are only sensitive to material properties if oblique incidence is considered. The current work shows that the oblique incidence is not a requirement. Transmitted beams, for normal incidence, are sensitive to composites’ internal structure.
Even though the current work reveals a proof of concept and is numerical, it is essential to note what type of measurements could be performed by experimentalists for further investigations. One may, for instance, apply a needle transducer [37] as a receiver and perform a pointbypoint scan of the transmitted sound beam. The needle may be piezoelectric or a metalcoated fiber tip paired with a heterodyne interferometer that senses an optical path alteration caused by the fiber tip’s oscillation in the sound field [38]. Alternatively, a singlemode fiberoptic probe hydrophone (FOPH) can be used for local sound measurements [39], or Acoustooptic diffraction can also be applied to measure a sound field [40]. The latter typically comprises a global measurement of the entire beam and is insufficiently resolved for pointbypoint information. Tomographic variations of acoustooptic Bragg diffraction [41–47] can approach an adaptation of acoustooptics to pointbypoint measurements [48, 49] and Schlieren imaging [50, 51–56].
2 Incident beam descriptions
One may hardly imagine what sort of beams will be manufactured in future developments in metamaterials. As proof of concept, three types of 3D beams are assumed in this research: Gaussian beams, pillar beams, and conicallike beams.
Pillar beams look like a filled cylinders. The conicallike beams studied here have beam flanks showing exponential decay. There is a particular reason for this choice. The beauty of conicallike beams is that they possess properties similar to inhomogeneous waves. It is well known [57] that the propagation features of inhomogeneous waves in crystals depend not only on the direction of propagation but also on the value of the inhomogeneity vector. Therefore, whereas classical homogeneous plane waves cannot reflect composites’ internal symmetry for normal incidence, inhomogeneous waves are susceptible to symmetry because rotating the inhomogeneity vector (for normal incidence) must result in different propagation properties. A conicallike beam contains inhomogeneity properties for each radial direction and may reflect the fundamental inplane symmetry at once; therefore, it is worthwhile considering.
A Gaussian beam G(x, y, z), is characterized by its profile at z = 0 as
(1)whereas a conicallike beam C(x, y, z) by
(2)and a pillar beam P(x, y, z) by
In order to compare the different beams, we enforce the same amount of energy, whence
Because of the symmetry of the considered profiles, i.e., each profile K(x, y, 0) = K(R) with the radius, equation (4) are equivalent to
Therefore
Modeling ultrasound beams is done by Fourier analysis, where each beam is a sum of plane waves. Therefore, the transmitted beam is the sum of all transmitted plane waves originating from the plane waves constituting the perpendicular incident beam. The interaction of plane waves with an anisotropic layer is discussed in Section 3.
3 Sound transmission
Coupling a transducer, with or without a metamaterial in between, to an anisotropic disk could be done in different ways. The current paper merely focuses on the effect caused by the anisotropic material and assumes liquid coupling. It is also assumed that the transmitted beam arises on the opposite side of the disk in a liquid. The liquid in both cases is water, with a speed of sound of 1480 m/s and a density of 1000 kg/m^{3}. We also assume that the disk is large enough so the reflections from its edges do not occur, an assumption that makes sense in the case of perpendicular incident bounded beams. The material is considered orthotropic, typical for fiberreinforced composites.
The interaction of sound with an anisotropic layer is welldescribed in the literature [58–61]. A summary follows.
3.1 Material description
With Einstein’s double suffix notation convention, the dynamics of an anisotropic material is described [20, 62] by
With i, j = 1, 2, 3, and σ_{ij} being the stress tensor, and u the particle displacement. r is the position and t is time.
Hooke’s law, is given by
Orthotropy considerations [20, 62] result in
The stiffness constants of orthotropic materials are related to the Engineering constants, with Young’s moduli E_{11}, E_{22} and E_{33}, the Poisson coefficients ν_{23}, ν_{13} and ν_{12}, and the shear moduli G_{23}, G_{13} and G_{12}. The connection is
The reverse Voigt procedure transforms the stiffness tensor C_{mn} of rank 2 to the stiffness tensor c_{ijkl} of rank 4 as (1→11), (2→22), (3→33), (4→23 = 32), (5→13 = 31) and (6→12 = 21), after which we can conveniently rotate the composite as needed for the different plane waves constituting a bounded beam. A rotation R transforms the intrinsic stiffness constants into stiffness constants in coordinates corresponding to a rotated (laboratory) system as follows:
(15)where R_{ij} are the entries of the rotation matrix for a rotation from the intrinsic lattice coordinate system to the laboratory coordinate system.
Equation (7) then becomes
A plane wave solution of (16) is of the form
(17)where n is the wave vector. If this is entered in (16), straightforward calculations result in
Equation (18) is called the Christoffel equation [20, 62]. It relates the slowness n/ω and the polarization U to the propagation direction and is solved by assuming nontrivial solutions, followed by the determination of the corresponding eigenvectors.
3.2 Plane wave interactions with interfaces
Apart from the above determined Christoffel’s equation which is comparable to a dispersion relation for isotropic media, we must account for Snell’s law to determine the wave vector components along the interface, as well as the continuity condition.
3.2.a Snell’s law
If sound inside the bulk of the composite laminate results from impinging plane waves (denoted by superscript “inc”), Snell’s law for interfaces perpendicular to n_{3}, states that
Then, requiring nontrivial solutions, (18) leads to a sixthdegree polynomial equation of the form
Furthermore, symmetry higher than or equal to monoclinic symmetry results in
(21)whence three independent solutions for n_{3} are found.
3.2.b Continuity of normal stress and displacement
For a plate, continuity of normal stress and normal displacement is required along the watersolid interface,
4 Numerical results
Each of the considered beams, i.e. Gaussian, pillar and conicallike, is decomposed into many harmonic homogeneous plane waves by a 2D Fast Fourier Transform (FFT). The FFT is taken on the interval [−7W, 7W]^{2} on 2^{2N} equally distributed samples, for N the smallest positive integer for which
(24)and W is expressed in [m]. The number of samples equals the number of waves in the plane wave decomposition of the considered beam, and the considered interval is sufficiently extended to ensure reliable simulations.
This paper’s numerical simulations are based on calculating the transmission coefficient for every plane wave constituting the incident beam.
The considered composite has a thickness of 3 mm and consists of one layer of unidirectional material – carbon fibers in epoxy. The elastic properties of such material are given in Table 1. The fibers are directed in the XYplane along the Xdirection. Therefore, the material is stiffer along the Xdirection than the Ydirection [10].
Properties of a Carbon/Epoxy unidirectional fiberreinforced composites plate assumed as the material used for the 3 mm disk through which each sound beam passes.
Figure 1 shows the profile of an incident Gaussian beam, having the same amount of energy as a conicallike beam, shown in Figure 2, for β = 100/m. The comparative pillar beam is given in Figure 3. We only discuss transmitted beam profiles because that is this paper’s purpose.
Fig. 1 Profile of an incident Gaussian beam, having the same amount of energy as a conicallike beam, shown in Figure 2, for β = 100/m. 
Fig. 2 Profile of a conicallike beam for β = 100/m. This beam shows remarkable resemblance with an inhomogeneous wave. Its amplitude diminishes exponentially from the center of the beam. 
Fig. 3 Profile of an incident pillar beam, having the same amount of energy as a conicallike beam, shown in Figure 2, for β = 100/m. 
Figure 4 shows the transmitted beam profile for an incident Gaussian beam with the same amount of energy as a conicallike beam for β = 400/m. Note that the profile along the Xdirection differs from the one along the Ydirection. A similar effect, caused by the anisotropic disk, was found for the reflected beam, though less outspoken and less useful for the applications the topic of this paper aims at. The Gaussian beam is not very much influenced by the anisotropic disk, though comparing results for the other profiles are appealing.
Fig. 4 Transmitted beam profile for an incident Gaussian beam with the same energy as a conicallike beam for β = 400/m. Note that the profile along the Xdirection differs from the one along the Ydirection. The anisotropic disk does not significantly alter the Gaussian beam in any case. 
Figure 5 depicts the transmitted beam profile for a perpendicular incident conicallike beam β = 100/m. A significant difference exists between the profile along the Xdirection and the Ydirection. Numerous simulations show that the conicallike beam profile is always more sensitive than the Gaussian beam. The reason must be that conicallike beams have an exponentially decaying amplitude, similar to inhomogeneous waves [63]. Indeed, the scattering and propagation of inhomogeneous waves depend also on the inhomogeneity vector’s direction relative to the fiber direction [14].
Fig. 5 Transmitted beam profile for a normal incident conicallike beam β = 1000/m. There is a significant difference between the profile along the Xdirection and the Ydirection. 
If we increase the inhomogeneity β, as in Figure 6, where the transmitted beam profile for a perpendicularly incident conicallike beam, β = 2000/m, is given, we see that the pattern is much different from the one in Figure 5. This is because the scattering and propagation properties of inhomogeneous waves also depend on the magnitude of the inhomogeneity β. Another reason for the latter phenomenon is that narrower beams diffract much more, whence more propagation angles are involved, resulting in better sensitivity to the anisotropic elastic properties of the encountered material.
Fig. 6 Transmitted beam profile for a normal incident conicallike beam β = 2000/m. The pattern is much different from the one in Figure 5. 
Because conicallike beams might be challenging to generate experimentally, we have also performed simulations for pillar beams, which are presumably easier to generate. We found that pillar beams are much more sensitive to material properties and fiber direction than classical Gaussian beams. Figure 7 depicts the transmitted beam profile for an incident pillar beam with the same energy as a conicallike beam for β = 100/m. The transmitted beam pattern is directiondependent for this type of beam, even for a relatively large diameter.
Fig. 7 Transmitted beam profile for an incident pillar beam having the same amount of energy as a conicallike beam for β = 100/m. The transmitted beam pattern is direction dependent for this type of beam, even for a relatively large diameter. 
If we decrease the beam diameter, such as in Figure 8, where the transmitted beam profile for an incident pillar beam having the same amount of energy as a conicallike beam for β = 400/m, is shown, the pattern becomes more outspoken, certainly if compared to Figure 4.
Fig. 8 Transmitted beam profile for an incident pillar beam with the same energy as a conicallike beam for β = 400/m. The pattern is more outspoken than in Figure 7. The difference with Figure 4 is significant. 
In extreme situations, such as in Figure 9 (to be compared with Fig. 6) where the transmitted beam profile for an incident pillar beam having the same amount of energy as a conicallike beam for β = 2000/m, is shown, the transmitted beam profile has an even more complicated shape and depends very strongly on the direction in the XYplane.
Fig. 9 Transmitted beam profile for an incident pillar beam having the same amount of energy as a conicallike beam for β = 2000/m. The transmitted beam profile has a complicated shape and depends on the direction. The difference with Figure 6 is significant. 
Finally, to show the capabilities of the applied simulations, Figure 10 exposes an example of a reflected field for normal incidence, matching the case of Figure 7, where the transmission was shown. A study of how this reflected field would interact with the metamaterial and transducer is beyond the scope of this feasibility study for transmission. Further investigations would require a complete investigation, including an emitting transducer, a metamaterial, and the anisotropic disk as one complex system.
Fig. 10 Reflected beam profile for an incident pillar beam having the same amount of energy as a conicallike beam for β = 100/m. The result can be compared with Figure 7 of the transmitted profile. 
4 Conclusions
The report shows that transmitted beam profiles change when transmitted through an anisotropic disk. Furthermore, conicallike and pillar beams are more influenced than Gaussian beams. Even though this feasibility study used one particular unidirectional fiberreinforced composite layer as a model, depending on the requirements, the material can be tuned to match specific applications in which an anisotropic disk is placed in front of a single element transducer or a transducer covered by a metamaterial, to change transmitted beam characteristics further. Even though tunability may not be as sophisticated or flexible as a metasurface, its simplicity may add to the practical realization of specific beam profiles in a less expensive and lower level of fabrication sophistication than when only metamaterials are considered.
References
 J. Kim, S. Kasoji, P.G. Durham, P.A. Dayton: Acoustic hologram lens made of nanoparticleepoxy composite molding for directing predefined therapeutic ultrasound beams, in 2022 IEEE International Ultrasonics Symposium (IUS), Venice, Italy, 2022, pp. 1–4. https://doi.org/10.1109/IUS54386.2022.9957379. [Google Scholar]
 S. JimenezGambin, N. Jimenez, JM. Benlloch, F. Camarena: Holograms to focus arbitrary ultrasonic fields through the skull. Physical Review Applied 12, 1 (2019) 014016. [CrossRef] [Google Scholar]
 H. Ahmed, S. Ghosh, T. Sain, S. Banerjee: Hybrid Bessel beam and metamaterial lenses for deep laparoscopic nondestructive evaluation. Journal of Applied Physics 129, 16 (2021) 165107. [CrossRef] [Google Scholar]
 T. Yang, Y. Jin, T.Y. Choi, N. Dahotre, A. Neogi; Mechanically tunable ultrasonic metamaterial lens with a subwavelength resolution at long working distances for bioimaging. Smart Materials and Structures 30, 1 (2021) 015022. [CrossRef] [Google Scholar]
 Y. Jin, R. Kumar, O. Poncelet, O. MondainMonval, T. Brunet: Flat acoustics with soft gradientindex metasurfaces. Nature Communications 10 (2019) 143. [CrossRef] [PubMed] [Google Scholar]
 H. Gao, X. Fang, Z. Gu, T. Liu, S. Liang, Y. Li, J. Zhu: Conformally mapped multifunctional acoustic metamaterial lens for spectral sound guiding and Talbot effect. Research (Wash DC) 2019 (2019) 1748537. https://doi.org/10.34133/2019/1748537. [Google Scholar]
 Kun Li, Bin Liang, Jing Yang, Jun Yang, Jianchun Cheng: Broadband transmissiontype coding metamaterial for wavefront manipulation for airborne sound. Applied Physics Express 11, 7 (2018) 077301. [CrossRef] [Google Scholar]
 C. Shen, X. Jun, N.X. Fang, Y. Jing: Anisotropic complementary acoustic metamaterial for canceling out aberrating layers. Physical Review X 4, 4 (2014) 041033. [CrossRef] [Google Scholar]
 F. Semperlotti, H. Zhu: Achieving selective interrogation and subwavelength resolution in thin plates with embedded metamaterial acoustic lenses. Journal of Applied Physics 116, 5 (2014) 054906. [CrossRef] [Google Scholar]
 S. Yaacoubi, P. McKeon, W. Ke, N.F. Declercq, F. Dahmene: Towards an ultrasonic guided wave procedure for health monitoring of composite vessels: application to hydrogenpowered aircraft. Materials 10, 9 (2017) 1097. https://doi.org/10.3390/ma10091097. [CrossRef] [PubMed] [Google Scholar]
 M. Veidt, W. Sachse: Ultrasonic evaluation of thin, fiberreinforced laminates. Journal of Composite Materials 28, 4 (1994) 329–342. [CrossRef] [Google Scholar]
 W.P. Rogers: Elastic property measurement using RayleighLamb waves. Research in Nondestructive Evaluation 6, 4 (1995) 185–208. [CrossRef] [Google Scholar]
 P.B. Nagy, A. Jungman, L. Adler: Measurements of backscattered leaky Lamb waves in composite plates. Materials Evaluation 46, 1 (1988) 97–100. [Google Scholar]
 S. Eckel, F. Meraghni, P. Pomarede, N.F. Declercq: Investigation of damage in composites using nondestructive nonlinear acoustic spectroscopy. Experimental Mechanics 57 (2017) 20717. https://doi.org/10.1007/s1134001602226. [CrossRef] [Google Scholar]
 D.E. Chimenti, J. Song: Performance of spherically focused aircoupled ultrasonictransducers. AIP Conference Proceedings 894 (2007) 862. [CrossRef] [Google Scholar]
 S.D. Holland, S.V. Teles, D.E. Chimenti: Quantitative aircoupled ultrasonic materials characterization with highly focussed acoustic beams. Review of Progress in Quantitative Nondestructive Evaluation 23a and 23b (2004) 1376–1381. [CrossRef] [Google Scholar]
 D. Fei, D.E. Chimenti, S.V. Teles: Material property estimation in thin plates using focused, syntheticaperture acoustic beams. Journal of the Acoustical Society of America 113, 5 (2003) 2599–2610. [CrossRef] [PubMed] [Google Scholar]
 D.E. Chimenti, S.D. Holland, D. Fei: Aircoupled ultrasound and rapid elastic property characterization using focused acoustic beams. 2003 IEEE Ultrasonics Symposium Proceedings 1 and 2 (2003) 266–275. [CrossRef] [Google Scholar]
 S.D. Holland, S.V. Teles, D.E. Chimenti: Aircoupled, focused ultrasonic dispersion spectrum reconstruction in plates. Journal of the Acoustical Society of America 115, 6 (2004) 2866–2872. [CrossRef] [Google Scholar]
 A.H. Nayfeh, D. E. Chimenti: Propagation of guided waves in fluidcoupled plates of fiberreinforced composite. Journal of the Acoustical Society of America 83, 5 (1988) 1736–1743. [CrossRef] [Google Scholar]
 W. Sachse, Y.H. Pao: Determination of phase and group velocities of dispersive waves in solids. Journal of Applied Physics 49, 8 (1978) 4320–4327. [CrossRef] [Google Scholar]
 M. Deschamps, B. Hosten: The effects of viscoelasticity on the reflection and transmission of ultrasonicwaves by an orthotropic plate. Journal of the Acoustical Society of America 91, 4 (1992) 2007–2015. [CrossRef] [PubMed] [Google Scholar]
 R.L. Weaver, W. Sachse, L. Niu: Transient ultrasonicwaves in a viscoelastic plate – applications to materials characterization. Journal of the Acoustical Society of America 85, 6 (1989) 2262–2267. [CrossRef] [Google Scholar]
 L. Satyanarayan, JM. Vander Weide, N.F. Declercq: Ultrasonic polar scan imaging of damaged fiberreinforced composites, Materials Evaluation 68, 6 (2010) 733–739. [Google Scholar]
 R. Raišutis, O. Tumšys: Application of dual focused ultrasonic phased array transducer in two orthogonal crosssections for inspection of multilayered composite components of the aircraft fuselage. Materials (Basel) 13, 7 (2020) 1689. https://doi.org/10.3390/ma13071689. [CrossRef] [PubMed] [Google Scholar]
 D. Hopkins, M. Datuin, M. Brassard: Challenges and solutions for ultrasonic phasedarray inspection of polymermatrix composites at production rates, in 45th Annual Review of Progress in Quantitative Nondestructive Evaluation vol. 38, AIP Conference Proceedings 2102, UNSP 1000021, 2019. [Google Scholar]
 D.W. Schindel: Ultrasonic imaging of solid surfaces using a focussed aircoupled capacitance transducer. Ultrasonics 35, 8 (1998) 587–594. [CrossRef] [Google Scholar]
 D.E. Chimenti, D. Fei: Scattering coefficient reconstruction in plates using focused acoustic beams. International Journal of Solids and Structures 39, 21–22 (2002) 5495–5513. [CrossRef] [Google Scholar]
 V.M. Levin, O.I. Lobkis, R.G. Maev: Investigation of the spatial structure of acoustic fields by a spherical focusing transducer. Soviet Physics AcousticsUSSR 36, 4 (1990) 391–395. [Google Scholar]
 B. Hosten, D.A. Hutchins, D.W. Schindel: Measurement of elastic constants in composite materials using aircoupled ultrasonic bulk waves. Journal of the Acoustical Society of America 99, 4 (1996) 2116–2123. [CrossRef] [Google Scholar]
 A. Safaeinili, O.I. Lobkis, D.E. Chimenti: Aircoupled ultrasonic estimation of viscoelastic stiffnesses in plates. IEEE Transactions on Ultrasonics Ferroelectrics and Frequency Control 43, 6 (1996) 1171–1180. [CrossRef] [Google Scholar]
 N. Miqoi, P. Pomarede, N.F. Declercq, L. Guillaumat, G. Le Coz, S. Delalande, F. Meraghni: Detection and evaluation of barely visible impact damage in woven glass fabric reinforced polyamide 6.6/6 using ultrasonic imaging, Xray tomography and optical profilometry. International Journal of Damage Mechanics 30 (2020) 323–348. https://doi.org/10.1177/1056789520957703. [Google Scholar]
 P. Pomarède, L. Chehami, N.F. Declercq, F. Meraghni, J. Dong, A. Locquet, D.S. Citrin: Application of ultrasonic coda wave interferometry for microcracks monitoring in woven fabric composites. Journal of Nondestructive Evaluation, Springer Verlag 38, 1 (2019) 26–34. [CrossRef] [Google Scholar]
 J. Dong, P. Pomarede, L. Chehami, A. Locquet, F. Meraghni, N.F. Declercq, D.S. Citrin: Visualization of subsurface damage in woven carbon fiberreinforced composites using polarizationsensitive terahertz imaging. NDT and E International 99 (2018) 72–79. [CrossRef] [Google Scholar]
 P. Pomarède, F. Meraghni, L. Peltier, S. Delalande, N.F. Declercq: Damage evaluation in woven glass reinforced polyamide 6.6/6 composites using ultrasound phaseshift analysis and XRay tomography. Journal of Nondestructive Evaluation 73, 12 (2018) 1–21. [Google Scholar]
 A.U. Rehman, C. Potel, J.F. de Belleval: Numerical modeling of the effects on reflected acoustic field for the changes in internal layer orientation of a composite. Ultrasonics 36 (1998) 343–348. [CrossRef] [Google Scholar]
 C. Koch: Sound field measurement in a double layer cavitation cluster by rugged miniature needle hydrophones. Ultrasonics Sonochemistry 29 (2016) 439–446. [CrossRef] [PubMed] [Google Scholar]
 C. Koch, K.V. Jenderka: Measurement of sound field in cavitating media by an optical fibretip hydrophone. Ultrasonics Sonochemistry 15, 4 (2008) 502–509. [CrossRef] [PubMed] [Google Scholar]
 J. Petelin, Z. Lokar, D. Horvat, R. Petkovsek: Localized measurement of a subnanosecond shockwave pressure rise time. IEEE Transactions on Ultrasonics Ferroelectrics and Frequency Control 69, 1 (2022) 369–376. [CrossRef] [PubMed] [Google Scholar]
 L. Jia, S. Chen, B. Xue, H. Wu, K. Zhang, X. Yang, Z. Zeng: Acoustic pressure measurement of pulsed ultrasound using acoustooptic diffraction. Proceedings of SPIE 10621 (2017) 75–84. [Google Scholar]
 M.G. Moharam, L. Young: Criterion for Bragg and RamanNath diffraction regimes, Applied Optics 17, 11 (1978) 1757–1759. [CrossRef] [PubMed] [Google Scholar]
 A. Korpel: Visualization of cross section of a sound beam by Bragg diffraction of light. Applied Physics Letters 9, 12 (1966) 425–427. [CrossRef] [Google Scholar]
 A. Korpel: Proceedings of the Second International Symposium on Acoustical Holography. Plenum, London, England, 1970, p. 39. [CrossRef] [Google Scholar]
 A. Korpel: AcoustoOptics, 2nd ed. (Marcel Dekker Inc, New York, 1997, p. 21–22, 206–219. [Google Scholar]
 K. Vandenabeele, M.A. Breazeale, O. Leroy, J.K. Na: Strong Interaction of arbitrary fields of sound and light – application to higher order Bragg Imaging. Journal of Applied Physics 75, 1 (1994) 84–95. [CrossRef] [Google Scholar]
 J.K. Na, M.A. Breazeale, O. Leroy: Ultrasonic Bragg imaging of flaws, Journal of the Acoustical Society of America 81, Suppl. 1 (1987) S43–S43. [Google Scholar]
 L.H.V. Wang: Ultrasoundmediated biophotonic imaging: A review of acoustooptical tomography and photoacoustic tomography. Disease Markers 19, 2–3 (2003) 123–138. [Google Scholar]
 A. Teklu, N.F. Declercq, M. McPherson: Acoustooptic Bragg imaging of biological tissue. Journal of the Acoustical Society of America 136, 2 (2014) 634–637. [CrossRef] [PubMed] [Google Scholar]
 N.F. Declercq, M.S. McPherson, M.A. Breazeale, A.A. Teklu: Optical Bragg imaging of acoustic fields after reflection. Journal of the Acoustical Society of America 127, 6 (2010) 3466–3469. [CrossRef] [PubMed] [Google Scholar]
 NF Declercq, A. Teklu, M.A. Breazeale, R.D. Hasse, J. Degrieck, O. Leroy: Detection of fiber direction in composites by means of high frequency wide bounded ultrasonic beam and Schlieren photography. Research in Nondestructive Evaluation 16, 2 (2005) 55–64. [CrossRef] [Google Scholar]
 G. Cammi, A. Spinelli, F. Cozzi, A. Guardone: Automatic detection of oblique shocks and simple waves in Schlieren images of twodimensional supersonic steady flows. Measurement 168 (2021) 108260. [CrossRef] [Google Scholar]
 H.D. Lim, X.F. Wei, B. Zang, U.S. Vevek, R. Mariani, T.H. New, Y.D. Cui: Shorttime proper orthogonal decomposition of timeresolved Schlieren images for transient jet screech characterization. Aerospace Science and Technology 107 (2020) 106276. [CrossRef] [Google Scholar]
 E. Lampsijärvi, J. Heikkilä, I. Kassamakov, A. Salmi, E. Hæggström: Calibrated quantitative stroboscopic Schlieren imaging of ultrasound in air, in IEEE International Ultrasonics Symposium (IUS), Glasgow, UK 2019 (2019) 1651–1654. https://doi.org/10.1109/ULTSYM.2019.8925916. [Google Scholar]
 Z. Xu, H. Chen, X. Yan, M.L. Qian, Q. Cheng: Threedimensional reconstruction of nonplanar ultrasound fields using Radon transform and the Schlieren imaging method. Journal of the Acoustical Society of America 142 (2017). EL82–EL88. https://doi.org/10.1121/1.4994282. [CrossRef] [PubMed] [Google Scholar]
 G. Caliano, A.S. Savoia, A. Iula: An automatic compact Schlieren imaging system for ultrasound transducer testing. IEEE Transactions on Ultrasonics Ferroelectric and Frequency Control 59, 9 (2012) 2102–2110. https://doi.org/10.1109/TUFFC.2012.2431. [Google Scholar]
 M. Ohno, N. Tanaka, Y. Matsuzaki: Schlieren imaging by the interference of two beams in RamanNath diffraction. Japanese Journal of Applied Physics 42, 5b (2003) 3067–3071. [CrossRef] [Google Scholar]
 N.F. Declercq, R. Briers, J. Degrieck, O. Leroy: The history and properties of ultrasonic inhomogeneous waves. IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control 52, 5 (2005) 776–791. https://doi.org/10.1109/TUFFC.2005.1503963. [CrossRef] [PubMed] [Google Scholar]
 A.H. Naefeh: Wave propagation in layered anisotropic media with applications to composites, in North Holland series in Applied Mathematics and Mechanics, 1995. [Google Scholar]
 S.I. Rokhlin, W. Wang: Double throughtransmission bulk wave method for ultrasonic phasevelocity measurement and determination of elasticconstants of compositematerials. Journal of the Acoustical Society of America 91, 6 (1992) 3303–3312. [CrossRef] [Google Scholar]
 B. Hosten, M. Deschamps, B.R. Tittmann: Inhomogeneous wave generation and propagation in lossy anisotropic solids – application to the characterization of viscoelastic compositematerials. Journal of the Acoustical Society of America 82, 5 (1987) 1763–1770. [CrossRef] [Google Scholar]
 O.I. Lobkis, D.E. Chimenti, H. Zhang: Inplane elastic property characterization in composite plates. Journal of the Acoustical Society of America 107, 4 (2000) 1852–1858. [CrossRef] [PubMed] [Google Scholar]
 M. Deschamps, B. Hosten: The effects of viscoelasticiy on the reflection and transmission of ultrasonic waves by an orthotropic plate. Journal of the Acoustical Society of America 91, 4 (1992) 2007–2015. [CrossRef] [PubMed] [Google Scholar]
 N.F. Declercq: Experimental study of ultrasonic beam sectors for energy conversion into Lamb waves and Rayleigh waves. Ultrasonics 54, 2 (2013) 609–613. [Google Scholar]
Cite this article as: Declercq NF. 2023. Numerical study of beam shape adaptation by anisotropic disk covering transducer or metamaterial. Acta Acustica, 7, 30.
All Tables
Properties of a Carbon/Epoxy unidirectional fiberreinforced composites plate assumed as the material used for the 3 mm disk through which each sound beam passes.
All Figures
Fig. 1 Profile of an incident Gaussian beam, having the same amount of energy as a conicallike beam, shown in Figure 2, for β = 100/m. 

In the text 
Fig. 2 Profile of a conicallike beam for β = 100/m. This beam shows remarkable resemblance with an inhomogeneous wave. Its amplitude diminishes exponentially from the center of the beam. 

In the text 
Fig. 3 Profile of an incident pillar beam, having the same amount of energy as a conicallike beam, shown in Figure 2, for β = 100/m. 

In the text 
Fig. 4 Transmitted beam profile for an incident Gaussian beam with the same energy as a conicallike beam for β = 400/m. Note that the profile along the Xdirection differs from the one along the Ydirection. The anisotropic disk does not significantly alter the Gaussian beam in any case. 

In the text 
Fig. 5 Transmitted beam profile for a normal incident conicallike beam β = 1000/m. There is a significant difference between the profile along the Xdirection and the Ydirection. 

In the text 
Fig. 6 Transmitted beam profile for a normal incident conicallike beam β = 2000/m. The pattern is much different from the one in Figure 5. 

In the text 
Fig. 7 Transmitted beam profile for an incident pillar beam having the same amount of energy as a conicallike beam for β = 100/m. The transmitted beam pattern is direction dependent for this type of beam, even for a relatively large diameter. 

In the text 
Fig. 8 Transmitted beam profile for an incident pillar beam with the same energy as a conicallike beam for β = 400/m. The pattern is more outspoken than in Figure 7. The difference with Figure 4 is significant. 

In the text 
Fig. 9 Transmitted beam profile for an incident pillar beam having the same amount of energy as a conicallike beam for β = 2000/m. The transmitted beam profile has a complicated shape and depends on the direction. The difference with Figure 6 is significant. 

In the text 
Fig. 10 Reflected beam profile for an incident pillar beam having the same amount of energy as a conicallike beam for β = 100/m. The result can be compared with Figure 7 of the transmitted profile. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.