• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.8
• /
• pp.1040-1047
• /
• 1999

The characteristics of flame propagation in inert particle-laden $H_2$/Air premixed gas are numerically investigated on this study. The 2nd order TVD scheme is applied to numerical analysis of governing equations and multi-step chemical reaction model and detailed transport properties are sued to solve chemical reaction terms. Radiation heat transfer is computed by applying the finite volume method to a radiative transfer equation. The burning velocities against the mole fractions of hydrogen agree well with results performed by different workers. The inert particles play significant roles in the flame propagation on account of momentum and heat transfer between gas and particles. Gas temperature, pressure and flame propagation speed are decreased as the loading ratio of particle is increased. Also the products behind flame zone contain lots of water vapor whose absorption coefficient is much larger than that of unburned gas. Thus, the radiation effect of gas and particles must be considered simultaneously for the flame propagation in a mixture of $H_2$/Air and inert particles. As a result, it is founded that because the water vapor emits much radiation and this emitted radiation is released at boundaries as radiant heat loss as well as reabsorbed by gas and particles, flame propagation speed and flame structure are altered with radiation effect.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.39 no.3
• /
• pp.283-289
• /
• 2015

In this study, the lightweight design of a brake bracket for a composite bogie was studied by considering two brake bracket models with thicknesses of 12t and 9t, respectively. For achieving this goal, finite element analysis and topology optimization were conducted. Firstly, the largest cross-sectional areas of the vertical and horizontal plates of the brake bracket were selected as the design variables. As the constraint, the Z-axis displacement of the brake bracket was increased by 2.5 units from the initial displacement value. The minimum volume fraction of the design regions was chosen as the objective function. The full model comprised a composite bogie frame and brackets attached together. However, to reduce the analysis time, 1D beam elements were used instead of the composite bogie frame by ensuring its equivalence with the full model. The result revealed that the weights of the 12t and 9t models of the brake bracket were reduced to 60 kg and 31 kg, respectively.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.41 no.6
• /
• pp.429-434
• /
• 2017

A molecular-continuum hybrid method was developed to simulate microscale and nanoscale fluids where continuum fluidics cannot be used to predict Couette flow. Molecular dynamics simulation is used near the solid surface where the flow cannot be predicted by continuum fluidics, and Navier-Stokes equations are used in the other regions. Numerical simulation of Couette flow was performed using the hybrid method to investigate the effect of solid-liquid interaction and surface roughness in a nanochannel. It was found that the solid-liquid interaction and surface roughness influence the boundary condition. When the surface energy is low, slippage occurs near the solid surface, and the magnitude of slippage decreases with increase in surface energy. When the surface energy is high, a locking boundary condition is formed. The roughness disturbs slippage near the solid surface and promotes the locking boundary condition.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.10
• /
• pp.1940-1954
• /
• 1992

Fine grid computations were attempted to analyze the turbulent flows in the near wall low Reynolds number region and the numerical analyses were incorporated by a finite-volume discretization with full find grid system and low Reynolds number k-.epsilon. model was employed in this region. For the improvement of low Reynolds number k-.epsilon. model, modification coefficient of eddy viscosity $f_{\mu}$ was derived as a function of turbulent Reynolds number $R_{+}$ and nondimensional length $y^{+}$ from the concept of two length scales of dissipation rate of turbulent kinetic energy. The modification coefficient $f_{\epsilon}$ in .epsilon. transport equation was also derived theoretically. In the turbulent kinetic energy equation, pressure diffusion term was added in order to consider low Reynolds number region effect. The main characteristics of this low Reynolds number k-.epsilon. model were founded as : (1) In high Reynolds number region, the present model has limiting behavior which approaches to the high Reynolds number model. (2) Present low Reynolds number k-.epsilon. model dose not need additional empirical constants for the transport equations of turbulent kinetic energy and dissipation of turbulent kinetic energy in order to consider wall effect. Present low Reynolds number turbulence model was tested in the pipe flow and obtained improved results in velocity profiles and Reynolds stress distributions compared with those from other k-.epsilon. models.s.s.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.2
• /
• pp.235-251
• /
• 1997

A numerical simulation on the combined natural convection and radiation is carried out in a partially open rectangular enclosure with a heater by using the finite volume and the S-8 discrete ordinate methods. The fluid inside the enclosure is considered as an absorbing, emitting and anisotropic scattering media. The heater causes a natural circulation of the fluid (10$^{5}$ $^{9}$ ) which results in significant in-flow of the ambient cold fluid through the partially open wall. Comparing the results of pure convection with those of the combined convection- radiation, the combined heat transfer results with small Planck numbers (P$_{l}$ <1.0) show much stronger circulation than those of the pure convection, and the fluid circulation is more evident for larger Rayleigh numbers. When one of three radiative properties - the medium absorption coefficient, the wall reflectivity, and the scattering albedo - increases, the fluid circulation and the heat transfer in the enclosure are reduced. The location of the heater and the open ratio of the right wall are also shown to affect the fluid circulation and heat transfer significantly. However, the anisotropy of the scattering phase function is shown to be unimportant for the fluid circulation and heat transfer within the enclosure considered in this study.

• Journal of Korea Water Resources Association
• /
• v.44 no.2
• /
• pp.145-156
• /
• 2011

Two dimensional numerical model of high-order accuracy is developed to analyze complex flow including transition flow, discontinuous flow, and wave propagation to dry bed emerging at natural river flow. The bed slope term of two dimensional shallow water equation consisting of integral conservation law is treated efficiently by applying quasi-steady wave propagation scheme. In order to apply Finite Volume Method using Fractional Step Method, MUSCL scheme is applied based on HLL Riemann solver, which is second-order accurate in time and space. The TVD method is applied to prevent numerical oscillations in the second-order accurate scheme. The developed model is verified by comparing observed data of two dimenstional levee breach experiment and dam breach experiment containing structure at lower section of channel. Also effect of the source term is verified by applying to dam breach experiment considering the adverse slope channel.

• Journal of Korea Water Resources Association
• /
• v.42 no.4
• /
• pp.271-279
• /
• 2009

The objective of this study is to develop the scheme to apply one-dimensional finite volume method (FVM) to natural river with complex geometry. In the previous study, FVM using the Riemann approximate solver was performed successfully in the various cases of dam-break, flood propagation, etc. with simple and rectangular cross-sections. We introduced the transform the natural into equivalent rectangular cross-sections. As a result of this way, the momentum equation was modified. The accuracy and applicability of newly developed scheme are demonstrated by means of a test example with exact solution, which uses triangular cross-sections. Secondly, this model is applied to natural river with irregular cross-sections and non-uniform lengths between cross-sections. The results shows that the aspect of flood propagation, location and height of hydraulic jump, and numerical solutions of maximum water level are in good agreement with the measured data. Using the developed scheme in this study, existing numerical schemes conducted in simple cross-sections can be directly applied to natural river without complicated numerical treatment.

• Membrane Journal
• /
• v.12 no.4
• /
• pp.226-239
• /
• 2002

A process simulation model of pervaporation process has been developed as a design tool to analyse and optimize the dehyhration of organic solvents through a commercial scale of pervaporative plate-and-frame modules that contain a stack of membrane sheets. In the simulation model, the mass balance, the heat balance and the concentration balance are integrated in a finite elements-in-succession method to simulate the overall process. In the integration method, a feed channel between membrane sheets in the modules was taken as differential unit element volume to simplify calculation procedure and shorten computing time. Some of permeation parameters used in the simulation model, were quantified directly from the dehydration experiment of ethanol through $AzeoSep^{TM}$-2002 membrane which is a commercial pervaporation membrane. The simulation model was verified by comparing the simulated values with experimental data. Using the model, continuous and batch pervaporation processes were simulated, respectively, to acquire basic data for analysing and optimizing in the dehydration of ethanol through the membrane. Based on the simulation results, a comparison between the continuous and the batch pervaporation processes would be discussed.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.7 no.4 s.19
• /
• pp.70-78
• /
• 2004

A study has been made to investigate the high temperature deformation behavior of Ti-6Al-4V alloyand to predict the final microstructure under given forming conditions. Equiaxed and $Widmanst\ddot{a}tten$ microstructures of Ti-6Al-4V alloys were prepared as initial microstructures. By performing the compression tests at high temperatures$(700\~1100^{\circ}C)$ and at a wide range of strain rates$(10^{-4}\~10^2/s)$, various parameters such as strain rate sensitivity(m) and activation energy(Q) were calculated and used to establish constitutive equations. When the specimens were deformed up to strain 0.6, equiaxed microstructure did not show any significant changes in microstructure, while $Widmanst\ddot{a}tten$ microstructure revealed considerable flow softening, which was attributed to the globularization of a platelet at the temperature range of $800\~970^{\circ}C$ and at the strain rate range of $10^{-4}\~10^{-2}/s$. To predict the final microstructure after forming, finite element analysis was performed considering the microstructural evolution during the deformation. The grain size and the volume fraction of second phase of deformed body were predicted and compared with the experimental results.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.146-153
• /
• 2004

The Numerical study has been carried out to investigate the effects of chemical reaction and thermal radiation on the rocket plume flow-field at various altitudes. The theoretical formulation is based on the Navier-Stokes equations for compressible flows along with the infinitely fast chemistry and thermal radiation. The governing equations were solved by a finite volume fully-implicit TVD(Total Variation Diminishing) code which uses Roe's approximate Riemann solver and MUSCL(Monotone Upstream-centered Schemes for Conservation Laws) scheme. LU-SGS (Lower Upper Symmetric Gauss Seidel) method is used for the implicit solution strategy. An equilibrium chemistry module for hydrocarbon mixture with detailed thermo-chemical properties and a thermal radiation module for optically thin media were incorporated with the fluid dynamics code. In this study, kerosene-fueled rocket was assumed operating at O/F ratio of 2.34 with a nozzle expansion ratio of 6.14. Flight conditions considered were Mach number zero at ground level, Mach number 1.16 at altitude 5.06km and Mach number 2.9 at altitude 17.34km. Numerical results gave the understandings on the detailed plume structures at different altitude conditions. The diffusive effect of the thermal radiation on temperature field and the effect of chemical recombination during the expansion process could be also understood. By comparing the results from frozen flow and infinitely fast chemistry assumptions, the excess temperature of the exhaust gas resulting from the chemical recombination seems to be significant and cannot be neglected in the view point of performance, thermal protection and flow physics.