The aim of this paper is to study higher-order Caginalp phase-field systems based on the Maxwell–Cattaneo law, instead of the classical Fourier law. More precisely, one obtains well-posedness results, as well as the existence of finite-dimensional attractors.
The Caginalp phase-field model
proposed in [4], has been extensively studied (see, e.g., [3,5,6,8,9,11,15–20,22] and the references therein). Here u denotes the order parameter, θ the relative temperature (defined as , where is the absolute temperature and the equilibrium melting temperature) and f the derivative of a double-well potential F (a typical choice of potential is , hence the usual cubic nonlinear term ). Furthermore, all physical constants have been set equal to one. This system models, e.g., melting-solidition phenomena in certain classes of materials.
These equations can be derived as follows: one introduces the (total Ginzburg-Landau) free energy
where Ω is the domain occupied by the system (one assumes here that it is a bounded and regular domain of , with boundary Γ). One then defines the enthalpy H as
where ∂ denotes a variational derivative, which gives
The governing equations for u and θ are then given by (see [5])
where q is the thermal flux vector. Assuming the classical Fourier law
one finds (1.1) and (1.2).
Now, a drawback of the Fourier law is the so-called “paradox of heat conduction”, namely, it predicts that thermal signals propagate with infinite speed, which, in particular, violates causality (see [12]). One possible modification, in order to correct this unrealistic feature, is the Maxwell–Cattaneo law
In that case, it follows from (1.7) that
hence the following equation for the temperature:
Integrating (1.11) between 0 and t (the terme in (1.11) would be very difficult to handle from a mathematical point of view) and setting
where α is called thermal displacement variable (here, is a priori fixed arbitrarily; see also below), one finally obtains
where
For simplification, one takes in (1.13): ; one notes that it follows from (1.14) that, if and are described, then, providing suitable assumptions on the initial data, one can always choose such that in (1.13), also one notes that can be obtained from (1.1), hence (1.13) becomes
Caginalp G. and Esenturk E. proposed in [7] higher-order phase-field models in order to account for anisotropic interfaces (see also [23] for other approaches which, however, do not provide an explicit way to compute the anisotropy). More precisely, these authors proposed the following modified (total) free energy
where, for ,
and, for ,
(one agrees that ). Noting that , this then yields the following evolution equation for the order parameter u:
In particular, for (anisotropic Caginalp phase-field system), one has an equation of the form
and, for (fourth-order anisotropic Caginalp phase-field system), one has an equation of the form
Cherfils L., Miranville A. and Peng S. studied in [10] the corresponding higner-order isotropic equation (without the coupling with the temperature), namely the equation
where
endowed with the Dirichlet/Navier boundary conditions
The aim of this paper is to study the model consisting of the higher-order anisotropic equation (1.17) and the temperature equation (1.15). In particular, one obtains well-posedness results, as well as the existence of finite-dimensional attractors.
Setting of the problem
One considers in this section the following initial and boundary value problem, for :
One assumes that
and one introduces the elliptic operator defined by
where is the topological dual of . Furthermore, denotes the usual -scalar product, with associated norm ; more generally, one denotes by the norm on the Banach space X. One can note that
is bilinear, symmetric, continuous and coercive, so that
is indeed well defined. It then follows from elliptic regularity results for linear elliptic operators of order (see [1,2]) that is a strictly positive, selfadjoint and unbounded linear operator with compact inverse, with domain
where, for ,
One further notes that and, for ,
One finally notes that (see, e.g., [24]) (resp., ) is equivalent to the usual -norm (resp., -norm) on (resp., ).
Similarly, one can define the linear operator
which is a strictly positive, selfadjoint and unbounded linear operator with compact inverse, with domain
where, for ,
Furthermore, and, for ,
Besides, (resp., ) is equivalent to the usual -norm (resp., -norm) on (resp., ).
Having this, one rewrites (2.1) as
where and, for ,
As far as the nonlinear term f is concerned, one assumes that
where . In particular, the usual cubic nonlinear term satisfies these asumptions.
Throughout the paper, the same letters c, and denote (generally positive) constants which may vary from line to line. Similarly, the same letter Q denotes (positive) monotone increasing (with respect to each argument) and continuous functions which may vary from line to line.
A priori estimates
The estimates derived in this section are formal, but they can easily be justified within a Galerkin approximation.
One multiplies (2.7) by , one has, integrating over Ω and by parts,
where and, for ,
(one notes that is not necessarily nonnegative). One can note that, owing to the interpolation inequality
there holds
This yields, employing (2.12),
whence
noting that, owing to Young’s inequality,
One then multiplies (2.7) by u and one has, owing to (2.11) and the interpolation inequality (3.3),
Summing finally (3.9) and times (3.12), where is chosen small enough so that
one obtains, owing to the interpolation inequality (3.3), a differential inequality of the form
where
satisfies
If follows from (3.13), (3.14) and Gronwall’s lemma that
and
for , given.
One now distinguishes between the cases and .
First case:
Setting , , and 3, one then considers the initial and boundary value problem:
Multipying (3.17) by , one has
which yields, owing to (2.9),
One now differentiates (3.17) with respect to time to find, owing to (3.18),
together with the boundary condition
Multiplying (3.22) by , one has, owing to (2.9),
Summing (3.13) (for ), times (3.21) and times (3.24), where are chosen small enough, one obtains a differential inequality of the form
where
satisfies
One finally rewrites (3.17) as an elliptic equation, for fixed,
Multiplying (3.27) by , one finds, owing to (2.9) and a classical elliptic regularity result,
In particular, it follows from (3.25) (for ) and Gronwall’s lemma that
where one has used the continuous embedding to deduce that
and the fact that
which yields
Combining (3.28) and (3.29), one finally finds, in view of (3.26),
Second case:
One multiplies (2.7) by and one obtains, owing to the interpolation inequality (3.3),
It follows from the continuity of f and F, the continuous embedding (one recalls that ) and (3.15) that
so that
One then differentiates (2.7) with respect to time and noting that
one obtains
One can note that, if and , then and
One multiplies (3.36) by and one finds, owing to (2.9) and the interpolation inequality (3.3),
Summing (3.13) (for ), (3.35) and times (3.40), where is chosen small enough, one obtains a differential inequality of the form
where
satisfies
In particular, it follows from (3.41)–(3.42) that
and
Next, one rewrites (2.7) as an elliptic equation, for fixed,
Multiplying (3.45) by , one has, owing to the interpolation inequality (3.3),
hence, owing to (3.34) and (3.43),
One now multiplies (2.2) by and by times , where is chosen small enough enough so that
one integrates over Ω and by parts, and one sums the resulting relations to obtain
In particular,
This yields, owing to (3.44) () and the uniform Gronwall’s lemma that
Furthermore, on , if one assumes that and , one deduces from (3.48), (3.44) and Gronwall’s lemma an -estimate on on which, combined with (3.49) and (3.47), gives an -estimate on , for all times, namely
One can note that one has estimate which is similar to (3.49) when .
The dissipative semigroup
This section is devoted to the existence and uniqueness of soltions to the problem (2.1)–(2.4), and of bounded absording sets to the semigroup.
One first has the following theorem.
One assumes that, withwhen, and. Then, (
2.1
)–(
2.4
) possesses a unique solutionsuch that,,,,,andwhereand, for,
If one further assumes thatand, thenand
Existence:
The proofs of existence and regulararity in (i) and (ii) follow from the a priori estimates derived in the previous section and, e.g., a standard Galerkin scheme.
Uniqueness:
Let now and be two solutions to (2.1)–(2.3) with initial data and , respectively. One sets
and
Then, satisfies
Multiplying (4.1) by u, one obtains, owing to (2.9) and the interpolation inequality (3.3),
Summing then (4.5) and times (4.6), where is small enough, one has a differential inequality of the form (one notes that )
where
satisfies
It follows from (4.7), (4.8) and Gronwall’s lemma that
hence the uniqueness, as well as the continuous dependence with respect to the initial data in the -norm. □
One sets . It follows from Theorem 4.1 that one can define the family of solving operators
(i.e., (identity operator) and , t, ) which are continuous with respect to the -topology. One then deduces from (3.15) the
The semigroupis dissipative in Φ, in the sense that it possesses a bounded absording set(i.e., for everybounded, there existssuch thatimplies).
It is easy to see that one can assume, without loss of generality, that is positively invariant by , i.e. , for every .
Existence of the global attractor
This section is devoted to the existence of the global attractor and one takes the coefficients .
One now states the following existence result of the global attractor.
One assumes that the assumptions of Theorem
4.1
, and (
2.8
) and (
2.9
) hold. Then, the semigrouppossesses the global attractorwhich is compact in Φ and bounded in.
The equality (2.2) is a damped wave equation, then one uses a semigroup decomposition argument consisting in splitting the semigroup , , into the sum of two operators families: , where operators go to zero as t tends to infinity while operators are compacts. In order to do that, one writes
where is the solution to
and satisfies
Multiplying (5.1) by and (5.2) by , and summing the resulting relations, one obtains
where
Summing (5.9) and times (5.10), where is chosen small enough so that
and one obtains adifferential inequality of the form
where
satisfies
It follows from (5.11), (5.12) and Gronwall’s lemma that
One can see that tends to zero when t tends to infinity.
Now, one considers (5.5)–(5.8). One multiplies (5.5) by and (5.6) by , one integrates over Ω and by parts, and summing the resulting relations, one then writes
One then has
Furthemore, one has, owing to Hölder’s inequality (),
Analogously, one finds
One deduce from (5.14) and (5.15) an inequality of the form
One then multiplies (5.6) by ψ and integrating over Ω and by parts, one has
Summing (5.16) and times (5.17), where is chosen small enough so that
one obtains a differential inequality of the form
where
satisfies
Multiplying (5.5) by and (5.6) by , integrating over Ω and by parts, and summing the two resulting relations to obtain, using Hölder’s inequality,
where
One multiplies then (5.6) by , to have
One deduces summing (5.20) and times (5.21), where is chosen small enough so that
Finally, summing (5.18) and (5.22), in particular, one finds
where
satisfies
It follows from the assumptions of the theorem that (see [14]) , for every , and
for some functions Q. One thus deduces from (5.23), (5.24) and (5.25) that
Hence, the operator is asymptotically compact in the sense of Kuratowski measure of noncompactness (see [20]), which completes the proof. □
One recalls that the global attractor is the smallest (for the inclusion) compact set of the phase space which is invariant by the flow (i.e., , for ) and attracts all bounded sets of initial data as time goes to infinity; it thus appears as a suitable object in view of the study of the asymptotic behavior of the system. Furthermore, the finite-dimensionality means, very roughly speaking, that, even though the initial phase space has infinite dimension, the reduced dynamics can be described by a finite number of parameters. One refers the interested reader to, e.g., [21,24] for more details and discussions on this.
Now that the existence of the global attractor is proved, a naturel question is to know it this attractor has finite dimension in terms of the Hausdorff or fractal dimension. That is the aim of the next section.
Existence of exponential attractors
In this section, one derives several estimates on the difference of solutions of problem (2.1)–(2.3) which are of fundamental significance for the study of exponential attracors.
One starts by stating an abstract result that will be useful in what follows (see [21]).
Let V and H be two Banach spaces such that V is compactly embedded into H andbe a semigroup acting on a closed subset X of H. One assumes that
withd is continuous,,as, and
is uniformly Hölder continuous on, for every, andbounded.
Thenpossesses an exponential attractor(contains the global attractor) on X.
In order to get the existence of exponential attractors in this case, Theorem 6.1 can be used. One has the following result.
Assume that the assumptions of Theorem
4.1
hold. Then, the semigroup,correspoding to equations (
2.1
)–(
2.4
) defined from Φ into itself satisfies a decomposition as in Theorem
6.1
.
Let now and be two solutions to (2.1)–(2.3) with initial data and , respectively. One sets
Then, satisfies
where
Now one decomposes the solution as follows:
where
are solutions to
and
respectively. Repeating for (6.5)–(6.8) estimates which led us to (5.11), one writes
where
satisfies as (5.12)
It follows from (6.13), (6.14) and Gronwall’s lemma that
Setting , , d is continuous and as .
One now considers (6.9)–(6.12). One multiplies (6.9) by and (6.10) by , one integrates over Ω and by parts, and summing the resulting relations, one then writes
From to (2.10) and owing to the Hölder’s inequality, one has ()
Analogously, one finds
One deduces from (6.16) and (6.17) that
One then multiplies (6.10) by ξ and integrating over Ω and by parts, one has
Summing (6.18) and times (6.19), where is chosen small enough so that
one obtains a differential inequality of the form
where
satisfies
Next, multiplying (6.9) by and (6.10) by , integrating over Ω and by parts, and summing the two resulting relations to obtain
Finally, summing (6.20) and (6.26), one obtains a differential inequality of the form
where
satisfies
From previous estimates, especially (3.10)–(3.11), one can see that
where Q is a monotonic function that depends only on time. Therefore, thanks to (6.27), (6.28) and the Gronwall’s lemma, one has the following inequality:
where
and
The function h is well continuous, which completes the proof of the proposition. □
Thus, the first condition of Theorem 6.1 is proved and it is sufficient to prove the second condition in order to obtain the existence of an exponential attractor.
Letbounded and. The semigroup,generated by the problem (
2.1
)–(
2.4
) is Hölder continuous on.
The lipschitz continuity in space being a consequence of (4.9), it just remains to prove the continuity in time (actually, a Hölder condition in time for the semigroup , ). One assumes that the initial data belong to . For every and , owing to the above estimates, one gets
where c depends on T. One multiplies (2.2) by to obtain
Integrating (6.29) between and , one deduces from the above estimates that
where c depends on T and . Finally, one has
where c depends on T and , which finishes the proof of the Hölder continuity with respect to t. □
One finally deduces from Proposition 6.1 and Proposition 6.2 the following result [13,14].
The semigrouppossesses an exponential attractor, i.e.,
has a finite fractal dimension,;
is positively invariant,, for every;
attracts exponentially fast the bounded subsets of Φ,wheredenotes the Hausdorff semidistance (for the norm of X) defined as
The exponential attractor being finite-dimensional and , this ensures the boundedness of the fractal dimension of the global attractor.
Compared to the global attractor, an exponential attractor is expected to be more robust under perturbations. Indeed, the rate of attraction of trajectories to the global attractor may be slow and it is very difficult, if not impossible, to estimate this rate of attraction with respect to the physical parameters of the problem in general. As a consequence, global attractors may change drastically under small perturbations. One refers the reader to [21,24] for discussions on this subject.
Footnotes
Acknowledgements
The authors wish to thank Alain Miranville for useful discussions. They also wish to thank an anonymous referees for their careful reading of the paper and useful comments.
References
1.
A.Agmon, Lectures on Elliptic Boundary Value Problems. Mathematical Studies, Van Nostrand, New York, 1965.
2.
S.Agmon, A.Douglis and L.Nirenberg, Estimates near the boundary for solutions of elliptic partial differential equations. I, Commun. Pure Appl. Math.12 (1959), 623–727. doi:10.1002/cpa.3160120405.
3.
G.Brochet, X.Chen and D.Hilhorst, Finite dimensional exponential attractors for the phase-field model, Appl. Anal.49 (1993), 197–212. doi:10.1080/00036819108840173.
4.
G.Caginalp, An analysis of a phase-field model of a free boundary, Arch. Rational Mech. Anal.92 (1986), 205–245. doi:10.1007/BF00254827.
5.
G.Caginalp, Conserved-phase field system: Implications for kinetic undercooling, Phys. Rev. B.38 (1988), 789–791. doi:10.1103/PhysRevB.38.789.
6.
G.Caginalp, The dynamics of conseved phase-field system: Stefan-Like, Hele-Shaw and Cahn-Hilliard models as asymptotic limits, IMA J. Appl. Math.44 (1990), 77–94. doi:10.1093/imamat/44.1.77.
7.
G.Caginalp and E.Esenturk, Anisotropic phase-field equations of arbitrary order, Discrete Contin. Dyn. Systems S.4 (2011), 311–350. doi:10.3934/dcdss.2011.4.311.
8.
L.Cherfils and A.Miranville, Some results on the asymptotic behavior of the Caginalp system with singular potentials, Adv. Math. Sci. Appl.17 (2007), 107–129.
9.
L.Cherfils and A.Miranville, On the Caginalp system with dynamic boundary conditions and syngular potential, Appl. Math.54 (2009), 89–115. doi:10.1007/s10492-009-0008-6.
10.
L.Cherfils, A.Miranville and S.Peng, Higher-order models in phase separation, J. Appl. Anal. Comput.7 (2017), 39–56.
11.
L.Cherfils, A.Miranville, S.Peng and W.Zhang, Higher-order generalized Cahn–Hilliard equation, Electron. J. Qualitative Theory Diff. Eqns.9 (2017), 1–22.
12.
C.I.Christov and P.M.Jordan, Heat conduction paradox involving second-sound propagation in moving media, Phys. Rev. Lett.94 (2005), 154.
13.
A.Eden, C.Foias, B.Nicolaenko and R.Temam, Exponential Attractors for Dissipative Evolution Equations, Research in Applied Mathematics, Vol. 37, John-Wiley, New York, 1994.
14.
M.Efendiev, A.Miranville and S.Zelik, Exponential attractors for a nonlinear reaction-diffusion system in , C.R. Acad. Sci. Paris Série I Math.330 (2000), 713–718. doi:10.1016/S0764-4442(00)00259-7.
15.
A.Miranville, Some mathematical models in phase transition, Discrete Contin. Dyn. Systems Seres S.7 (2014), 271–306.
A.Miranville and A.J.Ntsokongo, On anisotropic Caginalp phase-field type models with singular nonlinear terms, J. Appl. Anal. Comput.8 (2018), 655–674.
18.
A.Miranville and R.Quintanilla, Some generazalisations of the Caginalp phase-field system, Appl. Anal.88(6) (2009), 877–894. doi:10.1080/00036810903042182.
19.
A.Miranville and R.Quintanilla, A conserved phase-field system based on the Maxwell–Cattaneo law, Nonlinear Analysis: Real World Applications14 (2013), 1680–1692. doi:10.1016/j.nonrwa.2012.11.004.
20.
A.Miranville and S.Zelik, Robust exponential attractors for Cahn–Hilliard type equations with singular potentials, Math. Methods Appl. Sci.27 (2004), 545–582. doi:10.1002/mma.464.
21.
A.Miranville and S.Zelik, Attractors for dissipative partial differential equations in bounded and unbounded domains, in: Handbook of Differential Equations, Evolutionary Partial Differential Equations, C.M.Dafermos and M.Pokorny, eds, Vol. 4, Elsevier, Amsterdam, 2008, pp. 103–200.
22.
A.J.Ntsokongo, On higher-order anisotropic Caginalp phase-field systems with polynomial nonlinear terms, J. Appl. Anal. Comput.7 (2017), 992–1012.
23.
J.E.Taylor, Mean curvature and weighted mean curvature, Acta Metall. Mater.40 (1992), 1475–1485. doi:10.1016/0956-7151(92)90091-R.
24.
R.Temam, Infinite-Dimensional Dynamical Systems in Mechanics and Physics, 2nd edn, Applied Mathematical Sciences, Vol. 68, Springer-Verlag, New York, NY, 1997.