Schémas spatiaux hybrides pour les CFD: tout inconvénient du mélange par rapport à la commutation?

9

Mis à part le coût de calcul supplémentaire dû au fait de devoir calculer les deux flux sur une certaine région, y a-t-il un inconvénient à mélanger deux évaluations de flux pour un schéma hybride dans une méthode de volume fini? L'évaluation du flux ressemblerait à ceci:

Fje+12=Λje+12Fje+12c+(1-Λje+12)Fje+12u

L'interrupteur est basé sur un capteur de gradient de pression et / ou de densité selon votre application. est un schéma central (McCormack, compact, ...) et F u est un schéma au vent comme un découpage de différence de flux avec une reconstruction MUSCL. Y a-t-il des problèmes en termes de propriétés numériques et conservatrices si je mélange les deux schémas en utilisant une fonction continue pour Λ au lieu de simplement basculer entre les schémas avec Λ valant 0 ou 1?FcFuΛΛ

FrançaisKheldar
la source
L'exemple classique des interrupteurs qui ont mal tourné est le schéma de Murman-Cole de 1970 pour les équations du potentiel transsonique (fixé en 1973). Si votre interrupteur ou mélange ne donne pas de somme télescopique, vous violez la conservation. C'est simple à tester.
tpg2114
Comment allez-vous résoudre les équations? Discrétisation (pseudo-) temporelle implicite? La commutation pourrait alors avoir un effet négatif sur la convergence itérative.
chris
Ce serait utiliser un schéma de temps explicite, soit un prédicteur-correcteur, soit un Runge-Kutta de base.
FrenchKheldar

Réponses:

4

L'approche que vous utilisez maintiendra la conservation de toute façon. Il existe d'autres approches évidentes qui ne sont pas conservatrices et peuvent causer des problèmes.

Il est possible (et même probable) que vous perdiez un ordre de précision dans la région où vous basculez, si vous examinez l'erreur de troncature locale. Mais généralement, cette erreur est localisée de sorte que l'erreur globale est toujours de l'ordre attendu. Donc, selon mon expérience, vous verrez essentiellement le même comportement que vous utilisiez un commutateur dur ou une région de transition.

J'ai un manuscrit sur (plus ou moins) ce sujet: Analyse d'erreur des schémas explicites de Runge-Kutta partitionnés pour les lois de conservation .

Je serais très intéressé de savoir ce que vous voyez lorsque vous essayez les deux approches, si c'est différent de ce que je suggère.

David Ketcheson
la source
génial, j'espérais que vous seriez en mesure de commenter cela.
Aron Ahmadia
Existe-t-il un moyen théorique de vérifier, semblable aux sommes télescopiques?
tpg2114