Analysis of Control Problems for Stationary Magnetohydrodynamics Equations Under the Mixed Boundary Conditions for a Magnetic Field

: The optimal control problems for stationary magnetohydrodynamic equations under the inhomogeneous mixed boundary conditions for a magnetic ﬁeld and the Dirichlet condition for velocity are considered. The role of controls in the control problems under study is played by normal and tangential components of the magnetic ﬁeld given on different parts of the boundary and by the exterior current density. Quadratic tracking-type functionals for velocity, magnetic ﬁeld or pressure are taken as cost functionals. The global solvability of the control problems under consideration is proved, an optimality system is derived and, based on its analysis, a mathematical apparatus for studying the local uniqueness and stability of the optimal solutions is developed. On the basis of the developed apparatus, the local uniqueness of solutions of control problems for speciﬁc cost functionals is proved, and stability estimates of optimal solutions are established.


Introduction
Great attention has recently been paid to optimal control problems for the models of magnetohydrodynamics (MHD) for a viscous conducting incompressible fluid.There are a number of papers devoted to the theoretical study of such problems.Among these, we mention the papers [1][2][3][4][5][6][7][8], devoted to studying the control problems for a stationary MHD system, and papers [9][10][11][12], where the authors study control problems in a nonstationary case.In these papers, solvability of the control problems is proved, and optimality systems that describe the necessary conditions of the extremum are constructed and studied.In [5], the uniqueness and stability of the solutions of control problems are studied for some particular cases.
Along with optimal control problems, inverse or identification problems for the MHD models as well as for other hydrodynamic models play a key role in applications.In these problems, some parameters (constant or functional) that are included into a boundary value problem under study are unknown and are required to be determined together with the solution by using additional information about the state of the system.It is significant that the identification problems can be reduced to appropriate control problems by choosing a suitable tracking-type cost functional.As a result, control and identification problems can be analyzed using a common approach based on the theory of smooth-convex extremum problems in Hilbert or Banach spaces [13,14] (see also [15]).
In this paper, we study control problems for the stationary MHD model in which magnetohydrodynamic equations are considered under mixed boundary conditions for a magnetic field and under the Dirichlet boundary condition for velocity (see Section 2).
The role of controls in control problems under study is played by normal and tangential components of the magnetic field given on the different parts of the boundary and by exterior current density given within the flow region.Our goal is to construct a theory for studying local uniqueness and the stability of optimal solutions.Reaching of this goal will be based on generalization of the approach developed in [5,16] for studying the stability of optimal solutions for the stationary Navier-Stokes and MHD systems related to small perturbations of the cost functionals to be minimized.In this approach, there is no requirement to determine the second derivative of the cost functional under minimization, since it is based on the analysis of fundamental properties of the optimality system for the control problem under study and the use of special estimates for the difference of solutions of the original and perturbed control problems.The method is rather simple, natural and applicable for the models of hydrodynamics, heat convection, mass transfer and other hydrodynamic models based on the Navier-Stokes system [15].The use of this approach allows us to obtain stability estimates for optimal solutions with respect to small perturbations of the cost functional in an explicit and sufficiently easy-to-interpret form.
One of the features of this work consists in proving the solvability of optimal control problems under minimum requirements on normal and tangential components of the magnetic field given on the different parts of boundary and playing the role of controls in control problems.This result, in contrast to [5], allows the use of simple L 2 norms for the normal and tangential components of the magnetic field instead of the standard H 1/2 norms of boundary controls as Tikhonov regularizers while studying control problems for the MHD system.This regularization is needed to prove the local uniqueness and stability of the control problem solutions.The latter is the main goal of the study.
The structure of this paper is as follows.In Section 2, the boundary value problem is formulated for stationary MHD equations under inhomogeneous mixed boundary conditions for a magnetic field, and some notations that will be used throughout the paper are introduced along with the presentation of some additional facts that are necessary for studying optimal control problems.Presented in Section 3 are two lemmas on the existence of lifting of the velocity and of the magnetic component of the solution of the boundary value problem under consideration.In addition, the unique solvability of a generalized linear analogue of the original MHD boundary value problem is proved.In Section 4, two optimal control problems are formulated, and their solvability is proved, the optimality system describing the first-order necessary optimality conditions is derived and, on the basis of its analysis, the additional properties of optimal solutions are established.Finally, in Section 5, we will prove the local uniqueness and stability of the solutions for the control problems for "magnetic field-tracking" or "velocity-tracking" cost functionals.

Statement of the Boundary Value Problem and Notation
While studying the flows of electrically conducting fluids in real-life devices, the necessity often arises in flows of conducting fluid modeling in domains with boundaries consisting of parts with different electrophysical properties.Mathematical modeling of conducting flows in such types of domains gives rise to the study of boundary value problems for MHD equations under the mixed boundary conditions for a magnetic field.
Let Ω be a bounded domain of space R 3 with boundary Σ = ∂Ω consisting of two parts Σ τ and Σ ν .In this paper, we study control problems for the stationary magnetohydrodynamic equations of a viscous incompressible fluid considered in domain Ω under the following inhomogeneous boundary conditions: Here u is the velocity vector, H and E are magnetic and electric fields, respectively, p = P/ρ 0 , where P is the pressure, ρ 0 = const is a fluid density, κ = µ/ρ 0 , ν 1 = 1/ρ 0 σ = κν m , ν and ν m are constant kinematic and magnetic viscosity coefficients, σ is a constant electrical conductivity, µ is a constant magnetic permeability, n is the outer normal to ∂Ω, and j is the exterior current density.Further, we will refer to problem (1)-(3) for given f, j, g, q, q and k as Problem 1.It should be noted that all the quantities used in (1)-( 3) are dimensional and, moreover, their physical dimensions are defined in terms of SI units.
For the first time, global solvability of the homogeneous analogue of Problem 1 was proved in [23].In this work, mathematical tools from [24][25][26] were essentially used.Global solvability of the inhomogeneous mixed Problem 1 was proved in [27].In [28], the results obtained in [27] were generalized for the model of heat conducting magnetohydrodynamics.It should be noted that papers [29][30][31], magnetohydrodynamic equations are studied under mixed boundary conditions with respect to velocity and under the standard boundary conditions of the first type for an electromagnetic field.In [32], the author proves the existence of a very weak solution of the MHD boundary value problem using the Dirichlet boundary condition for a magnetic field.In [33], the MHD boundary value problem is studied in the case where a pressure and zero tangential velocity component are specified on an entire boundary.A local solvability of boundary value and boundary control problems for the model of magnetomicropolar flow is proved in paper [6], where the optimality systems are also derived and analyzed for the control problem.
Papers [4,28] and a number of others study the solvability of boundary value problems for stationary and nonstationary magnetohydrodynamics-Boussinesq systems.In [34], the authors study the properties of solutions (and, in particular, solvability) of the heat conducting magnetic hydrodynamic equations with the buoyancy effects due to temperature differences in the flow, Joule and viscous heating effects.In [35,36], the authors study the solvability of boundary value problems respectively for a steady or nonsteady MHD-Boussinesq system considered under mixed boundary conditions for temperature, magnetic field and velocity, in the general case when the thermal conductivity viscosity coefficient, electrical conductivity, magnetic permeability and specific heat of the fluid depend on the temperature.
We will assume below that the domain Ω and the boundary Σ partitioning into parts Σ τ and Σ ν satisfy (i) Ω is a bounded domain in the space R 3 , and its boundary Σ = ∂Ω consists of m+1 disjoint closed C 2 -surfaces Σ 0 , Σ 1 ,. . ., Σ m , each of which has a finite area where Σ 0 is the outer boundary of Ω.
(ii) Sets Σ τ and Σ ν are not empty, open and the following conditions take place: The conditions in (ii) mean that each of the parts Σ τ and Σ ν consists of a finite number of connected components of the boundary Σ.
It is assumed that, in the general case, Ω is a multi-connected domain and by n ≥ 1, we denote the number of handles of Σ.The case n = 0 corresponds to a simply connected domain.
The numbers n and m are respectively called the first and second Betti numbers (see, e.g., [15] (p.277) and [37]).Typical examples of domain Ω are shown in Figures 1 and 2, where a simply connected domain Ω with a disjoint boundary Σ = Σ 0 ∪ Σ 1 ∪ Σ 2 (with Betti numbers m = 2 and n = 0) and multi-connected toroidal domain Ω with a disjoint boundary Σ (with Betti numbers m = 2 and n = 1), respectively, are presented.Similarly to [37], we denote, by m + 1, the number of connected components of the part Σ ν of the boundary Σ, by Σ i , i = 1, . . ., m 0 ≤ m we denote the internal connected components of the boundary Let D(Ω) be the space of infinitely differentiable compactly supported functions in In addition to the spaces defined above, we will use the spaces and the space H DC (Ω) = H(div, Ω) ∩ H(curl, Ω), endowed with Hilbert norm defined by Here l is a dimensional factor of the dimension [l] = L 0 , and its value is equal to 1, L 0 denotes the SI dimension of the length (for more detail, see Section 4).
Any vector v defined on the boundary Σ (or on a part Σ 0 ⊂ Σ) can be decomposed to the sum of two vectors-the normal and tangential components v n and v 3 , the components v n and v T are given by formulas Here v n = v • n decribes the normal component of field v, v × n describes the tangential vector orthogonal to the normal n and to v T .It is clear that we denote the operator defined on H DC (Ω), which is placed in correspondence to every function h ∈ H DC (Ω) with the normal trace The following Green's formulae (see, e.g., [38]) will be used below: 3 the right-hand sides of ( 5) or (6) become Based on ( 5), ( 6) we say, following [25], that the function Similarly, we say that v × n = 0 weakly on Σ τ if Let H DCΣ τ (Ω) be the closure of C Σ τ 0 (Ω) 3 ∩H 1 (Ω) 3 with respect to norm • DC in (4).Set (7) are understood in the weak sense defined above.
Similarly, the basis of H Σ ν (Ω) consists of gradients ∇ ẑj of harmonic functions ẑj ∈ H 1 (Ω), j = m + 1, . . ., m satisfying the boundary conditions and harmonic vector fields y l ∈ H Σ ν (Ω) satisfying the condition Here ζ k denotes a cicle contained in Σ ν that is not homotopic to zero in Ω.
A number of important properties inherent in the function spaces defined above and proved in [25] are presented in the following lemma.Lemma 1.We assume that conditions (i), (ii) hold.Then: (1) H DCΣ τ (Ω) ⊂ H 1 (Ω) 3 and the norm • DC is equivalent to the norm • 1,Ω ; (2) there is a positive constant δ 1 depending on domain Ω and Σ τ such that the following holds: (3) the orthogonal decomposition of the space L 2 (Ω) 3 holds: Along with the spaces H DC (Ω) and H 0 (curl, Ω), their subspaces , will be used equipped, respectively, with norms The spaces H s+1/2 div (Ω) and H 0 Σ τ (curl, Ω) will be used below as the solution spaces for magnetic and electric components, respectively.In turn, the spaces will play the role of the solution spaces for velocity u.Besides Lemma 1, we also will use the following lemma (for details, see [25,26,38]).
Lemma 2. Under condition (i) there exist constants curlH Moreover, the following equality holds: The bilinear form defined by −(div•, •) satisfies Below, when formulating a result on the existence of magnetic lifting, we will need the space endowed with a norm and its subspace s ≥ 0, with the norm h Hs T (Σ ν ) = h s,Σ ν .Here div Σ ν is the linear surface divergence operator on the part Σ ν of the boundary Σ (see [15,37]).In the case s = 0, instead of HT (Σ ν ), we will write L2 T (Σ ν ).Let the following conditions for the data take place in addition to conditions (i), (ii): As usual, while studying control problems for the MHD system, we will deal with a weak form of Problem 1.This consists of finding a triple of functions (u, In order to obtain (19), one should multiply the first relation in (1) by v ∈ H 1 0 (Ω) 3 , the first relation in (2) by rot Ψ where Ψ ∈ V Σ τ (Ω), to integrate over Ω, to apply Green's formulas, to add the obtained results and to make use of the identity (see details in [27]) The identity (19) does not contain electric field E ∈ H 0 Σ τ (curl, Ω), which was eliminated with the help of ( 21).However, using a condition on a boundary vector k in (iii), vector E can be uniquely recovered from triple (u, 19) such that the first equation in relations (2) holds a.e. in Ω (see [27]).This allows us to refer, below, to the mentioned triple (u, H, p) satisfying relations (19), (20) as a weak solution to Problem 1.

Lifting of the Velocity and Magnetic Field, Solvability of Problem 1 and Its Linear Analogue
The proof of the existence of a solution of the inhomogeneous boundary value problem ( 1)-( 3) essentially uses the results for the existence of lifting of the velocity and magnetic field corresponding to the boundary conditions in (3).By velocity lifting, we mean the function u 0 ∈ H 1 div (Ω) satisfying the boundary condition u 0 | Σ = g.By magnetic field lifting, we mean the function existence of velocity lifting is ensured by the following lemma, proved in [3].Lemma 3.Under condition (i) for each function g ∈ H 1/2 T (Σ) and arbitrary number ε > 0, there is a function Here C ε is a constant depending on ε and Ω.
The existence of the lifting of the magnetic part of the solution is ensured by the following result, which is a generalization of Theorem 4.2 in [37], where it was proved in the case s = 0. Lemma 4. Let conditions (i), (ii) be satisfied.Then for any pair q ∈ H s (Σ τ ) and q ∈ Hs T (Σ ν ), where s ∈ [0, 1/2], there exists a unique function Here C Σ is a constant independent of q and q.
In what follows, the vector H 0 alone, defined in Lemma 4, will play the role of magnetic lifting for Problem 1.
Let us define two linear subspaces of H s+1/2 div (Ω): From Lemma 4, applied for the case q = 0 on Σ ν , it follows that for any function q ∈ H s (Σ τ ), there exists a unique solution h ∈ H s+1/2 div (Ω, Σ ν ) of problem (23) for q = 0, while for any function q ∈ Hs T (Σ ν ) there exists a unique solution h ∈ H s+1/2 div (Ω, Σ τ ) of problem (23) for the case q = 0.Moreover, when h runs through H s+1/2 div (Ω, Σ ν ), its normal component h • n| Σ τ runs through the space H s (Σ τ ).Similarly, when h runs through the space H s+1/2 div (Ω, Σ τ ), its tangential component h × n| Σ ν runs through the space Hs T (Σ ν ).Define the next products of spaces: Along with Problem 1, when studying the control problems below, an important role will be played by a linear analogue of Problem 1.This consists of finding a triple Here F ∈ W * is an arbitrary functional, "velocity" û ∈ H 1 div (Ω), "magnetic field" Ĥ ∈ H s+1/2 div (Ω) and χ ∈ L 2 0 (Ω) are given functions.In fact, we have somewhat generalized the linear analogue of Problem 1 by replacing the solenoidality condition divu = 0 by the more general condition divu = χ.Another generalization is that by the functional F in (26), we mean an arbitrary functional from W * , which does not necessarily coincide with the functional F defined in (19).
The following lemma about the unique solvability of problems ( 26), (27) holds.
Let us rewrite the problem ( 26), ( 27) in an equivalent operator form.To this end, we put in correspondence with the bilinear continuous forms (32), one can rewrite the identity (26) in the equivalent operator form By construction, Φ belongs to the space of continuous linear operators L(X, Y), and the original linear problems (26), ( 27) is equivalent to the operator equation From Lemma 5, it follows that (34) has a unique solution for any element y ∈ Y.This means that the operator Φ = (Φ 1 , Φ 2 , Φ 3 , Φ 4 , Φ 5 ) is surjective and invertible.Then it follows from Banach's inverse operator theorem that the operator Φ : X → Y is an isomorphism.Therefore, the following theorem holds.Theorem 1.Let assumptions (i), (ii) be satisfied.Then the operator Φ : X → Y, which is defined by ( 32), (33), is the isomorphism of the space X into Y.
We now formulate the following result concerning the sufficient conditions of the existence of the solution to Problem 1, which was proved in [27].

Statement of Control Problems, Optimality System and Additional Properties of Optimal Solutions
We note that problems (1)-( 3) contain constant parameters ν, ν 1 , ρ 0 , κ and functional parameters-boundary functions g, k, q, q and "volume" source densities f and j.To solve problems (1)-(3), one must specify values of respective parameters, boundary functions and sources.In practice, however, some of their specific elements may be unknown, and one should determine them and the solution (u, H, p) using certain information about the solution.
In this section, control problems for the MHD system (1)-(3) will be considered, for which we prove their global solvability.These problems are to minimize the so-called cost functionals, which depend on the variables (u, H, p) of the main state and other unknown functions (controls) satisfying the state Equations ( 1)-(3).We will choose one specific cost functional: Here the function 3 and p d ∈ L 2 (Q) have a similar sense for magnetic field or pressure.We note that functionals I 1 , I 2 , I 3 are used to solve inverse problems for the MHD system (1)-(3) using the optimization method.This method was developed by A.N. Tikhonov in the process of creating the famous Tikhonov regularization method [39].Currently, the optimization method is one of the fundamental methods for solving inverse problems arising in electromagnetism, acoustics, fluid mechanics, heat and mass transfer, design of complicated technical devices and in other fields of physics, natural science and engineering (for more detail, see [40][41][42][43][44][45][46][47][48][49][50][51]).
As controls in this paper, we choose three functions: q, q and j.The function j will play the role of distributed control, while q and q will play the role of boundary controls in the control problems stated below.We assume that the functions q, q and j change over sets K 1 , K 2 and K 3 satisfying the following conditions: (j) 3 are nonempty convex closed sets.Setting u = (q, q, j) and rewrite the weak form ( 19), (20) of Problem 1 as operator equation Let us introduce a functional I : X → R and nonnegative parameters µ 0 , µ 1 , µ 2 , µ 3 , and let us presuppose that the following conditions take place, in addition to (j): (jj) I : X → R is a weakly lower semicontinuous cost functional; (jjj) µ 0 > 0, µ 1 ≥ 0, µ 2 ≥ 0, µ 3 ≥ 0 and K 1 , K 2 , K 3 are limited sets, or µ 1 > 0, µ 2 > 0, µ 3 > 0 and the functional I : X → R is bounded below.
The following problem will be considered below: The parameters µ 0 , µ 1 , µ 2 , µ 3 serve in regulating the relative contribution of each of the terms in (40) and, moreover, to adjust their dimensions.Another purpose of using µ l is to ensure the uniqueness and stability for the control problems under study (see below).The last three terms in the structure of J are penalization terms.In what follows, we will refer to these terms as regularizers (strong in the case s > 0 and weak for s = 0).As will be shown in Section 5, just the presence of these regularizers in the structure of the functionals under minimization will allow us to prove the theorems concerning both the uniqueness and stability for optimal solutions.
By the conditions of Theorem 3 and in virtue of Theorem 2, we have the following estimates for q m , q m , j m , u m , H m , p m : Here c 1 , c 2 , . . .are some constants that do not depend on m.From the given estimates, it follows that there are weak limits 3 and strongly in L 4 (Ω) 3 , It is clear that To this end, we note that u m , H m , p m and j m satisfy the identity Let us pass to a limit in (43) as m → ∞.From (41), it follows that all linear terms in (43) pass to corresponding linear terms in (42) as m → ∞.Let us treat nonlinear terms beginning with ((u m • ∇)u m , v).Since u m ∈ H 1 T (Ω), then arguing as in [3], we easily derive that We now consider the second nonlinear term (rot Let us prove firstly that for the first term in (45), we have For this purpose, it is sufficient to prove that for an arbitrary pair, a number ε > 0 and a test function v ∈ H 1 0 (Ω) 3 , there exists a number As H 1 0 (Ω) is the closure of D(Ω) by norm • 1,Ω , then for mentioned function v ∈ H 1 0 (Ω) 3 there exists a sequence v n ∈ D(Ω) 3 , converging to v by norm • 1,Ω as n → ∞.It is clear that for all m, n ∈ N, we have From (41) follows that the norms rot H m Ω and H m − H * L 3 (Ω) 3 are uniformly bounded.Therefore for all m ∈ N there exists a number N=N(ε, v) such that for the second term in the right-hand side of (48) we have Let us consider the first term on the right-hand side of (48).Based on Hölder inequality for three functions, we have Since rot From the obtained inequality and from ( 48), (49), the estimate (47) follows.Since ε is an arbitrary positive number, the relation ( 46) is proved.Now we turn to the second term on the right-hand side of ( 45) and prove that From the weak convergence of the sequence H m in the space H s+1/2 div (Ω), the weak convergence of rot H m to rot H * in L 2 (Ω) 3 follows.As H * × v ∈ L 2 (Ω) 3 , we have Therefore, (50) is proved.From ( 46) and ( 50), it follows from ( 45) that Using the analogous scheme, one can show that for the last nonlinear term in (43), we have As a result, passing to the limit in (43) as m → ∞, we arrive by ( 44), ( 51) and ( 52) at (42).Finally, since J below is weakly semicontinuous on X × K functional, we have that According to Theorem 3, the solution of the control problem (40) exists for any value s ∈ [0, 1/2].The case s = 0, corresponding to weak regularizers when q ∈ L 2 (Σ τ ), q ∈ L 2 T (Σ ν ), is physically the most interesting.Another physically interesting case is the control problem with weak regularizers having the form provided that (x, u) ∈ X × K for s > 0 as in the situation of problem (40) with strong regularizers.Let the following additional condition apply: The proof of the following theorem is carried out analogously to the proof of Theorem 3. Theorem 4. Let the assumptions (i)-(iii) and (j), (jj), (jv) hold, and let the set Z ad be nonempty.Then the problem (53) has a solutuion (x, u) ∈ X × K for any s ∈ (0, 1/2]. Our next goal is to analyze the uniqueness and stability of the solutions of our control problems.For this purpose, we apply the technique developed in [5], which is based on using the additional properties of optimal solutions obtained by analyzing the optimality system.Therefore, the next stage of our study is to obtain the required optimality system.For concreteness, we will consider the case of control problem (40).Preliminarily, we define the spaces , which are dual of spaces X and Y defined in (25).
Below we will assume that multiplier λ 0 is dimensionless.Dimensions of adjoint state variables ξ, η, σ, ζ 1 , ζ 2 , ζ 3 clearly depend on dimension [µ 0 ] of parameter µ 0 .We assume that the dimension [µ 0 ] is chosen so that the dimensions of the adjoint state variables ξ, η and σ coincide with those of u, H and p in the main state x = (u, H, p), i.e., so that Here and below, L 0 , T 0 , I 0 and M 0 denote the SI dimensions of the length, time, electric current and mass units expressed in meters, seconds, amperes and kilograms, respectively.This allows us to refer to ξ, η and σ as "adjoint velocity", "adjoint magnetic field" and "adjoint pressure".64) that the adjoint velocity ξ and adjoint magnetic field η possess properties

Remark 2. It follows from conditions
Here χ Q is a characteristic function of the set Q. We emphasize that the adjoint velocity ξ, unlike η, is in a general case a nonsolenoidal vector function except for the case when the cost functional I is independent of pressure p.Only in this case div ξ = 0 and, moreover, ξ ∈ V.
To prove the uniqueness and also stability of the solution of (40), we need to introduce additional conditions for the data depending on the cost functional.Below, we preliminarily derive one important inequality with respect to the difference of a solution (x 1 , u 1 ) of problem (40) and a solution (x 2 , u 2 ) of the perturbed problem (40).In addition, we derive the estimates for the difference x 1 − x 2 via the difference u 1 − u 2 = (q 1 − q 2 , q 1 − q 2 , j 1 − j 2 ).
Let us denote by (x 1 , u 1 ) ≡ (u 1 , H 1 , p 1 , q 1 , q 1 , j 1 ) ∈ X × K an arbitrary solution to problem (40).By (x 2 , u 2 ) ≡ (u 2 , H 2 , p 2 , q 2 , q 2 , j 2 ) ∈ X × K, we denote the solution to problem Here Ĩ is a functional which is close to functional I.By virtue of Theorem 2, we have for triples (u i , H i , p i ), i = 1, 2: Let the values M 0 u and M 0 H be such that To make conditions (78) more illustrative and to simplify the subsequent presentation, we define the parameters They are analogues of the hydrodynamic dimensionless parameters [53], namely the Reynolds number Re, the magnetic Reynolds number Rm, the Hartman number Ha and the magnetic Prandtl number Pm.We emphasize that parameters Re, Rm, Ha and P m are dimensionless.To demonstrate this fact, one should know the dimensions of all parameters δ 0 , δ 1 , γ 0 , γ 1 and β defined in Lemma 2 and also of M 0 u and M 0 H entering into (78).In order to determine the dimensions of Re, Rm, Ha and P m, it will be assumed, below, the norms • Ω and • 1,Ω of a function u in L 2 (Ω) and in H 1 (Ω) and seminorm | • | 1,Ω in H 1 (Ω) are defined as follows: Here l is the dimensional factor having the dimension [l] = L 0 with value equal to 1. Using (80), one can verify that the dimensions of u Ω , |u| 1,Ω and u 1,Ω are connected with the dimension [u] of u by the formulas We recall also (see, e.g., [15] (p.272)) that Combining this with ( 10)-( 15), ( 17) and (74) yields Using (81), we have that i.e., all the parameters Re, Rm, Ha and P m defined in (79) are dimensionless.Since parameters Re and Rm are connected by relation Rm = (γ 1 /γ 0 )P mRe, we can rewrite conditions (78) in the following form containing only three nondimensional parameters Re, Ha and P m: 3 ), i = 1, 2, the Lagrange multipliers that correspond to solutions (x 1 , u 1 ) and (x 2 , u 2 ) of problems ( 40) and (76), respectively (these multipliers are determined uniquely under conditions (82)).By definition, they satisfy identities Here we renamed I = I 1 , Ĩ = I 2 .We define the following differences: We now deduce an important inequality for differences (85).While obtaining the inequality, we will use some ideas and results from [5], which we will sketch here for the reader's convenience.
We now consider problems (86), (87) with respect to differences u = u 1 − u 2 , H = H 1 − H 2 and p = p 1 − p 2 .In this problem, differences j = j 1 − j 2 , q = q 1 − q 2 and q = q 1 − q 2 play the role of the data together with functions u 1 , u 2 , H 1 and H 2 .Below, we will need estimates for norms of the differences u, H and p via norms of the differences q, q and j.In order to derive these estimates, we present the difference H ≡ H 1 − H 2 in the form H = H 0 + H where H 0 is a unique solution of ( 23) corresponding to differences q = q 1 − q 2 and q = q 1 − q 2 , while H ∈ V Σ τ (Ω) is a certain function.We set v = u and Ψ = H in (86).Taking into account (16) and condition rotH 0 = 0, we obtain ν(∇u, ∇u) Using estimates ( 10)-( 15), (77) and setting ν 1 = ν m κ, we obtain from (92) that Applying (79) and Young's inequality |ab| ≤ εa 2 /2 + b 2 /(2ε) or 2ab ≤ εa 2 /2 + 2b 2 /ε, ε = const > 0, at ε = δ 0 ν/2, ε = δ 1 ν m κ/4 or ε = 1, we consequently derive that Taking into account (94)-(97), we obtain from (93) that It follows from (78) that Using (99), we derive from (98) that where From (100), we conclude that Since H = H 0 + H, from (102), (103), using (11), we obtain the needed estimates for differences u and H: Here C u and C H are dimensionless constants defined by Based on (17), we now derive a similar estimate for the difference p = p 1 − p 2 .In view of (17) for function p and for any (small enough) number δ > 0, there exists a function Setting v = v 0 , Ψ = 0 in (86), we obtain Using the previous estimate (107) for −(divv 0 , p) and ( 10), ( 13)-( 15) from (108), we deduce that Dividing (109) by v 0 1,Ω = 0 and using (104), (105) yields Taking into account (104) and (105), we obtain the following estimate: Here C p is a dimensionless constant defined by Remark 3. We note two peculiarities of the mathematical apparatus used in this paper.On the one hand, we use the dimensional MHD system (1)-( 3) so that all formulae used or obtained in the paper are dimensional.On the other hand, we essentially use three dimensionless parameters Re, Ha, P m defined in (79) and three dimensionless constants C u , C H , C p defined in (106), ( 112).These constants depend on Re, Ha and P m snf contain important information on the MHD system (1), (2).In particular, using C u , C H and C p , we could write estimates of the norms of differences u, H and p via norms of differences q, q and j in the simple form (104), ( 105) and (111).We recall that these estimates hold under the condition (82).Based on these estimates and Theorem 6, in Section 5, we establish similar stability estimates for a number of specific control problems for the MHD system (1)-(3) under study.
Again, we note that we cannot prove the estimates for differences u = u 1 − u 2 and p = p 1 − p 2 , which are analogous to (134), but we can obtain coarser estimates for u, H and p even if Q is only a part of Ω.Indeed, using (134), we infer from (142) that εµ 1 q 2 s,Σ τ + εµ 2 q s,Σ ν + εµ 3 j 2 Ω ≤ µ 0 (− where ∆ is defined in (134).From (104), ( 105), ( 111) and (143), we come to the required stability estimates having the form (130). Thus, the following result was proved.Using an analogous scheme, similar theorems can be proved for control problems corresponding to weak regularizers (µ 1 /2) q 2 Σ τ and (µ 2 /2) q 2 Σ ν as in (53), or to the third cost functional I s (p) defined in (37).We leave the formulations and proofs of the corresponding theorems to the reader.

Conclusions
In this paper, the optimal control problems for stationary magnetohydrodynamic equations considered under the inhomogeneous mixed boundary conditions for a magnetic field and the Dirichlet condition for velocity were analyzed.The normal and tangential magnetic components at different parts of the boundary and the exterior current density within the flow region are the controls.A general theory of investigation of the control problems with special tracking-type cost functionals has been developed and applied for studying the uniqueness and stability for optimal solutions.The use of this theory made it possible to prove the existence of solutions of control problems in the wide class of cost functionals, to deduce an optimality system and, based on its analysis, to prove theorems on local uniqueness and stability of optimal solutions for a number of tracking-type functionals.In particular, for problems (113) and (131), it is shown that the uniqueness as well as stability of their solutions occurs, provided the conditions (115) or (133) for the data are accordingly fulfilled.This means that the presence of summands with squared norms of controls, with positive coefficients µ 1 , µ 2 , µ 3 in the expressions of functionals to be minimized in (113) and in (131), contributes a regularizing effect to the control problems under study.
Similarly, by n 1 (or n 2 ), we denote the number of handles of the part Σ ν (or Σ τ ).Clearly, n 1 + n 2 = n (where n is the number of handles of Σ).Below we will use the Sobolev spaces H s (D), s ∈ R, H 0 (D) ≡ L 2 (D), where D denotes Ω or the boundary Σ or a part Σ 0 ⊂ Σ.The corresponding spaces of vector functions are denoted by H s (D) 3 and L 2 (D) 3 .The inner products and norms in the spaces H s (D) and H s (D) 3 are denoted by (•, •) s,D and • s,D .The inner products and norms in L 2 (Ω) and L 2 (Ω) 3 are denoted by (•, •) and • Ω .By • 1,Ω and | • | 1,Ω we denote norm and seminorm in H 1 (Ω) or in H 1 (Ω) 3 .For arbitrary Hilbert space H by H * we denote the dual space of H.By H s T (Σ 0 ) we will denote the subspace in H s (Σ 0 ) 3 consisting of tangential on