Open AccessArticle

Theoretical Zero-Thickness Broadband Holograms Based on Acoustic Sieve Metasurfaces

Ye Tian

Shuyu Zuo

^2,*,

Qian Lv

Guanjun Yin

and

Jianzhong Guo

^1,*

Key Laboratory of Ultrasound of Shaanxi Province, School of Physics and Information Technology, Shaanxi Normal University, Xi’an 710062, China

College of Science, Henan University of Technology, Zhengzhou 450001, China

Authors to whom correspondence should be addressed.

Appl. Sci. 2022, 12(13), 6453; https://doi.org/10.3390/app12136453

Submission received: 29 April 2022 / Revised: 14 June 2022 / Accepted: 22 June 2022 / Published: 25 June 2022

(This article belongs to the Special Issue Recent Advance in Acoustic Metamaterials)

Download

Browse Figures

Review Reports Versions Notes

Abstract

Acoustic holography is an essential tool for controlling sound waves, generating highly complex and customizable sound fields, and enabling the visualization of sound fields. Based on acoustic sieve metasurfaces (ASMs), this paper proposes a theoretical design approach for zero-thickness broadband holograms. The ASM is a zero-thickness rigid screen with a large number of small holes that allow sound waves to pass through and produce the desired real image in the target plane. The hole arrangement rules are determined using a genetic algorithm and the Rayleigh–Sommerfeld theory. Because the wave from a hole has no extra phase or amplitude modulation, the intractable modulation dispersion can be physically avoided, allowing the proposed ASM-based hologram to potentially function in any frequency band as long as the condition of paraxial approximation is satisfied. Using a numerical simulation based on the combination of the finite element method (FEM) and the boundary element method (BEM), this research achieves broadband holographic imaging with a good effect. The proposed theoretical zero-thickness broadband hologram may provide new possibilities for acoustic holography applications.

Keywords:

acoustic holography; metasurface; broadband hologram; ultrathin hologram; sound field control

1. Introduction

Acoustic holography [1,2] is an essential tool for controlling sound waves, generating highly complex and customizable sound fields, and enabling the visualization of sound fields. Due to its powerful and flexible capability of arbitrary sound beam shaping and pre-designed sound field reconstruction, acoustic holography has recently received increased attention in many interdisciplinary fields related to acoustics, including medicine [3,4,5,6], engineering [7,8], and biology [9,10,11], in addition to the pure acoustic field. Many applications have been made possible by acoustic holography, including acoustic fabrication [7], beam shaping [12,13,14,15,16,17,18,19,20], acoustic holographic imaging [3,15,16,17,18,19,20,21,22,23,24,25,26,27,28], volumetric display [29,30], volumetric haptics [31,32], particle manipulation [11,21,26,33,34], 3D ultrasound imaging [8], cell manipulation [9,10], characterizing medical ultrasounds [4,5], neurostimulation and neuromodulation [6], and energy transfer [25].

Instead of acoustically recording a hologram from a physical scene with an exacting recording environment and a complex process, modern computer-generated holography directly calculates the required phase and/or amplitude maps in the acoustic holographic plane before rendering them for reconstruction [35]. The algorithms used include the iterative angular spectrum approach [21], the iterative backpropagation algorithm [33], the weighted Gerchberg–Saxton algorithm [20], the pseudo-inverse method [36], the conjugate field method [37], the deep learning method [35], and so on. Physically, these calculated maps are realized by either discrete and independently driven elements of active phased-array emitters [3,4,5,27,28,29,30,31,32,33,34] (such as loudspeakers in the air or piezoelectric transducers in water) or meta-pixels of passive acoustic metasurfaces [6,7,8,9,12,13,14,15,16,17,18,19,20,21,22,23,24,25,38], or a combination of both [26]. In the reconstruction process, these acoustic holographic devices with the required phase and/or amplitude distributions diffract sound waves to form the desired real image in the pre-designed image plane.

Metasurface-based acoustic holograms, which have become a research hotspot in the field of acoustic holography in recent years, can overcome the disadvantages of active phased-array holograms, such as a complex structure and circuit, high power consumption, and the need for careful calibration and tuning. That is, holography can be realized with low costs, a large aperture, on a large scale, with high precision, and with high frequency using metasurface-based acoustic holograms [18]. However, due to the structural complexity of metasurface units, their thickness is still not ultra-thin, especially in high-frequency bands, which makes expanding the application range of metasurface-based acoustic holograms difficult. To address the aforementioned issues, this paper proposes a design concept for an acoustic sieve metasurface (ASM) that can theoretically realize a zero-thickness broadband hologram. The ASM is a rigid screen with many small holes, through which the sound wave can pass to obtain the desired real image in the target plane. The arrangement rules of the holes on the ASM can be found using the Rayleigh–Sommerfeld theory [36] and optimization algorithms (a genetic algorithm [39,40] in this paper). There is no thickness limit as long as the ASM-based hologram is sufficiently rigid relative to the background medium, so it can theoretically achieve zero thickness. Furthermore, because the wave from a hole has no additional phase or amplitude modulations, the intractable modulation dispersion can be physically eliminated, allowing the proposed ASM-based hologram to potentially work in any frequency band as long as the condition of paraxial approximation is satisfied. Using the numerical simulation based on the combination of the finite element method (FEM) and the boundary element method (BEM), broadband holographic imaging be achieved, and the excellent imaging effects confirm the feasibility of our design scheme.

2. Theories and Methods

2.1. Rayleigh–Sommerfeld Theory

The complex acoustic pressure at a point in the field due to a time–harmonic radiating source, with

e^{i ω_{0} t}

time dependence, is given by the Rayleigh–Sommerfeld integral [36]

