• Journal of The Korean Astronomical Society
• /
• v.39 no.4
• /
• pp.95-105
• /
• 2006

We have calculated the cosmic ray(CR) acceleration at young remnants from Type Ia supernovae expanding into a uniform interstellar medium(ISM). Adopting quasi-parallel magnetic fields, gasdynamic equations and the diffusion convection equation for the particle distribution function are solved in a comoving spherical grid which expands with the shock. Bohm-type diffusion due to self-excited $Alfv\acute{e}n$ waves, drift and dissipation of these waves in the precursor and thermal leakage injection were included. With magnetic fields amplified by the CR streaming instability, the particle energy can reach up to $10^{16}Z$ eV at young supernova remnants(SNRs) of several thousand years old. The fraction of the explosion energy transferred to the CR component asymptotes to 40-50 % by that time. For a typical SNR in a warm ISM, the accelerated CR energy spectrum should exhibit a concave curvature with the power-law slope flattening from 2 to 1.6 at $E{\gtrsim}0.1$ TeV.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.6
• /
• pp.772-778
• /
• 2002

A numerical analysis of turbulent gas-particle two-phase flow is performed in conjunction with the experiments of Fackrell ＆ Robins and Raupach & Legg that considered ground-level source and/or elevated source flat plate flow. K-$\omega$ turbulence model is used in order to analyze fully turbulent flow field and the concentration equation with settling velocity is adopted for the concentration field. The model of Einstein and Chien is applied that couples the velocity field and the concentration field. Turbulent eddy viscosity is re-evaluated in this model. The present numerical results have good agreement between the simulation and the experimental data for the mean flow velocities and particle concentrations. While the previous study shows about 27% error in the vicinity of the source of particle concentration, the .present study allows about 14% error. A new turbulent gas-particle flow model developed by this study is able to cut down error by 13% at a near source.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.1
• /
• pp.110-115
• /
• 2002

The tunnels built in recent years are equipped with traffic counters and pollution sensors (mostly, CO and Vl sensors). Utilizing these built-in sensors, it is possible to develop an algorithm to estimate the amount of pollutants exhausted from the each class of cars passing through the tunnel. These estimated data can be effectively utilized not only for ventilation control but also for designing ventilation facilities. The diffusion of pollutants in a tunnel can be described with one-dimensional diffusion-convection equation. This equation is approximated with interpolation functions and weighted residual method converting to adequate form for standard state estimate algorithms. With this converted equations, a least square optimization based algorithm is developed, whose outputs are the estimated amounts of pollutants emitted from each class of cars. In order to verify the feasibility of the developed algorithms, simulations are performed with the real data acquisitioned from the Tunnae tunnel located in Young-Dong highway in Korea.

• Proceedings of the Korea Institute of Fire Science and Engineering Conference
• /
• 1996.11a
• /
• pp.5-9
• /
• 1996

In foreign country such as U.S.A and Japan, considerable research has been done regarding the spread of smoke in room of fire involvement by using computer. But, in our country it has not been. So, this paper presents a detailed qualitative description of phenomena which occures during typical fire scenarios through numerical analysis. This research, in the view of field model, is focused on finding out the smoke movement and temperature distribution. And it is planned to analyze governing equation including smoke diffusion equation by numerical analysis with finite volume method and non-staggered grid system. The SIMPLE method for pressure-velocity couple and power-law scheme for convection terms are used. It shows that a plume is formed, hot plume is formed, hot plume gases impinge on the ceiling and they spread across it. then, it eventually reaches the bounding walls of the enclosure. It takes 60s for smoke to fill the enclosure.

• International Journal of Ocean System Engineering
• /
• v.1 no.1
• /
• pp.37-45
• /
• 2011

The nonlinear free-surface motions interacting with a floating body were investigated using the Moving Particle Semi-implicit (MPS) method proposed by Koshizuka and Oka [6] for incompressible flow. In the numerical method, more realistic Lagrangian moving particles were used for solving the flow field instead of the Eulerian approach with a grid system. Therefore, the convection terms and time derivatives in the Navier-Stokes equation can be calculated more directly, without any numerical diffusion, instabilities, or topological failure. The MPS method was applied to a numerical simulation of predicting the efficiency of floating-type breakwater interacting with waves.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.12
• /
• pp.1715-1721
• /
• 2002

Numerical analysis was conducted to characterize the gas flow field and particle deposition on a horizontal freestanding semiconductor wafer under the laminar flow field at vacuum environment. In order to calculate the properties of gas, the gas was assumed to obey the ideal gas law. The particle transport mechanisms considered were convection, Brownian diffusion and gravitational settling. The averaged particle deposition velocities and their radial distributions fnr the upper surface of the wafer were calculated from the particle concentration equation in an Eulerian frame of reference for system pressures of 1 mbar~1 atm and particle sizes of 2nm~10$^4$ nm(10 ${\mu}{\textrm}{m}$). It was observed that as the system pressure decreases, the boundary layer of gas flow becomes thicker and the deposition velocities are increased over the whole range of particle size. One thing to be noted here is that the deposition velocities are increased in the diffusion dominant particle size range with decreasing system pressure, whereas the thickness of the boundary layer is larger. This contradiction is attributed to the increase of particle mechanical mobility and the consequent increase of Brownian diffusion with decreasing the system pressure. The present numerical results showed good agreement with the results of the approximate model and the available experimental data.

• 한국연소학회:학술대회논문집
• /
• 1998.10a
• /
• pp.123-129
• /
• 1998

A numerical study has been conducted to investigate the effect of shock waves on the mixing and the recirculation zone of a hydrogen jet diffusion flame in a supersonic combustor. The general trends are compared with the experimental results obtained from the supersonic combustor at the University of Michigan. For the numerical simulation of supersonic diffusion flames, multi-species Navier-Stokes equations and detailed chemistry reaction equations of $H_2$-Air are considered. The $K-{\omega}/k-{\varepsilon}$ blended two equation turbulent model is used. Roe's FDS method and MUSCL method are used for convection fluxes in governing equations. Numerical results show that when slender wedges are mounted at the combustor wall the mixing and the combustion are enhanced and the size of recirculation zone is increased . The flame shape of supersonic flames is different in the flame-tip; it is not closed but open. The flame shape is shown to be greatly affected by shock waves.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.4
• /
• pp.328-333
• /
• 2014

A hybrid particle-mesh method based on the vorticity-velocity formulation for solving the incompressible Navier-Stokes equations is a combination of the Vortex-In-Cell(VIC) method for convection and the penalization method for diffusion. The key feature of the numerical methods is to determine velocity and vorticity fields around a solid body on a temporary grid, and then the time evolution of the flow is computed by tracing the convection of each vortex element using the Lagrangian approach. Assuming that the vorticity and velocity fields are to be computed in time domain analysis, pressure fields are estimated through a complete set of solutions at present time step. It is possible to obtain vorticity and velocity fields prior to any pressure calculation since the pressure term is eliminated in the vorticity-velocity formulation. Therefore, pressure field is explicitly treated by solving a suitable Poisson equation. In this paper, we propose a simple way to numerically implement the vorticity-velocity-pressure formulation including a penalty term. For validation of the proposed numerical scheme, we illustrate the early development of viscous flows around an impulsive started circular cylinder for Reynolds number of 9500.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.10
• /
• pp.1284-1293
• /
• 1997

The electrostatic effect on particle deposition onto a heated, Horizontal free-standing wafer surface was investigated numerically. The deposition mechanisms considered were convection, Brownian and turbulent diffusion, sedimentation, thermophoresis and electrostatic force. The electric charge on particle needed to calculate the electrostatic migration velocity induced by the local electric field was assumed to be the Boltzmann equilibrium charge. The electrostatic forces acted upon the particle included the Coulombic, image, dielectrophoretic and dipole-dipole forces based on the assumption that the particle and wafer surface are conducting. The electric potential distribution needed to calculate the local electric field around the wafer was calculated from the Laplace equation. The averaged and local deposition velocities were obtained for a temperature difference of 0-10 K and an applied voltage of 0-1000 v.The numerical results were then compared with those of the present suggested approximate model and the available experimental data. The comparison showed relatively good agreement between them.

• Journal of applied mathematics & informatics
• /
• v.22 no.1_2
• /
• pp.49-65
• /
• 2006

A robust numerical method for a singularly perturbed second-order ordinary differential equation having two parameters with a discontinuous source term is presented in this article. Theoretical bounds are derived for the derivatives of the solution and its smooth and singular components. An appropriate piecewise uniform mesh is constructed, and classical upwind finite difference schemes are used on this mesh to obtain the discrete system of equations. Parameter-uniform error bounds for the numerical approximations are established. Numerical results are provided to illustrate the convergence of the numerical approximations.