WO2020203284A1 - 隆起構造および翼 - Google Patents

隆起構造および翼 Download PDF

Info

Publication number
WO2020203284A1
WO2020203284A1 PCT/JP2020/011911 JP2020011911W WO2020203284A1 WO 2020203284 A1 WO2020203284 A1 WO 2020203284A1 JP 2020011911 W JP2020011911 W JP 2020011911W WO 2020203284 A1 WO2020203284 A1 WO 2020203284A1
Authority
WO
WIPO (PCT)
Prior art keywords
raised structure
cross
raised
section
flow
Prior art date
Application number
PCT/JP2020/011911
Other languages
English (en)
French (fr)
Inventor
廣田　真
裕司 服部
優紀 井手
光 ▲高▼見
吉本　稔
Original Assignee
国立大学法人東北大学
国立研究開発法人宇宙航空研究開発機構
三菱重工業株式会社
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 国立大学法人東北大学, 国立研究開発法人宇宙航空研究開発機構, 三菱重工業株式会社 filed Critical 国立大学法人東北大学
Priority to JP2021511412A priority Critical patent/JP7465483B2/ja
Priority to US17/599,092 priority patent/US20220161923A1/en
Publication of WO2020203284A1 publication Critical patent/WO2020203284A1/ja

Links

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

