Dans tous les cas considérés jusqu'à présent, les équations dynamiques donnent naissance à une structure cellulaire dont les cellules sont symétriques non seulement dans la direction x mais aussi dans la direction perpendiculaire z, celle de la direction vicinale. Cependant, les expériences de croissance menées sur des surfaces vicinales de
Cu(1, 1, 17) [#!Maroutian99!#] montrent une structure de méandre fortement asymétrique dans la direction vicinale.
Des simulations Monté-Carlo effectuées par Pierre-Louis et al. [#!Pierre-Louis99!#] pour des surfaces vicinales en croissance mènent à des résultats similaires : l'instabilité de méandre observée conduit à la formation d'une structure cellulaire asymétrique entre le devant des marches et leur partie arrière. La morphologie
n'est pas symétrique sous la transformation
z - z.
Les équations d'évolution dynamiques dérivées jusqu'à présent des équations constitutives du modèle microscopiques ont été tronquées à l'ordre dominant en , le paramètre caractérisant l'écart par rapport au seuil de l'instabilité.
Ce paramètre
est aussi directement proportionnel à la vitesse moyenne du train de marches.
Or, à l'ordre dominant, les équations de la dynamique ne font pas apparaître ce paramètre (l'équation adimensionnée de la dynamique ne dépend, à l'ordre dominant que du paramètre
). Comme mentionné précédemment, l'ordre dominant correspond formellement à la limite
0 ; pour prendre en compte l'asymétrie du système par rapport à la direction z, le développement de l'équation d'évolution dynamique doit alors être poursuivi aux ordres supérieurs.
Nous nous restreindrons par soucis de simplicité au cas unilatéral.
Le calcul menant à la détermination de la contribution sous-dominante de la dynamique est développé en annexe ; donnons en ici le résultat. L'équation d'évolution du méandre est donnée par :
Afin de prendre en compte les termes supplémentaires avec une bonne précision, le schéma numérique précédent à été modifié.
Deux méthodes ont été utilisées qui ont conduit aux mêmes résultats.
Dans un premier temps, le schéma précédent à été simplement modifié en prenant une « grille numérique » adaptative. En adaptant cette grille numérique aux variations géométriques du méandre on a ainsi une meilleure évaluation des dérivées spatiales.
Ces dérivations spatiales se font à l'aide d'une méthode d'éléments finis par interpolation sur des bases locales de polynômes (voir annexe ).
Cependant la méthode à grille adaptative, bien que plus précise devient plus coûteuse en temps de calcul. En effet, la réadaptation de la grille ( « remeshing » ) introduit de nouvelles perturbations qu'il faut laisser relaxer.
Nous avons donc dans un second temps reformulé complètement les équations dynamiques en utilisant une représentation en coordonnées « intrinsèques » , en terme d'angle polaire
et d'abscisse curviligne s comme représenté sur la figure [
].
Cette reformulation a l'avantage, outre d'une écriture simplifiée des équations dynamiques, de permettre une reparamétrisation dynamique du méandre (voir annexe
).
La dépendance temporelle de la paramétrisation permet un « choix de jauge » dans la formulation de la vitesse tangentielle [#!Langer92!#], [#!Csahok99!#].
En utilisant cette représentation « intrinsèque », l'équation de la dynamique se formule de la façon suivante pour la vitesse normale vn :
vn = - ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
(3.100) |
où, comme précédemment, le temps a été redimensionné par
4/F
et les variables spatiales x et
, par
. La dynamique ne dépend alors plus que de deux paramètres adimensionnés :
= DLa/DS
qui mesure l'importance de la diffusion de ligne relativement à la diffusion de terrasse et
qui mesure l'écart par rapport au seuil de l'instabilité.
|
|
L'intégration numérique de cette équation donne des résultats similaires à ceux trouvés à l'ordre dominant : la dynamique donne naissance à une structure cellulaire dont la longueur d'onde caractéristique correspond à celle du mode le plus instable linéairement.
De plus, à temps long, la rugosité du méandre suit la même loi d'échelle :
w t1/2.
Cependant, la forme des cellules formées diffère qualitativement de celle trouvée en ne retenant que les termes d'ordre dominant : les cellules sont cette fois asymétriques par rapport à la direction z (voir figures [
] et [
]).
La forme asymétrique des cellules obtenue lors de nos simulations numériques de l'équation dynamique développée jusqu'à l'ordre sous-dominant ressemble qualitativement bien à celle observée dans les expériences de croissance de surfaces vicinales de cuivre menées par T. Maroutian et al. ainsi qu'à celle obtenue par simulation Monté-Carlo [#!Pierre-Louis99!#].
|
Cependant, un problème demeure : la rugosité du méandre croît au court du temps comme t1/2 (cf. fig.[]) et rien ne semble saturer (du moins dans la présente équation d'évolution) sa croissance.
Ce comportement singulier de la rugosité semble robuste car a priori indépendant du processus de relaxation envisagé.
Cela implique que la distance entre marches tend vers zéro en certains endroits, là où la pente du méandre devient importante.
En effet la distance locale entre marches est :