Vortex Generation in a Liquid Metal Duct Flow near a Magnetic Dipole

The generation of vortical structures by a strong magnetic dipole field in a liquid metal duct flow is studied by means of three-dimensional direct numerical simulations. The dipole is considered as the paradigm for a magnetic obstacle which will deviate the streamlines due to Lorentz forces which act on the fluid elements. Our model uses the quasi-static approximation applicable in the limit of small magnetic Reynolds numbers. The analysis covers the stationary flow regime at smaller flow Reynolds numbers Re as well as the fully time-dependent regimes at higher values with a turbulent flow in the wake of the magnetic obstacle. We present a systematic study of these two basic flow regimes on Re and the Hartmann number Ha, a measure of the strength of the magnetic dipole field. Furthermore, three orientations of the dipole are compared, the streamwise, spanwise and wall-normal ones. The most efficient generation of turbulence at a fixed distance above the duct occurs for the spanwise orientation in which we can observe the formation of Hartmann layers at the top plate. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Turbulent drag reduction at moderate Reynolds numbers via spanwise velocity waves

The question whether large turbulent drag reduction can be achieved at the high values of Re typical of applications is addressed. Answering such question, either by experiments or DNS, is obviously challenging. For DNS, the problem lies in the tremendous increase of the computational cost with Re, that has to be appreciated in view of the need of carrying out an entire parametric study at every Re, owing to the unknown location of the optimal forcing parameters. In this paper we limit ourselves to considering an open-loop technique based on spanwise forcing, the streamwise-traveling waves introduced by [1], and explore via Direct Numerical Simulations (DNS) how the drag reduction varies when the friction Reynolds number is increased from Re_τ = 200 to Re_τ = 2000. To achieve high Re while keeping the computational cost affordable, computational domains of reduced size are employed. We adopted special care to interpret results that are indeed still box-size dependent, as well as strategies to compute the random errors and give the results an error bar. Our results indicate that still R = 0.29 can be obtained at Re_τ = 2000 in the partial region of the parameter space studied. The maximum R is found to decrease as R ˜ Re_τ^−0.22 in the Reynolds range investigated. As most important outcome, we find that the sensitivity of R to Re becomes smaller when far from the low-Re optimum parameters: in this region, we suggest R ˜ Re_τ^−0.08. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Singularities on the boundaries of magnetorotational instabilities and scaling laws

In the theory of magnetorotational instability and its modern extensions such as the helical MRI, non-trivial scaling laws between the critical parameters are observed. In case of the standard MRI it is well known that the Reynolds and Hartmann numbers are scaled as Re ∼ Ha² while for the helical MRI Re ∼ Ha³ . What is less known is that the thresholds of SMRI and HMRI plotted as surfaces in the space of parameters, possess singularities that determine the scaling laws. Moreover, the two paradoxes of SMRI and HMRI in the limits of infinite and zero magnetic Prandtl number (Pm), respectively, sharply correspond to the singularities on the instability thresholds. In either case, it is the local Plücker conoid structure that explains the non-uniqueness of the critical Rossby number, and its crucial dependence on the Lundquist number. For HMRI, we have found an extension of the former Liu limit Ro_c ≃ −0.828 (valid for Lu = 0 ) to a somewhat higher value Ro ≃ −0.802 at Lu = 0.618 which is, however, still below the Kepler value. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Magnetohydrodynamic stationary and oscillatory convective stability in a mushy layer during binary alloy solidification

For the case of solidification of a bottom cooled binary alloy, the magnetohydrodynamic stationary and oscillatory convective stability in the mushy layer is investigated analytically using normal mode linear stability analysis. In the limit of large Stefan number (S_t), a near–eutectic approximation with large far field temperature is considered in the present research. To ascertain the instability in the mushy layer, the strength of the superimposed magnetic field is so chosen that it corresponds to a given mush Hartmann number (Ha_m) of the problem. The results are presented for various values of mush Hartmann numbers in the range, 0 ≤ Ha_m ≤ 50. The critical Rayleigh number for stationary convection shows a linear relationship with increasing Ha_m. The magnetohydrodynamic effect imparts a stabilizing influence during stationary convection. In comparison to that of the stationary convective mode, the oscillatory mode appears to be critically susceptible at higher values of β (β = S_t/℘² ϒ², ℘ is the compositional ratio, ϒ = 1 + S_t/℘), and vice versa for lower β values. Analogous to the behavior for stationary convection, the magnetic field also offers a stabilizing effect in oscillatory convection and thus influences global stability of the mushy layer. Increasing magnetic strength shows reduction in the wavenumber and in the number of rolls formed in the mushy layer.     

Curve Fitting, Differential Equations And The Riemann Hypothesis

It is known that the Riemann zeta function ζ (s) in the critical strip 0 < Re(s) < 1, may be represented as the Mellin transform of a certain function φ (x) which is related to one of the theta functions. The function φ (x) satisfies a well known functional equation, and guided by this property we deduce a family of approximating functions involving an arbitrary parameter α. The approximating function corresponding to the value of α = 2 gives rise to a particularly accurate numerical approximation to the function φ (x). Another approximation to φ (x), which is based upon the first one, is obtained by solving a certain differential equation. Yet another approximating function may be determined as a simple extension of the first. All three approximations, when used in conjunction with the Mellin transform expression for ζ (s) in the critical strip, give rise to an explicit expression from which it is clear that Re(s) = 1/2 is a necessary and sufficient condition for the vanishing of the imaginary part of the integral, the real part of which is non-zero. Accordingly, the analogy with the Riemann hypothesis is only partial, but nevertheless Re(s) = 1/2 emerges from the analysis in a fairly explicit manner. While it is generally known that the imaginary part of the Mellin transform must vanish along Re(s) = 1/2, the major contribution of this paper is the presentation of the actual calculation for three functions which approximate φ (x). The explicit nature of these calculation details may facilitate progress towards the corresponding calculation for the actual φ (x), which may be necessary in a resolution of the Riemann hypothesis.2000 Mathematics Subject Classification: Primary—11M06, 11M26     

Numerical simulation of magnetohydrodynamic buoyancy-induced flow in a non-isothermally heated square enclosure

Numerical simulation of magnetohydrodynamic (MHD) buoyancy-induced heat transfer and fluid flow has been analyzed in a non-isothermally heated square enclosure using finite volume method. The bottom wall of enclosure were heated and cooled with a sinusoidal function and top wall was cooled isothermally. Vertical walls of the enclosure were adiabatic. Effects of Rayleigh number (Ra = 10⁴, 10⁵ and 10⁶), Hartman number (Ha = 0, 50 and 100) and amplitude of sinusoidal function (n = 0.25, 0.5 and 1) on temperature and flow fields were analyzed. It was observed that heat transfer was decreased with increasing Hartmann number and decreasing value of amplitude of sinusoidal function.     

Analysis of a conducting fluid in a thin annulus with rotating insulated walls under radial magnetic effect

This work considers an electrically conducting fluid filled between two concentric cylindrical walls relatively close to each other. A theoretical solution for the steady Taylor–Couette flow between these two electrically insulated rotating cylinders under the influence of a radial magnetic field is provided in this work. By solving the appropriate set of governing equations simultaneously, the profiles of fluid tangential velocity component and induced magnetic field were obtained as complicated functions involving the modified Bessel functions of the first and second kinds of the first-order in terms of radial coordinates and Hartmann number. A computational study was also performed to validate the present theoretical solution. The analytical and computational results are identical when Ha = 1 while these results only slightly deviate from each other as Ha increases. Current results show that, the presence of the external magnetic field causes the flow close to the slower cylinder to accelerate while that close to the faster cylinder to decelerate. This has clearly implied the fact that an external magnetic field tends to make the velocity distribution across the inner and outer cylinders more uniform.     

Determining the limits of geometrical tortuosity from seepage flow calculations in porous media.

