Next Article in Journal
Blitz Vision: Development of a New Full-Electric Sports Sedan Using QFD, SDE and Virtual Prototyping
Next Article in Special Issue
Mass Transfer Effects on the Mucus Fluid with Pulsatile Flow Influence of the Electromagnetic Field
Previous Article in Journal
A Fingerprint Matching Algorithm Using the Combination of Edge Features and Convolution Neural Networks
Previous Article in Special Issue
An Investigation on the Vortex Effect of a CALM Buoy under Water Waves Using Computational Fluid Dynamics (CFD)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Solution of Steady Incompressible MHD Problems with Quasi-Least Square Method

1
Department of Mathematics, University of Baltistan, Skardu 16200, Pakistan
2
Department of Mathematical Sciences, Karakorm International University, Gilgit 15100, Pakistan
3
Ghizer Campus, Karakorm International University, Gahkuch 15200, Pakistan
*
Author to whom correspondence should be addressed.
Inventions 2022, 7(2), 40; https://doi.org/10.3390/inventions7020040
Submission received: 26 April 2022 / Revised: 30 May 2022 / Accepted: 31 May 2022 / Published: 6 June 2022
(This article belongs to the Special Issue Recent Trends in Nanofluids III)

Abstract

:
A quasi-least-squares (QLS) mixed finite element method (MFE) based on the L 2 -inner product is utilized to solve an incompressible magnetohydrodynamic (MHD) model. These models are associated with the three unknown terms, i.e., fluid velocity, fluid pressure, and magnetic field. For the MHD-based models, common theories and algorithms for approximation of the solutions are not always applicable because of the choice of the functional spaces during the utilization of the weak formulation. It is well known that the spaces used for the approximation of the different unknowns, e.g., the spaces for the unknowns, cannot be chosen independently for the variational formulation, and may have to satisfy strict stability conditions such as the inf-sup, or Ladyzhenskaya–Babuska–Brezzi (LBB) condition. The dependency of the selection of the spaces for the unknowns are critical and always not applicable for some pair of unknowns. Because of this, the numerical or theoretical solutions must have to face some stability issue. The proposed scheme ( L 2 -inner product) is introduced to circumvent this deficiency of the conditions (inf-sup or LBB) and obtained a well-posed solution theoretically. The model equations are nonlinear and highly coupled with the combination of Navier–Stokes and Maxwell relations. First, these nonlinear models are made linear around a specific state wherein the modified system represents an algebraic equation in a first-order symmetric form. Secondly, a direct iteration technique is applied to solve the nonlinearities and obtain a theoretical convergent rate for a general initial guess. Theoretical results show that only a single parameter with a single initial guess is sufficient to establish the well-posedness of the solution.

1. Introduction

Magnetohydrodynamics (MHD) is a combined study of hydrodynamic flows and electromagnetic fluid flow through coupling forces and their interactions. These equations consist of highly coupling equations with well-known Maxwell equations and the hydrodynamic equations for fluid flows. Alfven first proposed the field of MHD for single- and multi-fluid flows [1]. More recently, this field has become important because of its utilization and practical applications in geophysics, astrophysics, and many other engineering fields, like cooling metal, MHD propulsion [2], MHD pumps [3], process metallurgy [4], controlled thermonuclear fusion and seawater propulsion [5], electromagnetic casting of metals, MHD power generation, MHD ion propulsion [6,7] etc. Moreover, the hydrodynamical behavior of conducting fluids, e.g., electrolytes, liquid metal cooling in nuclear reactors [8], and plasmas are usually formulated via MHD models [9]. Furthermore, the theoretical and numerical investigations can be further seen in the following articles and are the given references [10,11,12].
Because these equations are not easily solvable for the analytical solution, the alternative way to solve the complex model equations are the numerical solutions where inf-sup or Ladyzhenskaya–Babuska–Brezzi (LBB) conditions need to satisfy for the well-posedness of the models. The complications caused by the well-known inf-sup condition have prompted the introduction of various stabilization techniques intended to circumvent these conditions [13], a Galerkin method, Galerkin least square method, penalty method, L 2 -inner products, different stabilization techniques and many modified methods are abundantly used to find the numerical solutions [14,15,16,17], and the well-posedness of the stationary equations with the inf-sup conditions have been given in [18]. MHD-modeled equations are highly non-linear coupling equations which must have to accept the complicated structure of the solution. Therefore, it does require an attractive, efficient numerical scheme which plays a key role to find the solution of such types of complicated structured models. In MFE literature, the least square methods are utilized for such types of complex and non-linear models to find the approximate solution of higher-order PDEs by defining some norms of functional spaces as Hilbert-type to solve [19,20,21,22,23], singular solutions of flows related to viscoelastic behaviours [24], heat transfer and flows [25], and fluids for which temperature-dependecy or viscosity-dependency are critical [17]. Standard least square methods are efficient, effective, and smart enough to give approximate solutions of the non-linear model equations in the linear algebraic form, which would be symmetric and positive-definite. During the last decade, many methods and theories of least square methods have been discussed in the literature. For further details on such applications, we refer interested readers to [26,27,28] among others.
Similarly, a quasi-least-square scheme (QLSFES) is developed [23] to solve the problems in the L 2 -inner product. Typically, stabilization relies on some form of modification of the discrete continuity equation. For example, in [29] the continuity equation is modified to penalize the stability issue with the addition of a single parameter to make the bilinear form coercive for some specific finite elements. The central point of the scheme is to study the problems in the L 2 -norm to the Oseen type (linearized) forms of given coupled nonlinear problems. It has several benefits; first, only L 2 -inner product with norms is used in these methods which are convenient for the computer programming sense. Second, by using the linearizing way (Oseen type), someone can get iterative methods with symmetric positive-definite coefficient metrics. However, this simplified iterative method is easily convergent in a complex domain of initial guess [23,30,31]. For the MHD-based models, common theories and algorithms for approximation of the solutions of such nonlinear problems are not always applicable because of the functional spaces, i.e., Hilbert space for the velocity and pressure. These are not convenient to utilize the Hilbert space, which will ultimately create some stability issue in the analysis. The QLSFES can be easily applied to circumvent this deficiency of the main conditions (inf-sup or LBB). Thus, the prime intention here is to develop a scheme for a coupled branch of non-linear problems to analyze the existence and convergence of model equations. The least-square MFE schemes are inadequate to address the local convergent properties of coupled non-linear problems. However, this method can be utilized as an effective and sufficient way for the well-posedness of the incompressible MHD models without the inf-sup conditons.
In this contribution, the practicality principle can always be met by transforming the given model equations into a transformed first-order system and forming least-squares functionals that use only L 2 -norms. We focus attention on quasi-least-squares methods for which the discretization step is invoked after the quasi-least-squares functional has been defined. The key point to utilize this setting is that it allows one to point out the variational interpretation of least-squares principles as projections in a Hilbert space with respect to problem-dependent inner products. From this point of view, the principal task in the formulation of the method becomes setting up a least-squares functional that is norm-equivalent ( L 2 -norms) in some Hilbert space. This in turn allows one to work in the variational setting as continuity and coercivity (well-posedness) without the LBB or inf-sup conditions. A quasi-least-squares (QLS) mixed finite element method (MFE) based on the L 2 -inner product is easily convinced to solve a coupled and nonlinear incompressible magnetohydrodynamic (MHD) model. Here, first, these nonlinear models are made linear around a particular state wherein linearized first-order equations represent an algebraic system of equations with symmetric matrices. Secondly, a direct iteration technique is applied to couple up the nonlinearities and obtain a theoretical convergent rate for general initial guess. As far as we know from the existing literature, this method has never been considered to find the solution of the MHD model equations without the LBB/or inf-sup conditions.
From a theoretical viewpoint, such a bona-fide (continuous) QLS finite element methods possess a number of significant and valuable properties, as follows.
  • The weak problems are, in general, coercive.
  • Conforming discretizations lead to stable and, ultimately, optimally accurate methods.
  • The resulting algebraic problems are symmetric and positive definite.
  • Essential boundary conditions may be imposed in a weak sense.
  • Finite element spaces of equal interpolation order, defined with respect to the same triangulation, can be used for all unknowns.
  • Algebraic problems can be solved by using standard and robust iterative methods, such as conjugate gradient methods.
  • Methods can be implemented without any matrix assemblies, even at the element level.
  • Only a single parameter with a single initial guess is sufficient to establish the well-posedness of the solution. In the existing literature this is the first time we apply for this model.
The rest of the work is arranged as follows. In Section 2, we demonstrate linear and non-linear incompressible stationary MHD equations with the proposed scheme (QLSFES). In Section 3, we investigate the well-posedness of the developed QLSFES system for investigation of the initial guess. In Section 4, we discuss convergence of the proposed scheme in the case of the nonsingular exact solution. In Section 5, the theoretical is proven and in conclusion the contribution is illustrated to end this work.

2. Model Introduction

This work illustrates the numerical resolution of the stationary-coupled magnetohydrodynamics system of equations [32,33]. The unknowns for this problem are supposed symbolically as velocity field v f , the magnetic field M B , and fluid pressure p under the connected two-dimensional domain ω .
The non-dimensional stationary MHD equations are as follows:
v f · v f 1 R e η v f + p + S ( M B × × M B ) = f i n ω ,
· v f = 0 i n ω ,
1 R m × ( × M B ) S ( × ( v f × M B ) ) = 0 i n ω ,
· M B = 0 i n ω ,
where R e , signifies the Reynolds number for hydrodynamic, R m the Reynolds number for magnetism, and S = M 2 R e R m , with M > 0 , coupling number respectively. In the literature of the industrial sense, we know that, the parameters are considered always R e 10 5 , R m 10 1 and S 1 . To acquire the values of the known velocity field v f = ( v 1 , v 2 ) and pressure field p and magnetic field M B = ( b 1 , b 2 ) are given in a bounded domain ω . We consider the load function f represents the external forces or inertial terms. Additionally, for the 2D form, the operator can be defined as × M B = b 2 x b 1 y , while the product of given vectors v f = ( v 1 , v 2 ) and M B = ( b 1 , b 2 ) is defined as v f × M B = v 1 b 2 v 2 b 1 .
The system of nonlinear stationary MHD can be solved under the set of boundary conditions within the connected domain ω R 2 on boundary ω is given from the references [33,34,35,36]:
v f = 0 on ω ,
M B · n = 0 on ω ,
× M B × n = 0 on ω ,
where n is a normal vector to the domain ω [17]. Here Equation (5) is representing the viscous nature fluids and is well-known as a no-slip boundary condition, whereas the Equations (6) and (7) are nominated for the perfect conduction wall.
Remark 1.
Alternately, some other boundary constraints are utilized in the literature [12,36,37,38,39] are M B · n = 0 in Equation (6) and × M B × n = 0 in Equation (7) for the unknown M B . We utilize M B × n = 0 and n × ( × M B ) = 0 , which is given for the nonlinear MHD for a single-fluid flow.

3. Quasi-Least-Square MFE

In this section, we discuss the QLES based on the L 2 -inner product with MFE for the stationary MHD models.

3.1. Definitions and Notations

Here, the following Sobolev spaces are introduced first. Let us consider C ( ω ) to be the set of infinitely differentiable continuous functions C 0 ( ω ) = { o , M B C ( ω ) } ; the vector-support functions o , M B are in domain ω . (In the coming sections, we suppose o, M B Oseen type iterative values utilizing for the conversion of the non-linear to linearization of the MHD model equations). Moreover, the Sobolev function spaces in standard form are given as:
W : = [ H 1 ( ω ) ] n × n , I : = H 0 1 ( ω ) n , P : = { q H 1 ( ω ) ; ω q d ω = 0 } , Y : = H 0 1 ( ω ) , N : = [ H 1 ( ω ) ] n .
Let V = ( V i j ) n × n = ( v f 1 , v f 2 , v f n ) be known as a vector function of degree n × n with column matrices v f k , 1 k n . For each n-dimensional matrix v f and vector U , introduce v f = ( v 1 , v 2 , v n ) and × U = ( × v f 1 , × v f 2 , × v f n ) , · U T = ( · v f 1 , · v f 2 , · v f n ) T and similarly B = × M B . For the conciseness and applicability of the QLS method, the first-order system of equations can be now equivalently written as:

3.2. The Stationary First Order MHD Model and QLSFE Scheme Algorithms

An algorithm is proposed to solve the incompressible stationary MHD Model (1)–(4) via the following first order system as following:
ν · U T + ( v f · ) v f + p + S ( M B × × M B ) = f in ω ,
U v f = 0 in ω ,
· v f = 0 in ω ,
ν m × ( B ) S ( × ( v f × M B ) ) = 0 in ω ,
B × M B = 0 in ω ,
· M B = 0 in ω .
In the further formulation for any o , M B I is considered from the Oseen type formulation [23,30] where o and M B are known functions usually chosen from the previous iterative step of the Picard iterations which are always known values for the next iteration and are also supposed regular for the unknown values. We represent the viscosity of the hydrodynamic ν = 1 R e , and the diffusivity of the magnetism ν m = 1 R m . Then the bilinear form of the functions are now introduced for the ( V , v f , p , B , M B ) and ( W , ρ , q , τ , c ) respectively ( W × I × P × Y × N ) as
L ( o , M B ; ( U , v f , p , B , M B ) , ( W , ρ , q , τ , σ ) ) = ( ν · V T + ( o · ) v f + p + S ( M B × × M B ) , ν · W T + ( o · ) ρ + q + S ( M B × × σ ) ) + ( V v f , W ρ ) + γ 2 ( · v f , · ρ ) + ( ν m × B × ) ( v f × M B ) , ν m ( × τ × ( ρ × M B ) ) + ( B × M B , τ × σ ) + η 2 ( · M B , · σ )
where the two positive constants γ and η are to be found in the upcoming sections. Suppose that ( v f , p , M B ) are the solutions of system Equations (1)–(4) which are in [ H 2 ( ω ) ] d × H 1 ( ω ) × H 2 ( ω ) . Let us consider V = v f and B = × M B , and now the solution of unknown variables ( V , v f , p , B , M B ) satisfies the following QLS variational formulation:
L ( v f , M B ; ( V , v f , p , B , M B ) , ( W , ρ , q , τ , σ ) ) = ( f , ν · W T + v f · ρ + q + S ( M B × × σ ) ) ( W , ρ , q , τ , σ ) W × I × P × Y × N .

3.3. Variational Form of the MHD

For the mathematical solutions and findings, few symbols utilized in functional spaces are introduced. Suppose that W m , p ( G ) is noted particularly written as H m ( ω ) with the norm · m and seminorm | · | m [40,41] for p = 2 and m = 0 then W 2 0 ( ω ) . We see that the inner product is always given by ( · , · ) , · = · 0 for the space L 2 ( ω ) respectively, and the general norm for this space can be written as L p ( ω ) norm, which is equivalent to · L p , with some notable cases of Hilbert space L 2 ( ω ) = · L 2 ( ω ) and the Banach space L ( ω ) = · . Here G is ignored for G = ω without any confusion for the notational simplicity.
However, with the above setting, the corresponding variational formulation, we introduce FE spaces ( W h × I h × P h × Y h × N h ) ( W × I × P × Y × N ) on the triangular elements T h where T h is a member of FE triangulations of the discretized domain ω and script h represents the largest element in the meshed domain or the triangles. By considering the QLS illustration, we can write:
QLSFES. Find ( V h , v f h , p h , B h , M B h ) ( W h × I h × P h × Y h × N h ) such that
L v f h , M B h ; ( V h , v f h , p h , B h , M B h ) , ( W h , ρ h , q h , τ h , σ h ) = f , ν · W h + v f h · ρ h + q h + S ( M B h × × σ h ) ( W h , ρ h , q h , τ h , σ h ) ( W h × I h × P h × Y h × N h ) .
Remark 2.
QLSFES is a simple symmetric form. For practical calculation, we define two parameters γ and η with the incompressibility conditions. These parameters are important to ensure the coercivity and continuity of the scheme. Indeed, it is not possible to ensure the coerciveness of the bilinear form of L ( o , M B : · , · ) for each o , M B I . But, by choosing suitable values for γ and η they became well-defined and positive symmetric in some bounded domain which holds the solution of (1)–(4). In the next sections, we will discuss the judicious guess of the parameters γ and η for the linear form. The ideal feature of this technique is that a solution of a nonlinear model equation can be obtain without any initial guess for the iterative procedure for a large domain. We confirm this feature in forthcoming formulations in our analysis part.

4. QLSFES Convergence and Existence of the Solutions

In this section, we obtain convergence and existence of the solution of proposed scheme with the illustration of four steps.
  • We demonstrate property of the L ( o , M B : · , · ) bilinear function. For any parameter γ and η are introduced in such a way that there exists a bounded function set L ( o , M B : · , · ) that is well defined and continuous as o and M B holds in this bounded domain.
  • In step two, we find out a large domain of bounded functions which holds all solutions of the System (8)–(13). We intend to solve the nonlinear coupled equations by the right choice of γ and η ,and L ( o , M B : · , · ) is well-defined as o and M B are in this large domain of functions.
  • To show the well-posedness of the proposed scheme QLSFES, we establish the nonlinear plan of the scheme in such a way that the solutions under some nominated domain are fixed points. In step one and step two, we summarize that for a specific value of η and γ , the system of the nonlinear model is distinctively executable in this specific bounded set (see detail in Lemmas 3 and 4). Moreover, in Theorem 1, by using the fixed point theory [23], we illustrate briefly the existence of solutions of the scheme.
  • In the fourth step, Theorem 2 is given to illustrate the convergence and existence of the proposed scheme QLSFES.
    Before proceeding to the actual contribution, we intend to understand several existing results which are utilized in the immanent sections. By the theory of embedding and the Poincaré’s inequality, the positive constants e, a 0 , and a 1 always depend on the fact that the domain can be recalled as
