• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.193-200
• /
• 2009

Recently, there have been efforts to construct hybrids among the existing methodologies for multiphase flow such as VOF, Level Set, and Front Tracking with the intention of facilitating simulations of general three-dimensional problems. As one of the hybrid method, we have developed the Level Contour Reconstruction Method (LCRM) for general three-dimensional multiphase flows including phase change. The main idea was focused on simplicity and a robust algorithm especially for the three-dimensional case. It combines characteristics of both Front Tracking and Level Set methods. While retaining an explicitly tracked interface using interfacial elements, the calculation of a vector distance function plays a crucial role in the periodic reconstruction of the interface elements in the LCRM method to maintain excellent mass conservation and interface fidelity. In addition, compact curvature formulation is incorporated for the calculation of the surface tension force thereby reducing parasitic currents to a negligible level.

• Proceedings of the Korean Statistical Society Conference
• /
• 2003.10a
• /
• pp.9-14
• /
• 2003

A Bayesian model-based clustering method is proposed for clustering objects on the basis of dissimilarites. This combines two basic ideas. The first is that tile objects have latent positions in a Euclidean space, and that the observed dissimilarities are measurements of the Euclidean distances with error. The second idea is that the latent positions are generated from a mixture of multivariate normal distributions, each one corresponding to a cluster. We estimate the resulting model in a Bayesian way using Markov chain Monte Carlo. The method carries out multidimensional scaling and model-based clustering simultaneously, and yields good object configurations and good clustering results with reasonable measures of clustering uncertainties. In the examples we studied, the clustering results based on low-dimensional configurations were almost as good as those based on high-dimensional ones. Thus tile method can be used as a tool for dimension reduction when clustering high-dimensional objects, which may be useful especially for visual inspection of clusters. We also propose a Bayesian criterion for choosing the dimension of the object configuration and the number of clusters simultaneously. This is easy to compute and works reasonably well in simulations and real examples.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.11
• /
• pp.827-833
• /
• 2009

The present study numerically investigates three-dimensional laminar flow past a rotating circular cylinder placed in a uniform stream. For the purpose of a careful analysis of the modification of flow by the effect of the rotation on the flow, numerical simulations are performed at a various range of rotational coefficients ($0{\leq}{\alpha}{\leq}2.5$) at one Reynolds number of 300. As ${\alpha}$ increases, flow becomes stabilized and finally a steady state beyond the critical rotational coefficient. The 3D (three dimensional) wake mode of the stationary cylinder defined at this Reynolds number has been disorganized according to ${\alpha}$, which were observed by the visualization of 3D vortical structures. The variation of the Strouhal number is very weak when the wake pattern is changed according to the rotational coefficient. As ${\alpha}$ increases, the lift increases, whereas the drag decreases.

• Journal of Information Display
• /
• v.3 no.2
• /
• pp.13-18
• /
• 2002

Numerical simulation is one of the most useful tools to study gas discharge phenomena that occur in alternating current plasma display panel (AC-PDP) cell. Most PDP cell simulations have been performed for two-dimensional cell, is cross-section along the address electrode. We developed a three-dimensional PDP simulator and applied it to a T-shaped electrode cell in order to show the effects of sustain electrode shape that cannot be included in two-dimensional simulation. The dependence of power consumption on electrode shape and area in the simulation showed the same trend as experiment.

• Proceedings of the KSME Conference
• /
• 2008.11b
• /
• pp.2733-2737
• /
• 2008

The present study numerically investigates three-dimensional laminar flow past a rotating circular cylinder placed in a uniform stream. For the purpose of a careful analysis of the modification of flow by the effect of the rotation on the flow, numerical simulations are performed at a various range of rotational speeds($0{\leq}{\alpha}{\leq}2.5$) at one Reynolds number of 300. As $\alpha$ increases, flow becomes stabilized and finally a steady state beyond the critical rotational speed. The 3D (three dimensional) wake mode of the stationary cylinder defined at this Reynolds number has been disorganized according to $\alpha$, which were observed by the visualization of 3D vortical structures. The variation of the Strouhal number is significant when the wake pattern is changed according to the rotational speed. As $\alpha$ increases, the lift increases, whereas the drag decreases.

• Bulletin of the Korean Chemical Society
• /
• v.22 no.8
• /
• pp.807-815
• /
• 2001

Two-dimensional (2D) correlation method, which generates the synchronous and the asynchronous 2D spectrum by complex cross correlation of the Fourier transformed spectra, is an analysis method for the changes of the sample spectrum induced by vari ous perturbations. In the present work, the 2D electronic correlation spectra have been simulated for the cases where the sample spectrum composed of two gaussian bands changes linearly. When only the band amplitudes of the sample spectrum change, the synchronous spectrum shows strong peaks at the band centers of the sample spectrum, but the asynchronous spectrum does not make peaks. When the sample spectrum shifts without changing intensity and width, the synchronous spectrum shows peaks around the initial and final positions of the band maximum and the asynchronous spectrum shows long peaks spanning the shifting range. The band width change produces the complex 2D correlation spectra. When the sample spectrum shifts with band broadening, the width change by 50% of full width at half maximum (FWHM) does not give so large an effect on the correlation spectrum as the spectral shift by one half of FWHM of the sample spectrum.

