https://ameblo.jp/einstein-eq/ 勉強にはコーヒーと音楽がつきもの。京都大学理学部、物理系4回生のとぅるーです。これ口癖です（笑）専門は非平衡理論です。興味ある分野は「粉体系」「散逸流子系」「漸近解析」「統計力学基礎」です。学部生もあと半年で終わると思うとしみじみします。 ja
とぅるーです。こんばんは。いやはや

卒業式も無事におわり
（自粛要請を振りだと解釈していろいろやりましたので、
聞きたい方は個人的にどうぞw）

卒研はポスター発表に一応参加できる程度にして
(もちろん研究としてはまだまだまだまだ)

今は流体ではなく粉体のシミュレーションをやっています。

大学4年間正直、いろいろな事がありすぎて振り返ったらキリがないので
まとめるということはしません。

ただ、楽しい4年間でした。

区切りということで、3年間ぐらいお世話になった

このblogに終止符を打とうと思います。

これまで、更新頻度がまばらで適当な僕のblogを

読んでくださってありがとうございました。

これからは大学院生という身分になり

またひと味もふた味も違った楽しい生活が待っています。

その生活の様子はここでは書きませんが、

http://blog.goo.ne.jp/tomo_1121

でもしかしたらチラ見できるかもしれないという噂です(笑)。

もしチラ見していただける方がいらしたらブックマークやリンクの変更をよろしくおねがいします。

それではだらだらと僕の醜態をさらすのは

ここらへんでやめて去る事にします。それでは＾＾／




]]> https://ameblo.jp/einstein-eq/entry-10847563445.html Fri, 01 Apr 2011 01:15:57 +0900
昨日引っ越しがありまして、ようやく平地に住むことになりました。

大学４年間は山の上で暮らしていたので、大学院では下界におりてくることにしました(どや)

少し広くなって快適です＾＾
もっとも、妹とルームシェアしている友人の家の広さにはかないませんでしたがw

卒業研究の方は、まあいつが終わりなののかよくわからんのですが、
2次元の流体シミュレーションは計算も終わって調べる量は調べたかなという感じで

来年のボスに報告したら、進んでよかったねと少しほめられた所です。
もちろん第一段階、第二段階を突破したにすぎないですが。

いま粘性流で計算を行っていますが１つのgeometryで計算させるのに数日かかり、
1回に2つずつやっているので論文提出までにぎりぎり間に合うか間に合わないか。

完全流体でのシミュレーションで研究会のポスターでは発表しようと思ってます。

あ、3/10にある京大ソフトマターという研究会でポスターやることになってまして
初めてポスター作ったんですが力加減がよくわからないので

明日担当教官にチェックしてもらいます。

調子に乗ってポスター全部英語で作ったけど

大丈夫かな？笑

それでは、がんばります。＾＾ノ
]]> https://ameblo.jp/einstein-eq/entry-10822926310.html Mon, 07 Mar 2011 00:54:48 +0900
最近はパソコンと会話する時間が多いかな。

＠セミナー室
とぅるー「おいおい、そう来たか...」

しーん. . .

元気ですよ笑

進み具合ですが、先日ボスに報告会があってパワポで
結果見せて、こうこういう方針で行こうといった矢先に

パワポで見せた結果が間違いである事がその日の夕方に発覚笑

修正したら見事に誤差が無くなっていい感じになりました。

３次元ダム決壊問題に移行したんですが
どうみても粒子数が足りない。

３次元ダム決壊のモチベーションとしては
３次元シミュレーションだとあるパラメータを下げないといけないらしいけど
どれくらいさげたらいいのかどこにも書いていない

→じゃあみつけるか。

と軽い気持ちでやったら実験値と全然合わない。
うむむ。

今日、担当教官に相談したら準２次元系でやって大きくしていくか
軸対称系でやるかどちらかになりましたが軸対称系だとかなりの量を書き換えないといけない。

