Hydrodynamic Behavior of a Submerged Spheroid in Close Proximity to the Sea Surface

: The principal objective of this investigation is to assess the hydrodynamic characteristics and the exciting forces induced by waves acting upon a shallowly submerged spheroid. This study focuses on an arbitrarily shaped spheroid body with a vertical axis, fully immersed beneath the free surface within waters of ﬁnite depth. The methodology outlined here necessitates solving the linear hydrodynamic diﬀraction and radiation problems, which entail discretizing the ﬂow ﬁeld around the body into ring-shaped ﬂuid regions. Within each region, expansions of axisymmetric eigenfunctions of the velocity potential are employed. Complementing the theoretical framework, numerical methodologies are employed utilizing panel models across the we�ed surface of the submerged body. Extensive numerical results concerning the exciting forces induced and the hydrodynamic coeﬃcients are presented in the framework of frequency domain formulations. Through the current analysis, the phenomenon of negative added mass and rapid variations in the added mass and damping coeﬃcients is conﬁrmed, a�ributed to the free surface eﬀect elucidated in terms of the presence of near-resonant standing waves above the submerged body.


Introduction
The examination of hydrodynamic characteristics of bodies submerged in close proximity to the sea surface assumes paramount importance within certain contexts of ocean engineering.Particularly, the maneuverability and control of unmanned underwater vehicles when they interact with the sea surface become pivotal for ensuring safe operations, necessitating precise prognostication of hydrodynamic loading and vessel motion characteristics.However, conventional analyses concentrate on surface vessels and submerged bodies situated at a considerable distance from the free surface.Conversely, limited analysis is available concerning scenarios wherein a vessel or a submerged body maneuvers near the free water surface.
Since the 1960s, Ogilvie [1] conducted a two-dimensional analysis focusing on horizontal cylindrical bodies with circular or square cross-sections submerged near the free water surface, revealing the occurrence of negative hydrodynamic (added) mass phenomena.Subsequently, in [2], experiments involving two-dimensional submerged circular or square cross-sections, subjected to oscillations in the sway and heave directions at various depths beneath the free water surface, were undertaken.The analysis elucidated negative added mass values for square sections when the submergence equaled to one-quarter of the semi-width.Newman et al. [3] further expanded on this topic, identifying negative added mass values, and rapid fluctuations in damping and added mass coefficients within specific frequency ranges, based on resonant free surface motion observed above a submerged body.In addition, McIver and Evans [4] conducted an investigation into the behavior of a submerged vertical cylinder undergoing heave motions.Their study revealed the potential of the added mass to a ain negative values under conditions of shallow submergence, a phenomenon arising when the potential energy of the fluid surpasses its kinetic energy [5].The occurrence of negative added mass values is a ributed to the presence of a near-resonant standing wave above the body at specific frequencies, in the vicinity of which added mass and damping coefficients undergo rapid variations [4,6].Additionally, negative added mass has been observed in the case of floating bodies encompassing a portion of the free surface.Examples in the realm of three-dimensional water wave theory include the floating torus [7], the toroidal body, referred as the "McIver toroid" [8][9][10], as well as non-axisymmetric trapping structures [11].These structures exhibit a trapped-mode frequency wherein no waves are radiated to infinity, corresponding to homogenous solutions of the linearized body-wave interaction problem, and representing unbounded solutions in the frequency domain.Similar physical arguments were employed by [12] to elucidate the negative added mass and sharp force coefficient peaks and hydrodynamic characteristics of groups of interacting axisymmetric submerged bodies near the free surface or seabed.In addition, Shao et al. [13] developed a self-regulating fuzzy depth control method to maintain a cylindrical submerged structure at a specific depth below the free water surface under free surface disturbances.Also, the hydrodynamic coefficients of twin connected circular cylinders were studied in [14].The analyses highlighted the discontinuity of the added mass coefficient at specific wave frequencies, wherein high positive values fell to negative values abruptly.These abrupt variations were a ributed to the resonance frequency of the piston mode and the even sloshing modes of the gap oscillation between the bodies.
The movement of fully submerged bodies near the free surface leads to the deformation of the water surface, resulting in the generation of propagating waves that dissipate energy, giving rise to a resisting force known as the wave drag force.This phenomenon was theoretically investigated by [15][16][17] for spheres and spheroids, while in [18], the influence of an immersed cylinder on wave propagation based on submerged sphere analyses was addressed.Additionally, in [19], fully submerged spheres were utilized to examine the effect of wave drag near the free surface, emphasizing its comparable magnitude to hydrodynamic drag when the top of the sphere was within one sphere radius from the surface.Submerged bodies have also been associated with wave energy conversion, as they surpass surface-piercing converters in enhanced survivability during storms and diminished visual impact [20].McCaulley et al. [21,22] considered shallowly submerged structures with application to ocean wave energy.They verified the negative added mass results and linked resonance of the fluid above the structures to the occurrence of negative added mass coefficients.
While analyses concerning submerged axisymmetric cylindrical/spherical bodies near the free surface are relatively limited, extensive research has been dedicated to submarine maneuverability in near-surface water environments.Submarines frequently descend to periscope or snorkeling depths, for target exploration or ba ery recharging, wherein they encounter environmental loads, including waves, currents, resistance and suction forces [23,24].Numerous studies have investigated the impact of the free water surface on underwater vehicles using both numerical and experimental methods [25][26][27][28][29][30][31][32][33][34][35].Given the imposed forces experienced by submerged structures in different water depths, submarines must maintain specific navigational poses (depth, roll and trim angle), necessitating extensive research on depth control strategies.Specifically, a two-step depth fuzzy controller for the stern angle to counteract the second-order exciting wave forces on submarines was developed by [36], whereas in [37,38], a mathematical model for evaluating wave exciting forces on submerged structures was proposed and depth control simulations using the proportional-integral-derivative method were conducted.In [39], the hydrodynamics of an autonomous underwater vehicle (AUV) during adverse weather conditions were examined.The analysis highlighted negative values of drift forces when the vehicle is close to the free surface.Recently, a novel adaptive trajectory tracking control of an AUV in six degrees of freedom was proposed by [40], and an adaptive control technique utilizing a neural network and a proportional-integral-derivative controller to regulate the depth of a submerged body close to the free water surface was elaborated in [41].The free surface effect on the resistance of an underwater vehicle based on various submergence depths was examined in [42].The analysis revealed that with increasing distance from the free surface, there is a corresponding decrease in the total resistance experienced by the vehicle.Furthermore, scenarios involving ice-covered water surfaces have also been considered [43,44].These studies theoretically and experimentally examined the motions of submerged bodies near the free water surface when the la er was covered by ice, observing significant variations in relative submergence for bodies moving at shallow depths.
The primary objective of this investigation is to assess the hydrodynamic characteristics and wave exciting forces acting on a shallowly submerged spheroid body in waters of finite depth.The present analysis employs a semi-theoretical formulation, integrating solutions to diffraction and motion radiation problems around the floating structure within the framework of linear potential theory, assuming an incompressible, inviscid fluid, and irrotational flow.Additionally, a numerical simulation tool, ANSYS AQWA 2021R1 software, is implemented to compare numerical results with theoretical outcomes.
The subsequent sections of this paper are structured as follows: Section 2 presents the developed theoretical formulation within the domain of linear potential theory, along with a focus on the applied numerical model.Section 3 is dedicated to presenting the outcomes of the two applied formulations, while conclusions are drawn in Section 4.

