• The Journal of the Acoustical Society of Korea
• /
• v.39 no.2
• /
• pp.77-86
• /
• 2020

In this study, the Eulerian/Lagrangian one-way coupling method is proposed to predict flow noise due to Blade-Tip Vortex Cavitation (BTVC). The proposed method consists of four sequential steps: flow field simulation using Computational Fluid Dynamics (CFD) techniques, reconstruction of wing-tip vortex using vortex model, generation of BTVC using bubble dynamics model and acoustic wave prediction using the acoustic analogy. Because the CFD prediction of tip vortex structure generally suffers from severe under-prediction of its strength along the steamwise direction due to the intrinsic numerical damping of CFD schemes and excessive turbulence intensity, the wing-tip vortex along the freestream direction is regenerated by using the vortex modeling. Then, the bubble dynamics model based on the Rayleigh-Plesset equation was employed to simulate the generation and variation of BTVC. Finally, the flow noise due to BTVC is predicted by modeling each of spherical bubbles as a monople source whose strength is proportional to the rate of time-variation of bubble volume. The validity of the proposed numerical methods is confirmed by comparing the predicted results with the measured data.

• Journal of the Korea Computer Graphics Society
• /
• v.25 no.5
• /
• pp.21-30
• /
• 2019

We propose a physical simulation method based on the shallow water equation(SWE) to represent water surface effectively. In this paper, the water which can be represented has a much larger width compared to the depth does not have a large vertical direction flow. In order to calculate the water flow efficiently, we start with the shallow water equation as the governing equation, which is a simplified version of the Navier-Stokes equation. In order to numerically calculate the advection term of the SWE, we introduce a new conservtive USCIP(CUSCIP) method which improves the Constrained Interpolation Profile (CIP) method to preserve the physical quantity while increasing the numerical accuracy. The proposed method is based on Kim et. al.'s Unsplit Semi-lagrangian CIP[9], and calculates advection term with additional constraints on term that consider integral values. The experimental results show that the CUSCIP method is robust to the loss of physical quantity due to numerical dissipation, which improves wave detail and persistence.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.27 no.1
• /
• pp.24-34
• /
• 1990

In this paper, the formulation of a new simplified finite element is made to analyze the ultimate strength of damaged tubular members subjected to combined axial force and end moment. A damaged tubular member that has the bending deformation and the local dent is modeled by beam elements. Tangent elastic stiffness matrix of a beam element which contains the effect of the geometric nonlinearity is derived by using the updated Lagrangian approach. Here the contribution of the stiffness in the dented area is neglected since its resistance against the external loads is considered to be small. A fully plastic interaction curve of the element under combined loads taking account of the local dent effect is selected as a yielding criterion at each nodal point. Also tangent elasto-plastic stiffness matrix of the element is formulated by plastic node method. Comparison with the present solution and the existing experimental results is made showing that the present method gives quite an accurate solution.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.8 no.1
• /
• pp.37-43
• /
• 1996

Breaking waves were generated in a 2-D flume. A piston-type wavemaker was operated in accordance with signals which consist of elementary harmonics with appropriate phase differences. These phase differences were estimated by using a linear wave theory so that wave crests were to be concentrated at the same position. The stroke of wavemaker was controlled to create plunging-type breaking waves. The signal with small amplitude could not generate breaking waves. In the case of moderate amplitudes, various breaking waves could be obtained. Most of breaking waves were spilling type. Only when the wavemaker was operated with appropriate amplitude, plunging-type breaking waves were generated. The parameters of breaking waves are the wave steepness and the frequency bandwidth. If the central frequency was low, breaking waves were not generated. Based on experimental data, we found that the wave height of breaking inception was H = 0.0113 gT$^2$. We also made computations by using a mixed Euler-Lagrangian scheme under the assumption of potential flow. The numerical results show good agreements with tank measurements.

• Journal of the Korean Geotechnical Society
• /
• v.35 no.2
• /
• pp.37-51
• /
• 2019

This study focuses on the numerical simulations of the cone penetration tests in a sand ground. The mechanical responses of sand were described using the modified Mohr Coulomb plasticity model based on the critical state soil mechanics. In the plasticity model, the dilatancy angle was not a constant, but a function of the distance to the critical state line from the current state of void ratio and mean effective stress. To simulate cone penetration tests numerically, this study relied on Lagrangian finite element method under the axisymmetric condition. To enable penetration of the cone penetrometer without tearing elements along the symmetric axis, the penetration guide concept was adopted in this study. The results of numerical simulations on the calibration chamber cone penetration tests had good agreement with the experimental results.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.26 no.5
• /
• pp.35-43
• /
• 2022

Precise modeling and analysis of closed bomb test(CBT) combustion using solid propellants was performed. The fluid structure interaction(FSI) method was implemented to analyze the gas and solid phases at the same time. The Eulerian analysis method was applied for the gas phase and grain combustion, and the Lagrangian analysis method was implemented for the grain movement. The interaction between the solid phase grains and the combustion gas was fully coupled through the source term. The volume of fluid(VOF) method was used to simulate the burning distance of the grain and the movement of the combustion surface. The force acting on the grain was comprised of the pressure and gravity acting on the grain burning surface, and the grain burning rate and grain movement speed were considered in the velocity term of the VOF. The combustion analysis was performed for both one and three grains, and fairly compared with the experiments. The acoustic field during grain combustion due to pressure fluctuations was also analyzed.

• Journal of the Korean Geotechnical Society
• /
• v.32 no.10
• /
• pp.17-29
• /
• 2016

A debris flow analysis model has been developed to simulate the erosion and entrainment of soil layer. Special attention is given to the model which represents strength softening behaviour of soil layer due to velocity of deformation. The 3D FE analysis by Coupled Eulerian-Lagrangian (CEL) model is conducted to simulate the debris flow. The model is validated using published data on laboratory experiment (Mangeny et al., 2010). It has been definitely shown that proposed model is in good agreement with the results of laboratory data. Futhermore, the FE analysis is conducted to ensure capability of simulating the real scale debris flow. The result of Ramian watershed, Korea shows that the debris flow has increased the volume and speed and it is in good agreement with field investigation. Based on this, it is confirmed that proposed model shows good agreement of the behavior of the actual and analytical debris flow.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.17 no.4
• /
• pp.375-387
• /
• 2004

A new finite element model will be presented to analyze the nonlinear behavior of RC beams and slabs strengthened by a patch repair. The numerical approach is based on the p-version degenerate shell element including theory of anisotropic laminated composites, theory of materially and geometrically nonlinear plates. In the nonlinear formulation of this model, the total Lagrangian formulation is adopted with large deflections and moderate rotations being accounted for in the sense of von Karman hypothesis. The material model is based on hardening rule, crushing condition, plate-end debonding strength model and so on. The Gauss-Lobatto numerical quadrature is applied to calculate the stresses at the nodal points instead of Gauss points. The validity of the proposed p-version nonlinear finite element model is demonstrated through the load-deflection curves, the ultimate loads, and the failure modes of RC beams or slabs bonded with steel plates or FRP plates compared with available result of experiment and other numerical methods.

• 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.

• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.2
• /
• pp.86-97
• /
• 1993

A numerical method is developed for the nonlinear motion of two-dimensional wedges and axisymmetric-forced-heaving motion using Semi-Largrangian scheme under assumption of potential flows. In two-dimensional-problem Cauchy's integral theorem is applied to calculate the complex potential and its time derivative along boundary. In three-dimensional-problem Rankine ring sources are used in a Green's theorem boundary integral formulation to salve the field equation. The solution is stepped forward numerically in time by integrating the exact kinematic and dynamic free-surface boundary condition. Numerical computations are made for the entry of a wedge with a constant velocity and for the forced harmonic heaving motion from rest. The problem of the entry of wedge compared with the calculated results of Champan[4] and Kim[11]. By Fourier transform of forces in time domain, added mass coefficient, damping coefficient, second harmonic forces are obtained and compared with Yamashita's experiment[5].