• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.10
• /
• pp.1089-1095
• /
• 2011

Numerical simulation is performed for nucleate boiling on a micro-finned surface, which has been widely used to enhance heat transfer, by solving the equations governing the conservation of mass, momentum, and energy in the liquid and vapor phases. The bubble motion is determined by a sharp-interface level-set method, which is modified to include the effect of phase change and to treat the no-slip and contact-angle conditions, as well as the evaporative heat flux from the liquid microlayer on immersed solid surfaces such as micro fins and cavities. The numerical results for bubble formation, growth, and departure on a microstructured surface including fins and cavities show that the bubble behavior during nucleate boiling is significantly influenced by the fin-cavity arrangement and the fin-fin spacing.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.23 no.4
• /
• pp.276-284
• /
• 2011

The rip current occurred at Haeundae beach was numerically investigated under directional random wave environment. The numerical simulation was performed using a fully nonlinear Boussinesq equation model, FUNWAVE which is capable of simulating nearshore circulation since it includes the effect of wave-induced momentum flux and horizontal turbulent mixing. The results of numerical simulation show the time-dependent evolution of the wave-induced nearshore circulation system (including rip current) that are caused by nonlinear transformation of directional irregular waves due to unique topography of Haeundae. From the results, it was found that rip current is well generated and developed where relatively lower wave height and relatively deeper water depth along the longshore direction, and sudden and strong events of rip current were observed.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2017.05a
• /
• pp.787-790
• /
• 2017

This study is a preliminary research on characterization of methane - oxygen combustion used in bipropellant thruster. The limit of combustion stability and flame shape of methane - oxygen non-premixed flame injected by shear coaxial injector in the model combustion chamber Experimental studies have been carried out. A direct image of the flame was photographed using a DSLR camera, and combustion characteristics and flame length were quantified through image post-processing. As a result, it was confirmed that the stabilized flame was generated at the stoichiometric ratio as the oxidizer Reynolds number ($Re_o$) was increased, and the length of the turbulent flame was increased under the same injector diameter condition.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.31 no.1
• /
• pp.97-103
• /
• 2003

Numerical research has been done for the transverse jet behind a rearward-facings step in turbulent supersonic flow without chemical reaction. Purpose of transverse jet is to enhance mixing of the fuel in the combustor. Two-dimensional unsteady flowfields generated by slot injection into supersonic flow are numerically simulated with the Navier-Stokes equations with two-equation k-$\varepsilon$ turbulence model. Numerical method is used high-order upwind TVD scheme. Eight cases are computed for different slot momentum flux ratios and slot position at downstream of the step. The flow is very similar to the cavity flow, because the jet acts as an obstacle. The numerical results thus show the periodic phenomenon.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.11 no.2
• /
• pp.262-270
• /
• 1987

A numerical method is presented which can solve the unsteady momentum and thermal boundary layers, coupled through the agency of buoyancy force, over a heated circular cylinder impulsively started from rest. By linearizing the nonlinear finite difference equations without sacrificing accuracy, numerical solutions are obtained at each time step without iteration. To get rid of the requirement of excessive number of grid points in the region of reversed flow, special form of transformed variables are used, by which the computational boundary layer thickness is maintained almost constant. These numerical properties enable the method to easily handle the region of reversed flow and how the singularity develops in the interior of the boundary layer. In order to investigated the thermal effects on the skin friction, heat flux, displacement thickness and on the separation, we have successfully solved three different cases of the buoyancy parameter .alpha.(Gr/Re$^{2}$).

• Nuclear Engineering and Technology
• /
• v.53 no.10
• /
• pp.3344-3351
• /
• 2021

The Rare Isotope Science Project was launched in 2011 in Korea toward constructing the Rare isotope Accelerator complex for ON line experiments (RAON). RAON will house several experimental systems, including the Muon Spin Rotation/Relaxation/Resonance (μSR) facility in High Energy Experimental Building B. This facility will use 600-MeV protons with a maximum current of 660 pμA and beam power of 400 kW. The key μSR features will facilitate projects related to condensed-matter and nuclear physics. Typical experiments require a few million surface muons fully spin-polarized opposite to their momentum for application to small samples. Here, we describe the design of a muon transport beam line for delivering the requisite muon numbers and the electromagnetic-component specifications in the μSR facility. We determine the beam-line configuration via beam-optics calculations and the transmission efficiency via single-particle tracking simulations. The electromagnet properties, including fringe field effects, are applied for each component in the calculations. The designed surface-muon beamline is 17.3 m long, consisting of 2 solenoids, 2 dipoles affording 70° deflection, 9 quadrupoles, and a Wien filter to eliminate contaminant positrons. The average incident-muon flux and spin rotation angle are estimated as 5.2 × 10⁶ μ⁺/s and 45°, respectively.

• The Bulletin of The Korean Astronomical Society
• /
• v.40 no.1
• /
• pp.73.2-74
• /
• 2015

The gas cloud G2 falling toward Sagittarius A* (Sgr A*), the supermassive black hole at the center of the Milky Way, is supposed to provide valuable information on the physics of accretion flows and the environment of the black hole. We observed Sgr A* with four European stations of the Global Millimeter Very Long Baseline Interferometry Array (GMVA) at 86 GHz on 1 October 2013 when parts of G2 had already passed the pericenter. We searched for possible transient asymmetric structure - such as jets or winds from hot accretion flows - around Sgr A* caused by accretion of material from G2. The interferometric closure phases (which are zero if the spatial brightness distribution of the target is symmetric, and deviate from zero otherwise) remained zero within errors during the observation time. We thus conclude that Sgr A* did not show significant asymmetric (in the observer frame) outflows in late 2013. Using simulations, we constrain the size of the outflows that we could have missed to ${\approx}2.5$ mas along the major axis, ${\approx}0.4$ mas along the minor axis of the beam, corresponding to approximately 232 and 35 Schwarzschild radii, respectively; we thus probe spatial scales on which the jets of radio galaxies are suspected to convert magnetic into kinetic energy. As probably less than 0.2 Jy of the flux from Sgr A* can be attributed to accretion from G2, one finds an effective accretion rate ${\eta}M{\leq}1.5{\times}10^9kg/s{\approx}7.7{\times}10^{-9}M_{earth}/yr$ for material from G2. Exploiting the kinetic jet power-accretion power relation of radio galaxies, one finds that the rate of accretion of matter that ends up in jets is limited to $M{\leq}10^{17}kg/s{\approx}0.5M_{Earth}/yr$ less than about 20% of the mass of G2. Accordingly, G2 appears to be largely stable against loss of angular momentum and subsequent (partial) accretion at least on time scales ${\leq}1$ year.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.29 no.3
• /
• pp.139-146
• /
• 2017

In order to understand the efficacy of submerged breakwater constructed for the beach protection, laboratory experiments were carried out by observing the characteristics of flow around a single shore-parallel submerged breakwater. The velocity field near the shoreline was measured by utilizing the LSPIV (Large-Scale Particle Image Velocimetry) technique, and mean surface and wave height distributions were observed around the submerged breakwater, according to various combinations of incident waves and submerged breakwaters. In this experiment, it was found that the mean flow pattern behind the submerged breakwater was determined by the balance among the gradients of mean water surface and excess wave-momentum flux (i.e., radiation stress tensors) which interact with the wave-induced current developed by the gradients on the rear and the side of the submerged breakwater. The divergent and convergent flow patterns behind the submerged breakwater (i.e., accretion and erosion response) of the numerical study of Ranasinghe et al.(2010) were observed in the measured velocity distributions, and their empirical formula mostly agreed with the experimental results. However, for some cases in this experiment, it was difficult to say that the flow pattern was one of them and was agreed with the empirical formula.

• Proceedings of the Korean Vacuum Society Conference
• /
• 2014.02a
• /
• pp.221-221
• /
• 2014

The present work proposes an improved numerical simulator for design and modification of large area capacitively coupled plasma (CCP) processing chamber. CCP, as notoriously well-known, demands the tremendously huge computational cost for carrying out transient analyses in realistic multi-dimensional models, because electron dissociations take place in a much smaller time scale (${\Delta}t{\approx}10-8{\sim}10-10$) than time scale of those happened between neutrals (${\Delta}t{\approx}10-1{\sim}10-3$), due to the rf drive frequencies of external electric field. And also, for spatial discretization of electron flux (Je), exponential scheme such as Scharfetter-Gummel method needs to be used in order to alleviate the numerical stiffness and resolve exponential change of spatial distribution of electron temperature (Te) and electron number density (Ne) in the vicinity of electrodes. Due to such computational intractability, it is prohibited to simulate CCP deposition in a three-dimension within acceptable calculation runtimes (<24 h). Under the situation where process conditions require thickness non-uniformity below 5%, however, detailed flow features of reactive gases induced from three-dimensional geometric effects such as gas distribution through the perforated plates (showerhead) should be considered. Without considering plasma chemistry, we therefore simulated flow, temperature and species fields in three-dimensional geometry first, and then, based on that data, boundary conditions of two-dimensional plasma discharge model are set. In the particular case of SiH4-NH3-N2-He CCP discharge to produce deposition of SiNxHy thin film, a cylindrical showerhead electrode reactor was studied by numerical modeling of mass, momentum and energy transports for charged particles in an axi-symmetric geometry. By solving transport equations of electron and radicals simultaneously, we observed that the way how source gases are consumed in the non-isothermal flow field and such consequences on active species production were outlined as playing the leading parts in the processes. As an example of application of the model for the prediction of the deposited thickness uniformity in a 300 mm wafer plasma processing chamber, the results were compared with the experimentally measured deposition profiles along the radius of the wafer varying inter-electrode gap. The simulation results were in good agreement with experimental data.

• Korean Journal of Agricultural and Forest Meteorology
• /
• v.12 no.1
• /
• pp.45-62
• /
• 2010

Zero plane displacement (d) is the elevated height of the apparent momentum sink exerted by the vegetation on the air. For a vegetative canopy, d depends on the roughness structure of a plant canopy such as leaf area index, canopy height and canopy density, and thus is critical for the analysis of canopy turbulence and the calculation of surface scalar fluxes. In this research note, we estimated d at the Gwangneung coniferous forest by employing two independent methods of Rotach (1994) and Martano (2000), which require only a single-level eddy-covariance measurement. In general, these two methods provided comparable estimates of $d/h_c$ (where $h_c$ is the canopy height, i.e., ~23m), which ranged from 0.51 to 0.97 depending on wind directions. These estimates of $d/h_c$ were within the ranges (i.e., 0.64~0.94) reported from other forests in the literature but were sensitive to the forms of the nondimensional functions for atmospheric stability. Our finding indicates that one should be careful in interepreation of zero plane displacement estimated from a single-level eddy covariance measurement that is conductaed within the roughness sublayer.