The method of variably scaled radial kernels for solving two-dimensional magnetohydrodynamic (MHD) equations using two discretizations: The CrankNicolson scheme and the method of lines (MOL)
MHD equations have many applications in physics and engineering. The model is coupled equations in velocity and magnetic field and has a parameter namely Hartmann. The value of Hartmann number plays an important role in the equations. When this parameter increases, using different meshless methods makes the oscillations in velocity near the boundary layers in the region of the problem. In the present paper a numerical meshless method based on radial basis functions (RBFs) is provided to solve MHD equations. For approximating the spatial variable, a new approach which is introduced by Bozzini etal. (2015) is applied. The method will be used here is based on the interpolation with variably scaled kernels. The methodology of the new technique is defining the scale function c c on the domain ??Rd ? ? R d . Then the interpolation problem from the data locations x j ? R d transforms to the new interpolation problem in the data locations ( x j , c ( x j ) ) ? R d + 1 (Bozzini etal., 2015). The radial kernels used in the current work are Multiquadrics (MQ), Inverse Quadric (IQ) and Wendlands function. Of course the latter one is based on compactly supported functions. To discretize the time variable, two techniques are applied. One of them is the CrankNicolson scheme and another one is based on MOL. The numerical simulations have been carried out on the square and elliptical ducts and the obtained numerical results show the ability of the new method for solving this problem. Also in appendix, we provide a computational algorithm for implementing the new technique in MATLAB software.
variably scaled radial kernel wendlands function two-dimensional unsteady mhd equations meshless collocation method cranknicolson scheme and mol differential algebraic equations (daes)