Nothing Special   »   [go: up one dir, main page]

CN117494323B - Design method of high-speed waverider with pressure-matched supersonic cooling air film - Google Patents

Design method of high-speed waverider with pressure-matched supersonic cooling air film Download PDF

Info

Publication number
CN117494323B
CN117494323B CN202410004438.0A CN202410004438A CN117494323B CN 117494323 B CN117494323 B CN 117494323B CN 202410004438 A CN202410004438 A CN 202410004438A CN 117494323 B CN117494323 B CN 117494323B
Authority
CN
China
Prior art keywords
air film
supersonic
nozzle
pressure
supersonic air
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
Application number
CN202410004438.0A
Other languages
Chinese (zh)
Other versions
CN117494323A (en
Inventor
易仕和
张博
陆小革
曾瑞童
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
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 National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN202410004438.0A priority Critical patent/CN117494323B/en
Publication of CN117494323A publication Critical patent/CN117494323A/en
Application granted granted Critical
Publication of CN117494323B publication Critical patent/CN117494323B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • Mathematical Optimization (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computational Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Physics (AREA)
  • Application Of Or Painting With Fluid Materials (AREA)

Abstract

The invention provides a high-speed waverider design method with a pressure-matched supersonic cooling air film, which comprises the following steps: generating a wave multiplier, giving a supersonic air film arrangement position on the wave multiplier, acquiring supersonic air film spray pipe outlet pressure under a pressure matching condition, determining supersonic air film mass flow rate, supersonic air film Mach number, supersonic air film spray pipe throat height and spray pipe outlet height of the supersonic air film spray pipe, obtaining supersonic air film spray pipe transonic solution, setting Mach number distribution of a supersonic air film spray pipe axis, using the supersonic air film spray pipe transonic solution, spray pipe outlet and Mach number distribution of a supersonic air film spray pipe axis as boundary conditions, determining a spray pipe non-sticking molded line, solving boundary layer displacement thickness, and obtaining a final spray pipe molded line, thereby completing the design of the high-speed wave multiplier with a pressure matching supersonic cooling air film. The invention solves the problem that the surface pressure of the high-speed waverider and the pressure of the supersonic air film are matched with each other.

Description

Design method of high-speed waverider with pressure-matched supersonic cooling air film
Technical Field
The invention mainly relates to the technical field of high-speed waverider aircrafts, in particular to a design method of a high-speed waverider with a pressure-matched supersonic cooling air film.
Background
Currently, high-speed wave-taking aircrafts are hot in research in all countries of the world due to the advantages of high flying speed, high speed, difficulty in interception of the 'water-float' flight track and the like. However, high speed flight gives severe aerodynamic heat to the waverider, causing ablation and destruction of the aircraft surface structure, optical window. Taking Mach 6 waverider aircraft as an example, 20km high altitude atmospheric temperature is about 216K, and the turbulent flow Plantain number is 0.9, then the aircraft surface adiabatic temperature is approximately 1720K;100km high air temperature of approximately 195K and laminar flow prandtl 0.7, the aircraft surface insulation temperature is approximately 1340K. The high-speed waverider aircraft is provided with a layer of cooling air film, which is an effective way for improving infrared stealth performance and is the only way for avoiding burning of an optical window.
At present, a supersonic air film of a high-speed wave-taking aircraft is generally arranged in front of an optical window, and a spray pipe is arranged in a wall surface slotting mode, so that the supersonic air film is generated. The supersonic air film for infrared stealth is arranged on the aircraft body, the width of the spray pipe is larger, and the effective cooling area is longer.
The current supersonic velocity air film spray pipe is generally designed based on a one-dimensional spray pipe flow basic theory, and the main design parameters are spray pipe Mach number, air film height, total temperature total pressure and mass flow rate, and the implementation method is as follows:
firstly, determining the supersonic speed air film Mach number and the air film mass flow rate required by the surface cooling length of an aircraft by means of a priori experiment or an empirical curve;
under the condition of total temperature and total pressure of an air source carried by the aircraft, the air film height and throat height of the nozzle outlet are determined along with the air film mass flow rate, and the static temperature and static pressure of the nozzle outlet are determined by an isentropic relation;
setting the length of the spray pipe, setting the Mach number distribution of the axis of the spray pipe, and determining the height of the spray pipe corresponding to the Mach number according to a Mach number and area ratio formula;
and smoothly connecting a series of spray pipe height profile points to obtain a supersonic spray pipe, and arranging the supersonic spray pipe on the surface of the aircraft so as to generate a supersonic air film.
However, current supersonic gas film nozzles suffer from several disadvantages. First of all, the design is theoretically flawed, and the jet flow is not one-dimensional in nature, but at least two-dimensional. In particular, for curved nozzles, the flow is three-dimensional or modified by two-dimensional flow theory in combination with yaw rate. The existing supersonic velocity air film spray pipe is designed by adopting a one-dimensional flow theory, the wall surface of the spray pipe cannot be perfectly wave-cut, in other words, the inner flow passage must have concentrated and converged compression waves, and the problems include: on one hand, the uniformity of the nozzle outlet is very low, the nozzle outlet is mixed with the main flow after being ejected, the heat-resistant performance is reduced, on the other hand, the existence of weak shock waves leads to unequal flow entropy, the static temperature and static pressure of the nozzle outlet are inconsistent with the theoretical results of one-dimensional flow, the air film pressure is difficult to effectively control to be matched with the main flow pressure of an aircraft, the pressure mismatch leads to a wave system structure, the mixing is enhanced, and meanwhile, the optical distortion is increased.
Disclosure of Invention
In the prior art, the supersonic air film cannot be matched with the surface pressure of an aircraft, and the air film and a main stream are mixed and aggravated and flow separated due to the fact that the pressure is not matched, the cooling efficiency is worsened, and the optical distortion is increased. Aiming at the technical problems in the prior art, the invention provides a high-speed waverider design method with a pressure matching supersonic cooling air film.
In order to achieve the above purpose, the technical scheme adopted by the invention is as follows:
on one hand, the invention provides a high-speed waverider design method with a pressure-matched supersonic cooling air film, which comprises the following steps:
generating a waverider, and giving a supersonic air film arrangement position on the waverider;
given the cooling length requirement of the supersonic air film, determining the mass flow rate and Mach number of the supersonic air film;
knowing the surface pressure of the waverider, and determining the outlet pressure of the supersonic velocity air film spray pipe by combining the Mach number of the supersonic velocity air film under the pressure matching condition;
determining the throat height of the supersonic velocity air film spray pipe according to the mass flow rate of the supersonic velocity air film and the outlet pressure of the supersonic velocity air film spray pipe;
determining the nozzle outlet height of the supersonic velocity air film nozzle according to the throat height of the supersonic velocity air film nozzle and the Mach number of the supersonic velocity air film;
solving a parabolic potential function equation for the transonic flow of the throat part of the supersonic gas film spray pipe by adopting a series expansion method to obtain a transonic solution of the supersonic gas film spray pipe;
setting Mach number distribution of the axis of the supersonic air film spray pipe by taking the supersonic air film spray pipe transonic velocity solution and the spray pipe outlet as a starting point and a finishing point respectively;
the supersonic velocity air film spray pipe transonic velocity solution, spray pipe outlet and Mach number distribution of the supersonic velocity air film spray pipe axis are used as boundary conditions, a supersonic velocity air film spray pipe characteristic line grid is constructed, a supersonic velocity flow field in the supersonic velocity air film spray pipe is solved, and then a spray pipe non-sticking line is determined according to a streamline control equation;
solving the displacement thickness of the boundary layer according to the Von Karman momentum integral relation;
performing viscosity correction on the non-sticky molded line of the spray pipe based on the displacement thickness of the boundary layer to obtain a final molded line of the spray pipe;
and connecting a contraction section molded line of the preset supersonic speed air film spray pipe with the spray pipe molded line at the throat to obtain the complete supersonic speed air film spray pipe, and completing the design of the high-speed waverider with the pressure-matched supersonic speed cooling air film.
Compared with the prior art, the invention has the technical effects that at least the following aspects are realized:
(1) The invention solves the problem that the surface pressure of the high-speed wave multiplier and the pressure of the supersonic air film are matched with each other, and the generated supersonic air film is suitable for pressure matching of the high-speed wave multiplier.
(2) The generated supersonic speed air film body is not separated, and the cooling efficiency is high.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and other drawings may be obtained according to the structures shown in these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of a method for designing a high-speed waverider with a pressure-matched supersonic cooling air film according to an embodiment of the present invention;
FIG. 2 is a Mach-number cloud of flow field cell 1 for 2 different incoming flow Mach numbers, wherein (a) is a Mach-number cloud of flow field cell 1 for an incoming flow Mach number of 3 and an object plane angle of 30 °; (b) Mach number cloud pictures of a 1 st flow field unit with an incoming flow Mach number of 5 and an object plane angle of 30 degrees;
FIG. 3 is a Mach-number cloud of a flow field obtained by solving a swirled eigenvector in one embodiment;
FIG. 4 is a view of a wave multiplier profile generated by an embodiment;
FIG. 5 is a transonic de-schematic diagram of a supersonic nozzle in accordance with an embodiment of the invention;
FIG. 6 is a graph of nozzle axis Mach number distribution in one embodiment;
FIG. 7 is a schematic diagram of a nozzle signature line network and a nozzle non-stick profile in one embodiment;
FIG. 8 is a graph of a high speed wavebody surface temperature profile with a pressure matched supersonic cooling film designed for one embodiment.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and fully with reference to the accompanying drawings, in which it is evident that the embodiments described are only some, but not all embodiments of the invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Referring to fig. 1, in one embodiment, a method for designing a high-speed wave multiplier with a pressure-matched supersonic cooling air film is provided, including:
generating a waverider, and giving a supersonic air film arrangement position on the waverider;
given the cooling length requirement of the supersonic air film, determining the mass flow rate and Mach number of the supersonic air film;
knowing the surface pressure of the waverider, and determining the outlet pressure of the supersonic velocity air film spray pipe by combining the Mach number of the supersonic velocity air film under the pressure matching condition;
determining the throat height of the supersonic velocity air film spray pipe according to the mass flow rate of the supersonic velocity air film and the outlet pressure of the supersonic velocity air film spray pipe;
determining the nozzle outlet height of the supersonic velocity air film nozzle according to the throat height of the supersonic velocity air film nozzle and the Mach number of the supersonic velocity air film;
solving a parabolic potential function equation for the transonic flow of the throat part of the supersonic gas film spray pipe by adopting a series expansion method to obtain a transonic solution of the supersonic gas film spray pipe;
setting Mach number distribution of the axis of the supersonic air film spray pipe by taking the supersonic air film spray pipe transonic velocity solution and the spray pipe outlet as a starting point and a finishing point respectively;
the supersonic velocity air film spray pipe transonic velocity solution, spray pipe outlet and Mach number distribution of the supersonic velocity air film spray pipe axis are used as boundary conditions, a supersonic velocity air film spray pipe characteristic line grid is constructed, a supersonic velocity flow field in the supersonic velocity air film spray pipe is solved, and then a spray pipe non-sticking line is determined according to a streamline control equation;
solving the displacement thickness of the boundary layer according to the Von Karman momentum integral relation;
performing viscosity correction on the non-sticky molded line of the spray pipe based on the displacement thickness of the boundary layer to obtain a final molded line of the spray pipe;
and connecting a contraction section molded line of the preset supersonic speed air film spray pipe with the spray pipe molded line at the throat to obtain the complete supersonic speed air film spray pipe, and completing the design of the high-speed waverider with the pressure-matched supersonic speed cooling air film.
The type of the waverider is not limited, the design method is not limited, a person skilled in the art can select a proper method to generate the waverider based on the prior art, the method comprises the following steps of, but not limited to, a oscillometric axisymmetric von Karman waverider design method disclosed in Chinese patent application with publication number CN109573092A, a oscillometric axisymmetric von Karman waverider design method fused with a low-speed airfoil disclosed in Chinese patent application with publication number CN109573093A, a variable wall pressure distribution law oscillometric flow field waverider design method disclosed in Chinese patent with publication number CN109598062B, and the like.
Without loss of generality, an embodiment provides a method for designing a waverider, including:
determining the cruising height and cruising Mach number of the high-speed waverider aircraft, and selecting a reference flow field bus, such as a von Karl curve with the minimum resistance characteristic, as an object plane;
solving a 1 st flow field unit by adopting a Taylor-Michael (Taylor-Maccoll) method;
solving a reference flow field by using a 1 st flow field unit and an object plane as boundary conditions through a swirled characteristic line method;
and according to the given waverider trailing edge line, determining the upper and lower surfaces of the waverider through streamline tracking, and generating the waverider.
Wherein the control equation for the 1 st flow field unit is as follows:
in the middle ofIs the velocity of the vertex ray direction, +.>Is the velocity of the ray normal, +.>Is sound velocity, & lt & gt>Is the angle between the ray emitted from the vertex and the x-axis. No value range exists. The right end of the formula is a known quantity, the left end +.>Also known, the unknown isAnd 2 equations solve 2 unknowns, and solve the ordinary differential equation set through a Dragon-Greek tower method to obtain the 1 st flow field unit. Figure 2 illustrates a 1 st flow field cell mach number cloud of 2 different incoming flow mach numbers in an embodiment wherein (a) is a 1 st flow field cell mach number cloud of 3 for an incoming flow mach number at an object plane angle of 30 °; (b) Mach number cloud for the 1 st flow field unit with incoming flow Mach number of 5 and object plane angle of 30 degrees.
The reference flow field is solved by a method of a rotation characteristic line, and the control equation is as follows:
where x represents the coordinates of the x-axis, y represents the coordinates of the y-axis,is the flow direction angle +.>Is the mach angle, subscript + denotes the left line feature line, subscript-denotes the right line feature line, and subscript 0 denotes the streamline. />Is the x-axis direction velocity, +.>Is the y-axis direction speed, +.>Is the closing speed->Is Mach number->Density (I)>Is a flow characteristic factor.
And solving the ordinary differential equation set by a second-order Euler method to obtain the reference flow field. As shown in fig. 3, a flow field mach number cloud obtained by solving a swirled characteristic line in one embodiment is shown.
Finally, according to the given waverider trailing edge line, the upper and lower surfaces of the waverider are determined through streamline tracking, and the waverider is generated, as shown in fig. 4, which is a waverider appearance diagram finally obtained by an embodiment.
Under the pressure matching condition, the outlet pressure of the supersonic velocity air film spray pipe is consistent with the surface pressure of the waverider, and according to the relation between the Mach number of the supersonic velocity air film and the isentropic flow, the method comprises the following steps:
obtaining the outlet pressure of the supersonic velocity air film spray pipeWherein->Is the specific heat ratio of the gas +.>Is supersonic speed air film Mach number->The total pressure of the air film, namely the total pressure of the air source provided for the supersonic air film spray pipe in the waverider body is a known quantity.
Throat height of supersonic velocity air film spray pipeDetermined by the following formula:
in the middle ofIs the specific heat ratio of the gas, +.>Is a gas constant->Is still warm and is filled with->Is the outlet pressure of the supersonic velocity air film spray pipe,is the supersonic gas film mach number.
Nozzle outlet height of supersonic velocity air film nozzleDetermined by the following formula:
in the middle ofIs the throat height of the supersonic velocity air film spray pipe, +.>Is the specific heat ratio of the gas, +.>Is the supersonic gas film mach number.
In the design method of the supersonic air film spray pipe provided by the embodiment, the adopted parabolic potential function equation is as follows:
in the middle ofIs critical sound velocity of sound velocity line of throat part of supersonic velocity air film spray pipe,/-for>Is the axial velocity of sonic velocity line at the throat of the supersonic velocity air film spray pipe,/->The radial velocity of the sonic velocity line at the throat of the supersonic velocity air film spray pipe is the radial velocity, and x and y respectively represent the x-axis coordinate and the y-axis coordinate on the sonic velocity line at the throat of the supersonic velocity air film spray pipe.
And solving the parabolic potential function equation by adopting a series expansion method for the transonic flow of the throat part of the supersonic gas film spray pipe to obtain a transonic solution of the supersonic gas film spray pipe, as shown in figure 5. Figure 5 is a throat portion of the supersonic velocity air film nozzle,the throat height is the throat height, the broken line is the sonic velocity line, and the sonic velocity line is obtained by solving the transonic velocity solution of the supersonic velocity air film spray pipe, and the Mach number of the sonic velocity line is 1. Upstream (left) of the sonic line is the subsonic region, with Mach numbers less than 1. Downstream (right) of the sonic line is the supersonic region, with Mach numbers greater than 1.
In the design method of the supersonic air film nozzle provided by the embodiment, the supersonic air film nozzle transonic solution and the nozzle outlet are respectively used as a starting point and a finishing point, mach number distribution of the supersonic air film nozzle axis is set, as shown in fig. 6, according to the supersonic air film nozzle transonic solution boundary condition and the supersonic air film nozzle transonic solution boundary condition, namely the supersonic air film nozzle transonic solution and the nozzle outlet are respectively used as the starting point and the finishing point, and position coordinates, mach number and Mach number 1-order derivatives of the starting point and the finishing point are substituted into a B spline curve expression to obtain Mach number distribution of the supersonic air film nozzle axis.
In the design method of the supersonic air film spray pipe provided by the embodiment, the supersonic air film spray pipe transonic solution, the spray pipe outlet and Mach number distribution of the supersonic air film spray pipe axis are used as boundary conditions, a supersonic air film spray pipe characteristic line grid is constructed based on a characteristic line method, a supersonic flow field in the supersonic air film spray pipe is solved, and further a spray pipe non-sticking line is determined according to a streamline control equation, as shown in fig. 7, in the embodiment, a spray pipe characteristic line network and a spray pipe non-sticking line schematic diagram are provided. The characteristic line method is a conventional technical means in the field, a two-dimensional flow hyperbolic equation is solved by adopting the characteristic line method to design the supersonic air film spray pipe, and the supersonic air film spray pipe with a two-dimensional structure has complete internal wave elimination, so that the shock wave is fundamentally eliminated, and the flow field parameters such as pressure, temperature, speed and the like are controllable in distribution. Two eigenvalue equations are given below:
characteristic line method equation 1:
characteristic line method equation 2:
in the middle ofThickness for loss of momentum>Mach angle->For local Mach number->Is a gasSpecific heat ratio of (c). />Is a flow characteristic factor.
The streamline control equation is:
in the middle ofIs axial speed, +.>Is the radial velocity.
In the design method of the supersonic velocity air film spray pipe provided by the embodiment, the boundary layer displacement thickness is solved through the following formulas
In the middle ofIs the thickness of the loss of momentum, +.>Is the boundary layer displacement thickness, +.>Is wall angle->Is boundary layer shape factor>Is a compressible friction coefficient>For local Mach number->Is the specific heat ratio of the gas.
Further, the boundary layer displacement thickness is obtained according to the methodThen, the non-sticking line of the spray pipe is increased in the normal directionIs used for achieving viscosity correction.
It is understood that, a person skilled in the art may use an existing method or use a direct given method to pre-define the designed contraction section profile of the supersonic gas film nozzle according to experience, and the specific design method of the contraction section profile of the supersonic gas film nozzle is not limited in this application.
In a preferred embodiment, a design method of a contraction section profile of a supersonic gas film spray pipe is provided, wherein the coordinates of each point on the contraction section profile of the supersonic gas film spray pipe are [ ]x,y) Determined from the cubic curve, as follows:
wherein the method comprises the steps ofIs the inlet height of the supersonic velocity air film jet pipe, < + >>Is the throat height of the supersonic velocity air film spray pipe, +.>Is a given shrink section length.
In a preferred embodiment, a design method of a contraction section profile of a supersonic gas film spray pipe is provided, wherein the coordinates of each point on the contraction section profile of the supersonic gas film spray pipe are [ ]x,y) Determined from the victims curve as follows:
wherein the method comprises the steps ofIs the inlet height of the supersonic velocity air film jet pipe, < + >>Is the throat height of the supersonic velocity air film spray pipe, +.>Is a given shrink section length.
In a preferred embodiment, a design method of a contraction section profile of a supersonic gas film spray pipe is provided, wherein the coordinates of each point on the contraction section profile of the supersonic gas film spray pipe are [ ]x,y) Determined from the hyperbola, as follows:
wherein the method comprises the steps ofIs the inlet height of the supersonic velocity air film jet pipe, < + >>Is the throat height of the supersonic velocity air film spray pipe, +.>Is given the length of the constriction->Is the relative position of the connecting point of the hyperbola, taking +.>
The supersonic air film jet pipe designed by the design method of the supersonic air film jet pipe has complete wave elimination, fundamentally eliminates shock waves, has controllable flow field parameter distribution of pressure, temperature, speed and the like, can carry out viscosity correction by depending on a von Karman momentum equation, and is suitable for generating a supersonic air film required by the surface heat reduction stealth of a high-speed wave-taking aircraft.
In order to verify the effectiveness of the above-mentioned design method of the high-speed waverider with the pressure-matching supersonic cooling air film, a simulation example is provided below:
as shown in fig. 8, the high-speed wave multiplier with the pressure-matched supersonic cooling air film designed by the embodiment isolates the high Wen Lailiu of the wave multiplier aircraft and avoids the direct heating of the wave multiplier and the optical window by the high Wen Lailiu. It can be seen that the surface temperature of the wave-taking aircraft downstream of the jet pipe is reduced from 1300K to 400K, and the heat reducing effect is obvious. Also shown in fig. 8 are supersonic film streamlines, and it can be seen that the film uniformly covers downstream along the waverider aircraft wall without flow separation, and thus the effective coverage length is up to the aircraft tail. In this process, although the film-main flow structure is continuously mixed, the film is continuously heated by the main flow, but the cooling effect can still control the temperature below 800K at the tail of the aircraft.
The invention is not a matter of the known technology.
The technical features of the above embodiments may be arbitrarily combined, and all possible combinations of the technical features in the above embodiments are not described for brevity of description, however, as long as there is no contradiction between the combinations of the technical features, they should be considered as the scope of the description.
The above examples merely represent a few embodiments of the present application, which are described in more detail and are not to be construed as limiting the scope of the invention. It should be noted that it would be apparent to those skilled in the art that various modifications and improvements could be made without departing from the spirit of the present application, which would be within the scope of the present application. Accordingly, the scope of protection of the present application is to be determined by the claims appended hereto.
The above description is only of the preferred embodiments of the present invention and is not intended to limit the present invention, but various modifications and variations can be made to the present invention by those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims (10)

