• Applied Science and Convergence Technology
• /
• v.23 no.4
• /
• pp.179-186
• /
• 2014

Fluid model based numerical simulation was carried out for an inductively coupled plasma assisted sputter deposition system. Power absorption, electron temperature and density distribution was modeled with drift diffusion approximation. Effect of an electrically conducting substrate was analyzed and showed confined plasma below the substrate. Part of the plasma was leaked around the substrate edge. Comparison between the quasi-neutrality based compact model and Poisson equation resolved model showed more broadened profile in inductively coupled plasma power absorption than quasi-neutrality case, but very similar Ar ion number density profile. Electric potential was calculated to be in the range of 50 V between a Cr rod source and a conductive substrate. A new model including Cr sputtering by Ar+was developed and used in simulating Cr deposition process. Cr was modeled to be ionized by direct electron impact and showed narrower distribution than Ar ions.

• Journal of the Korean Vacuum Society
• /
• v.8 no.2
• /
• pp.83-92
• /
• 1999

Pumping performance of a disk-type molecular drag pump for a hybrid molecular pump is numerically analyzed by the direct simulation onte-Carlo method. The flows in pumping channels are three-dimensional (3D) in a molecular transition regime. The main difficulty in modeling a 3D case comes from the rotating frame of reference. Thus, trajectories of particles ar no longer straight lines. In the present study, trajectories of particles are calculated by integrating a system of differential equations including the Coriolis and centrifugal forces. The null-collisions. The present numerical results molecular model is used for calculation of molecular collisions. The present numerical results significantly disagree with the previously known ones. This indicates that an actual pumping passage is very limited to a narrow region due to the significant backstreaming of molecules from the outlet.

• Journal of Ocean Engineering and Technology
• /
• v.24 no.3
• /
• pp.31-37
• /
• 2010

During typhoon periods, coastal regions are often directly flooded by typhoon-surges. There are also many cases where coastal regions are inundated by river inundations or dam breaks. However, most studies on coastal flooding by typhoons have been restricted to cases involving the sea. Flooding by river inundation has been excluded in those studies. Usually ocean numerical models are not applied to river flow because the governing equations for ocean flow and river flow are not the same. For a coastal flooding simulation with river inundation, POM, the three-dimensional numerical ocean model, was applied to the popular river flow problems, dam-break problem, and flows over a spillway. The simulated results showed good agreement with other numerical simulations and measured data, suggesting the possibility of using POM in coastal flooding simulations involving direct coastal surges and river inundations.

• Journal of the Korean Society of Safety
• /
• v.24 no.1
• /
• pp.14-20
• /
• 2009

Direct numerical simulations(DNS) were performed for the prediction of transient buoyant jet diffusion flames where the Froude numbers(Fr) are 5 and 160, respectively. The thermodynamic and transport properties were evaluated using CHEMKIN package to enhance the prediction performance of the DNS code. The simulated buoyant jet diffusion flame of Fr=5 and 160 showed the transient, dynamic motion well. It was identified that the buoyant jet flames were flickered periodically, and the simulated flickering frequency of the jet diffusion flame of Fr=5 was 12.5Hz, which was in good agreement with the experimental results. The flamelet structures of the buoyant jet diffusion flames could be well understood by comparing the scalar dissipation rates(SDR) and the heat release rates(HRR) of the flames. It was found that the SDR was strongly coupled with the HRR in the buoyant jet diffusion flames.

• Journal of Korean Society of Environmental Engineers
• /
• v.31 no.10
• /
• pp.933-940
• /
• 2009

To evaluate characteristics in spray flame, laminar counterflow is investigated on the effects of equivalence ratio and fuel by a two-dimensional DNS (direct numerical simulation). For the gaseous phase, Eulerian mass, momentum, energy, and species conservation equations are solved. For the disperse phase, all individual droplets are calculated by the Lagrangian method without the parcel model. n-Decane ($C_{10}H_{22}$) and n-heptane ($C_7H_{16}$) is used as a liquid spray fuel, and a one-step global reaction is employed for the combustion reaction model. As equivalence ratio increases, the fuel ignites early and the high temperature region spreads wider. The peak value of temperature, however, tends to once increase and then decreases with increasing equivalence ratio. The decrease in the peak value of temperature for the higher equivalence ratio condition is caused by the cooling effect associated with droplet group combustion. Since the evaporation of n-heptane is early, the high temperature region spreads wider than ndecane, but the peak values of temperature for both n-heptane and n-decane is almost same.

• Journal of the Korean Society of Combustion
• /
• v.8 no.3
• /
• pp.15-23
• /
• 2003

Partial quenching structure of turbulent diffusion flames in a turbulent mixing layer is investigated by the method of flame hole dynamics in order to develop a prediction model for turbulent flame lift off. The essence of flame hole dynamics is derivation of the random walk mapping, from the flame-edge theory, which governs expansion or contraction of flame holes initially created by local quenching events. The numerical simulation for flame hole dynamics is carried out in two stages. First, a direct numerical simulation is performed for constant-density fuel-air channel mixing layer to obtain the turbulent flow and mixing fields, from which a time series of two dimensional scalar dissipation rate array is extracted at a fixed virtual flame surface horizontally extending from the end of split plate to the downstream. Then, the Lagrangian simulation of the flame hole random walk mapping projected to the scalar dissipation rate array yields temporally evolving turbulent extinction process and its statistics on partial quenching characteristics. The statistical results exhibit that the chance of partial quenching is strongly influenced by the crossover scalar dissipation rate while almost unaffected by the iteration number of the mapping that can be regarded as a flame-edge speed.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.32 no.8
• /
• pp.636-644
• /
• 2008

The objective of this paper is to investigate characteristics of an autothermal reformer at various operating conditions. Numerical method has been used, and simulation model has been developed for the analysis. Pseudo-homogeneous model is incorporated because the reactor is filled with catalysts of a packed-bed type. Dominant chemical reactions are Full Combustion reaction, Steam Reforming(SR) reaction, Water-Gas Shift(WGS) reaction, and Direct Steam Reforming(DSR) reaction. Simulation results are compared with experimental results for code validation. Operating parameters of the autothermal reformer are inlet temperature, Oxygen to Carbon Ratio(OCR), Steam to Carbon Ratio(SCR), and Gas Hourly Space Velocity(GHSV). Temperature at the reactor center, fuel conversion, species at the reformer outlet, and reforming efficiency are shown as simulation results. SR reaction rate is improved by increased inlet temperature. Reforming efficiency and fuel conversion reached the maximum at 0.7 of OCR. SR reaction and WGS reaction are activated as SCR increases. When GHSV is increased, reforming efficiency increases but pressure drop from the increased GHSV may decrease the system efficiency.

• Computers and Concrete
• /
• v.29 no.6
• /
• pp.407-418
• /
• 2022

This paper numerically investigates the effect of changes in the mechanical properties (displacement, strain, and stress) of the ultra-high-performance concrete (UHPC) without rebar and the reinforced concrete (RC) using steel re-bars. This reinforced concrete is mostly used in the concrete bridge decks. A mixture of sand, gravel, cement, water, steel fiber, superplasticizer, and micro silica was used to fabricate UHPC specimens. The extended finite element method as used in the ABAQUS software is applied for considering the mechanical properties of UHPC, RC, and ordinary concrete specimens. To calibrate the ABAQUS, some experimental tests have been carried out in the laboratory to measure the direct tensile strength of UHPC by the compressive-to-tensile load converting (CTLC) device. This device contains a concrete specimen and is mounted on a universal tensile testing apparatus. In the experiments, three types of mixed concrete were used for UHPC specimens. The tensile strength of these specimens ranges from 9.24 to 11.4 MPa, which is relatively high compared with ordinary concrete specimens, which have a tensile strength ranging from 2 to 5 MPa. In the experimental tests, the UHPC specimen of size 150×60×190 mm with a central hole of 75 mm (in diameter)×60 mm (in thickness) was specially made in the laboratory, and its direct tensile strength was measured by the CTLC device. However, the numerical simulation results for the tensile strength and failure mechanism of the UHPC were very close to those measured experimentally. From comparing the numerical and experimental results obtained in this study, it has been concluded that UHPC can be effectively used for bridge decks.

• Journal of the Korean Society of Safety
• /
• v.10 no.3
• /
• pp.120-129
• /
• 1995

Since the rate and completeness of combustion in direct injection engines were controlled by the characteristics of gas flow fields and sprays, an understanding of those was essential to the design of the direct injection engines. In this study the numerical simulations of swirl effects on the characteristics of gas flow fields and sprays were performed using the spray model that could predict the interactions between gas fields and spray droplets. The governing equations were discretized by the finite volume method and the modified k- e model which included the compressibility effects due to the compression/expansion of piston was used. The results of numerical calculation of the spray characteristics in the quiescent environment were compared with the experimental data. There were good agreements between the results of calculation and the experimental data, except in the early stages of spray. In the motoring condition, the results showed that a substantial air entrainment into the spray volume was emerged and hence the squish motion was relatively unimportant during fuel injection periods. As the swirl ratio increased, the evaporation rate was increased due to the wide dispersion of the spray droplets and the strong interaction between spray droplets and gas fields.

• Structural Engineering and Mechanics
• /
• v.64 no.5
• /
• pp.621-640
• /
• 2017

This paper reports the numerical investigation conducted to study the influence of Local-Distortional (L-D) interaction mode buckling on post buckling strength erosion in fixed ended lipped channel cold formed steel columns. This investigation comprises of 81 column sections with various geometries and yield stresses that are carefully chosen to cover wide range of strength related parametric ratios like (i) distortional to local critical buckling stress ratio ($0.91{\leq}F_{CRD}/F_{CRL}{\leq}4.05$) (ii) non dimensional local slenderness ratio ($0.88{\leq}{\lambda}_L{\leq}3.54$) (iii) non-dimensional distortional slenderness ratio ($0.68{\leq}{\lambda}_D{\leq}3.23$) and (iv) yield to non-critical buckling stress ratio (0.45 to 10.4). The numerical investigation is carried out by conducting linear and non-linear shell finite element analysis (SFEA) using ABAQUS software. The non-linear SFEA includes both geometry and material non-linearity. The numerical results obtained are deeply analysed to understand the post buckling mechanics, failure modes and ultimate strength that are influenced by L-D interaction with respect to strength related parametric ratios. The ultimate strength data obtained from numerical analysis are compared with (i) the experimental tests data concerning L-D interaction mode buckling reported by other researchers (ii) column strength predicted by Direct Strength Method (DSM) column strength curves for local and distortional buckling specified in AISI S-100 (iii) strength predicted by available DSM based approaches that includes L-D interaction mode failure. The role of flange width to web depth ratio on post buckling strength erosion is reported. Then the paper concludes with merits and limitations of codified DSM and available DSM based approaches on accurate failure strength prediction.