University of Birmingham Dynamic Pressure Analysis of Hemispherical Shell Vibrating in Unbounded Compressible Fluid

: This paper is the ﬁrst to highlight the vibrations of a hemispherical shell structure interacting with both compressible and incompressible ﬂuids. To precisely calculate the pressure of the shell vibrating in the air, a novel analytical approach has been established that has existed in very few publications to date. An analytical formulation that calculates pressure was developed by integrating both the ‘small-density method’ and the ‘Bessel function method’. It was considered that the hemispherical shell vibrates as a simple harmonic function, and the ﬂuid is non-viscous. For comparison, the incompressible ﬂuid model has been analyzed. Surprisingly, it is the ﬁrst to report that the pressure of the shell surface is proportional to the vibration acceleration, and the velocity amplitude decreased at the rate of 1 r 2 when the ﬂuid was incompressible. Otherwise, the surface pressure of the hemispherical shell was proportional to the vibration velocity, and the velocity amplitude decreased with the rate of 1 r when the ﬂuid was compressible. The compressibility of ﬂuid played an important role in the dynamic pressure of the shell structure. Furthermore, the scale factor derived by the theoretical approach was the product of the density and the sound velocity of the ﬂuid ( ρ o c ) exactly. In this study, the analytical solutions were veriﬁed by the calibrated numerical simulations, and the analytical formulation were rigorously tested by extensive parametric studies. These new ﬁndings can be used to guide the optimal design of the spherical shell structure subjected to wind load, seismic load, etc. wind load of the spherical shell structure. The novel analytical formulation that was formed in this study can be used to predict the vibration of a spherical shell structure interacting with ﬂuids. This new formulation can be applied to the design and analysis of the shell structure in practice.


Introduction
A hemispherical shell is widely used in large span structures, which have been constructed in large public buildings such as gymnasiums, exhibition centers, and opera houses [1,2]. Such curved shapes are often used for steel shells and rod structure as they are naturally strong structures, allowing wide areas to be spanned without the use of internal supports, giving an open, unobstructed interior. Figure 1 shows a geodesic dome for digital artwork that was established in 2014 [3]. The use of steel as a building material reduces both material costs and construction costs, as steel is relatively inexpensive and easily made into compound curves. The shell structure tends to be immensely strong and safe. For example, modern monolithic dome houses have resisted hurricanes and fires. In recent literature, there has been renewed interest in the problem of a hemispherical shell vibrating in contact with air or other fluids [4]. The vibrations of a structure induce local air movement, while the moving air There are some publications in the literature devoted to the interaction of fluid and shell by theoretical approaches [5,6]. In the way of potential flow theory, Minami [7] derived the added mass ratio for a rectangular plate under the assumption that the vibration mode is simple harmonic. Actually, the researcher utilized the incompressible assumption so the obtained result was an added mass ratio. Wang [8] reanalyzed the same problem with the assumption that the structure vibrates as a quadratic function and raised another added mass formula. In both papers, the problem was actually treated as a one-dimensional model so that their method cannot be applied to twodimensional (2D) and three-dimensional (3D) problems. Yadykin [9] studied the effect of the added mass of a thin shell vibrating in a stationary fluid through using the energy method, and the pressure of the shell was calculated by the thin-airfoil theory. Ping and Wei [10] studied the problem of a beam vibrating in the fluid based on the classical Fourier series, and reported the effects of the added mass ratio on the one-dimensional beam vibrating in non-viscous fluid.
Most of the research about the vibrations of shell in fluids commonly studied them by a numerical method, especially by the computation fluid dynamics (CFD) method [11,12]. Sorokin and Terentiev [13] studied the effect of the generation and transmission of the vibro-acoustic energy in an elastic cylindrical shell filled with water. They focused more on vibro-acoustic energy than the pressure. Chen et al. [14] studied the Guangzhou International Convention and Exhibition Center with dimensions that were 210 m wide and 457 m long. They studied the pressure distribution of the exhibition according to the wind tunnel test. Stolarski [15] conducted some research studies on the pressure of compressible squeeze film. In the study, the CFD method was used. Jeong and Kim [16,17] developed a Rayleigh-Ritz method to study a proportional hydroelastic vibration of two annular plates coupled with a bounded fluid. They considered the effects of the compressibility of the fluid and obtained a theoretical result that could predict the fluid-coupled frequencies well. Kubenko et al. [18] conducted a number of studies on a cylindrical shell submerged in an unbounded medium. They also assumed the fluid to be an incompressible fluid. Eftekhari [19] developed a differential quadrature method (DQM) to investigate the problem of the free vibration of circular plates in contact with fluid. Hu and Tang [20,21] investigated a simple-based monolithic implicit method for the strong-couple FSI problem. Their method may be used to analyze the vibration of the shell. These methods consist of several numerical approaches, but they offered little findings on a basis of theoretical analysis.
In spite of previous research on shells in contact with an incompressible fluid by experimental or numerical methods, very few studies on shells vibrating in compressible fluid have been performed. There have also been no attempts to tackle the problem of hemispherical shell vibration, which is highlighted in this paper. Moreover, up to now, very few studies have reported how many errors occurred owing to the incompressible fluid assumption. According to previous studies, our research is the first to consider the compressibility of fluid, which cannot be neglected. In the present paper, the pressure of a hemispherical surface vibrating under both compressible and incompressible There are some publications in the literature devoted to the interaction of fluid and shell by theoretical approaches [5,6]. In the way of potential flow theory, Minami [7] derived the added mass ratio for a rectangular plate under the assumption that the vibration mode is simple harmonic. Actually, the researcher utilized the incompressible assumption so the obtained result was an added mass ratio. Wang [8] reanalyzed the same problem with the assumption that the structure vibrates as a quadratic function and raised another added mass formula. In both papers, the problem was actually treated as a one-dimensional model so that their method cannot be applied to two-dimensional (2D) and three-dimensional (3D) problems. Yadykin [9] studied the effect of the added mass of a thin shell vibrating in a stationary fluid through using the energy method, and the pressure of the shell was calculated by the thin-airfoil theory. Ping and Wei [10] studied the problem of a beam vibrating in the fluid based on the classical Fourier series, and reported the effects of the added mass ratio on the one-dimensional beam vibrating in non-viscous fluid.
Most of the research about the vibrations of shell in fluids commonly studied them by a numerical method, especially by the computation fluid dynamics (CFD) method [11,12]. Sorokin and Terentiev [13] studied the effect of the generation and transmission of the vibro-acoustic energy in an elastic cylindrical shell filled with water. They focused more on vibro-acoustic energy than the pressure. Chen et al. [14] studied the Guangzhou International Convention and Exhibition Center with dimensions that were 210 m wide and 457 m long. They studied the pressure distribution of the exhibition according to the wind tunnel test. Stolarski [15] conducted some research studies on the pressure of compressible squeeze film. In the study, the CFD method was used. Jeong and Kim [16,17] developed a Rayleigh-Ritz method to study a proportional hydroelastic vibration of two annular plates coupled with a bounded fluid. They considered the effects of the compressibility of the fluid and obtained a theoretical result that could predict the fluid-coupled frequencies well. Kubenko et al. [18] conducted a number of studies on a cylindrical shell submerged in an unbounded medium. They also assumed the fluid to be an incompressible fluid. Eftekhari [19] developed a differential quadrature method (DQM) to investigate the problem of the free vibration of circular plates in contact with fluid. Hu and Tang [20,21] investigated a simple-based monolithic implicit method for the strong-couple FSI problem. Their method may be used to analyze the vibration of the shell. These methods consist of several numerical approaches, but they offered little findings on a basis of theoretical analysis.
In spite of previous research on shells in contact with an incompressible fluid by experimental or numerical methods, very few studies on shells vibrating in compressible fluid have been performed. There have also been no attempts to tackle the problem of hemispherical shell vibration, which is highlighted in this paper. Moreover, up to now, very few studies have reported how many errors occurred owing to the incompressible fluid assumption. According to previous studies, our research is the first to consider the compressibility of fluid, which cannot be neglected. In the present paper, the pressure of a hemispherical surface vibrating under both compressible and incompressible fluid are analyzed using the Bessel function method. The principle of pressure on the shell surface is discussed, and some formulas about the velocity and pressure are highlighted. All of the results are verified using the calibrated numerical simulations. The outcome of this study will enable a novel analytical formulation that can be used to predict the vibration of a spherical shell structure interacting with fluids. This formulation can be used in the design and analysis of the structure in practice. Figure 2 shows a mechanical model of a hemispherical shell vibrating in a radial direction in the unbounded fluid. Obviously, structure vibration affected the boundary of flow, while the pressure of flow affected the structure simultaneously. With an inviscid fluid hypothesis, the assumption of p = c 2 ρ, and the character of axial symmetry, the FSI coupling equation can be simplified as follows [22], where m represents the areal density of the shell, U represents the displacement of vibration, k represents the bending stiffness, u represents the velocity of fluid in the radius direction, and ∆p represents the pressure differential between the upper surface and the bottom surface of the shell.