けど再現したい系は軸対称系。。。
ここはおいおいいじっていく事にする。

２次元ダム決壊問題でできたムービーをのせておこう
ググったらこんなものいっぱい出てくると思うんですが笑




次にjetだけど実験系と同じスケールになるように変換して
計算させたら壁でやっぱりうまく行ってないみたいなので
時間刻みを小さくしている最中。

もうちょいでこれもうまくいくのかもしれないです
ただしこの台詞は２週間ぐらい前から言っている笑

どっちにしても軸対称系でモデルを考えて行く必要ありますな。

あとはそろそろ粉体サイドから詰めて行く方針を少しずつ考えなければならないです
一応少しはアイディアありますが、少し勉強が必要です。

1、２回生の頃ってよく理論系の院生って何してるんだろう？って思ってたけど
常に研究かゼミしてるといつのまにか時間すぎてるんだなと実感。
]]> https://ameblo.jp/einstein-eq/entry-10782980066.html Sat, 29 Jan 2011 00:38:01 +0900 今年もこのblogをよろしくおねがいします。一応4月になったら院行くので気分転換で
blogの場所を変えると思います。ちゃんと告知するのでご心配なく笑

今年からタイトルは普通にしますw
曲紹介は面倒くさいのでやめますが、気が向いたらやりますw

今年の目標は、今やっている研究に何かしらの答えを出す事です。
最も理想的な形は今年中に論文が書くことです、もちろん難しさは重々承知しています。

去年は年末に粒子法のプログラムを0から書く事ができたので
わりとそこらへんに関して自信がついた気がしています、自信がついたというよりむしろ
シミュレーションへの恐れみたいなものはなくなったと思っています。
なので研究の手法が少し増えたような、そうでないような笑

いつもうだうだ言ってますが、なんだかんだで
結果が全てだと思っているので、結果を出して行きたい
しかし着実にしっかりと。
そう思う2011年の始まりです。

年末にバグが取れないとか言ってたのはとれました！
見事に！

8000粒子とかで計算が走るようになってパソコンを一時期いじめてましたが

効率悪すぎたので1000粒子以下で計算させています。

今日はベンチマーク問題、本実験の問題をパラメータをいじって
jobをパソコンに投げて帰宅しました。明日には計算が終わっている計算なので
楽しみ。

そいえば昨日はゼミのメンバーで鍋でした。
その際買ったお酒のラベルが

「静寂の時」

どさんこ　「　これでとぅるーが静かになる」
とぅるー　　「ぁん？」

飲んだり喋ったりした後
雪が降りしきる中帰りましたとさ。



昨日の迷言　「足がガンダム」

]]> https://ameblo.jp/einstein-eq/entry-10764528240.html Mon, 10 Jan 2011 22:13:03 +0900
年末という事ですがパソコンで作業してます。
となりのウィンドウで計算させてます、がチラ見したところやっぱり
だめそうですorz初期ではうまくいってるんですがねえ。
時間経つとおかしくなるみたいで、といってるまにダメである事は発覚。。。笑

good bye 2010ということで誰得な感じもありますが、
周囲の方々にならって今年を振り返ってみよう。

日記すくねえ。と思いました、ついったーのつぶやき
が増えてしまったからしかたないかねえ。

1月
くりこみ群ハァハァの時期。くりこめ詐欺という言葉があって
「もしもし、俺だけど発散振り込んでくんねえ？」
Wilson KogutやGoldenfeldを読んでいた。

2月
Phase - Orderingについて勉強していた、Onukiゼミ用に。
あとはセミナーに少し行って視点が広がった。
Kさんの生物のセミナーはおもろかったな、仁科賞おめでとうございます。

3月
くりこみ群による漸近解析の論文を本格的に読み込んでいた。
射影演算子法と組み合わせで何かできないか、と探っていると
１０年前にやられているオチ。これらの論文を読んで漸近解析の雰囲気を感じれたし、
Navier - Stokes eq.やBoltzmann eq.をいじるいい勉強になった。
今でもこの時期のノートを参考にしている。さきがけ数学塾も行ったなー

4月
漸近解析とかできない現象ってないん？と思いこの時期から粉体の勉強を始める。
手始めに粉体ガスを勉強してた、がかなりムズイと同時に、チャレンジングだと思った。
うっすらと専門を粉体にしようと思い始めた、院試ゼミもスタートした。

5月
MacBookが届く。この月からMac畑の住人となりました。
Onukiゼミでの僕の発表が荒れたのはいい思い出（笑）

6月
サークルで病む、院試勉強でなえる。
Jammingとかfluctuation theoremを勉強していたようだ。

7月
教授に凸。すごく丁寧に質問に答えてくださる方だと思った。
研究紹介もさくっとしていただけたし志望を出した後だったけど
決めてよかったと思った。あと院試ゼミが佳境、復習に入った。
粉体系の粗視化について考えて計算もしていたけど、失敗に終わる。
学部生が数週間でできるもんじゃない。

8月
同時にサークルのごたごたでダメージからの院試乙。
今だから言えるけど、筆記が終わったとき、うかったわと思いました笑。

9月
合格発表と奨学金もらえる事発覚。
教授に挨拶。来年一緒の研究室のT君とも会う。
教授「研究したい内容のスライド作ってきて」
...
しかし人前で話す事になると色々考え直すことになる
のでためになりました。

10月
卒業研究の方向性でぐらぐらしていた。グリッドでやると
いっていたのが月末に粒子法でしようとなった。うおーう。

11月
本格的に数値計算の指導を受ける。初歩の初歩から教えてもらったので
1ヶ月の指導でかなり書けるようになったし、感じもつかめた。
本当に感謝しています！あと研究会連発。
粉体の研究会は楽しかった！月末はこのため研究はほぼせず。

12月
バグ取りなうw
とれねえ。粒子数を2000でも計算できるように改良したけど
どうも変な感じ。バグ取る前に新年を越せそうにはないです。笑


ふー。激動でした。
どちらかというと運には恵まれた一年だったかと思います。
来年は運勢悪いみたいですが、勝負の年と思っています。
地道に着実に歩いて行こうと思っています。


それではみなさんよいお年を。



おやすみ/ ELLEGARDEN ]]> https://ameblo.jp/einstein-eq/entry-10753206865.html Fri, 31 Dec 2010 02:16:18 +0900
まあ近況ですね、ぼちぼちです。

ゼミの方はゆらぎのエネルギー論が割と面白い所です。
ゆらぎが支配的な系(Langevin equationの階層)で熱力学でやったような

カルノーサイクルを考えると等温サイクルはwell-definedですが
系を孤立させる断熱仮定に値するエネルギーが一意に決まらないため
マクロな熱力学とは全く異なることがおきる。
（まだ読み切ってないから自信はないけど、
たぶんマクロ系のカルノーサイクルより効率があがる）

所謂ナノ系の基礎物理への興味が出た一幕。
ナノクラスターとかあたりから勉強始めると面白いかもしれないし、

連続体近似が破れるあたり
(Kn ~ 0.01 Knはクヌーセン数で平均自由行程と流れの特徴的な長さの比)
の流体（マイクロ流体）とかも境界の影響ががもろに効いて面白そうだ。

そして肝心の卒業研究は、進み具合がいまいちだ。


自由界面の問題なのでバグがどこで生じるかわからないとボスに言われたように
どこがバグの起源なのかわかったはずでも、わかっていないことになったり

たぶん圧力の計算でおかしくなってるんだろうけど、？？？？
まず静水平衡を

(アルゴリズムの説明が書かれた本を読んで0から書いた)僕の書いたスキームで

再現して、その後にベンチマーク問題であるダム決壊問題を再現してという手順。

やりたい現象までは遠い。やりたい現象は流入(と流出)があるからたぶんさらに大変。
果たして2月までに終わるのか？

