Sensitivity Analysis of Underwater Structural-Acoustic Problems Based on Coupled Finite Element Method/Fast Multipole Boundary Element Method with Non-Uniform Rational B-Splines

: For the direct differentiation technique-based modeling of acoustic fluid–structure interaction and structural-acoustic sensitivity analysis, a coupling algorithm based on the finite element method (FEM) and the fast multipole boundary element method (FMBEM) is suggested. By bypassing the challenging volume parameterization process in isogeometric finite element techniques and the laborious meshing process in traditional FEM/BEM, non-uniform rational B-splines (NURBS) isogeometric analysis (IGA) is utilized to immediately perform numerical analysis on CAD models. The matrix-vector products in the boundary element analysis are accelerated using the fast multipole method (FMM). To hasten the solution of the linear system of equations, the iterative solver GMRES is used. The numerical prediction of the effects of arbitrarily shaped vibrating structures on the sound field is made feasible by the FEM/FMBEM technique. A number of numerical examples are provided to show the applicability and effectiveness of the suggested approach.


Introduction
Underwater acoustics frequently deal with the analysis of acoustic radiation or scattering from elastic objects in heavy fluid.It is only possible to find analytical solutions for acoustic fluid-structure interaction [1] issues when the structure has simple geometry and simple boundary conditions.More real-world problems with complicated geometries cannot be solved analytically; thus, efficient numerical approaches must be devised.
The dynamic behavior of problems involving fluid-structure interactions, acoustics, fracture mechanics, and electromagnetics has been widely studied using FEM.However, there are a number of problems with the FEM when modeling infinite domains.As is commonly known, the BEM has been effectively used to address acoustic concerns as it offers superior precision and straightforward mesh formation.Particularly for exterior acoustic problems, the Sommerfeld radiation condition at infinity is immediately satisfied [2].The boundary integral problem has been quantitatively solved using the Galerkin approach for BEM implementation [3,4].However, the engineering community has always favored the collocation technique.Therefore, the examination of problems with fluid-structure interaction is acceptable for the coupled FEM/BEM approach [5,6].However, due to the CBEM's generation of a dense and non-symmetric coefficient matrix that requires O(N 3 ) arithmetic operations to directly resolve the equation system, for example, when employing the Gauss elimination approach, coupling analysis of underwater structural-acoustic problems continues to be the bottleneck of high computational cost.The fast multipole method (FMM), the fast direct solver, and the adaptive cross approximation approach are just a few of the methods that have been utilized to expedite the resolution of the integral issue.The rapid direct solver was developed by Martinsson and Rokhlin [7,8].It quickly creates a reduced factorization of the matrix's inverse and performs well for problems requiring moderately ill-conditioned matrices.For problems needing a lot of iterations, the Bebendorf and Rjasanow-invented adaptive cross approximation approach [9] creates blockwise low-rank approximants from the BEM matrices.The CBEM system of equations may now be solved more quickly thanks to the development of FMM [10][11][12][13][14]. Therefore, utilizing a coupling strategy based on FEM/fast multipole boundary element method (FEM/FMBEM) [15], it may be possible to handle large-scale fluid-structure interaction problems.This study also recommends the coupling method FEM/FMBEM to handle the challenging problems involving fluid-structure interactions.
Design professionals are increasingly taking into account passive noise management by changing the geometry of constructions.This structural-acoustic optimization has great potential for reducing radiated noise [16,17], especially for thin-shell constructions.Because it may demonstrate how a geometry adjustment impacts the structure's acoustic performance, acoustic design sensitivity analysis is an essential stage in the procedures of acoustic design and optimization.The advancement of passive noise elimination through structural-acoustic optimization is summarized by Marburg [18].The finite difference method (FDM) has been utilized extensively for structural-acoustic optimization [19][20][21] because of how easy it is to apply.Yet, this approach has poor performance, especially when numerous design factors are given consideration at once.To solve this problem, we use the adjoint variable method (AVM) [22,23] or the direct differentiation method (DDM) [24].The sensitivity analysis for the fluid-structure interaction problem, as is well known, takes up the greatest time in the gradient-based optimization procedure.In order to speed up the analysis, the coupling method FEM/FMBEM is subjected to the structural-acoustic sensitivity analysis based on DDM in this work.
The FEM and BEM may be applied with the help of suitable software; this technique is called computer-aided engineering (CAE).Modern CAE requires, however, that the models produced by CAD software be converted into simulation-ready models as part of the preprocessing step.Geometric model data transfer in the CAE results in geometry errors.A proposed solution to this problem is the integration of boundary element method (IGABEM) [25,26] with geometric modeling and numerical simulation utilizing isogeometric analysis [27][28][29].Geometric errors and laborious preprocessing processes may be eliminated, and numerical simulation can be performed directly from the accurate models, all thanks to IGABEM.Elastic mechanics [30], potential difficulties [31,32], wave propagation [33,34], fracture mechanics [35], electromagnetics [36][37][38], and structural optimization [39][40][41][42][43][44] are just a few of the problems that IGABEM has been utilized to solve since its founding.The NURBS IGABEM is used in this investigation.
This work develops the sensitivity formulation for the coupled FEM/BEM analysis and adds FMBEM to the structural-acoustic coupling sensitivity analysis.Coupling FEM/FMBEM is recommended for issues involving fluid-structure interaction and structural acoustic sensitivity assessments.NURBS is applied for the first time to structuralacoustic coupling calculations and acoustic sensitivity analysis, which eliminates geometric errors and improves the calculation accuracy.The precision and effectiveness of this method are demonstrated using numerical examples.

