• Journal of the Korean Society of Safety
• /
• v.31 no.3
• /
• pp.8-15
• /
• 2016

The purpose of this study is to design a model with the structural stability so as not to lose the operational function due to structural plastic or fail of a sliding blast door by blast pressure to this aim, a numerical simulation was performed using full-size experiments and M&S (Modeling & Simulation) of the sliding blast door. The sliding blast door ($W3,000{\times}H2,500mm$) under the blast load is in the form of a sliding type 2-way metal grill, which was applied by a design blast pressure (reflected pressure $P_r$) of 17 bar. According to the experimental results of a real sliding blast door under blast load, the blast pressure reached the sliding blast door approximately 4.3 ms after the explosion and lasted about 4.0 ms thereafter. The maximum blast pressure($P_r$) was 347.7 psi (2,397.3 kPa), it is similar to the UFC 3-340-02 of Parameter(91 %). In addition, operation inspection that was conducted for the sliding blast door after real test showed a problem of losing the door opening function, which was because of the fail of the Reversal Bolt that was installed to prevent the shock due to rebound of the blast door from the blast pressure. According to the reproduction of the experiment through M&S by applying the blast pressure measurement value of the full-size experiments, the sliding blast door showed a similar result to the full-size experiment in that the reversal bolt part failed to lose the function. In addition, as the pressure is concentrated on the failed reversal bolt, the Principal Tensile Failure Stress was exceeded in only 1.25 ms after the explosion, and the reversal bolt completely failed after 5.4 ms. Based on the result of the failed reversal bolt through the full-size experiment and M&S, the shape and size of the bolts were changed to re-design the M&S and re-analyze the sliding blast door. According to the M&S re-analysis result when the reversal bolt was designed in a square of 25 mm ($625mm^2$), the maximum pressure that the reversal bolt receives showed 81% of the principal tensile failure stress of the material, in plastic stage before fail.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.24 no.5
• /
• pp.343-351
• /
• 2012

In this study, shallow-water design waves are calculated for the return period of 10, 20, 30, and 50 years, based on the extreme value analysis of the wave measurement data at Gangneung beach. These values are compared with the results of SWAN simulation with the boundary condition of the deep-water design waves of the corresponding return periods at the Gangneung sea area provided by the Fisheries Agency (FA, 1988) and Korea Ocean Research & Development Institute (KORDI, 2005). It is found that the shallow-water wave heights at Gangneung beach calculated by the deep-water design waves were significantly less than the observation data. As the return period becomes higher, the significant wave heights obtained by the extreme value analysis becomes higher than those computed by SWAN with the deep-water design waves of the corresponding return periods. KORDI computed the hindcast wave data from January 2004 to August 2008 by WAM with a finer-grid mesh system than those of previous studies. Comparisons of the wave hindcast results with the wave observation show that the reproducibility of the winter-season storm wave was considerably improved compared to the hindcast data from 1979 to 2003. Hereafter, it is necessary to carry out hindcast wave data for the years before 2004 using WAM with the finer-grid mesh system and to supplement the deep-water design wave.

• Journal of Korean Society of Steel Construction
• /
• v.20 no.1
• /
• pp.21-31
• /
• 2008

response surface method (RSM) is widely used to evaluate th e extremely smal probability of ocurence or toanalyze the reliability of very complicated structures. Althoug h Monte-Carlo Simulation (MCS) technique can evaluate any system, the procesing time of MCS dependson the reciprocal num ber of the probability of failure. The stochastic finite element method could solve thislimitation. However, it is limit ed to the specific program, in which the mean and coeficient o f random variables are programed by a perturbation or by a weigh ted integral method. Therefore, it is not aplicable when erequisite programing. In a few number of stage analyses, RSM can construct a regresion model from the response of the c omplicated structural system, thus, saving time and efort significantly. However, the acuracy of RSM depends on the dist ance of the axial points and on the linearity of the limit stat e functions. To improve the convergence in exact solution regardl es of the linearity limit of state functions, an improved adaptive response surface method is developed. The analyzed res ults have ben verified using linear and quadratic forms of response surface functions in two examples. As a result, the be st combination of the improved RSM techniques is determined and programed in a numerical code. The developed linear adapti ve weighted response surface method (LAW-RSM) shows the closest converged reliability indices, compared with quadratic form or non-adaptive or non-weighted RSMs.

• Journal of the Korean Institute of Gas
• /
• v.7 no.4 s.21
• /
• pp.44-52
• /
• 2003

In this work dynamic heat transfer in a CPFS (cable penetration fire stop) system built in the firewall of nuclear power plants is three-dimensionally investigated to develop a test-simulator that can be used to verify effectiveness of the sealant. Dynamic heat transfer in the fire stop system is formulated in a parabolic PDE (partial differential equation) subjected to a set of initial and boundary conditions. First, the PDE model is divided into two parts; one corresponding to heat transfer in the axial direction and the other corresponding to heat transfer on the vertical planes. The first PDE is converted to a series of ODEs (ordinary differential equations) at finite discrete axial points for applying the numerical method of SOR (successive over-relaxation) to the problem. The ODEs are solved by using an ODE solver In such manner, the axial heat flux can be calculated at least at the finite discrete points. After that, all the planes are separated into finite elements, where the time and spatial functions are assumed to be of orthogonal collocation state at each element. The initial condition of each finite element can be obtained from the above solution. The heat fluxes on the vertical planes are calculated by the Galerkin FEM (finite element method). The CPFS system was modeled, simulated, and analyzed here. The simulation results were illustrated in three-dimensional graphics. Through simulation, it was shown clearly that the temperature distribution was influenced very much by the number, position, and temperature of the cable stream, and that dynamic heat transfer through the cable stream was one of the most dominant factors, and that the feature of heat conduction could be understood as an unsteady-state process.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.20 no.6
• /
• pp.390-395
• /
• 2019

Hydrogymnastics so that sufficient exercise effect can be obtained using the resistance of water has a positive effect on patients who have to receive arthritis or rehabilitation treatment. However, the studies on the effect are insufficient, and the main cause of their effects has not been unclear yet. In this study, in order to identify the main cause of the effects of Hydrogymnastics, conducted Unsteady fluid flow simulation under the same conditions as the actual environment. The analysis model based on real hands, and the pressure fluctuation applied to the knuckle was analyzed by the computational fluid method. During the underwater movement of the hands, Various sizes of vortices were generated between fingers due to skin surface velocity and flow resistance. Pressure of about -500 Pa to +500 Pa is applied by the vortex flow. Also It was confirmed that the positive pressure and the negative pressure were continuously repeated up to maximum + 2000 Pa at the minimum of -2000 Pa at the portion where the direction was changed. Pressure fluctuations with a frequency of 20 Hz to 70 Hz were added continuously for each knuckle. These continuous pressure fluctuations provide a direct massage effect on the knuckles, an It is judged that the blood circulation at the relevant part is positively affected.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.3 no.3 s.10
• /
• pp.151-163
• /
• 2003

In this study, the algorithm of groundwater flow process was established for koreanized groundwater program development dealing with the geographic and geologic conditions of the aquifer have dynamic behaviour in groundwater flow system. All the input data settings of the 3-DFM model which is developed in this study are organized in Korean, and the model contains help function for each input data. Thus, it is designed to get detailed information about each input parameter when the mouse pointer is placed on the corresponding input parameter. This model also is designed to easily specify the geologic boundary condition for each stratum or initial head data in the work sheet. In addition, this model is designed to display boxes for input parameter writing for each analysis condition so that the setting for each parameter is not so complicated as existing MODFLOW is when steady and unsteady flow analysis are performed as well as the analysis for the characteristics of each stratum. Descriptions for input data are displayed on the right side of the window while the analysis results are displayed on the left side as well as the TXT file for this results is available to see. The model developed in this study is a numerical model using finite differential method, and the applicability of the model was examined by comparing and analyzing observed and simulated groundwater heads computed by the application of real recharge amount and the estimation of parameters. The 3-DFM model is applied in this study to Sehwa-ri, and Songdang-ri area, Jeju, Korea for analysis of groundwater flow system according to pumping, and obtained the results that the observed and computed groundwater head were almost in accordance with each other showing the range of 0.03 - 0.07 error percent. It is analyzed that the groundwater flow distributed evenly from Nopen-orum and Munseogi-orum to Wolang-bong, Yongnuni-orum, and Songja-bong through the computation of equipotentials and velocity vector using the analysis result of simulation which was performed before the pumping started in the study area. These analysis results show the accordance with MODFLOW's.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2003.05a
• /
• pp.91-93
• /
• 2003

A comprehensive numerical study is carried out to investigate for the understanding of the flow evolution and flame development in a supersonic combustor with normal injection of ncumally injecting hydrogen in airsupersonic flows. The formulation treats the complete conservation equations of mass, momentum, energy, and species concentration for a multi-component chemically reacting system. For the numerical simulation of supersonic combustion, multi-species Navier-Stokes equations and detailed chemistry of H2-Air is considered. It also accommodates a finite-rate chemical kinetics mechanism of hydrogen-air combustion GRI-Mech. 2.11[1], which consists of nine species and twenty-five reaction steps. Turbulence closure is achieved by means of a k-two-equation model (2). The governing equations are spatially discretized using a finite-volume approach, and temporally integrated by means of a second-order accurate implicit scheme (3-5).The supersonic combustor consists of a flat channel of 10 cm height and a fuel-injection slit of 0.1 cm width located at 10 cm downstream of the inlet. A cavity of 5 cm height and 20 cm width is installed at 15 cm downstream of the injection slit. A total of 936160 grids are used for the main-combustor flow passage, and 159161 grids for the cavity. The grids are clustered in the flow direction near the fuel injector and cavity, as well as in the vertical direction near the bottom wall. The no-slip and adiabatic conditions are assumed throughout the entire wall boundary. As a specific example, the inflow Mach number is assumed to be 3, and the temperature and pressure are 600 K and 0.1 MPa, respectively. Gaseous hydrogen at a temperature of 151.5 K is injected normal to the wall from a choked injector.A series of calculations were carried out by varying the fuel injection pressure from 0.5 to 1.5MPa. This amounts to changing the fuel mass flow rate or the overall equivalence ratio for different operating regimes. Figure 1 shows the instantaneous temperature fields in the supersonic combustor at four different conditions. The dark blue region represents the hot burned gases. At the fuel injection pressure of 0.5 MPa, the flame is stably anchored, but the flow field exhibits a high-amplitude oscillation. At the fuel injection pressure of 1.0 MPa, the Mach reflection occurs ahead of the injector. The interaction between the incoming air and the injection flow becomes much more complex, and the fuel/air mixing is strongly enhanced. The Mach reflection oscillates and results in a strong fluctuation in the combustor wall pressure. At the fuel injection pressure of 1.5MPa, the flow inside the combustor becomes nearly choked and the Mach reflection is displaced forward. The leading shock wave moves slowly toward the inlet, and eventually causes the combustor-upstart due to the thermal choking. The cavity appears to play a secondary role in driving the flow unsteadiness, in spite of its influence on the fuel/air mixing and flame evolution. Further investigation is necessary on this issue. The present study features detailed resolution of the flow and flame dynamics in the combustor, which was not typically available in most of the previous works. In particular, the oscillatory flow characteristics are captured at a scale sufficient to identify the underlying physical mechanisms. Much of the flow unsteadiness is not related to the cavity, but rather to the intrinsic unsteadiness in the flowfield, as also shown experimentally by Ben-Yakar et al. [6], The interactions between the unsteady flow and flame evolution may cause a large excursion of flow oscillation. The work appears to be the first of its kind in the numerical study of combustion oscillations in a supersonic combustor, although a similar phenomenon was previously reported experimentally. A more comprehensive discussion will be given in the final paper presented at the colloquium.

• Journal of Korean Society of Disaster and Security
• /
• v.14 no.2
• /
• pp.61-75
• /
• 2021

Recently, as the occurrence frequency of sudden floods due to climate change increased and the aging of the existing spillway, it is necessary to establish a plan to utilize an auxiliary spillway to minimize the flood damage of downstream rivers. Most studies have been conducted on the review of flow characteristics according to the operation of auxiliary spillway through the hydraulic experiments and numerical modeling. However, the studies on examination of flood damage in the downstream rivers and the stability of the revetment according to the operation of the auxiliary spillway were relatively insufficient in the literature. In this study, the stability of the revetment on the downstream river according to the outflow conditions of the existing and auxiliary spillway was examined by using 3D numerical model, FLOW-3D. The velocity, water surface elevation and shear stress results of FLOW-3D were compared with the permissible velocity and shear stress of design criteria. It was assumed the sluice gate was fully opened. As a result of numerical simulations of various auxiliary spillway operations during flood season, the single operation of the auxiliary spillway showed the reduction effect of maximum velocity and the water surface elevation compared with the single operation of the existing spillway. The stability of the revetment on downstream was satisfied under the condition of outflow less than 45% of the design flood discharge. However, the potential overtopping damage was confirmed in the case of exceeding the 45% of the design flood discharge. Therefore, the simultaneous operation with the existing spillway was important to ensure the stability on design flood discharge condition. As a result of examining the allocation ratio and the total allowable outflow, the reduction effect of maximum velocity was confirmed on the condition, where the amount of outflow on auxiliary spillway was more than that on existing spillway. It is because the flow of downstream rivers was concentrated in the center due to the outflow of existing spillway. The permissible velocity and shear stress were satisfied under the condition of less than 77% of the design flood discharge with simultaneous operation. It was found that the flood damage of downstream rivers can be minimized by setting the amount allocated to the auxiliary spillway to be larger than the amount allocated to the existing spillway for the total outflow with simultaneous operation condition. However, this study only reviewed the flow characteristics around the revetment according to the outflow of spillway under the full opening of the sluice gate condition. Therefore, the various sluice opening conditions and outflow scenarios will be asked to derive more efficient utilization of the auxiliary spillway in th future.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.26 no.4C
• /
• pp.247-254
• /
• 2006

Online testing method is one of the numerical experiment methods using experimental information for a numerical analysis directly. The method has an advantage in that analysis can be conducted without using an idealized mechanical model, because mechanical properties are updated from element test for a numerical analysis in real time. The online testing method has mainly been used for the geotechnical seismic engineering, whose major target is sand. A testing method that may be applied to a consolidation problem has recently been developed and laboratory and field verifications have been tried. Although related research thus far has mainly used a method to update average reaction for a numerical analysis by positioning an element tests at the center of a consolidation layer, a weakness that accuracy of the analysis can be impaired as the thickness of the consolidation layer becomes more thicker has been pointed out regarding the method. To clarify the effectiveness and possible analysis scope of the online testing method in relation to the consolidation problem, we need to review the results by applying experiment conditions that may completely exclude such a factor. This research reviewed the results of the online consolidation test in terms of reproduction of the consolidation settlement and the dissipation of excess pore water pressure of a clay specimen by comparing the results of an online consolidation test and a separated-type consolidation test carried out under the same conditions. As a result, the online consolidation test reproduced the change of compressibility according effective stress of clay without a huge contradiction. In terms of the dissipation rate of excess pore water pressure, however, the online consolidation test was a little faster. In conclusion, experiment procedure needs to improve in a direction that hydraulic conductivity can be updated in real time so as to more precisely predict the dissipation of excess pore water pressure. Further research or improvement should be carried out with regard to the consolidation settlement after the end of the dissipation of excess pore water pressure.

• Clean Technology
• /
• v.30 no.1
• /
• pp.37-46
• /
• 2024

Environmental problems caused by plastic waste have been continuously growing around the world, and plastic waste is increasing even faster after COVID-19. In particular, PP and PE account for more than half of all plastic production, and the amount of waste from these two materials is at a serious level. As a result, researchers are searching for an alternative method to plastic recycling, and plastic pyrolysis is one such alternative. In this paper, a numerical study was conducted on the pyrolysis behavior of non-condensable gas to predict the chemical reaction behavior of the pyrolysis gas. Based on gas products estimated from preceding literature, the behavior of non-condensable gas was analyzed according to temperature and residence time. Numerical analysis showed that as the temperature and residence time increased, the production of H₂ and heavy hydrocarbons increased through the conversion of the non-condensable gas, and at the same time, the CH₄ and C₆H₆ species decreased by participating in the reaction. In addition, analysis of the production rate showed that the decomposition reaction of C₂H₄ was the dominant reaction for H₂ generation. Also, it was found that more H₂ was produced by PE with higher C₂H₄ contents. As a future work, an experiment is needed to confirm how to increase the conversion rate of H₂ and carbon in plastics through the various operating conditions derived from this study's numerical analysis results.