• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.90-95
• /
• 1998

The two-dimensional incompressible and compressible Navier-Stokes codes are developed for the computation of the viscous turbulent flow over high-lift airfoils. Incompressible code using pseudo-compressibility and dual-time stepping method involves a conventional upwind differencing scheme for the convective terms and LU-SGS scheme for time integration. Compressible code also adopts an FDS scheme and LU-SGS scheme. Several two-equation turbulence models (the standard $k-{\varepsilon}$ model, the $k-{\omega}$ model. and $k-{\omega}$ SST model) are evaluated by computing the flow over single and multi-element airfoils. The compressible and incompressible codes are validated by computing the flow around the transonic RAE2822 airfoil and the NACA4412 airfoil, respectively. Both the results show a good agreement with experimental surface pressure coefficients and velocity profiles in the boundary layers. Also, the GA(W)-1 single airfoil and the NLR7301 airfoil with a flap are computed using the two-equation turbulence models. The grid systems around two- and three-element airfoil are efficiently generated using Chimera grid scheme, one of the overlapping grid generation methods.

• Journal of ILASS-Korea
• /
• v.13 no.1
• /
• pp.42-50
• /
• 2008

In the present work, atomization characteristics of shear-thinning liquid slurry discharging from fan spray nozzles were studied experimentally for spray painting applications. The effects of solid particle size and concentration on the properties (especially on the viscosity) of suspensions and mean drop size were examined by using model fluids. In the range of low particle concentration (below 3 wt%), the fluid viscosity was primarily determined by the particle size. On the other hand, in the range of high particle concentration (higher than 10 wt%), the agglomeration phenomenon and the oil absorption capability of solid particles played major roles in determining the fluid viscosity. In the high concentration region, which most of the paints belong to, the fluid became more viscous and the shear thinning behavior appeared more prominent as the particle concentration was increased. In this region, mean drop size(SMD) decreased more rapidly with the increase of the injection pressure. Also, SMD became larger with the higher particle concentration and the larger particle size.

• Economic and Environmental Geology
• /
• v.54 no.1
• /
• pp.105-115
• /
• 2021

The geological CO2 sequestration in underground geological formation such as deep saline aquifers and depleted hydrocarbon reservoirs is one of the most promising options for reducing the atmospheric CO2 emissions. The process in geological CO2 sequestration involves injection of supercritical CO2 (scCO2) into porous media saturated with pore water and initiates CO2 flooding with immiscible displacement. The CO2 migration and distribution, and, consequently, the displacement efficiency is governed by the interaction of fluids. Especially, the viscous force and capillary force are controlled by geological formation conditions and injection conditions. This study aimed to estimate the effects of surfactant on interfacial tension between the immiscible fluids, scCO2 and porewater, under high pressure and high temperature conditions by using a pair of proxy fluids under standard conditions through pendant drop method. It also aimed to observe migration and distribution patterns of the immiscible fluids and estimate the effects of surfactant concentrations on the displacement efficiency of scCO2. Micromodel experiments were conducted by applying n-hexane and deionized water as proxy fluids for scCO2 and porewater. In order to quantitatively analyze the immiscible displacement phenomena by n-hexane injection in pore network, the images of migration and distribution pattern of the two fluids are acquired through a imaging system. The experimental results revealed that the addition of surfactants sharply reduces the interfacial tension between hexane and deionized water at low concentrations and approaches a constant value as the concentration increases. Also it was found that, by directly affecting the flow path of the flooding fluid at the pore scale in the porous medium, the surfactant showed the identical effect on the displacement efficiency of n-hexane at equilibrium state. The experimental observation results could provide important fundamental information on immiscible displacement of fluids in porous media and suggest the potential to improve the displacement efficiency of scCO2 by using surfactants.

• The KSFM Journal of Fluid Machinery
• /
• v.10 no.1 s.40
• /
• pp.56-63
• /
• 2007

Because the mixing efficiency is influenced remarkably by varying the geometrical configurations, the study of flow characteristics inside the mechanical agitator is very important to improve the performances. The draught tube in the agitator makes intermixing between the screw and tube by interrupting radial flow, and it makes circulation region in a mixing chamber. In general, the helical screw agitator with a draught tube (HSA) is proved more efficient to mix than the others. Consequently, such as the shapes of helical screw, number of pitches and the variation of angular velocity are the main parameters for improving the capacity of HSA. And also the suspension of the solid particles in the agitator can be determined these parameters. The rate of solids suspension in the mixing chamber was quantified with a statistical average value, of. Numerical analyses were carried out, using a commercial CFD code, Fluent, to obtain the velocity, pressure and particle distributions under steady, laminar flow and no-slip conditions. Results are graphically depicted with various parameters.

• 한국전산유체공학회:학술대회논문집
• /
• 2000.05a
• /
• pp.39-44
• /
• 2000

Two-dimensional Navier-Stokes codes are used to simulate the shock tunnel in Konkuk university. In order to design experiments in impulse facilities properly and to interpret data from such facilities, it is necessary to understand how the flow approaches steady state. This is done by determining the transient flow field and flow establishment time around a given model. This will be accomplished by developing appropriate CFD codes which solve the Navier-Stokes equations, and simulating the starting process and resulting unsteady viscous flow phenomena. The starting process in a shock tunnel consists of multiple shock interactions and contact discontinuities, which are difficult to solve with the classical shock capturing schemes. A recently developed high resolution scheme is adapted for resolving the unsteady phenomena of those multiple shock interactions and contact surfaces during the starting process. The bifurcation phenomenon due to the interactions of the reflected shock from the end of the shock tube with the boundary layer generated by the incident shock becomes of particular interest. By comparing with the experiment results, the accuracy of the numerical analysis is validated and it is demonstrated that the properties which can hardly be obtained through the experiment can be estimated.

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.389-396
• /
• 2010

