Next Article in Journal
Tangled Cord of FTTM4
Next Article in Special Issue
Mathematical Analysis of Mixed Convective Peristaltic Flow for Chemically Reactive Casson Nanofluid
Previous Article in Journal
Omni-Domain Feature Extraction Method for Gait Recognition
Previous Article in Special Issue
Reconstruction of Unsteady Wind Field Based on CFD and Reduced-Order Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Gennadii Alekseev
1,2
1
Institute of Applied Mathematics FEB RAS, 7, Radio St., 690041 Vladivostok, Russia
2
Department of Mathematical and Computer Modelling, Far Eastern Federal University, 690922 Vladivostok, Russia
Mathematics 2023, 11(12), 2610; https://doi.org/10.3390/math11122610
Submission received: 29 May 2023 / Revised: 5 June 2023 / Accepted: 6 June 2023 / Published: 7 June 2023
(This article belongs to the Special Issue Theoretical Research and Computational Applications in Fluid Dynamics)

Abstract

:
The optimal control problems for stationary magnetohydrodynamic equations under the inhomogeneous mixed boundary conditions for a magnetic field 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 field given on different parts of the boundary and by the exterior current density. Quadratic tracking-type functionals for velocity, magnetic field 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 specific cost functionals is proved, and stability estimates of optimal solutions are established.

1. 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.

