La pression comme multiplicateur de Lagrange

12

Dans les équations de Navier-Stokes incompressibles,

ρ(ut+(u)u)=p+μΔu+fu=0
le terme de pression est souvent appelé multiplicateur de Lagrange imposant la condition d'incompressibilité.

Dans quel sens est-ce vrai? Existe-t-il une formulation des équations de Navier-Stokes incompressibles en tant que problème d'optimisation soumis à la contrainte d'incompressibilité? Dans l'affirmative, existe-t-il un analogue numérique dans lequel les équations du débit de fluide incompressible sont résolues dans un cadre d'optimisation?

Ben
la source

Réponses:

18

μΔu+p=fu=0
minuμ2u2(f,u)so thatu=0.

Cette équivalence entre problèmes n'est exploitée dans aucun schéma numérique (à ma connaissance) mais c'est un outil important dans l'analyse car elle montre que les équations de Stokes sont essentiellement l'équation de Poisson sur un sous-espace linéaire. Il en va de même pour les équations de Stokes dépendant du temps (qui correspond à l'équation de chaleur sur le sous-espace) et il peut être étendu aux équations de Navier-Stokes.

Wolfgang Bangerth
la source
Merci pour une excellente réponse. Savez-vous si cette formulation peut être étendue au cas dépendant du temps?
Ben
1
Oui, comme je le dis, cela conduit à une équation de chaleur sur le sous-espace des fonctions libres de divergence.
Wolfgang Bangerth
1
Désolé, j'aurais dû être plus clair. Existe-t-il un moyen de refondre les équations de Stokes dépendantes du temps (ou Navier-Stokes) en tant que problème d'optimisation, éventuellement d'une fonction intégrée dans le temps?
Ben
1
Pas comme un problème d'optimisation - la solution de l'équation de la chaleur ne minimise rien (bien que ce soit le point stationnaire d'une fonction lagrangienne). Mais vous pouvez formuler les équations de Stokes comme suit: Trouvez sorte que pour tous les sous la contrainte que . Notez que j'ai choisi l'espace d'essai plus petit que l'espace d'essai et donc les côtés gauche et droit de l'équation variationnelle ne seront pas égaux. La différence est la pression. uHdiv(ut,φ)+(u,φ)=(f,φ)φ{vHdiv:v=0}u=0
Wolfgang Bangerth