Codice Matlab Netz gestaffelt Harlow welsh

[marvelreveal]

Hallo, alle. Ich bin dabei, einen Matlab-Code für die Auflösung von inkompressiblen 2d Schiff stokes zu schreiben, mit einem Mesh gestaffelt zu dem harlow welsh in einer quadratischen Domäne. in meiner Anfrage benötigen wir einen bc zu dem Dirichlet, was ich glaube, die unord Wand, Funktion von x, also eine echte Funktion zu sein. Einfügen dieses bc erzeugt ein Zeilenfeld von nx Punkten.
das Problem ist, dass in der Netz gestaffelt u Matrix hat Größe (nx x ny-1).
Ich weiß daher nicht, wie zu lösen ist, weil die Größe meines bc mit der Matrix, in der es eingefügt werden sollte, unvereinbar ist. Hat jemand eine Idee? Dies ist, was mich erzeugt, ein Problem aus der dimensionalen Sicht, und in der Berechnung von diffusen Begriffen.

 

[Stan9411]

Nun, die Nachfrage ist extrem technisch, um in einem Forum platziert werden, wo wir hauptsächlich von mechanischen Design sprechen.
in dieser Hinsicht eine kleine Anmerkung: Nehmen Sie nicht an, dass, wenn Sie sich an jemanden wenden, er Ihren eigenen Vokabular teilt... "bc." sogar das Akronym eines Engländers... "Bedingung zur Kontur" kostet so viel?

Ich verstehe nichts über Fluiddynamik. Aber als ich an der Universität war, habe ich einige dynamische Modelle zu Teildriften gemacht und ich war immer noch leidenschaftlich darüber.
Und meiner Meinung nach ist Ihr Problem nicht fluidodynamische, sondern mathematisch.

Insbesondere denke ich, dass Sie das Konzept der Schmutzgrenze falsch interpretieren.
in der Praxis, wenn Sie eine Differentialgleichung lösen, erhalten Sie keine Zahlen als Lösung, sondern eine Gleichung, ein Verhalten von etwas innerhalb einer Domain.
die Schmutzbedingung verhängt, was der Wert Ihrer Funktion am Rand der Domain ist.

wenn Ihre Differentialgleichung auf einer 1d-Domain definiert ist:
= 0
der Zustand des Schmutzes wäre:
= b
x0 und xf sind die Extreme Ihrer Domain 1d (z.B.: Zeit)

dass Ihre Funktion auf einer Oberfläche definiert ist, ist es blinkend, dass die Kontur Ihrer Domäne kein Paar von Punkten ist, aber es ist eine geschlossene Linie (z.B.: die Domäne ist ein Kreis, die Kontur ist ihr Außenumfang).

so ist es richtig, dass Ihr Zustand ein 1xn-Array erzeugt. wenn die Domain 2d ist, ist die Randbedingung 1d. Punkt.
Was nicht zurückkommt, ist, dass du es nur auf einer Seite deines Platzes hast. theoretisch sollten Sie es für alle vier Seiten haben. aber vielleicht ist das eine Besonderheit des fluiddynamischen Problems und hier kann ich Ihnen nicht helfen.