Non-Uniform Rational B-Splines (NURBS)
In order to be thorough, this part addresses some of the fundamental NURBS ideas that are the basis of the isogeometric analysis.The readers are directed to [27,45] for further information.The knot vector made up by a collection of non-decreasing real numbers represented as in Equation ( 1) is a key idea in NURBS.
where n is the overall amount of basis functions, p is the polynomial order, and a is the knot index.One may think of a knot vector as a parametric space in one dimension, as shown in Figure 1.The Cox-de Boor recursion formula is used to express the B-spline basis functions for a given knot vector, as presented in Equations ( 2) and (3).When p = 0, we have and when p ≥ 1, Numerous advantageous characteristics of B-spline basis functions (as demonstrated in Figure 2, different colors representing different orders of NURBS basis functions), including linear independence and local support, make them suitable for numerical analysis.As demonstrated in Equation ( 4), the B-spline curve may be created by linearly combining B-spline basis functions and control points.x where the control point coordinates are represented by the coefficients P a,b .Because of this, a B-spline curve is fundamentally a conversion from a parametric one-dimensional space to actual space.There are the following knot vectors in every dimension of the following two-dimensional parametric spaces: The tensor product property may be used to build the B-spline surface, as illustrated in Equation (7).
where n and m are the corresponding numbers of basis functions in each dimension.It should be mentioned that B-spline control points are not usually on the surface (as demonstrated in Figure 3) due to the absence of the Kronecker delta property (as demonstrated in Figure 2).By linking a weight coefficient with every control point, NURBS extend B-splines.Using NURBS, designers may precisely represent various curves having conic segments, like ellipses and circles, and gain more control over the depicted curve without boosting the degree or the quantity of control points.The NURBS basis functions are constructed from the B-spline basis functions, which are represented in two dimensions as Equation (8).
where w is the weight coefficient.Similar to how B-spline surfaces are formed, NURBS basis functions and control points are used to define NURBS surfaces, as indicated in Equation (9).
Equation ( 9) may be recast as Equation (10) by dropping the notation p in the following and iterating between basis functions or control points using the global index A.
More control points can be added using the knot insertion operator without altering the structural form.By adopting the h-refinement method, this characteristic helps to increase the precision of estimating physical fields while maintaining geometric correctness.

