【本号の目次】
１． 抗力係数CDを求める
２． 運動方程式を解く
３． ネバー・ギブアップ
４． 臼尻水産実験場ふたたび
抗力係数
を求める ここでは，図2-1に示すスルメイカの数値模型の抗力係数を数値流体力学による手法により求める。ここで，模型の長さは21cm, 厚さ1.8cm, 幅2.1cm とする。
を数値流体力学による手法により求める。ここで，模型の長さは21cm, 厚さ1.8cm, 幅2.1cm とする。
図2-1：スルメイカの数値模型 (a)は上から見た図，(b)は横から見た図
まず，数値流体力学の手法により数値的に求まった，物体の表面の応力（垂直応力および剪断応力）を面積分することにより抗力
を求める。この力を式(2-1)のように無次元化することにより 値が得られる。
（2－1）ここで用いた各パラメーターは図2-2に示すとおり。Vは速度、Sは体表面積、Densityは海水密度
図2－2
図2-3にイカの泳ぐ速度に依存する抗力係数
を示す。
図2－3：抗力係数
運動方程式を解く
ここでは，第三十話「イカの泳ぐ速度をどのように予測するか」で求めたイカの運動方程式(1-5)の離散化を行う。得られた離散化方程式を数値的に解くことにより，イカの泳ぐ速度を予測することができる。また，離散化式に含まれる抗力係数は前節で求まった値（図2-3）を用いる。
（１－５）まず，式(1-5)を式(1-5)’の形に変形する。
（１－５）’式(1－5)’を一次精度で離散化すると式(3-1)になる。
（３－１）式(3-1)を予測速度 について纏めると式(3-2)のようになる。
（３－２）これがイカの泳ぐ速度を予測するために解くべき式である。この式にいままでに推定できた値（表3-1）をこの式にあてはめてグラフ化すると図７がえられる。したがって、この体長21㎝の数値模型の最大速度は1.4m/s となる事がわかる。赤で示した噴水時間と漏斗の直径は仮の値である。
図７：速度履歴
―――――――――――――――――――――
と、ここまでが徳山さんのレポートである。
面積分てなんだ？無次元化てなんだ？？離散化方程式てなんだ？？？。高校二年時で数学を放棄した私には、まったく理解の及ばないところではあるが、体長21㎝程のスルメイカが直径5㎜の漏斗から0.1秒で外套内の海水を一吹きすると、スピードは1.4 m/秒（時速約5㎞）に達するという推定には納得できた。ただし、噴水時間と漏斗の直径の値は、ダミーである。それにイカは一吹きした後、すぐさま外套膜に水を吸い込んで漏斗から水を吹きだすことを繰り返して加速していく。徳山さんに無理を言ってスルメイカの遊泳速度を推定してもらったが、このようなシミュレーションをダイオウイカの遊泳に拡張するのはかなり無理があるし、論文としても纏まらないだろうと徳山さんから指摘された。
ネバー・ギブアップ しかし、せっかくここまで追求してきた目標を「そうですか、それではこれで止めにしましょう」とあきらめるには、未練が残った。昨年、臼尻水産実験所でおこなったスルメイカの実験を思い返しながら、なにか他にやれることはないかと頭を絞った。一つアイデアが閃いた。氷温麻酔をつかえば、スルメイカの泳ぐ力を計測することが出来るかもしれない。スピードではなく、パワーを測定するのだ。仕掛けは簡単である。氷温麻酔をかけてスルメイカを動かなくして、その外套膜の背側前縁に細いテグスで釣り針を掛け常温の海水にもどし、麻酔がさめて泳ぎだした際にテグスを引っ張る力（張力）を測定すれば、泳ぐ力が分かるはずである。
早速、北海道大学水産学部の桜井教授、臼尻水産実験所の宗原准教授と連絡をとり、臼尻で飼育しているスルメイカを使って遊泳力を測定する実験計画について相談にのってもらった。今年はスルメイカが豊漁で、実験所の飼育水槽には十分な量のスルメイカが畜養されているとのことで、快く実験の許可を頂いた。昨年スルメイカの吸水量測定の際、力を貸してくれた大学院生の高原君と大島さんが今回も実験を手伝ってくれることになった。漁業測器講座から張力計一式を借り出す算段もついた。
漁業測器教室から借用した張力計一式。左円筒：センサー、中央2台：アンプ、右：モニター
臼尻水産実験所ふたたび
臼尻水産実験所の実験水槽に収納されたスルメイカ
