• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.10
• /
• pp.778-782
• /
• 2009

Recently, we developed a simple thermopneumatic micropump having neither a membrane nor a valve. This micropump discharges liquid by a thermopneumatic pressure and refill by a capillary attraction. In case of the micropump driven by the capillary attraction, the flow characteristic depends mainly on the geometry of the micropump. In this paper, we investigated the influence of the geometry of the micropump on the performance of the micropump to illustrate the properness of the micropump shape. We analyzed the micropump characteristics of six types having different geometries by FVM simulation with a commercial CFD tool. Also we fabricated the micropumps with PDMS and glass by micromachining, and tested the performances. The simulation and the test results illustrate that the discharge volume and the discharge time depend on the chamber volume. The expansion angle of the inlet channel location has influence on the refill time, while the front air channel direction has influence on the backward flow loss.

• Proceedings of the Korean Society of Marine Engineers Conference
• /
• 2005.06a
• /
• pp.1149-1154
• /
• 2005

The finite difference lattice Boltzmann method(FDLBM) is a quite recent approach for simulating fluid flow, which has been proven as a valid and efficient tool in a variety of complex flow problems. It is considered an attractive alternative to conventional FDM and FVM, because it recovers the Navier-Stokes equations and is computationally more stable, and easily parallelizable to simulate for various laminar flows and a direct simulation of aerodynamics sounds. However, the research of a numerical simulation of turbulent flow by FDLBM, which is important to analyze the structure of turbulent flow in engineering fields, is not carried out. In this research, the FDLBM built in the turbulent model is applied, and a flowfield around 2-dimensional square to validate the applied model with 2D9V is simulated. Besides, 2D computation of the cavity noise generated by flow over a cavity at a Mach number of 0.1 and a Reynolds number based on cavity depth of 5000 is calculated. The computation result is well presented a understanding of the physical phenomenon of tonal noise occurred primarily by well-jet shear layer and vortex shedding and an aeroacoustic feedback loop.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.916-921
• /
• 2003

In most industrial applications, the geometrical complexity is combined with the moving boundaries. These problems considerably increase the computational difficulties since they require, respectively, regeneration and deformation of the grid. As a result, engineering flow simulation is restricted. In order to solve this kind of problems the immersed boundary method was developed. In this study, the immersed boundary method is applied to the numerical simulation of stationary, rotating and oscillating cylinders in the 2-dimensional square cavity. No-slip velocity boundary conditions are given by imposing feedback forcing term to the momentum equation. Besides, this technique is used with a second-order accurate interpolation scheme in order to improve the accuracy of flow near the immersed boundaries. The governing equations for the mass and momentum using the immersed boundary method are discretized on the non-staggered grid by using the finite volume method(FVM). This study presents the possibility of the immersed boundary method to apply to the complex flow experienced in the industrial applications.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.177-180
• /
• 2008

The objective of this study is to get simulation data about pulsatile flow of a non-Newtonian fluid through a bifurcated tube. All the process was based on CFD method, with a commercial FVM code, SC/Tetra ver. 6.0 for solving, and with CATIA R16 for generating geometries. To define a non-Newtonian fluid, the following viscous models are used; the Powell-Eyring model, the modified Powell-Eyring model, the Cross model, the modified Cross model, the Carreau model, the Carreau-Yasuda model and the modified Power Law model. The flow calculation data using each model were compared with the other data of a existing paper. Finally, the Carreau model was recognized to give the best result with the SC/Tetra code, and the succeeding simulations are made with the model. For the pulsating flow condition, the sine wave type velocity profile is given as the inlet boundary condition. To investigate the effect of geometries and mesh, the pre-test is carried out with various curvature conditions of the bifurcated corner, and then with various mesh conditions. The final process is to calculate flow variables such as the wall shear stress (WSS) and the wall shear stress gradient (WSSG). To validate all the result, the simulation is compared with the existing data of the other papers. Generally speaking, there is a noticeable difference in the maximum and minimum value of WSS. It is not sure that the values in each data are on the exactly same location. However, the overall trend is similar. The next study needs to investigate the same situation by experimental method. Furthermore, if the flow is simulated with more pulsatile conditions, more data of flow field through a bifurcated tube could be achieved.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.177-180
• /
• 2008

The objective of this study is to get simulation data about pulsatile flow of a non-Newtonian fluid through a bifurcated tube. All the process was based on CFD method, with a commercial FVM code, SC/Tetra ver. 6.0 for solving, and with CATIA R16 for generating geometries. To define a non-Newtonian fluid, the following viscous models are used; the Powell-Eyring model, the modified Powell-Eyring model, the Cross model, the modified Cross model, the Carreau model, the Carreau-Yasuda model and the modified Power Law model. The flow calculation data using each model were compared with the other data of a existing paper. Finally, the Carreau model was recognized to give the best result with the SC/Tetra code, and the succeeding simulations are made with the model. For the pulsating flow condition, the sine wave type velocity profile is given as the inlet boundary condition. To investigate the effect of geometries and mesh, the pre-test is carried out with various curvature conditions of the bifurcated corner, and then with various mesh conditions. The final process is to calculate flow variables such as the wall shear stress (WSS) and the wall shear stress gradient (WSSG). To validate all the result, the simulation is compared with the existing data of the other papers. Generally speaking, there is a noticeable difference in the maximum and minimum value of WSS. It is not sure that the values in each data are on the exactly same location. However, the overall trend is similar. The next study needs to investigate the same situation by experimental method. Furthermore, if the flow is simulated with more pulsatile conditions, more data of flow field through a bifurcated tube could be achieved.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2004.05a
• /
• pp.222-228
• /
• 2004

The effects of rotation on the unsteady laminar flow past a circular cylinder is numerically investigated in the present study. We obtained the numerical solutions for unsteady two-dimensional governing equation for the flow using two different numerical schemes. One is an accurate spectral method and another is finite volume method. Above all, the flow around a stationary circular cylinder is investigated to understand the basic phenomenon of flow separation, bluff body wake. Also, the validation of our own codes, expecially based on FVM, is carried out by the comparison of results obtained from our simulations using two different schemes and previous numerical and experimental studies. By the effect of rotation, the mean lift increases and drag deceases, which well represent the previous study.

• Journal of the Korea Computer Graphics Society
• /
• v.15 no.1
• /
• pp.1-6
• /
• 2009

This paper presents a new CIP method for unstructured mesh to reduce the numerical dissipation. To reflect precise physical characteristics, CIP method updates both the physical quantity and the derivative information. The proposed method uses the Finite Volume Method(FVM) to solve the non-advection term of CIP equation. And we performed several experiments to improve the accuracy of third-order interpolation. Our result shows that our algorithm has less numerical dissipation than that of linear advection solver.

• Journal of the Korean Society for Aviation and Aeronautics
• /
• v.18 no.1
• /
• pp.33-38
• /
• 2010

Developments of numerical methods are very important to design and analysis for a high subsonic turbine blade. In general, Analysis by experimental investigation has needed a lot of human resources and required time, indispensably, and equipments still have a limit to measure in conditions of high temperature. Rapid technological developments of CPU and integration level of memory make it possible to advance computation with almost exactly simulation so, recent developments of numerical methods are in spotlight. In the present study, the panel method, which is well-known as relatively simplified numerical method, and 2-dimensional ordinary differential Falkner-Skan equation were computed in order to analyze the outer flow, and FVM-based solid heat transfer equation, was also computed to forecast the temperature distribution of the airfoil and the turbine blade. Unstructured grid was constructed in the turbine blade, which has double cooling holes, in order to analyze the internal heat transfer. Cooling fluid was assumed as fully-developed turbulent flow and that circulated in cooling holes.

• Korea-Australia Rheology Journal
• /
• v.16 no.1
• /
• pp.47-54
• /
• 2004

Low Reynolds number, dilute, and surfactant-free bubble suspensions are prepared by mechanical mixing after introducing carbon dioxide bubbles into a Newtonian liquid, polyol. The apparent shear viscosity is measured with a wide-gap parallel plate rheometer by imposing a simple shear flow of capillary numbers(Ca) of the order of $10^{-2}$ ~ $10^{-1}$ and for various gas volume fractions ($\phi$). Effects of capillary numbers and gas volume fractions on the viscosity of polyol foam are investigated. At high capillary number, viscosity of the suspension increases as the gas volume fraction increases, while at low capillary number, the viscosity decreases as the gas volume fraction increases. An empirical constitutive equation that is similar to the Frankel and Acrivos equation is proposed by fitting experimental data. A numerical simulation for deformation of a single bubble suspended in a Newtonian fluid is conducted by using a newly developed two-dimensional numerical code using a finite volume method (FVM). Although the bubble is treated by a circular cylinder in the two dimensional analysis, numerical results are in good agreement with experimental results.

• Journal of Korea Foundry Society
• /
• v.27 no.6
• /
• pp.237-242
• /
• 2007

To improve the quality of the product and the cost efficiency, the joining of A356 alloy to an Al-18wt%Si alloys has been performed by centrifugal casting. The influence of the mold preheating temperature, the pouring temperature and the rotational velocity of the mold on the microstructures of the shell in the centrifugal casting was investigated using the experimental and simulation methods. In the present study, the cellular automaton (CA) technique and the finite volume method (FVM) were adopted to simulate the evolution of the macro structures and to calculate the temperature profiles, respectively. The evolution of the microstructures was also simulated using a modified cellular automaton (MCA) model. The optimal rotational speed of the mold for obtaining the sound shape of the shell was estimated experimentally to be over 1200 rpm. For the uniform microstructure, the outer shell needs to be cast with higher preheated mold temperature and lower pouring temperature, and the melt was poured at lower temperature in the inner shell. In order to obtain the sound shape of the joining, the different materials were poured simultaneously.