Thermal instability in natural convection flow over a boundary layer subject to external electrical and magnetic fields

This paper presents a numerical study of external electrical and magnetic effects on the formation of longitudinal vortices in natural convection flow over a heated horizontal plate. The criterion on the position marking on the onset of longitudinal vortices is defined in the present paper. The onset position characterized by the Grashof number depends on the Prandtl number, the wave number, the electric field parameter, and the magnetic field parameter. The flow is found more stable as the value of the magnetic field parameter increases. Moreover, the stabilizing effect is also found on the flow when the positive electric field parameter Ec is applied. The results of the present numerical prediction show reasonable agreement with the experimental data with zero magnetic field and electric field parameters in literature.     

Effect of rotation on the formation of longitudinal vortices in mixed convection flow over a flat plate

This paper presents a numerical study of the effect of rotation on the formation of longitudinal vortices in mixed convection flow over a flat plate. The criterion on the position of marking the onset of longitudinal vortices is defined in this paper. The onset position characterized by the Goertler number G _δ depends on the Grashof number, the rotation number Ro, the Prandtl number Pr and the wave number. The results show that negative rotation stabilizes the boundary layer flow on the surface. On the contrary, positive rotation destabilizes the flow. The numerical data are compared with the experimental results.     

Effect of magnetic field on thermal instability in mixed convection flow over a rotating boundary layer

Effect of magnetic field on the formation of longitudinal vortices in mixed convection flow over a rotating heated flat plate is presented. The onset position is characterized by the Grashof number, the rotational number, the Prandtl number, the Eckert number, the magnetic field parameter, and the wave number. Negative rotation (clockwise) and external magnetic field stabilize the boundary layer flow. On the contrary, positive rotation (anti-clockwise), the Eckert number, and the Prandtl number destabilize the flow. The numerical data show agreement with the experimental data with the case of zero Hartmann number in the literature.     

Thermal instability of a boundary layer subjected to external electric and magnetic fields

This paper presents a numerical study of the electrical and magnetic fields on thermal instability in a boundary layer. The criterion on the position marking on the onset of longitudinal vortices is defined in the present paper. The results show that the onset position characterized by the Grashof number depends on the Prandtl number, the Reynolds number, the wave number, the electric field parameter, and the Hartmann number. The flow becomes more stable as the magnetic field increases. However, the destabilizing effect is found on the flow when the negative electric field parameter is applied. The results of the present numerical prediction show reasonable agreement with the experimental data in the case of zero Hartmann number and zero electric field parameter in the literature. Copyright © 2010 John Wiley & Sons, Ltd.     

Effect of ion‐slip current on the thermal instability of mixed connection in a boundary layer flow

Effect of ion‐slip current on the thermal instability in a boundary layer is studied. The criterion on the position marking the onset of longitudinal vortices is defined in the present paper. The results show that the onset position characterized by the Grashof number depends on the Prandtl number, the Reynolds number, the wave number, the Hall parameter, the ion‐slip parameter, and the Hartmann number. The flow becomes more stable as the magnetic field increases. However, the destabilizing effects are found on the flow when the Hall and ion‐slip currents are presented. The results of the present numerical prediction show reasonable agreement with the experimental data in the case of zero Hartmann number, ion‐slip parameter, and Hall parameter in the open literature. Copyright © 2011 John Wiley & Sons, Ltd.     

Numerical prediction of the formation of Goertler vortices on a concave surface with suction and blowing

This paper presents a numerical prediction of the formation of Goertler vortices on a concave surface with suction and blowing. Suction stabilizes the boundary layer flow on the surface, whereas blowing destabilizes the flow. The criterion on the position marking the onset of Goertler vortices is defined in the present paper. For facilitating the numerical study, the computation is carried out in the transformed x–η plane. The results show that the onset position characterized by the Goertler number depends on the local suction/blowing parameter, the Prandtl number and the wavenumber. The value of the critical Goertler number increases with the increase in suction, while the value of the Goertler number decreases with the increase in blowing. Both the experimental and the numerical data can be correlated by G_θ*=10.2(a′θ)*^3/2 without suction and blowing and by a simple relation G*_x=(G*_x)_γ=0 e^−γ with suction and blowing. The obtained critical Goertler number and wavenumber are in good agreement with the previous experimental data. Copyright © 1999 John Wiley & Sons, Ltd.     

Effect of rotation on Goertler vortices in the boundary layer flow on a curved surface

This paper presents a numerical prediction of the formation of Goertler vortices on a curved surface with effect of rotation. The criterion of flow visualization marking the onset position of Goertler vortices is employed in the present paper. For facilitating the numerical study, the computation is carried out in the transformed x and ηplane. The results show that the onset position characterized by the Goertler number, depends on the rotation number Ro, the Prandtl number and the wave number. The value of critical Goertler number increases with the increase in negative rotation, while the value of Goertler number decreases with the increase in positive rotation on a concave surface. On the contrary, the value of critical Goertler number decreases with the increase in negative rotation on a convex surface. The obtained critical Goertler number and wave number are compared with the previous theoretical and experimental data. Copyright © 2002 John Wiley & Sons, Ltd.     

Motion of a sphere in an electrically conducting rotating fluid

The combined effect of rotation and magnetic field is investigated for the axisymmetric flow due to the motion of a sphere in an inviscid, incompressible electrically conducting fluid having uniform rotation far upstream. The steady-state linearized equations contain a single parameter α=1/2βR _m, β being the magnetic pressure number and R _m the magnetic Reynolds number. The complete solution for the flow field and magnetic field is obtained and the distribution of vorticity and current density is found. The induced vorticity is O(α⁴) and the current density is O(R _m) on the sphere.     

Time resolved analysis of steady and oscillating flow in the upper human airways

In this experimental study a thorough analysis of the steady and unsteady flow field in a realistic transparent silicone lung model of the first bifurcation of the upper human airways will be presented. To determine the temporal evolution of the flow time-resolved particle-image velocimetry recordings were performed for a Womersley number range 3.3 ≤ α ≤ 5.8 and Reynolds numbers of Re _D = 1,050, 1,400, and 2,100. The results evidence a highly three-dimensional and asymmetric character of the velocity field in the upper human airways, in which the influence of the asymmetric geometry of the realistic lung model plays a significant role for the development of the flow field in the respiratory system. At steady inspiration, the flow shows independent of the Reynolds number a large zone with embedded counter-rotating vortices in the left bronchia ensuring a continuous streamwise transport into the lung. At unsteady flow the critical Reynolds number, which describes the onset of vortices in the first bifurcation, is increased at higher Womersley number and decreased at higher Reynolds number. At expiration the unsteady and steady flows are almost alike.     

The effects of local thermal nonequilibrium and MFD viscosity on the onset of Brinkman ferroconvection

The simultaneous effect of local thermal nonequilibrium (LTNE) and magnetic field dependent (MFD) viscosity on thermal convective instability in a horizontal ferrofluid saturated Brinkman porous layer in the presence of a uniform vertical magnetic field is studied analytically. The results indicate that the onset of Brinkman ferroconvection is delayed with increasing MFD viscosity parameter but the critical wave number is found to be independent of this parameter. When compared to the simultaneous presence of buoyancy and magnetic forces, it is observed that the onset of Brinkman ferroconvection is delayed more when the magnetic forces alone are present. Asymptotic solutions for both small and large values of scaled inter-phase heat transfer coefficient H _t are compared with those computed numerically and good agreement is found between them. Besides, the influence of magnetic and LTNE parameters on the stability characteristics of the system is also discussed. The available results in the literature are recovered as particular cases from the present study.     

On the formation of Taylor–G?rtler vortices in RMF-driven spin-up flows

This paper is concerned with a liquid metal flow driven by a rotating magnetic field inside a stationary cylinder. We consider especially the secondary meridional flow during the time when the fluid spins up from rest. The developing flow is investigated experimentally and by direct numerical simulations. The vertical profiles of the axial velocity are measured by means of the ultrasound Doppler velocimetry. Evolving instabilities in the form of Taylor–G?rtler vortices have been observed just above the instability threshold (Ta ≥ 1.5· Ta ^cr). The rotational symmetry may survive over a distinct time even if a first Taylor–G?rtler vortex pair has been formed as closed rings along the cylinder perimeter. The transition to a three-dimensional flow in the side layers results from the advection or a precession and splitting of the Taylor–G?rtler vortex rings. The predictable behaviour of the Taylor–G?rtler vortices disappears with increasing magnetic field strength. The numerical simulations agree very well with the flow measurements.     

Hall effects on MHD flow in a rotating system with heat transfer characteristics

Closed-form solutions are derived for the steady magnetohydrodynamic (MHD) viscous flow in a parallel plate channel system with perfectly conducting walls in a rotating frame of reference, in the presence of Hall currents, heat transfer and a transverse uniform magnetic field. A mathematical analysis is described to evaluate the velocity, induced magnetic field and mass flow rate distributions, for a wide range of the governing parameters. Asymptotic behavior of the solution is analyzed for large M ² (Hartmann number squared) and K ² (rotation parameter). The heat transfer aspect is considered also with Joule and viscous heating effects present. Boundary layers arise close to the channel walls for large K ², i.e. strong rotation of the channel. For slowly rotating systems (small K ²), Hall current parameter (m) reduces primary mass flow rate (Q _x /R ρ v). Heat transfer rate at the upper plate (d θ/d η) _η=1 decreases, while at the lower plate (d θ/d η) _η=−1 increases, with increase in either K ² or m. For constant values of the rotation parameter, K ², heat transfer rate at both plates exhibits an oscillatory pattern with an increase in Hall current parameter, m. The response of the primary and secondary velocity components and also the primary and secondary induced magnetic field components to the control parameters is also studied graphically. Applications of the study arise in rotating MHD induction machine energy generators, planetary and solar plasma fluid dynamics systems, magnetic field control of materials processing systems, hybrid magnetic propulsion systems for space travel etc.     