2. 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
ν Δ u + ( u · ) u + p κ rot H × H = f , div u = 0 ,
ν 1 rot H ρ 0 1 E + κ H × u = ν 1 j , div H = 0 , rot E = 0 ,
considered in domain Ω under the following inhomogeneous boundary conditions:
u | Σ = g , H · n | Σ τ = q , H × n | Σ ν = q , E × n | Σ τ = k .
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.
We emphasize that Equations (1) and (2) are considered under mixed boundary conditions with respect to electromagnetic field ( H , E ) in (3). These conditions generalize two types of previously used boundary conditions. The first type is described by relations H · n = 0 and E × n = 0 on Ω , corresponding to a perfectly conducting boundary (see, e.g., [2,3,4,5,17,18,19,20]). The second type is described by the condition H × n = 0 on Ω (see [8,21,22]) corresponding to a perfectly insulating boundary.
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 Figure 1 and Figure 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 Σ contained in Σ ν . It is clear that m ^ = m 0 if Σ 0 Σ ν and m ^ = m 0 1 if Σ 0 Σ τ . 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 T s ( Σ 0 ) we will denote the subspace in H s ( Σ 0 ) 3 consisting of tangential on Σ 0 Σ vector functions. Set H s ( Σ τ ) = H s ( Σ τ ) * , H T s ( Σ ν ) = H T s ( Σ ν ) * , H T 1 / 2 ( Σ ) = H T 1 / 2 ( Σ ) * for s 0 .
Let D ( Ω ) be the space of infinitely differentiable compactly supported functions in Ω , H 0 1 ( Ω ) be the closure of D ( Ω ) in H 1 ( Ω ) , V = { v H 0 1 ( Ω ) 3 : div v = 0 } , H 1 ( Ω ) 3 = ( H 0 1 ( Ω ) 3 ) * , L 0 2 ( Ω ) = { p L 2 ( Ω ) : ( p , 1 ) = 0 } , H 1 ( Ω , Σ τ ) = { φ H 1 ( Ω ) : φ | Σ τ = 0 } , C Σ τ 0 ( Ω ¯ ) 3 : = { v C 0 ( Ω ¯ ) 3 : v · n | Σ τ = 0 , v × n | Σ ν = 0 } . In addition to the spaces defined above, we will use the spaces H ( div , Ω ) = { v L 2 ( Ω ) 3 : div v L 2 ( Ω ) } , H ( curl , Ω ) = { v L 2 ( Ω ) 3 : curl v L 2 ( Ω ) 3 } , H 0 ( curl , Ω ) = { v H ( curl , Ω ) : curl v = 0 } and the space H D C ( Ω ) = H ( div , Ω ) H ( curl , Ω ) , endowed with Hilbert norm defined by
u D C 2 : = l 2 u Ω 2 + div u Ω 2 + curl u Ω 2 .
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 T : v = v n + v T . If v L 2 ( Σ ) 3 , the components v n and v T are given by formulas v n = ( v · n ) n v n n and v T = v v n ( n × v ) × n . 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 v T = 0 on Σ if v × n | Σ = 0 . By γ n | Σ 0 (or γ τ | Σ 0 ) we denote the operator defined on H D C ( Ω ) , which is placed in correspondence to every function h H D C ( Ω ) with the normal trace γ n | Σ 0 h = h · n | Σ 0 (or tangential trace γ τ | Σ 0 h = h × n | Σ 0 ).
The following Green’s formulae (see, e.g., [38]) will be used below:
Ω v · grad φ d x + Ω div v φ d x = Ω v · n φ d σ v H 1 ( Ω ) 3 , φ H 1 ( Ω ) ,
Ω ( v · curl w w · curl v ) d x = Ω ( v × n ) · w T d σ v , w H 1 ( Ω ) 3 .
If φ H 1 ( Ω , Σ τ ) or w C Σ τ 0 ( Ω ¯ ) 3 H 1 ( Ω ) 3 the right-hand sides of (5) or (6) become Σ ν v · n φ d σ or Σ τ ( v × n ) · w T d σ . Based on (5), (6) we say, following [25], that the function v H D C ( Ω ) weakly satisfies condition v · n = 0 on Σ ν if
Ω ( v · grad φ + div v φ ) d x = 0 φ H 1 ( Ω , Σ τ ) .
Similarly, we say that v × n = 0 weakly on Σ τ if
Ω ( v · curl w w · curl v ) d x = 0 w C Σ τ 0 ( Ω ¯ ) 3 H 1 ( Ω ) 3 .
Let H D C Σ τ ( Ω ) be the closure of C Σ τ 0 ( Ω ¯ ) 3 H 1 ( Ω ) 3 with respect to norm · D C in (4). Set
H Σ τ ( Ω ) = { h L 2 ( Ω ) 3 : div h = 0 , curl h = 0   in   Ω , h · n | Σ τ = 0 , h × n | Σ ν = 0 } ,
H Σ ν ( Ω ) = { h L 2 ( Ω ) 3 : div h = 0 , curl h = 0   in   Ω , h · n | Σ ν = 0 , h × n | Σ τ = 0 } ,
V Σ τ ( Ω ) = { v H D C Σ τ ( Ω ) : div v = 0   in   Ω } H Σ τ ( Ω ) .
Equalities h · n | Σ τ = 0 , h × n | Σ ν = 0 or h · n | Σ ν = 0 , h × n | Σ τ = 0 in (7) are understood in the weak sense defined above.
We remind that the spaces H Σ τ ( Ω ) and H Σ ν ( Ω ) are finite dimensional [37]. In particular, the dimension of H Σ τ ( Ω ) under conditions (i), (ii) is exactly m ^ + n 2 , and the basis of the space H Σ τ ( Ω ) consists of gradients z ^ j of harmonic functions z ^ j H 1 ( Ω ) , j = 1 , , m ^ satisfying boundary conditions
z ^ j / n | Σ τ = 0 , z ^ j | Σ ν Σ j = 0 , z ^ j | Σ j = 1 , j = 1 , , m ^ ,
and harmonic vector fields y l H Σ τ ( Ω ) , l = 1 , , n 2 satisfying the condition ς k y l · n d σ = δ l k for any cycle ς k , k = 1 , , n 2 , contained in Σ τ and not homotopic to zero in Ω . Here and below, δ l k is the Christoffel symbol equal to 1 at l = k and 0 at l k .
Similarly, the basis of H Σ ν ( Ω ) consists of gradients z ^ j of harmonic functions z ^ j H 1 ( Ω ) , j = m ^ + 1 , , m satisfying the boundary conditions
z ^ j / n | Σ ν = 0 , z ^ j | Σ τ Σ j = 0 , z ^ j | Σ j = 1 , j = m ^ + 1 , , m
and harmonic vector fields y l H Σ ν ( Ω ) satisfying the condition ζ k y l · n d σ = δ l k ζ k , k , l = n 2 + 1 , , n . 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 D C Σ τ ( Ω ) H 1 ( Ω ) 3 and the norm · D C is equivalent to the norm · 1 , Ω ;
(2) there is a positive constant δ 1 depending on domain Ω and Σ τ such that the following holds:
curl h 2 δ 1 h 1 , Ω 2 h V Σ τ ( Ω ) ;
(3) the orthogonal decomposition of the space L 2 ( Ω ) 3 holds:
L 2 ( Ω ) 3 = H 1 ( Ω , Σ τ ) curl H D C Σ τ ( Ω ) H Σ ν ( Ω ) .
Along with the spaces H D C ( Ω ) and H 0 ( curl , Ω ) , their subspaces
H div s + 1 / 2 ( Ω ) : = H s + 1 / 2 ( Ω ) 3 { h H D C ( Ω ) : div h = 0 } H Σ τ ( Ω ) ,
H Σ τ 0 ( curl , Ω ) : = { e H 0 ( curl , Ω ) : e × n | Σ τ L T 2 ( Σ τ ) } , 0 s 1 / 2 ,
will be used equipped, respectively, with norms
h H div s + 1 / 2 ( Ω ) = h s + 1 / 2 , Ω + curl h Ω , e H Σ τ 0 ( curl , Ω ) : = e Ω + e × n Σ τ .
The spaces H div s + 1 / 2 ( Ω ) and H Σ τ 0 ( curl , Ω ) will be used below as the solution spaces for magnetic and electric components, respectively. In turn, the spaces
H T 1 ( Ω ) = { v H 1 ( Ω ) : v · n | Σ = 0 } , H div 1 ( Ω ) : = { v H T 1 ( Ω ) : div v = 0 } ,
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 C 1 = C 1 ( Ω ) , δ i = δ i ( Ω ) > 0 , γ i ( Ω ) , γ i = γ i ( Ω ) > 0 , i = 0 , 1 , β = β ( Ω ) , depending on Ω such that
| ( u , w ) | u 1 , Ω w 1 , Ω u , w H 1 ( Ω ) 3 , ( v , v ) δ 0 v 1 , Ω 2 v H 0 1 ( Ω ) 3 ,
( curl Ψ , rot Ψ ) δ 1 Ψ 1 , Ω 2 Ψ V Σ τ ( Ω ) ,
curl H Ω H H div s + 1 / 2 ( Ω ) C 1 H 1 , Ω H H 1 ( Ω ) 3 .
| ( ( u · ) v , w ) | γ 0 u L 3 ( Ω ) 3 v 1 , Ω w 1 , Ω γ 0 u 1 , Ω v 1 , Ω w 1 , Ω u , v , w H 1 ( Ω ) 3 ,
| ( curl H 1 × H 2 , v ) | curl H 1 Ω H 2 L 3 ( Ω ) v L 6 ( Ω ) γ 1 H 1 1 , Ω H 2 s + 1 / 2 , Ω v 1 , Ω H 1 H ( curl , Ω ) , H 2 H s + 1 / 2 ( Ω ) 3 ,
| ( curl Ψ × h , v ) | rot Ψ Ω h L 3 ( Ω ) 3 u ^ L 6 ( Ω ) 3 γ 1 Ψ 1 , Ω H s + 1 / 2 , Ω v 1 , Ω H H div s + 1 / 2 ( Ω ) , v , Ψ H 1 ( Ω ) 3 .
Moreover, the following equality holds:
( ( u · ) v , v ) = 0 u H div 1 ( Ω ) , v H 0 1 ( Ω ) 3 .
The bilinear form defined by ( div · , · ) satisfies
sup v H 0 1 ( Ω ) 3 , v 0 ( div v , p ) / v 1 , Ω β p Ω p L 0 2 ( Ω ) .
Below, when formulating a result on the existence of magnetic lifting, we will need the space
H T 1 / 2 ( div Σ ν ; Σ ν ) = { h H T 1 / 2 ( Σ ν ) : div Σ ν h H 1 / 2 ( Σ ν ) } ,
endowed with a norm
h 1 / 2 , div , Σ ν 2 = h 1 / 2 , Σ ν 2 + div Σ ν h 1 / 2 , Σ ν 2 ,
and its subspace
H ˜ T s ( Σ ν ) = { h H T 1 / 2 ( div Σ ν ; Σ ν ) : div Σ ν h = 0   in   H 1 / 2 ( Σ ν ) ,
( h , m ) Σ ν = 0 m H Σ ν ( Ω ) , ( h , y l ) Σ ν = 0 , l = n 2 + 1 , , n } H T s ( Σ ν ) ,
s 0 , with the norm h H ˜ T s ( Σ ν ) = 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 H ˜ T ( Σ ν ) , we will write L ˜ T 2 ( Σ ν ) .
Let the following conditions for the data take place in addition to conditions (i), (ii):
(iii) f H 1 ( Ω ) 3 , g H T 1 / 2 ( Σ ) , k ( γ τ | Σ τ ) H Σ τ 0 ( curl , Ω ) ,
(iv) q H s ( Σ τ ) , q H ˜ T s ( Σ ν ) , j L 2 ( Ω ) 3 , s [ 0 , 1 / 2 ] .
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 , H , p ) H T 1 ( Ω ) × H div s + 1 / 2 ( Ω ) × L 0 2 ( Ω ) satisfying
ν ( u , v ) + ν 1 ( rot H , rot Ψ ) + ( ( u · ) u , v ) + κ [ ( rot Ψ × H , u ) ( rot H × H , v ) ] ( div v , p ) = F ( v , Ψ ) f , v + ν 1 ( j , rot Ψ ) + ρ 0 1 ( k , Ψ ) Σ τ ( v , Ψ ) H 0 1 ( Ω ) 3 × V Σ τ ( Ω ) ,
div u = 0   in   Ω , u = g   on   Σ , H · n = q   on   Σ τ , H × n = q   on   Σ ν .
In order to obtain (19), one should multiply the first relation in (1) by v H 0 1 ( Ω ) 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])
( E , curl Ψ ) = Σ τ ( E × n | Σ τ ) · Ψ T d σ = ( k , Ψ T ) Σ τ Ψ V Σ τ ( Ω ) .
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 , H , p ) H div 1 ( Ω ) × H div s + 1 / 2 ( Ω ) × L 0 2 ( Ω ) satisfying (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.

3. 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 div 1 ( Ω ) satisfying the boundary condition u 0 | Σ = g . By magnetic field lifting, we mean the function H 0 H div s + 1 / 2 ( Ω ) satisfying the mixed boundary conditions H 0 · n | Σ τ = q and H 0 × n | Σ ν = q . The existence of velocity lifting is ensured by the following lemma, proved in [3].
Lemma 3.
Under condition (i) for each function g H T 1 / 2 ( Σ ) and arbitrary number ε > 0 , there is a function u ε H T 1 ( Ω ) such that
div u ε = 0   in   Ω , u ε = g   on   Γ , u ε 1 , Ω C ε g 1 / 2 , Σ , u ε L 4 ( Ω ) 3 ε g 1 / 2 , Σ .
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 H ˜ T s ( Σ ν ) , where s [ 0 , 1 / 2 ] , there exists a unique function H 0 H div s + 1 / 2 ( Ω ) such that
rot H 0 = 0 , div H 0 = 0   in   Ω , H 0 · n = q   on   Σ τ , H 0 × n = q   on   Σ ν ,
H 0 H div s + 1 / 2 ( Ω ) = H 0 s + 1 / 2 , Ω C Σ ( q s , Σ τ + q s , Σ ν ) .
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 div s + 1 / 2 ( Ω ) :
H div s + 1 / 2 ( Ω , Σ ν ) = { h H div s + 1 / 2 ( Ω ) : h × n = 0   on   Σ ν } ,
H div s + 1 / 2 ( Ω , Σ τ ) = { h H div s + 1 / 2 ( Ω ) : h · n = 0   on   Σ τ } .
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 div s + 1 / 2 ( Ω , Σ ν ) of problem (23) for q = 0 , while for any function q H ˜ T s ( Σ ν ) there exists a unique solution h H div s + 1 / 2 ( Ω , Σ τ ) of problem (23) for the case q = 0 . Moreover, when h runs through H div s + 1 / 2 ( Ω , Σ ν ) , its normal component h · n | Σ τ runs through the space H s ( Σ τ ) . Similarly, when h runs through the space H div s + 1 / 2 ( Ω , Σ τ ) , its tangential component h × n | Σ ν runs through the space H ˜ T s ( Σ ν ) .
Define the next products of spaces:
X = H T 1 ( Ω ) × H div s + 1 / 2 ( Ω ) × L 0 2 ( Ω ) , H = H T 1 ( Ω ) × H div s + 1 / 2 ( Ω ) ,
W = H 0 1 ( Ω ) 3 × V Σ τ ( Ω ) , Y = W * × L 0 2 ( Ω ) × H T 1 / 2 ( Σ ) × H s ( Σ τ ) × H ˜ T s ( Σ ν ) .
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 ( u , H , p ) X H T 1 ( Ω ) × H div s + 1 / 2 ( Ω ) × L 0 2 ( Ω ) satisfying the functions
a u ^ H ^ ( ( u , H ) , ( v , Ψ ) ) ν ( u , v ) + ν 1 ( rot H , rot Ψ ) + ( ( u ^ · ) u , v ) + κ ( rot Ψ × H ^ , w )
( rot H × H ^ , v ) = F , ( v , Ψ ) + ( div v , p ) ( v , Ψ ) W ,
div u = χ   in   Ω , u | Σ = g , H · n | Σ τ = q   and   H × n | Σ ν = q .
Here F W * is an arbitrary functional, “velocity” u ^ H div 1 ( Ω ) , “magnetic field” H ^ H div s + 1 / 2 ( Ω ) and χ L 0 2 ( Ω ) are given functions. In fact, we have somewhat generalized the linear analogue of Problem 1 by replacing the solenoidality condition div u = 0 by the more general condition div u = χ . 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.
Lemma 5.
Let conditions (i), (ii) be satisfied. Then for any quintuple ( F , χ , g , q , q ) Y for arbitrary s [ 0 , 1 / 2 ] , there exists a unique solution ( u , H , p ) X of problems (26), (27), and for the solution the following estimates hold:
u 1 , Ω M ^ u 0 , H H div s + 1 / 2 ( Ω ) M ^ H 0 , p Ω M ^ p 0 .
Here M ^ u 0 , M ^ H 0 , M ^ p 0 are nondecreasing continuous functions of χ Ω , F W * , g 1 / 2 , Σ , q s , Σ τ , q s , Σ ν .
Proof of Lemma 5.
The existence of the solution of (26), (27) and estimate (28) are proved using the scheme proposed in [3] (see also [15], Chapter 6). For proving the uniqueness, let us assume that there exist two solutions ( u i , H i , p i ) X , i = 1 , 2 , of problem (26), (27). Then the quantities u = u 1 u 2 , H = H 1 H 2 and p = p 1 p 2 belong to V × V Σ τ ( Ω ) × L 0 2 ( Ω ) and satisfy
ν ( u , v ) + ν 1 ( curl H , curl Ψ ) + ( ( u ^ · ) u , v ) +
+ κ [ ( curl Ψ × H ^ , w ) ( curl H × H ^ , v ) ] ( div v , p ) = 0 ( v , Ψ ) W ,
div u = 0   in   Ω , u | Σ = 0 , H · n | Σ τ = 0 , H × n | Σ ν = 0 .
Setting v = u , Ψ = H in (29) and using (16), we arrive at the relation
ν ( u , u ) + ν 1 ( curl H , curl H ) = 0 .
where (31) takes place by (10), (11) if u = 0 and H = 0 or u 1 = u 2 and H 1 = H 2 in Ω . Then from (29), it follows that ( div v , p ) = 0 for all v H 0 1 ( Ω ) 3 . This means by inf-sup condition (17), that p = 0 or p 1 = p 2 in Ω . □
Let us rewrite the problem (26), (27) in an equivalent operator form. To this end, we put in correspondence with the bilinear continuous forms
a u ^ , H ^ : H × W ( H T 1 ( Ω ) × H div s + 1 / 2 ( Ω ) ) × ( H 0 1 ( Ω ) 3 × V Σ τ ( Ω ) ) R , b : W × L 0 2 ( Ω ) R ,
linear continuous operators A : H W * , B : W L 0 2 ( Ω ) * L 0 2 ( Ω ) , B * : L 0 2 ( Ω ) W * , acting by
A ( u , H ) , ( v , Ψ ) = a u ^ , H ^ ( ( u , H ) , ( v , Ψ ) ) ( u , H ) H   and   ( v , Ψ ) W ,
B ( v , Ψ ) , r = b ˜ ( ( v , Ψ ) , r ) b ( v , r ) = B * r , ( v , Ψ ) r L 0 2 ( Ω ) , ( v , Ψ ) W .
Using (32), one can rewrite the identity (26) in the equivalent operator form
A ( u , H ) + B * p = F .
Setting x = ( u , H , p ) , we define the operator Φ = ( Φ 1 , Φ 2 , Φ 3 , Φ 4 , Φ 5 ) : X Y by
Φ 1 ( x ) = A ( u , H ) + B * p , Φ 2 ( x ) = div u , Φ 3 ( x ) = u | Σ ,
Φ 4 ( x ) = H · n | Σ τ , Φ 5 ( x ) = H × n | Σ ν .
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
Φ ( x ) = y ( F , χ , g , q , q ) Y .
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].
Theorem 2.
Under assumptions (i)–(iv) there exists a weak solution ( u , H , p ) to Problem 1, and for this solution the following estimates hold:
u 1 , Ω M u , H H div s + 1 / 2 ( Ω ) M H , p Ω M p .
Here M u , M p and M H are nondecreasing continuous functions of norms f 1 , Ω , j Ω , k Σ τ , g 1 / 2 , Σ , q s , Σ τ , q s , Σ ν . If, besides, elements f , j , k , g , q and q are small (or “viscosity coefficients” ν, ν m are, vice versa, great) in the sense
γ 0 M u + γ 1 ( κ / 2 ) M H < δ 0 ν , γ 1 M u + γ 1 ( κ / 2 ) M H < δ 1 ν m ,
where constants δ 0 , δ 1 , γ 0 , γ 1 were defined in Lemmas 1 and 2, then the weak solution ( u , H , p ) is unique.
Remark 1.
Theorem 2 was proved in [27] under an additional condition on function g , namely that g is a tangential vector on Σ. If this condition does not hold, i.e., g · n | Σ ¬ 0 , one can prove only the local solvability of Problem 1. Thus, the problem of global solvability of inhomogeneous problem (1)–(3) in the general case when g · n | Σ 0 is still an open problem (see related discussion in [15,22]).

4. 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:
I 1 ( u ) = u u d Q 2 , I 2 ( H ) = H H d Q 2 , I 3 ( p ) = p p d Q 2 .
Here the function u d L 2 ( Q ) 3 describes, in a subdomain Q Ω velocity field, given functions H d L 2 ( Q ) 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) K 1 H s ( Σ τ ) , K 2 H ˜ T s ( Σ ν ) , K 3 L 2 ( Ω ) 3 are nonempty convex closed sets.
Setting u = ( q , q , j ) K = K 1 × K 2 × K 3 , x = ( u , H , p ) X H T 1 ( Ω ) × H div s + 1 / 2 ( Ω ) × L 0 2 ( Ω ) , we define an operator F ( F 1 , F 2 , F 3 , F 4 , F 5 ) : X × K Y ( H 1 ( Ω ) 3 × V Σ τ ( Ω ) * ) × L 0 2 ( Ω ) × H T 1 / 2 ( Σ ) × H s ( Σ τ ) × H ˜ T s ( Σ ν ) by
F 1 ( x , u ) , ( v , Ψ ) = ν ( u , v ) + ν 1 ( rot H , rot Ψ ) + ( ( u · ) u , v ) ( div v , p ) +
+ κ [ ( rot Ψ × H , u ) ( rot H × H , v ) ] f , v ν 1 ( j , rot Ψ )
ρ 0 1 ( k , Ψ T ) Σ τ ( v , Ψ ) H 0 1 ( Ω ) 3 × V Σ τ ( Ω ) ,
F 2 ( x ) , r = ( div u , r ) r L 0 2 ( Ω ) , F 3 ( x ) = u | Σ g H T 1 / 2 ( Σ ) ,
F 4 ( x , u ) = H · n | Σ τ q H s ( Σ τ ) , F 5 ( x , u ) = H × n | Σ ν q H ˜ T s ( Σ ν )
and rewrite the weak form (19), (20) of Problem 1 as operator equation
F ( x , u ) F ( u , H , p , q , q , j ) = 0 .
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:
J ( x , u ) = μ 0 2 I ( x ) + μ 1 2 q s , Σ τ 2 + μ 2 2 q s , Σ ν 2 + μ 3 2 j Ω 2 inf , F ( x , u ) = 0 , ( x , u ) X × K .
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.
Let Z a d = { ( x , u ) X × K , F ( x , u ) = 0 , J ( x , u ) < } be a set of possible pairs for problem (40).
Theorem 3.
Let the assumptions (i)–(iii), (j)–(jjj) hold, and let the set Z a d be nonempty. Then for arbitrary s [ 0 , 1 / 2 ] , the problem (40) has at least one solution ( x , u ) X × K .
Proof of Theorem 3.
Let us denote by ( x m , u m ) ( u m , H m , p m , q m , q m , j m ) Z a d , m N = { 1 , 2 , } a minimizing sequence, for which
lim m J ( x m , q m , q m , j m ) = inf ( x , q , q , j ) Z a d J ( x , q , q , j ) J * .
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 :
q m s , Σ τ c 1 , q m s , Σ ν c 2 , j m Ω c 3 , u m 1 , Ω c 4 ,
H m H div s + 1 / 2 ( Ω ) c 5 , p m Ω c 6 .
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 q * H s ( Σ τ ) , q * H T s ( Σ ν ) , j * L 2 ( Ω ) 3 , u * H T 1 ( Ω ) , H * H div s + 1 / 2 ( Ω ) , p * L 0 2 ( Ω ) of some subsequences of sequences { q m } , { q m } , { j m } , { u m } , { H m } , { p m } . As usual, one should consider that as m
q m q *   weakly   in   H s ( Σ τ ) , q m q *   weakly   in   H T s ( Σ ν ) , j m j *   weakly   in   L 2 ( Ω ) 3 ,
p m p *   weakly   in   L 2 ( Ω ) , u m u *   weakly   in   H 1 ( Ω ) 3   and   strongly   in   L 4 ( Ω ) 3 ,
H m H *   weakly   in   H div s + 1 / 2 ( Ω ) L 3 ( Ω ) 3   and   strongly   in   L 2 ( Ω ) 3 .
It is clear that u * ( q * , q * , j * ) K , F 2 ( x * ) = 0 , F 3 ( x * ) = 0 , F 4 ( x * , u * ) = 0 and F 5 ( x * , u * ) = 0 where x * ( u * , H * , p * ) . Let us show that F 1 ( x * , u * ) = 0 , i.e., that
ν ( u * , v ) + ν 1 ( rot H * , rot Ψ ) + ( ( u * · ) u * , v ) ( div v , p * ) +
+ κ [ ( rot Ψ × H * , u * ) ( rot H * × H * , v ) ]
= f , v + ν ( j * , rot Ψ ) + ρ 0 1 ( k , Ψ T ) Σ τ ( v , Ψ ) H 0 1 ( Ω ) 3 × V Σ τ ( Ω ) .
To this end, we note that u m , H m , p m and j m satisfy the identity
ν ( u m , v ) + ν 1 ( rot H m , rot Ψ ) + ( ( u m · ) u m , v ) ( div v , p m ) +
+ κ [ ( rot Ψ × H m , u m ) ( rot H m × H m , v ) ] =
= f , v + ν 1 ( j m , rot Ψ ) + ρ 0 1 ( k , Ψ T ) Σ τ ( v , Ψ ) H 0 1 ( Ω ) 3 × V Σ τ ( Ω ) .
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 T 1 ( Ω ) , then arguing as in [3], we easily derive that
( ( u m · ) u m , v ) ( ( u * · ) u * , v ) v H 0 1 ( Ω ) 3   as   m .
We now consider the second nonlinear term ( rot H m × H m , v ) in (43). It is clear that
( rot H m × H m , v ) ( rot H * × H * , v ) =
= ( rot H m × ( H m H * ) , v ) + ( rot ( H m H * ) × H * , v ) .
Let us prove firstly that for the first term in (45), we have
( rot H m × ( H m H * ) , v ) 0 v H 0 1 ( Ω ) 3   as   m .
For this purpose, it is sufficient to prove that for an arbitrary pair, a number ε > 0 and a test function v H 0 1 ( Ω ) 3 , there exists a number M = M ( ε , v ) such that
| ( rot H m × ( H m H * ) , v ) | ε v H 0 1 ( Ω ) 3 , m M .
As H 0 1 ( Ω ) is the closure of D ( Ω ) by norm · 1 , Ω , then for mentioned function v H 0 1 ( Ω ) 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
| ( rot H m × ( H m H * ) , v ) | | ( rot H m × ( H m H * ) , v n ) | +
+ | ( rot H m × ( H m H * ) , v n v ) | m , n N .
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
| ( rot H m × ( H m H * ) , v n v ) | ε / 2 n N , m N .
Let us consider the first term on the right-hand side of (48). Based on Hölder inequality for three functions, we have
| ( rot H m × ( H m H * ) , v n ) | rot H m Ω H m H * Ω v n L ( Ω ) 3 .
Since rot H m Ω H m H div s + 1 / 2 ( Ω ) c 5 and H m H * in L 2 ( Ω ) 3 as m by (41), then there exists a number M = M ( ε , v , N ) such that the estimate | ( rot H m × ( H m H * ) , v n ) | ε / 2 holds for all m M , n = N . 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
| ( ( rot ( H m H * ) × H * ) , v ) | 0 v H 0 1 ( Ω ) 3   as   m .
From the weak convergence of the sequence H m in the space H div s + 1 / 2 ( Ω ) , the weak convergence of rot H m to rot H * in L 2 ( Ω ) 3 follows. As H * × v L 2 ( Ω ) 3 , we have
( rot ( H m H * ) × H * , v ) = ( rot H m rot H * ) , H * × v ) 0 v H 0 1 ( Ω ) 3   as   m 0 .
Therefore, (50) is proved. From (46) and (50), it follows from (45) that
( rot H m × H m , v ) ( rot H * × H * , v ) v H 0 1 ( Ω ) 3   as   m .
Using the analogous scheme, one can show that for the last nonlinear term in (43), we have
( rot Ψ × H m , u m ) ( rot Ψ × H * , u * ) Ψ V Σ τ ( Ω )   as   m .
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 J ( x * , u * ) = J * . □
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 T 2 ( Σ ν ) , is physically the most interesting. Another physically interesting case is the control problem with weak regularizers having the form
J ( x , q ) = ( μ 0 / 2 ) I ( x ) + ( μ 1 / 2 ) q Σ τ 2 + ( μ 2 / 2 ) q Σ ν 2 +
+ ( μ 3 / 2 ) j Ω 2 inf , F ( x , u ) = 0 , ( x , u ) X × K
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:
(jv) μ 0 > 0 ; μ 1 0 , μ 2 0 , μ 3 0 and K is a bounded set in H s ( Σ τ ) × H ˜ T s ( Σ ν ) × L 2 ( Ω ) 3 f o r f i x e d s ( 0 , 1 / 2 ] .
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 a d 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
X * = H T 1 ( Ω ) * × H div s + 1 / 2 ( Ω ) * × L 0 2 ( Ω ) ,
Y * = ( H 0 1 ( Ω ) 3 × V Σ τ ( Ω ) ) × L 0 2 ( Ω ) × H T 1 / 2 ( Σ ) * × H s ( Σ τ ) * × H ˜ T s ( Σ ν ) * ,
which are dual of spaces X and Y defined in (25).
By [13], the derivation of the optimality system uses finding of the Fréchet partial derivative with respect to x of F : ( F 1 , F 2 , F 3 , F 4 , F 5 ) : X × K Y defined in (38). The simple analysis shows that the mentioned partial derivative at any point ( x ^ , u ^ ) ( u ^ , H ^ , p ^ , q ^ , q ^ , j ^ ) X × K is a linear continuous operator F x ( x ^ , u ^ ) : X Y , which associates every element ( w , h , r ) X with the element F x ( x ^ , u ^ ) ( w , h , r ) = ( y ^ 1 , y ^ 2 , y ^ 3 , y ^ 4 , y ^ 5 ) Y . Here the elements y ^ 1 W * = H 1 ( Ω ) 3 × V Σ τ ( Ω ) * , y ^ 2 L 0 2 ( Ω ) , y ^ 3 H T 1 / 2 ( Σ ) , y ^ 4 H s ( Σ τ ) , y ^ 5 H ˜ T s ( Σ ν ) are defined by the triples ( u ^ , H ^ , p ^ ) and ( w , h , r ) from relations
y ^ 1 , ( v , Ψ ) = ν ( w , v ) + ν 1 ( rot h , rot Ψ ) + ( ( w · ) u ^ , v ) +
+ ( ( u ^ · ) w , v ) + κ [ ( rot Ψ × h , u ^ ) + ( rot Ψ × H ^ , w ) ( rot h × H ^ , v )
( rot H ^ × h , v ) ] + b ( v , r ) ( v , Ψ ) H 0 1 ( Ω ) 3 × V Σ τ ( Ω ) ,
y ^ 2 , r = b ( w , r ) r L 0 2 ( Ω ) , y ^ 3 = w | Σ , y ^ 4 = h · n | Σ τ , y ^ 5 = h × n | Σ ν .
By F x ( x ^ , u ^ ) * : Y * X * , we denote the operator adjoint to F x ( x ^ , u ^ ) , which is defined by
F x ( x ^ , u ^ ) * y * , x X * × X = y * , F x ( x ^ , u ^ ) x Y * × Y ( x , y * ) X × Y * .
Following to [13] (Chapter 1), let us introduce an element y * = ( ( ξ , η ) , σ , ζ 1 , ζ 2 , ζ 3 ) Y * having the sense of an adjoint state and define the Lagrangian L : X × K × R + × Y * R where R + = { x R : x 0 } by
L ( x , u , λ 0 , y * ) = λ 0 J ( x , u ) + y * , F ( x , u ) Y * × Y λ 0 J ( x , u ) + F 1 ( x , u ) , ( ξ , η ) +
+ ( F 2 ( x ) , σ ) + ζ 1 , F 3 ( x ) Σ + ζ 2 , F 4 ( x , u ) Σ τ + ζ 3 , F 5 ( x , u ) Σ ν .
Here and below, · , · Σ , · , · Σ τ or · , · Σ ν denote the duality pairings between H T 1 / 2 ( Σ ) and H T 1 / 2 ( Σ ) , between H s ( Σ τ ) and H s ( Σ τ ) or between H T s ( Σ ν ) and H T s ( Σ ν ) , respectively. Based on the results of [13] (Chapter 1), one can prove the following theorem about the justification of the Lagrange principle for problem (40) and the regularity for the Lagrange multiplier ( λ 0 , y * ) .
Theorem 5.
Let, under assumptions (i)–(iii) and (j), (jjj), the element ( x ^ , u ^ ) ( u ^ , H ^ , p ^ , q ^ , q ^ , j ^ ) X × K be a local minimizer for problem (40), and let the cost functional I have the continuous Fréchet derivative with respect to x at x ^ . Then:
(1) there is a nonvanishing Lagrange multiplier ( λ 0 , y * ) = ( λ 0 , ( ξ , η ) , σ , ζ 1 , ζ 2 , ζ 3 ) R + × Y * for which the Euler–Lagrange equation
F x ( x ^ , u ^ ) * y * = λ 0 ( μ 0 / 2 ) I ( x ^ )   i n   X *
for adjoint state y * is satisfied and the minimum principle holds, having the form
L ( x ^ , u ^ , λ 0 , y * ) L ( x ^ , u , λ 0 , y * ) u K .
(2) If, besides, condition (36) holds for each triple u ( q , q , j ) K , then any nonvanishing Lagrange multiplier ( λ 0 , y * ) satisfying (56) is regular, i.e., λ 0 = 1 . Moreover, it is determined uniquely if the value is given.
Proof of Theorem 5.
To prove statement 1 of Theorem 5 by [13] Chapter 1 (see also [15], Chapter 6), it is sufficient to prove that the operator F x ( x ^ , u ^ ) : X Y is a Fredholm operator. By virtue of (54), the operator F ^ F x ( x ^ , u ^ ) : X Y can be written in the form
F ^ = Φ + Φ ^ ( Φ 1 , Φ 2 , Φ 3 , Φ 4 , Φ 5 ) + ( Φ ^ 1 , 0 , 0 , 0 , 0 ) .
Here the operators Φ 1 , Φ 2 , Φ 3 , Φ 4 , Φ 5 : X Y are defined by relations (32) and (33), while the operator Φ ^ 1 : X W * = H 1 ( Ω ) 3 × V Σ τ ( Ω ) * is defined by
Φ ^ 1 ( w , h , r ) , ( v , Ψ ) = ( ( w · ) u ^ , v ) + κ [ ( rot Ψ × h , u ^ ) ( rot H ^ × h , v ) ] .
Using estimates (13)–(15), we deduce that
| ( ( w · ) u ^ , v ) | γ 0 w L 3 ( Ω ) 3 u ^ 1 , Ω v 1 , Ω v , w H T 1 ( Ω ) ,
| ( rot Ψ × h , u ^ ) | rot Ψ Ω h L 3 ( Ω ) 3 u ^ L 6 ( Ω ) 3 Ψ H div s + 1 / 2 ( Ω ) , h H div s + 1 / 2 ( Ω ) ,
| ( rot H ^ × h , v ) | rot H ^ Ω h L 3 ( Ω ) 3 v L 6 ( Ω ) 3 v H T 1 ( Ω ) , h H div s + 1 / 2 ( Ω ) .
Since the space H 1 ( Ω ) 3 and the space H div s + 1 / 2 ( Ω ) for any s > 0 are continuously and compactly embedded into the space L 3 ( Ω ) 3 , then estimates (60)–(62) imply that the operator Φ ^ = ( Φ ^ 1 , 0 , 0 , 0 , 0 ) : X Y where Φ ^ 1 is defined by (59) is continuous and compact. In addition, it follows from Theorem 1 that the operator Φ : X Y , defined by relations (32) and (33) is an isomorphism. This means that the operator F x ( x ^ , u ^ ) is Fredholm as the sum of the isomorphism Φ and the continuous compact operator Φ ^ . Therefore, statement 1 about existence of nonzero Lagrange multiplier ( λ 0 , y * ) is proved.
It remains to prove the regularity of the multiplier ( λ 0 , y * ) , i.e., that λ 0 0 under conditions (36). Arguing as in paper [3], one can verify that (56) is equivalent to three identities with respect to the adjoint state ( ξ , η , σ , ζ 1 , ζ 2 , ζ 3 ) having the form:
ν ( w , ξ ) + ν 1 ( rot h , rot η ) + ( ( u ^ · ) w , ξ ) + ( ( w · ) u ^ , ξ ) +
+ κ [ ( rot η × H ^ , w ) + ( rot η × h , u ^ ) ]
κ [ ( rot H ^ × h , ξ ) + ( rot h × H ^ , ξ ) ] ( div w , σ ) + ζ 1 , w Σ + ζ 2 , h · n Σ τ +
+ ζ 3 , h × n Σ ν =
= λ 0 ( μ 0 / 2 ) I u ( x ^ ) , w λ 0 ( μ 0 / 2 ) I H ( x ^ ) , h ( w , h ) H T 1 ( Ω ) × H div s + 1 / 2 ( Ω ) ,
( div ξ , r ) = λ 0 ( μ 0 / 2 ) ( I p ( x ^ ) , r ) r L 0 2 ( Ω ) .
The mentioned statement about the regularity of the multiplier ( λ 0 , y * ) is equivalent to the statement about the nonexistence of nonvanishing solutions of the systems (63), (64) at λ 0 = 0 , in which the elements x ^ = ( u ^ , H ^ , p ^ ) and u ^ = ( q ^ , q ^ , j ^ ) are connected by the relation F ( x ^ , u ^ ) = 0 . To prove this, one should denote by y * = ( ξ , η , σ , ζ 1 , ζ 2 , ζ 3 ) Y * an arbitrary solution of problems (63), (64) at λ 0 = 0 . Setting w = ξ , h = η and λ 0 = 0 in (63), (64) and using (16), we have
ν ( ξ , ξ ) + ν 1 ( rot η , rot η ) + ( ( ξ · ) u ^ , ξ ) + κ [ ( rot η × η , u ^ ) ( rot H ^ × η , ξ ) ] = 0 ,
div ξ = 0   in   Ω .
Arguing as in [8], i.e., applying inequalities (8), (10), (13)–(15) for estimating all terms in (65) and using estimates (35) at u = u ^ , H = H ^ , we derive the following inequality from (65):
( δ 0 ν γ 0 M u ^ γ 1 ( κ / 2 ) M H ^ ) ξ 1 , Ω 2 + ( δ 1 ν m γ 1 M u ^ γ 1 ( κ / 2 ) M H ^ ) κ η 1 , Ω 2 0 .
From (66), it follows under smallness conditions (36) at u = u ^ , H = H ^ that ξ = 0 and η = 0 in Ω .
Setting ξ = 0 , η = 0 in (63) at λ 0 = 0 , we have
( div w , σ ) + ζ 1 , w Σ + ζ 2 , h · n Σ τ + ζ 3 , h × n Σ ν = 0 ( w , h ) H T 1 ( Ω ) × H div s + 1 / 2 ( Ω )
which in the case h = 0 , transforms to
( div w , σ ) + ζ 1 , w Σ = 0 w H T 1 ( Ω ) .
From (68), it follows that
( div w , σ ) = 0 w H 0 1 ( Ω ) 3 .
By inf-sup condition (17), this identity holds if σ = 0 . Setting σ = 0 in (68), we obtain ζ 1 , w Σ = 0 for all w H T 1 ( Ω ) . This means that ζ 1 = 0 in H T 1 / 2 ( Σ ) .
Setting σ = 0 , ζ 1 = 0 in (67), we arrive at
ζ 2 , h · n Σ τ + ζ 3 , h × n Σ ν = 0 h H div s + 1 / 2 ( Ω ) .
Choosing h H div s + 1 / 2 ( Ω , Σ ν ) in (70), where H div s + 1 / 2 ( Ω , Σ ν ) is defined in (24), we obtain ζ 2 , h · n Σ τ = 0 for all h H div s + 1 / 2 ( Ω , Σ ν ) . This means that ζ 2 = 0 in H s ( Σ τ ) . Analogously, we deduce that ζ 3 = 0 in H T s ( Σ ν ) . Thus, we obtained y * = 0 , and the regularity of the Lagrage multiplier is proved. Concerning the uniqueness property of the Lagrange multiplier ( 1 , y i * ) , it is a consequence of Fredholm property of the linear operator F x ( x ^ , u ^ ) : X Y . □
We note that the Lagrangian L defined in (55) is a continuously differentiable function of controls q K 1 , q K 2 and j K 3 , and its partial derivatives L q ( x ^ , u ^ , λ 0 , y * ) with respect to q, L q ( x ^ , u ^ , λ 0 , y * ) with respect to q and L j ( x ^ , u ^ , λ 0 , y * ) with respect to j in any point ( x ^ , u ^ , λ 0 , y * ) are determined by
L q ( x ^ , u ^ , λ 0 , y * ) , q = λ 0 μ 1 ( q ^ , q ) s , Σ τ ζ 2 , q Σ τ q K 1 ,
L q ( x ^ , u ^ , λ 0 , y * ) , q = λ 0 μ 2 ( q ^ , q ) s , Σ ν ζ 3 , q Σ ν q K 2 ,
L j ( x ^ , u ^ , λ 0 , y * ) , j = λ 0 μ 3 ( j ^ , j ) ν 1 ( j , rot η ) j K 3 .
As the triple u ^ ( q ^ , q ^ , j ^ ) is a minimizer of the function L ( x ^ , · , λ 0 , y * ) on a closed convex set K = K 1 × K 2 × K 3 by (57), the following inequalities hold (see, e.g., [52]):
L q ( x ^ , u ^ , λ 0 , y * ) , q q ^ λ 0 μ 1 ( q ^ , q q ^ ) s , Σ τ ζ 2 , q q ^ Σ τ 0 q K 1 ,
L q ( x ^ , u ^ , λ 0 , y * ) , q q ^ λ 0 μ 2 ( q ^ , q q ^ ) s , Σ ν ζ 3 , q q ^ Σ ν 0 q K 2 ,
L j ( x ^ , u ^ , λ 0 , y * ) , j j ^ λ 0 μ 3 ( j ^ , j j ^ ) ν 1 ( j j ^ , rot η ) 0 j K 3 .
Identities (63) and (64), the variational inequalities (71)–(73) and the operator constraint (state equation) (39), which is equivalent to the weak forms (19) and (20) of Problem 1, constitute the optimality system for our control problem (40). It describes the first-order necessary conditions of optimality.
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
[ ξ ] = [ u ] = L 0 / T 0 , [ η ] = [ H ] = I 0 / L 0 , [ σ ] = [ p ] = L 0 2 / T 0 2 .
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”.
Remark 2.
It follows from conditions ξ H 0 1 ( Ω ) 3 , η V Σ τ ( Ω ) and (64) that the adjoint velocity ξ and adjoint magnetic field η possess properties
ξ | Σ = 0 , η · n | Σ τ = 0 , η × n | Σ ν = 0 , div η = 0   and   div ξ = λ 0 ( μ 0 / 2 ) χ Q I p ( x ^ )   in   Ω .
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
J ˜ ( x , u ) = μ 0 2 I ˜ ( x ) + μ 1 2 q s , Σ τ 2 + μ 2 2 q s , Σ ν 2 + μ 3 2 j Ω 2 inf , F ( x , u ) = 0 , ( x , u ) X × K .
Here I ˜ 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 :
u i 1 , Ω M u 0 = sup u K M u ( u ) , H i H div s + 1 / 2 ( Ω ) M H 0 = sup u K M H ( u ) ,
p i Ω M p 0 = sup u K M p ( u ) .
Let the values M u 0 and M H 0 be such that
( γ 0 / δ 0 ν ) M u 0 + ( γ 1 / δ 0 ν ) ( κ / 2 ) M H 0 < 1 / 2 ,
( γ 1 / δ 1 ν m ) M u 0 + ( γ 1 / δ 1 ν m ) ( κ / 2 ) M H 0 < 1 / 2 .
To make conditions (78) more illustrative and to simplify the subsequent presentation, we define the parameters
R e = γ 0 M u 0 / ( δ 0 ν ) , R m = γ 1 M u 0 / ( δ 1 ν m ) , H a = γ 1 κ M H 0 / ( δ 0 ν ) , P m = δ 0 ν / ( δ 1 ν m ) .
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 R e , R m , H a 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 u 0 and M H 0 entering into (78). In order to determine the dimensions of R e , R m , H a 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:
u Ω 2 = Ω u 2 d Ω , | u | 1 , Ω 2 = Ω | u | 2 d Ω , u 1 , Ω 2 = l 2 u Ω 2 + | u | 1 , Ω 2 , [ l ] = L 0 .
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
[ u Ω ] = [ u ] L 0 3 / 2   and   [ | u | 1 , Ω ] = [ u 1 , Ω ] = [ u ] L 0 1 / 2 .
We recall also (see, e.g., [15] (p. 272)) that
[ g ] = [ u ] = L 0 / T 0 , [ q ] = [ q ] = [ H ] = I 0 / L 0 , [ j ] = I 0 / L 0 2 ,
[ ν m ] = [ ν ] = L 0 2 / T 0 , [ κ ] = L 0 4 / T 0 2 I 0 2 .
Combining this with (10)–(15), (17) and (74) yields
[ δ i ] = 1 , [ γ i ] = L 0 1 / 2 , [ β ] = 1 ,
[ M u 0 ] = [ u 1 , Ω ] = [ ξ 1 , Ω ] = [ g 1 / 2 , Σ ] = L 0 3 / 2 / T 0 , [ M p 0 ] = L 0 7 / 2 / T 0 2 ,
[ M H 0 ] = [ H H div s + 1 / 2 ( Ω ) ] = [ η 1 , Ω ] = [ q Σ τ ] = [ q ] Σ ν = [ j Ω ] = I 0 / L 0 1 / 2 ,
[ κ H H div s + 1 / 2 ( Ω ) ] = L 0 3 / 2 / T 0 = [ M u 0 ] .
Using (81), we have that
[ R e ] = [ R m ] = [ H a ] = [ P m ] = 1 ,
i.e., all the parameters R e , R m , H a and P m defined in (79) are dimensionless. Since parameters R e and R m are connected by relation R m = ( γ 1 / γ 0 ) P m R e , we can rewrite conditions (78) in the following form containing only three nondimensional parameters R e , H a and P m :
R e + ( 1 / 2 ) H a < 1 / 2 , ( γ 1 / γ 0 ) P m R e + ( 1 / 2 ) P m H a < 1 / 2 .
Denote, by ( 1 , y i * ) ( 1 , ξ i , η i , σ i , ζ 1 ( i ) , ζ 2 ( i ) , ζ 3 ( i ) ) , 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
ν ( w , ξ i ) + ν 1 ( rot h , rot η i ) + ( ( u i · ) w , ξ i ) + ( ( w · ) u i , ξ i ) + κ [ ( rot η i × H i , w )
+ ( rot η i × h , u i ) ] κ [ ( rot H i × h , ξ i ) + ( rot h × H i , ξ i ) ] ( div w , σ i ) + ζ 1 ( i ) , w Σ +
+ ζ 2 ( i ) , h · n Σ τ + ζ 3 ( i ) , h × n Σ ν =
= ( μ 0 / 2 ) ( I i ) u ( x i ) , w ( μ 0 / 2 ) ( I i ) H ( x i ) , h ( w , h ) H T 1 ( Ω ) × H div s + 1 / 2 ( Ω ) , i = 1 , 2 ,
( div ξ i , r ) = ( μ 0 / 2 ) ( ( I i ) p ( x i ) , r ) r L 0 2 ( Ω ) , i = 1 , 2 .
Here we renamed I = I 1 , I ˜ = I 2 . We define the following differences:
u = u 1 u 2 , H = H 1 H 2 , p = p 1 p 2 , q = q 1 q 2 , q = q 1 q 2 , j = j 1 j 2 ,
ξ = ξ 1 ξ 2 , σ = σ 1 σ 2 , ζ 1 = ζ 1 ( 1 ) ζ 1 ( 2 ) , ζ 2 = ζ 2 ( 1 ) ζ 2 ( 2 ) , ζ 3 = ζ 3 ( 1 ) ζ 3 ( 2 ) .
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.
Subtract relations (19), (20), written for u 2 , H 2 , p 2 , u 2 , from (19), (20) for u 1 , H 1 , p 1 , u 1 . Using (85), we obtain
ν ( u , v ) + ν 1 ( rot H , rot Ψ ) + [ ( ( u · ) u 1 , v ) + ( ( u 2 · ) u , v ) ] + κ [ ( rot Ψ × H , u 1 ) +
+ ( rot Ψ × H 2 , u ) ] κ [ ( rot H 1 × H , v ) + ( rot H × H 2 , v ) ] ( div v , p ) =
= ν 1 ( j , rot Ψ ) ( v , Ψ ) H 0 1 ( Ω ) 3 × V Σ τ ( Ω ) ,
div u = 0   in   Ω , u | Σ = 0 , H · n | Σ τ = q , H × n | Σ ν = q .
Set j = j 2 in the inequality (73) under λ 0 = 1 , written at j ^ = j 1 , η = η 1 , and then we set j = j 1 in (73) under λ 0 = 1 , written at j ^ = j 2 , η = η 2 . We obtain μ 3 ( j 1 , j ) + ν 1 ( j , rot η 1 ) 0 and μ 3 ( j 2 , j ) ν 1 ( j , rot η 2 ) 0 . Adding these inequalities yields the relation
ν 1 ( j , rot η ) μ 3 j Ω 2 .
In the same manner, we derive from (71) and (72) the inequalities
ζ 2 , q Σ τ μ 1 q s , Σ τ 2 , ζ 3 , q Σ ν μ 2 q s , Σ ν 2 .
Let us subtract identities (83), (84) for i = 2 from (83), (84) for i = 1 and set w = u , h = H , r = p . Adding the results and using conditions u | Σ = 0 , H · n | Σ τ = q , H × n | Σ ν = q , we obtain
ν ( u , ξ ) + ν 1 ( rot H , rot η ) + ( ( u 1 · ) u , ξ ) + 2 ( ( u · ) u , ξ 2 ) + ( ( u · ) u 1 , ξ ) +
+ κ [ ( rot η × H 1 , u ) + 2 ( rot η 2 × H , u ) + ( rot η × H , u 1 ) ]
κ [ ( rot H 1 × H , ξ ) + 2 ( rot H × H , ξ 2 ) + ( rot H × H 1 , ξ ) ] ( div ξ , p ) =
= ζ 2 , q Σ τ ζ 3 , q Σ ν ( μ 0 / 2 ) I u ( x 1 ) I ˜ u ( x 2 ) , u
( μ 0 / 2 ) I H ( x 1 ) I ˜ H ( x 2 ) , H ( μ 0 / 2 ) I p ( x 1 ) I ˜ p ( x 2 ) , p .
We now set v = ξ , Ψ = η in (86) and subtract from (90). Using (88), (89) and relations
2 ( ( u · ) u , ξ 2 ) + ( ( u 1 · ) u , ξ ) ( ( u 2 · ) u , ξ ) = 2 ( ( u · ) u , ξ 2 ) + ( ( u · ) u , ξ ) =
= ( ( u · ) u , ξ 1 + ξ 2 ) ,
2 ( rot H × H , ξ 2 ) + ( rot H × H 2 , ξ ) ( rot H × H 1 , ξ ) = ( rot H × H , ξ 1 + ξ 2 ) ,
2 ( rot η 2 × H , u ) + ( rot η × H 1 , u ) ( rot η × H 2 , u ) = ( rot ( η 1 + η 2 ) × H , u ) ,
we arrive at
( ( u · ) u , ξ 1 + ξ 2 ) κ ( rot H × H , ξ 1 + ξ 2 ) + κ ( rot ( η 1 + η 2 ) × H , u ) +
+ ( μ 0 / 2 ) I u ( x 1 ) I ˜ u ( x 2 ) , u + ( μ 0 / 2 ) I H ( x 1 ) I ˜ H ( x 2 ) , H + ( μ 0 / 2 ) ( I p ( x 1 ) I ˜ p ( x 2 ) , p )
μ 1 q s , Σ τ 2 μ 2 q s , Σ ν 2 μ 3 j Ω 2 .
This inequality alone will take the important role in Section 5 when studying the local uniqueness and stability for optimal solutions. It is appropriate to formulate this result as the next theorem.
Theorem 6.
Let, under assumptions of Theorem 3 for cost functionals I and I ˜ and under (82), pairs ( x 1 , u 1 ) = ( u 1 , H 1 , p 1 , q 1 , q 1 , j 1 ) X × K and ( x 2 , u 2 ) = ( u 2 , H 2 , p 2 , q 2 , q 2 , j 2 ) X × K be solutions to (40) and (76), respectively. Let y i * = ( ξ i , η i , σ i , ζ 1 ( i ) , ζ 2 ( i ) , ζ 3 ( i ) ) Y * , i = 1 , 2 , be adjoint states which correspond to these solutions. Then the relation (91) for differences u , H , p , q , q , j defined in (85) holds.
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 rot H 0 = 0 , we obtain
ν ( u , u ) + ν 1 ( rot H ˜ , rot H ˜ ) = ( ( u · ) u 1 , u ) κ ( rot H ˜ × H ˜ , u 1 ) κ ( rot H ˜ × H 0 , u 1 )
κ ( rot H 1 × H ˜ , u ) κ ( rot H 1 × H 0 , u ) + ν 1 ( j , rot H ˜ ) .
Using estimates (10)–(15), (77) and setting ν 1 = ν m κ , we obtain from (92) that
δ 0 ν u 1 , Ω 2 + δ 1 ν m κ H ˜ 1 , Ω 2 γ 0 M u 0 u 1 , Ω 2 + κ γ 1 M u 0 H ˜ 1 , Ω 2
+ κ γ 1 M H 0 u 1 , Ω H ˜ 1 , Ω + κ γ 1 M H 0 u 1 , Ω H 0 s + 1 / 2 , Ω
+ κ γ 1 M u 0 H ˜ 1 , Ω H 0 s + 1 / 2 , Ω + C 1 ν m κ j Ω H ˜ 1 , Ω .
Applying (79) and Young’s inequality | a b | ε a 2 / 2 + b 2 / ( 2 ε ) or 2 a b ε a 2 / 2 + 2 b 2 / ε , ε = const > 0 , at ε = δ 0 ν / 2 , ε = δ 1 ν m κ / 4 or ε = 1 , we consequently derive that
κ γ 1 M H 0 u 1 , Ω H ˜ 1 , Ω ( κ / 2 ) γ 1 ( u 1 , Ω 2 + κ H ˜ 1 , Ω 2 ) ,
κ γ 1 M H 0 H 0 s + 1 / 2 , Ω u 1 , Ω ε 2 u 1 , Ω 2 + ( γ 1 κ M H 0 H 0 s + 1 / 2 , Ω ) 2 2 ε =
= δ 0 ν 4 u 1 , Ω 2 + κ δ 0 ν H a 2 H 0 s + 1 / 2 , Ω 2 ,
κ γ 1 M H 0 H 0 s + 1 / 2 , Ω H ˜ 1 , Ω ε 2 H ˜ 1 , Ω 2 + ( γ 1 κ M u 0 ) 2 2 ε H 0 s + 1 / 2 , Ω 2 ) =
= ( κ δ 1 ν m / 8 ) H ˜ 1 , Ω 2 + 2 κ δ 0 ν P m ( γ 1 / γ 0 ) 2 R e H 0 s + 1 / 2 , Ω 2 ,
C 1 κ ν m j Ω H ˜ 1 , Ω 2 κ ν m C 1 2 δ 1 1 j Ω 2 + δ 1 κ ν m 8 H ˜ 1 , Ω .
Taking into account (94)–(97), we obtain from (93) that
( δ 0 ν γ 0 M u 0 + γ 1 ( κ / 2 ) M H 0 ) u 1 , Ω 2 + ( δ 1 ν m γ 1 M H 0 + γ 1 ( κ / 2 ) M H 0 ) κ H ˜ 1 , Ω 2
( δ 0 ν / 4 ) u 1 , Ω 2 + ( δ 1 ν m κ / 4 ) H ˜ 1 , Ω 2 + κ δ 0 ν ( H a 2 +
+ 2 P m ( γ 1 / γ 0 ) 2 R 2 ) H 0 s + 1 / 2 , Ω 2 + 2 κ ν m C 1 2 δ 1 1 j Ω 2 .
It follows from (78) that
( δ 0 ν / 2 ) < δ 0 ν γ 0 M u 0 γ 1 ( κ / 2 ) M H 0 , ( δ 1 ν m / 2 ) < δ 1 ν m γ 1 M u 0 γ 1 ( κ / 2 ) M H 0 .
Using (99), we derive from (98) that
δ 0 ν u 1 , Ω 2 + δ 1 ν m κ H ˜ 1 , Ω 2 2 κ δ 0 ν R 2 H 0 s + 1 / 2 , Ω 2 + 8 κ ν m C 1 2 δ 1 1 j Ω 2 )
where
R = 2 H a 2 + 2 P m ( γ 1 / γ 0 ) 2 R e 2 .
From (100), we conclude that
u 1 , Ω ( 2 κ R 2 H 0 s + 1 / 2 , Ω 2 + 8 κ ν m C 1 2 δ 1 1 δ 0 ν j Ω 2 ) 1 / 2
κ R H 0 s + 1 / 2 , Ω + 2 C 1 δ 1 κ P m j Ω ,
H ˜ 1 , Ω ( 2 δ 0 ν δ 1 ν m R 2 H 0 s + 1 / 2 , Ω 2 + 8 C 1 2 δ 1 2 j Ω 2 ) 1 / 2 P m R H 0 s + 1 / 2 , Ω + 2 C 1 δ 1 j Ω .
Since H = H 0 + H ˜ , from (102), (103), using (11), we obtain the needed estimates for differences u and H :
u 1 , Ω C u κ ( q Σ τ + q Σ ν + j Ω ) ,
H H div s + 1 / 2 ( Ω ) C H ( q Σ τ + q Σ ν + j Ω ) .
Here C u and C H are dimensionless constants defined by
C u = max { C Σ R , 2 C 1 δ 1 P m } , C H = max { C Σ ( 1 + C 1 P m R ) , 2 C 1 2 δ 1 } .
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 v 0 H 0 1 ( Ω ) 3 , v 0 0 such that
( div v 0 , p ) β 0 v 0 1 , Ω p Ω , β 0 = ( β δ ) > 0 .
Setting v = v 0 , Ψ = 0 in (86), we obtain
ν ( u , v 0 ) + [ ( ( u · ) u 1 , v 0 ) + ( ( u 2 · ) u , v 0 ) ] κ [ ( rot H 1 × H , v 0 ) +
+ ( rot H × H 2 , v 0 ) ] ( div v 0 , p ) = 0 .
Using the previous estimate (107) for ( div v 0 , p ) and (10), (13)–(15) from (108), we deduce that
β 0 v 0 1 , Ω p Ω ( div v 0 , p ) ν v 0 1 , Ω u 1 , Ω +
+ 2 γ 0 M u 0 v 0 1 , Ω u 1 , Ω + 2 κ γ 1 M H 0 v 0 1 , Ω H H div s + 1 / 2 ( Ω ) .
Dividing (109) by v 0 1 , Ω 0 and using (104), (105) yields
p Ω β 0 1 [ ( ν + 2 γ 0 M u 0 ) u 1 , Ω + 2 κ γ 1 M H 0 H H div s + 1 / 2 ( Ω ) ] =
= β 0 1 δ 0 ν [ ( δ 0 1 + 2 R e ) u 1 , Ω + 2 κ H a H H div s + 1 / 2 ( Ω ) ] .
Taking into account (104) and (105), we obtain the following estimate:
p Ω C p ν κ ( q Σ τ + q Σ ν + j Ω ) .
Here C p is a dimensionless constant defined by
C p δ 0 β 0 1 [ ( δ 0 1 + 2 R e ) C u + 2 H a C H ] .
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 R e , H a , P m defined in (79) and three dimensionless constants C u , C H , C p defined in (106), (112). These constants depend on R e , H a 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.

5. Analysis of Uniqueness and Stability for Solutions to Control Problems

We begin with consideration of the case when I = I 1 in (40), i.e., we consider control problem
J ( x , u ) μ 0 2 u u d Q 2 + μ 1 2 q s , Σ τ 2 + μ 2 2 q s , Σ ν 2 + μ 3 2 j Ω 2 inf ,
F ( x , u ) = 0 , ( x , u ) X × K .
Here x = ( u , H , p ) , u = ( q , q , j ) . Denote by ( x 1 , u 1 ) ( u 1 , H 1 , p 1 , q 1 , q 1 , j 1 ) a solution to problem (113) that corresponds to the function u d u d ( 1 ) L 2 ( Q ) 3 . By ( x 2 , u 2 ) ( u 2 , H 2 , p 2 , q 2 , q 2 , j 2 ) , we denote a solution to problem (113) that corresponds to another function u ˜ d u d ( 2 ) L 2 ( Q ) 3 .
We define dimensionless parameter R e 0 (Reynolds number for the data u d ( i ) ) by
R e 0 = ( γ 0 / δ 0 ν l ) max ( u d ( 1 ) Q , u d ( 2 ) Q )
where l is a dimensional constant defined in Section 2. We assume that the data for problem (113) or parameters μ 0 , μ 1 , μ 2 and μ 3 are such that the following condition with some sufficiently small ε > 0 takes place:
min { ( 1 ε ) μ 1 , ( 1 ε ) μ 2 , ( 1 ε ) μ 3 } 12 μ 0 κ γ ( R e + R e 0 ) ( γ 0 C u 2 + γ 1 C H 2 + γ 1 P m C u C H ) .
Here C u , C H are dimensionless constants defined in (106).
Lemma 6.
Let, under assumptions (i)–(iii), (j) and (82), a pair ( x i , u i ) = ( u i , H i , p i , q i , q i , j i ) be a solution to problem (113) corresponding to function u d ( i ) L 2 ( Q ) 3 , i = 1 , 2 , where Q Ω is an arbitrary nonempty open subset. Assume that condition (115) is satisfied. Then the following estimate for u u 1 u 2 holds:
u 1 u 2 Q Δ u d ( 1 ) u d ( 2 ) Q .
Proof of Lemma 6.
Setting u d = u d ( 1 ) u d ( 2 ) , in addition to (85), we have that ( I 1 ) p = 0 , ( I 1 ) H = 0 ,
( I 1 ) u ( x i ) , w = 2 ( u i u d ( i ) , w ) Q , ( I 1 ) u ( x 1 ) ( I ˜ 1 ) u ( x 2 ) , u =
= 2 ( u u d , u ) Q = 2 ( u Q 2 ( u , u d ) Q ) .
In view of (117), identities (83), (84) for adjoint states y i * = ( ξ i , η i , σ i , ζ 1 ( i ) , ζ 2 ( i ) , ζ 3 ( i ) ) , i = 1 , 2 , corresponding to solutions ( u i , H i , p i , u i ) and the main inequality (91) for differences u , H , p, q, q , j defined in (85), take by Remark 2 the form
ν ( w , ξ i ) + ν 1 ( rot h , rot η i ) + ( ( u i · ) w , ξ i ) + ( ( w · ) u i , ξ i ) + κ [ ( rot η i × H i , w ) +
+ ( rot η i × h , u i ) ] κ [ ( rot H i × h , ξ i ) + ( rot h × H i , ξ i ) ] ( div w , σ i ) + ζ 1 ( i ) , w Σ +
+ ζ 2 ( i ) , h · n Σ τ + ζ 3 ( i ) , h × n Σ ν =
= μ 0 ( u i u d ( i ) , w ) Q ( w , h ) H T 1 ( Ω ) × H div s + 1 / 2 ( Ω ) ; ξ i V , i = 1 , 2 ,
( ( u · ) u , ξ 1 + ξ 2 ) κ ( rot H × H , ξ 1 + ξ 2 ) + κ ( rot ( η 1 + η 2 ) × H , u ) +
+ μ 0 ( u Q 2 ( u , u d ) Q ) μ 1 q s , Σ τ 2 μ 2 q s , Σ ν 2 μ 3 j Ω 2 .
Our nearest purpose is to estimate adjoint state variables ξ i and η i via ( R e + R e 0 ) . To this end, we set w = ξ i , h = η i in (118). Using (16) and conditions div ξ i = 0 in Ω , ξ i | Σ = 0 , η i · n | Σ τ = 0 and η i × n | Σ ν = 0 , we obtain
ν ( ξ i , ξ i ) + ν 1 ( rot η i , rot η i ) + ( ( ξ i · ) u i , ξ i ) + κ ( rot η i × η i , u i )
κ ( rot H i × η i , ξ i ) = μ 0 ( u i u d ( i ) , ξ i ) Q .
Taking into account estimates (10)–(15), (77) and (79), (80), (114), we have
( ξ i , ξ i ) δ 0 ξ i 1 , Ω 2 , | ( ( ξ i · ) u i , ξ i ) | γ 0 u i 1 , Ω ξ i 1 , Ω 2 γ 0 M u 0 ξ i 1 , Ω 2 ,
κ | ( rot H i × η i , ξ i ) | γ 1 κ M H 0 η i 1 , Ω ξ i 1 , Ω γ 1 ( κ / 2 ) M H 0 ( ξ i 1 , Ω 2 + κ η i 1 , Ω 2 ) ,
( rot η i , rot η i ) δ 1 η i 1 , Ω 2 , | ( rot η i × η i , u i ) | γ 1 M u 0 η i 1 , Ω 2 ,
u i u d ( i ) Q u i Q + u d ( i ) Q l M u 0 + u d ( i ) Q δ 0 ν l γ 0 1 ( R e + R e 0 ) .
In virtue of (121)–(124) and (99), we infer from (120) that
( δ 0 ν / 2 ) ξ i 1 , Ω 2 + ( δ 1 ν m / 2 ) κ η i 1 , Ω 2 [ δ 0 ν γ 0 M u 0 γ 1 ( κ / 2 ) M H 0 ] ξ i 1 , Ω 2 +
+ [ δ 1 ν m γ 1 M u 0 γ 1 ( κ / 2 ) M H 0 ] κ η i 1 , Ω 2
μ 0 l ( l M u 0 + u d ( i ) Q ) ξ i 1 , Ω μ 0 δ 0 ν γ ( R e + R e 0 ) ξ i 1 , Ω
where γ = l 2 γ 0 1 . Using (79) and (124), we conclude from the last inequality that
ξ i 1 , Ω 2 μ 0 γ ( R e + R e 0 ) , κ η i 1 , Ω 2 μ 0 γ P m ( R e + R e 0 ) , i = 1 , 2 .
Using (13)–(15), (104), (105), (115) and (125), we have
| ( ( u · ) u , ξ 1 + ξ 2 ) κ ( rot H × H , ξ 1 + ξ 2 ) + κ ( rot ( η 1 + η 2 ) × H , u ) |
γ 0 u 1 , Ω 2 ( ξ 1 1 , Ω + ξ 2 1 , Ω ) + γ 1 κ H H div s + 1 / 2 ( Ω ) 2 ( ξ 1 1 , Ω + ξ 2 1 , Ω ) +
+ γ 1 κ H H div s + 1 / 2 ( Ω ) u 1 , Ω ( η 1 1 , Ω + η 2 1 , Ω )
12 μ 0 κ γ ( R e + R e 0 ) ( γ 0 C u 2 + γ 1 C H 2 + γ 1 P m C u C H ) ( q s , Σ τ 2 + q s , Σ ν 2 + j Ω 2
( 1 ε ) μ 1 q s , Σ τ 2 + ( 1 ε ) μ 2 q s , Σ ν 2 + ( 1 ε ) μ 3 j Ω 2 .
Taking into account (126) from (119), we arrive at
μ 0 ( u Q 2 ( u , u d ) Q ) ( ( u · ) u , ξ 1 + ξ 2 ) + κ ( rot H × H , ξ 1 + ξ 2 )
κ ( rot ( η 1 + η 2 ) × H , u ) μ 1 q Σ τ 2 μ 2 q Σ ν 2 μ 3 j Ω 2
ε μ 1 q s , Σ τ 2 ε μ 2 q s , Σ ν 2 ε μ 3 j Ω 2 .
Omitting the nonpositive term ε μ 1 q s , Σ τ 2 ε μ 2 q s , Σ ν 2 ε μ 3 j Ω 2 , we derive from (127) that u Q 2 u Q u d Q 0 or
u Q u d Q .
Since u = u 1 u 2 , u d = u d ( 1 ) u d ( 2 ) , Lemma 6 is proved. □
If Q = Ω , the estimate (116) is the stability estimate in L 2 ( Ω ) 3 norm for the component u ^ of the solution ( u ^ , H ^ , p ^ , q ^ , q ^ , j ^ ) of problem (113) with respect to small disturbances of function u d L 2 ( Ω ) 3 in the norm of L 2 ( Ω ) 3 . The same estimate was obtained in viscous hydrodynamics [16]. Additionally, if u d ( 1 ) = u d ( 2 ) we conclude from (116) that u 1 = u 2 in Q. This yields together with (127) and (86) that q 1 = q 2 , q 1 = q 2 , j 1 = j 2 . In turn, it follows from this fact and (104), (105), (111) that u 1 = u 2 , H 1 = H 2 , p 1 = p 2 in Ω . The latter is equivalent to the uniqueness for the solution of (113).
We cannot prove the estimates for differences H = H 1 H 2 and p = p 1 p 2 , which are analogous to (116). Based on (128), however, one can obtain coarser stability estimates for all differences u , H and p even in situations in which Q Ω , i.e., Q is only a part of Ω . Indeed, let us consider inequality (127). Using (128), we deduce from (127) that
ε μ 1 q s , Σ τ 2 + ε μ 2 q s , Σ ν 2 + ε μ 3 j Ω 2 μ 0 ( u Q 2 + u Q u d Q ) μ 0 u d Q 2 = μ 0 Δ 2
where Δ is defined in (116). From (104), (105), (111) and (129), we arrive at
q 1 q 2 s , Σ τ μ 0 / ε μ 1 Δ , q 1 q 2 s , Σ ν μ 0 / ε μ 2 Δ , j 1 j 2 Ω μ 0 / ε μ 3 Δ ,
u 1 u 2 1 , Ω C u κ ( μ 0 / ε μ 1 + μ 0 / ε μ 2 + μ 0 / ε μ 3 ) Δ ,
H 1 H 2 H div s + 1 / 2 ( Ω ) C H ( μ 0 / ε μ 1 + μ 0 / ε μ 2 + μ 0 / ε μ 3 ) Δ ,
p 1 p 2 Ω C p ν κ ( μ 0 / ε μ 1 + μ 0 / ε μ 2 + μ 0 / ε μ 3 ) Δ .
Let us describe the obtained result as the theorem:
Theorem 7.
Let parameter R e 0 be defined by relation (114) and let assumptions of Lemma 6 be fulfilled. Then the stability estimates (130) for problem (113) hold where Δ = u d ( 1 ) u d ( 2 ) Q .
We now study the uniqueness and stability for solutions to the control problem (40) in the case when I = I 2 , i.e., we consider control problem
J ( x , u ) μ 0 2 H H d Q 2 + μ 1 2 q s , Σ τ 2 + μ 2 2 q s , Σ ν 2 + μ 2 2 j Ω 2 inf ,
F ( x , u ) = 0 , ( x , u ) X × K .
Here, as usual, x = ( u , H , p ) , u = ( q , q , j ) . Denote by ( x 1 , u 1 ) ( u 1 , H 1 , p 1 , q 1 , q 1 , j 1 ) a solution to the problem (131) that corresponds to the function H d H d ( 1 ) L 2 ( Q ) 3 . By ( x 2 , u 2 ) ( u 2 , H 2 , p 2 , q 2 , q 2 , j 2 ) , we denote the solution to problem (131) that corresponds to another function H ˜ d H d ( 2 ) L 2 ( Q ) 3 .
Let us define a parameter γ and dimensionless Hartman number H a 0 for the data H d ( i ) by
γ = l 2 γ 1 1 , H a 0 = ( γ 1 κ / δ 0 ν l ) max ( H d ( 1 ) Q , H d ( 2 ) Q ) .
We assume that the data for problem (131) or parameters μ 0 , μ 1 , μ 2 and μ 3 are such that the following condition with some sufficiently small ε > 0 takes place:
min { ( 1 ε ) μ 1 , ( 1 ε ) μ 2 , ( 1 ε ) μ 3 }
12 μ 0 γ P m ( H a + H a 0 ) ( γ 0 C u 2 + γ 1 C H 2 + γ 1 P m C u C H ) .
Lemma 7.
Let, under assumptions (i)–(iii), (j) and (82), a pair ( x i , u i ) = ( u i , H i , p i , q i , q i , j i ) be a solution to problem (131) that corresponds to the function H d ( i ) L 2 ( Q ) 3 , i = 1 , 2 , where Q Ω is an arbitrary nonempty open subset. Assume that the condition (133) is satisfied. Then the following estimate for difference H H 1 H 2 holds:
H 1 H 2 Q Δ , Δ H d ( 1 ) H d ( 2 ) Q .
Proof of Lemma 7.
Setting H d = H d ( 1 ) H d ( 2 ) , in addition to (85), we have that ( I 2 ) u = 0 , ( I 2 ) p = 0 ,
( I 2 ) H ( x i ) , h = 2 ( H i H d ( i ) , h ) Q , ( I 2 ) H ( x 1 ) ( I ˜ 2 ) H ( x 2 ) , H = 2 ( H Q 2 ( H , H d ) Q ) .
In view of (135), identities (83), (84) and the main inequality (91) become
ν ( w , ξ i ) + ν 1 ( rot h , rot η i ) + ( ( u i · ) w , ξ i ) + ( ( w · ) u i , ξ i ) + κ [ ( rot η i × H i , w ) +
+ ( rot η i × h , u i ) ] κ [ ( rot H i × h , ξ i ) + ( rot h × H i , ξ i ) ] ( div w , σ i ) + ζ 1 ( i ) , w Σ +
+ ζ 2 ( i ) , h · n Σ τ + ζ 3 ( i ) , h × n Σ ν =
= μ 0 ( H i H d ( i ) , h ) Q ( w , h ) H T 1 ( Ω ) × H div s + 1 / 2 ( Ω ) ; ξ i V , i = 1 , 2 ,
( ( u · ) u , ξ 1 + ξ 2 ) κ ( rot H × H , ξ 1 + ξ 2 ) + κ ( rot ( η 1 + η 2 ) × H , u ) +
+ μ 0 ( H Q 2 ( H , H d ) Q ) μ 1 q s , Σ τ 2 μ 2 q s , Σ ν 2 μ 3 j Ω 2 .
We firstly estimate adjoint state variables ξ i and η i via ( H a + H a 0 ) . To this end, we set w = ξ i , h = η i in (136). In virtue of (16) and conditions ξ i | Σ = 0 , η i · n | Σ τ = 0 and η i × n | Σ ν = 0 , we obtain
ν ( ξ i , ξ i ) + ν 1 ( rot η i , rot η i ) + ( ( ξ i · ) u i , ξ i ) + κ ( rot η i × η i , u i )
κ ( rot H i × η i , ξ i ) = μ 0 ( H i H d ( i ) , η i ) Q .
Using (79), (80), (132), we derive, in addition to (121)–(123), that
H i H d ( i ) Q l H i H div s + 1 / 2 ( Ω ) + H d ( i ) Q l ( M H 0 + l 1 H d ( i ) Q )
δ 0 ν l ( H a + H a 0 ) / ( γ 1 κ ) .
Taking into account (99), (121)–(123) and (139), we infer from (138) that
( δ 0 ν / 2 ) ξ i 1 , Ω 2 + ( δ 1 ν m / 2 ) κ η i 1 , Ω 2 μ 0 ( H i Q ) + H d ( i ) Q ) η i Q
μ 0 l 2 ( M H 0 + l 1 H d ( i ) Q ) η i 1 , Ω μ 0 δ 0 ν γ κ 1 / 2 ( H a + H a 0 ) η i 1 , Ω
where γ = l 2 γ 1 1 . Using (79) and (139), we conclude from the last inequality that
ξ i 1 , Ω 2 μ 0 γ κ 1 P m ( H a + H a 0 ) , κ η i 1 , Ω 2 μ 0 γ κ 1 P m ( H a + H a 0 ) .
Combining (15), (104), (105), (115) and (140) gives the estimate
| ( ( u · ) u , ξ 1 + ξ 2 ) κ ( rot H × H , ξ 1 + ξ 2 ) + κ ( rot ( η 1 + η 2 ) × H , u ) |
12 μ 0 γ P m ( H a + H a 0 ) ( γ 0 C u 2 + γ 1 C H 2 + γ 1 P m C u C H ) ( q s , Σ τ 2 + q s , Σ ν 2 + j Ω 2 )
( 1 ε ) μ 1 q s , Σ τ 2 + ( 1 ε ) μ 2 q s , Σ ν 2 + ( 1 ε ) μ 3 j Ω 2 .
Using (141) from (137), we arrive at
μ 0 ( H Q 2 ( H , H d ) Q ) ( ( u · ) u , ξ 1 + ξ 2 ) +
+ κ ( rot H × H , ξ 1 + ξ 2 ) κ ( rot ( η 1 + η 2 ) × H , u ) μ 1 q s , Σ τ 2 μ 2 q s , Σ ν 2 μ 3 j Ω 2
ε μ 1 q s , Σ τ 2 ε μ 2 q s , Σ ν 2 ε μ 3 j Ω 2 .
Omitting the nonpositive term ε μ 1 q s , Σ τ 2 ε μ 2 q s , Σ ν 2 ε μ 3 j Ω 2 , we infer from (142) that H Q 2 H Q H d Q 0 or H Q H d Q . □
If Q = Ω , the estimate (134) is the stability estimate in the L 2 ( Ω ) 3 norm for the magnetic component H ^ of the solution ( u ^ , H ^ , p ^ , q ^ , q ^ , j ^ ) of problem (131) with respect to small disturbances of function H d L 2 ( Ω ) 3 in (131). If, besides, H d ( 1 ) = H d ( 2 ) in Q it follows from (134) that H 1 = H 2 in Q. This, together with (104), (105), (111) and (142), yields that q 1 = q 2 , q 1 = q 2 , j 1 = j 2 , u 1 = u 2 , H 1 = H 2 , p 1 = p 2 . The latter is equivalent to the uniqueness for the solution of (131).
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 s , Σ τ 2 + ε μ 2 q s , Σ ν + ε μ 3 j Ω 2 μ 0 ( H Q 2 + H Q H d Q )
μ 0 ( H Q 2 + H Q H d Q ) μ 0 H d Q 2 = μ 0 Δ 2
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.
Theorem 8.
Let parameters γ and H a 0 be defined by relations (132) and let conditions of Lemma 7 be fulfilled. Then the stability estimates (130) for control problem (131) hold where Δ = H d ( 1 ) H d ( 2 ) Q .
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.

6. 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.

Funding

This work was supported by the Russian Science Foundation Grant No. 22-21-00271.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Hou, L.S.; Meir, A.J. Boundary optimal control of MHD flows. Appl. Math. Optim. 1995, 32, 143–162. [Google Scholar] [CrossRef]
  2. Alekseev, G.V. Control problems for the steady–state equations of magnetohydrodynamics of a viscous incompressible fluid. J. Appl. Mech. Phys. 2003, 44, 890–898. [Google Scholar] [CrossRef]
  3. Alekseev, G.V. Solvability of control problems for stationary equations of magnetohydrodynamics of a viscous fluid. Siberian Math. J. 2004, 45, 197–213. [Google Scholar] [CrossRef]
  4. Alekseev, G.V.; Brizitskii, R.V. Control problems for stationary magnetohydrodynamic equations of a viscous heat-conducting fluid under mixed boundary conditions. Comput. Math. Math. Phys. 2005, 45, 2049–2065. [Google Scholar]
  5. Alekseev, G.V.; Brizitskii, R.V. Stability estimates of solutions of control problems for the stationary equations of magnetic hydrodynamics. Differ. Equ. 2012, 48, 397–409. [Google Scholar] [CrossRef]
  6. Mallea-Zepeda, E.; Ortega-Torres, E. Control problem for a magneto-micropolar flow with mixed boundary conditions for the velocity field. J. Dyn. Control Syst. 2019, 25, 599–618. [Google Scholar] [CrossRef]
  7. Griesse, R.; Kunish, K. Optimal control for stationary MHD system in velocity-current formulation. SIAM J. Control Opt. 2006, 45, 1822–1845. [Google Scholar] [CrossRef] [Green Version]
  8. Alekseev, G.V.; Brizitskii, R.V. Boundary control problems for the stationary magnetic hydrodynamic equations in the domain with non-ideal boundary. J. Dyn. Control Syst. 2020, 26, 641–661. [Google Scholar] [CrossRef]
  9. Gunzburger, M.; Trenchea, C. Analysis and discretization of an optimal control problem for the time–periodic MHD equations. J. Math. Anal. Appl. 2005, 308, 440–466. [Google Scholar] [CrossRef] [Green Version]
  10. Gunzburger, M.; Trenchea, C. Analysis of an optimal control problem for the three-dimensional coupled modified Navier-Stokes and Maxwell equations. J. Math. Anal. Appl. 2007, 333, 295–310. [Google Scholar] [CrossRef] [Green Version]
  11. Gunzburger, M.; Peterson, J.; Trenchea, C. The velocity tracking problem for MHD flows with distributed magnetic field controls. Int. J. Pure Appl. Math. 2008, 42, 289–296. [Google Scholar]
  12. Ravindran, S. On the dynamics of controlled magnetohydrodynamic systems. Nonlinear Anal. Model. Control 2008, 13, 351–377. [Google Scholar] [CrossRef]
  13. Ioffe, A.; Tikhomirov, V. Theory of Extremal Problems; North-Holland Publishing Co.: Amsterdam, The Netherlands, 1979. [Google Scholar]
  14. Fursikov, A.V. Optimal Control of Distributed Systems. Theory and Applications; AMS: Providence, RI, USA, 1999. [Google Scholar]
  15. Alekseev, G.V. Optimization in Stationary Problems of Heat and Mass Transfer and Magnetic Hydrodynamics; Nauchniy Mir: Moscow, Russia, 2010. (In Russian) [Google Scholar]
  16. Alekseev, G.V.; Brizitskii, R.V. On uniqueness and stability of solutions of extremum problems for stationary Navier–Stokes equation. Differ. Equ. 2010, 46, 70–82. [Google Scholar] [CrossRef]
  17. Solonnikov, V. On some stationary boundary value problems in magnetohydrodynamics. Trudy Mat. Inst. Steklov. 1960, 59, 174–187. (In Russian) [Google Scholar]
  18. Gunzburger, M.D.; Meir, A.J.; Peterson, J.S. On the existence, uniqueness, and finite element approximation of solution of the equation of stationary, incompressible magnetohydrodynamics. Math. Comput. 1991, 56, 523–563. [Google Scholar] [CrossRef]
  19. Schotzau, D. Mixed finite element methods for stationary incompressible magneto-hydrodynamics. Numer. Math. 2004, 96, 771–800. [Google Scholar] [CrossRef]
  20. Consiglieri, L.; Necasova, S.; Sokolowski, J. Incompressible Maxwell–Boussinesq approximation: Existence, uniqueness and shape sensitivity. Cont. Cybern. 2009, 38, 1193–1215. [Google Scholar] [CrossRef] [Green Version]
  21. Alekseev, G.V. Solvability of an inhomogeneous boundary value problem for the stationary magnetohydrodynamic equations for a viscous incompressible fluid. Differ. Equ. 2016, 52, 739–748. [Google Scholar] [CrossRef]
  22. Alekseev, G.V.; Brizitskii, R.V. Solvability analysis of a mixed boundary value problem for stationary magnetohydrodynamic equations of a viscous incompressible fluid. Symmetry 2021, 13, 2088. [Google Scholar] [CrossRef]
  23. Alekseev, G.V.; Brizitskii, R.V. Solvability of the boundary value problem for stationary magnetohydrodynamic equations under mixed boundary conditions for the magnetic field. Appl. Math. Lett. 2014, 32, 13–18. [Google Scholar] [CrossRef]
  24. Fernandes, P.; Gilardi, G. Magnetostatic and electrostatic problems in inhomogeneous anisotropic media with irregular boundary and mixed boundary conditions. Math. Model. Methods Appl. Sci. 1997, 7, 957–991. [Google Scholar] [CrossRef]
  25. Auchmuty, G. The main inequality of vector field theory. Math. Model. Methods Appl. Sci. 2004, 14, 79–103. [Google Scholar] [CrossRef]
  26. Auchmuty, G.; Alexander, J.S. Finite energy solutions of mixed 3D div-curl systems. Quart. Appl. Math. 2006, 64, 335–357. [Google Scholar] [CrossRef] [Green Version]
  27. Alekseev, G.V. Mixed boundary value problems for steady-state magnetohydrodynamic equations of viscous incompressible fluid. Comput. Math. Math. Phys. 2016, 56, 1426–1439. [Google Scholar] [CrossRef]
  28. Alekseev, G. Mixed boundary value problems for stationary magnetohydrodynamic equations of a viscous heat-conducting fluid. J. Math. Fluid Mech. 2016, 18, 591–607. [Google Scholar] [CrossRef]
  29. Alekseev, G.; Brizitskii, R. Solvability of the mixed boundary value problem for the stationary magnetohydrodynamic equations of viscous fluid. Daln. Math. J. 2002, 3, 285–301. (In Russian) [Google Scholar]
  30. Brizitskii, R.V.; Tereshko, D.A. On the solvability of boundary value problems for the stationary magnetohydrodynamic equations with inhomogeneous mixed boundary conditions. Differ. Equ. 2007, 43, 246–258. [Google Scholar] [CrossRef]
  31. Meir, A.J. The equation of stationary, incompressible magnetohydrodynamics with mixed boundary conditions. Comput. Math. Appl. 1993, 25, 13–29. [Google Scholar] [CrossRef] [Green Version]
  32. Villamizar-Roa, E.J.; Lamos-Diaz, H.; Arenas-Dias, G. Very weak solutions for the magnetohydrodynamic type equations. Discret. Contin. Dyn. Syst. B 2008, 10, 957–972. [Google Scholar] [CrossRef]
  33. Poirier, J.; Seloula, N. Regularity results for a model in magnetohydrodynamics with imposed pressure. C. R. Math. 2020, 58, 1033–1043. [Google Scholar] [CrossRef]
  34. Bermudez, A.; Munoz-Sola, R.; Vazquez, R. Analysis of two stationary magnetohydrodynamics systems of equations including Joule heating. J. Math. Anal. Appl. 2010, 368, 444–468. [Google Scholar] [CrossRef] [Green Version]
  35. Kim, T. Existence of a solution to the steady magnetohydrodynamics-Boussinesq system with mixed boundary conditions. Math. Meth. Appl. Sci. 2022, 45, 9152–9193. [Google Scholar] [CrossRef]
  36. Kim, T. Existence of a solution to the non-steady magnetohydrodynamics-Boussinesq system with mixed boundary conditions. J. Math. Anal. Appl. 2023, 525, 127183. [Google Scholar] [CrossRef]
  37. Alonso, A.; Valli, A. Some remarks on the characterization of the space of tangential traces of H(rot; Ω) and the construction of the extension operator. Manuscr. Math. 1996, 89, 159–178. [Google Scholar] [CrossRef]
  38. Girault, V.; Raviart, P. Finite Element Methods for Navier-Stokes Equations: Theory and Algorithms; Springer: Berlin/Heidelberg, Germany, 1986. [Google Scholar]
  39. Tikhonov, A.N.; Goncharsky, A.V.; Stepanov, V.V.; Yagola, A.G. Numerical Methods for the Solution of Ill-Posed Problems; Springer: Amsterdam, The Netherlands, 2013. [Google Scholar]
  40. Klibanov, M.V.; Kolesov, A.E. Convexification of a 3-D coefficient inverse scattering problem. Comp. Math. Appl. 2019, 77, 1681–1702. [Google Scholar] [CrossRef] [Green Version]
  41. Klibanov, M.V.; Li, J. Inverse Problems and Carleman Estimates; Walter de Gruyter GmbH: Berlin, Germany, 2021; Volume 63. [Google Scholar]
  42. Beilina, L.; Hosseinzadegan, S. An adaptive finite element method in reconstruction of coefficients in Maxwell’s equations from limited observations. Appl. Math. 2016, 61, 253–286. [Google Scholar] [CrossRef] [Green Version]
  43. Beilina, L.; Smolkin, E. Computational design of acoustic materials using an adaptive optimization algorithm. Appl. Math. Inf. Sci. 2018, 12, 33–43. [Google Scholar] [CrossRef] [Green Version]
  44. Cakoni, F.; Kovtunenko, V.A. Topological optimality condition for the identification of the center of an inhomogeneity. Inverse Probl. 2018, 34, 035009. [Google Scholar] [CrossRef]
  45. Kovtunenko, V.A.; Kunisch, K. High precision identification of an object: Optimality conditions based concept of imaging. SIAM J. Control Optim. 2014, 52, 773–796. [Google Scholar] [CrossRef]
  46. Fershalov, Y.Y.; Fershalov, M.Y.; Fershalov, A.Y. Energy efficiency of nozzles for axial microturbines. Proc. Eng. 2017, 206, 499–504. [Google Scholar] [CrossRef]
  47. Fershalov, Y.Y.; Fershalov, A.Y.; Fershalov, M.Y. Microturbinne with new design of nozzles. Energy 2018, 157, 615–624. [Google Scholar] [CrossRef]
  48. Mallea-Zepeda, E.; Ortega-Torres, E.; Villamizar-Roa, E.J. A boundary control problem for micropolar fluids. J. Optim. Theory Appl. 2016, 169, 349–369. [Google Scholar] [CrossRef]
  49. Almeida, A.; Chemetov, N.V.; Cipriano, F. Uniqueness for optimal control problems of two-dimensional second grade fluids. Electron. J. Differ. Equ. 2022, 22, 1–12. [Google Scholar]
  50. Baranovskii, E.S. Optimal boundary control of the Boussinesq approximation for polymeric fluids. J. Optim. Theory Appl. 2021, 189, 623–645. [Google Scholar] [CrossRef]
  51. Shestopalov, Y.V.; Smirnov, Y.G. Determination of permittivity of an inhomogeneous dielectric body in a waveguide. Inverse Probl. 2011, 27, 095010. [Google Scholar] [CrossRef]
  52. Cea, J. Optimization: Theory and Algorithms; Springer: New York, NY, USA, 1978. [Google Scholar]
  53. Shercliff, J. A Textbook of Magnetohydrodynamics; Pergamon Press: Oxford, UK, 1965. [Google Scholar]
Figure 1. Example of simply connected domain with n = 0 , m = 2 .
Figure 1. Example of simply connected domain with n = 0 , m = 2 .
Mathematics 11 02610 g001
Figure 2. Example of multi-connected domain with n = 1 , m = 2 .
Figure 2. Example of multi-connected domain with n = 1 , m = 2 .
Mathematics 11 02610 g002
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Alekseev, G. Analysis of Control Problems for Stationary Magnetohydrodynamics Equations under the Mixed Boundary Conditions for a Magnetic Field. Mathematics 2023, 11, 2610. https://doi.org/10.3390/math11122610

AMA Style

Alekseev G. Analysis of Control Problems for Stationary Magnetohydrodynamics Equations under the Mixed Boundary Conditions for a Magnetic Field. Mathematics. 2023; 11(12):2610. https://doi.org/10.3390/math11122610

Chicago/Turabian Style

Alekseev, Gennadii. 2023. "Analysis of Control Problems for Stationary Magnetohydrodynamics Equations under the Mixed Boundary Conditions for a Magnetic Field" Mathematics 11, no. 12: 2610. https://doi.org/10.3390/math11122610

APA Style

Alekseev, G. (2023). Analysis of Control Problems for Stationary Magnetohydrodynamics Equations under the Mixed Boundary Conditions for a Magnetic Field. Mathematics, 11(12), 2610. https://doi.org/10.3390/math11122610

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop