Advances in Mathematical Physics
Volume 2018 (2018), Article ID 2028657, 11 pages
https://doi.org/10.1155/2018/2028657

Research Article

Strong Solutions for the Fluid-Particle Interaction Model with Non-Newtonian Potential

Yukun Song, Heyuan Wang, and Fengming Liu

College of Science, Liaoning University of Technology, Jinzhou 121001, China

Correspondence should be addressed to Yukun Song; songyukun8@163.com

Received 19 January 2018; Revised 2 May 2018; Accepted 11 June 2018; Published 17 July 2018

Academic Editor: Manuel De León

Copyright © 2018 Yukun Song et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This paper deals with a mathematical fluid-particle interaction model used to describing the evolution of particles dispersed in a viscous compressible non-Newtonian fluid. It is proved that the initial boundary value problems with vacuum admits a unique local strong solution in the dimensional case. The strong nonlinearity of the system brings us difficulties due to the fact that the viscosity term and non-Newtonian gravitational potential term are fully nonlinear.

1. Introduction

Fluid-particle interaction model arises in many practical applications in science and engineering [1–4] and is of primarily importance in the sedimentation analysis of disperse suspensions of particles in fluids. We focus on the fluid-particle interaction model that describes the evolution of particles dispersed in a viscous fluid.

Carrillo and Goudon first derived a fluid-particle interaction system by formal asymptotics from a mesoscopic description (see [5]). In the microscopic description, the cloud of the particles is related to its distribution function , which is the solution to a dimensionless Vlasov-Fokker-Planck equation. On the other hand, the fluid is described by its density and its velocity field . We assume that the fluid is compressible and isentropic, then solves the compressible Euler equations for the inviscid case or the Navier-Stokes equations for the viscous case, respectively.

In [5], for the inviscid case, the coupling between the kinetic and the fluid equations is obtained through the friction forces that the fluid and the particles exert mutually. The friction force is assumed to follow the Stokes law and thus is proportional to the relative velocity of the fluid and the particles given by . Furthermore, both phases are affected by external forces, which are supposed to derive from a time independent potential . The system was given as follows:

where

are some related dimensionless parameters, are the pressure, the mass density of particles and fluid, respectively. Setting , , , with , then (1) becomes

Finally, letting , then (3) converges to system

For the viscous case, with the dynamic viscosity terms taken into consideration, then will be have an additional term of . In [6], the viscous stress tensor is assumed to satisfy Newton’s Law for viscosity which requires that

where and are constant viscosity coefficients satisfying

Then

and thus the system turns into the following equation

Moreover, if the influence of gravitational potential was taken into consideration, there will be a Poisson equation coupled to the above system, as for the Navier-Stokes-Poisson equation in [7].

Carrillo et al. obtained the global existence and asymptotic behavior of the weak solutions and stability properties to (8). Subsequently, Fang et al. [8] studied the existence of global classical solutions in dimension one. In [9, 10], Balew and Trivisa obtained the existence of global weak solutions and weakly dissipative solutions by entropy method in dimension three. The two-phase flow hydrodynamic models have been proposed in [3]. For some mathematical results on the Navier-Stokes coupled equations such as the nematic liquid crystal flows models where viscous effects are included, for more details, we refer to [11–13] and references therein.

On the other hand, as we know, the viscous stress tensor is depends on the rate of strain , where

If the stress and rate of strain satisfy the following linear relation

then the fluid is called Newtonian. The coefficient of proportionality is called the viscosity coefficient, and it is a characteristic material quantity for the fluid concerned, which in general depends on density, temperature, and pressure. The governing equations of motions of them will be the Navier-Stokes equations. If the relation is not linear, the fluid is called non-Newtonian. Examples of non-Newtonian fluids are molten plastics, polymer solutions, dyes, varnishes, suspensions, adhesives, paints, greases, paper pulp, and biological fluids like blood. The simplest model of the stress-strain relation for such fluids given by the power laws, which states that

for (see [14]). Ladyzhenskaya (see [15]) proposed a special form for on the incompressible model:

These models are called

For , if then it is a pseudo-plastic fluid, and if then it is a dilatant fluid (see [14]). In the view of physics, the model captures the shear thinning fluid for the case of , and captures the shear thickening fluid for the case of .

