lzpel/fluid-simulation

Hi there 👋

🔭 I’m currently working on ...

• 2022~ cloudmama：日時に合わせて適切なアドバイスをスピーカーから再生する私的サービス

• 2021~ liquid simulation : Implementation of Moving-Particle Semi-Implicit Method

🎨 I made ...

• 2024 あと何日？：今日から選択した日付まで、あと何日かを計算します。
• 2021 日経500風見鶏 : 投資家のツイートを重回帰分析・ベイズ推定し明日の株の値動きを予測するシステム
• 2017 gorogoro動画 : A video sharring website

🌱 I'm learning ...

• 2021 learning-gl : OpenGL practice

fluid-simulation's Issues

陰解法問題

正しくない

Q. そもそも収束しているのか？
収束していないならばガウスザイデル法の実装が誤っている
収束しているならば係数やベクトルに問題がある可能性が高い
t=0,100の推移を比べてみる

Q. nは連続的か
nが連続的かつガウスザイデル法が正しいならPも連続的になる

Q. 静水圧を壁面に作用させるのは正しいのか。
二章のコードで同じことをして弾け飛ぶか調べる

疑問解消

Q 平均粒子間距離で格子状に初期配置したのち、計算開始と同時に圧縮（例えば六方最密構造へ）、膨張などが生じないのはなぜか

粒子間距離で格子状に配置した場合の粒子数密度n_0を最初に求めている。同様の格子状初期配置の粒子数密度nは必ずn_0に等しくなるため圧縮も膨張も生じない。

Q なんでここにクーラン条件が出てくるの？DTならともかく DB = r * (1.0 + CRT_NUM); //バケット1辺の長さ

Q SNDをもっとわかりやすい値で表現できないか

段取りが良い実装

• 同時実行
  • 0.1秒後の水柱崩壊で答え合わせ 2D, 3D
  • 20cm2D静水圧
• C4のコードの圧力をC2の手法で求める

廃案

• 構造体を使っても使わなくても良い実装、MPS, CMPSを取り替えたい
  • offsetバイト数と基本型ポインタを引数に持つ、何にも依存しない関数の集合として粒子法を構成する
• MPS法のソースには近傍粒子検索が使われていない

images

自作コードは卒論に使えない

朝メール見てもうダムブレークできたんかと驚いたんだけど←嬉しい