p (r) = \frac{i ρ_{0} c_{0} k_{0}}{2 π} \int_{S^{'}}^{} u (r^{'}) \frac{e^{- i k_{0} | r - r^{'} |}}{| r - r^{'} |} d S^{'}

where

i = \sqrt{- 1}

ρ_{0}

and

c_{0}

are, respectively, the density and the speed of sound in the background medium,

k_{0} = ω_{0} / c_{0}

is the wave number,

S^{'}

is the surface of the source,

u

is the particle velocity normal to the surface of the source, and

r

and

r^{'}

are the observation and source point, respectively.

2.2. Modified Genetic Algorithm

A modified genetic algorithm [39] that only contains the mutation and evolution operations for a chromosome encoding the value (0 or 1) of every grid node is used to implement the design. The hologram information is encoded into one chromosome containing

N \times N

genes, each of which is labeled by the value of a grid node. This algorithm has only mutation operations and no crossover operations. The mutation operation changes only one gene in the parent chromosome at each iteration. Offspring chromosomes are obtained after altering the gene. When calculating the offspring sound field, the sound field of the entire offspring chromosome is not calculated once, but only the influence of the changed gene on the sound field is calculated by using the Rayleigh–Sommerfeld theory, and then the influence is superimposed on the parent sound field. This superposition calculation method that reduces the amount of calculation significantly expedites an evolution operation by four orders of magnitude [39], which makes it more feasible in carrying out the design of a large-scale hologram.

2.3. Numerical Simulation

A combination of the FEM and the BEM is used for the numerical simulations and the software used is COMSOL Multiphysics. During the simulation of holographic imaging, a zero-thickness rigid boundary with a special arrangement of small holes is used to simulate the proposed hologram. The space below the rigid boundary is simulated by the FEM, and the space above the rigid boundary is simulated by the BEM. The plane where the rigid boundary is located is the coupling plane between the two methods. The plane 1/6 wavelength below the rigid boundary is the emission plane. The emitted sound wave is a plane wave with a sound pressure amplitude of 1 Pa and its propagation direction (i.e., z-direction) is perpendicular to the rigid boundary (the rigid boundary is located on the plane of z = 0 and its center is at the point of (0, 0, 0)). All boundaries in the calculation part of the FEM are set as plane-wave radiation boundaries (they can also be set as perfectly matched layers, and the calculation results of the two settings are basically the same) to achieve an anechoic environment, and the whole space in the calculation part of the BEM is set as free open space. For domains, the mesh is set to free tetrahedral, and for boundaries, the mesh is set to free triangular. The maximum element size of meshes is set to be no larger than 1/6 wavelength.

3. Results

3.1. Hologram Realized by the Array of Identical Point Sources

First, the holographic principle based on an array of identical point sources (PSs) is discussed. As shown in Figure 1, the PSs (represented by red dots) with source strength

Q_{0}

are located at the grid nodes (represented by

G_{m, n}

) of an N × N array with a distance D between two neighboring nodes, where m, n = 1, 2, 3, …, N. Every grid node is assigned a value of 1 or 0, with 1 indicating the presence of a PS and 0 indicating the absence of a PS. According to the Rayleigh–Sommerfeld theory, the sound pressure of an array of PSs at any point

(x, y, z)

in space can be described as a superposition of sound fields from all of these PSs

p (x, y, z) = \frac{i ρ_{0} c_{0} k_{0} Q_{0}}{2 π} \sum_{m = 1}^{N} \sum_{n = 1}^{N} A_{m, n} \frac{e^{- i k_{0} R_{m, n}}}{R_{m, n}}

(1)

where

A_{m, n} = 0 or 1

is the value of

G_{m, n}

R_{m, n} = \sqrt{{(x - x_{m, n})}^{2} + {(y - y_{m, n})}^{2} + z^{2}}

is the distance from

G_{m, n}

to the point

(x, y, z)

, and

(x_{m, n}, y_{m, n}, 0)

is the coordinate of

G_{m, n}

Combined with Equation (1), a desired real image (e.g., the four-pointed star pattern in this paper, as shown in Figure 2a; the black parts in the figure represent the presence of sound energy with equal magnitudes and the sound energy at these locations can be considered to be 1; the white parts represent the absence of sound energy, and the sound energy at these locations can be considered to be 0) in the target plane can be designed by using a modified genetic algorithm that only contains the mutation and evolution operations for a chromosome encoding the value (0 or 1) of every grid node. The specific operation process of the modified genetic algorithm is as follows:

(1): Randomly generate a chromosome $A_{1}$ (each subsequent generation of the chromosome is named $A_{2}$ , $A_{3}$ , …, $A_{n}$ , …). $A_{1}$ contains all $A_{m, n}$ , each $A_{m, n}$ as a gene, and the value of each gene is only 0 or 1. Moreover, let $A = A_{1}$ .
(2): Calculate the sound field $p_{1} (x, y, z_{0})$ (each subsequent generation of the sound field is named $p_{2}$ , $p_{3}$ , …, $p_{n}$ , …) generated by chromosome $A$ in the desired image plane (the distance from the hologram is $z_{0}$ ) based on Equation (1), and let $p (x, y, z_{0}) = p_{1} (x, y, z_{0})$ .
(3): Using $p (x, y, z_{0})$ and the desired sound field $p_{desired} (x, y, z_{0})$ in the image plane, the root-mean-square error (RMSE) is calculated as the value of the fitness function $F_{1}$ (each subsequent generation of the fitness function is represented by $F_{2}$ , $F_{3}$ , …, $F_{n}$ , …), and let $F = F_{1}$ .

$RMSE = \sqrt{\sum_{i = 1}^{I} [\sum_{j = 1}^{J} {[\frac{E (x_{i}, y_{j}, z_{0})}{MAX (E)} - \frac{E_{desired} (x_{i}, y_{j}, z_{0})}{MAX (E_{desired})}]}^{2}] / (I \cdot J - 1)}$

where $E (x_{i}, y_{j}, z_{0}) \propto p^{2} (x_{i}, y_{j}, z_{0})$ and $E_{desired} (x_{i}, y_{j}, z_{0}) \propto p_{desired}^{2} (x_{i}, y_{j}, z_{0}),$ respectively, represent the sound energy generated by the hologram and the sound energy of the desired image at $(x_{i}, y_{j}, z_{0})$ in the image plane; $MAX (E)$ and $MAX (E_{desired}),$ respectively, represent the maximum values of $E$ and $E_{desired}$ ; I and J, respectively, represent the numbers of samples along the x-direction and the y-direction. The RMSE is a sensitive method to show the profile differences between the sound energy generated by the hologram and the sound energy of the desired image and can be used to characterize the imaging quality in the image plane. The smaller the value of the RMSE, the closer the sound field of the image plane to the desired sound field. When the RMSE = 0, it means that the image generated by the hologram is the same as the desired image.
(4): Perform a mutation operation ( $1 \to 0$ or $0 \to 1$ ) for a random gene $A_{a, b}$ whose coordinate is $(x_{m, n}, y_{m, n}, 0)$ in the parent chromosome A to generate the offspring chromosome $A_{n}$ .
(5): Calculate the effect of this mutation for the sound field in the desired image plane using the Rayleigh–Sommerfeld theory, represented by $p_{mutated} (x, y, z_{0})$ . When the mutation operation is $1 \to 0$ , according to Equation (1), $p_{mutated} (x, y, z_{0})$ satisfies

$p_{mutated} (x, y, z_{0}) = - \frac{i ρ_{0} c_{0} k_{0} Q_{0}}{2 π} \frac{e^{- i k_{0} \sqrt{{(x_{a, b} - x)}^{2} + {(y_{a, b} - y)}^{2} + z_{0}^{2}}}}{\sqrt{{(x_{a, b} - x)}^{2} + {(y_{a, b} - y)}^{2} + z_{0}^{2}}}$

When the mutation operation is $0 \to 1$ , $p_{mutated} (x, y, z_{0})$ satisfies

$p_{mutated} (x, y, z_{0}) = \frac{i ρ_{0} c_{0} k_{0} Q_{0}}{2 π} \frac{e^{- i k_{0} \sqrt{{(x_{a, b} - x)}^{2} + {(y_{a, b} - y)}^{2} + z_{0}^{2}}}}{\sqrt{{(x_{a, b} - x)}^{2} + {(y_{a, b} - y)}^{2} + z_{0}^{2}}}$
(6): Superimpose this effect on the parent sound field $p (x, y, z_{0})$ to get the offspring sound field $p_{n} (x, y, z_{0})$ in the desired image plane, as follow

$p_{n} (x, y, z_{0}) = p (x, y, z_{0}) + p_{mutated} (x, y, z_{0})$
(7): Calculate the value (RMSE) of the offspring fitness function $F_{n}$ by $p_{n} (x, y, z_{0})$ and $p_{desired} (x, y, z_{0})$ .
(8): Compare the magnitudes of $F_{n}$ and $F$ , and take the smaller one as the new parent fitness function. The chromosome with the smaller fitness function is taken as the new parent chromosome and the sound field with the smaller fitness function is taken as the new parent sound field. That is, if $F_{n} < F$ , then let $F = F_{n}$ , $A = A_{n}$ , $p (x, y, z_{0}) = p_{n} (x, y, z_{0})$ ; if $F < F_{n}$ , then let $F = F$ , $A = A$ , $p (x, y, z_{0}) = p (x, y, z_{0})$ . If $F$ is always less than $F_{n}$ for a certain number (such as 100,000 times) of iterations, then record chromosome $A$ and end the iteration.
(9): If the iteration does not end, go back to step 4.

Based on the desired image shown in Figure 2a, the overall map (as shown in Figure 2b) of the encoded PSs can be designed by the modified genetic algorithm at

z_{0} = 120 λ_{0}

, where

λ_{0} = c_{0} / f_{0} = 10 mm

is the wavelength (the background medium is chosen as air with a sound velocity of

c_{0} = 340 m / s

and a density of

ρ_{0} = 1.2 kg / m^{3}

, and the radiation frequency is set at

f_{0} = 34 kHz

). The distance between grid nodes is

D = λ_{0}

, and there are 50 grid nodes in the horizontal and vertical directions, so the size of the hologram is

50 λ_{0} \times 50 λ_{0}

, which is the same as the size of the desired image. The BEM-based numerical simulation is used to confirm the performance of the PS array. Figure 2c shows the simulation results, which show a clear image that closely matches the desired image. The RMSE between this holographic image and the desired image shown in Figure 2a is 0.1737.

3.2. Theoretical Zero-Thickness Holograms Based on Acoustic Sieve Metasurfaces

To approximate the PS array, an acoustic sieve metasurface (ASM) that is a zero-thickness rigid screen with many small holes of subwavelength size is introduced. A sound field similar to that of a PS can be created by sound waves passing through a small hole. Despite the fact that a small hole in a rigid screen cannot be considered an ideal PS because it must radiate an ideal spherical wave, it is found that the sound field diffracted from a small hole is nearly identical to the sound field diffracted from an ideal PS (see the next paragraph for specific evidence), which is important because the reconstruction of a holographic image is largely dependent on the fields of diffraction units (such as holes or PSs) in the image plane. Furthermore, the transmitted wave from the small hole propagates forward, whereas the reflected wave has no bearing on image reconstruction. As a result, the concept of ASMs is introduced for the realization of holograms.

First, as shown in Figure 3a, the lateral profile of the sound field from a PS in the middle of a circular rigid boundary with a diameter of

30 λ_{0} is investigated

. On the red dashed line in Figure 3a, the lateral (x-y plane) energy profiles in the z-cut plane of

z = 4 λ_{0}

are calculated. This distance (

4 λ_{0}

) is beyond the evanescent region while remaining small enough that the simulation model does not consume excessive computational resources or memory space. Then, in the same model, the PS is replaced with holes of various diameters (

d = 0.01 λ_{0}, 0.02 λ_{0}, 0.03 λ_{0}, \dots, λ_{0}

) to allow plane sound waves to pass through. The lateral energy profiles of the same location (the red dashed line) generated by plane sound waves passing through holes of various diameters are calculated using the FEM-based numerical simulation. The lateral normalized energy profiles

E_{PS / hole} (x_{k}, 0, 4 λ_{0}) / E_{PS / hole} (0, 0, 4 λ_{0})

of the ideal PS and holes of various diameters are shown in Figure 3b. As can be seen, the smaller the hole’s diameter, the closer the energy profile approaches to the ideal PS. By calculating their RMSEs (shown by blue circles in Figure 3c), the degree of approximation between PS and holes of various diameters can be quantitatively verified

RMSE = \sqrt{\sum_{k = 1}^{K} {[\frac{E_{hole} (x_{k}, 0, 4 λ_{0})}{E_{hole} (0, 0, 4 λ_{0})} - \frac{E_{PS} (x_{k}, 0, 4 λ_{0})}{E_{PS} (0, 0, 4 λ_{0})}]}^{2} / (K - 1)}

where

E_{hole} (x_{k}, 0, 4 λ_{0}) \propto p_{hole}^{2} (x_{k}, 0, 4 λ_{0})

and

E_{PS} (x_{k}, 0, 4 λ_{0}) \propto p_{PS}^{2} (x_{k}, 0, 4 λ_{0})

are the energy profiles of the hole and ideal PS, respectively, and

K = 1500

is the sampling number in the x-direction. The calculated RMSE values (blue circles) and relative transmitted energies (red dots) for holes of various diameters are shown in Figure 2c. Each blue circle represents an RMSE between

E_{hole} (x_{k}, 0, 4 λ_{0}) / MAX (E_{hole})

(color curves in Figure 2b) from a d-diameter hole and

E_{PS} (x_{k}, 0, 4 λ_{0}) / MAX (E_{PS})

(black curve in Figure 2b) from the ideal PS. Each red dot represents the logarithm of the simulated

E_{hole} (0, 0, 4 λ_{0}) / {E_{hole} (0, 0, 4 λ_{0}) |}_{d = λ_{0}}

from a d-diameter hole. The RMSE decreases as the hole size decreases, as shown in the figure, implying that the smaller the hole size, the better the effect of the hole simulating the PS. Furthermore, the figure shows that as the hole size decreases, the transmitted energy also decreases, and when the hole diameter is less than 0.2

λ_{0}

, the energy decay rate is significantly accelerated. As a result,

d = 0.2 λ_{0}

(i.e., 2 mm) is chosen as the diameter value of the ASM holes to balance the effect and efficiency.

Then, instead of the PSs in Figure 2b, the holes with a diameter of

d = 0.2 λ_{0}

are employed to construct the ASM-based hologram shown in Figure 4a. The light gray area is the zero-thickness rigid screen and the small circles are the holes. Figure 4b shows the image of the ASM-based hologram calculated using numerical simulation based on a combination of the FEM and the BEM in the image plane. In the numerical simulation, all of the acoustic and geometric parameters of the holographic imaging system based on the ASM are the same as those of the holographic imaging system based on the PS array. The RMSE between this holographic image and the desired image shown in Figure 2a is 0.1755. As can be seen, this image closely resembles the desired image.

3.3. Theoretical Zero-Thickness Broadband Holograms

It is worth noting that the wave from a hole has no additional phase or amplitude modulations, which can physically remove the intractable modulation dispersion. As a result, even though it is designed for a specific frequency (

f_{0} = 34 kHz

), the proposed ASM should be feasible for a broadband hologram. Although the modulation dispersion has been removed, the intrinsic propagation dispersion [40] is provided by the propagation kernel

e^{- i k_{0} R_{m, n}} / R_{m, n}

in Equation (1) describing the diffraction of a PS located at

G_{m, n}

. In principle, propagation dispersion cannot be eliminated, but it primarily affects a hologram’s imaging position and does not result in significant distortions or aberrations in the reconstructed holographic image [40,41,42]. The imaging distance (

z

) from a hologram varies with the radiation frequency (

f

) and follows the mathematical relationship [42]

\frac{2 π f}{c_{0}} (\sqrt{{(x_{m, n} - x)}^{2} + {(y_{m, n} - y)}^{2} + z^{2}} - z) = \frac{2 π f_{0}}{c_{0}} (\sqrt{{(x_{m, n} - x)}^{2} + {(y_{m, n} - y)}^{2} + z_{0}^{2}} - z_{0})

Under the paraxial approximation, the mathematical relationship can be simply expressed as [42]

\frac{z}{f} = \frac{z_{0}}{f_{0}}

(2)

indicating that as the frequency increases, the imaging distance increases linearly.

The sound energy distributions in the image planes (according to Equation (2), the imaging distances (

z

) can be calculated as

60 λ_{0}

80 λ_{0}

100 λ_{0}

144 λ_{0}

180 λ_{0}

and

240 λ_{0}

, respectively) of the ASM-based hologram are calculated using numerical simulation based on the combination of the FEM and the BEM when the incident frequencies (

f

) are, respectively,

f_{0} / 2

f_{0} / 1.5

f_{0} / 1.2

1.2 f_{0}

1.5 f_{0},

and

2 f_{0}

, as shown in Figure 5a–f. The figures show that the imaging effects in the six image planes are quite good. The RMSEs between these holographic images and the desired image shown in Figure 2a are 0.1942, 0.1813, 0.1771, 0.1743, 0.1849, and 0.1839, respectively. The graph (Figure 6) of the RMSE against frequency further shows that in the studied frequency range, all the RMSE values between the holographic images and the desired image shown in Figure 2a are modest and nearly constant.

4. Discussion

4.1. Selection of Geometric Parameters

Regarding the setting of the distance between the grid nodes, according to the Nyquist sampling theorem, the sampling condition in the calculation of the Rayleigh–Sommerfeld formula should satisfy [43]

Δ x \leq \frac{λ \sqrt{z^{2} + Δ L^{2} / 2}}{Δ L}

where,

Δ x

is the sampling interval, which is the distance between the grid nodes (D), and

Δ L

is the width of the diffraction field, which can be approximated as the diagonal length (

50 \sqrt{2} λ_{0}

) of the proposed hologram. Substitute the relevant data into the above equation to obtain

D = Δ x \leq 23.9 mm = 2.39 λ_{0}

Regarding the setting of the image size, if the desired image size is set too small, the image detail size may be smaller than the minimum resolution distance (according to the Rayleigh criterion, the minimum resolution distance should be

1.22 λ z / a

, where

a

can be approximately regarded as the diagonal length of the proposed hologram), resulting in a loss of image details. If the desired image size is set too large, the paraxial approximation will not be satisfied or the information content [21] of the hologram will be less than that of the image. Therefore, in general, for acoustic holographic imaging, the desired image size is set to be equal to or similar to the hologram size in most research [3,15,16,17,18,19,21,22,23,24,25,26,27,28].

4.2. Robustness

To verify the robustness of the proposed hologram, three cases that may cause errors are calculated using numerical simulation based on the combination of the FEM and the BEM. These three cases are, respectively, small displacements of holes from their ideal places, a dispersion of hole size, and an incidence wave with a wrong angle. In the case of small displacements of holes from their ideal places, every hole position is added a random error within

\pm 0.1 λ_{0}

in the x-direction and y-direction, respectively. Figure 7a–c show the calculated images for the three frequencies of

f_{0} / 2

f_{0},

and

2 f_{0}

in this case. It can be seen from the figures that the images are almost identical to Figure 4c and Figure 5a,f but without the error. The RMSEs between the images and the desired image (Figure 2a) are 0.2012, 0.1788, and 0.1841, respectively. The RMSEs increase only slightly (rates of increase are less than 3.7%) compared to the non-error case. In the case of the dispersion of hole size, a random error within

\pm 0.1 d

is added to the diameter of every hole. Figure 7d–f show the calculated images for the three frequencies of

f_{0} / 2

f_{0},

and

2 f_{0}

in this case. It can be seen from the figures that the images are almost identical to Figure 4c and Figure 5a,f but without the error. The RMSEs between the images and the desired image (Figure 2a) are 0.2009, 0.1780, and 0.1868, respectively. The RMSEs increase only slightly (rates of increase are less than 3.5%) compared to the non-error case. In the case of the incident wave with a wrong angle, the plane wave that was originally incident parallel to the z-axis is inclined to the y-direction by 5°. Figure 7g–i show the calculated images for the three frequencies of

f_{0} / 2

f_{0},

and

2 f_{0}

in this case. Since the incident direction of the plane sound wave is inclined, the position of the image is shifted by

z \tan 5 °

. The RMSEs between the images and the desired image (Figure 2a) are 0.2062, 0.1873, and 0.1868, respectively. The increases in RMSEs are also not significant (rates of increase are less than 6.8%) compared to the non-error case. To sum up, the proposed hologram offers a high level of robustness.

4.3. Actual Situation

To verify the effectiveness of the proposed hologram design method in practical applications, based on the combination of the FEM and the BEM, holographic imaging for the three frequencies of

f_{0} / 2

f_{0},

and

2 f_{0}

is simulated in a case where the hologram is fabricated with a

0.1 λ_{0}

-thick steel plate with holes, as shown in Figure 8a–c. It can be seen from the figures that the images are almost identical to Figure 4c and Figure 5a,f of the ideal rigid screen. The RMSEs between the images and the desired image (Figure 2a) are 0.2008, 0.1780, and 0.1853, respectively. The RMSEs increase only slightly (rates of increase are less than 3.4%) compared to the case of the ideal rigid screen.

It is worth noting that, the proposed hologram design method may not be well-suited for the background medium of water in practical applications, because most materials are not rigid enough relative to water.

5. Conclusions

In conclusion, a theoretical zero-thickness broadband hologram is proposed and comprehensively investigated. A design approach for an ASM capable of achieving the zero-thickness broadband hologram is offered in theory. The ASM is a hole-filled zero-thickness rigid screen. A genetic algorithm and the Rayleigh–Sommerfeld theory can be used to obtain the arrangement rules of holes in the ASM. It should have no thickness restrictions as long as the ASM-based hologram is sufficiently rigid in comparison to the background medium, and it can theoretically approach zero thickness. The intractable modulation dispersion could also be physically avoided because the wave from a hole has no extra phase or amplitude modulations, allowing the proposed ASM-based hologram to possibly function in any frequency band as long as the condition of paraxial approximation is satisfied. However, because propagation dispersion cannot be eliminated, one disadvantage is that the imaging distance from a hologram varies with the radiation frequency. Broadband holographic imaging with an excellent effect is achieved by using numerical simulation based on a combination of the FEM and the BEM. The proposed theoretical zero-thickness broadband hologram could open up new possibilities for expanding the application scenarios of acoustic holography.

Author Contributions

Conceptualization, S.Z.; Data curation, Y.T.; Formal analysis, Y.T.; Funding acquisition, Y.T. and J.G.; Investigation, Y.T. and S.Z.; Project administration, S.Z. and J.G.; Resources, J.G.; Software, Y.T.; Supervision, J.G.; Validation, Q.L. and G.Y.; Visualization, Q.L. and G.Y.; Writing—original draft, Y.T. and S.Z.; Writing—review and editing, Y.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (grant numbers 11904221, 12104284, 12004237, 12034005, and 11727813), the China National Postdoctoral Program for Innovative Talents (grant number BX20190193), the China Postdoctoral Science Foundation (grant numbers 2019M663612 and 2020M683416), Youth Talent Promotion Project of Xi’an Association for Science Technology (grant number 095920221302), the Fundamental Research Funds for the Central Universities (grant numbers GK202103017, TD2020009Z, and 2020TS023), and the High-level Talent Fund of Henan University of Technology (grant number 31401242).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to extend our deep gratitude to all those who have offered cordial and selfless support in writing this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

Mueller, R.K. Acoustic holography. Proc. IEEE 1971, 59, 1319–1335. [Google Scholar] [CrossRef]
Ahmed, M.; Wang, K.Y.; Metherell, A.F. Holography and its application to acoustic imaging. Proc. IEEE 1979, 67, 466–483. [Google Scholar] [CrossRef]
Hertzberg, Y.; Navon, G. Bypassing absorbing objects in focused ultrasound using computer generated holographic technique. Med. Phys. 2011, 38, 6407–6415. [Google Scholar] [CrossRef]
Sapozhnikov, O.A.; Tsysar, S.A.; Khokhlova, V.A.; Kreider, W. Acoustic holography as a metrological tool for characterizing medical ultrasound sources and fields. J. Acoust. Soc. Am. 2015, 138, 1515–1532. [Google Scholar] [CrossRef] [Green Version]
Kreider, W.; Yuldashev, P.V.; Sapozhnikov, O.A.; Farr, N.; Partanen, A.; Bailey, M.R.; Khokhlova, V.A. Characterization of a multi-element clinical HIFU system using acoustic holography and nonlinear modeling. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2013, 60, 1683–1698. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Jimenez-Gambin, S.; Jimenez, N.; Benlloch, J.M.; Camarena, F. Holograms to Focus Arbitrary Ultrasonic Fields through the Skull. Phys. Rev. Appl. 2019, 12, 014016. [Google Scholar] [CrossRef]
Melde, K.; Choi, E.; Wu, Z.G.; Palagi, S.; Qiu, T.; Fischer, P. Acoustic Fabrication via the Assembly and Fusion of Particles. Adv. Mater. 2018, 30, 1704507. [Google Scholar] [CrossRef]
Kruizinga, P.; van der Meulen, P.; Fedjajevs, A.; Mastik, F.; Springeling, G.; de Jong, N.; Bosch, J.G.; Leus, G. Compressive 3D ultrasound imaging using a single sensor. Sci. Adv. 2017, 3, e1701423. [Google Scholar] [CrossRef] [Green Version]
Ma, Z.C.; Holle, A.W.; Melde, K.; Qiu, T.; Poeppel, K.; Kadiri, V.M.; Fischer, P. Acoustic Holographic Cell Patterning in a Biocompatible Hydrogel. Adv. Mater. 2020, 32, 1904181. [Google Scholar] [CrossRef] [Green Version]
Baudoin, M.; Thomas, J.L.; Al Sahely, R.; Gerbedoen, J.C.; Gong, Z.X.; Sivery, A.; Matar, O.B.; Smagin, N.; Favreau, P.; Vlandas, A. Spatially selective manipulation of cells with single-beam acoustical tweezers. Nat. Commun. 2020, 11, 4244. [Google Scholar] [CrossRef]
Baudoin, M.; Gerbedoen, J.C.; Riaud, A.; Matar, O.B.; Smagin, N.; Thomas, J.L. Folding a focalized acoustical vortex on a flat holographic transducer: Miniaturized selective acoustical tweezers. Sci. Adv. 2019, 5, eaav1967. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Wang, H.; Gao, W.; Zhu, R.; Wang, Z.; Xu, Z.; Zheng, B. Ultrathin Acoustic Metasurface Holograms with Arbitrary Phase Control. Appl. Sci. 2019, 9, 3585. [Google Scholar] [CrossRef] [Green Version]
Jiménez, N.; Groby, J.-P.; Romero-García, V. Spiral sound-diffusing metasurfaces based on holographic vortices. Sci. Rep. 2021, 11, 10217. [Google Scholar] [CrossRef] [PubMed]
Zhang, H.; Zhang, W.; Liao, Y.; Zhou, X.; Li, J.; Hu, G.; Zhang, X. Creation of acoustic vortex knots. Nat. Commun. 2020, 11, 3956. [Google Scholar] [CrossRef]
Zhang, J.; Yang, Y.; Zhu, B.; Li, X.; Jin, J.; Chen, Z.; Chen, Y.; Zhou, Q. Multifocal point beam forming by a single ultrasonic transducer with 3D printed holograms. Appl. Phys. Lett. 2018, 113, 243502. [Google Scholar] [CrossRef]
Tian, Y.; Zuo, S.; Cheng, Y.; Liu, X. Phase-controlled acoustic metasurface with high efficiency and its applications. J. Appl. Acoust. 2018, 37, 691–700. [Google Scholar] [CrossRef]
Tian, Y.; Wei, Q.; Cheng, Y.; Liu, X. Acoustic holography based on composite metasurface with decoupled modulation of phase and amplitude. Appl. Phys. Lett. 2017, 110, 191901. [Google Scholar] [CrossRef]
Brown, M.D.; Cox, B.; Treeby, B.E. Stackable acoustic holograms. Appl. Phys. Lett. 2020, 116, 261901. [Google Scholar] [CrossRef]
Zhu, Y.; Hu, J.; Fan, X.; Yang, J.; Liang, B.; Zhu, X.; Cheng, J. Fine manipulation of sound via lossy metamaterials with independent and arbitrary reflection amplitude and phase. Nat. Commun. 2018, 9, 1632. [Google Scholar] [CrossRef] [Green Version]
Xie, Y.; Shen, C.; Wang, W.; Li, J.; Suo, D.; Popa, B.-I.; Jing, Y.; Cummer, S.A. Acoustic Holographic Rendering with Two-dimensional Metamaterial-based Passive Phased Array. Sci. Rep. 2016, 6, 35437. [Google Scholar] [CrossRef]
Melde, K.; Mark, A.G.; Qiu, T.; Fischer, P. Holograms for acoustics. Nature 2016, 537, 518–522. [Google Scholar] [CrossRef] [PubMed]
Zhu, Y.; Assouar, B. Systematic design of multiplexed-acoustic-metasurface hologram with simultaneous amplitude and phase modulations. Phys. Rev. Mater. 2019, 3, 045201. [Google Scholar] [CrossRef]
Fan, S.W.; Zhu, Y.F.; Cao, L.Y.; Wang, Y.F.; Chen, A.L.; Merkel, A.; Wang, Y.S.; Assouar, B. Broadband tunable lossy metasurface with independent amplitude and phase modulations for acoustic holography. Smart Mater. Struct. 2020, 29, 105038. [Google Scholar] [CrossRef]
Brown, M.D. Phase and amplitude modulation with acoustic holograms. Appl. Phys. Lett. 2019, 115, 053701. [Google Scholar] [CrossRef] [Green Version]
Bakhtiari-Nejad, M.; Elnahhas, A.; Hajj, M.R.; Shahab, S. Acoustic holograms in contactless ultrasonic power transfer systems: Modeling and experiment. J. Appl. Phys. 2018, 124, 244901. [Google Scholar] [CrossRef]
Cox, L.; Melde, K.; Croxford, A.; Fischer, P.; Drinkwater, B.W. Acoustic Hologram Enhanced Phased Arrays for Ultrasonic Particle Manipulation. Phys. Rev. Appl. 2019, 12, 064055. [Google Scholar] [CrossRef] [Green Version]
Fushimi, T.; Yamamoto, K.; Ochiai, Y. Acoustic hologram optimisation using automatic differentiation. Sci. Rep. 2021, 11, 12678. [Google Scholar] [CrossRef]
Morales, R.; Ezcurdia, I.; Irisarri, J.; Andrade, M.A.B.; Marzo, A. Generating Airborne Ultrasonic Amplitude Patterns Using an Open Hardware Phased Array. Appl. Sci. 2021, 11, 2981. [Google Scholar] [CrossRef]
Hirayama, R.; Martinez Plasencia, D.; Masuda, N.; Subramanian, S. A volumetric display for visual, tactile and audio presentation using acoustic trapping. Nature 2019, 575, 320–323. [Google Scholar] [CrossRef]
Fushimi, T.; Marzo, A.; Drinkwater, B.W.; Hill, T.L. Acoustophoretic volumetric displays using a fast-moving levitated particle. Appl. Phys. Lett. 2019, 115, 064101. [Google Scholar] [CrossRef]
Hoshi, T.; Takahashi, M.; Iwamoto, T.; Shinoda, H. Noncontact Tactile Display Based on Radiation Pressure of Airborne Ultrasound. IEEE Trans. Haptics 2010, 3, 155–165. [Google Scholar] [CrossRef] [PubMed]
Long, B.; Seah, S.A.; Carter, T.; Subramanian, S. Rendering Volumetric Haptic Shapes in Mid-Air using Ultrasound. Acm Trans. Graph. 2014, 33, 181. [Google Scholar] [CrossRef] [Green Version]
Marzo, A.; Drinkwater, B.W. Holographic acoustic tweezers. Proc. Natl. Acad. Sci. USA 2019, 116, 84–89. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Marzo, A.; Seah, S.A.; Drinkwater, B.W.; Sahoo, D.R.; Long, B.; Subramanian, S. Holographic acoustic elements for manipulation of levitated objects. Nat. Commun. 2015, 6, 8661. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Lin, Q.; Wang, J.; Cai, F.; Zhang, R.; Zhao, D.; Xia, X.; Wang, J.; Zheng, H. A deep learning approach for the fast generation of acoustic holograms. J. Acoust. Soc. Am. 2021, 149, 2312–2322. [Google Scholar] [CrossRef]
Ebbini, E.S.; Cain, C.A. Multiple-focus ultrasound phased-array pattern synthesis: Optimal driving-signal distributions for hyperthermia. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 1989, 36, 540–548. [Google Scholar] [CrossRef] [PubMed]
Ibbini, M.S.; Cain, C.A. A field conjugation method for direct synthesis of hyperthermia phases-array heating patterns. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 1989, 36, 3–9. [Google Scholar] [CrossRef]
Zhang, J.; Tian, Y.; Cheng, Y.; Liu, X. Acoustic holography using composite metasurfaces. Appl. Phys. Lett. 2020, 116, 030501. [Google Scholar] [CrossRef]
Huang, K.; Liu, H.; Garcia-Vidal, F.J.; Hong, M.H.; Luk’yanchuk, B.; Teng, J.H.; Qiu, C.W. Ultrahigh-capacity non-periodic photon sieves operating in visible light. Nat. Commun. 2015, 6, 7059. [Google Scholar] [CrossRef] [Green Version]
Huang, K.; Liu, H.; Si, G.Y.; Wang, Q.; Lin, J.; Teng, J.H. Photon-nanosieve for ultrabroadband and large-angle-of-view holograms. Laser Photonics Rev. 2017, 11, 1700025. [Google Scholar] [CrossRef]
Huang, K.; Dong, Z.G.; Mei, S.T.; Zhang, L.; Liu, Y.J.; Liu, H.; Zhu, H.B.; Teng, J.H.; Luk’yanchuk, B.; Yang, J.K.W.; et al. Silicon multi-meta-holograms for the broadband visible light. Laser Photonics Rev. 2016, 10, 500–509. [Google Scholar] [CrossRef]
Huang, L.L.; Chen, X.Z.; Muhlenbernd, H.; Zhang, H.; Chen, S.M.; Bai, B.F.; Tan, Q.F.; Jin, G.F.; Cheah, K.W.; Qiu, C.W.; et al. Three-dimensional optical holography using a plasmonic metasurface. Nat. Commun. 2013, 4, 2808. [Google Scholar] [CrossRef] [Green Version]
Li, J.C.; Peng, Z.J.; Fu, Y.C. Diffraction transfer function and its calculation of classic diffraction formula. Opt. Commun. 2007, 280, 243–248. [Google Scholar] [CrossRef]

Figure 1. Schematic diagram of a hologram based on an array of identical PSs. The PSs (represented by red dots) are at the grid nodes (represented by

G_{m, n}

) of an

N \times N

array with a distance d between two neighboring nodes, where

m, n = 1, 2, 3, \dots, N

. The grid node

G_{m, n}

takes the value

A_{m, n} = 0 or 1

, where 1 means a PS exists in this grid node and 0 means no PS. According to the Rayleigh–Sommerfeld theory, the diffraction pattern in the image plane is the superposition of the sound waves generated by these PSs.

Figure 1. Schematic diagram of a hologram based on an array of identical PSs. The PSs (represented by red dots) are at the grid nodes (represented by

G_{m, n}

) of an

N \times N

array with a distance d between two neighboring nodes, where

m, n = 1, 2, 3, \dots, N

. The grid node

G_{m, n}

takes the value

A_{m, n} = 0 or 1

Figure 2. Holographic imaging based on the PS array. The background medium is chosen as air with a sound velocity

c_{0} = 340 m / s

and a density

ρ_{0} = 1.2 kg / m^{3}

, and the radiation frequency is set at

f_{0} = 34 kHz

, so the wavelength is

λ_{0} = c_{0} / f_{0} = 10 mm

. (a) Desired image (a four-pointed star). The black parts in the figure represent the presence of sound energy with equal magnitudes (the sound energy at these locations can be considered to be 1), and the white parts represent the absence of sound energy (the sound energy at these locations can be considered to be 0); (b) The overall map of the hologram based on the PS array designed by the modified genetic algorithm. Each black dot represents a PS. The distance between grid nodes is

D = λ_{0}

and the distance of the desired image plane from the hologram is

z_{0} = 120 λ_{0}

; (c) Normalized energy distribution in the image plane.

Figure 2. Holographic imaging based on the PS array. The background medium is chosen as air with a sound velocity

c_{0} = 340 m / s

and a density

ρ_{0} = 1.2 kg / m^{3}

, and the radiation frequency is set at

f_{0} = 34 kHz

, so the wavelength is

λ_{0} = c_{0} / f_{0} = 10 mm

D = λ_{0}

and the distance of the desired image plane from the hologram is

z_{0} = 120 λ_{0}

; (c) Normalized energy distribution in the image plane.

Figure 3. Mimicking of a PS by a small hole. (a) Illustration for a simulation model of PS (or holes of various diameters) radiation. The green part is a rigid circular boundary with a diameter of

30 λ_{0}

; (b) Normalized energy profiles

E_{hole} (x_{k}, 0, 4 λ_{0}) / E_{hole} (0, 0, 4 λ_{0})

of holes of various diameters calculated using numerical simulation based on the model of (a). The various pseudocolors of the curves represent the various hole diameters d. The black curve is the normalized energy profile

E_{PS} (x_{k}, 0, 4 λ_{0}) / E_{PS} (0, 0, 4 λ_{0})

radiating from an ideal PS; (c) A quantitative evaluation of radiations from holes of various diameters. Each blue circle represents an RMSE between

E_{hole} (x_{k}, 0, 4 λ_{0}) / MAX (E_{hole})

(color curves in (b)) from a d-diameter hole and

E_{PS} (x_{k}, 0, 4 λ_{0}) / MAX (E_{PS})

(black curve in (b)) from the ideal PS. Each red dot represents the logarithm of the simulated

E_{hole} (0, 0, 4 λ_{0}) / {E_{hole} (0, 0, 4 λ_{0}) |}_{d = λ_{0}}

from a d-diameter hole.

30 λ_{0}

; (b) Normalized energy profiles

E_{hole} (x_{k}, 0, 4 λ_{0}) / E_{hole} (0, 0, 4 λ_{0})

E_{PS} (x_{k}, 0, 4 λ_{0}) / E_{PS} (0, 0, 4 λ_{0})

radiating from an ideal PS; (c) A quantitative evaluation of radiations from holes of various diameters. Each blue circle represents an RMSE between

E_{hole} (x_{k}, 0, 4 λ_{0}) / MAX (E_{hole})

(color curves in (b)) from a d-diameter hole and

E_{PS} (x_{k}, 0, 4 λ_{0}) / MAX (E_{PS})

(black curve in (b)) from the ideal PS. Each red dot represents the logarithm of the simulated

E_{hole} (0, 0, 4 λ_{0}) / {E_{hole} (0, 0, 4 λ_{0}) |}_{d = λ_{0}}

from a d-diameter hole.

Figure 4. Holographic imaging of an ASM. (a) The overall map of the ASM-based hologram based on small holes. The light gray area is the zero-thickness rigid screen and the small circles are the holes with a diameter of

d = 0.2 λ_{0}

; (b) Normalized energy distribution in the image plane calculated by numerical simulation from the hologram shown in (a). In the numerical simulation, all of the acoustic and geometric parameters of the holographic imaging system based on the ASM are the same as those of the holographic imaging system based on the PS array.

d = 0.2 λ_{0}

Figure 5. Acoustic broadband holographic imaging. (a–f) The normalized sound energy distributions in the image planes of the ASM-based hologram are calculated using the numerical simulation when the incident frequencies (

f

) are

f_{0} / 2

f_{0} / 1.5

f_{0} / 1.2

1.2 f_{0}

1.5 f_{0,}

and

2 f_{0}

, respectively. The imaging distances (

z

) are

60 λ_{0}

80 λ_{0}

100 λ_{0}

144 λ_{0}

180 λ_{0},

and

240 λ_{0}

, respectively.

f

) are

f_{0} / 2

f_{0} / 1.5

f_{0} / 1.2

1.2 f_{0}

1.5 f_{0,}

and

2 f_{0}

, respectively. The imaging distances (

z

) are

60 λ_{0}

80 λ_{0}

100 λ_{0}

144 λ_{0}

180 λ_{0},

and

240 λ_{0}

, respectively.

Figure 6. Graph of RMSE against frequency. The RMSE values between the holographic images and the desired image shown in Figure 2a are calculated using numerical simulation based on the combination of the FEM and the BEM. The studied frequency range is from

f_{0} / 2

2 f_{0}

f_{0} / 2

2 f_{0}

Figure 7. Holographic images of the cases that may cause errors. (a–c) The normalized sound energy distributions in the image planes of the ASM-based hologram with a small displacement of holes from their ideal places are calculated using numerical simulation when the incident frequencies are

f_{0} / 2

f_{0},

and

2 f_{0}

, respectively. (d–f) The normalized sound energy distributions in the image planes of the ASM-based hologram with a dispersion of hole size are calculated using numerical simulation when the incident frequencies are

f_{0} / 2

f_{0},

and

2 f_{0}

, respectively. (g–i) The normalized sound energy distributions in the image planes of the ASM-based hologram for incidence wave with a 5° declination are calculated using numerical simulation when the incident frequencies are

f_{0} / 2

f_{0},

and

2 f_{0}

, respectively. Since the incident direction of the plane sound wave is inclined, the position of the image is shifted by

z \tan 5 °

f_{0} / 2

f_{0},

and

2 f_{0}

f_{0} / 2

f_{0},

and

2 f_{0}

f_{0} / 2

f_{0},

and

2 f_{0}

, respectively. Since the incident direction of the plane sound wave is inclined, the position of the image is shifted by

z \tan 5 °

Figure 8. Holographic images of

0.1 λ_{0}

-thick steel plate with holes. (a–c) The normalized sound energy distributions in the image planes of the hologram fabricated with

0.1 λ_{0}

-thick steel plate with holes are calculated using numerical simulation when the incident frequencies are

f_{0} / 2

f_{0},

and

2 f_{0}

, respectively.

Figure 8. Holographic images of

0.1 λ_{0}

-thick steel plate with holes. (a–c) The normalized sound energy distributions in the image planes of the hologram fabricated with

0.1 λ_{0}

-thick steel plate with holes are calculated using numerical simulation when the incident frequencies are

f_{0} / 2

f_{0},

and

2 f_{0}

, respectively.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Tian, Y.; Zuo, S.; Lv, Q.; Yin, G.; Guo, J. Theoretical Zero-Thickness Broadband Holograms Based on Acoustic Sieve Metasurfaces. Appl. Sci. 2022, 12, 6453. https://doi.org/10.3390/app12136453

AMA Style

Tian Y, Zuo S, Lv Q, Yin G, Guo J. Theoretical Zero-Thickness Broadband Holograms Based on Acoustic Sieve Metasurfaces. Applied Sciences. 2022; 12(13):6453. https://doi.org/10.3390/app12136453

Chicago/Turabian Style

Tian, Ye, Shuyu Zuo, Qian Lv, Guanjun Yin, and Jianzhong Guo. 2022. "Theoretical Zero-Thickness Broadband Holograms Based on Acoustic Sieve Metasurfaces" Applied Sciences 12, no. 13: 6453. https://doi.org/10.3390/app12136453

APA Style

Tian, Y., Zuo, S., Lv, Q., Yin, G., & Guo, J. (2022). Theoretical Zero-Thickness Broadband Holograms Based on Acoustic Sieve Metasurfaces. Applied Sciences, 12(13), 6453. https://doi.org/10.3390/app12136453

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu

Theoretical Zero-Thickness Broadband Holograms Based on Acoustic Sieve Metasurfaces

Abstract

1. Introduction

2. Theories and Methods

2.1. Rayleigh–Sommerfeld Theory

2.2. Modified Genetic Algorithm

2.3. Numerical Simulation

3. Results

3.1. Hologram Realized by the Array of Identical Point Sources

3.2. Theoretical Zero-Thickness Holograms Based on Acoustic Sieve Metasurfaces

3.3. Theoretical Zero-Thickness Broadband Holograms

4. Discussion

4.1. Selection of Geometric Parameters

4.2. Robustness

4.3. Actual Situation

5. Conclusions

Author Contributions

Funding

Institutional Review Board Statement

Informed Consent Statement

Data Availability Statement

Acknowledgments

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI