Inhomogeneous Boundary Value Problems for the Generalized Boussinesq Model of Mass Transfer

: We consider boundary value problems for a nonlinear mass transfer model, which generalizes the classical Boussinesq approximation, under inhomogeneous Dirichlet boundary conditions for the velocity and the substance’s concentration. It is assumed that the leading coefficients of viscosity and diffusion and the buoyancy force in the model equations depend on concentration. We develop a mathematical apparatus for studying the inhomogeneous boundary value problems under consideration. It is based on using a weak solution of the boundary value problem and on the construction of liftings of the inhomogeneous boundary data. They remove the inhomogeneity of the data and reduce initial problems to equivalent homogeneous boundary value problems. Based on this apparatus we will prove the theorem of the global existence of a weak solution to the boundary value problem under study and establish important properties of the solution. In particular, we will prove the validity of the maximum principle for the substance’s concentration. We will also establish sufficient conditions for the problem data, ensuring the local uniqueness of weak solutions.


Introduction and Statement of the Boundary Value Problem
A large number of works are devoted to the issues concerning mathematical modeling of the processes of fluid flows taking into account convective effects.These papers explore both applied questions, concerning the derivation and justification of basic convection models, and purely mathematical questions.These include methods for constructing exact and approximate solutions, the theoretic-group analysis (or Li-Ovsyannikov symmetry) method of investigating the qualitative properties of solutions to differential equations underlying the convection models under study, and theoretical analysis of the solvability and uniqueness of solutions of boundary value problems for basic convection models.
Among the various convection models for binary and/or thermally conducting liquids, an important role is played by so-called Oberbeck-Boussinesq models [1,2], which are derived from the exact Navier-Stokes equations of fluid dynamics, taking into account the following assumptions (see [1] Section 54): 1. Convective motion is similar to that of an incompressible fluid with constant density ρ 0 , but the possible deviation of the true density ρ from ρ 0 is taken into account in the momentum conservation equation in the form of a term describing the additional volumetric force -the buoyancy (Archimedes) force.
2. The change in density is caused by changes in the temperature and concentration of the dissolved substance but not changes in the pressure.
3. The velocity gradients are small enough so that the process of transition of work to heat during movement does not lead to a change in the temperature of the medium.
The model that results from these simplifications is called the Oberbeck-Boussinesq model.In turn, the Oberbeck-Boussinesq model allows for further simplification, often called the simplified or classical Boussinesq model.
When modeling heat transfer (HT) processes in a viscous heat-conducting liquid, by the (classical) Boussinesq model one usually understands the model in which the bulk buoyancy force included in the momentum conservation equation linearly depends on temperature, and, in addition, the main parameters of the fluid; namely, the viscosity coefficient and the thermal conductivity coefficient are positive constants.When modeling mass transfer (MT) processes in a binary fluid, by the (classical) Boussinesq model one usually understands the model in which the bulk buoyancy force (included in the momentum conservation equation) linearly depends on the concentration of the dissolved substance while the main parameters, namely the viscosity coefficient, diffusion coefficient, as well as another parameter, called the reaction coefficient, are the positive constants.The latter parameter is responsible for the possible decay of the dissolved substance in the main medium due to the chemical reaction.When all or some of the above conditions are not fulfilled, the corresponding model is often referred to as a generalized Boussinesq model.
Let us emphasize that for the classical Boussinesq model, many theoretical issues are quite fully studied.This, in particular, holds for the study of the correctness of boundary or initial-value problems for stationary or non-stationary models of heat and mass transfer (HMT).Among many works in this area, we note the works [3,4] by the first author and his coauthors on the study of the correctness of boundary value problems for stationary equations of HT and/or MT.We emphasize that in these works, in addition to the study of correctness of a boundary value problems, a theoretical analysis of control problems for the models of HMT was performed.The analysis of the results obtained in [3,4] made it possible to identify interesting regularities related to the interaction of hydrodynamic and thermal fields in binary and/or heat-conducting media and, in particular, to establish the most effective mechanisms for controlling thermohydrodynamic processes in viscous liquids.The close problems of boundary or distributed control for the HT equations in the Boussinesq approximation have also been investigated in the works [5][6][7][8][9][10].
Theoretical questions for the generalized Boussinesq model have been studied to a much lesser extent.However, significant progress has been made in recent years in this area as well.Over the past decades, a large number of papers have been published regarding the study of HMT equations with variable transfer coefficients and with variable buoyancy force depending on the temperature and/or concentration of the dissolved substance.These works can be divided into several groups.The first group contains papers that develop methods for finding exact solutions to these equations (see, for example, [11][12][13][14][15] and reviews [16]).The second group contains works devoted to application of the Li-Ovsyannikov symmetry method to study qualitative properties of solutions of equations of HMT in viscous binary and/or heat-conducting liquids.This group includes a very large quantity of works (see, e.g., [17][18][19][20][21], monographs [22] and reviews [23,24]).Another group of works is that in which mathematical modeling of fluid motion processes takes into account thermodiffusion effects (or Sorét effects) and/or concentration diffusion effects (or Dufort effects).A detailed list and analysis of these works can be found in the reviews [23,24].
At the same time, the authors know of only a few papers in which the solvability of boundary value problems for equations of HMT with variable coefficients is investigated.The mentioned works can be divided into several groups.The first group includes works [25][26][27][28], devoted to the study of the solvability of boundary value problems for stationary equations of HMT.The second group includes works [29][30][31][32][33][34], devoted to the study of the solvability of boundary value and control problems for non-stationary Boussinesq equations of HT (or MT).The works [35][36][37] form one more group in which the solvability of boundary value and control problems is studied for the stationary MT model in the case where the reaction coefficient can depend on the concentration of matter and spatial variables.
Close questions on the study of the correctness of boundary value or control problems for stationary equations of magnetic hydrodynamics of viscous incompressible or heatconducting liquid in the Boussinesq approximation were investigated in [38][39][40][41][42][43].In [44], the solvability of the initial-boundary problem for the non-stationary MHD-Boussinesq system is considered under mixed boundary conditions for velocity, magnetic field, and temperature, in the case when the viscosity coefficient, magnetic permeability, electrical conductivity, thermal conductivity, and specific heat of the fluid depend on the temperature.
The purpose of this work is to analyze the global solvability and local uniqueness of solutions of the boundary value problem for a generalized Boussinesq MT model describing the flow of binary fluid in which the diffusion, viscosity, and reaction coefficients and the buoyancy force depend on the substance concentration.
The paper is organized as follows.In Section 2, we will formulate the main boundary value problem, to which we will refer below as Problem 1.In addition, we introduce functional spaces and formulate a number of auxiliary results in the form of Lemmas 1-3, which will be used when studying the solvability and uniqueness of Problem 1.In Section 3, we will formulate and prove the theorem on the global existence of a weak solution to Problem 1 and establish the maximum principle for substance concentration φ.In Section 4, we will establish sufficient conditions on the data of Problem 1 that provide the conditional uniqueness of the weak solution having an additional property of smoothness for concentration.Section 5 contains a discussion on the relationship of results of this paper with mathematical modeling method.The last Section 6 contains a brief summary of the results obtained in our paper.

Statement of the Main Problem, Functional Spaces
Let Ω be a bounded domain in the space R 3 with a Lipschitz boundary Γ.Below, we will consider the following boundary value problem describing the motion of binary fluid within the framework of the generalized Boussinesq model of MT: −div(λ(φ)∇φ) Here, u is the velocity vector, φ is the concentration of dissolved substance, p = P/ρ 0 , where is the gravitational acceleration, and f or f is the bulk density of external forces or of the external sources of substance, respectively.One can read about the derivation of the model ( 1), (2) in the general case of convective approach and of boundary conditions (3) in ( [1] Sections 54 and 55).Below, the problem (1)-(3) for the given functions ν(φ), λ(φ), b(φ), k(φ), f, and f will be referred to as Problem 1.
Unlike [25,26,28], we consider inhomogeneous boundary conditions for velocity and for concentration.This allows us to apply the theory, developed in our work, for a wide class of physically interesting problems of flowing binary fluid through a limited region.Typical examples of this kind of fluid flows are presented in Figure 1, showing possible fluid flow situations in flat channels (Figure 1a,b) and in a cylindrical channel (Figure 1c).Clearly, physical statements of fluid flow problems in such channels require the use of inhomogeneous boundary conditions for both velocity and concentration at the inflow Γ in and outflow Γ out sites.An important role in our analysis will be played by the following functional spaces: It is well known (see, for example, [54,55]) that for the domain Ω with the Lipschitz boundary, the spaces H and V are characterized as follows: and denote by X * the space (H −1 (Ω) 3 × H −1 (Ω), which is the dual of X.
Recall that our goal is to prove the global solvability of the boundary value problem (1)-(3) and the local uniqueness of its solution.To achieve this goal, we will need to introduce a number of assumptions on the problem's data.More specifically, we will assume in what follows that the following conditions are met.
2.1.Ω is a bounded domain in R 3 with a boundary Γ ∈ C 0,1 consisting of N component

2.4.
For any function φ ∈ H 1 (Ω), the embedding b(φ) ∈ L p (Ω) is valid where p ≥ 3/2 is a fixed number independent of φ, and the vector function b Here βp is a positive constant dependent on p.In addition, for any pair of functions φ 1 , φ 2 ∈ H 1 (Ω) belonging to sphere B r = {φ ∈ H 1 (Ω) : ∥φ∥ 1,Ω ≤ r} of radius r, the following inequality is true: Here, L b is a constant that depends on b and on r but does not depend on φ 1 , φ 2 ∈ B r .

2.5.
For any function φ ∈ H 1 (Ω), the embedding k(φ) ∈ L p + (Ω) is valid, where p ≥ 3/2 is a fixed number independent of φ, and the following estimate for ∥k(φ) Here, γp is a positive constant dependent on p. Also, for any sphere B r = {φ ∈ H 1 (Ω) : ∥φ∥ 1,Ω ≤ r}, the following inequality holds.Here, L k is a constant that depends on k and on r but does not depend on φ 1 , , and there are positive constants ν min , ν max , λ min , and Consider the function µ ∈ C 0 (R), satisfying the following condition: Clearly, µ(h) ∈ L ∞ (Ω) for any h ∈ H 1 (Ω), and the following estimates in Ω, the Lebesgue's theorem on majorant convergence implies that The property (11) will play an important role when proving the solvability of Problem 1, which contains the variable leading coefficients ν(φ) and λ(φ).
Below, we will often use the following inequalities: Here, C s is a constant dependent on Ω and s ∈ [1,6], and C p is a constant dependent on Ω and p at p ≥ 3/2.The inequality ( 12) is a consequence of Sobolev's embedding theorem, according to which the space H 1 (Ω) is imbedded into L s (Ω) continuously at s ≤ 6 and compactly at s < 6.The inequality ( 14) is a consequence of the Hölder inequality for the three functions.In turn, the following inequality is the consequence of inequalities ( 12) and ( 14), where A similar inequality holds for the scalar functions q, φ, and h.It has the form Along with inequalities ( 12)-( 16), we will use a number of other important inequalities and properties of bilinear and trilinear forms, which we will write as the following Lemma.Lemma 1.Let conditions 2.1, 2.4, 2.5, and 2.6 be met, and let u ∈ H 1 div (Ω) be a given function.Then, there are the positive constants δ 0 , δ 1 , γ 1 , γ ′ 1 , γ 2 , γ ′ 2 , γ 3 , and β, which depend on Ω, and the constants β p and γ p depending on Ω and p, such that the following relationships are fulfilled: sup Here , where the constants C P , βp , γp and Ĉ6 were defined in (4), ( 5), (7), and ( 15) respectively; ν min , ν max , λ min and λ max -constants defined in (9).We prove, for example, (27) and establish a relationship between the constant β p and the constant βp defined in (5).To do this, we use the inequality (15) at q = b(φ) and property (5).Using these relations successively, we have : Proof of the remaining statements constituting Lemma 1, we leave to the reader.They can also be found in [54][55][56][57].
Similarly, the following estimate for the difference k(φ 1 ) − k(φ 2 ) follows from the estimate (25) and the property (8) of the function k(φ): The following lemmas about the existence of liftings for the velocity and concentration will play an important role below.Lemma 2. Let under the assumption 2.1 the boundary vector g satisfy the condition 2.3..Then, for any arbitrary number ε > 0, there exists a function (velocity lifting) Here, C ε is a constant dependent on ε and Ω.
Lemma 3. Let assumption 2.1 be fulfilled.Then, there is a family of continuous non-decreasing functions M ϵ : R + ≡ (0, ∞) → R with M ϵ (0) = 0 depending on the parameter ϵ ∈ (0.1] as well as from Ω, such that for any function ψ ∈ H 1/2 (Γ) that is not equal to zero, there exists a function φ ϵ ∈ H 1 (Ω) that satisfies the conditions Proof of Lemma 3. Detailed proof of Lemmas 2 and 3 can be found in book [56] (Appendix 2).

Global Solvability of Problem 1
Our nearest goal is to derive the weak formulation of Problem 1 and to prove the existence of its solution.To this end, we multiply the first equation in (1) by the function v ∈ H 1 0 (Ω) 3 , Equation (2) by the function h ∈ H 1 0 (Ω), and integrate the result over Ω using Green's formulas.As a result, we will get a weak formulation of Problem 1.It consists in finding a trio of functions (u, φ, p) ∈ H 1 (Ω) 3 × H 1 (Ω) × L 2 0 (Ω) that satisfy the relations: The specified three functions (u, φ, p) satisfying ( 31)- (33), will be called below a weak solution to Problem 1.Consider the restriction of the identity (31) by the space V, which, taking into account the condition divv = 0 for v ∈ V, becomes : It is well known that for the proof of existence of a weak solution of Problem 1, it is enough to prove the existence of the solution (u, φ) ∈ H 1 (Ω) 3 × H 1 (Ω) of problem ( 32)- (34), and then, using the standard scheme, to restore the pressure p ∈ L 2 0 (Ω) so that the identity (31) is fulfilled.One can read more about pressure recovery in ( [55,56], p. 134, Section.3.2).
We will look for the solution (u, Here, u 0 ≡ u ε 0 and φ 0 = φ ϵ 0 are the velocity and concentration liftings defined above, and ũ ∈ V and φ ∈ H 1 0 (Ω) are the new unknown functions that we are looking for.The value ε 0 of the parameter ε will be chosen so that the following conditions are fulfilled: The corresponding value ϵ 0 of the parameter ϵ will be selected below.Substituting (35) in ( 32), (34), we come to the following relations with respect to the pair ( ũ, φ): 37), ( 38), we apply Schauder's fixed point theorem according to the scheme proposed in ( [56], Section 4.2).To this end, we define the pairs z = (w, η) ∈ W, y = ( ũ, φ) ∈ W and construe the operator F : W → W acting by: F(z) = y, where y = ( ũ, φ) ∈ W is the solution of the linear problem Here, (w, η) ∈ V × H 1 0 (Ω) is a given pair of functions, and functionals f 1 : V → R and f 1 : Using ( 17)-( 19), ( 22)-( 25), ( 27), ( 28), ( 35), (36), and Lemmas 2 and 3, we deduce that From ( 43), ( 46), ( 48), ( 49), (51), and ( 52), it follows that f 1 ∈ V * , f 1 ∈ H −1 (Ω), and the following estimates are valid: ∥ In addition, for each fixed pair (w, η) : 40) is continuous and coercive with the constant λ * defined in (23), while the right-hand side of the identity (40) defines a linear continuous functional over H 1 0 (Ω), and the following estimates hold: |⟨ In this case, it follows from the Lax-Milgram theorem that for any pair (w, η) 40) exists and is unique.In addition, by ( 55) and ( 56) and Lemma 3, the following estimates are performed for φ and φ ≡ φ 0 + φ: Let us turn to the problem (39) where we put φ = φw,η .From ( 51), (53), and (58), it follows that the right-hand side in (39) is a value on the element v ∈ V of the linear continuous functional of V * and we have In turn, from the estimates (43), (44), and (45) and from the second identity in (20), it follows that the bilinear form a w,η 1 : V × V → R defined in (39) is continuous and coercive with constant (ν * /2), where ν * is defined in (18).From Lax-Milgram's theorem, it then follows that for any pair (w, η) ∈ V × H 1 0 (Ω), there is a single solution ũ ∈ V of the problem (39), and the following estimate is performed: Let us assume that ν * ≤ 2β p c * γ ′ 2 and choose the parameter ϵ from the condition From ( 59), we infer that Considering (60), from (57), we conclude that An important feature of the estimates (61) and ( 62) is the fact that their right-hand parts depend only on ∥w∥ 1,Ω but do not depend on ∥η∥ 1,Ω .This will allow us to select below the convex closed set K ⊂ W, which the operator F maps into itself.Still, it was assumed that (w, η) is an arbitrary pair of functions from W. Suppose now that η is still arbitrary, while w satisfies the condition By this assumption, it follows from (61) that According to (63), we have that In this case, from (64), we conclude that Similarly, considering (63), from (62), we get Thus, we showed that the operator F : W → W maps a pair (w, η) ∈ W, where w satisfies (63) and η is an arbitrary function, into the pair ( ũ, φ) ∈ W, for which the next estimates are fulfilled: Therefore, if we choose a convex closed set in the space W, then the obtained estimates (67) mean that the operator F maps the set K to itself.
To apply Schauder's theorem, we have to prove that the operator F is continuous and compact on the set K defined in (68).To this end, denote by z n = (w n , η n ), n = 1, 2, . . .an arbitrary sequence from K. Let us put y n ≡ ( ũn , φn ) = F(z n ), n = 1.2, . .., and show that from the sequence y n , we can extract a subsequence converging in the norm defined in X to some element y ∈ K.
From ( 78) and (79), it follows that the operator F : W → W is continuous and compact on the set K. In this case, Schauder's fixed point theorem implies that the operator F has a fixed point y ≡ ( ũ, φ) = F(z) ∈ K, which by construction is the desired solution to the problem (37) and (38) and satisfies estimates in (67).From this important fact, in turn, it follows that the pair (u, φ) ∈ H 1 (Ω) 3 × H 1 (Ω) defined in (35) is the desired solution to the problem ( 32)-( 34) and the following estimates are valid: We formulate the results obtained in the form of the following theorem.
Remark 2. The theorem corresponding to the case when the condition 2 can be proved in a similar way.
The existence of pressure p ∈ L 2 0 (Ω), which together with the specified pair (u, φ) satisfies the relation (31), is proved as in ( [56], Section 3.2).It remains to derive the estimate for p.For this purpose, we will use the inf-sup condition (21), according to which for the above function p and any (arbitrarily small) number δ > 0, there exists a function Assuming v = v 0 in (31) and using the last inequality and estimates ( 17), ( 19), ( 27), we deduce that Dividing by ∥v 0 ∥ 1,Ω ̸ = 0 and taking into account the estimate (84), we derive from (86) that We formulate the result in the form of the following theorem.
In conclusion of this Section, we will establish sufficient conditions on the data of Problem 1, under which the maximum (or minimum) principle is valid for the concentration component φ of the solution (u, φ, p) of Problem 1.
Let ψ min , ψ max be nonnegative numbers, f max be a positive number and, in addition to 2.1-2.6, the following conditions hold: in Ω; ψ min ≤ ψ ≤ ψ max a.e. on Γ.

3.2.
The nonlinearity k(φ)φ is monotonic in the following sense:

3.3.
It is assumed that every of functional for M 1 or m 1 equations has at least one (positive) solution M 1 or m 1 .
We set Theorem 3. Let conditions 2.1-2.6 and 3.1-3.3be satisfied.Then, for the component φ of the weak solution (u, φ, p) ∈ H 1 (Ω) 3 × H 1 (Ω) × L 2 0 (Ω) of Problem 1, the maximum and minimum principle hold having the form: m ≤ φ ≤ M everywhere in Ω. (90) Here, m and M are constants defined in relation ( 89) in which M 1 is chosen as a minimum root of the first equation in (88) while m 1 is chosen as a maximum root of the second equation in (88).
Proof of Theorem 3. Let M be constant as defined in (89).First, we will prove that φ ≤ M a.e. in Ω.To this end, we will define the function φ = max{φ − M, 0}.It is clear that the principle of maximum or estimate φ ≤ M holds a.e. in Ω if and only if φ = 0 everywhere in Ω.
Denote by Ω M ⊂ Ω an open measurable subset of Ω in which φ > M. From ([57], p. 152), it follows that ∇ φ = ∇φ everywhere in Ω M and, in addition, φ ∈ H 1 0 (Ω).So, the following equalities take place: Taking into account (91), let us set h = φ in (32).We obtain that The following equations follow from the properties of function φ: Due to the property 3.2 for the functions φ 1 = φ + M and φ 2 = M, which belong to H 1 (Ω), and by condition φ = 0 in Ω \ Ω M , the following relation takes place: Now, we subtract the equality (k(M)M, φ) ≡ (k(M)M, φ) Ω M from both parts of (92).As a result, we arrive at Using ( 23) and (93), we deduce from (94) that By definition of M in (89), it follows from (95) that (if M is selected from the first condition in (88)) then φ = 0 and therefore φ ≤ M in Ω.
The principle of minimum is proved in a similar way using the non-positive function φ = min{φ − m, 0} (see in more detail in [58]).

Conditional Uniqueness of Solution of Problem
In this Section, we will prove the conditional uniqueness of the weak solution (u, φ, p) of Problem 1 provided that the component φ has an additional property of smoothness, namely: φ ∈ H 2 (Ω).We will assume, in addition to the property 2.6, that the leading coefficients ν(•) and λ(•) have the following properties: 4.1.λ(•) belongs to the space C 1 (R), and where λ ′ min and λ ′ max are positive constants.4.2.Functions ν, λ and λ ′ are Lipschitz continuous, i.e., |ν(s Here, L ν , L λ and L ′ λ are some positive constants.To achieve our goal, we will use as in [25,26,28] the equivalence between the standard ∥ • ∥ 2,Ω and L 2 -norm ∥∆ • ∥ Ω of the Laplace operator in the space H 2 (Ω) ∩ H 1 0 (Ω) for the domain with the boundary Γ ∈ C 2 (see [54,55]).The mentioned equivalence is described by the following inequalities: Here and below, Ci , i = 1, 2, . . .are positive constants that depend on Ω and may be also some indexes .
Recall that in [28] the proof of the local uniqueness of the solution to the boundary value problem under study was based precisely on the equivalence property of the norms ∥ • ∥ 2,Ω and ∥∆ • ∥ Ω , which is valid for the components , where the equations of the Boussinesq model are considered under homogeneous Dirichlet conditions, in our case, the specified equivalence property for components φ 1 and φ 2 does not work because they satisfy the inhomogeneous boundary condition φ i | Γ = ψ, i =1.2.However, this property remains valid for the difference φ = φ 1 − φ 2 , which belongs to H 1 0 (Ω).And, as will be shown below, this property is quite enough to prove the local uniqueness of the weak solution to Problem 1, which has an additional smoothness of the form φ ∈ H 2 (Ω).
Below, we will also use the following estimates: which result from the continuity of the embedding The latter is a consequence of estimates (98) and (101).In (99), C3 is a positive constant that depends on Ω and r ∈ [1, ∞).
The following result is valid regarding the uniqueness of a weak solution to Problem 1, which has a certain property of smallness: Theorem 4. Let, in addition to conditions 2.1-2.6 and 4.1, 4.2, the following conditions be met: Γ ∈ C 2 , f ∈ L 2 (Ω), while the conditions 2.1 and 2.4 are met for some p ≥ 2. If there exists such a number ε > 0 that there is a weak solution (u, φ, p) then this solution is unique.
It remains to prove that the components p 1 and p 2 introduced at the beginning of the proof of Theorem 4 coincide.To this end, subtract the identity (31) for (u 2 , φ 2 , p 2 ) from ( 31) for (u 1 , φ 1 , p 1 ) and take into account that u 1 = u 2 and φ 1 = φ 2 .As a result, we obtain that the difference p = p 1 − p 2 satisfies the identity From (133), it follows due to the inf-sup condition (21) that p = 0 or p 1 = p 2 .Thus, the local uniqueness property of the weak solution of Problem 1 from the space H 1 (Ω) 3 × H 2 (Ω) × L 2 0 (Ω) is proven.

Discussion
It is well known that the method of mathematical modeling (MMM) plays a fundamental role in the study of various kinds of processes, phenomena, and systems.This method allows one to investigate processes, phenomena, objects, and systems using mathematical methods, analytical or numerical, oriented to the use of computers and computational experiment (see in more detail in [59,60]).
The use of MMM for the study of a specific physical (biological, chemical, or other) process can be divided into several stages.The first stage consists in the fact that on the basis of natural laws describing the behavior of the process, a system of mathematical relations is derived with respect to the values characterizing the process in question.This system, called a mathematical model, usually consists of differential equations to which boundary (and initial in the case of a non-stationary process) conditions can be applied to distinguish a unique solution with respect to the resulting boundary value problem of mathematical physics.
The next stage, known as "theoretical analysis of boundary value problems", is to investigate the important properties of the solution to the boundary value problem obtained at the first stage without knowing (i.e., without finding) the solution itself, but only using information about the mathematical model under consideration and about the boundary and initial conditions.During such a stage, the correctness of the boundary value problem is investigated, i.e., conditions are established for the initial data of the considered problem, under which the solution exists, that are unique and stable with respect to small disturbances in the problem's data.Finally, in the third stage, the solution to the problem under consideration is supposed to be found: preferably, an exact solution, but, as a rule, an approximate one by using approximate analytical or numerical methods focused on the use of computers.
We emphasize that during the third stage, a specific boundary value problem of mathematical physics is usually solved, corresponding to a specific domain and specific coefficients of the model under consideration, while in the second stage, a whole class of problems under study can be investigated, including different domains and different coefficients of the model.
Let us explain the aforementioned on the example to be used in our manuscript, which examines a whole class of boundary value problems for a stationary model of MT within the generalized Boussinesq approximation under inhomogeneous boundary conditions for the velocity and concentration.The role of the problem's data is played by the flow domain Ω, which might have an arbitrary shape and be piecewise smooth, by the basic coefficients of viscosity, diffusion, reaction, and buoyancy, depending on the concentration, as well as by the functions included in the boundary conditions.Just for this large class of boundary value problems, we have established sufficient conditions on the data (they have the form of conditions 2.1-2.6),under which the solution of the boundary value problem (1)-( 3) exists (Theorems 1 and 2), and it has the property of local uniqueness (Theorem 4) and the maximum principle is valid for it (Theorem 3).We emphasize that the results obtained are valid for any data satisfying conditions 2.1-2.6.
But, if one chooses a specific problem from our class, i.e., one considers a specific flow domain Ω, specific coefficients of viscosity, diffusion, reaction, buoyancy force, and boundary functions, then one can, using a suitable numerical method, find the approximate solution we are looking for and plot the necessary graphs to visualize the properties of this solution.However, the latter is no longer the goal of our work, which is devoted precisely to the theoretical analysis of the boundary value problem (1)- (3).But this can be the goal of another work devoted to the numerical study of the particular boundary value problem under consideration.

Conclusions
In this paper, the inhomogeneous boundary value problem for the stationary model of MT within the generalized Boussinesq approximation was studied.It is assumed that the coefficients ν, λ, k, and b of the model depend on the concentration of the substance.We defined a concept of a weak solution of the boundary value problem under study and established the conditions on the coefficients ν, λ, k, b and on other data that provide the global solvability of the problem.We also found additional important properties of the weak solution and in particular the maximum principle for the substance's concentration.Also, we have proved conditional uniqueness of the weak solution having an additional smoothness property for the concentration.
Teh main results of our work were obtained by using a number of the author's findings.The first finding is related to the possibility of spreading the mathematical theory developed in [28] in the study of homogeneous boundary value problems for the MT model to the case of the inhomogeneous boundary value problem considered in the manuscript.This possibility was realized by using special extensions of inhomogeneous boundary data entered in (3) inside the domain Ω.The existence of these extensions, called liftings of inhomogeneous boundary data, is ensured by Lemmas 2 and 3 given in Section 2.
The second finding is related to the proof of the maximum principle for concentration, which makes it possible to estimate the values of concentration within the flow region by the values for the concentration at the boundary of the domain.The use of the maximum principle plays an important role in the development and analysis of effective numerical algorithms for solving the boundary value problem under consideration.Finally, the third finding is related to the method proposed in the work to prove the local uniqueness of a weak solution.It is based on using thin estimates of bilinear and trilinear forms included in the definition of the weak solution.These estimates are formulated in Lemma 1.

Figure 1 .
Figure 1.Typical examples of domains Ω having the form of the channel when it is required to set inhomogeneous boundary conditions for both velocity and concentration at the inflow Γ in and outflow Γ out .When studying Problem 1, we will use the Sobolev functional spaces H s (D), s ∈ R. Here, D means either domain Ω or some subset Q ⊂ Ω, or the boundary Γ.By ∥ • ∥ s,Q , | • | s,Q , and (•, •) s,Q the norm, half-norm, and inner product in H s (Q) will be defended, respectively.The norms and scalar product in L 2 (Q) or in L 2 (Ω) are denoted by ∥ • ∥ Q , (•, •) Q or by ∥ • ∥ Ω and (•, •), respectively.X * denotes the dual space of Hilbert space X, while the duality relation for the pair of dual spaces X and X * is written as ⟨•, •⟩ X * ×X , or simply as ⟨•, •⟩.