Escape over a saddle by coloured noise: theory and numerics

Résumé Technique : Évasion par-dessus une selle par un bruit coloré : Théorie et Numérique

Énoncé du Problème
Les systèmes dynamiques stochastiques sont essentiels pour modéliser les transitions induites par des perturbations aléatoires, particulièrement les évasions depuis des équilibres attracteurs à travers la variété stable d'un point de selle. Bien que la limite de faible bruit soit théoriquement bien comprise via le principe de grande déviation de Freidlin-Wentzell (LDP), les méthodes de calcul existantes font face à des limitations significatives. Les approches actuelles peinent généralement avec :

1. Le bruit dégénéré : Des situations où le bruit n'agit pas dans toutes les directions de l'espace des phases (par exemple, le bruit n'entrant que dans les équations de quantité de mouvement dans les systèmes mécaniques).
2. Le bruit coloré : Un forçage réaliste avec des temps de corrélation finis, souvent modélisé via des filtres linéaires sur un bruit blanc, plutôt que par un bruit blanc idéalisé.
3. Les horizons temporels non bornés : La réalité physique selon laquelle les événements d'évasion se produisent sur un temps infini, alors que les méthodes numériques reposent souvent sur des intervalles finis, ce qui entraîne des dégénérescences de translation temporelle.

Les méthodes existantes pour calculer les hauteurs de barrière sur des horizons temporels infinis dépendent généralement de l'invertibilité de la matrice de covariance du bruit, ce qui les rend inapplicables aux cas de bruit dégénéré. De plus, l'échantillonnage direct d'événements rares est prohibitivement coûteux en calcul, et la discrétisation standard des trajectoires à temps infini gaspille des ressources dans les régions où la dynamique est presque linéaire.

Méthodologie
Les auteurs proposent un cadre unifié combinant trois innovations méthodologiques distinctes pour relever ces défis :

1. Contrôle Optimal Hamiltonien (HOC) pour le Bruit Dégénéré et Coloré :

   • Pour gérer le b bruit dégénéré (où la matrice de covariance $C = \sigma\sigma^T$ est singulière), l'article adopte la méthode de contrôle optimal hamiltonien. Au lieu de minimiser la fonctionnelle d'action impliquant $C^{-1}$ (qui est mal définie pour le bruit dégénéré), le problème est reformulé comme une tâche variationnelle déterministe dans un espace des phases étendu.
   • En introduisant des variables conjuguées (adjointes) pour chaque composante d'état, la méthode dérive les équations du mouvement de Hamilton. Cette approche évite totalement l'inversion de la matrice de covariance du bruit, accommodant naturellement la dégénérescence sans nécessiter de pré-conditionnement explicite ou d'hypothèses strictes d'hypoellipticité.
   • Le bruit coloré est modélisé en augmentant l'espace des états avec des variables de filtre régies par des équations différentielles stochastiques linéaires (processus d'Ornstein-Uhlenbeck). Cela transforme le problème du bruit coloré en un cadre d'SDE standard avec un vecteur d'état de dimension supérieure, où le bruit n'agit que sur les variables de filtre, résultant automatiquement en une structure de bruit dégénéré pour le système d'origine.

2. La Méthode de Division (MOD) pour les Horizons Temporels Infinis :

   • Pour remédier à l'infaisabilité computationnelle de la simulation d'intervalles de temps infinis et à l'instabilité numérique causée par la dégénérescence de translation temporelle, les auteurs introduisent la Méthode de Division.
   • La trajectoire de transition est décomposée en trois phases distinctes :
     • Évasion Initiale : Le système s'échappe de l'équilibre stable (puits) sur $(-\infty, 0]$. Cette phase est traitée analytiquement en utilisant la dynamique linéarisée autour du puits. L'action est calculée via le Gramien de contrôlabilité, résolvant une équation de Lyapunov.
     • Phase de Transition : La transition non linéaire centrale se produit sur un intervalle fini $[0, T]$. Cette partie est la seule résolue numériquement via la méthode HOC.
     • Approche Finale : Le système approche la selle sur $[T, \infty)$. Cette phase est également traitée analytiquement en linéarisant la dynamique autour de la selle. La direction instable est gérée explicitement, tandis que les directions stables convergent naturellement vers la selle.
   • Briser la Dégénérescence de Translation Temporelle : Pour fixer l'ambiguïté numérique où les trajectoires optimales peuvent glisser librement le long de l'axe temporel, la méthode impose une contrainte terminale. La trajectoire doit traverser un ellipsoïde transverse prédéfini autour de la selle au temps $T$. Cette contrainte brise la symétrie, assurant un problème numérique bien conditionné.