Velocity Potential
We consider a cylindrical submerged body with a vertical axis of symmetry, which is exposed to the action of regular waves propagating in water depth d with frequency ω.A cylindrical co-ordinate system (r, θ, z) is defined with its origin at the sea bo om and its vertical axis Oz directed upwards (Figure 1).Let the linear translation and rotational vector of the body's motion be denoted by = ( , , ) and = ( , , ) , respectively.Their components correspond to the surge (j = 1), sway (j = 2), heave (j = 3), roll (j = 4), pitch (j = 5) and yaw (j = 6) modes of motion.Under the assumption of a symmetrical mass distribution, a vertical cylindrical body undergoes three-degree of freedom motion in the wave propagation plane under the influence of a regular wave train, i.e., two translations (surge, ; heave, ) and one rotation (pitch, ).
Assuming further that the viscous effects are negligible, the fluid is incompressible, and the displacements of the body and the wave elevation are small, classical linearized water wave theory can be employed.The fluid flow can be expressed by the potential function: where the complex potential function ( , , ) can be articulated, on the basis of linear modeling, as a superposition of the incident, , sca ered, , and radiated wave fields due to the motion of the body, i.e., In Equation ( 2), stands for the velocity potential of the wave field resulting from the forced oscillation of the body in the jth mode of motion with unit velocity amplitude ̇ .
Furthermore, all contributions to the total potential are required to solve Laplace's equation across the entirety of the fluid domain.Additionally, contributions must satisfy the prescribed boundary conditions: = 0, and = , = 1, 3, 5 on The term ( ) ⁄ in Equation ( 5) denotes the derivative in the direction of the outward unit normal vector n to the mean we ed surface S of the body, and are its generalized normal components specified by: = ( , , ), × = ( , , ) Here, represents the position vector of a point on a body's we ed surface with respect to its reference point of motion.The vector is expressed in the co-ordinate system (r, θ, z).Finally, a radiation condition, specifying that propagating disturbances must exhibit outward behavior should be fulfilled (see Equation (13) in Ref. [45]).
The velocity potential of the undisturbed incident wave system, denoted as , which propagates along the positive x-axis can be articulated within the cylindrical coordinate system (r, θ, z) in the following manner: where denotes the wave amplitude, is the Neumann's symbol (i.e., = 1; = 2 for > 0 ), is the th order Bessel function of the first kind, and is the orthonormal function which forms: Here, ( ) denotes the derivative at the free water surface, = .The wave number in Equations ( 7) and ( 8) is linked to the wave frequency by the dispersion equation: = tanh( ).To solve the diffraction and radiation problems, we will employ the method of matched axisymmetric eigenfunction expansions.Specifically, the flow field surrounding the submerged body is divided into coaxial ring-shaped fluid regions, designated by the numerals , , , as illustrated in Figure 1.Within each of these regions, distinct series expansions of the velocity potentials are employed.In alignment with the series representation of the undisturbed incident wave potential, denoted by Equation (7), the diffraction potential within each fluid region ℓ = , , is expressed as: For the radiation velocity potential ℓ at each fluid domain, ℓ, yields: In the unknown functions ℓ (see Equations ( 9) and ( 10)), the initial subscript = , 1, 3, 5 indicates the boundary value problem under consideration, while the subsequent subscript indicates the values of to be accounted for.It is noteworthy that the fluid flow resulting from the imposed oscillation of the cylindrical body in otherwise calm water exhibits symmetry about the = 0-plane and antisymmetry about the = -plane for the surge (j = 1) and pitch (j = 5) mode of motions, whereas it is symmetrical with respect to both planes for the heave mode, (j = 3).Consequently, Equation (10) can be reformulated as follows: Furthermore, it is imperative that both the velocity potential =1, 3, 5, D, and its derivative , = 1, 3, 5, D, maintain continuity at the vertical boundaries of adjacent macro-elements.This requirement yields: In the notation above, the superscripts and denote quantities pertaining to the th and th macro-elements of types and , respectively, with and being the notation for the outmost ring macro-elements of the type and , respectively.Conversely, the superscript is linked to the infinite ring element.
The unknown functions ℓ , ℓ = , , , = 1, 3, 5, represent solutions to the Laplace equation within each distinct fluid domain.These are properly selected to ensure compliance with the kinematic boundary condition at the horizontal and vertical walls of the submerged body, the linearized condition on the free surface, the kinematic condition on the sea bo om, and the radiation condition at infinity.The various potential solutions are then matched by Equations ( 13)-( 20).This formulation yields linear systems of equations, essential for determining the unknown functions ℓ within each fluid domain.The methodology has been extensively detailed in the prior literature, including works by [45] for an arbitrary shaped vertical axisymmetric body, by [46] for a bo omless cylindrical body, by [47] for two coaxial surface piercing truncated cylinders, and by [48] for a permeable vertical cylindrical body, to name a few.Therefore, further elaboration is deemed unnecessary here.Nevertheless, to ensure comprehensiveness, we include appropriate series representations of the functions ℓ , = , 1, 3, 5; ℓ = , in each fluid domain surrounding the submerged body in Appendix A.
As the submerged body is floating near the free water surface, the fluid domain is mainly affected by the submergence of the body.The velocity potential approximation , = 1, 3, 5, for the th ring element equals: In Equation ( 21), equals to = , = 1,3, , = 5 , whereas the function equals to: The functions , * , are expressed as: Here, , are the th order modified Bessel function of the first and second kind, respectively.Also, a stand for the roots of the transcendental equation: Equation ( 24) has two imaginary and an infinite number of real roots.Assuming a = − as the imaginary root, the well-known dispersion equation is obtained.
In Equation ( 21), are orthonormal functions for ∈ [ − , ] obtained by: Finally, the coefficients , , , * in Equation ( 21) denote the unknown Fourier coefficients of the -th ring element of type .
Considering the extreme positions of the submerged body near the free surface, i.e., → 0, special a ention must be paid as the respective arguments of the Bessel functions involved in the series representation (see Equation ( 21)) become too large.To circumvent the difficulty, asymptotic expressions for the Bessel functions for large arguments are introduced [49] and the corresponding expression for the velocity potential representation in the second fluid domain, II, is recast as follows: where: Here, = 4 ; the rest of the symbols and functions involved in Equation ( 26) have been defined previously.

