Further Application of Surface Capturing Method and Cartesian Cut Cell Mesh on Hydroelastic Water-Entry Problems of Free-Falling Elastic Wedge

In order to study the interactions between fluid and elastic structure (such as marine lifeboat falling down and ship), this paper presents a newCFDmethod on hydroelastic water-entry problem of free-falling elastic wedge, which canmore 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 lath-beam 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 fluid-structure interaction. Finally, by comparing numerical results with experimental data, the present method is validated to be available and feasible for hydroelastic water-entry problem and further successfully adopted to analyze the motion characteristics of free-falling elastic wedge.


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, green-water 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 fluid-structure interaction.
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][43][44] and so on), the surface capturing method [45][46][47][48] based on Cartesian cut cell mesh [47][48][49][50][51][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 2 Mathematical Problems in Engineering 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 water-entry problem of free-falling 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 dual-time stepping technique with artificial compressibility method is used for time advancing.On the other hand, in the structural solver finite element method and lath-beam structural model are used to analyze structural response of elastic wedge.Particularly based on the geometry of Cartesian cut cell mesh and lath-beam 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 free-falling 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 water-entry phase.

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, V  and V  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   = 0 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 no-penetration 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 v  = V  i + V  j is fluid velocity at the solid boundary, n =   i +   j is normal unit vector of solid edge, and V  is normal projection of average velocity of solid edge.Furthermore, a detailed description about how to achieve V  and ∇V  based on the fluid-structure interaction on solid boundary should be demonstrated in Section 4.1.2 of data transfer method.

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).
Based on the Cartesian cut cell mesh, here a cell-central 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 least-square 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, dual-time 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].

Introduction of Lath-Beam Structural Model
For the bottom plate with cylindrical bending of elastic wedge structure, a lath-beam model can be used to analyze structural deformation.Compared with Euler-Bernoulli beam, due to bilateral constraints, the cross-section of lath-beam model cannot freely deform and keep the original shape.On this basis, only if by introducing  1 = /(1 −  2 ) ( is the material elastic modulus and  is Poisson ratio) as new elastic modulus into Euler-Bernoulli beam theory can the elastic deformation of lath-beam 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 Euler-Bernoulli beam theory should be applied to calculate structural response on solid boundaries of free-falling elastic wedge.Initially, the global (fixed) and local (moving with body) coordinate systems are created.Then, lath-beam structure is divided into many beam elements with equal length, which can support bending and tension-compression effects simultaneously.Here in the local coordinate system, the nodal displacement vector {}  of element  can be demonstrated as where  1 ,  1 ,  2 , and  2 are translational displacements, and  1 and  2 are rotation angles on two nodes.Next, in local coordinate system the corresponding mass and stiffness matrix of lath-beam 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 water-entry 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 lath-beam 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 lath-beam structural model,   is the structural density, ℎ is the element thickness,  is the element length,  = 1.0 * ℎ is the cross-section area,  = 1.0 * ℎ 3 /12 is the moment of inertia,  1 = /(1 −  2 ),  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.

Analysis of Fluid-Structure Interaction
For the water-entry model of free-falling elastic wedge, the global motion and local deformation of elastic body should be calculated from the fluid-structure interaction, so here the data transfer method on elastic boundary and fluid-structure coupling strategy are discussed, respectively.

Data Transfer Method on Elastic Boundary.
In order to study and analyze the fluid-structure 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 lath-beam 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.

Transfer of Distributional Fluid Pressure to Structural
Nodes.Based on respective geometries of Cartesian cut cell mesh and lath-beam element, there should be many kinds of position relations between the solid boundary  in cut cell and structural element  1  2 to be considered (three classic cases as shown in Figure 1), which are full inclusion of  in  1  2 (Figure 1(a)) or  1  2 in  (Figure 1(b)), partial intersection between  1  2 with  (Figure 1(c)), and some other situations (such as that  coincides with  1  2 ).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 Firstly, local coordinate system       should be created on structural  1  2 .Therein coordinateorigin   coincides with node  1 , coordinates axis     is clockwise and, along with element  1  2 , axis     points from fluid to structure.Thus unit vectors of     and     can be written: where ( 1 ,  1 ) and ( 2 ,  2 ) are coordinates of nodes  1 and  2 , respectively.Next, by using the coordinates (  ,   ) of centre  in cut cell, position vector OL 1 can be achieved: Then by combining ( 9) with (11), position vector Δr (from centre  to arbitrary point  on ) can be written: where   = |L 1 C| 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  1 and  2 on nodes  1 and  2 of structural lath-beam element  1  2 , which can be described as follows: where   and Δp =   i +   j 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.

Feedback of Structural Response on CFD Boundary
Condition.Here among many kinds of position relations of Cartesian cut cell mesh and lath-beam 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.
Firstly based on the shape function of structural element  1  2 , the velocities of point  in local coordinate system can be obtained: where  1 ∼  4 and  5 ∼  6 are the shape function versus bending and tension-compression effects, respectively.For the lath-beam element,  1 Then by the same method, the velocities u  and ẇ  of point  also can be calculated.
Here it should be noted that V  and ∇V  terms in ( 5) are the sum velocity, which includes local elastic deformation and global rigid translation.Subsequently, the average normal velocity V  and tangential velocity V  of solid boundary  can be calculated: where n and t are normal and tangential unit vectors of solid boundary and v  = V  i + V  j 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 ∇V  = (V  /)i + (V  /)j 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.

Fluid-Structure Coupled Solution Strategy.
In the present study, the fluid-structure interaction of water-entry model for free-falling 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 V   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 body is obtained by direct integration of the new pressure over body surface.When  = (V +1  ) +1 − (V +1  )  < 10 −3 , the subiteration is convergence and (V +1  ) +1 is taken as the new velocity boundary condition V +1  of the solid body for the next physical time step  + 1.
On the other hand, compared with military missile at high speed, ship and ocean engineering structures usually enter water with medium-low 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 two-way "weak" coupling method as shown in Figure 2 is used for the coupled solution of fluid and local structural deformation.

Results of Test Cases and Discussions
5.1.Validation.In this section, the current numerical method for hydroelastic analysis on water-entry problems of freefalling wedge should be tested by comparing calculated results with experimental data of Sun et al. [55] about three different water-entry 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 ( = 0) 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 Δ = Δ = 0.01 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.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.

Mathematical Problems in Engineering
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 two-way "weak" coupling strategy between fluid and local structure artificially increases the numerical stiffness of water-entry model.

Description of Free Surface Profile and Velocity Vectors.
By using the present method to study water-entry 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 water-entry moments ( = 0.00 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 water-entry phenomenon of the free-falling elastic wedge.At  = 0.00 s the first graph shows the initial location of elastic wedge and calm free surface.Then the elastic wedge freely falls down ( = 0.12 s) and the air near solid surface moves with the wedge.As the body pieces the free surface ( = 0.18 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 ( = 0.24 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 ( = 0.42 s).After the free surface is reconnected ( = 0.48 s), a jet flow is created upwards above the wedge ( = 0.60 s).

Analysis on Hydroelasticity of Free-Falling 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 water-entry models versus rigid and elastic wedge are created simultaneously.By comparing with the numerical results of different models, the hydroelasticity of free-falling wedge during the water-entry phase is studied and analyzed.

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 high-frequency 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 Savitzky-Golay method with 5 points and 4 degrees was used to filter high-frequency 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 water-entry phase.The other is that there are many high-frequency 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 free-falling 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 free-falling 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 water-entry.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 water-entry moments ( = 0.15 s, 0.18 s, 0.21 s, and 0.24 s).In the figure, at the initial stage of water-entry ( = 0.15 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 air-flow to press on the wedge edges.With the wedge continuously entering water ( = 0.18 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 ( = 0.18 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.

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 time-frequency 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 ( 1 −  2 )/4 and high-frequency vibration ( 1 +  2 )/4 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 free-falling 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 two-way "weak" coupling method in Figures 10 and 11.
At initial water-entry moments (from  = 0.15 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 two-way "weak" coupling method, which consequently causes the strain of "pure" postprocessing treatment to irrationally and impractically enlarge during water-entry phase as shown in 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 rigid-point and so the simplified half-wedge model can be efficiently applied in the numerical simulation to enhance computational efficiency.

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 free-falling 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 fluid-structure 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 free-falling elastic wedge are numerically  Mathematical Problems in Engineering simulated.By comparing with experimental data, the results show that the present method has a good ability to study hydroelastic water-entry problems of free-falling elastic wedge.
(4) Through further analysis, the global hydrodynamic performance and physical phenomenon of water entry problem (such as free surface rising, pressure distribution, jet-flow, 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 water-entry model of wedge, it has trouble in water-entry 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 fluid-structure 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.

Figure 1 :
Figure 1: Different position relations of cut cell and structural elements; (a) full inclusion of  in  1  2 ; (b) full inclusion of  1  2 in ; (c) partial intersection between  1  2 and .

Figure 2 :
Figure 2: Flow chart of iteration procedure for the fluid-structure coupled solution strategy.

Figure 3 :Figure 4 :
Figure 3: Sketch of elastic model and description of geometric parameters.

Figure 7 :
Figure 7: Time history of global hydrodynamic acceleration and velocity versus elastic and rigid wedge model with 7.787 kg mass and 45 ∘ deadrise angle.

Figure 8 :
Figure 8: Pressure contour in fluid field for elastic wedge with 7.787 kg mass and 45 ∘ deadrise angle;  = 0.15 s, 0.18 s, 0.21 s, and 0.24 s (from left to right and from up to bottom).

Figure 9 :
Figure 9: Time history and frequency analysis of local pressure and strain on the midpoint of wedge edge versus elastic and rigid model with 7.787 kg mass and 45 ∘ deadrise angle.

Figure 12 :
Figure 12: Time history of local structural deformation and global displacement on the vertex of wedge with 1.5 mm thickness, 7.787 kg mass and 45 ∘ deadrise angle.

Table 1 :
Geometric parameters of three different elastic wedges.
(19)Equation(19)demonstrates a subiterative process of fluidglobal body motion coupling at each physical moment.At physical time step , the iterative initial value (V +1  ) 0 is taken as V   .For the subiterative  step, based on the vertical body velocity (V +1  )  , a new pressure field is calculated and then