Definitions

  • the present invention relates to raised structures and blades that reduce viscous frictional resistance.
  • the present application claims priority based on Japanese Patent Application No. 2019-069153 filed in Japan on March 29, 2019, the contents of which are incorporated herein by reference.
  • Objects placed in the flow of viscous fluid are subject to viscous frictional resistance.
  • an aircraft is subject to frictional resistance due to the viscosity of air during flight. It is said that about half of the total resistance that an aircraft receives during flight is frictional resistance. If the frictional resistance is reduced, the cruising range of the aircraft can be extended and the fuel consumption can be reduced. Therefore, various means for reducing the viscous frictional resistance have been proposed.
  • Patent Documents 1 and 2 propose to reduce viscous frictional resistance by providing a plurality of protrusions on a wing.
  • a seamless film is provided on the surface of the wing, which extends long in the chord direction and repeats a peak valley substantially parallel to the span direction orthogonal to the chord direction.
  • the cross-sectional shape of the convex portion in the span direction in the seamless film tapers toward the surface of the wing, and a sharply pointed top is formed at the tip.
  • an aerodynamic device having a central protrusion extending in the chord direction and two outer protrusions arranged substantially parallel to the central protrusion with the central protrusion sandwiched in the span direction. However, they are densely or dispersedly arranged on the surface of the wing.
  • the recess between the central protrusion and the outer protrusion consists of an inclined surface that gradually separates from the surface of the wing from the front end side to the rear end side in the chord direction.
  • the cross-sectional shape of this aerodynamic device in the span direction is a substantially triangular shape having a top at the tip.
  • Patent Documents 1 and 2 are considered to be techniques for reducing friction by forming a stable turbulent boundary layer on the surface of the blade by means of protrusions.
  • the techniques described in Patent Documents 1 and 2 cannot reduce the frictional resistance due to the turbulent flow transition caused by the cross flow instability in the swept wing widely used in aircraft.
  • Non-Patent Document 1 reports from the results of wind tunnel experiments that cross-flow instability can be suppressed by providing a minute isolated roughness (Discrete Roughness Elements) in the vicinity of the stagnation point (adhesion line) on the wing. ..
  • the shape of each isolated roughness is substantially cylindrical.
  • Non-Patent Document 2 proposes an improved shape with isolated roughness based on consideration for the analysis result, and shows the numerical verification result. Since the improved isolated roughness shape in Non-Patent Document 2 has a sinusoidal curve in the cross section in the span direction and a Gaussian distribution in the cross section in the chord direction, Sinusoidal Roughness Elements (SRE, sinusoidal roughness) Degree).
  • SRE sinusoidal roughness
  • Patent Documents 1 and 2 have a problem that the frictional resistance in the swept wing cannot be reduced so much because the turbulent transition caused by the cross flow instability generated in the swept wing cannot be suppressed.
  • the frictional resistance may increase as a result of increasing the turbulent transition of the cross flow.
  • the lower the wave number of the disturbance the more easily the amplitude of the disturbance increases, and the more the saturation state is reached, the easier it is to transition to the turbulent flow.
  • Non-Patent Document 2 the isolated roughness of Non-Patent Document 1 excites a plurality of high wavenumber disturbances at the same time, and the excitation ability of the amplitude of a specific single wavenumber is low. There is also a problem that the suppressive effect of is low.
  • the improved sinusoidal roughness according to Non-Patent Document 2 has a higher excitation ability of the amplitude of a specific wave number, and therefore has an improved effect of suppressing turbulent flow transition.
  • turbulent transition occurs immediately when the height of isolated roughness exceeds a certain value, so from the viewpoint of stability of turbulent transition control, greater suppression is performed within the range not exceeding this constant value. It is preferable that the effect is obtained. In this sense, there is still room for further improvement in the isolated roughness shape.
  • the present invention has been made in view of the above problems, and an object of the present invention is to provide a raised structure and a blade capable of reducing frictional resistance due to cross-flow instability.
  • the raised structure according to the first aspect of the present invention is a raised structure that reduces frictional resistance due to mainstream viscosity flowing toward the object in a direction forming a sharp angle with respect to the end of the object.
  • a plurality of ridges provided on the surface of the object downstream of the stagnation point of the mainstream at the end and extending in the flow direction of the mainstream on the surface, and each of the plurality of ridges.
  • the height of the object gradually increases from both ends toward the intermediate portion and changes along a convex smooth curve that is maximum at the intermediate portion.
  • the plurality of raised bodies form a first uneven shape that changes periodically in a first cross section having a constant distance from the stagnation point and orthogonal to the surface, and a line (adhesion line) connecting the stagnation points and a line (adhesion line) connecting the stagnation points.
  • a second uneven shape having a concave portion and a convex portion that change periodically and the height of each apex of the convex portion changes along a convex wrapping line They are arranged to form.
  • the first uneven shape is a sinusoidal shape having a constant wave number, and the envelope is a bell-shaped curve. Good.
  • the maximum height of each of the plurality of raised bodies is 1. the boundary layer exclusion height formed by the mainstream. It may be 7 times or more and 2.3 times or less.
  • the wing of the fourth aspect of the present invention is provided with a raised structure of any one of the first to third aspects on the surface.
  • the frictional resistance caused by the cross flow instability can be reduced.
  • the position x of the chordwise is a graph showing the boundary layer edge speed U e, W e, boundary layer displacement thickness of [delta], and the relationship between the Reynolds number Re [delta].
  • FIG. 1 is a schematic perspective view showing an example of a raised structure and a wing according to an embodiment of the present invention.
  • the wing 1 (object) of the present embodiment shown in FIG. 1 is used, for example, as a swept wing of an aircraft.
  • FIG. 1 illustrates only a part of the blade 1 in the span direction.
  • the blade 1 has an appropriate airfoil in which the blade thickness increases from the front edge f to the trailing edge b, reaches a maximum, and then gradually decreases.
  • the specific airfoil is not particularly limited.
  • the bulge of the upper surface 1a (surface) is larger than the bulge of the lower surface 1b with respect to the chord line connecting the front edge f and the trailing edge b.
  • the wing 1 is attached to a fuselage (not shown) at an appropriate attachment angle at which the front edge f is higher than the trailing edge b.
  • the reference line of the wing is inclined diagonally backward with respect to the front-rear axis of the fuselage. Therefore, in a plan view, the front edge f is also inclined diagonally backward with respect to the front-rear axis of the fuselage.
  • the flow F represents the flow of air received by the wing 1 when the airframe to which the wing 1 is attached goes straight.
  • a stagnation point St (see the alternate long and short dash line) of the flow F is formed at a fixed position at the tip of the wing 1 (the end of the object). The stagnation point St extends linearly along the front edge f at the tip of the wing 1.
  • the raised structure 2 of the present embodiment is formed on the upper surface 1a on the downstream side of the stagnation point St at the tip of the wing 1.
  • the raised structure 2 is provided to reduce frictional resistance due to the viscosity of the flow F (mainstream) flowing toward the blade 1.
  • the position of the stagnation point St for determining the arrangement of the raised structure 2 changes depending on the flight condition.
  • a stagnation point St formed in a flight state in which frictional resistance is particularly desired to be reduced is used.
  • the stagnation point St when the aircraft flies horizontally may be used.
  • the raised structure 2 includes a plurality of unit raised 2a (raised bodies). Each unit ridge 2a extends in the flow direction of the flow F on the upper surface 1a (see streamline S schematically shown in FIG. 1).
  • the streamline S at an arbitrary position on the upper surface 1a can be obtained by numerical analysis given the shape of the blade 1 from which the raised structure 2 has been removed, the velocity of the flow F, the viscosity of the flow F, and the like.
  • the height of each unit ridge 2a gradually increases from both ends toward the middle part and changes along a convex smooth curve that becomes maximum in the middle part in both the longitudinal direction and the lateral direction. ing.
  • each unit ridge 2a is, in the longitudinal cross section, from the upstream tip 2f side and the downstream rear end 2b side toward the intermediate portion 2m between the tip 2f and the rear end 2b. It gradually increases and changes along a convex smooth curve that becomes maximum at the middle part 2 m.
  • the intermediate portion 2m may be the center of the unit ridge 2a in the longitudinal direction, but may be deviated from the center in the longitudinal direction.
  • the cross section in the lateral direction is the same.
  • the tip surface in the ridge direction is a rounded and smooth curved surface.
  • the extending direction of the unit ridge 2a is represented by a ridge line connecting the tops in the cross section in the lateral direction. It is more preferable that the ridge line extends in the flow direction of the flow F. However, since the ridge line of the unit ridge 2a is a stop point of a smooth convex portion in the cross section in the lateral direction, the ridge line may extend in the substantially flow direction in consideration of measurement error, manufacturing error, and the like. ..
  • the unit ridges 2a of the ridge structure 2 are arranged at equal intervals on a line having a constant distance from the stagnation point St.
  • the raised structure 2 forms a first uneven shape that changes periodically in the cross section (first cross section) along the arrangement direction.
  • the raised structure 2 forms a second uneven shape in which the height of the apex of each convex portion changes along the convex envelope in the line connecting the stagnation points St and the second cross section orthogonal to the first cross section. ing.
  • a specific example of the raised structure 2 will be described later.
  • FIG. 2 is a schematic diagram showing an example of the flow velocity distribution of the three-dimensional boundary layer on the blade.
  • Figure 3 is a graph showing the position x of the chordwise, boundary layer edge speed U e, W e, boundary layer displacement thickness of [delta], and the relationship between the Reynolds number Re [delta].
  • the horizontal axis x represents the length [mm] that advances in the positive direction in the x direction from the stagnation point St.
  • the vertical axis represents the boundary layer edge speed U e, W e, boundary layer displacement thickness of [delta], and the Reynolds number Re [delta].
  • the boundary layer exclusion thickness ⁇ is indicated by a dimensionless numerical value of 1 mm.
  • FIG. 4 is a graph showing a neutral stability curve of cross flow instability.
  • FIG. 5 is a graph showing the N value of cross flow instability.
  • the direction along the straight line connecting the stagnation points St is referred to as the z direction.
  • the z direction is also the span direction of the blade 1.
  • the positive direction in the z direction is a direction from the fuselage (not shown) toward the tip of the wing 1 in the longitudinal direction.
  • the direction orthogonal to the straight line connecting the stagnation points St of the wing 1 and along the upper surface 1a is referred to as the x direction.
  • the x direction may be simply referred to as the chord direction.
  • the positive direction in the x direction is the direction from the stagnation point St toward the trailing edge b along the upper surface 1a.
  • the direction orthogonal to the z direction and the x direction is referred to as the y direction.
  • the y direction is the uplift direction of the unit uplift 2a.
  • the length and arrangement pitch of the unit ridge 2a of the ridge structure 2 are sufficiently shorter than the width in the chord direction and the length in the span direction of the wing 1, as will be described later. Therefore, in the formation range of the raised structure 2, the upper surface 1a can be approximated by a plane in analysis.
  • FIG. 2 shows an example of the flow velocity distribution of the three-dimensional boundary layer formed on the upper surface 1a at the tip of the blade 1.
  • the U-axis represents the flow velocity in the chord direction
  • the y / ⁇ axis is the position in the y direction standardized by the boundary layer exclusion thickness ⁇ .
  • the main flow of the boundary layer edge rate respectively chordwise and spanwise velocity component of U e, W e Has. Outside the boundary layer, no U e, the velocity component other than the W e.
  • the chord direction which is the direction of the pressure gradient, and the mainstream direction are not parallel, so that a cross flow (CF) occurs in a direction intersecting the mainstream (MF) flow direction.
  • the velocity distribution of the cross flow is 0 at the surface of the upper surface 1a and the outer edge of the boundary layer, and is a chevron distribution having a maximum value in the middle portion thereof. Therefore, the synthetic flow velocity distribution inside the boundary layer is twisted with respect to the cross section along the mainstream flow.
  • FIG. 3 shows an analysis example of the change in the boundary layer edge speed U e, W e, boundary layer displacement thickness of [delta], and the Reynolds number Re [delta] in the x direction when having cross flow.
  • ⁇ 0 is the kinematic viscosity coefficient of air (fluid) at the stagnation point St.
  • the detailed calculation method of the graph of FIG. 3 is as described in Non-Patent Document 2.
  • Boundary layer edge speed W e in the span direction is constant regardless of the x.
  • the boundary layer exclusion thickness ⁇ increases substantially linearly after an initial rapid increase. Therefore, it can be seen that the thickness of the boundary layer increases toward the rear end in the chord direction.
  • the Reynolds number Re ⁇ increases substantially linearly as x increases. Therefore, it can be seen that the laminar boundary layer is likely to cause a turbulent flow transition on the rear end side of the blade 1.
  • the turbulent flow transition is considered to be excited by deposits on the upper surface 1a, minute irregularities on the blade surface itself, and the like.
  • a neutral stability curve of cross-flow instability as shown by the curve 40 shown in FIG. 4 can be obtained. That is, considering the standing mode of the disturbance, the disturbance is locally approximated by a plane wave ( ⁇ exp (i ⁇ x + i ⁇ z)).
  • i is an imaginary unit
  • ⁇ is a wave number in the chord direction (complex number)
  • ⁇ is a wave number in the span direction (real number).
  • the chord-direction wavenumber ⁇ Given the span-direction wavenumber ⁇ , the chord-direction wavenumber ⁇ can be obtained as a complex eigenvalue.
  • the horizontal axis is the position x in the chord direction
  • the vertical axis is the span-direction wavenumber ⁇ in the disturbance approximated by the local plane wave.
  • the stable region when the wave number ⁇ is high, the stable region is reached when the wave number ⁇ is far from the front edge f in the chord direction, whereas the unstable region extends to the rear end as the wave number ⁇ is low, which is not good.
  • the stable area is expanding.
  • FIG. 5 shows a graph in which the N value (N factor) with respect to the wave number ⁇ (mm -1 ) in the span direction is plotted.
  • N values is determined by the amplification factor of disturbance (gain) and e N. The larger the N value, the more likely the turbulent transition will occur.
  • curves 51, 52, and 53 in FIG. 5 it can be seen that the smaller x is, the more the disturbance is amplified by the high wavenumber component, and the larger the x is, the more the disturbance is amplified by the low wavenumber component.
  • the magnitude of the amplification factor increases as the wave number component increases.
  • Non-Patent Documents 1 and 2 if the isolated roughness is arranged at intervals that easily excite the disturbance of the span direction wave number ⁇ , which is about 1.5 times the span direction wave number ⁇ of the disturbance to be suppressed, the turbulence transition position is further arranged. It has been reported that it can be moved backwards. The present inventor has made diligent research, thinking that if the excitation of the disturbance due to the isolated roughness can be made more efficient, the disturbance to be suppressed can be suppressed more effectively, and the present invention has been reached.
  • FIG. 6 is a perspective view showing an example of a raised structure according to an embodiment of the present invention.
  • FIG. 7 is a plan view showing an example of a raised structure according to an embodiment of the present invention.
  • FIGS. 8A and 8B are graphs showing cross-sectional shapes in the span direction and the chord direction of an example of the raised structure according to the embodiment of the present invention, respectively.
  • x, y, and z are dimensionless at 1 mm.
  • the shape y y R (x, z) of the raised structure 2 shown in FIGS. 6 and 7 is represented by the following equations (1) to (5).
  • the domain of z in the equation (1) is appropriately set within the range of the wingspan of the blade 1.
  • it is more preferable that the edge of the domain of the raised structure 2 in the z direction is further connected to an end where y gradually decreases toward 0.
  • h is the maximum height of the raised structure 2
  • ⁇ R is the wave number in the span direction (wavelength is 2 ⁇ / ⁇ R )
  • x R is the position of the raised structure 2 in the chord direction
  • w is the raised structure in the chord direction. It is one width of 2.
  • y R is only z. It is a function and is a sine wave (first uneven shape) having a wave number ⁇ R.
  • ⁇ R wave number
  • y R is a sine wave that oscillates at a wave number ⁇ R from 0 to h.
  • the sine wave is the same except that the maximum value of y changes to less than h according to (xx R ).
  • y R is a function of x only in the cross section in the span direction obtained by fixing the position z in the span direction, that is, in the line connecting the stagnation points St and the second cross section orthogonal to the first cross section.
  • y R in the second section, wavenumber ⁇ R ⁇ (W e ( x R) / U e (x R)) ⁇ recesses and protrusions periodically change respectively (sinusoidal uneven portion ), And a second uneven shape in which the height of the apex of each convex portion changes along a convex envelope (Gaussian function) is formed.
  • Gaussian function a second uneven shape in which the height of the apex of each convex portion changes along a convex envelope
  • the parameters ⁇ R , x R , h, and w in the equation (1) are set for the purpose of increasing the length of the laminar boundary layer on the upper surface 1a by shifting the turbulent transition position in the chord direction to the downstream side. Will be done. An example of how to set each parameter will be briefly described below. In the cross-flow instability analysis, items and methods well known to those skilled in the art will be omitted as appropriate.
  • the wave number be ⁇ , assuming that the disturbance is Fourier-decomposed in the span direction.
  • the spatial growth rate ⁇ i (x, ⁇ ) of the stationary cross-flow instability can be evaluated by the following equation (2).
  • c 1 , c 2 , c 3 , and c 4 generally take values on the order of these even if the shape of the blade surface is different.
  • a more specific and accurate spatial growth rate- ⁇ i (x, ⁇ ) can be calculated by numerically performing a linear stability analysis of the boundary layer flow.
  • the disturbance amplification factor (gain) can be evaluated by the N value obtained by integrating the above-mentioned spatial growth rate with x (see the following equation (6)).
  • the N value is a function of x and ⁇ .
  • the lower the wavenumber of ⁇ the larger the N value tends to be, which tends to cause a turbulent flow transition.
  • a disturbance with an N of 8 or more may be regarded as an environment that causes a turbulent flow transition.
  • the wave number ⁇ R in the uplift structure 2 is the wave number that excites the disturbance. Therefore, a wave number that gives an N value smaller than the N value selected as the threshold value is selected as the wave number ⁇ R so that the turbulent flow transition does not occur due to the excitation by the raised structure 2.
  • the arrangement position of the raised structure 2 represented by the formula (1) is defined at the center of each unit raised 2a and at a position where the height of the unit raised 2a is equal to the maximum height h.
  • h becomes too high the amplitude of the disturbance excited by the raised structure 2 becomes too large. As a result, a turbulent flow transition may occur immediately.
  • h should be 3 times or less of the boundary layer exclusion thickness ⁇ (x R ) at position x R.
  • h is more preferably about twice as much as ⁇ (x R ), for example, 1.7 times or more and 2.3 times or less.
  • ⁇ (x R ) 0.087 (mm)
  • h / ⁇ (x R ) is 2.874.
  • the width w in the chord direction of the raised structure 2 is not particularly limited as long as it is sufficiently longer than the maximum height h of the raised structure 2.
  • w may be 10 times or more of h, and more preferably 30 times or more.
  • w 4 (mm), so w is 16 times h.
  • the shape of the raised structure 2 is not limited to this.
  • the uneven shape in the chord direction in the raised structure 2 may be represented by a periodic function other than the sine function as long as the dominant wave number is ⁇ R.
  • the periodic function is preferable because the change in the gradient in the chord direction is smooth.
  • the uneven shape in the chord direction in the raised structure 2 may be a shape that is not represented by a specific periodic function as long as it is a periodic uneven shape having a dominant wave number of ⁇ R.
  • the envelope of the concave-convex shape in the span direction in the raised structure 2 is a convex shape in which the height of the concave-convex shape changes gently so as not to excite the disturbance of the wave number other than ⁇ R
  • the Gaussian function can be used.
  • the convex shape is not particularly limited as long as it is a shape that does not easily disturb the flow, such as a streamlined shape or a bell shape, and may be a shape that is not represented by a function.
  • the method of forming the raised structure 2 on the blade 1 is not particularly limited as long as the uneven shape required for the raised structure 2 can be formed.
  • the raised structure 2 may be formed on the upper surface 1a at the time of molding the upper surface 1a.
  • the raised structure 2 may be manufactured as a sheet member having a required uneven shape formed on the surface.
  • the seat member may be attached to the wing 1 on the upper surface 1a by adhesion or the like.
  • Examples of the method for forming the raised structure 2 in the sheet member include mold molding, press molding, and removal processing.
  • the raised structure 2 may be formed on the upper surface 1a by curing the raw material applied on the upper surface 1a into a required uneven shape.
  • the raw material for example, ultraviolet irradiation, heating, or the like may be used.
  • a method for forming the uneven shape a molding die may be used, or an unnecessary portion may be removed after forming the layered portion.
  • the removal processing includes mechanical removal processing, chemical removal processing, laser processing and the like.
  • the raised structure 2 may be formed by adding or depositing a raw material on a sheet member to be attached to the upper surface 1a or the upper surface 1a.
  • the addition / deposition method is not particularly limited, and examples thereof include a 3D printer and an inkjet printing.
  • the material of the raised structure 2 is not particularly limited as long as it is a material having durability under the usage conditions of the wing 1.
  • the material of the raised structure 2 one or more kinds of materials selected from the group consisting of metal, non-metal inorganic material, and organic material may be used.
  • Each unit ridge 2a constituting the ridge structure 2 extends along the direction of the streamline (S in FIG. 1) at the outer edge of the boundary layer.
  • the longitudinal and lateral ends of each unit ridge 2a are smooth curved surfaces that gently increase in height from the top surface 1a.
  • the intermediate portion 2m of each unit ridge 2a is a smooth convex curved surface whose height gradually increases to the top in the longitudinal direction and the lateral direction. Since it has such a shape, the unit ridge 2a is a streamlined type that does not easily obstruct the flow in the mainstream (MF) direction of FIG. 2, and is a ridge that allows the mainstream to pass smoothly.
  • the raised structure 2 extends long in the span direction as a whole in a dome-shaped range in which the envelope is translated in the span direction, and the cross-sectional shape in the span direction is a periodic uneven shape having a wave number ⁇ R. Therefore, by cross flow having a vertical directional component in the main flow (MF) direction (CF) passes through a periodic uneven shape of wavenumber beta R, disturbance of wavenumber beta R is excited. As described above, the raised structure 2 forms a roughness structure that selectively excites the disturbance of the wave number ⁇ R only in the transverse flow component inside the boundary layer at the position of the distance x R from the stagnation point St.
  • FIG. 9 is a graph showing the action of the raised structure according to the embodiment of the present invention.
  • the horizontal axis is the position x in the chord direction, and the vertical axis is max y
  • FIG. 9 is a semi-logarithmic graph whose vertical axis is a logarithmic scale.
  • is a numerical value obtained by discretely Fourier transforming the wall vertical velocity v in the span direction and taking the maximum value of the amplitude in the wall vertical direction (y direction). In the following, max y
  • the target mode the mode of disturbance in which the wave number that artificially excites the disturbance by the raised structure 2 is dominant is called the killer mode.
  • the curve 110 represents the amplitude of the disturbance on the upper surface 1a when the raised structure 2 is not provided.
  • x 20 (mm)
  • the flow transitions to almost turbulent flow near x 250 (mm).
  • the amplitude is reduced over the lower wavenumber amplitude. Therefore, as shown in FIG. 9, the effect of the raised structure 2 can be sufficiently predicted only by examining the amplitude growth of the killer mode, the target mode, and the mode having a lower wave number than the target mode.
  • FIG. 10 is a perspective view showing a raised structure of a comparative example.
  • FIG. 11 is a plan view showing a raised structure of a comparative example.
  • 12 (a) and 12 (b) are graphs showing cross-sectional shapes of the raised structure of the comparative example in the span direction and the chord direction, respectively.
  • x, y, and z are dimensionless at 1 mm.
  • the raised structure 20 of the comparative example includes a substantially columnar unit raised 20a instead of the unit raised 2a in the raised structure 2 of the embodiment.
  • the raised structure 20 has a structure substantially similar to the isolated roughness disclosed in Non-Patent Document 1. However, the specific dimensions are appropriately changed so as to be easily compared with the raised structure 2.
  • the unit ridge 20a is a substantially columnar shape having an average diameter d of 1 mm and a height H of 0.12 mm.
  • the side surface of the unit ridge 20a has a tapered inclination in the projecting direction, and the outer edge of the tip portion is rounded.
  • FIG. 13 shows the result of numerically calculating the amplitude of the disturbance in the blade having the raised structure 20 in the same manner as in the embodiment.
  • FIG. 13 is a graph showing the action of the raised structure of the comparative example.
  • the vertical and horizontal axes of FIG. 13 are the same as the vertical and horizontal axes of FIG.
  • FIG. 13 a curve 110 similar to that in FIG. 9 is plotted for comparison.
  • FIG. 14 is a graph showing the relationship between the height of the raised structure and the killer amplitude.
  • the maximum height h was changed and the same calculation as described above was performed to find the relationship between the maximum height h and the amplitude of the disturbance (killer amplitude) in the killer mode.
  • the horizontal axis is the maximum height h of the raised structure
  • the polygonal lines 120 and 220 show the results of the raised structures 2 and 20, respectively.
  • the killer amplitude can be controlled according to the magnitude of h in the range of the maximum height h of 0 mm or more and 0.16 mm. Further, in the range where h is 0.16 mm or more and 0.25 mm or less, the killer amplitude can be kept constant and the turbulent flow transition can be suppressed.
  • the reason why the turbulent flow transition cannot be suppressed when h exceeds 0.25 mm is that h exceeds 3 times the boundary layer exclusion thickness ⁇ , so that the convex portion of the raised structure 2 extends to the high-speed flow region outside the boundary layer. It is probable that it protruded and immediately caused a turbulent flow transition due to flow separation at the convex portion in the high-speed flow region.
  • h is more preferably 1.7 times or more and 2.3 times or less of ⁇ .
  • the excitation ability of the killer amplitude was lower than that of the raised structure 2, and the effect of suppressing the turbulent transition was also lower.
  • h exceeds 0.10 mm, it has no excitation ability. Since the substantially cylindrical shape is far from the streamlined shape, it is considered that turbulent flow transition due to peeling is caused even if h / ⁇ is less than 2.
  • the turbulent transition position can be shifted to the downstream side in the boundary layer where the cross flow occurs.
  • the turbulent boundary layer in which the frictional resistance becomes significantly larger can be reduced and the region of the laminar boundary layer in which the frictional resistance is smaller can be expanded as compared with the case where the raised structure 2 is not provided.
  • the frictional resistance caused by the lateral flow instability can be reduced.
  • the object provided with the raised structure is not limited to the swept wing of the aircraft as long as it is an object in which cross-flow instability can occur.
  • Examples of an object to which a raised structure can be provided include a wing used for other than an aircraft, a retreating cylinder in flow, a conical rotating body whose diameter increases from the upstream side to the downstream side in the flow direction, and the like.
  • the type of the viscous fluid is not limited as long as the fluid in which the object provided with the raised structure is arranged is a viscous fluid.
  • the viscous fluid may be a compressible fluid or an incompressible fluid.

Abstract

隆起構造は、物体の端部に対して鋭角をなす方向に物体に向かって流れる主流の粘性による摩擦抵抗を低減する。隆起構造は、端部における主流の淀み点よりも下流側における物体の表面に設けられ、表面における主流の流れ方向に延びる複数の隆起体を含む。複数の隆起体のそれぞれの高さは、長手方向および短手方向のいずれの断面においても、両端部から中間部に向かって漸増して前記中間部で最大となる凸状の滑らかな曲線に沿って変化している。複数の隆起体は、淀み点からの距離が一定かつ表面に直交する第１断面において、周期的に変化する第１凹凸形状を形成する。複数の隆起体は、淀み点を連ねる線および第１断面に直交する第２断面において、周期的に変化する凹部および凸部を有し、凸部の各頂点の高さが凸状の包絡線に沿って変化する第２凹凸形状を形成するように、配列されている。

Description

隆起構造および翼
　本発明は、粘性摩擦抵抗を低減する隆起構造および翼に関する。
　本願は、２０１９年３月２９日に日本に出願された特願２０１９－０６９１５３号に基づき優先権を主張し、その内容をここに援用する。
　粘性流体の流れに配置された物体は粘性摩擦抵抗を受ける。例えば、航空機は飛行中に空気の粘性による摩擦抵抗を受ける。航空機が飛行中に受ける全抵抗のうち、約半分は摩擦抵抗であると言われている。摩擦抵抗を低減すれば、航空機における航続距離の伸長、燃費の低減などが可能になるので、粘性摩擦抵抗の低減手段が種々提案されている。
　例えば、特許文献１、２には、翼に複数の突起物を設けることによって、粘性摩擦抵抗を低減することが提案されている。
　特許文献１に記載の技術では、翼の表面に、翼弦方向に長く延び、翼弦方向と直交するスパン方向には実質的に平行なピークヴァレーが繰り返されるシームレスフィルムが設けられている。シームレスフィルムにおける凸部のスパン方向の断面形状は、翼の表面から離れるにつれて先細になっており、先端には鋭くとがった頂部が形成されている。
　特許文献２に記載の技術では、翼弦方向に延びる中央部突起と、スパン方向において中央部突起を挟んで中央部突起と実質的に平行に配置された２つの外側突起と、を有する空力装置が、翼の表面に密集もしくは分散して配置されている。中央部突起と外側突起の間の凹部は、翼弦方向において先端側から後端側に向かって翼の表面から漸次離間する傾斜面からなる。この空力装置のスパン方向における断面形状は、先端に頂部を有する略三角形状である。
　特許文献１、２に記載の技術は、突起物によって、翼の表面に安定した乱流境界層を形成することによって、摩擦低減を図る技術であると考えられる。しかしながら、特許文献１、２に記載の技術は、航空機に広く用いられる後退翼における横流れ不安定性に起因する乱流遷移による摩擦抵抗を低減することはできない。
　非特許文献１には、風洞実験の結果から、翼における淀み点（付着線）の近傍に、微小な孤立粗度（Discrete Roughness Elements）を設けると、横流れ不安定性を抑制できることが報告されている。各孤立粗度の形状は、略円柱状である。非特許文献１に記載の技術では、孤立粗度の配列間隔を適正に設定して横流れ不安定性に起因する乱流遷移の位置を後方にずらすことによって、摩擦低減を図る。
　非特許文献２では、非特許文献１における孤立粗度の作用が数値解析によって解析（評価）されている。さらに、非特許文献２では、解析結果に対する考察に基づいて孤立粗度の改良形状が提案され、数値的な検証結果が示されている。
　非特許文献２における改良型の孤立粗度の形状は、スパン方向の断面においては、正弦曲線を有し、翼弦方向における断面においてはガウス分布を有するので、Sinusoidal Roughness Elements（ＳＲＥ、正弦波状粗度）と称される。
国際公開第２０１２／０８２６６７号 国際公開第２０１８／１６５３１３号
サーリック　外２名、「アメリカン・インスティテュート・オブ・エアロノーティックス・アンド・アストロノーティックス・ペーパー(American Institute of Aeronautics and Astronautics Paper)　Ｎｏ．１９９８－０７８１」（米国）、１９９８、ｐｐ．１－１３ 廣田、井手、林田、服部、「孤立粗度による横流れ不安定性抑制効果の数値的検証」、「数値流体力学シンポジウム講演論文集　ＣＦＤ　２０１８　第３２回」ＤＶＤ－ＲＯＭ、日本流体力学会、２０１８
　しかしながら、上記のような従来技術には、以下のような問題がある。
　特許文献１、２に記載の技術は、上述したように、後退翼に発生する横流れ不安定性に起因する乱流遷移を抑制できないので、後退翼における摩擦抵抗をあまり低減できないという問題がある。特に、特許文献１、２における構成は、凸部が翼弦方向に長く延びているので、横流れの乱流遷移を増長する結果、摩擦抵抗が増大するおそれがある。
　非特許文献２の解析によれば、擾乱の波数が低いほど擾乱の振幅が増大しやすく、飽和状態に達して乱流に遷移しやすい。これに対して、擾乱の波数が高いと飽和時の振幅が低くかつ一旦飽和しても直ちに乱流に遷移することはないという特徴がある。孤立粗度の大きさ、配置ピッチ等を調整することは、横流れにおいて特定の波数成分の擾乱を選択的に励起することによって、より低波数の擾乱の成長を抑制することに相当している。
　このため、非特許文献２に記載されたように、孤立粗度を設ける場合には、孤立粗度によって励起する特定波数の「ｋｉｌｌｅｒ　ｍｏｄｅ（成長率（またはＮ値）の比較的低い横流れモード）によって乱流遷移が起きてしまわないように、細心の注意を払わねばならない。」という問題がある。
　さらに、非特許文献２の解析によれば、非特許文献１の孤立粗度は複数の高波数擾乱を同時に励起してしまい、特定の単一波数の振幅の励起能力が低いため、乱流遷移の抑制効果が低いという問題もある。
　非特許文献２による改良型の正弦波状粗度は、特定波数の振幅の励起能力がより高いため、乱流遷移の抑制効果が向上している。
　ただし、いずれの形状にしても、孤立粗度の高さが一定値を超えるとただちに乱流遷移が起こるので、乱流遷移制御の安定性の観点ではこの一定値を超えない範囲でより大きな抑制効果が得られることが好ましい。この意味で孤立粗度の形状には依然としてさらなる改良の余地がある。
　本発明は、上記のような問題に鑑みてなされたものであり、横流れ不安定性に起因する摩擦抵抗を低減することができる隆起構造および翼を提供することを目的とする。
　上記の課題を解決するために、本発明の第１の態様の隆起構造は、物体の端部に対して鋭角をなす方向に前記物体に向かって流れる主流の粘性による摩擦抵抗を低減する隆起構造であって、前記端部における前記主流の淀み点よりも下流側における前記物体の表面に設けられ、前記表面における前記主流の流れ方向に延びる複数の隆起体を備え、前記複数の隆起体のそれぞれの高さは、長手方向および短手方向のいずれの断面においても、両端部から中間部に向かって漸増して前記中間部で最大となる凸状の滑らかな曲線に沿って変化しており、前記複数の隆起体は、前記淀み点からの距離が一定かつ前記表面に直交する第１断面において、周期的に変化する第１凹凸形状を形成し、前記淀み点を連ねる線（付着線）および前記第１断面に直交する第２断面において、周期的に変化する凹部および凸部を有し、前記凸部の各頂点の高さが凸状の包絡線に沿って変化する第２凹凸形状を形成するように、配列されている。
　本発明の第２の態様の隆起構造によれば、第１の態様において、前記第１凹凸形状は、一定の波数を有する正弦波状であり、前記包絡線は、釣り鐘型の曲線であってもよい。
　本発明の第３の態様の隆起構造によれば、第１または第２の態様において、前記複数の隆起体のそれぞれにおける最大高さは、前記主流によって形成される境界層排除高さの１．７倍以上２．３倍以下であってもよい。
　本発明の第４の態様の翼は、第１～第３のいずれか１つの態様の隆起構造を表面に備える。
　本発明の第１～第３の態様のいずれか１つの態様の隆起構造および第４の態様の翼によれば、横流れ不安定性に起因する摩擦抵抗を低減することができる。
本発明の実施形態の隆起構造および翼の一例を示す模式的な斜視図である。 翼上の三次元境界層の流速分布の一例を示す模式図である。 翼弦方向の位置ｘと、境界層外縁速度Ｕｅ、Ｗｅ、境界層排除厚さδ、およびレイノルズ数Ｒｅδとの関係を示すグラフである。 横流れ不安定性の中立安定曲線を示すグラフである。 横流れ不安定性のＮ値を示すグラフである。 本発明の実施形態の隆起構造の一例を示す斜視図である。 本発明の実施形態の隆起構造の一例を示す平面図である。 本発明の実施形態の隆起構造の一例のスパン方向および翼弦方向の断面形状を示すグラフである。 本発明の実施形態の隆起構造の作用を示すグラフである。 比較例の隆起構造を示す斜視図である。 比較例の隆起構造を示す平面図である。 比較例の隆起構造のスパン方向および翼弦方向の断面形状を示すグラフである。 比較例の隆起構造の作用を示すグラフである。 隆起構造の高さとキラー振幅との関係を示すグラフである。
　以下では、本発明の実施形態の隆起構造および翼について添付図面を参照して説明する。
　図１は、本発明の実施形態の隆起構造および翼の一例を示す模式的な斜視図である。
　図１に示す本実施形態の翼１（物体）は、例えば、航空機の後退翼として用いられる。ただし、図１は、翼１におけるスパン方向の一部のみを図示している。
　翼１は、前縁ｆから後縁ｂに向かって、翼厚が増大し、最大になった後、漸次減少する適宜の翼型を有する。具体的な翼型については特に限定されない。図１に示す例では、前縁ｆと後縁ｂとを結ぶ翼弦線に対して、上面１ａ（表面）の膨らみが下面１ｂの膨らみよりも大きな翼型を有する。
　例えば、翼１は、前縁ｆが後縁ｂよりも高くなる適宜の取り付け角で、図示略の胴体に対して取り付けられている。翼の基準線は、平面視では、胴体の前後軸に対して斜め後方向に傾斜している。このため、平面視では、前縁ｆも、胴体の前後軸に対して斜め後方向に傾斜している。
　図１において、流れＦは、翼１が取り付けられた機体が直進したときに翼１が受ける空気の流れを表す。
　航空機の飛行時には、翼１の先端部（物体の端部）の定位置に、流れＦの淀み点Ｓｔ（二点鎖線参照）が形成される。淀み点Ｓｔは、翼１の先端部において前縁ｆに沿って線状に延びている。
　翼１の先端部において淀み点Ｓｔの下流側の上面１ａ上には、本実施形態の隆起構造２が形成されている。
　隆起構造２は、翼１に向かって流れる流れＦ（主流）の粘性による摩擦抵抗を低減するために設けられている。隆起構造２の配置を決定するための淀み点Ｓｔの位置は、飛行状態によって変化する。隆起構造２の設計にあたっては、特に摩擦抵抗を低減したい飛行状態に形成される淀み点Ｓｔが用いられる。例えば、航空機が水平飛行する際の淀み点Ｓｔが用いられてもよい。
　隆起構造２は、複数の単位隆起２ａ（隆起体）を備える。
　各単位隆起２ａは、上面１ａ上の流れＦの流れ方向（図１に模式的に示す流線Ｓ参照）に延びている。上面１ａ上の任意の位置の流線Ｓは、隆起構造２が除去された翼１の形状、流れＦの速度、流れＦの粘性等が与えられれば、数値解析によって求められる。
　各単位隆起２ａの高さは、長手方向および短手方向のいずれの断面においても、両端部から中間部に向かって漸増して中間部で最大となる凸状の滑らかな曲線に沿って変化している。例えば、各単位隆起２ａの高さは、長手方向の断面においては、上流側の先端２ｆ側および下流側の後端２ｂ側から、先端２ｆおよび後端２ｂの間の中間部２ｍに向かって、漸増して中間部２ｍで最大となる凸状の滑らかな曲線に沿って変化している。ここで、中間部２ｍは、単位隆起２ａの長手方向の中央であってもよいが、長手方向の中央からずれていてもよい。図示は省略するが、短手方向の断面も同様である。
　このように、単位隆起２ａにおいて隆起方向の先端面は丸みを帯びた滑らかな曲面である。単位隆起２ａの延在方向は、短手方向の断面における頂部を連ねた尾根線で表される。尾根線は、流れＦの流れ方向に延びることがより好ましい。しかし単位隆起２ａの尾根線は短手方向の断面においては滑らかな凸部の停留点になっているので、尾根線は測定誤差、製作誤差等を考慮して略流れ方向に延びていればよい。
　図１に示す例では、隆起構造２の各単位隆起２ａは、淀み点Ｓｔからの距離が一定の線上に等間隔に配列されている。これにより、隆起構造２は、配列方向に沿う断面（第１断面）では、周期的に変化する第１凹凸形状を形成している。隆起構造２は、淀み点Ｓｔを連ねる線および第１断面に直交する第２断面においては、各凸部の頂点の高さが凸状の包絡線に沿って変化する第２凹凸形状を形成している。
　隆起構造２の具体例については後述する。
　隆起構造２の具体例について説明する前に、翼１上の三次元境界層と横流れ不安定性について説明する。
　図２は、翼上の三次元境界層の流速分布の一例を示す模式図である。図３は、翼弦方向の位置ｘと、境界層外縁速度Ｕｅ、Ｗｅ、境界層排除厚さδ、およびレイノルズ数Ｒｅδとの関係を示すグラフである。図３において、横軸ｘは、淀み点Ｓｔからｘ方向の正方向に進む長さ［ｍｍ］を表す。縦軸は、境界層外縁速度Ｕｅ、Ｗｅ、境界層排除厚さδ、およびレイノルズ数Ｒｅδを示す。ただし、境界層排除厚さδは、１ｍｍで無次元した数値で示されている。図４は、横流れ不安定性の中立安定曲線を示すグラフである。図５は、横流れ不安定性のＮ値を示すグラフである。
　図１に示すように、以下の説明では、淀み点Ｓｔを連ねた直線に沿う方向をｚ方向と称する。ｚ方向は翼１のスパン方向でもある。ｚ方向の正方向は、図示略の胴体から翼１の長手方向の先端に向かう方向である。
　翼１の淀み点Ｓｔを連ねた直線に直交し、上面１ａに沿う方向をｘ方向と称する。特に誤解の生じるおそれがない場合には、ｘ方向を単に翼弦方向と称する場合もある。ｘ方向の正方向は、淀み点Ｓｔから上面１ａに沿って後縁ｂに向かう方向である。
　ｚ方向およびｘ方向に直交する方向をｙ方向と称する。ｙ方向は、単位隆起２ａの隆起方向である。
　なお、隆起構造２の単位隆起２ａの長さおよび配列ピッチは、後述するように、翼１の翼弦方向の幅およびスパン方向の長さに比べて充分短い。このため、隆起構造２の形成範囲において、上面１ａは、解析上、平面で近似できる。
　図２に、翼１の先端部における上面１ａ上に形成される三次元境界層の流速分布の一例を示す。図２において、Ｕ軸は、翼弦方向の流速、Ｗ軸はスパン方向の流速を表す。いずれの流速も淀み点Ｓｔ（ｘ＝０）における音速３１１ｍ／ｓで規格化（無次元化）されている。図２においてｙ／δ軸は、境界層排除厚さδで規格化されたｙ方向の位置である。
　流れＦは、翼１の後退角に応じて淀み点Ｓｔに対して斜め方向に入射するので、主流の境界層外縁速度は、翼弦方向およびスパン方向にそれぞれ、Ｕｅ、Ｗｅの速度成分を有する。境界層の外部では、Ｕｅ、Ｗｅ以外の速度成分を有しない。
　これに対して、境界層の内部では、圧力勾配の方向である翼弦方向と主流の方向が平行でないため、主流（ＭＦ）の流れ方向に交差する方向に横流れ（ＣＦ）が発生する。横流れの速度分布は、上面１ａの表面と、境界層外縁と、において０であり、その中間部で最大値を取る、山形の分布である。
　このため、境界層の内部における合成流速分布は、主流の流れに沿う断面に対してねじれている。
　流れＦは、上面１ａに沿って進むにつれて、翼１の翼型の形状に応じて、翼弦方向に増速する。図３に、横流れを有する場合のｘ方向における境界層外縁速度Ｕｅ、Ｗｅ、境界層排除厚さδ、およびレイノルズ数Ｒｅδの変化の解析例を示す。ここで、レイノルズ数Ｒｅδは、Ｒｅδ＝｛√（Ｕｅ ２＋Ｗｅ ２）｝δ／ν０で定義される。ただし、ν０は、淀み点Ｓｔにおける空気（流体）の動粘性係数である。
　図３のグラフの詳しい計算方法については、非特許文献２に記載された通りである。
　翼弦方向の境界層外縁速度Ｕｅは、淀み点Ｓｔ（ｘ＝０）からｘが増大するにつれて、急峻に増大して、次第に飽和に向かう。スパン方向の境界層外縁速度Ｗｅは、ｘによらず一定である。
　境界層排除厚さδは、初期に急増した後、略直線的に増大する。このため、翼弦方向における後端ほど、境界層の厚さが増大することが分かる。一方、レイノルズ数Ｒｅδはｘが増大するにつれて、略直線的に増大する。このため、翼１の後端側では層流境界層が乱流遷移を起こしやすいことが分かる。
　乱流遷移は、上面１ａにおける付着物、翼面自体の微小な凹凸などによって励起されると考えられる。擾乱を局所的に平面波で近似することによって、図４に示す曲線４０で示すような横流れ不安定性の中立安定曲線が求められる。
　すなわち、擾乱の定在モードを考えて、擾乱を局所的に平面波（∝ｅｘｐ（ｉαｘ＋ｉβｚ））で近似する。ここで、ｉは虚数単位、αは翼弦方向波数（複素数）、βはスパン方向波数（実数）である。スパン方向波数βを与えると、翼弦方向波数αが複素固有値として求まる。
　図４において、横軸は、翼弦方向の位置ｘであり、縦軸は、局所平面波近似した擾乱におけるスパン方向波数βである。
　曲線４０は、擾乱において複素数で表される翼弦方向波数α（＝αｒ＋ｉαｉ）の複素数成分αｉが０となる曲線である。このため、曲線４０の外側（「Ｓｔａｂｌｅ（αｉ＞０）」と記載の領域）は、αｉ＞０になるため、安定領域である。これに対して、曲線４０の内側（「Ｕｎｓｔａｂｌｅ（αｉ＜０）」と記載の領域）は、αｉ＜０になるため、不安定領域である。
　図４によれば、波数βが高い場合には、翼弦方向において前縁ｆから遠ざかれば安定領域になるのに対して、波数βが低くなるほど、不安定領域が後端に延び、不安定領域が拡大している。
　図５に、スパン方向波数β（ｍｍ－１）に対するＮ値（Ｎ　ｆａｃｔｏｒ）をプロットしたグラフを示す。曲線５１、５２、５３は、それぞれ、ｘ＝１００、２０２、３０３（ｍｍ）におけるＮ値を示す。Ｎ値は、擾乱の増幅率（ゲイン）をｅＮとすることによって求められる。Ｎ値が大きいほど乱流遷移が起こりやすい。
　図５に曲線５１、５２、５３で示されるように、ｘが小さいほど高波数成分によって擾乱が増幅され、ｘが大きいほど低波数成分によって擾乱が増幅されることが分かる。さらに、増幅率の大きさは、低波数成分ほど大きくなっている。例えば、ｘ＝２０２（ｍｍ）ではβ＝３．３３３（ｍｍ－１）による増幅が支配的なモードであり、Ｎ値は８を超えている。
　したがって、擾乱を形成するエネルギーが高波数成分に分配されれば、Ｎ値を増大させることなく、低波数成分の擾乱の成長を抑制できると考えられる。この場合、高波数の擾乱の振幅は増大する。しかし、高波数の擾乱は、成長しても低波数の擾乱より乱流遷移を起こしにくいので、境界層全体としては乱流遷移が起こりにくくなる。
　非特許文献１、２には、特に抑制したい擾乱のスパン方向波数βの１．５倍程度のスパン方向波数βの擾乱を励起しやすい間隔で孤立粗度を配列すると、乱流遷移位置をより後方にずらすことができることが報告されている。
　本発明者は、孤立粗度による擾乱の励起をより効率化できれば、抑制対象の擾乱をさらに効果的に抑制できるのではないかと考えて、鋭意研究し、本発明に到った。
　次に、隆起構造２の具体例について詳細に説明する。
　図６は、本発明の実施形態の隆起構造の一例を示す斜視図である。図７は、本発明の実施形態の隆起構造の一例を示す平面図である。図８における（ａ）、（ｂ）は、それぞれ本発明の実施形態の隆起構造の一例のスパン方向および翼弦方向の断面形状を示すグラフである。図６～８において、ｘ、ｙ、ｚは１ｍｍで無次元化されている。
　図６、７に示す隆起構造２の形状ｙ＝ｙＲ（ｘ，ｚ）は、下記式（１）～（５）で表される。ただし、式（１）におけるｚの定義域は、翼１の翼長の範囲で適宜設定される。この場合、隆起構造２のｚ方向における定義域の端縁には、さらにｙが漸次０に向かって滑らかに減少する端部が接続されることがより好ましい。
Figure JPOXMLDOC01-appb-M000001
　ここで、ｈは隆起構造２の最大高さ、βＲはスパン方向波数（波長は、２π／βＲ）、ｘＲは翼弦方向における隆起構造２の位置、ｗは翼弦方向における隆起構造２の片幅である。Ｕｅ（ｘＲ）およびＷｅ（ｘＲ）は、それぞれＵｅ（ｘ）およびＷｅ（ｘ）におけるｘ＝ｘＲの値である。
　式（１）において、翼弦方向の位置ｘを固定して得られる翼弦方向の断面、すなわち淀み点Ｓｔからの距離が一定かつ上面１ａに直交する第１断面では、ｙＲはｚのみの関数であって、波数βＲの正弦波（第１凹凸形状）である。
　例えば、図８の（ａ）に曲線１０１で示すように、ｘ＝ｘＲでは、ｙＲは、０～ｈまで波数βＲで振動する正弦波である。ｘ≠ｘＲでは、（ｘ－ｘＲ）に応じてｙの最大値がｈ未満に変化する以外は、同様の正弦波である。
　式（１）において、スパン方向の位置ｚを固定して得られるスパン方向の断面、すなわち淀み点Ｓｔ連ねる線および第１断面に直交する第２断面では、ｙＲはｘのみの関数であって、波数｛βＲ×（Ｗｅ（ｘＲ）／Ｕｅ（ｘＲ））｝の正弦波成分と、中心がｘＲ、翼弦方向における片幅ｗのガウス関数成分と、の積である。ただし、片幅ｗは、高さｅｘｐ（－π）（＝０．０４３２）における片幅を表す。
　このため、ｙＲは、第２断面において、波数｛βＲ×（Ｗｅ（ｘＲ）／Ｕｅ（ｘＲ））｝でそれぞれ周期的に変化する凹部および凸部（正弦波状の凹凸部）を有し、各凸部の頂点の高さが凸状の包絡線（ガウス関数）に沿って変化する第２凹凸形状を形成している。
　例えば、図８の（ｂ）に曲線１０２で示すように、ｚ＝０とすると、ｙＲは、ｘ＝ｘＲから幅方向の両端部に向かって、曲線１０３で表されるガウス関数に沿って振幅が減衰し、波数｛βＲ×（Ｗｅ（ｘＲ）／Ｕｅ（ｘＲ））｝でそれぞれ周期的に変化する凹部および凸部を形成している。図８の（ｂ）に示す例では、ガウス関数の片幅はｗ＝４（ｍｍ）である。
　式（１）におけるパラメータβＲ、ｘＲ、ｈ、ｗは、翼弦方向における乱流遷移位置を下流側にずらすことによって、上面１ａ上の層流境界層の長さを長くする目的で設定される。
　以下に、各パラメータの設定方法の一例について簡単に説明する。なお、横流れ不安定性解析において当業者に周知の事項、手法については適宜説明を省略する。
　擾乱をスパン方向にフーリエ分解したとし、波数をβとする。各βに応じて、定在横流れ不安定性の空間的成長率－αｉ（ｘ，β）が次式（２）のように評価できる。
Figure JPOXMLDOC01-appb-M000002
　ここで、ｃ１、ｃ２、ｃ３、ｃ４は、翼面の形状に依存した無次元定数であり、例えば、ｃ１＝０．４５、ｃ２＝１．５、ｃ３＝１５、ｃ４＝０．８３である。ただし、ｃ１、ｃ２、ｃ３、ｃ４は、翼面の形状が異なっていても、概ねこのようなオーダーの値をとる。
　より具体的で正確な空間的成長率－αｉ（ｘ，β）は、境界層流れの線形安定性解析を数値的に行えば計算することができる。
　擾乱の増幅率(ゲイン) は、上述の空間成長率をｘで積分したＮ値（下記式（６）参照）で評価できる。Ｎ値は、ｘ、βの関数である。
Figure JPOXMLDOC01-appb-M000003
　一般に低波数のβほどＮ値が大きい傾向があり、乱流遷移を引き起こしやすい。
　環境のノイズレベルによっても異なるが、例えば、Ｎが８以上の擾乱が乱流遷移を引き起こす環境であると見なしてもよい。例えば、航空機の場合、乱流遷移を引き起こすＮ値の閾値は、８～１０であると言われている。したがって、航空機の場合、Ｎ＝８は、安全側の見積もりである。
　隆起構造２における波数βＲは、擾乱を励起する波数である。このため、隆起構造２による励起によって乱流遷移が起こらないように、波数βＲとしては、閾値に選定されたＮ値よりも小さいＮ値を与える波数が選ばれる。例えば、閾値がＮ＝８の場合、波数βＲの選定用のＮ値は、例えば、Ｎ＝６が用いられてもよい。
　この場合、式（６）の左辺にＮ＝６を代入すれば、βＲが求められる。
　次に、β＝βＲとした空間的成長率－αｉ（ｘ，βＲ）が最大になる位置ｘを、隆起構造２の配置位置ｘＲとする。式（１）で表される隆起構造２の配置位置は、各単位隆起２ａの中心であり、かつ単位隆起２ａの高さが最大高さｈに等しくなる位置で規定される。
　隆起構造２の最大高さｈは、高いほど隆起構造２による擾乱を励起しやすくなるので、低波数の擾乱を抑制しやすくなる。しかし、ｈが高くなりすぎると、隆起構造２によって励起される擾乱の振幅が大きくなりすぎる。この結果、ただちに乱流遷移が発生するおそれがある。そのため、ｈは、位置ｘＲにおける境界層排除厚さδ（ｘＲ）の３倍以下とするべきである。安全側に見積もると、ｈはδ（ｘＲ）の２倍程度、例えば、１．７倍以上２．３倍以下であるとより好ましい。
　例えば、図８の（ａ）、（ｂ）に示す例では、δ（ｘＲ）＝０．０８７（ｍｍ）とすると、ｈ／δ（ｘＲ）は、２．８７４である。
　隆起構造２における翼弦方向の片幅ｗは、隆起構造２の最大高さｈよりも充分に長ければ特に限定されない。例えば、ｗは、ｈの１０倍以上あればよく、３０倍以上だとより好ましい。
　例えば、図８の（ｂ）に示す例では、ｗ＝４（ｍｍ）であるので、ｗはｈの１６倍である。
　以上、隆起構造２の具体例について説明したが、隆起構造２の形状は、これには限定されない。
　例えば、隆起構造２における翼弦方向の凹凸形状は、卓越波数がβＲであれば、正弦関数以外の周期関数で表されてもよい。ただし、周期関数は、翼弦方向における勾配の変化が滑らかであることより好ましい。
　さらに、隆起構造２における翼弦方向の凹凸形状は、卓越波数がβＲである周期的な凹凸形状であれば、特定の周期関数で表されない形状であってもよい。
　例えば、隆起構造２におけるスパン方向の凹凸形状の包絡線は、βＲ以外の波数の擾乱を励起しにくいように、凹凸形状の高さがなだらかに変化する凸形状であれば、ガウス関数には限定されない。例えば、凸形状は、流線型、釣り鐘型など、流れを擾乱しにくい形状であれば、具体的な関数形は特に限定されず、関数で表されない形状であってもよい。
　翼１上における隆起構造２の形成方法は、隆起構造２として必要な凹凸形状を形成できれば特に限定されない。
　例えば、隆起構造２は、上面１ａの成形時に上面１ａに形成されてもよい。
　例えば、隆起構造２は、必要な凹凸形状が表面に形成されたシート部材として製造されてもよい。シート部材は、上面１ａ上に接着などによって、翼１に取り付けられてもよい。シート部材における隆起構造２の形成方法としては、モールド成形、プレス成形、除去加工などが挙げられる。
　例えば、隆起構造２は、上面１ａ上に塗布された原料を必要な凹凸形状に硬化させることによって、上面１ａ上に形成されてもよい。原料の硬化方法としては、例えば、紫外線照射、加熱などが用いられてもよい。凹凸形状の形成方法としては、成形型が用いられてもよいし、層状部を形成した後に不要部が除去加工されてもよい。例えば、除去加工としては、機械的除去加工、化学的除去加工、レーザ加工などが挙げられる。
　例えば、隆起構造２は、原料を上面１ａまたは上面１ａに貼り付けるシート部材上に付加または堆積することによって形成されてもよい。付加・堆積方法としては、特に限定されないが、例えば、３Ｄプリンタ、インクジェット印刷などが挙げられる。
　隆起構造２の材料としては、翼１の使用条件において耐久性を有する材料であれば、特に限定されない。例えば、隆起構造２の材料としては、金属、非金属無機材料、および有機材料からなる群から選ばれた１種類以上の材料が用いられてもよい。
　次に、隆起構造２の作用について説明する。
　隆起構造２を構成する各単位隆起２ａは、境界層外縁における流線（図１のＳ）の方向に沿って延びている。各単位隆起２ａの長手方向および短手方向の端部は、上面１ａからなだらかに高さが増大する滑らかな曲面である。さらに、各単位隆起２ａの中間部２ｍは、長手方向および短手方向において、頂部まで漸次高さが増大する滑らかな凸曲面である。
　このような形状を有するため、単位隆起２ａは、図２の主流（ＭＦ）方向の流れを阻害しにくい流線型であり、主流を円滑に通過させることができる隆起である。
　隆起構造２を全体として見ても、単位隆起２ａがスパン方向に周期的に整列しているため、同じく主流方向に通過する流れを阻害しにくい。
　隆起構造２は、包絡線をスパン方向に平行移動したドーム状の範囲で全体としてスパン方向に長く延びており、スパン方向の断面形状は、波数βＲの周期的な凹凸形状からなる。このため、主流（ＭＦ）方向に垂直な方向成分を有する横流れ（ＣＦ）が波数βＲの周期的な凹凸形状を通過することによって、波数βＲの擾乱が励起される。
　このように、隆起構造２は、淀み点Ｓｔからの距離ｘＲの位置において、境界層内部の専ら横流れ成分に波数βＲの擾乱を選択的に励起する粗度構造を形成している。
　次に、波数βＲの擾乱の成長が乱流遷移を抑制する作用について、数値解析結果に基づいて説明する。
　図９は、本発明の実施形態の隆起構造の作用を示すグラフである。図９において、横軸は翼弦方向の位置ｘ、縦軸は擾乱の振幅に対応するｍａｘｙ｜ｖｎ｜（ただし、ｖにおけるハット記号（サーカムフレックスアクセント）は省略。以下も同じ。）を示す。なお、図９は縦軸が対数目盛とされた片対数グラフである。
　図９に示すグラフは、ｘ＝２０（ｍｍ）の位置において湧き出し/吸い込みによる擾乱源を配置し、適当な初期振幅で、擾乱の波数のモードを励起した結果を示す。縦軸のｍａｘｙ｜ｖｎ｜は、壁垂直方向速度ｖをスパン方向に離散フーリエ変換し、壁垂直方向（ｙ方向）に関する振幅の最大値を取った数値である。以下では、ｍａｘｙ｜ｖｎ｜を擾乱の振幅と称する。
　スパン方向波数βは、β＝１．６６６ｎ（ｍｍ－１）（ただし、ｎ＝１，２，…）のように整数ｎでラベル付けして表している。また、グラフの見やすさを考慮して、図９には、ｎ＝１，２，３のみの結果を示している。
　計算に用いた隆起構造２の形状は、上記式（１）において、βＲ＝５．０００（ｍｍ－１）（ｎ＝３）、Ｗｅ（ｘＲ）／Ｕｅ（ｘＲ）＝０．５７、ｈ＝０．２５（ｍｍ）、ｗ＝４（ｍｍ）、ｘＲ＝３０（ｍｍ）とされた形状である。
　ここで、βＲ＝５．０００（ｍｍ－１）としているのは、Ｎ値が大きいβ＝３．３３３（ｍｍ－１）の擾乱を抑制するためである。本例では、スパン方向波数βＲは、擾乱の抑制対象のモード（ターゲットモードと称する）の波数β＝３．３３３（ｍｍ－１）の１．５倍である。
　ターゲットモードを抑制するために、隆起構造２によって人為的に擾乱を励起する波数が支配的な擾乱のモードをキラーモードと称する。
　図９において曲線１１０は、隆起構造２を設けない場合の上面１ａにおける擾乱の振幅を表す。
　曲線１１０に示すように、隆起構造２を設けない場合、擾乱の振幅は、ｘ＝２０（ｍｍ）から、横流れ不安定性によって指数関数的に増大し、ｘ＝１８０（ｍｍ）の付近で飽和している。さらにｘが増加すると、ｘ＝２００（ｍｍ）の付近から二次不安定性が発生し始め、ｘ＝２５０（ｍｍ）の付近ではほぼ乱流に遷移していることが分かる。ただし、ここでは二次不安定性の引き金として、ｘ＝１８０（ｍｍ）の位置にて微小な振動擾乱を数値的に与えている。
　図９において、曲線１１１、１１２、１１３は、それぞれｎ＝１、２、３のスパン方向波数βに対応する擾乱の振幅のフーリエ成分の変化を表す。
　曲線１１３に示すように、ｎ＝３の場合、ｘ＝３０（ｍｍ）を中心として配置された隆起構造２の形成領域で、振幅が急峻に増大して最大値を取り、この後、ｘが増加するにつれて、振幅が減少している。特に、ｘ＝２５０（ｍｍ）の付近では、振幅は略０である。さらにｘが増加すると振幅は増大するが、ｘ＝４００（ｍｍ）でも乱流遷移は発生していない。
　曲線１１３の最大値は、曲線１１０における乱流遷移時の振幅を上回っているにもかかわらず乱流遷移が発生していない。この理由は、曲線１１０における乱流遷移はβ＝５．０００（ｍｍ－１）よりも低い波数の擾乱の成長によって形成されているためであると考えられる。
　曲線１１２に示すターゲットモード（ｎ＝２）の場合、ｘ＝２４０（ｍｍ）の付近までは、擾乱の振幅が略０に抑制されている。さらにｘが増加すると、振幅は指数関数的に増大するが、ｘ＝４００（ｍｍ）でも乱流遷移は発生していない。
　曲線１１１に示すターゲットモードよりも低波数（ｎ＝１）の場合、ｘが増加するにつれて指数関数的に増大している。ターゲットモードに比べると振幅の抑制効果は少ないが、ｘ＝４００（ｍｍ）でも乱流遷移は発生していない。
　図９における図示は省略するが、ｎが４以上のモードの場合、３の倍数のモードではｎ＝３と同様に擾乱が励起され、その他のモードでは振幅が抑制される。いずれのモードでもより低波数の振幅よりは振幅が低減される。このため、図９に示すように、キラーモード、ターゲットモード、ターゲットモードよりも低波数のモードの振幅成長を検討するだけでも充分に隆起構造２の効果を予測できる。
　次に本実施形態の隆起構造２の作用について比較例と対比して説明する。
　図１０は、比較例の隆起構造を示す斜視図である。図１１は、比較例の隆起構造を示す平面図である。図１２における（ａ）、（ｂ）は、それぞれ比較例の隆起構造のスパン方向および翼弦方向の断面形状を示すグラフである。
　図１０～１２において、ｘ、ｙ、ｚは１ｍｍで無次元化されている。
　図１０、１１に示すように、比較例の隆起構造２０は、実施形態の隆起構造２における単位隆起２ａに代えて、略円柱状の単位隆起２０ａを備える。以下、上述の隆起構造２と異なる点を中心に説明する。
　隆起構造２０は、非特許文献１に開示の孤立粗度と略同様の構成を有する。ただし、具体的な寸法は、隆起構造２と対比しやすいように適宜変更している。
　隆起構造２０によるターゲットモードは、実施形態と同様、β＝３．３３３（ｍｍ－１）である。
　図１２の（ａ）に曲線２０１で示すように、単位隆起２０ａは、平均直径ｄが１ｍｍ、高さＨが０．１２ｍｍの略円柱状である。ただし、単位隆起２０ａの側面は突出方向に先細の傾斜を有しており、先端部の外縁には丸みが付けられている。
　図１２の（ｂ）に示すように、単位隆起２０ａは、ｘ＝ｘＲにおいて、スパン方向に等ピッチの１列をなして配列されている。単位隆起２０ａのスパン方向の配列ピッチｐ（図１２の（ａ）参照）は、β＝５．０００（ｍｍ－１）の擾乱を励起するため、２π／５．０００ｍｍである。
　実施形態と同様にして、隆起構造２０を有する翼における擾乱の振幅を数値計算した結果を図１３に示す。
　図１３は、比較例の隆起構造の作用を示すグラフである。図１３の縦軸および横軸は、図９の縦軸および横軸と同様である。
　図１３には、比較のために、図９におけると同様の曲線１１０がプロットされている。
　曲線２１１、２１２、２１３は、それぞれｎ＝１、２、３のスパン方向波数βに対応する擾乱の振幅のフーリエ成分の変化を表す。
　曲線２１３に示すように、ｎ＝３のキラーモードの場合、ｘ＝３０（ｍｍ）を中心として配置された隆起構造２０の形成領域で、振幅が急峻に増大して最大値を取り、この後、急峻に振幅が低下した後、振幅は増加に転じた。ｎ＝３のキラーモードの場合、振幅は、ｘ＝１００（ｍｍ）の近辺で極大、ｘ＝２００（ｍｍ）の近辺で極小、ｘ＝２６０（ｍｍ）の近辺で極大となる２峰性の変化を示した。この後、振幅は、ｘ＝３００（ｍｍ）の近辺で不安定になった後、急峻に増大した。ｘ＝３３０（ｍｍ）の近辺では、乱流に遷移したことが分かる。
　曲線２１２に示すターゲットモード（ｎ＝２）の場合、ｘ＝１６０（ｍｍ）付近までは、擾乱の振幅が指数関数的に増大し、極大となった後、ｘ＝２００（ｍｍ）の付近で極小となった。ｘがさらに増加すると、振幅が指数関数的に増大して、ｘ＝２５０（ｍｍ）付近で飽和した。この後、キラーモードと同様、振幅は、ｘ＝３００（ｍｍ）の近辺で不安定になった後、急峻に増大した。ｘ＝３３０（ｍｍ）の近辺では、乱流に遷移したことが分かる。
　曲線２１１に示す、ターゲットモードよりも低波数（ｎ＝１）の場合、ターゲットモードと同様に、ｘが増加するにつれて指数関数的に増大した。ただし、曲線２１１に示すように、振幅は、ｘ＝１６０（ｍｍ）を超えても増大し、ｘ＝３００（ｍｍ）の付近で飽和した後、乱流に遷移したことが分かる。
　比較例では、ｎ＝１、２、３のいずれもモードでも、ｘ＝３３０（ｍｍ）の付近では乱流遷移していた。このため、隆起構造２０の効果は、乱流遷移位置をｘ＝２５０（ｍｍ）からｘ＝３３０（ｍｍ）まで、８０ｍｍ後退させることができる程度である。
　これに対して、本実施形態の隆起構造２を有する計算結果によれば、少なくとも計算範囲のｘ＝４００（ｍｍ）まで乱流遷移が発生しておらず、発生する兆候もなかった。したがって、隆起構造２は、乱流遷移位置は１５０ｍｍ以上後退させることができている。
　次に、隆起構造２の最大高さｈの検討結果について説明する。
　図１４は、隆起構造の高さとキラー振幅との関係を示すグラフである。
　隆起構造２、２０において、最大高さｈを変えて、上述と同様の計算を行い、最大高さｈと、キラーモードにおける擾乱の振幅（キラー振幅）との関係を求めた結果を図１４に示す。
　図１４における横軸は隆起構造の最大高さｈ、縦軸はｘ＝４０（ｍｍ）におけるキラーモードにおける擾乱の振幅である。なお、本解析において、ｘＲ＝３０（ｍｍ）における境界層排除厚さδは０．０８７ｍｍであった。
　折れ線１２０、２２０は、それぞれ隆起構造２、２０の結果を示す。
　折れ線１２０に示すように、隆起構造２では、最大高さｈが増大するにつれてキラー振幅は、増大し、ｈ＝０．１６（ｍｍ）の近辺で略一定になった。しかし、ｈ＝０．２５（ｍｍ）の近辺からキラー振幅はさらに増大した。ただし、ｈが０．３０（ｍｍ）以上では、ｘＲ＝３０（ｍｍ）より後方からただちに乱流領域Ｔに入っていた。
　折れ線２２０に示すように、隆起構造２０では、ｈ＝０．０７（ｍｍ）までは、最大高さｈの増加とともにキラー振幅が単調に増加した。しかし、変化率は、隆起構造２に比べると格段に小さかった。さらに、キラー振幅は、ｈ＝０．０７（ｍｍ）で略飽和し、ｈ＝０．１０（ｍｍ）を超えると急峻に減少した。
　隆起構造２０を有する場合、ｈが０．１０（ｍｍ）を超えると、キラーモードはほとんど励起されなかった。さらに、ｈ＝０．１５（ｍｍ）以上になると、やはりｘＲ＝３０（ｍｍ）より後方からただちに乱流領域Ｔに入っていた。
　実施形態の隆起構造２の場合、最大高さｈが０ｍｍ以上０．１６ｍｍの範囲では、ｈの大きさに応じてキラー振幅を制御できる。さらにｈが０．１６ｍｍ以上０．２５ｍｍ以下の範囲では、キラー振幅を一定に保ち、かつ乱流遷移を抑制できる。
　ｈが０．２５ｍｍを超えると、乱流遷移が抑制できない理由としては、ｈが境界層排除厚さδの３倍を超えたため、隆起構造２の凸部が境界層外部の高速流領域にまで突出し、高速流領域における凸部において流れの剥離による乱流遷移をただちに引き起こしたことなどが考えられる。流れの剥離が起きると乱流遷移を逆に促進してしまうので、安全側に見積もるならば、ｈは、δの１．７倍以上２．３倍以下であることがより好ましい。
　比較例の隆起構造２０の場合、隆起構造２に比べるとキラー振幅の励起能力が低く、乱流遷移の抑制効果も低かった。例えば、ｈが０．１０ｍｍを超えると励起能力を有していなかった。略円柱形状が流線型とはかけ離れているために、ｈ／δが２に満たない程度でも剥離による乱流遷移を引き起こすようになると考えられる。
　以上説明したように、隆起構造２および翼１によれば、横流れが発生する境界層において、乱流遷移位置を下流側にずらすことができる。これにより、隆起構造２を有しない場合に比べて、摩擦抵抗が格段大きくなる乱流境界層を縮小し、より摩擦抵抗が小さい層流境界層の領域を拡大できる。
　隆起構造２および翼１によれば、横流れ不安定性に起因する摩擦抵抗を低減することができる。
　なお、上記実施形態では、隆起構造が航空機の後退翼に設けられる場合の例で説明した。しかし、隆起構造が設けられる物体は、横流れ不安定性が発生しうる物体であれば、航空機の後退翼には限定されない。
　隆起構造を設けることができる物体の例としては、航空機以外に用いられる翼、流れ中の後退円柱、流れ方向の上流側から下流側に向かって拡径する円錐状の回転体などが挙げられる。
　さらに、隆起構造が設けられる物体が配置される流体は、粘性流体であれば、粘性流体の種類は限定されない。粘性流体は圧縮性流体であってもよいし、非圧縮性流体でもよい。
　以上、本発明の好ましい実施形態を説明したが、本発明はこの実施形態に限定されることはない。本発明の趣旨を逸脱しない範囲で、構成の付加、省略、置換、およびその他の変更が可能である。
　また、本発明は前述した説明によって限定されることはなく、添付の特許請求の範囲によってのみ限定される。
　上記実施形態によれば、横流れ不安定性に起因する摩擦抵抗を低減することができる隆起構造および翼を提供できる。
１　翼（物体）
１ａ　上面（表面）
１ｂ　下面
２　隆起構造
２ａ　単位隆起（隆起体）
ｆ　前縁
Ｆ　流れ
Ｓｔ　淀み点

Claims (4)

  1. 　物体の端部に対して鋭角をなす方向に前記物体に向かって流れる主流の粘性による摩擦抵抗を低減する隆起構造であって、
    　前記端部における前記主流の淀み点よりも下流側における前記物体の表面に設けられ、前記表面における前記主流の流れ方向に延びる複数の隆起体を備え、
    　前記複数の隆起体のそれぞれの高さは、
    　長手方向および短手方向のいずれの断面においても、両端部から中間部に向かって漸増して前記中間部で最大となる凸状の滑らかな曲線に沿って変化しており、
    　前記複数の隆起体は、
    　　前記淀み点からの距離が一定かつ前記表面に直交する第１断面において、周期的に変化する第１凹凸形状を形成し、
    　　前記淀み点を連ねる線および前記第１断面に直交する第２断面において、周期的に変化する凹部および凸部を有し、前記凸部の各頂点の高さが凸状の包絡線に沿って変化する第２凹凸形状を形成するように、
    　配列されている、
    隆起構造。
  2. 　前記第１凹凸形状は、
    　一定の波数を有する正弦波状であり、
    　前記包絡線は、
    　釣り鐘型の曲線である、
    請求項１に記載の隆起構造。
  3. 　前記複数の隆起体のそれぞれにおける最大高さは、
    　前記主流によって形成される境界層排除高さの１．７倍以上２．３倍以下である、
    請求項１または２に記載の隆起構造。
  4. 　請求項１～３のいずれか１項に記載の隆起構造を表面に備える、翼。
PCT/JP2020/011911 2019-03-29 2020-03-18 隆起構造および翼 WO2020203284A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2021511412A JP7465483B2 (ja) 2019-03-29 2020-03-18 隆起構造および物体
US17/599,092 US20220161923A1 (en) 2019-03-29 2020-03-18 Raised structure and wing

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2019-069153 2019-03-29
JP2019069153 2019-03-29

Publications (1)

Publication Number Publication Date
WO2020203284A1 true WO2020203284A1 (ja) 2020-10-08

Family

ID=72668750

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2020/011911 WO2020203284A1 (ja) 2019-03-29 2020-03-18 隆起構造および翼

Country Status (3)

Country Link
US (1) US20220161923A1 (ja)
JP (1) JP7465483B2 (ja)
WO (1) WO2020203284A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7369386B1 (ja) 2022-05-23 2023-10-26 三菱重工業株式会社 隆起構造、翼、隆起構造の設計方法及びその設計プログラム

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116552777B (zh) * 2023-07-05 2023-09-12 中国空气动力研究与发展中心计算空气动力研究所 一种涡流调控器以及一种飞行器

Citations (2)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (2)

* Cited by examiner, † Cited by third party
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

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7369386B1 (ja) 2022-05-23 2023-10-26 三菱重工業株式会社 隆起構造、翼、隆起構造の設計方法及びその設計プログラム

Also Published As

Publication number Publication date
JPWO2020203284A1 (ja) 2020-10-08
US20220161923A1 (en) 2022-05-26
JP7465483B2 (ja) 2024-04-11

Similar Documents

Publication Publication Date Title
US6431498B1 (en) Scalloped wing leading edge
Hansen Effect of leading edge tubercles on airfoil performance.
CN102666275B (zh) 提供有机身和机翼的飞机的方法和由该方法制造的飞机
WO2020203284A1 (ja) 隆起構造および翼
US20060060722A1 (en) Controlling bondary layer fluid flow
US20070194178A1 (en) Passive micro-roughness array for drag modification
US8460779B2 (en) Microstructures for reducing noise of a fluid dynamic structure
US5167387A (en) Porous airfoil and process
EP0532093B1 (en) Perforated wing panel with variable porosity
EP1578663B1 (en) Perforated skin structure for laminar-flow systems
SE464861B (sv) Kropp med reducerat basbromsmotstaand
EP0273851A2 (en) Gas turbine engine airfoil
EP0273850A2 (en) Bodies with reduced base drag
Magill et al. Exploring the feasibility of pulsed jet separation control for aircraft configurations
Dwivedi et al. Effect of surface roughness on aerodynamic performance of the wing with NACA 4412 airfoil at Reynolds number 1.7× 105
Jackson et al. Control of afterbody vortices by blowing
Vasile et al. Interaction of a finite-span synthetic jet near the tip of a sweptback wing
Weihs et al. Comb wings for flapping flight at extremely low Reynolds numbers
Genç et al. Unsteady flow over flexible wings at different low Reynolds numbers
Kaakkunen et al. Fabrication of surfaces with reduced friction using nanosecond laser
de Paula et al. The Thickness Effect on Symmetrical Airfoil Flow Characteristics at low Reynolds number
JP6762604B2 (ja) 翼及び航空機
US20210231142A1 (en) Submerged periodic riblets
JP7369386B1 (ja) 隆起構造、翼、隆起構造の設計方法及びその設計プログラム
JP7015582B2 (ja) 翼及び航空機

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 20782430

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2021511412

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 20782430

Country of ref document: EP

Kind code of ref document: A1