周囲は割と結果が出始めたらしいけど、まだ僕はテスト問題。。
問題設定自体は結構速かったけど時間かかってるな。

もちろん担当教官に、それは難しい。粉体の方が簡単でしょ笑
と言われてたので覚悟はしていたけど、辛抱しかない。。

辛抱だ。



Centerfold / J. Geils Band
歌詞がせつないです、「高校時代のあの娘がピンナップに...」という歌。笑
]]> https://ameblo.jp/einstein-eq/entry-10740078833.html Sat, 18 Dec 2010 00:55:26 +0900
生きてます、12月ですね。。先生も走るといいいますが

僕も塾講師を始めましたが、かなり自由な感じで
こんなであんな額もらっていいのかという気分になります笑

近況ですか？別に対した事ないです笑



卒業研究　：　思った以上に進まないorz



につきます。２週間ぐらい前には
「もうできたに等しいぜ！！！！」と周囲に謳っていた僕でしたが

即　発　散

原因を調べた所圧力のPoisson方程式を解く所で詰んでいる様子。

境界条件から求まっている解があって、それを一定に保ちながら
計算が進むようなプログラムにしたはずが

どこか数直線の遥か彼方へ飛んで行きます。

逃げちゃダメだとプログラムの調査に入った所少しずつ間違いがわかってきましたが
未だに完了しません。。

担当教官いわく
「流体シミュレーションを始めてやるときは、皆圧力項の計算で泣くから
決して悪い所でミスしている訳ではないよ。」

ということなので、負けないようにがんばっていこうと思います。

土日でバグ取りきれたらいいんだけどなあ。

ていうか月曜日は

ボスに進捗状況報告なので

そろそろ、
「あれれーまだ実装してないのおお」
的な感じにつっこまれるかも。。。

おおおおおおおおおお。

たすけてーーw

あ、来週の土日はサークルの用事で東京に行きます。

traveling / 宇多田ヒカル
]]> https://ameblo.jp/einstein-eq/entry-10726015406.html Fri, 03 Dec 2010 22:11:48 +0900 引き受けてしまった仕事は全うせねばならないので、仕方ありません。

そのサークル関係の仕事で少し腹が立った事がありましたが
それはついったーのTLで書いたり喋ったりでだいぶすっきりしました。

簡単に言うと自分に酔っている不条理な事を言うおっさんを
怒らしてきたという感じです笑
沸点が低すぎて手応えがありませんでしたがw

今日でまあ、だいたいの用事は済んだので
ようやく復帰って感じがしますが。

１１月後半は研究会に参加したり、サークル関係で
研究や勉強をする時間がありませんでした。ただのいい訳ですが。

研究会はひさしぶりの先輩と喋ったりなかなか有意義でした。
特に粉体系の研究会ですが実験系の話が結構面白くてテンションがあがりました。
理論的に説明できるかわかりませんが。。

しかし粉体の研究をしようと思えば
やっぱり静力学(全ての粉体粒子が止まっていてぎちぎちに詰まったときのレオロジーなど)の知識
が必要な予感がしました。本は購入して少し読んでますが。

一部の基本知識(塑性など)は工学部では勉強するみたいで
M0の友人は聞いた事があるみたいでした。おそろしや工学部。。。


研究の方は数日パソコンと会話していなかったのですが、
バグはとれなさそうです笑

また担当教官に聞くしかないですかね。。



Automatic / 宇多田ヒカル
]]> https://ameblo.jp/einstein-eq/entry-10721121490.html Sun, 28 Nov 2010 21:27:46 +0900
気づけば1週間なんてあっというまにすぎている。

自分は何をしているのだろう、何をしたいのだろう、
妙な気持ちになる

秋はいやな季節だ。変にセンチメンタルになるし
変に自暴自棄になる。もうすぐ誕生日でまた１つ年をとるのに
成長していないと感じるし、

頭が固くなって行くように感じる。自分を少しでも大きく見せようとしたり
自分に嘘ついたり、高校のときから成長してない。

別に最近人間関係に悩んでいるとかじゃなくて、ただ
時の流れに思いをよせてるというか

風にでもなってしまいたいと思う今日この頃なんです。


卒研の調子ですか？
粒子法による流体シミュレーションの骨格はほぼかんせいしたので
ベンチマーク問題にチャレンジします。

初期条件と境界条件をうまく決めるだけで、
今日の夕方にだいたいできて、バグ取りをするだけで実装できそうです。
まあまだ2次元だし、粉体まで持って行くのはきっと大変です。

担当教官のおかげで数値シミュレーションのコードの書き方にだいぶ慣れました。
この点本当に感謝してます。わからないところは丁寧に教えて頂けるし
もっと教えてくれって言ったらがんがん教えて頂けるし、本当に感謝です。

たぶん粉体系のシミュレーションをするのであれば
時間発展に衝突をいれるだけだし、割と拡張は簡単にできるのかな？

純理論的にどこまでできるかわからないけど。

あ、そういえば今、研究会の期間です。
明日も朝早くからだ。

今日は最近ちまたを騒がした所謂Maxwell demonの話があった。
実験家ぱねえと思ったにつきます。Jarzynski等式の話は勉強してても


ふーん

って感じになってまうんよな。


あと、名前だけ知ってる人とかもいらっしゃって
ほーこんな人かと思った。ほー

明日はポスターセッションがある。


FINAL DISTANCE/宇多田ヒカル
]]> https://ameblo.jp/einstein-eq/entry-10711286204.html Thu, 18 Nov 2010 23:05:16 +0900
近況ですか、そうですね、来年のボスとのゼミを隔週やっています。
（通称M0ゼミor M0ミーティング）
ゼミとは行っても卒研の調子を報告するという感じのゼミです。
まだ１回しかありませんが来週の月曜にあります。報告が終わった後に
今度はD3の方がMDの講義をしてくださるそうです。

先週あってジェットのレビューをしたと書きました、その２週間後にあたります。
２週間でそんなに研究進まなくて、何報告すんべーとなっていまして、

そんなおりに先日M0の友人T君にばったりルネで出会い
最近研究どう的な探り合いをした所お互いに同じ感じだったらしいので

きっと今度のミーティングきっとはMDのレクチャーが主だよねという確認をしました笑

そして昨日の朝、寝ぼけながらメールの確認
「今度の月曜やる事はとぅるー、Tの研究の進捗状況。時間があったらMDのレクチャー。」


φ(_ _)zzZ → ちょ m9(^ _ ^;)


となりました笑
時間にして1時間30分しかないのでそんなもんかという気もしないでもないですが
少なくとも粒子法(SPH法とMPS法どちらも？)について解説をしようかと思っています。


といってもまだ両者の違いが微妙なんですよね。


自分の中では
MPS法：流体力学方程式を粒子の運動方程式にモデル化（勾配やLaplacianを）して時間発展させる。
SPH法：粒子の運動方程式を流体力学方程式にスムーズな関数で書き換える。


どちらもグリッドいらずでおもろいと思いますが
自分が扱っている問題が、自由界面の問題なので時間発展の際に
"どこが界面か？"を定義しなければならなくて前者の方法の方が
おそらく簡単なんですな。


これができるなと思ったので自分でうりゃーーーと本気出して
プログラム書き始めたのですが

圧力を決めるのに

陰的にLaplace方程式を解きますが、
MPSであればこれがただの線形の方程式になりますが

添字が粒子の番号なので粒子数2000の系をやろうとすると
2000×2000の行列の逆行列を求める事が必要で

＼(^o^)／

自明に0とおける成分がないので、ただの対称行列。。
と詰んでいる今日この頃。


Rock and Roll Music / The Beatles
]]> https://ameblo.jp/einstein-eq/entry-10696275388.html Wed, 03 Nov 2010 20:41:47 +0900