• Proceedings of the Korean Society of Soil and Groundwater Environment Conference
• /
• 2004.04a
• /
• pp.57-62
• /
• 2004

Streamline simulation researches have been extensively accomplished due to the swiftness of computation and the reduction of numerical dispersion. In this study, we developed a streamline simulation model using a semianalytical solution of ID transport equation. To validate accuracy of the developed model, we compared simulation results of contaminant transport, which were acquired by streamline simulation models using an analytical solution, a numerical solution, and a semianalytical solution. The developed model using the semianalytical solution matched well with the model using an analytical solution. However, streamline simulation model using a numerical solution showed numerical dispersion. For an advection-dominant flow, there was little difference in the simulation results between the developed model and tile analytical model, but the differences between the analytical model and the numerical model were cleary shown. From the comparison of computing time we know that the streamline simulation using the semianalytical solution is 2-60 times as fast as the streamline simulation using the numerical solution.

• Journal of Korea Water Resources Association
• /
• v.30 no.6
• /
• pp.753-759
• /
• 1997

A new numerical model was doveloped to simulate density-dependent ground water flow and solute transport. Accuracy of a numerical model depends upon how well it simulates advection dominant situations because numerical oscillations can spoil solutions for these situations. Nonlinear oscillation-absorption finite element method. based on the variational principle, was employed. Unlike previous numerical models, this model can easily be expanded for more complex situations. Accuracy of the model is evaluated by comparing with analytical solutions and results of other numerical model.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.2
• /
• pp.467-474
• /
• 1993

Viscous flow and heat transfer phenomena in an axisymmetric cylinder which models a diesel engine have been numerically studied. In order to search for a way to minimize numerical diffusion, the effectiveness and the appropriateness of two selected numerical schemes for convective terms in the governing equations have been tested. They are Linear Upwind Difference Scheme and Hybrid Scheme. Using a standard k-.epsilon. turbulence model, the calculation has been carried out basically up to 180.deg. of crank angle. As a result, it was shown from comparison with previous experimental data that Linear Upwind Difference Scheme is less influenced than Hybrid Scheme by the numerical diffusion and it was suggested that these effects of numerical diffusion can be more significant than those due to turbulence modeling.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.14 no.4
• /
• pp.14-20
• /
• 2015

Molding technologies for plastic microstructures have been extensively investigated during the last two decades, and theoretical and numerical studies on the micro molding process have provided efficient tools for the development of such molding technologies. In this paper, we present a review of numerical simulation methods for the micro molding process. Basic models for a description of the material property, governing equations of the flow and heat transfer during the molding process, and numerical methods will be described. Particularly, numerical simulations for micro injection molding and hot embossing processes will be presented, and their main features noted and compared to those for conventional molding processes.

• Structural Engineering and Mechanics
• /
• v.67 no.5
• /
• pp.505-516
• /
• 2018

In this paper, the effects of particle size and model scale of concrete has been investigated on the failure mechanism of PFC2D numerical models under uniaxial compressive test. For this purpose, rectangular models with same particle sizes and different model dimensions, i.e., $3mm{\times}6mm$, $6mm{\times}12mm$, $12mm{\times}24mm$, $25mm{\times}50mm$ and $54mm{\times}108mm$, were prepared. Also rectangular models with dimension of $54mm{\times}108mm$ and different particle sizes, i.e., 0.27 mm, 0.47 mm, 0.67 mm, 0.87 mm, 1.07 mm, 1.87 mm and 2.27 mm were simulated using PFC2D and tested under uniaxial compressive test. Concurrent with uniaxial test, direct shear test was performed on the numerical models. Dimension of the models were $75{\times}100mm$. Two narrow bands of particles with dimension of $37.5mm{\times}20mm$ were removed from upper and lower of the model to supply the shear test condition. The particle sizes in the models were 0.47 mm, 0.57 mm, 0.67 mm and 0.77 mm. The result shows that failure pattern was affected by model scale and particle size. The uniaxial compressive strength and shear strength were increased by increasing the model scale and particle size.

• Earthquakes and Structures
• /
• v.12 no.1
• /
• pp.67-77
• /
• 2017

In this paper, it is aimed to present a comparative study about the structural behavior of tall buildings consisting of different type of materials such as concrete, steel or timber using finite element analyses and experimental measurements on shaking table. For this purpose, two 1/60 scaled 28 and 30-stories wooden building models with $40{\times}40cm$ and $35{\times}35cm$ ground/floor area and 1.45 m-1.55 m total height are built in laboratory condition. Considering the frequency range, mode shapes, maximum displacements and relative story drifts for structural models as well as acceleration, displacement and weight limits for shaking table, to obtain the typical building response as soon as possible, balsa is selected as a material property, and additional masses are bonded to some floors. Finite element models of the building models are constituted in SAP2000 program. According to the main purposes of earthquake resistant design, three different earthquake records are used to simulate the weak, medium and strong ground motions. The displacement and acceleration time-histories are obtained for all earthquake records at the top of building models. To validate the numerical results, shaking table tests are performed. The selected earthquake records are applied to first mode (lateral) direction, and the responses are recorded by sensitive accelerometers. Comparisons between the numerical and experimental results show that shaking table tests are enough to identify the structural response of wooden buildings. Considering 20%, 10% and 5% damping rations, differences are obtained within the range 4.03-26.16%, 3.91-65.51% and 6.31-66.49% for acceleration, velocity and displacements in Model-1, respectively. Also, these differences are obtained as 0.49-31.15%, 6.03-6.66% and 16.97-66.41% for Model-2, respectively. It is thought that these differences are caused by anisotropic structural characteristic of the material due to changes in directions parallel and perpendicular to fibers, and should be minimized using the model updating procedure.

• Geomechanics and Engineering
• /
• v.29 no.5
• /
• pp.567-582
• /
• 2022

In the present study, a series of physical experiments and numerical simulations were conducted to investigate the effects of mode I and mixed-mode I/II cracks on the fracture modes and stability of roadway tunnel models. The experiments and simulations incorporated different inclination angle flaws under both static and dynamic loads. The quasi-static and dynamic testing were conducted by using an electro-hydraulic servo control device and drop weight impact system (DWIS), and the failure process was simulated by using rock failure process analysis (RFPA) and AUTODYN software. The stress intensity factor was also calculated to evaluate the stability of the flawed roadway tunnel models by using ABAQUS software. According to comparisons between the test and numerical results, it is observed that for flawed roadways with a single radical crack and inclination angle of 45°, the static and dynamic stability are the lowest relative to other angles of fractured rock masses. For mixed-mode I/II cracks in flawed roadway tunnel models under dynamic loading, a wing crack is produced and the pre-existing cracks increase the stress concentration factor in the right part of the specimen, but this factor will not be larger than the maximum principal stress region in the roadway tunnel models. Additionally, damage to the sidewalls will be involved in the flawed roadway tunnel models under static loads.

• Journal of Ocean Engineering and Technology
• /
• v.38 no.2
• /
• pp.64-73
• /
• 2024

Oil spills pose significant threats to marine ecosystems, human health, socioeconomic aspects, and coastal communities. Accurate real-time predictions of oil slick transport along coastlines are paramount for quick preparedness and response efforts. This study used an open-source OpenOil numerical model to simulate the fate and trajectories of oil slicks released during the 2007 Hebei Spirit accident along the Korean coasts. Six combinations of input parameters, derived from a five-day met-ocean dataset incorporating various hydrodynamic, meteorological, and wave models, were investigated to determine the input variables that lead to the most reasonable results. The predictive performance of each combination was evaluated quantitatively by comparing the dimensions and matching rates between the simulated and observed oil slicks extracted from synthetic aperture radar (SAR) data on the ocean surface. The results show that the combination incorporating the Hybrid Coordinate Ocean Model (HYCOM) for hydrodynamic parameters exhibited more substantial agreement with the observed spill areas than Copernicus Marine Environment Monitoring Service (CMEMS), yielding up to 88% and 53% similarity, respectively, during a more than four-day oil transportation near Taean coasts. This study underscores the importance of integrating high-resolution met-ocean models into oil spill modeling efforts to enhance the predictive accuracy regarding oil spill dynamics and weathering processes.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.20 no.4
• /
• pp.363-371
• /
• 2008

The aim of this study is to compare the numerical results obtained by 2-D and 3-D models which are used to examine the wave field around a permeable submerged breakwater. At first, the numerical model, which is able to consider the flow through a porous medium with inertial, laminar and turbulent resistance terms and determine the eddy viscosity with LES turbulent model, is used and validated by comparing with existing experimental data. And then, the numerical test on the wave field around a permeable submerged breakwater is performed. It is revealed from the numerical results that, at the onshore side of the submerged breakwater, the wave height by 2-D analysis is higher than that by 3-D analysis. Also, the time-averaged mean flow around a submerged breakwater is discussed in detail.

• Geomechanics and Engineering
• /
• v.3 no.4
• /
• pp.291-321
• /
• 2011

The paper describes a simple numerical FLAC model that was developed to simulate the dynamic response of two instrumented reduced-scale model reinforced soil walls constructed on a 1-g shaking table. The models were 1 m high by 1.4 m wide by 2.4 m long and were constructed with a uniform size sand backfill, a polymeric geogrid reinforcement material with appropriately scaled stiffness, and a structural full-height rigid panel facing. The wall toe was constructed to simulate a perfectly hinged toe (i.e. toe allowed to rotate only) in one model and an idealized sliding toe (i.e. toe allowed to rotate and slide horizontally) in the other. Physical and numerical models were subjected to the same stepped amplitude sinusoidal base acceleration record. The material properties of the component materials (e.g. backfill and reinforcement) were determined from independent laboratory testing (reinforcement) and by back-fitting results of a numerical FLAC model for direct shear box testing to the corresponding physical test results. A simple elastic-plastic model with Mohr-Coulomb failure criterion for the sand was judged to give satisfactory agreement with measured wall results. The numerical results are also compared to closed-form solutions for reinforcement loads. In most cases predicted and closed-form solutions fall within the accuracy of measured loads based on ${\pm}1$ standard deviation applied to physical measurements. The paper summarizes important lessons learned and implications to the seismic design and performance of geosynthetic reinforced soil walls.