a 0 u L 4 2 u L 2 ( ω ) 2 , a 1 u 1 2 u L 2 ( ω ) 2 , u I . v f L 4 2 e v f 1 2 v f I .
M B L 2 ( ω ) 2 a 2 ( × M B L 2 ( ω ) 2 + · M B L 2 ( ω ) 2 ) , M B N M B L 4 2 C M B 1 2
× M B L 2 ( ω ) 2 a 3 M B 1 2
M B L 2 ( ω ) 2 e × M B L 2 ( ω ) 2
e M B L 2 ( ω ) 2
According to Green’s formula and integration by parts [42], we can state
( v f , u ) ( v f , u ) = 0 , v f , u I ( o · u , u ) + 1 2 ( d i v o , u 2 ) = 0 , u , o I ( M B × × M B , u ) ( v f × M B , × M B ) = 0 , u I , M B , M B N o × M B L 2 ( ω ) 2 e o L 2 M B L 2 ( ω ) 2 o , M B ( L 2 ( ω ) n × n ) .
The bilinear form ( · , · ) is well known as an inner product in space L 2 ( ω ) or L 2 ( ω ) n . It is always known to be a fact that q L 2 ( ω ) holds the boundary condition ( q , 1 ) = 0 if ϕ [ H 0 1 ( ω ) ] n such that
d i v · ϕ = q , ϕ 1 C e q L 2 ( ω ) .
Lemma 1.
Let constants K o > 0 , 0 < η < 2 a 0 m i n ( ν , ν m 2 ) and γ > m a x ( K o , 1 ) . There exists two constants β and β such that for each o I , b M N satisfying o L 2 ( ω ) K o and · o L 2 ( ω ) η , M B L π and for each ( W h , u h , q h , τ h , σ h ) ( W , I , P , Y , N ) ,
( i ) β u L 2 ( ω ) 2 + × σ L 2 ( ω ) 2 L o , M B ; ( W , u f , q , τ , σ ) , ( W , u f , q , τ , σ ) , ( i i ) Λ Λ * W L 2 ( ω ) 2 + τ L 2 ( ω ) 2 + q L 2 ( ω ) 2 L o , M B ; ( W , u f , q , τ , σ ) , ( W , u f , q , τ , σ ) ,
where the constants Λ and Λ * are defined as follows
Λ = m i n ( Λ 1 , Λ 2 , Λ 3 , Λ 4 , Λ 5 ) , w h e r e Λ j , 1 j 5 ,
and Λ * = m i n 1 2 ( 1 + Λ ) , 1 2 ( C e ) 2 Λ , 1 4 ν 2 C e 2 ( 1 + Λ ) , a 0 2 a 1 2 C e 2 ( η + K ) 2 a 1 .
Proof. 
From Young’s inequality and the given conditions, we can demonstrate
u f L 2 ( ω ) 2 = ( u f W , u f ) + ν 1 ( ν · W T + ( o · ) u f + q + S M B × × σ , u f ) + 1 2 ν ( · o u f , u f ) + ν 1 ( q , · u f ) + ν 1 [ ( S M B × × σ , u f ) ] 1 2 a 1 2 ϵ 0 ν 2 ν · W T + ( o · ) u f + q + S M B × × σ L 2 ( ω ) 2 + 1 2 ϵ 0 u f W L 2 ( ω ) 2 + γ 2 4 ϵ 1 · u f L 2 ( ω ) 2 + ϵ 1 γ 2 ν 2 q L 2 ( ω ) 2
+ ( 2 ϵ 0 + η 2 ν 2 a 0 ) u f L 2 ( ω ) 2 + S 2 π 2 4 ν 2 a 0 ϵ 0 × σ L 2 ( ω ) 2
and
W L 2 ( ω ) 2 2 u f L 2 ( ω ) 2 + u f W L 2 ( ω ) 2 .
Let us consider ϕ holds (23), and we obtain the following estimation
q L 2 2 = ( ν · W T + ( o · ) u f + q + S M B × × σ , ϕ ) + ν ( W , ϕ ) + ( o · ) u f + S M B × × σ , ϕ [ ( ν · W T + ( o · ) u f + q + S M B × × σ ) L 2 ( ω ) + ν W L 2 ( ω ) + 1 a 0 ( · o L 2 ( ω ) + o L 2 ( ω ) ) u f L 2 ( ω ) + S π × σ L 2 ( ω ) ] ϕ H 1 .
By utilizing Equation (23), we can state
q L 2 2 2 C L 2 ( ω ) 2 [ ν · W T + ( o · ) u f + q + S M B × × σ L 2 ( ω ) 2 + ν 2 W L 2 ( ω ) 2 + ( η + K ) 2 a 0 2 u f L 2 ( ω ) 2 + S 2 π 2 × σ L 2 ( ω ) 2 ] ,
in a similar way,
× σ L 2 ( ω ) 2 = ( × σ , × σ τ ) + ν m 1 ( σ , ν m × τ × ( u f × M B ) ) + ν m 1 ( σ , × ( u f × M B ) )
× σ L 2 ( ω ) 2 4 × σ τ L 2 ( ω ) 2 + 4 ν m a 2 ν m 2 a 2 η ν m × τ × ( u f × M B ) L 2 ( ω ) 2 + 2 π 2 ν m 2 u f L 2 ( ω ) 2 + 1 4 · σ L 2 ( ω ) 2 .
By substituting Equations (30) and (32) in Equation (26), we have
u f L 2 2 1 2 a 1 2 ϵ 0 ν 2 ( ν · W T + ( o · ) u f + q + S M B × × σ ) L 2 2 + 1 2 ϵ 0 u f W L 2 ( ω ) 2 + γ 2 4 ϵ 1 · u f L 2 ( ω ) 2 + ( ϵ 0 + η 2 ν a 0 ) u f L 2 ( ω ) 2 + a 2 S 2 4 ν 2 ϵ 1 ν m × τ × ( u f × M B ) L 2 ( ω ) 2 + 2 ϵ 1 × σ L 2 ( ω ) 2 + ϵ 1 · σ L 2 ( ω ) 2 + S 2 ν m 2 4 ν 2 ϵ 1 × σ τ L 2 ( ω ) 2 + 2 ϵ 1 C L 2 ( ω ) 2 γ 2 ν 2 [ ν · W T + ( o · ) u f + q + S M B × × σ 2 + ν 2 W L 2 ( ω ) 2 + ( η + K ) 2 a 0 2 u f L 2 ( ω ) 2 + S 2 π 2 × σ L 2 ( ω ) 2 ] 1 2 a 1 2 ϵ 0 ν 2 + 2 ϵ 1 C L 2 ( ω ) 2 γ 2 ν 2 ( ν · W T + ( o · ) u f + q + S M B × × σ ) L 2 2 + 4 ϵ 1 C 0 2 γ 2 + 1 2 ϵ 0 u f W L 2 ( ω ) 2 + γ 2 4 ϵ 1 · u f 2 + ϵ 0 + η 2 ν a 0 + 2 ϵ 1 C 0 2 ( η + K o ) 2 γ 2 ν 2 a 0 2 + 4 ϵ 1 C 0 2 γ 2 u f L 2 ( ω ) 2 + 2 ϵ 1 C 0 2 S 2 π 2 γ 2 ν 2 + 2 ϵ 1 × σ L 2 ( ω ) 2 + a 2 S 2 4 ν 2 ϵ 1 ν m × τ × ( u f × M B ) L 2 ( ω ) 2 + ϵ 1 · σ L 2 ( ω ) 2 + S 2 ν m 2 4 ν 2 ϵ 1 × σ τ L 2 ( ω ) 2 .
Considering the constants
ϵ 0 = 1 4 ( 1 η 2 ν a 0 ) , ϵ 1 = a 0 ν ν m 2 ( 2 a 0 ν η ) 16 ( C 0 2 ( η + K ) 2 ν m 2 + 2 C 0 2 ν 2 ν m 2 a 0 2 + 2 a 0 2 C 0 2 S 2 π 4 + 2 a 0 2 ν 2 π 2 ) ,
we can state
u f L 2 2 1 4 a 1 2 ϵ 0 ν 2 + C L 2 ( ω ) 2 ϵ 1 ϵ 0 ν 2 ( ν · W T + ( o · ) u f + q + S M B × × σ ) L 2 2 + 2 C 0 2 ϵ 1 ϵ 0 + 1 4 ϵ 0 2 u f W L 2 ( ω ) 2 + γ 2 8 ϵ 0 ϵ 1 · u f L 2 ( ω ) 2 + ϵ 1 C 0 2 S 2 π 2 4 ϵ 0 ν 2 + 3 ϵ 1 4 ϵ 0 · σ L 2 ( ω ) 2 + a 2 S 2 8 ϵ 0 ϵ 1 ν 2 + 4 a 2 ϵ 1 C 0 2 S 2 π 2 ϵ 0 ν 2 ν m 2 + 4 a 2 ϵ 1 ϵ 0 ν m 2 ν m × τ × ( u f × M B ) L 2 ( ω ) 2 + S 2 ν m 2 8 ϵ 0 ϵ 1 ν 2 + 4 ϵ 1 C 0 2 S 2 π 2 ϵ 0 ν 2 + 4 ϵ 1 ϵ 0 × σ τ L 2 ( ω ) 2 .
It is further estimated that τ L 2 ( ω ) 2 is:
τ L 2 ( ω ) 2 2 × σ τ L 2 ( ω ) 2 + × σ L 2 ( ω ) 2 2 × σ τ L 2 ( ω ) 2 + 4 ν m Λ 0 ν m Λ 0 η × σ τ L 2 ( ω ) 2 + 1 Λ 1 ν m 2 ν m × ( M B × u f ) L 2 ( ω ) 2 .
From the Equations (32) and (34), we understand that the conclusion of the condition (24)(i) holds for the Λ defined by (24) and (25), whereas by (24)( i i ), (28), (30) and (35) also satisfies. □
Remark 3.
Lemma 1 deals three positive constants K o , η, and γ where K o and η are nominated constraints for the gradient and divergence of function o illustrated below. We would be able to indicate that K o is playing a key role to control the side terms in right-hand. By employing the constant K o , we may figure out the bounded function set which restrains all solutions of the given MHD model equations. Otherwise, it could not be easy to handle the divergence-free conditions in the theoretical formulation. Here, we manage the · o L 2 ( ω ) η (divergence-free) and · M B L 2 ( ω ) π . The Λ is a coefficient which is actually free from the constat K o and γ but it depends strongly on η and π. It is worth noting that Λ becomes smaller if η is larger. Hence, Lemma 1 illustrates the way by which the parameter Λ can be formulated. This will be utilized to execute γ in practical usage. The value of γ is a penalty factor which is used to control the boundary condition for the velocity field divergence. It is important to note that the parameter Λ * has acting no such type of usage in the mathematical formulation at all.
Remark 4.
Because a magnetic field is a solenoidal field, it may not be considered as compressible or incompressible [11]. To penalize the magnetic effect as divergence or × of the field there is no coefficient so far discussed in the literature. It means the well-psedness of the MHD is still open for the researchers. This will be a challenging problem for our future work.
In active usage, the well-defined bilinear form L ( o , M B ; · , · ) for all o I is not compulsory and a similar remedy is considered for the M B is applied. However, one can find the approximate numerical values of the unknown functions in the particular domain, which have all the exact solutions. We are suppose to seek this bounded function set. To this end, some functions are supposed in W × I × P × Y × N as:
L ˜ ( η , K o , ρ ) = [ ( W , u f , q , τ , σ ) ; · u f L 2 ( ω ) η , [ u f L 2 ( ω ) 2 + × σ L 2 ( ω ) 2 ] 1 / 2 K o , [ W L 2 2 + τ L 2 ( ω ) 2 + q L 2 ( ω ) 2 ] 1 / 2 K o ρ ] ,
L h ˜ ( η , K o , ρ ) = L ˜ ( η , K , ρ ) [ W h × I h × P h × Y h × N h ] .
The given lemma illustrates that for some positive parameters η , K o , ρ , Equations (8)–(13) holds solutions in the bounded set L ˜ .
Lemma 2.
Consider 0 < η < 2 a 0 m i n ( ν , ν m 2 ) . Suppose that Λ and Λ * holds ((24)i,ii) and K o satisfies
f L 2 ( ω ) m i n ( & # x 0039 B ; , 1 ) K o .
All the possible solutions of (8)–(13) are in the function set L ˜ ( η , K o , Λ * ) .
Proof. 
Let ( V , v f , p , B , M B ) be the only one solution of the given system Equations (8)–(13). For all γ 0 , it can be seen that
L ( v f , M B ; ( V , v f , p , B , M B ) , ( V , v f , p , B , M B ) ) = ( f , ν · V T + v f · v f + p + S M B × × M B ) .
This states:
L ( v f , M B ; ( V , v f , p , B , M B ) , ( V , v f , p , B , M B ) ) f L 2 ( ω ) 2 .
We know the term for the incompressible conditon · v f = 0 and · M B = 0 i n ω , so
v f L 2 2 = ν 1 ( ν · V T + ( v f · ) v f + p + S M B × × M B , v f ) + ( v f V , v f ) ν 1 ( S M B × × M B , v f ) ν 1 ( ν · V T + ( v f · ) v f + p + S M B × × M B , v f ) + ( v f V , v f ) + S ν ( ν m × B × ( v f × M B ) , M B ) + S ν m ν ( × M B B , × M B ) 2 a 1 2 ν 2 ν · V T + ( v f · ) v f + p + S M B × × M B L 2 2 + 2 v f V L 2 ( ω ) 2 + 1 4 v f L 2 ( ω ) 2 + 2 a 2 S 2 ν 2 ν m × B × ( v f × M B ) L 2 ( ω ) 2 + 1 4 × M B L 2 ( ω ) 2 + 2 S 2 ν m 2 ν 2 × M B B L 2 ( ω ) 2 ,
including
× M B L 2 ( ω ) 2 2 × M B B L 2 ( ω ) 2 + 2 a 2 ν m 2 ν m × B × ( v f × M B ) L 2 ( ω ) 2 + 2 a 1 2 S 2 ν m 2 ν · V T + ( v f · ) v f + p + S M B × × M B L 2 2 + 2 ν 2 S 2 ν m 2 v f V L 2 ( ω ) 2 + 1 4 v f L 2 ( ω ) 2 + 1 4 × M B L 2 ( ω ) 2 .
Therefore,
Λ ˜ [ v f L 2 ( ω ) 2 + × M B L 2 ( ω ) 2 ] L ( v f , M B ; ( V , v f , p , B , M B ) , ( V , v f , p , B , M B ) ) f L 2 ( ω ) 2 ,
where Λ ˜ = m i n ( a 1 2 S 2 ν 2 ν m 2 4 ( S 2 ν m 2 + ν 2 ) , S 2 ν m 2 4 ( S 2 ν m 2 + ν 2 ) , ν 2 4 a 2 S 2 , ν 2 4 S 2 ν m 2 ) . Accordingly, it can be easily seen that Λ Λ ˜ . Because v f L 2 ( ω ) 2 K o . By taking γ m a x ( 1 , K o ) and by using Lemma 1, one can see that ( V , v f , p , B , M B ) is in L ˜ ( η , K o , Λ * ) . □
We acquire the approximate solution of the System (8)–(13). Let us consider the nonlinear map based upon the defined spaces W h × I h × P h × Y h × N h into W h × I h × P h × Y h × N h as
L ( W h , u f h , q h , τ h , σ h ) = ( W ^ h , u f ^ h , q ^ h , τ ^ h , σ ^ h )
such that for each ( W ¯ h , u f ¯ h , q ¯ h , τ ¯ h , σ ¯ h ) W h × I h × P h × Y h × N h ,
L ( u f h , σ h ; ( W ^ h , u f ^ h , q ^ h , τ ^ h , σ ^ h ) , ( W ¯ h , u f ¯ h , q ¯ h , τ ¯ h , σ ¯ h ) ) = ( f , ν · W ¯ h + u f h · u f ¯ h + q ¯ h + S M B ¯ h × τ ¯ h ) .
Now, it is certain that the System (44) is linear with respect to ( W ^ h , u f ^ h , q ^ h , τ ^ h , σ ^ h ) . Now the non-linear MHD can be estimated through the following results given below.
Lemma 3.
Suppose the Lemmas 1 and 2 satisfy. Then the bilinear form of the system L from L ˜ h ( η , K o , Λ * ) to W h × I h × P h × Y h × N h is distinctly defined.
The key point is that the Lemma 3 is the straightforward result of Lemma 1.
Lemma 4.
Let us consider that the results of Lemmas 1 and 2 satisfies all the conditions for γ and it satisfies
γ f L 2 ( ω ) max ( 1 , 1 / η , 1 / Λ ) .
Now the bilinear operator L shows L ˜ h ( η , K o , Λ * ) to itself.
Proof. 
Hence from Equation (44)
L ( u f h , σ h ; ( W ^ h , u f ^ h , q ^ h , τ ^ h , σ ^ h ) , ( W ^ h , u f ^ h , q ^ h , τ ^ h , σ ^ h ) ) f L 2 ( ω ) 2 ,
because we have u f h L 2 ( ω ) K o and · u f h L 2 ( ω ) η so Equation (24) tends to
  • u f ^ h L 2 2 + × σ ^ h L 2 ( ω ) 2 1 Λ f L 2 ( ω ) 2 K o 2 ;
  • W ^ h L 2 2 + τ ^ h L 2 ( ω ) 2 + q ^ h L 2 ( ω ) 2 1 Λ Λ * f L 2 ( ω ) 2 K o 2 Λ * .
Furthermore,
· u f ^ h L 2 γ 1 f L 2 ( ω ) η .
Therefore, ( W ^ h , u f ^ h , q ^ h , τ ^ h , σ ^ h ) L ˜ h ( η , K o , Λ * ) . This proves completed. □
Theorem 1.
Let 0 < η < 2 a 0 m i n ( ν , 1 ν m ) , K holds (38), γ holds (45), Λ and Λ * holds ((24)i,ii). Therefore the scheme QLSFES satisfies and holds a minimum one solution in bilinear form L ˜ h ( η , K o , Λ * ) . However, all the illustrated unknown solutions of the nonlinear model Equation (16) are in L ˜ h ( η , K , Λ * ) .
Proof. 
Lemma 4 states that, the operator L relates the bounded domain L ˜ h ( η , K , Λ * ) into itself under the specific conditions of Theorem 1. It satisfies the fixed point Browners theory that the nonlinear System (16) has minimum one solution in L ˜ h ( η , K o , Λ * ) . Indeed it is true that all the solutions of the unknowns of model nonlinear Equation (16) are in L ˜ h ( η , K o , Λ * ) . Hence the proof of Theorem 1 is completed. □
Moreover, the solutions of QLSFES convergence is demonstrated below as:
Theorem 2.
Suppose that Theorem 1 holds and ( V h , v f h , p h , B h , M B h ) is a unique sequence of the solutions of the proposed scheme as h 0 . Then the sequence of solution ( V h , v f h , p h , B h , M B h ) can be further subdivided into many subsequences which will converge weakly to the different solutions of the first-order MHD (8)–(13). Particularly, components of the MHD v f h , M B h are weakly convergent into the subsequences and are strongly convergent to the similar limit components in [ H s ( ω ) d ] × H s ( ω ) for specific condition 1 > s 0 .
Before proceeding to the prove of the Theorem 1 we entail some lemmas for the understanding, which are given as
Lemma 5.
We suppose the inequality 0 < η < 2 a 0 ν . For a specific Oseen-type function o ¯ [ H 1 ( ω ) ] d holds the relation · o ¯ L 2 ( ω ) η and a particular load function f ¯ [ L 2 ( ω ) ] d . At present the problem with boundary values can be stated as
ν η v f * + o ¯ · v f * + p = f ¯ , in ω ;
· v f * = 0 , ( p * , 1 ) = 0 , v f * = 0 , on ω ;
which holds a unique solution ( v f * , p * ) in H 2 ( ω ) × H 1 ( ω ) .
The reader can see [23] appendix 1. Hence, the embedding theory between Sobolev spaces and some results reported in [23,40] are directly utilized here as in the following lemmas given below.
Lemma 6.
Assume that R is a Hilbert space in which F might be a bounded function set within R , i.e., there is a domain-dependent constant K o > 0 in such a way that f R K for every f F . If T is a function set which is weakly compact in the given Hilbert space R [23]:
lim n f n , g = f , g , g R .
Here · , · denotes the inner product of two functions in space R .
Lemma 7.
Let T be a set of bounded functions in Hilbert space H 1 ( ω ) , such that there exists K o such that v H 1 ( ω ) K o for each v T . Then the T set is compact strongly in space H s ( ω ) for every interval 0 s < 1 , i.e., the { f n } n = 1 , which is strongly convergent in the Hilbert space H s ( ω ) as n , which may be deduced from T.
Lemma 8.
Suppose a positive constant C exists and holds for interval 0 s 1 and every v H s ( ω )
e v H s ( ω ) v L q ( ω ) , 1 q 2 n n 2 s .
Hence Theorem 2 might be resumed as:
Proof. 
It comes from Theorem 1 that the solutions of QLSFES ( V h , v f h , p h , B h , M B h ) are bounded [23]. Lemmas 6 and 7 concludes that the results are divided into many other subsequences which are weakly convergent in the L 2 × H 1 ( ω ) × L 2 × H 1 ( ω ) . Moreover, the weak convergence of the subsequence of the proposed scheme, might still represent it by ( V h , v f h , p h , B h , M B h ) and its weak limitation by ( V ¯ , v f ¯ , p ¯ , B ¯ , M B ¯ ) in [ L 2 ( ω ) × H 1 ( ω ) × L 2 ( ω ) × H 1 ( ω ) ]. However, from Lemma 7, someone may know that components ( v f h , M B h ) are strongly convergent to ( v f ¯ , M B ¯ ) in [ H s ( ω ) × H s ( ω ) ] for condition 0 s < 1 .
We are supposed to justify that ( V ¯ , v f ¯ , p ¯ , B ¯ , M B ¯ ) is the first-order linear system Solution (8)–(13). As a result, ( v f ¯ , p ¯ , M B ¯ ) is a unique solution of the MHD. To this end, tentative functions ( v f * , p * ) are introduced as:
v f ¯ · v f * η v f * + p * = f S ( M B ¯ × × M B ¯ ) , i n ω ;
· v f * = 0 , i n ω ;
v f * = 0 , o n ω ;
and similarly for M B *
1 μ σ × ( × M B * ) × ( v f ¯ × M B * ) = 0 , i n ω ;
· M B * = 0 , i n ω ;
M B * = 0 , o n ω ;
Right now, Systems (49a) and (50a) represent two different linear systems, which are independent Lemma 5 concludes that for the ( v f * , p * ) [ H 2 ( ω ) ] d × [ H 1 ( ω ) ] d these two free system of equations are distinctly executable. Let V * = v f * and B * = × M B * . It is true that ( V ¯ , v f ¯ , p ¯ , B ¯ , M B ¯ ) is a unique solution for the Equations (8)–(13) if the relation holds as ( V ¯ , v f ¯ , p ¯ , B ¯ , M B ¯ ) = ( V * , v f * , p * , B * , M B * ) . We intend to justify this in three steps.
In Step 1, we execute
ν · ( V ¯ V * ) T + ( v f ¯ · ) ( v f ¯ v f * ) + ( p ¯ p * ) L 2 ( ω ) = 0 .
In Step 2, the proof is
( V ¯ V * ) ( v f ¯ v f * ) = 0 , i n ω ;
· ( v f ¯ v f * ) = 0 , i n ω .
In Step 3, from (51), (52a) and (52b) it is satisfied that
ν η ( v f ¯ v f * ) + ( v f ¯ · ) ( v f ¯ v f * ) + ( p ¯ p * ) = 0 , i n ω ;
· ( v f ¯ v f * ) = 0 , i n ω .
v f ¯ v f * = 0 , o n ω .
Equation (53a) represents that ( V ¯ , v f ¯ , p ¯ ) = ( V * , v f * , p * ) because the Equation (53a) holds a unique but trivial solution. Similarly, one can prove ( B ¯ , M B ¯ ) = ( B * . M B * ) . □
The proof of (51) can be written and it is obvious that (51) is identical to
( ν · ( V ¯ V * ) T + ( v f ¯ · ) ( v f ¯ v f * ) + ( p ¯ p * ) , φ ) = 0 , φ [ C 0 ] d .
For every φ [ C 0 ( ω ) ] n , one might have a test functions for the variables velocity and pressure as (v,q).
ν η v + ( v f ¯ · ) v + q = φ , i n ω ,
· v = 0 , i n ω ,
v = 0 , o n ω .
From Lemma 5, we conclude that ( v , q ) [ H 2 ( ω ) ] d × H 1 ( ω ) . By assuming V = v , one can have
( ν · ( V ¯ V * ) T + ( v f ¯ · ) ( v f ¯ v f * ) + ( p ¯ p * ) , φ ) = ν · ( V ¯ V h ) T + ( v f ¯ · ) ( v f ¯ v f h ) + ( p ¯ p h ) + S M B ¯ × × M B ¯ M B h , φ + ( ( v f ¯ v f h · ) v f h , φ ) + ( ν · ( V h V * ) T + ( v f h · ) v f h ( v f ¯ · ) v f * + ( p h p * ) + S M B ¯ × × ( M B ¯ M B h ) , ( ( v f ¯ v f h ) · ) v ) + [ ( ν · ( V h V * ) T + ( v f h · ) v f h ( v f ¯ · ) v f * + ( p h p * ) + S M B ¯ × × ( M B ¯ M B h ) , ν · V T + ( v f h · ) v + q + ( V h v f h , V v ) + γ 2 ( · ( v f h v f * ) , · v ) ) ] .
The weak convergence property of ( V h , v f h , p h , τ h , M B h ) is the unique solution of the MHD System (8)–(13) for the nominated positive parameter γ . This might be further seen equivalently as
( ν · ( V ¯ V h ) T + ( v f ¯ · ) ( v f ¯ v f h ) + ( p ¯ p h ) + S M B ¯ × × ( M B ¯ M B h ) , φ ) ( v f ¯ · ) ( v f ¯ v f h ) + S M B ¯ × × ( M B ¯ M B h ) , φ + ν · ( V ¯ V h , φ ) + ( p p h , d i v · φ ) 0 , a s h 0 .
Therefore Lemma 8 and the ( v f h , M B h ) concludes that
v f ¯ v f h L 4 C v f ¯ v f h H d / 4 0 , a s h 0 ,
together with
( ( ( v f ¯ v f h ) · ) v f h , φ ) + ν · ( V h V * ) T + ( v f h · ) v f h ( v f ¯ · ) v f * + ( p h p * ) + M B ¯ × × ( M B ¯ M B h ) , ( ( v f ¯ v f h ) · v ) C v f ¯ v f h L 4 [ v f h 0 φ L 4 + f 0 v L 4 ] .
From Equations (8)–(13), it can be further estimated as
( ν · ( V h V * ) T + ( v f h · ) v f h ( v f ¯ · ) v f * + ( p h p * ) + M B ¯ × × ( M B ¯ M B h ) , ν · V T + ( v f h · ) v + q ) + ( V h v f h , V v ) + ( · ( v f h v f * ) , · v ) = inf ( V h , v h , q h ) ( W h × X h × Q h ) [ ( ν · ( V h V * ) T + ( v f h · ) v f h ( v f ¯ · ) v f * + ( p h p * ) + M B ¯ × × ( M B ¯ M B h ) , ν · ( V V h ) T + ( v f h · ) ( v v h ) + ( q q h ) ) + ( V h v f h , V V h ( v v h ) ) + γ 2 ( ( v f h v f * ) , · ( v v h ) ) ] 0 , h 0 .
If we replace Equations (57)–(59) with (56), it yields (54). Now we may proceed the proof of the second step (52a) and (52b). With the help of weak convergence solution sequence of the FE space approximations, it can be further seen for each ( W , u f ) [ C 0 ( ω ) ] d × d × [ C 0 ( ω ) ] d , and we have
( V ¯ V * ( v f ¯ v f * ) , W u f ) + γ 2 ( · ( v f ¯ v f * ) , · u f ) = ( ν · ( V ¯ V h ) T + ( v f ¯ · ) ( v f ¯ v f h ) + ( p ¯ p h ) + M B ¯ × × ( M B ¯ M B h ) , ν · W T + ( v f ¯ · ) u f ) + ( V ¯ V h ( v f ¯ v f h ) , W u f ) + γ 2 ( · ( v f ¯ v f h ) , · u f ) + ( ( ( v f ¯ v f h ) · ) v f h , ν · W T + ( v f ¯ · ) u f ) + ( ν · ( V h V * ) T + ( v f h · ) v f h ( v f ¯ · ) v f * + ( p h p * ) + M B ¯ × × ( M B ¯ M B h ) , ( ( v f ¯ v f h ) · ) u f ) + inf ( W h , u f h ) ( W h × X h ) [ ( ν · ( V h V * ) T + ( v f h · ) v f h ( v f ¯ · ) v f * + ( p h p * ) + M B ¯ × × ( M B ¯ M B h ) , ν · ( W W h ) T + ( v f h · ) ( u f u f h ) ) + ( V h V * ( v f h u f * ) , W W h ( u f u f h ) ) + γ 2 · ( v f h v f * ) , · ( u f u f h ) ] 0 , as h 0 .
Now, this shows a complete proof of the Theorem 2 which is indeed a solution of QLSFES in general cases. Moreover, close to the approximate solutions of a singular solution, someone can only execute approximate solutions of weak convergence of the subsequence. However, in the next section, the proof of the strong convergence of the MHD holding with the uniform convergence rate in non-singular solutions are briefly illustrated.

5. Convergent Rate of Non-Singular Solution of QLSFES

In this section, the non-singular solutions and convergence of the system Equations (8)–(13) are demonstrated. A solution of (8)–(13) is considered as the non-singular solution in some special cases i.e., if this solution does not depend on any other solution and the first-order differential approximation of the system is non-singular at this solution [23]. We consider a linear system ( F 0 * , F 1 * , f * ) [ H 1 ( ω ) ] d × [ L 2 ( ω ) ] n × n × [ L 2 ( ω ) ] , which states that if
ν · W T + ( v f · ) u f + q + S ( M B × × σ ) = F 0 * in ω ,
W u f = F 1 * in ω ,
· u f = f * in ω ,
u f = 0 on ω , ( q , 1 ) = 0
ν m × τ S ( × ( v f × σ ) ) = 0 in ω ,
τ × σ = 0 in ω ,
σ = 0 on ω ,
· σ = 0 in ω ,
it holds a unique solution ( W , u f , q , τ , σ ) [ L 2 ( ω ) ] n × n × [ H 0 1 ( ω ) ] n × [ L 2 ( ω ) ] × [ H 0 1 ( ω ) ] . For all of the above, the constant e holds the relation
W T 0 + u f 1 + q 0 + σ 0 + τ 0 C F 0 * H 1 F 1 * 0 f * 0 in ω .
We may assume the Stokes equation for further analysis as
ν η u f + q = f in ω ,
u f = 0 on ω , ( q , 1 ) = 0 ,
· u f = f * in ω ,
and
ν η w = g in ω ,
w = 0 on ω .
Here, H 2 + t is a regular solution for any t ( 0 , 1 ] so the solution ( u f , q , w ) [ H 2 + t ( ω ) ] n × [ H 1 + t ( ω ) ] × [ H 2 + t ( ω ) ] ,
u f H 2 + s + q H 1 + t f H t ,
and
w H 2 + t C g H t .
Furthermore, we assume that the MFE spaces W h × I h × P h × Y h × N h satisfies the approximation properties such that there exists an approximation order r of the MFE space which is defined as ( r 1 ) and C such that
inf ( W h W h ) [ W W h 0 + h W W h 1 ] C h r + 1 W r + 1 W [ H r + 1 ] n × n ; inf ( u h I h ) [ u u h 0 + h u u h 1 ] C h r + 1 u r + 1 u [ H r + 1 ] n × n , u I [ H r + 1 ] n ; inf ( q h P h ) [ q q h 0 + h q q h 1 ] C h r + 1 q r + 1 q [ H r + 1 ] n × n , q P [ H r + 1 ] n ; inf ( τ h N h ) [ τ τ h 0 + h τ τ 1 ] C h r + 1 τ r + 1 τ [ H r + 1 ] d × n , τ N [ H r + 1 ] n . inf ( σ h Y h ) [ σ σ h 0 + h σ σ h 1 ] C h r + 1 σ r + 1 σ [ H r + 1 ] n × n , σ Y [ H r + 1 ] n ;
here, r is the order of approximation for spaces defined.
Theorem 3.
Let us suppose the conditions of Theorem 2 holds. Assume that ( V h , v f h , p h , B h , M B h ) is the unique approximate solution which weakly converges to ( V , v f , p , B , M B ) (8)–(13). Suppose that ( V , v f , p , B , M B ) is a non-singular solution of the given model, and then the sequence of solution ( V h , v f h , p h , B h , M B h ) converges strongly in [ L 2 ( ω ) ] n × [ L 2 ( ω ) ] × [ L 2 ( ω ) ] n × n × [ H 1 ( ω ) ] as h 0 . Thus the estimation results priori satisfies the following relation:
V V h 0 + v f v f h 1 + p p h 0 + B B h 0 + M B M B h 1 C h r
Proof. 
Suppose ( V h , v f h , p h , B h , M B h ) is a solution sequence for a discrete scheme, which is weakly convergent to a unique continuous solution (non-singular) ( V , v f , p , B , M B ) of its lowest order Model (8)–(13) and belongs to [ L 2 ( ω ) ] n × [ L 2 ( ω ) ] × [ L 2 ( ω ) ] n × n × [ H 1 ( ω ) ] as h 0 .
Furthermore, we can represent × as:
× ( v f × M B ) = ( · M B ) v f + ( M B · ) v f ( · v f ) M B + ( v f · ) M B = ( M B · ) v f ( v f · ) M B .
Then, we note here
( a ) ν · ( V V h ) T + ( v f · ) ( v f v f h ) + ( ( v f v f h ) · ) v f + ( p p h ) + S ( M B × × ( M B M B h ) ) = ν · ( V V h ) T + ( v f · ) v f ( v f h · ) v f h + ( p p h ) + S ( M B × × ( M B M B h ) + ( ( v f v f h ) · ) ( v f v f h ) , i n ω ; ( b ) ν m × ( B B h ) × ( ( v f v f h ) × ( M B M B h ) ) = ν m × ( B B h ) + { ( M B · ) ( v f v f h ) + ( M B M B h ) · ) v f } + { ( v f · ) ( M B M B h ) + ( v f v f h ) · ) M B } = ν m × ( B B h ) + { ( M B · ) M B ( M B h · ) v f h + ( M B M B h ) · ) ( v f v f h ) } + { ( v f · ) v f ( v f h · ) M B h + ( v f v f h ) · ) ( M B M B h ) } + ( ( v f v f h ) · ) ( v f v f h ) ( ( M B M B h ) · ) ( M B M B h )
Equation (69) gives the inequality as
V V h 0 + v f v f h 1 + p p h 0 + B B h 0 + M B M B h 1 C { ν · ( V V h ) T + ( v f · ) v f ( v f h · ) v f h + ( p p h ) + S M B × × ( M B M B h ) 0 2 } + V V h ( v f v f h ) 0 2 + γ 2 · ( v f v f h ) L 2 ( ω ) 2 + ν m × ( B B h ) { ( M B · ) M B ( M B h · ) v f h ( M B M B h ) · ) ( v f v f h ) } { ( v f · ) v f ( v f h · ) M B h ( v f v f h ) · ) ( M B M B h ) } 1 2 + B B h × ( M B M B h ) 0 2 + η 2 · ( M B M B h ) 0 2 + ( v f v f h ) · ) ( v f v f h ) H 1 2 + ( M B M B h ) · ) ( v f v f h ) H 1 + ( v f v f h ) · ) ( M B M B h ) H 1 .
For the error estimations, we may bound right-hand terms of Equation (79) as
( v f v f h ) · ) ( v f v f h ) H 1 2 + ( M B M B h ) · ) ( v f v f h ) H 1 + ( v f v f h ) · ) ( M B M B h ) H 1 C ( v f v f h ) L 4 2 [ ( v f v f h ) + ( M B M B h ) L 2 ( ω ) 2 ] + ( M B M B h ) · ) ( v f v f h ) H 1 .
It follows from (8)–(13) and (16) we can get
ν · ( V V h ) T + ( v f · ) v f ( v f h · ) v f h + ( p p h ) + S M B × × ( M B M B h ) 0 2 } + V v f ( V h v f h ) 0 2 + · ( v f v f h ) L 2 ( ω ) 2 + ν m × ( B B h ) { ( M B · ) M B ( M B h · ) v f h ( M B M B h ) · ) ( v f v f h ) } { ( v f · ) v f ( v f h · ) M B h ( v f v f h ) · ) ( M B M B h ) } 1 2 = inf ( W h , u f h , q h , τ h , σ h ) ( W h × X h × P h × Y h × N h ) [ ( ν · ( V V h ) T + ( v f · ) v f ( v f h · ) v f h + ( p p h ) + S M B × × ( M B M B h ) , ν · ( V W h ) T + ( v f h · ( v f u f h ) ) + ( p q h ) + S M B × × ( M B σ h ) ) ) + ( V v f ( V h v f h ) , V W h ( v f u f h ) ) + ( · ( v f v f h ) , · ( v f u f h ) ) + ( ν m × ( B B h ) { ( M B · ) M B ( M B h · ) v f h ( M B M B h ) · ) ( v f v f h ) } { ( v f · ) v f ( v f h · ) M B h ( v f v f h ) · ) ( M B M B h ) } ) , ν m × ( B τ h ) { ( M B · ) M B ( M B h · ) v f h ( M B σ h ) · ) ( v f v f h ) } { ( v f · ) v f ( v f h · ) M B h ( v f u f h ) · ) ( M B σ h ) } ) + ( · ( M B M B h ) , · ( M B σ h ) ) ) ] ν · ( V V h ) T + ( v f · ) v f ( v f h · ) v f h + ( p p h ) + S M B × × ( M B M B h ) , ( ( v f v f h ) · ) v f ) + ( ν m × ( B B h ) { ( M B · ) M B ( M B h · ) v f h ( M B M B h ) · ) ( v f v f h ) } { ( v f · ) v f ( v f h · ) M B h ( v f v f h ) · ) ( M B M B h ) ) , ( ( v f v f h ) · ) ( v f v f h ) + ( M B M B h ) · ( M B M B h ) .
To bind the right-hand side final terms of Equation (81), we can define the alternative equation as
ν η v + ( v f · ) v + q + S M B × × ( σ ) = ( ( v f v f h ) · ) v f
· v = 0 , in ω ; ( q , 1 ) = 0 ;
ν m × × ( σ ) ( M B · ) σ ( v f · ) v = ( ( v f v f h ) · ) v f + ( ( M B M B h ) · ) M B
Hence this System (82)–(84) possesses one unique solution ( v , q , σ ) in [ H 2 ( ω ) ] d × H 1 ( ω ) × H 2 ( ω ) . Let us consider V = v and τ = × σ . Then by auxiliary ( V , v , q , τ , σ ) we have the following relation:
( ν · ( V V h ) T + ( v f · ) v f ( v f h · ) v f h + ( p p h ) + S ( M B M B h ) × × M B , ( ( v f v f h ) · ) v f ) + ( ν m × ( B B h ) ( M B · ) M B ( M B h · ) v f h ( v f · ) v f ( v f h · ) M B h , ( ( v f v f h ) · ) v f + ( ( M B M B h ) · ) M B ) = inf ( V h , v h , q h , τ h , σ h ) ( W h × X h × P h × Y h × N h ) [ ( ν · ( V V h ) T + ( v f · ) v f ( v f h · ) v f h + ( p p h ) + S ( M B M B h ) × × M B , ν · ( V V h ) T + ( v f h · ( v v h ) ) + ( p q h ) + S ( σ σ h ) × × M B ) ) + ( V V h ( v f v f h ) , V V h ( v v h ) ) + γ 2 ( · ( v f v f h ) , · ( v v h ) ) ( B B h + × ( M B M B h ) , τ τ h + × ( σ σ h ) ) + ν m × ( B B h ) ( M B · ) M B ( M B h · ) v f h ( v f · ) v f ( v f h · ) M B h , ν m × ( τ τ h ) ( M B h · ) ( σ σ h ) ( v f v f h · ) ( v v h ) ] ( ν · ( V V h ) T + ( v f · ) v f ( v f h · ) v f h + ( p p h ) + S ( M B M B h ) × × M B , ( ( v f v f h ) · ) v ) + ( ν m × ( τ τ h ) ( M B · ) M B ( M B h · ) v f h ( v f · ) v f ( v f h · ) M B h , ( ( v f v f h ) · ) v + ( ( M B M B h ) · ) σ ) [ ν · ( V V h ) T + ( v f · ) v f ( v f h · ) v f h + ( p p h ) + S ( M B M B h ) × × M B 0 2 + V V h ( v f v f h ) 0 2 + γ 2 · ( v f v f h ) 0 2 + B B h + × ( M B M B h ) 0 2 + ν m × ( τ τ h ) ( M B · ) M B ( M B h · ) v f h ( v f · ) v f ( v f h · ) M B h 0 2 ] + C h 2 r v f H r + 1 2 + v f H r + 1 2 + p H r + 1 2 + M B H r + 1 2 + h 2 s v H s + 2 2 + q H s + 1 2 + v f v f h L 4 2 ( v H 2 2 + τ H 2 2 )
By Substituting (75), (76), and (85) into (81), we get
ν · ( V V h ) T + ( v f · ) v f ( v f h · ) v f h + ( p p h ) + S ( M B M B h ) × × ( M B ) 0 2 + V v f ( V h v f h ) 0 2 + · ( v f v f h ) L 2 ( ω ) 2 + ν m × ( B B h ) { ( M B · ) M B ( M B h · ) v f h ( M B M B h ) · ) ( v f v f h ) } ( v f · ) v f ( v f h · ) M B h ( v f v f h ) · ( M B M B h ) } 1 2 C { h 2 r v f H r + 2 2 + p H r + 1 2 + M B H r + 2 2 }
+ ( h 2 t + v f v f h L 4 2 ) ( v f v f h ) · v f H s 2 + ( v f v f h ) · M B H t 2 } .
Having the Equations (80) and (86) into (79), and the assignment of t = 0 leads to the following relation:
V V h 0 + v f v f h 1 + p p h 0 + B B h 0 + M B M B h 1 C { h r v f H r + 2 2 + p H r + 1 2 + M B H r + 2 2
+ ( 1 + v f v f h L 4 2 ) v f v f h L 4 2 ) v f H 2 + M B H 2 }
where 1 > t > 0
V V h 0 + v f v f h 1 + p p h 0 + B B h 0 + M B M B h 1 C { h r v f H r + 2 2 + p H r + 1 2 + M B H r + 2 2
+ ( h s + v f v f h L 4 2 ) v f v f h H 1 2 ) v f H 3 + M B H 3 } .
Consequently, if v f v f h L 4 is a relation which converges to 0 as h 0 , Equation (86) leads us to
V V h 0 + v f v f h 1 + p p h 0 + B B h 0 + M B M B h 1 0 , as h 0 .
For t = 0 , Equation (88) leads to (78) for interval 1 > t > 0 . This is a complete proof of the theorem. □

Limitation and Future Work

This scheme can be applied in many problems, such as
  • an error analysis of uasi-least squares finite element method of velocity-pressure-magnetic field formulation for MHD problem [22,43];
  • quasi-least-square method to solve the MHD with four unknowns, i.e., velocity of the fluid , velocity of the magnetic field, pressure of the fluid and pressure for the magnetic field;
  • quasi-least square method for the Maxwell equations [44]; and
  • quasi-least square method for the second order MHD model equations [45].

6. Numerical Examples

In this section, some numerical tests by using the public domain software Freefem++ [46] are performed to confirm the correctness of the theoretical prediction. The well-known Taylor–Hood elements ( P 2 P 2 P 1 ) are utilized in computational code under UMFPACK solver to approximate the velocity, magnetic field, and pressure respectively.
Consider the following two-dimensional problem by setting domain ω = [ 0 , 1 ] × [ 0 , 1 ] with a homogeneous Dirichlet boundary conditions on all the boundaries. This domain has been discretized or meshed with 2592 linear triangular elements, 1390 vertices, and 144 boundary edges. The components of velocity v f and magnetic field B M that correspond to ( v 1 , v 2 ) and ( B 1 , B 2 ) are supposed. The computational code has been tested for five iterations by assuming physical parameters as R e = 1 , R m = 200 , S = 1 , ν = 1 , iterative tolerance ϵ = 1.0 × e 6 . The initial guesses for all variables are set to be zero. The initial condition, Dirichlet boundary condition, and the source term f are posed in such a way that it follows the exact solution [47] of the model.
v 1 = x 2 ( x 1 ) 2 y ( y 1 ) ( 2 y 1 ) ; v 2 = x ( x 1 ) ( 2 x 1 ) y 2 ( y 1 ) 2 ; P = ( 2 x 1 ) ( 2 y 1 ) ; B 1 = s i n ( π x ) c o s ( π y ) ; B 2 = s i n ( π y ) c o s ( π x ) .

6.1. Example 1

This example is demonstrated to illustrate the convergence rate concerning the L 2 -error norm with respect to first and second order errors. In Table 1, the approximate errors and convergence rate of the hydrodynamic variable and magnetic field variables are displayed and observed that v f v f h L 2 , p p h L 2 , M B M B h L 2 converges to zero with the varying spacing h = 1 / 4 , 1 / 8 , 1 / 16 , 1 / 32 , a n d 1 / 64 , for the well-known MINI elements P 1 b P 1 P 1 .
Mathematically, L 2 -error estimation between the exact solution v f and approximate solution v f h is designated as
v f v f h L 2 = ω v f v f h 2 d x d y 1 2 .
In Table 2, the well-known Taylor–Hood elements ( P 2 P 2 P 1 ) with the varying spacing h = 1 / 4 , 1 / 9 , 1 / 16 , 1 / 25 , and 1/36, are utilized in computational code to approximate the velocity, magnetic field, and pressure respectively. Moreover, the computational code for the convergence rate or convergence order R is delivered by the following formulations:
E 1 = C h 1 R
E 2 = C h 2 R
E 1 E 2 = C h 1 h 2 R
log E 1 E 2 = R log h 1 h 2
R = log E 1 E 2 log h 1 h 2 .

6.2. Example 2

In this example, the software figures are directly pasted to display the geometry: In Figure 1, the domain with its mesh generation is exemplified. In Figure 2, we see flow lines known as a velocity contour. Figure 3, encapsulates the pressure lines similarly, and magnetic field lines are instantiated in Figure 4. All these figure are decorated to verify for the proposed scheme.

7. Conclusions

A QLS mixed finite element method (MFE) based on L 2 -inner product is discussed to solve a coupled and nonlinear incompressible magnetohydrodynamic (MHD) model. The model equations are highly coupled and nonlinear. Some functional spaces were not convenient to transform model equations into the weak formulation due to which the stability issue occurs. Therefore the L 2 -inner product is utilized as a stabilization technique to circumvent this deficiency. First, we utilized the Oseen-type technique for the conversion of nonlinearity into linearity. Secondly, a direct iteration technique is applied to overcome the nonlinearities and obtained a theoretical convergence rate for the general initial guesswork. Theoretical results show that only a single parameter with a single initial guess is sufficient to establish the well-posedness of the solution.

Author Contributions

S.H. has formulated the main work and performed analysis, S.u.R. has performed the composing and sequencing. S.A. has done the numerical formulation and computational work for the error and convergence rate. M.A.A. supervised all the work thoroughly. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors thank the editor and referees for their valuable comments and suggestions which helped us to improve the results of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Alfvén, H. Existence of electromagnetic-hydrodynamic waves. Nature 1942, 150, 405. [Google Scholar] [CrossRef]
  2. Mitchell, D.L.; Gubser, D.U. Magnetohydrodynamic ship propulsion with superconducting magnets. J. Supercond. 1988, 1, 349–364. [Google Scholar] [CrossRef]
  3. Al-Habahbeh, O.; Al-Saqqa, M.; Safi, M.; Khater, T.A. Review of magnetohydrodynamic pump applications. Alex. Eng. J. 2016, 55, 1347–1358. [Google Scholar] [CrossRef] [Green Version]
  4. Davidson, P.A. Magnetohydrodynamics in material processing. Annu. Rev. Fluid Mech. 1999, 31, 273–300. [Google Scholar] [CrossRef]
  5. Yuksel, G.; Ingram, R. Numerical Analysis of a Finite Element, Crank-Nicolson Discretization for MHD Flow at Small Magnetic Reynolds Number; Tech. Report; University of Pittsburgh: Pittsburgh, PA, USA, 2011. [Google Scholar]
  6. Shercliff, J.A. Steady motion of conducting fluids in pipes under transverse magnetic fields. Proc. Camb. Philos. Soc. 1953, 49, 136–144. [Google Scholar] [CrossRef]
  7. Gerbeau, J.-F.; Bris, C.L.; Lelievre, T. Mathematical Methods for the Magnetohydrodynamics of Liquid Metals, Numerical Mathematics and Scientific Computation; Oxford University Press: New York, NY, USA, 2006. [Google Scholar]
  8. Barleon, L.; Casal, V.; Lenhart, L. MHD flow in liquid-metal-cooled blankets. Fusion Eng. Des. 1991, 14, 401–412. [Google Scholar] [CrossRef]
  9. Davidson, P.A. An Introduction to Magnetohydrodynamics; Cambridge Texts in Applied Mathematics; Cambridge University Press: Cambridge, UK, 2001. [Google Scholar]
  10. Roberts, P.H. An Introduction to Magnetohydrodynamics; Elsevier: Alpharetta, GA, USA, 1967. [Google Scholar]
  11. Codina, R.; Hernadez-Silva, N. Stabilized finite element approximation of the stationary magneto-hydrodynamics equations. Comput. Mech. 2006, 38, 344–355. [Google Scholar] [CrossRef]
  12. Gunzburger, M.D.; Meir, A.J.; Peterson, J.S. On the existence and uniqueness and finite element approximation of solutions of the equations of stationary incompressible magnetohydrodynamics. Math. Comput. 1991, 56, 523–563. [Google Scholar] [CrossRef]
  13. Barbosa, H.; Hughes, T. The finite element method with Lagrange multipliers on the boundary: Circumventing the Babuska-Brezzi condition. Comput. Meth. Appl. Mech. Engrgy 1991, 85, 109–128. [Google Scholar] [CrossRef]
  14. Greif, C.; Li, D.; Schotzau, D.; Wei, X. A mixed finite element method with exactly divergence-free velocities for incompressible magnetohydrodynamics. Comput. Methods Appl. Mech. Engrgy 2010, 45–48, 2840–2855. [Google Scholar] [CrossRef]
  15. Hasler, U.; Schneebeli, A.; Schözau, D. Mixed finite element approximation of incompressible MHD problems based on weighted regularization. Appl. Numer. Math. 2004, 51, 19–45. [Google Scholar] [CrossRef]
  16. Schneebeli, A.; Schötzau, D. Mixed finite elements for incompressible magnetohydrodynamics. C. R. Acad. Sci. Paris Ser. I 2003, 337, 71–74. [Google Scholar] [CrossRef]
  17. Gerbeau, J.F. A stabilized finite element method for the incompressible magnetohydrodynamic equations. Numer. Math. 2000, 87, 83–111. [Google Scholar] [CrossRef]
  18. Shi, D.; Yu, Z. Nonconforming mixed finite element methods for stationary in-compressible Magnetohydrodynamics. Int. J. Numer. Anal. Model. 2013, 10, 904–919. [Google Scholar]
  19. Bochev, P.B.; Gunzburger, M.D. Least-squares method for the velocity pressure-stress formulation of the Stokes equations. Comput. Meth. Appl. Mech. Eng. 1995, 126, 267–287. [Google Scholar] [CrossRef] [Green Version]
  20. Bramble, J.H.; Pasciark, J.E. Least-squares methods for the Stokes equations based on a discrete minus one inner project. J. Comput. Appl. Math. 1996, 74, 155–173. [Google Scholar] [CrossRef] [Green Version]
  21. Bramble, J.H.; Schatz, A.H. Least-squares methods for 2m-th order elliptic boundary-value problems. Math. Comput. 1971, 25, 1–32. [Google Scholar]
  22. Chang, C.L.; Jiang, B.N. An error analysis of least-squares finite element method of velocity- pressure- vorticity formulation for the Stokes problem. Comput. Methods Appl. Mech. Eng. 1990, 84, 247–255. [Google Scholar] [CrossRef]
  23. Yang, D.; Wang, L. Quasi-least-squares finite element method for steady flow and heat transfer with system rotation. Numer. Math. 2006, 104, 377–411. [Google Scholar] [CrossRef]
  24. Qureshi, I.H.; Nawaz, M.; Rana, S.; Nazir, U.; Chamkha, A.J. Investigation of variable thermo-physical properties of viscoelastic rheology: A Galerkin finite element approach. AIP Adv. 2018, 8, 075027. [Google Scholar] [CrossRef] [Green Version]
  25. Wang, L.; Pang, S.Y.; Cheng, L. Bifurcation and stability of forced convection in tightly coiled ducts: Stability. Chaos Solitons Fractals 2006, 26, 991–1005. [Google Scholar] [CrossRef]
  26. Bochev, P.B. Analysis of least-squares finite element method for Navier-Stokes equations. SIAM J. Numer. Anal. 1997, 34, 1817–1844. [Google Scholar] [CrossRef] [Green Version]
  27. Bochev, P.B. Negative norm least-squares finite element mfethod for velocity vorticity pressure Navier Stokes equations. Numer. Meth. Part. Diff. Eqs. 1999, 15, 237–256. [Google Scholar] [CrossRef]
  28. Cai, Z.; Lazarov, R.; Manteuffel, T.A.; McCormick, S.F. First-order least squares for second-order partial differential equations: Part I. SIAM J. Numer. Anal. 1994, 6, 1785–1799. [Google Scholar] [CrossRef] [Green Version]
  29. Chang, C.L.; Yang, S.Y. Analysis of the L2 least-squares finite method for the velocity-vorticity-pressure Stokes equations with velocity boundary conditions. Appl. Math. Comput. 2002, 130, 121–144. [Google Scholar] [CrossRef]
  30. Hussain, S.; Mahbub, M.A.A.; Nasu, N.J.; Zheng, H. stabilized lowest equal-order mixed finite element method for the oseen viscoelastic fluid flow. Adv. Differ. Eq. 2018, 2018, 461. [Google Scholar] [CrossRef] [Green Version]
  31. Selmi, M.; Naudakumar, K. Bifurcation study of flow through a rotating curved duet. Phys. Fluids 1999, 11, 2030–2043. [Google Scholar] [CrossRef]
  32. Codina, R.; Hernadez, N. approximation of the thermally coupled MHD problemusing a stabilized finite element method. J. Comput. Phys. 2010, 230, 1281–1303. [Google Scholar] [CrossRef] [Green Version]
  33. Schötzau, D. Mixed finite element methods for stationary incompressible magneto-hydrodynamics. Numer. Math. 2004, 96, 771–800. [Google Scholar] [CrossRef]
  34. Gunzburger, M.D.; Ladyzhenskaya, O.A.; Peterson, J.S. On the global unique solvability of initial-boundary value problems for the coupled modified Navier-Stokes Maxwell equations. J. Math. Fluid Mech. 2004, 6, 462–482. [Google Scholar] [CrossRef]
  35. Prohl, A. Convergent finite element discretizations of the nonstationary incompressible Magnetohydrodynamics system. Esaim Math. Model. Numer. Anal. 2008, 42, 1065–1087. [Google Scholar] [CrossRef] [Green Version]
  36. Ladyzhenskaya, O.A.; Solonnikov, V. Solution of Some Nonstationary Magnethydrodynamical Problems for Incompressible Fluid. Trusy Steklov Math. Inst. 1960, 59, 115–173. [Google Scholar]
  37. Hughes, W.F.; Young, F.J. The Electromagneto-Hydrodynamics of Fluids; Wiley: New York, NY, USA, 1966. [Google Scholar]
  38. Jackson, J.D. Classical Electrodynamics; Wiley: New York, NY, USA, 1975. [Google Scholar]
  39. Monk, P. Finite Element Methods for Maxwell’s Equations; Oxford University Press: New York, NY, USA, 2003. [Google Scholar]
  40. Adams, R.; Fournier, J. Sobolev Spaces. Pure and Applied Mathematics; Elsevier Science: Amsterdam, The Netherlands, 2003. [Google Scholar]
  41. Adams, R.A. Sobolev Space. Pure and Applied Mathematics; Academic Press: New York, NY, USA, 1975; Volume 65. [Google Scholar]
  42. Girault, V.; Raviart, P.A. Finite Element Methods for Navier-Stockes Equations, Theory and Algorithms; Springer: Berlin, Germany, 1986. [Google Scholar]
  43. Chang, C. An Error Analysis of Least Squares Finite Element Method of Velocity-Pressure-Vorticity Formulation for Stokes Problem: Correction; Mathematics Research Report; Department of Mathematics, Cleveland State University: Cleveland, OH, USA, 1995. [Google Scholar]
  44. Chang, C.; Li, J.; Xiang, X.; Yu, Y.; Ni, W. Least squares finite element method for electromagnetic fields in 2-D. Appl. Math. Comput. 1993, 58, 143–167. [Google Scholar]
  45. Carey, G.F.; Jiang, B.N. Least-squares finite elements for first-order hyperbolic systems. Int. J. Num. Meth. Eng. 1988, 26, 81–93. [Google Scholar] [CrossRef]
  46. Hecht, F. FreeFEM++. J. Numer. Math. 2012, 20, 251–265. [Google Scholar] [CrossRef]
  47. Yang, J.; Mao, S.; He, X.; Yang, X.; He, Y. A diffuse interface model and semi-implicit energy stable fnite element method for two-phase magnetohydrodynamics flows. Comput. Methods Appl. Mech. Eng. 2019, 356, 435–464. [Google Scholar] [CrossRef]
Figure 1. Mesh configuration.
Figure 1. Mesh configuration.
Inventions 07 00040 g001
Figure 2. Fluid path or velocity contours.
Figure 2. Fluid path or velocity contours.
Inventions 07 00040 g002
Figure 3. The pressure field lines.
Figure 3. The pressure field lines.
Inventions 07 00040 g003
Figure 4. The magnetic field lines.
Figure 4. The magnetic field lines.
Inventions 07 00040 g004
Table 1. The error estimate for MHD with standard FE P 1 b P 2 P 1 pair.
Table 1. The error estimate for MHD with standard FE P 1 b P 2 P 1 pair.
h | | v f v f h | | L 2 | | p p h | | L 2 | | M B M B h | | L 2
1 / 4 0.221220.205620.12221
1 / 8 0.102670.064300.04332
1 / 16 0.048830.021870.01579
1 / 32 0.023900.007240.00621
1 / 64 0.011850.002470.00253
R 333
Table 2. The error estimate for MHD with standard FE P 2 P 2 P 1 pair.
Table 2. The error estimate for MHD with standard FE P 2 P 2 P 1 pair.
h | | v f v f h | | L 2 | | p p h | | L 2 | | M B M B h | | L 2
1 / 4 0.04185640.04861370.00786763
1 / 9 0.003500270.009566520.000672493
1 / 16 0.0006084290.003025950.000119331
1 / 25 0.0001598060.001239373.12644 × 10 5
1 / 36 5.70465 × 10 5 0.000597691.04688 × 10 5
R 3.0022.9993.00
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hussain, S.; Rahman, S.u.; Abbas, S.; Abbas, M.A. Solution of Steady Incompressible MHD Problems with Quasi-Least Square Method. Inventions 2022, 7, 40. https://doi.org/10.3390/inventions7020040

AMA Style

Hussain S, Rahman Su, Abbas S, Abbas MA. Solution of Steady Incompressible MHD Problems with Quasi-Least Square Method. Inventions. 2022; 7(2):40. https://doi.org/10.3390/inventions7020040

Chicago/Turabian Style

Hussain, Shahid, Shams ur Rahman, Suhail Abbas, and Munawwar Ali Abbas. 2022. "Solution of Steady Incompressible MHD Problems with Quasi-Least Square Method" Inventions 7, no. 2: 40. https://doi.org/10.3390/inventions7020040

APA Style

Hussain, S., Rahman, S. u., Abbas, S., & Abbas, M. A. (2022). Solution of Steady Incompressible MHD Problems with Quasi-Least Square Method. Inventions, 7(2), 40. https://doi.org/10.3390/inventions7020040

Article Metrics

Back to TopTop