1.具有压力匹配超声速冷却气膜的高速乘波体设计方法,其特征在于,包括:1. The design method of high-speed waverider with pressure-matched supersonic cooling air film, which is characterized by: 生成乘波体,在所述乘波体上给定超声速气膜布置位置;Generate a waverider body, and specify the supersonic air film arrangement position on the waverider body; 给定超声速气膜的冷却长度要求,确定超声速气膜质量流率与超声速气膜马赫数;Given the cooling length requirement of the supersonic air film, determine the supersonic air film mass flow rate and supersonic air film Mach number; 已知所述乘波体的表面压力,在压力匹配条件下结合超声速气膜马赫数,确定超声速气膜喷管出口压力,其中在压力匹配条件下,超声速气膜喷管出口压力与乘波体的表面压力一致;The surface pressure of the waverider is known, and under pressure matching conditions combined with the supersonic air film Mach number, the supersonic air film nozzle outlet pressure is determined. Under the pressure matching conditions, the supersonic air film nozzle outlet pressure is the same as the waverider. The surface pressure is consistent; 根据超声速气膜质量流率以及超声速气膜喷管出口压力,确定超声速气膜喷管喉部高度;According to the supersonic air film mass flow rate and the supersonic air film nozzle outlet pressure, determine the throat height of the supersonic air film nozzle; 根据超声速气膜喷管喉部高度与超声速气膜马赫数,确定超声速气膜喷管的喷管出口高度;According to the throat height of the supersonic air film nozzle and the supersonic air film Mach number, determine the nozzle exit height of the supersonic air film nozzle; 对超声速气膜喷管喉部的跨声速流动,采用级数展开方法求解抛物型势函数方程,得到超声速气膜喷管跨声速解;For the transonic flow in the throat of the supersonic air film nozzle, the series expansion method is used to solve the parabolic potential function equation, and the transonic solution of the supersonic air film nozzle is obtained; 以超声速气膜喷管跨声速解与喷管出口分别作为起点与终点,设置超声速气膜喷管轴线的马赫数分布;Taking the transonic solution of the supersonic air film nozzle and the nozzle outlet as the starting point and end point respectively, set the Mach number distribution of the axis of the supersonic air film nozzle; 以超声速气膜喷管跨声速解、喷管出口和超声速气膜喷管轴线的马赫数分布作为边界条件,构建超声速气膜喷管特征线网格,求解超声速气膜喷管内部的超声速流场,进而根据流线控制方程,确定喷管无粘型线;Using the transonic solution of the supersonic film nozzle, the Mach number distribution of the nozzle outlet and the axis of the supersonic film nozzle as boundary conditions, a grid of characteristic lines of the supersonic film nozzle is constructed to solve the supersonic flow field inside the supersonic film nozzle. , and then determine the inviscid profile of the nozzle according to the streamline control equation; 根据冯卡门动量积分关系式,求解边界层位移厚度;According to the von Karman momentum integral relationship, the boundary layer displacement thickness is solved; 基于边界层位移厚度对喷管无粘型线进行粘性修正,得到最终的喷管型线;Based on the displacement thickness of the boundary layer, the inviscid profile of the nozzle is modified with viscosity to obtain the final nozzle profile; 将预设的超声速气膜喷管的收缩段型线与喷管型线在喉部相连,得到完整的超声速气膜喷管,完成具有压力匹配超声速冷却气膜的高速乘波体的设计。Connect the preset contraction section profile of the supersonic air film nozzle to the nozzle profile at the throat to obtain a complete supersonic air film nozzle and complete the design of a high-speed waverider with a pressure-matched supersonic cooling air film. 2.根据权利要求1所述的具有压力匹配超声速冷却气膜的高速乘波体设计方法,其特征在于,还包括,根据超声速气膜马赫数与等熵流关系式,如下:2. The design method of a high-speed waverider with a pressure-matched supersonic cooling air film according to claim 1, further comprising: according to the relationship between the Mach number of the supersonic air film and the isentropic flow, as follows: ; 得到气膜总压,其中p为超声速气膜喷管出口压力,为气体的比热比,/>是超声速气膜马赫数,/>为气膜总压,也就是乘波体内部为超声速气膜喷管提供气源的总压,为待求量。The total air film pressure is obtained, where p is the exit pressure of the supersonic air film nozzle, is the specific heat ratio of the gas,/> is the supersonic air film Mach number,/> is the total pressure of the air film, that is, the total pressure inside the waverider that provides the air source for the supersonic air film nozzle, and is the quantity to be found. 3.根据权利要求1所述的具有压力匹配超声速冷却气膜的高速乘波体设计方法,其特征在于,超声速气膜喷管喉部高度通过下式确定:3. The design method of a high-speed waverider with pressure-matched supersonic cooling air film according to claim 1, characterized in that the height of the throat of the supersonic air film nozzle Determined by the following formula: ; 式中是气体的比热比,/>是气体常数,/>是静温,/>是超声速气膜喷管出口压力,/>是超声速气膜马赫数。in the formula is the specific heat ratio of the gas,/> is the gas constant,/> It’s static temperature,/> is the supersonic air film nozzle outlet pressure,/> is the supersonic air film Mach number. 4.根据权利要求1所述的具有压力匹配超声速冷却气膜的高速乘波体设计方法,其特征在于,超声速气膜喷管的喷管出口高度通过下式确定:4. The design method of a high-speed waverider with pressure-matched supersonic cooling air film according to claim 1, characterized in that the height of the nozzle outlet of the supersonic air film nozzle Determined by the following formula: ; 式中是超声速气膜喷管喉部高度,/>是气体的比热比,/>是超声速气膜马赫数。in the formula is the throat height of the supersonic air film nozzle,/> is the specific heat ratio of the gas,/> is the supersonic air film Mach number. 5.根据权利要求1或权利要求2或权利要求3或权利要求4所述的具有压力匹配超声速冷却气膜的高速乘波体设计方法,其特征在于,所述抛物型势函数方程为:5. The design method of a high-speed waverider with a pressure-matched supersonic cooling air film according to claim 1 or claim 2 or claim 3 or claim 4, characterized in that the parabolic potential function equation is: ; 式中是超声速气膜喷管喉部声速线的临界声速,/>是超声速气膜喷管喉部声速线的轴向速度,/>是超声速气膜喷管喉部声速线的径向速度,x、y分别代表超声速气膜喷管喉部声速线上的x轴坐标、y轴坐标。in the formula is the critical sound speed of the sound speed line at the throat of the supersonic air film nozzle,/> is the axial velocity of the sonic line at the throat of the supersonic air film nozzle,/> is the radial velocity of the sound speed line at the throat of the supersonic air film nozzle, x and y respectively represent the x-axis coordinate and y-axis coordinate on the sound speed line at the throat of the supersonic air film nozzle. 6.根据权利要求5所述的具有压力匹配超声速冷却气膜的高速乘波体设计方法,其特征在于,边界层位移厚度通过联立以下公式求解得到:6. The design method of a high-speed waverider with pressure-matched supersonic cooling air film according to claim 5, characterized in that the boundary layer displacement thickness It is obtained by combining the following formulas: ; ; ; 式中是动量损失厚度,/>是边界层位移厚度,/>为壁面角度,/>为边界层形状因子,为可压缩摩擦系数,/>为当地马赫数,/>为气体的比热比。in the formula is the momentum loss thickness,/> is the boundary layer displacement thickness,/> is the wall angle,/> is the boundary layer shape factor, is the compressible friction coefficient,/> is the local Mach number,/> is the specific heat ratio of the gas. 7.根据权利要求5所述的具有压力匹配超声速冷却气膜的高速乘波体设计方法,其特征在于,对喷管无粘型线进行粘性修正的方法是:喷管无粘型线在法向增加的距离实现粘性修正。7. The design method of a high-speed waverider with pressure-matched supersonic cooling air film according to claim 5, characterized in that the method for viscosity correction of the nozzle inviscid profile is: nozzle inviscid profile method to increase distance to achieve viscous correction. 8.根据权利要求1或权利要求2或权利要求3或权利要求4或权利要求6或权利要求7所述的具有压力匹配超声速冷却气膜的高速乘波体设计方法,其特征在于,超声速气膜喷管的收缩段型线上的各点坐标(x,y)由维托辛斯基曲线确定,如下:8. The design method of a high-speed waverider with a pressure-matched supersonic cooling air film according to claim 1 or claim 2 or claim 3 or claim 4 or claim 6 or claim 7, characterized in that the supersonic air The coordinates ( x , y ) of each point on the shrinkage section line of the membrane nozzle are determined by the Witosinski curve, as follows: ; 其中是超声速气膜喷管的入口高度,/>是超声速气膜喷管喉部高度,/>是给定的收缩段长度。in is the entrance height of the supersonic air film nozzle,/> is the throat height of the supersonic air film nozzle,/> is the given contraction segment length. 9.根据权利要求1或权利要求2或权利要求3或权利要求4或权利要求6或权利要求7所述的具有压力匹配超声速冷却气膜的高速乘波体设计方法,其特征在于,超声速气膜喷管的收缩段型线上的各点坐标(x,y)由五次方曲线确定,如下:9. The design method of a high-speed waverider with a pressure-matched supersonic cooling air film according to claim 1 or claim 2 or claim 3 or claim 4 or claim 6 or claim 7, characterized in that the supersonic air The coordinates ( x , y ) of each point on the shrinkage section line of the membrane nozzle are determined by the fifth power curve, as follows: ; 其中是超声速气膜喷管的入口高度,/>是超声速气膜喷管喉部高度,/>是给定的收缩段长度。in is the entrance height of the supersonic air film nozzle,/> is the throat height of the supersonic air film nozzle,/> is the given contraction segment length. 10.根据权利要求1或权利要求2或权利要求3或权利要求4或权利要求6或权利要求7所述的具有压力匹配超声速冷却气膜的高速乘波体设计方法,其特征在于,超声速气膜喷管的收缩段型线上的各点坐标由双三次曲线确定,如下:10. The design method of a high-speed waverider with a pressure-matched supersonic cooling air film according to claim 1 or claim 2 or claim 3 or claim 4 or claim 6 or claim 7, characterized in that the supersonic air Coordinates of each point on the shrinkage section line of the membrane nozzle It is determined by the bicubic curve as follows: ; 其中是超声速气膜喷管的入口高度,/>是超声速气膜喷管喉部高度,/>是给定的收缩段长度,/>是双三次曲线衔接点相对位置,取/>in is the entrance height of the supersonic air film nozzle,/> is the throat height of the supersonic air film nozzle,/> is the given shrinkage segment length,/> is the relative position of the connecting point of the bicubic curve, take/> .
CN202410004438.0A 2024-01-03 2024-01-03 Design method of high-speed waverider with pressure-matched supersonic cooling air film Active CN117494323B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202410004438.0A CN117494323B (en) 2024-01-03 2024-01-03 Design method of high-speed waverider with pressure-matched supersonic cooling air film

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202410004438.0A CN117494323B (en) 2024-01-03 2024-01-03 Design method of high-speed waverider with pressure-matched supersonic cooling air film

Publications (2)

Publication Number Publication Date
CN117494323A CN117494323A (en) 2024-02-02
CN117494323B true CN117494323B (en) 2024-03-26

Family

ID=89683428

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202410004438.0A Active CN117494323B (en) 2024-01-03 2024-01-03 Design method of high-speed waverider with pressure-matched supersonic cooling air film

Country Status (1)

Country Link
CN (1) CN117494323B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114608784A (en) * 2022-05-10 2022-06-10 中国空气动力研究与发展中心高速空气动力研究所 Method for obtaining dynamic running pressure matching point of jet flow in jet wind tunnel through ultrasonic velocity jet flow
CN116894303A (en) * 2023-06-06 2023-10-17 中国人民解放军国防科技大学 Design method of supersonic nozzle, supersonic nozzle and high temperature supersonic wind tunnel test platform
WO2023213196A1 (en) * 2022-05-06 2023-11-09 北京航空航天大学 Forward jet drag reduction and heat shielding method for hypersonic pointed-cone aircraft

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20170137116A1 (en) * 2009-07-10 2017-05-18 Peter Ireland Efficiency improvements for flow control body and system shocks
US20180285497A1 (en) * 2017-03-31 2018-10-04 The Government Of The United States Of America, As Represented By The Secretary Of The Navy Numerical Modeling and Performance Analysis of a Scramjet Engine with a Controllable Waverider Inlet Design

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023213196A1 (en) * 2022-05-06 2023-11-09 北京航空航天大学 Forward jet drag reduction and heat shielding method for hypersonic pointed-cone aircraft
CN114608784A (en) * 2022-05-10 2022-06-10 中国空气动力研究与发展中心高速空气动力研究所 Method for obtaining dynamic running pressure matching point of jet flow in jet wind tunnel through ultrasonic velocity jet flow
CN116894303A (en) * 2023-06-06 2023-10-17 中国人民解放军国防科技大学 Design method of supersonic nozzle, supersonic nozzle and high temperature supersonic wind tunnel test platform

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Experimental investigation on supersonic film cooling of hypersonic optical dome under different nozzle pressure ratios;Xiaobin Sun等;Aerospace Science and Technology;20230629;第1-9页 *
Hypervelocity imperfect gas nozzle design with shared wave-elimination contour;Bo Zhang等;Physics of Fluids;20230808;第1-13页 *
超声速气膜冷却时的光学性能优化设计;易司琪等;应用光学;20170731;第549-554页 *
适用高超声速飞行环境的超声速气膜冷却光学窗口研究进展;易仕和等;空天防御;20211231;第1-13页 *
高超声速飞行器气动设计中的若干关键问题;安复兴;中国科学:物理学 力学 天文学;20210918;第6-25页 *

