• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.177-185
• /
• 1999

The airfoil design optimization procedures based on the Navier-Stokes equations were developed, This procedure enables more realistic and practical transonic airfoil designs. The modified Hicks-Henne functions were used to generate the shape of airfoils. Five Hick-Henne functions were used to design upper surface of airfoil only. To enhance the ability of Hick-Henne function to generate various airfoil shape with limited number of functions, the positions of control points were adjusted through optimization procedure. The design procedure was applied to the single-point design for the drag minimization problem with lift and area constraints. The result shows the capability of the procedure to generate much realistic airfoils with very small drag-creep in the low transonic regime. This is mainly due to the viscosity effect of Navier-Stokes flow analysis. However, in the higher transonic range tile drag-creep appears. The multi-point design is shown to be an effective way to avoid the drag-creep and improve off-design performance which is very similar in the Euler design.

• Journal of computational fluids engineering
• /
• v.12 no.4
• /
• pp.14-19
• /
• 2007

We present a direct numerical simulation technique and some preliminary results of the capillary spreading of a droplet containing particles on the solid substrate. We used the level-set method with the continuous surface stress for description of droplet spreading with interfacial tension and employed the discontinuous Galerkin method for the stabilization of the interface advection equation. The distributed Lagrangian-multipliers method has been combined for the implicit treatment of rigid particles. We investigated the droplet spreading by the capillary force and discussed effects of the presence of particles on the spreading behavior. It has been observed that a particulate drop spreads less than the pure liquid drop. The amount of spread of a particulate drop has been found smaller than that of the liquid with effectively the same viscosity as the particulate drop.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2010.11a
• /
• pp.19-23
• /
• 2010

Injection characteristics of a kerosene swirl injector of liquid rocket engine operating at supercritical environment have been investigated. Kerosene surrogate models are proposed to model the kerosene properties. Turbulent numerical model is based on large eddy simulation and contains Soave modification of Redlich-Kwong equation of state and Chung's model. Numerical analysis results at supercritical environment are compared with the one at transcritical condition. Differences of density and viscosity are analyzed at both liquid film and core gas in the swirl injector.

• 한국항공운항학회:학술대회논문집
• /
• 2004.11a
• /
• pp.240-246
• /
• 2004

For the prediction on the onset of oblique shock wave-induced vortex breakdown, computational studies on the Oblique Shock wave/Vortex Interaction (OSVI) are conducted and compared with both experimental results and analytic model. A Shock-stable numerical scheme, the Roe scheme with Mach number-based function (RoeM), and a two-equation eddy viscosity-transport approach are used for three-dimensional turbulent flow computations. The computational configuration is identical to available experiment, and we attempt to ascertain the effect of parameters such as a vertex strength, streamwise velocity deficit, and shock strength at a freestream Mach number of 2.49. Numerical simulations using the ${\kappa}-{\omega}SST$ turbulence model and suitably modeled vortex profiles are able to accurately reproduce many fine features through a direct comparison with experimental observations. The present computational approach to determine the criterion on the onset of oblique shock wave-induced vortex breakdown is found to be in good agreement with both the experimental result and the analytic prediction.

• Korea-Australia Rheology Journal
• /
• v.17 no.4
• /
• pp.199-205
• /
• 2005

Carbon dioxide was absorbed into aqueous polyethylene oxide (PEO) solution containing AMP in a flat-stirred vessel to investigate the effect of non-Newtonian rheological behavior of PEO on the rate of chemical absorption of $CO_2$, where the reaction between $CO_2$ and AMP was assumed to be a first-order reaction with respect to the molar concentration of $CO_2$ and AMP respectively. The liquid-side mass transfer coefficient, which was obtained from the dimensionless empirical equation containing the properties of viscoelasticity of the non-Newtonian liquid, was used to estimate the enhancement factor due to chemical reaction. PEO with elastic property of non-Newtonian liquid made the rate of chemical absorption of $CO_2$ accelerate compared with Newtonian liquid based on the same viscosity of the solution.

• Proceedings of the Korean Society of Marine Engineers Conference
• /
• 2001.05a
• /
• pp.52-57
• /
• 2001

Laminar two-dimensional time-dependent flow past a circular cylinder is numerically investigated using direct numerical simulation for the low Reynolds number (Re=164∼280). The higher-order finite difference scheme is employed for the spatial distributions along with the second order Adams-Bashforth and the first order backward-Euler time integration. The convection term is applied by the 7th order up wind scheme and the pressure and viscosity terms are applied by the 4th order central difference. The grid system makes use of the regular grid system and it is generated by an equation. The calculated results of drag coefficients, lift coefficients, pressure distributions, and vorticity contours and other information are compared with experimental and numerical ones. These results obtained by the present DNS show good agreement with the previous studies.