On Temporal Stability Analysis for Hydromagnetic Flow in a Channel Filled with a Saturated Porous Medium

In this paper, a linear stability analysis is presented to trace the time evolution of an infinitesimal, two-dimensional disturbance imposed on the base flow of an electrically conducting fluid in a channel filled with a saturated porous medium under the influence of a transversely imposed magnetic field. An eigenvalue problem is obtained and solved numerically using the Chebyshev collocation spectral method. The critical Reynolds number Re _c, the critical wave number α _c and the critical wave speed c _c are obtained for a wide range of the porous medium shape factor parameter S and Hartmann number H. It is found that an increase in the magnetic field intensity and a decrease in porous medium permeability have a stabilizing effect on the fluid flow.     

Maximum density effects on convective instability of horizontal plane poiseuille flows in the thermal entrance region

A linear stability analysis is used to study the conditions marking the onset of secondary flow in the form of longitudinal vortices for plane Poiseuille flow of water in the thermal entrance region of a horizontal parallel-plate channel by a numerical method. The water temperature range under consideration is 0∼30°C and the maximum density effect at 4°C is of primary interest. The basic flow solution for temperature includes axial heat conduction effect and the entrance temperature is taken to be uniform at far upstream location jackie=−∞ to allow for the upstream heat penetration through thermal entrance jackie=0. Numerical results for critical Rayleigh number are obtained for Peclet numbers 1, 10, 50 and thermal condition parameters (λ ₁, λ ₂) in the range of −2.0≤λ ₁≤−0.5 and −1.0≤λ ₂≤1.4. The analysis is motivated by a desire to determine the free convection effect on freezing or thawing in channel flow of water.     

Understanding coherent vortices through computational fluid dynamics

A new definition of coherent vortices in turbulence is proposed, where the vorticity equation reduces to a cyclostrophic balance. Afterward, we describe five fundamental vortex interactions, the sheet, the spiral, the pairing, the even longitudinal, and the odd longitudinal modes. Numerous examples of these interactions are provided from direct numerical or large-eddy simulations. The resulting vortices are responsible for the internal intermittent character of turbulence, with highly nongaussian tails for the probability density functions of vorticity, passive scalar, and low pressure. In a mixing layer, the combination of the odd longitudinal and the pairing modes (helical pairing) is inhibited by compressibility, above a convective Mach number of 0.7. When turbulence is submitted to a solid-body rotation, anticyclonic vortices of local Rossby number of the order of 1 transform into intense perpendicular Görtler-type alternate longitudinal vortices.The support of CCVR, CEA, CNRS, Dassault/CNES, DRET, LHF, and Région Rhône-Alpes is acknowledged. This paper is the text of an invited lecture given at the IUTAM Symposium on Eddy Structure Identification in Free Turbulent Shear Flows, Poitiers, 12–14 October 1992.     

Effect of Thermal Non-Equilibrium on Convective Instability in a Ferromagnetic Fluid-Saturated Porous Medium

The effect of local thermal non-equilibrium (LTNE) on the onset of thermomagnetic convection in a ferromagnetic fluid-saturated horizontal porous layer in the presence of a uniform vertical magnetic field is investigated. A modified Forchheimer-extended Darcy equation is employed to describe the flow in the porous medium, and a two-field model is used for temperature representing the solid and fluid phases separately. It is found that both the critical Darcy–Rayleigh number and the corresponding wave number are modified by the LTNE effects. Asymptotic solutions for both small and large values of scaled interphase heat transfer coefficient H _t are presented and compared with those computed numerically. An excellent agreement is obtained between the asymptotic and the numerical results. Besides, the influence of magnetic parameters on the instability of the system is also discussed. The available results in the literature are recovered as particular cases from the present study.     

Effect of magnetic Reynolds number on the two‐dimensional hydromagnetic flow around a cylinder