Formulation for Hemispherical Shell
Equation (1) is the structural control equation of vibration. The first term of Equation (1) presents the inertial force (pressure) derived from the acceleration of motion, and the second term of Equation (1) shows the elastic force contributed from the deformation of the shell. According to Newton's second law, the combination of inertial force, elastic force, and pressure without damping force are in the condition of balance. The boundary conditions can be written as Equation (4), where the first term condition means the displacement is zero at the initial time, and the second term shows that the initial velocity of the shell is in compliance to the u s (0) at the initial time.
Equations (2) and (3) are fluid control equations, which are well known as Naiver-Stokes equations (3D axial symmetry), where the density, static pressure, sound speed, and velocity of the fluid is denoted by ρ, p, c, and u, respectively. Equation (2) was derived from the conservation of mass, and Equation (3) was deduced from the conservation of momentum. Let: The velocity of fluid boundary vibration is, u = u s (t). Here, it assumes that |u s | = c. Therefore, the boundary conditions can be rewritten as Equation (6), where, the first term shows that the boundary of fluid is just the position of the shell, and both the second and third terms present that on the boundary of the far field, the velocity and pressure are approximate to the initial condition 0, p 0 , respectively. u| r=r s = u s (t), u| r=R−>∞ = 0, p| r=R−>∞ = p 0 (6) According to the physical conditions, the boundary of the fluid is the structural shell itself actually; thus, the initial value of the radius is the radius of the hemispherical shell, r s , therefore: where r s (t) is the radius of the fluid boundary, which varies with time. Since U is taken into account in a small value, r s (t) can be denoted as r s . The initial condition can be expressed by: According to the physical conditions, the boundary of the fluid is the structural shell itself actually; thus, the initial value of the radius is the radius of the hemispherical shell, s r , therefore: where ( ) s r t is the radius of the fluid boundary, which varies with time. Since U is taken into account in a small value, ( ) s r t can be denoted as s r . The initial condition can be expressed by:

Solution on the Compressible Fluid Case
Using the Fourier series method, the velocity can be written in a harmonic function with a frequency of ω . This implies that: For the physical property of the fluid, the density changes slightly, but the change of density over time and space cannot be neglected. Thus, in this paper, the density is assumed to be a constant as ρ o , but the changes of density over time and space have to be considered in the solution. From Equation (2), pressure can be obtained by the integration between limits from r to the unbounded position: From Equation (10), it is found that the pressure comes from two parts: the local acceleration, which is the velocity changes with respect to time, and the convective acceleration, which is the change in velocity with respect to space. It should be noted that there is only one boundary condition with respect to the pressure as prescribed in Equation (6). Therefore, the pressure at the position of r is an unknown value. This implies that there is an arbitrary constant that needs to be determined in Equation (10). For the assumption of s u c  , the second term of Equation (10) is very small. After neglecting the convective acceleration (the second term) and substituting the local acceleration (the first term) to the continuous equation of Equation (2), it can be obtained as follows:

