transformations between different solid states of a material are essential for many technological and scientific applications , and understanding their kinetics is therefore not only interesting because of the scientific variety and beauty , but also important for tailoring of new materials with specific properties . often , these processes are accompanied by elastic deformations , which can be e.g.  due to density differences or structural changes , which provoke stresses at interfaces between adjacent phases .
the influence on the thermodynamics of the transitions between different phases has been thoroughly discussed in the literature @xcite , whereas the kinetics of these processes are still far less understood .    facing this lack of understanding
, we would like to mention that very fruitful concepts have been developed to describe the dynamics of solidification or melting processes as moving boundary problems .
it means that the evolution of the different phases is described on the very detailed level of the motion of the interfaces between them , for which suitable equations of motion have to be posed .
a particular application is the famous example of dendritic growth , where heat or component diffusion is the rate limiting process . at the interface between the solid and melt phase the temperature or concentration is close to its phase equilibrium value .
the concentration jump or the release of latent heat due to partitioning leads to a jump in the corresponding fluxes at the interfaces , and the magnitude of this jump is directly related to the local interface velocity . altogether , this provides a closed description for the motion of each interface point , and many efforts have been done towards analytical and numerical solutions of this problem .
it turns out that the questions concerning existence , shape and growth velocity of steady state patterns in solidification crucially depend on selection mechanisms . in this spirit
the influence of many different physical effects on solidification or melting processes has been studied throughout the years @xcite .
for example , the effect of isotropic surface tension was proven not to serve as a selection mechanism for a solution , whereas anisotropic surface tension leads to a unique solution .
also , the appearance of triple junctions can provide a selection mechanism , and can therefore lead to a solution of the free boundary problem @xcite .
of course , the situation considering solid - solid transformation kinetics differs significantly from the solidification or melting problem , where elastic effects often play only a minor role .
still , we follow here the same successful concepts to formulate a free boundary description for the transformation kinetics in solids , now including the important effect of elasticity . in real systems
many other effects apart from elastic effects come into play  like inhomogeneous compositions , the mullins - sekerka instability , crystal anisotropies , polycrystalline structures , etc .  and also
quite successful attempts have been made to model this complex behavior ( see for example in @xcite and references therein ) . here , in contrast , we focus on gaining a fundamental understanding of the elastic influence on pattern formation processes , and we restrict our investigations therefore to simple and  clean " model systems .
correspondingly , we consider only a single diffusion field , namely the thermal diffusion in pure materials .
the formulation of such a model system requires to write down explicit expressions for the boundary conditions and suitable equations of motion at the propagating interfaces , which therefore have to be tracked during the entire dynamical process . finally , it provides a detailed description of the microstructure evolution on the level of individual grains in a polycrystalline structure . as for the dendritic growth , we focus here on the regime of diffusion limited growth .
then , propagation of a front generates or consumes latent heat , which has to diffuse through the system , and we assume that this is the rate - limiting process .
this is in contrast to isothermal solid - solid transformations that were discussed in @xcite , which are driven by local interface kinetics , therefore allowing fast front propagation also on the scale of the sound speed .
many important properties result from the influence of elastic effects due to structural transformations and density differences .
apart from sharp interface methods , that are based on green s function techniques , we also use phase field descriptions , which are nowadays widely used .
more importantly , however , is that we supplement these numerical investigations by analytical calculations which allow to understand these processes and the relevant parameters on a deeper basis .
specifically , we consider the growth of an equilibrium solid phase @xmath0 from another metastable solid phase @xmath1 , as depicted in fig .
[ fig : geometry ] . at a temperature @xmath2 ,
both stress - free phases are in bulk equilibrium with each other , whereas below this value the new @xmath0 phase is thermodynamically favorable .
far ahead of the propagating front , the system is therefore exposed to a lower temperature @xmath3 , and we define a dimensionless temperature @xmath4 , where @xmath5 denotes the latent heat of the transition and @xmath6 the heat capacity . within the solid phases the evolution of the temperature field obeys a diffusion equation , @xmath7 where @xmath8 is the thermal diffusivity , which we assume to be equal in both phases , as well as the heat capacity ( symmetrical model ) .
the motion of the interfaces follows from energy conservation principles , which account here for the net heat flux and the generation or absorption of latent heat due to the local motion of the interface with normal velocity @xmath9 , @xmath10}-\nabla w_{int}^{[\alpha ] } ) , \label{eq : heat - conservation - equation}\end{aligned}\ ] ] where @xmath11 is the interface normal and the superscripts denote the phases .
elastic effects enter into the problem via a modification of the local phase equilibrium condition at the moving interfaces .
in contrast to previous investigations @xcite , we now assume the phases to be coherently connected , i.e.  we neglect the appearance of any slips , detachments or defect formation at the interfaces , which is reasonable for small misfits .
it is important to mention here , that with the assumption of coherency we formulate a completely different model being capable to describe qualitatively different types of elastic effects .
for example the well - known effect of elastic hysteresis @xcite , which was proven to disappear without coherency @xcite , is now included in the model .    in a lagrangian description of linear elasticity , the displacement field
is denoted by @xmath12 , and the coherency condition reads , @xmath13}=\mathbf{u}^{[\beta]}$ ] at the interface . the strain field
is then given by the symmetrical spatial derivative , @xmath14 .
this implies that tangential and shear strains @xmath15 are continuous at the interface with @xmath16 being the two tangential directions .
force balance demands the continuity of the normal and shear stresses , @xmath17 , where @xmath18 is the normal direction . here , the stresses are defined as the derivative of the free energy with respect to the strains , @xmath19 .    , of a thermodynamically favored phase @xmath0 into a metastable phase @xmath1 , which is initially at temperature @xmath20 .
the motion of the interface is locally described by the interface - normal velocity @xmath21 , which is given by the local heat conservation eq .
( [ eq : heat - conservation - equation ] ) .
the dashed lines at the two side borders indicate that we will discuss two different types of outer boundary conditions : free growth , where the strains decay for @xmath22 , and channel growth , with fixed displacements and thermal insulation at the walls . in the growth direction
@xmath23 the system is assumed to be infinite in both cases . [
fig : geometry],scaledwidth=55.0% ]    the interface temperature for a planar , stress free interface is given by @xmath24 , with @xmath25 .
capillary corrections induce a curvature dependent term , and for simplicity we discuss only the case of isotropic surface tension @xmath26 .
finally , elastic deformations lead to an additional change of the interface temperature . altogether , this leads to the local equilibrium boundary condition @xcite @xmath27 with @xmath28 being the capillary length and @xmath29 the local curvature of the interface , counted positive for a convex phase @xmath0 ( see fig .
[ fig : geometry ] ) ; notice that due to the coherency constraint @xcite , the elastic shift of the equilibrium temperature is proportional to the difference of a new potential @xmath30 } - \tilde{f}_{el}^{[\beta]}$ ] , which is related to the elastic free energy @xmath31 via the legendre transformation @xmath32 .
here , we use the sum convention for repeated indices for the sake of brevity .
from now on , we restrict our considerations to linear isotropic elasticity . taking the relaxed state of the @xmath1 phase as reference , the elastic contribution to the free energy there reads @xmath33 } & = \frac{e}{2(1+\nu)}\left(\frac{\nu}{1 - 2\nu}\epsilon_{ii}^{2}+\epsilon_{ik}^{2}\right),\label{eq : free - energy - reference}\end{aligned}\ ] ] where @xmath34 and @xmath35 are young s elastic modulus and poisson ratio , respectively .
in contrast , the elastic contribution to the free energy density of the new phase @xmath0 reads @xmath36 } & = \frac{e}{2(1+\nu)}\left(\frac{\nu}{1 - 2\nu}\left(\epsilon_{ii}-\epsilon_{ii}^{0}\right)^{2}+\left(\epsilon_{ik}-\epsilon_{ik}^{0}\right)^{2}\right),\label{eq : free - energy - new}\end{aligned}\ ] ] which has an other state of zero elastic energy due to the lattice strain @xmath37 assigned to the phase transformation ( see for example @xcite ) . in the following
we will call @xmath37 the eigenstrain of the transformation .
particular cases of eigenstrains and their interpretation will be discussed in the next section .
the mechanical equilibrium and coherency conditions provide expressions for the discontinuous jumps of the strains at the interface @xmath38}-\epsilon_{nn}^{[\alpha]}= & \epsilon_{nn}^{0}+\frac{\nu}{1-\nu}(\epsilon_{ss}^{0}+\epsilon_{\tau\tau}^{0}),\label{eq : strain - jump - n}\\ \epsilon_{n\tau}^{[\beta]}-\epsilon_{n\tau}^{[\alpha]}= & \epsilon_{n\tau}^{0},\quad\epsilon_{ns}^{[\beta]}-\epsilon_{ns}^{[\alpha]}=\epsilon_{ns}^{0},\label{eq : strain - jump - tau / s}\end{aligned}\ ] ] which are only related to the imposed eigenstrain .
defining an eigenstress tensor via hooke s law @xmath39 we obtain after a few straightforward algebraic manipulations for the elastic contribution to the local equilibrium condition eq .
( [ eq : phase - equilibrium - condition ] ) @xmath40 } - \frac{e}{2(1-\nu^{2})}\left((\epsilon_{\tau\tau}^{0})^{2}+(\epsilon_{ss}^{0})^{2}+2\nu\epsilon_{\tau\tau}^{0}\epsilon_{ss}^{0}+2(1-\nu)(\epsilon_{\tau s}^{0})^2\right ) .
\label{eq : phase - equilibrium - elast - cont}\end{aligned}\ ] ] apparently , this expression depends on the strain state of the @xmath1 phase at the interface .    as already implicitly assumed above
, the elastic degrees of freedom relax fast on diffusive timescales , and therefore the application of static elasticity is legitimate .
therefore , newton s second law becomes @xmath41    the above set of equations now has to be supplemented by boundary conditions at the external system boundaries , and we will discuss different scenarios below .
they complete the self - sonsistent description of the moving boundary problem , which require the simultaneous solution of the elastic and diffusion equations , coupled via the boundary conditions at the propagating interfaces , leading to a complicated nonlinear and nonlocal problem .
we note that we will later also discuss more complicated situations with more than one growing phase .
then , also similar boundary conditions and equations have to be set up for the interfaces between them .
typically , then also triple junctions will appear , where the contact angles are given by young s law .
the simplest case of a characteristic lattice strain is to assume the bond length of the new phase @xmath0 to be uniformly longer or shorter in all directions in comparison to the reference phase , i.e.  @xmath42 , which implies a density change .
notice , that this kind of transformations , to which we will refer to as _ dilatational eigenstrain _ , only changes the volume of an elementary cell , and not its relative dimensions or its scaled shape
.    the opposite case are transformations , where the volume of the unit cell is conserved and only the shape is changed .
these _ shear transformations _ are characterized by a traceless eigenstrain tensor . a particular transition involving shear
strain can occur e.g.  in hexagonal crystals , where the transition leads to a lowering of the symmetry from @xmath43 to @xmath44 .
this is the case for example in hexagonal - orthorombic transitions in ferroelastics ( see @xcite and references therein ) .
let the principal axis @xmath45 be oriented in @xmath46 direction . for simplicity
we neglect all other possible strains with higher ( axial ) symmetry . by proper choice of the crystal orientation around the main axis in the initial phase ,
we obtain the new phase in three possible states in the @xmath47 plane due to the original hexagonal symmetry ( see fig .  [
fig : hexagonaltoothorombic ] ) .
visualisation of the hexagonal to orthorhombic transformation . from the symmetry
its obvious , that the orientation angle @xmath48 of the new orthorhombic phase can have only the three different values @xmath49.,scaledwidth=60.0% ]    then the nonvanishing components of the eigenstrain tensor @xmath50 are @xmath51 and @xmath52 , where the orientation angle @xmath48 is either @xmath53 or @xmath54 .
we also consider a more generic mixed mode case , where the self - strain of the transformation is a linear combination of the above two different transformations .
we will discuss the following two mixing cases @xmath55 and refer to them as positive or negative mixing .
it is important to mention that switching from positive to negative mixing corresponds to a rotation of the eigenstrain tensor by an angle @xmath56 . since our model is isotropic in all respects except the choice of the eigenstrain tensor , the only predefined direction is given by the growth velocity
therefore , we can interpret the negative mixing case to represent physically the same system as the positive mixing case but with a growth direction perpendicular to it .
as a first application , we discuss the growth of the stable @xmath0 phase from the metastable @xmath1 phase in an infinitely large system , which is exposed to the temperature @xmath3 far away from the interface . for simplicity , we discuss an effectively two - dimensional infinite system , assuming translational invariance in @xmath46 direction . in particular , we assume a plane strain situation from point of view of elasticity . also , we choose the elastic boundary conditions such that all stresses decay far away from the interface .    seeking for steady state solutions in free space ,
we use green s function methods for solving both the elastic and the diffusional problem , to finally obtain self - consistently the shape and the velocity of the growing front . while the derivation of the integral representation of the thermal field is well - known , see e.g.  @xcite , we present the calculation of @xmath57 in more detail @xcite . to this end , we express the elastic problem using the eigenstrain - independent stress tensor @xmath58 ; for the @xmath1 phase it is the same as the nomimal stress , @xmath59}=\sigma_{ik}^{[\alpha]}$ ] , whereas in the @xmath0 phase it is @xmath60}=\sigma_{ik}^{[\beta]}+\sigma_{ik}^{0}$ ] . in terms of the new stress tensor @xmath58 , the mechanical equilibrium equation
@xmath61 introduces a force density which is localized on the interface and vanishes in the bulk , @xmath62}-\tilde{\sigma}_{in}^{[\alpha]}=\sigma_{in}^{0}. \label{eq : artificial - interface - forces}\end{aligned}\ ] ] this means that the problem of two coherently connected solids with different reference state is equivalent to a monolithic material with a distribution of point forces along the interface .
we can therefore use an integral representation of the displacement field , @xmath63 where the integration is performed along the interface , and the free - space green s tensor @xmath64 of the elastic problem is given by @xcite @xmath65 this function is the displacement response at @xmath66 to a point force located in the origin . as a result , the mismatch provokes long - ranged elastic deformations , and therefore induces a nonlocal character to the problem .
the above expressions can then be used to calculate the elastic influence on the interface temperature @xmath67 .
far behind the tip the interface profile becomes parabolic and is described by the asymptotic ivantsov parabola @xcite . in this region ,
the elastic strains have decayed , and only the constant local contribution to @xmath68 ( second term in eq .
[ eq : phase - equilibrium - elast - cont ] ) remains .
consequently , the effective driving force @xmath69 consists not only of the thermal undercooling @xmath70 , which is controlled by the applied temperature @xmath3 , but also a contribution @xmath71 from the elasticity , @xmath72}{2(1-\nu^{2})l^2}. \label{eq : definition - deltael}\end{aligned}\ ] ] we can consider @xmath71 as a parameter quantifying the strength of the elastic effects .
the elastic effects therefore induce a shift of the equilibrium temperature , and the effective undercooling @xmath69 is lower than the actually applied thermal driving force , @xmath73 which reflects the elastic hysteresis @xcite .    by elimination of the thermal field @xcite we can derive the steady state equation for the shape of the solid - solid interface in a closed dimensionless representation , which reads in the co - moving frame of reference @xcite , @xmath74k_{0}(p\mu(x , x ' ) ) .
\label{mainpartclosedequation}\end{aligned}\ ] ] here , @xmath75 denotes the peclet number ( @xmath76 is the steady state velocity ) , and it is related to the driving force by the implicit relation @xmath77 .
furthermore , @xmath78 is the modified bessel function of third kind in zeroth order , and @xmath79^{1/2}$ ] . in the above expression we rescaled all lengths by the radius of curvature @xmath80 of the asymptotic ivantsov parabola .
consequently , the central equation ( [ mainpartclosedequation ] ) , together with the expression for the elastic contribution ( [ eq : phase - equilibrium - elast - cont ] ) , where the strains are expressed via the integral representation ( [ formal_solution_displacements ] ) results in a complicated integro - differential equation for the interface shape .
it differs from conventional dendritic growth by the appearance of the elastic term . to better understand its influence
, we briefly review the solvability conditions for classical dendritic growth @xcite , for situations with pure dilatational and shear strain @xcite , and cases with a triple junction @xcite : in classical dendritic growth , any solution for isotropic surface tension with nonzero solvability parameter @xmath81 must exhibit a finite negative tip kink angle @xmath82 .
for this reason , dendritic growth is not possible for isotropic surface energy , since no smooth dendrite profiles exist .
a dilatational eigenstrain simply shifts the equilibrium interface temperature , and therefore suffers from the same missing selection as the dendrite problem without surface tension anisotropy .
the presence of a shear eigenstrain , which will also be discussed in more detail below , leads to the propagation of a bicrystal . in this case
, the elastic effects operate as new selection mechanism .
independent of this elastically induced selection mechanisms , the presence of a triple junction alone also works as a selection mechanism .
we solve the eigenvalue problem eq .
( [ mainpartclosedequation ] ) numerically for the steady state shape and the stability parameter @xmath83 in the spirit of ref .
@xcite , and the results are presented in section [ results ] .
the growth of the @xmath0 phase in a channel of finite width @xmath84 differs significantly from the free growth scenario discussed above .
the specific behavior depends of course on the boundary conditions , and we assume fixed grip conditions on the long sides of the channel , i.e. the displacement is constant .     phases growing with constant velocity @xmath76 into the metastable @xmath1 phase in a finite channel of width @xmath84 .
the shapes in @xmath85 and @xmath86 correspond to the dilatational case ( @xmath87 ) , with @xmath88 and an undercooling , which is chosen such that according to eq .
( [ eq : fraction - of - beta - phase ] ) it would correspond to an asymptotic phase fraction @xmath89 . while @xmath85 shows the shape of the symmetrical finger , @xmath86 shows the shape of the sightly faster parity broken finger ( see fig .
[ fig : diagonaleigenstrain]a ) . in panel
@xmath90 we have plotted the shape in the case of negative mixing eigenstrain , @xmath91 , which in that case corresponds to growth in the privileged direction .
the parameters are @xmath92 and @xmath93 ( for the resulting velocity see fig .  [
fig : diagonaleigenstrain]b ) .
, scaledwidth=65.0% ]    again , we discuss an effectively two - dimensional system by the assumption of translational invariance in @xmath46 direction ( see fig .  [ dilatationalsetup ] for the setup ) .
furthermore , we assume that the @xmath1 phase far in front of the growing tip is fully relaxed , which is in agreement with the boundary conditions @xmath94 and @xmath95 along the boundaries of the whole , infinitely long strip . from a thermal point of view , we assume thermal insulation along these boundaries , so no temperature flux through the surfaces , @xmath96 on the outer boundaries .
far ahead of the tip , the temperature is given by the dimensionless undercooling @xmath70 .
due to the phase transition , the release of the latent heat leads to a heating of the system , therefore the temperature far behind the tip @xmath97 is higher than in front . under the given circumstances it is quite obvious that in this tail region a homogeneous state is reached . in particular ,
since the interfaces do not move there , the temperature becomes indeed constant .    for simplicity
, we discuss here only the specific case of a purely diagonal eigenstrain , i.e.  @xmath98 for @xmath99 , to demonstrate the analytical procedure .
in contrast , in the case of an eigenstrain tensor with nonzero off - diagonal elements , bicrystal patterns emerge , and we will briefly discuss this situation below ; it is quite obvious that the analysis can be performed in an analogous way .    due to the confinement of the system , phase coexistence is possible in a finite range of undercoolings , and
we can calculate the asymptotic volume fraction @xmath100 of the new phase ( see fig .
[ dilatationalsetup ] ) .
notice , that in contrast to the isothermal situation discussed in @xcite , the temperature in the tail region is not a control parameter but has to be found self - consistently , due to the thermally insulating boundaries . in terms of the asymptotic phase fraction @xmath100
, the width of the @xmath0 phase in the tail region is written as @xmath101 .
mechanical compatibility requires @xmath102}+{\lambda}\epsilon_{xx}^{[\beta ] } & = 0,\label{eq : mixedmode - coherency - cond_yy}\end{aligned}\ ] ] since the average lateral strain must vanish due to the boundary conditions for the displacement .
together with the mechanical equilibrium condition of force balance , we can then calculate the elastic fields in the tail region , and consequently the elastic energy .
then , we calculate the energy excess @xmath103 , due to the appearance of a finite fraction of the new phase @xmath100 .
it is defined as the difference between the energy for @xmath104 without any phase @xmath0 and the energy with an arbitrary but finite fraction , @xmath105}+(1-{\lambda})f_{el}^{[\alpha]}\right)-2\gamma+w\int\limits _ { 0}^{\lambda}\left(\frac{l(t_{eq}-t(\lambda'))}{t_{eq}}\right)d\lambda'.\end{aligned}\ ] ] initially , there is only the elastically relaxed but metastable phase @xmath1 at the constant temperature @xmath70 .
as soon as we have a finite phase fraction @xmath100 , the energy increases by elastic and capillary contributions .
finally , the integral appears due to thermal insulation , since an increase of the amount of the phase @xmath0 by phase transformation processes generates latent heat , and therefore causes an increase of the temperature . using the heat conservation condition we find the relation between this temperature and the phase fraction to be @xmath106 therefore , we obtain for the energy excess as a function of the asymptotic phase fraction @xmath100 @xmath107 where the parameter @xmath71 , as already introduced above in eq .
( [ eq : definition - deltael ] ) , defines the strength of elastic effects , and @xmath108 is a parameter depending on the type of eigenstrain .
the maximum of this energy excess , @xmath109 , provides the equilibrium fraction , which expresses the thermodynamic balance between elastic deformation and free energy release due to the phase transition , @xmath110 if the fraction @xmath111 is in the range @xmath112 , we obtain phase coexistence in the asymptotic regime far behind the tip .
we would like to point out that one can derive this equation also by using the phase equilibrium condition eq .
( [ eq : phase - equilibrium - condition ] ) . then growth of the @xmath0 phase demands that the energy excess ( [ abc123 ] ) for the equilibrium fraction ( [ eq : fraction - of - beta - phase ] ) has to be positive , @xmath113 this is equivalent to the condition @xmath114 , where @xmath115 is given by @xmath116 we note that for the shear transformations discussed above , the most favorable process is the growth of a bicrystal , where two twinned phases @xmath0 and @xmath117 with @xmath118 grow together , because this is energetically most favorable @xcite . for the bicrystal configuration , the analysis can be performed in a similar way , and then also the @xmath119 grain boundary energy appears in the energy balance analogous to eq .
( [ qwert ] ) .
considering the channel geometry , we use a phase field formulation to solve the full free boundary problem numerically .
this method is very versatile and allows to study even more complicated scenarios , e.g.  the full dynamical evolution and not only the steady state regime , three - dimensionsal situations with full anisotropy of elasticity and surface tension or simulations beyond the symmetrical model , where the two phases have different diffusion constants .
nevertheless , for brevity of presentation , we confine our investigations to the simple , isotropic , two - dimensional and thermally symmetric case that we also used for the green s function approach in the free space geometry and the above analytical calculations .
the phase field method introduces an additional numerical lengthscale @xmath120 , the width of the transition region between the different solid phases .
this means that an order parameter , which discriminates between the different phases , has a profile @xmath121/2 $ ] in normal direction through an interface .
we operate here in the sharp interface limit , which means that this lengthscale is significantly smaller than all other physical lengthscales present in the problem .
then the description recovers the desired sharp interface equations that were proposed above .
hence , we introduce the phase field @xmath82 having a value @xmath122 in the metastable initial phase , and a value @xmath123 in the new phase .
we start from a free energy functional @xmath124 & = \int\limits _ { v}\left(f_{s}+f_{dw}+f\right)dv,\label{eq : phasefieldmodel - free - energy - functional}\end{aligned}\ ] ] where @xmath125 is the gradient energy density and @xmath126 is the double well potential , guaranteeing that the free energy functional has two local minima at @xmath123 and @xmath122 corresponding to the two distinct phases of the system .
the form of the double well potential and the gradient energy density are chosen such that the phase field parameter @xmath120 defines the interface width and the parameter @xmath26 corresponds to the interface energy of the sharp interface description @xcite . the free energy density @xmath127 of our phase field description is chosen to be an interpolation between the two distinct free energies of the model , @xmath128}_{el}h(\phi)+f^{[\beta]}_{el}(1-h(\phi ) ) + l\frac{t - t_{eq}}{t_{eq } } ( 1-h(\phi)),\label{eq : phasefieldmodel - free - energy - interpolation}\end{aligned}\ ] ] where we choose the interpolation function @xmath129 to be @xmath130 .
it is the simplest polynom satisfying the necessary interpolation conditions @xmath131 and @xmath132 , and having a vanishing slope at @xmath123 and @xmath122 , in order not to shift the bulk states .
the evolution equation of the phase field is given by the variational expression @xmath133 for large values of the mobility @xmath134 we recover the desired case of diffusion limited growth .
the elastodynamic equations read @xmath135 which recover static elasticity for slowly propagating interfaces in comparison to the sound speed @xmath136 , with @xmath137 being the mass density .
for the temperature field we have the usual thermal diffusion equation with the motion of the phase field or interface being a source of latent heat , @xmath138 where the prime denotes the derivative with respect to @xmath82 .
the phase field model presented here is very similar to the model in @xcite .
however , there only the very special case of melting or crystallization of a confined sphere with elastic effects was discussed .    for the simulations we typically use parameters @xmath139 and @xmath140 .
apart from the theoretical investigations above , we also performed numerical simulations for both the free growth in an infinite system and in a channel or narrow strip .
as already mentioned before , we observe the growth of a single new phase if the eigenstrain tensor is purely diagonal , and twin structures for shear transformations with @xmath118 . also , we consider mixed mode situations which are between these two extreme cases .
we first consider the situation of free growth . here
, the simplest case is that of a dilatational eigenstrain where all diagonal elements have the same value , i.e.  @xmath141 ; this case is obtained from eq .
( [ mixedmodeeigenstraindefinition ] ) by setting @xmath142 . as already mentioned before , selection of a steady solution in the absence of elastic effects
is not possible for isotropic surface tension .
it is quite remarkable that also the inclusion of elastic effects due to the dilatational eigenstrain of the transition does not change this situation ; the reason is that the elastic effects can be mapped to a change of the undercooling only , and therefore we recover again the missing selection of conventional free dendritic growth @xcite .
it turns out that also for the volume - preserving shear strain transitions with @xmath143 , where all other components @xmath144 vanish , i.e. @xmath145 and @xmath53 , steady state solution do not exist either .
( see eq .
( [ eq : stabilityparameterdefinition ] ) ) as function of the mixing parameter @xmath146 for the privileged growth direction with an eigenstrain tensor @xmath147 when @xmath148 . on the right panel
we show the stability parameter @xmath149 for two different values of mixing @xmath146.,title="fig:",scaledwidth=45.0% ]   ( see eq .
( [ eq : stabilityparameterdefinition ] ) ) as function of the mixing parameter @xmath146 for the privileged growth direction with an eigenstrain tensor @xmath147 when @xmath148 . on the right panel
we show the stability parameter @xmath149 for two different values of mixing @xmath146.,title="fig:",scaledwidth=45.0% ]    under these circumstances it is very suprising , that for values of the mixing parameter in between zero and one @xmath150 , @xmath53 we do find steady state solutions for free growth and negative mixing .
this is shown in fig .
[ mixedmoderesultsinglecrystal2]a , where the selected `` stability parameter '' @xmath151 , which is a dimensionless measure for the growth velocity , is plotted as function of the mixing parameter @xmath146 . in general , the eigenvalue @xmath152 depends on @xmath153 and @xmath154 , but in the regime of small peclet numbers the elastic driving force and @xmath154 appear only in the combination @xmath155 .
[ mixedmoderesultsinglecrystal2]b shows the dependence of the selected stability parameter on @xmath155 for two different values of @xmath146 . for the other type of mixing @xmath156 corresponding to the unprivileged direction of growth , we do not find steady state solutions in free space .
we note that from the selected stability parameter @xmath157 one can extract the steady state growth velocity @xmath76 via the relation @xmath158    , where we choose a dilatational eigenstrain @xmath159 . for small driving forces
we obtain velocities and shapes corresponding to a symmetric solution ( see fig .  [
dilatationalsetup]@xmath160 ) . for higher driving forces
a first order transition towards asymmetrical solutions takes place ( see parity broken finger morphology in fig .  [
dilatationalsetup]@xmath161 ) . in the right panel
the comparison of growth velocities between phase field and the free space calculations is shown , for the case of negative mixing corresponding to the privileged growth direction @xmath162 .
the agreement is reasonably good , especially for driving forces @xmath100 slightly above the critical phase fraction @xmath163 of zero growth velocity.,title="fig:",scaledwidth=45.0% ] , where we choose a dilatational eigenstrain @xmath159 .
for small driving forces we obtain velocities and shapes corresponding to a symmetric solution ( see fig .  [
dilatationalsetup]@xmath160 ) . for higher driving forces
a first order transition towards asymmetrical solutions takes place ( see parity broken finger morphology in fig .  [ dilatationalsetup]@xmath161 ) . in the right panel
the comparison of growth velocities between phase field and the free space calculations is shown , for the case of negative mixing corresponding to the privileged growth direction @xmath162 .
the agreement is reasonably good , especially for driving forces @xmath100 slightly above the critical phase fraction @xmath163 of zero growth velocity.,title="fig:",scaledwidth=45.0% ]    in contrast , for the growth in a channel , steady state solutions do also exist for the case of a dilatational eigenstrain , @xmath87 ( see fig .
[ fig : diagonaleigenstrain]a ) . as expressed through eqs .
( [ eq : free - energy - gain ] ) and ( [ eq : delta - equilibrium ] ) , the phase @xmath0 can not grow for phase fractions @xmath100 below the threshold @xmath163 . beyond this value
the growth of a symmetrical finger ( dendrite ) is the most favorable solution .
we note that in situations where free growth solutions do not exist , the still existing solutions in the channel must have a scale that is determined by the size @xmath84 of the channel . for higher driving forces
an asymmetric finger grows with higher velocities , and is therefore more likely to be observed ; typical finger shapes are shown in fig .
[ dilatationalsetup ] .
for the situation of @xmath164 , where steady state solutions exist also in the free space , we performed a quantitative comparison of the results of the two methods in fig .
[ fig : diagonaleigenstrain]b . for this comparison
we directly relate the fraction @xmath165 to the driving force @xmath69 for free growth via eq .
( [ eq : fraction - of - beta - phase ] ) and eq .
( [ eq : shifted - undercooling ] ) .
then , in the limit of small peclet numbers , the latter is related to @xmath69 according to @xmath166 . finally , via the calculated function @xmath149 and eq .
( [ eq : stabilityparameterdefinition ] ) we obtain the growth velocity as a function of @xmath100 , which is then compared to the phase field results .
the results for the growth velocity agree well in a range of intermediate fractions of the new phase @xmath0 .
in particular , we confirmed that the steady state solutions which were obtained by the green s function method are indeed stable , as the dynamical phase field method tracks only these stable branches . for low phase fractions , growth is not possible below the threshold @xmath163 , which tends to zero for infinitely wide channels .
for very high driving forces , the growing finger is influenced by the geometrical confinement of the channel width and also affected by kinetic corrections to the sharp interface equations of motion , which therefore lead to a deviation from the free growth result which was obtained by the green s function method .
it is interesting to compare the situation to a case without elasticity , where the problem reduces to classical growth in a channel @xcite .
it is known that there even in wide channels symmetric solutions do not exist below a fraction @xmath167 , whereas we found that the eigenstrain leads to the existence of solutions at already significantly lower driving forces @xmath100 , see fig .
[ fig : diagonaleigenstrain]a .       shows the single tip solution and @xmath86 the twinned finger morphology , which corresponds to the most favorable solution branch for higher driving forces .
the positive mixing chosen for @xmath85 and @xmath86 leads to growth along the privileged direction and the parameters are @xmath168 .
the panel @xmath90 shows the shape of the slow growth direction for @xmath169 .
[ bicrystalshapes],scaledwidth=65.0% ]    in contrast to the cases studied above , the presence of off - diagonal eigenstrain entries lead to the preferential growth of bicrystal patterns .
this behavior can easily be understood intuitively , since the eigenstrain in the @xmath0 phase favors the presence of off - diagonal strain .
then the boundary conditions of fixed displacements require that the surrounding @xmath1 phase is sheared in the opposite direction .
this implies a finite elastic energy which is stored far behind the tip .
if , on the other hand , a second finger @xmath117 of the same size with an off - diagonal eigenstrain with opposite sign appears , the elastic energy contribution that is associated with the off - diagonal elements of the strain tensor vanishes .
consequently , the total elastic energy in the tail region is reduced , and the growth of the bicrystal therefore thermodynamically favored .
an additional energy contribution emerges from this bicrystal configuration , since the boundary between @xmath0 and @xmath117 generates a grain boundary energy @xmath170 .
we consider in particular the case of vanishingly small twin boundary energy , @xmath171 .
this affects also the geometrical situation at the tip , which is now a triple junction , and the contact angle @xmath82 becomes zero , so the overall finger has a smooth contour .
typical shapes of such bicrystal pattern are shown in fig .
[ bicrystalshapes ] , which are here related to the growth in a channel , but they look similar for free growth . we note that we simulated only half of the channel using the symmetry of the growing structures @xcite .    :
comparison of growth velocities obtained from the boundary integral technique for the free space geometry @xcite with the velocities from phase field simulations for the channel geometry ( @xmath172 ) .
[ fig : bicrystalpureshearcomp],scaledwidth=60.0% ]    considering the free space geometry , the first important observation is that in contrast to the single crystal scenario , steady state solution do exist also for pure shear , @xmath145 , which is already discussed in @xcite .
new is here that we found the same morphology also in the channel geometry , showing surprisingly good agreement in the comparison of growth velocities , done in the same way as was done for the single crystal described above ( see fig .  [
fig : bicrystalpureshearcomp ] ) .
again , we would like to point out that the stability of the bicrystal solution obtained from the green s function technique is confirmed by the phase field simulations .
we also did investigations of the bicrystal free growth considering mixed eigenstrains , and the results are summarized in fig .
[ bicrystalsigmaofeta ] and [ bicrystalpositiveandnegativedilatinshearsigmadeltael ] . in the first graph ,
the stability parameter @xmath157 , which is again a measure for the dimensionless velocity , is shown as function of the admixture parameter @xmath146 for fixed ratio @xmath155 .
as before , the elastic driving force and the @xmath154 appear only the combination @xmath155 for small peclet numbers .
as function of the admixture parameter @xmath146 , for the two different types of mixing and @xmath173 .
the solid line corresponds to the privileged direction of growth , i.e.  positive mixing , while the dotted line shows the dependence for the unprivileged negative mixing configuration.,scaledwidth=60.0% ]     as function of the control parameter @xmath155 , for positive mixing , whereas the in the right panel @xmath149 is plotted for negative mixing .
additionally , both graphs contain the @xmath149 dependence for the pure shear configuration @xmath145 , as a reference to compare both types of mixings , which was taken from @xcite .
, title="fig:",scaledwidth=45.0% ]   as function of the control parameter @xmath155 , for positive mixing , whereas the in the right panel @xmath149 is plotted for negative mixing .
additionally , both graphs contain the @xmath149 dependence for the pure shear configuration @xmath145 , as a reference to compare both types of mixings , which was taken from @xcite .
, title="fig:",scaledwidth=45.0% ]    as already mentioned above it is quite interesting that the appearance of an additional dilatational strain leads to an effective anisotropy of the system : whereas in the case of a positive admixture @xmath174 the velocity grows with increasing @xmath146 , the opposite happens for negative admixture @xmath175 .
this suggests that in free growth the bicrystal structure will eventually grow in the direction with the higher propagation velocity . for the fast branch
, we did not find solutions beyond the point @xmath176 .    in fig .
[ fig : bicrystallmixedmodecompchannelandfreegrowth ] we show the comparison of the two geometries for the bicrystal case .
again , the phase field results confirm the stability of the solution branch for low driving forces . in this region
, we obtain a convincing quantitative agreement between the growth velocities obtained for the free growth scenario using the green s function method , and the channel geometry with the phase field method . in the left panel
the favored direction is shown ( positive admixture ) for @xmath177 . in this case
we again observe dendritic solutions for small driving forces @xmath100 ( see also fig . [
bicrystalshapes]a ) . for the chosen parameters in the channel geometry , all solutions on the upper branch
are dendritic . on the lower branch , however ,
the steady state shapes are double - fingers , which means here the curvature directly at the point of symmetry has changed its sign and we obtain two tips ( see fig . [ bicrystalshapes]b ) .
the velocity of these branched solutions is here always smaller than for the single - tip structures . for the other direction , i.e. negative admixture , and a mixing parameter @xmath91
, we show the comparison of growth velocities in the right panel of fig .
[ fig : bicrystallmixedmodecompchannelandfreegrowth ] . here , we obtain a convincing agreement between the different methods and associated geometries at least for small driving forces , indicating dynamical stability also for this type of solution . obviously , the growth velocities are here substantially lower than for the previous case , and for fully free growth an initial bicrystal seed would therefore preferentially grow in the  fast " direction .     for bicrystal growth with mixed eigenstrain : comparison of free growth results obtained by green s function methods to phase field simulations for the growth in a channel . left panel : positive admixture with @xmath178 .
right panel : negative admixture , @xmath164.,title="fig:",scaledwidth=45.0% ]   for bicrystal growth with mixed eigenstrain : comparison of free growth results obtained by green s function methods to phase field simulations for the growth in a channel . left panel : positive admixture with @xmath178 .
right panel : negative admixture , @xmath164.,title="fig:",scaledwidth=45.0% ]
in summary , we investigated the diffusion - limited kinetics of solid - solid transformations in the presence of lattice strain .
we studied in particular the free growth with green s function methods and the growth in a channel by phase field techniques , exhibiting a comparable behavior .
it turns out that the elastic effects have a strong influence on the selection of the arising patterns and the corresponding steady state velocities , which differ crucially from conventional dendritic growth .
in particular , we find that the selected stability parameter and the corresponding growth velocity can be much higher compared to the classical case of dendritic growth , where the solutions are selected by the tiny effect of the anisotropy of surface tension .
different structural transitions , which are described by dilatational and shear eigenstrains , lead to a rich behavior , and we find very different structures already for the simple cases that were discussed here .      in comparison to the green s function method , which is computationally very efficient , the phase field technique has the advantage of higher flexibility . in particular , it can easily be extended to three - dimensional cases , it can be used also for non - steady state situations or anisotropic or non - symmetrical situations with different material parameters in the phases .