Issue 
Acta Acust.
Volume 5, 2021



Article Number  39  
Number of page(s)  12  
Section  Computational and Numerical Acoustics  
DOI  https://doi.org/10.1051/aacus/2021030  
Published online  01 September 2021 
Scientific Article
Feasibility analysis for active near/far field acoustic pattern synthesis in free space and shallow water environments
^{1}
Department of Electrical and Computer Engineering, University of Houston, 4800 Calhoun Rd, Houston, TX 77004, USA
^{2}
Institute of Mathematical Sciences and Physics, University of the Philippines Los Baños, Victoria M. Ela Ave., Los Baños, Laguna 4030, Philippines
^{3}
Department of Mathematics, University of Houston, 4800 Calhoun Rd, Houston, TX 77004, USA
^{*} Corresponding author: jchen82@central.uh.edu
Received:
19
January
2021
Accepted:
4
August
2021
In this paper, a detailed sensitivity and feasibility analysis of the active manipulation scheme for scalar Helmholtz fields proposed in our previous works, in both free space and constantdepth homogeneous ocean environments, is presented. We apply the method of moments (MoM) together with Tikhonov regularization with the Morozov discrepancy principle to investigate the effects of varying the problem parameters to the accuracy and feasibility of the proposed active field control strategy. We discuss the feasibility of the active scheme (with respect to power budget, control accuracy and process error) as a function of the frequency, the distance between the control region and the source, the mutual distance between the control regions, and the size of the control region. Process error is considered as well to investigate the possibility of an accurate active control in the presence of manufacturing or feeding noise. The numerical simulations show the accuracy of the active field control scheme and indicate some challenges and limitations for its physical implementation.
Key words: Active field control / Feasibility analysis / Inverse source problem / Integral equation method
© C. Qi et al., Published by EDP Sciences, 2021
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
The active control of acoustic fields has been extensively explored in the past decades and is an emerging research area in modern acoustics. Compared with passive control schemes, active control techniques are more suited for low frequencies, as the barriers used in passive control schemes are relatively large and ineffective at low frequencies [1]. Instead of using the interaction between sound and specific materials to control the acoustic field, active control strategies characterize a source so that it is capable of approximating given field patterns in prescribed exterior regions. The current literature has significantly addressed the idea of the active control of Helmholtz scalar fields in broad applications. These include, but are not limited to, active noise control [2–4], personal sound zones or multizone sound reproduction [5–7], active acoustic cloaking [2, 8–10], remote sensing [11, 12] and metamaterial design [13–15]. Active sound control techniques are becoming increasingly ubiquitous to enhance soundbased systems [16–21]. Forward and inverse problems for sound in underwater environments have been widely studied in the literature (see monographs [22–24] and references therein). Several comprehensive reviews are available in [1, 25, 26], discussing the technicalities and recent advances in various active control schemes.
Directional farfield control can potentially play an important role in some applications with farfield manipulation involved. In this context, in [27, 28], the authors investigated the directional acoustic manipulation via multiphase forehead structure of the porpoises to realize beamforming. In the same paradigm, the works [29–31] made use of loudspeaker arrays to achieve a directional source.
The majority of the active control strategies have been focused on the method of pressure matching (PM) [25]. The PM approach aims to match the target field pattern in the given region with minimum error. In this approach, the active control problem is cast as an inverse source problem (ISP).
Unlike in free space, the active sound control in underwater environments is more challenging. It makes use of additional boundary conditions and the approximated Green’s function, which adds a layer of complexity than free space scenario. We mention here the works [32, 33] where the authors develop a singlemode excitation with a feedback control algorithm to realize both near and farfield sound control.
A general ocean environment can be modeled as a horizontally stratified waveguide [23, 24, 34]. In general, it is fairly difficult to find the analytical fundamental solution and thus, in this paper we follow the paradigm proposed in [22] and consider a simpler marine environment modeled as a shallow water or a homogeneous finitedepth ocean. Our sensitivity analysis builds up on the numerical framework developed in [20, 21, 35]. We use the associated Green’s function to represent the solution of the Helmholtz equation and employ the integral equation (IE) method to formulate the forward propagator. The method of moments (MoM) approach [36] is used to reduce the original integral equation to a discrete linear system. Then a Tikhonov regularization scheme with the Morozov discrepancy principle is applied to solve the resulting system of equations. In the underwater environment, additional boundary conditions, including pressure release boundary on airwater interface and sound hard boundary on the sea floor, were imposed. Consequently, the Green’s function was modified in the formulation of the forward propagator. We used the normal mode representation to formulate the Green’s function in the homogeneous finitedepth ocean [22, 23].
In this paper, we present a detailed sensitivity study for the problem of controlling threedimensional scalar Helmholtz fields in several prescribed exterior regions while maintaining desired farfield pattern values in given fixed directions. Moreover, we also look at the process error, i.e., small noise that may be present in the feeding current when the source is physically implemented. We discuss the feasibility of the active scheme (in terms of power budget, control accuracy and process error) with respect to variations in frequency, the distance between the control region and the source, the mutual distance between the control regions, and the control region size. By a feasible source, we mean a source that requires low power (within the dB levels of common sound sources such as loudspeakers), achieves the control effects with high accuracy (errors of less than 10%) and with high tolerance to feeding noise.
The rest of this paper is organized as follows. In Section 2, we formally describe the problem and provide relevant theoretical results obtained in [35]. Section 3 shows the numerical results and sensitivity analysis in free space. In Section 4, the numerical results and sensitivity analysis in the shallow water environment are presented. Finally, we conclude the paper with some remarks in Section 5.
2 Theory
2.1 Problem formulation
In this section, we present a general description of the active manipulation scheme for Helmholtz fields proposed in our previous works. The unified functional and numerical framework have already been discussed in [20, 21, 35]. We shall briefly recall several essential theoretical results and describe some geometric configurations of interest.
The problem is to characterize a source (modeled as surface pressure or surface normal velocity) so that its generated field approximate some prescribed fields in several exterior regions of interest while maintaining desired patterns in several given far field directions. In this paper, the active field manipulation scheme is explored in both free space and homogeneous ocean with a constant depth. The problem geometry in these two environments are sketched in Figures 1 and 2. Although the theoretical discussion in [20, 21, 35] indicates that an arbitrary number of source regions, exterior control regions, and farfield directions can be considered in the active control scheme, here we shall only consider a single source D_{a}, two control regions D_{1}, D_{2}, and two farfield directions x_{1}, x_{2} for illustrative purposes. A single source D_{a} ⋐ is modeled as a compact region in both free space and homogeneous ocean. The control regions D_{1} and D_{2} are mutually disjoint smooth domains, i.e., D_{1} ∩ D_{2} = ∅. We also assume that the control regions are wellseparated from the source region, i.e., (D_{1} ∪ D_{2}) ∩ D_{a} = ∅. Furthermore, we consider two distinct directions x_{1} and x_{2} representing the farfield directions of interest, which can be arbitrary.
Figure 1 Sketch of the problem geometry showing the near controls D_{1}, D_{2} and the far field directions x_{1} and x_{2} in free space. 
Figure 2 Sketch of the problem geometry showing the near controls D_{1}, D_{2} and the far field directions x_{1} and x_{2} in homogeneous ocean with constant depth. 
Mathematically the problem is to find the boundary input on the source, either a Dirichlet input data p (pressure) or a Neumann input data v_{n} (normal velocity) such that for any desired field f = (f_{1}, f_{2}) on the control regions D_{1}, D_{2} and prescribed far field pattern values f_{∞} = (f_{∞,1}, f_{∞,2}), the solution u of the following exterior Helmholtz problem:
(1)satisfies the control constraints,
where 0 < μ ≪ 1 is the desired control accuracy threshold and u_{∞} denotes the the far field pattern of u. In (1) and (2), n denotes the outwardpointing normal vector to ∂D_{a} and is the unit vector along the direction x. The subscript C^{2} denotes that the norm is defined in the space of smooth functions with continuous derivatives up to the second order and is computed as a sum of the L^{2} norms of all the partial derivatives up to the second order. Recall that the L^{2} norm on the space of square integrable functions on a given domain D is defined as . The boundary and radiation conditions in (1) depend on the medium (environment) and will be given in Sections 2.2 and 2.3.
In [20, 21], it was demonstrated that problem (1) together with (2) admits a solution if the wavenumber k is not a resonance, i.e., −k^{2} is neither a Neumann eigenvalue for the Laplacian in the source region nor a Dirichlet eigenvalue in the control regions. In order to ease the analysis and integral computations, our scheme makes use of a “fictitious source”, i.e., an arbitrary sphere compactly embedded in the actual source region D_{a}. In general, the physical source D_{a} can have any shape as long as it has a Lipschitz boundary, compactly includes the fictitious source and is well separated from the control regions. Meanwhile, our scheme uses slightly larger mutually disjoint regions W_{1} and W_{2} such that D_{1} ⋐ W_{1}, D_{2} ⋐ W_{2}, and because, as shown in [20], an accurate control in the sense of the L^{2}norm on ∂W_{1} and ∂W_{2} implies, via regularity and uniqueness results for the solution of the interior Helmholtz equation, the smooth control required in (2). As pointed out in [21, 35], within the framework mentioned above, the boundary input data, either normal velocity v_{n} or pressure p on the surface of the source can be characterized using a smooth density function such that,
for x ∈ ∂D_{a} and where ρ is the density of the surrounding environment, c is the speed of sound in the given medium and ϕ(x, y) is the fundamental solution of the 3D Helmholtz equation. In general, the solution of problem (1) can be represented by a linear combination of a single and a double layer potential illustrated in [37, 38]. For simplicity of the computations, we only use the singlelayer potential operator throughout this paper.
2.2 Free space environment
In this section, the active manipulation of Helmholtz fields in a freespace environment is investigated. In this medium, problem (1) reads,
The fundamental solution in this case is given by . Then, the control problem is to characterize p or v_{n} such that the solution u of (5) satisfies (2). In free space, the solution u can be written as the image of a density function w under a forward propagator operator . To define , we make use of the near field operator on the control region ∂W_{l} defined as,
where for each , z_{l} ∈ ∂W_{l} and . Moreover, following the derivation in [38], the far field pattern operator can be defined as,
where , is the unit vector pointing in the farfield direction of interest. Hence, the overall propagator operator in free space is defined as,
2.3 Homogeneous ocean environment
Compared with the free space regime, the active control scheme in the homogeneous ocean with a constant depth is much more complicated. As shown in Figure 2, two additional boundary conditions must be specified. More explicitly, problem (1) now reads,
where following the framework proposed in [22], we employ cylindrical coordinates in our analysis and the functions u_{p}’s represent the normal modes in the representation of u. The main control problem is then to characterize v_{n} or p such that the solution u of (9) satisfies (2). We employ the following Green’s representation for u,
where w is the density function defined on the fictitious surface source and G is the associated Green’s function in the medium. For any observation point x = (r, θ, z) = (ξ, z) and source point y = (r′, θ′, z′) = (ξ′, z′), the Green’s function has the following normal mode representation [24]:
where is the zero order Hankel function of the first kind, ϕ_{p} is the pth modal solution with associated eigenvalue a_{p} [22–24]. These eigenvalues are,
while the separated modal solutions ϕ_{p} are given by,
In the farfield region, the field u has an asymptotic form given by [22],
where N is the number of propagating modes (i.e., the largest integer so that ), g_{p} is given by,
and,
where ϵ_{0} = 1, and ϵ_{q} = 2 for q ≥ 1 with,
where J_{p}(x) is the Bessel function of the first kind of order p. Therefore, the farfield pattern in a given direction x = (1, θ, z) can be defined as [22],
Remark 2.1
In (5) and (9), we assume that the free space and homogeneous environment are idealized. In free space, there is no additional boundary condition considered. In the shallow water environment, only airwater surface and the sea floor are applied as two boundary conditions. These boundary conditions are used to demonstrate the sensitivity study. Adaptive or random boundary conditions and the corresponding Green’s function have to be imposed when a practical environment is considered.
Similar to the free space regime, we define a propagator operator that calculates the generated field on the control regions and the farfield pattern in the given directions. For each and , define,
Then the overall propagator operator is given by,
2.4 Optimization scheme
In Sections 2.2 and 2.3, we already defined the propagator operator that evaluates the field in the exterior control regions and the farfield pattern values. The problem (1), (2), formulated in the respective context of free space model (5) or constantdepth homogeneous ocean model (9) can be summarized as,
where the operator is defined in (8) and (12), w is the unknown density function defined on the fictitious source and f denotes the prescribed fields in the control regions. Following the approach in [17] the density function w in (23) is determined using the method of moments (MoM) by discretizing the control regions into a discrete mesh of collocation points and w being expressed as a linear combination (with unknown coefficients) of local basis functions spanning the space of square integrable functions on . Thus, the integral equation in (23) is reduced to a linear system,
where w_{d} represents the vector of unknown coefficients in the local basis representation of w, A represents the matrix of moments computed from the propagator and b is the vector of values of f in the mesh of evaluation points distributed within the control regions together with the prescribed two far field directions. The matrix A is not invertible in most cases, thus the linear system (24) is solved using a regularization routine to minimize the sum of squared residuals. Following the strategy in [16, 20, 21], the unknown coefficients in w_{d} are obtained by using Tikhonov regularization, and can be compactly written as
where l = 1 or 2 denoting the index of the control region. α is the regularization parameter representing the penalty weight for the power required by the solution. The optimal α is determined by the Morozov discrepancy principle [39, 40]. The unknown discrete coefficients in w_{d} are taken to be the Tikhonov solution,
where I is the identity matrix and A^{*} is the complex conjugate transpose of A.
Remark 2.2
In (25), the weighting coefficient η_{l} can be applied to change the contributions of two data misfit in near region and far direction. The weighting coefficient determines the control effort made by the source in each region. For example, if η_{1} → ∞, the data misfit in the near region 1 will be minimized to approach 0 such that the total misfit is less than a desired threshold. Note that the weighting coefficient also plays an important role to the complexity of the boundary input.
To estimate the power on the actual source D_{a}, the average radiated power P_{ave} and the stored energy P_{stor} are computed as,
where u^{*} denotes the complex conjugate and B_{R} is some sphere of radius R containing the actual source D_{a}. In our analyses, the calculated power is expressed in dB relative to a reference level of 10^{−12} W. Finally, the sought boundary input: either normal velocity v_{n} or pressure p on the actual source is obtained from (3) and (4), respectively.
3 Numerical results in free space
In this section, we present several relevant numerical simulations to support the above mentioned theoretical framework. We start from a simplified geometric configuration as shown in Figure 3a with one near control region D_{1} and one farfield direction x_{1} which is exactly behind the near control. Then, we extend our numerical study into a multipleregion regime with two near field control regions D_{1} and D_{2} and two far field directions x_{1} behind D_{1} and x_{2} behind D_{2} as sketched in Figure 3b. The source and control regions are in free space (with medium parameters c = 343 m/s and ρ = 1.225 kg/m^{3}). Throughout this section, the fictitious source region is the sphere of radius 0.2 m centered at the origin. In general, the actual source D_{a} can be arbitrarily shaped as long as it is Lipschitz and compactly embeds . In our simulations, for exemplification, we assume the actual source to be the sphere of radius 0.22 m centered at the origin. In the spirit of [35], the section starts with Sections 3.1 and 3.2 discussing the performance of our strategy, which is quantified by the relative or absolute error in control regions, in each of the above mentioned configurations and then continues with Section 3.3 where we present a detailed sensitivity analysis for free space.
Figure 3 Sketch of the top view of the problem geometry showing the near control(s) and the far field direction(s). (a) One near control and one far field direction. (b) Two near controls and two far field directions. 
3.1 A null in the near control region and a nonzero far field pattern
In this subsection we show the performance of our scheme in creating a null in D_{1} and a given pattern f_{∞,1} in the prescribed far field direction x_{1} (see Fig. 3a for one possible configuration). In our simulation, we set the wavenumber to k = 10 and consider D_{1} ⋐ W_{1} where W_{1} is an annular sector given in the spherical coordinates (with respect to the origin) by,
The farfield direction is exactly behind the near control, i.e., , for large r. The desired field in region D_{1} is f_{1} = 0 while the desired far field pattern in direction x_{1} is given by f_{∞,1} = 0.01 + i · 0.02 and . The simulation results are shown in Figures 4 and 5.
Figure 4 Generated field on the plane z = 0. 
Figure 5 Boundary input v_{n} on the actual source ∂D_{a}. 
The generated field on a 2D cross section through z = 0 around the control region is shown in Figure 4. This shows that the field inside the control region has amplitude not exceeding −120 dB. In the linear scale, the generated field’s maximum pointwise absolute value is less than 1.6 × 10^{−8} while its L^{2} norm is just about 6.69 × 10^{−8}. As a numerical stability check, these values are computed using points slightly off from the mesh points used in the collocation scheme. To describe the accuracy of the scheme, we use the pointwise measure of error e_{i} defined as,
where u = Aw_{d} is the generated field and f_{i} is the prescribed value in the ith evaluation point. In the far field direction x_{1}, the maximum pointwise error of the real and imaginary parts are 4.79 × 10^{−8} and 1.68 × 10^{−8}, respectively suggesting good control accuracy.
Remark 3.1
In Figure 4, once a scatterer besides its reflection appears in the energy radiation direction, two approaches can be applied to treat the reflections. The first method modifies the surface input of the source accounting for nulls around the scatterer, i.e., creating nulls around the scatterer such that the energy would be by design radiated away from the scatterer. Theoretically, it is possible in far field or for small scaterrers in near field. If the scatterer is large or there is a scattering boundary, the second method which incorporates the scattered field in the total propagator has to be used. This requires a different mathematical formulation and is not considered in this paper.
Figure 5 shows the normal velocity v_{n} on the source ∂D_{a} in a rectangular (θ, ϕ)plot. The pointwise amplitudes of v_{n} are quite small and within order 10^{−2}. The average power and stored energy in the actual source are 7.1 dB and 7.8 dB, respectively, which implies that it is feasible for physical implementation, insofar as power budget is concerned.
Remark 3.2
The boundary input data shown in Figure 5 can approximate the field patterns in the exterior regions and far directions. In further numerical simulations [41], our continuous sources are proven to be able to achieve controls at regions that are much (10 wavelengths) further away from the source and that are much (15 times) larger than the source. This suggests that our approach may be useful to push the current understanding about the physical limits.
3.2 Two near control regions and two far field directions
In this subsection we show the performance of our scheme in creating null fields in D_{1} and D_{2} while approximating two distinct prescribed patterns f_{∞,1}, f_{∞,2} in farfield directions x_{1} and x_{2}, respectively (see Fig. 3b for one possible configuration). One possible application of this configuration would be in developing a strategy to establish and maintain communication in several given farfield directions while avoiding near field interference located in D_{1} and D_{2}. In this test, k = 10, D_{1} ⋐ W_{1} with W_{1} defined in (29) while D_{2} ⋐ W_{2} with W_{2} given in spherical coordinates as,
In this geometry, and , r ≫ 1 while the prescribed farfield pattern values are, f_{∞,1} = 0.01 + i·0.02 in direction x_{1} and f_{∞,2} = 0.05 + i · 0.03 in direction x_{2}. The simulation results are shown in Figures 6 and 7. In Figure 6, we see that the absolute errors in the two near controls are within order 10^{−7}. In the far directions x_{1} and x_{2}, the supremum relative errors of the real parts are 3.29 × 10^{−7} and 6.39 × 10^{−8} while for the imaginary parts, the supremum errors are 8.02 × 10^{−8} and 1.61 × 10^{−7}, respectively. In Figure 7 we present the normal velocity required on the actual physical source D_{a}. The average radiated power is 8.4 dB and the stored energy is 9.7 dB. The power budget is slightly larger than that in Section 3.1 as the source is projecting two nonzero farfield patterns in two directions.
Figure 6 Pointwise magnitudes of the generated fields on (a) D_{1} and (b) D_{2}, both approximating a null field. 
Figure 7 Normal velocity v_{n} on the actual source ∂D_{a}. 
3.3 Sensitivity analysis
The aforementioned results support the analysis of [35] and show that our strategy works for each of the two configurations depicted in Figure 3. In the following tests, we aim to study the sensitivity of our strategy with respect to variations in several physically relevant parameters, such as wavenumber k, the distance between the control region and the source, the control region size and, mutual distance between the control regions (in the case of more control regions and far field directions Fig. 3b). The feasibility of the active control scheme is also discussed by looking at the overall control accuracy, power budget and its performance against process error. The geometry in the sensitivity analysis is depicted in Figure 8. To distinguish the original region (dark), the regions with modified parameters are shown in a light color. For instance, denotes the near control which is shifted away from the source, in which the superscript `2’ corresponds to the experiment number in Section 3.3.2.
Figure 8 Sketch of the geometry in sensitivity analysis. D_{1}, D_{2} and x_{1}, x_{2} are original near controls and far field directions, respectively. They are shown in dark color. The lightcolored regions with , where n = 2, 3 and 4, are corresponding to the experiments in Sections 3.3.2 to 3.3.4. 
3.3.1 Varying the wavenumber k
We start with the initial geometry in Figure 3a, i.e., only one near control and a farfield direction. The prescribed field in D_{1} is zero and the farfield pattern is a given nonzero complex number. In this simulation we let the wavenumber vary from 1 to 31. In every single simulation, we keep the geometry fixed and only change the wavenumber. The simulation results are shown in Figure 9. The first plot in Figure 9 shows the supremum pointwise absolute error in D_{1} and the relative error in the far field direction x_{1}. The behavior is less optimal in the low frequency case since the condition number of the associated matrix grows large with decrease in frequency while our numerics suggest the performance of the strategy gets worse for larger wavenumbers (k > 200). The second plot is the process error, which is the supremum pointwise error when the computed feeding current w_{d} is contaminated with small noise. Mathematically, the process error can be obtained by replacing u with u_{δ} in (30), where u_{δ} = Aw_{dδ} and δ = 0.0001 is the noise threshold. In this paper, we consider a random Gaussian noise such that w_{dδ} = w_{d}·(1 + δ·_{2}·) and (0,1). The plot of the process error in Figure 9 indicates that the system is capable of overcoming feedingcurrent noise of a certain level. The power budget and L^{2} norm of v_{n} show that the control effort decreases as the operating frequency increases.
Figure 9 Results showing the control accuracy and the power budget varying with k. From left to right, (1) Supremum error. (2) Process error. (3) Power budget on D_{a}. (4) L^{2} norm of normal velocity v_{n} on D_{a}. 
3.3.2 Varying the distance between the near control region and the source
The following results show the effect of variations in the distance between the near control region and the source. We still use the initial model in Figure 3a. The near control D_{1} is shifted further away from the source ( in Fig. 8a) while all other parameters are fixed. The results are depicted in Figure 10. We notice that the control accuracy in the near region keeps on improving when the near control is moved further away from the source, while the control accuracy in the far direction is converging to a certain value. This indicates that the generated field on D_{1} becomes closer to the null field as it is moved further away from the source. This may be desired for scenarios when D_{1} is an acoustic scatterer. The power budget on the source follows a similar trend as the control accuracy.
Figure 10 Results showing the control accuracy and the power budget varying with mutual distance between D_{1} and . From left to right, (1) Supremum error. (2) Process error. (3) Power budget on D_{a}. L^{2} norm of normal velocity v_{n} on D_{a}. 
3.3.3 Varying the near control region size
Now we consider the behavior of the control accuracy and the power budget with respect to incremental increase in the outer radius of the near control region D_{1} ( in Fig. 8a) with all the other parameters kept fixed. The results are shown in Figure 11. Notice that the control accuracy and power budget are slightly oscillating in the entire range of near control size. The results indicate a good performance for larger obstacles.
Figure 11 Results showing the control accuracy and the power budget varying with the size of the near control D_{1}. The size of D_{1} is the difference between the outer and inner radii, i.e., the thickness of the sectorial region. From left to right, (1) Supremum error. (2) Process error. (3) Power budget on D_{a}. (4) L^{2} norm of normal velocity v_{n} on D_{a}. 
3.3.4 Varying the mutual distance between the near control regions
In this sensitivity test, we consider two near control regions and vary the mutual distance between them. In Figure 8b, D_{1} is fixed and we rotate D_{2} around D_{a} to obtain a new secondary control region . For a fair comparison, the farfield directions are kept exactly behind the controls, i.e., the second far field direction is also rotated (see ). The mutual distance between the two near controls is determined by the mutual angle ϕ shown in Figure 8b. The control accuracy and power budget are recorded as the angle ϕ is varied from 3.6° to 356.4°.The results are shown in Figure 12. We find that the relative error is less than 10^{−5} and the power budget are within 12 dB if ϕ is in the range of [30°, 330°]. If ϕ is out of this range, i.e., the near controls are too close to each other, the relative error and source power become excessive. In the lower and upper bounds of the ϕrange, the relative error in the farfield patterns can be as high as 15%. This indicates that the control regions cannot be too close, otherwise the accurate control effects are not guaranteed. This makes sense as it is hard to generate two very different patterns in two far field direction which are very close due to the continuity of Helmholtz fields.
Figure 12 Results showing the control accuracy and the power budget varying with the mutual distance between near controls D_{1} and D_{2} (given as a function of the angle ϕ between them). From left to right, (1) Supremum error. (2) Process error. (3) Power budget on D_{a}. (4) L^{2} norm of normal velocity v_{n} on D_{a}. 
4 Numerical results in the homogeneous ocean
In this section, we extend our sensitivity analysis to the homogeneous ocean regime. We follow a similar procedure as in the previous section that deals with the problem in free space. We only show part of the numerical results due to the page limitation. The rest of the numerical examples are available in the supplementary document [41]. In the entire section, the fictitious source is the sphere of radius 0.2 m with center at (0, 0, −50) while, for exemplification, the physical source is chosen to be the concentric sphere of radius 0.22 m. The control region D_{1} and the farfield direction x_{1} are given by the control region and farfield direction described in the free space simulations but this time shifted 50 m downward (see Fig. 2). The depth of the ocean environment is h = 100 m, the speed of sound c is assumed to be 1515 m/s and the density ρ is 1020 kg/m^{3}.
4.1 A null in the near control region and a nonzero far field pattern
We start from the initial geometry in Figure 3a but with k = 1 and recall, in the spirit of [35], the performance of our strategy. In D_{1}, we prescribe a null field, and the desired farfield pattern value is f_{∞,1} = 0.01 + i·0.02. The results are shown in Figures 13 and 14. Figure 13 shows the generated fields on two cross sections of the ocean. The left plot (a) shows that the generated field on W_{1} has pointwise amplitude less than −60 dB. In the linear scale, this corresponds to a maximum magnitude of the generated field within order 10^{−5}. The right plot (b) shows the field on the plane y = 0, where the reflections from the ocean’s surface and bottom are observable.
Figure 13 Generated fields on cut planes: (a) z = −50 and (b) y = 0 
Figure 14 Normal velocity v_{n} on the actual source ∂D_{a}. 
In the exact direction x_{1}, the relative errors of the real part and imaginary part are 5.24 × 10^{−4} and 1.83 × 10^{−4}, respectively. In Figure 9, we show the boundary input v_{n} on the actual source ∂D_{a}. The average power and stored energy on the actual source are 11.2 dB and 11.6 dB, respectively. This low power budget suggests the possibility of physically instantiating this source.
4.2 Sensitivity analysis
In this subsection we repeat the same experiments performed in the free space environment to explore the limitations or challenges of the active control in the homogeneous ocean environment. We only present here the effect of frequency and of the mutual distance between the near control and the source on the control accuracy and power budget. Some other experiment results are available online via the clickable https://drive.google.com/file/d/1E1ApCAUwBokeGuqiumu1GKzUIKR3u3N1/view?usp=sharing; Google Drive link [41].
4.2.1 Varying the wavenumber k
In this simulation, the wavenumber varies from 1 to 31 while the other parameters are fixed. The simulation results are shown in Figure 15. In contrast to the free space results shown in Figure 9, the relative error increases as the frequency increases. This is due to the complex model (propagator, boundary condition, and so on) used in the homogeneous ocean environment. In practical underwater acoustic applications, the underwater environment poses serious challenges which are much more complicated than that in the free space. Hence, the control effort required on the source to accurately manipulate the near field and directional far field pattern is more than that in free space. Furthermore, the active scheme in the underwater environment is more sensitive to the feeding noise. For the smaller noise threshold (δ = 10^{−8}) the process error in the far direction is kept smaller than 10%. This implies that even a low noise level in the feeding network may have a major effect on the radiated pattern in a far field direction. However, the system still maintains a good performance in the near control. The overall trend of the power budget is also increasing as the wavenumber increases. It is noted that the spikes on the curves are due to the wavenumber being close to resonance frequencies.
Figure 15 Results showing the control accuracy and the power budget varying with wavenumber k. From left to right, (1) Supremum error. (2) Process error. (3) Power budget on D_{a}. (4) L^{2} norm of normal velocity v_{n} on D_{a}. 
4.2.2 Varying the distance between the near control region and the source
In this simulation, we test the sensitivity of our strategy as D_{1} is moved away from the source as shown in Figure 8. The initial near control region is D_{1} and it is incrementally pushed away from the source to obtain a new control region . The results are shown in Figure 16. We observe that the overall performance of control accuracy and power budget is better when the near control is further away from the source, save some spikes due to resonances. The reason accounting for this trend is that the near control, where we want a field with low amplitude (as it may be interpreted as an obstacle or scatterer for some applications), is pushed away and hence the source effort is eased.
Figure 16 Results showing the control accuracy and the power budget varying with the distance between the source and the control region. From left to right, (1) Supremum error. (2) Process error. (3) Power budget on D_{a}. (4) L^{2} norm of normal velocity v_{n} on D_{a}. 
5 Conclusions
In this paper, the feasibility of the active manipulation of Helmholtz fields both in free space and in a homogeneous ocean of constant depth is presented. We build up on our previous works and demonstrated the possibility of of having a feasible characterization of an active source (modeled as surface pressure or surface normal velocity) such that it is capable of approximating a priori given field some near control regions while simultaneously projecting desired patterns in several farfield directions. By a feasible source, we mean a source with low power requirement, for instance 30 dB, which can radiate the prescribed fields within some accuracy thresholds. We showed the control accuracy and the power budget of the proposed active control mechanism for each environment. Then we explored the behavior of physically relevant parameters (power budget and control accuracy) with respect to variations in the frequency, outward shift, the outer radius of the near control and the mutual distance between near controls.
In our simulations, we considered the two initial models shown in Figure 3. The first one contains one near control and one farfield direction and the second one has two near control regions and two farfield directions. The farfield directions are placed exactly behind the near control regions in each of the two models. In this paper, we only show the cases in which the near control regions are prescribed to be have a null field and the given farfield directions to have a nonzero pattern.
In free space, the control accuracy is within order 10^{−8} both in the near control and in the far field direction. In the geometry shown in Figure 3a, the operating frequency is first varied while all other parameters are fixed. The frequency is swept from 54.59 Hz to 1.69 kHz (k is from 1 to 31). Moreover, the power requirement is kept at low levels which suggests the feasibility of a physical implementation of the calculated source. In the second simulation, the near control region is moved outward. We noticed that the control accuracy in the near region is continuously improving as the near control region is pushed further away from the active source. However, accuracy converges to a certain level in the farfield direction. Next, we vary the outer radius of the near control region to explore the effect of the near region’s size on the control accuracy and power budget. The simulation results show that the active control scheme works in the entire range of values of the outer radii. We also considered the geometry in Figure 3b. Here, we rotate the second near control region together with the farfield direction behind it. We find that the control accuracy and power budget don’t change significantly if the two near control regions are separated by an angle ϕ ∈ [30°, 300°]. However, outside this range the control performance is gradually degrading.
We then extended our sensitivity analysis into the case of a homogeneous ocean of constant depth. The Green’s function, corresponding to a pressurerelease surface and a totally reflecting bottom, is expressed using the normal mode representation. Accordingly, the farfield pattern propagator is defined. In this case, the control accuracy is within orders 10^{−5} and 10^{−4} in the near control and in the far direction, respectively. In the first sensitivity test we varied the wavenumber from 1 to 31. The results show that the overall performance of the scheme, with respect to the control accuracy and the power budget decreases with increase in frequency. These suggest that the control scheme is more suitable to low frequencies in the homogeneous ocean environment. In the second test, we move the near control further away from the source. The results are similar to that in the free space. More sensitivity test results are available in the supplementary material [41].
Conflict of interest
Author declared no conflict of interests.
Acknowledgments
This research has been supported in part by the Army Research Office under award number W911NF1710478 and in part by the National Science Foundation under award 1801925.
References
 S.M. Kuo, D.R. Morgan: Active noise control: A tutorial review. Proceedings of the IEEE 87, 6 (1999) 943–973. [Google Scholar]
 D. Eggler, H. Chung, F. Montiel, J. Pan, N. Kessissoglou: Active noise cloaking of 2d cylindrical shells. Wave Motion 87 (2019) 106–122. [Google Scholar]
 J. Pan, R. Paurobally, X. Qiu: Active noise control in workplaces. Acoustics Australia 44, 1 (2016) 45–50. [Google Scholar]
 S. Wang, H. Sun, J. Pan, X. Qiu: Nearfield error sensing for active directivity control of radiated sound. The Journal of the Acoustical Society of America 144, 2 (2018) 598–607. [CrossRef] [PubMed] [Google Scholar]
 T. Betlehem, P.D. Teal: A constrained optimization approach for multizone surround sound, in 2011 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), IEEE, 2011, pp. 437–440. [Google Scholar]
 L.J. Brännmark, A. Bahne, A. Ahlen: Compensation of loudspeaker–room responses in a robust mimo control framework. IEEE Transactions on Audio, Speech, and Language Processing 21, 6 (2013) 1201–1216. [Google Scholar]
 J. Cheer, S.J. Elliott, M.F. Simón Gálvez: Design and implementation of a car cabin personal audio system. Journal of the Audio Engineering Society 61, 6 (2013) 412–424. [Google Scholar]
 D. Eggler, M. Karimi, N. Kessissoglou: Active acoustic cloaking in a convected flow field. The Journal of the Acoustical Society of America 146, 1 (2019) 586–594. [CrossRef] [PubMed] [Google Scholar]
 F.G. Vasquez, G.W. Milton, D. Onofrei: Active exterior cloaking for the 2d laplace and helmholtz equations. Physical Review Letters 103, 7 (2009) 073901. [CrossRef] [PubMed] [Google Scholar]
 F.G. Vasquez, G.W. Milton, D. Onofrei: Exterior cloaking with active sources in two dimensional acoustics. Wave Motion 48 (2011) 515–524. [Google Scholar]
 S.J. Elliott, J. Cheer: Modeling local active sound control with remote sensors in spatially random pressure fields. The Journal of the Acoustical Society of America 137, 4 (2015) 1936–1946. [CrossRef] [PubMed] [Google Scholar]
 W. Jung, S.J. Elliott, J. Cheer: Combining the remote microphone technique with headtracking for local active sound control. The Journal of the Acoustical Society of America 142, 1 (2017) 298–307. [CrossRef] [PubMed] [Google Scholar]
 H. Chen, C.T. Chan: Acoustic cloaking in three dimensions using acoustic metamaterials. Applied Physics Letters 91, 18 (2007), 183518. [Google Scholar]
 S.A. Cummer, J. Christensen, A. Alù: Controlling sound with acoustic metamaterials. Nature Reviews Materials 1, 3 (2016) 16001. [Google Scholar]
 R.V. Kohn, J. Lu, B. Schweizer, M.I. Weinstein: A variational perspective on cloaking by anomalous localized resonance. 2012. arXiv:1210.4823 [math.AP]. [Google Scholar]
 N.J.A. Egarguin, D. Onofrei, E. Platt: Sensitivity analysis for the active manipulation of helmholtz fields in 3d. Inverse Problems in Science and Engineering 28, 3 (2020) 314–339. [Google Scholar]
 N.J.A. Egarguin, S. Zeng, D. Onofrei, J. Chen: Active control of helmholtz fields in 3d using an array of sources. Wave Motion 94, 102523 (2020) 1–27. [Google Scholar]
 C. House, J. Cheer, S. Daley: On the use of virtual sensing for the realtime detection and active control of a scattered acoustic field, in International Congress of Sound and Vibration. 2019. [Google Scholar]
 P.A. Nelson, S.J. Elliott: Active Control of Sound. Academic London, 1992. [Google Scholar]
 D. Onofrei: Active manipulation of fields modeled by the helmholtz equation. Journal Of Integral Equations and Applications 26, 4 (2014) 553–579. [Google Scholar]
 D. Onofrei, E. Platt: On the synthesis of acoustic sources with controllable near fields. Wave Motion 77 (2018) 12–27. [Google Scholar]
 J.L. Buchanan, R.P. Gilbert, A. Wirgin, Y.S. Xu: Marine Acoustics: Direct and Inverse Problems. SIAM, 2004. [Google Scholar]
 F.B. Jensen, W.A. Kuperman, M.B. Porter, H. Schmidt: Computational Ocean Acoustics. Springer, 2011. [Google Scholar]
 J.B. Keller, J.S. Papadakis: Wave propagation and underwater acoustics, in Number 70 in Lecture Notes in Physics, SpringerVerlag. 1977. [Google Scholar]
 T. Betlehem, W. Zhang, M.A. Poletti, T.D. Abhayapala: Personal sound zones: Delivering interfacefree audio to multiple listeners. IEEE Signal Processing Magazine 32, 2 (2015) 81–91. [Google Scholar]
 J. Jiang, Y. Li: Review of active noise control techniques with emphasis on sound quality enhancement. Applied Acoustics 136 (2018) 139–148. [Google Scholar]
 E. Dong, Y. Zhang, Z. Song, T. Zhang, C. Cai, N.X. Fang: Physical modeling and validation of porpoises’ directional emission via hybrid metamaterials. National Science Review 6, 5 (2019) 921–928. [Google Scholar]
 Y. Zhang, Z. Song, X. Wang, W. Cao, W.W.L. Au: Directional acoustic wave manipulation by a porpoise via multiphase forehead structure. Physical Review Applied 8, 6 (2017) 064002. [Google Scholar]
 E. Fisher, B. Rafaely: Nearfield spherical microphone array processing with radial filtering. IEEE Transactions on Audio, Speech, and Language Processing 19, 2 (2010) 256–265. [Google Scholar]
 M. Poletti, T. Betlehem: Design of a prototype variable directivity loudspeaker for improved surround sound reproduction in rooms, in Audio Engineering Society Conference: 52nd International Conference: Sound Field ControlEngineering and Perception, Audio Engineering Society. 2013. [Google Scholar]
 B. Rafaely, D. Khaykin: Optimal modelbased beamforming and independent steering for spherical loudspeaker arrays. IEEE Transactions on Audio, Speech, and Language Processing 19, 7 (2011) 2234–2238. [Google Scholar]
 J.R. Buck, J.C. Preisig, M. Johnson, J. Catipovic: Singlemode excitation in the shallowwater acoustic channel using feedback control. IEEE Journal of Oceanic Engineering 22, 2 (1997) 281–291. [Google Scholar]
 D. Peng, T. Gao, J. Zeng: Study on singlemode excitation in timevariant shallow water environment. Journal of Computational Acoustics 22, 01 (2014) 1440001. [Google Scholar]
 N.C. Makris, P. Ratilal: A unified model for reverberation and submerged object scattering in a stratified ocean waveguide. The Journal of the Acoustical Society of America 109, 3 (2001) 909–941. [CrossRef] [PubMed] [Google Scholar]
 N.J.A. Egarguin, D. Onofrei, C. Qi, J. Chen: Active manipulation of helmholtz scalar fields: Near field synthesis with directional far field control. Inverse Problems 36, 9 (2020) 095005. [Google Scholar]
 R.F. Harrington: Field computation by moment methods. WileyIEEE Press, 1993. [Google Scholar]
 J. Ahrens: The singlelayer potential approach applied to sound field synthesis including cases of nonenclosing distributions of secondary sources. 2010. [Google Scholar]
 D. Colton, R. Kress: Inverse Acoustic and Electromagnetic Scattering Theory. 3rd ed., SpringerVerlag, 2013, pp. 83–94. [Google Scholar]
 T. Bonesky: Morozov’s discrepancy principle and tikhonovtype functionals. Inverse Problems 25, 1 (2008) 015015. [Google Scholar]
 D. Colton, R. Kress: Inverse Acoustic and Electromagnetic Scattering Theory. SpringerVerlag, 3rd ed, 2013. [Google Scholar]
 C. Qi, N.J.A. Egarguin, D. Onofrei, J. Chen: Supplementary material. https://drive.google.com/file/d/1E1ApCAUwBokeGuqiumu1GKzUIKR3u3N1/view?usp=sharing, 2021. [Google Scholar]
Cite this article as: Qi C. Egarguin N.J.A. Onofrei D. & Chen J. 2021. Feasibility analysis for active near/far field acoustic pattern synthesis in free space and shallow water environments. Acta Acustica, 5, 39.
All Figures
Figure 1 Sketch of the problem geometry showing the near controls D_{1}, D_{2} and the far field directions x_{1} and x_{2} in free space. 

In the text 
Figure 2 Sketch of the problem geometry showing the near controls D_{1}, D_{2} and the far field directions x_{1} and x_{2} in homogeneous ocean with constant depth. 

In the text 
Figure 3 Sketch of the top view of the problem geometry showing the near control(s) and the far field direction(s). (a) One near control and one far field direction. (b) Two near controls and two far field directions. 

In the text 
Figure 4 Generated field on the plane z = 0. 

In the text 
Figure 5 Boundary input v_{n} on the actual source ∂D_{a}. 

In the text 
Figure 6 Pointwise magnitudes of the generated fields on (a) D_{1} and (b) D_{2}, both approximating a null field. 

In the text 
Figure 7 Normal velocity v_{n} on the actual source ∂D_{a}. 

In the text 
Figure 8 Sketch of the geometry in sensitivity analysis. D_{1}, D_{2} and x_{1}, x_{2} are original near controls and far field directions, respectively. They are shown in dark color. The lightcolored regions with , where n = 2, 3 and 4, are corresponding to the experiments in Sections 3.3.2 to 3.3.4. 

In the text 
Figure 9 Results showing the control accuracy and the power budget varying with k. From left to right, (1) Supremum error. (2) Process error. (3) Power budget on D_{a}. (4) L^{2} norm of normal velocity v_{n} on D_{a}. 

In the text 
Figure 10 Results showing the control accuracy and the power budget varying with mutual distance between D_{1} and . From left to right, (1) Supremum error. (2) Process error. (3) Power budget on D_{a}. L^{2} norm of normal velocity v_{n} on D_{a}. 

In the text 
Figure 11 Results showing the control accuracy and the power budget varying with the size of the near control D_{1}. The size of D_{1} is the difference between the outer and inner radii, i.e., the thickness of the sectorial region. From left to right, (1) Supremum error. (2) Process error. (3) Power budget on D_{a}. (4) L^{2} norm of normal velocity v_{n} on D_{a}. 

In the text 
Figure 12 Results showing the control accuracy and the power budget varying with the mutual distance between near controls D_{1} and D_{2} (given as a function of the angle ϕ between them). From left to right, (1) Supremum error. (2) Process error. (3) Power budget on D_{a}. (4) L^{2} norm of normal velocity v_{n} on D_{a}. 

In the text 
Figure 13 Generated fields on cut planes: (a) z = −50 and (b) y = 0 

In the text 
Figure 14 Normal velocity v_{n} on the actual source ∂D_{a}. 

In the text 
Figure 15 Results showing the control accuracy and the power budget varying with wavenumber k. From left to right, (1) Supremum error. (2) Process error. (3) Power budget on D_{a}. (4) L^{2} norm of normal velocity v_{n} on D_{a}. 

In the text 
Figure 16 Results showing the control accuracy and the power budget varying with the distance between the source and the control region. From left to right, (1) Supremum error. (2) Process error. (3) Power budget on D_{a}. (4) L^{2} norm of normal velocity v_{n} on D_{a}. 

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.