Solution on the Compressible Fluid Case
Using the Fourier series method, the velocity can be written in a harmonic function with a frequency of ω. This implies that: For the physical property of the fluid, the density changes slightly, but the change of density over time and space cannot be neglected. Thus, in this paper, the density is assumed to be a constant as ρ o , but the changes of density over time and space have to be considered in the solution. From Equation (2), pressure can be obtained by the integration between limits from r to the unbounded position: From Equation (10), it is found that the pressure comes from two parts: the local acceleration, which is the velocity changes with respect to time, and the convective acceleration, which is the change in velocity with respect to space. It should be noted that there is only one boundary condition with respect to the pressure as prescribed in Equation (6). Therefore, the pressure at the position of r is an unknown value. This implies that there is an arbitrary constant that needs to be determined in Equation (10). For the assumption of u s c, the second term of Equation (10) is very small. After neglecting the convective acceleration (the second term) and substituting the local acceleration (the first term) to the continuous equation of Equation (2), it can be obtained as follows: Equation (11) can be simplified as Equation (12) below after a derivative by r. Having an arbitrary constant in pressure discussed above, two boundary conditions are necessary to solve Equation (11), although it is a first order differential equation. The two boundary conditions are the same as Equation (12): Let s = ω c r, and the above equation can be expressed in terms of the classical Bessel equation as follows: Equation (13) is a classical Bessel equation, so the solution can be derived as: where A and B are the arbitrary functions of the time, and j 1 (s), y 1 (s) are the first and the second-order Bessel function, respectively. Applying the boundary conditions into Equation (14), it can be rewritten as Equation (15) shows: Thus, to get the exact solution, there are only A and B left to be confirmed. It should be noted that when r = R → ∞ , both j 1 (s R ) and y 1 (s R ) converge to zero. Therefore, the second term of the boundary condition is satisfied naturally. In such cases, taking into account that u is a harmonic function, it can be assumed that: According to the first term of the boundary conditions, the amplitude and phase angle will be equal to that of initial vibration, so that A and B can be confirmed. In addition, it is found that the solution can be determined totally by the boundary condition u s (t).

The Analysis of the Velocity
Due to the physical boundary conditions, the vibration at the boundary can be expressed as: where D represents the amplitude of vibration. So, the first term of the boundary condition can be expressed as Equation (18): Afterward, it may be simplified as follows: where, s r = ω c r r , and tan(φ) = j 1 (s r ) y 1 (s r ) . Furthermore, it may be obtained that: By the nature of the Bessel function [23,24], if s is very large, then there are: Finally, the solution can be derived as Equation (22): Accordingly, it can be concluded that the characteristic of vibration at time t 0 will transfer to the position r = ct with the same frequency between time t, and the amplitude will be decreased with the rate of 1 r .

The Analysis of the Pressure
By the characters of the Bessel function, it also can be written that: Thus, the formula about pressure can be obtained: where s, ϕ, D are the same as the Section 2.3. if s is large, there are: Finally, the solution of the pressure can be induced: It can be seen from Equation (26) that the pressure is proportional to the velocity of vibration and the ratio is the product of the density and the sound speed. Furthermore, according to Equation (22), the convective acceleration is 1 2 ρ o u 2 , which is the commonly used pressure formula derived by Bernoulli law [25]. Besides, compared to Equation (26), the value of the convective acceleration is very small. This means that neglecting the value of convective acceleration cannot cause a large error.

Incompressible Fluid Solution
As for a comparison, it is necessary to consider the model of a shell vibrating in the incompressible fluid. The control equation can be expressed as follows: ∂u ∂r ∂u ∂t + u ∂u ∂r Equation (27) is the structural control equation, which is the same as Equation (1). Equation (28) is the continuity equation with the constant density. Equation (29) is the equation for the conservation of momentum. The symbols and the boundary conditions of structure are the same as Equation (4), whilst the boundary of the fluid domain is as follows: According to the continuity equation, it may be concluded as Equation (31): While according to the momentum equation, it can be written as: Substituting the boundary condition of the far field to Equation (32), it can be found that: From the solution of Equation (31), in the condition of incompressible fluid, it is found that the motion of boundary will transfer to the far field simultaneously, and the amplitude decreases with the rate of 1 r 2 . According to Equation (33), it can be observed that the pressure is contributed from two parts [26]. The first term of pressure can be seen as added mass, and the second term is also convective acceleration, which is the same as Equation (5). The pressure under the shell is ∆p = ρ 0 r s ∂u s ∂t + 1 2 ρ 0 u 2 s , which is very different to the solution of compressible fluid condition.

Discussion
For the incompressible fluid condition, Equation (33) shows that it has an important influence, as the incompressible assumption was adopted, and the pressure is contributed from two parts. The first part is the local acceleration, which is proportional to both the boundary acceleration and the reciprocal of the radius of the fluid domain. Another part is proportional to the square of the fluid velocity. As for compressible fluid, according to the solution, ∆p = ρ o c · u, it can be concluded that the pressure is proportional to the velocity of fluid. The ratio is the product of the sound speed and the density. In such case, it can be obtained from Equation (1) that the pressure effect of fluid on the structure is equivalent to the viscous damping whose value is ρ o c, exactly.
In the physical sense, it can be seen that between times ∆t, the amount of fluid that can be squeezed by the structure is 4 3 π[(r s + u s ∆t) 3 − r 3 s )], and can be spread to the range of 4 3 When ∆t → 0 , the fluid compression ratio is: Moreover, the fluid compression modulus is ρ o c 2 , Therefore, according to Hooke's law, the pressure of compressible fluid is: The pressure of compressible fluid is identical to Equation (26) exactly. The discussion above shows that in the condition of compressible fluid, the motion of boundary propagates to the far field at the speed of sound. At the beginning, the motion of the shell compresses the air (fluid), which will be transferred to the far field afterward. The compression of air produces the local acceleration pressure, which is proportional to the fluid velocity.

CFD Method
It is well-known that the second-order method is the most stable discretization scheme, which is derived by Taylor's series: Equations (36) and (37) are one and two-dimensional schemes, respectively, in which u is the value of physical quantities such as velocity and pressure, and x, y are the coordinates. Usually, the grid will not be the rectangular shape. Let the transformation of domain be as Equation (36), shown in Figure 3a.
In that case, the differential can be obtained by the chain rule, and the second-order differential cab be rewritten as Equations (39) and (40): On the basis of the preceding analysis, in order to find and check the validity of the theory, a hemispherical shell submerged in unbounded fluid was examined using a CFD analysis using a commercial computer code, ANSYS (release 17.0, Ansys Company, Canonsburg, PA, USA) [27]. The CFD model has the identical shell geometry, boundary conditions, and material properties that will be used in the theoretical calculation. Similar applications can be found in some papers [28,29]. The radius of the hemispherical shell is set to 10 m, and the unbounded fluid is simulated as a far field boundary whose radius is 500 m (see Figure 3b). The physical properties of the fluid are as follows: density = 1.225 kg/m 3 , sound velocity c = 340 m/s. The viscosity of the fluid is neglected in both the theoretical analysis and the CFD calculation. The purpose of the present paper is to verify the principle of pressure under the vibration of the shell. The movement of the hemispherical shell was set as uniform vibration in the radial direction, which was also the boundary of the fluid. As for the numerical simulation parameters, the far field boundary is set to be the pressure far field, reference pressure to zero, movement boundary to U(t) = sin(10t), calculating time to 1.0 s, and the delta-time to 0.001 s. A simple scheme second-order method was accepted for the CFD simulation, and the initial condition of fluid was set to be stationary. For comparison, both the compressible fluid case and the incompressible fluid case were carried out. To verify the solution, the line from r s to R, which data was obtained from, is shown as LINE-X in Figure 3b.  In order to demonstrate the influence of the discretization, the numerical result with different meshing sizes is displayed. Table 1 shows the comparison of the total energy of all over the fluid at time = 0.2 s with different density of shell mesh. Since the results of size = 0.2 m and size = 0.1 m are the same, the size of 0.2 m has been adopted in this paper. After meshing, there are 363,923 structural hexahedron elements.  In order to demonstrate the influence of the discretization, the numerical result with different meshing sizes is displayed. Table 1 shows the comparison of the total energy of all over the fluid at time = 0.2 s with different density of shell mesh. Since the results of size = 0.2 m and size = 0.1 m are the same, the size of 0.2 m has been adopted in this paper. After meshing, there are 363,923 structural hexahedron elements.     Compared with both left and right plot in Figure 5, it can be found that the velocity of fluid in the incompressible fluid is a function of radii, which decreases with the function of 1 r 2 . Therefore, it is found that the velocity far away from boundary shell trends to converge to zero quickly. On the other hand, the velocity in the compressible fluid is a function of both radii and time and the periodicity of radius from formula is 2πc ω = 213, which is verified by the result in Figure 4 (left). Besides, the velocity amplitude decreases with the function of 1 r c, which is slower than that in the incompressible fluid condition.

Parametric Analyses
To extend the applications of the analytical formulations, two more hemispherical shells rotated by the generatrix of the elliptical curve were conducted, which was inspired by the previous studies [30][31][32][33]. The length of major axis A was set to 10 m, which is the same as the conditions of hemisphere shell above, and the length of minor axis B was set to 8 m and 6 m, respectively. To compare the difference in the major axis and minor axis, the results of two lines, including X ranging from 10 m to 300 m with Y = 0, which is called LINE-X, and LINE-Y ranging from 8/6 to 300 m with X = 0, which is called LINE-Y, were investigated below.
4.1. Ellipse Geometry, A = 10, B = 8 Figure 9 shows the plot of velocity versus time on the incompressible and compressible fluid condition along LINE-X and LINE-Y, respectively, when time = 1.0 s. Figure 10 presents the contour slice of velocity ranging from 10 m to 50 m, as the values of the velocity that are far from 50 m are too small to distinguish it. It is noted that there is a small deviation at about 55 m. This may be caused by the error of the numerical calculation, since the mesh size at the place far from the boundary is not small enough. It can be found that there is no evident difference regarding velocity. Also, the data is compliant with the rules that this paper presents.

Parametric Analyses
To extend the applications of the analytical formulations, two more hemispherical shells rotated by the generatrix of the elliptical curve were conducted, which was inspired by the previous studies [30][31][32][33]. The length of major axis A was set to 10 m, which is the same as the conditions of hemisphere shell above, and the length of minor axis B was set to 8 m and 6 m, respectively. To compare the difference in the major axis and minor axis, the results of two lines, including X ranging from 10 m to 300 m with Y = 0, which is called LINE-X, and LINE-Y ranging from 8/6 to 300 m with X = 0, which is called LINE-Y, were investigated below. Figure 9 shows the plot of velocity versus time on the incompressible and compressible fluid condition along LINE-X and LINE-Y, respectively, when time = 1.0 s. Figure 10 presents the contour slice of velocity ranging from 10 m to 50 m, as the values of the velocity that are far from 50 m are too small to distinguish it. It is noted that there is a small deviation at about 55 m. This may be caused by the error of the numerical calculation, since the mesh size at the place far from the boundary is not small enough. It can be found that there is no evident difference regarding velocity. Also, the data is compliant with the rules that this paper presents.      Figure 11 shows the plot of pressure versus time of the incompressible and the compressible fluid along LINE-X and LINE-Y, respectively. The results of the numerical solution and theoretical prediction agree well. In addition, Figure 12 presents the contour slice of velocity ranging from 10 m to 50 m, since the data values of the velocity that are far from 50 m are too small to be distinguished.    Figure 11 shows the plot of pressure versus time of the incompressible and the compressible fluid along LINE-X and LINE-Y, respectively. The results of the numerical solution and theoretical prediction agree well. In addition, Figure 12 presents the contour slice of velocity ranging from 10 m to 50 m, since the data values of the velocity that are far from 50 m are too small to be distinguished.     Figure 13 shows the curve of the pressure versus velocity in the numerical solution and the analytical solution. In this situation, the slope of the fitting curve is 427 with 95% confidence bounds, and the predicted value is 407, which is consistent with the expected results. It is noticeable that the smaller the radius, the larger the gap between the analytical and numerical result. That is because when s , which is linear to R, is smaller, is the approximations that are shown in Equations (33) and (15) are worse.

Elliptic Geometry A = 10, B = 6
In the situation of A = 10, B = 6, the simulation with the identical parameter is carried out. To compare with the theoretical result, the data of t = 0.65 s is shown below. Figure 14 plots the figure of velocity versus time in incompressible and compressible fluid conditions along LINE-X and LINE-Y, respectively. To a certain degree, the flow field will be a nonuniform flow, so it can be seen from Figure 14a that there is some difference between the data values obtained from LINE-X and LINE-Y.  Figure 13 shows the curve of the pressure versus velocity in the numerical solution and the analytical solution. In this situation, the slope of the fitting curve is 427 with 95% confidence bounds, and the predicted value is 407, which is consistent with the expected results. It is noticeable that the smaller the radius, the larger the gap between the analytical and numerical result. That is because when s, which is linear to R, is smaller, is the approximations that are shown in Equations (33) and (15) are worse.  Figure 13 shows the curve of the pressure versus velocity in the numerical solution and the analytical solution. In this situation, the slope of the fitting curve is 427 with 95% confidence bounds, and the predicted value is 407, which is consistent with the expected results. It is noticeable that the smaller the radius, the larger the gap between the analytical and numerical result. That is because when s , which is linear to R, is smaller, is the approximations that are shown in Equations (33) and (15) are worse.

Elliptic Geometry A = 10, B = 6
In the situation of A = 10, B = 6, the simulation with the identical parameter is carried out. To compare with the theoretical result, the data of t = 0.65 s is shown below. Figure 14 plots the figure of velocity versus time in incompressible and compressible fluid conditions along LINE-X and LINE-Y, respectively. To a certain degree, the flow field will be a nonuniform flow, so it can be seen from Figure 14a that there is some difference between the data values obtained from LINE-X and LINE-Y.

Elliptic Geometry A = 10, B = 6
In the situation of A = 10, B = 6, the simulation with the identical parameter is carried out. To compare with the theoretical result, the data of t = 0.65 s is shown below. Figure 14 plots the figure of velocity versus time in incompressible and compressible fluid conditions along LINE-X and LINE-Y, respectively. To a certain degree, the flow field will be a non-uniform flow, so it can be seen from Figure 14a that there is some difference between the data values obtained from LINE-X and LINE-Y.  Figure 13 shows the curve of the pressure versus velocity in the numerical solution and the analytical solution. In this situation, the slope of the fitting curve is 427 with 95% confidence bounds, and the predicted value is 407, which is consistent with the expected results. It is noticeable that the smaller the radius, the larger the gap between the analytical and numerical result. That is because when s , which is linear to R, is smaller, is the approximations that are shown in Equations (33) and (15) are worse.

Elliptic Geometry A = 10, B = 6
In the situation of A = 10, B = 6, the simulation with the identical parameter is carried out. To compare with the theoretical result, the data of t = 0.65 s is shown below. Figure 14 plots the figure of velocity versus time in incompressible and compressible fluid conditions along LINE-X and LINE-Y, respectively. To a certain degree, the flow field will be a nonuniform flow, so it can be seen from Figure 14a that there is some difference between the data values obtained from LINE-X and LINE-Y.  Figure 15 illustrates the contour slice of the pressure in the incompressible and compressible fluid conditions, respectively. From Figure 15, it is clear that in the middle of the shell, the velocity is larger than that of the side. Since in this situation, the slope at the side is larger than that at the middle, the fluid can distribute within the other domain close to the side boundary. In this case, the velocity of fluid decreases faster than that of the middle position.   Figure 17 depicts the contour slice of the pressure in the incompressible and compressible fluid conditions. To a certain degree, the flow will be a non-uniform flow, so it can be seen from Figure 16a that the values in LINE-X and LINE-Y are obviously different.  Figure 15 illustrates the contour slice of the pressure in the incompressible and compressible fluid conditions, respectively. From Figure 15, it is clear that in the middle of the shell, the velocity is larger than that of the side. Since in this situation, the slope at the side is larger than that at the middle, the fluid can distribute within the other domain close to the side boundary. In this case, the velocity of fluid decreases faster than that of the middle position.  Figure 15 illustrates the contour slice of the pressure in the incompressible and compressible fluid conditions, respectively. From Figure 15, it is clear that in the middle of the shell, the velocity is larger than that of the side. Since in this situation, the slope at the side is larger than that at the middle, the fluid can distribute within the other domain close to the side boundary. In this case, the velocity of fluid decreases faster than that of the middle position.   Figure 17 depicts the contour slice of the pressure in the incompressible and compressible fluid conditions. To a certain degree, the flow will be a non-uniform flow, so it can be seen from Figure 16a that the values in LINE-X and LINE-Y are obviously different.   Figure 17 depicts the contour slice of the pressure in the incompressible and compressible fluid conditions. To a certain degree, the flow will be a non-uniform flow, so it can be seen from Figure 16a that the values in LINE-X and LINE-Y are obviously different.  Figure 15 illustrates the contour slice of the pressure in the incompressible and compressible fluid conditions, respectively. From Figure 15, it is clear that in the middle of the shell, the velocity is larger than that of the side. Since in this situation, the slope at the side is larger than that at the middle, the fluid can distribute within the other domain close to the side boundary. In this case, the velocity of fluid decreases faster than that of the middle position.   Figure 17 depicts the contour slice of the pressure in the incompressible and compressible fluid conditions. To a certain degree, the flow will be a non-uniform flow, so it can be seen from Figure 16a that the values in LINE-X and LINE-Y are obviously different.  Figure 18 shows the relationship of pressure versus velocity and the correspondent fitting results. In the conditions of LINE-Y and LINE-X, the slope of fitting is 430 and 477, respectively. Based on the theoretical results, the slope is 417 without any change. This implies that the more irregular the shell, the larger the error between numerical and theoretical values will be. Besides, it is also found that the smaller the R, the worse the approximation will be.

