1998-11-29[n年前へ]
■鴨川カップルの謎 
そうだ、京都、行こう
京都の風物詩の一つに「鴨川カップル」がある。京都を流れる鴨川の川縁に、カップルが等間隔に並ぶ現象である。鴨川の三条大橋から四条大橋までがその舞台である。この鴨川カップルについては、いくつかの性質が知られている。
- 基本的にはカップルらは等間隔に座る。
- 暗くなるに従い、カップル間の間隔が狭まる。
- 電灯の近く、すなわち、比較的明るいところではカップル間の間隔は広がる。
- 夏はカップルの間隔が狭まり、冬はカップルの間隔が広がる。
- 夏はカップル数が多く、冬は少ない。
- 男性が連れ立って、しかも何人かで座っていると、その周りにはカップルはなかなか座らない。
 (四条大橋の上から撮影) |
 |
 |
(この後、1999.05.04午前中に撮影した画像も示す。午前でもあるにも関わらず、上の時より密度が高い。もちろん、夜の密度の高さはこんなものではない)
 |
今回の目的は、「鴨川カップル」のこのような性質はなぜ生じ、そこに効いているパラメータを実測することである。
まず、性質2,3(暗くなるに従い、カップル間の間隔が狭まる、電灯の近く、すなわち、比較的明るいところではカップル間の間隔は広がる)より、明るいと隣のカップルと離れようとする力が大きくなると予想される。また、カップル同士が離れようとする力、斥力、のバランスにより性質1(基本的にカップルらは等間隔に座る)が形成されていると考えられる。
性質3(夏はカップルの間隔が狭まり、冬はカップルの間隔が広がる。)の季節による間隔の違いは性質4のカップル数の違いから来ているのかもしれないし、それ以外の何らかのパラメータがあるのかもしれない。性質4(夏はカップル数が多く、冬は少ない)はやはり「寒さ」のせいであろう。京都は本当に寒い。いくらアツアツのカップルとは言えども、寒いものは寒い(多分、想像だが)。当然、建物の中でのデートとなるだろう。自ずと、鴨川カップルは少なくなる。
性質5(男性が連れ立って、しかも何人かで座っていると、その周りにはカップルはなかなか座らない)については、「男性の群れ」というのは普通のカップルに比べて非常に強い斥力が働いていると考えられる。この性質5に関しては、私の実体験を持って語ることができる。私と友人連中が鴨川の川岸に近づくと、座っている鴨川カップルが離れていく、という経験を何度もしたことがある。やはり、強い斥力が働くのだ。
計算は全てMathematica3.0を使った。計算の例をMathematicaのNotebookで示す。また、計算の考え方を以下に示す。
 | 鴨川の川縁を真上から見た所。画面上部が鴨川。画面下部が歩道。画面中央に街灯がある。そのため、画面中央が一番明るい。 |
 | カップル達はこの辺りに座る。 |
 | 対称性を考えて、画面の左半分のみを考える。また、赤い線上にカップルが座るものとする。 |
 | その線上の明るさを模したもの。X=100の所が街灯の直下とする。 このような明かりに照らされた状態で、カップルが照らされているとする。 |
 | 100m内にカップルが5組として方程式を解くとこうなる。棒グラフの先端の位置がカップルの座っている位置を示している。1組目は0m地点、5組目は100m地点、すなわち、街灯の直下である。 明るい所ではカップルの間隔が離れているのが判ると思う。 カップルの居心地(隣のカップルの存在感=隣のカップルの明るさ/距離)はいずれも、8という値である。カップルの居心地指数は少ないほど快適であることを示す。 |
それでは、条件を振ってやってみたい。上からだんだん暗くなっていく。鴨川の半日である。また、カップル数はいずれも100m辺り5カップルである。
 |  |  | 20 |
 |  |  | 16 |
 |  |  | 10 |
|  |  | 8 |
 |  |  | 6.4 |
