Embodiment one
Referring to Fig. 1-5, the embodiment of the present invention provides the hypersonic inlet stream of a kind of boundary layer guidance and suction combinationFlowing control method, comprising the following steps:
Step 1, boundary layer guide profile is formed on hypersonic inlet precursor, thus in the guide profile of boundary layerTransverse-pressure gradient is generated, is oriented by the flow direction that transverse-pressure gradient flows through the air-flow of boundary layer guide profileChange;
By the working principle of supersonic speed Bump, the flow direction for flowing through the air-flow of boundary layer guide profile is orientedChange, precursor type face generates transverse-pressure gradient, to carry out active guidance to boundary layer airflow, the direction of guidance can be pre-Fixed, theoretically any direction is ok, and is also possible to being set according to actual design condition or that analysis obtains.Draw in boundary layerConductivity type face is a part of hypersonic aircraft precursor, and its role is to flow through the flowing of boundary layer guide profile overdraughtDirected change occurs for direction, and directed change here refers to flowing towards scheduled direction;And then boundary layer is influenced on boundaryThe distribution of thickness, is adjusted boundary layer thickness under the action of guide profile in boundary layer in layer guide profile;
Step 2, boundary layer is directed to by the transverse-pressure gradient that acts in the boundary layer guide profile describedThe part of boundary layer guide profile;
It can be any part in this local theory, be by being distributed to boundary layer guide profile pressure in the present embodimentIt is well-designed, so that boundary layer is directed separately to the middle part and bottom of boundary layer guide profile under the action of transverse-pressure gradientThen portion is respectively arranged suction hole in the two regions, realizes the suction to boundary layer.It is optimal for being arranged in the two positionsBecause being typically all only to be arranged in bottom, but hypersonic condition boundary layer is difficult bottom, place at middle part andBottom has carried out row respectively and has moved.This is also invention distinguishes in other importances.
In boundary layer, the distribution of boundary layer thickness can change under the action of guide profile.According to boundary layer in boundary layerCharacteristic distributions in guide profile determine puff profile band by given threshold, then on each puff profile band respectivelySuction hole is set.
Step 3, in the local setting aspirating hole of the boundary layer guide profile.In boundary layer, the part of guide profile canThicker boundary layer is formed, in the thicker place setting aspirating hole in boundary layer, the quantity of aspirating hole can be reduced, while improving pumpingAbsorption energy.
Preferably, include: the step of formation boundary layer guide profile in the step 1
Step 11, the pressure on each osculating plane is defined by functional relation using variable θ and x as independent variable;θ be exhibition toVariable, x are to flow to variable, the i.e. axial deflection of air intake duct.The ginseng on a two-dimensional surface is assured that by the two variablesNumber, flow field symmetrical for outer shaft, the osculating plane, which refers to, flows to section with the plane of symmetry angle o degree;Referring to Fig. 1, region 1 is sideInterlayer guide profile;Boundary layer guide profile is a part of hypersonic aircraft precursor, and its role is to flow through boundaryDirected change occurs for the flow direction of layer guide profile overdraught, and the plane where face 2 is the plane of symmetry, and the plane where face 3 isIt is in the osculating plane of θ with the plane of symmetry, point 4 is coordinate origin, while being also boundary layer guide profile starting point, and point 5 is at x=LPoint, while being also boundary layer guide profile terminating point, each osculating plane is defined with variable θ, for Tidal Flow, closelyFace refers to parallel with the plane of symmetry, and distance is that θ flows to section;Axisymmetric flow field is turned for inside/outside, osculating plane refers to is in the plane of symmetryθ angle flows to section.
Step 12, the functional relation of the pressure on each osculating plane and independent variable obtains each and opens up to scale factorPressure distribution curve on osculating plane;
Step 13, the corresponding flow field of pressure distribution curve is obtained based on the method for characteristic curves in each osculating plane;
Step 14, the flow field wall surface line in all osculating planes is subjected to curved surface setting-out, the type face of acquisition is drawn as boundary layerConductivity type face.
Preferably, the step 11 includes:
The intersection point 4 that coordinate origin is boundary layer guide profile start line and the plane of symmetry is defined, following functional relation is passed throughDefine the pressure on each osculating plane:
P (θ, x)=kθy(x)+P0 (1)
Wherein, P0It is given value for the pressure value of coordinate origin;Independent variable x is abscissa, and value range is 0 < x < L,Middle L is boundary layer guide profile terminated line abscissa, that is, puts 5 abscissa;Y (x) is the function of x, meets the following conditions: with xIt is increased monotonically;It is 0 in coordinate origin functional value, i.e. y (0)=0;The symmetrical flow field of the outer shaft of θ indicates osculating plane locating for the flow fieldAngular range with the plane of symmetry is-θc<θ<θc, wherein θcFor given value;kθIt is the function of θ to open up to scale factor;Due to heightSupersonic speed precursor needs to undertake a part of compression duty, therefore boundary layer guide profile flow field edge flows to pressure and constantly increases(boundary layer guide profile plays the role of compression to air-flow, so wall surface constantly increases along stroke pressure).
The step 12 includes:
By changing kθ, controlling each and opening up to scale factor is kθOsculating plane in pressure distribution curve, realize to pressurePower be distributed in exhibition to control, in boundary layer, guide profile forms transverse-pressure gradient, realizes guidance to boundary layer.
Preferably, the step 12 includes:
kθAs θ increases monotone decreasing, by second dervative by kθCurve is divided into following three sections:
θ is located at [0, θ1] section when, kθ" (θ) > 0, and kθ(0)=1;
θ is located at [θ1,θ2] section when, kθ"(θ)<0;
θ is located at [θ2,θc] section when, kθ" (θ) > 0, and kθ(θc)=kmin, wherein kminFor specified value, range is 0 < kmin<1;
kθCurve first derivative at θ=0 is zero, θ=θ1With θ=θ2It is continuous to locate first derivative.
Preferably, k in each sectionθSpecific functional relation includes the letters such as given polynomial type, trigonometric function type, exponential typeThe specific analytic expression of at least one of number relationship.
As a preferred embodiment of the present invention: for example given y (x) is the functional relation of parabolic, a in formulaThe compression of air-flow can be required to give according to practical precursor, value range is generally [0.01,1].
P (θ, x)=kθ(ax2)+P0 (2)
kθIt is the function of θ to open up to scale factor.The present invention is by changing kθ, realize to pressure be distributed in exhibition to controlSystem to form transverse-pressure gradient in boundary layer guide profile, and then realizes the guidance to boundary layer.Below to kθWith θ'sFunctional relation is described.kθAlong the direction θ monotone decreasing, curve can be divided into three sections by second dervative is positive and negative.First segment is[0,θ1] section, there is k in this sectionθ" (θ) > 0, meets kθ(0)=1.Second segment is [θ1,θ2] section, k in this sectionθ"(θ)<0.Third section is [θ2,θc] section, there is k in this sectionθ" (θ) > 0, meets kθ(θc)=kmin, wherein kminFor specified value, rangeFor 0 < kmin<1.In addition, kθCurve first derivative at θ=0 is zero, θ=θ1With θ=θ2It is continuous to locate first derivative, to guaranteeLine smoothing.Concrete functional form can give polynomial type, trigonometric function type, exponential type etc. and specifically parse in each sectionFormula.
One specific implementation case: the given quadratic function distribution form in three sections takes θ1=10, θ2=25, θc=60;θ1Corresponding kθTake 0.8, kmin=0.6.Simultaneously curve meet it is above-mentioned required, equation number etc. under these conditionsIn variable number, the quadratic function expression formula in three sections can be thus solved.(directly solve equation or be converted into lineProperty algebraic solution can obtain totally 9 variables and 9 equations)
By k in each sectionθValue is discrete along the direction θ, and discrete accuracy rating is [1 °, 10 °].It is obtained according to the solution of formula 1It is k that each, which is opened up to scale factor,θOsculating plane in pressure distribution curve.
To osculating plane pressure controllable flow field calculation, given pressure distribution is solved based on the method for characteristic curves in each osculating planeThe corresponding flow field of curve.This method is this field publicly-owned technology, solution procedure can refer to application No. is: 201710784957.3Chinese patent literature and Beijing Higher Education Publishing House -2012 " aerodynamics [M] ".
It is that known two phases are solved using the method estimated-corrected based on there is the rotation method of characteristic curves in specific solution procedureThe downstream wall millet cake these two types unit process that the downstream inner of adjacent internal point selects, known upstream wall millet cake and adjacent inner are selected.SoThe corresponding entire flow field of pressure distribution curve and wall surface molded line are solved according to space step-by-step system afterwards.
The corresponding wall surface of pressure distribution curve in osculating plane is thus obtained.Then it solves according to the method all closeThe corresponding wall surface of pressure distribution curve in section.
The flow field wall surface line in all osculating planes is finally subjected to curved surface setting-out, obtained type face is boundary layer leading typeFace.
Preferably, the step 2 includes:
Step 21, numerical simulation is carried out to boundary layer guide profile, obtains boundary layer guide profile top thickness of boundary layer pointCloth.In k of the present inventionθSetting boundary layer close to θ1<θ<θ2And θ=θcRegion, boundary layer thickness extend to formed pairAnswer the first protrusion and the second protrusion.
The step 3 includes:
Step 31, the characteristic distributions according to boundary layer on boundary layer type face determine puff profile band by threshold value;
Step 32, suction hole is respectively set on each puff profile band.
Preferably, the step 31 includes:
Step 311, in the case that boundary layer leading type is not added in boundary layer, x=x0The boundary layer thickness conduct in sectionT0;
Step 312, for the first elevated regions of boundary layer, give proportional factor r 1, when boundary layer thickness meet T > r1 ×T0When, assert the position on the first puff profile band;
For the second elevated regions of boundary layer, proportional factor r 2 is given, the boundary layer thickness at the position meets T > r2×T0When, assert the position on the second puff profile band;
R1 and r2 is given value, and r1, r2 are all larger than zero, and r1 < r2;
Step 313, puff profile band position is carried out really perpendicular to the section in the direction x to each according to step 312Recognize, obtains the first puff profile band and the second puff profile band in entire boundary layer guide profile.
As a specific embodiment of the invention, (axial direction of air intake duct) arranges one in the guide profile of boundary layer in the x-directionSerial equidistantly to open up to section, section spacing value range is [10,100].
This, which sentences x=x, is confirmed to puff profile band position into section in each open up0It is illustrated for section.With no addition boundary layer guide profile, identical x=x in the case of the Free Development of boundary layer0The boundary layer thickness in section is as baseQuasi- thickness, is denoted as To。
For the first elevated regions of boundary layer, proportional factor r is given1, when boundary layer thickness meets T > r1·To, it is assumed thatThe position is on the first puff profile band.For the second elevated regions of boundary layer, proportional factor r is given2, side at the positionThickness of boundary layer meets T > r2·To, it is assumed that the position is on the second puff profile band.r1And r2For given value, passed through according to designIt tests, the two value range takes [0.5,2] proper, in general meets r1<r2。
Each is opened up to section according to the method and carries out confirmation of the puff profile with position, available entire sidePuff profile band in interlayer guide profile.
The first puff profile band and the second puff profile band are obtained in the guide profile of boundary layer as a result,.
One specific implementation case: x=x is provided0Section inner boundary layer extend to distribution illustrate puff profileThe confirmation process of band.Boundary layer thickness extends apparent raised to forming two in figure.Numerical simulation obtains T in the sectiono=0.58.Take r1=0.9, take r2=1.04, the puff profile band range for finally obtaining the first elevated regions is [10.2,28.8], is obtainedPuff profile band range to the second elevated regions is [48.0,58.6].
Preferably, the step 32 specifically includes:
Step 321, according to the area of given shape, suction area ratio and single suction hole, the sum of aspirating hole is obtainedAmount;Suction area ratio refers to that the area of the total suction hole percentage with area that accounts for puff profile, total area that aspirates are taken out divided by singleTotal hole count can be obtained in sucker area.
Step 322, the distribution form in hole can be by the way of array equidistant on puff profile band, it is of course possible to rootIt is arranged according to actual demand.Suction hole is laid respectively on each puff profile band according to the distribution form in hole.
The arrangement of suction hole is the publicly-owned technology in this field, and suction hole arrangement can be divided into two: determining suction hole shape,Quantity;The distribution form of given suction hole.
The shape of given bore first, suction hole can be circle, triangle, rectangle etc..Then give suction area ratio andThe area of single suction hole, suction area ratio refers to the area of the total suction hole percentage with area that accounts for puff profile, total to aspirateTotal hole count can be obtained divided by single suction hole area in area.
The distribution form in hole can be by the way of array equidistant on puff profile band, it is of course possible to according to practical needIt asks and is arranged.
The working principle of supersonic speed Bump is introduced hypersonic precursor design by the present invention, and precursor type face is made to generate laterally pressureForce gradient is distributed in boundary layer thickness in the guide profile of boundary layer and forms two to carry out active guidance to boundary layer airflowA protrusion, then in the place arrangement suction hole of protrusion, to realize that a kind of boundary layer suitable for hypersonic inlet is drawnLead the boundary layer flow flowing control method of cooperation suction.
Compared with traditional liposuction technique, under identical suction effect, suction hole of the present invention is less, suction chamberAlso smaller, be conducive to the damaged condition for reducing fuselage material structural strength, while also saving the space of aircraft interior;WithBump air intake duct is compared, and the present invention proposes boundary layer aiming at the problem that boundary layer row moves relative difficulty in hypersonic flowingThe boundary layer controllable way that guidance and suction match, solves the problems, such as this substantially;Cooperate with existing Bump and suctionControl mode is compared, and the present invention guides boundary layer, then targetedly taken out by precursor type face special designingIt inhales, the effect of suction hole can be played to a greater extent.By numerical simulation, it was demonstrated that the program is feasible.
The above description is only a preferred embodiment of the present invention, is not intended to limit the scope of the invention, all at thisUnder the inventive concept of invention, using equivalent structure transformation made by description of the invention and accompanying drawing content, or directly/use indirectlyIt is included in other related technical areas in scope of patent protection of the invention.