• Bulletin of the Society of Naval Architects of Korea
• /
• v.27 no.2
• /
• pp.21-30
• /
• 1990

Theoretical and experimental techniques to analyze the two-dimensional liquid motion in a tank are discussed. A Lagrangian FEM with a velocity correction procedure is introduced to describe incompressible free surface fluid flow. A mesh rezoning technique is used to prevent strong distortion of finite elements in the Lagrangian description. Model test technique for sloshing tank is developed using a hydraulic type bench tester. The influence of the variation in the exciting frequency and amplitude are observed for various fill depths. The results of theoretical calculations are compared with those of experiments.

• Journal of Korea Water Resources Association
• /
• v.52 no.12
• /
• pp.1047-1055
• /
• 2019

Various types of numerical modeling techniques are used to predict the behavior of pollutants under various water environmental conditions in the event of a water pollutant accident. Among them, a hydraulic model that can consider water flow characteristics is the most basic and very important. The process of evaluating whether the hydraulic model accurately predicts the applied river characteristics is very important. In the verification of the modeling result, the measuring data are often used in the river. Currently, ADCP and FlowTrackers are widely used to measure the flow velocity of rivers. However, ADCP is not accurate when the depth is less than 0.6 m and also when the ratio of irreversibility near the surface is more than 40%. Futhermore, FlowTracker has a limitation in measuring at high depth and high velocity due to the direct measurement method in rivers. Simuation results, which are validated by these methods, are not reliable for low depth conditions of low flowrate and high velocity conditions of high flowrate. In this study, Lagrangian GPS floaters which measures physical quantity of water according to particle movement is used without the conventional method measured by Eulerian technique. The verification method of the model results was studied by comparing the simulation results of the hydraulic model with the velocities measured using the GPS floaters. When comparing the traveling distance of the GPS floaters with the traveling distance of the LPT simulations, the average error rate was 13.6% on distances, and the average error rate was 3.2% on velocities except for the stagnant section. Therefore, GPS floaters can be used for a correction and verification method of hydraulic model simulations.

• Coupled systems mechanics
• /
• v.7 no.5
• /
• pp.583-610
• /
• 2018

This paper discusses the issues associated with modeling frictional contact between solid bodies undergoing large deformations. The most common model for friction on contact interfaces in solid mechanics is the Coulomb friction model, in which two distinct responses are possible: stick and slip. Handling the transition between these two phases computationally has been a source of algorithmic instability, lack of convergence and non-unique solutions, particularly in the presence of large deformations. Most computational models for frictional contact have used penalty or updated Lagrangian approaches to enforce frictional contact conditions. These two approaches, however, present some computational challenges due to conditioning issues in penalty-type implementations and the iterative nature of the updated Lagrangian formulation, which, particularly in large simulations, may lead to relatively slow convergence. Alternatively, a plasticity-inspired implementation of frictional contact has been shown to handle the stick-slip conditions in a local, algorithmically efficient manner that substantially reduces computational cost and successfully avoids the issues of instability and lack of convergence often reported with other methods (Laursen and Simo 1993). The formulation of this approach, however, has been limited to the small deformations realm, a fact that severely limited its application to contact problems where large deformations are expected. In this paper, we present an algorithmically consistent formulation of this method that preserves its key advantages, while extending its application to the realm of large-deformation contact problems. We show that the method produces results similar to the augmented Lagrangian formulation at a reduced computational cost.

• Journal of the Korean Geotechnical Society
• /
• v.16 no.4
• /
• pp.191-199
• /
• 2000

