Abstract
We consider the sharp interface limit for the scalar-valued and vector-valued Allen–Cahn equation with homogeneous Neumann boundary condition in a bounded smooth domain Ω of arbitrary dimension
Keywords
Introduction
In the following we introduce the Allen–Cahn equation and the vector-valued variant considered in this paper. Moreover, we motivate and review results concerning the corresponding sharp interface limits.
Let us begin with the scalar Allen–Cahn equation. Let

Typical form of the double-well potential,
The Allen–Cahn equation (similar to (AC1)) was originally introduced by Allen and Cahn [11] to describe the evolution of antiphase boundaries in certain polycrystalline materials. Moreover, one can directly verify that equation (AC1)–(AC3) is the
The Allen–Cahn equation is a diffuse interface model: the

Diffuse interface and sharp interface limit.
Next, let us introduce the vector-valued Allen–Cahn equation considered in the paper. Let N, Ω,
In analogy to the scalar case one can compute directly that equation (vAC1)–(vAC3) is the
The potential W. Here we allow two types of potentials W. On the one hand, we consider W with exactly two distinct minima and symmetry with respect to the hyperplane in the middle of these. On the other hand, we consider

Typical triple-well potential W. The image is taken from Kusche [54].
Let W has exactly two global minima W is a symmetric triple well-potential for
An example for a typical triple-well potential that fulfils the conditions in Definition 1.1 can be found in [54, Section 3.4]. See Fig. 3 above. The example stems from Haas [45], see also Garcke, Haas [43]. Compared to the scalar case, we always require symmetry properties for the potential W. The assumption is used e.g. in Section 4.3.1 below, but it might be possible to relax this.
Analogously to the scalar case one can argue with formal arguments (system of fast reaction and slow diffusion; gradient flow to the energy (1.4)) that diffuse interfaces for solutions of (vAC) should develop after short time. Note that the transition between minima of W runs in
Sharp/diffuse interface models and sharp interface limits. In general, sharp interface models and diffuse interface models are important model categories for the description of interfaces and moving boundaries in a large variety of applications. Some prominent examples are the melting of ice, the motion of an oil droplet in water, crystal growth, biological membranes, porous media, tumour evolution, spinodal decomposition of polymers and grain boundaries, see e.g. [11,14,15,19,33,34,42,63,68,70] and the references therein.
It is an important task to connect diffuse interface models and sharp interface models via their sharp interface limits for the following reasons (see also the partly universal introduction in Caginalp, Chen [24] and general comments in Caginalp, Chen, Eck [25]):
Modelling and Analysis: both types of models can usually be derived or motivated with physical principles, phenomenological observations or geometrical arguments etc., but one always incorporates some constitutive assumptions. Often the derivation for the sharp interface models is more transparent and these models appear simpler and more qualitative. On the other hand, diffuse interface models are usually advantageous in more complicated situations (see e.g. [14], p. 141) and solutions typically have better analytical properties. Especially topology changes do not impose any difficulties. By identifying the sharp interface limit one confirms that the assumptions in the derivations are appropriate as well as that the models are compatible with each other and can be used to describe the same situation. Another motivation is the concept of using the diffuse interface model to extend solutions of the corresponding sharp interface model past singularities. Numerics: diffuse interface models are often simpler to solve numerically. By considering the sharp interface limit one justifies that the numerical solution to the diffuse interface model can be used to approximate the solution to the sharp interface model.
Concerning results for sharp interface limits: in general there are formal results and rigorous proofs for convergence. The formal sharp interface limits are typically based on formal asymptotic expansions or numerical experiments. However, see also Anderson, McFadden, Wheeler [14], p. 156ff for a “pillbox argument”, i.e. reasoning with a small test volume. Regarding rigorous sharp interface limits, one can basically group such results into two types:
Local time results that are applicable before singularities appear, i.e. as long as the interface does not develop singularities and stays smooth. Relatively “strong” results are obtained, e.g. norm estimates.
Global time results using some kind of weak notion for the sharp interface system, e.g. viscosity solutions for mean curvature flow, varifold solutions, distributional solutions, etc.
Formal and rigorous sharp interface limit for the Allen–Cahn-equation. In the situation of the scalar Allen–Cahn equation (AC) formal asymptotic analysis by Rubinstein, Sternberg, Keller [72] yields that the limit sharp interface
For the vector-valued Allen–Cahn equation (vAC) formal asymptotic calculations in Bronsard, Reitich [23] yield that for a triple-well potential W in the sharp interface limit
A 90°-contact angle if a transition of two phases meets the boundary.
A 120°-triple junction if the three phases meet at an interior point.
For the sake of completeness, note that the limit
Of course this “Γ” has a different meaning than the Γ in Fig. 2.
For results in the direction of Γ-convergence with respect to
There are many rigorous results on the sharp interface limit for the Allen–Cahn equation ((AC1) on
We start with the local time results. Via a comparison principle and the construction of sub- and supersolutions, Chen [27] proves local in time convergence as long as the interface stays smooth. Moreover, de Mottoni and Schatzman [30] consider the
For global time results one has to use some weak formulation of (MCF). There is the notion of viscosity solutions used by Evans, Soner, Souganidis [37] for
Altogether, there is a large variety of results. At the time of writing, [27] and [5] were the only results of local type that allow boundary contact for the diffuse interfaces. Moreover, at the time of writing there was only the conditional result [56] in the vector-valued case on the convergence of the vector-valued Allen–Cahn equation to multiphase mean curvature flow in a BV-setting. But very recently, the relative entropy approach from [40] has been generalized to the vector-valued case by Fischer, Marveggio [41], even for potentials with more than three wells, but at the moment for a slightly different class of potentials compared to Definition 1.1 above, e.g. not allowing
Note that there is a work by Sáez [74] on the vector-valued problem, but unfortunately there is the same gap in the proof as in [73], cf. [5] for details. It is difficult to generalize the methods in [27] because comparison principles are used. On the other hand, the method by de Mottoni and Schatzman [30] has also proven to be versatile and was applied to many other diffuse interface models as well, see the comments below.
This motivates to extend the result in [5] to more complicated geometrical situations and equations. In this paper we generalize the latter in two directions. On the one hand, we look at the higher-dimensional setting and on the other hand, we also consider the situation of a two-phase transition for the vector-valued Allen–Cahn equation. The following results are obtained:
Convergence of (solutions to) the scalar-valued Allen–Cahn equation (AC) to (MCF) with 90°-contact angle in any dimension Convergence of (solutions to) the vector-valued Allen–Cahn equation (vAC) to (MCF) with 90°-contact angle in any dimension
The results are part of the PhD thesis of the author, cf. Moser [66]. In a separate paper by Abels, Moser [6] the work [5] is extended for the case of a non-linear Robin boundary condition which is designed in such a way that in the sharp interface limit (MCF) with α-contact angle is attained, where α is fixed but close to 90°.
The method of de Mottoni and Schatzman. One assumes that there exists a local smooth solution to the limit sharp interface problem. This can typically be shown for small times. Then
One rigorously constructs an approximate solution to the diffuse interface model using asymptotic expansions based on the evolving surface that is (part of) the solution to the limit problem. To this end one has to solve model problems for the series coefficients. Then one estimates the difference between exact and approximate solutions with a Gronwall-type argument. This typically involves a spectral estimate for a linear operator associated to the diffuse interface equation and the approximate solution.
In the method comparison principles are not needed in a fundamental way compared to many other approaches. However, the application can be difficult for complicated equations or too singular situations. On the other hand, one gets a strong result: norm estimates and the typical profile of the solution are obtained.
Therefore the method was used for many other diffuse interface models as well. These results are based on general spectrum estimates in Chen [28] for Allen–Cahn, Cahn–Hilliard and phase-field-type operators. For the subtle variations in the rigorous asymptotic expansions and the spectral estimates used in applications of the method by de Mottoni and Schatzman see the inceptions to Sections 5 and 6 in [66]. There are results for the Cahn–Hilliard equation by Alikakos, Bates, Chen [9], the phase-field equations by Caginalp, Chen [24], the mass-conserving Allen–Cahn equation by Chen, Hilhorst, Logak [29], the Cahn–Larché system by Abels, Schaubeck [7] and a Stokes/Allen–Cahn system by Abels, Liu [2]. Moreover, there is the paper [1] on the sharp interface limit for a Navier–Stokes/Allen–Cahn system with different viscosities. Moreover, Marquardt [62] (see also [3,4]) studied the sharp interface limit for a Stokes/Cahn–Hilliard system. See also Schaubeck [75] for a result on a convective Cahn–Hilliard equation. Furthermore, there is the result by Fei, Liu [38], where a phase field approximation for the Willmore flow is considered. Finally, [5] is the first result obtained with the method of de Mottoni and Schatzman that allows boundary contact for the diffuse interfaces.
Structure of the paper. In Sections 1.1–1.2 we formulate the main results for the scalar-valued and vector-valued Allen–Cahn equation, (AC) and (vAC) respectively. In Section 2 we fix some notation and introduce function spaces. In Section 3 suitable curvilinear coordinates are constructed. In Section 4 we solve the model problems appearing in the asymptotic expansions in Section 5. The difference estimates and the proofs of the convergence theorems are carried out in Section 6. Therefore suitable spectral estimates are needed, see [67] for those.
Note that Section 4 on the model problems can be read independently since it solely relies on Section 2 on function spaces. Moreover, Section 3 on the existence of curvilinear coordinates is also self-contained, and does only need some notation from the beginning of Section 2. Section 5 on the asymptotic expansions then needs the results of Sections 3–4, but just the definitions from Section 2. The main point is then the short Section 5.1.3 and 5.2.3, respectively. The latter can be understood without precisely reading the remainder of Section 5 and is the only part that is needed from Section 5 for Section 6. Finally, apart from that, Section 6 relies on Section 3 and the spectral estimates from [67].
Comparison to [
5
]. Compared to [5] the geometrical situation is more difficult because the dimension is arbitrary. However, with similar but more technical arguments we also obtain suitable curvilinear coordinates in Section 3. The model problems in the scalar case are the same as in [5], cf. Sections 4.1–4.2, but for the half space problem we use a new setting with exponentially weighted Sobolev spaces which allows for solution operators. With the latter the iteration of the asymptotic expansion in Section 5.1 is more efficient and one can solve the appearing model problems in [6] via the Implicit Function Theorem for angles α close to
In the vector-valued case, new ideas are required for the ODE model problems on the real line, cf. Section 4.3. For the nonlinear vector-valued ODE in Section 4.3.1 an energy minimization argument is used, which was done before by [22] and [54]. In order to solve the linearized vector-valued ODE in Section 4.3.3 we make the assumption that the kernel of the corresponding linearized operator is 1-dimensional, cf. Remark 4.11. This then gives rise to a solution operator for the linearized vector-valued ODE, also in exponentially weighted spaces via perturbation arguments. Moreover, the assumption in Remark 4.11 yields a spectral gap as in the scalar case, cf. Lemma 4.12. Based on this, one can solve the vector-valued model problem on the half space in the analogous way as in the scalar case, see Section 4.4. The asymptotic expansion in Section 5.2 and the difference estimate in Section 6.3 is then analogous to the scalar case. Finally, note that the spectral estimates needed for Section 6.3 also had to be adjusted suitably from the scalar case to the vector-valued case, cf. [67]. It was unclear a priori how to achieve the latter, because in the scalar case there was frequent use of theorems such as the maximum principle for scalar 1D-operators. This is another reason why the assumption on the kernel from Remark 4.11 is used.
Mean curvature flow ( MCF ) with 90°-contact angle, coordinates and notation. “Mean curvature flow” for evolving hypersurfaces means that the normal velocity equals mean curvature, where we define for convenience “mean curvature” as the sum of the principal curvatures.
For the convergence result below we will assume that (MCF) together with a 90°-contact angle condition at
The local well-posedness and existence of a smooth solution for small time starting from suitable initial sharp interfaces is basically well-known. At this point let us give some references in this direction. In Katsoulakis, Kossioris, Reitich [52, Section 2], a parametric approach is used to show local existence and uniqueness of classical solutions for (MCF) in arbitrary dimension and with fixed contact angle. In principle, it is also possible to reduce the evolution to a parabolic PDE by writing it over a reference hypersurface via suitable coordinates. For the typical procedure in the case of a closed interface see Prüss, Simonett [70]. For curvilinear coordinates in the situation of boundary contact see Vogel [79] and Section 3 below. Moreover, note that in Huisken [48] the special case of (MCF) with 90°-contact angle in the graph case for cylindrical domains is considered and global existence and uniqueness of smooth solutions as well as convergence to a constant graph is obtained.
We need some notation in the context of the curvilinear coordinates in order to formulate the main theorems below.
For the details see Section 3.
Domain. Let Sharp Interface. Consider Coordinates. We construct appropriate curvilinear coordinates
Result in scalar case
Finally, we state the new convergence result for (AC) obtained with the method of de Mottoni and Schatzman [30]. Let
Then there are
If
If
If
Interpretation of Theorem
1.4
. One can interpret Layout of the Proof. Required model problems, some ODEs on Well-Posedness of (AC). In Theorem 1.4 the existence of solutions The approximate solution The level sets In principle also estimates of stronger norms (e.g. Hölder norms of arbitrary order) are possible in the situation of Theorem 1.4, but better estimates for the initial values could be required. The basic idea is to interpolate the already controlled norms with stronger norms that can be estimated for exact solutions by some negative ε-orders. Cf. Alikakos, Bates, Chen [9, Theorem 2.3] for a similar idea in the case of the Cahn–Hilliard equation. However, note that this does not improve the approximation of Γ in the sense of 5. Theorem 1.4 and the above comments hold analogously for closed Γ moving by (MCF) and compactly contained in Ω. The proof is basically contained since the constructions are localizable. (Convergence of (AC) to (MCF) with 90°-Contact Angle).
(Convergence of (vAC) to (MCF) with 90°-Contact Angle).
Let
Then there are
If
If
If
The interpretation of Theorem 1.6 is analogous to the one of Theorem 1.4, where convergence of (AC) to (MCF) with 90°-contact angle is obtained, cf. Remark 1.5, 1. Layout of the Proof. The new model problems, some vector-valued ODEs on Well-Posedness of (vAC). In general the analysis of systems is more challenging than that of single equations. However, the derivatives in (vAC1)–(vAC3) are decoupled and hence some methods from the scalar case can also be used for (vAC), e.g. regularity theory. Equation (vAC1)–(vAC3) matches the general setting of Lunardi [60, Section 7.3.1], where a semigroup method and Hölder-spaces are used. Moreover, by reduction to a scalar equation and maximum principle arguments, one can obtain a priori boundedness of classical solutions, see Section 6.1.2 below. Hence global well-posedness for regular, bounded initial data follows. Higher regularity can be obtained with linear theory for scalar equations, cf. Lunardi, Sinestrari, von Wahl [61]. The comments for Theorem 1.4 in Remark 1.5, 4.–7. hold similarly. Here
Let
Consider
Let
Next we introduce the used Sobolev spaces on domains in
Let Let
For the definitions and properties of scalar-valued function spaces, in particular embeddings, interpolation results and trace theorems see Adams, Fournier [8], Alt [12], Leoni [58] and Triebel [77,78]. Many properties can be generalized to vector-valued function spaces over domains, see e.g. Kreuter [53] and the references therein.
Furthermore, we need Sobolev spaces on submanifolds of
Let Then for For
Consider the situation of Definition
2.2
. Let
The proof is technical but straightforward and uses density and localization arguments, cf. [66, Lemma 2.15]. □
Later on
Finally, we define the used spaces with exponential weight.
Let Then we introduce with canonical norms
In a similar way we define Let
The following lemma summarizes all the needed properties for these spaces:
The spaces in Definition
2.4
are Banach spaces. Equivalent norms: Let Density of smooth functions with compact support: For all Embeddings: It holds Traces of weighted functions on Reverse Fundamental Theorem for weighted
The proofs of 1.–5. are straightforward and use the properties of the unweighted spaces. Cf. [66, Lemma 2.22]. By density it is enough to prove the estimate for Let
Let
For convenience. The considerations can be adapted for the case of finitely many connected components.
For the definition of an evolving hypersurface cf. Depner [32, Definition 2.31].
Let
Such
Since continuous bijections of compact topological spaces into Hausdorff spaces are homeo-morphisms, we know that
Additionally, we require
Let the assumptions in Section 3.1 hold. We consider the outer unit conormal
(Curvilinear Coordinates ).
Let the above assumptions in Section
3
up to this point hold. There exist
Let
Let
Let Let For transformation arguments we define
The proof of Theorem 3.2 is similar to [5, Theorem 2.1], where the case
There is an
The assertions are compatible with shrinking η for small
Let
Let
We choose
Now we define X. Let
First of all, X is well-defined due to the cutoff-function. Moreover, X is smooth and
Next we prove that
Consider
Let
Finally, we show the properties of b. Let
Finally, we show relations of
Let
It holds
Set
We only consider The second equivalence estimate is evident since Consequence of 1. and The chain rule yields
Analogous assertions hold for ψ defined on
Unless otherwise stated we use real-valued function spaces in this section.
Some scalar-valued ODE problems on
In this section we recall existence and regularity results needed for ODEs appearing in the inner asymptotic expansion for (AC), cf. Section 5.1.1 below. For the potential
The ODE for the optimal profile
The ODE system for the lowest order (cf. Section 5.1.1 below) is
([75, Lemma 2.6.1]).
Let f be as in (
1.1
). Then (
4.1
) has a unique solution
The idea is to solve the equivalent first order ODE
We call

Typical optimal profile
The following theorem is concerned with the solvability of the equation which is obtained by linearization of (4.1) at
([75, Lemma 2.6.2]).
Let
Let
For our purpose The result follows from the proof of [75, Lemma 2.6.2]. The idea is to reduce to a first order ODE for the derivative of
Let
In Section 4.2.1 we recall existence and uniqueness assertions for weak solutions from [5]. For the proofs see [5, Section 2.4.1]. Moreover, in [5, Section 2.4.2] exponential decay estimates were proven via ordinary differential inequality arguments. Here we proceed differently. In Section 4.2.2 we introduce a functional analytic setting with several types of exponentially weighted Sobolev spaces in order to have solution operators for (4.3)–(4.4). The rough idea is always to multiply the equation with the weights, use the product rule and known isomorphisms.
The framework with exponentially weighted Sobolev spaces helps with additional independent variables and simplifies the induction procedure needed in the asymptotic expansion in Section 5.1 below. Note that in the 2D-case an induction is not necessary and hence was not carried out in [5]. Moreover, the isomorphism property can be used to solve via the Implicit Function Theorem a corresponding model problem appearing in the sharp interface limit for an Allen–Cahn equation with a nonlinear Robin boundary condition designed to approximate (MCF) with a constant α-contact angle for α close to
Weak solutions and regularity
The definition of a weak solution is canonical:
([5, Definition 2.7]).
Let
Regarding weak solutions we have the following theorem:
([5, Theorem 2.9]).
Let
If
Weak solutions are unique.
If
If
If
Here the spectral properties of the linear operator
In Theorem 4.4, 6. weaker conditions on G are enough, cf. Section 4.2.2. The point is included for aesthetic reasons. Altogether the following existence theorem for weak solutions is obtained:
([5, Corollary 2.10]).
Solution operators in exponentially weighted spaces
In the following the superscript “⊥” always means
For
For all
For
Let Via (4.6) we have isomorphic splitting operators from Let The idea is similar as in the proof of Theorem 4.6, 1. In order to solve (4.10), we observe that Now let Let Now we show The remaining assertions and estimates will be shown by differentiating and rearranging the first equation in (
The structure of this section is similar to Section 4.1 which is the analogue in the scalar case. We consider vector-valued ODEs appearing in the inner asymptotic expansion of (vAC) at a two-phase transition (cf. Section 5.2.1 below) and also the linear operator belonging to a linearized ODE. In order to solve the linearized ODE, we will assume that the kernel of the linearization is 1-dimensional. The latter is fulfilled for a typical potential, cf. the example in Remark 4.11 below.
Let
The nonlinear ODE
The nonlinear ODE problem in the lowest order (cf. Section 5.2.1 below) is the following: Find Let m, W,
From now on, we fix a The proof relies on minimizing an energy over an appropriate set (see below). Similar to Bronsard, Gui, Schatzman [22, Section 2], where transitions between two minima in the triple-well case is considered, it should be possible to determine the qualitative behaviour of the set of 1D-minimizers for both types of W in Definition 1.1 precisely. E.g. in the triple-well case the 1D-minimizers are trapped in the smaller sector between In Bronsard, Gui, Schatzman [22, Section 4] the existence of a solution of the two-dimensional problem
Let ([54, Section 2.1], [22, Section 2]).
We look at the operator obtained by linearization of the left hand side of the ODE (4.11) at
Let
That (Assumption
We have to consider the ODE that arises from the linearization of (4.11) at
The additional condition with Let
Let
There is a
Dependence on parameteres. In Theorem 4.13 we obtained linear solution operators in suitable spaces with exponential decay. Therefore, if the right hand side in the first differential equation of (4.14) depends on additional parameters and satisfies such exponential decay estimates, the regularity and decay carries over to the solution. Using Theorem 4.13 one can directly obtain a result for right hand sides Kusche [54, Proposition 1.6 and Corollary 1.2] yield pointwise exponential decay estimates. The latter would be enough for our purpose. But the downside is that the exponent shrinks. Moreover, the proof of Theorem 4.13 is self-contained and simpler. However, in [54] there are also uniform estimates for finite large intervals which are important for the spectral estimates, cf. [67]. Consider Now let We prove this with similar ideas as in Section 4.2.2, i.e. for
This section is analogous to Section 4.2, where the scalar case was done. Let
Weak solutions and regularity
In the vector-valued case weak solutions are defined as follows:
Let
We obtain the analogue of Theorem 4.4:
Let
If
Weak solutions are unique.
If
If
If
One can essentially copy the proof of Theorem 4.4. All multiplications of functions that are now vector-valued have to be interpreted as scalar products and
Again this yields an existence theorem for weak solutions analogously to Corollary 4.5:
Let
The first part is a direct consequence of Theorem 4.16. The second assertion can be shown similar to the proof of [5, Corollary 2.10, 2.] using iteratively scalar regularity theory for every component of the elliptic equation, where
With analogous adjustments as above one obtains solution operators in exponentially weighted Sobolev spaces similar to Section 4.2.2, cf. Theorems 4.6–4.8. We will just need the analogue of Theorem 4.8, hence we only formulate the latter in the vector-valued setting:
Let
(Solution Operators for the Vector-valued Case).
Asymptotic expansions
In this section we carry out the rigorous asymptotic expansions for (AC) and (vAC) in the situations mentioned in the introduction. The expansions are based on the curvilinear coordinates from Section 3 and use the solutions for the model problems in Section 4.
Asymptotic expansion of (AC) in ND
Let
Let
In Section 5.1.1 we construct the inner expansion and in Section 5.1.2 the contact point expansion. Finally, in Section 5.1.3 we show that the construction yields a suitable approximate solution
Inner expansion of (AC) in ND
For the inner expansion we consider the following ansatz: Let
In the following the surface gradient
Let
The differential operator Note the similarity to [5, Lemma 3.1]. The terms on the right hand side are more sophisticated, but are always of the same type as in the 2D-case. Therefore in the expansion they will contribute in the analogous way.
Basically, this follows from the chain and product rules as well as from the properties of
To expand the Allen–Cahn equation
Moreover, we expand functions of
For the higher orders in the expansion the following definition is useful:
We call Let Let For Note that this set only appears in the contact point expansion later.
Now we expand (AC1) for
Inner expansion:
Inner expansion:
Inner expansion:
The norm of the entirety of derivatives in local coordinates up to any fixed order
If additionally f is even, then
Inner expansion:
Then with the notation in Definition 5.3 it holds for all
With these identities one can compute the
In the contact point expansion we proceed similarly as in the 2D-case, cf. [5, Section 3.2], but the computations are more technical. Here we have more contact points, namely all
For the expansion we calculate the action of the differential operators on
Let
This can be shown in a similar manner as in the proof of Lemma 5.1. □
The formulas in Lemma 5.6 without the
Contact point expansion: The bulk equation. We expand the
Moreover, we expand terms arising from Lemma 5.6 in (5.11) that are functions of
Let us introduce some notation for the higher orders in the expansion:
We denote with Let Let The set
In the following we expand the bulk equation (5.11) with the above formulas into ε-series with coefficients in
Bulk Equation:
Bulk Equation:
Then with the notation from Definition 5.8 it holds for all those
With these identities, one can compute that the
Contact point expansion: The Neumann boundary condition. The boundary conditions complementing (5.15) and (5.16) will be obtained from the expansion of the Neumann boundary condition (AC2) for
For
In the following we expand the Neumann boundary condition into ε-series with coefficients in
Neumann Boundary Condition:
Neumann Boundary Condition:
If f is even, then so is
Now we solve (5.5) together with (5.22). We show that the principal part in (5.5) satisfies a suitable ellipticity condition and that (5.22) fulfils a non-tangentiality condition in local coordinates. Based on this one can show maximal regularity results in Hölder spaces (similar to Lunardi, Sinestrari, von Wahl [61]) and Sobolev spaces (similar to Prüss, Simonett [70, Chapters 6.1–6.4] and Denk, Hieber, Prüss [31]) with typical localization procedures. This always involves compatibility conditions for the initial value. In our case these can be avoided via extension arguments similar to [5, Section 3.2.2]. All these arguments involve many technical computations, but are in principle well-known. Therefore we refrain from going into details.
The ellipticity condition. Let
The non-tangentiality condition. Let
Finally, we obtain a smooth solution
Neumann Boundary Condition:
With the notation as in Definition 5.8 it holds for
As in the last subparagraph, the equations (5.16), (5.23) are equivalent to
Because of the computations in the last subparagraph one can solve (5.10) from Section 5.1.1 together with (5.26) and get a smooth solution
Finally, by induction the jth order is determined for all
Let
The function
The estimate also holds without the The third and the last equation are evident from the construction. Moreover, the rigorous Taylor expansions (5.2)–(5.3), (5.12)–(5.14) and (5.17) together with the remarks for the remainders and Sections 5.1.1–5.1.2 yield
It is left to prove the remaining assertion for
Let
Let
In Section 5.2.1 we construct the inner expansion and in Section 5.2.2 the contact point expansion. Finally, in Section 5.2.3 the result on the approximation error of
Inner expansion of (vAC) in ND
For the inner expansion we consider the following ansatz: Let
In the following we use the same conventions as in Lemma 5.1. Moreover, for a sufficiently smooth
Let
This follows from Lemma 5.1 applied to every component. □
For the expansion of (vAC1) for
Moreover, we expand functions of
For the higher orders in the expansion we use analogous definitions as in the scalar case:
We call Let For
Now we expand (vAC1) for
Inner expansion:
Inner expansion:
Inner expansion:
If
Inner expansion:
Contact point expansion of (vAC) in ND
This is analogous to the scalar case, cf. Section 5.1.2. We make the ansatz
For the expansions we calculate the action of the differential operators on
Let
This can be shown by applying Lemma 5.6 to every component. □
Contact point expansion: The bulk equation. We expand the
Functions of
As in the scalar case we use some notation for the higher orders in the expansion:
We call Let The set
In the following we expand (5.37) into ε-series with coefficients in
Bulk Equation:
Bulk Equation:
Contact point expansion: The Neumann boundary condition. As in the scalar case, the boundary conditions complementing (5.38)–(5.39) will be obtained from the expansion of the Neumann boundary condition (vAC2) for
For
In the following we expand the Neumann boundary condition into ε-series with coefficients in
Neumann Boundary Condition:
Neumann Boundary Condition:
Hence we get
Neumann Boundary Condition:
Analogously as in the scalar case, the
As in the last subparagraph, the equations (5.39), (5.45) are equivalent to
Because of the remarks and computations in Section 5.1.2 we can solve (5.36) from Section 5.2.1 together with (5.48) and obtain a smooth solution
Finally, by induction the jth order is constructed for all
The approximate solution for (vAC) in ND
Let
The function
The analogous assertions as in Remark 5.11 hold.
In this section we estimate the difference of the exact and approximate solutions with a Gronwall-type argument in both cases. This is the second step in the method by de Mottoni and Schatzman [30]. Here the major ingredients are the spectral estimates from [67]. Moreover, we have to control certain nonlinear terms stemming from differences of potential terms. We will estimate these with interpolation inequalities. Therefore as preparation we prove uniform a priori bounds for exact classical solutions of (AC) and (vAC) in Sections 6.1.1 and 6.1.2, respectively. Moreover we recall some Gagliardo–Nirenberg estimates in Section 6.1.3. For Allen–Cahn type models such uniform bounds in ε for the exact solutions is typical, see de Mottoni, Schatzman [30, Section 6] for the standard Allen–Cahn equation as well as [2, Remark 1.2 and Section 5.2] for the Allen–Cahn equation coupled with the Stokes system. In Section 6.2 we prove the difference estimate in the scalar case and Theorem 1.4. This is similar to [5, Sections 5–6], but more technical. Finally, in Section 6.3 we show the difference estimate in the vector-valued case and Theorem 1.6. The latter works analogously to the scalar case.
Preliminaries
Uniform a priori bound for classical solutions of (AC)
Let N, Ω, The proof is done via contradiction with ideas from the maximum principle. The argument works similarly in the vector-valued case which is slightly more complicated, see Lemma 6.2 below. Therefore we leave out the proof. See [66, Lemma 7.1] for the details. □
Let N, Ω, Let We use a contradiction argument and ideas from the proof of the weak maximum principle for parabolic equations, cf. Renardy, Rogers [71, Theorem 4.25]. We make the assumption that Finally, let
Let us recall some Gagliardo–Nirenberg inequalities.
Let
See Leoni [58, Theorem 12.83]. □ With suitable extension operators the estimate in Lemma 6.3 carries over to domains Ω with uniform Lipschitz boundary if the (Gagliardo–Nirenberg Inequality).
We prove in Section 6.2.1 a rather abstract estimate for the difference of exact solutions and suitable approximate solutions for the Allen–Cahn equation (AC1)–(AC3) in ND. Then in Section 6.2.2 we show the Theorem 1.4 about convergence by verifying the requirements for the difference estimate applied to the approximate solution from Section 5.1.3.
Difference estimate
(Difference Estimate for (AC)).
Let
Moreover, let
For some
The parameter M corresponds to the order of the approximate solution in Section 5.1. The parameter β was introduced in order to obtain a result valid for all times Note that weaker requirements in the theorem also work, e.g. when one does not have the two additional terms on the right hand side of the spectral estimate or only an estimate with the full That the parameter
The continuity of the objects on the left hand side in (6.5) yields that
We carry out a general computation first and return back to the different cases later. The difference of the left hand sides in (AC1) for
In the following we estimate the
Under the assumptions in Theorem
6.5
it holds
Note that one could apply the same standard Gagliardo–Nirenberg Inequality used for
It remains to estimate
Now we consider the cases in the theorem.
If
Let
Finally, let
The proof of Theorem 6.5 is completed. □ First, let us estimate (6.10) for Next we estimate (6.10) for Finally, we collect the above estimates. To reduce the number of
Let
Note that
First we estimate the boundary term in (6.3). Lemma 5.10 yields
Finally, we estimate the
Difference estimate and proof of the convergence theorem for (vAC) in ND
We show in Section 6.3.1 the difference estimate for exact and suitable approximate solutions for the vector-valued Allen–Cahn equation (vAC1)–(vAC3). Then in Section 6.3.2 we prove the Theorem 1.6 about convergence by checking the requirements for the difference estimate applied to the approximate solution from Section 5.2.3. All computations are analogous to the scalar case in the last Section 6.2.
Difference estimate
(Difference Estimate for (vAC)).
Let
Moreover, let
For some
The parameter M corresponds to the order of the approximate solution constructed in Section 5.2. The comments for the scalar case in Remark 6.6, 2.–4., on the role of the parameters β,
The continuity of the objects on the left hand side in (6.20) yields that
We do a general computation first and consider the specific cases later. The difference of the left hand sides in (vAC1) for
Now we estimate the
In order to control
Now we consider the different cases in the theorem.
If
Let
Finally, let
The proof of Theorem 6.9 is completed. □
Let
Theorem 1.6 follows directly from Theorem 6.9 if we prove the conditions 1.–4. in Theorem 6.9. The requirement 1. (“uniform boundedness”) is satisfied because of Lemma 5.16 for
Footnotes
Acknowledgements
The author gratefully acknowledges support through DFG, GRK 1692 “Curvature, Cycles and Cohomology” during parts of the work.
