ボルツマン法、数理最適化、C++コード付き)とか勉強したことをまとめます。右のカテゴリーから興味のある記事を探してください。最近はクラシックの名演も紹介しています。noteにも書いています。Galerkin)法で解いていきます。今までずっとやろうと思いながらも放置してきましたが、やっとやりました。1次元移流方程式とはのような一階の偏微分方程式のことをいいます。風などによる移流速度です。SUPG法とは?SUPG法とは、有限要素法における風上化の一種です。Galerkin法とは、形状関数とは異なる重み関数を用いる有限要素法のことを言います。風上方向に重み関数を歪ませたのがSUPG法です。移流方程式にたいして通常のGalerkin法を適用すると数値振動が発生することから開発されました。簡単に離散化を説明します(詳細な離散化は現在執筆中です)。行列で書くととなります。がにたいするSUPG項です。それぞれ、のような形をしています。SUPG項は移流行列と拡散行列に適当な係数をかけた形になっています。キン法で一度離散化していれば楽に計算できます。あとは時間項に対して適当な時間進行を入れればOKです。今回はオイラー法を使います(陰解法にすると数値振動を抑えられます)。初期条件は領域の左端の矩形波(長方形の波)です。