Recent investigations have found a distinct correlation of effective properties of porous media to sigmoidal functions, where one axis is the Reynolds number Re and the other is the effective property dependent of Re, Θ_i = S_i(Re). One of these properties is tortuosity. At very low Re (seepage flow), there is a characteristic value of tortuosity, and it is the upper horizontal asymptote of the sigmoidal function. With higher values of Re (transient flow) the tortuosity value decreases, until a lower asymptote is reached (turbulent flow). Estimations of this parameter have been limited to the low Reynolds regime in the study of porous media. The current state of the art presents different numerical measurements of tortuosity, such as skeletization, centroid binding, and arc length of streamlines. These are solutions for the low Re regime. So far, for high Re, only the arc length of stream lines has been used to calculate tortuosity. The present approach involves the simulation of fluid flow in large domains and high Re, which requires numerous resources, and often presents convergence problems. In response to this, we propose a geometrical method to estimate the limit of tortuosity of porous media at Re → ∞, from the streamlines calculated at low Re limit. We test our method by calculating the tortuosity limits in a fibrous porous media, and comparing the estimated values with numerical benchmark results. Ongoing work includes the geometric estimation of different intrinsic properties of porous media. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Investigating Asymptotic Suction Boundary Layers using a One-Dimensional Stochastic Turbulence Model

The one-dimensional turbulence model (ODT) is applied to study turbulent asymptotic suction boundary layers for a Reynolds number of Re = u_∞/v₀ = 333, where u_∞ and v₀ are the free stream and suction velocity, respectively. In here we will demonstrate that a large eddy suppression mechanism may reduce the influence of ODT model parameters, such as the viscous cut-off parameter Z. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On normal forms of singular Levi-flat real analytic hypersurfaces

Let F(z) = Re(P(z)) + h.o.t be such that M = (F = 0) defines a germ of real analytic Levi-flat at 0 ∈ ℂ ⁿ , n ≥ 2, where P (z) is a homogeneous polynomial of degree k with an isolated singularity at 0 ∈ ℂ ⁿ and Milnor number μ. We prove that there exists a holomorphic change of coordinate ϕ such that ϕ(M) = (Re(h) = 0), where h(z) is a polynomial of degree μ + 1 and j ₀ ^k (h) = P.     

The Trailing Edge Problem For Mixed Convection Flow Past a Horizontal Plate