Followed by the Ladyzhenskaya model, in this paper, we investigate the compressible non-Newtonian fluid-particle interaction model in one-dimensional case, then system (8) changes to be

with the initial and boundary conditions

and the no-flux condition for the density of particles

where is a bounded interval, . denote the fluid density, velocity, and the density of particle in the mixture, respectively. is the pressure where , stands for the non-Newtonian gravitational potential, and the given function denotes the external potential. is the viscosity coefficient and , , , are constants.

As to the non-Newtonian fluids, there has been much research both theoretically and experimentally, see ([15–23]). Indeed, we have investigated the existence results of solutions for and in the absence of term for (14) in [22, 23]. However, the influence of non-Newtonian gravitational potential for a practical model was not taken into consideration there. To our knowledge, there seems very few mathematical results for the case of the fluid-particle interaction model systems with non-Newtonian gravitational potential, even in dimension one. The existence results to problem (14)-(16) when which describes the motion of the compressible viscous isentropic gas flow is driven by a non-Newtonian gravitational force is still open up to now. We are interested in the existence and uniqueness of strong solutions on a one-dimensional bounded domain. In fact, the strong nonlinearity of (14) brings us new difficulties in getting the upper bound of and the method used in [8] does not suitable for us. Motivated by Cho etal’s [24, 25] work on Navier-Stokes equations, we establish local existence and uniqueness of strong solutions by the iteration techniques.

Throughout the paper we assume that . In the following sections, we will use simplified notations for standard Sobolev spaces and Bochner spaces, such as , , .

1.1. Main Results

Theorem 1. Let be a positive constant and , and assume that the initial data satisfy the following conditions:

and

for some , where is a constant. Then there exist a time and a unique strong solution to (14)-(16) such that

2. A Priori Estimates for Smooth Solutions

In this section, we will prove the local existence of strong solutions. Provided that is a smooth solution of (14)-(16) and , where is a positive number. We denote by and introduce an auxiliary function

Then we estimate each term of in terms of some integrals of and apply arguments of Gronwall-type and thus prove that is locally bounded.

2.1. Estimate for

First we need to do the following estimates. Using , we rewritten as

By virtue of

Then

Taking the above inequality by norm, we get

We deal with the term of .

Multiplying by and integrating over , we obtain

then we have

Hence, we deduce that

From , taking it by norm, we get

Multiplying by , integrating over , we have

Integrating by parts, using Sobolev inequality, we deduce that

Differentiating with respect to , and multiplying it by , integrating over , using Sobolev inequality, we have

From (30) and (31), by Gronwall’s inequality, it follows that

Besides, we can also get the following estimates. By using ,

From , we have

then

Differentiating with respect to time , multiplying it by , integrating over to , and using Young’s inequality, we have

thus, we get

where is a positive constant, depending only on .

2.2. Estimate for

Multiplying (21) by , integrating (by parts) over , we have

We deal with each term as follows:

From we get

Substituting the above into (38), and using Young’s inequality, we obtain

Using , we have

Combining (42)-(43), yields

where is a positive constant, depending only on .

2.3. Estimate for and

Multiplying by , integrating the resulting equation over , and using the boundary conditions (16), Young’s inequality, we have

Multiplying by , integrating (by parts) over , and using the boundary conditions (16), Young’s inequality, we have

Differentiating with respect to , multiplying the resulting equation by , and integrating (by parts) over , we get

Combining (45)-(47), we get

2.4. Estimate for

Differentiating (21) with respect to , multiplying the result equation by , and integrating it over with respect to , we have

Note that

Combining (40), (49) can be rewritten into

By using Sobolev inequality, Hölder inequality and Young’s inequality, (35), (37), we estimate each term of as follows:

Substituting into (51), and integrating over on the time variable, we have

To obtain the estimate of , we need to estimate . Multiplying (21) by and integrating over , we get

According to the smoothness of , we obtain

Therefore, taking a limit on in (53), as , we conclude that

where is a positive constant, depending only on .

Combining the estimates of (27), (28), (32), (33), (44), (48), (56) and the definition of , we conclude that

where are positive constants, depending only on . This means that there exist a time and a constant , such that

3. Proof of the Main Theorem

In this section, our proof will be based on the usual iteration argument and some ideas developed in [24, 25]. We construct the approximate solutions, by using the iterative scheme, inductively, as follows: first define and assuming that was defined for , let be the unique smooth solution to the following problems:

with the initial and boundary conditions

where

We directly construct approximate solutions of the problem (59)–(63). More precisely, we first find from (59) and (63) with smooth function , i.e.,

and

It follows from the classical linear hyperbolic theory that there is a unique solution on this above initial problem. Using the method of characteristics, we have

By (68) and (69), we have

Using (67), then

which means

Next, combining the classical stableness results of the elliptic equation, the existence of can be obtained by (61) and (63), then by (62) and (63) we get . The last, with being given, by virtue of (72), from (60) and (63), according to the classical theorem of quasi-linear parabolic equation (see [17], Chapter VI, Theorem 5.2), there exists a unique smooth solution . With the process, the nonlinear coupled system has been deduced into a sequence of decoupled problems and each problem admits a smooth solution. And the following estimates hold:

where is a generic constant depending only on , but independent of .

Next, we have to prove that the approximate solution converges to a solution to the original problem (14) in a strong sense. To this end, let us define

then we can verify that the functions satisfy the system of equations

Multiplying (75) by , integrating over , and using Young’s inequality, we obtain

where is a positive constant, depending on and for all and .

Multiplying (76) by , integrating over , and using Young’s inequality, we obtain

Let

and then

We estimate the second term of (80) as follows:

Similarly, we have

and then we have

Substituting (83) into (80) and using Young’s inequality, we have

where , for all and . Using (73) we derive

Multiplying (78) by , integrating over , and using (73) and Young’s inequality, we have

Collecting (79), (86), and (88), we obtain

where depends only on and , for all and . Using (73), we have

Integrating (89) over with respect to t, using Gronwall’s inequality, we have

from the above recursive relation, choose and such that , using Gronwall’s inequality, we deduce that

Since all of the constants do not depend on , as , we conclude that sequence converges to a limit in the following convergence:

and there also holds

For each small , let , is a mollifier on , and is a smooth solution of the boundary value problem:

where and satisfies , .

We deduce that is a solution of the following initial boundary value problem:

where .

By the proof of Lemma 2.3 in [20], there exists a subsequence of , as , in , in , Hence, satisfies (18) of Theorem 1. By virtue of the lower semi-continuity of various norms, we deduce that satisfies the following uniform estimate:

where is a positive constant, depending only on .

The uniqueness of solution can be obtained by the same method as the above proof of convergence; we omit the details here. This completes the proof.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (nos. 11526105; 11572146) and the Funds of Education Department of Liaoning Province (JQL201715411).

References

A. A. Anthony, “Kiva-3V Release2, Improvements to Kiva-3V,” Tech. Rep., Los Alamos National Laboratory, 1999.

S. Berres, R. Bürger, K. H. Karlsen, and E. M. Tory, “Strongly degenerate parabolic-hyperbolic systems modeling polydisperse sedimentation with compression,” SIAM Journal on Applied Mathematics, vol. 64, no. 1, pp. 41–80, 2003.

L. Boudin, L. Desvillettes, and R. Motte, “A modeling of compressible droplets in a fluid,” Communications in Mathematical Sciences, vol. 1, no. 4, pp. 657–669, 2003.

I. Vinkovic, C. Aguirre, S. Simoëns, and M. Gorokhovski, “Large eddy simulation of droplet dispersion for inhomogeneous turbulent wall flow,” International Journal of Multiphase Flow, vol. 32, no. 3, pp. 344–364, 2006.

J. A. Carrillo and T. Goudon, “Stability and asymptotic analysis of a fluid-particle interaction model,” Communications in Partial Differential Equations, vol. 31, no. 7-9, pp. 1349–1379, 2006.

J. A. Carrillo, T. Karper, and K. Trivisa, “On the dynamics of a fluid-particle interaction model: the bubbling regime,” Nonlinear Analysis. Theory, Methods & Applications. An International Multidisciplinary Journal, vol. 74, no. 8, pp. 2778–2801, 2011.

J. Yin and Z. Tan, “Global existence of strong solutions of Navier-Stokes-Poisson equations for one-dimensional isentropic compressible fluids,” Chinese Annals of Mathematics, Series B, vol. 29, no. 4, pp. 441–458, 2008.

D. Fang, R. Zi, and T. Zhang, “Global classical large solutions to a 1D fluid-particle interaction model: the bubbling regime,” Journal of Mathematical Physics, vol. 53, no. 3, 033706, 21 pages, 2012.

J. Ballew and K. Trivisa, “Suitable weak solutions and low stratification singular limit for a fluid particle interaction model,” Quarterly of Applied Mathematics, vol. 70, no. 3, pp. 469–494, 2012.

J. Ballew and K. Trivisa, “Weakly dissipative solutions and weak-strong uniqueness for the Navier-Stokes-SMOluchowski system,” Nonlinear Analysis. Theory, Methods & Applications. An International Multidisciplinary Journal, vol. 91, pp. 1–19, 2013.

T. Huang, C. Wang, and H. Wen, “Strong solutions of the compressible nematic liquid crystal flow,” Journal of Differential Equations, vol. 252, no. 3, pp. 2222–2265, 2012.

H. Wen and C. Zhu, “Blow-up criterions of strong solutions to 3D compressible Navier-Stokes equations with vacuum,” Advances in Mathematics, vol. 248, pp. 534–572, 2013.

S. Jiang and Y. Ou, “A blow-up criterion for compressible viscous heat-conductive flows,” Acta Mathematica Scientia B, vol. 30, no. 6, pp. 1851–1864, 2010.

G. Bohme, Non-Newtonian Fluid Mechanics, vol. 31 of North-Holland Series in Applied Mathematics and Mechanics, North-Holland, Amsterdam, The Netherlands, 1987.

O. A. Ladyzhenskaya, “New equations for the description of viscous incompressible fluids and solvability in the large of the boundary value problems for them,” in Boundary Value Problems of Mathematical Physics, vol. V, American Mathematical Society, Providence, RI, USA, 1970.

O. A. Ladyzhenskaya, The Mathematical Theory of Viscous Incompressible Flow, Gordon and Breach, New York, NY, USA, 1969.

O. A. Ladyzenskaja, V. A. Solonnikov, and N. N. Uralceva, Linear and Quasilinear Equations of Parabolic Type, Translations of Mathematical Monographs, American Mathematical Society, Providence, RI, USA, 1968.

J. Málek, J. Necas, Rokyta. M., and M. Ružicka, Weak and Measure-v alued Solutions to Evolutionary PDEs, Chapman-Hall, New York, NY, USA, 1996.

A. E. Mamontov, “Global regularity estimates for multidimensional equations of compressible non-Newtonian fluids,” Annali dell'Università di Ferrara. Nuova Serie. Sezione VII. Scienze Matematiche, vol. 46, pp. 139–160, 2000.

H. Yuan and X. Xu, “Existence and uniqueness of solutions for a class of non-Newtonian fluids with singularity and vacuum,” Journal of Differential Equations, vol. 245, no. 10, pp. 2871–2916, 2008.

O. Rozanova, “Nonexistence results for a compressible non-Newtonian fluid with magnetic effects in the whole space,” Journal of Mathematical Analysis and Applications, vol. 371, no. 1, pp. 190–194, 2010.

Y. Song, H. Yuan, Y. Chen, and Z. Guo, “Strong solutions for a 1D fluid-particle interaction non-newtonian model: the bubbling regime,” Journal of Mathematical Physics, vol. 54, no. 9, Article ID 091501, 12 pages, 2013.

Y. Song, H. Wang, Y. Chen, and Y. Zhang, “The strong solutions for a class of fluid-particle interaction non-Newtonian models,” Boundary Value Problems, Paper No. 108, 17 pages, 2016.

Y. Cho, H. J. Choe, and H. Kim, “Unique solvability of the initial boundary value problems for compressible viscous fluids,” Journal de Mathématiques Pures et Appliquées, vol. 83, no. 2, pp. 243–275, 2004.

Y. Cho and H. Kim, “Existence results for viscous polytropic fluids with vacuum,” Journal of Differential Equations, vol. 228, no. 2, pp. 377–411, 2006.