• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.21 no.2
• /
• pp.89-107
• /
• 2017

In this paper, we have proposed a modified Marker-And-Cell (MAC) method to investigate the problem of an unsteady 2-D incompressible flow with heat and mass transfer at low, moderate, and high Reynolds numbers with no-slip and slip boundary conditions. We have used this method to solve the governing equations along with the boundary conditions and thereby to compute the flow variables, viz. u-velocity, v-velocity, P, T, and C. We have used the staggered grid approach of this method to discretize the governing equations of the problem. A modified MAC algorithm was proposed and used to compute the numerical solutions of the flow variables for Reynolds numbers Re = 10, 500, and 50000 in consonance with low, moderate, and high Reynolds numbers. We have also used appropriate Prandtl (Pr) and Schmidt (Sc) numbers in consistence with relevancy of the physical problem considered. We have executed this modified MAC algorithm with the aid of a computer program developed and run in C compiler. We have also computed numerical solutions of local Nusselt (Nu) and Sherwood (Sh) numbers along the horizontal line through the geometric center at low, moderate, and high Reynolds numbers for fixed Pr = 6.62 and Sc = 340 for two grid systems at time t = 0.0001s. Our numerical solutions for u and v velocities along the vertical and horizontal line through the geometric center of the square cavity for Re = 100 has been compared with benchmark solutions available in the literature and it has been found that they are in good agreement. The present numerical results indicate that, as we move along the horizontal line through the geometric center of the domain, we observed that, the heat and mass transfer decreases up to the geometric center. It, then, increases symmetrically.

• Journal of computational fluids engineering
• /
• v.18 no.3
• /
• pp.34-41
• /
• 2013

We develop a novel immersed boundary (IB) method based on implicit direct forcing scheme for incompressible flows. The proposed IB method is based on an iterative procedure for calculating the direct forcing coupled with the momentum equations in order to satisfy no-slip boundary conditions on IB surfaces. We perform simulations of two-dimensional flows over a circular cylinder for low and moderate Reynolds numbers. The present method shows that the errors for estimated velocities on IB surfaces are significantly reduced even for low Reynolds number with a fairly large time step while the previous methods based on direct forcing failed to provide no-slip boundary conditions on IB surfaces.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.9
• /
• pp.1263-1271
• /
• 2001

An experimental study is carried out for an aluminum foam heat sink attached to an isolated heat source to evaluate high potential of aluminum foam as a heat sink with impinging jets. The effects of the pore density and the height of the aluminum foam heat sink, the jet Reynolds number, and the nozzle diameter are delineated in comparison with a conventional pin type heat sink. It is found that the aluminum foam with small pores is inefficient for the heat transfer enhancement due to the large flow friction at the given porosity. In the parameter ranges of the present study, the change in the nozzle diameter shows no significant effects on the surface temperature of the aluminum foam heat sink at a given Reynolds number. The heat transfer enhancement is strongly dependent on the jet Reynolds number and shows a maximum value at a moderate Reynolds number.

• Wind and Structures
• /
• v.34 no.1
• /
• pp.15-27
• /
• 2022

This work comes within the framework of the "Benchmark on the Aerodynamics of a Rectangular Cylinder" that investigates a rectangular cylinder of length-to-depth ratio equal to 5. The present study reports and discusses velocity fields acquired using planar Particle Image Velocitmetry for several angles of attack and Reynolds numbers. In particular, for a cylinder depth-based Reynolds number of 2 × 10⁴ and zero incidence angle, the flow features along the lateral (parallel to the freestream) upper and lower surfaces of the cylinder are reported. Using first and second order statistics of the velocity field, the main flow features are discussed, especially the size and location of the time-averaged flow structures and the distribution of the Reynolds stresses. The variation of the flow features with the incidence is also studied considering angles of attack up to 6°. It is shown that the time-averaged flow is fully detached for incidence higher than 2°. For an angle of attack of 0°, the effects of the Reynolds number varying between 5 × 10³ and 2 × 10⁴ are investigated looking at flow statistics. It is shown that the time-averaged location of the reattachment point and the shape and position of the time-averaged main vortex are mostly constant with the Reynolds number. However, the size of the inner region located below the time-averaged shear layer and just downstream the leading edge corner appears to be strongly dependent on the Reynolds number.

• Advances in nano research
• /
• v.15 no.2
• /
• pp.91-104
• /
• 2023

Development of thin carbon nanotubes (CNTs) nanoliquid film over the rough surface of a horizontal rotating disk is investigated by considering symmetric roughness either along the azimuthal or radial directions. The disk surface is either heated or cooled axisymmetrically from below. The effects of single-walled carbon nanotubes (SWCNTs), multi-walled carbon nanotubes (MWCNTs) are analyzed on the film thinning process with different types of base liquids. Closed form solutions for velocity and temperature field are obtained for small values of Reynolds number whereas the numerical solution is derived for moderate values of Reynolds number. It is found that fluid retention / depletion takes place when the roughness is symmetric along the azimuthal / radial directions. It is also seen that the film thinning rate enhances for MWCNTs compare to SWCNTs. Further it is found that two different heat transfer regions exits within the flow domain depending on the fact that heat is transferred from disk to liquid film and vice-versa.

• Journal of Advanced Marine Engineering and Technology
• /
• v.36 no.8
• /
• pp.1091-1096
• /
• 2012

The aerodynamic characteristics of circular cylinder in a channel are studied to make clear the flow feature by solving the Navier-Stokes equation based on the finite volume method with unstructured grids. Reviews are made on with the vorticity, velocity, dynamic pressure, residual and drag, where the Reynolds numbers are 50 and 100. The flows for $Re{\succeq}50$ shows the vortex shedding in the wake, and the result is the same as the case of moving cylinder. The ground effect of flat bottom results in the growth of vortex, being generated in the upper side of the cylinder and elongated in the rear. As the cylinder approaches to wall, for example 0.6, the cylinder plays as a role of blockage to obstruct the flow between the cylinder and wall. The drag coefficients are compared with others' results to confirm the validity of the present numerical simulation.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.10
• /
• pp.2738-2750
• /
• 1994

The Effects of radiative heat transfer on turbulent flow in a partitioned enclosure is studied numerically. The enclosure is partially divided by a thin, poorly conducting vertical divider projecting from the ceiling of the enclosure. The low Reynolds number $k-{\epsilon}$ model is adopted to calculate the turbulent flow field. The solutions to the radiative transfer equations are obtained by the discrete ordinates method(DOM). This method is based on control volume method and is compatible with the SIMPLER algorithm used to solve the momentum and energy equations. The effects of optical thickness and Planck number on the flow, temperature fields and heat transfer rates are investigated for a moderate Rayleigh number($=10^9$). The changes in buoyant flow fields and temperature distributions due to the variation of baffle length are also analyzed. From the predictions, radiant heat exchange between the baffle and the sidewalls strongly influences the temperature distribution in the baffle and its vicinity and total heat transfer increases as the optical thickness and the baffle length decrease. It is possible to neglect the radiative heat transfer effect when Planck number is over one.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.30 no.5 s.248
• /
• pp.438-445
• /
• 2006

We performed the simulation of the unsteady three dimensional flow over a square cylinder in a wind tunnel in moderate Reynolds number range, $100{\sim}2500$ by using LBM. SGS model was applied for the turbulent flow. Frist of all we compared LBM(Lattice Boltzmann Method) solution of Poiseuille flow applied Farout and bounce back boundary conditions with the analytical and FOAM solutions to verify the applicability of the boundary conditions. For LBM simulation the calculation domain was formed by structured grids and prescribed uniform velocity and density inlet and Farout boundary conditions were imposed on the in-out boundaries. Bounceback and wind tunnel boundary conditions were applied to the cylinder walls and the boundaries of calculation domain respectively. The maximum Strouhal number of the vortex shedding is 0.2025 at Re = 750. and the number maintains the constant value of 0.18 when Re>1000. We also predicted that the critical reynolds number of the turbulent flow is in the range of $250{\sim}500$.

• Wind and Structures
• /
• v.16 no.5
• /
• pp.411-431
• /
• 2013

A computational study of vortex-induced transverse vibrations of a cylinder with low mass-damping is presented. An Arbitrary Lagrangian-Eulerian (ALE) formulation of the Unsteady Reynolds-Averaged Navier-Stokes equations (URANS), along with the Spalart-Allmaras (SA) one-equation turbulence model, are coupled conservatively with rigid body motion equations of the cylinder mounted on elastic supports in order to study the amplitude and frequency response of a freely vibrating cylinder, its flow-induced motion, Vortex Street, near-wake flow structure, and unsteady loading in a moderate range of Reynolds numbers. The time accurate response of the cylinder from rest to its limit cycle is studied to explore the effects of Reynolds number on the start of large displacements, motion amplitude, and frequency. The computational results are compared with published physical experiments and numerical studies. The maximum amplitudes of displacements computed for various Reynolds numbers are smaller than the experimental values; however, the overall agreement of the results is quite satisfactory, and the upper branch of the limit-cycle displacement amplitude vs. reduced velocity response is captured, a feature that was missed by other studies. Vortex shedding modes, lock-in phenomena, frequency response, and phase angles are also in agreement with experiments.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2009.11a
• /
• pp.492-495
• /
• 2009

Supersonic microjets acquire considerable research interest from a fundamental fluid dynamics perspective, in part because the combination of highly compressible flow at low-to-moderate Reynolds number is not very common, and in part due to the complex nature of the flow itself. In addition, microjets have a great variety engineering applications such as micro-propulsion, MEMS (Micro-Electro Mechanical Systems) components, microjet actuators and fine particle deposition and removal. Numerical simulations have been carried out at moderate nozzle pressure ratios and for different nozzle exit diameters to investigate and to understand in-depth of aerodynamic characteristics of supersonic microjets.