3. Implémentation Numérique :

   • Le problème d'optimisation global minimise la somme des actions des trois phases par rapport à la trajectoire de transition et aux points de terminaison optimaux ($\phi_0$ près du puits et $\phi_T$ sur l'ellipsoïde).
   • La méthode du Lagrangien Augmenté est employée pour résoudre le problème d'optimisation sous contrainte. Cette approche convexifie la fonction objectif et met à jour les multiplicateurs de Lagrange et les paramètres de pénalité de manière itérative pour satisfaire la contrainte de l'ellipsoïde terminal sans nécessiter de paramètres de pénalité infinis, évitant ainsi le mauvais conditionnement.
   • L'algorithme utilise des méthodes de gradient conjugué non linéaires pour une convergence efficace.

Contributions Clés

• Extension au Bruit Dégénéré : L'article fournit un cadre de calcul robuste pour calculer les taux de transition et les chemins les plus probables dans des systèmes avec bruit dégénéré, un scénario auparavant inaccessible aux méthodes reposant sur l'inversion de la matrice de covariance.
• Gestion du Bruit Coloré : En intégrant des filtres linéaires dans l'espace d'état et en utilisant HOC, la méthode gère efficacement le bruit coloré physiquement réaliste sans recourir à des approximations de bruit blanc.
• Solution d'Horizon Temporel Infini : La Méthode de Division offre une stratégie de calcul novatrice pour les événements rares sur des temps infinis. Elle isole la région de transition non linéaire coûteuse en calcul tout en résolvant analytiquement le comportement asymptotique près des équilibres, améliorant considérablement l'efficacité de calcul et la convergence.
• Suppression de la Dégénérescence de Translation Temporelle : L'introduction d'une contrainte terminale transverse fournit un mécanisme rigoureux pour stabiliser les solutions numériques des trajectoires à temps infini.

Résultats et Applications
L'efficacité du cadre proposé est démontrée à travers deux études de cas spécifiques :

1. Potentiel à Double Puits Inversé : Un modèle simplifié utilisé pour illustrer les capacités théoriques de la méthode pour gérer le bruit dégénéré et coloré.
2. Modèle de Chavirage de Navire : Un modèle simplifié de roulis-pilonnement représentant une application critique en ingénierie maritime. La méthode calcule avec succès les chemins d'évasion les plus probables et les lois d'échelle exponentielles des taux de transition pour un navire soumis à un forçage de vagues stochastiques, capturant la contribution exponentielle dominante au taux de chavirage.

Signification
Cet article positionne ce travail comme un pont entre l'analyse stochastique et la théorie des systèmes dynamiques non autonomes, répondant spécifiquement au défi du chavirage des navires. Il affirme offrir une approche physiquement réaliste et numériquement robuste pour l'ingénierie de la sécurité et de la fiabilité, où l'estimation du taux de transitions très improbables (telles que le chavirage) est cruciale. En permettant le calcul des hauteurs de barrière et des chemins les plus probables pour le bruit dégénéré et coloré sur des horizons temporels infinis, la méthode fournit une base rigoureuse pour les évaluations de sécurité à long terme qui étaient auparavant hors de portée des techniques de calcul existantes. Ce travail sert de compagnon stochastique aux approches déterministes précédentes dans la dynamique non autonome, intégrant ces domaines pour mieux comprendre les transitions d'état dans des systèmes physiques complexes.