今月の１６日から京都に研究会で行くのだけれど、どうやら僕は京都大学の宿舎を申し込んでなかったらしい。 それをキレた感じで京都大学の事務の人に今日知らされて急いでホテルを探した。１１月１６日ってのは紅葉の超ハイシーズンらしくて色々ホテルを見た…

三次元ダイナモの数値計算をおこなうためには表面での電流なしの磁場の境界条件が必要である。 そのサブルーチンを作っていたのだけれど今週なかなかうまくいかなかった。ついさっきすごいミスに気づいてようやく電流なし磁場の境界条件のサブルーチンが出来…

昨日の話だけど国立天文台の一般公開の手伝いに行ってきた。 僕は磁場の実験の手伝いをした。具体的にはペットボトルの中にサラダ油とスチールウールを入れて、さらに筒を通す。 そこに磁石を入れると双極子磁場の形が立体的に見えるというわけだ。こどもた…

今日のセミナーは僕の順番だった。 2時間半の長丁場でかなり疲れたが皆さんに考え方がわかってもらえてよかった。それからみんなで新宿にいすを買いにいったが高くて買えなかった。

来週のセミナーは自分の番だ。 太陽の内部角速度分布を実現するための研究について話したいが思ってたより量が多くなってしまいそうだ。 頑張って準備したい。

今日は駒場にて乱流セミナーに参加してきた。リサーチキャンパスなるところに初めて行ったがすごくきれいでうらやましかった。セミナー自体まず太陽の磁場活動について南極の氷床からとった10 Beを使って過去の活動周期を探るという事をやっていた。それなり…

三次元ダイナモコードはある意味でほとんど出来上がった。あとテストすることは二次元でやった時と周期などは同じになるかとdiv B=0がどの程度満たされているかである。

前回の日記の続きである。 赤道の境界を内側にとると負の動径方向速度は発生しなくなることがわかった。 しかし、これでは別の問題が発生する。

今はRempel(2005)の太陽の差動回転と子午面環流を実現するシミュレーションの再現をやっているのだけどどうにもうまくいかない。今の問題をとりあえず、メモしておく。 低緯度のtachoclineからそれより下のところでなぜか負の速度が発生してそれによって、そ…

太陽の内部構造を計算しているがなかなかうまくいかない。 かなり時間をかけているがうーん

今日は横山研セミナーに国立天文台の藤村さんが来た。色々とおもしろい話が聞けてよかった。それはそうと並列計算プログラミングの課題が大変だ。 明日が期限だけどなんとか終えたい。

http://www.solarmonitor.org/index.php?date=20090923そろそろ極大期に向かって走り出すか？

今日は三次補間をYin-Yangに実装した。 補間する点を原点として三次関数で補間すると切片の値のみがわかればいいという事に気づいた。並列計算プログラミングの課題もそろそろけりを付けなければいけない。最初の課題のみをやったので どうにか金曜までには…

並列計算プログラミングのためにcenjuにログインして色々やろうとしたがxemacsにてfortranモードにならないとかBackspaceがきかないとかいろいろあって進まなかったのだけど解決したのでメモをしておく。とにかくxemacsの設定ファイルは~/.xemacs/init.elと…

今日はhinode-3という国際会議への参加の許しがでたので申し込みをした。 英語がちゃんとしゃべれるかわからないのでそれまでにどうにか練習しなければいけないがどうしたものか？

やたら時間がかかってしまったが(書きだしたのは5月13日)今日やっと論文をThe Astrophysical Journalに投稿した。いちおうメールは来たがうまく投稿できたのかは微妙。取り合えずeditorからの連絡をまつことにする。なんとかacceptされるといいがさて…。

今日は天文学会最終日。今日東京に帰るのだけど、記事が書けないかもしれないので今のうちにいちおう書いておく。

今日は天文学会二日目。おもしろかった研究をまとめておく。MRIによる磁場生成太陽の表面近くの動径方向の勾配は磁気回転不安定が起きて磁場が増幅される。 先生に話したら、これは対流が起きたら起きないんじゃないかと言っていた。 うーん。grand minimaで…

今日は天文学会の発表。午前中は太陽系の発表を見た。 彗星などの観測の発表が多くて、あまりよくわからなかったが、天体が地球に衝突しようとするとき、どのような衝撃を与えれば地球からそらせるかという研究をしていたのが個人的におもしろかった。 結局…

天文学会のために山口に到着した。 朝4時半に起きて、東京駅に行き昼ごろに山口についた。 それから、秋芳洞と秋吉台に行った。かなり歩いたのでつらかったのだけど色々おもしろいものが見えてよかった。明日はとうとう発表。学会自体は二回目だけど今回は太…

明日は天文学会のために山口に行ってくる。 そのために今日は狛江でバッグを買った。安くていいものが買えたと思う。 あとは山口について色々調べていた。明日は観光をするために７時に東京集合なので早起きしなければいけない。がんばろう。とにかく発表が…

天文学会は来週の月曜日から。僕の発表順番は一番最初なので緊張するがしっかりやりたい。 今日のコロキウムは僕の研究室の先輩の発表だった。なぜかいつも質問する某先生が静かだった。 三次元ダイナモは項に分けて、やろうかと思ったがやはり無理だったよ…

今日は三次元のダイナモコードを引き続き作りつつ、新着セミナーと天文学会の資料を作った。新着セミナーは磁場が強くなると乱流が抑えられて拡散が減るという立場からダイナモの二次元のシミュレーションをやったもので僕はそれなりに期待していたがあまり…

今日は少し２流体不安定について勉強した。 電子と陽子が速度差をもっているとき、いろいろな不安定が起こるようだ。 今日見たブーネマン不安定は二つの粒子によってつくられる電場がうまいことなり不安定となると ただその不安定の条件は二つの速度差が小さ…

三次元太陽ダイナモコードにまた問題発生。 なにやら表面でおかしな挙動を示している。α効果がおかしいのだと思うけどどうすればいいのかさっぱりわからない。今日は昼頃にやってきて新着セミナー。不作の日でほとんど新着論文はなかった。 α効果についての…

今日は東洋大学にて、流体力学会年会2009。初めての学会発表だ。 事前登録を忘れていて、お金を払わなければいけなかったのでまず受付に行き 3000払った。この時大学の住所を書く欄があって郵便番号が分からなかったのだが黙って出した。さっそく宇宙惑星セ…

どうやら３次元太陽ダイナモコードはうまくまわり始めたようだ。これで一安心。 やはり、軸対称性を維持するにはそれなりのresolutionが必要な用でやはり大変な計算になりそうだ。うーむ。今日はそれから明日の流体力学会の発表の練習をした。 説明が多いと…

今日は非常に長い一日だった。最近体調が悪いのでインフルエンザかどうかチェックしに病院に行ったら、健康診断の時期だったらしく 二時間近く待たされた。そのおかげで「太陽」が二章くらい読めた。実際診察してもらうと一瞬ですぐにインフルエンザで ない…

今日も並列計算プログラミングの集中講義。 ただ風邪気味なため1時間遅れて行った。 ベクトル化のところで色々と有用なことを聞けてそれを自分のプログラムで試してみたら。 1.5倍くらいの速さになったのでよかった。朝はランチパック(たまご) 昼はメトロで…

今日は並列計算プログラミングの4日目。 午前中の前半は課題の解説。後半は線形ソルバーの解説。 前に国立天文台にRadiative Magnetohydrodynamicsのworkshopで行ったときに聞いた リーマンソルバーの話を思い出した。 午後は東大オープンスパコンT2Kの話を…