392dc54ce0a39523a52d6e7d12d2a2ff5b0cfdbb

【本号の目次】
１． 強力な助っ人
２． 徳山さんのレポート
３． 推進力
４． 抗力

強力な助っ人

　海洋生物学一辺倒で過ごしてきた私の周りには、そのような遊泳力のシミュレーションを考えてくれそうな研究者は見当たらなかった。そこで、姉貴の長男で東京大学大学院の複雑理工学専攻博士課程を修了した甥の小林徹也（現：東京大学生産技術研究所准教授）に相談をもちかけた。彼は東大で日本学術振興会特別研究員として「数理・情報で解き明かす生命現象」を研究テーマに、私には理解を超えた融合科学の最先端の研究をしていた。自分は忙しいので手伝いは出来ないが、興味をもってくれそうな研究者を知っているので話をしてみてはと、徳山佳央さん（M&T株式会社）を紹介してくれた。
　当時、徳山さんはシミュレーションを通じたコンピュータソフトウエアの企画・開発を目指す新しい会社を立ち上げたばかりで忙しい中、私の話を聞いてくれた。ダイオウイカの遊泳速度をシミュレーションするなんて、とても面白そうな企画ですねと頷いて、仕事の合間でよければちょっと考えてみますと請け負ってくれた。博物館からは、正式な業務として「ダイオウイカ数値模型作成受託業務」の発注依頼書を出してもらった。

徳山さんのレポート

　一か月ほど経ってから、徳山さんからレポートが届いた。思っていたより大変なことになってきた。以下に紹介するが、理解するには何度も読み返す必要があった。

１． イカの泳ぐ速度をどのように予測するのか？　

　一般に運動する物体の速度は、運動方程式を解くことにより求めることが可能である。運動方程式をたてるとは、力のつり合い式を求めることである。イカが泳いでいる時は推進力と抗力の力のバランスが成り立っている。

「F (受ける力）＝　6240fcecac93eec509108de34c79796646f17928（推進力）－　cf184cc149fb95a8156c7f0a5531f25935cbcf80（抗力）」

c60c0ff5a4fdfa2825ff4e45d7431da5609a21ff
　　図1-1：泳ぐイカの力の釣り合い

一般に運動方程式は次式のように書くことができる。
　　　　　　　f1fa63429a74fee831e6a6a125c3413b695487fd　　　　　（1-1）

この式をイカに当てはめて考えると，fbb96fb30da1f0746eedfa36b3c14d45144857c5 がイカに働く推進力 6240fcecac93eec509108de34c79796646f17928と抗力 cf184cc149fb95a8156c7f0a5531f25935cbcf80の合力で， e1151bc0042d909bb9dbe560a81fa0a61dea7041はイカの質量， fa716a189317e033a18384332ee04ca56f4efd6eはイカの泳ぐ速度になる。イカの泳ぐ速度は，この合力 Fを求め，運動方程式(1-1)を解くことにより求まる。

 推進力　6240fcecac93eec509108de34c79796646f17928

　まずは、推進力　6240fcecac93eec509108de34c79796646f17928について考える。
　推進力は漏斗から噴出した水から受ける力であり，作用反作用の法則により水がイカから受ける力に等しくなる。

a237d46281a456f75a2be70140e3939d9cc56057
　　　　　　　　　　　　　　　図1-2：推進力

この力は水の単位時間あたりの運動量の変化に等しく次式のように求まる。

　　　　　 fc6316a29be015d360fc190661fae71e91022874　　　（1-2）

ここで9ce7a5d8600387da5755adbcde11881c088f9f35 はイカの漏斗から噴出する水の流量 。また，a8135a6c32d1fe30f6efb1641b51ebae53f57450は水の密度、 977b432cc5071c4ad0bb9ca9def384807dc16ba9は漏斗の断面積で、 u は漏斗から噴出する水の速度だ。よって推進力 は式(1-3)のように求まる。

　　　　a8f6a7ef7c7f33ee0bf5c1f8ec464dbbddcfbe18　　　　　　　　（1-3）

水の密度a8135a6c32d1fe30f6efb1641b51ebae53f57450 および，漏斗の断面積977b432cc5071c4ad0bb9ca9def384807dc16ba9 は既知の量であるので，単位時間あたりに噴出する水量 Qを与えることにより，漏斗から噴出する流速u が求まり，6240fcecac93eec509108de34c79796646f17928推進力 が定まる．

※ Q を求める為には，イカが何秒で水を全て噴出しているかが重要になる。


 抗力 cf184cc149fb95a8156c7f0a5531f25935cbcf80

次に，抗力 cf184cc149fb95a8156c7f0a5531f25935cbcf80 について考えてみる。一般的に流体中にある物体に働く力 cf184cc149fb95a8156c7f0a5531f25935cbcf80は式(1-4)のようにモデル化される。

　　　　  7b57da7b291cd972357b84dd38633dd8e2b64273　　　　　（1-4）

ここで a24ef041039894ea9cf35153ca1553d5d9b88adcは抗力係数， a8135a6c32d1fe30f6efb1641b51ebae53f57450は流体の密度， V は物体と流体の相対速度 S は物体の代表面積。
残念ながら，この抗力係数 a24ef041039894ea9cf35153ca1553d5d9b88adcは物体の形状・流体の粘性・流れの速さなどにより変化し理論的に求めることは困難で，実験または数値流体解析により求めるのが一般的である。近年では、コンピューターや数値流体力学の発展により、抗力係数 a24ef041039894ea9cf35153ca1553d5d9b88adcなどを数値解析により求める手法が可能となってきた。

以上をまとめると，F = 6240fcecac93eec509108de34c79796646f17928-cf184cc149fb95a8156c7f0a5531f25935cbcf80 よりイカの運動方程式は式(1-5)のように書くことができる。

　　　　 ef1c7d9870e09985d0d8702f0956b4434e4f6838　　　（1-5）

ただし， a24ef041039894ea9cf35153ca1553d5d9b88adcは速度 V の関数値（実験値）である。この運動方程式(1-5)を解くことにより，イカの泳ぐ速度 を予測することができることになる。

と、ここまではなんとか付いて行けたが、問題は次のステップである。