We study the well-posedness and asymptotic behavior of solutions to the Cauchy problem of a three-dimensional sixth-order Cahn–Hilliard equation arising in oil-water-surfactant mixtures. First, by using the pure energy method and a standard continuity argument, we prove that there exists a unique global strong solution provided that the -norm of the initial data is sufficiently small. Moreover, we establish suitable negative Sobolev norm estimates and obtain the optimal decay rates of the higher-order spatial derivatives of the strong solution.
In 1990s, Gompper et. al. [7,8] introduced the following free energy functional:
with density given by
to describe the dynamics of phase transitions in ternary oil-water-surfactant systems, where describes a scalar parameter that is proportional to the local difference between oil and water concentrations; δ denotes the mobility and the second gradient energy coefficient; is the first gradient energy coefficient, which may be of arbitrary sign; and denotes the multiwell volumetric free energy density [18,19,24].
Let M be the mobility, and μ the chemical potential difference between the oil and water phases. Applying mass conservation, we obtain
with the mass flux j given by
Moreover, the chemical potential can be defined by the constitutive equation
where is the first variation of the function . Let the mobility , and we end up with the following sixth-order Cahn–Hilliard type equation:
Many papers have studied the initial boundary value problem of equation (1.2) from the point of view of global well-posedness. Pawlow and Zajaczkowski [18] assumed that the considered space is a bounded domain with a boundary of class , and proved the existence of a unique global smooth solution which depends continuously on the initial data. Moreover, by using the Bäcklund transformation and the Leray–Schauder fixed point theorem, Pawlow and Zajaczkowski [20] proved the global unique solvability of equation (1.2) in the Sobolev space under the assumption that the initial datum is in . Schimperna and Pawlow [24] discussed the existence, uniqueness, and parabolic regularization of a weak solution to the initial boundary value problem of equation (1.2) with a singular potential. The authors also supposed that the parameter , investigated a fourth-order system, and considered the existence of weak solutions under very general conditions by means of a fixed-point argument. In 2013, Schimperna and Pawlow [23] continued to study the fourth-order system with singular diffusion. The authors proved that for any final time T, the system admits a unique energy type weak solution, and for any , such a solution is classical. Later, based on Leray–Schauder’s fixed-point theorem and Campanato spaces, Liu and Wang [14] proved the existence of time-periodic solutions to the initial-boundary value problem of equation (1.2) in two spatial dimensions.
There are also some results on equation (1.2) with a viscous term. Pawlow and Zajaczkowski [19] applied Leray–Schauder’s fixed-point theorem and suitable estimates to establish the existence and uniqueness of a global in-time regular solution. Very recently, based on the estimates of weighted Sobolev spaces, Duan and Zhao [6] considered the existence of a global attractor for equation (1.2) with a viscous term in a 2D unbounded domain.
A Cauchy problem in mathematics asks for the solution of a partial differential equation that satisfies certain conditions that are given on a hypersurface in the domain. It is worth pointing out that the study of the Cauchy problem on higher-order nonlinear diffusion equations is also interesting. Many classical results are related to this topic (e.g., Caffarelli and Muler [1], Dlotko, Kania, and Sun [4], Cholewa and Rodriguez-Bernal [2] for the global existence of higher-order diffusion equations; Dlotko and Sun [5] and Savostianov and Zelik [22] for the global dynamics of higher-order diffusion equations). As far as we know, there is no reference concerning this aspect of the sixth-order Cahn–Hilliard equation arising in oil-water-surfactant mixtures. Can we establish some well-posedness results for the Cauchy problem of the sixth-order Cahn–Hilliard equation (1.2)? We will answer that question in this paper.
Due to its important applications in everyday life, the model of oil-water-surfactant mixtures is increasingly used in many fields. The Cauchy problem of equation (1.2) may be used to describe the handling process of offshore oil spills, where the oil pollution region may be large and irregular. To handle a spill, we must draw support from the surfactant, and consider the relations between oil, water, and surfactant.
Consider the Cauchy problem of the sixth-order viscous Cahn–Hilliard equation (1.2) in . The problem is stated as follows:
There are some different choices of the functions and . For example, on the basis of Landau–Ginzburg free energy, Gompper et al. [7,8] chose and as
where , , and are constants. It is easy to see that the free energy is nonconvex, and it has at least two ground states, 1 and . Based on the fourth-order gradient free energy, Pawlow and Zajaczkowski [20] supposed that
A specific free energy with composition-dependent gradient energy coefficient also will arise in the modeling of phase separation in polymers [3,19]. Known as Flory–Huggins–de Gennes energy, it has the form (1.1) with ,
and the singular coefficient
The sixth-order Cahn–Hilliard equation model is only a phenomenological model. Hence, various modifications, e.g., [12,13,15,16] and references therein, have been proposed to better capture the dynamical picture of the phase transition phenomena.
In this paper, we suppose that and satisfy (1.4). There are two cases to be considered: First, if the parameter , then equation (1.3)1 can be rewritten as
Second, if , then we would like to rewrite problem (1.3) as
For simplicity, we consider the following Cauchy problem in :
where , , , , and are constants. If , , and , then equation (1.8)1 is equivalent to (1.6). Moreover, if , , and , then we obtain equation (1.7).
Since we consider problem (1.8) in , the Laplacian () can be defined through the Fourier transform, i.e.,
where is the Fourier transform of f. Moreover, with the integer stands for the usual spatial derivatives of order l. If or l is not a positive integer, then stands for , as defined by (1.9). We use the notation to mean that for a universal constant that only depends on the parameters coming from the problem, and the indices N and s come from the regularity of the data. We also employ C as a positive constant, depending again on the initial data.
The first purpose of this paper is to consider the global well-posedness of solutions to problem (1.8). Let . Then, we have the following theorem.
Let. Suppose thatwhereis a positive constant to be determined in (
3.28
). Then there exists a unique global solutionsuch that for all,Moreover, iffor any integer, then
The asymptotic behavior of solutions is also an interesting topic in the study of the Cauchy problem of dissipative equations. In this paper, we show that the solutions of problem (1.8) satisfy some negative exponent decay rate.
Under the assumptions of Theorem
1.5
, and assuming thatfor some, then for all,and
Note that the Hardy–Littlewood–Sobolev theorem implies that for , with . Therefore, based on Theorem 1.5, we obtain the following corollary on the optimal decay estimates.
Under the assumptions of Theorem
1.5
, if we replace theassumption by, then the following decay estimate holds:where
The temporary decay rate of problem (1.8) is optimal because it is equivalent to the decay rate of the fourth-order generalized heat equation,
The main difficulties regarding the Cauchy problem of the sixth-order Cahn–Hilliard equation arising in oil-water-surfactant mixtures in are dealing with the second order nonlinear term and obtaining a negative Sobolev estimate to study the decay rate of solutions, since the principle part of problem (1.8) is a sixth-order linear term, and the nonlinear term is only second-order. Due to Sobolev’s embedding , we cannot control through . To overcome this difficulty, we borrow a fourth-order term from , and rewrite (1.3)1 as (1.8). Moreover, to consider the temporal decay rate of solutions of dissipative equations, one of the main tools is the standard Fourier splitting method. In this paper, since there exists the lower-order linear term on the right side of problem (1.8), it is difficult to use the Fourier splitting method to study the decay rate of solutions. By using the pure energy method [10,26,27], i.e., using a family of scaled energy estimates with minimum derivative counts and interpolations among them, we overcome the difficulty caused by the lower-order linear term, obtain suitable a priori estimates in the Sobolev space and negative Sobolev space (), and establish the optimal decay rate of problem (1.8) in .
The structure of this paper is organized as follows. In Section 2, we introduce some preliminary results, which are useful to prove our main results. Section 3 is devoted to prove the small data global well-posedness of problem (1.8). In the last section, we establish the time decay rate of solutions.
Preliminaries
First of all, we introduce the Kato–Ponce inequality which is of great importance in this paper.
We also give the following inequality on the -norm of , which can be seen as a corollary of Lemma 2.5. In fact, if we set , , , and in Lemma 2.5, the following inequality is easily to proved.
If, then u is (almost everywhere equal to) a continuous function and
The purpose of this subsection is to establish the a priori nonlinear energy estimates for problem (1.8).
We begin with the energy estimates including itself.
Suppose that all assumptions in Theorem
1.5
hold. Then for, we have
Applying to (1.8)1, multiplying the resulting identity by , and then integrating over by parts, we arrive at
We can estimate the first term of the right hand side of (3.2) as
For the second term of the right hand side of (3.2), we have
Moreover, the last term of the right hand side of (3.2) satisfies
Note that
Similarly, we also have
Plugging the estimates (3.3)–(3.7) into (3.2), we then obtain (3.1) and complete the proof. □
Local well-posedness
In this subsection, we prove the local well-posedness of solution in -norm.
We first construct the solution by solving iteratively the Cauchy problem:
for , where holds. One denote in short hand by and denote as . In the following, we shall show that is a Cauchy sequence in Banach space with suitable small. Then, by take the limit and a continuous argument, one propose to prove that is a global solution to Cauchy problem (1.8).
Suppose that all assumptions in Theorem
1.5
hold. There are constantsandsuch that if, then for each,is well defined andMoreover,is a Cauchy sequence in Banach space, the corresponding limit function denoted bybelongs towithandis a solution overto problem (
1.8
). Finally, for the Cauchy problem (
1.8
), there exists at most one solutioninsatisfying (
3.10
).
The inequality(3.9) will be proved by induction. By using the assumption at the initial step, we get , which means holds. Next, suppose that (3.9) holds for with small enough to be determined later, we are going to prove it also holds for . Hence, we need some energy estimates on . On the basis of (3.8)1, we obtain for and ,
Since it is trivial for the case , thus, we only need to consider . For the term , if , we estimate as follows
For the case , we can estimate as
For the term , if , we have
For the case of , we estimate as
Note that
and
Combining (3.15)–(3.17) together, if , we obtain
Next, for , if , we have
By Lemma 2.1, we consider the case of :
Now, summing up the estimates (3.11), (3.12), (3.13), (3.14), (3.18), (3.19) and (3.20), we arrive at
By taking time integration, we have
which from the inductive assumption implies
for any . Take and sufficiently small, such that
Then, by (3.22), we find that
for any . Therefore, (3.9) is true for if so for j, which implies (3.9) is proved for all .
Next, by using (3.21), we deduce that
for any . Therefore, due to (3.23), the time integral in the last inequality is finite, and hence is continuous in t for each . On the other hand, we also need to consider the convergence of the sequence . Taking the difference of (3.8)1 for j and , it yields that
Appealing to the same energy estimate as before, we get
which is equivalent to
Based on (3.23), by taking time integration, it holds that
Let be sufficiently small, then there exists a constant such that
for any . Hence, is a Cauchy sequence in the Banach space , the limit function
exists in , satisfies
Hence, (3.10) is proved. Finally, suppose that and are two solutions in satisfying (3.10). By using the same process as in (3.27) to prove the convergence of , we find that
for , which implies that . The proof of uniqueness is complete and thus the proof of Lemma 3.2 is complete too. □
Global well-posedness
In this subsection, we shall combine all the energy estimates that we have derived in the previous sections and the Sobolev interpolation to prove Theorem 1.5.
Summing up the estimates (3.2) of Lemma 3.1 from to 2, we easily obtain
where is a positive constant. Keeping in mind the condition (1.10), that is
one may deduce that there exists a time such that for all ,
Indeed, we argue (3.29) by contradiction. Suppose that (3.29) is not true. Then there exists a time (the first time) such that
Hence,
Since is one of the solutions for equation . It follows from (3.28) that
which implies
This is a contradiction with (3.30). In other words, (3.29) is true. It follows from (3.28) that for all , if ,
Therefore, we obtain the estimate (1.11). Summing up the estimates (3.2) of Lemma 3.1 from to N, we easily obtain
On the basis of the condition (1.10), we easily obtain
Then integrating it directly in time, we obtain (1.11). The proof of Theorem 1.5 is complete. □
In this subsection, we derive the evolution of the negative Sobolev norms of the solution.
Suppose that all assumptions in Theorem
1.5
hold. Then, for, we have
Applying to (1.8), multiplying the resulting identity by , integrating over by parts, we deduce that
We will estimate the three terms of the right hand side of (4.2) one by one. For the first term, we have
We estimate as
where we have used the fact , which was proved in Theorem 1.5. can be estimated as
where we have used the fact , which was proved in Theorem 1.5. We also have
and
where we have used the fact , which was proved in Theorem 1.5. Combining (4.3)–(4.7) together gives
Next, on the basis of the estimates (4.4)–(4.6), we estimate as
For , we have
Note that
For and , we have
and
where we have used the facts that and Sobolev’s embedding in (4.12) and (4.13). Similarly, we can easily obtain
Combining (4.11)–(4.14) together gives
Similarly, we obtain
It then follows from (4.15)–(4.16) that
Plugging (4.8), (4.9) and (4.17) into (4.2), we deduce (4.1). □
Define
For inequality (4.1), integrating in time, by the bound (1.11), we have
which implies (1.13) for , that is
Moreover, if , we may use Lemma 2.5 to have
Then, by this facts and (4.19), we get
Hence, for ,
Thus, we deduce from (3.34) the following inequality
Solving this inequality directly gives
which means (1.14) holds. Hence, we complete the proof of Theorem 1.6. □
Footnotes
Acknowledgements
This paper was supported by the Fundamental Research Funds for the Central Universities (grant No. N2005031). The author would like to thank the anonymous referee for his/her useful comments and suggestions, which led to the final version of this paper. The author also thank Prof. Hao Wu from Fudan University for his encourage and useful suggestions.
References
1.
L.A.Caffarelli and N.E.Muler, An bound for solutions of the Cahn–Hilliard equation, Arch. Ration. Mech. Anal.133 (1995), 129–144. doi:10.1007/BF00376814.
2.
J.W.Cholewa and A.Rodriguez-Bernal, On the Cahn–Hilliard equation in , J. Differential Equations253 (2012), 3678–3726. doi:10.1016/j.jde.2012.08.033.
3.
P.G.de Gennes, Dynamics of fluctuations and spinodal decomposition in polymer blends, J. Chem. Phys.72 (1980), 4756–4763. doi:10.1063/1.439809.
4.
T.Dlotko, M.B.Kania and C.Sun, Analysis of the viscous Cahn–Hilliard equation in , J. Differential Equations252 (2012), 2771–2791. doi:10.1016/j.jde.2011.08.052.
5.
T.Dlotko and C.Sun, Dynamics of the modified viscous Cahn–Hilliard equation in , Topol. Methods Nonlinear Anal.35 (2010), 277–294.
6.
N.Duan and X.Zhao, Global attractor for a class of sixth-order viscous Cahn–Hilliard equation in an unbounded domain, J. Dyn. Control Syst.25 (2019), 95–108. doi:10.1007/s10883-018-9403-1.
7.
G.Gompper and J.Goos, Fluctuating interfaces in microemulsion and sponge phase, Phys. Rev. E50 (1994), 1325–1335. doi:10.1103/PhysRevE.50.1325.
8.
G.Gompper and S.Zschocke, Ginzburg–Landau theory of oil-water-surfactant mixtures, Phys. Rev. A46 (1992), 4836–4851. doi:10.1103/PhysRevA.46.4836.
9.
L.Grafakos, Classical and Modern Fourier Analysis, Pearson Education, Inc., Prentice-Hall, 2004.
10.
Y.Guo and Y.Wang, Decay of dissipative equations and negative Sobolev spaces, Commu. Partial Differential Equations37 (2012), 2165–2208. doi:10.1080/03605302.2012.696296.
11.
T.Kato and G.Ponce, Commutator estimates and the Euler and Navier–Stokes equations, Comm. Pure Appl. Math.41 (1988), 891–907. doi:10.1002/cpa.3160410704.
12.
M.D.Korzec, P.L.Evans, A.Münch and B.Wagner, Stationary solutions of driven fourth- and sixth-order Cahn–Hilliard-type equations, SIAM J. Appl. Math.69 (2008), 348–374. doi:10.1137/070710949.
13.
M.D.Korzec, P.Nayar and P.Rybka, Global weak solutions to a sixth order Cahn–Hilliard type equation, SIAM J. Math. Anal.44 (2012), 3369–3387. doi:10.1137/100817590.
14.
C.Liu and Z.Wang, Time periodic solutions for a sixth order nonlinear parabolic equation in two space dimensions, Commun. Pure Appl. Anal.13 (2014), 1087–1104. doi:10.3934/cpaa.2014.13.1087.
15.
A.Miranville, Asymptotic behavior of a sixth-order Cahn–Hilliard system, Cent. Eur. J. Math.12 (2014), 141–154.
16.
A.Miranville, Sixth-order Cahn–Hilliard systems with dynamic boundary conditions, Math. Methods Appl. Sci.38 (2015), 1127–1145. doi:10.1002/mma.3134.
17.
L.Nirenberg, On elliptic partial differential equations, Annali della Scuola Normale Superiore di Pisa13 (1959), 115–162.
18.
I.Pawlow and W.M.Zajaczkowski, A sixth order Cahn–Hilliard type equation arising in oil-water-surfactant mixtures, Comm. Pure Appl. Anal.10 (2011), 1823–1847. doi:10.3934/cpaa.2011.10.1823.
19.
I.Pawlow and W.M.Zajaczkowski, On a class of sixth order viscous Cahn–Hilliard type equations, Discrete Contin. Dyn. Syst. Ser6 (2013), 517–546. doi:10.3934/dcdss.2013.6.517.
20.
I.Pawlow and W.M.Zajaczkowski, The global solvability of a sixth order Cahn–Hilliard type equation via the Bäcklund transformation, Commun. Pure Appl. Anal.13 (2014), 859–880. doi:10.3934/cpaa.2014.13.859.
21.
J.C.Robinson, J.L.Rodrigo and W.Sadowski, The Three-Dimensional Navier–Stokes Equations, Cambridge University Press, Cambridge, 2016.
22.
A.Savostianov and S.Zelik, Global well-posedness and attractors for the hyperbolic Cahn–Hilliard–Oono equation in the whole space, Math. Models Methods Appl. Sci.26 (2016), 1357–1384. doi:10.1142/S0218202516500329.
23.
G.Schimperna and I.Pawlow, A Cahn–Hilliard equation with singular diffusion, J. Differential Equations254 (2013), 779–803. doi:10.1016/j.jde.2012.09.018.
24.
G.Schimperna and I.Pawlow, On a class of Cahn–Hilliard models with nonlinear diffusion, SIAM J. Math. Anal.45 (2013), 31–63. doi:10.1137/110835608.
25.
E.M.Stein, Singular Integrals and Differentiability Properties of Functions, Princeton Unversity Press, Princeton, NJ, 1970.
26.
Z.Tan, W.Wu and J.Zhou, Global existence and decay estimate of solutions to magneto-micropolar fluid equations, J. Differential Equations266 (2019), 4137–4169. doi:10.1016/j.jde.2018.09.027.
27.
Y.Wang, Decay of the Navier–Stokes–Poisson equations, J. Differential Equations253 (2012), 273–297. doi:10.1016/j.jde.2012.03.006.