昼間はカップルが等間隔であり、夕方になり、街灯で照度分布ができると、カップルの間隔も分布ができている。といっても、そういう条件で解いているのだから、当たり前だが。
また、カップルの居心地は暗くなった夜の方が快適であるのがわかると思う。それは、カップルの様子を見ていても、その通りであると思う。しかし、快適であるからといって、カップルが何をやってもいいという意味ではない。そこは、はっきりしておきたい。
今回は、カップル数を全て同じにしたが、逆に同じ居心地指数であるという条件下で解けば、夜の方がカップルが多いという性質も再現できる。これらの計算モデルと実験を比較していくことにより、鴨川カップルの性質を実証していくことができるだろう。なお、今回は計算の簡単のため、男性連れの効果は組み入れていない。また、いつかもう少しまともな計算をしてみたい。
計算を行った感想だが、実に不毛な計算であった。気が向けば、モンテカルロシミュレーションによる鴨川カップルの検証も行いたい、と思う。が、気が向く日はきっと来ない。
2000-08-16[n年前へ]
■エアコンの風は心地よく吹くか？
真夏の夜の夢 流体力学入門編
私は長野県の野辺山という高原で幼い時期を過ごしたせいか、暑さにとても弱い体である。なので、真夏の夜はエアコンが欠かせない体と根性になってしまった。エアコン無しではろくな夢が見られ無いどころか、眠れなかったりするのである。気持ちの良い「真夏の夜の夢」を見るためには、エアコンがとっても重要なのである。
つい先日、そんなエアコンの話題が「今日の必ずトクする一言」に載っていた。それが、
- エアコン用流体素子(PATPEND.)のナゾ
- ( http://www.bekkoame.ne.jp/~jh6bha/higawari.html#000819 )
- 復活する冷蔵庫脱臭剤のナゾ(1/f揺らぎと活性炭編)
- ( http://www.bekkoame.ne.jp/~jh6bha/higawari.html#000901)
（おそらくhiraxさんあたりがやってくれるのではなかろうか）と話題のパスまでされている。いきなり、ボーとしているところを授業で当てられた気分である。
とはいえ、「エアコンの流体力学」というのもちょっと興味のある話題でもあるし、パスされたからにはやってみなければなるまい。そういうわけで、真夏の眠れない夜のパズル代わりに挑戦してみることにした。
さて、部屋の中の「エアコンから送られる風の様子」を計算するということは、流体力学の計算をするということになる。流体力学の運動方程式（ナビエ・ストークス式）に代表されるような方程式群を解かなければならないのである。そこで、以前
の時にいじりかけたMichael Griebel氏らによる非圧縮性流体のNast2Dのコードをもう一度引っ張り出して使ってみることにした。このNast2Dは非圧縮性二次元流れを計算する教育用のソースコードである。詳しくは- TUM INFO V- Homepage - ( http://www5.informatik.tu-muenchen.de/home_us.html )
- Numerical Simulation in Fluid Dynamics A Practical IntroductionISBN 0-89871-398-6
計算するモデルは次の図に示すような部屋である。四畳半一間であるかもしれないし、100畳以上の大きな広間かもしれないが、とにかく正方形の部屋だ。向かって左の青い丸部分にはエアコンがあり、そこから冷たい空気が送られてくるのである。そして、この部屋には何故かタンスがおいてある。「やっぱり、この部屋は四畳半一間じゃないの」というツッコミは言ってはいけない約束である。とりあえず、このタンスが向かって下側にある紫の部分である。このタンスが、エアコンから送られてくる冷た〜い空気流の障害物となるのである。
 |
暑い真夏の夜に、こんな部屋をエアコンで涼しくする時のことを考えてみよう。気持ちよく眠るために、何はともあれエアコンのスイッチを入れるわけだ。そうしないと、暑くて眠れないから当然である。
そして、そのエアコンには送風モードが何故か三つあるのだ。次のような三種類の送風、
- 真っ直ぐ送風するモード
- 単純な首振り送風をするモード
- 山本式エアコン用流体素子を用いた送風をするモード
空気流の速度分布の計算結果を次に示してみよう。
|
| を用いた送風の場合  |
また、これらの場合の計算結果を動画で示したものをMPEG4形式のAVIファイルとReal形式のファイルにしたものを以下に置いておく。エアコンから空気が送られる様子を知るには、何はともあれこの動画を見て頂きたい。なお、手元のRealProducerの制限のために、このReal形式のファイルは古いバージョンのRealPlayerだとアップデートが必要になってしまった。また、MPEG4のCodecが導入されていない場合には、DivXなどをインストールする必要がある。
とりあえず
- 真っ直ぐ送風している場合 MPEG4 AVI ( 289KB ) Real形式( 73KB )
- 単純な首振り送風の場合 MPEG4 AVI ( 577KB ) Real形式( 89KB )
- 山本式エアコン用流体素子を用いた送風の場合 MPEG4 AVI( 693KB ) Real形式 ( 92KB )
また、それぞれのMPEG4形式の動画のファイルサイズを見れば、それぞれの場合の送風による空気流の速度分布の複雑さは一目瞭然である、と思う。真っ直ぐ送風している場合はとにかく単純な速度分布であって動画ファイルサイズも結果的に小さくなっているのに対して、山本式エアコン用流体素子を用いた送風の場合はかなり複雑な速度分布になっていて動画ファイルも結果的に大きくなっているのである。
動画を見ることができない人のために、それぞれの場合の連続した3つの瞬間における部屋内の空気流の速度分布の静止画も示しておく。
1 | 1 | を用いた送風の場合 1  |
 |  | |
 |  |  |
これを見ると、「真っ直ぐ送風している場合」にはエアコンの正面は強烈に風が当たっている(つまり冷えまくっている)ことがわかるが、タンスの陰になっている部屋の隅などはほとんど空気が動いていないことがわかる。また、「単純な首振り送風の場合」は送風方向を動かしているとはいえ、その送風方向の変化はかなりゆっくりであって部屋の中の空気流の速度分布はそれほど急激には変化していないことも判ると思う(そういう風に計算しただけではあるが...)。いずれにせよ、もしもこの部屋の中にあなたがいたとしたら、体の決まった部分にのみ冷た〜い風があたることになるわけだ。それは体にはちょっとよろしくなさそうである。
一方、「山本式エアコン用流体素子を用いた送風の場合」は時々刻々と送風方向が変化しており、まるで部屋を舐め回すかのように、冷た〜い空気が送られていることがわかると思う。この部屋の中にもしあなたがいたとしても、体のごく一部分だけに冷たい風が当たるようなことはなく、それほど体に悪くないことが予想されるわけである。
さて、このページもかなり重いページになってきた。本来はこのタンスの裏側の空気が澱みやすい場所の「空気の入れ替わり」を、送風の具合を変えた上で調べてみたいわけであるが、それは次回のお楽しみ、ということにしておきたい。とりあえず、「真夏の夜の夢 流体力学入門編」はここら辺で終わりにしたい。
さて、夏休というわけでこんな景色のところで気持ちの良い風に吹かれてみたりするわけである。この写真の先の方の海の向こうに見えているのは左が伊豆半島で、右が三保の松原の辺である。
暑い夏の夜はエアコンが欠かせない体ではあるけれど、気持ちの上で言えばエアコンよりはおんぼろの扇風機の方が好きだし、扇風機よりも山の上の風の方がずっと好きだ。いつか、山の上を吹き抜ける風の音を録音して、その1/f揺らぎでも調べてみようかなと思うのである... それとも、そんなことをしてもツマラナイだけかな？
2002-10-19[n年前へ]
■Photoshop 鳥瞰図（立体グラフ）プラグイン
「Photoshop 鳥瞰図（立体グラフ）プラグイン」を作った。プリントダイアログがでないけれど、まぁ良しとしておく。表計算フィルタと統合したバージョンを明日作ってアップデートの予定。
2002-11-25[n年前へ]
■SDKライセンス変更？
SDKライセンス変更？ASNに入ると会員用のサイトで日本語化されたドキュメントなんかも見ることができるのだけれど、あぁいうのだってフリーで閲覧できるようにして欲しいなぁ、と思う。それに、ASNに加入する時の書類だってちょっとめんどくさすぎだと思う。少なくとも日本版はそうだった。
とりあえず、ASNでは11/18 Photoshop 7.0 SDK API Guide日本語版がアップデートされています。
2002-12-09[n年前へ]
■数式処理プラグイン+ATOK15+Windows98
の環境下で「入力しているウィンドウのフォーカスをいったん他に移してから戻してやらないと表示されない」という症状が起こるとのメール。こちらでもその症状を確認しましたので、適当に直してみました。というわけで、アップデート版に差し替えました。
具体的に直したところは、「AmetStartで他のWindowを作成してからすぐに消去していた」ところを「AmetResultでその処理を行う」ようにしました。
また、Fさんのリクエストに応えて、"list"コマンドで自分で定義した変数定義一覧を出力するようにしてみました。
Windows2000+ATOK13でのインストールが上手くいかない、というものについてはまた後ほど…（アドバイス募集中）。
」など。