Preprint
Article

On the Stabilization of the Solution to the Initial-Boundary Value Problem for One-Dimensional Isothermal Equations of Viscous Compressible Multicomponent Media Dynamics

This version is not peer-reviewed.

Submitted:

16 June 2023

Posted:

19 June 2023

You are already at the latest version

A peer-reviewed article of this preprint also exists.

Abstract
An initial-boundary value problem is considered for one-dimensional isothermal equations of the dynamics of viscous compressible multicomponent media, which are a generalization of the Navier--Stokes equations. We prove the stabilization of the solution to the initial-boundary value problem while the time tends to infinity, without simplifying assumptions for the structure of the viscosity matrix, except for the standard physical requirements of symmetry and positive definitenes.
Keywords: 
viscous compressible medium; multicomponent flows; stabilization of solution
Subject: 
Computer Science and Mathematics  -   Applied Mathematics

1. Statement of the initial-boundary value problem and formulation of the main result

Modeling the motion of multicomponent media, and hence the solution of the resulting mathematical problems, is an interesting and relatively little-studied problem of both physics and mathematics. Until now, a unified approach to this problem has not been developed, and there is no advanced mathematical theory about the existence, uniqueness, and properties of solutions to initial-boundary value problems that arise in this modeling. The purpose of this article does not include a detailed review of this problem as a whole; to a certain extent, an idea of it can be obtained from monographs [1,2], as well as from the reviews given in the articles [3,4]. The aim of the article is a mathematical study of the asymptotic behavior of the solution at t + of the following initial-boundary value problem for one-dimensional isothermal equations of viscous compressible multicomponent media:
ρ t + ( ρ v ) x = 0 ,
ρ u i t + ρ v u i x + K ρ x = j = 1 N ν i j 2 u j x 2 , i = 1 , , N ,
ρ | t = 0 = ρ 0 ( x ) , u i | t = 0 = u 0 i ( x ) , i = 1 , , N ,
u i | x = 0 = u i | x = 1 = 0 , i = 1 , , N .
Here ρ , u i , i = 1 , , N are respectively the density and velocity of the medium components, which are sought functions of time t [ 0 , T ] , 0 < T < + , and of the point x of the flow domain Ω = { x R | 0 < x < 1 } , v = 1 N i = 1 N u i is and average velocity of the medium, N 2 is the number of components, K = const > 0 , and constant viscosity coefficients ν i j , i , j = 1 , , N from a symmetric matrix N > 0 .
As can be seen, equations (1), (2) are a generalization of the one-dimensional isothermal Navier—Stokes equations for the multicomponent case. Equations (2) contain higher order derivatives of the velocities of all components, since, unlike the Navier—Stokes equations, in which viscosity is a scalar, in the multicomponent case, due to the composite structure of viscous stress tensors, the viscosities form a matrix whose entries are responsible for viscous friction. Diagonal elements are responsible for viscous friction within each component, while off-diagonal elements are responsible for friction between components. This does not allow us to automatically extend the known results for the Navier—Stokes equations to the multicomponent case. In the case of a diagonal viscosity matrix, the equations are connected possibly via the lower order terms only. In this paper, we consider a more complex case of an off-diagonal viscosity matrix. Stabilization of solutions to one-dimensional Navier—Stokes equations is studied in [5,6,7,8,9]. Unique solvability and asymptotic behavior (as t + ) of the solution to the considered one-dimensional equations of viscous compressible multicomponent media in the polytropic case is studied in [10,11,12]. Similar questions for related one-dimensional models of multicomponent media are discussed in [13,14,15,16,17,18,19,20].
It will be useful to use mass Lagrangian coordinates in the study of the equations (1), (2). Let us consider t and y ( x , t ) = 0 x ρ ( s , t ) d s as new independent variables. Then the system (1), (2) takes the form
ρ t + ρ 2 v y = 0 , v = 1 N j = 1 N u j ,
u i t + K ρ y = j = 1 N ν i j y ρ u j y , i = 1 , , N .
After that, the flow domain Ω transforms to the domain Ω ˜ = { y R | 0 < y < d } , where d = 0 1 ρ 0 ( x ) d x > 0 , and the initial and boundary conditions arrive at the form
ρ | t = 0 = ρ ˜ 0 ( y ) , u i | t = 0 = u ˜ 0 i ( y ) , i = 1 , , N ,
u i | y = 0 = u i | y = d = 0 , i = 1 , , N .
The main result of the article is the following assertion.
Theorem 1. 
Assume that ρ 0 W 2 1 ( 0 , 1 ) , ρ 0 > 0 , u 0 i W 2 1 ( 0 , 1 ) , i = 1 , , N . Then ρ d , u i 0 , i = 1 , , N , as t + , in the norm of the space W 2 1 ( 0 , 1 ) .
Proof. 
Let us obtain a priori estimates for the solution to the problem (1)–(4), which would be uniform in t, and from which the stabilization of the solution follows.

2. A priori estimates

Let us multiply the equations (2) by u i , summarize in i and integrate in y, thus, using (1), we obtain
1 2 d d t i = 1 N 0 1 ρ u i 2 d x + i , j = 1 N ν i j 0 1 u i x u j x d x = K i = 1 N 0 1 u i ρ x d x .
Since N > 0 , the second summand in the left-hand side of (9) satisfies the inequality 1
i , j = 1 N ν i j 0 1 u i x u j x d x C 1 ( N ) i = 1 N 0 1 u i x 2 d x .
Multiplying (1) by ln ρ ln d and integrating in x, we obtain for the right-hand side of (9)
K i = 1 N 0 1 u i ρ x d x = K N i = 1 N 0 1 ρ v x d x = K N d d t 0 1 ( ρ ln ρ ( ln d + 1 ) ρ + d ) d x .
Thus, using (10) and (11), from (9) we derive
1 2 d d t i = 1 N 0 1 ρ u i 2 d x + K N d d t 0 1 ( ρ ln ρ ( ln d + 1 ) ρ + d ) d x + C 1 i = 1 N 0 1 u i x 2 d x 0 .
In the Lagrangian coordinates ( t , y ) the formula (12) looks like
1 2 d d t i = 1 N 0 d u i 2 d y + K N d d t 0 d 1 ρ ( ρ ln ρ ( ln d + 1 ) ρ + d ) d y + C 1 i = 1 N 0 d ρ u i y 2 d y 0 .
From (13), after integration in t, we obtain the estimate
i = 1 N 0 d u i 2 d y + 0 d ρ ln ρ ( ln d + 1 ) ρ + d ρ d y + i = 1 N 0 t 0 d ρ u i y 2 d y d τ C 2 ,
where C 2 = C 2 ( C 1 , K , N , d , inf ( 0 , d ) ρ ˜ 0 , sup ( 0 , d ) ρ ˜ 0 , u ˜ 0 i L 2 ( 0 , d ) ) . Note that all terms in the left-hand side of (14) are nonnegative.
Let us rewrite the equations (6) as
1 N j = 1 N ν ˜ i j u j t + K N j = 1 N ν ˜ i j ρ y = 1 N y ρ u i y , i = 1 , , N ,
where ν ˜ i j , i , j = 1 , , N are the entries of the matrix N ˜ = N 1 . Summarizing (15) in i, we obtain
1 N i , j = 1 N ν ˜ i j u j t + K ˜ ρ y = y ρ v y ,
where K ˜ = K N i , j = 1 N ν ˜ i j > 0 . The equations (5) imply that
ρ v y = ln ρ t .
Substituting ρ v y from (17) into the equality (16), we arrive at the relation
2 ln ρ t y + K ˜ ρ y = 1 N i , j = 1 N ν ˜ i j u j t .
Let us multiply (18) by ln ρ y and integrate in y, then we obtain the equality
1 2 d d t 0 d ln ρ y 2 d y + K ˜ 0 d ρ ln ρ y 2 d y = 1 N i , j = 1 N ν ˜ i j 0 d u j t ln ρ y d y .
The right-hand side of (19) can be converted via integration by parts and using (8), (17):
1 N i , j = 1 N ν ˜ i j 0 d u j t ln ρ y d y = d d t 1 N i , j = 1 N ν ˜ i j 0 d u j ln ρ y d y + + 1 N i , j = 1 N ν ˜ i j 0 d ρ u j y v y d y .
Hence from (19) it follows that
1 2 d d t 0 d ln ρ y 2 d y + K ˜ 0 d ρ ln ρ y 2 d y = d d t 1 N i , j = 1 N ν ˜ i j 0 d u j ln ρ y d y + + 1 N i , j = 1 N ν ˜ i j 0 d ρ u j y v y d y .
Let us integrate (21) in t:
1 2 0 d ln ρ y 2 d y + K ˜ 0 t 0 d ρ ln ρ y 2 d y d τ = = 1 2 0 d 1 ρ ˜ 0 2 ρ ˜ 0 y 2 d y 1 N i , j = 1 N ν ˜ i j 0 d u j ln ρ y d y + + 1 N i , j = 1 N ν ˜ i j 0 d u ˜ 0 j ρ ˜ 0 ρ ˜ 0 y d y + 1 N i , j = 1 N ν ˜ i j 0 t 0 d ρ u j y v y d y d τ .
The last equality and (14) lead to the estimate
0 d ln ρ y 2 d y + 0 t 0 d ρ ln ρ y 2 d y d τ C 3 ,
where C 3 = C 3 ( C 2 , N ˜ , K ˜ , N , inf ( 0 , d ) ρ ˜ 0 , ρ ˜ 0 W 2 1 ( 0 , d ) , { u ˜ 0 i L 2 ( 0 , d ) } ) .
Let note that the equation (5) implies, for every t [ 0 , T ] , the existence of a point y ˜ ( t ) [ 0 , d ] such that
ρ ( t , y ˜ ( t ) ) = d .
Since
ln ρ ( t , y ) = ln ρ ( t , y ˜ ( t ) ) + y ˜ ( t ) y ( ln ρ ( t , s ) ) s d s ,
then by the Hölder inequality, taking into account (22) and (23), we deduce that
| ln ρ ( y , t ) | | ln d | + d ln ρ y L 2 ( 0 , d ) C 4 ( C 3 , d ) ,
and we immediately conclude that
0 < 1 C 5 ρ ( t , y ) C 5 < + , C 5 = C 5 ( C 4 ) .
Hence, from (14), (22) and (24) we obtain the estimate
i = 1 N 0 t 0 d u i y 2 d y d τ + 0 d ρ y 2 d y + 0 t 0 d ρ y 2 d y d τ C 6 ( C 2 , C 3 , C 5 ) .
In order to derive the next group of estimates, we multiply (6) by 2 u i y 2 and integrate in y, thus we obtain
1 2 d d t 0 d u i y 2 d y + j = 1 N ν i j 0 d ρ 2 u i y 2 2 u j y 2 d y = j = 1 N ν i j 0 d ρ y 2 u i y 2 u j y d y + + K 0 d ρ y 2 u i y 2 d y .
Let us summarize the equation (26) over i and integrate over t, then we obtain
1 2 i = 1 N 0 d u i y 2 d y + i , j = 1 N ν i j 0 t 0 d ρ 2 u i y 2 2 u j y 2 d y d τ = 1 2 i = 1 N 0 d u ˜ 0 i y 2 d y i , j = 1 N ν i j 0 t 0 d ρ y 2 u i y 2 u j y d y d τ + K i = 1 N 0 t 0 d ρ y 2 u i y 2 d y d τ .
The left-hand side of (27) satisfies the estimate
1 2 i = 1 N 0 d u i y 2 d y + i , j = 1 N ν i j 0 t 0 d ρ 2 u i y 2 2 u j y 2 d y d τ C 7 ( C 1 , C 5 ) i = 1 N 0 d u i y 2 d y + i = 1 N 0 t 0 d 2 u i y 2 2 d y d τ .
Due to (25) and the inequalities ( i = 1 , , N )
u i y C [ 0 , d ] 2 2 u i y L 2 ( 0 , d ) 2 u i y 2 L 2 ( 0 , d ) ,
the second summand in the right-hand side of (27) satisfies the estimate
i , j = 1 N ν i j 0 t 0 d ρ y 2 u i y 2 u j y d y d τ C 7 4 i = 1 N 0 t 0 d 2 u i y 2 2 d y d τ + C 8 ,
where C 8 = C 8 ( C 6 , C 7 , N , N ) . The third summand in the right-hand side of (27) can be estimated as follows:
K i = 1 N 0 t 0 d ρ y 2 u i y 2 d y d τ C 7 4 i = 1 N 0 t 0 d 2 u i y 2 2 d y d τ + C 9 ( C 6 , C 7 , K , N ) .
Thus, from (27), due to (28), (30) and (31), the estimate
i = 1 N 0 d u i y 2 d y + i = 1 N 0 t 0 d 2 u i y 2 2 d y d τ C 10 C 7 , C 8 , C 9 , u ˜ 0 i W 2 1 ( 0 , d )
follows. Finally, integrating (26) in t, we arrive at the inequalities
0 t d d τ 0 d u i y 2 d y d τ C 11 C 5 , C 6 , C 10 , N , K , N , i = 1 , , N .

3. Stabilization of the solution with an unlimited increase in time

From (25) and (33), the following convergence follow
u i y L 2 ( 0 , d ) 0 , i = 1 , , N
as t + . Differentiating (5) with respect to y and multiplying by ρ y , using (29), we arrive at the estimate
0 t d d τ 0 d ρ y 2 d y d τ C 12 C 5 , C 6 , C 10 , N .
Hence, as t + , we have
ρ y L 2 ( 0 , d ) 0 .
Thus, it is proved that (see (23)) in the norm of W 2 1 ( 0 , d )
u i 0 , i = 1 , , N , ρ d
as t + . It is easy to verify now that the same convergence takes place in the Euler variables in the norm of the space W 2 1 ( 0 , 1 ) . The theorem is proved. □

4. Conclusions

For the system of differential equations of isothermal viscous compressible multicomponent media with an non-diagonal, symmetric, and positive-definite viscosity matrix, an analysis is made of the asymptotic behavior (as t + ) of the solution to the initial-boundary value problem in the case of one spatial variable. The main difficulty was to obtain a priori estimates. To overcome it, the mass Lagrangian coordinates have been used. As a result, new a priori estimates are obtained and stabilization of the solution of the initial-boundary value problem is proved.

Funding

The research was supported by Russian Science Foundation (project 23–21–00381).

Conflicts of Interest

The author report no conflict of interest.

References

  1. Nigmatulin, R.I. Dynamics of multiphase media: v. 1–2; Hemisphere: New York, 1990. [Google Scholar]
  2. Rajagopal, K.L.; Tao, L. Mechanics of mixtures; World Scientific Publishing: New Jersey, 1995. [Google Scholar]
  3. Mamontov, A.E.; Prokudin, D.A. Viscous compressible multi–fluids: modeling and multi–D existence. Methods and Applications of Analysis 2013, 20, 179–195. [Google Scholar] [CrossRef]
  4. Mamontov, A.E.; Prokudin, D.A. Viscous compressible homogeneous multi–fluids with multiple velocities: barotropic existence theory. Siberian Electronic Mathematical Reports 2017, 14, 388–397. [Google Scholar]
  5. Zlotnik, A.A.; Ducomet, B. Stabilization rate and stability for viscous compressible barotropic symmetric flows with free boundary for a general mass force. Sbornik: Mathematics 2005, 196, 1745–1799. [Google Scholar] [CrossRef]
  6. Straskraba, I.; Zlotnik, A. Global properties of solutions to 1D–viscous compressible barotropic fluid equations with density dependent viscosity. Zeitschrift fur angewandte Mathematik und Physik 2003, 54, 593–607. [Google Scholar] [CrossRef]
  7. Zlotnik, A.A. ; Nguen, Zha Bao Properties and asymptotic behavior of solutions of some problems of one-dimensional motion of a viscous barotropic gas. Mathematical Notes 1994, 55, 471–482. [Google Scholar] [CrossRef]
  8. Antontsev, S.N.; Kazhikhov, A.V.; Monakhov, V.N. Boundary value problems in mechanics of nonhomogeneous fluids; North–Holland Publishing Co.: Amsterdam, 1990. [Google Scholar]
  9. Kazhikhov, A.V. Stabilization of solutions of the initial–boundary value problem for barotropic viscous fluid equations. Differentsial’nye Uravneniya 1979, 15, 662–667. [Google Scholar]
  10. Prokudin, D.A. Unique solvability of initial–boundary value problem for a model system of equations for the polytropic motion of a mixture of viscous compressible fluids. Siberian Electronic Mathematical Reports 2017, 14, 568–585. [Google Scholar]
  11. Prokudin, D.A. Global solvability of the initial boundary value problem for a model system of one-dimensional equations of polytropic flows of viscous compressible fluid mixtures. Journal of Physics: Conference Series 2017, 894, 012076. [Google Scholar] [CrossRef]
  12. Prokudin, D.A. On the stabilization of solutions to the initial–boundary value problem for the equations of dynamics of viscous compressible multicomponent media. Siberian Electronic Mathematical Report 2021, 18, 1278–1285. [Google Scholar] [CrossRef]
  13. Li, S. On one–dimensional compressible Navier–Stokes equations for a reacting mixture in unbounded domains. Zeitschrift fur angewandte Mathematik und Physik 2017, 68, 106. [Google Scholar] [CrossRef]
  14. Bresh, D.; Huang, X.; Li, J. Global weak solutions to one–dimensional non–conservative viscous compressible two–phase system. Communications in Mathematical Physics 2012, 309, 737–755. [Google Scholar] [CrossRef]
  15. Papin, A.A. On the uniqueness of the solutions of an initial boundary–value problem for the system of a heat–conducting two–phase mixture. Mathematical Notes 2010, 309, 594–598. [Google Scholar] [CrossRef]
  16. Papin, A.A.; Akhmerova, I.G. Solvability of the System of Equations of One-Dimensional Motion of a Heat-Conducting Two-Phase Mixture. Mathematical Notes 2010, 87, 230–243. [Google Scholar] [CrossRef]
  17. Zlotnik, A.A. Uniform estimates and stabilization of solutions to equations of one-dimensional motion of a multicomponent barotropic mixture. Mathematical Notes 1995, 58, 885–889. [Google Scholar] [CrossRef]
  18. Petrov, A.N. Well-posedness of initial-boundary value problems for one-dimensional equations of mutually penetrating flows of ideal gases. Dinamika Sploshnoi Sredy 1982, 56, 105–121. [Google Scholar]
  19. Antontsev, S.N.; Papin, A.A. Approximation methods of solving problems of two–phase filtration. Doklady Akademii Nauk SSSR 1979, 247, 521–525. [Google Scholar]
  20. Kazhikov, A.V.; Petrov, A.N. Well-posedness of the initial-boundary value problem for a model system of equations of a multicomponent mixture. Dinamika Sploshnoi Sredy 1978, 35, 61–73. [Google Scholar]

Notes

1
Everywhere, C k , k N stand for various positive constants depending on the values given in brackets, but independent of t.
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.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.

Altmetrics

Downloads

77

Views

21

Comments

0

Subscription

Notify me about updates to this article or when a peer-reviewed version is published.

Email

Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2025 MDPI (Basel, Switzerland) unless otherwise stated