Hydrodynamic Forces
Upon establishing the diffraction and radiation velocity potentials surrounding the submerged body, the exciting wave forces and hydrodynamic reaction forces exerted on the body can be derived through the following equations.
The horizontal and vertical exciting forces are determined by: Here, stands for the generalized normal components, which are defined by = ( , , ) and × = ( , , ), with the position vector of a point on the body's we ed surface .
The horizontal moments are given by the sum of components, i.e., ( ) = ( ) + ( ), where: Similarly, the hydrodynamic reaction forces and moments can be wri en as: Here is defined by Equation ( 6). Thus: and ( ) = ( ) + ( ) given by: Following [50] the hydrodynamic (added) mass, , , and damping coefficients, , , of the submerged body in the th direction due to forced oscillation in the th direction can be wri en as: Consequently, from Equations ( 33) and ( 30), it can be derived that:

Validation of the Results
Initially, the theoretical outcomes obtained from the aforementioned analysis are compared with those existing in the literature.The present analytical model is applied to the case of an oblate spheroid at finite water depth in order to compare the results with those of [51].In addition, the utilization of commercial software ANSYS AQWA 2021R1 is incorporated.ANSYS AQWA, a component of the ANSYS Mechanical Enterprise suite, facilitates diffraction and radiation analyses, based on potential theory [52].The physics of AQWA are applicable for finite depth waters and are solved within the frequency domain framework.The computational demands, expressed in terms of CPU time, for conducting the numerical simulations employing 5612 we ed elements, amount to approximately twenty seconds for each wave frequency.With regard to the precision of the theoretical modeling, its integrity is influenced by the method employed for evaluating the Fourier coefficients within each fluid domain surrounding the body.In the current computations, 80 terms are employed for the series expansions of the velocity potential in the outer and upper fluid domains, while 150 terms are preserved for the velocity representation in the lower domains.The theoretical outcomes are generated utilizing proprietary inhouse computational software [53], developed in the Fortran programming language, with the CPU time constraints set to less than one second per each analyzed wave frequency.
The examined oblate spheroid of a semi-axis α is subjected to an incident wave, at a water depth 20α.The body is fully submerged below the free water surface, and the distance between the spheroid's center and the free surface is 1.5α.Also, the distance from the body's center to the pole along its symmetry axis is α/1.25 (see Figure 2).The comparison is made in terms of the dimensionless exciting forces and moments (see Equations ( 28) and ( 29)), i.e.,: /( /2), /( /2), and /( /2), where is the water density, is the gravity acceleration, and /2 is the wave amplitude.The results are plo ed against the dimensionless value of kα, where k is the wave number.Figure 3 illustrates the horizontal and vertical exciting forces and horizontal moments on the oblate spheroid.An excellent correlation between the theoretical results (present formulation), the numerical outcomes (ANSYS AQWA), and those reported in [51] can be observed.Consequently, it may be inferred that the current theoretical model aptly simulates the hydrodynamic behavior of a submerged spheroid body under wave interactions.

Oblate Spheroid
Subsequently, an investigation is conducted on an oblate and a prolate spheroid under various levels of submergence.Initially, an oblate spheroid of a semi-axis , and distance from the body's center to pole along its symmetry axis = /4, is considered.Here, the water depth equals to = 20 .The body is assumed to be floating at ℎ = /10, /5, /2 and below the free water surface (see Figure 4).In Appendix B, the dimensions of the coaxial ring elements for the oblate spheroid utilized in the theoretical formulation are presented.Figures 5-7 present the hydrodynamic forces on the oblate spheroid evaluated by the aforementioned theoretical formulation, under various submergence values versus kα.Also, the AQWA software is applied for comparative purposes for the shallowest investigated submersion, i.e., ℎ = /10.Specifically, Figure 5 depicts the exciting forces and moments on the body for various considered submergences.Here, the results are nondimensionalized by the term ( /2) for the forces, and ( /2) for the moments.It is evident that the employed methodology involving the division of the flow field surrounding the submerged body into coaxial ring-shaped fluid regions yields precise results in terms of the exciting forces and moments on the body.The theoretical predictions align with the numerical findings from AQWA (see Figure 5a,c,e).As far as the influence of the submergence value is concerned, it is apparent that both the forces and moments exhibit a progressive increase up to a maximum occurring at a relatively small wave frequency, followed by a subsequent decrease to zero.For greater submergence depths, the maxima of the surge and heave exciting forces shift towards lower frequencies compared to the ones for shallower submergence values.Furthermore, the submergence parameter notably impacts the frequencies of the minimum loading components, as evidenced in Figure 5b,d,f.As the submergence value increases, the zeroing of the forces and moments is shifted towards higher wave frequencies.It should be noted that the moments delineated herein and in the subsequent section are assessed with respect to the water surface, whereas for the comparisons with AQWA, they have been assessed in relation to the body's center of gravity.Concerning the magnitudes of the exciting forces and moments, it is demonstrated that a reduction in the distance between the free surface and the body leads to an increase in the exciting wave loads.Figure 6 examines the hydrodynamic added mass coefficients of the immersed oblate spheroid in the surge, heave, pitch and surge-pitch directions.Here the results are normalized by the terms ( ) for the surge and heave, ( ) for the pitch, and ( ) for the surge-pitch coeffients, respectively.Also, denotes the submerged volume of the body.It is demonstrated that the hydrodynamic added mass coefficients can also be accurately evaluated utilizing the present theoretical formulation.The few a ained discrepancies at /( ) between the theoretical and numerical outcomes (see Figure 6a), in the viscinity of (1.5,4), can be considered negligible.Conversely, it is immediately apparent that the comparisons of the heave, pitch and surge-pitch added mass coefficients, between the present method and AQWA, are quite favorable (Figure 6c,e,g).Further, the impact of the submergence value on the body's hydrodynamic characteristics is presented in Figure 6b,d,f,h.Clearly, as the distance between the body and the free water surface decreases, the hydrodynamic coefficients a ain a more tense oscillatory behavior.Specifically, for ℎ = /10, the surge, heave and pitch added mass, there is a sharp increase in their values, followed by an abrupt decrease to negative values.However, this is not the case for higher submergencies, i.e., ℎ = /5, /2 and , where /( ), /( ) and /( ) generally a ain positive values.In contrast, the surgepitch added mass initially a ains negative values, whereas as increases, and the submergence decreases, /( ) presents positive values.Furthermore, it can be noted that there is a smooth variation pa ern of the body's hydrodynamic coefficients for large submergence, i.e., ℎ = .Figure 7 presents the hydrodynamic damping coefficients of the immersed oblate spheroid in the surge, heave, pitch and surge-pitch directions.Here, the results are normalized by the terms ( ) for the surge and heave, ( ) for the pitch, and ( ) for the surge-pitch coefficients, respectively, where is the wave frequency.An excellent correlation between the theoretical and numerical results is evident (see Figure 7a,c,e,g).Consequently, the accuracy of the formulation remains unaltered.As far as the effect of the submergence values on the body's hydrodynamic damping coefficients are concerned, it is apparent that the ℎ value impacts the variation pa ern of these coefficients.Specifically, as ℎ decreases, the values of the damping term increase.In addition, a sharp increase in , = 1, 3, 5, followed by an intense decrease to zero is observed as the submergence value decreases.Nonetheless, this oscillatory behavior diminishes, and the wave frequency at which , =1, 3, 5, zero, shifts towards higher ω with decreasing ℎ values.On the contrary, the surge-pitch damping coefficient exhibits oscillatory behavior in the vicinity of negative values.These behaviors tend to subside as the submergence values increase, leading to a more gradual fluctuation pa ern.

