空間積分

アンチアワグラスと接触力のある並進速度に対する全体の運動方程式は次のように書き表されます:(1)
M v t = F ext F int + F bod + F hgr + F cont + F trm MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWHnbWaaS aaaeaacqGHciITcaWG2baabaGaeyOaIyRaamiDaaaacqGH9aqpcaWH gbWaaSbaaSqaaiaadwgacaWG4bGaamiDaaqabaGccqGHsislcaWHgb WaaSbaaSqaaiGacMgacaGGUbGaaiiDaaqabaGccqGHRaWkcaWHgbWa aSbaaSqaaiaadkgacaWGVbGaamizaaqabaGccqGHRaWkcaWHgbWaaS baaSqaaiaadIgacaWGNbGaamOCaaqabaGccqGHRaWkcaWHgbWaaSba aSqaaiaadogacaWGVbGaamOBaiaadshaaeqaaOGaey4kaSIaaCOram aaBaaaleaacaWG0bGaamOCaiaad2gaaeqaaaaa@598B@
ここで、
F t r m MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWHgbWaaS baaSqaaiaadshacaWGYbGaamyBaaqabaaaaa@3A3C@
次のように与えられる輸送運動量ベクトル
(2)
F t r m = e l e m e n t s f t r m MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWHgbWaaS baaSqaaiaadshacaWGYbGaamyBaaqabaGccqGH9aqpdaaeqbqaaiaa hAgadaahaaWcbeqaaiaadshacaWGYbGaamyBaaaaaeaacaWGLbGaam iBaiaadwgacaWGTbGaamyzaiaad6gacaWG0bGaam4Caaqab0Gaeyye Iuoaaaa@48E7@
(3)
f iI trm =( 1+ η I ) v ρ Φ I ( w j v j ) v i x j dV MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWHMbWaa0 baaSqaaiaadMgacaWGjbaabaGaamiDaiaadkhacaWGTbaaaOGaeyyp a0ZaaeWaaeaacaaIXaGaey4kaSIaeq4TdG2aaSbaaSqaaiaadMeaae qaaaGccaGLOaGaayzkaaGaeyyXIC9aa8quaeaacqaHbpGCcqGHflY1 cqqHMoGrdaWgaaWcbaGaamysaaqabaGcdaqadaqaaiaadEhadaWgaa WcbaGaamOAaaqabaGccqGHsislcaWG2bWaaSbaaSqaaiaadQgaaeqa aaGccaGLOaGaayzkaaGaeyyXIC9aaSaaaeaacqGHciITcaWG2bWaaS baaSqaaiaadMgaaeqaaaGcbaGaeyOaIyRaamiEamaaBaaaleaacaWG QbaabeaaaaGccaWGKbGaamOvaaWcbaGaamODaaqab0Gaey4kIipaaa a@60F8@

上の式で定義された全てのマトリックスとベクトルは空間領域で積分されます。質量マトリックスは密度と領域が時間で変化するため一定ではありません。Peclet数として知られたメッシュパラメータが限界値を超えると、方程式の解は振動することが示されています。この空間の安定性は数値的な拡散をスキームに付加することで回避可能で、そうでない場合は一般的に拡散が足りず、このように不安定になります。材料入力に運動量風上係数が定義されます。完全な風上の場合 (coef=1)がRadiossのデフォルト値で、一般的に用いられます。

1つの風上テクニックがRadiossで可能です:
  • Petrov-Galerkin法による流線風上(SUPG):
    SUPGは形状関数の運動量の対流項を考慮した修正から成ります。バージョン2018より、SUPG定式化はデフォルトで有効とされています。SUPGは、Engineファイルオプションによりオフにすることができます:
    /UPWM/SUPG