Fundamental solution for coupled magnetohydrodynamic flow equations

In this paper, a fundamental solution for the coupled convection–diffusion type equations is derived. The boundary element method (BEM) application then, is established with this fundamental solution, for solving the coupled equations of steady magnetohydrodynamic (MHD) duct flow in the presence of an external oblique magnetic field. Thus, it is possible to solve MHD duct flow problems with the most general form of wall conductivities and for large values of Hartmann number. The results for velocity and induced magnetic field is visualized in terms of graphics for values of Hartmann number M?300$M?300$.     

A stabilized finite element method for convection–diffusion problems

A stabilized finite element method (FEM) is presented for solving the convection–diffusion equation. We enrich the linear finite element space with local functions chosen according to the guidelines of the residual‐free bubble (RFB) FEM. In our approach, the bubble part of the solution (the microscales) is approximated via an adequate choice of discontinuous bubbles allowing static condensation. This leads to a streamline‐diffusion FEM with an explicit formula for the stability parameter τ_K that incorporates the flow direction, has the capability to deal with problems where there is substantial variation of the Péclet number, and gives the same limit as the RFB method. The method produces the same a priori error estimates that are typically obtained with streamline‐upwind Petrov/Galerkin and RFB. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2011     

A stable boundary elements method for magnetohydrodynamic channel flows at high Hartmann numbers

The article is devoted to extension of boundary element method (BEM) for solving coupled equations in velocity and induced magnetic field for time dependent magnetohydrodynamic (MHD) flows through a rectangular pipe. The BEM is equipped with finite difference approach to solve MHD problem at high Hartmann numbers up to 10⁶. In fact, the finite difference approach is used to approximate partial derivatives of unknown functions at boundary points respect to outward normal vector. It yields a numerical method with no singular boundary integrals. Besides, a new approach is suggested in this article where transforms 2D singular BEM's integrals to 1D nonsingular ones. The new approach reduces computational cost, significantly. Note that the stability of the numerical scheme is proved mathematically when computational domain is discretized uniformly and Hartmann number is 40 times bigger than length of boundary elements. Numerical examples show behavior of velocity and induced magnetic field across the sections.     

Local radial point interpolation method for the fully developed magnetohydrodynamic flow

In this paper, a local radial point interpolation method (LRPIM) is presented to obtain the numerical solutions of the coupled equations in velocity and magnetic field for the fully developed magnetohydrodynamic (MHD) flow through a straight duct of rectangular section with arbitrary wall conductivity and orientation of applied magnetic field. Local weak forms are developed using weighted residual method locally for the governing equations of fully developed MHD flow. The shape functions from LRPIM possess the delta function property. Therefore, essential boundary conditions can be applied as easily as that in the finite-element method. The implementation procedure of LRPIM method is node based, and it doesn’t need any “mesh” or “element”. Computations have been carried out for different Hartmann numbers, wall conductivities and orientations of applied magnetic field.     

High Hartmann number flow through a rectangular duct with unequally and finitely conducting walls

Nonsymmetric Hartmann flow through a rectangular duct is investigated for thin duct walls with, generally, unequal but finite conductivities. A high Hartmann number is adopted. Consistent with known phenomena, both Hartmann layers transverse to the applied magnetic field are assumed to be separated from the two side boundary layers by four corner regions plus four inner corner regions. The method of singular perturbations and matched asymptotic expansions is applied to the coupled system. The equations governing the core and Hartmann layers are first partially resolved for leading terms. This is then followed by tackling equations governing one side layer and two adjacent corner regions. The latters' incorporation secures, for the former, only those boundary conditions that are compatible along the transverse walls. Both corner regions are denied access to non-required boundary conditions along the neighbouring side wall by the adjoining inner corner regions. However, the latters' boundary value problems need not be tackled for the acquirement of only dominant terms beyond all four inner corner regions. The complementary side layer and associated corners are accounted for by a non-symmetric reflection principle. Results reveal that a difference between conductivities in the transverse walls together with at least one finitely conducting side wall impart to disturbances within the core and Hartmann layers (i) a nontrivial dependence on the transverse coordinate relative to the magnetic field and flow in addition to the (usual) dependence on the field aligned coordinate, (ii) a dependence on side wall parameters in addition to the dependence on transverse wall parameters. Applications to related situations are considered. These include the case for a perfectly conducting lower wall, a finitely conducting upper wall, and equally and finitely conducting side walls.     

