可変速度の移流方程式は保守的ですか?


13

速度係数が可変の移流方程式をもう少し良く理解しようとしています。特に、方程式がどのように保守的であるかはわかりません。

移流方程式、

∂u∂t+∂∂x(vu)=0

u(x,t)を、物理種(cm−3)の濃度または作成または破壊できない物理量として解釈してみましょう。ドメイン上でu(x、t)を統合する場合、u(x,t)定数を取得する必要があります。

∫xminxmaxu(x,t)dx=constant

(これは保守的であることの意味です。)

ここで、速度を空間(および時間)の関数v(x,t)にすると、チェーンルールを適用して、

∂u∂t+v∂u∂x+u∂v∂x⏟?=0

最後の用語はソース用語のように「見え」ますが、これは紛らわしいものです。速度場の発散に応じて、量u増減します。

この質問に続いて、保存境界条件を課す方法を理解しています。ただし、可変速度移流方程式については、チェーンルールを適用することによって導入される追加の「ソース項」のために、保存境界条件をどのように導出できるかがわかりません。この方程式は保守的ですか?もしそうなら、どのように正しい境界条件を適用できますか?

回答:


15

輸送の基本的な量であるフラックス、移流のため。発散定理は、vu

∫Ω∇⋅(vu)=∫∂Ω(vu)⋅n.

この等式が保存されている場合、方程式は保守的です。 1Dにドロップし、方程式、Ω=(a,b)ut+(vu)x=0

(∫abu)t=∫abut=−∫ab(vu)x=−vu|ab

ここで、右側の用語は、左右の境界間のフラックスの差にすぎません。

2番目の観察に関して、非保守的(非発散)の形式は誤解を招きます(そして、スムーズなソリューションに対してのみ正当化されます)。製品あるない場合は保守的な輸送発散フリー(1Dにおけるすなわち、定数)ではありません。保存的な形式を固守し、保存プロパティを評価するときにチェーンルールを適用する衝動に抵抗する必要があります。v⋅∇uv


本当に明確な答えをありがとう、ジェド!これについてフォローアップの質問をすることになると思いますが、最初に提案を実行する必要があります。
— boyfarrell
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.