Development of indirect EFBEM for radiating noise analysis including underwater problems
- Author: Kwon Hyun-Wung, Hong Suk-Yoon, Song Jee-Hun
- Organization: Kwon Hyun-Wung; Hong Suk-Yoon; Song Jee-Hun
- Publish: International Journal of Naval Architecture and Ocean Engineering Volume 5, Issue3, p392~403, 30 Sep 2013
For the analysis of radiating noise problems in medium-to-high frequency ranges, the Energy Flow Boundary Element Method (EFBEM) was developed. EFBEM is the analysis technique that applies the Boundary Element Method (BEM) to Energy Flow Analysis (EFA). The fundamental solutions representing spherical wave property for radiating noise problems in open field and considering the free surface effect in underwater are developed. Also the directivity factor is developed to express wave’s directivity patterns in medium-to-high frequency ranges. Indirect EFBEM by using fundamental solutions and fictitious source was applied to open field and underwater noise problems successfully. Through numerical applications, the acoustic energy density distributions due to vibration of a simple plate model and a sphere model were compared with those of commercial code, and the comparison showed good agreement in the level and pattern of the energy density distributions.
Energy flow boundary element method (EFBEM) , Radiating noise , Free surface effect , Directivity factor.
When ships or underwater objects move through water, whose impedance is much higher than that of air, small vibrating motions of the inner ship body can cause remarkable noises in wide frequency ranges. No single analysis method of the noise phenomenon can be effectively applied to all ranges of noise problems. At low frequency ranges, the analysis of vibration problems is analyzed by the conventional methods such as the Finite Element Method (FEM) and Boundary Element Method (BEM). But, at high frequency ranges, those methods require more computation time and costs, and thus alternative methods are needed.
Among many alternative methods, EFA has received much attention. This method was introduced by Belov et al. (1997) in 1997. Nefske and Sung (1989) implemented Energy Flow Finite Element Method (EFFEM) to solve the transverse vibration of a beam. Wholever and Bernhard (1992) derived the energy governing equation for diverse vibrating waves of a beam. Bouthier and Bernhard (1992) derived the energy governing equation of a flexural wave, considering only the far-field component of the transverse vibration of a membrane and thin plate. They expanded the research of EFA to vibration problems of twodimensional structures. Cho (1993) researched the energy boundary condition in the connection point between elements of a beam by applying the wave transmission approach to a complex structure. Park (1999) and Park et al. (2001) derived the energy governing equation for an in-plane wave of a thin plate and studied the spatial distribution and transmission path of vibration energy for a plate structure, which is connected at a non-determined angle. Seo et al. (2003) expanded the application of EFA to beam-plate coupled structures. Lee et al. (2008) applied EFBEM to the vibration analysis of beam and plate problems. Also Wang et al. (2004) applied energy boundary element formulation to the sound radiation problems.
In this paper, the energy governing equation having spherical wave property is developed in open field. And the directivity effect which is represented in high frequency range problems but can’t express in EFA is studied. And the fundamental solution and energy governing equation for underwater problems are developed. The developed equation and directivity effect are applied to the simple case and the results are compared with commercial noise analysis program, SYSNOISE and reliable results are obtained.
In a linearly elastic medium, the energy balance equation is derived from the following Eq. (1). The amount of incoming and out coming power through the surface of an object and the rate of change of the total energy in the object are same. From this fact, energy balance equation is represented as follows:
eis the total energy density in the volume of an object, CV.
is the displacement vector at the boundary of an object, CS.
is a small area vector perpendicular to the surface of the boundary.
σrepresents the stress on the surface of the boundary. πinand πdissare expressed as input power and loss power acting on a unit volume at a unit time, respectively. Intensity, defined as the power per unit area flowing out of an object, can be obtained from the stress acting on the surface of the boundary and the velocity of the surface of the boundary as follows:
is the intensity, the power per unit area. In Eq. (1), the first integral on the right hand side is rewritten with the application of Gauss’s theorem as follows:
Therefore, the energy balance equation in the volume of an object is obtained by
From Eq. (4), the energy balance equation for a small volume is expressed by
Eq. (5) is the energy balance equation for all elastic mediums, steady state and transient state. For a steady state elastic medium, the rate of energy density with respect to the time is zero; the left term in Eq. (5) is zero. Therefore, the steady-state energy balance equation is represented as follows:
From Eq. (6), the input power due to the exterior exciting force is expressed as the sum of the power lost in the object and the power flowing out the adjacent medium.
The loss of vibration energy can be represented by a damping structure model. Cremer et al. (2004) showed that energy density loss,
ediss, due to the damping structure per one period at any point in an elastic medium vibrating at frequency ωis proportional to reversible vibration energy density, eR. Therefore, Eq. (7) is obtained by
If the inner damping coefficient
ηis very small, the reversible energy density, eR, can be substituted for the time average value of the total energy density, < e> , which is the sum of the kinetic energy density and the potential energy density. Therefore, Eq. (8) is given by
where, < > means the time average. The time average of power loss during one period, <
πdiss> can be obtained by dividing the energy loss during one period, edissby period, Tbeing 2 π/ ω. And it can be obtained approximately by using the time average of the total energy density, < e>, as follows:
Eq. (9) is the energy loss equation, which is derived from the application of structural damping to any point on an elastic medium vibrating at frequency
ω. We assumed that the kinetic energy is the same as the potential energy.
For analyzing the radiating noise problems underwater, we must consider the effect of free surface. According to Fig. 1 and reflection coefficient of surface, spherical wave in underwater is represented as follows:
r1 is the direct distance between a source point and field point. r2 is the distance arrived from the source point to the field point which is reflected at the free surface of the water.
is the wave number considering the damping effect.
ηis damping loss factor and kis the wave number. From the Euler’s equation, the velocity of the wave can be obtained.
Therefore, energy density and intensity can be expressed as follows for a far-field with low damping:
From Eqs. (11) and (12), the following energy transmission equation can be obtained as follows:
D= 1 + e-2 ηkhcos θ？2cos(2 khcos θ) e- ηkhcos θ.
From the energy balance equation, Eq. (6), the energy loss equation, Eq. (9), and the energy transmission equation, Eq. (13), the energy governing equation can be obtained as follows:
Eq. (14) is the energy governing equation representing the property of a spherical wave and including the effect of the water surface.
The fundamental solution
Gis the exact solution for input power, as well as the fundamental solution representing the relation between the virtual source and energy density. Therefore, the shape of the fundamental solution Gcorresponds to that of the energy density of Eq. (12). According to this equation, the fundamental solution Gis expressed as follows:
Hrepresenting the relation between the virtual source and intensity is obtained from Eq. (14) as follows:
When intensity boundary condition is applied to Eq. (16), the value of
C1 is obtained as follows:
In the left term of Eq. (17), the integral of
is equal to one, because of the property of the Dirac delta function,
δ( r). According to this property and Eq. (16), Eq. (17) is expressed as follows:
In the right term of Eq. (18),
Dremains. The e-2 ηkhcos θ？2cos(2 khcos θ) e- ηkhcos θterm of Drepresents the effect of the reflected wave. If radius approaches zero, this effect diminishes; that is, D=1. Therefore, C1 can be obtained from Eq. (18) as follows:
From Eq. (19), we know that the non-determined constant
Therefore, the fundamental solution representing the relation between the virtual source and energy density or intensity is as follows:
Eq. (20) is the fundamental solutions having the property of a spherical wave and incorporating the effect of the water surface.
BEM which is one of traditional noise analysis methods shows the noise’s directivity pattern at high frequency ranges. This is caused by the difference of wave’s phase. But EFA doesn’t show the noise’s directivity pattern, because EFA doesn’t have the phase information. So representing the noise’s directivity, the directivity factor is developed.
When n sources are existed, the energy is as follows:
is newly added. This equation can be divided to two terms as follows:
From Eqs. (21) and (22), the following equation can be obtained as follows:
is directivity factor.
In the indirect method of the boundary element, the real system of finite size is extended to the infinite field and the virtual source is assumed to exist on the boundary of the real system. And energy and intensity in the concerned field are obtained by using fundamental solution as the sum of the effects due to virtual sources. Based on this concept, in three-dimensional problems, equations of energy density and intensity are obtained in the concerned field as follows:
Vis assumed as the concerned field in three dimensions and Sis assumed as the boundary surface encompassing the concerned field. The fundamental solution, Gand Hhas the same shape as that of Eq. (20).
indicates a field point of the concerned field.
means the position of the virtual source on the boundary.
is the location of input power.
means the virtual source on the boundary surface. The virtual source on the boundary surface,
can be obtained by using Eqs. (24) and (25) according to the property of the boundary condition. If
approaches the boundary surface, the boundary integral of the right term in Eqs. (24) and (25) has a singular integral whose point is centered on the boundary surface. This point is the identity of
Therefore, the integral term of the Eq. (24) is divided into the following terms as follows:
where, if the boundary of the concerned field is a smooth curve surface, the boundary surface
Sεis a hemisphere with the small radius ε, centering on the boundary point. And the boundary surface S？ Sεindicates the part which is subtracted from the boundary surface Sεfrom total boundary surface S. The first right term of the Eq. (26) is as follows:
where, if the radius
εapproaches zero, the area of the boundary surface Sεbecomes 2 πε2. So Eq. (28) is obtained by taking the limit of Eq. (27) as follows:
And if radius
εapproaches zero, the boundary integral term with respect to the boundary surface S？ Sε, is replaced with the boundary integral with respect to the total boundary surface, S. Eq. (27) is expressed about one point on the smooth boundary surface as follows:
EFA is developed by assumptions which are high damping value and medium-to-high frequencies. According to these assumptions, the nearfield evanescent waves can be neglected. So EFA’ s result is not good agreement with a classical solution at low damping values, because the effect of nearfield waves will be affected on the vibration phenomena in low damping system. Therefore EFBEM will be applied to the high damping systems. To verify the accuracy of the developed works, the boundary integral of indirect EFBEM was applied to the radiating noise problems of structures with the vibration of a simple plate model and a simple sphere model in open field and underwater. And the result of the indirect EFBEM was compared with that of SYSNOISE. Figs. 3 and 4 show energy density distribution, respectively, in the z-direction by SYSNOISE and indirect EFBEM when the plate which size is 1
m× 1 m, is located at x-y plane when fis 1 kHz, vis 0.1 m/sand ηis 0. Field points which size 7 m× 3 mis located x-z plane. Though EFBEM can consider the damping loss factor in a medium, a damping loss factor of zero was used so that EFBEM could be compared with SYSNOISE. Fig. 5 shows a simple sphere structure which radius is 0.5 m, with a uniform vibration in open field. When fis 1 kHz, vis 0.1 m/sand ηis 0, Figs. 6 and 7 show the energy density distributions, respectively, in the z-direction by SYSNOISE and indirect EFBEM when the sphere is located under the field points. Field points’ size is 2 m× 2 m. Figs. 8 and 9 show the energy density distribution in the x-z plane predicted by SYSNOISE and indirect EFBEM when the sphere is set at the center of the field points. From Figs. 3,4, and 6-9, we can see that the results of indirect EFBEM agree with those of SYSNOISE in open field. Fig. 10 shows the energy density distribution in the zdirection predicted by SYSNOISE and indirect EFBEM when the sphere is located near the water surface. Though the difference of the results between SYSNOISE and indirect EFBEM is about 1 dB, the patterns of the results are the same. Fig. 11 shows a simple barge type structure with a uniform vibration in open field. When fis 63 Hz, vis 0.1 m/sand ηis 0, Fig. 12 shows the energy density distributions, respectively, in the z-direction by SYSNOISE and indirect EFBEM when the barge type structure is located above the field points. Though the difference of the results between SYSNOISE and indirect EFBEM is about 5 dB, the patterns of the results are the same. And Fig. 13 shows the energy density distribution when fis 1 kHz, vis 0.1 m/sand ηis 0. Because the barge type structure is large, SYSNOISE can’t analyze noise analyses but EFBEM can analyze those over 1 kHz.
This paper analyzed radiating noise problems in the medium-to-high frequency range. The energy governing equation and fundamental solutions having the spherical wave property were developed. And the directivity which was the property of high frequency analysis but couldn’t represent in existing EFA, was developed. Also, the fundamental solutions representing the effect of the water surface were developed for the underwater radiating noise problems. Last indirect EFBEM applying BEM to EFA was developed for this analysis. To verify the developed works, the results of the indirect EFBEM were compared with those of commercial code, SYSNOISE in open field. For plate, the directivity was represented similar shapes in SYSNOISE and developed indirect EFBEM. And the directivity was vanished for sphere model. Also the effect of the water surface was examined by the comparison of the results between indirect EFBEM and SYSNOISE for an underwater case. These comparisons showed satisfactory results and developed indirect EFBEM is applied to barge type structure at high frequency.
[Fig. 1] Wave propagation from medium 1 to medium 2.
[Fig. 2] The distance between field point and source point.
[Fig. 3] The energy density distribution in z-direction when the plate is vibrating: (a) SYSNOISE, (b) Indirect EFBEM.
[Fig. 4] The energy density distribution in z-direction when the plate is vibrating: *, SYSNOISE; o, Indirect EFBEM.
[Fig. 5] Uniformly vibrating 3-dimensional spherical structure.
[Fig. 6] The energy density distribution in z-direction when the sphere is vibrating: (a) SYSNOISE, (b) Indirect EFBEM.
[Fig. 7] The energy density distribution in z-direction when the sphere is vibrating: *, SYSNOISE; o, Indirect EFBEM.
[Fig. 8] The energy density distribution in x-z plane when the sphere is vibrating: (a) SYSNOISE, (b) Indirect EFBEM.
[Fig. 9] The energy density distribution in x-direction when the sphere is vibrating: *, SYSNOISE; o, Indirect EFBEM.
[Fig. 10] The energy density distribution in z-direction when the sphere is vibrating: (a) SYSNOISE, (b) Indirect EFBEM.
[Fig. 11] Uniformly vibrating barge type structure.
[Fig. 12] The energy density distribution in z-direction when barge type structure is vibrating: (a) SYSNOISE, (b) Indirect EFBEM.
[Fig. 13] The energy density distribution in z-direction when barge type structure is vibrating.