FEM Analysis
Fritze et al. [6,46] detailed the whole structural-acoustic analysis technique, and associated expressions are provided in this section.If a harmonic load is given to the structure, it is possible to infer from the frequency response analysis what the structure's steady-state reaction will be.The structural-acoustic equation's linear system is derived in Equation (11).
where K is the stiffness matrix, i = √ −1 the imaginary unit, ω the harmonic load's excitation frequency, C the damping matrix, M the mass matrix, u the nodal displacement vector, and f the full excitation.
Note that, although the steady-state response maintains the same frequency with the applied load, damping may cause it to have a distinct phase angle.By breaking down the time-dependent impulses into the frequency domain, Equation (11) can still be applied in the event that the load being imposed is not harmonic.A coupling matrix is added to shift the structural nodal load from the fluid effect to fluid nodal pressure to consider how the acoustic pressure exerted on structural surfaces affects aspects.The full excitation, which combines the acoustic load and the structural load, might then be expressed using Equation (13).
where C sf is the coupling matrix, p the fluid nodal pressure, C sf p the acoustic load denoting the effect of the acoustic pressure applied to the structural surfaces, f s the structural load, N s the interpolation function for structural domain, n the external normal direction of the structural surface, N f the interpolation function for the fluid domain, and Γ the interaction surface situated between the areas of fluidity and structure.The coupling matrix C sf directs the structural nodal load from the fluid effect to the fluid nodal pressure.Equation ( 15) may then be used to calculate the nodal displacement.
The fluid loads applied on underwater bodies are very important.Rehman et al. [47] looked at a fast and less expensive semi-empirical method to determine the hydrodynamic coefficients for a convoluted transportation system with two UUVs.As for the details of the acoustic pressure applied to the structural surfaces mentioned above, please refer to Rehman et al. [47].

BEM Analysis
The time-harmonic wave field of sound in the Helmholtz equation is described by Equation ( 16), while the boundary conditions are expressed by Equation (17).
where p is the sound pressure, k the wave number, () the known function given on the border, Γ p the Dirichlet boundary condition, ρ the structural density, ω the frequency of the incoming force, v the normal velocity, Γ q the Neumann boundary condition, z the acoustic impedance, and Γ z the Robin boundary condition.
From Equations ( 16) to (18), one may construct a boundary integral equation that is specific to the structural boundary Γ.
where x is the source point, y the field point, c(x) = 1/2 if the boundary Γ is smooth in the vicinity of the source point x, p the intensity of the incoming wave, G(x, y) the Green's function, q the normal derivative of p, and F the normal derivative of G.
For acoustic issues in three-dimensional problems, Equation ( 19) offers the expression of Green's function G(x, y).
In cases when the source point x has a smooth border Γ, the integral representation's derivative in Equation ( 18) with reference to the outer normal at point x may be expressed as Equation (21).
It is well known that nonuniqueness makes it difficult to apply a single Helmholtz boundary integral equation to problems requiring external boundary values.In this work, the method known as the Burton-Miller strategy [48] combining the linear Equations ( 18) and (21), is employed to successfully address the issue of nonuniqueness.The Hadamard finite part integral method and the Cauchy principal value may also be used to directly and efficiently compute the singular boundary integrals [24] brought about by Equations ( 18) and (21).
The system of linear algebraic equations shown in Equation ( 22) can be obtained [49] if the border Γ is divided into elements by putting all of the collocation point equations for each element's center together and presenting them using matrix representations.
where p i is the nodal pressure brought on by the incoming wave.

FEM-BEM Coupling Analysis
The precise FEM/BEM modeling formulas were published by Fritze et al. [6], and associated expressions are provided in this section.The governing equations, as described in the preceding section, are connected by the continuity constraint across the interaction surface, as indicated in Equation (23).Then, in accordance with Equation (24), it is reasonable to represent the normal velocity v as a function of the displacement u.
As a result of the sluggish convergence of the direct iterations on Equation ( 28), the direct system equation solution would require significantly greater processing power and storage capacity.It is also difficult to obtain highly precise numerical results.We deliver the subsequent method to solve the above non-symmetric linear system without utilizing an iterative solution.You may obtain the coupled boundary element equation [6] in Equation ( 29) by substituting Equation (15) into Equation (27).
A sparse direct solver might be used to obtain the linear equations' solution in Equation (29).The generalized minimum residual (GMRES) iterative solver and FMM are employed to hasten the answer.

