OE Letters

Fractional Fourier analysis of elastic wave scattering in inhomogeneous materials

[+] Author Affiliations
Yuanzhang Fan

Beijing Institute of Technology, School of Information and Electronics, 5 South Zhongguancun Street, Haidian District, Beijing 100081, China

Zhongyuan University of Technology, School of Electronic and Information Engineering, 41 Zhongyuan Road (M), Zhengzhou 450007, China

Xiangyang Lu

Zhongyuan University of Technology, School of Electronic and Information Engineering, 41 Zhongyuan Road (M), Zhengzhou 450007, China

Opt. Eng. 55(11), 110501 (Nov 07, 2016). doi:10.1117/1.OE.55.11.110501
History: Received August 30, 2016; Accepted October 7, 2016
Text Size: A A A

Open Access Open Access

Abstract.  To describe the elastic wave scattering, which reflects the performance of propagation control materials, the approximate directional cloaks of an elastic wave are designed using the zero- and first-order approximation coordinate transformation method. Because the convergence features of fractional Fourier transform (FRFT) are more acute and sensitive to the frequency change than those of short time Fourier transform, the spatial signals in the designed materials are transformed in the FRFT domain. The spatial frequency changes of elastic waves through inhomogeneous materials are quantitatively analyzed under several circumstances. The provided time-frequency analysis method with FRFT can support the design evaluation of the material parameters.

Figures in this Article

The coordinate transformation method, proposed for electromagnetic waves controlling and manipulating metamaterials by Leonhardt and Pendry et al., has been extended to elastic waves in recent years.14 Due to zero- and first-order approximation method itself and the impedance mismatch in the materials, the scattering phenomena, which refers to the frequency changes as waves pass through the media, is caused in the material design of wave propagation control.5 Generally, the operating mechanism of elastic waves is more complicated and difficult to control than electromagnetic waves.3 Furthermore, according to the existing control theory of wave propagation, some assumptions about the characteristics of the material, such as the anisotropy of the medium parameters, the continuity in the radial direction, and lossless dielectric are necessary but difficult to implement in practice. It is common that scattering exists in material designs of wave propagation control, including the directional cloak design of elastic waves. To reduce the scattering effects, many methods that simplify the sets of medium parameters or lower the complexity of the parameters are utilized.68 However, as a precondition, these scattering phenomena should be described adequately and grasped accurately, which, especially on the material boundary, only relies on macroscopic observation of the time domain waveform for qualitative conclusions without quantitative analysis.9

In analyzing the dynamic behavior and characteristics of the elastic media, the complexity and computational burden can be reduced by the models of homogeneous materials,10,11 which are inhomogeneous in practice.12 The deformation mechanisms of waves are complex in media, often elastic, plastic, nonlinear, and anisotropic. Meanwhile, the resulting models, based on material composition, structure, temperature, stress paths, history, and so on, are generally applied under strict conditions.

Since scattering phenomena lead to spatial frequency changes, the fractional Fourier transform (FRFT)13,14 is a mathematically suitable tool to analyze scattering phenomena in both time and frequency domains and even in more complex scenarios of the propagation medium. Based on the FRFT analysis of scattering phenomena, we established the correspondence between the maximum matching order of transformation and the scattering, and provided the basis for assessing the material design of elastic wave propagation control.

It is complex and difficult to give a completely accurate description of the elastic waves’ control under normal conditions. But in some applications, the first-order approximate equations are sufficient. The governing equations of elastodynamics are Display Formula

σ=ω2ρu,σ=Cu,(1)
where σ is the second-order stress tensor in a certain direction, ω is the frequency, ρ is the density, u is the displacement vector, and C is the fourth-order elasticity tensor. Removing the tensor item and then obtaining Navier’s equation about the displacement of elastic wave control Display Formula
i(Cijklkul)=ω2ρijui.(2)

Its first-order approximation is written as (without summation) Display Formula

σij+12σijxidxi=σji+12σjixjdxj,i,j=1,2,3,ij,(3)
where dxi is the size of an element. As Eq. (3) is difficult to solve, it can be simplified to Display Formula
σij=σji.(4)
by zero-order approximation with symmetric elasticity tensor.

In fact, Eq. (4) is proper only for high-frequency waves or slowly varying materials. That is, this approximation method can be used to design corresponding materials, which also results in scattering phenomena. The coordinate transformation method to design elastic wave propagation control devices is shown in Refs. 15 and 16. We provide the design method of linear change materials as follows.

To ensure the materials change linearly in the design, we use an undetermined coefficients method to solve the transformation relationship, which is expressed as Display Formula

α=ax+b.(5)

When the materials change from 1 to N multiples of transformations, e.g., N=2, the equation Display Formula

