• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.3
• /
• pp.439-447
• /
• 2002

The microchannel flow in miniature TCDs (thermal conductivity detectors) is investigated numerically. The solutions based on the boundary layer approximation are not very accurate in the region of the duct inlet for low Reynolds numbers. In this study, two-dimensional Navier-Stokes equations are considered to analyze the gas flow in a miniature TCD. Effects of channel size, inlet and boundary conditions on the heat transfer rate are examined. When the gas stream is not preheated, the distances for a miniature TCD to reach the conduction-dominant region for duct flow are found to be approximately two and three times the thermal entry length for duct flow with constant properties, respectively, leer constant wall temperature and constant wall heat flux boundary conditions. If the gas temperature at the channel inlet is close to the mean gas temperature in the conduction-dominant region, the entrance region is much shorter compared to other cases considered in this study.

• Applied Chemistry for Engineering
• /
• v.20 no.6
• /
• pp.646-652
• /
• 2009

Present study deals with a three phase flow problem of determining drag acting on spheres wetted by liquid flow by gas flow through the spheres in simple cubic (SC), body-center cubic (BCC) and face-centered cubic (FCC) array, respectively, when the inertia of gas is negligibly small. The liquid flow driven by gravity on the spheres is assumed to be unaffected by the countercurrent gas flow. A perturbation method coupled with a multipole expansion method is used to calculate the hydrodynamic interactions between spheres and hence determine the effect of liquid film and flow on the gas flow for each periodic array of spheres. An approximate method for evaluating the effect of the liquid film is also presented for simple estimations. It is found that the approximation results are in a reasonable agreement with the numerical calculations.

• Journal of the Society of Naval Architects of Korea
• /
• v.42 no.3
• /
• pp.205-211
• /
• 2005

In this study, a new numerical procedure for the generation of a nonlinear tailored group of waves is presented. The procedure is based on the transient wave group technique. In order to integrate the nonlinearity during the wave propagation in the computational method, the Navier-Stokes equations are applied as governing equations. The governing equations are discretized by finite volume approximation. The deformation of the free water surface in each time step is pursued with a moving grid. A two-dimensional, numerical wave tank for the simulation of the wave propagation is developed and tested in detail. The numeric results are compared first with analytical wave theories and with measurements, in order to examine the correctness of the numerical wave tank. Wave surface elevation and associated fields of velocity and pressure are numerically computed and compared with measurements. Very good agreements show up.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.32 no.10
• /
• pp.754-760
• /
• 2008

Computation of moving interface by the level set method typically requires the reinitialization of level set function. An inaccurate estimation of level set function $\phi$ results in incorrect free-surface capturing and thus errors such as mass gain/loss. Therefore, an accurate and robust reinitialization process is essential to the simulation of free-surface flows. In the present paper, we pursue further development of the reinitialization process, which evaluates level set function directly using a normal vector on the interface without solving there-distancing equation of hyperbolic type. The Taylor-Galerkin approximation and P1P1 splitting/SUPG (Streamline Upwind Petrov-Galerkin) FEM are adopted to discretize advection equation of the level set function and the incompressible Navier-Stokes equation, respectively. Advection equation and re-initialization process of free surface capturing are validated with benchmark problems, i.e., a broken dam flow and timereversed single vortex flow. The simulation results are in good agreement with the existing results.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.4 no.4
• /
• pp.219-224
• /
• 1992

A numerical method using the truncated Fourier series is presented to predict the wave potential and water surface profile for two dimensional nonlinear standing waves. The unknown coefficients of the series are to be determined through the Newton solution of nonlinear simultaneous equations given by the governing equation and boundary conditions of the problem. In order to prove the effectiveness of the present method. an existing Stokes-like perturbation method is considered together, and a hydraulic experiment for measuring water surface profile and wave pressure is performed as well. The results are such that the present method can generally give exact solutions even for relatively big wave stiffness regardless of the water depth condition. It also demonstrates its validity by showing double humps in the crest of temporal wave pressure profile which normally appear in strongly nonlinear standing waves.

• Journal of The Korean Astronomical Society
• /
• v.32 no.1
• /
• pp.65-73
• /
• 1999

In this study we present a new improved nonlinear calibration method for vector magnetograms made by the Solar Flare Telescope of BOAO. To identify Fe I 6302.5 line, we have scanned monochromatic images of the line integrated over filter passband, changing the location of the central transmission wavelength of a Lyot filter. Then we obtained a filter-convolved line profile, which is in good agreement with spectral atlas data provided by the Sacramento Peak Solar Observatory. The line profile has been used to derive calibration coefficients of longitudinal and transverse fields, employing the conventional line slope method under the weak field approximation. Our improved nonlinear calibration method has also been used to calculate theoretical Stokes polarization signals with various angles of inclination of magnetic fields. For its numerical test, we have compared input magnetic fields with the calibrated ones, which have been derived from the new improved non-linear method and the conventional method respectively. The numerical test shows that the calibrated fields obtained from the improved method are consistent with the input fields, but not with those from the conventional method. Finally, we applied our new improved method to a dipole model which characterizes a typical field configuration of a single, round sunspot. It is noted that the conventional method remarkably underestimates the transverse field component near the inner penumbra.

• International Journal of Aeronautical and Space Sciences
• /
• v.3 no.1
• /
• pp.95-104
• /
• 2002

Numerical solution of the full Navier-Stokes equation as well as the energy equation has been obtained for the unsteady natural convection in a rectangular enclosure. One side wall was maintained at very high temperature simulating fires. Especially the effect of surface radiation was taken into account. While the enclosed air was assumed to be transparent, the internal walls directly interacted one another through the surface radiation. Due to a significant temperature difference in the flow field, the equation of state was used instead of the Boussinesq approximation. It was found that the rapid heating of the adiabatic ceiling and floor by the incoming radiation from the hot wall made the evolution at thermo-fluid field highly unstable in the initial period. Therefore, the secondary cells brought about at the floor region greatly affected the heat transfer mechanism inside the enclosure. The heat transfer rate was augmented by the radiation, resulting in requiring less time for the flow to reach the steady state. At the steady state neglecting radiation two internal hydraulic jumps were clearly observed in upper/left as well as in lower/right comer. However, the hydraulic jump in the lower/right comer could not be observed for the case including radiation due to its high momentum flow over the bottom wall. Radiation resulted in a faster establishment of the steady state phenomena.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.2
• /
• pp.355-365
• /
• 1994

An analytical model of return flow is presented outside the surf zone. The governing equation is derived from the Navier-Stokes equation and the continuity. Each term of the governing equation is evaluated by the ordering analysis. Then the infinitesimal terms, i.e. the turbulent normal stress, the squared vertical velocity of water particle and the streaming velocity, are neglected. The driving forces of return flow are calculated using the linear wave theory for the shallow water approximation. Especially, the space derivative of local wave heights is described considering a shoaling coefficient. The vertical distribution of eddy viscosity is discussed to the customary types which are the constant, the linear function and the exponential function. Each coefficient of the eddy viscosities which sensitively affect the precision of solutions is uniquely decided from the additional boundary condition which the velocity becomes zero at the wave trough level. Also the boundary conditions at the bottom and the continuity relation are used in the integration of the governing equation. The theoretical solutions of present model are compared with the various experimental results. The solutions show a good agreement with the experimental results in the case of constant or exponential function type eddy viscosity.

• Journal of Hydrogen and New Energy
• /
• v.22 no.3
• /
• pp.386-401
• /
• 2011

Mathematical models for various steps in coal gasification reactions were developed and applied to investigate the effects of operation parameters on dynamic behavior of gasification process. Chemical reactions considered in these models were pyrolysis, volatile combustion, water shift reaction, steam-methane reformation, and char gasification. Kinetics of heterogeneous reactions between char and gaseous agents was based on Random pore model. Momentum balance and Stokes' law were used to estimate the residence time of solid particles (char) in an up-flow reactor. The effects of operation parameters on syngas composition, reaction temperature, carbon conversion were verified. Parameters considered here for this purpose were $O_2$-to-coal mass ratio, pressure of reactor, composition of coal, diameter of char particle. On the basis of these parametric studies some quantitative parameter-response relationships were established from both dynamic and steady-state point of view. Without depending on steady state approximation, the present model can describe both transient and long-time limit behavior of the gasification system and accordingly serve as a proto-type dynamic simulator of coal gasification process. Incorporation of heat transfer through heterogenous boundaries, slag formation and steam generation is under progress and additional refinement of mathematical models to reflect the actual design of commercial gasifiers will be made in the near futureK.

• Nuclear Engineering and Technology
• /
• v.56 no.4
• /
• pp.1310-1319
• /
• 2024

The reactor cavity cooling system (RCCS) is a passive reactor safety system commonly present in the designs of High-Temperature Gas-cooled Reactors (HTGR) that removes heat from the reactor pressure vessel by means of natural convection and radiation. It is one of the factors responsible for ensuring that the reactor does not melt down under any plausible accident scenario. For the simulation of accident scenarios, which are transient phenomena unfolding over a span of up to several days, intermediate fidelity methods and system codes must be employed to limit the models' execution time. These models can quantify radiation heat transfer well, but heat transfer caused by natural convection must be quantified with the use of correlations for the heat transfer coefficient. It is difficult to obtain reliable correlations for HTGR RCCS heat transfer coefficients experimentally due to such a system's size. They could, however, be obtained from high-fidelity steady-state simulations of RCCSs. The Rayleigh number in RCCSs is too high for using a Direct Numerical Simulation (DNS) technique; thus, a Reynolds-Averaged Navier-Stokes (RANS) approach must be employed. There are many RANS models, each performing best under different geometry and fluid flow conditions. To find the most suitable one for simulating an RCCS, the RANS models need to be validated. This work benchmarks various RANS models against three experiments performed on the HTTR RCCS Mockup by the Japanese Atomic Energy Agency (JAEA) in 1993. This facility is a 1/6 scale model of a vessel cooling system (VCS) for the High Temperature Engineering Test Reactor (HTTR), which is operated by JAEA. Multiple RANS models were evaluated on a simplified 2d-axisymmetric geometry. They were found to reproduce the experimental temperature profiles with errors of up to 22% for the lowest temperature benchmark and 15% for the higher temperature benchmarks. The results highlight that the pragmatic turbulence models need to be validated for high Rayleigh natural convection-driven flows and improved accordingly, more publicly available experimental data of RCCS resembling experiments is needed and indicate that a 2d-axisymmetric geometry approximation is likely insufficient to capture all the relevant phenomena in RCCS simulations.