---

Hallo alle zusammen,

ich bin gerade dabei eine neue Randbedingung für Epsilon am Eintritt (inlet) in OpenFOAM zu implementieren. Mein Problem dabei ist, dass ich für diese Randbedingung unteranderem den Geschwindigkeitsgradient haben muss.

Ich weiß dass ich das Geschwindigkeitsfeld des Patches durch folgenden Code bekomme:

const fvPatchField<vector>& U = patch().lookupPatchField<volVectorField, vector>("U");

Allerdings ist dies dann ein PatchField. Leider kann ich aber nicht einfach den Gradient durch fvc::grad(U) berechnen, diese Funktion ein volVectorField bzw. volScalarField haben will.

Gibt es eine Möglichkeit auf das gesamte Geschwindigkeitsfeld zuzugreifen, um dann fvc::grad(U) auszurechnen und dann den Wert auf dem gewünschten Patch zu bestimmen?

Gruß

Christian

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP