• Journal of applied mathematics & informatics
• /
• v.35 no.3_4
• /
• pp.241-260
• /
• 2017

An analysis is made to study the solute transport in a Casson fluid flow through an annulus in presence of oscillatory flow field and determine how this flow influence the solute dispersion along the annular region. Axial dispersion coefficient and the mean concentration expressions are calculated using the generalized dispersion model. Dispersion coefficient in oscillatory flow is found to be a function of frequency parameter, Schmidt number, and the pressure fluctuation component besides its dependency on yield stress of the fluid, annular gap and time in the case of steady flow. Due to the oscillatory nature of the flow, the dispersion coefficient changes cyclically and the amplitude and magnitude of the dispersion increases initially with time and reaches a non - transient state after a certain critical time. This critical value varies with frequency parameter and independent of the other parameters. It is found that the presence of inner cylinder and increase in the size of the inner cylinder inhibits the dispersion process. This model may be used in understanding the dispersion phenomenon in cardiovascular flows and in particular in catheterized arteries.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.2
• /
• pp.456-466
• /
• 1994

Experimental investigation was made to study the mechanism of fluid and thermal oscillation phenomena of surface-tension driven flow in a cylindrical liquid column heated from above which is the low-gravity floating zone simulated on earth. Hexadecane, octadecane, silicon oil (10cs), FC-40 and water are used as the test liquids. The onset of the oscillatory thermocapillary convection appears when Marangoni number exceeds its criteria value and is found to be due to the coupling among velocity and temperature field with the free surface deformation. The frequency of temperature oscillation decreases with increasing aspect ratio for a given diameter and Marangoni number and the oscillation level increases with Marangoni number. The flow pattern in the liquid column appears either as symmetric or asymmetric 3-D flow due to the oscillatory flow in the azimuthal direction. The free surface deformation also occurs either as symmetric or asymmetric mode and its frequency is consistent with those of flow and temperature oscillations. The amplitude of surface deformation also increases with Marangoni number.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.3
• /
• pp.433-441
• /
• 2001

Natural convection of a fluid with intermediate Prand시 number of Pr=0.2 in a horizontal annulus is considered, and the bifurcation phenomena and chaotic flows are numerically investigated. The unsteady two-dimensional streamfunction-vorticity equation is solved with finite difference method. The steady downward flow with two counter-rotating eddies bifurcates to a simple periodic flow with a fundamental frequency. And afterwards, second Hopf bifurcation occurs, and a quasi-periodic flow with two incommensurable frequencies appears. However, a new time-periodic flow is established after experiencing quasi-periodic states. As Rayleigh number is increased further, the chaotic flow regime is reached after a sequence of successive Hopf bifurcation to quasi-periodic and chaotic flow regimes. A scenario similar to the Ruelle-Takens-Newhouse scenario of the onset of chaos is observed.

• Journal of computational fluids engineering
• /
• v.16 no.3
• /
• pp.29-36
• /
• 2011

The natural convection in a horizontal enclosure heated from the bottom wall, cooled at the top wall, and having a square adiabatic body at its centered area was studied. Three different Prandtl numbers (0.01, 0.7 and 7) were considered for an effect of the Prandtl number on natural convection. A two-dimensional solution for unsteady natural convection was obtained, using Chebyshev spectral methodology for different Rayleigh numbers varying over the range of $10^4$ to $10^6$. It had been experimentally and numerically reported [1,2] that the heat transfer mode becomes oscillatory when Pr is out of a specific Pr band beyond the critical Ra. In this study, we reproduced this phenomenon numerically. The variation of time- and surface-averaged Nusselt numbers on the hot and cold walls for different Rayleigh numbers and Prandtl numbers was presented to show the overall heat transfer characteristics in the system. And also, the isotherms and streamline distributions were presented in detail to compare the physics related to their thermal behavior.

• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• pp.311-315
• /
• 2000

A numerical investigation is made of time-dependent buoyant convection in a square of a non-Boussinesq fluid. The density-temperature$({\rho}-T)$ relation is modeled by a quadratic function, with the maximum density ${\rho}_M$ at temperature $T_M$. The horizontal walls of the square are insulated, and a pulsating temperature $T_H=T_M+{\Delta}T'\;sin({\omega}{\tau})$ is imposed on the hot vertical sidewall. The temperature at the cold wall $T_c$ is constant. Extensive numerical solutions to the governing Navier-Stokes equations are portrayed. Resonance is identified by monitoring the amplitude of the mid-plane Nusselt number, $A(Nu^*)$. The primary resonance frequency is found by matching ${\omega}$ to the nondimensional basic mode $N_1$ of internal gravity oscillations. Due to the quadratic$({\rho}-T)$ relationship, the effective pulsation frequency for density, $2{\omega}$, is meaningful, which brings forth the secondary resonance frequency, i.e., $2{\omega}=N_1$

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.3
• /
• pp.381-388
• /
• 2001

Presented are heat data which describe the effect of interaction between bulk flow pulsations and a vortex embedded in a turbulent boundary layer. The pulsation frequencies are 3 Hz, 15 Hz and 30 Hz. A half delta wing with the same height as the boundary layer thickness is used to generate the vortex flow. The convection heat transfer coefficients on a constant heat-flux surface are measured by embedded 77 T-type thermocouples. Spanwise profiles of convection heat transfer coefficients show that upwash region of vortex flow is influenced by bulk flow pulsations. The local heat transfer coefficient increases approximately by 7 percent. The increase in the local change of convection heat transfer coefficient is attributed to the spanwise oscillatory motion of vortex flow especially at the low Strouhal number and to the periodic change of vortex size.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.9
• /
• pp.2953-2964
• /
• 1996

Heat Transfer with periodic fluctuation of fluid temperature caused by oscillatory flow or compression expansion can be out of phase with balk fluid-wall temperature difference. Newton's law of convection is inadequate to describe this phenomenon. In order to solve this problem the concept of the complex Nusselt number has been introduced by severla researchers. The complex Nusselt number expresses out of phase excellently while the first harmonic is dominant in the variations of both fluid-wall temperature difference and heat flux. However, in the case of oscillatory flow with non-linear wall temperature distribution, the complex Nusselt number is not appropriate to predict the heat transfer phenomena since the higher order harmonic components appear in periodic temperature variation. Analytic solutions to the heat transfer with an sinusoidal well temperature distribution were obtained to investagate the effect of non-linear wall temperature distribution. A new formula considering the thermal boundary layer was suggested based on the solutions. A comparison was also made with the complex Nusselt number. It was verified that the new formula describes well the heat transfer of oscillating flow even if the first harmonic component is not dominant in the fluid-wall temperature difference.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.7 no.3
• /
• pp.528-537
• /
• 1995

Natural convective flow and heat transfer characteristics in a partially opened enclosure fitted with a horizontal divider are investigated numerically. The enclosure is composed of a lower hot and a upper cold horizontal walls and adiabatic vertical walls. A divider is attached perpendicularly to the vertical insulated wall. The governing equations are solved by using the finite element method with Galerkin method. The computations have been carried out by varying the length of divider, the opening size, and the Rayleigh number based on the temperature difference between two horizontal walls and the enclosure height for air(Pr=0.71). As result, when the opening size is fixed, the intensity of the secondary flow is weaken as the length of divider increases. The maximum heat transfer rate over the upper cold wall occurs at a position bounded on the opening. However, when the length of divider is increased considerably, its maximum occurs at the right wall. The stability and frequency of oscillation are affected by the Rayleigh number and length of divider. The Nusselt number is increased with the increase of the opening size and the increase of Rayleigh number.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.12 no.2
• /
• pp.150-160
• /
• 2000

Turbulent natural convective flow and heat transfer in a square enclosure with horizontal partition are investigated numerically. The enclosure is composed of a lower hot and a upper cold horizontal walls and adiabatic vertical walls. Partitions carried with the upward, downward, and both control plates are attached perpendicularly to the one of the vertical insulated walls, respectively. The low Reynolds number $k-\varepsilon$ model is adopted to calculate the turbulent thermal convection. The governing equations are solved by using the finite element method with Galerkin method. The computations have been carried out by varying the length of partition, the position of control plates, and the Rayleigh number based on the temperature difference between two horizontal walls and the enclosure height for water(Pr=4.95). When the control plates are attached at the edge of partition, the stability of oscillating flow grows wrose with the increase of Rayleigh number and the partition length. The heat transfer rate has been reducer than that of no control plate due to the restraint of control plates with the increase of Rayleigh number.

• Journal of applied mathematics & informatics
• /
• v.25 no.1_2
• /
• pp.51-66
• /
• 2007

A layer of compressible, rotating, elastica-viscous fluid heated & soluted from below is considered in the presence of vertical magnetic field to include the effect of Hall currents. Dispersion relation governing the effect of viscoelasticity, salinity gradient, rotation, magnetic field and Hall currents is derived. For the case of stationary convection, the Rivlin-Erickson fluid behaves like an ordinary Newtonian fluid. The compressibility, stable solute gradient, rotation and magnetic field postpone the onset of thermosolutal instability whereas Hall currents are found to hasten the onset of thermosolutal instability in the absence of rotation. In the presence of rotation, Hall currents postpone/hasten the onset of instability depending upon the value of wavenumbers. Again, the dispersion relation is analyzed numerically & the results depicted graphically. The stable solute gradient and magnetic field (and corresponding Hall currents) introduce oscillatory modes in the system which were non-existent in their absence. The case of overstability is discussed & sufficient conditions for non-existence of overstability are derived.