Comparison study on the numerical stability of dual reciprocity boundary element method for the MHD slip flow problem

ŞENEL P.

Engineering Analysis with Boundary Elements, cilt.151, ss.370-386, 2023 (SCI-Expanded) identifier identifier

  • Yayın Türü: Makale / Tam Makale
  • Cilt numarası: 151
  • Basım Tarihi: 2023
  • Doi Numarası: 10.1016/j.enganabound.2023.03.010
  • 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.370-386
  • Anahtar Kelimeler: MHD, Wall slip, Unsteady flow, Dual reciprocity boundary element method, Numerical stability, MAGNETOHYDRODYNAMIC FLOW, DRBEM SOLUTION, PIPE, SIMULATION, CAVITY, WALLS, DUCT
  • Karadeniz Teknik Üniversitesi Adresli: Evet

Özet

In this paper, unsteady Magnetohydrodynamics (MHD) flow in a cavity with insulated and slipping walls is studied. The nonlinear coupled MHD equations containing the velocity and the induced magnetic field are solved by the dual reciprocity boundary element method (DRBEM). The time derivative is discretized by the explicit Euler or the central difference methods. The grid independence is tested. The numerical stability is analyzed through the spectral radius of the coefficient matrix and the performances of the two time discretization methods are compared. In the studied ranges, the proper choices of relaxation parameter and time increment are found. The influences of Hartmann number (Ha), Reynolds number (Re), magnetic Reynolds number (Rm), and the slip length on the numerical stability are investigated. The flow behaviors in the steady-state and transient time levels are presented. It is found that relaxation parameter and time step close to one guarantee the numerical stability in the studied range. Increasing Re alters the stability condition more than Rm increment. DRBEM with the central difference method is the appropriate numerical technique considering the computational cost and the stability for solving MHD slip flow problems. The flow elongation time is postponed for large Rm and Re.