JPH10311773A - 衝撃波位置推定方法および衝撃波位置推定装置 - Google Patents
衝撃波位置推定方法および衝撃波位置推定装置Info
- Publication number
- JPH10311773A JPH10311773A JP13593997A JP13593997A JPH10311773A JP H10311773 A JPH10311773 A JP H10311773A JP 13593997 A JP13593997 A JP 13593997A JP 13593997 A JP13593997 A JP 13593997A JP H10311773 A JPH10311773 A JP H10311773A
- Authority
- JP
- Japan
- Prior art keywords
- shock wave
- pressure
- estimating
- measurement
- wave position
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Withdrawn
Links
Landscapes
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
(57)【要約】
【課題】 簡易なシステムにより衝撃波の移動量や位置
をリアルタイムに推定がなし得る衝撃波位置推定方法お
よび衝撃波位置推定装置を提供する。 【解決手段】 ラム燃焼器などへの空気導入路における
衝撃波の位置を空気導入路に一定間隔で設けられた所要
数の圧力検出手段１１を有する圧力測定部１０により圧
力分布を測定し、衝撃波位置推定部２０によりその測定
値の中で最低圧力が測定された位置近傍に衝撃波位置が
あると推定するものである。
をリアルタイムに推定がなし得る衝撃波位置推定方法お
よび衝撃波位置推定装置を提供する。 【解決手段】 ラム燃焼器などへの空気導入路における
衝撃波の位置を空気導入路に一定間隔で設けられた所要
数の圧力検出手段１１を有する圧力測定部１０により圧
力分布を測定し、衝撃波位置推定部２０によりその測定
値の中で最低圧力が測定された位置近傍に衝撃波位置が
あると推定するものである。
Description
【０００１】
【発明の属する技術分野】本発明は衝撃波位置推定方法
および衝撃波位置推定装置に関する。さらに詳しくは、
空気導入路内の圧力変化により衝撃波の位置を推定ある
いは検出する衝撃波位置推定方法および衝撃波位置推定
装置に関する。
および衝撃波位置推定装置に関する。さらに詳しくは、
空気導入路内の圧力変化により衝撃波の位置を推定ある
いは検出する衝撃波位置推定方法および衝撃波位置推定
装置に関する。
【０００２】
【従来の技術】従来より、衝撃波の位置を検出する方法
としてシュリーレン写真あるいはシュリーレン画像を利
用した方法が知られている。そして、この方法によれば
衝撃波の発生状況を画像として把握できるという利点が
ある反面、その移動量や位置を簡易に数値として把握す
ることはできないという欠点がある。
としてシュリーレン写真あるいはシュリーレン画像を利
用した方法が知られている。そして、この方法によれば
衝撃波の発生状況を画像として把握できるという利点が
ある反面、その移動量や位置を簡易に数値として把握す
ることはできないという欠点がある。
【０００３】そこで、シュリーレン画像に画像処理など
の再加工を施して衝撃波の移動量や位置を検出すること
が検討あるいは計画されている。
の再加工を施して衝撃波の移動量や位置を検出すること
が検討あるいは計画されている。
【０００４】しかしながら、前記画像処理などを用いた
方法には、画像処理に伴う時間遅れがあるためにリアル
タイム性に欠けるという問題、画像処理装置のためにシ
ステムが大掛かりとなるとともに複雑化するという問
題、画像処理に伴うデータの再加工により検出精度が低
下するという問題などがある。
方法には、画像処理に伴う時間遅れがあるためにリアル
タイム性に欠けるという問題、画像処理装置のためにシ
ステムが大掛かりとなるとともに複雑化するという問
題、画像処理に伴うデータの再加工により検出精度が低
下するという問題などがある。
【０００５】
【発明が解決しようとする課題】本発明はかかる従来技
術の課題に鑑みなされたものであって、簡易なシステム
により衝撃波の移動量や位置をリアルタイムに推定ある
いは検出がなし得る衝撃波位置推定方法および衝撃波位
置推定装置を提供することを目的としている。
術の課題に鑑みなされたものであって、簡易なシステム
により衝撃波の移動量や位置をリアルタイムに推定ある
いは検出がなし得る衝撃波位置推定方法および衝撃波位
置推定装置を提供することを目的としている。
【０００６】
【課題を解決するための手段】本発明の衝撃波位置推定
方法は、ラム燃焼器などへの空気導入路における衝撃波
の位置を推定する衝撃波位置推定方法であって、空気導
入路における圧力を適宜間隔または一定間隔で測定し、
最低圧力が得られた近傍、例えば圧力測定間隔の１／２
以内に衝撃波位置があると推定することを特徴とする。
方法は、ラム燃焼器などへの空気導入路における衝撃波
の位置を推定する衝撃波位置推定方法であって、空気導
入路における圧力を適宜間隔または一定間隔で測定し、
最低圧力が得られた近傍、例えば圧力測定間隔の１／２
以内に衝撃波位置があると推定することを特徴とする。
【０００７】さらに、本発明の衝撃波位置推定方法にお
いて圧力測定を一定間隔で行った場合は、具体的には、
下記のようにして衝撃波位置が推定される。
いて圧力測定を一定間隔で行った場合は、具体的には、
下記のようにして衝撃波位置が推定される。
【０００８】（１）第１推定方法 この第１推定方法は下記式により衝撃波位置を推定する
ものである。
ものである。
【０００９】X＝Xi-X0／２+X0＊（Pi-1-Pi）／（（Pi-1
-Pi）+（Pi+1-Ｐi））
-Pi）+（Pi+1-Ｐi））
【００１０】ここに、 Ｘ：推定衝撃波位置 Ｘi：最低圧力測定位置 Ｘ0：圧力測定間隔 Ｐi：測定最低圧力 Ｐi-1：最低圧力測定位置直前の測定圧力 Ｐi+1：最低圧力測定位置直後の測定圧力
【００１１】（２）第２推定方法 この第２推定方法は下記式により衝撃波位置を推定する
ものである。
ものである。
【００１２】 Ｘ＝Ｘi-1＋Ｘ0＊（Ｐi2−Ｐi）／（Ｐi2−Ｐi1）
【００１３】ここに、 Ｘ：推定衝撃波位置 Ｘi-1 ：最低圧力測定位置より一つ上流側の圧力測定位
置 Ｘ0：圧力測定間隔 Ｐi：測定最低圧力 Ｐi1：測定範囲の下流に衝撃波があるときの測定圧力 Ｐi2：測定位置Ｘiの直前の測定位置に衝撃波があると
きの測定圧力
置 Ｘ0：圧力測定間隔 Ｐi：測定最低圧力 Ｐi1：測定範囲の下流に衝撃波があるときの測定圧力 Ｐi2：測定位置Ｘiの直前の測定位置に衝撃波があると
きの測定圧力
【００１４】（３）第３推定方法 この第３推定方法は、空気導入路における圧力を一定間
隔で測定し、最低圧力が得られた位置の前後近傍におけ
る圧力分布を直線近似することにより、降下直線および
上昇直線を求め、その降下直線および上昇直線の交点に
衝撃波位置があるとするものである。
隔で測定し、最低圧力が得られた位置の前後近傍におけ
る圧力分布を直線近似することにより、降下直線および
上昇直線を求め、その降下直線および上昇直線の交点に
衝撃波位置があるとするものである。
【００１５】この第３推定方法においては、予め降下直
線を求めておいてもよく、また予め上昇直線の勾配を求
めておいてもよい。
線を求めておいてもよく、また予め上昇直線の勾配を求
めておいてもよい。
【００１６】一方、本発明の衝撃波位置推定装置の第１
形態は、ラム燃焼器などへの空気導入路における衝撃波
の位置を推定する衝撃波位置推定装置であって、空気導
入路に適宜間隔または一定間隔で設けられた所要数の圧
力検出手段を有する圧力測定部と、前記圧力測定部から
の測定値により衝撃波位置を推定する衝撃波位置推定部
とを備え、前記衝撃波位置推定部により最低測定圧力位
置近傍に衝撃波位置があると推定することを特徴とす
る。
形態は、ラム燃焼器などへの空気導入路における衝撃波
の位置を推定する衝撃波位置推定装置であって、空気導
入路に適宜間隔または一定間隔で設けられた所要数の圧
力検出手段を有する圧力測定部と、前記圧力測定部から
の測定値により衝撃波位置を推定する衝撃波位置推定部
とを備え、前記衝撃波位置推定部により最低測定圧力位
置近傍に衝撃波位置があると推定することを特徴とす
る。
【００１７】本発明の衝撃波位置推定装置の第２形態
は、ラム燃焼器などへの空気導入路における衝撃波の位
置を推定する衝撃波位置推定装置であって、空気導入路
に適宜間隔または一定間隔で設けられた所要数の圧力検
出手段を有する圧力測定部と、前記圧力測定部からの測
定値により衝撃波位置を推定する衝撃波位置推定部とを
備え、前記衝撃波位置推定部により圧力測定間隔の１／
２以内に衝撃波位置があると推定することを特徴とす
る。
は、ラム燃焼器などへの空気導入路における衝撃波の位
置を推定する衝撃波位置推定装置であって、空気導入路
に適宜間隔または一定間隔で設けられた所要数の圧力検
出手段を有する圧力測定部と、前記圧力測定部からの測
定値により衝撃波位置を推定する衝撃波位置推定部とを
備え、前記衝撃波位置推定部により圧力測定間隔の１／
２以内に衝撃波位置があると推定することを特徴とす
る。
【００１８】本発明の衝撃波位置推定装置の第３形態
は、ラム燃焼器などへの空気導入路における衝撃波の位
置を推定する衝撃波位置推定装置であって、空気導入路
に一定間隔で設けられた所要数の圧力検出手段を有する
圧力測定部と、前記圧力測定部からの測定値により衝撃
波位置を推定する衝撃波位置推定部とを備え、前記衝撃
波位置推定部により下記式により衝撃波位置を推定する
ことを特徴とする。
は、ラム燃焼器などへの空気導入路における衝撃波の位
置を推定する衝撃波位置推定装置であって、空気導入路
に一定間隔で設けられた所要数の圧力検出手段を有する
圧力測定部と、前記圧力測定部からの測定値により衝撃
波位置を推定する衝撃波位置推定部とを備え、前記衝撃
波位置推定部により下記式により衝撃波位置を推定する
ことを特徴とする。
【００１９】X＝Xi-X0／２+X0＊（Pi-1-Pi）／（（Pi-1
-Pi）+（Pi+1-Ｐi））
-Pi）+（Pi+1-Ｐi））
【００２０】ここに、 Ｘ：推定衝撃波位置 Ｘi ：最低圧力測定位置 Ｘ0：圧力測定間隔 Ｐi：測定最低圧力 Ｐi-1：最低圧力測定位置直前の測定圧力 Ｐi+1：最低圧力測定位置直後の測定圧力
【００２１】本発明の衝撃波位置推定装置の第４形態
は、ラム燃焼器などへの空気導入路における衝撃波の位
置を推定する衝撃波位置推定装置であって、空気導入路
に一定間隔で設けられた所要数の圧力検出手段を有する
圧力測定部と、前記圧力測定部からの測定値により衝撃
波位置を推定する衝撃波位置推定部とを備え、前記衝撃
波位置推定部により下記式により衝撃波位置を推定する
ことを特徴とする。
は、ラム燃焼器などへの空気導入路における衝撃波の位
置を推定する衝撃波位置推定装置であって、空気導入路
に一定間隔で設けられた所要数の圧力検出手段を有する
圧力測定部と、前記圧力測定部からの測定値により衝撃
波位置を推定する衝撃波位置推定部とを備え、前記衝撃
波位置推定部により下記式により衝撃波位置を推定する
ことを特徴とする。
【００２２】 Ｘ＝Ｘi-1＋Ｘ0＊（Ｐi2−Ｐi）／（Ｐi2−Ｐi1）
【００２３】ここに、 Ｘ：推定衝撃波位置 Ｘi-1 ：最低圧力測定位置より一つ上流側の圧力測定位
置 Ｘ0：圧力測定間隔 Ｐi：測定最低圧力 Ｐi1：測定範囲の下流に衝撃波があるときの測定圧力 Ｐi2：測定位置Ｘiの直前の測定位置に衝撃波があると
きの測定圧力
置 Ｘ0：圧力測定間隔 Ｐi：測定最低圧力 Ｐi1：測定範囲の下流に衝撃波があるときの測定圧力 Ｐi2：測定位置Ｘiの直前の測定位置に衝撃波があると
きの測定圧力
【００２４】本発明の衝撃波位置推定装置の第５形態
は、ラム燃焼器などへの空気導入路における衝撃波の位
置を推定する衝撃波位置推定装置であって、空気導入路
に一定間隔で設けられた所要数の圧力検出手段を有する
圧力測定部と、前記圧力測定部からの測定値により衝撃
波位置を推定する衝撃波位置推定部とを備え、前記衝撃
波位置推定部により、最低圧力が得られた位置の前後近
傍おける圧力分布を直線近似して降下直線および上昇直
線を算出し、その算出した降下直線および上昇直線の交
点を衝撃波位置であると推定することを特徴とする。
は、ラム燃焼器などへの空気導入路における衝撃波の位
置を推定する衝撃波位置推定装置であって、空気導入路
に一定間隔で設けられた所要数の圧力検出手段を有する
圧力測定部と、前記圧力測定部からの測定値により衝撃
波位置を推定する衝撃波位置推定部とを備え、前記衝撃
波位置推定部により、最低圧力が得られた位置の前後近
傍おける圧力分布を直線近似して降下直線および上昇直
線を算出し、その算出した降下直線および上昇直線の交
点を衝撃波位置であると推定することを特徴とする。
【００２５】本発明の衝撃波位置推定装置の第５形態に
おいては、予め降下直線を求めておいてもよく、また予
め上昇直線の勾配を求めておいてもよい。
おいては、予め降下直線を求めておいてもよく、また予
め上昇直線の勾配を求めておいてもよい。
【００２６】
【作用】本発明においては、ラム燃焼器などへの空気導
入路における圧力を、例えば一定間隔で測定し、そして
最低圧力が得られた測定位置の近傍に衝撃波位置がある
としているので、衝撃波位置を推定あるいは検出するた
めの装置構成が簡素化されるともに、衝撃波位置の推定
あるいは検出が短時間でなし得る。つまり、リアルタイ
ムで衝撃波位置の推定あるいは検出がなし得る。また、
この推定あるいは検出を連続して行うことにより、衝撃
波位置の移動も把握できる。
入路における圧力を、例えば一定間隔で測定し、そして
最低圧力が得られた測定位置の近傍に衝撃波位置がある
としているので、衝撃波位置を推定あるいは検出するた
めの装置構成が簡素化されるともに、衝撃波位置の推定
あるいは検出が短時間でなし得る。つまり、リアルタイ
ムで衝撃波位置の推定あるいは検出がなし得る。また、
この推定あるいは検出を連続して行うことにより、衝撃
波位置の移動も把握できる。
【００２７】
【発明の実施の形態】以下、添付図面を参照しながら本
発明を実施の形態に基づいて説明するが、本発明はかか
る実施の形態のみに限定されるものではない。
発明を実施の形態に基づいて説明するが、本発明はかか
る実施の形態のみに限定されるものではない。
【００２８】実施の形態１ 本発明の実施の形態１に係る衝撃波位置推定方法は、ラ
ムジェットのインテーク（以下、特に誤解の生じない場
合には単にインテークという）Ｉにおける衝撃波位置推
定方法であって、このインテークＩにおいては衝撃波の
前まで圧力が降下する一方、衝撃波の後では圧力が上昇
するという現像を利用するもので、図１に示すような所
要数の圧力センサ（以下、単にセンサという）１１をイ
ンテークＩ側壁に所定間隔で配置してなる圧力測定部１
０と、この圧力測定部１０からの測定値に基づいて衝撃
波位置を推定する衝撃波位置推定部２０と、出力部３０
とを備えた衝撃波位置推定装置Ａを用いて衝撃波位置を
推定するものであって、衝撃波位置推定部２０におい
て、圧力測定部１０の最低圧力を検出したセンサ１１L
の前後のセンサ１１B，１１Aの圧力と最低圧力との各差
圧により、そのセンサ１１Lの前後の１／２間を比例配
分して衝撃波の位置を推定あるいは検出するものであ
る。ここで、衝撃波位置推定部２０は、例えばマイクロ
コンピュータにより構成されており、そのＲＯＭには後
述する各種演算あるいは処理をするためのプログラムが
格納されている。また、そのＲＡＭには圧力測定値が一
時的に格納される。
ムジェットのインテーク（以下、特に誤解の生じない場
合には単にインテークという）Ｉにおける衝撃波位置推
定方法であって、このインテークＩにおいては衝撃波の
前まで圧力が降下する一方、衝撃波の後では圧力が上昇
するという現像を利用するもので、図１に示すような所
要数の圧力センサ（以下、単にセンサという）１１をイ
ンテークＩ側壁に所定間隔で配置してなる圧力測定部１
０と、この圧力測定部１０からの測定値に基づいて衝撃
波位置を推定する衝撃波位置推定部２０と、出力部３０
とを備えた衝撃波位置推定装置Ａを用いて衝撃波位置を
推定するものであって、衝撃波位置推定部２０におい
て、圧力測定部１０の最低圧力を検出したセンサ１１L
の前後のセンサ１１B，１１Aの圧力と最低圧力との各差
圧により、そのセンサ１１Lの前後の１／２間を比例配
分して衝撃波の位置を推定あるいは検出するものであ
る。ここで、衝撃波位置推定部２０は、例えばマイクロ
コンピュータにより構成されており、そのＲＯＭには後
述する各種演算あるいは処理をするためのプログラムが
格納されている。また、そのＲＡＭには圧力測定値が一
時的に格納される。
【００２９】例えば、インテークＩ側壁に配置されたセ
ンサ群により、図２に示すような圧力測定結果が得られ
たとすると、衝撃波位置推定部２０は、衝撃波位置は最
低圧力（図示例ではＰ3）が測定されたセンサ（図示例
ではＫ3）位置（図示例ではＸ3）の前後近傍にあると仮
定する。つまり、衝撃波位置Ｘはセンサ間隔をＸ0とす
ると、Ｘ3±Ｘ0／２の範囲にあると推定する。そして、
この±Ｘ0／２の範囲を最低圧力が検出されたセンサ１
１Lの前後のセンサ（図示例ではＫ2，Ｋ4）の測定圧力
（図示例ではＰ2，Ｐ4）とセンサＫ3の測定圧力Ｐ3との
差圧を用いて比例配分する。これを式で示せば下記式
（１）となる。
ンサ群により、図２に示すような圧力測定結果が得られ
たとすると、衝撃波位置推定部２０は、衝撃波位置は最
低圧力（図示例ではＰ3）が測定されたセンサ（図示例
ではＫ3）位置（図示例ではＸ3）の前後近傍にあると仮
定する。つまり、衝撃波位置Ｘはセンサ間隔をＸ0とす
ると、Ｘ3±Ｘ0／２の範囲にあると推定する。そして、
この±Ｘ0／２の範囲を最低圧力が検出されたセンサ１
１Lの前後のセンサ（図示例ではＫ2，Ｋ4）の測定圧力
（図示例ではＰ2，Ｐ4）とセンサＫ3の測定圧力Ｐ3との
差圧を用いて比例配分する。これを式で示せば下記式
（１）となる。
【００３０】 Ｘ＝Ｘ3−Ｘ0／２＋Ｘ0＊Ａ／（Ａ＋Ｂ） （１）
【００３１】ここに、 Ａ＝Ｐ2−Ｐ3 Ｂ＝Ｐ4−Ｐ3
【００３２】このように、この実施の形態１によれば、
インテークＩ側壁に所定間隔でセンサ群を設け、そのセ
ンサ群の中で最低圧力を測定したセンサ１１Lの前後間
をその前後のセンサ、つまり直前にあるセンサ（直前セ
ンサ）１１Bおよび直後にあるセンサ（直後センサ）１
１Aによって測定された圧力と最低圧力との差圧により
比例配分するという簡単な処理により衝撃波位置が推定
あるいは検出される。そのため、衝撃波位置をリアルタ
イムに推定あるいは検出できる。また、この推定あるい
は検出を連続的に行えば衝撃波位置の移動量を推定ある
いは検出することもできる。
インテークＩ側壁に所定間隔でセンサ群を設け、そのセ
ンサ群の中で最低圧力を測定したセンサ１１Lの前後間
をその前後のセンサ、つまり直前にあるセンサ（直前セ
ンサ）１１Bおよび直後にあるセンサ（直後センサ）１
１Aによって測定された圧力と最低圧力との差圧により
比例配分するという簡単な処理により衝撃波位置が推定
あるいは検出される。そのため、衝撃波位置をリアルタ
イムに推定あるいは検出できる。また、この推定あるい
は検出を連続的に行えば衝撃波位置の移動量を推定ある
いは検出することもできる。
【００３３】実施の形態２ 本発明の実施の形態２に係る衝撃波位置推定方法は、イ
ンテークＩにおける衝撃波位置推定方法であって、この
インテークＩにおいては衝撃波後の圧力勾配が衝撃波の
位置により異なるという現像を利用するもので、まずラ
ム燃焼器後方の排気ノズル（図示省略）を開閉し、イン
テークＩにおける衝撃波の発生位置を移動させてその移
動に伴う圧力分布を測定する。そして、衝撃波位置の推
定あるいは検出に際し、センサ群の測定圧力からどのセ
ンサ間に衝撃波があるかを仮定し、このときの圧力と先
に得た圧力とを差を用いて衝撃波位置のセンサ位置から
のずれを算出して衝撃波位置を推定あるいは検出するも
のである。なお、インテークＩの圧力の測定は、実施の
形態１と同様に、図１に示す衝撃波位置推定装置Ａによ
るものとされている。
ンテークＩにおける衝撃波位置推定方法であって、この
インテークＩにおいては衝撃波後の圧力勾配が衝撃波の
位置により異なるという現像を利用するもので、まずラ
ム燃焼器後方の排気ノズル（図示省略）を開閉し、イン
テークＩにおける衝撃波の発生位置を移動させてその移
動に伴う圧力分布を測定する。そして、衝撃波位置の推
定あるいは検出に際し、センサ群の測定圧力からどのセ
ンサ間に衝撃波があるかを仮定し、このときの圧力と先
に得た圧力とを差を用いて衝撃波位置のセンサ位置から
のずれを算出して衝撃波位置を推定あるいは検出するも
のである。なお、インテークＩの圧力の測定は、実施の
形態１と同様に、図１に示す衝撃波位置推定装置Ａによ
るものとされている。
【００３４】衝撃波が圧力測定部１０の位置検出範囲
（センサ設置範囲）の後方にある場合には、例えば図３
に示すような圧力分布が得られ、またラム燃焼器の排気
ノズル角をごくゆっくりと、つまり静止状態に近い状態
（準静的状態）で変化させて衝撃波を下流から上流に向
けてセンサ設置範囲全域にわたって移動させると、例え
ば図４に示すような圧力分布が得られる。そして、この
一連の準備が完了すると衝撃波位置の推定あるいは検出
を行う。この推定あるいは検出において、例えば図５に
示すような圧力分布が得られたとする。なお、図５には
説明の都合上、準備段階で得られた測定圧力も示されて
いる。
（センサ設置範囲）の後方にある場合には、例えば図３
に示すような圧力分布が得られ、またラム燃焼器の排気
ノズル角をごくゆっくりと、つまり静止状態に近い状態
（準静的状態）で変化させて衝撃波を下流から上流に向
けてセンサ設置範囲全域にわたって移動させると、例え
ば図４に示すような圧力分布が得られる。そして、この
一連の準備が完了すると衝撃波位置の推定あるいは検出
を行う。この推定あるいは検出において、例えば図５に
示すような圧力分布が得られたとする。なお、図５には
説明の都合上、準備段階で得られた測定圧力も示されて
いる。
【００３５】しかして、衝撃波位置推定部２０は、図５
に示す圧力分布から衝撃波の位置は、センサＫ2とセン
サＫ3との間にあると仮定する。そこで、衝撃波位置の
推定あるいは検出時におけるセンサＫ3の測定圧力Ｐ3と
衝撃波がセンサＫ2の位置Ｘ2にあるときのセンサＫ3の
測定圧力Ｐ32との差圧および衝撃波がセンサ設置範囲の
下流にあるときのセンサＫ3の測定圧力Ｐ31と前記圧力
Ｐ32との差圧の比により、センサＫ2とセンサＫ3との間
隔Ｘ0を比例配分して衝撃波位置Ｘを推定あるいは検出
する。これを式で示せば下記式（２）となる。
に示す圧力分布から衝撃波の位置は、センサＫ2とセン
サＫ3との間にあると仮定する。そこで、衝撃波位置の
推定あるいは検出時におけるセンサＫ3の測定圧力Ｐ3と
衝撃波がセンサＫ2の位置Ｘ2にあるときのセンサＫ3の
測定圧力Ｐ32との差圧および衝撃波がセンサ設置範囲の
下流にあるときのセンサＫ3の測定圧力Ｐ31と前記圧力
Ｐ32との差圧の比により、センサＫ2とセンサＫ3との間
隔Ｘ0を比例配分して衝撃波位置Ｘを推定あるいは検出
する。これを式で示せば下記式（２）となる。
【００３６】 Ｘ＝Ｘ2＋Ｘ0＊Ｄ／Ｃ （２）
【００３７】ここに、 Ｃ＝Ｐ32−Ｐ3 Ｄ＝Ｐ32−Ｐ31
【００３８】このように、この実施の形態２によれば、
インテークＩ側壁に所定間隔でセンサ群を設け、それに
より衝撃波の下流の圧力勾配を求め、その得られた圧力
勾配を利用して衝撃波位置を推定あるいは検出している
ので、実施の形態１と同様に簡易に衝撃波位置をリアル
タイムに推定あるいは検出できる。また、この推定ある
いは検出を連続的に行えば衝撃波位置の移動量を推定あ
るいは検出することもできる。
インテークＩ側壁に所定間隔でセンサ群を設け、それに
より衝撃波の下流の圧力勾配を求め、その得られた圧力
勾配を利用して衝撃波位置を推定あるいは検出している
ので、実施の形態１と同様に簡易に衝撃波位置をリアル
タイムに推定あるいは検出できる。また、この推定ある
いは検出を連続的に行えば衝撃波位置の移動量を推定あ
るいは検出することもできる。
【００３９】実施の形態３ 本発明の実施の形態３に係る衝撃波位置推定方法は、イ
ンテークＩのセンサ群により測定された衝撃波前後の圧
力分布を局部的に見れば直線近似が可能であるところか
ら、まずセンサ群の中で最低圧力を示したセンサ１１L
の近傍に衝撃波があると仮定し、ついでその仮定された
衝撃波前後の近傍における圧力勾配を直線近似し、その
直線近似により得られた二直線の交点を求めてその交点
を衝撃波位置とするものである。ここで、圧力勾配を直
線近似するために採用する測定圧力は、仮定された衝撃
波位置の前の数個あるいは後の数個のセンサ１１の測定
圧力とされる。なお、インテークＩの圧力分布の測定
は、実施の形態１と同様に、図１に示す衝撃波位置推定
装置Ａによるものとされている。
ンテークＩのセンサ群により測定された衝撃波前後の圧
力分布を局部的に見れば直線近似が可能であるところか
ら、まずセンサ群の中で最低圧力を示したセンサ１１L
の近傍に衝撃波があると仮定し、ついでその仮定された
衝撃波前後の近傍における圧力勾配を直線近似し、その
直線近似により得られた二直線の交点を求めてその交点
を衝撃波位置とするものである。ここで、圧力勾配を直
線近似するために採用する測定圧力は、仮定された衝撃
波位置の前の数個あるいは後の数個のセンサ１１の測定
圧力とされる。なお、インテークＩの圧力分布の測定
は、実施の形態１と同様に、図１に示す衝撃波位置推定
装置Ａによるものとされている。
【００４０】例えば、衝撃波位置推定部２０において、
図６に示すように仮定された衝撃波位置よりも上流側近
傍にあるセンサ１１により得られた圧力分布を直線近似
することにより直線Ｌ1が得られたとし、また前記衝撃
波位置よりも下流側近傍にあるセンサ１１により得られ
た圧力分布を直線近似することにより直線Ｌ2が得られ
たとすると、直線Ｌ1と直線Ｌ2との交点が衝撃波位置を
与えるとするものである。
図６に示すように仮定された衝撃波位置よりも上流側近
傍にあるセンサ１１により得られた圧力分布を直線近似
することにより直線Ｌ1が得られたとし、また前記衝撃
波位置よりも下流側近傍にあるセンサ１１により得られ
た圧力分布を直線近似することにより直線Ｌ2が得られ
たとすると、直線Ｌ1と直線Ｌ2との交点が衝撃波位置を
与えるとするものである。
【００４１】このように、この実施の形態３によれば、
インテークＩ側壁に所定間隔でセンサ群を設け、それら
のセンサ群により測定された圧力分布より衝撃波位置を
仮定するとともに、その仮定された衝撃波の前後数個の
センサ１１の位置および測定値により、衝撃波前後近傍
における圧力勾配を近似する直線を算出し、それら２直
線の交点を衝撃波位置としているので、簡易かつ迅速、
つまりリアルタイムに衝撃波位置が推定あるいは検出さ
れる。また、この推定あるいは検出を連続的に行えば衝
撃波位置の移動量を推定あるいは検出することもでき
る。
インテークＩ側壁に所定間隔でセンサ群を設け、それら
のセンサ群により測定された圧力分布より衝撃波位置を
仮定するとともに、その仮定された衝撃波の前後数個の
センサ１１の位置および測定値により、衝撃波前後近傍
における圧力勾配を近似する直線を算出し、それら２直
線の交点を衝撃波位置としているので、簡易かつ迅速、
つまりリアルタイムに衝撃波位置が推定あるいは検出さ
れる。また、この推定あるいは検出を連続的に行えば衝
撃波位置の移動量を推定あるいは検出することもでき
る。
【００４２】実施の形態４ 本発明の実施の形態４に係る衝撃波位置推定方法は、実
施の形態３を改変したものであって、事前にセンサ１１
が設置されている全域にわたって衝撃波前の圧力分布を
測定しておいて、その圧力分布より衝撃波前の圧力勾配
を近似する直線を算出しておき、衝撃波位置の検出時に
仮定された衝撃波位置の下流近傍の測定圧力より圧力勾
配を求めるとともに、その圧力勾配を近似する直線を実
施の形態３と同様にして算出し、その算出された直線と
前記衝撃波前の直線との交点を算出し、その交点を衝撃
波位置とするものである。
施の形態３を改変したものであって、事前にセンサ１１
が設置されている全域にわたって衝撃波前の圧力分布を
測定しておいて、その圧力分布より衝撃波前の圧力勾配
を近似する直線を算出しておき、衝撃波位置の検出時に
仮定された衝撃波位置の下流近傍の測定圧力より圧力勾
配を求めるとともに、その圧力勾配を近似する直線を実
施の形態３と同様にして算出し、その算出された直線と
前記衝撃波前の直線との交点を算出し、その交点を衝撃
波位置とするものである。
【００４３】例えば、衝撃波位置推定部２０において、
図７に示すように事前に衝撃波前の圧力勾配を近似する
直線Ｌ1を求めておいて、衝撃波位置の検出時に最低圧
力を測定しているセンサＫ3の位置Ｘ3と圧力Ｐ3および
センサＫ4の位置Ｘ4と圧力Ｐ4により衝撃波後の圧力勾
配を近似する直線Ｌ2を算出し、そしてこれらの２直線
Ｌ1，Ｌ2の交点を算出してその交点を衝撃波位置とする
ものである。
図７に示すように事前に衝撃波前の圧力勾配を近似する
直線Ｌ1を求めておいて、衝撃波位置の検出時に最低圧
力を測定しているセンサＫ3の位置Ｘ3と圧力Ｐ3および
センサＫ4の位置Ｘ4と圧力Ｐ4により衝撃波後の圧力勾
配を近似する直線Ｌ2を算出し、そしてこれらの２直線
Ｌ1，Ｌ2の交点を算出してその交点を衝撃波位置とする
ものである。
【００４４】このように、この実施の形態４によれば、
衝撃波前の圧力勾配を近似する直線が事前に算出されて
いるので、より短時間で衝撃波位置の推定あるいは検出
がなし得る。
衝撃波前の圧力勾配を近似する直線が事前に算出されて
いるので、より短時間で衝撃波位置の推定あるいは検出
がなし得る。
【００４５】実施の形態５ 本発明の実施の形態５に係る衝撃波位置推定方法は、実
施の形態４を改変したものであって、事前にセンサ１１
が設置されている全域にわたって衝撃波前の圧力分布を
測定しておいて、その圧力分布より衝撃波前の圧力勾配
を近似する直線を算出するとともに、排気ノズル角を準
静的に変化させて衝撃波の位置を上流側あるいは下流側
に移動させながら、そのときの衝撃波後の圧力勾配を測
定し、その圧力勾配を近似する直線の勾配を算出してお
き、衝撃波位置の検出時にその直線の勾配と最低圧力を
測定したセンサ１１Lの一つ後のセンサ１１Aの位置と測
定圧力により、衝撃波後の圧力勾配を近似する直線を算
出し、その算出された直線と前記衝撃波前の直線との交
点を求めて、その交点を衝撃波位置とするものである。
施の形態４を改変したものであって、事前にセンサ１１
が設置されている全域にわたって衝撃波前の圧力分布を
測定しておいて、その圧力分布より衝撃波前の圧力勾配
を近似する直線を算出するとともに、排気ノズル角を準
静的に変化させて衝撃波の位置を上流側あるいは下流側
に移動させながら、そのときの衝撃波後の圧力勾配を測
定し、その圧力勾配を近似する直線の勾配を算出してお
き、衝撃波位置の検出時にその直線の勾配と最低圧力を
測定したセンサ１１Lの一つ後のセンサ１１Aの位置と測
定圧力により、衝撃波後の圧力勾配を近似する直線を算
出し、その算出された直線と前記衝撃波前の直線との交
点を求めて、その交点を衝撃波位置とするものである。
【００４６】例えば、衝撃波位置推定部２０において、
図８に示すように事前に衝撃波前の圧力勾配を近似する
直線Ｌ1を求めるとともに、衝撃波後の圧力勾配を近似
する直線の勾配を求めておいて、衝撃波位置の検出時に
その直線の勾配と最低圧力を測定しているセンサＫ3の
次のセンサＫ4の位置Ｘ4と圧力Ｐ4により衝撃波後の圧
力勾配を近似する直線Ｌ2を算出し、そしてこれらの２
直線Ｌ1，Ｌ2の交点を算出してその交点を衝撃波位置と
するものである。
図８に示すように事前に衝撃波前の圧力勾配を近似する
直線Ｌ1を求めるとともに、衝撃波後の圧力勾配を近似
する直線の勾配を求めておいて、衝撃波位置の検出時に
その直線の勾配と最低圧力を測定しているセンサＫ3の
次のセンサＫ4の位置Ｘ4と圧力Ｐ4により衝撃波後の圧
力勾配を近似する直線Ｌ2を算出し、そしてこれらの２
直線Ｌ1，Ｌ2の交点を算出してその交点を衝撃波位置と
するものである。
【００４７】このように、この実施の形態５によれば、
衝撃波前の圧力勾配を近似する直線および衝撃波後の圧
力勾配を近似する直線の勾配が事前に算出されているの
で、より短時間で衝撃波位置の推定あるいは検出がなし
得る。
衝撃波前の圧力勾配を近似する直線および衝撃波後の圧
力勾配を近似する直線の勾配が事前に算出されているの
で、より短時間で衝撃波位置の推定あるいは検出がなし
得る。
【００４８】
【実施例】以下、より具体的な実施例により本発明をよ
り具体的に説明する。
り具体的に説明する。
【００４９】実施例１〜２および比較例１〜２ 図９に示すダミーインテークＩD、ラム燃焼器ＲD、可変
排気ノズルＮDとを備えた試験装置Ｔに、図１に示す衝
撃波位置推定装置Ａを装着し、そして実施の形態１によ
り衝撃波位置の推定あるいは検出を行った。衝撃波位置
の推定あるいは検出は、計測開始後１９．１秒後（実施
例１）および計測開始後２３．１後（実施例２）につい
て行い、その結果を表１に示す。また、比較のために実
施例１および実施例２と同一時刻におけるシュリーレン
画像により衝撃波位置を計測し（比較例１および比較例
２）、その結果を表１に併せて示した。なお、このシュ
リーレン画像による衝撃波位置の計測に当たっては、圧
力センサが試験装置ＴのダミーインテークＩDの上側に
設けられているところから、衝撃波がダミーインテーク
ＩDの上側に到達している個所を衝撃波位置とした。
排気ノズルＮDとを備えた試験装置Ｔに、図１に示す衝
撃波位置推定装置Ａを装着し、そして実施の形態１によ
り衝撃波位置の推定あるいは検出を行った。衝撃波位置
の推定あるいは検出は、計測開始後１９．１秒後（実施
例１）および計測開始後２３．１後（実施例２）につい
て行い、その結果を表１に示す。また、比較のために実
施例１および実施例２と同一時刻におけるシュリーレン
画像により衝撃波位置を計測し（比較例１および比較例
２）、その結果を表１に併せて示した。なお、このシュ
リーレン画像による衝撃波位置の計測に当たっては、圧
力センサが試験装置ＴのダミーインテークＩDの上側に
設けられているところから、衝撃波がダミーインテーク
ＩDの上側に到達している個所を衝撃波位置とした。
【００５０】
【表１】
【００５１】表１より実施例１および実施例２における
衝撃波の推定あるいは検出位置は、それぞれ比較例１お
よび比較例２の検出位置によく一致しているのがわか
る。とりわけ、衝撃波がスロート近傍に位置するときに
は、シュリーレン画像による計測値と本発明の衝撃波位
置推定方法による推定値とは完全に一致している。した
がって、従来のシュリーレン画像による検出方法を本発
明の推定方法に置き換えても問題のないのがわかる。
衝撃波の推定あるいは検出位置は、それぞれ比較例１お
よび比較例２の検出位置によく一致しているのがわか
る。とりわけ、衝撃波がスロート近傍に位置するときに
は、シュリーレン画像による計測値と本発明の衝撃波位
置推定方法による推定値とは完全に一致している。した
がって、従来のシュリーレン画像による検出方法を本発
明の推定方法に置き換えても問題のないのがわかる。
【００５２】以上、本発明を実施の形態および実施例に
基づいて説明してきたが、本発明はかかる実施の形態お
よび実施例のみに限定されるものではなく、種々改変が
可能である。例えば、本発明はラムジェットエンジン、
ラム燃焼器試験設備、ショックダウン試験設備、衝撃波
位置検出をなすための設備などに適用できる。また、実
施の形態においては圧力センサは所定間隔で配置されて
いるが、適宜間隔とされてもよい。例えば、精度よく推
定したい部分には密に配置し、残りの部分については粗
く配置するようにしてもよい。
基づいて説明してきたが、本発明はかかる実施の形態お
よび実施例のみに限定されるものではなく、種々改変が
可能である。例えば、本発明はラムジェットエンジン、
ラム燃焼器試験設備、ショックダウン試験設備、衝撃波
位置検出をなすための設備などに適用できる。また、実
施の形態においては圧力センサは所定間隔で配置されて
いるが、適宜間隔とされてもよい。例えば、精度よく推
定したい部分には密に配置し、残りの部分については粗
く配置するようにしてもよい。
【００５３】
【発明の効果】以上詳述したように本発明によれば、ラ
ム燃焼器などの空気導入路の側壁に適宜間隔または所定
間隔で圧力検出手段を設けて空気導入路内の圧力分布を
測定し、そして最低測定圧力が得られた位置の近傍に衝
撃波位置があるとしているので、衝撃波位置を簡易な構
成により迅速に推定あるいは検出できるという優れた効
果が得られる。つまり、衝撃波位置をリアルタイムに推
定あるいは検出できるという優れた効果が得られる。ま
た、本発明による衝撃波位置の推定あるいは検出を連続
して行えば簡易に衝撃波位置の移動も把握できるという
効果も得られる。
ム燃焼器などの空気導入路の側壁に適宜間隔または所定
間隔で圧力検出手段を設けて空気導入路内の圧力分布を
測定し、そして最低測定圧力が得られた位置の近傍に衝
撃波位置があるとしているので、衝撃波位置を簡易な構
成により迅速に推定あるいは検出できるという優れた効
果が得られる。つまり、衝撃波位置をリアルタイムに推
定あるいは検出できるという優れた効果が得られる。ま
た、本発明による衝撃波位置の推定あるいは検出を連続
して行えば簡易に衝撃波位置の移動も把握できるという
効果も得られる。
【図１】本発明の衝撃波位置推定方法に用いられる衝撃
波位置推定装置の一例の概略図である。
波位置推定装置の一例の概略図である。
【図２】本発明の実施の形態１に係る衝撃波位置推定方
法の説明図である。
法の説明図である。
【図３】本発明の実施の形態２の推定方法の説明図であ
って、圧力測定範囲の下流に衝撃波があるときの圧力分
布を模式図的に示すグラフである。
って、圧力測定範囲の下流に衝撃波があるときの圧力分
布を模式図的に示すグラフである。
【図４】本発明の実施の形態２の推定方法の説明図であ
って、センサＫ2の位置に衝撃波があるときの圧力分布
を模式図的に示すグラフである。
って、センサＫ2の位置に衝撃波があるときの圧力分布
を模式図的に示すグラフである。
【図５】本発明の実施の形態２の衝撃波位置推定方法の
説明図である。
説明図である。
【図６】本発明の実施の形態３の衝撃波位置推定方法の
説明図である。
説明図である。
【図７】本発明の実施の形態４の衝撃波位置推定方法の
説明図である。
説明図である。
【図８】本発明の実施の形態５の衝撃波位置推定方法の
説明図である。
説明図である。
【図９】本発明の実施例および比較例において用いられ
た試験設備の側面図である。
た試験設備の側面図である。
【図１０】比較例１のシュリーレン画像である。
【図１１】比較例２のシュリーレン画像である。
１０ 圧力測定部 １１ 圧力センサ ２０ 衝撃波位置推定部 ３０ 出力部 Ａ 衝撃波位置推定装置 Ｉ インテーク Ｔ 試験装置 ＩD ダミーインテーク ＲD ラム燃焼器 ＮD 可変排気ノズル
───────────────────────────────────────────────────── フロントページの続き (72)発明者 神戸 勝啓 明石市川崎町１番１号 川崎重工業株式会 社明石工場内 (72)発明者 二村 尚夫 明石市川崎町１番１号 川崎重工業株式会 社明石工場内 (72)発明者 大嶋 孝夫 明石市川崎町１番１号 川崎重工業株式会 社明石工場内 (72)発明者 榎本 吉也 明石市川崎町１番１号 川崎重工業株式会 社明石工場内 (72)発明者 唄野 真 明石市川崎町１番１号 川崎重工業株式会 社明石工場内
Claims (15)
- 【請求項１】 ラム燃焼器などへの空気導入路における
衝撃波の位置を推定する衝撃波位置推定方法であって、 空気導入路における圧力を適宜間隔で測定し、最低圧力
が得られた近傍に衝撃波位置があると推定することを特
徴とする衝撃波位置推定方法。 - 【請求項２】 ラム燃焼器などへの空気導入路における
衝撃波の位置を推定する衝撃波位置推定方法であって、 空気導入路における圧力を一定間隔で測定し、最低圧力
が得られた近傍に衝撃波位置があると推定することを特
徴とする衝撃波位置推定方法。 - 【請求項３】 前記近傍が圧力測定間隔の１／２以内と
されていることを特徴とする請求項１または２記載の衝
撃波位置推定方法。 - 【請求項４】 衝撃波位置を下記式により推定すること
を特徴とする請求項２記載の衝撃波位置推定方法。 X＝Xi-X0／２+X0＊（Pi-1-Pi）／（（Pi-1-Pi）+（Pi+1
-Ｐi）） ここに、 Ｘ：推定衝撃波位置 Ｘi：最低圧力測定位置 Ｘ0：圧力測定間隔 Ｐi：測定最低圧力 Ｐi-1：最低圧力測定位置直前の測定圧力 Ｐi+1：最低圧力測定位置直後の測定圧力 - 【請求項５】 衝撃波位置を下記式により推定すること
を特徴とする請求項２記載の衝撃波位置推定方法。 Ｘ＝Ｘi-1＋Ｘ0＊（Ｐi2−Ｐi）／（Ｐi2−Ｐi1） ここに、 Ｘ：推定衝撃波位置 Ｘi-1 ：最低圧力測定位置より一つ上流側の圧力測定位
置 Ｘ0：圧力測定間隔 Ｐi：測定最低圧力 Ｐi1：測定範囲の下流に衝撃波があるときの測定圧力 Ｐi2：測定位置Ｘiの直前の測定位置に衝撃波があると
きの測定圧力 - 【請求項６】 ラム燃焼器などへの空気導入路における
衝撃波の位置を推定する衝撃波位置推定方法であって、 空気導入路における圧力を一定間隔で測定し、最低圧力
が得られた位置の前後近傍における圧力分布を直線近似
することにより、降下直線および上昇直線を求め、その
降下直線および上昇直線の交点を衝撃波位置とすること
を特徴とする請求項２記載の衝撃波位置推定方法。 - 【請求項７】 予め降下直線を求めておくことを特徴と
する請求項６記載の衝撃波位置推定方法。 - 【請求項８】 予め上昇直線の勾配を求めておくことを
特徴とする請求項６または７記載の衝撃波位置推定方
法。 - 【請求項９】 ラム燃焼器などへの空気導入路における
衝撃波の位置を推定する衝撃波位置推定装置であって、
空気導入路に適宜間隔または一定間隔で設けられた所要
数の圧力検出手段を有する圧力測定部と、前記圧力測定
部からの測定値により衝撃波位置を推定する衝撃波位置
推定部とを備え、 前記衝撃波位置推定部により最低測定圧力位置近傍に衝
撃波位置があると推定することを特徴とする衝撃波位置
推定装置。 - 【請求項１０】 ラム燃焼器などへの空気導入路におけ
る衝撃波の位置を推定する衝撃波位置推定装置であっ
て、空気導入路に適宜間隔または一定間隔で設けられた
所要数の圧力検出手段を有する圧力測定部と、前記圧力
測定部からの測定値により衝撃波位置を推定する衝撃波
位置推定部とを備え、 前記衝撃波位置推定部により圧力測定間隔の１／２以内
に衝撃波位置があると推定することを特徴とする衝撃波
位置推定装置。 - 【請求項１１】 ラム燃焼器などへの空気導入路におけ
る衝撃波の位置を推定する衝撃波位置推定装置であっ
て、空気導入路に一定間隔で設けられた所要数の圧力検
出手段を有する圧力測定部と、前記圧力測定部からの測
定値により衝撃波位置を推定する衝撃波位置推定部とを
備え、 前記衝撃波位置推定部により下記式により衝撃波位置を
推定することを特徴とする衝撃波位置推定装置。 X＝Xi-X0／２+X0＊（Pi-1-Pi）／（（Pi-1-Pi）+（Pi+1
-Ｐi）） ここに、 Ｘ：推定衝撃波位置 Ｘ ：最低圧力測定位置 Ｘ0：圧力測定間隔 Ｐi：測定最低圧力 Ｐi-1：最低圧力測定位置直前の測定圧力 Ｐi+1：最低圧力測定位置直後の測定圧力 - 【請求項１２】 ラム燃焼器などへの空気導入路におけ
る衝撃波の位置を推定する衝撃波位置推定装置であっ
て、空気導入路に一定間隔で設けられた所要数の圧力検
出手段を有する圧力測定部と、前記圧力測定部からの測
定値により衝撃波位置を推定する衝撃波位置推定部とを
備え、 前記衝撃波位置推定部により下記式により衝撃波位置を
推定することを特徴とする衝撃波位置推定装置。 Ｘ＝Ｘi-1＋Ｘ0＊（Ｐi2−Ｐi）／（Ｐi2−Ｐi1） ここに、 Ｘ：推定衝撃波位置 Ｘi-1 ：最低圧力測定位置より一つ上流側の圧力測定位
置 Ｘ0：圧力測定間隔 Ｐi：測定最低圧力 Ｐi1：測定範囲の下流に衝撃波があるときの測定圧力 Ｐi2：測定位置Ｘiの直前の測定位置に衝撃波があると
きの測定圧力 - 【請求項１３】 ラム燃焼器などへの空気導入路におけ
る衝撃波の位置を推定する衝撃波位置推定装置であっ
て、空気導入路に一定間隔で設けられた所要数の圧力検
出手段を有する圧力測定部と、前記圧力測定部からの測
定値により衝撃波位置を推定する衝撃波位置推定部とを
備え、 前記衝撃波位置推定部により、最低圧力が得られた位置
の前後近傍おける圧力分布を直線近似して降下直線およ
び上昇直線を算出し、その算出した降下直線および上昇
直線の交点を衝撃波位置であると推定することを特徴と
する衝撃波位置推定装置。 - 【請求項１４】 予め降下直線を求めておくことを特徴
とする請求項１３記載の衝撃波位置推定装置。 - 【請求項１５】 予め上昇直線の勾配を求めておくこと
を特徴とする請求項１３または１４記載の衝撃波位置推
定装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP13593997A JPH10311773A (ja) | 1997-05-09 | 1997-05-09 | 衝撃波位置推定方法および衝撃波位置推定装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP13593997A JPH10311773A (ja) | 1997-05-09 | 1997-05-09 | 衝撃波位置推定方法および衝撃波位置推定装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH10311773A true JPH10311773A (ja) | 1998-11-24 |
Family
ID=15163385
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP13593997A Withdrawn JPH10311773A (ja) | 1997-05-09 | 1997-05-09 | 衝撃波位置推定方法および衝撃波位置推定装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH10311773A (ja) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8041957B2 (en) | 2003-04-08 | 2011-10-18 | Qualcomm Incorporated | Associating software with hardware using cryptography |
| US8499171B2 (en) | 2005-11-18 | 2013-07-30 | Qualcomm Incorporated | Mobile security system and method |
| CN105738067A (zh) * | 2016-02-01 | 2016-07-06 | 南京航空航天大学 | 一种两道同侧斜激波相交后参数的快速确定方法 |
| CN105841882A (zh) * | 2016-03-17 | 2016-08-10 | 中国人民解放军国防科学技术大学 | 一种隔离段内激波串前缘探测装置 |
| CN107576446A (zh) * | 2017-08-02 | 2018-01-12 | 南京航空航天大学 | 一种超燃冲压发动机隔离段激波串前缘的主动探测装置及方法 |
-
1997
- 1997-05-09 JP JP13593997A patent/JPH10311773A/ja not_active Withdrawn
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8041957B2 (en) | 2003-04-08 | 2011-10-18 | Qualcomm Incorporated | Associating software with hardware using cryptography |
| US8499171B2 (en) | 2005-11-18 | 2013-07-30 | Qualcomm Incorporated | Mobile security system and method |
| CN105738067A (zh) * | 2016-02-01 | 2016-07-06 | 南京航空航天大学 | 一种两道同侧斜激波相交后参数的快速确定方法 |
| CN105738067B (zh) * | 2016-02-01 | 2018-04-06 | 南京航空航天大学 | 一种两道同侧斜激波相交后参数的快速确定方法 |
| CN105841882A (zh) * | 2016-03-17 | 2016-08-10 | 中国人民解放军国防科学技术大学 | 一种隔离段内激波串前缘探测装置 |
| CN105841882B (zh) * | 2016-03-17 | 2018-08-03 | 中国人民解放军国防科学技术大学 | 一种隔离段内激波串前缘探测装置 |
| CN107576446A (zh) * | 2017-08-02 | 2018-01-12 | 南京航空航天大学 | 一种超燃冲压发动机隔离段激波串前缘的主动探测装置及方法 |
| CN107576446B (zh) * | 2017-08-02 | 2020-10-20 | 南京航空航天大学 | 一种超燃冲压发动机隔离段激波串前缘的主动探测方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104021676B (zh) | 基于车辆动态视频特征的车辆定位及车速测量方法 | |
| US7519483B2 (en) | Unsteady flow meter | |
| KR970012113A (ko) | 좌표 입력 장치 및 방법 | |
| Costantini et al. | Feasibility of skin-friction field measurements in a transonic wind tunnel using a global luminescent oil film | |
| ATE211817T1 (de) | Vorrichtung zur volumenmessung strömender medien sowie entsprechendes verfahren | |
| JPH10311773A (ja) | 衝撃波位置推定方法および衝撃波位置推定装置 | |
| Schwieger | Sensitivity analysis as a general tool for model optimisation–examples for trajectory estimation | |
| JPS63132107A (ja) | 光切断線法による物体測定装置 | |
| Robbins | Detailed Characterization of Flowfields and Uncertainty in a Speed-Bump Turbulent Separated Flow Validation Experiment | |
| Naughton et al. | Skin-friction distribution near a cylinder mounted on a flat plate | |
| Beltrame | Accuracy assessment of the SU2 flow solver for non-ideal organic vapor supersonic expansions using experimental data | |
| JP2830810B2 (ja) | 追尾処理方法及び装置 | |
| JPH07214124A (ja) | 圧延機制御装置 | |
| RU1805293C (ru) | Способ измерени геометрических размеров объекта и устройство дл его осуществлени | |
| JPS6011172A (ja) | 測定出力のピ−ク検出回路 | |
| JPH0882663A (ja) | 目標運動解析方法 | |
| Kashinskii et al. | Phase velocity in a gas-liquid bubble flow | |
| JPS60165555A (ja) | 物体まわりの風圧測定装置 | |
| JPS63266318A (ja) | 異音発生箇所探査装置 | |
| JPH09319404A (ja) | 自己調整型制御装置及び方法 | |
| JPH06213700A (ja) | 流量計器差試験装置 | |
| JPH07270147A (ja) | 放射線厚さ計の演算処理装置 | |
| JPS63249069A (ja) | 異音発生箇所探査装置 | |
| JPH05342486A (ja) | ガス漏洩源の推定方法及びその装置 | |
| JPH06281429A (ja) | 形鋼材の寸法測定装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A300 | Withdrawal of application because of no request for examination |
Free format text: JAPANESE INTERMEDIATE CODE: A300 Effective date: 20040803 |