• Journal of the Korean earth science society
• /
• v.40 no.4
• /
• pp.392-405
• /
• 2019

Biogeochemical processes play an important role in ocean environments and can affect the entire Earth's climate system. Using an ocean-biogeochemistry model (NEMO-TOPAZ), we investigated the effects of changes in albedo and wind stress caused by phytoplankton in the equatorial Pacific. The simulated ocean temperature showed a slight decrease when the solar reflectance of the regions where phytoplankton were present increased. Phytoplankton also decreased the El $Ni{\tilde{n}}o$-Southern Oscillation (ENSO) amplitude by decreasing the influence of trade winds due to their biological enhancement of upper-ocean turbulent viscosity. Consequently, the cold sea surface temperature bias in the equatorial Pacific and overestimation of the ENSO amplitude were slightly reduced in our model simulations. Further sensitivity tests suggested the necessity of improving the phytoplankton-related equation and optimal coefficients. Our results highlight the effects of altered albedo and wind stress due to phytoplankton on the climate system.

• The Bulletin of The Korean Astronomical Society
• /
• v.44 no.2
• /
• pp.47.3-47.3
• /
• 2019

In this talk, we present numerical simulations of driven hydrodynamic and magnetohydrodynamic (MHD) turbulence with weak/strong imposed magnetic fields. We mainly focus on turbulence driven compressively (∇ × f = 0). Our main goal is to examine how magnetic fields play a role in generating solenoidal modes in compressive turbulence. From our simulation analysis, we find that solenoidal energy densities in hydrodynamic and weak magnetic field cases are generated up to ~ 30% of total ones. On the other hand, in the case of strong magnetic fields, solenoidal energy densities are excited up to ~ 70%. To interpret the results, we further analyze vorticity (w = ∇ × u) equation and find that magnetic fields directly create solenoidal motions, and magnetic tension is most effective in this sense. In hydrodynamic simulations, however, we find that viscous dissipation provides vorticity seeds at the very early stage and they are amplified via stretching process. Lastly, in weak magnetic fields cases, we find that solenoidal motions are created by the effects of magnetic fields, viscosity, and stretching in conjunction.

• Journal of the Korean Chemical Society
• /
• v.19 no.5
• /
• pp.289-293
• /
• 1975

According to the viscosity and density studies of an aqueous solution of dodecyl pyridinium chloride (DPC) it was shown that the 2nd critical micelle concentration (2nd cmc) existed in the solution in addition to the 1st cmc. The volume fraction ${\phi}_m$ of the hydrated micelle was calculated by using the following equation: ${\Pi}_{rm}=1+2.5{\phi}_m+14.1{\phi}_m^2$ It has been found that the increment of ${\phi}_m$ in the ${\phi}_m$ vs. $C_m$ curve decreased at near the 2nd cmc. And the partial specific volume ($\={v}$) of DPC obtained from the density measurement also decreased rapidly at near the 2nd cmc and remains constant value above the 2nd cmc. This may be attributed to a change in the micelle structure at the 2nd cmc caused by a variation in the type of aggregation and by a decrease in the counterion binding by the micelle.

• Journal of Biosystems Engineering
• /
• v.17 no.3
• /
• pp.229-236
• /
• 1992

In an effort to localize the spray-dryer which is markedly used for drying food materials, a experiment was carried out with a wheel type atomizer locally designed and manufactured to evaluate the effect of rotational speed, feed rate and physical properties of liquid food material on the droplet size, and to develop a model to predict the droplet size sprayed at various operational conditions. The result are summarized as follows. 1. The frequency of droplet size sprayed from the atomizer at every treatment were similar to normal distributions. 2. Under the test conditions adopted in this study, that is, rotational speed of the atomizer ranging from 15,000 to 20,000 rpm (55.0 m/sec - 73.3 m/sec), feed rate from 14 to 37 kg/hr and viscosity of the material from 1.14 to 350 cP, the mean volume-surface dia. of droplets was decreased as increase in rotational speed and was not affected significantly by the feed rate and viscosity. 3. Through the dimensional analysis, a prediction model was developed as follows : $$\frac{Dvs}{r}=K[\frac{Q}{{\mu}r}]^a[\frac{rN^2}{g}]^b[\frac{{\rho}^2r^3g}{{\mu}^2}]^c[\frac{L}{r}]^d$$ and it was proved that the above model was better in degree of fitness than other models reported. 4. A prediction equation for the droplet size sprayed from the atomizer under the test was expressed as follows : $$\frac{Dvs}{r}=0.0215[\frac{Q}{{\mu}r}]^{0.06}[\frac{rN^2}{g}]^{0.3314}[\frac{{\mu}^2}{{\rho}^2r^3g}]^{0.0158}$$.