The present disclosure claims the priority to Chinese Patent Application No. 201711100044.1, titled “METHOD FOR DESIGNING FIXED-PLANFORM WAVERIDER BASED ON OSCULATING CONE THEORY”, filed on Nov. 9, 2017 with the China National Intellectual Property Administration, the content of which is incorporated herein by reference.
The disclosure relates to the field of aerodynamic design of configurations of hypersonic vehicles, and in particular, to configurations of waveriders.
High-lift supersonic or hypersonic configurations have always been an unremitting pursuit of human beings. Aerodynamic performances of a vehicle can be greatly improved according to a hyperbolic characteristic of a hypersonic inviscid flow. A waverider is a typical configuration in utilizing such characteristic. The waverider constrains the high-pressure aerodynamics under a lower surface of the vehicle through attachment of shock waves, so as to prevent flow leakage. A lift-drag barrier of a hypersonic vehicle is effectively broken with a high lift-drag ratio. The waverider has gradually developed in recent decades from a single configuration to complex configurations with different characteristics. In particular, an osculating-cone method is proposed, in which the waverider can be designed under a designated inlet capture curve, and a configuration of the wave rider is provided with various characteristics.
There are still many limitations in engineering application of the waverider. Main problems include poor aerodynamic performance at a low speed, and difficulty in ensuring longitudinal stability. The configuration of the waverider is generally obtained through flow capture based on a hypersonic flow field. A generated curved surface has a unique characteristic, and thereby is difficult to design freely. A planform of the waverider may be modified through a design curve, and it is effective to improve performances of the waverider by controlling the planform of the waverider.
A technical issue addressed by the present disclosure is as follows. A relationship between a planform contour of a configuration and a design curve is established in an osculating-cone method for designing a waverider. Through differential equations, it is capable to designate a planform of the waverider in configuration design of the waverider. Flexibility is improved in designing the waverider, and novel concepts are provided in addressing performance defects such as poor low-speed performances and poor longitudinal stability.
A technical solution of the present disclosure is as follows. A method for designing a fixed-planform waverider based on an osculating cone theory is provided, including following steps:
Further, the equation in the step (1) is established in a following manner:
An method for designing a fixed-planform waverider based on an osculating cone theory is provided, including following steps:
Further, the equation set in the third step is as follows:
Further, a boundary condition for the solving in the fourth step is: values of the three functions are equal at a half yK of a spanwise length, that is, f(y)=c(y)=p(y)|y=y
Further, a process of the solving c(y) in the case that f(y) and p(y) are predetermined in the fourth step includes following steps:
Further, a process of the solving f(y) in the case that c(y) and p(y) are predetermined in the fourth step includes following steps:
Further, the step Δy ranges from yK/2000 to yK/100.
Further, the step Δy is Δy=yK/1000 as optimum.
Further, the PLF may be, but not limited to, a configuration of a delta-wing waverider, a configuration of a double-sweepback waverider, or a configuration of an S-shaped leading edge waverider.
Following beneficial effects are provided according to the present disclosure in comparison with conventional technology.
(1) In conventional design of a waverider, the planform is derived from other design curves and cannot be freely designated. In an osculating-cone method, known design variables are a flow capture start line and an inlet capture curve. It is necessary to apply continuous trial-and-error approximation when designing a configuration with certain planar characteristics, and thereby design flexibility is poor. The present disclosure establishes the relationship between the planform of the waverider and design parameters, gives relational equations, and allows the planform to be freely designated in design. Flexibility is improved in designing the waverider. Since the planform has a great influence on performances of a vehicle, such method for designing the fixed-planform waverider improves performances, such as a low-speed performance and longitudinal stability of the waverider.
(2) In the present disclosure, the relationship between the planform contour of the waverider and the two design parameters is described by a differential equation set, which can be solved by numerical recursion. The boundary condition is that the three curves intersect with each other at a middle of the spanwise length, and the solution is continuously recurred from the middle of the spanwise length toward inside, ensuring of the recursion process.
(3) It is necessary that the step of the recursive solution is reasonably selected. The range provided in the present disclosure can ensure both efficiency of the numerical solution and a reasonable distribution of obtained points on the curve. Smoothness of the solved curve is ensured.
(4) The configurations of the waverider, such as single-sweepback, double-sweepback, and S-shaped leading edge, can be generated based the group of equations and the solution according to embodiments of the present disclosure, providing a basis for improving the low-speed performance and the longitudinal stability.
A design principle of the present disclosure is as follows. A corresponding relationship among the inlet capture curve, the flow capture start line and the planform contour of the waverider is derived, and a manner of numerical solution is determined, according to several elements and assumptions in an osculating-cone method. A design configuration of a planform contour of a waverider can be set, that is, a planform can be determined, according to such relationship. A fixed-planform waverider that is reasonably designed, for example, waveriders with a double-sweepback configuration, S-shaped leading edge configuration, or the like, is advantageous in performances such as a low-speed performance and longitudinal stability.
First, a design principle of the osculating cone waverider is briefly introduced. As shown in
An osculating-cone waverider can be treated as a combination of configurations of sub-waveriders within a series of osculating plane. Herein an arbitrary one of the osculating planes is taken as an example to derive a geometric relationship. Reference is made to
A length and a width of the sub-waverider corresponding to the osculating plane FG are Llocal and Wlocal, respectively. According to the definition of the leading-edge sweepback angle, there is a following equation.
tan λ=Llocal/Wlocal.
In each osculating plane, a quasi-two-dimensional conical flow of a corresponding scale is selected according to a local radius of curvature. In a case that the radius of curvature is infinite, a two-dimensional wedge flow is selected. In the osculating-cone method, a shock wave angle β of a flow in each osculating plane is same, and a shock wave in each conical flow and each wedge flow is along a straight line. Therefore, there is a following equation.
local tan(β).
It is noted that signs of δ1 and δ2 are same as sings of slopes of local tangent of the ICC and the FCT. Therefore, ∠FHG=δ1−δ2. There is a following geometric relationship.
Another equation is obtained based on the above three equations.
A configuration of a waverider with a fixed-sweepback can be generated based on the equation (1). Generally, the leading edge of the waverider is designated as a straight line with a fixed tangent angle λ. One of the ICC or the FCT is given, that is, δ1 or δ2 serves as a basis, to solve distribution of δ2 or δ1, respectively. A configuration of the waverider is generated through a conventional osculating-cone method.
According to the definitions of f(y), c(y) and p(y), there are three equations as follows.
tan(δ1)=c(1)(yG)
tan(δ2)=f(1)(yF)
tan(λ)=p(1)(yF) (2)
According to the definition of the osculating plane, there is an equation as follows.
yF and yG are y-coordinates of point F and point G, respectively. The superscript ‘(1)’ represents calculating a first-order derivative.
The equation set formed by equations (1), (2), and (3) is a geometric relationship according to an embodiment of the present disclosure, and can be solved through a numerical method. A boundary conditions need to be set in solving. In an embodiment of the present disclosure, the boundary condition is an intersection point K of the three curves, that is, f(y)=c(y)=p(y)|y=y
As an alternative of the boundary condition being an intersection at point K, the boundary condition may be set at a symmetry axis of the vehicle.
It should be noted that although the equation set (4) is derived from a single osculating plane, the relationship fits within the whole spanwise length of the waverider. In the equation set (4), yG, yF, δ1, δ2 and λ are unknown variables, and β is a known variable as the shock wave angle of a conical flow. A quantity of equations is 5. In a case that any two of the functions f(y),c(y),p(y) is known, a quantity of the unknowns is 6, including five unknown variables and an unknown equation. Hence, a third of the functions f(y),c(y),p(y) can be solved, according to an ordinary differential equation theory.
According to the equation set (4), any of the functions f(y),c(y) and p(y) can be solve in a case that the other two are predetermined. A configuration of the waverider can be designed based on a planform, that is, the contour p(y), determined through the equation set (4). Specific steps of implementation are as follows. The curve p(y) is given, which is generally a quadratic differentiable curve. One of c(y) and f(y) is given to solve another curve. The c(y) or f(y) obtained from the above method serves as an input of the method for designing the osculating-cone waverider, and a planform contour of a waverider configuration generated by the method is the p(y). There may be two cases as follows.
{circle around (1)} c(y) is solved based on f(y) and p(y). In a case that an upper surface of the waverider is generated from a free-flow surface, the f(y) curve is a contour curve of an upper surface at an inlet, and is a projection of the waverider contour on the y-z plane. p(y) is a projection of the waverider contour on the x-y plane. In such case, the waverider is generated by a three-dimensional contour of the configuration.
{circle around (2)} f(y) is solved based on c(y) and p(y). c(y) represents a contour line of a shock wave at an inlet. In the osculating-cone method, c(y) determines a reference flow field generated by the waverider. In such case, the method is for a design when the planform and the reference flow field of the waverider are predetermined.
In the above two cases, the c(y) or the f(y) can be solved through numerical recursion. In a case that f(y) is solved based on c(y) and p(y) that are predetermined, a process of solving is as follows.
{circle around (1)} Processing is recurred from a boundary at yK toward yF=0, and there are (yG)0=(yF)0=yK and c((yG)0)=f((yF)0) at the boundary;
{circle around (2)} f((yF)i+1) is solved based on a previous processing point ((yG)i,c((yG)i)) in c(y) and (yF)i, where a processing step is Δy, (yF)i+1=(yF)i−Δy. (δ2)i+1, λi+1, (δ1)i+1 are solved based on f(y) and p(y), according to the equation set. A relationship between c(y) and δ1 is discretized to be:
according to a differential rule. ((yG)i+1,c((yG)i+1)) are solved based on the above discretized relationship in combination with
{circle around (3)} Recursion is performed by repeating the step {circle around (2)}, until (yF)i+1=0, such that all coordinate points on c(y) can be obtained. That is, a form of c(y)|y=[0,yK] is obtained.
In a case that c(y) is solved based on f(y) and p(y) that are predetermined, a process of solving is as follows.
{circle around (1)} Processing is recurred from a boundary at yK toward yG=0, and there are (yF)0=(yG)0=yK and f((yF)0)=c((yG)0) at the boundary.
{circle around (2)} (δ1)i, c((yG)i+1) and c(1)((yG)i+1) are solved based on a previous proceeding point ((yF)i,f(yF)i)) in f(y) and (yG)i, where a processing step is Δy, (yG)i+1=(yG)i−Δy. λi is solved based on p(y). (δ2), is solved based on (δ1)i and λi. A relationship between f(y) and δ2 is discretized to be
according to a differential rule. ((yF)i+1,c((yF)i+1)) is solved based on the above discretized relationship in combination with
{circle around (3)} Recursion is performed by repeating the step {circle around (2)} until (yG)i+1=0, such that all coordinate points on f(y) can be obtained. That is, a form of f(y)|y=[0,yK] is obtained.
c(y) and f(y) are obtained in the above two cases, respectively. Then, the waverider configuration can be generated through a conventional osculating-cone method for designing a waverider. In such case, a planform contour of a configuration of the waverider is the designated curve p(y). Thereby, the method allows customizing the planform of the waverider.
Detailed description which is not included herein belongs to common knowledge of those skilled in the art.
Number | Date | Country | Kind |
---|---|---|---|
201711100044.1 | Nov 2017 | CN | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/CN2018/085426 | 5/3/2018 | WO | 00 |