Also Published As

Publication number Publication date
CN117494323A (en) 2024-02-02

Similar Documents

Publication Publication Date Title
CN105667812B (en) Integrated design method of hypersonic vehicle precursor, inlet and wing waveriders
CN109927917B (en) Integrated design method for internal rotation type wave-rider forebody air inlet channel of supersonic aircraft
CN106005475B (en) Hypersonic inside and outside flow integrated full Waverider aircraft method for designing
CN107672821B (en) Variable Mach number parallel wide-speed-range waverider aircraft design method based on osculating cone theory
CN107514311A (en) Based on rotatable air intake duct/waverider forebody derived integrated design method in precursor shock wave
Chen et al. A novel approach for design and analysis of volume-improved osculating-cone waveriders
CN110001974A (en) Use the propulsion system and supersonic plane of large-scale vortex generator
Otto et al. Inward-turning streamline-traced inlet design method for low-boom, low-drag applications
CN103678774B (en) Designing method for supersonic velocity thrust exhaust nozzle considering inlet parameter unevenness
CN111003196B (en) A full waverider aircraft and its design method and system
CN107336842B (en) A kind of aerodynamic layout method of hypersonic waverider canard
CN108999845B (en) Three-dimensional variable cross-section curved flow channel design method and device based on geometric fusion
Yu et al. 3D inverse method of characteristics for hypersonic bump-inlet integration
CN106741976A (en) A kind of mimetic design method of waverider forebody derived air intake duct integration configuration
Slater Methodology for the design of streamline-traced external-compression supersonic inlets
CN110633522A (en) A Reverse Design Method of Supersonic Thrust Nozzle Based on Maximum Thrust Theory
CN117494323B (en) Design method of high-speed waverider with pressure-matched supersonic cooling air film
Musa et al. New parent flowfield for streamline-traced intakes
Xie et al. Effects of forebody boundary layer on the performance of a submerged inlet
Maru et al. Multi-row disk arrangement concept for spike of axisymmetric air inlet
CN109815564B (en) An inverse design method for supersonic thrust nozzles that can simulate the distribution of aerodynamic parameters at the outlet and determine the shape of the outlet
CN109677630A (en) Waverider design method under the controllable strong geometrical constraint of benchmark flow field profile of shock wave
CN113148102B (en) Mach number variable waverider inverse design method based on local deflection osculating theory
CN103470400B (en) Design method of an air-breathing hypersonic vehicle exhaust nozzle with controllable inlet and outlet shapes
CN119305721A (en) Method for reducing heat and hiding waverider in near space

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant