JP7465483B2 - 隆起構造および物体 - Google Patents
隆起構造および物体 Download PDFInfo
- Publication number
- JP7465483B2 JP7465483B2 JP2021511412A JP2021511412A JP7465483B2 JP 7465483 B2 JP7465483 B2 JP 7465483B2 JP 2021511412 A JP2021511412 A JP 2021511412A JP 2021511412 A JP2021511412 A JP 2021511412A JP 7465483 B2 JP7465483 B2 JP 7465483B2
- Authority
- JP
- Japan
- Prior art keywords
- flow
- raised structure
- protruding
- cross
- height
- 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.)
- Active
Links
- 230000007704 transition Effects 0.000 claims description 46
- 239000012530 fluid Substances 0.000 claims description 11
- 230000008859 change Effects 0.000 claims description 7
- 238000011144 upstream manufacturing Methods 0.000 claims description 6
- 230000001154 acute effect Effects 0.000 claims description 4
- 238000000034 method Methods 0.000 description 14
- 230000000052 comparative effect Effects 0.000 description 12
- 230000000694 effects Effects 0.000 description 10
- 238000004458 analytical method Methods 0.000 description 9
- 238000009826 distribution Methods 0.000 description 6
- 230000000737 periodic effect Effects 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 230000003321 amplification Effects 0.000 description 4
- 230000007423 decrease Effects 0.000 description 4
- 230000005284 excitation Effects 0.000 description 4
- 230000001976 improved effect Effects 0.000 description 4
- 239000000463 material Substances 0.000 description 4
- 238000000465 moulding Methods 0.000 description 4
- 238000003199 nucleic acid amplification method Methods 0.000 description 4
- 229920006395 saturated elastomer Polymers 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 230000007935 neutral effect Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 239000002994 raw material Substances 0.000 description 3
- 238000000926 separation method Methods 0.000 description 3
- 238000000151 deposition Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000007717 exclusion Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000004026 adhesive bonding Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000002902 bimodal effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000007641 inkjet printing Methods 0.000 description 1
- 229910010272 inorganic material Inorganic materials 0.000 description 1
- 239000011147 inorganic material Substances 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 150000002739 metals Chemical class 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- JHJNPOSPVGRIAN-SFHVURJKSA-N n-[3-[(1s)-1-[[6-(3,4-dimethoxyphenyl)pyrazin-2-yl]amino]ethyl]phenyl]-5-methylpyridine-3-carboxamide Chemical group C1=C(OC)C(OC)=CC=C1C1=CN=CC(N[C@@H](C)C=2C=C(NC(=O)C=3C=C(C)C=NC=3)C=CC=2)=N1 JHJNPOSPVGRIAN-SFHVURJKSA-N 0.000 description 1
- 239000011368 organic material Substances 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64C—AEROPLANES; HELICOPTERS
- B64C21/00—Influencing air flow over aircraft surfaces by affecting boundary layer flow
- B64C21/10—Influencing air flow over aircraft surfaces by affecting boundary layer flow using other surface properties, e.g. roughness
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64C—AEROPLANES; HELICOPTERS
- B64C2230/00—Boundary layer controls
- B64C2230/26—Boundary layer controls by using rib lets or hydrophobic surfaces
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F15—FLUID-PRESSURE ACTUATORS; HYDRAULICS OR PNEUMATICS IN GENERAL
- F15D—FLUID DYNAMICS, i.e. METHODS OR MEANS FOR INFLUENCING THE FLOW OF GASES OR LIQUIDS
- F15D1/00—Influencing flow of fluids
- F15D1/002—Influencing flow of fluids by influencing the boundary layer
- F15D1/0025—Influencing flow of fluids by influencing the boundary layer using passive means, i.e. without external energy supply
- F15D1/003—Influencing flow of fluids by influencing the boundary layer using passive means, i.e. without external energy supply comprising surface features, e.g. indentations or protrusions
- F15D1/0035—Influencing flow of fluids by influencing the boundary layer using passive means, i.e. without external energy supply comprising surface features, e.g. indentations or protrusions in the form of riblets
- F15D1/004—Influencing flow of fluids by influencing the boundary layer using passive means, i.e. without external energy supply comprising surface features, e.g. indentations or protrusions in the form of riblets oriented essentially parallel to the direction of flow
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F15—FLUID-PRESSURE ACTUATORS; HYDRAULICS OR PNEUMATICS IN GENERAL
- F15D—FLUID DYNAMICS, i.e. METHODS OR MEANS FOR INFLUENCING THE FLOW OF GASES OR LIQUIDS
- F15D1/00—Influencing flow of fluids
- F15D1/10—Influencing flow of fluids around bodies of solid material
- F15D1/12—Influencing flow of fluids around bodies of solid material by influencing the boundary layer
Landscapes
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Structures Of Non-Positive Displacement Pumps (AREA)
Description
本願は、2019年3月29日に日本に出願された特願2019-069153号に基づき優先権を主張し、その内容をここに援用する。
例えば、特許文献1、2には、翼に複数の突起物を設けることによって、粘性摩擦抵抗を低減することが提案されている。
特許文献1に記載の技術では、翼の表面に、翼弦方向に長く延び、翼弦方向と直交するスパン方向には実質的に平行なピークヴァレーが繰り返されるシームレスフィルムが設けられている。シームレスフィルムにおける凸部のスパン方向の断面形状は、翼の表面から離れるにつれて先細になっており、先端には鋭くとがった頂部が形成されている。
特許文献2に記載の技術では、翼弦方向に延びる中央部突起と、スパン方向において中央部突起を挟んで中央部突起と実質的に平行に配置された2つの外側突起と、を有する空力装置が、翼の表面に密集もしくは分散して配置されている。中央部突起と外側突起の間の凹部は、翼弦方向において先端側から後端側に向かって翼の表面から漸次離間する傾斜面からなる。この空力装置のスパン方向における断面形状は、先端に頂部を有する略三角形状である。
非特許文献1には、風洞実験の結果から、翼における淀み点(付着線)の近傍に、微小な孤立粗度(Discrete Roughness Elements)を設けると、横流れ不安定性を抑制できることが報告されている。各孤立粗度の形状は、略円柱状である。非特許文献1に記載の技術では、孤立粗度の配列間隔を適正に設定して横流れ不安定性に起因する乱流遷移の位置を後方にずらすことによって、摩擦低減を図る。
非特許文献2では、非特許文献1における孤立粗度の作用が数値解析によって解析(評価)されている。さらに、非特許文献2では、解析結果に対する考察に基づいて孤立粗度の改良形状が提案され、数値的な検証結果が示されている。
非特許文献2における改良型の孤立粗度の形状は、スパン方向の断面においては、正弦曲線を有し、翼弦方向における断面においてはガウス分布を有するので、Sinusoidal Roughness Elements(SRE、正弦波状粗度)と称される。
特許文献1、2に記載の技術は、上述したように、後退翼に発生する横流れ不安定性に起因する乱流遷移を抑制できないので、後退翼における摩擦抵抗をあまり低減できないという問題がある。特に、特許文献1、2における構成は、凸部が翼弦方向に長く延びているので、横流れの乱流遷移を増長する結果、摩擦抵抗が増大するおそれがある。
非特許文献2の解析によれば、擾乱の波数が低いほど擾乱の振幅が増大しやすく、飽和状態に達して乱流に遷移しやすい。これに対して、擾乱の波数が高いと飽和時の振幅が低くかつ一旦飽和しても直ちに乱流に遷移することはないという特徴がある。孤立粗度の大きさ、配置ピッチ等を調整することは、横流れにおいて特定の波数成分の擾乱を選択的に励起することによって、より低波数の擾乱の成長を抑制することに相当している。
このため、非特許文献2に記載されたように、孤立粗度を設ける場合には、孤立粗度によって励起する特定波数の「killer mode(成長率(またはN値)の比較的低い横流れモード)によって乱流遷移が起きてしまわないように、細心の注意を払わねばならない。」という問題がある。
さらに、非特許文献2の解析によれば、非特許文献1の孤立粗度は複数の高波数擾乱を同時に励起してしまい、特定の単一波数の振幅の励起能力が低いため、乱流遷移の抑制効果が低いという問題もある。
非特許文献2による改良型の正弦波状粗度は、特定波数の振幅の励起能力がより高いため、乱流遷移の抑制効果が向上している。
ただし、いずれの形状にしても、孤立粗度の高さが一定値を超えるとただちに乱流遷移が起こるので、乱流遷移制御の安定性の観点ではこの一定値を超えない範囲でより大きな抑制効果が得られることが好ましい。この意味で孤立粗度の形状には依然としてさらなる改良の余地がある。
本発明の第3の態様の隆起構造によれば、第1または第2の態様において、前記複数の隆起体は、前記淀み点を連ねる線および前記第1断面に直交する第2断面において、周期的に変化する凹部および凸部を有し、前記凸部の各頂点の高さが凸状の包絡線に沿って変化する第2凹凸形状を形成するように、配列されてもよい。
本発明の第4の態様の隆起構造によれば、第3の態様において、前記第1凹凸形状は、一定の波数を有する正弦波状であってもよい。前記包絡線は、釣り鐘型の曲線であってもよい。
本発明の第5の態様の隆起構造によれば、第1の態様において、前記複数の隆起体のそれぞれにおける最大高さは、前記流れによって形成される境界層排除高さの1.7倍以上2.3倍以下であってもよい。
上記の課題を解決するために、本発明の第6の態様の隆起構造は、物体の前縁に対して鋭角をなす方向に前記物体に沿って流れる流れの層流から乱流への遷移を抑制し、粘性による摩擦抵抗を低減する隆起構造であって、前記前縁における前記流れの淀み点よりも下流側における前記物体の表面に設けられ、前記表面における前記流れに沿った方向に延びる複数の隆起体を備える。前記複数の隆起体の各々は、先端と、後端と、中間部とを有する。前記先端、前記後端、及び前記中間部は、前記表面から隆起する凸部を形成する。前記先端は、前記表面において前記流れに沿った方向における上流側に位置しており、前記凸部と前記表面との間の境界を形成する。前記後端は、前記表面において前記流れに沿った方向における下流側に位置しており、前記凸部と前記表面との間の境界を形成する。前記表面から隆起する前記凸部の隆起方向に見て、前記中間部は、前記先端と前記後端との間に位置する。前記凸部においては、前記中間部は、前記表面からの前記隆起方向における高さが最大となる頂点に位置する。前記複数の隆起体は、前記淀み点からの距離が一定かつ前記表面に直交する第1断面において、周期的に変化する第1凹凸形状を形成する。
本発明の第8の態様の物体によれば、第7の態様において、前記表面において、前記隆起構造は、前記前縁の近くに設けられてもよい。
図1は、本発明の実施形態の隆起構造および翼の一例を示す模式的な斜視図である。
翼1は、前縁fから後縁bに向かって、翼厚が増大し、最大になった後、漸次減少する適宜の翼型を有する。具体的な翼型については特に限定されない。図1に示す例では、前縁fと後縁bとを結ぶ翼弦線に対して、上面1a(表面)の膨らみが下面1bの膨らみよりも大きな翼型を有する。
例えば、翼1は、前縁fが後縁bよりも高くなる適宜の取り付け角で、図示略の胴体に対して取り付けられている。翼の基準線は、平面視では、胴体の前後軸に対して斜め後方向に傾斜している。このため、平面視では、前縁fも、胴体の前後軸に対して斜め後方向に傾斜している。
図1において、流れFは、翼1が取り付けられた機体が直進したときに翼1が受ける空気の流れを表す。
航空機の飛行時には、翼1の先端部(物体の端部)の定位置に、流れFの淀み点St(二点鎖線参照)が形成される。淀み点Stは、翼1の先端部において前縁fに沿って線状に延びている。
隆起構造2は、翼1に向かって流れる流れF(主流)の粘性による摩擦抵抗を低減するために設けられている。隆起構造2の配置を決定するための淀み点Stの位置は、飛行状態によって変化する。隆起構造2の設計にあたっては、特に摩擦抵抗を低減したい飛行状態に形成される淀み点Stが用いられる。例えば、航空機が水平飛行する際の淀み点Stが用いられてもよい。
各単位隆起2aは、上面1a上の流れFの流れ方向(図1に模式的に示す流線S参照)に延びている。上面1a上の任意の位置の流線Sは、隆起構造2が除去された翼1の形状、流れFの速度、流れFの粘性等が与えられれば、数値解析によって求められる。
各単位隆起2aの高さは、長手方向および短手方向のいずれの断面においても、両端部から中間部に向かって漸増して中間部で最大となる凸状の滑らかな曲線に沿って変化している。例えば、各単位隆起2aの高さは、長手方向の断面においては、上流側の先端2f側および下流側の後端2b側から、先端2fおよび後端2bの間の中間部2mに向かって、漸増して中間部2mで最大となる凸状の滑らかな曲線に沿って変化している。ここで、中間部2mは、単位隆起2aの長手方向の中央であってもよいが、長手方向の中央からずれていてもよい。図示は省略するが、短手方向の断面も同様である。
このように、単位隆起2aにおいて隆起方向の先端面は丸みを帯びた滑らかな曲面である。単位隆起2aの延在方向は、短手方向の断面における頂部を連ねた尾根線で表される。尾根線は、流れFの流れ方向に延びることがより好ましい。しかし単位隆起2aの尾根線は短手方向の断面においては滑らかな凸部の停留点になっているので、尾根線は測定誤差、製作誤差等を考慮して略流れ方向に延びていればよい。
隆起構造2の具体例については後述する。
図2は、翼上の三次元境界層の流速分布の一例を示す模式図である。図3は、翼弦方向の位置xと、境界層外縁速度Ue、We、境界層排除厚さδ、およびレイノルズ数Reδとの関係を示すグラフである。図3において、横軸xは、淀み点Stからx方向の正方向に進む長さ[mm]を表す。縦軸は、境界層外縁速度Ue、We、境界層排除厚さδ、およびレイノルズ数Reδを示す。ただし、境界層排除厚さδは、1mmで無次元した数値で示されている。図4は、横流れ不安定性の中立安定曲線を示すグラフである。図5は、横流れ不安定性のN値を示すグラフである。
翼1の淀み点Stを連ねた直線に直交し、上面1aに沿う方向をx方向と称する。特に誤解の生じるおそれがない場合には、x方向を単に翼弦方向と称する場合もある。x方向の正方向は、淀み点Stから上面1aに沿って後縁bに向かう方向である。
z方向およびx方向に直交する方向をy方向と称する。y方向は、単位隆起2aの隆起方向である。
なお、隆起構造2の単位隆起2aの長さおよび配列ピッチは、後述するように、翼1の翼弦方向の幅およびスパン方向の長さに比べて充分短い。このため、隆起構造2の形成範囲において、上面1aは、解析上、平面で近似できる。
流れFは、翼1の後退角に応じて淀み点Stに対して斜め方向に入射するので、主流の境界層外縁速度は、翼弦方向およびスパン方向にそれぞれ、Ue、Weの速度成分を有する。境界層の外部では、Ue、We以外の速度成分を有しない。
これに対して、境界層の内部では、圧力勾配の方向である翼弦方向と主流の方向が平行でないため、主流(MF)の流れ方向に交差する方向に横流れ(CF)が発生する。横流れの速度分布は、上面1aの表面と、境界層外縁と、において0であり、その中間部で最大値を取る、山形の分布である。
このため、境界層の内部における合成流速分布は、主流の流れに沿う断面に対してねじれている。
図3のグラフの詳しい計算方法については、非特許文献2に記載された通りである。
翼弦方向の境界層外縁速度Ueは、淀み点St(x=0)からxが増大するにつれて、急峻に増大して、次第に飽和に向かう。スパン方向の境界層外縁速度Weは、xによらず一定である。
境界層排除厚さδは、初期に急増した後、略直線的に増大する。このため、翼弦方向における後端ほど、境界層の厚さが増大することが分かる。一方、レイノルズ数Reδはxが増大するにつれて、略直線的に増大する。このため、翼1の後端側では層流境界層が乱流遷移を起こしやすいことが分かる。
すなわち、擾乱の定在モードを考えて、擾乱を局所的に平面波(∝exp(iαx+iβz))で近似する。ここで、iは虚数単位、αは翼弦方向波数(複素数)、βはスパン方向波数(実数)である。スパン方向波数βを与えると、翼弦方向波数αが複素固有値として求まる。
図4において、横軸は、翼弦方向の位置xであり、縦軸は、局所平面波近似した擾乱におけるスパン方向波数βである。
曲線40は、擾乱において複素数で表される翼弦方向波数α(=αr+iαi)の複素数成分αiが0となる曲線である。このため、曲線40の外側(「Stable(αi>0)」と記載の領域)は、αi>0になるため、安定領域である。これに対して、曲線40の内側(「Unstable(αi<0)」と記載の領域)は、αi<0になるため、不安定領域である。
図4によれば、波数βが高い場合には、翼弦方向において前縁fから遠ざかれば安定領域になるのに対して、波数βが低くなるほど、不安定領域が後端に延び、不安定領域が拡大している。
図5に曲線51、52、53で示されるように、xが小さいほど高波数成分によって擾乱が増幅され、xが大きいほど低波数成分によって擾乱が増幅されることが分かる。さらに、増幅率の大きさは、低波数成分ほど大きくなっている。例えば、x=202(mm)ではβ=3.333(mm-1)による増幅が支配的なモードであり、N値は8を超えている。
したがって、擾乱を形成するエネルギーが高波数成分に分配されれば、N値を増大させることなく、低波数成分の擾乱の成長を抑制できると考えられる。この場合、高波数の擾乱の振幅は増大する。しかし、高波数の擾乱は、成長しても低波数の擾乱より乱流遷移を起こしにくいので、境界層全体としては乱流遷移が起こりにくくなる。
非特許文献1、2には、特に抑制したい擾乱のスパン方向波数βの1.5倍程度のスパン方向波数βの擾乱を励起しやすい間隔で孤立粗度を配列すると、乱流遷移位置をより後方にずらすことができることが報告されている。
本発明者は、孤立粗度による擾乱の励起をより効率化できれば、抑制対象の擾乱をさらに効果的に抑制できるのではないかと考えて、鋭意研究し、本発明に到った。
図6は、本発明の実施形態の隆起構造の一例を示す斜視図である。図7は、本発明の実施形態の隆起構造の一例を示す平面図である。図8における(a)、(b)は、それぞれ本発明の実施形態の隆起構造の一例のスパン方向および翼弦方向の断面形状を示すグラフである。図6~8において、x、y、zは1mmで無次元化されている。
例えば、図8の(a)に曲線101で示すように、x=xRでは、yRは、0~hまで波数βRで振動する正弦波である。x≠xRでは、(x-xR)に応じてyの最大値がh未満に変化する以外は、同様の正弦波である。
このため、yRは、第2断面において、波数{βR×(We(xR)/Ue(xR))}でそれぞれ周期的に変化する凹部および凸部(正弦波状の凹凸部)を有し、各凸部の頂点の高さが凸状の包絡線(ガウス関数)に沿って変化する第2凹凸形状を形成している。
例えば、図8の(b)に曲線102で示すように、z=0とすると、yRは、x=xRから幅方向の両端部に向かって、曲線103で表されるガウス関数に沿って振幅が減衰し、波数{βR×(We(xR)/Ue(xR))}でそれぞれ周期的に変化する凹部および凸部を形成している。図8の(b)に示す例では、ガウス関数の片幅はw=4(mm)である。
以下に、各パラメータの設定方法の一例について簡単に説明する。なお、横流れ不安定性解析において当業者に周知の事項、手法については適宜説明を省略する。
より具体的で正確な空間的成長率-αi(x,β)は、境界層流れの線形安定性解析を数値的に行えば計算することができる。
環境のノイズレベルによっても異なるが、例えば、Nが8以上の擾乱が乱流遷移を引き起こす環境であると見なしてもよい。例えば、航空機の場合、乱流遷移を引き起こすN値の閾値は、8~10であると言われている。したがって、航空機の場合、N=8は、安全側の見積もりである。
隆起構造2における波数βRは、擾乱を励起する波数である。このため、隆起構造2による励起によって乱流遷移が起こらないように、波数βRとしては、閾値に選定されたN値よりも小さいN値を与える波数が選ばれる。例えば、閾値がN=8の場合、波数βRの選定用のN値は、例えば、N=6が用いられてもよい。
この場合、式(6)の左辺にN=6を代入すれば、βRが求められる。
隆起構造2の最大高さhは、高いほど隆起構造2による擾乱を励起しやすくなるので、低波数の擾乱を抑制しやすくなる。しかし、hが高くなりすぎると、隆起構造2によって励起される擾乱の振幅が大きくなりすぎる。この結果、ただちに乱流遷移が発生するおそれがある。そのため、hは、位置xRにおける境界層排除厚さδ(xR)の3倍以下とするべきである。安全側に見積もると、hはδ(xR)の2倍程度、例えば、1.7倍以上2.3倍以下であるとより好ましい。
例えば、図8の(a)、(b)に示す例では、δ(xR)=0.087(mm)とすると、h/δ(xR)は、2.874である。
例えば、図8の(b)に示す例では、w=4(mm)であるので、wはhの16倍である。
例えば、隆起構造2における翼弦方向の凹凸形状は、卓越波数がβRであれば、正弦関数以外の周期関数で表されてもよい。ただし、周期関数は、翼弦方向における勾配の変化が滑らかであることより好ましい。
さらに、隆起構造2における翼弦方向の凹凸形状は、卓越波数がβRである周期的な凹凸形状であれば、特定の周期関数で表されない形状であってもよい。
例えば、隆起構造2におけるスパン方向の凹凸形状の包絡線は、βR以外の波数の擾乱を励起しにくいように、凹凸形状の高さがなだらかに変化する凸形状であれば、ガウス関数には限定されない。例えば、凸形状は、流線型、釣り鐘型など、流れを擾乱しにくい形状であれば、具体的な関数形は特に限定されず、関数で表されない形状であってもよい。
例えば、隆起構造2は、上面1aの成形時に上面1aに形成されてもよい。
例えば、隆起構造2は、必要な凹凸形状が表面に形成されたシート部材として製造されてもよい。シート部材は、上面1a上に接着などによって、翼1に取り付けられてもよい。シート部材における隆起構造2の形成方法としては、モールド成形、プレス成形、除去加工などが挙げられる。
例えば、隆起構造2は、上面1a上に塗布された原料を必要な凹凸形状に硬化させることによって、上面1a上に形成されてもよい。原料の硬化方法としては、例えば、紫外線照射、加熱などが用いられてもよい。凹凸形状の形成方法としては、成形型が用いられてもよいし、層状部を形成した後に不要部が除去加工されてもよい。例えば、除去加工としては、機械的除去加工、化学的除去加工、レーザ加工などが挙げられる。
例えば、隆起構造2は、原料を上面1aまたは上面1aに貼り付けるシート部材上に付加または堆積することによって形成されてもよい。付加・堆積方法としては、特に限定されないが、例えば、3Dプリンタ、インクジェット印刷などが挙げられる。
隆起構造2を構成する各単位隆起2aは、境界層外縁における流線(図1のS)の方向に沿って延びている。各単位隆起2aの長手方向および短手方向の端部は、上面1aからなだらかに高さが増大する滑らかな曲面である。さらに、各単位隆起2aの中間部2mは、長手方向および短手方向において、頂部まで漸次高さが増大する滑らかな凸曲面である。
このような形状を有するため、単位隆起2aは、図2の主流(MF)方向の流れを阻害しにくい流線型であり、主流を円滑に通過させることができる隆起である。
隆起構造2を全体として見ても、単位隆起2aがスパン方向に周期的に整列しているため、同じく主流方向に通過する流れを阻害しにくい。
隆起構造2は、包絡線をスパン方向に平行移動したドーム状の範囲で全体としてスパン方向に長く延びており、スパン方向の断面形状は、波数βRの周期的な凹凸形状からなる。このため、主流(MF)方向に垂直な方向成分を有する横流れ(CF)が波数βRの周期的な凹凸形状を通過することによって、波数βRの擾乱が励起される。
このように、隆起構造2は、淀み点Stからの距離xRの位置において、境界層内部の専ら横流れ成分に波数βRの擾乱を選択的に励起する粗度構造を形成している。
図9は、本発明の実施形態の隆起構造の作用を示すグラフである。図9において、横軸は翼弦方向の位置x、縦軸は擾乱の振幅に対応するmaxy|vn|(ただし、vにおけるハット記号(サーカムフレックスアクセント)は省略。以下も同じ。)を示す。なお、図9は縦軸が対数目盛とされた片対数グラフである。
スパン方向波数βは、β=1.666n(mm-1)(ただし、n=1,2,…)のように整数nでラベル付けして表している。また、グラフの見やすさを考慮して、図9には、n=1,2,3のみの結果を示している。
計算に用いた隆起構造2の形状は、上記式(1)において、βR=5.000(mm-1)(n=3)、We(xR)/Ue(xR)=0.57、h=0.25(mm)、w=4(mm)、xR=30(mm)とされた形状である。
ここで、βR=5.000(mm-1)としているのは、N値が大きいβ=3.333(mm-1)の擾乱を抑制するためである。本例では、スパン方向波数βRは、擾乱の抑制対象のモード(ターゲットモードと称する)の波数β=3.333(mm-1)の1.5倍である。
ターゲットモードを抑制するために、隆起構造2によって人為的に擾乱を励起する波数が支配的な擾乱のモードをキラーモードと称する。
曲線110に示すように、隆起構造2を設けない場合、擾乱の振幅は、x=20(mm)から、横流れ不安定性によって指数関数的に増大し、x=180(mm)の付近で飽和している。さらにxが増加すると、x=200(mm)の付近から二次不安定性が発生し始め、x=250(mm)の付近ではほぼ乱流に遷移していることが分かる。ただし、ここでは二次不安定性の引き金として、x=180(mm)の位置にて微小な振動擾乱を数値的に与えている。
図9において、曲線111、112、113は、それぞれn=1、2、3のスパン方向波数βに対応する擾乱の振幅のフーリエ成分の変化を表す。
曲線113に示すように、n=3の場合、x=30(mm)を中心として配置された隆起構造2の形成領域で、振幅が急峻に増大して最大値を取り、この後、xが増加するにつれて、振幅が減少している。特に、x=250(mm)の付近では、振幅は略0である。さらにxが増加すると振幅は増大するが、x=400(mm)でも乱流遷移は発生していない。
曲線113の最大値は、曲線110における乱流遷移時の振幅を上回っているにもかかわらず乱流遷移が発生していない。この理由は、曲線110における乱流遷移はβ=5.000(mm-1)よりも低い波数の擾乱の成長によって形成されているためであると考えられる。
曲線111に示すターゲットモードよりも低波数(n=1)の場合、xが増加するにつれて指数関数的に増大している。ターゲットモードに比べると振幅の抑制効果は少ないが、x=400(mm)でも乱流遷移は発生していない。
図9における図示は省略するが、nが4以上のモードの場合、3の倍数のモードではn=3と同様に擾乱が励起され、その他のモードでは振幅が抑制される。いずれのモードでもより低波数の振幅よりは振幅が低減される。このため、図9に示すように、キラーモード、ターゲットモード、ターゲットモードよりも低波数のモードの振幅成長を検討するだけでも充分に隆起構造2の効果を予測できる。
図10は、比較例の隆起構造を示す斜視図である。図11は、比較例の隆起構造を示す平面図である。図12における(a)、(b)は、それぞれ比較例の隆起構造のスパン方向および翼弦方向の断面形状を示すグラフである。
図10~12において、x、y、zは1mmで無次元化されている。
隆起構造20は、非特許文献1に開示の孤立粗度と略同様の構成を有する。ただし、具体的な寸法は、隆起構造2と対比しやすいように適宜変更している。
隆起構造20によるターゲットモードは、実施形態と同様、β=3.333(mm-1)である。
図12の(a)に曲線201で示すように、単位隆起20aは、平均直径dが1mm、高さHが0.12mmの略円柱状である。ただし、単位隆起20aの側面は突出方向に先細の傾斜を有しており、先端部の外縁には丸みが付けられている。
図12の(b)に示すように、単位隆起20aは、x=xRにおいて、スパン方向に等ピッチの1列をなして配列されている。単位隆起20aのスパン方向の配列ピッチp(図12の(a)参照)は、β=5.000(mm-1)の擾乱を励起するため、2π/5.000mmである。
図13は、比較例の隆起構造の作用を示すグラフである。図13の縦軸および横軸は、図9の縦軸および横軸と同様である。
曲線211、212、213は、それぞれn=1、2、3のスパン方向波数βに対応する擾乱の振幅のフーリエ成分の変化を表す。
曲線213に示すように、n=3のキラーモードの場合、x=30(mm)を中心として配置された隆起構造20の形成領域で、振幅が急峻に増大して最大値を取り、この後、急峻に振幅が低下した後、振幅は増加に転じた。n=3のキラーモードの場合、振幅は、x=100(mm)の近辺で極大、x=200(mm)の近辺で極小、x=260(mm)の近辺で極大となる2峰性の変化を示した。この後、振幅は、x=300(mm)の近辺で不安定になった後、急峻に増大した。x=330(mm)の近辺では、乱流に遷移したことが分かる。
曲線211に示す、ターゲットモードよりも低波数(n=1)の場合、ターゲットモードと同様に、xが増加するにつれて指数関数的に増大した。ただし、曲線211に示すように、振幅は、x=160(mm)を超えても増大し、x=300(mm)の付近で飽和した後、乱流に遷移したことが分かる。
比較例では、n=1、2、3のいずれもモードでも、x=330(mm)の付近では乱流遷移していた。このため、隆起構造20の効果は、乱流遷移位置をx=250(mm)からx=330(mm)まで、80mm後退させることができる程度である。
これに対して、本実施形態の隆起構造2を有する計算結果によれば、少なくとも計算範囲のx=400(mm)まで乱流遷移が発生しておらず、発生する兆候もなかった。したがって、隆起構造2は、乱流遷移位置は150mm以上後退させることができている。
図14は、隆起構造の高さとキラー振幅との関係を示すグラフである。
図14における横軸は隆起構造の最大高さh、縦軸はx=40(mm)におけるキラーモードにおける擾乱の振幅である。なお、本解析において、xR=30(mm)における境界層排除厚さδは0.087mmであった。
折れ線120、220は、それぞれ隆起構造2、20の結果を示す。
折れ線120に示すように、隆起構造2では、最大高さhが増大するにつれてキラー振幅は、増大し、h=0.16(mm)の近辺で略一定になった。しかし、h=0.25(mm)の近辺からキラー振幅はさらに増大した。ただし、hが0.30(mm)以上では、xR=30(mm)より後方からただちに乱流領域Tに入っていた。
折れ線220に示すように、隆起構造20では、h=0.07(mm)までは、最大高さhの増加とともにキラー振幅が単調に増加した。しかし、変化率は、隆起構造2に比べると格段に小さかった。さらに、キラー振幅は、h=0.07(mm)で略飽和し、h=0.10(mm)を超えると急峻に減少した。
隆起構造20を有する場合、hが0.10(mm)を超えると、キラーモードはほとんど励起されなかった。さらに、h=0.15(mm)以上になると、やはりxR=30(mm)より後方からただちに乱流領域Tに入っていた。
hが0.25mmを超えると、乱流遷移が抑制できない理由としては、hが境界層排除厚さδの3倍を超えたため、隆起構造2の凸部が境界層外部の高速流領域にまで突出し、高速流領域における凸部において流れの剥離による乱流遷移をただちに引き起こしたことなどが考えられる。流れの剥離が起きると乱流遷移を逆に促進してしまうので、安全側に見積もるならば、hは、δの1.7倍以上2.3倍以下であることがより好ましい。
比較例の隆起構造20の場合、隆起構造2に比べるとキラー振幅の励起能力が低く、乱流遷移の抑制効果も低かった。例えば、hが0.10mmを超えると励起能力を有していなかった。略円柱形状が流線型とはかけ離れているために、h/δが2に満たない程度でも剥離による乱流遷移を引き起こすようになると考えられる。
隆起構造2および翼1によれば、横流れ不安定性に起因する摩擦抵抗を低減することができる。
隆起構造を設けることができる物体の例としては、航空機以外に用いられる翼、流れ中の後退円柱、流れ方向の上流側から下流側に向かって拡径する円錐状の回転体などが挙げられる。
さらに、隆起構造が設けられる物体が配置される流体は、粘性流体であれば、粘性流体の種類は限定されない。粘性流体は圧縮性流体であってもよいし、非圧縮性流体でもよい。
また、本発明は前述した説明によって限定されることはなく、添付の特許請求の範囲によってのみ限定される。
1a 上面(表面)
1b 下面
2 隆起構造
2a 単位隆起(隆起体)
f 前縁
F 流れ
St 淀み点
Claims (8)
- 物体の前縁に対して鋭角をなす方向に前記物体に沿って流れる流れの粘性による摩擦抵抗を低減する隆起構造であって、
前記前縁における前記流れの淀み点よりも下流側における前記物体の表面に設けられ、前記表面における前記流れに沿った方向に延びる複数の隆起体を備え、
前記複数の隆起体の各々は、先端と、後端と、中間部とを有し、
前記先端、前記後端、及び前記中間部は、前記表面から隆起する凸部を形成し、
前記先端は、前記表面において前記流れに沿った方向における上流側に位置しており、前記凸部と前記表面との間の境界を形成し、
前記後端は、前記表面において前記流れに沿った方向における下流側に位置しており、前記凸部と前記表面との間の境界を形成し、
前記表面から隆起する前記凸部の隆起方向に見て、前記中間部は、前記先端と前記後端との間に位置し、
前記凸部においては、前記中間部は、前記表面からの前記隆起方向における高さが最大となる頂点に位置し、
前記複数の隆起体は、
前記淀み点からの距離が一定かつ前記表面に直交する第1断面において、周期的に変化する第1凹凸形状を形成する、
隆起構造。 - 前記凸部においては、前記先端から前記中間部に向けて、前記表面からの高さが漸増するように滑らかな曲線が形成されており、
前記凸部においては、前記後端から前記中間部に向けて、前記表面からの高さが漸増するように滑らかな曲線が形成されている、
請求項1に記載の隆起構造。 - 前記複数の隆起体は、前記淀み点を連ねる線および前記第1断面に直交する第2断面において、周期的に変化する凹部および凸部を有し、前記凸部の各頂点の高さが凸状の包絡線に沿って変化する第2凹凸形状を形成するように、配列されている、
請求項1または2に記載の隆起構造。 - 前記第1凹凸形状は、
一定の波数を有する正弦波状であり、
前記包絡線は、
釣り鐘型の曲線である、
請求項3に記載の隆起構造。 - 前記複数の隆起体のそれぞれにおける最大高さは、
前記流れによって形成される境界層排除高さの1.7倍以上2.3倍以下である、
請求項1に記載の隆起構造。 - 物体の前縁に対して鋭角をなす方向に前記物体に沿って流れる流れの層流から乱流への遷移を抑制し、粘性による摩擦抵抗を低減する隆起構造であって、
前記前縁における前記流れの淀み点よりも下流側における前記物体の表面に設けられ、前記表面における前記流れに沿った方向に延びる複数の隆起体を備え、
前記複数の隆起体の各々は、先端と、後端と、中間部とを有し、
前記先端、前記後端、及び前記中間部は、前記表面から隆起する凸部を形成し、
前記先端は、前記表面において前記流れに沿った方向における上流側に位置しており、前記凸部と前記表面との間の境界を形成し、
前記後端は、前記表面において前記流れに沿った方向における下流側に位置しており、前記凸部と前記表面との間の境界を形成し、
前記表面から隆起する前記凸部の隆起方向に見て、前記中間部は、前記先端と前記後端との間に位置し、
前記凸部においては、前記中間部は、前記表面からの前記隆起方向における高さが最大となる頂点に位置し、
前記複数の隆起体は、
前記淀み点からの距離が一定かつ前記表面に直交する第1断面において、周期的に変化する第1凹凸形状を形成する、
隆起構造。 - 請求項1に記載の隆起構造を表面に備え、流体に配置される物体。
- 前記表面において、前記隆起構造は、前記前縁の近くに設けられている、
請求項7に記載の物体。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019069153 | 2019-03-29 | ||
JP2019069153 | 2019-03-29 | ||
PCT/JP2020/011911 WO2020203284A1 (ja) | 2019-03-29 | 2020-03-18 | 隆起構造および翼 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPWO2020203284A1 JPWO2020203284A1 (ja) | 2020-10-08 |
JP7465483B2 true JP7465483B2 (ja) | 2024-04-11 |
Family
ID=72668750
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2021511412A Active JP7465483B2 (ja) | 2019-03-29 | 2020-03-18 | 隆起構造および物体 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20220161923A1 (ja) |
JP (1) | JP7465483B2 (ja) |
WO (1) | WO2020203284A1 (ja) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP7369386B1 (ja) | 2022-05-23 | 2023-10-26 | 三菱重工業株式会社 | 隆起構造、翼、隆起構造の設計方法及びその設計プログラム |
CN116552777B (zh) * | 2023-07-05 | 2023-09-12 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种涡流调控器以及一种飞行器 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009000703A1 (en) | 2007-06-25 | 2008-12-31 | Politecnico Di Milano | Method for reducing the viscous friction between a fluid and an object |
WO2014026246A1 (en) | 2012-08-16 | 2014-02-20 | Adelaide Research & Innovation Pty Ltd | Improved wing configuration |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3129908A (en) * | 1961-08-25 | 1964-04-21 | Richard G Harper | Device for selectively altering lift characteristics of an airfoil |
US6345791B1 (en) * | 2000-04-13 | 2002-02-12 | Lockheed Martin Corporation | Streamwise variable height riblets for reducing skin friction drag of surfaces |
-
2020
- 2020-03-18 JP JP2021511412A patent/JP7465483B2/ja active Active
- 2020-03-18 US US17/599,092 patent/US20220161923A1/en active Pending
- 2020-03-18 WO PCT/JP2020/011911 patent/WO2020203284A1/ja active Application Filing
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009000703A1 (en) | 2007-06-25 | 2008-12-31 | Politecnico Di Milano | Method for reducing the viscous friction between a fluid and an object |
WO2014026246A1 (en) | 2012-08-16 | 2014-02-20 | Adelaide Research & Innovation Pty Ltd | Improved wing configuration |
Also Published As
Publication number | Publication date |
---|---|
JPWO2020203284A1 (ja) | 2020-10-08 |
US20220161923A1 (en) | 2022-05-26 |
WO2020203284A1 (ja) | 2020-10-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20070194178A1 (en) | Passive micro-roughness array for drag modification | |
EP3363733B1 (en) | Passive flow control mechanism for reducing and/or suppressing tollmien-schlichting waves, delaying transition to turbulence and reducing drag | |
US5167387A (en) | Porous airfoil and process | |
US8460779B2 (en) | Microstructures for reducing noise of a fluid dynamic structure | |
Glezer | Some aspects of aerodynamic flow control using synthetic-jet actuation | |
US6345791B1 (en) | Streamwise variable height riblets for reducing skin friction drag of surfaces | |
JP7465483B2 (ja) | 隆起構造および物体 | |
US20150017385A1 (en) | Passive drag modification system | |
Hrynuk et al. | The effects of leading-edge tubercles on dynamic stall | |
JPH01122784A (ja) | 非流線形のベース部を有する物体及び陸上車輌 | |
Greenblatt et al. | Influence of finite span and sweep on active flow control efficacy | |
Klän et al. | Experimental analysis of the flow field over a novel owl based airfoil | |
JP2005529298A (ja) | 境界層の流体の流れの制御 | |
US11577823B2 (en) | Noise reduction apparatus, aircraft, and noise reduction method | |
US11614106B2 (en) | Partially submerged periodic riblets | |
Arivoli et al. | Rudimentary emulation of covert feathers on low-AR wings for poststall lift enhancement | |
Ashraf et al. | Thrust generation and wake structure for flow across a pitching airfoil at low Reynolds number | |
US20210231141A1 (en) | Staggered periodic riblets | |
Vasile et al. | Interaction of a finite-span synthetic jet near the tip of a sweptback wing | |
Genç et al. | Unsteady flow over flexible wings at different low Reynolds numbers | |
Gong et al. | On the aerodynamic loads and flow statistics of airfoil with deformable vortex generators | |
WO2019116697A1 (ja) | 遷音速翼型、翼及び航空機 | |
Tregidgo et al. | Fluid-structure interactions for a low aspect-ratio membrane wing at low Reynolds numbers | |
Kaakkunen et al. | Fabrication of surfaces with reduced friction using nanosecond laser | |
US10988233B2 (en) | Wing and aircraft |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
RD01 | Notification of change of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7426 Effective date: 20220207 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20220207 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20230208 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20231017 |
|
A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20231218 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20240216 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20240305 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20240321 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 7465483 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |