• Coupled systems mechanics
• /
• v.7 no.4
• /
• pp.421-440
• /
• 2018

In this study, a computational method for wall shear stress combined with an implicit direct-forcing immersed boundary method is presented. Near the immersed boundaries, the sub-grid stress is determined by a wall model in which the wall shear stress is directly calculated from the Lagrangian force on the immersed boundary. A coupling mathematical model of the transition process for a model Francis turbine comprising turbulent flow and rotating rigid guide vanes is established. The spatiotemporal distributions of pressure, velocity, vorticity and turbulent quantity are gained with the transient process; the drag and lift coefficients as well as other forces (moments) are also obtained as functions of the attack angle. At the same time, analysis is conducted of the characteristics of pressure pulsation, velocity stripes and vortex structure at some key parts of flowing passage. The coupling relations among the turbulent flow, the dynamical force (moment) response of blade and the rotating of guide vane are also obtained.

• 한국전산유체공학회:학술대회논문집
• /
• 2002.10a
• /
• pp.84-87
• /
• 2002

When computing the flow around complex three dimensional configurations, the generation of grid is the sunt time consuming part of any calculation. The object of this study is to develop the grid duster techniques capable of resolving complex flows with shock waves, expansion waves, shear layers, and cursive shapes, The Dot insert method of Non-Uniform Rational B-Splines is described as a id control method.

• Journal of computational fluids engineering
• /
• v.19 no.3
• /
• pp.37-43
• /
• 2014

In this paper, hydraulic & thermal developing and fully developed laminar forced convection flow of a water-$Al_2O_3$ nanofluid in a circular horizontal tube with uniform heat flux at the wall, are investigated numerically. A single phase model employed with temperature independent properties. The thermal entrance length is presented in this paper. The variations of the convective heat transfer coefficient and shear stress are shown in the entrance region and fully developed region along different nanoparticles concentration and Reynolds numbers. Convective heat transfer coefficient for nanofluids is larger than that of the base fluid. It is shown that heat transfer is enhanced and shear stress is increased as the particle volume concentration increases. The heat transfer improves, as Reynolds number increases.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.17-22
• /
• 2002

Non-uniform and unsteady inflow into a Horizontal Axis Wind Turbine (HAWT) brings about an asymmetric flow field on the rotor plane and an unsteady aerodynamic load on the blades. In the present paper effects of yawed inflow and wind shear are analyzed by an inviscid aerodynamic model based on the asymptotic acceleration potential method. In the analysis the rotor blades are represented by spanwise and chordwise pressure distribution composed of analytical first-order asymptotic solutions for the Laplace equation. As the actual wind field experienced by a HAWT is turbulent, the effects of the turbulence are also examined using the Veers' model.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.15-28
• /
• 1998

As an alternative for solving the incompressible Navier-Stokes equations, we present a vorticity-based integro-differential formulation for vorticity, velocity and pressure variables. One of the most difficult problems encountered in the vorticity-based methods is the introduction of the proper value-value of vorticity or vorticity flux at the solid surface. A practical computational technique toward solving this problem is presented in connection with the coupling between the vorticity and the pressure boundary conditions. Numerical schemes based on an iterative procedure are employed to solve the governing equations with the boundary conditions for the three variables. A finite volume method is implemented to integrate the vorticity transport equation with the dynamic vorticity boundary condition . The velocity field is obtained by using the Biot-Savart integral derived from the mathematical vector identity. Green's scalar identity is used to solve the total pressure in an integral approach similar to the surface panel methods which have been well-established for potential flow analysis. The calculated results with the present mettled for two test problems are compared with data from the literature in order for its validation. The first test problem is one for the two-dimensional square cavity flow driven by shear on the top lid. Two cases are considered here: (i) one driven both by the specified non-uniform shear on the top lid and by the specified body forces acting through the cavity region, for which we find the exact solution, and (ii) one of the classical type (i.e., driven only by uniform shear). Secondly, the present mettled is applied to deal with the early development of the flow around an impulsively started circular cylinder.

• Korea-Australia Rheology Journal
• /
• v.21 no.3
• /
• pp.193-200
• /
• 2009

We present a direct simulation technique for two-dimensional mold-filling simulations of fluids filled with a large number of circular disk-like rigid particles. It is a direct simulation in that the hydrodynamic interaction between particles and fluid is fully considered. We employ a pseudo-concentration method for the evolution of the flow front and the DLM (distributed Lagrangian multipliers)-like fictitious domain method for the implicit treatment of the hydrodynamic interaction. Both methods allow the use of a fixed regular discretization during the entire computation. The discontinuous Galerkin method has been used to solve the concentration evolution equation and the rigid-ring description has been introduced for freely suspended particles. A buffer zone, the gate region of a finite area subject to the uniform velocity profile, has been introduced to put discrete particles into the computational domain avoiding any artificial discontinuity. From example problems of 450 particles, we investigated the particle motion and effects of particles on the flow for both Newtonian and shear-thinning fluid media. We report the prolonged particle movement toward the wall in case of a shear-thinning fluid, which has been interpreted with the shear rate distribution.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.1
• /
• pp.274-284
• /
• 1991

Mean flows and Reynolds stresses through circular and elliptic wire-mesh screens in the wind tunnel are measured by using the hot-wire system, and flow structures are investigated. Flow in the core of the wake are nearly uniform and the shear layer is developed along the edge of the screen The turbulent kinetic energy in the core decreases at the fast rate. However turbulence components are not in local equilibrium in the shear layer. The shear layer of the circular screen develops outward according to the radial mean motion. On the other hand, 3-dimensional transverse mean motion was turned to the main mechanism for the elliptic shape of the wake to be circular at the downstream.

• Journal of computational fluids engineering
• /
• v.12 no.3
• /
• pp.20-28
• /
• 2007

A fast and robust method of grid generation to multiple functions has been developed for flow analysis in three dimensional space. It is based on the Non-Uniform Rational B-Spline(NURBS) of an approximation method. Many of NURBS intrinsic properties are introduced and much more easily understood. The grid generation method, details of numerical implementation. examples of application, and potential extensions of the current method are illustrated in this paper. The object of this study is to develop the surface grid generation and the grid cluster techniques capable of resolving complex flows with shock waves, expansion waves, shear layers. The knot insert method of Non-Uniform Rational B-Spline seems well worked. In addition, NURBS has been widely utilized to generate grids in the computational fluid dynamics community. Computational examples associated with practical configurations have shown the utilization of the algorithm.

• Journal of Biomedical Engineering Research
• /
• v.28 no.3
• /
• pp.377-386
• /
• 2007

Mechanical loading to bone cells using simple sine wave or constant wave fluid flow has been widely used for in vitro experiments. Human gait is characterized by a complex loading to bones of lower extremities which results from a series of events consisting of heel strike, foot flat and push-off during the stance phase of the gait cycle. Telemetric force analyses have shown that human femora are subject to multiphasic loading. Therefore, it would be ideal if the physiologic loading conditions during human walking can be used for in vitro mechanotransduction studies. Here, for a mechanotransduction study, we develop it fluid flow system (FFS) in order to simulate human physiologic mechanicalloading on bone cells. The development methods of the FFS including the COR (Center for Orthopedic Research), monitor program are presented. The FFS could generate various multiphasic loading conditions of human gaits with output flow. Wall shear distribution was very uniform, with 81 % of the effective loading area of the culture on a glass slide. Our results demonstrated that the FFS, provide a new translational approach for unveiling molecular mechanotransduction pathways in bone cells.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.10 no.1
• /
• pp.50-55
• /
• 1986

Characteristics of the flow and heat transfer in a channel of 37 rods are investigated numerically. The flow is taken to be a fully developed incompressible laminar flow and it has an uniform temperature profile at the inlet and flows down through the channel of constant wall temperature. A boundary-fitted coordinate system is used for the complex geometry. Calculation is initiated by calculating the developed flow profile and then proceeds to temperature development. Through the calculation the details of the flow and temperature distribution characteristics are found, and discussion is made on the mechanism of the transport phenomena in the complex geometry in terms of wall shear stress distribution, non-dimensionalized velocity, friction factor, Nusselt number distribution, Reynolds number, and porosity. Also the effects of the eccentricity in rod configuration are analyzed and its importance is emphasized.