Numerical Study of Creeping Liquid Metal Flow in the Presence of a Permanent Magnet

In this paper, we present the numerical study of interaction of a liquid metal flow with a small permanent magnet. Such an analysis is of fundamental interest for applications involving Lorentz Force Velocimetry (LFV). As a canonical problem, we consider the flow of liquid metal in the creeping flow regime (Re=0.01) through a square duct with electrically insulating walls. For this setup, we perform magnetohydrodynamic (MHD) simulations by coupling the commercial finite volume solver FLUENT and finite element solver COMSOL Multiphysics. Parametric analyses are performed at different relative positions of the permanent magnet with respect to the duct and with different magnetisation directions. The analyses provide good reference results for qualitative understanding of MHD flows exposed to non-uniform magnetic fields. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

ON THE STABILITY OF THE RESIDUAL-FREE BUBBLES FOR THE NAVIER-STOKES EQUATIONS

This paper considers the Calerkin finite element method for the incompressible Navier-Stokes equations in two dimensions, where the finite-dimensional spaces employed consist of piecewise polynomials enriched with residual-free bubble （RFB） functions. The stability features of the residual-free bubble functions for the linearized Navier-Stokes equations are analyzed in this work. It is shown that the enrichment of the velocity space by bubble functions stabilizes the numerical method for any value of the viscosity parameter for triangular elements and for values of the viscosity parameter in the vanishing limit case for quadrilateral elements.     

A finite volume spectral element method for solving magnetohydrodynamic (MHD) equations

In this paper, the coupled equations in velocity and magnetic field for unsteady magnetohydrodynamic (MHD) flow through a pipe of rectangular section are solved using combined finite volume method and spectral element technique, improved by means of Hermit interpolation. The transverse applied magnetic field may have an arbitrary orientation relative to the section of the pipe. The velocity and induced magnetic field are studied for various values of Hartmann number, wall conductivity and orientation of the applied magnetic field. Comparisons with the exact solution and also some other numerical methods are made in the special cases where the exact solution exists. The numerical results for these sample problems compare very well to analytical results.     

The boundary elements method for magneto-hydrodynamic (MHD) channel flows at high Hartmann numbers

In this article the constant and the continuous linear boundary elements methods (BEMs) are given to obtain the numerical solution of the coupled equations in velocity and induced magnetic field for the steady magneto-hydrodynamic (MHD) flow through a pipe of rectangular and circular sections having arbitrary conducting walls. In recent decades, the MHD problem has been solved using some variations of BEM for some special boundary conditions at moderate Hartmann numbers up to 300. In this paper we develop this technique for a general boundary condition (arbitrary wall conductivity) at Hartmann numbers up to 10⁵${10}^5$ by applying some new ideas. Numerical examples show the behavior of velocity and induced magnetic field across the sections. Results are also compared with the exact values and the results of some other numerical methods.     

Simulation and analysis of turbulent MHD channel flow with a streamwise magnetic field

We perform large-eddy simulations of turbulent MHD channel flow with a streamwise magnetic field using a pseudo spectral method. The streamwise magnetic field leads to turbulent drag reduction due to the selective Joule damping of certain flow structures. Near the walls, the turbulent mean velocity profile retains the logarithmic layer but the von Karman constant decreases with increasing magnetic field strength. In the outer region, the flow is characterized by persistent streaky structures of large streamwise extent, which lead to a rather flat mean velocity profile. In addition, the streamwise velocity fluctuations develop a pronounced second peak upon increasing the magnetic induction as well as a second logarithmic layer that increases in steepness. We find that Prandtl's classical mixing-length model with a variable Kármán constant can describe the modified logarithmic layer reasonably accurately in a wide range of Reynolds and Hartmann numbers. However, the flow modification near the center of the channel is not properly captured by this approach. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Application of Adapted-Bubbles to the Helmholtz Equation with Large Wavenumbers in 2D

