倾斜非均匀磁场下导电方管磁流体管流的数值模拟研究

热核聚变反应堆液态金属包层应用中的一个重要问题是液态金属在导电管中流动和强磁场相互作用产生的额外的磁流体动力学压降.这种磁流体动力学压降远远大于普通水力学压降.美国阿贡国家实验室ALEX研究小组,对非均匀磁场下导电管中液态金属磁流体动力学效应进行了实验研究,其实验结果成为液态金属包层数值验证的标准模型之一.液态金属包层在应用中会受到不同方向的磁场作用,本文以ALEX的非均匀磁场下导电方管中液态金属管流实验中的一组参数为基础,保持哈特曼数、雷诺数和壁面电导率不变,采用三维直接数值模拟的方法,研究了外加磁场与侧壁之间的倾角对导电方管内液态金属流动的速度、电流和压降分布的影响.研究结果表明:沿流向相同横截面上的速度、电流以及压力分布均随磁场的倾斜而同向旋转.倾斜磁场均匀段,横截面上的高速区位于平行磁场方向的哈特曼层和平行层交叉位置,压力梯度随磁场倾角的增大先增大后减小.倾斜磁场递减段,在三维磁流体动力学效应作用下,横截面上的高速射流位置向垂直磁场方向偏移.磁场递减段的三维磁流体动力学压降随磁场倾角的增大而增大.随磁场倾斜,截面上的射流峰值逐渐减小,二次流增强,引发层流向湍流的转捩.      

Boundary element method solution of MHD flow in a rectangular duct

The magnetohydrodynamic (MHD) flow of an incompressible, viscous, electrically conducting fluid in a rectangular duct with an external magnetic field applied transverse to the flow has been investigated. The walls parallel to the applied magnetic field are conducting while the other two walls which are perpendicular to the field are insulators. The boundary element method (BEM) with constant elements has been used to cast the problem into the form of an integral equation over the boundary and to obtain a system of algebraic equations for the boundary unknown values only. The solution of this integral equation presents no problem as encountered in the solution of the singular integral equations for interior methods. Computations have been carried out for several values of the Hartmann number (1 ? M ? 10). It is found that as M increases, boundary layers are formed close to the insulated boundaries for both the velocity and the induced magnetic field and in the central part their behaviours are uniform. Selected graphs are given showing the behaviours of the velocity and the induced magnetic field.     

Combined effects of rotation and translation on a MHD flow due to compressible viscous fluid

Summary The modification of an axi-symmetric viscous flow due to a relative rotation of a disk or fluid by a translation of the boundary are studied. The fluid is taken to be compressible and electrically conducting. The equations governing the motion are solved iteratively through a central-difference scheme. The effect of an axial magnetic field and disk temperature on the flow and heat transfer are included in the present analysis. The translation of the disk or fluid generates a velocity field at each plane parallel to the disk (secondary flow). The cartesian components of the velocity due to the secondary flow are oscillatory in nature when a rigid body rotation of the free stream along with a translation of the disk is considered. The magnetic field damps out the velocity field, and reduces the thickness of the boundary layer. The cross component of wall shear due to secondary flow acts in a direction opposite to the rotation of the disk or fluid for all cases of the motion. The rise in disk temperature produces an increment in the magnitude of the wall shear associated with the secondary flow.     

The structure of parallel layers in steady two-dimensional magnetohydrodynamic flows in sudden duct expansions and contractions

The structure of non-linear, steady, two-dimensional parallel layers at high values of the Hartmann number, Ha, the Reynolds number, Re, and the interaction parameter, N, for duct expansions and contractions has been investigated. The magnetic field is transverse to the flow. For the flow regime, in which the electromagnetic force balances the inertial force in the layer, a viscous sublayer at the solid wall parallel to the magnetic field has been obtained, in which the flow is driven by the pressure gradient induced in the outer, inviscid layer.     

MHD forced convection boundary layer flow with a flat plate and porous substrate

The flow and heat transfer for an electrically conducting fluid with a porous substrate and a flat plate under the influence of magnetic field is considered. The magnetic field is assumed to be uniform and also along normal to the surface. The momentum and energy equations are transformed to ordinary differential equations by using suitable similarity transformation and are solved by standard techniques. But the energy equation is solved by considering two boundary layers, one in the porous substrate and the other above the porous substrate. Numerical results are presented through graphs with various values of magnetic parameter for both velocity and thermal boundary layers along with Nusselt number and for various values of Prandtl number and Eckert number in thermal boundary layer.     

Linear approximation for the flow of a conducting gas with a high magnetic reynolds number through a channel

In this paper, we discuss the flow of a nonviscous and non-heat-conducting gas through a channel of variable cross section under the influence of a transverse magnetic field. For high magnetic Reynolds numbers, the flow is shown to consist of a core and current layers at the electrodes and at the fixed channel walls. The distributions of currents and other parameters in the core and in the current layers are found analytically, in a linear approximation. The Joule dissipation in the current layers may be more intense than that in the core. The longitudinal currents and Joule dissipation increase with increasing Hall parameter in the electrode layers. Zhigulev [1] has shown that magnetic boundary layers may form in the flow of a conducting gas when there is a high magnetic Reynolds number (R_m«1). He illustrated this situation by the shielding of a plasma flow from the magnetic fields produced near a plate which is electrically isolated from the plasma and through which a current is flowing. In an incompressible fluid, the layer thickness is proportional to R_m ^?1/2. Morozov and Shubin [2] have offered a linear-approximation treatment of the structure of the electromagnetic near-electrode layers which arise during the flow of a nonviscous plasma with a high R_m and a small “exchange” parameter ξ≈H/R_m, for flow transverse to a magnetic field and near a corrugated wall. They pointed out the possible formation of “dissipationless” near-electrode layers with thicknesses on the order of the Debye or electron Larmor radii, and a “dissipative” layer whose thickness increases along the length of the electrodes and is proportional to (R_mc_B ²/c_T ²)^?1/2, where c_B and c_T are the magnetic and thermal sound velocities. Morozov and Shubin studied the properties of dissipationless and dissipative electromagnetic layers at segmented accelerator electrodes through which a current is passing, for an arbitrary “exchange” parameter, in [2] and [3], respectively. The exchange parameter ξ was found in [4]. Such layers should also exist at solid electrodes and at the nonconducting walls of an accelerator channel. Study of the two-dimensional flow in a channel is significantly simplified when such layers are present.     

Approximate solution of the problem of the flow of a conducting liquid in a duct with rectangular cross section

An approximate analytic solution is constructed for the problem of the fully developed stationary flow of a viscous incompressible liquid with a finite isotropic conductivity in a duct of rectangular cross section in the presence of an external magnetic field at right angles to nonconducting walls. An investigation is made of the extent to which flows in ducts with two electrodes parallel to the field resemble flows with four nonconducting walls. Theoretical and experimental investigations devoted to this problem are reviewed.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 33–40, September–October, 1970.     

The effects of the side walls on the flow of a second grade fluid in ducts with suction and injection

The effects of the side walls on the flow in ducts with suction and injection are examined. Three illustrative examples are given. The first example considers the effect of the side walls on the flow over a porous plate. The second example considers the flow between two parallel porous plates and the third example is devoted to the investigation of the flow in a rectangular duct with two porous walls. Exact solution of the governing equation using the no-slip boundary condition and an additional condition are obtained. The expression of the velocity, the volume flux and the vorticity are given. It is found that for large values of the cross-Reynolds number near the suction region the flow for a Newtonian fluid does not satisfy the boundary condition, but it does not behave in the same way for a second grade fluid. Three examples considered show that there are pronounced effects of the side walls on the flows of a second grade fluid in ducts with suction and injection.     

Higher order stability effects in a natural convection boundary layer over a vertical heated wall

 The stability of a laminar boundary layer flow under natural convection on a vertical isothermally heated wall is studied analytically. The analysis is performed by using two different two-dimensional linear models: (1) The non-parallel flow model in which the steady mean flow as well as the disturbance amplitude functions can change in the streamwise direction; (2) The parallel flow model in which the effects of the mean flow and disturbance changes in the streamwise direction are neglected. The linear non-parallel stability analysis is based on the so-called parabolised stability equations (PSEs) which have been successfully applied to the stability analysis of forced convection boundary layers. In this study the PSE equations are applied to natural convection boundary layers in order to show the difference between parallel and non-parallel stability analysis. A second part of this study deals with the effects of variable properties, which are always present in natural convection flows. They are analysed by an extended version of the Orr–Sommerfeld equation (EOSE). Received on 31 May 2000     

Hydromagnetic flow in a rotating channel

The steady flow in a parallel plate channel rotating with an angular velocity Ω and subjected to a constant transverse magnetic field is analysed. An exact solution of the governing equations is obtained. The solution in the dimensionless form contains two parameters: the Hartmann number, M ², and K ² which is the reciprocal of the Ekman number. The effects of these parameters on the velocity and magnetic field distributions are studied. For large values of the parameters, there arise thin boundary layers on the walls of the channel.     

An application of the method of multiple scales to magnetohydrodynamic flow

In this paper we study the flow of an incompressible conducting fluid along an elliptic duct imposed in an uniform magnetic field. In case Hartmann number of the flow is sufficiently large, the method of multiple scales is adopted for constructing the asymptotic approximation of solution up to any order. Our method can also be applied to study the magnetohydrodynamic flow along any duct whose cross section has smooth boundary.     

Magnetohydrodynamic flow on a half-plane

We investigate the magnetohydrodynamic flow (MHD) on the upper, half of a non-conducting plane for the case when the flow is driven by the current produced by an electrode placed in the middle of the plane. The applied magnetic field is perpendicular to the plane, the flow is laminar, uniform, steady and incompressible. An analytical solution has been developed for the velocity field and the induced magnetic field by reducing the problem to the solution of a Fredholm's integral equation of the second kind, which has been solved numerically. Infinite integrals occurring in the kernel of the integral equation and in the velocity and magnetic field were approximated for large Hartmann numbers by using Bessel functions. As the Hartmann number M increases, boundary layers are formed near the non-conducting boundaries and a parabolic boundary layer is developed in the interface region. Some graphs are given to show examples of this behaviour.     

Magnetohydrodynamic flow through a rectangular duct due to a periodic pressure gradient-II

An analysis is made of the flow of an incompressible, viscous, electrically conducting fluid in a long channel of rectangular cross section due to a periodic pressure gradient, in the presence of a uniform transverse magnetic field. Exact solutions are obtained and asymptotic forms valid for large Hartmann numbers in the boundary layers parallel to the field are discussed.     

Solution of the laminar duct flow problem by a boundary integral equation method

Summary A boundary integral equation method is proposed for approximate numerical and exact analytical solutions to fully developed incompressible laminar flow in straight ducts of multiply or simply connected cross-section. It is based on a direct reduction of the problem to the solution of a singular integral equation for the vorticity field in the cross section of the duct. For the numerical solution of the singular integral equation, a simple discretization of it along the cross-section boundary is used. It leads to satisfactory rapid convergency and to accurate results. The concept of hydrodynamic moment of inertia is introduced in order to easily calculate the flow rate, the main velocity, and the fRe-factor. As an example, the exact analytical and, comparatively, the approximate numerical solution of the problem of a circular pipe with two circular rods are presented. In the literature, this is the first non-trivial exact analytical solution of the problem for triply connected cross section domains. The solution to the Saint-Venant torsion problem, as a special case of the laminar duct-flow problem, is herein entirely incorporated.     

Magnetohydrodynamic flow in a rectangular duct

The magnetohydrodynamic flow of an incompressible, viscous, electrically conducting fluid in a rectangular duct, with an external magnetic field applied transverse to the flow, has been investigated. One of the duct's boundaries which is perpendicular to the magnetic field is taken partly insulated, partly conducting. An analytical solution has been developed for the velocity field and magnetic field by reducing the problem to the solution of a Fredholm integral equation of the second kind, which has been solved numerically. Solutions have been obtained for Hartmann numbers M up to 100. All the infinite series obtained are transformed to infinite integrals first and then to finite integrals which contain modified Bessel functions of the second kind. In this way, the difficulties associated with the computation of infinite integrals with oscillating integrands and slowly converging infinite series, the convergence of which is further affected for large values of M, have been avoided. It is found that, as M increases, boundary layers are formed near the non-conducting boundaries and in the interface region, and a stagnant region is developed in front of the conducting boundary for velocity field. The maximm value of magnetic field takes place on the conducting part. These behaviours are shown on some graphs.     

Diffusion and Electrical Processes in the Turbulent Boundary Layer and in the Neighborhood of the Stagnation Point

The problem of electric current (engine current) formation in aircraft jet engine ducts as a result of the development of electrical diffusion boundary layers on the surfaces of the duct and internal engine components is investigated. It is assumed that the outer flow containing electrons and positive ions is quasi-neutral and that the electrical quasi-neutrality is violated (and the electric engine current develops) in the wall flow zone as a result of the difference between the electron and ion diffusion coefficients. The problem of the development of an electrical diffusion boundary layer inside the turbulent gasdynamic boundary layer on a plane surface is formulated and solved. The engine current distribution along the duct is found for various values of a turbulent viscosity on the boundary of the gasdynamic boundary layer which affect the laminar-turbulent transition point.The electrical diffusion processes that occurs when an electrically quasi-neutral hydrodynamic stream impinges on a plane surface (simulation of the flow in the neighborhood of a stagnation point) is studied. In this case the Navier-Stokes equations have a self-similar solution. It is shown that the system of electrohydrodynamic equations also has a self-similar solution. The electrical parameter fields are determined and the engine current is found on the basis of this solution.     

PIV analysis of merging flow in a simplified model of a rotary kiln

Rotary kilns are used in a variety of industrial applications. The focus in this work is on characterizing the non-reacting, isothermal flow field in a rotary kiln used for iron ore pelletization. A downscaled, simplified model of the kiln is experimentally investigated using particle image velocimetry. Five different momentum flux ratios of the two inlet ducts to the kiln are investigated in order to evaluate its effect on the flow field in general and the recirculation zone in particular. Time-averaged and phase-averaged analyses are reported, and it is found that the flow field resembles that of two parallel merging jets, with the same characteristic flow zones. The back plate separating the inlet ducts acts as a bluff body to the flow and creates a region of reversed flow behind it. Due to the semicircular cross-section of the jets, the wake is elongated along the walls. Conclusions are that the flow field shows a dependence on momentum flux ratio of the jets; as the momentum flux ratio approaches unity, there is an increasing presence of von Kármán-type coherent structures with a Strouhal number of between 0.16 and 0.18. These large-scale structures enhance the mixing of the jets and also affect the size of the recirculation zone. It is also shown that the inclination of the upper inlet duct leads to a decrease in length of the recirculation zone in certain cases.     

Two-phase gas-liquid flow regimes for smooth and ribbed rectangular ducts

In this research, adiabatic two-phase air-water flow was investigated, and results for smooth and ribbed rectangular ducts are presented here. The test fluids were air and water at approximately atmospheric conditions. Three ribs of different heights were used; the rib width and pitch were held constant. The ribs were positioned in the duct at three different locations to establish three different conditions: on the bottom wall (water side), on the top wall (air side) and on both the top and bottom walls. The flow regimes in the smooth and ribbed ducts, which were recorded with a video camera, were classified as plug, stratified, slug and wavy flow. The location of the ribs in the duct did not alter the shape of the flow regimes, but the regime boundaries were considerably changed (repositioned). The effects of using ribs of different heights on regime boundaries are represented with flow map diagrams and discussed in detail. Compared to the smooth duct, the ribbed duct had different regime boundary positions. Increasing the rib height initiated hydrodynamical instability at lower fluid velocities. These findings are relevant for the operation and design of pipes, boilers and heat exchangers.     

Pulsed wire velocity anemometry near walls

A new pulsed wire probe for making velocity and turbulence measurements in the near wall region of incompressible, isothermal boundary layers of all kinds is described. Results of careful calibrations of the probe response in both laminar and turbulent flows are presented, with particular emphasis on the effects of diffusion in the very near wall region. Analytic results for the motion and distortion of a heat puff in linear shear flow near a wall are developed and these are shown to validate a very simple approximate theory that accounts for the diffusional effects. It is demonstrated that correction procedures based on the theory can be successfully implemented. Examples of the use of the probe in highly turbulent, separated flows, as well as more standard boundary layers, are given and its response near the wall is contrasted with that of the corresponding (parallel wire) probe used for surface shear stress measurements.