• 한국연소학회:학술대회논문집
• /
• 2002.11a
• /
• pp.193-207
• /
• 2002

The present study is mainly motivated to investigate the vaporization, autoignition, and combustion of liquid fuel spray injected into high pressure environment. In order to represent these phenomena realistically, discrete droplet model (DDM) which simulates the spray using finite number of representative droplets was adopted for detailed consideration of the finite rate of uansport between liquid and gas phases. The Eulerian-Lagrangian formulation was used to analyze the two-phase interactions. The high pressure vaporization model was applied using the thermodynamic and phase equilibrium at droplet surface. The high pressure effect as well as high temperature effect was considered in the calculation of liquid and gas properties. The characteristics of spray in high pressure environment were explained by comparison with normal pressure case.

• Nuclear Engineering and Technology
• /
• v.51 no.5
• /
• pp.1261-1271
• /
• 2019

Droplet size and distribution are important parameters determining venturi scrubber performance. In this paper, we proposed physical models for a maximum stable droplet size prediction and upper limit log-normal (ULLN) distribution parameters. For the proposed maximum stable droplet size prediction model, a Eulerian-Lagrangian framework and a Reitz-Diwakar breakup model are solved simultaneously using CFD calculations to reflect the effect of multistage breakup and droplet acceleration. Then, two ULLN distribution parameters are suggested through best fitting the previously published experimental data. Results show that the proposed approach provides better predictions of maximum stable droplet diameter and Sauter mean diameter compared to existing simple empirical correlations including Boll, Nukiyama and Tanasawa. For more practical purpose, we developed the simple, one dimensional (1-D) calculation of Sauter mean diameter.

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

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.35 no.12
• /
• pp.1082-1088
• /
• 2007

Demand for high speed sea transportation modes has been increased dramatically last few decades. The WIG(Wing-in-ground effect) is considered as next generation maritime transportation system. In the structural design of high speed marine vessels, an estimation of water impact loads is essential. The dynamic structural responses of the WIG excited by the water impact loads may bring an important contribution to their damage process. The work presented in this paper is focused on the numerical simulation of the water impact on the WIG craft when it lands. It is aimed to study the structural responses of the WIG craft subjected to the water impact loads. The Arbitrary Lagrangian-Eulerian (ALE) finite element method is used to simulate the water impact of the WIG craft during a landing phase. A full 3D shell element is used to model the WIG craft in carbon composites, and a developed FE model is used to investigate the effect of the water impact loads on the structural responses of the WIG craft. In the analysis, two different landing scenarios are considered and their effects on the structural responses are investigated.

• Water for future
• /
• v.26 no.3
• /
• pp.137-148
• /
• 1993

A hybrid finite difference method for the longitudinal dispersion equation was developed. The method is based on combining the Holly-Preissmann scheme with the fifth-degree Hermite interpolating polynomial and the generalized Crank-Nicholson scheme. Longitudinal dispersion of an instantaneously-loaded pollutant source was simulated by the model and other characteristics-based numerical methods. Computational results were compared with the exact solution. The present method was free from wiggles regardless of the Courant number, and exactly reproduced the location of the peak concentration. Overall accuracy of the computation increased for smaller value of the weighting factor, $\theta$ of the model. Larger values of $\theta$ overestimated the peak concentration. Smaller Courant number gave better accuracy, in general, but the sensitivity was very low, especially when the value of $\theta$ was small. From comparisons with the hybrid method using the third-degree interpolating polynomial and with split-operator methods, the present method showed the best performance in reproducing the exact solution as the advection becomes more dominant.

• Journal of Korean Society of Water and Wastewater
• /
• v.27 no.2
• /
• pp.261-272
• /
• 2013

Algae boom (Red tide) in south coastal area of Korea has been appeared several times during a decade. If algae boom appears in the desalination plant, media filter and UF filter are clogged quickly, and the plant should be shutdown. In general, Algae can be removed from water by flotation better than by sedimentation, because of the low density of algal cell. The purpose of this study conducts the CFD simulation of DAF flotation basin to apply the design of the dissolved air flotation with ball filter in the Test Bed for SWRO desalination plant. In this study, Eulerian-Eulerian multiphase model was applied to simulate the behavior of air bubbles and seawater. Density difference model and gravity were used. But de-sludge process and mass transfer between air bubbles and seawater were ignored. Main parameter is hydraulic loading rate which is varied from 20 m/hr to 27.5 m/hr. Geometry of flotation basin were changed to improve the DAF performance. According to the result of this study, the increase of hydraulic loading rate causes that the flow in the separation basin is widely affected and the concentration of air is increased. The flow pattern in the contact zone of flotation basin is greatly affected by the location of nozzle header. When the nozzle header was installed not the bottom of the contact zone but the above, the opportunity of contact between influent and recycle flow was increased.

• Journal of Korean Society of Water and Wastewater
• /
• v.18 no.2
• /
• pp.201-207
• /
• 2004

Dissolved air flotation (DAF) is a solid-liquid separation process that uses fine rising bubbles to remove particles in water. Most of particle-bubble collision occurs in the DAF contact zone. This initial contact considered by the researchers to play a important role for DAF performance. It is hard to make up conceptual model through simple mass balance for estimating collision efficiency in the contact zone because coupled behavior of the solid-liquid-gas phase in DAF system is 90 complicate. In this study, 2-phase(gas-liquid) flow equations for the conservation of mass, momentum and turbulence quantities were solved using an Eulerian-Eulerian approach based on the assumption that very small particle is applied in the DAF system. For the modeling of turbulent 2-phase flow in the reactor, the standard $k-{\varepsilon}$ mode I(liquid phase) and zero-equation(gas phase) were used in CFD code because it is widely accepted and the coefficients for the model are well established. Particle-bubble collision efficiency was calculated using predicted turbulent energy dissipation rate and gas volume fraction. As the result of this study, the authors concluded that bubble size and recycle ratio play important role for flow pattern change in the reactor. Predicted collision efficiency using CFD showed good agreement with measured removal efficiency in the contact zone. Also, simulation results indicated that collision efficiency at 15% recycle ratio is higher than that of 10% and showed increasing tendency of the collision efficiency according to the decrease of the bubble size.

• Steel and Composite Structures
• /
• v.45 no.3
• /
• pp.389-408
• /
• 2022

Residual capacity is defined as the load carrying capacity of an RC column after undergoing severe damage. Evaluation of residual capacity of RC columns is necessary to avoid damage initiation in RC structures. The central aspect of the current research is to propose an empirical formula to estimate the residual capacity of RC columns after undergoing severe damage. This formula facilitates decision making of whether a replacement or a repair of the damaged column is adequate for further use. Available literature mainly focused on the simulation of explosion loads by using simplified pressure time histories to develop residual capacity of RC columns and rarely simulated the actual explosive. Therefore, there is a gap in the literature concerning general relation between blast damage of columns with different explosive loading conditions for a reliable and quick evaluation of column behavior subjected to blast loading. In this paper, the Arbitrary Lagrangian Eulerian (ALE) technique is implemented to simulate high fidelity blast pressure propagations. LS-DYNA software is utilized to solve the finite element (FE) model. The FE model is validated against the practical blast tests, and outcomes are in good agreement with test results. Multivariate linear regression (MLR) method is utilized to derive an analytical formula. The analytical formula predicts the residual capacity of RC columns as functions of structural element parameters. Based on intensive numerical simulation data, it is found that column depth, longitudinal reinforcement ratio, concrete strength and column width have significant effects on the residual axial load carrying capacity of reinforced concrete column under blast loads. Increasing column depth and longitudinal reinforcement ratio that provides better confinement to concrete are very effective in the residual capacity of RC column subjected to blast loads. Data obtained with this study can broaden the knowledge of structural response to blast and improve FE models to simulate the blast performance of concrete structures.

• Journal of Advanced Marine Engineering and Technology
• /
• v.25 no.6
• /
• pp.1308-1316
• /
• 2001

This paper described a numerical investigation performed to understand better the effects of flow parameters in an entrained flow combustor on the flow characteristics. The computational model was based on the gas phase Eulerian equations of mass, momentum and energy. The code was formulated with RNG $k-\varepsilon$ model for turbulent flow. The calculation parameters were the ratio of primary and secondary jet velocity and the height difference between primary and secondary jet As the secondary jet velocity increased, the upper recirculation 3one of the primary jet was strengthened. It was found that as the primary jet velocity increased, there was a critical jet Velocity at which the size of upper and lower recirculation zone was reversed.

• 한국연소학회:학술대회논문집
• /
• 2002.06a
• /
• pp.24-32
• /
• 2002

The autoignition and subsequent flame propagation of initially nonpremixed turbulent system have been numerically analyzed. The unsteady flamelet modeling based on the RIF (Representative Interactive Flamelet) concept has been employed to account for the influences of turbulence on these essentially transient combustion processes. In this RIF approach, the partially premixed burning, diffusive combustion and formation of pollutants(NOx, soot) can be consistently modeled by utilizing the comprehensive chemical mechanism. To treat the spatially distributed inhomogeneity of scalar dissipation rate, the multiple RIFs are employed in the framework of EPFM(Eulerian Particle Flamelet Model) approach. Computations are made for the various initial conditions of pressure, temperature, and fuel composition. The present turbulent combustion model reasonably well predicts the essential features of autoignition process in the transient gaseous fuel jets injected into high pressure and temperature environment.