KABIRA https://blog.hatena.ne.jp/KABIRA/ KABIRAの日記 https://kabira.hatenadiary.org/ 拡散 シミュレーション 拡散方程式のシミュレーション 昨日の記事と同じベクトル場を用意する 境界からの流出入をなくす ソースのなかで (i != 1)や(i != Nx)の論理値 をかけることで境界では流出入が0になる 総量の保存について 初期条件によることであるが t=0 で36、 t=300 で 34.2 になった 端で流出している 今回は真ん中まで集まってくる やはりまだスケーリングを気にしないといけない 移動量の計算を３回繰り返しているの部分があるので、計算を工夫したほうがいいかもしれない X<-30 #空間 Y<-30 T<-100 #時間 dx<-1 dy<-1 dt<-0.5 #時間間隔 Nx<-X/d… 190 <iframe src="https://hatenablog-parts.com/embed?url=https%3A%2F%2Fkabira.hatenadiary.org%2Fentry%2F20110111%2F1293865211" title="移動項 境界 ベクトル場 - KABIRAの日記" class="embed-card embed-blogcard" scrolling="no" frameborder="0" style="display: block; width: 100%; height: 190px; max-width: 500px; margin: 10px 0px;"></iframe> https://cdn-ak.d.st-hatena.com/diary/KABIRA/2011-01-11.gif Hatena Blog https://hatena.blog 2011-01-11 16:00:11 rich https://kabira.hatenadiary.org/entry/20110111/1293865211 1.0 100%