• Journal of computational fluids engineering
• /
• v.11 no.4 s.35
• /
• pp.75-80
• /
• 2006

Navier-Stokes simulation of the flow in a micro viscous pump is carried out. The micro viscous pump consists of a rotating circular rotor placed in a two dimensional channel. All simulation is carried out by using a finite volume approach, at the Reynolds number of 0.5, to study the performance of the micro viscous pump. Length of channel of the pump is varied to simulate the effects of the pumping load. Numerical solutions show that the net flow of the pump is realized by two counter rotating vortices formed on both sides of the rotor. The volume flow rate of the pump is decreased as length of the channel is increased, while the static pressure difference across the rotor is increased. The static pressure difference across the rotor is observed to be inversely proportional to the volume flow rate as inertia effects are negligibly small. The efficiency of the pump is found to reach a maximum when two counter rotating vortices on both sides of the rotor becomes to merge forming an outer enveloping vortex.

• 연구논문집
• /
• s.33
• /
• pp.17-25
• /
• 2003

In this work, thermal design of a PCR chip for LOC is systematically conducted. From the numerical simulation of a PCR chip based on the finite volume method, how to control the average temperature of a PCR chip and the temperature difference between the denaturation zone and the annealing zone is presented. The average temperature is shown to be controlled by adjusting heat input and a cooler as well as a heater is shown to be necessary to obtain three individual temperature zones for polymerase chain reaction. To reduce the time required, a heat sink for the cooler is not included in the calculation domain for the PCR chip and heat sink design is conducted separately by using a compact modeling method, the porous medium approach.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.24 no.3
• /
• pp.9-16
• /
• 1987

In this paper, the fluid flow in the two-dimensional tank is analyzed by the Finite Difference Method. The Navier-Stokes equation is modified for the tank fixed coordinate system. For the treatment of the free surface, the Volume of Fluid Method by Hirt and Nichols is adopted. The continuity equation and the Poisson equation which is derived from the Navier-Stokes equation to find the pressure are solved by the Successive-Line-Overrelaxation Method. The comparison of the calculated results with experimental data show a favorable agreement. The fluid flow in the two-dimensional tank can be predicted reasonably before the free surface reaches breaking by this numerical method.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.05b
• /
• pp.245-249
• /
• 2005

The numerical simulation of a two-phase flow In a discrete fracture network model is presented in this paper, The purpose of this work is to consider density-driven flows induced by the density difference between hot autochthonous heavy brines and injected cold water. Mechanical consequences of high pressure waves on the fracture permeability and heat exchanges between fluids and rock matrix are neglected in this study. The finite volume method is employed to discretize spatially and the system is solved by using an IMPES(Implicit Pressure-Explicit Saturation) scheme. In order to solve the strong non-linearity of the system, the Newton-Raphson algorithm is used. The well-known Buckeley-Leverett problem is adapted to validate results calculated from the model and a relatively good agreement is obtained.

• Wind and Structures
• /
• v.5 no.2_3_4
• /
• pp.317-328
• /
• 2002

A numerical investigation on the turbulent flows over a three-dimensional steep hill is presented. The numerical model developed for the present work is based on the finite volume method and the SIMPLE algorithm with a non-staggered grid system. Standard $k-{\varepsilon}$ model and Shih's non-linear model are tested for the validation of the prediction accuracy in the 3D separated flow. Comparisons of the mean velocity and turbulence profiles between the numerical predictions and the measurements show good agreement. The Shih's non-linear model is found to predict mean flow and turbulence better than the Standard $k-{\varepsilon}$. Flow patterns have also been examined to explain the difference in the cavity zone between 2D and 3D hills.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.11
• /
• pp.3031-3040
• /
• 1995

Two-dimensional incompressible Navier-Stokes equations have been solved by the node-centered finite volume method with the unstructured triangular meshes. The pressure-velocity coupling is handled by the artificial compressibility algorithm due to its computational efficiency associated with the hyperbolic nature of the resulting equations. The convective fluxes are obtained by the Roe's flux difference splitting scheme using edge-based connectivities and higher-order differences are achieved by a reconstruction procedure. The time integration is based on an explicit four-stage Runge-Kutta scheme. Numerical procedures with local time stepping and implicit residual smoothing have been implemented to accelerate the convergence for the steady-state solutions. Comparisons with experimental data and other numerical results have proven accuracy and efficiency of the present unstructured approach.

• Steel and Composite Structures
• /
• v.28 no.6
• /
• pp.735-748
• /
• 2018

In this study, the semi-analytical finite strip method is adopted to examine the free vibration of cylindrical shells made up of functionally graded material. The properties of functionally graded shells are assumed to be temperature-dependent and vary continuously in the thickness direction according to a simple power law distribution in terms of the volume fraction of ceramic and metal. The material properties of the shells and stiffeners are assumed to be continuously graded in the thickness direction. Theoretical formulations based on the smeared stiffeners technique and the classical shell theory with first-order shear deformation theory which accounts for through thickness shear flexibility are employed. The finite strip method is applied to five different shell theories, namely, Donnell, Reissner, Sanders, Novozhilov, and Teng. The approximate procedure is compared favorably with three-dimensional finite elements. Finally, a detailed numerical study is carried out to bring out the effects of power-law index of the functional graded material, stiffeners, and geometry of the shells on the difference between various shell theories. Finally, the importance of choosing the shell theory in simulating the functionally graded cylindrical shells is addressed.

• Journal of Korea Water Resources Association
• /
• v.44 no.6
• /
• pp.429-438
• /
• 2011

In order to simulate a free surface flow in a trench channel, a three-dimensional incompressible unsteady Reynolds-averaged Navier-Stokes (RANS) equations are closed with the ${\kappa}-{\epsilon}$ model. The artificial compressibility (AC) method is used. Because the pressure fields can be coupled directly with the velocity fields, the incompressible Navier-Stokes (INS) equations can be solved for the unknown variables such as velocity components and pressure. The governing equations are discretized in a conservation form using a second order accurate finite volume method on non-staggered grids. In order to prevent the oscillatory behavior of computed solutions known as odd-even decoupling, an artificial dissipation using the flux-difference splitting upwind scheme is applied. To enhance the efficiency and robustness of the numerical algorithm, the implicit method of the Beam and Warming method is employed. The treatment of the free surface, so-called interface-tracking method, is proposed using the free surface evolution equation and the kinematic free surface boundary conditions at the free surface instead of the dynamic free surface boundary condition. AC method in this paper can be applied only to the hydrodynamic pressure using the decomposition into hydrostatic pressure and hydrodynamic pressure components. In this study, the boundary-fitted grids are used and advanced each time the free surface moved. The accuracy of our RANS solver is compared with the laboratory experimental and numerical data for a fully turbulent shallow-water trench flow. The algorithm yields practically identical velocity profiles that are in good overall agreement with the laboratory experimental measurement for the turbulent flow.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.214-219
• /
• 2007

Both the bubble rising in a fully filled container and the droplet splash are simulated by a solution code(PowerCFD). This code employs an unstructured cell-centered method based on a conservative pressure-based finite-volume method with interface capturing method (CICSAM) in a volume of fluid(VOF) scheme for phase interface capturing. The present results are compared with other numerical solutions found in the literature. It is found that the present code simulate complex free surface flows such as multi phase flows due to large density difference efficiently and accurately.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2009.05a
• /
• pp.458-461
• /
• 2009

Thermo-mechanical simulation of the Friction Stir Spot Welding (FSSW) processes was performed for the AA5083-H18 sheets, utilizing commercial Finite Element Method (FEM) and Finite Volume Method (FVM) which are based on Lagrangian and Eulerian formulations, respectively. The Lagrangian explicit dynamic FEM code, PAM-CRASH, and the Eulerian Computational Fluid Dynamics (CFD) FVM code, STAR-CD, were utilized to understand the effect of pin geometry on weld strength and material flow under the unsteady state condition. Using FVM code, material flow pattern near the tool boundary was analyzed to explain the weld strength difference between the weld by cylindrical pin and the weld by triangular pin, while the frictional energy concept using the FEM code had limitation to explain the weld strength difference.