• Bulletin of the Society of Naval Architects of Korea
• /
• v.27 no.3
• /
• pp.107-116
• /
• 1990

The liquid sloshing in an elastic tank is a fluid-structure interaction problem. It requires nonlinear analysis to solve the complicated physics involved in the large interaction of fluid-structure, the variation of dynamic characteristics of structure due to hydrodynamic loading, and the distorsion of fluid flow due to structural vibration. In this paper a Lagrangian FEM is introduced to analyze the liquid sloshing in an elastic tank assuming that the elastic wall is one degree of freedom rigid wall. Numerical integration is performed using an implicit-explicit algorithm, which is formed by mixing the predictor-corrector method and the Runge-Kutta 4th order method. The influence of dynamic characteristics of the sloshing tank on the fluid flow is discussed. The numerical method is also applied for the simulation of the wall generated wave in the tank.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.6 no.3
• /
• pp.139-149
• /
• 2002

An explicit direct time integration method based solution algorithm is proposed to predict dynamic postbuckling response of thin plates. Based on the von Karman's plate equations and Marquerre's shallow shell theory, a rectangular plate finite element is formulated and utilized in this study. The element formulation takes into account geometrical nonlinearity and initial deflection of plates. The solution algorithm employs the central difference method. Using the computer program developed by the authors, dynamic postbuckling behavior of elastic thin plates under impact loading is investigated by considering the time variation of load and load duration. The efficiency of the proposed solution algorithm is examined through illustrative numerical examples.

• Journal of Korean Society of Water and Wastewater
• /
• v.31 no.6
• /
• pp.511-519
• /
• 2017

Multi-regional water supply system, which installed for supplying multiple water demands, is characterized by large-sized, long-distance, tree-type layout. This system is vulnerable to long-standing service interruption when a pipe breaks is occurred. In this study, a numerical method is proposed to calculate drainage time that directly affects time of service interruption. To begin with, governing equations are formulated to embed the delayed drainage effect by the friction loss, and to resolve complicated connection of pipelines, which are derived from the continuity and energy equations. The nonlinear hydraulic equations are solved by using explicit time integration method and the Newton-Raphson method. The developed model is verified by comparing the result with analytical solution. Furthermore, the model's applicability is validated by the examples of pipelines in serial, in parallel, and complex layout. Finally, the model is utilized to suggest an appropriate actions to reduce the deviation of draining time in the C transmission line of the B multi-regional water supply system.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 1995.10a
• /
• pp.223-226
• /
• 1995

A three-dimensional elastic-plastic finite element analysis using the explicit time integration method has been performed for the characterization of theimpact forming machines. The block upsetting using a forging hammer has been analyzed. The effects of machine type, work capacity of equipment and the mass ratio in an anvil-type hammer have been studied through the analysis.

• Wind and Structures
• /
• v.1 no.2
• /
• pp.127-144
• /
• 1998

In this paper, the wake structures behind a square cylinder at the Reynolds number of 22,000 are simulated using the large eddy simulation, and the main features of the wake structure associated with unsteady vortex-shedding are investigated. The Smagorinsky model is used for parametrization of the subgrid scales. The finite element method with isoparametric linear elements is employed in the computations. Unsteady computations are performed using the explicit method with streamline upwind scheme for the advection term. The time integration incorporates a subcycling strategy. No-slip condition is enforced on the wall surface. A comparative study between two-and three-dimensional computations puts a stress on the three-dimensional effects in turbulent flow simulations. Simulated three-dimensional wake structures are compared with numerical and experimental results reported by other researchers. The results include time-averaged, phase-averaged flow fields and numerically visualized vortex-shedding pattern using streaklines. The results show that dynamics of the vortex-shedding phenomenon are numerically well reproduced using the present method of finite element implementation of large eddy simulation.

• Journal of Ocean Engineering and Technology
• /
• v.33 no.5
• /
• pp.400-410
• /
• 2019

