Comment les coefficients non constants doivent-ils être traités avec un schéma au vent de premier ordre à volume fini?

11

En commençant par l'équation d'advection sous forme de conservation.

u_{t} = (une (X) u)_{X}

$u_t = (a(x)u)_x$

où $a(x)$ est une vitesse qui dépend de l'espace, et $u$ est une concentration d'une espèce qui est conservée.

La discrétisation du flux (où le flux $f=a(x)u$ , est défini sur les bords des cellules entre les points de maillage) donne,

u_{t} = \frac{1}{h} (F_{j - \frac{1}{2}} - F_{j + \frac{1}{2}})

$u_t = \frac{1}{h}\left( f_{j-{\frac{1}{2}}} - f_{j+{\frac{1}{2}}} \right)$

En utilisant un vent de premier ordre, nous approchons les flux comme,

F_{j - \frac{1}{2}} = une (X_{j - \frac{1}{2}}) u_{j - 1} F_{j + \frac{1}{2}} = une (X_{j + \frac{1}{2}}) u_{j}

$f_{j-{\frac{1}{2}}} = a(x_{j-\frac{1}{2}})u_{j-1} \\ f_{j+{\frac{1}{2}}} = a(x_{j+\frac{1}{2}})u_{j}$ Ce qui donne,

u_{t} = \frac{1}{h} (une (X_{j - \frac{1}{2}}) u_{j - 1} - une (X_{j + \frac{1}{2}}) u_{j})

$u_t = \frac{1}{h} \left( a(x_{j-\frac{1}{2}})u_{j-1} - a(x_{j+\frac{1}{2}})u_{j} \right)$

Si $a(x)$ était constant, cela se ramènera au schéma au vent familier, c'est-à-dire $u_t = \frac{a}{h} \left( u_{j-1} - u_{j} \right)$ .

Ma question est, comment pouvons-nous traiter les coefficients non constants de l'équation d'advection? La vitesse est définie au centre des cellules, donc une approche simple serait la suivante,

une (X_{j - \frac{1}{2}}) \to une (X_{j - 1}) une (X_{j + \frac{1}{2}}) \to une (X_{j})

$a(x_{j-\frac{1}{2}}) \rightarrow a(x_{j-1}) \\ a(x_{j+\frac{1}{2}}) \rightarrow a(x_{j})$

C'est mon approche préférée car elle est très simple à mettre en œuvre.

Cependant, nous pourrions également utiliser (je suppose) un schéma de moyenne pour définir la vitesse aux bords des cellules,

une (X_{j - \frac{1}{2}}) \to \frac{1}{2} une (X_{j - 1}) + \frac{1}{2} une (X_{j}) une (X_{j + \frac{1}{2}}) \to \frac{1}{2} une (X_{j}) + \frac{1}{2} une (X_{j + 1})

$a(x_{j-\frac{1}{2}}) \rightarrow \frac{1}{2}a(x_{j-1} ) + \frac{1}{2}a(x_{j} ) \\ a(x_{j+\frac{1}{2}}) \rightarrow \frac{1}{2}a(x_{j} ) + \frac{1}{2}a(x_{j+1} ) \\$

Dans le livre de LeVeque, il dit:

Jusqu'à présent, nous avons supposé que la vitesse variable est spécifiée par une valeur constante dans la j-ième cellule de la grille. Dans certains cas, il est plus naturel de supposer à la place qu'une vitesse $a(x)$ $a_j$ est spécifié à chaque interface de cellule. $a_{j−{\frac{1}{2}}}$

Mais il n'élabore pas vraiment trop après cela. Qu'est-ce qu'une approche commune?

