• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.1
• /
• pp.60-68
• /
• 1991

The updated Lagrangian Finite Element Method is introduced to analyse rigid body-fluid impact problem which is characterized by incompressible Navier-Stokes equations and impact-contact conditions between free surface and rigid body. For the convenience of numerical computation, velocity fields are splinted into vicous and pressure parts, and then the governing equations and boundary conditions are decomposed in accordance with the decomposition. However, Viscous stresses acting an the solid boundaries are neglected on the assumption that very small velocity gradients may occur during extremely small time interval of the impact. Four coded quadrilateral elements are used to discretize the space domain and the fully explicit time-marching algorithm is employed with a reasonably small time step. At the beginning of each time step, contact velocity of the rigid body is computed from the momentum balance between the body and the fluid. The velocity field is then computed to satisfy the discretized equations of motions and incompressibility and contact constraints as well as an exact free surface boundary condition. At the end of each time step, the fluid domain is updated from the velocity field. In the present time stepping numerical analysis, behaviour of the free surface near the body can be observed without any difficulty which is very important in the water impact problem. The applicability of the algorithm is illustrated by a wedge type falling body problem. The numerical solutions for time-varying pressure distributions and impact loadings acting ion the surface are obtained.

• Journal of the Korean Geotechnical Society
• /
• v.29 no.12
• /
• pp.35-44
• /
• 2013

In this study, a particle method without using grid was applied for analysing large deformation problems in soil flows instead of using ordinary finite element or finite difference methods. In the particle method, a continuum equation was discretized by various particle interaction models corresponding to differential operators such as gradient, divergence, and Laplacian. Soil behavior changes from solid to liquid state with increasing water content or external load. The Mohr-Coulomb failure criterion was incorporated into the particle method to analyze such three-dimensional soil behavior. The yielding and hardening behavior of soil before failure was analyzed by treating soil as a viscous liquid. First of all, a sand column test without confining pressure and strength was carried out and then a self-standing clay column test with cohesion was carried out. Large deformation from such column tests due to soil yielding or failure was used for verifying the developed particle method. The developed particle method was able to simulate the three-dimensional plastic deformation of soils due to yielding before failure and calculate the variation of normal and shear stresses both in sand and clay columns.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.5B
• /
• pp.559-573
• /
• 2008

In this study, we newly proposed 3-D nonlinear wave driver utilizing the Navier-Stokes Eq. the numerical integration of which is carried out using SPH (Smoothed Particle Hydrodynamics), an internal wave generation with the source function of Gaussian distribution and an energy absorbing layer. For the verification of new 3-D nonlinear wave driver, we numerically simulate the sloshing problem within a parabolic water basin triggered by a Gaussian hump and uniformly inclined water surface by Thacker (1981). It turns out that the qualitative behavior of sloshing caused by relaxing the external force which makes a free surface convex or uniformly inclined is successfully simulated even though phase error is visible and an inundation height shrinks as numerical simulation more proceeds. For the more severe test, we also simulate the nonlinear shoaling and refraction over uniform beach of wedge shape. It is shown that numerically simulated waves are less refracted than the linear counterpart by Hamiltonian ray theory due to nonlinearity, energy dissipation at the bottom and side walls, energy loss induced by breaking, and the hydraulic jump occurring when breaking waves encounter a down-rush by the preceding wave.