Codice Matlab maille décalée Harlow Welsh

[marvelreveal]

Bonjour à tous. Je suis sur le point d'écrire un code matlab pour la résolution d'incompressible 2d ship stokes, avec un maillage décalé vers le harlow Welsh dans un domaine carré. Dans ma demande, nous avons besoin d'un bc au dirichlet sur ce que je crois être le mur non-bord, fonction de x, donc une vraie fonction. l'insertion de ce bc génère un tableau de lignes de points nx.
le problème est que dans la matrice de maille décalée u a la taille (nx x ny-1).
Je ne sais donc pas comment résoudre, car la taille de mon bc est incompatible avec la matrice dans laquelle il doit être inséré. Quelqu'un a une idée ? c'est ce qui me génère est un problème du point de vue dimensionnel, et dans le calcul des termes diffus.

 

[Stan9411]

Eh bien la demande est extrêmement technique à placer dans un forum où nous parlons principalement de conception mécanique.
à cet égard une petite note: ne présumez pas que lorsque vous vous adressez à quelqu'un il partage votre propre vocabulaire... "bc." même l'acronyme d'un Anglais... "condition au contour" coûte tellement cher ? :

Je ne comprends rien à la dynamique des fluides. Mais quand j'étais à l'université, j'ai fait plusieurs modèles dynamiques à des dérives partielles et j'étais toujours passionné par ça.
Et à mon avis, votre problème n'est pas fluidodynamique mais mathématique.

En particulier, je pense que vous interprétez mal le concept de l'état de la frontière du dirichlet.
en pratique, lorsque vous résolvez une équation différentielle, vous n'obtenez pas de nombres comme solution, mais une équation, un comportement de quelque chose dans un domaine.
la condition de dirichlet impose la valeur de votre fonction au bord du domaine.

si votre équation différentielle est définie sur un domaine 1d:
= 0
l'état du dirichlet serait:
= b
où x0 et xf sont les extrêmes de votre domaine 1d (par exemple : temps)

que votre fonction est définie sur une surface, il est clignotant que le contour de votre domaine n'est plus une paire de points, mais une ligne fermée (par exemple: le domaine est un cercle, le contour est sa circonférence extérieure).

Il est donc juste que votre état génère un tableau 1xn. si le domaine est 2d, la condition limite est 1d. point.
Ce qui ne revient pas, c'est que tu ne l'as que d'un côté de ta place. En théorie, vous devriez l'avoir pour les quatre côtés. Mais c'est peut-être une particularité du problème des fluides et je ne peux pas vous aider.