Prolate Spheroid
The methodology outlined herein can be effectively employed for a prolate spheroid under various levels of submergence.Specifically, a prolate spheroid of a semi-axis , and distance from the body's center to the pole along its symmetry axis = 2 , is considered.The water depth equals to = 40 , and the body is assumed to be floating at ℎ = 2 /10, 2 /5, , 2 below the free water surface (see Figure 8).In Appendix B, the dimensions of the coaxial ring elements for the prolate spheroid utilized in the theoretical formulation are presented.Figures 9-11 present the hydrodynamic forces on the prolate spheroid evaluated by the presented theoretical formulation, under various submergence values versus kα.Also, the AQWA software is applied herein for comparative purposes for the shallowest investigated submersion, i.e., ℎ = 2 /10. Figure 9 depicts the dimensionless exciting forces and moments acting on the body.The normalized factor has been defined in Section 3.2.It is evident that the theoretical formulation adeptly captures the diffraction hydrodynamics of the prolate spheroid, as evidenced by the favorable comparisons with the AQWA outcomes (see Figure 9a,c,e).However, there are some discrepancies between the theoretical and numerical outcomes regarding the pitch moment at > 3.2 (see Figure 9e).Notably, this wave frequency surpasses the typical frequency bandwidth for gravity waves.Concerning the influence of the submergence values on the exciting forces and moments, a progressive increase in forces is observed, accompanied by a corresponding smooth decrease irrespective of the ℎ value.Nevertheless, the upper bound of the forces a ains higher values as ℎ decreases.Also, the maximum of the exciting forces shifts to a smaller frequency as ℎ increases.Noteworthy is the absence of a similar oscillatory pa ern in the exciting forces as observed in an oblate spheroid (see Section 3.2).As far as the pitch exciting moment is concerned, Figure 9f  Figure 10 illustrates the hydrodynamic added mass coefficients of the immersed prolate spheroid in the surge, heave, pitch and surge-pitch directions.The results presented here are normalized using the same factor as detailed in Section 3.2.It is evident that while the theoretical formulation accurately predicts the surge hydrodynamic added mass coefficients (see Figure 10a), notable discrepancies arise when compared to AQWA results for the heave, pitch and surge-pitch added mass (Figure 10c,e,g).These discrepancies are more pronounced at higher values of .Such deviations between the theoretical and the numerical results can be a ributed to the number of the considered fluid domains surrounding the body, which in this study amounts to 18 (i.e., 9 above the semi-axis α, and 9 below α).Regarding the effect of the submergence values on the hydrodynamic added mass of the spheroid, it is apparent that larger immersion depths result in a smoother variation in the added mass for both the surge and heave added mass.On the other hand, the pitch and surge-pitch added mass exhibit smooth fluctuations regardless of the ℎ value.It is worthwhile to note that although the spheroid is submerged near the free surface, the surge, heave and pitch added mass consistently remain positive, without reaching negative values.
Figure 11 illustrates the hydrodynamic damping coefficients of the immersed prolate spheroid in the surge, heave, pitch and surge-pitch directions.These results are normalized by the factor previously mentioned (refer to Section 3.2).The theoretical formulation aligns well with the numerical outcomes, particularly noticeable in Figure 11a,c.However, some discrepancies are notable concerning the pitch and surge-pitch damping coefficients, which are more pronounced in the vicinity of ∈ (1,2).Nevertheless, these can be considered negligible given that both methodologies exhibit a similar trend and variation pa ern.As far as the influence of the submergence values on the body's hydrodynamic damping coefficients are concerned, it is apparent that the ℎ value impacts the variation pa ern of these coefficients.Larger immersion depths result in a smoother variation in the damping coefficients.Also, it can be observed that as ℎ increases, the wave frequencies at which the damping coefficients reach their maximum values shift to lower frequencies.

