Hallo Lucas,
da ich mir auch nicht so ganz sicher war, musste ich bei einem OF-Meister mal nachfragen Ich möchte dir an dieser Stelle hier seine Meinung zitieren:
Zitat:
UEqn
(
fvm::Div(phi, U)
- fvm::Sp(fvc::Div(phi), U)
+ turbulence->divDevRhoReff(U)
) == -fvc::grad(p) Es gilt allgemein div(phi, U) = phi & grad U + U div(phi) und folglich phi & grad U = div(phi, U) - U div(phi) (also was auf der linken Seite der UEqn steht). Dies könnte man als schwache Form der Erhaltungsgleichung bezeichnen, und man kann sagen, im stationären Fall muss der Term div(phi) und damit - fvm::Sp(fvc::Div(phi), U) soundso Null sein. Ich addiere folglich lediglich eine Null hinzu (Massenerhaltung).
Im instationären, kompressiblen Fall ist der Term nicht a priori null, weil er nur zusammen mit der Zeitableitung Null wird. Deswegen steht er evtl. nicht in rhoPimpleFoam.
Der Term bewirkt somit nur etwas im nicht auskonvergierten Zustand.
Es scheint, es handelt sich hier bei also um einen numerischen Trick.
Ich hoffe, du kannst damit etwas anfangen? Wenn du es genau weißt, würden ich mich freuen, wenn du Info gibst.
Bis dann
Thomas
Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP