• Ocean Systems Engineering
• /
• v.4 no.2
• /
• pp.83-97
• /
• 2014

Wave propagation in a three-dimensional (3D) fully nonlinear numerical wave tank (NWT) is studied based on velocity potential theory. The governing Laplace equation with fully nonlinear boundary conditions on the moving free surface is solved using the indirect desingularized boundary integral equation method (DBIEM). The fourth-order predictor-corrector Adams-Bashforth-Moulton scheme (ABM4) and mixed Eulerian-Lagrangian (MEL) method are used for the time-stepping integration of the free surface boundary conditions. A smoothing algorithm, B-spline, is applied to eliminate the possible saw-tooth instabilities. The artificial wave speed employed in MTF (multi-transmitting formula) approach is investigated for fully nonlinear wave problem. The numerical results from incorporating the damping zone (DZ), MTF and MTF coupled DZ (MTF+DZ) methods as radiation condition are compared with analytical solution. An effective MTF+DZ method is finally adopted to simulate the 3D linear wave, second-order wave and irregular wave propagation. It is shown that the MTF+DZ method can be used for simulating fully nonlinear wave propagation very efficiently.

• Journal of the Korean Society for Precision Engineering
• /
• v.29 no.6
• /
• pp.632-639
• /
• 2012

A fully coupled thermo-mechanical model is adopted to study the temperature distribution and the material deformation in friction stir welding(FSW) process. Rotational speed is most important parameters in this research. Three dimension results under different process parameters were presented. Result indicate that the maximum temperature is lower than the melting point of the welding material. The higher temperature gradient occurs in the leading side of the workpiece. The maximum temperature can be increased with increasing the tool angular velocity, rpm in the current numerical modeling. In this research ABAQUS Ver.6.7 is to analyze a fully coupled thermo-mechanical model. ALE(Arbitrary Lagrangian-Eulerian) finite element formulation is used for the large deformation in FSW process and using the Mass scaling for the analysis time efficiency.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.26 no.1
• /
• pp.1-17
• /
• 2024

Analyzing the tunnel excavation behavior and its effect on the surrounding ground involves large deformation behavior. Therefore, in order to properly simulate the tunnel excavation process and rigorously investigate the actual effect of excavation on surrounding ground and tunnel structure large deformation analysis method is required. In this study, two major numerical approaches capable of considering large deformations behavior were applied to investigate the effect of tunnel boring machine excavation on the surrounding ground: coupled Eulerian-Lagrangian (CEL) and the automatic remeshing (AR) method. Relative performance of both approaches was evaluated through the ground response due to TBM excavation. The ground response will be quantified by estimating the range of the excavation damaged zone (EDZ). By comparing the results, the range of the EDZ will be suggested on the vertical and horizontal direction along the TBM excavation surface. Based on the computed results, it was found that the size of EDZ around the excavation surface and the tendencies was in good agreement among the two approaches. Numerical results clearly show that the size of the EDZ around the tunnel tends to be larger for rock with higher RMR rating. The size of the EDZ is found to be direct proportional to the tunnel diameter, whereas the depth of the tunnel is inversely proportional due to higher confinement stress around the excavation surface.

• Coupled systems mechanics
• /
• v.1 no.4
• /
• pp.361-379
• /
• 2012

This paper presents a fully coupled three-dimensional solver for the analysis of interaction between pulsatile flow and large deformation structure. A partitioned time marching algorithm is employed for the solution of the time dependent coupled discretised problem, enabling the use of highly developed, robust and well-tested solvers for each field. Conservative transfer of information at the fluid-structure interface is combined with an effective multi-predict-correct iterative scheme to enable implicit coupling of the interacting fields at each time increment. The three-dimensional unsteady incompressible fluid is solved using a powerful implicit time stepping technique and an ALE formulation for moving boundaries with second-order time accurate is used. A full spectrum of total variational diminishing (TVD) schemes in unstructured grids is allowed implementation for the advection terms and finite element shape functions are used to evaluate the solution and its variation within mesh elements. A finite element dynamic analysis of the highly deformable structure is carried out with a numerical strategy combining the implicit Newmark time integration algorithm with a Newton-Raphson second-order optimisation method. The proposed model is used to predict the wave flow fields of a particular flow-induced vibrational phenomenon, and comparison of the numerical results with available experimental data validates the methodology and assesses its accuracy. Another test case about three-dimensional biomedical model with pulsatile inflow is presented to benchmark the algorithm and to demonstrate the potential applications of this method.

• Journal of Ocean Engineering and Technology
• /
• v.33 no.6
• /
• pp.556-563
• /
• 2019

The ice keel gouge and seabed interaction is one of the major considerations in the design of an Arctic pipeline system. Ice keel and seabed interaction engineering models based on experimental data, which give an explicit equation for estimating the ice gouging depth, have been suggested. The suggested equations usually overestimate the ice keel gouging depth. In addition, various types of numerical analyses have been carried out to verify the suggested engineering model equations in comparison to the experimental data. However, most of numerical analysis results were also overestimated compared with the laboratory experimental data. In this study, a numerical analysis considering the contact condition and geostatic stress was carried out to predict the ice keel gouging depth and compared with the previous studies. Considering the previously mentioned conditions, more accurate results were produced compared with the laboratory experiment results and the error rate was reduced compared to previous numerical analysis studies.

• Journal of Advanced Research in Ocean Engineering
• /
• v.1 no.4
• /
• pp.198-210
• /
• 2015

This paper presents the damage assessment of a free-fall dropped object on the seabed. The damage to a dropped object totally depends on the relationship between the impact energy and the soil strength at the mudline. In this study, unexpected dropping scenarios were first assumed by varying the relevant range of the impact velocity, structure geometry at the moment of impact, and soil strength profile along the penetration depth. Theoretical damage assessments were then undertaken for a free-fall dropping event with a fixed final embedment depth for the structure. This paper also describes the results of a three-dimensional large deformation finite element analysis undertaken for the purpose of validation. The analyses were carried out using the coupled Eulerian-Lagrangian approach, modifying the simple elastic-perfectly plastic Tresca soil model. The validation exercises for each dropping scenario showed good agreement, and the present numerical approach was capable of predicting the behavior of a free-fall dropped object.

• Journal of Ocean Engineering and Technology
• /
• v.31 no.6
• /
• pp.397-404
• /
• 2017

In this study, an anchor dragging simulation was performed using the CEL method to design a rock-berm, which is a protection method for submarine cables. In order to simulate an anchor drag, preliminary simulations were first performed to determine the initial anchor penetration depth, anchor drag velocity, drag angle, and distance between the anchor and rock-berm. Based on the preceding simulation results, a safe rock-berm design for protecting the submarine cables was simulated to calculate the anchor penetration depth by the anchor dragging. As a result, the penetration depth of the anchor was found to be shallower in a hard seabed, and the penetration depth was deeper in a soft seabed, the height of the rock-berm was determined according to the physical properties of the seabed.

• Wind and Structures
• /
• v.15 no.3
• /
• pp.189-208
• /
• 2012

Lock-in and drag amplification phenomena are studied for a flexible cantilever using a simplified fluid-structure interaction approach. Instead of solving the 3D domain, a simplified setup is devised, in which 2D flow problems are solved on a number of planes parallel to the wind direction and transversal to the structure. On such planes, the incompressible Navier-Stokes equations are solved to estimate the fluid action at different positions of the line-like structure. The fluid flow on each plane is coupled with the structural deformation at the corresponding position, affecting the dynamic behaviour of the system. An Arbitrary Lagrangian-Eulerian (ALE) approach is used to take in account the deformation of the domain, and a fractional-step scheme is used to solve the fluid field. The stabilization of incompressibility and convection is achieved through orthogonal quasi-static subscales, an approach that is believed to provide a first step towards turbulence modelling. In order to model the structural problem, a special one-dimensional element for thin walled cross-section beam is implemented. The standard second-order Bossak method is used for the time integration of the structural problem.

• Journal of Ocean Engineering and Technology
• /
• v.32 no.3
• /
• pp.177-183
• /
• 2018

This study simulated ice load and the motion response of a moored semi-submersible rig in pack-ice conditions using a finite element method. Ice flows of random size and shape were modeled, and interactions for ice-sea, ice-structure, ice-ice were simulated using a simplified method. Parameters for the simplified method such as drag force coefficient and the pressure-penetration relation were obtained based on the result of detailed analysis using the coupled Eulerian-Lagrangian method. The mooring lines were modeled by spring elements based on their stiffness. As a result of the simulation over 1,400 seconds, the force and motion response of the rig were obtained and validated using discrete elements and compared with the results found by the Krylov State Research Centre.

• Journal of the Society of Naval Architects of Korea
• /
• v.43 no.2 s.146
• /
• pp.177-185
• /
• 2006

Cavitation is the dominant noise source of the marine vehicle. Of the various types of cavitation , tip vortex cavitation is the first appearance type of marine propeller cavitation and it generates high frequency noise. In this study, tip vortex cavitation behavior and noise are numerically investigated. A numerical scheme using Eulerian flow field computation and Lagrangian particle trace approach is applied to simulate the tip vortex cavitation on the hydrofoil. Vortex flow field is simulated by combined Moore and Saffman's vortex core radius equation and Sculley vortex model. Tip vortex cavitation behavior is analyzed by coupled Rayleigh-Plesset equation and trajectory equation. The cavitation nuclei are distributed and released in the vortex flow result. Vortex cavitation trajectories and radius variations are computed according to nuclei initial size. Noise is analyzed using time dependent cavitation bubble position and radius data. This study may lay the foundation for future work on vortex cavitation study and it will provide a basis for proper underwater propeller noise control strategies.