Conclusions
A new analytical approach to calculate the pressure of the hemispherical shell vibrating in compressible fluid has been firstly developed. It is rigorously verified by the theoretical approach  Figure 18 shows the relationship of pressure versus velocity and the correspondent fitting results. In the conditions of LINE-Y and LINE-X, the slope of fitting is 430 and 477, respectively. Based on the theoretical results, the slope is 417 without any change. This implies that the more irregular the shell, the larger the error between numerical and theoretical values will be. Besides, it is also found that the smaller the R, the worse the approximation will be.

Conclusions
A new analytical approach to calculate the pressure of the hemispherical shell vibrating in compressible fluid has been firstly developed. It is rigorously verified by the theoretical approach  Figure 18 shows the relationship of pressure versus velocity and the correspondent fitting results. In the conditions of LINE-Y and LINE-X, the slope of fitting is 430 and 477, respectively. Based on the theoretical results, the slope is 417 without any change. This implies that the more irregular the shell, the larger the error between numerical and theoretical values will be. Besides, it is also found that the smaller the R, the worse the approximation will be.  Figure 18 shows the relationship of pressure versus velocity and the correspondent fitting results. In the conditions of LINE-Y and LINE-X, the slope of fitting is 430 and 477, respectively. Based on the theoretical results, the slope is 417 without any change. This implies that the more irregular the shell, the larger the error between numerical and theoretical values will be. Besides, it is also found that the smaller the R, the worse the approximation will be.

Conclusions
A new analytical approach to calculate the pressure of the hemispherical shell vibrating in compressible fluid has been firstly developed. It is rigorously verified by the theoretical approach

Conclusions
A new analytical approach to calculate the pressure of the hemispherical shell vibrating in compressible fluid has been firstly developed. It is rigorously verified by the theoretical approach based on the compressible fluid assumption, which can predict the values of the velocity and pressure very well. It is found that the assumption of the compressibility of fluid plays a main role in the velocity and pressure of the FSI system. In addition, the exact solution and theoretical formulation is obtained by an integrated approach. According to the solutions, it is found that the pressure on the shell surface is proportional to the acceleration of boundary vibration when the fluid is assumed to be incompressible. Meanwhile, under the assumption of compressible fluids, the boundary pressure is proportional to the velocity of the boundary vibration. Considering the velocity of fluids, the amplitude of velocity decreases with the rate of 1/r in the compressible condition, whilst it decreases with the rate of 1/r 2 in the incompressible condition.
In addition, the physical implication of the developed formula has been discussed. It is found that the motion of the boundary will propagate to the far field at the speed of sound velocity with the same frequency, but will decrease with distant r. Initially, the motion of the shell compresses the air locally, which will transfer to the far field afterward. The compression of air produces the local acceleration pressure, which is proportional to fluid velocity. The results reveal that the compression force can be exactly expressed by the new formula.
Furthermore, it can be clarified that incompressible fluid assumption for fluids can induce a large numerical error in various cases. It is also found that under the incompressible fluid assumption, the effect of air pressure is equivalent to the effect of the added mass, which confirms the scientific findings by many previous research studies and common theories that are applied today. While under the compressible fluid consideration, the effect of air pressure is equivalent to the effect of viscous damping, whose ratio is exactly ρ o c. Therefore, our present findings and established formula can be used in the design and application of the resisting wind load of the spherical shell structure. The novel analytical formulation that was formed in this study can be used to predict the vibration of a spherical shell structure interacting with fluids. This new formulation can be applied to the design and analysis of the shell structure in practice.