• Advances in aircraft and spacecraft science
• /
• v.6 no.4
• /
• pp.333-348
• /
• 2019

Low fuel regression rate is the main drawback of hybrid rocket which should be overcome. One of the improvement techniques to this problem is usage of a solid fuel grain with a complicated geometry port, which has been promoted owing to the recent development of additive manufacturing technologies. In the design of a hybrid rocket fuel grain with a complicated geometry port, the understanding of fuel regression behavior is very important. Numerical investigations of fuel regression behavior requires a capturing method of solid fuel surface, i.e. gas-solid interface. In this study, level set method is employed as such a method and the preliminary numerical tool for capturing a hybrid rocket solid fuel surface is developed. At first, to test the adequacy of the numerical modeling, the simulation results for circular port are compared to the experimental results in open literature. The regression rates and oxidizer to fuel ratios show good agreements between the simulations and the experiments, after passing enough time. However, during the early period of combustion, there are the discrepancies between the simulations and the experiments, owing to transient phenomena. Second, the simulations of complicated geometry ports are demonstrated. In this preliminary step, a star shape is employed as complicated geometry of port. The slot number effect in star port is investigated. The regression rate decreases with increasing the slot number, except for the star port with many slots (8 slots) in the latter half of combustion. The oxidizer to fuel ratio increases with increasing the slot number.

• Nuclear Engineering and Technology
• /
• v.44 no.7
• /
• pp.709-718
• /
• 2012

The pressurized water reactor APR1400 adopts DVI (Direct Vessel Injection) for the emergency cooling water in the upper downcomer annulus. The International Standard Problem number 50 (ISP-50) was launched with the aim to investigate thermal hydraulic phenomena during a 50% DVI line break scenario with best estimate codes making use of the experimental data available from the ATLAS facility located at KAERI. The present work describes the calculation results obtained for the ISP-50 using the RELAP5/MOD3.3 system code. The work aims at validation and assessment of the code to reproduce the observed phenomena and investigate about its limitations to predict complicated mixing phenomena between the subcooled emergency cooling water and the two-phase flow in the downcomer. The obtained results show that the overall trends of the main test variables are well reproduced by the calculations. In particular, the pressure in the primary system show excellent agreement with the experiment. The loop seal clearance phenomenon was observed in the calculation and it was found to have an important influence on the transient progression. Moreover, the collapsed water levels in the core are accurately reproduced in the simulations. However, the drop in the downcomer level before the activation of the DVI from safety injection tanks was underestimated due to multi-dimensional phenomena in the downcomer that are not properly captured by one-dimensional simulations.

• Korea-Australia Rheology Journal
• /
• v.19 no.3
• /
• pp.141-146
• /
• 2007

We examine the wall contact of a $3\;{\mu}m$ tethered DNA chain's free end under shear with a focus on developing schemes for double-tethering in the application of making scaffolds for molecular wires. At this scale our results are found to be highly dependent on small length scale rigidity. Chain-end-wall contact frequency, mean fractional extension deficit upon contact, and standard deviation in extension upon contact are examined for scaling with dimensionless flow strength, Wi. Predictions made using a one dimensional approximation to the Smoluchowski equation for a dumbbell and three dimensional dumbbell simulations produce extension deficit, standard deviation, and frequency scaling exponents of -1/3, -1/3, and 2/3, respectively whereas more fine-grained Kratky-Porod (KP) simulations produce scaling exponents of -0.48, -0.42, and 0.76. The contact frequency scaling of 2/3 is derived from the known results regarding cyclic dynamics Analytical scaling predictions are in agreement with those previously proposed for ${\lambda}-DNA$. [Ladoux and Doyle, 2000, Doyle et al., 2000]. Our results suggest that the differences between the dumbbell and the KP model are associated with the addition of chain discretization and the correct bending potential in the latter. These scaling results will aide future exploration in double tethering of DNA to a surface.

• Wind and Structures
• /
• v.19 no.2
• /
• pp.219-232
• /
• 2014

The characteristics of turbulent boundary layers over hilly terrain depend strongly on the hill slope and upstream condition, especially inflow turbulence. Numerical simulations are carried out to investigate the neutrally stratified turbulent boundary layer over two-dimensional hills. Two kinds of hill shape, a steep one with stable separation and a low one without stable separation, two kinds of inflow condition, laminar turbulent, are considered. An auxiliary simulation, based on the local differential quadrature method and recycling technique, is performed to simulate the inflow turbulence be imposed at inlet boundary of the turbulent inflow, which preserves very well in the computational domain. A large separation bubble is established on the leeside of the steep hill with laminar inflow, while reattachment point moves upstream under turbulent inflow condition. There is stable separation on the side of low hill with laminar inflow, whilw not turbulent inflow. Besides increase of turbulence intensity, inflow can efficiently enhance the speedup around hills. So in practice, it is unreasonable to study wind flow over hilly terrain without considering inflow turbulence.