in this article we study the allen - cahn / cahn - hilliard model ( ac - ch model for short ) , which combines and extends two famous diffuse interface models : the allen - cahn equation and the cahn - hilliard equation .
both have been applied successfully to model segregation , precipitation and phase change phenomena in alloys and liquid mixtures in materials science , geology , physics , and biology , among others .
the ac - ch system was first introduced in @xcite , and its mathematical properties have been studied extensively , see , e.g. , @xcite and references therein .
we study here an extension of this model with a particular ansatz for the elastic energy .
the allen - cahn equation was first introduced without elasticity in @xcite ; and the cahn - hilliard equation was first introduced without elasticity in @xcite .
the allen - cahn equation is a second order partial differential equation of ginzburg - landau type for an _ unconserved _ order - parameter @xmath3 and thus can be used to model segregation and precipitation in solids , or other more general situations where a reordering of the crystal lattice occurs .
conversely , the cahn - hilliard equation is a fourth - order partial differential equation for a _ conserved _ order parameter @xmath4 striving to model phase change phenomena where one physical quantity like the volume of the phases is preserved . throughout the paper we set @xmath5 as this plays a special role .
the allen - cahn system with linear elasticity was studied before in @xcite , the cahn - hilliard system with linear elasticity in @xcite , @xcite and @xcite .
an extension of the cahn - hilliard system with geometrically linear elasticity valid for single crystals was recently found in @xcite for @xmath6 .    except for the latter work ,
elastic effects due to small scale microstructures within the phases have been neglected . here
we treat the combined ac - ch model with elasticity in @xmath7 dimensions .
we state the new extended model in  in section  [ secacch ] .
this model provides a basis for the generalization of further existing isothermal diffuse interface models .
our definition of the extended model does not require any special assumptions ( except regarding the regularity ) on the stored energy functional @xmath8 .
however , for the proof of existence of weak solutions , we require the assumption  ( @xmath9 ) phrased in section  [ secexist ] below .
we show that assumption  ( @xmath9 ) is satisfied in the following cases :    * materials which follow the linear theory of elasticity developed by eshelby @xcite in the context of elastic inclusions and inhomogeneities ( see a remark after assumption  ( @xmath9 ) ) ; * materials in @xmath6 which are well - described by a geometrically linear theory of elasticity ( theorem  [ thmcor1 ] , for @xmath10 , cf .
also @xcite ) ; * materials in @xmath1 which are well - described by a geometrically linear theory of elasticity , where the deformations are assumed to be scalar - valued functions .
this corresponds to the situation of anti - plane shear , see below for details and cf .
theorem  [ thmcor1 ] .
the geometrically linear theory of elasticity allows for fine microstructures within the phases of the elastic materials taken into account .
this becomes apparent in the explicit formulas for the relaxed energy functionals , which are , together with formulas for their derivatives , the basis of our proof of theorem  [ thmcor1 ] . in section  [ secexpl ]
we recall an explicit formula in @xmath10 derived by chenchiah and bhattacharya @xcite .
moreover we prove explicit formulas for @xmath2 , and for @xmath1 in the setting of scalar - valued deformations .
( in @xcite partial results for @xmath1 in the vectorial , i.e. the non - scalar , setting are shown . )    the coupling to elasticity changes significantly the morphology of the precipitates and the coarsening patterns , see , e.g. , the classification in @xcite .
for an intuitive picture of the microstructures taken into account we refer the reader to figures  [ fig : intuitivemacro ] and [ fig : intuitivemicro ] .
macro.pstex_t    micro2.pstex_t micro1.pstex_t    figure  [ fig : intuitivemacro ] displays exemplary macroscopic phases of a segregated material with phases @xmath11 being almost one ( dark gray ) and almost zero ( light gray ) . in figure
[ fig : intuitivemicro ] we show the corresponding microscopic length scale featuring microstructures ; the left is a blow - up of a small region within the phase with @xmath11 being almost one and the right is a blow - up of a very small region within the phase with @xmath11 being almost zero . on this microscopic scale where we assume that it is sufficient to treat the elastic energies of single crystals , fine microstructures in the form of laminates occur , see , e.g. , @xcite .
the methods developed here apply generally to any established phase change and segregation model provided the temperature is conserved .
( for non - isothermal settings , the validity of the second law of thermodynamics requires additional corrections which are not studied here . ) for further discussions of our model and the analytical results we refer to the conclusions in section  [ secconclude ] .
throughout this paper , let @xmath12 for @xmath13 be a bounded domain with lipschitz boundary which serves as an ( unstressed ) reference configuration .
for a stop time @xmath14 , let @xmath15 denote the space - time cylinder . to the allen - cahn / cahn - hilliard system , first derived in @xcite
, we add elasticity , possibly respecting the lamination microstructure , by introducing the system @xmath16 where the function @xmath17 is a _ conserved _ order parameter , typically a concentration , @xmath18 is an _
unconserved _ order - parameter , specifying the reordering of the underlying lattice , @xmath19 denotes the positive semi - definite mobility tensor , @xmath20 is a small constant determining the interfacial thickness .
the choice of @xmath8 determines whether lamination microstructure occurs in the system , see .    by @xmath21
we describe the displacement field , such that a material point @xmath22 in the undeformed body @xmath23 is at @xmath24 after the deformation . then the (
linearized ) strain tensor is defined by @xmath25 where @xmath26 denotes the transpose of a matrix @xmath27 . as usual
, @xmath28 stands for the inner product in @xmath29 , that is @xmath30 , and for @xmath31 we denote the inner product in @xmath32 by @xmath33 moreover , @xmath34 for @xmath35 is the frobenius norm .
we denote the symmetric matrices in @xmath32 by @xmath36 .
the functional @xmath37 represents the stored elastic energy density .
we will choose it either according to the linear ansatz by eshelby , @xmath38 , or as a single crystal composite lamination energy , @xmath39 .
the definition of @xmath40 is given in the coming section [ ssectelast ] .
the linear theory by eshelby @xcite developed in the context of elastic inclusions and inhomogeneities , can be summarized in the following ansatz for the elastic energy @xmath41 for all @xmath42 , @xmath43 , and @xmath44 with a constant @xmath45 .    by @xmath46
we denote the symmetric , positive definite and concentration - dependent elasticity tensor of the system that maps symmetric tensors in @xmath32 to themselves .    the system ( [ ach1])([ach3 ] )
is completed with the definition of the free energy @xmath47 see @xcite , where @xmath48 is the free energy density assumed to be @xmath49 for scalars @xmath50 .
the term @xmath51 $ ] in ( [ psidef ] ) defines the entropic part of the free energy , given in the canonical bernoulli form for perfect mixing , and @xmath52 is the constant temperature .    the functional @xmath53 in ( [ fdef ] )
represents energy effects due to applied forces . in the absence of body forces ,
the work necessary to transform the undeformed body @xmath23 into a state with displacement @xmath54 is then @xmath55 where we use that the applied stress @xmath56 is constant and symmetric .
consequently , @xmath57 is the energy density of the applied outer forces .
the system ( [ ach1])([ach3 ] ) has to be solved in @xmath58 subject to the initial conditions @xmath59 for given functions @xmath60 subject to the neumann boundary conditions for @xmath4 , the no - flux boundary conditions , and the equilibrium condition for applied forces @xmath61 here , @xmath62 defines the stress .    in ( [ bc ] ) , @xmath63 is the unit outer normal to @xmath64 . for simplicity ,
body forces are neglected and it is assumed that the boundary tractions are dead loads given by a constant symmetric tensor @xmath56 . by @xmath65
we denote the mass flux , given by @xmath66 with @xmath67 the chemical potential .
the valid parameter range of @xmath4 and @xmath3 is , see theorem  [ theo1 ] , @xmath68 the inequalities are strict unless @xmath69 or @xmath70 .    for @xmath71 in @xmath23
, we obtain the pathological solution @xmath72 in @xmath58 . for @xmath73 , @xmath74 in @xmath23
we obtain a cahn - hilliard equation in @xmath3 and a pathological equation for @xmath4 .
the system ( [ ach1])([ach3 ] ) includes as special case the _
elastic cahn - hilliard system _ ( setting @xmath75 , @xcite ) @xmath76 with @xmath77 and the boundary and initial conditions correspondingly to above .
moreover the system ( [ ach1])([ach3 ] ) includes as special case the _ elastic allen - cahn equations _ ( setting @xmath78 , @xcite ) which for rescaled @xmath3 with @xmath79 read @xmath80 with @xmath81 and @xmath82 as well as the boundary and initial conditions are as above .
the system ( [ ach1])([ach3 ] ) is exemplary for an isothermal model that exhibits simultaneous ordering and phase transitions .
equation  ( [ ach1 ] ) is a diffusion law for @xmath4 governed by the flux @xmath65 and states the conservation of mass in @xmath23 .
equation  ( [ ach2 ] ) is a simple gradient flow in the descent direction @xmath83 . equation  ( [ ach3 ] ) is a consequence of newton s second law under the additional assumption that the acceleration @xmath84 originally appearing on the left hand side can be neglected ( this can be proved formally by a scaling argument and formally matched asymptotics ) . the vector equation ( [ ach3 ] ) serves to determine the unknown displacement @xmath54 .
[ rem1 ] the equations  ( [ ach1])([ach3 ] ) can be generalized to vector - valued mappings @xmath4 , @xmath3 .
this allows to study situations with more than two phases present . to fix ideas and for the sake of a clear presentation
, we restrict ourselves throughout this paper to scalars @xmath4 and @xmath3 .
[ rem2 ] the equations  ( [ ach1])([ach3 ] ) with boundary conditions ( [ bc ] ) and and comply with the second law of thermodynamics , which in case of isothermal conditions reads for a closed system @xmath85 this inequality can be verified by direct inspection similar to the calculations in @xcite .    in the proof of theorem  [ theo2 ]
we apply the following explicit formulation of ( [ ach1])([ach3 ] ) with constant mobility @xmath86    @xmath87\!\!,\\ \label{ach2e } { \partial_t}b & = & \lambda\triangle b+\frac{{\vartheta}}{2}\left[g'(a - b ) -g'(a+b)\right]+2\kappa_2b-\frac{\partial{{\cal w}}}{\partial d}(a+b,{\varepsilon}({\mathbf{u}})),\\ \label{ach3e } \mathbf{0 } & = & { \mathrm{div}}\left({\partial_\varepsilon}{{\cal w}}(a+b,{\varepsilon}({\mathbf{u}}))\right).\end{aligned}\ ] ]    the subsequent section is devoted to the geometrically linear elastic energy for single crystals .
this is a prerequisite to the discussion of existence and uniqueness results for the extensions of the allen - cahn / cahn - hilliard ( ac - ch ) model studied in section  [ secexist ] .
our main objective in this subsection is to study a geometrically linear theory of elasticity in the context of isothermal phase transitions . for systematic reasons ,
we first recall the linear ansatz dating back to eshelby , @xcite . as a byproduct of the existence theory proved in section  [ secexist ]
we obtain a new existence result for the ac - ch equations with linear elasticity .
subsequently , we introduce the geometrically linear elasticity theory that takes the laminates of the material into account .    in the following
we assume that two phases are present in the considered material which may form microstructures as displayed , e.g. , in figures  [ fig1 ] and [ fig2 ] .
we refer to the energy @xmath88 , @xmath89 of each of the phases as _ microscopic energy _ , cf .  , and to the energy @xmath90 in , which reflects the effective behavior of the system with microstructures , as the _ macroscopic energy_.    to determine the energy @xmath90 in the geometrically linear theory we need to solve a local minimization problem , eqn .
( [ wmin ] ) below , which we shall outline now .
we assume that the volumes occupied by each of the two phases in @xmath23 are measurable sets .
in particular , if @xmath91 , @xmath92 characterize the two phases on the microscale , we have @xmath93 and @xmath94 a.e . in @xmath23 .
the symbol @xmath95 denotes the space of functions of bounded variation , see , e.g. , @xcite . by @xmath96
we denote the average of a function @xmath97 in @xmath23 , where @xmath98 is the @xmath99-dimensional lebesgue measure of a set @xmath100 .
let @xmath101 , @xmath89 , be the stress - free strain ( or eigenstrain ) of the @xmath102-th phase relative to the chosen reference configuration and @xmath103 be its positive definite elasticity tensor .
then the elastic energy density of phase @xmath102 subject to a strain @xmath104 is given by @xmath105 for @xmath106 .
in ( [ widef ] ) we assumed for simplicity that @xmath103 , @xmath107 and @xmath108 are constants , independent of @xmath109 and the order parameter @xmath11 .    under the assumption that the elastic energy adapts infinitely fast and that the surface energy between laminates of the microstructure can be neglected
, the effective elastic energy is , @xcite , @xmath110,\ ] ] where we used @xmath111 the definition ( [ wmin ] ) requires further clarification .
firstly , @xmath112 , and instead of integrating over @xmath23 , one may integrate over @xmath113 , the open ball of radius @xmath114 around @xmath115 , where the mean @xmath116 is now taken over @xmath113 . by homotopy arguments or by results in @xcite ,
any @xmath117 yields the same value of @xmath118 , as long as @xmath119 . taking
the union of such balls then leads to ( [ wmin ] ) .
secondly , the infimum over @xmath120 is the result of homogenization subject to the constraint that the volume fraction of the selected phase is preset by @xmath121 , see ( * ? ? ?
* chapter  10 ) .
this infimum is taken over functions @xmath122 as explained above ensuring @xmath123 .
this is why ( [ wmin ] ) is only meaningful for @xmath124 $ ] .
the second infimum is taken over functions @xmath125 where the condition @xmath126 has to be read as @xmath127 for a.e .
this originates from the requirement that the functional @xmath40 thus defined must be quasi - convex , see @xcite , and is the result of relaxation theory , @xcite , @xcite , as follows . if for prescribed @xmath5 the microscopic elastic energy density is denoted by @xmath129 , then @xmath130 is the elastic energy density of the material with macroscopic strain @xmath131 after microstructure has formed . for @xmath10 ,
explicit analytic formulas for @xmath40 are known , see @xcite as well as formulas of the partial derivatives of @xmath40 .
these will be recalled in section  [ secexpl ] .
the existence of solutions to the allen - cahn / cahn - hilliard equation without elasticity was studied in @xcite with the help of a semigroup calculus . existence and uniqueness of weak solutions to the cahn - hilliard equation with linear elasticity
is proved in @xcite , with geometrically linear elasticity in @xcite .
existence and uniqueness of weak solutions to the allen - cahn equation with linear elasticity is shown in @xcite .
subsequently we provide existence and uniqueness results for ( [ ach1])([ach3 ] ) , where @xmath8 is an elastic energy density satisfying the following assumption ( @xmath9 ) .
( @xmath9 ) the elastic energy density @xmath132 satisfies the conditions    ( @xmath91 ) @xmath133 is strongly monotone uniformly in @xmath11 , i.e. , there exists a constant + @xmath134 such that for all @xmath135 and all @xmath136 @xmath137 ( @xmath92 ) there exists a constant @xmath138 such that for all @xmath136 and all @xmath42 @xmath139 , @xmath140 , @xmath141 and @xmath142 , but are independent of @xmath11 and @xmath131 .
condition  ( @xmath91 ) states that @xmath8 is convex in @xmath131 .
the problem becomes non - convex through the dependence on @xmath11 .
one prominent example satisfying assumption ( @xmath9 ) is the elastic energy @xmath143 in . in section  [ secexam ]
we prove that also the relaxed energy functional @xmath40 defined in ( [ wmin ] ) satisfies the assumption ( @xmath9 ) .
we require the condition @xmath144 where @xmath145 is the solution of ( [ ach3 ] ) for @xmath146 , @xmath147 .
[ theo1 ] let the mobility tensor @xmath148 be positive definite and continuous for all @xmath4 , @xmath3 satisfying ( [ abrange ] ) , @xmath8 fulfill ( @xmath9 ) , @xmath149 be given by ( [ psidef ] ) and the initial data @xmath150 satisfy ( [ abrange ] ) and ( [ w0bound ] ) .
then there exists a weak solution @xmath151 to ( [ ach1])([ach3 ] ) that satisfies    \(i ) @xmath152;\,l^2(\omega)\right)$ ] , + ( ii ) @xmath153 , @xmath154 , + ( iii ) @xmath155 , + ( iv ) the feasible parameter range of @xmath156 is given by ( [ abrange ] ) .
: the statements of the theorem can be proved with the methods developed in @xcite for an allen - cahn system with linear elasticity and in @xcite for a cahn - hilliard system with linear elasticity .
we sketch the main steps .
first we introduce the operator @xmath157 associated to @xmath158 as a mapping from @xmath159 to its dual by @xmath160 from the poincar inequality and the lax - milgram theorem ( which can be applied since @xmath148 is assumed to be positive definite ) we know that @xmath157 is invertible and we denote its inverse by @xmath161 , the green s function .
we have @xmath162 for @xmath163 , we define the inner product @xmath164 with the corresponding norm @xmath165    for a small discrete step size @xmath166 , chosen such that @xmath167 , for time steps @xmath168 with @xmath169 , and given values @xmath170 , @xmath171 , we introduce the discrete free energy functional @xmath172 where ( in case of @xmath173 ) it holds @xmath174 , @xmath175 , the initial values of @xmath4 and @xmath3 . by the direct method in the calculus of variations and assumption  ( @xmath9 ) , it is possible to show that for @xmath176 sufficiently small
, @xmath177 possesses a minimizer @xmath178 .
this minimizer solves the fully implicit time discretization of ( [ ach1])([ach3 ] ) .
next the discrete solution is extended affine linearly to @xmath179 by setting for @xmath180 with suitable @xmath181 $ ] @xmath182 the validity of the second law of thermodynamics ( cf.remark  [ rem2 ] ) together with ( [ w0bound ] ) implies that @xmath183 is non - increasing in time . in combination with a higher integrability condition on @xmath54 , @xcite , this allows to derive uniform estimates for @xmath179 .
compactness arguments then allow to pass to the limit @xmath184 and the limit solves ( [ ach1])([ach3 ] ) . in general ,
the uniqueness of solutions to  is open .
however , we prove it in a special case for the linear elastic energy density @xmath38 .
[ theo2]@xmath185let @xmath186 be given by ( [ wlin ] ) , the material be homogeneous , i.e. the elasticity tensor @xmath187 be independent of @xmath11 , and let @xmath86 . then the solution @xmath151 of theorem  [ theo1 ] is unique in the spaces stated there .    : the proof is very similar to the proof of uniqueness for the cahn - hilliard system @xcite but we repeat it here because we later need to modify it .
fix @xmath188 .
let @xmath189 , @xmath190 be two pairs of solutions to ( [ ach1e])([ach3e ] ) and ( [ wlin ] ) .
the differences @xmath191 , @xmath192 , @xmath193 with corresponding difference of the chemical potentials @xmath194 solve the weak equations @xmath195{\,\mathrm{d}x}{\,\mathrm{d}t}=0,\\ & { \int\limits_{\omega_{\cal t}}}\left[{\partial_t}b\eta+\lambda\nabla b\cdot\nabla\eta-{\underline{\varepsilon}}:c\left({\varepsilon}({\mathbf{u } } ) -(a+b){\underline{\varepsilon}}\right)\eta\right]{\,\mathrm{d}x}{\,\mathrm{d}t}{\nonumber}\\ \label{test2 } & \ ; = \!{\int\limits_{\omega_{\cal t}}}\left[\frac{{\vartheta}}{2}\left(g'(a^2\!+\!b^2)-g'(a^1\!+\!b^1 ) + g'(a^1\!-\!b^1)-g'(a^2\!-\!b^2)\right)\eta+2\kappa_2b\eta\right]{\,\mathrm{d}x}{\,\mathrm{d}t},\\ \label{test3 } & { \int\limits_{\omega_{t_0}}}c\left({\varepsilon}({\mathbf{u}})-{\underline{\varepsilon}}(a+b)\right):{\varepsilon}({\mathbf{u}}){\,\mathrm{d}x}{\,\mathrm{d}t}=0\end{aligned}\ ] ] for every @xmath196 with @xmath197 , @xmath198 , where in order to get ( [ test3 ] ) we plugged in @xmath199 as a test function and integrated by parts . as a test function in ( [ test1 ] )
we pick @xmath200 this shows @xmath201 the difference of the chemical potentials fulfils , with the help of ( [ fdef ] ) , @xmath202{\,\mathrm{d}x}{\,\mathrm{d}t}.\end{aligned}\ ] ] we pick @xmath203 . with ( [ step1 ] )
we obtain @xmath204|a|{\,\mathrm{d}x}{\,\mathrm{d}t}.\end{aligned}\ ] ] in ( [ test2 ] ) we choose @xmath205 as a test function and add the resulting equation to ( [ step2 ] ) and use ( [ test3 ] ) .
we end up with @xmath206{\,\mathrm{d}x}{\,\mathrm{d}t}\\ & & \le{\int\limits_{\omega_{t_0}}}2(\kappa_1|a|^2+\kappa_2|b|^2){\,\mathrm{d}x}{\,\mathrm{d}t}\\ & & \quad+\!\!{\int\limits_{\omega_{t_0}}}\!\frac{{\vartheta}}{2}\big[\big|g'(a^1\!+\!b^1 ) -g'(a^2\!+\!b^2)\big|+\big|g'(a^1\!-\!b^1)-g'(a^2\!-\!b^2)\big|\big ] ( |a|\!+\!|b|){\,\mathrm{d}x}{\,\mathrm{d}t}.\end{aligned}\ ] ] from theorem  [ theo1 ] we know that the terms @xmath207 , @xmath89 are finite , and @xmath208 is lipschitz continuous . applying first young s inequality , then gronwall s inequality , as @xmath188 was arbitrary , we find @xmath209 in @xmath58 .
this finally yields @xmath210 with korn s inequality this proves @xmath211 in @xmath58 .
in many situations like the numerical implementation of the extended models , the above definition ( [ wmin ] ) of @xmath40 is not practical since it is indirect and based on a local minimization . for these applications and for direct later use ,
we collect here some explicit formulas of the relaxed energy @xmath40 for @xmath6 and for the scalar case in @xmath1 .
as shown in @xcite , it holds @xmath212 where @xmath213 , @xmath214 and @xmath215 are defined below .
first we need to fix further notations following @xcite .
let @xmath216 be given by @xmath217 where @xmath218 is the reciprocal of the largest eigenvalue of @xmath219 , @xmath103 is the elastic modulus of laminate @xmath102 , and the operator @xmath220 is given by @xmath221 in @xcite a recipe is given for the practical computation of @xmath222 . here , we only remark that if the space groups of the two existing laminates are cubic , it holds @xmath223 the first subscript of @xmath187 denotes here the phase , the other two indices are the coefficients of the reduced elasticity tensor in voigt notation , @xcite .    as shown in @xcite , the scalar @xmath224 $ ] determines the amount of translation of the laminates defined by @xmath225 in this definition , @xmath226 is the unique solution of @xmath227 with @xmath228 defined by @xmath229,\\ \triangle{\varepsilon}^ * & = & \triangle{\varepsilon}^*(\beta , d,{\varepsilon } ) \,\,:=\,\,{\varepsilon}_2^*(\beta , d,{\varepsilon})-{\varepsilon}_1^*(\beta , d,{\varepsilon}),{\nonumber}\end{aligned}\ ] ] and the yet undefined functions are specified below .
the four regimes have the following crystallographic interpretation , which follows from the construction of the optimal microstructure in the calculation of @xmath40 .
* regime  0 * : the material is homogeneous and the energy does not depend on the microstructure .
this occurs when @xmath230 .
* regime  i * : there exist two optimal rank - i laminates .
this is characterized by @xmath231<0.\ ] ]    * regime  ii * : the unique optimal microstructure is a rank - i laminate .
this regime occurs when the function @xmath232\ni\beta\mapsto\det\left[\left(d_2\alpha_1+d_1\alpha_2 -\beta t\right)^{-1}\left(\alpha_2({\varepsilon}_2^t-{\varepsilon})- \alpha_1({\varepsilon}_1^t-{\varepsilon})\right)\right]\ ] ] has a unique root ( which we denote by @xmath233 ) .
* regime  iii * : there exist two optimal rank - ii laminates .
this regime is present if the operator @xmath234 is invertible and @xmath235>0.\ ] ]    for illustration , we visualize prototypes of rank - i and rank - ii laminates in figures  [ fig1 ] and [ fig2 ] .
( 13.,3.2 ) yy ( + 4.6,0.0)(9.57,1.4)(1,0)1    ( 13.,3.2 ) yy ( + 3.45,-6.1 )    to complete the definition ( [ phidef ] ) , set @xmath236 hence @xmath237^{-1}e({\varepsilon})$ ] .
explicit computations of @xmath238 and @xmath239 are lengthy .
we recall the following results for the partial derivatives of @xmath40 in @xmath10 which are proved in @xcite . set @xmath240    [ lem1 ] let @xmath10 and let @xmath103 and @xmath241 commute .
then @xmath242 alternatively , in regime  iii , @xmath243 moreover it holds @xmath244 @xmath245 \frac{\partial\beta^*}{\partial{\varepsilon } } & = \begin{cases } \frac{1}{\left ( ( d_2\alpha_1+d_1\alpha_2-\beta t)^{-1}(t\triangle{\varepsilon}^*)\right ) : t\triangle{\varepsilon}^*}(\alpha_1-\alpha_2)\alpha^{-1}t\triangle{\varepsilon}^ * & \text{in regime ii},\\ 0 & \text{otherwise}. \end{cases}\label{eq : beps}\end{aligned}\ ] ]      the scalar setting in three dimensions is characterized by the ansatz @xmath246 for the deformation , where @xmath247 is a scalar function ( hence the name scalar theory , see @xcite , @xcite ) .
physically , ( [ uscal ] ) corresponds to anti - plane shear in the @xmath248-plane . since @xmath249 @xmath250 determines the strain .
this justifies to work with vectors @xmath251 , not with matrices @xmath252 .
we replace ( [ wmin ] ) by @xmath253,\ ] ] where the domain of integration @xmath254 is now two - dimensional .
therefore , @xmath255 , and it is well - understood from the context that @xmath256 denotes the two - dimensional average . in addition , in ( [ wmins ] ) we adapted the common notation for the scalar case and wrote @xmath257 instead of @xmath131 for the strain .
finally , we set @xmath258 with @xmath259 the transformation strains in the scalar setting and @xmath103 positive definite matrices .    in @xcite , partial results for @xmath40 in the non - scalar ,
i.e. vectorial setting are available for @xmath1 .
yet , in general , the computation of @xmath238 , @xmath239 for @xmath1 as required for assumption ( @xmath9 ) remains currently open due to its complexity .
this is the reason why we restrict ourselves to a special three - dimensional case .
our main result in this subsection is the following theorem that provides us with an explicit formula for @xmath40 in the scalar case .
[ theo5 ] the functional @xmath40 , given by ( [ wmins ] ) , satisfies the explicit representation formula @xmath260 with @xmath261,\\ \label{etasdef2 } \nabla\eta_2^*(d , f ) & : = ( d_2\alpha_1+d_1\alpha_2)^{-1}\big[\alpha_1f+d_1 ( \alpha_2f_2^t-\alpha_1f_1^t)\big].\end{aligned}\ ] ]    theorem  [ theo5 ] states in particular that the scalar elastic theory in 3d is directly related to the non - scalar geometrically linear elasticity theory in 2d which was discussed in subsection  [ subexpl2d ] .    : we apply the translation method , see for instance @xcite , @xcite for an overview . in general , the function @xmath262 which describes the translation is only quasi - convex . here , as a benefit of the scalar theory
, we can pick @xmath262 as a convex function , cf .  .
starting from the inequality @xmath263 which is satisfied by any quasi - convex function @xmath262 , we obtain , similar to the reasoning in @xcite for the non - scalar case in @xmath10 dimensions , @xmath264.\end{aligned}\ ] ] consider @xmath265 then @xmath266 next we apply jensen s inequality which is possible since @xmath267 is convex for @xmath89 .
furthermore we possibly enlarge the set of admissible functions in the minimization over @xmath247 in by considering the set of admissible functions @xmath268 with the constraint .
then becomes @xmath269    if @xmath10 , @xmath262 is chosen to be @xmath270 .
here we set analogously @xmath271 since @xmath262 is quadratic , there exists a unique linear operator @xmath272 such that @xmath273 , here simply @xmath274 .
it holds @xmath275.\end{aligned}\ ] ] by a direct computation the remainder term @xmath276 can be rewritten , @xmath277\\ & = & -d_1d_2t(\nabla\eta_2-\nabla\eta_1)\!:\!(\nabla\eta_2-\nabla\eta_1)\\ & = & -2d_1d_2{\varphi}(\nabla\eta_2-\nabla\eta_1).\end{aligned}\ ] ] so we have found @xmath278 next we compute the optimal strains @xmath279 , @xmath280 . after differentiating the argument on the right in ( [ h2 ] ) for fixed @xmath281
, we obtain @xmath282 using the constraint @xmath283 , after rearrangement this yields the formulas @xmath284 setting @xmath285 ( see below ) , this coincides with ( [ etasdef1 ] ) , ( [ etasdef2 ] ) .    the maximum over @xmath281 in ( [ h2 ] )
is attained at @xmath286 as @xmath287 .
hence , the translation is trivial in this setting .
we obtain @xmath288 it remains to estimate @xmath40 from above which amounts to showing that @xmath279 , @xmath280 yield an optimal microstructure .
plugging in any microstructure @xmath120 with @xmath289 , from the definition ( [ wmins ] ) of @xmath40 , we get the upper bound @xmath290 as the domain of integration @xmath254 in the scalar case with @xmath1 is two - dimensional , finding the optimal microstructure subsumes to the non - scalar case in @xmath10 dimensions . hence ,
as can be seen from ( [ bsdef ] ) , @xmath286 occurs either for regime  0 ( no microstructure ) or regime  i. the latter is equivalent to @xmath291 . from a well - known argument ,
see ( * ? ? ?
* lemma  4.1 ) , this implies that @xmath214 , @xmath215 are compatible .
besides this compatibility condition , any optimal microstructure must also satisfy the equilibrium condition @xmath292\vec{n}=0\ ] ] which states that the jump of the stress in the normal direction of the laminate must vanish .
equation  ( [ state ] ) originates from the euler - lagrange equation of the variational problem . here , it is automatically satisfied since , due to ( [ nosjump ] ) , @xmath293\vec{n}=\sigma_2^*-\sigma_1^*=0 $ ] .
so there exists a unique rank - i lamination microstructure ( unique up to a sign @xmath294 of direction of the laminate ) that connects @xmath214 and @xmath215 , and the strain of phase @xmath102 is @xmath295 , @xmath89 .
this microstructure is optimal , @xmath296 , which proves ( [ whwsexp ] ) .
the methods of the previous subsection also apply to derive rigorously an explicit expression for @xmath40 in @xmath2 . here , quasi - convex functions are always convex , and the vector and scalar settings coincide .
however , we use the notation in analogy to the vector case .
[ theoa ] for @xmath2 , the functional @xmath40 , given by ( [ wmin ] ) , satisfies the explicit representation formula @xmath297 with @xmath298 the partial derivatives of @xmath40 are given by @xmath299}{d_2\alpha_1+d_1\alpha_2},\\ \frac{\partial{\widehat w}}{\partial d}(d,{\varepsilon } ) & = w_1({\varepsilon}_1^*)-w_2({\varepsilon}_2^ * ) + \frac{\alpha_1\alpha_2}{(d_2\alpha_1+d_1\alpha_2)^2}\big [ ( \alpha_1-\alpha_2){\varepsilon}^2+d_1\alpha_1({\varepsilon}_1^t)^2-d_2\alpha_2({\varepsilon}_2^t)^2{\nonumber}\\ & \hspace*{170pt } + \big(\alpha_2{\varepsilon}_2^t\!-\!\alpha_1{\varepsilon}_1^t \!+\!(\alpha_2\!-\!\alpha_1)(d_1{\varepsilon}_1^t\!+\!d_2{\varepsilon}_2^t)\big){\varepsilon}{\nonumber}\\ & \hspace*{170pt } \label{dwdd1d } + ( d_2\alpha_1-d_1\alpha_2){\varepsilon}_1^t{\varepsilon}_2^t\big].\end{aligned}\ ] ]    : we apply the same methods as in the proof of theorem  [ theo5 ] .
first , we estimate @xmath300 , @xmath301 from below by their convexification .
similar to ( [ h2 ] ) , we have for any convex function @xmath228 and @xmath302 , @xmath303    secondly , we choose the optimal @xmath281 and @xmath228 . the later is an art when @xmath304 . with arguments identical to the scalar case in @xmath1 , it holds @xmath305 , leaving us with the estimate @xmath306    next we compute the optimal strains @xmath307 , @xmath308 on the right . for @xmath309 , only @xmath310
is admissible and hence optimal .
now , let @xmath311 . we can then resolve the constraint by setting @xmath312 and need to calculate @xmath313 after differentiation with respect to @xmath307 and using ( [ widef ] ) for the optimal @xmath214 in ( [ star2 ] ) , we are left with @xmath314 rearrangement and resolution of @xmath315 gives the formulas ( [ e1db ] ) , ( [ e1dc ] ) of the optimal strains . since @xmath316 , these formulas also state the correct solution when @xmath309 .    in the final step , due to the definition of @xmath40 , it holds @xmath317 and @xmath122 may represent any microstructure with @xmath289 .
the explicit construction of the optimal microstructure for which the upper bound and the lower bound coincide , @xmath318 is again as in the proof of theorem  [ theo5 ] .
this leads to ( [ wh1d ] ) .
finally , the verification of ( [ dwdeps1d ] ) , ( [ dwdd1d ] ) follows from @xmath319 and @xmath320}{(d_2\alpha_1+d_1\alpha_2)^2},\\ \label{depss1d2b } \frac{\partial{\varepsilon}_2^*}{\partial d } & = \frac{\alpha_1\big[\alpha_2{\varepsilon}_2^t -\alpha_1{\varepsilon}_1^t+(\alpha_1-\alpha_2){\varepsilon}\big]}{(d_2\alpha_1+d_1\alpha_2)^2}\end{aligned}\ ] ] which all can be verified by elementary computations .
using the relationship @xmath321 we find ( [ dwdeps1d ] ) and finally @xmath322 } { ( d_2\alpha_1+d_1\alpha_2)^2}.{\nonumber}\end{aligned}\ ] ] using @xmath323 this simplifies to identity ( [ dwdd1d ] ) .
above we have already remarked that the energy @xmath324 satisfies assumption ( @xmath9 ) and thus theorem  [ theo1 ] applies . in this section
we investigate whether the results transfer to the geometrically linear theory of elasticity stated in the previous section .
the following statement is a consequence of theorem  [ theo1 ] .
[ thmcor1 ] let @xmath39 be a function on @xmath325\times{\mathbb{r}}^{d\times d}_\mathrm{sym}$ ] .
assume that the mobility tensor @xmath148 be positive definite and continuous for all @xmath4 , @xmath3 satisfying ( [ abrange ] ) .
moreover , let @xmath149 be given by ( [ psidef ] ) and the initial data @xmath150 satisfy ( [ abrange ] ) and ( [ w0bound ] ) .    in @xmath326 ,
let @xmath40 be given by ( [ whwsexp ] ) corresponding to the scalar setting .
if @xmath10 , let    1 .
@xmath327 be independent of @xmath131 , and 2 .
@xmath103 and @xmath241 commute whenever @xmath328 .
if @xmath2 , let @xmath40 be given by ( [ wh1d ] ) . then there exists a solution @xmath151 to ( [ ach1])([ach3 ] ) that satisfies    \(i ) @xmath152;\,l^2(\omega)\right)$ ] , + ( ii ) @xmath329 , @xmath330 , + ( iii ) @xmath155 , + ( iv ) the feasible parameter range of @xmath156 is given by ( [ abrange ] ) .
: we first consider the case @xmath0 .
the idea is to apply theorem  [ theo1 ] .
therefore we    1 .
verify that @xmath40 satisfies assumption ( @xmath9 ) for @xmath124 $ ] , 2 .
extend @xmath40 to @xmath331 in such a way that ( @xmath9 ) is fulfilled for all @xmath136 .
then theorem  [ theo1 ] can be applied and we obtain by part ( iv ) therein that effectively @xmath332 $ ] in the system ( [ ach1])([ach3 ] ) . hence theorem  [ thmcor1 ] follows with @xmath39 as asserted .
we start with the second step and give the proofs of the former in dimensions @xmath2 , @xmath10 , and the scalar case in @xmath1 thereafter .
we introduce the following extension @xmath333 of @xmath334 by @xmath335 for functions @xmath336 , @xmath337 that we may choose such that @xmath338 and that @xmath339 fulfils ( @xmath91)(@xmath92 ) for @xmath340 and @xmath341 .
we further mention that the function @xmath339 defined in ( [ wext ] ) is positive and coercive for @xmath342 .
we note that this construction is required for the first part of the proof of theorem  [ theo1 ] , where polynomial approximations of @xmath149 are constructed , where @xmath343 $ ] may occur .
it remains to show that @xmath40 fulfils ( @xmath9 ) for @xmath332 $ ] , which we prove independently for @xmath2 , @xmath10 , and the scalar case in @xmath1 .
\(i ) _ verification of ( @xmath9 ) for @xmath2 _
:    from  , ( [ wh1d ] ) , and , it follows @xmath344 .    from ( [ dwdeps1d ] )
, we obtain @xmath345 the condition ( @xmath91 ) , restricted to @xmath124 $ ] , follows , since by , for @xmath89 , @xmath346 where @xmath347 by assumption .
from , , for @xmath89 and @xmath124 $ ] we have by similar arguments latexmath:[\[\label{visest }    this leads to ( @xmath92)@xmath349 , since by ( [ wh1d ] ) , , @xmath350    for @xmath124 $ ] we easily compute , for @xmath89 , @xmath351    we use this estimate in ( [ star3 ] ) , to obtain by @xmath352 which is ( @xmath92)@xmath353 restricted to @xmath124 $ ] .    for @xmath124 $ ]
we derive with the help of ( [ dwdeps1d ] ) , ( [ star4 ] ) , @xmath354 showing the validity of ( @xmath92)@xmath355 restricted to @xmath124 $ ] .
\(ii ) _ verification of ( @xmath9 ) for @xmath10 _ :    the regularity of @xmath40 required for ( @xmath9 ) follows directly from and the related definitions .
the constant @xmath222 solely depends on @xmath356 , @xmath140 since @xmath357 and @xmath218 is the reciprocal of the largest eigenvalue of @xmath219 , see . since @xmath224 $ ] , latexmath:[\[\label{bsest }    independent of @xmath11 and @xmath131 .    from ( [ bsest ] ) , ( [ epsistardef ] ) we obtain for @xmath124 $ ] latexmath:[\[\label{e1d2 }     with   and
this shows for @xmath124 $ ] @xmath360 from this we easily verify ( @xmath92)@xmath349 restricted to @xmath124 $ ] , as the first two terms in ( [ wdef ] ) can be estimated as in the one - dimensional case .    in order to show ( @xmath92)@xmath353 , we apply and ( [ visest ] ) to find @xmath361 finally , from ( [ betrus ] ) and the fact that @xmath241 is an isometry w.r.t .
the frobenius norm , @xmath362 the terms @xmath363 , @xmath89 remaining in can be estimated as in the one - dimensional case .
this verifies ( @xmath92)@xmath353 restricted to @xmath124 $ ] .
now we prove ( @xmath92)@xmath355 . for the regimes 0 , i and ii this follows directly from ( [ visest ] ) . when in regime  iii , we start from   and find with   for @xmath124 $ ] @xmath364    the validation of ( @xmath91 ) relies on . by assumption ( i ) of the theorem , @xmath365 and , as @xmath124
$ ] is fixed in ( @xmath91 ) , @xmath366 is a constant tensor .
assumption  ( @xmath91 ) then follows from @xmath367 a similar equality for @xmath368 , and the positive definiteness of @xmath366 and @xmath369 , @xmath89 .
indeed , the latter is equivalent to @xmath370 positive definite , and since @xmath371 , this is implied by all eigenvalues of @xmath372 being positive .
however , from ( [ gsdef ] ) and the definition of @xmath218 , this is true .    for the scalar case in @xmath1 , comparing ( [ wdef ] ) and ( [ whwsexp ] ) and
since the partial derivatives coincide , the theory subsumes to the geometrically linear theory for @xmath10 and @xmath286 .
hence we can proceed as above and show that @xmath40 given by ( [ whwsexp ] ) satisfies ( @xmath9 ) .
the assumptions ( @xmath3731 ) , ( @xmath3732 ) are not needed here because @xmath286 .
the following theorem states the uniqueness of weak solutions for the system with geometrically linear theory .
[ theo4 ] let @xmath39 be given by ( [ whwsexp ] ) for @xmath1 ( scalar setting ) , by ( [ wdef ] ) for @xmath10 , or by ( [ wh1d ] ) for @xmath2 . assume further that @xmath86 and the elastic moduli of the two phases are equal , i.e. @xmath374 . for @xmath10 , let @xmath375 . then the solution @xmath151 in theorem  [ thmcor1 ] is unique in the spaces stated there .    :
( i ) _ the case @xmath2 . _    for @xmath374 , the equations ( [ e1db ] ) , ( [ e1dc ] ) read @xmath376 which implies @xmath377 leading to @xmath378 .
this yields @xmath379 in this special case , the partial derivatives of @xmath40 are @xmath380 now we pass through the steps in the proof of theorem  [ theo2 ] and make the necessary modifications .
let again @xmath381 , @xmath190 be two solutions to ( [ ach1])([ach3 ] ) and set @xmath191 , @xmath192 , @xmath382 , @xmath383 for @xmath190 and @xmath384 , @xmath385 .
the equations ( [ test2 ] ) , ( [ test3 ] ) now read @xmath386{\,\mathrm{d}x}{\,\mathrm{d}t}{\nonumber}\\ \label{2test2 } & \quad = \!{\int\limits_{\omega_{\cal t}}}\left[\frac{{\vartheta}}{2}\left(g'(a^2\!+\!b^2)-g'(a^1\!+\!b^1 ) + g'(a^1\!-\!b^1)-g'(a^2\!-\!b^2)\right)\eta+2\kappa_2b\eta\right]{\,\mathrm{d}x}{\,\mathrm{d}t},\\ \label{2test3 } & { \int\limits_{\omega_{t_0}}}\alpha_1\big({\varepsilon}+d({\varepsilon}_2^t-{\varepsilon}_1^t)\big){\varepsilon}(u){\,\mathrm{d}x}{\,\mathrm{d}t}=0.\end{aligned}\ ] ] analogous to ( [ step2 ] ) we find @xmath387|a|{\,\mathrm{d}x}{\,\mathrm{d}t}.\end{aligned}\ ] ] choosing @xmath205 in ( [ 2test2 ] ) and adding the resulting equation to ( [ 2step4 ] ) , we obtain @xmath388{\,\mathrm{d}x}{\,\mathrm{d}t}\\ & \le\!{\int\limits_{\omega_{t_0}}}\!2(\kappa_1|a|^2\!+\!\kappa_2|b|^2 ) + \!\!{\int\limits_{\omega_{t_0}}}\!\!\frac{{\vartheta}}{2}\big[\big|g'(a^1\!\!+\!\!b^1 ) \!-\!g'(a^2\!\!+\!\!b^2)\big|\!+\!\big|g'(a^1\!\!-\!\!b^1)\!-\!g'(a^2\!\!-\!\ ! b^2)\big|\big](|a|\!+\!|b|){\,\mathrm{d}x}{\,\mathrm{d}t}{\nonumber}.\end{aligned}\ ] ] it is noteworthy that we are only able to estimate @xmath389 which differs from the total mechanical energy @xmath390 .    with the lipschitz continuity of @xmath208 , and by applying the inequalities of young and gronwall to ( [ wnote ] ) , we obtain as in the proof of theorem  [ theo2 ] @xmath391 with korn s inequality , @xmath392 in @xmath58 .
\(ii ) _ the scalar case in @xmath1 . _    by theorem  [ theo5 ] , for @xmath374 , @xmath393 consequently , @xmath394 leading to @xmath395 . since @xmath286 ,
the formula for @xmath40 has now the same structure as for @xmath2 and the proof follows exactly as in the first part ( i ) .
\(iii ) _ the case @xmath10 . _    by assumption , @xmath396 , which results in @xmath397 in @xmath58 . in addition , for @xmath374 , by eqn .
( [ epsistardef ] ) , @xmath398 consequently , @xmath399    so again @xmath400 , @xmath378 , and the proof can be carried out as before .
in this article we derived extensions of the allen - cahn / cahn - hilliard system to elastic materials .
in particular we included ( i ) the linear elastic energy derived by eshelby , ( ii ) a geometrically linear theory of elasticity for @xmath6 , and ( iii ) a geometrically linear theory of elasticity for @xmath1 in the scalar setting .
for all three cases we showed in a mathematically rigorous way the existence and the uniqueness of weak solutions , asserting the correctness of our approach . as a future goal , it is desirable to generalize the proof to the non - scalar , i.e. vectorial case in @xmath1 and to @xmath401 .
the generalized ac - ch models contain as special cases both the allen - cahn @xcite and the cahn - hilliard equation @xcite with linear elasticity . in our work
the existence and uniqueness of weak solutions to the ac - ch model with linear elasticity follows in a straightforward way since assumption  ( @xmath9 ) required in the existence result theorem  [ theo1 ] can easily be checked to hold true for the linear energy .
this had not been shown before .
we point out that for the new cases ( ii ) and ( iii ) , the formulas collected in section  [ secexpl ] for the specific energy @xmath40 are essential for any numerical studies of the ac - ch models extended to microstructure . for related investigations and numerical methods
we refer the interested reader to @xcite and @xcite .    besides the limiting assumption of constant temperature ( cf .
introduction ) , the most important pending restriction is the postulation of small strain , included in ( [ vedef ] ) .
that is , it would be desirable to combine the ac - ch system with energies within the framework of the geometrically non - linear theory of elasticity .
the technical problems of a large strain theory are striking , see , e.g. , @xcite , and it is in many cases not known how to compute explicit formulas for the relaxed energy functionals .
similar mathematical difficulties arise when one wants to extend the ac - ch systems to elastic materials which are of polycrystalline structure .
this is an open problem for future research . for related modelling aspects we refer to @xcite , for related mathematical aspects to @xcite
we have worked on this project while we were employed or were visitors at the max planck institute for mathematics in the sciences , leipzig ; the institute of applied mathematics , university of bonn ; and the institute for mathematics , university of wrzburg .
some part of this work was performed while tb was visiting the hausdorff research institute for mathematics , bonn .
we acknowledge the hospitality of all these institutions .
99 allen , s.  m. ; cahn , j.  w. microscopic theory for antiphase boundary motion and its application to antiphase domain coarsening , _ acta metallurgica _ 1979 , * 27 * , 10851095 .
ambrosio , l. ; fusco , n. ; pallara , d. _ functions of bounded variation and free discontinuity problems , _ clarendon press , new york , 2000 .
barrett , j.  w. ; garcke , h. ; nrnberg , r. a phase field model for the electromigration of intergranular voids .
_ interfaces and free boundaries _ 2007 , * 9 * , 171210 .
bhattacharya , k. comparison of the geometrically nonlinear and linear theories of martensitic transformation .
_ contin .
_ 1993 , * 5 * , 205242 .
bhattacharya , k. ; kohn , r.  v. elastic energy minimization and the recoverable strains of polycrystalline shape - memory materials .
.  rational mech .
_ 1997 , * 139 * , 99180 .
bhattacharya , k. ; schlmerkemper , a. stess - induced phase transformations in shape - memory polycrystals , _ arch .  rational mech .
_ 2010 , * 196 * , 715751 .
blesgen , t. ; chenchiah , i.  v. a generalized cahn - hilliard equation based on geometrically linear elasticity , _ interfaces and free boundaries _ 2011 , * 13 * , 127 .
blesgen , t. the elastic properties of single crystals with microstructure and applications to diffusion induced segregation . _
crystal research and technology _ 2008 , * 43 * , 905913 .
blesgen , t. ; weikard , u. on the multicomponent allen - cahn equation for elastically stressed solids .
_ electr .
j.  diff .
_ 2005 , * 89 * , 117 .
blesgen , t. ; schlmerkemper , a. towards diffuse interface models with a nonlinear polycrystalline elastic energy . in : a.
doughett and p. asnarez ( eds . ) , composite laminates : properties , performance and applications , nova science publishing 2010 , pp .  465 - 489 .
brochet , d. ; hilhorst , d. ; novick - cohen , a. finite - dimensional exponential attractor for a model for order - disorder and phase separation .
_ 1994 , * 7 * , 8387 .
cahn , j.  w. ; hilliard , j.  e. free energy of a non - uniform system i. interfacial free energy .
_ j. chem .
_ 1958 , * 28 * , 258267 .
cahn , j.  w. ; larch , f.  c. the effect of self - stress on diffusion in solids .
_ acta metall .
_ 1982 , * 30 * , 18351845 .
cahn , j.  w. ; novick - cohen , a. evolution equations for phase separation and ordering in binary alloys .
_ j.  stat .
_ 1994 , * 76 * , 877909 .
chenchiah , i.  v. ; bhattacharya , k. the relaxation of two - well energies with possibly unequal moduli .
rational mech .
_ 2008 , * 187 * , 409479 .
cherkaev , a. v. variational methods for structural optimization , _ appl .
_ 2000 , * 140 * , springer publishing .
cioranescu , d. ; donato , p. _ an introduction to homogenization , _ oxford university press , 1999 .
dacorogna , b. _ direct methods in the calculus of variations , _ springer , berlin , 1989 .
eshelby , j.  d. elastic inclusions and inhomogeneities .
solid mech .
_ 1961 , * 2 * , 89140 .
fratzl , p. ; penrose , o. ; lebowitz , j.  l. modelling of phase separation in alloys with coherent elastic misfit .
_ j.  stat .
_ 1999 , * 95 * , 14291503 .
garcke , h. on cahn - hilliard systems with elasticity .
edinburgh sec .
a _ 2003 , * 133 * , 307331 .
gottstein , g. physical foundations of materials science .
springer berlin , 2004 .
kohn , r.  v. ; vogelius , m. relaxation of a variational method for impedance computed tomography .
pure appl .
_ 1987 , * 60 * , 745777 .
kohn , r.  v. the relaxation of a double - well energy .
_ continuum mechanics and thermodynamics _ 1991 , * 3(3 ) * , 193236 .
kohn , r.  v. ; niethammer , b. geometrically nonlinear shape - memory polycrystals made from a two - variant material .
_ 2000 , * 34 * , 377 - 398 .
milton , g.w . a brief review of the translation method for bounding effective elastic tensors of composites , _ cont . models and discrete systems
maugin , ed . ) , longman scientific and technical , 1990 , 6074 .
novick - cohen , a. triple - junction motion for an allen - cahn / cahn - hilliard system .
_ physica d _ 2000 , * 137 * , 124 .
nye , j.  f. _ physical properties of crystals : their representation by tensors and matrices , _ clarendon press ; new york : oxford university press , 1984 .
onuki , a. ginzburg - landau approach to elastic effects in the phase separation of solids .
_ j.  phys .
japan _ 1989 , * 58 * , 30653068 .
salje , e. phase transitions in ferroelastic and co - elastic crystals .
cambridge university press , 1990 .
ziemer , w. _ weakly differentiable functions : sobolev spaces and functions of bounded variation , _ graduate texts in mathematics , springer , new york , 1989 .