Je résout un problème de conservation (j'utilise l'équation d'advection comme équation de continuité), donc je veux m'assurer qu'après avoir appliqué la discrétisation, la propriété de conservation est préservée. Je voudrais éviter toute surprise cachée concernant ces coefficients variables! Quelqu'un a-t-il des commentaires et des conseils généraux?

Mise à jour Il y a deux très bonnes réponses ci-dessous et je n'ai pu en choisir qu'une :(

discretization finite-volume numerical-analysis boyfarrell
la source

4

Selon le type de système que vous regardez, il peut être plus pratique de considérer la vitesse comme constante par morceaux dans chaque cellule, ou qu'elle est définie aux interfaces des cellules. Par exemple, en météorologie, les grilles décalées sont assez courantes, où la pression peut être définie à l'intérieur des cellules et la vitesse aux interfaces des cellules. Vous pourriez tout aussi bien penser à la vitesse définie dans les cellules. Tout compte fait: le choix de la représentation ne doit pas affecter la convergence de votre méthode *, à condition que votre discrétisation soit stable et cohérente. $a$

Ce qui est le plus important (et vous l'avez déjà évoqué dans votre question), c'est que le système discrétisé est toujours conservateur. A condition que votre schéma puisse être écrit sous la forme

$\frac{\partial u_j}{\partial t} = F_{j-\frac{1}{2}}(u_{j-1},u_j) - F_{j+\frac{1}{2}}(u_j,u_{j+1})$

alors il devrait être conservateur, car

. $\frac{\partial}{\partial t}\int u dx = \sum_j\frac{\partial u_j}{\partial t}\delta x = \sum_j( F_{j-\frac{1}{2}}- F_{j+\frac{1}{2}} )\delta x = (F_{-\frac{1}{2}}-F_{N+\frac{1}{2}})\delta x$

Votre approche simple devrait fonctionner correctement, tout comme la moyenne de la vitesse entre les cellules pour la définir sur les interfaces cellulaires, à condition que la vitesse soit toujours positive. De plus, je ne pense pas que la moyenne vous apportera une précision supérieure, vous avez donc raison de préférer la méthode simple.

Si vous résolvez également la vitesse et que vous avez un système d'équations, vous devrez peut-être être plus prudent. De même, si vous résolvez un PDE hyperbolique non linéaire et utilisez des limiteurs de flux, vous devez être encore plus prudent.

* Cependant, pour un système de PDE hyperboliques, l'utilisation de grilles décalées peut considérablement améliorer la dispersion / diffusion artificielle. Si vous voulez en savoir plus, consultez les grilles C d'Arakawa ou consultez le chapitre 4 de ce livre .

Daniel Shapero
la source

Merci d'avoir expliqué. Et votre intuition est correcte; Je résous un système d'équations où l'une des équations est la vitesse (un PDE des autres variables). Le système d'équations est uniquement 1D, je prévois d'utiliser une méthode adaptative au premier ordre au vent (peut basculer entre le deuxième ordre au centre et au vent), peut-être avec un ajustement exponentiel. Je n'utilise pas de limiteurs de flux, mais le système n'est pas linéaire. Dois-je être "plus prudent" dans cette situation?

boyfarrell

Tout dépend si vous vous attendez à ce que des ondes de choc et similaires se forment, s'il y a une possibilité que la vitesse passe en dessous de zéro dans certaines régions, ou si la vitesse peut devenir suffisamment élevée pour que vous exécutiez la condition Courant-Friedrichs-Lewy à un moment donné. Cela dit, j'essaierais d'abord l'approche simple pour voir si cela fonctionne, ce qui pourrait bien le faire. Si cela échoue, cela se produira de manière spectaculaire et sans ambiguïté, donc je ne pense pas que vous ayez à vous soucier d'avoir quelque chose de mal à glisser sous votre radar.

Daniel Shapero

Oui, je m'attends à ce que la vitesse soit uniquement non nulle uniquement au centre de mon domaine, puis j'approche rapidement de zéro lorsque l'on s'éloigne du centre. Je choisis le pas de temps pour que la condition CFL soit satisfaite (en utilisant la vitesse max.), Le maillage est fixe. Quels sont les critères d'une onde de choc? Je ne prévois pas voir ça (mais on ne sait jamais).

boyfarrell

5

$a(x)$

Ce que je veux dire par cohérent, c'est que la seule condition que l'interpolation doit satisfaire est

{une}_{je + 1 / 2}^{+} = {une}_{je + 1 / 2}^{-}

$a_{i+1/2}^+ = a_{i+1/2}^-$

En d'autres termes, tant que votre méthode d'interpolation est continue à travers les limites des cellules, votre discrétisation est garantie de rester conservatrice.

Cela peut ne pas sembler être un gros problème ici dans 1D (et cela ne devrait pas), mais peut provoquer des problèmes aux interfaces grossières et fines sur les grilles AMR à plusieurs niveaux.

GradGuy
la source

u_{j + \frac{1}{2}}

$u_{j+\frac{1}{2}}$

a (x_{j + \frac{1}{2}})

$a(x_{j+\frac{1}{2}})$

a (x_{j + \frac{1}{2}}) \to a (x_{j + 1})

$a(x_{j+\frac{1}{2}}) \rightarrow a(x_{j+1})$

u_{j + \frac{1}{2}}

$u_{j+\frac{1}{2}}$

@boyfarrell Ce serait bien dans le sens où la méthode reste conservatrice. Cela affecte cependant la précision de la solution. Souvent, par exemple dans les schémas ENO, on se rapproche de l'ensemble de la fonction de flux et non de la vitesse et de la solution séparément.

GradGuy

4

$a(x_{j-\frac{1}{2}})$

Pour voir pourquoi il en est ainsi, considérons que la définition analytique du conservateur est que

\frac{\partial}{\partial t} \int_{ré} u (X) ré X = \int_{\partial ré} une (X) u (X) ré S,

$\frac{\partial}{\partial t}\int_D u(x)\, dx = \int_{\partial D} a(x)u(x) dS,$

$D$

Si notre discrétisation est de la forme

u_{t} (X_{j}) = \frac{1}{h} (une (X_{j - \frac{1}{2}}) u_{j - \frac{1}{2}} - une (X_{j + \frac{1}{2}}) u_{j + \frac{1}{2}})

$u_t(x_j) = \frac{1}{h}\left(a(x_{j-\frac{1}{2}})u_{j-\frac{1}{2}} - a(x_{j+\frac{1}{2}})u_{j+\frac{1}{2}}\right)$

$x_1,\ldots, x_n$ $D = [c,d]$ $c = x_{\frac{1}{2}}$ $d = x_{n+\frac{1}{2}}$

\frac{1}{h} \sum_{j = 1}^{n} (une (X_{j - \frac{1}{2}}) u_{j - \frac{1}{2}} - une (X_{j + \frac{1}{2}}) u_{j + \frac{1}{2}}) = une (X_{\frac{1}{2}}) u_{\frac{1}{2}} - une (X_{n + \frac{1}{2}}) u_{n + \frac{1}{2}},

$\frac{1}{h}\sum_{j=1}^n\left(a(x_{j-\frac{1}{2}})u_{j-\frac{1}{2}} - a(x_{j+\frac{1}{2}})u_{j+\frac{1}{2}}\right) = a(x_{\frac{1}{2}})u_{\frac{1}{2}} - a(x_{n+\frac{1}{2}})u_{n+\frac{1}{2}},$

$u_{j-\frac{1}{2}} = u_{j-1}$ $u_{j+\frac{1}{2}} = u_j$ $a(x)u$

$a(x)$ $a(x_{j-r}), \ldots, a(x_{j+s})$ $a(x_{j-\frac{1}{2}})$

Ben
la source

Comment les coefficients non constants doivent-ils être traités avec un schéma au vent de premier ordre à volume fini?

Réponses: