(57)【要約】【課題】  推力偏向性能と高推力とを両立させた流体式推力方向制御装置を提供する。【解決手段】  流体式推力方向制御装置におけるノズル１０を、ノズル入口１１の下流側に設けられ、その下流側端部の断面積がその上流側端部の断面積よりも狭く形成され、その下流側端部が第一スロート１３とされた第一絞り部１２と、第一絞り部１２の下流側に設けられ、その下流側端部の断面積がその上流側端部の断面積よりも広く形成された第一広がり部１４と、第一広がり部１４の下流側に設けられ、その下流側端部の断面積がその上流側端部の断面積よりも狭く形成され、その下流側端部が第二スロート１６とされた第二絞り部１５と、第二絞り部１５の下流側に設けられ、ノズル出口１８に向かって断面積が広がるように形成された第二広がり部１７とを有するものとした。【選択図】  図２
2018178970.tif 000002
【特許請求の範囲】
【請求項１】
ノズル内をノズル入口からノズル出口に向かって流れる主流に対し、ノズルの中途部分に設けた二次噴流注入口から二次噴流を注入することによって、主流の流れに変化を起こさせ、全体としての噴流の向きを変更させる流体式推力方向制御装置であって、
ノズルが、
前記ノズル入口の下流側に設けられ、その下流側端部の断面積がその上流側端部の断面積よりも狭く形成され、その下流側端部が第一スロートとされた第一絞り部と、
第一絞り部の下流側に設けられ、その下流側端部の断面積がその上流側端部の断面積よりも広く形成された第一広がり部と、
第一広がり部の下流側に設けられ、その下流側端部の断面積がその上流側端部の断面積よりも狭く形成され、その下流側端部が第二スロートとされた第二絞り部と、
第二絞り部の下流側に設けられ、前記ノズル出口に向かって断面積が広がるように形成された第二広がり部と
を有することを特徴とする流体式推力方向制御装置。
【請求項２】
ノズルの中心線に対する第二広がり部の内壁の傾斜角度が３０〜７０°とされた請求項１記載の流体式推力方向制御装置。
【請求項３】
ノズルの前記ノズル入口の開口径（Ｄ１とする。）に対する、第二広がり部のノズルの中心線に沿った方向の長さ（Ｌ１とする。）の比Ｌ１／Ｄ１が、０．２〜２とされた請求項１又は２記載の流体式推力方向制御装置。
【請求項４】
ノズルとして、その中心線を含む一の平面に平行な断面の形状が同一のものを用いた請求項１〜３いずれか記載の流体式推力方向制御装置。
【請求項５】
ノズルとして、その中心線を軸とした回転体形状を為すものを用いた請求項１〜３いずれか記載の流体式推力方向制御装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ロケットエンジンやジェットエンジン等の推力の方向を制御するための流体式推力方向制御装置に関する。
【背景技術】
【０００２】
航空機や各種の飛翔体の多くは、空力操舵により運動の制御を行うが、動圧が低い状況下では、空力操舵が有効に働かない場合がある。例えば、宇宙空間や高高度において運用される宇宙ロケットのように、機体の周りが真空若しくは低密度の場合、又は、地上から発射された各種飛翔体の発射直後のように機体速度が非常に遅い場合に、空気力を用いた運動制御が著しく困難となる。このような条件下で機体を制御するために、推力方向制御（ＴＶＣ：  Ｔｈｒｕｓｔ  Ｖｅｃｔｏｒ  Ｃｏｎｔｒｏｌ）が用いられる。推力方向制御とは、ロケットエンジンあるいはジェットエンジンから排出される噴流の向きを何らかの方法で変えることにより、推力の向きを制御する技術である。これにより、空力操舵が有効でない状況でも、機体の姿勢や進行方向の制御が可能となる。また、近年では、航空機における高機動性を確保するための手段として、空力操舵とＴＶＣとが併用されるケースもある。
【０００３】
推力方向制御装置の方式は、機械式と流体式の２つに大別できる。機械式の推力方向制御装置（機械式推力方向制御装置）は、噴流の流路の形状を機械的に変化させるものとなっており、機械式推力方向制御装置としては、エンジンノズル内に傾斜板を設けるものや、ノズル形状を変化させるものや、エンジン自体の向きを変えるもの等がある。一方、流体式の推力方向制御装置（流体式推力方向制御装置）は、流れ自体の作用を用いるものであり、その代表例として、ノズル内を流れる主流に対して二次噴流を注入することにより、主流の流れに変化を起こさせ、全体としての噴流の向きを偏向させるものが挙げられる。流体式推力方向制御は、軽量化を図りやすく、応答性能を高くしやすいという点等で、機械式推力方向制御よりも優れている。
【０００４】
二次噴流を用いる流体式推力方向制御装置としては、デュアルスロートノズルを用いたものが提案されている（例えば、非特許文献１及び２を参照。）。ここで、「デュアルスロートノズル」とは、ロケットエンジン等のノズルとして通常用いられるラバールノズルの出口を絞って、第二スロートを設けた形状のノズルである。図１に、従来から提案されているデュアルスロートノズル１０’’の概念図を示す。同図において、符号「１０’’」はデュアルスロートノズルを、符号「１０ａ」はノズルの中心線を、符号「１１」はノズル入口を、符号「１２」は第一絞り部を、符号「１３」は第一スロートを、符号「１４」は第一広がり部を、符号「１５」は第二絞り部を、符号「１６」は第二スロートを、符号「１８」はノズル出口を、符号「１９」は二次噴流注入口を、符号「ＦＩＮ１」はエンジン燃焼室からノズル入口に流入する主流を、符号「ＦＩＮ２」は二次噴流注入口から注入する二次噴流を、符号「ＦＯＵＴ」はノズル出口から噴出する噴流を、符号「Ｆａ」、「Ｆｂ」、「Ｆｃ」及び「Ｆｄ」は代表的な流線を、符号「α１」は剥離領域（同図において網掛けハッチングで示した領域）をそれぞれ示している。
【０００５】
図１に示すように、二次噴流ＦＩＮ２により流線Ｆｄの近傍の流れは、デュアルスロートノズル１０’’における下側の内壁面から剥離し、二次噴流注入口１９と当該下側の内壁面における第二スロート１６との間に剥離領域α１を生じる。剥離領域α１の広がりは、デュアルスロートノズル１０’’の下側の内壁面における第二絞り部１５の存在により、第二絞り部１５が存在しないラバールノズルの場合よりも広範囲にわたり、その結果、主流ＦＩＮ１は上側に大きく押しやられる。しかし、デュアルスロートノズル１０’’の上側の内壁面付近においては、ノズル出口１８の直前の第二絞り部１５で逆方向に押し戻される。このため、ノズル出口１８から噴出する噴流ＦＯＵＴは、二次噴流ＦＩＮ２がない場合における右向きから、右下を向くように偏向し、その反作用として生じる推力が、二次噴流ＦＩＮ２のない場合における左向きから、左上を向くように偏向する。このようなメカニズムにより、ラバールノズルを用いた場合よりも、大きな推力偏向角を得ることが可能となっている。
【先行技術文献】
【非特許文献】
【０００６】
【非特許文献１】Ｆｌａｍｍ，Ｊ．Ｄ．，Ｄｅｅｒｅ,，Ｋ．Ａ．，Ｍａｓｏｎ，Ｍ．Ｌ．，Ｂｅｒｒｉｅｒ，Ｂ．Ｌ．，ａｎｄ  Ｊｏｈｎｓｏｎ，Ｓ． Ｋ．，“Ｅｘｐｅｒｉｍｅｎｔａｌ  Ｓｔｕｄｙ  ｏｆ  ａｎ  Ａｘｉｓｙｍｍｅｔｒｉｃ  Ｄｕａｌ  Ｔｈｒｏａｔ  Ｆｌｕｉｄｉｃ  Ｔｈｒｕｓｔ  Ｖｅｃｔｏｒｉｎｇ  Ｎｏｚｚｌｅ  ｆｏｒ  Ｓｕｐｅｒｓｏｎｉｃ  Ａｉｒｃｒａｆｔ  Ａｐｐｌｉｃａｔｉｏｎ”，ＡＩＡＡ  Ｐａｐｅｒ，２００７−５０８４，２００７．
【非特許文献２】Ｓｈｉｎ，Ｃ．Ｓ．，Ｓｕｒｙａｎ，Ａ．，Ｋｉｍ，Ｈ．Ｄ．，ａｎｄ  Ｓｅｔｏｇｕｃｈｉ，Ｔ．，“Ａ  Ｃｏｍｐｕｔａｔｉｏｎａｌ  Ｓｔｕｄｙ  ｏｎ  ｔｈｅ  Ｓｕｐｅｒｓｏｎｉｃ  Ｆｌｏｗ  ｉｎ  ａ  Ｄｕａｌ  Ｔｈｒｏａｔ  Ｎｏｚｚｌｅ”，Ｐｒｏｃ． １３ｔｈ  Ａｓｉａｎ  Ｃｏｎｇｒｅｓｓ  ｏｆ  Ｆｌｕｉｄ  Ｍｅｃｈａｎｉｃｓ，ｐｐ．３２９−３３３，２０１０．
【発明の概要】
【発明が解決しようとする課題】
【０００７】
上記のデュアルスロートノズルは、推力方向の偏向を効果的に行うよう工夫されたものであるが、その一方で、ノズル出口の断面積が絞られているため、推力自体は低下するという欠点を有している。というのも、超音速流においては、流体の圧縮性の影響により、断面積がより広い所で流速がより大きくなるため、ラバールノズルでは、この効果を利用してノズル出口から高速で流体を噴出させ高推力を得ているところ、デュアルスロートノズルの出口に設けた第二スロートが、この効果を一部打ち消してしまうからである。一方、ラバールノズルに二次噴流を注入した場合、推力偏向効果は、デュアルスロートノズルよりも著しく劣ることが判明している。したがって、これらのノズル形状を用いる限り、推力偏向性能と高推力とを両立させることはできず、一方を重視すれば他方の性能を犠牲にせざるを得なかった。
【０００８】
本発明は、上記課題を解決するために為されたものであり、推力偏向性能と高推力とを両立させた流体式推力方向制御装置を提供するものである。
【課題を解決するための手段】
【０００９】
上記課題は、
ノズル内をノズル入口からノズル出口に向かって流れる主流に対し、ノズルの中途部分に設けた二次噴流注入口から二次噴流を注入することによって、主流の流れに変化を起こさせ、全体としての噴流の向きを変更させる流体式推力方向制御装置であって、
ノズルが、
前記ノズル入口の下流側に設けられ、その下流側端部の断面積がその上流側端部の断面積よりも狭く形成され、その下流側端部が第一スロートとされた第一絞り部と、
第一絞り部の下流側に設けられ、その下流側端部の断面積がその上流側端部の断面積よりも広く形成された第一広がり部と、
第一広がり部の下流側に設けられ、その下流側端部の断面積がその上流側端部の断面積よりも狭く形成され、その下流側端部が第二スロートとされた第二絞り部と、
第二絞り部の下流側に設けられ、前記ノズル出口に向かって断面積が広がるように形成された第二広がり部と
を有することを特徴とする流体式推力方向制御装置
を提供することによって解決される。
【００１０】
このように、ノズル（デュアルスロートノズル）の出口直前に、第二広がり部を付加することにより、第二スロートで速度が低下した噴流を第二広がり部で再加速させることができる。このため、本発明の流体式推力方向制御装置は、第二広がり部を有さない従来から提案されているデュアルスロートノズルを用いた場合と比較して、推力を増加させることが可能なものとなっている。また、本発明の流体式推力方向制御装置におけるノズルには、従来から提案されているデュアルスロートノズルと同様、ノズル出口付近に第二スロートが存在している。このため、本発明の流体式推力方向制御装置は、第二スロートを有さない従来のラバールノズルを用いた場合と比較して、優れた推力偏向性能を得ることができるものとなっている。
【００１１】
本発明の流体式推力方向制御装置において、ノズルの中心線に対する第二広がり部の内壁の傾斜角度（図２における傾斜角度θ１を参照。傾斜角度θ１が場所によって異なる場合にはその平均値。）は、０°よりも大きく、９０°よりも小さければ、特に限定されない。しかし、この傾斜角度θ１を小さくしすぎる（第二広がり部の広がりが弱い）と、推力偏向角度を大きくしにくくなる虞がある。このため、傾斜角度θ１は、３０°以上とすると好ましい。後述するように、傾斜角度θ１が３０°以上であれば、従来のラバールノズルよりも推力偏向角度を大きくできることが確認できている。一方、この傾斜角度θ１を大きくしすぎると、推力を大きくしにくくなる虞がある。このため、傾斜角度θ１は、７０°以下とすると好ましい。後述するように、傾斜角度θ１が７０°以下であれば、従来から提案されているデュアルスロートノズルよりも推力を大きくできることが確認できている。
【００１２】
また、本発明の流体式推力方向制御装置において、ノズルの前記ノズル入口の開口径（図２における開口径Ｄ１を参照。前記ノズル入口が非円形である場合には、等価円直径で定義する。）に対する、第二広がり部のノズルの中心線に沿った方向の長さ（図２における長さＬ１を参照。）の比Ｌ１／Ｄ１は、特に限定されないが、比Ｌ１／Ｄ１を小さくしすぎると、第二広がり部による高推力化が限定的になる虞がある。このため、比Ｌ１／Ｄ１は、０．２以上とすると好ましく、０．３以上とするとより好ましく、０．４以上とするとさらに好ましい。一方、比Ｌ１／Ｄ１を大きくしすぎると、偏向性能が低下する虞がある。このため、比Ｌ１／Ｄ１は、２以下とすると好ましい。
【００１３】
さらに、本発明の流体式推力方向制御装置において、ノズルは、ノズル入口からノズル出口に向かって、第一絞り部、第一スロート、第一広がり部、第二絞り部、第二スロート、第二広がり部の順で各部が設けられたものであれば特に限定されない。ノズルとしては、その中心線を含む一の平面に平行な断面の形状が同一のもの（いわゆる二次元ノズル）や、その中心線を軸とした回転体形状を為すもの等が挙げられる。
【発明の効果】
【００１４】
以上のように、本発明によって、推力偏向性能と高推力とを両立させた流体式推力方向制御装置を提供することが可能になる。本発明の流体式推力方向制御装置を用いると、従来から提案されているデュアルスロートノズルを用いたエンジンよりも推進性能を高めるとともに、従来のラバールノズルを用いたエンジンよりも推力方向の制御性能を改善することも可能である。したがって、航空機や飛翔体の飛行性能や機動性を向上させることも可能となる。
【図面の簡単な説明】
【００１５】
【図１】従来から提案されているデュアルスロートノズルの概念図である。
【図２】本発明の流体式推力方向制御装置で用いるノズルの概念図である。
【図３】本発明の流体式推力方向制御装置で用いるノズルの形態例を示した斜視図である。
【図４】本発明に係るノズルの計算モデルを示した図である。
【図５】従来のラバールノズルの計算モデルを示した図である。
【図６】従来から提案されているデュアルスロートノズルの計算モデルを示した図である。
【図７】テストシミュレーションにより得られたラバールノズルの中心線上のマッハ数Ｍの分布を示したグラフである。
【図８】テストシミュレーションにより得られたラバールノズルの中心線上の流体密度ρの分布を示したグラフである。
【図９】テストシミュレーションにより得られたラバールノズルの中心線上の絶対温度Ｔの分布を示したグラフである。
【図１０】本シミュレーションにより得られた結果を示したグラフである。
【発明を実施するための形態】
【００１６】
１．本発明の流体式推力方向制御装置の概要
本発明の流体式推力方向制御装置の好適な実施態様について、図面を用いてより具体的に説明する。図２は、本発明の流体式推力方向制御装置で用いるノズル１０の概念図である。同図において、符号「１０」はノズルを、符号「１０ａ」はノズルの中心線を、符号「１１」はノズル入口を、符号「１２」は第一絞り部を、符号「１３」は第一スロートを、符号「１４」は第一広がり部を、符号「１５」は第二絞り部を、符号「１６」は第二スロートを、符号「１７」は第二広がり部を、符号「１８」はノズル出口を、符号「１９」は二次噴流注入口を、符号「Ｄ１」はノズル入口の開口径を、符号「Ｌ１」は第二広がり部のノズルの中心線に沿った方向の長さを、符号「θ１」はノズルの中心線に対する第二広がり部の内壁の傾斜角度を、符号「ＦＩＮ１」はエンジン燃焼室からノズル入口に流入する主流を、符号「ＦＩＮ２」は二次噴流注入口から注入する二次噴流を、符号「ＦＯＵＴ」はノズル出口から噴出する噴流をそれぞれ示している。
【００１７】
上記の図２は、ノズル１０を、その中心線１０ａを含む平面（ｘｚ平面に平行な平面）で切断した断面図として描いている。ノズル１０の立体的形状としては、図３（ａ）に示すように、その中心線１０ａを含む一の平面（図中のｘｚ平面）に平行な断面の形状が同一のもの（いわゆる二次元ノズル）や、図３（ｂ）に示すように、その中心線１０ａを軸とした回転体形状を為すもの等が挙げられる。図３は、本発明の流体式推力方向制御装置で用いるノズル１０の形態例を示した斜視図である。
【００１８】
本発明の流体式推力方向制御装置は、図２に示すように、ノズル１０内をノズル入口１１からノズル出口１８に向かって流れる主流ＦＩＮ１に対し、ノズル１０の中途部分に設けた二次噴流注入口１９から二次噴流ＦＩＮ２を注入することによって、主流ＦＩＮ１の流れに変化を起こさせ、全体としての噴流ＦＯＵＴの向きを変更することが可能なものとなっている。
【００１９】
ノズル１０において、ノズル入口１１の下流側には、第一絞り部１２が形成されている。この第一絞り部１２は、上流側（ｘ軸方向負側）から下流側（ｘ軸方向正側）に向かって絞られて形成された部分となっており、第一絞り部１２の下流側端部の断面積（開口面積）は、第一絞り部１２の上流側端部の断面積（開口面積）よりも狭くなっている。第一絞り部１２の下流側端部は、第一スロート１３となっている。
【００２０】
また、ノズル１０における第一絞り部１２の下流側には、第一広がり部１４が形成されている。この第一広がり部１４は、上流側から下流側に向かって広げられて形成された部分となっており、第一広がり部１４の下流側端部の断面積（開口面積）は、第一広がり部１４の上流側端部の断面積（開口面積）よりも広くなっている。
【００２１】
さらに、ノズル１０における第一広がり部１４の下流側には、第二絞り部１５が形成されている。この第二絞り部１５は、上記の第一絞り部１２と同様、上流側から下流側に向かって絞られて形成された部分となっており、第二絞り部１５の下流側端部の断面積（開口面積）は、第二絞り部１５の上流側端部の断面積（開口面積）よりも狭くなっている。第二絞り部１５の下流側端部は、第二スロート１６となっている。
【００２２】
ここまでは、図１に示したデュアルスロートノズル１０’’と同様である。しかし、本発明の流体式推力方向制御装置に係るノズル１０では、図２に示すように、第二絞り部１５のさらに下流側に、第二広がり部１７が形成されている。この第二広がり部１７は、上流側から下流側のノズル出口１８に向かって広げられて形成された部分となっており、第二広がり部１７の下流側端部の断面積（開口面積）は、第二広がり部１７の上流側端部の断面積（開口面積）よりも広くなっている。
【００２３】
また、ノズル１０の周壁部には、ノズル１０内に二次噴流ＦＩＮ２を注入するための二次噴流注入口１９が設けられている。この二次噴流注入口１９は、通常、第二広がり部１７よりも上流側に設けられ、好ましくは、第二絞り部１５よりも上流側に設けられる。本実施態様のノズル１０において、二次噴流注入口１９は、第一スロート１３付近に設けている。
【００２４】
二次噴流注入口１９は、図示省略のガス移送手段に接続され、当該ガス移送手段と二次噴流注入口１９とを接続するガス移送路には、当該ガス移送路の開閉を行うための開閉弁が設けられる。二次噴流注入口１９に接続するガス移送手段は、主流ＦＩＮ１を供給するエンジン燃焼室として、このエンジン燃焼室からのガス移送路を、主流ＦＩＮ１を移送する主流移送用のものと、二次噴流ＦＩＮ２を移送する二次噴流移送用のものとに分岐させてもよい。
【００２５】
二次噴流注入口１９を設ける個数は、特に限定されないが、１つしか設けていないと、推力を同じ側にしか偏向させることができない。このため、二次噴流注入口１９を複数箇所に設け、それらのうちの一部の二次噴流注入口１９に接続されたガス移送路の開閉弁のみを開くことにより、目的とする方向に推力を偏向させることが可能となる。例えば、図３（ａ）に示した二次元ノズルの場合には、二次噴流注入口１９は図２における下側の二次噴流注入口１９ａと上側の二次噴流注入口１９ｂとの二つを設け、推力を偏向させようとする向きに応じて、いずれか一方のみの二次噴流注入口１９から二次噴流ＦＩＮ２を注入させる。これにより、推力の向きを図３（ａ）におけるｙ軸に垂直な面内で、左上または左下の方向に向けることが可能となる。推力を偏向させる角度の調節は、ガス移送路の開閉弁の開度を調節して二次噴流ＦＩＮ２の流量を調節することにより、行うことが可能である。一方、図３（ｂ）に示した回転体形状のノズルの場合は、二次噴流注入口１９を、ノズル１０の中心線１０ａ回りの１８０°を超える範囲に少なくとも３箇所に設け、そのうち一部の二次噴流注入口１９のみから二次噴流ＦＩＮ２を注入させ、それぞれの二次噴流注入口１９からの二次噴流ＦＩＮ２の流量を調節することにより、推力を、中心線１０ａ回りの３６０°の範囲のいずれの向きとすることも可能となる。
【００２６】
上記のノズル１０を備えた流体式推力方向制御装置による、推力の偏向原理は、以下の通りである。
【００２７】
図示省略のエンジン燃焼室（ロケットエンジンやジェットエンジン等の燃焼室）から送出された燃焼ガスが主流ＦＩＮ１として、ノズル入口１１からノズル１０内に流入し、噴流ＦＯＵＴとしてノズル出口１８から噴出するようになる。この噴流ＦＯＵＴの反作用によって、噴流ＦＯＵＴの持つ総運動量と逆向きの力をエンジンが受けるようになり、推力が発生する。
【００２８】
ノズル１０は、その内部に主流ＦＩＮ１のみが流入している場合（二次噴流ＦＩＮ２が注入されていない場合）において、噴流ＦＯＵＴ全体が持つ運動量がノズル１０の中心線１０ａに平行な方向の一側（ｘ軸方向正側）を向き、その他側（ｘ軸方向負側）を向く推力が生じるように設計されている。
【００２９】
ここで、例えば推力の向きを、ｚ軸方向正側に傾けたい場合（図２の紙面における左上方向に偏向させようとする場合）には、同図中の下側の二次噴流注入口１９ａから、ノズル１０内に二次噴流ＦＩＮ２を注入する。すると、主流ＦＩＮ１は、図１に示したデュアルスロートノズル１０’’の場合と同様に、ノズル１０内で上側（ｚ軸方向正側）に押しやられるが、第二絞り部１５によって右下方向を向くように下側（ｚ軸方向負側）に曲げられるので、第二スロート１６を通過する噴流全体の運動量は、下向きの成分を持つようになる。その結果、エンジンの推力は、図２の紙面における左上方向に偏向する。このとき、第二広がり部１７によって、ノズル出口１８の断面積（開口面積）は、第二スロート１６の断面積（開口面積）よりも広がるため、第二広がり部１７内を通過するガスが加速され、第二広がり部１７を設けていない場合よりも、噴流ＦＯＵＴの総運動量の絶対値が大きくなり、推力が増大する。
【００３０】
一方、推力の向きを、ｚ軸方向負側に傾けたい場合（図２における左下方向に偏向させようとする場合）には、図中の上側の二次噴流注入口１９ｂより、ノズル１０内に二次噴流ＦＩＮ２（図示省略）を注入する。これにより、下側の二次噴流注入口１９ａから二次噴流ＦＩＮ２を注入した場合と同様の原理によって、推力は、図２の紙面における左下方向に偏向する。
【００３１】
２．シミュレーション（本シミュレーション及びテストシミュレーション）
本発明の流体式推力方向制御装置に係るノズルが、推力偏向性能と高推力性とに関し、どの程度の効果を発揮できるものであるかを確認するため、ノズル内の流れのシミュレーションを行った。
【００３２】
２．１  シミュレーション方法
２．１．１  本シミュレーションの方法
このシミュレーションにおいて、本発明のノズル１０は、図３（ａ）に示したものと同様に、ｙ軸方向に垂直な任意断面の形状が同一となる二次元ノズルとして定義した。本発明のノズル１０の断面（ｙ軸方向に垂直な断面）における各部の寸法（ノズル入口１１の開口径Ｄ１に対する各部の寸法）は、図４に示す値に設定した。図４は、本発明に係るノズル１０の計算モデルを示した図である。
【００３３】
このシミュレーションでは、ノズル入口１１からノズル１０内に主流ＦＩＮ１（図２を参照）を流入させるとともに、二次噴流注入口１９ａから二次噴流ＦＩＮ２（図２を参照）も注入した。二次噴流ＦＩＮ２の流量は、主流ＦＩＮ１の流量（「ＱＩＮ１」とする。）と二次噴流ＦＩＮ２の流量（「ＱＩＮ２」とする。）との比（ＱＩＮ１：ＱＩＮ２）が９５：５となるように設定した。また、ノズル１０の中心線１０ａに対する二次噴流ＦＩＮ２の注入角度は、３０°（デュアルスロートノズルでの予備計算において最大の推力偏向角度が得られた値）に設定した。
【００３４】
また、このシミュレーションにおいて、本発明に係るノズル１０は、第二広がり部１７の長さＬ１及び傾斜角度θ１の異なる計１４パターン（実施例１．１〜１．７及び実施例２．１〜２．７）についてシミュレーションを行った。実施例１．１〜１．７及び実施例２．１〜２．７における第二広がり部１７の長さＬ１及び傾斜角度θ１の値は、下記表１に示す通りである。
【００３５】
【表１】
2018178970.tif 000003
【００３６】
また、このシミュレーションでは、上記の本発明に係るノズルの範疇に属するもの（実施例１．１〜１．７及び２．１〜２．７）のほか、従来から提案されているラバールノズルに属するもの（比較例１）と、従来から提案されているデュアルスロートノズルに属するもの（比較例２）とについても行った。図５は、比較例１として採用した、従来のラバールノズル１０’’の計算モデルを示した図である。図６は、比較例２として採用した、従来から提案されているデュアルスロートノズル１０’の計算モデルを示した図である。
【００３７】
比較例１，２のノズル１０’，１０’’も、本発明のノズル１０と同様、ｙ軸方向に垂直な任意断面の形状が同一となる二次元ノズルとして定義した。比較例１，２のノズル１０’，１０’’の断面（ｙ軸方向に垂直な断面）における各部の寸法（ノズル入口１１の開口径Ｄ１に対する各部の寸法）は、それぞれ図５，６に示す値に設定した。図６のデュアルスロートノズル１０’’は、図４の本発明に係るノズル１０から第二広がり部１７を取り除いたこと以外は、図４の本発明に係るノズル１０と同じ寸法形状を有している。また、図５のラバールノズル１０’は、図６のデュアルスロートノズル１０’’から第二絞り部１５を取り除いたこと以外は、図６のデュアルスロートノズル１０’’と同じ寸法形状を有している（図４の本発明に係るノズル１０から第二広がり部１７と第二絞り部１５を取り除いたこと以外は、図４の本発明に係るノズル１０と同じ寸法形状を有している）。
【００３８】
比較例１，２におけるシミュレーションでも、上記の実施例１．１〜１．７及び２．１〜２．７と同様、ノズル入口１１からノズル１０’，１０’’内に主流ＦＩＮ１（図２を参照）を流入させるとともに、二次噴流注入口１９ａから二次噴流ＦＩＮ２（図２を参照）も注入した。主流ＦＩＮ１の流量ＱＩＮ１と二次噴流ＦＩＮ２の流量ＱＩＮ２との比ＱＩＮ１：ＱＩＮ２や、ノズルの中心線に対する二次噴流ＦＩＮ２の注入角度も、上記の実施例１．１〜１．７及び２．１〜２．７と同一に設定した。
【００３９】
上記のシミュレーションにおける計算は、モンテカルロ直接シミュレーション（ＤＳＭＣ:  Ｄｉｒｅｃｔ  Ｓｉｍｕｌａｔｉｏｎ  Ｍｏｎｔｅ  Ｃａｒｌｏ）法において、分子モデルを剛体球とし、ノズルの内周壁での反射条件を鏡面反射として行った。推力（ベクトル量「ｆ」とする。）は、微小時間（「Δｔ」とする。）の間にノズル出口から噴出した分子の総運動量（ベクトル量「Δｐ」とする。）の値から、下記式１を用いて求めた。また、推力偏向角度（「δ」とする。）は、下記式２を用いて求めた。
【００４０】
【数１】
2018178970.tif 000004
【数２】
2018178970.tif 000005
【００４１】
２．１．２  テストシミュレーションの方法
また、上記のシミュレーションで得られた結果の妥当性を確認するため、上記のシミュレーション（以下、「本シミュレーション」と呼ぶ。）に先立って、モンテカルロ直接シミュレーション（ＤＳＭＣ法）によるテストシミュレーションも行った。このテストシミュレーションでは、上記の比較例１のラバールノズル１０’（図５）において、二次噴流注入口１９から二次噴流を注入しない状態で、ノズル入口１１から主流のみを流入させて行った。このテストシミュレーションでは、ノズル１０’の中心線上におけるマッハ数Ｍ、流体密度ρ及び絶対温度Ｔの値の分布を算出し、これらの算出されたそれぞれの値を、ノズル入口１１、第一スロート１３及びノズル出口１８における、等エントロピー流理論による予測値と比較した。
【００４２】
２．２  シミュレーションの結果
２．２．１  テストシミュレーションの結果
まず、テストシミュレーションの結果について説明する。
【００４３】
図７は、テストシミュレーションにより得られたラバールノズル１０’の中心線上のマッハ数Ｍの分布を示したグラフである。図７における「Ｍａ」、「Ｍｂ」及び「Ｍｃ」は、それぞれ、ノズル入口１１、第一スロート１３及びノズル出口１８における、等エントロピー流理論によるマッハ数Ｍの予測値である。図８は、テストシミュレーションにより得られたラバールノズル１０’の中心線上の流体密度ρの分布を示したグラフである。図８における「ρａ」、「ρｂ」及び「ρｃ」は、それぞれ、ノズル入口１１、第一スロート１３及びノズル出口１８における、等エントロピー流理論による流体密度ρの予測値である。図９は、テストシミュレーションにより得られたラバールノズル１０’の中心線上の絶対温度Ｔの分布を示したグラフである。図９における「Ｔａ」、「Ｔｂ」及び「Ｔｃ」は、それぞれ、ノズル入口１１、第一スロート１３及びノズル出口１８における、等エントロピー流理論による絶対温度Ｔの予測値である。図７〜９において、横軸（ｘ軸）は、ノズル１０’の中心線に沿う座標である。図８における流体密度ρ及び図９における絶対温度Ｔは、それぞれ、ノズル入口１１の流体密度ρ及び絶対温度Ｔを「１」として規格化した値で示している。
【００４４】
図７〜９を見ると、テストシミュレーションによる結果（グラフ曲線）は、揺らぎ（上記のＤＳＭＣ法に伴って必然的に現れる揺らぎ）を含んでいるものの、マッハ数Ｍ、流体密度ρ及び絶対温度Ｔのいずれの物理量においても、理論値（等エントロピー流理論による予測値）とよく一致していることが分かる。このことから、モンテカルロ直接シミュレーション（ＤＳＭＣ法）は、本発明のノズル１０の有用性を検証するシミュレーション方法として、適切なものであることが確認された。
【００４５】
２．２．２  本シミュレーションの結果
続いて、本シミュレーションの結果について説明する。
【００４６】
図１０は、本シミュレーションにより得られた結果を示したグラフである。図１０における、マーカー「Ｅ１．１」、「Ｅ１．２」、「Ｅ１．３」、「Ｅ１．４」、「Ｅ１．５」、「Ｅ１．６」及び「Ｅ１．７」、並びに、「Ｅ２．１」、「Ｅ２．２」、「Ｅ２．３」、「Ｅ２．４」、「Ｅ２．５」、「Ｅ２．６」及び「Ｅ２．７」は、それぞれ、本発明に係る、実施例１．１、実施例１．２、実施例１．３、実施例１．４、実施例１．５、実施例１．６及び実施例１．７、並びに、実施例２．１、実施例２．２、実施例２．３、実施例２．４、実施例２．５、実施例２．６及び実施例２．７の結果を示したものであり、マーカー「Ｃ１」は、ラバールノズルに係る比較例１の結果を示したものであり、マーカー「Ｃ２」は、デュアルスロートノズルに係る比較例２の結果を示したものである。
【００４７】
また、図１０における横軸は、実施例１．１〜１．７及び実施例２．１〜２．７、並びに、比較例１及び比較例２のノズルにおける推力偏向角度δを、デュアルスロートノズルである比較例２（同図中の「Ｃ２」）の値を「１」として規格化して示したものである。さらに、図１０における縦軸は、実施例１．１〜１．７及び実施例２．１〜２．７、並びに、比較例１及び比較例２のノズルにおける推力の絶対値を、ラバールノズルである比較例１（同図中の「Ｃ１」）の値を「１」として規格化して示したものである。
【００４８】
図１０を見ると、本発明に係る実施例１．１〜１．７及び実施例２．１〜２．７の全てのノズルにおいて、推力ｆの絶対値は、１よりも大きくなっており、デュアルスロートノズルに係る比較例２（マーカーＣ２）の推力ｆの絶対値を上回っている。このことから、本発明の流体式推力方向制御装置のノズルのように、第二広がり部１７を設けることによって、従来から提案されているデュアルスロートノズルよりも、推力ｆを増大させることが可能であることが確認できた。
【００４９】
また、図１０を見ると、本発明に係る実施例１．１〜１．７及び実施例２．１〜２．７の多くのノズルにおいて、推力偏向角度δは、１よりも大きくなっており、ラバールノズルに係る比較例１（マーカーＣ１）の推力偏向角度δを上回っている。具体的には、第二広がり部１７の傾斜角度θ１が３０°未満の実施例（実施例１．１、実施例１．２、実施例２．１及び実施例２．２）を除いた全ての実施例（傾斜角度θ１が３０°以上の実施例１．３〜１．７及び実施例２．３〜２．７）において、推力偏向角度δが、比較例１（マーカーＣ１）の推力偏向角度δを上回っている。このことから、本発明の流体式推力方向制御装置のノズルのように、第二広がり部１７を設け、その傾斜角度θ１を適切な範囲（例えば、３０〜７０°の範囲）に設定することにより、従来のラバールノズルよりも、推力偏向角度δを増大させることが可能であることが確認できた。
【００５０】
ところで、図１０の結果からは、第二広がり部１７の長さＬ１が０．５９×Ｄ１と比較的短い実施例１．１〜１．７よりも、第二広がり部１７の長さＬ１が１．１８×Ｄ１と比較的長い実施例２．１〜２．７の方が、推力ｆの増大効果が高くなる傾向にあることも読み取れる。また、図１０の結果からは、第二広がり部１７の傾斜角度θ１が２０°を超える場合には、第二広がり部１７の傾斜角度θ１が小さいほど、推力ｆの増大効果が大きくなる傾向があることも読み取れる。さらに、図１０の結果からは、第二広がり部１７の傾斜角度θ１が大きいほど、推力偏向角度δが大きくなり、殆どのケースでは、第二広がり部１７の長さＬ１が短い場合に、推力偏向角度δがやや大きくなることも読み取れる。
【符号の説明】
【００５１】
１０      ノズル
１０ａ    ノズルの中心線
１０’    従来から提案されているラバールノズル
１０’’  従来から提案されているデュアルスロートノズル
１１      ノズル入口
１２      第一絞り部
１３      第一スロート
１４      第一広がり部
１５      第二絞り部
１６      第二スロート
１７      第二広がり部
１８      ノズル出口
１９      二次噴流注入口
１９ａ    下側の二次噴流注入口
１９ｂ    上側の二次噴流注入口
Ｄ１    ノズル入口の開口径
Ｌ１    第二広がり部のノズルの中心線に沿った方向の長さ
Ｆａ    代表的な流線
Ｆｂ    代表的な流線
Ｆｃ    代表的な流線
Ｆｄ    代表的な流線
ＦＩＮ１    エンジン燃焼室からノズル入口に流入する主流
ＦＩＮ２    二次噴流
ＦＯＵＴ    ノズル出口から噴出する噴流
α１    剥離領域
θ１    ノズルの中心線に対する第二広がり部の内壁の傾斜角度
【図１】
2018178970.tif 000006
【図２】
2018178970.tif 000007
【図３】
2018178970.tif 000008
【図４】
2018178970.tif 000009
【図５】
2018178970.tif 000010
【図６】
2018178970.tif 000011
【図７】
2018178970.tif 000012
【図８】
2018178970.tif 000013
【図９】
2018178970.tif 000014
【図１０】
2018178970.tif 000015