Vertical Cylinder
Based on the analysis, it has been concluded that negative added mass manifests in a submerged oblate spheroid body when floating near the free surface.Conversely, this phenomenon does not occur with a prolate spheroid of equivalent mass.Consequently, an investigation is undertaken to examine whether negative added mass is similarly observed in a fully submerged vertical cylinder, possessing the same mass as the aforementioned spheroids, when situated near the free surface.The objective is to examine the validity of this phenomenon in the context of the considered vertical cylinder characterized by a radius α, and a height = 4.89 .It is assumed that the water depth is 48.9 and the distance between the free surface and the body's upper surface is ℎ = α/4.09(see Figure 12).Figure 13 depicts the hydrodynamic coefficients of the cylinder, with Figure 13a,b presenting the surge, heave, pitch and surge-pitch hydrodynamic added mass and damping coefficients as functions of , respectively.These results are normalized by the factors previously mentioned (refer to Section 3.2).It is evident from Figure 13 that despite the vertical cylinder possessing a mass equivalent to that of the two examine spheroids, and its height matching the vertical axis of the prolate spheroid, the heave added mass exhibits negative values in the vicinity of ∈ (1,2).Conversely, the prolate spheroid demonstrates positive added mass values.Furthermore, the surge, pitch and surge-pitch added masses a ain a similar oscillation pa ern, consistently obtaining positive values across the entire spectrum (with the exception of the surge-pitch added mass, which a ains negative values within the examined ).As for the surge, heave and pitch damping coefficients, these are always positive as in the spheroid cases, while the surge-pitch damping coefficients assume negative values within the examined range.

Sphere
In the subsequent analysis, a ention is directed towards the particular case of a spheroid characterized by equal semi-axes, effectively constituting a sphere.This sphere under examination shares the identical mass as the previously mentioned spheroids and possesses a radius denoted as α.The water depth equals to = 31.75and the sphere is assumed to be floating at ℎ = /6.3,/3.15, /1.26, /0.63, as illustrated in Figure 14.In Appendix B, the dimensions of the coaxial ring elements for the sphere utilized to the theoretical formulation are presented.Given that the exciting forces and moments on a submerged sphere have been presented in [54], and the theoretical approach has been validated against numerical formulations and high fidelity CFD simulations with very good correlation, further elaboration is omi ed herein.Figures 15 and 16 depict the hydrodynamic coefficients of the sphere (added mass and damping coefficients) normalized by the same parameters outlined in Section 3.2 for varying degrees of submersion.Here, the outcomes are derived by the theoretical formulation, whereas AQWA software is utilized for comparative purposes for a submergence ℎ = /6.3.From Figure 15, it is evident that the added mass in both the surge and heave directions a ains positive values for ℎ = /6.3, as well as for ℎ = /3.15,/1.26, /0.63.This trend was consistent for the prolate spheroid (refer to Section 3.3), whereas the oblate spheroid exhibited negative added mass values for the lowest considered submergence (refer to Section 3.2).Additionally, the vertical cylinder only demonstrated a negative added mass in the heave direction for the lowest submergence value (refer to Section 3.3).Regarding the comparisons between the theoretical and the numerical results, a very good agreement is evident.However, some discrepancies, particularly regarding the heave added mass coefficients (see Figure 16c), are observed, similar to those found for the prolate spheroid.As far as the influence of the submergence depth on the added mass is concerned, it can be seen that as the immerse depth decreases, , exhibit an oscillatory variation pa ern.
Figure 16 illustrates the damping coefficients of the sphere in the surge and heave directions.A similar variation pa ern to the aforementioned spheroids is observed.Specifically, it is evident that greater immersion depths result in a smoother variation in the damping coefficients, whereas as ℎ increases, the wave frequencies at which the damping coefficients reach their maximum values shift to lower frequencies (see Figure 16b,d).Regarding the accuracy of the theoretical formulation, a very good correlation between theoretical and numerical outcomes can be observed.Consequently, the theoretical formulation can be confidently extended to a submerged sphere.
In regard to the number of region discretizations, the present work selects eight discretizations for the sphere for each of the and fluid domains.Figure 17 depicts the hydrodynamic coefficients for a submergence depth of ℎ = /6.3, under four different values of region discretization for the upper and lower fluid domains, namely, 5, 8, 10 and 11.It can be observed that the number of region discretizations has a minor effect on the surge hydrodynamic coefficients.However, it affects the damping coefficients, particularly near peak frequencies, as well as the heave added mass across the entire wave frequency band.The la er appears to be principally affected by the number of region discretizations.Nevertheless, as the number of fluid regions increases, the values of the heave added mass converge.Figure 18 presents the surge and heave added mass coefficients for the three examined bodies (i.e., oblate sphereoid, prolate spheroid, and sphere) for the same submergence depth of the three bodies consolidated into one diagram.The results are presented non -dimensionally as a function of the body's volume, which remains constant, and are plo ed as a function of kα.It is observed that the surge added mass of the prolate spheroid a ains higher values compared to the other geometries examined.Regarding the heave added mass, the coefficient exhibit a more pronounced oscillatory behavior as the maximun radius of the body increases (this is the case of the oblate spheroid compared to its sphere and prolate spheroid equal-volumed counterparts).

Conclusions
In this work, an examination is conducted into the influence of submergence on the hydrodynamics of a spheroid body possessing constant mass, namely, oblate, prolate and spherical shapes, while floating beneath the free water surface.The investigation employs a theoretical formulation, with its accuracy assessed against numerical results obtained from ANSYS AQWA software.The principal conclusions derived from this study are outlined as follows:


The hydrodynamics of submerged bodies undergo significant alteration depending on their depth below the free water surface.Notably, when the distance between the body and the free water surface is minimal, pronounced peaks in the body's hydrodynamics are observed.