The influence of buoyancy onto the boundary‐layer flow past a horizontal plate aligned parallel to a uniform free stream is characterized by the buoyancy parameter K = Gr/Re^5/2 where Gr and Re are the Grashof and Reynolds number, respectively. An asymptotiy analysis of the complete flow field including potential flow, boundary layer, wake and interaction region is given for small buoyancy parameters and large Reynolds numbers in the distinguished limit KRe^1/4 = O(1). (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Singular perturbations for the exterior three-dimensional slow viscous flow problem

In this paper the rigorous justification of the formal asymptotic expansions constructed by the method of matched inner and outer expansions is established for the three-dimensional steady flow of a viscous, incompressible fluid past an arbitrary obstacle. The justification is based on the series representation of the solution to the Navier-Stokes equations due to Finn, and it involves the reductions of various exterior boundary value problems for the Stokes and Oseen equations to boundary integral equations of the first kind from which existence as well as asymptotic error estimates for the solutions are deduced. In particular, it is shown that the force exerted on the obstacle by the fluid admits the asymptotic representation F = A₀ + A₁Re + O(Re² ln Re⁻¹) as the Reynolds number Re → 0⁺, where the vectors A₀ and A₁ can be obtained from the method of matched inner and outer expansions.     

DNS of a Turbulent Channel Flow with Streamwise Rotation - Study of the Reverse Effect of the Cross Flow

Modeling of rotating turbulent flows is a major issue in engineering applications. In this work a turbulent channel flow rotating about the streamwise direction is presented. The theory is based on the investigations of [3] employing Lie group analysis. It was found that a cross flow in spanwise direction is induced. A series of direct numerical simulations (DNS) has been conducted for both different rotation rates and different Reynolds numbers to validate the cross flow. In addition some new interesting effects were observed. The averaged profile ū₃ of the cross flow is formed like a ‘S’ that means it exhibits a triple zero-crossing which denotes regions of reverse flow. Alaso a reverse effect is seen which means that for small rotation rates up to Ro=10 the spanwise mean velocity profiles increase and at rotation number Ro=14 this effect appears to reverse. Both effects were observed at two different Reynolds numbers Re = 180 and Re = 270. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Longitudinal and transversal two-point correlations in a turbulent flow

Based on measurements with hot-wire anemometers in a fully developed turbulent flow in the wake of a cylinder we investigate different aspects of longitudinal and transversal two-point correlations. In particular we use Taylor's hypothesis as well as simultaneous measurements with two X-hot-wire probes, and investigate how these correlations change with the relative orientation of their separation vector with respect to the mean flow direction. For the wake at a distance of eighty cylinder diameters at a Reynolds number Re = 28 000 (Taylor-based Reynolds number Re_λ = 290), we find evidence of strong anisotropy. We analyze the data using increment statistics, in particular increment distributions and structure functions. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Unsteady magnetohydrodynamic Hartmann–Couette flow and heat transfer in a Darcian channel with Hall current,ionslip, viscous and Joule heating effects: Network numerical solutions

A theoretical study of unsteady magnetohydrodynamic viscous Hartmann–Couette laminar flow and heat transfer in a Darcian porous medium intercalated between parallel plates, under a constant pressure gradient is presented. Viscous dissipation, Joule heating, Hall current and ionslip current effects are included as is lateral mass flux at both plates. The dimensionless conservation equations for the primary (x¹-direction), secondary (z¹-direction) momentum and also energy conservation equation are derived and solved using a computational technique known as Network Simulation Methodology (NSM). Velocity distributions (u¹, w¹) and temperature distribution (T¹) at the channel centre (y¹ = 0) over time (t¹) are studied graphically for the effects of Darcy number (Da), Hartmann number (Ha), transpiration (N_t), Hall current parameter (Be), ionslip parameter (Bi), pressure gradient parameter (dP/dx¹) with Prandtl number prescribed at 7.0 (electrically conducting water), Eckert number held constant at 0.25 (heat convection from the plates to the fluid) and Reynolds number (Re) fixed at 5.0 (for Re < 10, Darcian model is generally valid). Increasing Darcy number causes an increase in temperature, T¹; values are however significantly reduced for the higher Hartmann number case (Ha = 10). For the case of low transpiration (i.e. N_t = 1 which corresponds to weak suction at the upper plate and weak injection at the lower plate), both primary velocity (u¹) and secondary velocity (w¹) are increased with a rise in Darcy number (owing to a simultaneous decrease in Darcian porous drag); temperature T¹ is also increased considerably with increasing Da. However, for stronger transpiration (N_t = 10), magnitudes of u¹, w¹ and T¹ are significantly reduced and also significant overshoots are detected prior to the establishment of steady state flow. With increasing Hall current parameter, Be, (for the purely fluid regime i.e. Da → ∞), primary velocity is considerably increased, whereas secondary velocity is reduced; temperatures are decreased in the early stages of flow but effectively increased in the steady state with increasing Be. With strong Darcian drag present (Da = 0.01 i.e. very low permeability), magnitudes of u¹, w¹ and T¹ are considerably reduced and temperatures are found to be reduced for all t¹, with increasing Hall current effect (Be). Increasing ionslip current parameter (Bi) increases primary velocity (u¹), decreases secondary velocity (w¹) and also temperature (T¹) for all time (t¹), in the infinite permeability case (Da → ∞). For weakly Darcian flow, ionslip parameter (Bi) has a much reduced effect on the velocity distributions. Temperature, T¹ is strongly increased with a rise in pressure gradient parameter, dP¹/dx¹, as is primary velocity (u¹); however, secondary velocity (w¹) is reduced. The present study has applications in hybrid magnetohydrodynamic (MHD) energy generators, materials processing, geophysical hydromagnetics, etc.     

Structure of bounded solutions for a class of abstract differential equations

Summary The paper is concerned with bounded solutions of an equationu′(t)=Bu(t) in Hilbert spaces, . A representation formula is obtained depending on the zeros of Rez(θ).

---

Riassunto è studiata la struttura delle soluzioni limitate di una equazione differenzialeu′(t)=Bu(t) in uno spazio di Hilbert, ove , in funzione dei zeri di Rez(θ).

---

---

This research is supported through a grant of the National Research Council Canada.     

Numerical study of magnetohydrodynamic viscous plasma flow in rotating porous media with Hall currents and inclined magnetic field influence

A numerical solution is developed for the viscous, incompressible, magnetohydrodynamic flow in a rotating channel comprising two infinite parallel plates and containing a Darcian porous medium, the plates lying in the x–z plane, under constant pressure gradient. The system is subjected to a strong, inclined magnetic field orientated to the positive direction of the y-axis (rotational axis, normal to the x–z plane). The Navier–Stokes flow equations for a general rotating hydromagnetic flow are reduced to a pair of linear, viscous partial differential equations neglecting convective acceleration terms, for primary velocity (u′) and secondary velocity (v′) where these velocities are directed along the x and y axes. Only viscous terms are retained in the momenta equations. The model is non-dimensionalized and shown to be controlled by a number of dimensionless parameters. The resulting dimensionless ordinary differential equations are solved using a robust numerical method, Network Simulation Methodology. Full details of the numerics are provided. The present solutions are also benchmarked against the analytical solutions presented recently by Ghosh and Pop [Ghosh SK, Pop I. An analytical approach to MHD plasma behaviour of a rotating environment in the presence of an inclined magnetic field as compared to excitation frequency. Int J Appl Mech Eng 2006;11(4):845–856] for the case of a purely fluid medium (infinite permeability). We study graphically the influence of Hartmann number (Ha, magnetic field parameter), Ekman number (Ek, rotation parameter), Hall current parameter (Nh), Darcy number (Da, permeability parameter), pressure gradient (Np) and also magnetic field inclination (θ) on primary and secondary velocity fields. Additionally we investigate the effects of these multiphysical parameters on the dimensionless shear stresses at the plates. Both primary and secondary velocity are seen to be increased with a rise in Darcy number, owing to a simultaneous reduction in Darcian drag force. Primary velocity is seen to decrease with an increase in Hall current parameter (Nh) but there is a decrease in secondary velocity. The study finds important applications in magnetic materials processing, hydromagnetic plasma energy generators, magneto-geophysics and planetary astrophysics.     

Characteristic features of the dynamics of the Ginzburg-Landau equation in a plane domain

We study the boundary value problem w_t=ℵ₀Δw+ℵ₁w-ℵ₂w|w|²,w|∂Ω₀=0 in the domain Ω₀={(x,y):0 ≤ x ≤ l₁,0 ≤ y ≤ l₂}. Here, w is a complex-valued function, Δ is the laplace operator, and ℵ_j, j=0,1,2, are complex constants withRe ℵ_j > 0. We show that under a rather general choice of the parameters l₁ and l₂, the number of stable invariant tori in the problem, as well as their dimensions, grows infinitely asRe ℵ₀ → 0 andRe ℵ₀ → 0. Translated from Teoreticheskaya i Matematicheskaya Fizika, Vol. 125, No. 2, pp. 205–220, November, 2000.     

Effect of aligned magnetic field on the steady viscous flow past a circular cylinder

The steady viscous incompressible and slightly conducting fluid flow around a circular cylinder with an aligned magnetic field is simulated for the range of Reynolds numbers 100 ? Re ? 500 using the Hartmann number, M. The multigrid method with defect correction technique is used to achieve the second order accurate solution of complete non-linear Navier–Stokes equations. The magnetic Reynolds number is assumed to be small. It is observed that volume of the separation bubble decreases and drag coefficient increases as M is increased. We noticed that the upstream base pressure increases slightly with increase of M whereas downstream base pressure decreases with increase of M. The effect of the magnetic field on the flow is discussed with contours of streamlines, vorticity, plots of surface pressure and surface vorticity.