DRBEM solution of singularly perturbed coupled MHD flow equations

ARSLAN ÖLÇER S., TEZER M.

Engineering Analysis with Boundary Elements, cilt.155, ss.696-706, 2023 (SCI-Expanded) identifier identifier

  • Yayın Türü: Makale / Tam Makale
  • Cilt numarası: 155
  • Basım Tarihi: 2023
  • Doi Numarası: 10.1016/j.enganabound.2023.06.038
  • Dergi Adı: Engineering Analysis with Boundary Elements
  • Derginin Tarandığı İndeksler: Science Citation Index Expanded (SCI-EXPANDED), Scopus, Academic Search Premier, Aerospace Database, Aquatic Science & Fisheries Abstracts (ASFA), Communication Abstracts, INSPEC, Metadex, zbMATH, DIALNET, Civil Engineering Abstracts
  • Sayfa Sayıları: ss.696-706
  • Anahtar Kelimeler: Boundary layer, DRBEM, Hartmann number, MHD flow, Singular perturbation, Transition point
  • Orta Doğu Teknik Üniversitesi Adresli: Evet

Özet

In this study, the numerical solution of singularly perturbed magnetohydrodynamic (MHD) flow in a square duct with no-slip, and insulated or perfectly conducting walls is investigated by using the Dual Reciprocity Boundary Element Method (DRBEM). The steady, laminar and fully-developed MHD flow of an incompressible, viscous, and electrically conducting fluid in a long channel of square cross-section (duct) is driven by a pressure gradient. The governing MHD flow equations are convection–diffusion type and coupled in terms of the velocity V(x,y) and induced magnetic field B(x,y). When the intensity of horizontally applied external magnetic field is high, the Hartmann number (Ha) which is the coefficient of convection terms becomes large, that is, the coupled MHD flow equations become convection dominated. In other words, the coefficient of the diffusion terms is very small giving singularly perturbed MHD duct flow which exhibits thin boundary layers. The numerical scheme uses Shishkin mesh which depends on the number of points on one side and Ha. Numerical results obtained by DRBEM reveal that the well-known characteristics of the MHD flow problem are found as Ha increases and it is possible to obtain V(x,y) and B(x,y) for large values of Hartmann number up to Ha=1000.