Negative added mass coefficients manifest at specific wave frequencies and for shallow submergence depths.However, this trend diminishes as the distance between the body and the free surface increases.Nonetheless, damping coefficients consistently exhibit positive values across the examined frequency range, regardless of the submergence depth. Oblate spheroids exhibit negative added mass coefficients in the surge, heave and pitch directions at the lowest examined submergence depth, while prolate spheroids demonstrate positive added mass coefficients in these directions.Similarly, the examined sphere exhibits positive added mass coefficients in the surge and heave directions.Conversely, cylinders only exhibit negative added mass coefficients in the heave direction.


With respect to the methodologies employed, the results derived from the developed theoretical formulation closely align with the variation pa ern observed from the numerical method outcomes.However, deviations are evident at wave frequencies corresponding to peaks in the body's hydrodynamics.Furthermore, it is pertinent to note that the present formulation is not applicable to ellipsoid-shaped bodies.


The analysis findings indicate that in near-surface operations (i.e., for submergence depths less than the radius of the submerged body), wave motion primarily comprises a standing wave, resulting in rapid fluctuations in both added mass and damping coefficients.Hence, when designing submerged bodies, meticulous consideration of the geometric a ributes in relation to the wave characteristics is imperative, as unexpected forces on the hull may arise, augmenting the body's inertia and influencing its acceleration.The functions are derived from Equation (A5).In Equations (A4) and (A8), is the Newmann's symbol, whereas , , , * are the unknown Fourier coefficients to be determined through the solution process.

Appendix B
In Tables A1-A3, the dimensions of the coaxial ring elements for the fully submerged bodies are presented.

Figure 1 .
Figure 1.Discretization of the flow field around the considered submerged body with vertical axis of symmetry.

Figure 2 .
Figure 2. 2D representation of the examined oblate spheroid (left image); 3D image of an axisymmetric oblate spheroid from AQWA (center image); 3D image of an axisymmetric oblate spheroid using the discretization of the flow field into coaxial ring-shaped fluid regions (right image).

Figure 3 .
Figure 3. Dimensionless exciting forces and moments on an oblate spheroid calculated by theoretical and numerical formulations.The results are also compared with the outcomes of [51].

Figure 4 .
Figure 4. 3D representation (left figure) and plane view (right figure) of the examined oblate spheroid.

Figure 5 .
Figure 5. Dimensionless exciting forces and moments on an oblate spheroid calculated by theoretical and numerical formulations: (a,c,e) Dimensionless surge, heave force and pitch moment at the shallowest investigated submersion, ℎ = /10, respectively.The results are compared with the outcomes from AQWA; (b,d,f) Dimensionless surge, heave force and pitch moment at various levels of submergence ℎ = /10, /5, /2, .

Figure 8 .
Figure 8. 3D representation (left figure) and plane view (right figure) of the examined prolate spheroid.

Figure 9 .
present the hydrodynamic forces on the prolate spheroid evaluated by the presented theoretical formulation, under various submergence values versus kα.Also, the AQWA software is applied herein for comparative purposes for the shallowest investigated submersion, i.e., ℎ = 2 /10.Figure9depicts the dimensionless exciting forces and moments acting on the body.The normalized factor has been defined in Section 3.2.It is evident that the theoretical formulation adeptly captures the diffraction hydrodynamics of the prolate spheroid, as evidenced by the favorable comparisons with the AQWA outcomes (see Figure9a,c,e).However, there are some discrepancies between the theoretical and numerical outcomes regarding the pitch moment at > 3.2 (see Figure9e).Notably, this wave frequency surpasses the typical frequency bandwidth for gravity waves.Concerning the influence of the submergence values on the exciting forces and moments, a progressive increase in forces is observed, accompanied by a corresponding smooth decrease irrespective of the ℎ value.Nevertheless, the upper bound of the forces a ains higher values as ℎ decreases.Also, the maximum of the exciting forces shifts to a smaller frequency as ℎ increases.Noteworthy is the absence of a similar oscillatory pa ern in the exciting forces as observed in an oblate spheroid (see Section 3.2).As far as the pitch exciting moment is concerned, Figure9freveals a similar fluctuation pa ern across varying ℎ values.Nevertheless, the wave frequency maximizing the acting moment shifts to lower values as ℎ increases.

Figure 12 .Figure 13 .
Figure 12. 3D representation (left figure) and plane view (right figure) of the examined vertical cylinder.

Figure 17 .
Figure 17.Dimensionless hydrodynamic added mass and damping coefficients of a sphere calculated for different values of region discretization: (a,b) Surge, , and heave, , added mass coefficients; (c,d) Surge and heave , , damping coefficients.