Numerical experiments have been conducted to study the effect of magnetic Reynolds number on the steady, two‐dimensional, viscous, incompressible and electrically conducting flow around a circular cylinder. Besides usual Reynolds number Re, the flow is governed by the magnetic Reynolds number R_m and Alfvén number β. The flow and magnetic field are uniform and parallel at large distances from the cylinder. The pressure Poisson equation is solved to find the pressure fields in the entire flow region. The effects of the magnetic field and electrical conductivity on the recirculation bubble, drag coefficient, standing vortex and pressure are presented and discussed. For low interaction parameter (N<1), the suppression of the flow‐separation is nearly independent of the conductivity of the fluid, whereas for large interaction parameters, the conductivity of the fluid strongly influences the control of flow‐separation. Copyright © 2008 John Wiley & Sons, Ltd.     

PIV–LES analysis of channel flow rotating about the streamwise axis

The flow field of a channel rotating about the streamwise axis is analyzed experimentally and numerically. The current investigations were carried out at a bulk velocity based Reynolds number of Re_m = 2850 and a friction velocity based Reynolds number of Re_τ = 180, respectively. Particle-image velocimetry (PIV) measurements are compared with large-eddy simulation data to show earlier direct numerical simulation findings to generate too large a reverse flow region in the center region of the spanwise flow. The development of the mean spanwise velocity distribution and the influence of the rotation on the turbulent properties, i.e., the Reynolds stresses and the two-point correlations of the flow, are confirmed in both investigations. The rotation primarily influences those components of the Reynolds shear stresses, which contain the spanwise velocity component. The size of the correlation areas and thus the length scales of the flow generally grow in all three coordinate directions leading to longer structures. Furthermore, experimental results of the same channel flow at a significantly lower bulk Reynolds number of Re_{m, l} = 665, i.e., a laminar flow in a non-rotating channel, are introduced. The experiments show the low Reynolds number flow to become turbulent under rotation and to develop the same characteristics as the high Reynolds number flow.     

Three-dimensional flow of Newtonian and Boger fluids in square–square contractions

The flow of a Newtonian fluid and a Boger fluid through sudden square–square contractions was investigated experimentally aiming to characterize the flow and provide quantitative data for benchmarking in a complex three-dimensional flow. Visualizations of the flow patterns were undertaken using streak-line photography, detailed velocity field measurements were conducted using particle image velocimetry (PIV) and pressure drop measurements were performed in various geometries with different contraction ratios. For the Newtonian fluid, the experimental results are compared with numerical simulations performed using a finite volume method, and excellent agreement is found for the range of Reynolds number tested (Re₂ ≤ 23). For the viscoelastic case, recirculations are still present upstream of the contraction but we also observe other complex flow patterns that are dependent on contraction ratio (CR) and Deborah number (De₂) for the range of conditions studied: CR = 2.4, 4, 8, 12 and De₂ ≤ 150. For low contraction ratios strong divergent flow is observed upstream of the contraction, whereas for high contraction ratios there is no upstream divergent flow, except in the vicinity of the re-entrant corner where a localized atypical divergent flow is observed. For all contraction ratios studied, at sufficiently high Deborah numbers, strong elastic vortex enhancement upstream of the contraction is observed, which leads to the onset of a periodic complex flow at higher flow rates. The vortices observed under steady flow are not closed, and fluid elasticity was found to modify the flow direction within the recirculations as compared to that found for Newtonian fluids. The entry pressure drop, quantified using a Couette correction, was found to increase with the Deborah number for the higher contraction ratios.     

流向磁场作用下圆柱绕流的直接数值模拟

绕流是托卡马克装置中液态包层内常见的流动形态,对流场与热量分布有着重要的影响.本文通过直接数值模拟(DNS),研究了不同磁场强度下$Re=3900$的圆柱绕流,分析了磁场强度对于湍流尾迹的影响.无磁场情况下,直接数值模拟的结果与前人的实验及模拟结果吻合很好.圆柱下游的尾迹中,随着流向距离的增大, 流向速度剖面逐渐从U型进化呈V型, 并慢慢趋于平缓,这表明尾迹中的流动结构受圆柱影响逐渐减小.圆柱后方两侧的剪切层中,由于Kelvin-Helmholtz不稳定性的影响,可以清晰地看到小尺度剪切层涡的脱落.通过对无磁场的计算结果施加流向磁场,本文计算了哈特曼数($Ha$)分别为20, 40和80的工况,以研究磁场效应对于湍流的影响.结果表明磁场较弱时,流动依然呈三维湍流状态.随着磁场增强, 近圆柱尾流区受磁场抑制明显,回流区被拉长,剪切层失稳位置向下游转移.圆柱后方的涡结构由于受到竖直方向洛伦兹力的挤压作用,随着哈特曼数的增加尾迹区域逐渐变窄.相比于无磁场情况的涡结构,由于磁场的耗散作用,相应的涡结构尺度变小.该研究不仅扩展了现有磁场下湍流运动的参数范围,对于液态包层的设计及安全运行同样具有重要的理论指导意义和工程应用价值.