Sensitivity Analysis for Shape Design
Finding the optimal design parameters within preset restrictions that define the desired form of the provided structure is the aim of shape optimization.Marjan and Huang [50] demonstrated a topology optimization technique that finds the best load route on the structure to provide a revolutionary jacket foundation design.Huang et al. [51] gave an introduction of using machine learning techniques to improve the sustainability of ships.The study addresses the foundations of machine learning as well as its applications to ship design, operational performance, and trip planning.
Utilizing shape design sensitivity analysis, one may determine the gradients of specified cost functions.The direction in which to look for the ideal ranges of the design variables may then be determined using the gradients that were acquired.Thus, in the process of designing and optimizing acoustic forms, the acoustic shape sensitivity research [24,52] is often the initial and most crucial step.The direct approach, which first determines the sensitivity of the variables, is employed to calculate the function's sensitivity using the chain rule of differentiation.This approach is highly effective since it is so closely related to the analytical process.
By differentiating Equation (11) with respect to the design variable in the sensitivity evaluation for shape design using FEM, Equation (31) may be produced.
Equations ( 18) and ( 21) are differentiated with regard to any chosen design variable to arrive at Equations ( 32) and (33) when the border Γ is smooth surrounding the source point x.
We have Equation (34) The Cauchy principal value and the Hadamard finite part integral approach [24] may be used to directly and efficiently compute the singular boundary integrals introduced by Equations ( 32) and (33).
Equation ( 35) may be created for the sensitivity analysis for shape design utilizing coupling FEM-BEM by applying Equation ( 27) and differentiating Equation ( 29) in relation to the design variable.
It takes a lot of computation effort to explicitly solve Equation ( 35) using the standard BEM since the matrices are complete and asymmetric.However, the FMM and GMRES can be used to accelerate the computing process.In order to speed up the matrix-vector combinations in Equations ( 29) and (35), the FMM is utilized, while the corresponding sensitivity equation and the formula for the FEM-BEM coupling could be solved using GMRES.

Fast Multipole Boundary Element Method (FMBEM)
This section introduces the FMM to quicken the matrix-vector product in Equations ( 29) and (35).The iterative solution GMRES is used to solve the coupling BEM problem and related sensitivity equation.

FMM Formulations for Acoustic State Analysis
The Green's function (Equation ( 19)) is enlarged into Equation (39) for 3D issues.
where y c is one expansion point near point y, () the complex conjugates, j n the first-kind, n-th order spherical Bessel function, h n the n-th order spherical Hankel function of the first type, Y m n the spherical harmonics, defined as Equation (42).
where P m n denotes the related Legendre functions.r, θ, and ϕ are spherical coordinates of vector ⃗ a, such as − → y c x or − → y c y, for example.
The symbol Γ 0 represents a section of the border Γ, which is distant from the source point x.The integrals in Equation ( 18) can be rewritten as Equation (44).
where M m n is the multipole moment, and y c is close to Γ 0 .Introducing the M2M, M2L, and L2L translation operations, we can obtain Equation (48).
where L m n is the low-frequency FMM's local expansion coefficient, and x 1 l an expansion point close to point x.Zheng et al. [53] provides extensive details on the M2M, M2L, and L2L translation processes as well as the formulation of L m n .

Numerical Examples
The efficiency of the suggested approach is demonstrated by the numerical examples for real-world engineering issues in this section.Fortran 95 with OpenMP parallelization is used to build the technique for finding solutions for the numerical analysis.The component linked to incident wave is maintained in boundary integral equations since the examples are exterior acoustic scattering issues.