An adapted-bubbles approach which is a modification of the residual-free bubbles (RFB) method, is proposed for the Helmholtz problem in 2D. A new two-level finite element method is introduced for the approximations of the bubble functions. Unlike the other equations such as the advection-diffusion equation, RFB method when applied to the Helmholtz equation, does not depend on another stabilized method to obtain approximations to the solutions of the sub-problems. Adapted-bubbles (AB) are obtained by a simple modification of the sub-problems. This modification increases the accuracy of the numerical solution impressively. We provide numerical experiments with the AB method up to $ch = 5$ where $c$ is the wavenumber and $h$ is the mesh size. Numerical tests show that the AB method is better by far than higher order methods available in the literature.     

Enhanced RFB method

The residual-free bubble method (RFB) is a parameter-free stable finite element method that has been successfully applied to a wide range of boundary-value problems exhibiting multiple-scale behaviour. If some local features of the solution are known a priori, the approximation properties of the RFB finite element space can be improved by enriching it on selected edges of the partition by edge-bubbles that are supported on pairs of neighbouring elements. Motivated by this idea, we define and analyse an enhanced residual-free bubble method for the solution of convection-dominated convection-diffusion problems in 2-D. Our a priori analysis highlights the limitations of the RFB method and the improved global approximation properties of the new method. The theoretical results are supported by detailed numerical experiments.     

Magneto-hydrodynamic stability of plane Poiseuille flow with flexible boundaries

MHD stability of plane Poiseuille flow between parallel flexible walls with coplaner magnetic field is analysed. The study is restricted to sufficiently low values of magnetic Reynolds number. The eigen-value problem so posed is then solved graphically and neutral curves are obtained for various sets of magnetic parameter and flexible wall parameters. The nature of influence on flow stability depends on the values of these parameters.     

可渗透壁面上Falkner-Skan磁流体动力学流动的近似解

研究了可渗透壁面上Falkner-Skan磁流体动力学（MHD）边界层流动问题．利用结合了微分变换法（DTM）和Padé近似的DTM-Padé方法,得到了边界层问题的近似解和壁摩擦因数值．通过建立一个迭代程序,边界层问题的近似解被表示为幂级数的形式,而且以图和表形式对不同参数下的近似解结果与打靶法得到的数值结果进行了对比,近似解和数值解结果高度吻合,从而验证了所得问题近似解和结论的可靠性和有效性．并且,对求得的边界层问题近似解结果进行了讨论,分析了不同物理参数对边界层流动的影响．     

Subcritical instability of liquid metal channel flow in the presence of a spanwise magnetic field

The linear and nonlinear evolution of perturbations is investigated in a magnetohydrodynamic channel flow with electrically insulating walls. The applied magnetic field is parallel to the walls and orthogonal to the stream. Linear optimal perturbations and their maximum amplifications over finite time intervals are computed using a scheme based on the direct and adjoint governing equations. It is shown that dominant optimal perturbations are no more the classical streamwise modes and how the flow is two-dimenzionalized for high enough Hartmann numbers. For fixed Reynolds and Hartmann numbers, direct numerical simulations are applied to investigate how the transition to turbulence is affected by the magnetic field. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Comparison of the simplified and full MHD models for laminar incompressible flow past a circular cylinder

In the majority of research on incompressible magnetohydrodynamic (MHD) flows, the simplified model with the low magnetic Reynolds number assumption has been adopted because it reduces the number of equations to be solved. However, because the effect of flow on magnetic field is also neglected, the solutions of the simplified model may be different from those of the full model. As an example, the flow of an electrically conducting fluid past a circular cylinder under a magnetic field is investigated numerically using the simplified and full models in this paper. To solve the problems, two second-order compact finite difference algorithms based on the streamfunction-velocity formulation of the simplified model and the quasi-streamfunction-velocity formulation of the full model are developed respectively.Numerical simulations are carried out over a wide range of Hartmann number for steady-state laminar problems with both models. For the full model, magnetic Reynolds number (Re_m) is chosen from 0.01 to 10. The computed results show that solutions of the simplified MHD model are not exactly the same as those of the full MHD model for this flow problem in most cases even if Re_m in the full model is very low. Only in the special case that a strong external magnetic field is exerted perpendicular to the dominant flow direction, can the simplified MHD model be regarded as an approximation of the full MHD model with low Re_m.     