In this paper, a new computational code was developed using Chorin's artificial compressibility method to solve the two-dimensional incompressible Navier-Stokes equations. In spatial derivatives, Roe's flux difference splitting was used for the inviscid flux, while central differencing was used for the viscous flux. Furthermore, AF-ADI with dual time stepping method was implemented for accurate unsteady computations. Two-equation turbulence models, Menter's $k-{\omega}$ SST model and Coakley's $q-{\omega}$ model, hae been adopted to solve high-Reynolds number flows. A number of numerical simulations were carried out for steady laminar and turbulent flow problems as well as unsteady flow problem. The code was verified and validated by comparing the results with other computational results and experimental results. The results of numerical simulations showed that the present developed code with the artificial compressibility method can be applied to slve steady and unsteady incompressible flows.

• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.604-609
• /
• 2001

This numerical analysis uses the lifting surface method and frequency-domain panel method based on the linear compressible aerodynamic theory. Increased knowledge of flow conditions within mixed-flow fan should indicates means of improving performance of these turbomachines. Thus, only an approximate solution is obtained whose prime intent is to recognize the most significant characteristics of the "ideal" geometry. For a given set of operating condition, the flow conditions within mixed-flow fan depend on the geometry of the machine (three-dimensional flow effects) and on the properties of the fluid. But most treatments of the problem have been concerned with the two-dimensional flow effects for incompressible, non-viscous fluids. Interest in the field of mixed-flow fan resulted in the undertaking of a program to develop reliable design procedures that would avoid the need for lengthy development work.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03b
• /
• pp.354-357
• /
• 2008

Drag reduction is one of main concerns for commercial aircraft companies than ever because fuel price has been tripled in ten years. In this research, Natural Laminar Flow airfoil is designed and tested to reduce drag at cruise condition, $c_l$=0.3, Re=3.4${\times}$10$^6$ and M=0.6. NLF airfoil is characterized by delayed transition from laminar to turbulent flow, which comes from maintaining favorable pressure gradient to downstream. Transition is predicted by solving Boundary Layer equations in viscous boundary layer and by solving Euler Equation outside the boundary layer. Once boundary layer thickness and momentum thickness are obtained, $e^N$-method is used for transition point prediction. As results, KARI's NLF airfoil is designed and shows better characteristics than NLF-0115. The characteristics are tested and verified at low Reynolds numbers, but at high Reynolds numbers, laminar flow characteristics are not obtainable because of fully turbulent flow over airfoil surfaces. Precious experiences, however, relating NLF airfoil design, subsonic and transonic tests are acquired.

• Journal of Magnetics
• /
• v.22 no.2
• /
• pp.181-187
• /
• 2017

Flow of electrically conducting nanofluids is of pivotal importance in countless industrial and medical appliances. Fluctuations in thermophysical properties of such fluids due to variations in temperature have not received due attention in the available literature. Present investigation aims to fill this void by analyzing the flow of copper-water nanofluid with temperature dependent viscosity past a Riga plate. Strong wall suction and viscous dissipation have also been taken into account. Numerical solutions for the resulting nonlinear system have been obtained. Results are presented in the graphical and tabular format in order to facilitate the physical analysis. An estimated expression for skin friction coefficient and Nusselt number are obtained by performing linear regression on numerical data for embedded parameters. Results indicate that the temperature dependent viscosity alters the velocity as well as the temperature of the nanofluid and is of considerable importance in the processes where high accuracy is desired. Addition of copper nanoparticles makes the momentum boundary layer thinner whereas viscosity parameter does not affect the boundary layer thickness. Moreover, the regression expressions indicate that magnitude of rate of change in effective skin friction coefficient and Nusselt number with respect to nanoparticles volume fraction is prominent when compared with the rate of change with variable viscosity parameter and modified Hartmann number.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.57-65
• /
• 2009

An efficient and high-fidelity design approach for wing-body shape optimization is presented. Depending on the size of design space and the number of design of variable, aerodynamic shape optimization process is carried out via different optimization strategies at each design stage. In the first stage, global optimization techniques are applied to planform design with a few geometric design variables. In the second stage, local optimization techniques are used for wing surface design with a lot of design variables to maintain a sufficient design space with a high DOF (Degree of Freedom) geometric change. For global optimization, Kriging method in conjunction with Genetic Algorithm (GA) is used. Asearching algorithm of EI (Expected Improvement) points is introduced to enhance the quality of global optimization for the wing-planform design. For local optimization, a discrete adjoint method is adopted. By the successive combination of global and local optimization techniques, drag minimization is performed for a multi-body aircraft configuration while maintaining the baseline lift and the wing weight at the same time. Through the design process, performances of the test models are remarkably improved in comparison with the single stage design approach. The performance of the proposed design framework including wing planform design variables can be efficiently evaluated by the drag decomposition method, which can examine the improvement of various drag components, such as induced drag, wave drag, viscous drag and profile drag.