Spherical Models
Two different models of underwater, thin-pulsed shell are used in this subsection.The first is a simulation of a spherical outer shell with a 1.0 m radius, while the second is a similar shell with the NURBS control point (−1, 0, 0) shifted to (−1.5, 0, 0).The NURBS surface descriptions and accompanying control grids of these two models are shown in Figure 4.The fluid is water with density ρ f = 1.0 × 10 3 kg/m 3 , and the sound wave's velocity in water is c = 1482 m/s.The presented technique's acoustic evaluation and shape sensitivity assessment are verified using the underwater, thin-pulsing, spherical surface example, as displayed in Figure 4a.While taking into consideration an incoming sound wave with an amplitude of 1.0 in the positive x direction, the position (0, 0, 0) is the center of the spherical shell.
Figure 5 gives the relative error of sound pressure between the numerical and analytical solutions, at point (2, 0, 0) with a frequency at 200 Hz, for the spherical shell model shown in Figure 4a.As the total number of elements rises, the computational results converge quickly, as seen in Figure 5. Based on that, the total number of elements in the discretized thin-shell model is set to be 6144.
The outcomes are shown at location (2, 0, 0) in Figure 6. Figure 7a displays the numerical and analytical results for the sound pressure at position (2, 0, 0), respectively, presented in terms of frequencies.The FMM technique is used by GMRES implementation to speed up the linear solution.The conventional BEM's exceptional accuracy is maintained by the FMM technique, as evidenced by the substantial concordance between the analytical and numerical outcomes in this figure.Now, we go on to the investigation of shape sensitivity, which is crucial to shape optimization.In this case, at location (2, 0, 0), the sound pressure serves as the objective function, and the x-coordinate at the control point (−1, 0, 0) is the design variable.The sound pressure sensitivity is depicted in Figure 6b.Figures 6a,b illustrates how the sensitivity and pressure of the sound grows significantly at resonance peaks.The frequency range [100 Hz, 500 Hz] is crucial for this spherical shell model because there the sound pressure is noticeably greater and more variable to changes in the location of the previously mentioned NURBS point.In light of this finding, the frequency range used in the study of the spheroid-like shell problem is taken as [100 Hz, 500 Hz].A spherical shell-like model may be obtained by moving the NURBS control point from its original location (−1, 0, 0) to a new one (−1.5, 0, 0), as illustrated in Figure 4b.The sound pressure at location (2, 0, 0) for the spherical (as in Figure 4a) and spherical-like (NURBS control point (−1, 0, 0) moved to (−1.5, 0, 0), as in Figure 4b) thin-shell models are compared in Figure 7a.In Figure 7b, the sound pressure's sensitivity at (2, 0, 0) to the x-location of the NURBS control point (−1.5, 0, 0) for the Figure 4b model is contrasted with that of the Figure 4a model.Basically, the results for these two models show a similar pattern.It should be observed that in Figure 7b, the extreme of sensitivity (the red line) for the model of Figure 4b occurs at the frequency of about 180 Hz, where the peak of the sound pressure resonance is not as high, as in Figure 7a.According to this phenomena, the highest resonance peak of sound pressure is not always where the bigger value of frequency sensitivity occurs.As a result, each sound pressure resonance peak should be taken into account when calculating sensitivity, necessitating a computation in terms of a frequency range.
Figure 8 shows the sound pressure on the spherical shell's boundary surface at the frequencies of 100 Hz, 300 Hz, and 500 Hz.When the NURBS control point is moved from (−1, 0, 0) to (−1.5, 0, 0), Figure 9 gives the sound pressure on the shell's boundary surface at the frequencies of 100 Hz, 300 Hz, and 500 Hz.Four ways of views are given for each calculation frequency.These figures show symmetrical features for both the x − y and x − z planes.Given that the plane wave occurs along with the x axis, these results make sense.

Submarine Model
This section uses two separate underwater, simplified, thin-walled submarine models, where the amplitude of the plane wave is 1.0 and it propagates down the x axis in a positive direction, just as in the above numerical model.The left end control point of the first model, which is axially oriented along the x axis from left to right, as shown in Figure 10a, is situated at (0, 0, 0).The model consists of three sections and has a total length of approximately 1.35 m.A conical surface is on the left, a cylindrical surface is in the middle, and a hemispherical surface is on the right.By extending the radius of a few of the cylinder surface's NURBS control points in Figure 10a, a partly expanded simple submarine model may be produced, as seen in Figure 10b, where a circle of control points at x = 0.5 m is inflated to twice the radius.Figure 11a, Figure 11c, and Figure 11e, respectively, give the sound pressure comparison at location (5, 0, 0), (20, 0, 0) and (100, 0, 0), for the models shown in Figure 10.The sound pressure's sensitivity in relation to the cylinder's radius at x = 0.5 m for these two models are given in Figure 11b (at point (5, 0, 0), the sound pressure is the objective function), Figure 11d (at point (20, 0, 0), the sound pressure is the objective function), and Figure 11f (at point (100, 0, 0), the sound pressure is the objective function).In these figures, sound pressure and sensitivity show a similar pattern of increasing and then decreasing.Additionally, the value of sound pressure and sensitivity decrease as the separation between the model and the computation point increases (calculation point from (5, 0, 0) to (100, 0, 0)).This outcome seems sensible given how energy decays.
The sound pressure on the model's boundary surface in Figure 10a is displayed in Figure 12 at the frequency of 100 Hz, 300 Hz, and 500 Hz.For the model in Figure 10b, the sound pressure on the boundary surface is provided by Figure 13 at the frequencies of 100 Hz, 300 Hz, and 500 Hz.Three ways of view are given for each frequency.Considering that the plane wave moves in the direction of the x axis, the x − y and x − z planes' symmetrical properties are depicted in these figures, just like in Figures 12 and 13   (h) 3D view in the direction of (−1, 1, −1).(i) 3D view in the direction of (1, −1, −1).

Fish Model
A simplified thin-shell fish model is given in this section to conduct sensitivity analysis.Similar to the submarine models, the outside load is the incident wave, which has an amplitude of 1.0 and propagates down the x axis in a positive direction.Figure 14 shows the prototype example of the fish and its calculation model.The total length of the model is no more than 0.5 m and it keeps the basic shape of a manta ray.Due to the shape of the model, the calculation point is fixed at (5, 0, 0).Furthermore, the variable thickness is chosen for sensitivity analysis and the value is set at 0.011 m.
Before giving the frequency scan figure, we calculate the sound pressure at a fixed frequency 200 Hz to verify the numerical method.Figure 15 gives the sound pressure distribution on the surface of the fish model, with a subfigure of real part (left) and a subfigure of imaginary part (right) of the sound pressure.From the figure, it can be seen that the sound pressure distributes in strips and shows a slowly decreasing trend along the direction of incidence; obviously, this result is from the model structure, as that the model itself is symmetric.Moreover, the maximum sound pressure of both the real and imaginary parts is around 1, which coincides with the magnitude of the external action.Thus, the figure proves the correctness of the numerical method intuitively.
Then, the sensitivity analysis was conducted using thickness as independent variable.A frequency scan with 0.011 m thickness and sensitivity to thickness are given in Figure 16.In Figure 16b, sound pressure and sensitivity shows a similar pattern of increasing and then decreasing.Here, we find that, in some frequency points, the sound pressure changes sharply and this also results in the sharp change in sensitivity.This kind of error is mainly caused by resonance phenomenon, and statistics of these points are not taken into consideration in this work.Thus, ignoring the resonance phenomenon, the outcome also seems reasonable and once again proves the effectiveness of our method.The numerical simulations clearly show that the fluid effect needs to be considered while analyzing the vibro-acoustic radiation problem for underwater thin-shell constructions.Thus, it is necessary to perform the fluid-structure coupling analysis.Delineating high-quality meshes is crucial since the mesh quality has a direct impact on the computational correctness of the coupled analysis.This means that using IGA like NURBS to increase computational accuracy has significant benefits for both engineering and academia.

Conclusions
Sensitivity analysis and the modeling of the acoustic-structure interaction are conducted utilizing a coupling method grounded in the BEM and FEM.The structural elements of the problem are modeled using FEM.To avoid the necessity of meshing the acoustic space, the border of the structure being studied, which is also the acoustic domain's border, is discretized using the BEM.To speed up the matrix-vector output, the FMM is employed.The utilisation of NURBS IGABEM enables the direct examination of the sensitivity of the structural-acoustic interaction using CAD models, without the necessity for meshing, thereby removing geometric mistakes.In coupling structural-acoustic systems, sound pressure sensitivity equations are created.Numerical examples are shown to demonstrate the accuracy and practicality of the proposed method.In real-world scenarios, the proposed method might be used to quantitatively predict how design elements would affect the sound field.

Figure 5 .Figure 6 .
Figure 5.The relative error of sound pressure between the numerical and analytical solutions, at point (2, 0, 0) with frequency at 200 Hz, for the spherical shell model shown in Figure 4a.

Figure 8 .
Figure 8. Sound pressure on the spherical shell's boundary surface at a frequency of 100 Hz (the top row), 300 Hz (the middle row), and 500 Hz (the bottom row).(a) View along the −z axis.(b) View along the −x axis.(c) View along the x axis.(d) View in the direction of (1, 1, −1).(e) View along the −z axis.(f) View along the −x axis.(g) View along the x axis.(h) View in the direction of (1, 1, −1).(i) View along the −z axis.(j) View along the −x axis.(k) View along the x axis.(l) View in the direction of (1, 1, −1).

Figure 9 .
Figure 9. Sound pressure on the spherical shell's boundary surface at a frequency of 100 Hz (the top row), 300 Hz (the middle row), and 500 Hz (the bottom row), with the NURBS control point (−1, 0, 0) moved to (−1.5, 0, 0).(a) View along the −z axis.(b) View along the −x axis.(c) View along the x axis.(d) View in the direction of (1, 1, −1).(e) View along the −z axis.(f) View along the −x axis.(g) View along the x axis.(h) View in the direction of (1, 1, −1).(i) View along the −z axis.(j) View along the −x axis.(k) View along the x axis.(l) View in the direction of (1, 1, −1).

Figure 10 .
Figure 10.Control points distribution in the NURBS models of simple submarine.The left end control point is situated at (0, 0, 0).The positive direction of the x axis is from the model's left end to its right endpoint.(a) The NURBS model of a simple submarine.(b) The NURBS model of a partly expanded simple submarine.
Figure11a, Figure11c, and Figure11e, respectively, give the sound pressure comparison at location (5, 0, 0), (20, 0, 0) and (100, 0, 0), for the models shown in Figure10.The sound pressure's sensitivity in relation to the cylinder's radius at x = 0.5 m for these two models are given in Figure11b(at point (5, 0, 0), the sound pressure is the objective function), Figure11d(at point (20, 0, 0), the sound pressure is the objective function), and Figure11f(at point (100, 0, 0), the sound pressure is the objective function).In these figures, sound pressure and sensitivity show a similar pattern of increasing and then decreasing.Additionally, the value of sound pressure and sensitivity decrease as the separation between the model and the computation point increases (calculation point from (5, 0, 0) to (100, 0, 0)).This outcome seems sensible given how energy decays.The sound pressure on the model's boundary surface in Figure10ais displayed in Figure12at the frequency of 100 Hz, 300 Hz, and 500 Hz.For the model in Figure10b, the sound pressure on the boundary surface is provided by Figure13at the frequencies of 100 Hz, 300 Hz, and 500 Hz.Three ways of view are given for each frequency.Considering that the plane wave moves in the direction of the x axis, the x − y and x − z planes' symmetrical properties are depicted in these figures, just like in Figures12 and 13.