Semi‐analytical solution of MHD flow through boundary integrals on the pipe wall

A mathematical model is given for the magnetohydrodynamic (MHD) pipe flow as an inner Dirichlet problem in a 2D circular cross section of the pipe, coupled with an outer Dirichlet or Neumann magnetic problem. Inner Dirichlet problem is given as the coupled convection‐diffusion equations for the velocity and the induced current of the fluid coupling also to the outer problem, which is defined with the Laplace equation for the induced magnetic field of the exterior region with either Dirichlet or Neumann boundary condition. Unique solution of inner Dirichlet problem is obtained theoretically reducing it into two boundary integral equations defined on the boundary by using the corresponding fundamental solutions. Exterior solution is also given theoretically on the pipe wall with Poisson integral, and it is unique with Dirichlet boundary condition but exists with an additive constant obtained through coupled boundary and solvability conditions in Neumann wall condition. The collocation method is used to discretize these boundary integrals on the pipe wall. Thus, the proposed procedure is an improved theoretical analysis for combining the solution methods for the interior and exterior regions, which are consolidated numerically showing the flow behavior. The solution is simulated for several values of problem parameters, and the well‐known MHD characteristics are observed inside the pipe for increasing values of Hartmann number maintaining the continuity of induced currents on the pipe wall.     

Direct numerical simulation of transition in MHD duct flow

Transition in the flow of electrically conducting fluid in a square duct with insulating walls is studied by direct numerical simulations. A uniform magnetic field is applied in the transverse direction. Moderate values of the Reynolds (Re = 5000 ) and Hartmann (Ha = 0 … 30 ) numbers are considered that correspond to the classical Hartmann & Lazarus [1] experiments. It is shown that the laminarization begins in the Hartmann layers, whereas the sidewall layers remain turbulent. Complete re-laminarization occurs in the range of R = Re/Ha ≈︁ 220 , which is in agreement with the H. & L. experiments. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Solitary fluid transients in rectangular ducts with transverse magnetic fields,II

Summary This paper treats a fluid hammer wave which is propagating into a fully devloped MHD duct flow. The wave is produced by suddenly closing a valve at some cross section of a rectangular duct with a uniform, transverse, applied magnetic field, with perfectly conducting walls parallel to the field and with either insulating or perfectly conducting walls perpendicular to the field. The Mach and magnetic Reynolds numbers are assumed to be small, while the Hartmann number is assumed to be large. The jump in velocity and pressure across the wave decreases exponentially in time. The fully developed flow ahead of the wave is undisturbed, and solutions for the velocity and pressure between the valve and the wave are presented.

---

Résumé Cet article traite de l'onde de choc produite par la fermeture subite d'une vanne et qui se propage dans un fluide conducteur d'électricité s'écoulant dans une conduite de section rectangulaire à laquelle est appliquée un champ magnétique transversal et homogène. Deux parois de la conduite sont parfaitement conductrices et parallèles au champ magnétique, les deux autres étant parfaitement conductrices ou isolantes. On admet que le nombre de Mach et le nombre magnétique de Reynolds sont petits et que le nombre de Hartmann est grand. Les changements de la vitesse et de la pression à travers l'onde diminuent exponentiellement à temps. L'écoulement à la tête de l'onde n'est pas perturbée. Des solutions pour la vitesse et la pression entre la vanne et l'onde sont presentées.

     

Bubble—函数稳定化混合有限元分析

本文针对混合结构抽象问题，基于「９」的非标准稳定化有限元方法的一般框架研究了ｂｕｂｂｌｅ－函数稳定化方法，该逼近代格式使得Ｂａｂｕｓｋａ－Ｂｒｅｚｚｉ条件是不必要的。