In this research, the finite element analysis of piezocone penetration and dissipation tests has been conducted using the anisotropic elastoplastic-viscoplastic bounding surface model, virtual work equation, and theory of mixtures formulated in the Up[dated Lagrangian reference frame for the large deformation and finite strain nature of piezocone penetration. The formulated equations have been implemented into a finite element program. The cone resistance, excess pore water pressure, and dissipation of excess pore water pressure from the finite element analysis have been compared and investigated. An effective simulation could be performed with the use of the anisotropic and viscous soil model. The finite element formulations and the results are described in part 'I' and part 'II' respectively.

• Journal of Korean Society for Atmospheric Environment
• /
• v.15 no.6
• /
• pp.757-765
• /
• 1999

Lagrangian particle dispersion model(LPDM) is an effective tool to calculate the dispersion from a point source since it dose not induce numerical diffusion errors in solving the pollutant dispersion equation. Fictitious particles are released to the atmosphere from the emission source and they are then transported by the mean velocity and diffused by the turbulent eddy motion in the LPDM. The concentration distribution from the dispersed particles in the calculation domain are finally estimated by applying a particle count method or a Gaussian kernel method. The two methods for calculating concentration profiles were compared each other and tested against the analytic solution and the tracer experiment to find the strength and weakness of each method and to choose computationally time saving method for the LPDM. The calculated concentrations from the particle count method was heavily dependent on the number of the particles released at the emission source. It requires lots fo particle emission to reach the converged concentration field. And resulting concentrations were also dependent on the size of numerical grid. The concentration field by the Gaussian kernel method, however, converged with a low particle emission rate at the source and was in good agreement with the analytic solution and the tracer experiment. The results showed that Gaussian kernel method was more effective method to calculate the concentrations in the LPDM.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.19 no.1
• /
• pp.81-96
• /
• 2007

Nonlinear shoaling characteristics over surf zone are numerically investigated based on spatially averaged NavierStokes equation. We also test the validity of gradient model for turbulent stresses due to wave breaking using the data acquainted during SUPERTANK LABORATORY DATA COLLECTION PROJECT(Krauss et al., 1992). It turns out that the characteristics length scale of breaking induced current is not negligible, which firmly stands against ever popular gradient model, ${\kappa}-{\varepsilon}$ model, but favors Large Eddy Simulation with finer grid. Based on these observations, we model the residual stress of spatially averaged NavierStokes equation after Lagrangian Dynamic Smagorinsky(Meneveau et al., 1996). We numerically integrate newly proposed wave equations using SPH with Gaussian kernel function. Severely deformed water surface profile, free falling water particle, queuing splash after landing of water particle on the free surface and wave finger due to structured vortex on rear side of wave crest(Narayanaswamy and Dalrymple, 2002) are successfully duplicated in the numerical simulation of wave propagation over uniform slope beach, which so far have been regarded very difficult features to mimic in the computational fluid mechanics.

• Journal of ILASS-Korea
• /
• v.15 no.4
• /
• pp.170-176
• /
• 2010

This paper describes the grid-size dependency of the conventional Eulerian-Lagrangian method to spray characteristics such as spray penetration and SMD in modeling DME sprays. In addition, the reduction of the grid-size dependency of the present Gas-jet model was investigated. The calculations were performed using the KIVA code and the calculated results were compared to those of experimental result. The results showed that the conventional Eulerian-Laglangian model predicts shorter spray penetration for large cell because of inaccurate calculation of momentum exchange between liquid and gas phase. However, it was shown that the gas-jet model reduced grid-size dependency to spray penetration by calculating relative velocity between liquid and ambient gas based on gas jet velocity.

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.20 no.3
• /
• pp.321-338
• /
• 2022

The transport of radionuclides at oceanic scales can be assessed using a Lagrangian model. In this review an application of such a model to the Atlantic, Indian and Pacific oceans is described. The transport model, which is fed with water currents provided by global ocean circulation models, includes advection by three-dimensional currents, turbulent mixing, radioactive decay and adsorption/release of radionuclides between water and bed sediments. Adsorption/release processes are described by means of a dynamic model based upon kinetic transfer coefficients. A stochastic method is used to solve turbulent mixing, decay and water/sediment interactions. The main results of these oceanic radionuclide transport studies are summarized in this paper. Particularly, the potential leakage of ¹³⁷Cs from dumped nuclear wastes in the north Atlantic region was studied. Furthermore, hypothetical accidents, similar in magnitude to the Fukushima accident, were simulated for nuclear power plants located around the Indian Ocean coastlines. Finally, the transport of radionuclides resulting from the release of stored water, which was used to cool reactors after the Fukushima accident, was analyzed in the Pacific Ocean.

• Journal of Soil and Groundwater Environment
• /
• v.11 no.1
• /
• pp.37-44
• /
• 2006

This study summarizes advantages and disadvantages of numerical methods and compares ELLAM and LEZOOMPC to develop an efficient numerical modeling technique on contaminant transport. Eulerian-Lagrangian method and Eulerian method are commonly used numerical techniques. However Eulerian-Lagrangian method does not conserve mass globally and fails to treat boundary in a straightforward manner. Also, Eulerian method has restrictions on the size of Courant number and mesh Peclet number because of time truncation error. ELLAM (Eulerian Lagrangian Localized Adjoint Method) which has been popularly used for past 10 years in numerical modeling, is known for overcoming these numerical problems of Eulerian-Lagrangian method and Eulerian method. However, this study investigates advantages and disadvantages of ELLAM and suggests a change for the better. To figure out the disadvantages of ELLAM, the results of ELLAM, LEZOOMPC (Lagrangian-Eulerian ZOOMing Peak and valley Capturing), and visual MODFLOW are compared for four examples having different mesh Peclet numbers. The result of ELLAM generates numerical oscillation at infinite of mesh Peclet number, but that of LEZOOMPC yields accurate simulations. The simulation results suggest that the numerical error of ELLAM could be alleviated by adopting some schemes in LEZOOMPC. In other words, the numerical model which combines ELLAM with backward particle tracking, forward particle tracking, adaptively local zooming, and peak/valley capturing of LEZOOMPC can be developed for not only overcoming the numerical error of ELLAM, but also keeping the numerical advantage of ELLAM.

• Journal of the Institute of Electronics Engineers of Korea SP
• /
• v.42 no.6
• /
• pp.63-72
• /
• 2005

In this paper, we propose the bit-rate control algorithm for the block based image compression like H.263, H.263+ or MPEG-4. The proposed algorithm is designed to identify the quantization parameter set through the Lagrangian optimization technique based on the well-known linear source model. We set the Lagrangian cost function with the rates and distortion calculated from the linear source model. We calculate the quantization parameter set using the Vitervi algorithm to solve the Lagrangian optimization problem considering the Dquant method of H.263 and MPEG-4. The proposed algorithm improves the video quality by up to 1.5 dB compared with the TMN8 scheme, and is more effective in the video sources with dynamic activities than the consistent quality approaches.