• Journal of computational fluids engineering
• /
• v.11 no.1 s.32
• /
• pp.29-35
• /
• 2006

A code is developed to analyze a spherically symmetric underwater explosion. The arbitrary Lagrangian-Eulerian(ALE) Godunov scheme for two-phase flow is used to calculate numerical fluxes through moving control surfaces. For detonation gas of TNT and liquid water, the Jones-Wilkins-Lee(JWL) equation of states and the isentropic Tait relation are used respectively. It is suggested to use the Godunov variable to estimate the velocity of a material interface. The code is validated through comparisons with other results on the gas-water shock tube problem. It is shown that the code can handle generation of discontinuity and recovering of continuity in the normal velocity near the material interface during shock waves interact with the material interface. The developed code is applied to analyze a spherically symmetric underwater explosion. Repeated transmissions of shock waves are clearly captured. The calculated period and maximum radius of detonation gas bubble show good agreements with experimental and other numerical results.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.15 no.4
• /
• pp.65-72
• /
• 2011

Numerical simulation has been performed to investigate the bubbly flow in the liquid with vertical temperature gradient. The objective of this study is to establish an accurate numerical prediction program of gas-liquid two-phase flows in a vertical temperature gradient condition, whose mathematical model is formulated by the Eulerian-Lagrangian model. The present numerical results reveal the temperature mixing mechanism and the fluid dynamical characteristics induced by the bubbly flow in the liquid with stratified temperature. The effects of bubble radius, void fraction, and gas flow rate on bubble-driven convective flow are considered, too.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.1
• /
• pp.60-68
• /
• 2009

The steam reformer for hydrogen production from methane is studied by a numerical method. Langmuir- Hinshelwood model is incorporated for catalytic surface reactions, and the pseudo-homogeneous model is used to take into account local equilibrium phenomena between a catalyst and bulk gas. Dominant chemical reactions are Steam Reforming (SR) reaction, Water-Gas Shift (WGS) reaction, and Direct Steam Reforming (DSR) reaction. The numerical results are validated with experimental results at the same operating conditions. Using the validated code, parametric study has been numerically performed in view of the steam reformer performance. As increasing a wall temperature, the fuel conversion increases due to the high heat transfer rate. When Steam to Carbon Ratio (SCR) increases, the concentration of carbon monoxide decreases since WGS reaction becomes more active. When increasing Gas Hourly Space Velocity (GHSV), the fuel conversion decreases due to the heat transfer limitation and the low residence time. The reactor shape effects are also investigated. The length and radius of cylindrical reactors are changed at the same catalyst volume. The longer steam reformer is, the better steam reformer performs. However, system energy efficiency decreases due to the large pressure drop.

• Journal of Power System Engineering
• /
• v.10 no.4
• /
• pp.25-30
• /
• 2006

This paper represents the numerical analysis in unsteady state on the variation of diameter and environmental velocity of carbon heating source. In general heating system, the oil and sheath heater is widely used, but these systems have many problems. So, the heating source with carbon ingredient has been researched in many country about manufacture, thermal and electrical properties. In this research, the carbon heating source was studied through numerical analysis on several conditions of unsteady state, heat generation, diameter and environmental velocity. The temperature distributions at steady state are appeared as a non-proportional linear pattern with variations of environmental velocity due to the Nesselt number with convective heat flux is proportioned to 0.805 of Reynolds number. As the radius is increasing, the temperature distributions is appeared the minus tilt because of the environmental condition is cooling by constant temperature. So, the correlation equation between temperature at steady state and environmental velocity was obtained.

• International Journal of Highway Engineering
• /
• v.18 no.2
• /
• pp.51-60
• /
• 2016

PURPOSES : The applicability of the mechanics-based similarity concept (suggested by Feng et al.) for determining scaled variables, including length and load, via laboratory-scale tests and discrete element analysis, was evaluated. METHODS: Several studies on the similarity concept were reviewed. The exact scaling approach, a similarity concept described by Feng, was applied in order to determine an analytical solution of a free-falling ball. This solution can be considered one of the simplest conditions for discrete element analysis. RESULTS : The results revealed that 1) the exact scaling approach can be used to determine the scale of variables in laboratory tests and numerical analysis, 2) applying only a scale factor, via the exact scaling approach, is inadequate for the error-free replacement of small particles by large ones during discrete element analysis, 3) the level of continuity of flowable materials such as SCC and cement mortar seems to be an important criterion for evaluating the applicability of the similarity concept, and 4) additional conditions, such as the kinetics of particle, contact model, and geometry, must be taken into consideration to achieve the maximum radius of replacement particles during discrete element analysis. CONCLUSIONS : The concept of similarity is a convenient tool to evaluate the correspondence of scaled laboratory test or numerical analysis to physical condition. However, to achieve excellent correspondence, additional factors, such as the kinetics of particles, contact model, and geometry, must be taken into consideration.

• Journal of the Society of Naval Architects of Korea
• /
• v.58 no.2
• /
• pp.73-83
• /
• 2021

In this study, the numerical investigation of the flow around the SUBOFF submarine model is performed by using the Detached Eddy Simulation (DES) method which is developed based on the SST k-ω turbulence model. At the DES analysis level, complex vortical flows around the submarine model are caused mainly by the vortices due to the appendages and their interactions with the flows from the hull boundary layer and other appendages. The complexity and scale of the vortical flow obtained from the numerical simulations are highly dependent on the grid. The computed local flow properties of the submarine model are compared with the available experimental data showing a good agreement. The DES analysis more reasonably estimates the physical phenomena inherent in the experimental result in a low radius of the propeller plane where vortical flows smaller than the RANS scale are dominant.

• Journal of the Korean Geotechnical Society
• /
• v.36 no.11
• /
• pp.149-156
• /
• 2020

Permeation grouting effectively enhances soil strength and decreases permeability of soil; however, the flow of grout is heavily affected by anisotropy of hydraulic conductivity in layers. Therefore, this study investigates the effect of permeability anisotropy on the effective radius of horizontal permeation grout using computational fluid dynamics (CFD). We modeled the horizontal permeation grout flow as a two-phase viscous fluid flow in porous media, and the model incorporated the chemical diffusion and the viscosity variation due to hardening. The numerical simulation reveals that the permeability anisotropy shapes the grout bulb to be elliptic and the dissolution-driven diffusion causes a gradual change in grout pore saturation at the edge of the grout bulb. For the grout pore saturations of 10%, 50% and 90%, the horizontal and vertical radii of grout bulb are estimated when the horizontal-to-vertical permeability ratio varies from 0.01 to 100, and the predictive model equations are suggested. This result contributes to more efficient design of injection strategy in formation layers with permeability anisotropy.

• Computers and Concrete
• /
• v.10 no.2
• /
• pp.135-147
• /
• 2012

At mesoscale, concrete is considered as a three-phase composite material consisting of the aggregate particles, the cement matrix and the interfacial transition zone (ITZ). The reconstruction of the internal structures for concrete composites requires the identification of the boundary of the aggregate particles and the cement matrix using digital imaging technology followed by post-processing through MATLAB. A parameter study covers the subsection transformation, median filter, and open and close operation of the digital image sample to obtain the optimal parameter for performing the image processing technology. The subsection transformation is performed using a grey histogram of the digital image samples with a threshold value of [120, 210] followed by median filtering with a $16{\times}16$ square module based on the dimensions of the aggregate particles and their internal impurity. We then select a "disk" tectonic structure with a specific radius, which performs open and close operations on the images. The edges of the aggregate particles (similar to the original digital images) are obtained using the canny edge detection method. The finite element model at mesoscale can be established using the proposed image processing technology. The location of the crack determined through the numerical method is identical to the experimental result, and the load-displacement curve determined through the numerical method is in close agreement with the experimental results. Comparisons of the numerical and experimental results show that the proposed image processing technology is highly effective in reconstructing the internal structures of concrete composites.

• Journal of the Korean Institute of Intelligent Systems
• /
• v.16 no.6
• /
• pp.683-690
• /
• 2006

Reducing power consumption to extend network lifetime is one of the most important challenges in designing wireless sensor networks. One promising approach to conserving system energy is to keep only a minimal number of sensors active and put others into low-powered sleep mode, while the active sensors can maintain a connected covet set for the target area. The problem of computing such minimum working sensor set is NP-hard. In this paper, a centralized Voronoi tessellation (CVT) based approximate algorithm is proposed to construct the near optimal cover set. When sensor's communication radius is at least twice of its sensing radius, the covet set is connected at the same time; In case of sensor's communication radius is smaller than twice of its sensing radius, a connection scheme is proposed to calculate the assistant nodes needed for constructing the connectivity of the cover set. Finally, the performance of the proposed algorithm is evaluated through theoretical analysis and extensive numerical experiments. Experimental results show that the proposed algorithm outperforms the greedy algorithm in terms of the runtime and the size of the constructed connected cover set.

• Tunnel and Underground Space
• /
• v.32 no.5
• /
• pp.312-326
• /
• 2022

The generalized Hoek-Brown (GHB) criterion, which is recognized as one of the standard failure conditions for rock mass, is specialized for rock engineering applications and covers a wide range of rock mass conditions. Accordingly, many research efforts have been devoted to the incorporation of this criterion into the stability analysis of rock structures. In this study, the slip-line analysis method, which is a kind of elastoplastic analysis method, is combined with the GHB failure criterion to derive analytical equations that can easily calculate the plastic radius and stress distribution in the vicinity of the circular tunnel. In the process of derivation of related formulas, it is assumed that the behavior of rock mass after failure is perfectly plastic and the in-situ stress condition is hydrostatic. In the formulation, it is revealed that the plastic radius can be calculated analytically using the two respective tangential friction angles corresponding to the stress conditions at tunnel wall and elastic-plastic boundary. It is also shown that the plastic radius and stress distribution calculated using the derived analytical equations coincide with the results of Lee & Pietruszczak's numerical method published in 2008. In the latter part of this paper, the influence of the quality of the rock mass on the size of the plastic zone, the stress distribution, and the change of the tangential friction angle was investigated using the derived analytical equations.