{2.5x+b=12.5x+b=2(6)
can be obtained.

The factor of change materials can then be obtained by solving this equation. The material parameters of elastic waves vary in accordance with the above equations, and then the corresponding linear change materials can be obtained by the transformation method.

In linear change materials, the spatial frequency of an elastic wave also varies linearly as it passes through the medium. The simulation results about the signal propagation with linear change media are shown in Fig. 1, where the constant scaling factor λ=3 and the background media is structural steel with material parameters E=200  Gpa, υ=0.33, and ρ=7850  kg/m3. The obstacle has material parameters Eobs=10E, ρobs=ρ/10, and υobs=υ.

Graphic Jump Location
Fig. 1
F1 :

Displacement vector simulation in inhomogeneous materials. (a) Original waves without obstacle, (b) waves without cloak, and (c) waves with the designed cloak.

Fractional Fourier Transform

The scattering effect leads to the change of signal spatial frequency in dielectric materials, i.e., instantaneous frequency changes as the space varies. To analyze the nonstationary signal spatial frequency change, we use the FRFT method, a signal analysis tool with better convergence that is more focused on energy than the Fourier transform.13 The FRFT of a signal is defined by Display Formula

Xp(u)={Fp[x(t)]}(u)=x(t)·Kp(t,u)dt,(7)
where α=pπ/2 is the rotation angle, p is the fractional order of transformation, and Display Formula
K(t,u)={1j·cotα2π·exp(j·t2+u22·cotαj·u·tcscα),αnπδ(tu),α=2nπδ(t+u),α=(2n+1)π(8)
is the transformation kernel.

This means that the signal x(t) is expanded into a set of orthogonal basis of the transform kernel Kp(t,u) function space. Then the u domain is the so-called fractional Fourier domain. The FRFT degenerates to a classical Fourier transform when p=1.

Estimation Performance Comparison between Fractional Fourier Transform and Short Time Fourier Transform

An example illustrates that the FRFT has better convergence than the STFT with the following general quintic polynomial (the signal to match), as shown in Fig. 2Display Formula

x(t)=exp[j2π*(100*t5+100*t4+160*t3+10*t2+30*t)].(9)

Graphic Jump Location
Fig. 2
F2 :

Fitting results of FRFT and STFT on a general quintic polynomial. (a) The original signal, (b) STFT estimation, and (c) FRFT estimation.

Figure 2(a) is the original polynomials to fit. Figures 2(b) and 2(c) are the fitting results with STFT and FRFT, respectively. As can be seen from the figures, FRFT has better focused features and a better estimation performance than STFT.

Fractional Fourier Analysis and Results

To transform the above spatial domain waveforms into the fractional Fourier domain, substitute the signals acquired from the materials into Eq. (7), and then get the spatial frequency diagram, as shown in Fig. 3. The fractional order change in the transformation is a search process to find the greatest undulation in frequency. We use a method, the so-called “first-rough and after-accurate scale,” to improve the frequency accuracy: first, scan the whole rotation angle range to find the rough maximum in the transform domain with a big fractional order step, then search the accurate maximum near the rough maximum with a small fractional order step.

Graphic Jump Location
Fig. 3
F3 :

Analysis of scattering effects on directional cloaks. Horizontal: (a) without obstacles, (b) with obstacles but without cloaks, and (c) with both obstacles and cloaks. Vertical: (d) waveforms of intermediate position in the designed inhomogeneous materials; (e) 3-D diagrams of FRFT signals in (d) column, where x-label is the transformation angle, y-label is the original signal, and z-label is the FRFT result; and (f) matching order of FRFT results.

Figure 3(d) shows the signal waveform of the elastic wave propagation in the middle of the linear change materials. In the scenario of linear design changes of 1 to 2 times and without obstacles, the spatial frequency of the elastic wave in the middle of materials changed linearly.

It is found from the waveforms in the middle of linear change elastic materials that the waveform is changed linearly. The frequency modulation rate is changed from 37.2  Hz/s [the basic waveform without obstacles as shown in Figs. 3(d)3(a)] to 3.1  Hz/s [with the obstacle in the materials as shown in Figs. 3(d)3(b)], and is improved to 58.2  Hz/s with the designed obstacle and cloak [Figs. 3(d)3(c)]. The signal modulation frequency rate is greatly changed when there is no cloak. And when the cloak design is added, the signal is improved and is closer to the original signal.

The three-dimensional (3-D) diagram of above signals FRFT is shown in Fig. 3(e). When there are no obstacles in the medium, the waveform is a linear frequency modulated signal. When there is an obstacle, the frequency modulation rate is changed. But after adding the cloak, it is better compared with the design without cloaks.

Three cases of FRFT results with matching order are shown in Fig. 3(f). This figure shows that it is smooth in the spatial FRFT domain in the absence of obstacles [Figs. 3(f)3(a)], and new frequency components are increased in the case of with obstacles [Figs. 3(f)3(b)]. Nevertheless, the increased parts become focused with both obstacle and cloak scenarios [Figs. 3(f)3(c)].

From the above results, whether in spatial domain diagrams or in FRFT transform domain diagrams, it can be seen that the signals are significantly improved after the cloak is added.

Based on the method of the continuum theory, an elastic wave approximation cloak can be designed to improve the signals interrupted by obstacles. But scattering phenomena occur in the propagation process. FRFT is a suitable tool used in the signal processing with scattering variable parameters. FRFT clearly plays a very important role in the variable parameters signal processing, as it can not only analyze the change of signals and evaluate the design of material parameters, but also quantitatively illustrate the results of the transformation, thus providing support for the evaluation of the design methods. This paper studies the design of a simple device and the theoretical analysis using the zero- and first-order approximation method. In practical applications, it is required to further consider the impact of material factors, such as the size and shape effects on the scattering.

This work was supported by the National Natural Science Foundation of China (61301089, 61471045, 61421001 and 61575022), Natural Science Research Fund of Education Department of Henan Province (15A140044), and Basic & Advanced Technology Research Fund of Henan Province (162300410200).

Leonhardt  U., “Optical conformal mapping and dielectric invisibility devices,” Science. 56, (12 ), 1 –4 (2006). 0036-8075 
Pendry  J. B., , Schurig  D., and Smith  D. R., “Controlling electromagnetic fields,” Science. 312, (5781 ), 1780 –1782 (2006). 0036-8075 CrossRef
Chang  Z., , Hu  J., and Hu  G. K., “Transformation method and wave control,” Acta Mech. Sin.. 26, (6 ), 889 –898 (2010). 0567-7718 CrossRef
Hai  X. Z., “The form-invariance of wave equations without requiring a priori relations between field variables,” Sci. China Phys. Mech. Astron.. 57, (12 ), 2285 –2296 (2014).CrossRef
Hu  J., , Liu  X., and Hu  G., “Constraint condition on transformed relation for generalized acoustics,” Wave Motion. 50, (2 ), 170 –179 (2009). 0165-2125 CrossRef
AlMuhaidib  A. M.  et al., “Finite difference elastic wave modeling including surface topography,” in SEG Technical Program Expanded. , pp. 2941 –2946 (2011).
Brun  M., , Guenneau  S., and Movchan  A. B., “Achieving control of in-plane elastic waves,” Appl. Phys. Lett.. 94, (6 ), 061903  (2009). 0003-6951 CrossRef
Mei  J.  et al., “Theory for elastic wave scattering by a two-dimensional periodical array of cylinders: an ideal approach for band-structure calculations,” Phys. Rev. B. 67, (24 ), 841 –845 (2003). 1098-0121 CrossRef
Lee  J., “Analysis of elastic wave scattering problems using volume integral equation method,” J. Mech. Sci. Technol.. 22, (22 ), 228 –239 (2008).CrossRef
Labate  G., , Matekovits  L., and Orefice  M., “Inverse scattering homogenization method for conformal metamaterial structures,” in  2015 IEEE Int. Symp. on Antennas and Propagation and USNC/URSI National Radio Science Meeting , pp. 63 –64 (2015).CrossRef
Lu  X., , Hu  J., and Tao  R., “Enhanced fractional Fourier lens with isotropic transformation media,” Opt. Eng.. 52, (6 ), 060501  (2013).CrossRef
Tsang  M., and Psaltis  D., “Magnifying perfect lens and superlens design by coordinate transformation,” Phys. Rev. B. 77, (3 ), 035122  (2008). 1098-0121 CrossRef
Tao  R., , Zhang  F., and Wang  Y., “Fractional power spectrum,” IEEE Trans. Signal Process.. 56, (9 ), 4199 –4206 (2008).CrossRef
Sohn  A., “Fractional Fourier transforms and their optical implementation. II,” J. Opt. Soc. Am. A Opt. Image Sci.. 37, (11 ), 2130 –2141 (1993). CrossRef
Hu  J., , Chang  Z., and Hu  G., “An approximate method for controlling solid elastic waves by transformation media,” Physics. 84, (20 ), 201101  (2010). 1943-2879 CrossRef
Zheng  C.  et al., “Controlling elastic waves with isotropic materials,” Appl. Phys. Lett.. 98, (12 ), 121904  (2011). 0003-6951 CrossRef
© The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.

Citation

Yuanzhang Fan and Xiangyang Lu
"Fractional Fourier analysis of elastic wave scattering in inhomogeneous materials", Opt. Eng. 55(11), 110501 (Nov 07, 2016). ; http://dx.doi.org/10.1117/1.OE.55.11.110501


Figures

Graphic Jump Location
Fig. 1
F1 :

Displacement vector simulation in inhomogeneous materials. (a) Original waves without obstacle, (b) waves without cloak, and (c) waves with the designed cloak.

Graphic Jump Location
Fig. 2
F2 :

Fitting results of FRFT and STFT on a general quintic polynomial. (a) The original signal, (b) STFT estimation, and (c) FRFT estimation.

Graphic Jump Location
Fig. 3
F3 :

Analysis of scattering effects on directional cloaks. Horizontal: (a) without obstacles, (b) with obstacles but without cloaks, and (c) with both obstacles and cloaks. Vertical: (d) waveforms of intermediate position in the designed inhomogeneous materials; (e) 3-D diagrams of FRFT signals in (d) column, where x-label is the transformation angle, y-label is the original signal, and z-label is the FRFT result; and (f) matching order of FRFT results.

Tables

References

Leonhardt  U., “Optical conformal mapping and dielectric invisibility devices,” Science. 56, (12 ), 1 –4 (2006). 0036-8075 
Pendry  J. B., , Schurig  D., and Smith  D. R., “Controlling electromagnetic fields,” Science. 312, (5781 ), 1780 –1782 (2006). 0036-8075 CrossRef
Chang  Z., , Hu  J., and Hu  G. K., “Transformation method and wave control,” Acta Mech. Sin.. 26, (6 ), 889 –898 (2010). 0567-7718 CrossRef
Hai  X. Z., “The form-invariance of wave equations without requiring a priori relations between field variables,” Sci. China Phys. Mech. Astron.. 57, (12 ), 2285 –2296 (2014).CrossRef
Hu  J., , Liu  X., and Hu  G., “Constraint condition on transformed relation for generalized acoustics,” Wave Motion. 50, (2 ), 170 –179 (2009). 0165-2125 CrossRef
AlMuhaidib  A. M.  et al., “Finite difference elastic wave modeling including surface topography,” in SEG Technical Program Expanded. , pp. 2941 –2946 (2011).
Brun  M., , Guenneau  S., and Movchan  A. B., “Achieving control of in-plane elastic waves,” Appl. Phys. Lett.. 94, (6 ), 061903  (2009). 0003-6951 CrossRef
Mei  J.  et al., “Theory for elastic wave scattering by a two-dimensional periodical array of cylinders: an ideal approach for band-structure calculations,” Phys. Rev. B. 67, (24 ), 841 –845 (2003). 1098-0121 CrossRef
Lee  J., “Analysis of elastic wave scattering problems using volume integral equation method,” J. Mech. Sci. Technol.. 22, (22 ), 228 –239 (2008).CrossRef
Labate  G., , Matekovits  L., and Orefice  M., “Inverse scattering homogenization method for conformal metamaterial structures,” in  2015 IEEE Int. Symp. on Antennas and Propagation and USNC/URSI National Radio Science Meeting , pp. 63 –64 (2015).CrossRef
Lu  X., , Hu  J., and Tao  R., “Enhanced fractional Fourier lens with isotropic transformation media,” Opt. Eng.. 52, (6 ), 060501  (2013).CrossRef
Tsang  M., and Psaltis  D., “Magnifying perfect lens and superlens design by coordinate transformation,” Phys. Rev. B. 77, (3 ), 035122  (2008). 1098-0121 CrossRef
Tao  R., , Zhang  F., and Wang  Y., “Fractional power spectrum,” IEEE Trans. Signal Process.. 56, (9 ), 4199 –4206 (2008).CrossRef
Sohn  A., “Fractional Fourier transforms and their optical implementation. II,” J. Opt. Soc. Am. A Opt. Image Sci.. 37, (11 ), 2130 –2141 (1993). CrossRef
Hu  J., , Chang  Z., and Hu  G., “An approximate method for controlling solid elastic waves by transformation media,” Physics. 84, (20 ), 201101  (2010). 1943-2879 CrossRef
Zheng  C.  et al., “Controlling elastic waves with isotropic materials,” Appl. Phys. Lett.. 98, (12 ), 121904  (2011). 0003-6951 CrossRef

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging & repositioning the boxes below.

Related Book Chapters

Topic Collections

Advertisement
  • Don't have an account?
  • Subscribe to the SPIE Digital Library
  • Create a FREE account to sign up for Digital Library content alerts and gain access to institutional subscriptions remotely.
Access This Article
Sign in or Create a personal account to Buy this article ($20 for members, $25 for non-members).
Access This Proceeding
Sign in or Create a personal account to Buy this article ($15 for members, $18 for non-members).
Access This Chapter

Access to SPIE eBooks is limited to subscribing institutions and is not available as part of a personal subscription. Print or electronic versions of individual SPIE books may be purchased via SPIE.org.