Research Article  Open Access
Further Application of Surface Capturing Method and Cartesian Cut Cell Mesh on Hydroelastic WaterEntry Problems of FreeFalling Elastic Wedge
Abstract
In order to study the interactions between fluid and elastic structure (such as marine lifeboat falling down and ship), this paper presents a new CFD method on hydroelastic waterentry problem of freefalling elastic wedge, which can more conveniently handle moving solid boundaries. In the CFD solver, a surface capturing method and the Cartesian cut cell mesh are employed to deal with the moving free surface and solid boundaries, respectively. On the other hand, in structural analysis, the finite element method and lathbeam structural model are introduced to calculate the elastic response. Furthermore, based on the current CFD and structural solver, a particular data transfer method and coupling strategy are presented for the fluidstructure interaction. Finally, by comparing numerical results with experimental data, the present method is validated to be available and feasible for hydroelastic waterentry problem and further successfully adopted to analyze the motion characteristics of freefalling elastic wedge.
1. Introduction
Interactions between fluid and elastic structure, especially water entry problem of elastic body, are of concern and practical importance in many ocean and marine engineering applications, such as seaplane landing on sea, ship slamming, sloshing, greenwater impacting on deck, and marine lifeboat falling down. For the water entry problem, it is complicated and difficult to numerically study moving solid bodies and free surface and the challenge is further higher if global motion and elastic deformation of solid body need to be calculated from the fluidstructure interaction.
The water entry problem has been widely studied by many scientists and engineers using various methods since the initial study of von Karman [1], such as Wagner [2], Yamamoto et al. [3], Worthington [4], Bottomley [5], Mayo [6], Stoffmacher [7], Chuang [8], Huges [9], Ochi and Motter [10], Belyschko and Mullen [11], Korobkin and Pukhnachov [12–15], Greenhow and Yanbao [16], Cointe and Armand [17], Wilson [18], Miloh [19], Zhao and Faltinsen [20, 21], Takagi [22], Fraenkel and McLeod [23], Wu et al. [24–26], Greenhow and Moyo [27], Sames et al. [28], Mei et al. [29], Lu et al. [30], Scolan and Korobkin [31], Aquelet and Souli [32–34], Battistin and Iafrati [35], Kleefsman et al. [36], Sun [37], Faltinsen and Semenov [38], Tveitnes et al. [39], Truscott and Techet [40], and Xu et al. [41].
Therein boundary element method is based on the theory assumption of potential flow, so the rotational motion near free surface cannot be well simulated. Thus in order to better predict the force and motion of body and achieve fluid information near free surface simultaneously, a CFD method should be applied in water entry problem. Furthermore, compared with the traditional CFD method ([36, 42–44] and so on), the surface capturing method [45–48] based on Cartesian cut cell mesh [47–52] as a novel method can not only well simulate complicated cases of wave breaking and trapping, but also has higher computational efficiency to treat moving boundaries by updating a few cut cells locally rather than remeshing the whole flow domain. On the other hand, so far there are only few documents which consider the interaction among fluid field, global motion, and local deformation of elastic body simultaneously for the water entry problem. Therefore, a new CFD method of surface capturing method and Cartesian cut cell mesh should be taken and developed to handle hydroelastic waterentry problem of freefalling wedge.
This paper presents surface capturing method and Cartesian cut cell mesh to treat moving free surface and solid boundaries. In the CFD solver, incompressible Euler equations are presented as governing equations for a variable density fluid and the location of free surface can be captured as a contact discontinuity in the density field. Furthermore, finite volume method is applied in numerical spatial discretization. Therein Roe’s approximate Riemann solver is adopted to evaluate numerical flow flux and dualtime stepping technique with artificial compressibility method is used for time advancing. On the other hand, in the structural solver finite element method and lathbeam structural model are used to analyze structural response of elastic wedge. Particularly based on the geometry of Cartesian cut cell mesh and lathbeam element here a particular data transfer method on solid boundary is deduced and presented. Furthermore, the flow field, structural response, and global motion are calculated in a coupled manner, which not only computes hydrodynamic force on solid boundary, but also allows the feedback of elastic deformation and global motion into the CFD solver. Finally, some test cases of water entry for various kinds of freefalling elastic wedge are numerically simulated and the calculated results show the feasibility and availability of the present method. Furthermore, by this method the global motion and local response of elastic wedge are studied and analyzed during waterentry phase.
2. Description of CFD Solver
2.1. Governing Equation and Boundary Condition
For the 2D, incompressible, unsteady, inviscid fluid system with a variable density field, the  and coordinate axis are along and vertical to free surface, respectively; thus the governing equations (mass conservation, , directional momentum, and incompressibility constraint equations) for surface capturing method can be written in conservation form, which have been verified and described by Kelecy and Pletcher [45], Pan and Chang [46], Qian et al. [47], and W. Wang and Y. Wang [48]: where, is the fluid density, and represent the and directional fluid velocities, is the fluid pressure, and are the and directional accelerations of body force, and, if assume the only body force is gravity, the source term and , respectively.
In the present study, various kinds of boundary conditions can be classified.(1)Outlet or open boundary: a zero gradient condition is applied in the velocity and density and the pressure at this boundary is fixed to be static pressure, which allows fluid to enter or leave the computational domain freely according to the local flow velocity and direction.(2)Solid body boundary: the nopenetration condition can be applied in velocity and the density is assumed to have a zero normal gradient. For pressure boundary condition, here the momentum equations (2) and (3) are projected on the normal direction of solid boundary and can be rewritten as where is fluid velocity at the solid boundary, is normal unit vector of solid edge, and is normal projection of average velocity of solid edge. Furthermore, a detailed description about how to achieve and based on the fluidstructure interaction on solid boundary should be demonstrated in Section 4.1.2 of data transfer method.
2.2. Numerical Solver
In the present study, Cartesian cut cell mesh is employed for spatial discretization, which can be generated by cutting solid bodies out of a background Cartesian mesh. Thus fluid cell, solid cell and cut cell are accordingly created, which include various storage data (such as identification of cell type, geometry, and fluid variables). For more information, refer to Coirier and Powell [49], Yang et al. [50, 51], Causon et al. [52], and Qian et al. [47].
Based on the Cartesian cut cell mesh, here a cellcentral finite volume method is applied in numerical discretization. In this numerical scheme, Roe’s approximate Riemann solver is adopted to calculate the numerical flux on each edge of fluid cell where fluid variables are reconstructed by using a piecewise linear upwind scheme. In the linear representation, a leastsquare method [48] is taken to achieve variable gradients and Superbee limiter [47] is used to maintain monotony of numerical scheme and control spurious oscillations. Then the flux on solid boundary can be achieved by means of exact Riemann discontinuous solution [48]. Furthermore, dualtime stepping technique with artificial compressibility method is used for time advancing. Finally, an approximate LU factorization (ALU) scheme [53] is applied to solve linear equations. The details about numerical scheme can be further found from the study of W. Wang and Y. Wang [48].
3. Introduction of LathBeam Structural Model
For the bottom plate with cylindrical bending of elastic wedge structure, a lathbeam model can be used to analyze structural deformation. Compared with EulerBernoulli beam, due to bilateral constraints, the crosssection of lathbeam model cannot freely deform and keep the original shape. On this basis, only if by introducing ( is the material elastic modulus and is Poisson ratio) as new elastic modulus into EulerBernoulli beam theory can the elastic deformation of lathbeam model be studied. For more information, please refer to Shu and Tan [54]. Thus, in structural solver, here finite element method and lathbeam model based on EulerBernoulli beam theory should be applied to calculate structural response on solid boundaries of freefalling elastic wedge.
Initially, the global (fixed) and local (moving with body) coordinate systems are created. Then, lathbeam structure is divided into many beam elements with equal length, which can support bending and tensioncompression effects simultaneously. Here in the local coordinate system, the nodal displacement vector of element can be demonstrated as where , , , and are translational displacements, and and are rotation angles on two nodes.
Next, in local coordinate system the corresponding mass and stiffness matrix of lathbeam elements are generated and then load vectors of element nodes are calculated from the fluid pressure on solid boundary. Subsequently, the mass matrix, stiffness matrix, and load vector are transformed into global coordinate system. For the waterentry model of elastic wedge, the elastic deformations basically take place on sloping sides; therefore here on the corner of sloping and straight sides the rigid connection is selected as boundary constraint of lathbeam structure. Then the global dynamical equation without structural damping can be written as where is the nodal load vector and the detail calculation method can be referred to in Section 4.1.1. In matrix and , for the structural element , the and can be expressed as where, for lathbeam structural model, is the structural density, is the element thickness, is the element length, is the crosssection area, is the moment of inertia, , is the Poisson ratio, and is the material elastic modulus.
Finally, by means of numerical Newmark method, the nodal vector of displacement, velocity, and acceleration are computed and fed back to CFD solver and meanwhile local structural response on elastic boundary can be obtained.
4. Analysis of FluidStructure Interaction
For the waterentry model of freefalling elastic wedge, the global motion and local deformation of elastic body should be calculated from the fluidstructure interaction, so here the data transfer method on elastic boundary and fluidstructure coupling strategy are discussed, respectively.
4.1. Data Transfer Method on Elastic Boundary
In order to study and analyze the fluidstructure interaction, it is necessary to transfer the fluid pressure to the element node and then feed the global motion and structural deformation back to CFD solver.
In this paper, Cartesian cut cell mesh in CFD solver and lathbeam element in structural solver are different from each other. So here a particular data transfer method on solid boundary is deduced and demonstrated, which includes transfer of distributional fluid pressure to structural nodes and feedback of structural response on CFD boundary condition.
4.1.1. Transfer of Distributional Fluid Pressure to Structural Nodes
Based on respective geometries of Cartesian cut cell mesh and lathbeam element, there should be many kinds of position relations between the solid boundary in cut cell and structural element to be considered (three classic cases as shown in Figure 1), which are full inclusion of in (Figure 1(a)) or in (Figure 1(b)), partial intersection between with (Figure 1(c)), and some other situations (such as that coincides with ). Nevertheless, the principle of transferring the pressure distribution to the FEM nodes is the same in all cases. Thus, the important issue is how to estimate the pressure distribution along the beam elements. Here take the structural element in Figure 1(a), for example, to demonstrate how to calculate the nodal forces and of the nodes and from fluid pressure distributed on .
(a)
(b)
(c)
Firstly, local coordinate system should be created on structural . Therein coordinateorigin coincides with node , coordinates axis is clockwise and, along with element , axis points from fluid to structure. Thus unit vectors of and can be written: where and are coordinates of nodes and , respectively.
Next, by using the coordinates of centre in cut cell, position vector can be achieved:
Then by combining (9) with (11), position vector (from centre to arbitrary point on ) can be written: where is the coordinate value along coordinates axis .
Finally, according to the equivalence principle about total force and moment, the distributional fluid pressure on can be transferred to the equivalent concentrated forces and on nodes and of structural lathbeam element , which can be described as follows: where and are the fluid pressure and its gradient stored in centre of cut cell, which can be obtained by CFD solver. Similarly, the equivalent concentrated forces on element nodes can be calculated for the other situations of position relations.
4.1.2. Feedback of Structural Response on CFD Boundary Condition
Here among many kinds of position relations of Cartesian cut cell mesh and lathbeam element as mentioned in Section 4.1.1 (Figure 1), take one (as shown in Figure 1(a)), for example, to describe how to adopt structural response to update boundary condition of CFD solver.
In Figure 1(a), structural element fully includes the solid boundary . Here the velocities of nodes and in local coordinate system have been calculated by finite element analysis.
Firstly based on the shape function of structural element , the velocities of point in local coordinate system can be obtained: where and are the shape function versus bending and tensioncompression effects, respectively. For the lathbeam element, , , , , , and , and are natural coordinates of structural element . Then by the same method, the velocities and of point also can be calculated.
Here it should be noted that and terms in (5) are the sum velocity, which includes local elastic deformation and global rigid translation. Subsequently, the average normal velocity and tangential velocity of solid boundary can be calculated: where and are normal and tangential unit vectors of solid boundary and is the global rigid velocity of solid boundary, which can be calculated by a subiteration coupling approach in Section 4.2.
Meanwhile, gradient of normal velocity on solid boundary can be achieved:
Finally, the tangential velocity, average normal velocity and its gradient can be obtained from (15)~(18), and then should be applied in the new boundary condition (5) of cut cell.
4.2. FluidStructure Coupled Solution Strategy
In the present study, the fluidstructure interaction of waterentry model for freefalling elastic wedge includes the coupling calculation of fluid with global body motion and local structural deformation.
On the one hand, in order to ensure enough numerical stability, a subiteration coupling approach for the global velocity of wedge proposed by Kleefsman et al. [36] is taken. Here assume that the elastic wedge moves with one degree (up and down) with an acceleration; the direction velocity is zero and the directional global rigid velocity of solid boundary can be computed by where is the physical time step, is the wedge mass, and is the directional global rigid velocity of solid edge at the physical time . is a relaxation parameter to control the stability of subiteration, which is in the range from 0 to 1. The larger the is, the faster the convergence rate of subiteration is. Equation (19) demonstrates a subiterative process of fluidglobal body motion coupling at each physical moment. At physical time step , the iterative initial value is taken as . For the subiterative step, based on the vertical body velocity , a new pressure field is calculated and then is obtained by direct integration of the new pressure over body surface. When , the subiteration is convergence and is taken as the new velocity boundary condition of the solid body for the next physical time step .
On the other hand, compared with military missile at high speed, ship and ocean engineering structures usually enter water with mediumlow velocity and therefore the deformation of elastic boundary is relatively small and its reaction to the fluid field can be neglected. Thus in order to increase the computational efficiency, here a twoway “weak” coupling method as shown in Figure 2 is used for the coupled solution of fluid and local structural deformation.
5. Results of Test Cases and Discussions
5.1. Validation
In this section, the current numerical method for hydroelastic analysis on waterentry problems of freefalling wedge should be tested by comparing calculated results with experimental data of Sun et al. [55] about three different waterentry test cases.
Based on the parameters of water tank in Sun’s experiment, 2D CFD computational domain was a square of 0.8 m × 1.45 m and the water depth was 1.1 m. The material properties of elastic wedges were 1.059 × 10^{3} kg/m^{3} density, 2.67 × 10^{9} kg/(ms^{2}) elastic module, and 0.357 passion ratio. Furthermore, geometric parameters of three different wedges were in Table 1 and the diagrammatic sketch is shown in Figure 3, where represents length of elastic wedge, is width, and is deadrise angle.

In the numerical calculation, the zero physical time () was selected as the moment in which the bottom of the wedge is 0.1 m up from the initial free surface and then the elastic wedge began to fall freely. The physical time was 0.0001 s, fictitious time was 0.01 s, the artificial compressibility coefficient was 500, the gravity acceleration was 9.81 m/s^{2}, and the length of structural element on solid boundary was 0.01 m. Furthermore, the inhomogeneous meshes were generated in fluid field. Near solid boundary the local meshes with equal spacing m were refined and regular and then the mesh size gradually enlarged away from boundary. By taking the model 1, for example, an inhomogeneous mesh system 66 × 90 is shown in Figure 4.
(a)
(b)
In Figure 5, the time histories of global acceleration and local strain of observation point on the sloping edge versus numerical solution and experiment data for three different elastic wedges are compared. Therein Figures 5(a), 5(b), and 5(c) represent various wedge models (number 1~3 in Table 1), respectively. From the figure, it can be found that the numerical solutions by this paper and experimental data are in good agreement with each other, except acceleration peak and negative strain of elastic wedge. On the one hand, for the sake of resistance of experimental equipment, wedge accelerations in experiment before water entry are less than the numerical simulation, which should cause a little difference about acceleration peak. On the other hand, by means of numerical method, the absolute value of negative strain on solid boundary is a bit underpredicted. The reason may be that the twoway “weak” coupling strategy between fluid and local structure artificially increases the numerical stiffness of waterentry model.
(a)
(b)
(c)
5.2. Description of Free Surface Profile and Velocity Vectors
By using the present method to study waterentry model 1 of Table 1, the developing process of free surface shape and velocity vectors in flow field can be successfully numerically simulated and the results at various waterentry moments ( s, 0.12 s, 0,18 s, 0.24 s, 0.30 s, 0.42 s, 0.48 s, and 0.60 s) are shown in Figure 6 to describe the waterentry phenomenon of the freefalling elastic wedge.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
At s the first graph shows the initial location of elastic wedge and calm free surface. Then the elastic wedge freely falls down ( s) and the air near solid surface moves with the wedge. As the body pieces the free surface ( s), the water rises up along each sloping side and the displacement waves are formed originating at the vertex of wedge. Meanwhile, the air between the solid boundary and free surface is forced around wedge edge and separates at the corner to form a pair of symmetric vortices behind the body in air. As the body travels further into the water ( s and 0.30 s), the two intersections between wedge and free surface move up along both straight sides and subsequently overturn towards each other ( s). After the free surface is reconnected ( s), a jet flow is created upwards above the wedge ( s).
5.3. Analysis on Hydroelasticity of FreeFalling Wedge
In order to study the global dynamic behaviors and local structural responses of elastic wedge, here, based on the geometric parameters of model 1 in Table 1, two kinds of waterentry models versus rigid and elastic wedge are created simultaneously. By comparing with the numerical results of different models, the hydroelasticity of freefalling wedge during the waterentry phase is studied and analyzed.
5.3.1. Global Dynamic Behaviors of Elastic Wedge
Firstly, by means of the current numerical method with or without elastic deformation, the global hydrodynamic acceleration and velocity versus elastic and rigid wedge model were calculated and compared as shown in Figure 7.
In Figure 7, for the sake of structural deformation, some highfrequency vibrations exist in the variation curves of global hydrodynamic and acceleration for the elastic model. In order to directly and clearly study the dynamic behaviors of elastic wedge, here a SavitzkyGolay method with 5 points and 4 degrees was used to filter highfrequency oscillation and smooth variable curves. For the rigid and elastic model, the time history of hydrodynamic and acceleration are in substantial agreement except two little differences. One is the peak value of hydrodynamic and acceleration curves and the results of rigid model are a little larger than the elastic model. The reason may be that the elastic deformation can provide buffer action and weaken fluid hydrodynamic during the waterentry phase. The other is that there are many highfrequency vibrations in hydrodynamic and acceleration curves, which can be caused by local deformation of elastic structure. However, the amplitudes of these vibrations are too small to influence the global motion of wedge. Thus, the velocity curves versus two models are basically in agreement with each other.
On this basis, in order to further study the long global motion of freefalling wedge during the water entry phase, here the water entry model was numerically simulated for a long physical time and the results of wedge velocity were shown in the last picture of Figure 7. From the picture, the motion process of freefalling wedge during water entry phase can be described as follows: firstly the wedge in air free falls straight down and then slows down due to large hydrodynamic at the initial moment of waterentry. With the continuous hydrodynamic acting on the wedge, the velocity gradually decreases until zero. Next, because the buoyancy of wedge is larger than its gravity, along with the closure of free surface, the wedge begins to move upwards. Subsequently, the jet successively grows upwards, falls downwards, and finally disappears on free surface, which will influence the global motion of wedge. Furthermore, by combining the effect of fluid force and gravity, wedge should vibrate up and down with less and less amplitude.
Figure 8 shows the pressure contour in fluid field at classic waterentry moments ( s, 0.18 s, 0.21 s, and 0.24 s). In the figure, at the initial stage of waterentry ( s), wedge has just begun to pierce free surface and the pressure peak appears near the vertex of wedge. Furthermore, the air escapes with high velocity from the gap of free surface and solid boundary to both sides and hence forms a strong airflow to press on the wedge edges. With the wedge continuously entering water ( s), free surface rises up along each sloping side and the amplitude and influence sphere of fluid pressure on wedge edge gradually enlarge. Furthermore, it can be found that the variation rate of pressure distribution on the wedge edge is related to the location of free surface. Therein the pressure slowly changes in water and rapidly reduces in air. As the body travels further into the water ( s and 0.24 s), the two intersections between wedge edge and free surface move up through corners and along both straight sides. Because of the continuous effect of fluid hydrodynamic, the wedge velocity becomes smaller and smaller, which results in the reduction of pressure distribution on wedge edge. Furthermore, it should be noted that, because near the unsmooth corners the fluid velocity become very large, the negative pressure on solid boundary should take place over there.
5.3.2. Local Structural Response on Elastic Edge
The time history and frequency analysis of local pressure and strain on the midpoint of wedge edge are shown in Figure 9. In order to better discuss the characteristic of local structural deformation; here fast Fourier transform method is applied in timefrequency analysis and the local pressure and strain in frequency domain are achieved and shown. From the frequency analysis, it can be found that the variations of local pressure and stain on the midpoint mainly include two characteristics.
One is the pressure peak with high amplitude and narrow range in low frequency domain. Based on the effect of max pressure (frequency tends to zero) in Figure 9 and global hydrodynamic in Figure 7, elastic boundaries should be deformed with 1st mode of vibration. Furthermore, the corresponding frequency of max strain is 1.3021 Hz and approximately equal to the first order natural frequency (1.2241 Hz) of the elastic model which is calculated by ANSYS software.
The other is the pressure oscillation with approximately 72~85 Hz range and low amplitude in high frequency domain which includes two crests and one trough. From the modal analysis, the pressure oscillation possibly causes the midpoint to locally vibrate with 11th and 12th natural frequencies 73.13 and 86.822 Hz which are calculated by ANSYS software. Furthermore, because the two natural frequencies are very close, the vibration phenomenon of “beats” should appear on the midpoint of wedge edge. Here according to the properties of “beats” phenomenon, the theoretical frequencies versus amplitude profile and highfrequency vibration are 6.846 and 79.976 Hz, respectively. Correspondingly, from the variation of local strains in frequency and time domain as shown in Figure 9, it can be found that the calculated frequencies (7.33 and 75.52 Hz) are close to the theoretical values, which theoretically and numerically verifies the “beats” phenomenon of elastic structural vibration.
Thus the local strain of midpoint on wedge edge can be divided into main deformation versus 1st principal mode and local deformation of “beat” phenomenon versus 11th and 12th modes.
In order to evaluate the hydroelastic effect of freefalling wedge, here case 1 is numerically simulated without feedback loop of elastic deformation, which means that the local structural response is handled as just a “pure” postprocessing step. Then the results are shown and compared with twoway “weak” coupling method in Figures 10 and 11.
(a)
(b)
(c)
(d)
(e)
(f)
(a)
(b)
At initial waterentry moments (from s to 0.20 s), for the pressure distributed on elastic boundaries (as shown in Figure 10) and local pressure on the midpoint (from Figure 11(a)), the results of “pure” postprocessing treatment without feedback loop are mostly larger than the twoway “weak” coupling method, which consequently causes the strain of “pure” postprocessing treatment to irrationally and impractically enlarge during waterentry phase as shown in Figure 11(b). Hence, from Figures 10 and 11, it can be found that the feedback of structural response on CFD boundary condition is very important to ensure the calculation precision. If the feedback loop is neglected, the local structural deformation on elastic boundaries should be overvalued.
In Figure 12, the time history of local structural deformation and global displacement of the wedge vertex are calculated and drawn. From the figure, during the water entry phase, the  and structural deformation and rotationangle are under 0.08 μm, 10 μm, and 1.5 μrad, respectively. By comparing with the wedge thickness 0.0015 m, the deformations of wedge vertex are very small and can be neglected. Furthermore, from the last picture, it also can be found that the two time histories of global displacement of the wedge vertex versus rigid and elastic model are almost the same. Thus, a conclusion can be made that, for some waterentry cases, the vertex of elastic wedge can be handled as rigidpoint and so the simplified halfwedge model can be efficiently applied in the numerical simulation to enhance computational efficiency.
6. Conclusion
(1)In this paper, surface capturing method and Cartesian cut cell mesh are successfully developed to handle the interaction among fluid field, global motion, and local deformation of elastic body simultaneously for water entry model of freefalling wedge.(2)In finite element analysis, based on the different characteristic of fluid and structural mesh, a particular data transfer method on solid boundary is successfully deduced and presented. Furthermore, an effective fluidstructure coupled solution strategy is taken to solve the interaction among fluid field, global motion, and local deformation of elastic body.(3)Based on the present method, some test cases of water entry for freefalling elastic wedge are numerically simulated. By comparing with experimental data, the results show that the present method has a good ability to study hydroelastic waterentry problems of freefalling elastic wedge.(4)Through further analysis, the global hydrodynamic performance and physical phenomenon of water entry problem (such as free surface rising, pressure distribution, jetflow, and negative pressure) for elastic wedge are discussed in detail.(5)For the local structural response on elastic boundaries, especial “beats” phenomenon of elastic vibration versus certain model without structural damping are discussed. Furthermore, based on the results of vertex for elastic wedge, a conclusion has been made that, for some test cases of elastic wedge, simplified half model can be applied in the numerical simulation to enhance computational efficiency.
Furthermore, the method should be developed further on the following aspects.(1)Although the present method is available and accurate for the waterentry model of wedge, it has trouble in waterentry problems of solid body with complex shape because of the invalidation of lathbeam structural model and a more finite element model (shell or plate) should be applied. Furthermore, the current fluidstructure coupling strategy should be improved and gradually transformed into other advanced coupling approaches (such as twoway “tight” coupling method).(2)Because here, the structural finite element analysis is based on linear elastic theory, the present method is not suitable for water entry problems of elastic body with large deformation. Hence in future the geometric nonlinearity should be taken into account to extend the application fields.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors are grateful to Dr Qian of the Manchester Metropolitan University in the UK for providing very useful theoretical help and many suggestions about the numerical method. This work was financially supported by the National Innovation Team Foundation under Grant no. 50921001 (China) and the National Natural Science Foundation of China (Grant no. 11202047).
References
 T. von Karman, “The impact of seaplane floats during landing,” Technical Note 321, NACA, Washington, Wash, USA, 1929. View at: Google Scholar
 H. Wagner, “Über Stoss und Gletivorgänge an der Oberfläche von Flüssigkeiten,” Zeitschriff für Angewandte Mathematik und Mechanik, vol. 12, no. 4, pp. 193–235, 1932. View at: Google Scholar
 Y. Yamamoto, K. Iida, T. Fukasawa, T. Murakami, M. Arai, and A. Ando, “Coupled analysis of nonlinear interaction between fluid and structure during impact,” International Shipbuilding Progress, vol. 32, no. 369, pp. 124–136, 1985. View at: Google Scholar
 A. M. Worthington, “Impact with a liquid surface studied with aid of instantaneous photography,” Philosophical Transactions of the Royal Society A, vol. 189, pp. 137–148, 1897. View at: Google Scholar
 G. H. Bottomley, “The impact of a model seaplane floats on water: reports and Memoranda,” Tech. Rep. 583, NACA, Washington, DC, USA, 1919. View at: Google Scholar
 W. L. Mayo, Hydrodynamic Impact of a System with a Single Elastic Mode, NACA, Washington, DC, USA, 1947.
 G. Stoffmacher, “Hydrodynamic impact of conicalnosed vehicles during vertical water entry,” Journal of Spacecraft and Rockets, vol. 1, no. 2, pp. 222–223, 1964. View at: Google Scholar
 S. L. Chuang, “Investigation of impact of rigid and elastic bodies with water,” Naval Warfare and Marine Engineering, 1970. View at: Google Scholar
 O. F. Hughes, “Solution of wedge entry problem by numerical conformal mapping,” Journal of Fluid Mechanics, vol. 56, no. 1, pp. 173–192, 1972. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. K. Ochi and L. E. Motter, “Prediction of slamming characteristics and hull response for ship design,” Transactions SNAME, vol. 81, pp. 144–176, 1973. View at: Google Scholar
 T. Belytschko and R. Mullen, “Twodimensional fluidstructure impact computations with regularization,” Computer Methods in Applied Mechanics and Engineering, vol. 27, no. 2, pp. 139–154, 1981. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A. A. Korobkin and V. V. Pukhnachov, “Initial asymptotic in problem of blunt body entrance into liquid,” in Proceedings of the 3rd International Conference on Numerical Ship Hydrodynamics, pp. 579–591, 1981. View at: Google Scholar
 A. Korobkin and V. V. Pukhnachov, “Initial asymptotic in contact hydrodynamics problems,” in Proceedings of the 4th International Conference on Numerical Ship Hydrodynamics, pp. 138–151, 1985. View at: Google Scholar
 A. Korobkin and V. V. Pukhnachov, “Initial stage of water impact,” Annual Review of Fluid Mechanics, vol. 20, pp. 159–185, 1988. View at: Google Scholar
 A. A. Korobkin and V. V. Pukhnachov, “Initial stage of water impact,” Annual Review of Fluid Mechanics, vol. 20, pp. 159–185, 1998. View at: Google Scholar
 M. Greenhow and L. Yanbao, “Added masses for circular cylinders near or penetrating fluid boundariesreview, extension and application to waterentry, exit and slamming,” Ocean Engineering, vol. 14, no. 4, pp. 325–348, 1987. View at: Publisher Site  Google Scholar
 R. Cointe and J. L. Armand, “Hydrodynamic impact analysis of a cylinder,” Journal of Offshore Mechanics and Arctic Engineering, vol. 109, pp. 237–243, 1989. View at: Google Scholar
 S. K. Wilson, “A mathematical model for the initial stages of fluid impact in the presence of a cushioning fluid layer,” Journal of Engineering Mathematics, vol. 25, no. 3, pp. 265–285, 1991. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 T. Miloh, “On the oblique waterentry problem of a rigid sphere,” Journal of Engineering Mathematics, vol. 25, no. 1, pp. 77–92, 1991. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 R. Zhao and O. Faltinsen, “Water entry of twodimensional bodies,” Journal of Fluid Mechanics, vol. 246, pp. 593–612, 1993. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 O. M. Faltinsen and R. Zhao, “Water entry of ship sections and axisymmetric bodies,” in AGARD FDP Workshop on High Speed Body Motion in Water, Kiev, Ukraine, 1997. View at: Google Scholar
 K. Takagi, “Influence of elasticity on hydrodynamic impact Problem,” Journal of the Kansai Society of Naval Architects, vol. 222, pp. 97–106, 1994. View at: Google Scholar
 L. E. Fraenkel and J. B. McLeod, “Some results for the entry of a blunt wedge into water,” Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, vol. 355, no. 1724, pp. 523–535, 1997. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 G. X. Wu, “Hydrodynamic force on a rigid body during impact with liquid,” Journal of Fluids and Structures, vol. 12, no. 5, pp. 549–559, 1998. View at: Publisher Site  Google Scholar
 G. X. Wu, H. Sun, and Y. S. He, “Numerical simulation and experimental study of water entry of a wedge in free fall motion,” Journal of Fluids and Structures, vol. 19, no. 3, pp. 277–289, 2004. View at: Publisher Site  Google Scholar
 G. X. Wu, “Numerical simulation of water entry of twin wedges,” Journal of Fluids and Structures, vol. 22, no. 1, pp. 99–108, 2006. View at: Publisher Site  Google Scholar
 M. Greenhow and S. Moyo, “Water entry and exit of horizontal circular cylinders,” Philosophical Transactions of the Royal Society of London A, vol. 355, no. 1724, pp. 551–563, 1997. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 P. C. Sames, T. E. Schellin, S. Muzaferija, and M. Peric, “Application of a twofluid finite volume method to ship slamming,” Journal of Offshore Mechanics and Arctic Engineering, vol. 121, no. 1, pp. 47–52, 1999. View at: Publisher Site  Google Scholar
 X. M. Mei, Y. M. Liu, and D. K. P. Yue, “On the water impact of general twodimensional sections,” Applied Ocean Research, vol. 21, no. 1, pp. 1–15, 1999. View at: Publisher Site  Google Scholar
 C. H. Lu, Y. S. He, and G. X. Wu, “Coupled analysis of nonlinear interaction between fluid and structure during impact,” Journal of Fluids and Structures, vol. 14, no. 1, pp. 127–146, 2000. View at: Publisher Site  Google Scholar
 Y.M. Scolan and A. A. Korobkin, “Energy distribution from vertical impact of a threedimensional solid body onto the flat free surface of an ideal fluid,” Journal of Fluids and Structures, vol. 17, no. 2, pp. 275–286, 2003. View at: Publisher Site  Google Scholar
 N. Aquelet and M. Souli, “Damping effect in fluidstructure interaction: application to slamming problem,” in Proceedings of the ASME Pressure Vessels and Piping Conference, pp. 233–242, Cleveland, Ohio, USA, July 2003. View at: Google Scholar
 N. Aquelet and M. Souli, “Fluidstructure coupling in a waterwedge impact problem,” in Proceedings of the ASME/JSME Pressure Vessels and Piping Conference, pp. 91–99, San Diego, Calif, USA, July 2004. View at: Google Scholar
 N. Aquelet, M. Souli, and L. Olovsson, “EulerLagrange coupling with damping effects: application to slamming problems,” Computer Methods in Applied Mechanics and Engineering, vol. 195, no. 1–3, pp. 110–132, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 D. Battistin and A. Iafrati, “A numerical model for the jet flow generated by water impact,” Journal of Engineering Mathematics, vol. 48, no. 34, pp. 353–374, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 K. M. T. Kleefsman, G. Fekken, A. E. P. Veldman, and B. Buchner, “A volumeoffluid based simulation method for wave impact problems,” Journal of Computational Physics, vol. 206, no. 1, pp. 363–393, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. Sun, A Boundary Element Method Applied to strongly Nonlinear WaveBody Interaction Problems, Norwegian University of Science and Technology, Trondheim, Norway, 2007.
 O. M. Faltinsen and Y. A. Semenov, “Nonlinear problem of flatplate entry into an incompressible liquid,” Journal of Fluid Mechanics, vol. 611, pp. 151–173, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 T. Tveitnes, A. C. FairlieClarke, and K. Varyani, “An experimental investigation into the constant velocity water entry of wedgeshaped sections,” Ocean Engineering, vol. 35, no. 1415, pp. 1463–1478, 2008. View at: Publisher Site  Google Scholar
 T. T. Truscott and A. H. Techet, “Water entry of spinning spheres,” Journal of Fluid Mechanics, vol. 625, pp. 135–165, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 G. D. Xu, W. Y. Duan, and G. X. Wu, “Simulation of water entry of a wedge through free fall in three degrees of freedom,” Proceedings of The Royal Society of London A: Mathematical, Physical and Engineering Sciences, vol. 466, no. 2120, pp. 2219–2239, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 C. W. Hirt and B. D. Nichols, “Volume of fluid (VOF) method for the dynamics of free boundaries,” Journal of Computational Physics, vol. 39, no. 1, pp. 201–225, 1981. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 D. L. Youngs, “Time dependent multimaterial flow with large fluid distortion,” in Numerical Methods for Fluid Dynamics, K. W. Morton and M. J. Baines, Eds., pp. 273–285, Academic Press, London, UK, 1982. View at: Google Scholar
 J. M. Hill, Y. H. Wu, and B. Wiwatanapataphee, “Analysis of flux flow and the formation of oscillation marks in the continuous caster,” Journal of Engineering Mathematics, vol. 36, no. 4, pp. 311–326, 1999. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 F. J. Kelecy and R. H. Pletcher, “The development of a free surface capturing approach for multidimensional free surface flows in closed containers,” Journal of Computational Physics, vol. 138, no. 2, pp. 939–980, 1997. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 D. Pan and C. H. Chang, “The capturing of fee surfaces in incompressible multifluid flows,” International Journal of Numerical Methods, vol. 33, pp. 203–222, 2000. View at: Google Scholar
 L. Qian, D. M. Causon, D. M. Ingram, and C. G. Mingham, “A freesurface capturing method for two fluid flows with moving bodies,” Proceedings of the Royal Society A, vol. 462, pp. 21–42, 2006. View at: Google Scholar
 W. Wang and Y. Wang, “An improved free surface capturing method based on Cartesian cut cell mesh for waterentry and exit problems,” Proceedings of The Royal Society of London A, vol. 465, no. 2106, pp. 1843–1868, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 W. J. Coirier and K. G. Powell, “An accuracy assessment of Cartesianmesh approaches for the Euler equations,” Journal of Computational Physics, vol. 117, no. 1, pp. 121–131, 1995. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 G. Yang, D. M. Causon, D. M. Ingram, R. Saunders, and P. Batten, “A cartesian cut cell method for compressible flows part A: atatic body problems,” Aeronautical Journal, vol. 101, no. 1002, pp. 47–56, 1997. View at: Google Scholar
 G. Yang, D. M. Causon, D. M. Ingram, R. Saunders, and P. Battent, “A cartesian cut cell method for compressible flows part B: moving body problems,” Aeronautical Journal, vol. 101, no. 1002, pp. 57–65, 1997. View at: Google Scholar
 D. M. Causon, D. M. Ingram, and C. G. Mingham, “A Cartesian cut cell method for shallow water flows with moving boundaries,” Advances in Water Resources, vol. 24, no. 8, pp. 899–911, 2001. View at: Publisher Site  Google Scholar
 D. Pan and H. Lomax, “A new approximate $LU$ factorization scheme for the Reynoldsaveraged NavierStokes equations,” The American Institute of Aeronautics and Astronautics, vol. 26, no. 2, pp. 163–171, 1988. View at: Publisher Site  Google Scholar  MathSciNet
 H. Y. Shu and L. S. Tan, Structural Mechanics for Ships, Huazhong University of Science and Technology Press, Wuhan, China, 1993 (Chinese).
 H. Sun, Z. H. Lu, and Y. S. He, “Experimental research on the fluidstructure interaction in water entry of 2D elastic wedge,” Journal of Hydrodynamics, vol. 18, no. 1, pp. 104–109, 2003 (Chinese). View at: Google Scholar
Copyright
Copyright © 2014 Wenhua Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.