In this study, the numerical code for the 3D nonlinear dynamic analysis of an SLWR (Steel Lazy Wave Riser) was developed using the lumped mass line model in a FORTRAN environment. Because the lumped mass line model is an explicit method, there is no matrix operation. Thus, the numerical algorithm is simple and fast. In the lumped mass line model, the equations of motion for the riser were derived by applying the various forces acting on each node of the line. The applied forces at the node of the riser consisted of the tension, shear force due to the bending moment, gravitational force, buoyancy force, riser/ground contact force, and hydrodynamic force based on the Morison equation. Time integration was carried out using a Runge-Kutta fourth-order method, which is known to be stable and accurate. To validate the accuracy of the developed numerical code, simulations using the commercial software OrcaFlex were carried out simultaneously and compared with the results of the developed numerical code. To understand the nonlinear dynamic characteristics of an SLWR, dynamic simulations of SLWRs excited at the hang-off point and of SLWRs in regular waves were carried out. From the results of these dynamic simulations, the displacements at the maximum bending moments at important points of the design, like the hang-off point, sagging point, hogging points, and touch-down point, were observed and analyzed.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.1
• /
• pp.104-112
• /
• 1998

The mold filling process has been a central issue in the development of numerical methods to solve the casting processes. A mold filling which is inherently transient free surface fluid flow, is important because the quality of casting highly depends on such phenomenon, Most of the existing numerical schemes to solve mold filling process have severe limitations in time step restrictions or Courant criteria since explicit time integration is used. Therefore, a large computation time is required to analyze casting processes. In this study, the well known SOLA-VOF method has been modified implicitly to simulate the mold filling process. Solutions to example filling problems show that the proposed method is more efficient in computation time than the original SOLA -VOF method.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.3
• /
• pp.576-585
• /
• 2002

A three-dimensional Lagrangian explicit time-integration finite element code for analyzing the dynamic impact phenomena was developed. It uses four node tetrahedral elements. In order to consider the effects of strain rate hardening, strain hardening and thermal softening, which are frequently observed in high-velocity deformation phenomena, Johnson-Cook model is used as constitutive model. For more accurate and robust contact force computation, the defense node contact algorithm was adopted and implemented. In order to evaluate the performance of the newly developed three-dimensional hydrocode NET3D, numerical simulations of the oblique impact of mild steel plate by mild steel sphere were carried out. Ballistic limit about various oblique angle between 0 degree and 80 degree was estimated through a series of simulations with different initial velocities of sphere. Element eroding by equivalent plastic strain was applied to mild steel spheres and targets. Ballistic limits and fracture characteristics obtained from simulation were compared with experimental results conducted by Finnegan et al. From numerical studies, the following conclusions were reached. (1) Simulations could successfully reproduce the key features observed in experiment such as tensile failure termed "disking"at normal impacts and outwards bending of partially formed plus segments termed "hinge-mode"at oblique impacts. (2) Simulation results fur 60 degrees oblique impact at 0.70 km/s and 0.91 km/s were compared with experimental results and Eulerian hydrocode CTH simulation results. The Lagrangian code NET3D is superior to Eulerian code CTH in the computational accuracy. Agreement with the experimentally obtained final deformed cross-sections of the projectile is excellent. (3) Agreement with the experimental ballistic limit data, particularly at the high-obliquity impacts, is reasonably good. (4) The simulation result is not very sensitive to eroding condition but slightly influenced by friction coefficient.

• 한국전산유체공학회:학술대회논문집
• /
• 1995.10a
• /
• pp.64-69
• /
• 1995

A computer code for solving the Reynolds averaged full Navier-Stokes equations has bent developed for analysis of gas and steam turbine cascade flows with the option of using one of two types of turbulence model. One is the Baldwin-Lomax model and the other is standard $k-{\varepsilon}$ model. The numerical integration is based on the explicit four stage Runge-Kutta scheme and finite volume method. To be verified, the resulting code is applied to VKI turbine cascade and compared with the previous experimental results. Finally, the flowfield around a steam turbine cascade is analyzed. Comparisons with experimental data show that present numerical scheme is an accurate Navier-Stokes solver and can give very good predictions for both gas and steam turbine cascade flow.

• Structural Engineering and Mechanics
• /
• v.22 no.1
• /
• pp.107-130
• /
• 2006

The paper deals with the numerical solution of the dynamic problem of masonry structures. Masonry is modelled as a non-linear elastic material with zero tensile strength and infinite compressive strength. Due to the non-linearity of the adopted constitutive equation, the equations of the motion must be integrated directly. In particular, we apply the Newmark or the Hilber-Hughes-Taylor methods implemented in code NOSA to perform the time integration of the system of ordinary differential equations obtained from discretising the structure into finite elements. Moreover, with the aim of evaluating the effectiveness of these two methods, some dynamic problems, whose explicit solutions are known, have been solved numerically. Comparisons between the exact solutions and the corresponding approximate solutions obtained via the Newmark and Hilber-Hughes-Taylor methods show that in the cases under consideration both numerical methods yield satisfactory results.