CN108647419A - One kind is with height change low latitude bright eruption infrared signature predictor method and device - Google Patents
One kind is with height change low latitude bright eruption infrared signature predictor method and device Download PDFInfo
- Publication number
- CN108647419A CN108647419A CN201810413157.5A CN201810413157A CN108647419A CN 108647419 A CN108647419 A CN 108647419A CN 201810413157 A CN201810413157 A CN 201810413157A CN 108647419 A CN108647419 A CN 108647419A
- Authority
- CN
- China
- Prior art keywords
- height
- bright eruption
- template
- flow field
- evaluated
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
The present invention relates to technical field of data processing, provide one kind includes with height change low latitude bright eruption infrared signature predictor method and device, this method:The bright eruption characteristic parameter of the flow field under the first height and the second height is obtained based on emulation, and calculates the infrared intensity of bright eruption under the first height and the second height;Bright eruption characteristic parameter of the flow field under the first height and the second height that are obtained according to emulation calculates the bright eruption flow field scale with height change to be evaluated, and obtains the infrared intensity function with height change to be evaluated according to the bright eruption flow field scale of height change to be evaluated and the spectral absorptance of each subregion;It corrects to obtain bright eruption infrared intensity by formula.The low latitude bright eruption infrared signature rapid Estimation with height change can be achieved in the present invention, solve the problems, such as in the past low based on detailed modeling method bright eruption flow field and radiation characteristic computational efficiency.
Description
Technical field
The present invention relates to technical field of data processing, more particularly to one kind is with height change low latitude bright eruption infrared signature
Predictor method and device.
Background technology
Rocket engine design when need to ensure that it be in flowing full state, at this time engine expansion ratio (nozzle exit pressure and
Environmental stress ratio) it is about 1, and when being gone up to the air by ground launch with rocket, as environmental stress reduces in flight course, expansion
Than becoming larger, engine is substantially at slight deficient expansion and deficient swelling state, and another rocket is in gradually acceleration mode.Resume combustion effect is
The important parameter of bright eruption infrared signature is influenced, resume combustion effect is determined by recovery temperature, concentration of component and blending parameter,
With the increase of height, expansion ratio is bigger, undergo mach disk after be lost it is bigger, and then expand after recovery temperature reduce, oxygen
Gas concentration reduces with height, another aspect bright eruption muzzle velocity (about 2000m/s to 3000m/s), as rocket speed increases,
Speed difference is smaller, and blending is weaker, and under this assumed condition, resume combustion effect dies down with the increase of height, is characterized as maximum temperature
Decline.
Therefore, with the variation of height, low latitude bright eruption infrared signature is different.And a certain allusion quotation is obtained based on emulation mode
Flow field and radiation characteristic calculating time cost under type height is higher, less efficient, thus is unsuitable for the low latitude spray of each height
Flame infrared signature is estimated.
Invention content
The technical problem to be solved in the present invention is, in the prior art based on detailed modeling method bright eruption flow field and spoke
The low defect of characteristic computational efficiency is penetrated, provides one kind with height change low latitude bright eruption infrared signature predictor method and dress
It sets, can fast and effeciently realize and estimate.
In order to solve the above technical problem, the present invention provides a kind of pre- with height change low latitude bright eruption infrared signature
Estimate method, including:
The bright eruption characteristic parameter of the flow field under the first height and the second height is obtained based on emulation, and calculates the first height and the
The infrared intensity I of bright eruption under two heighttemplate(H1) and Itemplate(H2);
Bright eruption characteristic parameter of the flow field under the first height and the second height that are obtained according to emulation is calculated with height to be evaluated
The bright eruption flow field scale of variation, and according to the spectral absorption system of bright eruption flow field scale and each subregion with height change to be evaluated
Number obtains the infrared intensity function I changed with height H to be evaluatednew(H);
It corrects to obtain bright eruption infrared intensity I by following formula:
Wherein, Inew(H1) and Inew(H2) it is respectively by the first height H1With the second height H2Substitute into infrared intensity letter
Number Inew(H) infrared intensity obtained after.
Optionally, the bright eruption characteristic parameter of the flow field under obtained first height and the second height according to emulation calculate with
The bright eruption flow field scale of height change to be evaluated, including:
Corresponding expansion ratio N is calculated according to height to be evaluatedPR, select the expansion ratio N under the first height and the second heightPR1With
NPR2In immediate height as basic template, the spray of foundation forms is calculated according to the bright eruption characteristic parameter of the flow field of foundation forms
Flame maximum equivalent radius RtemplateWith bright eruption equivalent length Ltemplate;
According to the bright eruption maximum equivalent radius R of the foundation formstemplateWith bright eruption equivalent length Ltemplate, by with
Lower formula calculates the bright eruption flow field scale changed with height H to be evaluated:
Wherein, RnewFor the bright eruption maximum equivalent radius of height H to be evaluated variations, LnewFor what is changed with height H to be evaluated
Bright eruption equivalent length, NPRAnd NPR_templateFor the expansion ratio for expanding when foundation forms of height to be evaluated, f (U∞template) and
g(U∞template) and for by the adjoint flow velocity degree U of foundation forms∞templateAs with flow velocity degree U∞High speed is substituted into stream shadow
Ring function f (U∞) and g (U∞) obtain, kLAnd kRIt is fitting constant.
Optionally, the high speed is with stream influence function f (U∞) and g (U∞) be respectively:
Wherein, U∞For with flow velocity degree, U*For bright eruption limit speed of expansion, and
γexitFor the specific heat ratio of engine export, R is gas constant, TexitFor engine jet pipe outlet temperature.
Optionally, the basis is with the bright eruption flow field scale of height change to be evaluated and the spectral absorptance of each subregion
Obtain the infrared intensity function I changed with height H to be evaluatednew(H), including:
According to the spectral absorption system of each bright eruption flow field subregion of bright eruption flow field dimension calculation changed with height H to be evaluated
Number;
The infrared intensity changed with height H to be evaluated is obtained according to the spectral absorptance of each bright eruption flow field subregion
Function Inew(H)。
Optionally, the basis is with flow velocity degree U∞The each bright eruption flow field subregion of bright eruption flow field dimension calculation of variation
Spectral absorptance, including:
The spectral absorptance in the area is calculated according to the jet pipe engine export parameter in flow field exits area;
The spectral absorptance in the area is calculated according to the flow field parameter of the incoming zone of influence;
According to the first height H1With the second height H2Flow field parameter calculate height H to be evaluated resume combustion area flow field parameter,
And the spectral absorptance in the area is calculated according to the flow field parameter in resume combustion area.
Optionally, described according to the first height H1With the second height H2Flow field parameter calculate the resume combustion area of height H to be evaluated
Flow field parameter include:
(1) it is calculated by the following formula the temperature T in resume combustion areaafterburning:
Tafterburning=T1_max+(T2_max-T1_max)(H-H1)/(H2-H1);
Wherein T1_maxFor the first height H1Template in temperature maximum, T2_maxFor the second height H2Template in temperature most
Big value;
(2) it is calculated by the following formula the pressure P in resume combustion areaafterburning:
Pafterburning=0.5* (P1(x″,0)+P2(x″′,0));
Wherein P1(x ", 0) is the first height H1Template in temperature maximum corresponding position pressure, P2(x " ', 0) is the
Two height H2Template in temperature maximum corresponding position pressure;
(3) it is calculated by the following formula the density p in resume combustion areaafterburning:
ρafterburning=0.5* (ρ1(x″,0)+ρ2(x″′,0));
Wherein ρ1(x ", 0) is the first height H1Template in temperature maximum corresponding position density, ρ2(x " ', 0) is the
Two height H2Template in temperature maximum corresponding position density;
(4) it is calculated by the following formula the speed U in resume combustion areaafterburning:
Uafterburning=0.5* (U1(x″,0)+U2(x″′,0));
Wherein U1(x ", 0) is the first height H1Template in temperature maximum corresponding position speed, U2(x " ', 0) is the
Two height H2Template in temperature maximum corresponding position speed;
It is calculated by the following formula the CO in resume combustion area2Constituent mass concentration XCO2_afterburning:
XCO2_afterburning=(X1_co2(x″,0)-X1_co2(x ', 0)) * αafterburning/α1+X1_co2(x′,0);
Wherein, X1_co2(x ", 0) is the first height H1Template in CO at the position (x ", 0)2Constituent mass concentration, X1_co2
(x ', 0) is the first height H1Template in CO at the position (x ', 0)2Constituent mass concentration, wherein (x ', 0) is the first height H1's
OH constituent mass concentration X in template1_OHPosition when (x, 0) >=0, αafterburningFor according to the speed U in resume combustion areaafterburning
With the temperature T in resume combustion areaafterburningThe volume of calculating inhales coefficient, α1According to the first height H1Template in temperature maximum correspond to position
The speed U set1(x ", 0) and temperature T1_maxThe volume of calculating inhales coefficient;The CO constituent mass in resume combustion area is calculated using same method
Concentration XCO_afterburningAnd H2O constituent mass concentration XH2O_afterburning。
The present invention also provides one kind with height change low latitude bright eruption infrared signature estimating device, including:Emulation is single
Unit and intensity amending unit are estimated in member, variation;
The simulation unit obtains the bright eruption characteristic parameter of the flow field under the first height and the second height based on emulation, and counts
Calculate the infrared intensity I of bright eruption under the first height and the second heighttemplate(H1) and Itemplate(H2);;
Unit is estimated in the variation, the bright eruption flow field characteristic under the first height and the second height for being obtained according to emulation
Parameter calculate with height change to be evaluated bright eruption flow field scale, and according to height change to be evaluated bright eruption flow field scale and
The spectral absorptance of each subregion obtains the infrared intensity function I changed with height H to be evaluatednew(H);
The intensity amending unit, is used for
It corrects to obtain bright eruption infrared intensity I by following formula:
Wherein, Inew(H1) and Inew(H2) it is respectively by the first height H1With the second height H2Substitute into infrared intensity letter
Number Inew(H) infrared intensity obtained after.
Optionally, the variation estimates unit for performing the following operations to calculate the bright eruption stream with height change to be evaluated
Field scale:
Corresponding expansion ratio N is calculated according to height to be evaluatedPR, select the expansion ratio N under the first height and the second heightPR1With
NPR2In immediate height as basic template, the spray of foundation forms is calculated according to the bright eruption characteristic parameter of the flow field of foundation forms
Flame maximum equivalent radius RtemplateWith bright eruption equivalent length Ltemplate;
According to the bright eruption maximum equivalent radius R of the foundation formstemplateWith bright eruption equivalent length Ltemplate, by with
Lower formula calculates the bright eruption flow field scale changed with height H to be evaluated:
Wherein, RnewFor the bright eruption maximum equivalent radius of height H to be evaluated variations, LnewFor what is changed with height H to be evaluated
Bright eruption equivalent length, NPRAnd NPR_templateFor the expansion ratio for expanding when foundation forms of height to be evaluated, f (U∞template) and
g(U∞template) and for by the adjoint flow velocity degree U of foundation forms∞templateAs with flow velocity degree U∞High speed is substituted into stream shadow
Ring function f (U∞) and g (U∞) obtain, kLAnd kRIt is fitting constant.
Optionally, the high speed is with stream influence function f (U∞) and g (U∞) be respectively:
Wherein, U∞For with flow velocity degree, U*For bright eruption limit speed of expansion, and
γexitFor the specific heat ratio of engine export, R is gas constant, TexitFor engine jet pipe outlet temperature.
Optionally, the variation the estimate unit infra-red radiation to change with height H to be evaluated for performing the following operations
Intensity function Inew(H):
According to the spectral absorption system of each bright eruption flow field subregion of bright eruption flow field dimension calculation changed with height H to be evaluated
Number;
The infrared intensity changed with height H to be evaluated is obtained according to the spectral absorptance of each bright eruption flow field subregion
Function Inew(U∞)。
Implement it is provided in an embodiment of the present invention with height change low latitude bright eruption infrared signature predictor method and device, until
It has the advantages that less:
1, the present invention can based on emulation obtain first height and the second height under bright eruption infrared intensity, and
Based on the infrared intensity function with height change to be evaluated that own transmission principle obtains, different altitude height condition is solved
Fluid field and radiation variation rule, can each height of rapid Estimation bright eruption infrared signature, solve in the past based in detail modeling
Method bright eruption flow field and the low problem of radiation characteristic computational efficiency.
2, the present invention constructs bright eruption flow field scale correction formula, can be counted according to the bright eruption flow field scale of foundation forms
The bright eruption flow field scale changed with height H to be evaluated is calculated, to which the simulation result under level altitude is transformed into Different Altitude
Height H, convenient for being solved to different height above sea level infrared signatures.
3, the present invention is also analyzed by the statistical property to bright eruption flow field, fits high speed with stream influence function f
(U∞) and g (U∞) formula, obtain low latitude bright eruption flow field bright eruption flow field scale with flow velocity degree U∞Changing rule.
4, the present invention is in the infrared intensity I for estimating new height Hnew(HH) after, in order to obtain the spectrum of higher precision
Radiation intensity, the I also further acquired using emulationnew(H1) and Itemplate(H2) be modified, to calculate higher precision
Spectral radiance under different height.
Description of the drawings
Fig. 1 is that the embodiment of the present invention one provides flow with height change low latitude bright eruption infrared signature predictor method
Figure;
Fig. 2 is low latitude bright eruption flow field schematic diagram;
Fig. 3 a and 3b are bright eruption infrared signature matched curve figures in low latitude according to the present invention;
Fig. 4 is the distribution character figure that low latitude bright eruption according to the present invention flow field is estimated;
Fig. 5 is that the embodiment of the present invention five provides signal with height change low latitude bright eruption infrared signature estimating device
Figure;
In figure:501:Simulation unit;502:Unit is estimated in variation;503:Intensity amending unit.
Specific implementation mode
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention
In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is
A part of the embodiment of the present invention, instead of all the embodiments.Based on the embodiments of the present invention, ordinary skill people
The every other embodiment that member is obtained without making creative work, shall fall within the protection scope of the present invention.
Embodiment one
As shown in Figure 1, it is provided in an embodiment of the present invention with height change low latitude bright eruption infrared signature predictor method, it can
To include the following steps:
Step S101:The first height H is obtained based on emulation1With the second height H2Under bright eruption characteristic parameter of the flow field, and calculate
First height H1With the second height H2The infrared intensity I of lower bright eruptiontemplate(H1) and Itemplate(H2);
Step S102:It is calculated according to the bright eruption characteristic parameter of the flow field under emulation obtained the first height and the second height with waiting for
Estimate the bright eruption flow field scale of height change, and according to the light of bright eruption flow field scale and each subregion with height change to be evaluated
Spectral absorption coefficient obtains the infrared intensity function I changed with height H to be evaluatednew(H);H1≤H≤H2;
Step S103:It corrects to obtain bright eruption infrared intensity I by following formula:
Wherein, Inew(H1) and Inew(H2) it is respectively by the first height H1With the second height H2Substitute into infrared intensity letter
Number Inew(H) infrared intensity obtained after.
It is provided in an embodiment of the present invention with height change low latitude bright eruption infrared signature predictor method, emulation can be based on
The infrared intensity of bright eruption under the first height and the second height that obtain, and based on own transmission principle obtain with waiting estimating
The infrared intensity function of calculated altitude variation, solves different altitude height condition Fluid field and radiation variation rule, can be quick
Estimate the bright eruption infrared signature of each height.
Embodiment two
On the basis of embodiment one is provided with height change low latitude bright eruption infrared signature predictor method, step
The first height H is obtained based on emulation in S1011With the second height H2Under bright eruption characteristic parameter of the flow field, and calculate the first height H1
With the second height H2The infrared intensity I of lower bright eruptiontemplate(H1) and Itemplate(H2) process, specifically can be by as follows
Mode is realized:
A1, CFD++ methods are based on using engine jet pipe outlet parameter and environmental parameter as input, emulation obtains two not
With height (the first height H1, the second height H2) under the conditions of bright eruption characteristic parameter of the flow field.First height H in the present invention1, second
Height H2Height above sea level is referred both to height H to be evaluated.Wherein, engine jet pipe outlet parameter includes:Nozzle exit radius
Rexit_template, back pressure Pexit_template, outlet temperature Texit_template, outlet density ρexit_template, muzzle velocity
Uexit_template, outlet specific heat ratio γexit_templateWith constituent mass concentration (such as co2Mass fractionCo mass
Score xcoexit_template、h2O mass fractions).The environmental parameter includes:With stream pressure P∞_template, companion
With stream temperature T∞_templateWith adjoint flow velocity degree U∞_template.The bright eruption characteristic parameter of the flow field includes:Temperature T (x, y), pressure
P (x, y), density p (x, y), constituent mass concentration distribution X (x, y), wherein x and y respectively represent position in flow field, above-mentioned parameter
Characterize the parameter with change in location.Subscript template represents template in the present invention, and exit represents outlet, and it is adjoint that ∞ represents environment
Stream.
Wherein height above sea level is the first height H1When, corresponding environmental stress is P∞1, environment temperature T∞1, speed of incoming flow U∞1,
Expansion ratio N at this timePR1=Pexit/P∞1;When height above sea level is the second height H2When, corresponding environmental stress is P∞2, environment temperature T∞2、
Speed of incoming flow U∞2, expansion ratio N at this timePR2=Pexit/P∞2。
A2, using bright eruption characteristic parameter of the flow field as input, the absorption of gas in bright eruption flow field is solved based on line-by-line integration method
Coefficient solves the radiation transfer equation in bright eruption flow field based on apparent light method (LOS), obtains two different heights (the first height
H1With the second height H2) bright eruption infrared intensity (i.e. spectral radiance).CFD++ and by-line product may be used in the present invention
Divide the template with bright eruption flow field and infrared signature under LOS methods emulation exemplary height, extract bright eruption characteristic parameter of the flow field,
Be conducive to subsequently obtain the infrared intensity function I changed with height H to be evaluatednew(H).The present invention does not limit template
Concrete form, emulates under typical rate, typical two height or flow field that actual measurement obtains and radiation characteristic can be used as template.
Embodiment three
On the basis of embodiment two is provided with height change low latitude bright eruption infrared signature predictor method, step
Bright eruption characteristic parameter of the flow field under the first height and the second height that are obtained according to emulation in S102 is calculated to be become with height to be evaluated
The process of the bright eruption flow field scale of change, can specifically be achieved by the steps of:
B1, corresponding expansion ratio N is calculated according to height to be evaluatedPR, select the expansion ratio N under the first height and the second heightPR1
And NPR2In immediate height as basic template, foundation forms is calculated according to the bright eruption characteristic parameter of the flow field of foundation forms
Bright eruption maximum equivalent radius RtemplateWith bright eruption equivalent length Ltemplate。
Bright eruption flow field is determined by nozzle exit flow parameter and with stream (i.e. ambient air incoming), can be attributed to swollen
It is swollen to compare NPRThe specific heat ratio γ of (nozzle exit pressure and environmental stress ratio), nozzle exitexitWith adjoint flow velocity degree U∞。
When height to be evaluated is H, it is P to correspond to environmental stress at this time∞, environment temperature T∞, speed of incoming flow U∞, at this time this wait for
Estimate the corresponding expansion ratio N of height HPR=Pexit/P∞, since height H to be evaluated is between the first height H1With the second height H2It
Between, i.e. NPRBetween NPR1And NPR2Between, with wherein NPRIt is immediate that be highly basic template, it is assumed that | NPR-NPR1| < |
NPR-NPR2|, then it is the first height H to take foundation forms1, conversely, it is the second height H then to take foundation forms2。
Parameter R is determined by following relationship in templatetemplateAnd Ltemplate, generation method is:The temperature on bright eruption boundary is taken to drop
As low as environment temperature of incoming flow T∞_template1.05 times, as bright eruption maximum equivalent radius Rtemplate, it is from nozzle exit to ring
Oxygen x in the gas of bordero2∞_templateMass fraction is reduced to distance when 0.1, as bright eruption equivalent length Ltemplate。
Characteristic distributions of the present invention according to low latitude bright eruption flow field characteristic, determine bright eruption most from foundation forms by following relationship
Big equivalent redius RtemplateWith bright eruption equivalent length Ltemplate:The temperature T on bright eruption boundary is taken to be reduced to environment temperature of incoming flow
T∞_template1.05 times of radius, as bright eruption maximum equivalent radius Rtemplate;Oxygen in from nozzle exit to environmental gas
xo2∞_templateMass fraction YaDistance when being reduced to 0.1, as bright eruption equivalent length Ltemplate.As shown in Fig. 2, being low latitude
Bright eruption flow field schematic diagram.Therefore, parameter RtemplateAnd LtemplateIt can as stated above be handled and be obtained by template.
B2, the bright eruption maximum equivalent radius R according to the foundation formstemplateWith bright eruption equivalent length Ltemplate, pass through
Following formula calculates the bright eruption flow field scale changed with height H to be evaluated:
Wherein, RnewFor the bright eruption maximum equivalent radius of height H to be evaluated variations, LnewFor what is changed with height H to be evaluated
Bright eruption equivalent length, NPRAnd NPR_templateFor the expansion ratio for expanding when foundation forms of height to be evaluated, f (U∞template) and
g(U∞template) and for by the adjoint flow velocity degree U of foundation forms∞templateAs with flow velocity degree U∞High speed is substituted into stream shadow
Ring function f (U∞) and g (U∞) obtain, kLAnd kRIt is fitting constant.For liquid-propellant rocket engine, it is preferable that kL=1.54, kR
=1.2;For solid propellant rocket, it is preferable that kL=1.32, kR=1.1.
Present invention data through a large number of experiments, analyze the statistical property in bright eruption flow field, find for given
Specific heat ratio γexitAnd U∞< U*/ 2, with flow velocity degree U∞Increase, corresponding bright eruption equivalent length increases therewith, works as companion
Bright eruption limit speed of expansion U is leveled off to flow velocity degree*When, bright eruption equivalent length reaches its maximum value, after with stream increase,
Bright eruption equivalent length reduces.As best shown in figures 3 a and 3b, it is according to low latitude bright eruption infrared signature matched curve figure of the invention.
The curve is the L counted in the corresponding flow field of 100 groups of difference speed of incoming flow based on numerical simulationtemplateIt is fitted based on data
It obtains, wherein γexit=1.2.
Fig. 3 a are bright eruption equivalent length LtemplateWith expansion ratio NPRMatched curve, Fig. 3 b be bright eruption equivalent length
LtemplateWith with flow velocity degree U∞Matched curve.
According to above-mentioned matched curve, high speed can be obtained with stream influence function f (U∞) and g (U∞) be respectively:
Wherein, U∞For with flow velocity degree, U*For bright eruption limit speed of expansion, and
γexitFor the specific heat ratio of engine export, R is gas constant, TexitFor engine jet pipe outlet temperature.
Obtaining the bright eruption maximum equivalent radius R for calculating foundation formstemplateWith bright eruption equivalent length Ltemplate, and
High speed is with stream influence function f (U∞) and g (U∞) after, so that it may to acquire height H variations to be evaluated by formula (2) and (3)
Bright eruption flow field scale RnewAnd Lnew。
The present invention constructs bright eruption flow field scale correction formula, can be calculated according to the bright eruption flow field scale of foundation forms
Go out the bright eruption flow field scale changed with height H to be evaluated, to which the simulation result under level altitude is transformed into Different Altitude height
H is spent, convenient for being solved to different height above sea level infrared signatures.
Further, the present invention is also analyzed by the statistical property to bright eruption flow field, fits high speed with stream shadow
Ring function f (U∞) and g (U∞) formula, obtain low latitude bright eruption flow field bright eruption flow field scale with flow velocity degree U∞Variation rule
Rule.
Example IV
On the basis of embodiment three is provided with height change low latitude bright eruption infrared signature predictor method, step
Basis is obtained with the bright eruption flow field scale of height change to be evaluated and the spectral absorptance of each subregion with to be evaluated in S102
The infrared intensity function I of height H variationsnew(H) process can be specifically achieved by the steps of:
C1, according to height H to be evaluated variation each bright eruption flow field subregion of bright eruption flow field dimension calculation spectral absorption
Coefficient km,η, wherein m is partition number, and η is wavelength.The bright eruption flow field subregion of the present invention includes resume combustion area, the incoming zone of influence and stream
Field outlet area, corresponding spectral absorptance is respectively k1,η、k2,ηAnd k3,η。
C2, according to the spectral absorptance of each bright eruption flow field subregion obtain it is strong with the infra-red radiation of height change to be evaluated
Spend function Inew(H)。
Referring to Fig. 4, for the distribution character figure estimated according to the low latitude bright eruption flow field of the present invention.As shown in figure 4, of the invention
Bright eruption flow field be divided into 3 bright eruption flow field subregions, the region of each subregion and corresponding characterisitic parameter are respectively:
1, flow field exits area:0≤l≤k1*Lnew, 0≤R≤Rexit_template;
The parameter in flow field exits area has outlet temperature Texit_template, back pressure Pexit_template, outlet density
ρexit_template、co2Mass fractionh2O mass fractionsCo mass fractions xcoexit_templateFor
Input parameter in parameter, that is, foundation forms of jet pipe engine export.Above-mentioned parameter is abbreviated as T in Fig. 4e、Pe、ρe、XCO、XH2O、
XCO2。
2, the incoming zone of influence:0≤l≤k1*Lnew, Rexit_template≤R≤Rnew;
The parameter of the incoming zone of influence has an impact the temperature in area, pressure, density, co mass fractions, h2O mass fractions and co2Matter
Measure score, i.e. Ti、Pi、ρi、XCO_i、XH2O_i、XCO2_i.These parameters are running parameter of the nozzle exit to incoming environment, with radial direction
Distance R is changed linearly,
Such as:When radial distance R is equal to Rexit_template, T at this timei=Texit_template, when radial distance R is equal to Rnew,
T at this timei=1.05T∞_template;
When distance R is between Rexit_templateAnd RnewBetween when:
Ti=Texit_template+(R-Rexit_template)*(1.05T∞_template-Texit_template)/(Rnew-
Rexit_template) (6)
Method equally solves remaining parameter according to this.
3, resume combustion area:k1*Lnew≤l≤Lnew, 0≤R≤Rnew。
The temperature in resume combustion area, pressure, density, co mass fractions, h2O mass fractions and co2Mass fraction Tafterburning、
Pafterburning、ρafterburning、XCO_afterburning、XH2O_afterburning、XCO2_afterburningTo consider the flow field of resume combustion effect
Parameter.
L in above-mentioned subregion division is the length across flow field, i.e., along aircraft flight direction, R is radial distance (i.e. half
Diameter).Rexit_templateFor the nozzle exit radius in known engine jet pipe outlet parameter.k1Computational methods it is as follows:From step
Expansion ratio N under the first height and the second height that are selected in rapid S101PR1And NPR2In take X in immediate templateoh(x, y), table
Show that hydroxyl (oh) constituent mass concentration changes with coordinate (x, y), enables y=0, work as Xoh(x,0)>=0 when occurring, and is obtained by template
The x of this position, then k1=x/Rtemplate。
The present invention is according to the first height H1With the second height H2Flow field parameter calculate height H to be evaluated resume combustion area stream
Field parameters include the following steps:
1) it is calculated by the following formula the temperature T in resume combustion areaafterburning:
Tafterburning=T1_max+(T2_max-T1_max)(H-H1)/(H2-H1) (7)
Wherein T1_maxFor the first height H1Template in temperature maximum, T2_maxFor the second height H2Template in temperature most
Big value;
In the first height H1Known bright eruption flow field Temperature Distribution T in template1(x, y) and OH constituent mass concentration distributions X1_OH
(x, y), enables y=0, then has temperature profile function T of the temperature along axis (i.e. the directions L)1(x, 0) and OH constituent mass concentration distributions
Function X1_OH(x, 0), works as X1_OHWhen (x ', 0) >=0, position x ' is acquired, enables x '≤x≤LtemplateIn section, its maximum value is obtained
T1(x ", 0) (is expressed as T1_max), similarly obtain the second height H2Middle temperature maximum of T2(x " ', 0) (is expressed as T2_max).Due to
With the increase of height, expansion ratio (back pressure is than upper environmental stress) is bigger, undergo mach disk after be lost it is bigger, and then expand
Recovery temperature later reduces, and oxygen concentration reduces with height, and (about 2000m/s is arrived another aspect bright eruption muzzle velocity
3000m/s), as rocket speed increases, speed difference is smaller, and blending is weaker, and under this assumed condition, resume combustion effect is with height
Increase die down, be characterized as the decline of maximum temperature.So height above sea level H is between H1And H2Between so that TafterburningBetween
T1_maxAnd T2_maxBetween, to highly making linear interpolation, obtain the T of above-mentioned formula (7) calculatingafterburning:
2) it is calculated by the following formula the pressure P in resume combustion areaafterburning:
Pafterburning=0.5* (P1(x″,0)+P2(x″′,0)) (8)
Wherein P1(x ", 0) is the first height H1Template in temperature maximum corresponding position pressure, P2(x " ', 0) is the
Two height H2Template in temperature maximum corresponding position pressure;
3) it is calculated by the following formula the density p in resume combustion areaafterburning:
ρafterburning=0.5* (ρ1(x″,0)+ρ2(x″′,0)) (9)
Wherein ρ1(x ", 0) is the first height H1Template in temperature maximum corresponding position density, ρ2(x " ', 0) is the
Two height H2Template in temperature maximum corresponding position density;
Since the variation of pressure and density is little, the pressure in resume combustion area can be calculated using above-mentioned (8) and (9)
And density.
(4) it is calculated by the following formula the speed U in resume combustion areaafterburning:
Uafterburning=0.5* (U1(x″,0)+U2(x″′,0)) (10)
Wherein U1(x ", 0) is the first height H1Template in temperature maximum corresponding position speed, U2(x " ', 0) is the
Two height H2Template in temperature maximum corresponding position speed;
It is calculated by the following formula the CO in resume combustion area2Constituent mass concentration XCO2_afterburning:
XCO2_afterburning=(X1_co2(x″,0)-X1_co2(x′,0))*αafterburning/α1+X1_co2(x′,0) (11)
Wherein, X1_co2(x ", 0) is the first height H1Template in CO at position2Constituent mass concentration, X1_co2(x ', 0) is
First height H1Template in CO at the position (x ', 0)2Constituent mass concentration, wherein (x ', 0) is the first height H1Template in OH
Constituent mass concentration X1_OHPosition when (x, 0) >=0, αafterburningFor according to the speed U in resume combustion areaafterburningWith resume combustion area
Temperature TafterburningThe volume of calculating inhales coefficient, α1According to the first height H1Template in temperature maximum corresponding position speed
U1(x ", 0) and temperature T1_maxThe volume of calculating inhales coefficient.
Bright eruption is described with volume suction factor alpha in order to characterize resume combustion effect, in the present invention and rolls up suction gas mixing from ambient air
Empirical coefficient, it is the parameter changed with Mach number.Volume inhales coefficient can be calculated by following formula (12):
M in formulaexitFor outside nozzle Mach number, ρexitThat is outside nozzle density, u are the speed at volume suction position, T
The temperature at position is inhaled for volume.
Therefore, by u=U1(x ", 0), T=T1_maxIt substitutes into formula (12) and acquires volume suction factor alpha1;By u=Uafterburning, T
=TafterburningIt substitutes into formula (12) and acquires volume suction factor alphaafterburning。
The CO constituent mass concentration X in resume combustion areaCO_afterburningAnd H2O constituent mass concentration XH2O_afterburningMethod is similar
With above-mentioned CO2The method for solving of constituent mass concentration is identical, is no longer repeated herein.
Preferably, the basis is with flow velocity degree U∞The each bright eruption flow field subregion of bright eruption flow field dimension calculation of variation
The step of spectral absorptance, including:
(1) spectral absorptance in the area is calculated according to the jet pipe engine export parameter in flow field exits area;
(2) spectral absorptance in the area is calculated according to the flow field parameter of the incoming zone of influence;
(3) according to the first height H1With the second height H2Flow field parameter calculate height H to be evaluated resume combustion area flow field ginseng
It counts, and calculates the spectral absorptance in the area according to the flow field parameter in resume combustion area.
After having divided 3 bright eruption flow field subregions, so that it may to use HITRAN/HITEMP spectra databases, calculate above-mentioned
The absorption coefficient k of each bright eruption flow field subregionm,η.Circular is as follows:
For same gas, the spectral absorptance k at wave number ηηEqual to line of each overlapped spectral line at wave number η
Absorption coefficient kη,iThe sum of, i.e.,:
In formula, kηIt is absorption coefficient, F (η-η0,i) it is spectral line linear function, η0,iFor in computational domain at i-th core
Wave number, SiTo be standardized as the spectral line integrated intensity of individual molecule, N is the number density of molecule of one-component.
Standard state (P is given in HITRAN/HITEMP spectra databases0=1.01325 × 105pa、T0=
Under 296K) in air each component spectral line integrated intensity S*(T0) (being standardized as 296K individual molecules).Individually divide at other temperature
Sub- Radiation Instensity Expression In Atomic Emission Spectrometry:
In formula:η0It is low state spectral term, Q for core position, E "V(T) it is vibrational partition function, QR(T) it is rotation partition
Function, S*(T) when be temperature being T individual molecule spectral line integrated intensity, c is the light velocity, and h is Planck's constant, and k is Boltzmann
Constant.
Wherein spectral line linear function has lorentzian curve (Lorentz), Doppler's line style (Doppler), the present invention to use
Fo Aote (Viogt) line style, calculation formula are as follows:
In formula, WLFor Lorentz spectral line whole-line width, WDFor Doppler's spectral line whole-line width, WVFor Fo Aote line style whole-line widths,
IV,maxFor the value of Fo Aote linear functions at core.
Wherein
Obtain each bright eruption flow field subregion temperature, pressure, density, co mass fractions, h2O mass fractions and co2Quality
After score, so that it may calculate the spectral line integrated intensity of individual molecule to substitute into formula (14) respective temperature as temperature T
S*(T) S i.e. in formula (13)i.By Si、F(η-η0,i) and N substitute into the spectral absorption that single gas is obtained in formula (13)
Coefficient kη.Wherein number density of molecule N is converted to by the density of bright eruption flow field subregion.
For there is the total absorption coefficient of the mixture of n ' kinds of components (i.e. the absorption coefficient of bright eruption flow field subregion)Wherein kη,jThe spectral absorptance for indicating jth kind gas in the subregion is obtained by above-mentioned formula (13).
Preferably, mainly consider that gas component is co in the present invention2、co、h2Tri- kinds of components of o, wherein n '=3.
When due to being observed along different directions, spectral radiance is different, and this field basic technology personnel can basis
The difference of observation angle using the spectral absorptance of each bright eruption flow field subregion come calculate change with height H to be evaluated it is infrared
Radiation intensity function Inew(H).For example, illustrating by for radial direction (the above-mentioned directions R), then it is:
K in formula1,ηFor the spectral absorptance in resume combustion area, k2,ηFor the spectral absorptance of the incoming zone of influence, k3,ηFor flow field
The spectral absorptance of outlet area, k in formula4,η=k2,η。c1For first radiation constant, c2For second radiation constant, λ is radiated wave
The η of long λ=10000/.
According to the method described above, it can calculate and obtain Inew(H), for the first height H1, H=H1When can obtain Inew(H1), it is right
In the second height H2, H=H2When, I can also be calculatednew(H2), in order to obtain the spectral radiance of higher precision, utilize
The I that step S101 is acquiredtemplate(H1) and Itemplate(H2) be modified, the light under different height to calculate higher precision
Spectrum intensity, method are as follows:
By the first height H1Template Itemplate(H1) correct obtain, Inew(H)*Itemplate(H1)/Inew(H1);It is high by second
Spend H2Template Itemplate(H2) correct obtain, Inew(H)*Itemplate(H2)/Inew(H2)。
Thus mean value is taken to can get with bright eruption infrared intensity I in height section:
Wherein I is the spectral radiance that the height that this method estimates acquisition is H, Itemplate(H1) and Itemplate(H2) be
The first height H that the first step is generated based on fine emulation1With the second height H2Bright eruption infrared intensity, Inew(H1) and Inew
(H2) it is respectively by the first height H1With the second height H2Substitute into the infrared intensity function I of estimationnew(H) what is obtained after is infrared
Radiation intensity.
To verify with the efficiency with height change low latitude bright eruption infrared signature predictor method under the influence of stream, the present invention
Emulation platform simulation calculation based on 780 software and hardware configurations of Dell OptipleX from 1km to 40km under interval 1km height altogether
The bright eruption flow field of 40 exemplary heights and infrared signature are emulated using CFD approach (CFD++ softwares) under each height condition
Bright eruption flow field and infrared signature about 11hour are generated, under conditions of not using parallel algorithm, take about 40 in total ×
11hour.Using this predictor method, calculates typicalness (1km and 40km) total time-consuming about 2 × 11hour and prepare template, generate it
He takes 15min by each height bright eruption infrared signature altogether, and time cost greatly reduces, therefore with height change low latitude bright eruption
The quick calculating of low latitude bright eruption flow field and infrared signature may be implemented in infrared signature predictor method.
Embodiment five
As shown in figure 5, it is provided in an embodiment of the present invention with height change low latitude bright eruption infrared signature estimating device, it can
To include:Unit 502 and intensity amending unit 503 are estimated in simulation unit 501, variation;
Simulation unit 501, for obtaining the bright eruption characteristic parameter of the flow field under the first height and the second height based on emulation, and
Calculate the infrared intensity I of bright eruption under the first height and the second heighttemplate(H1) and Itemplate(H2).The simulation unit
The operation of 501 execution is identical as step S101 in preceding method.
Unit 502 is estimated in variation, for the bright eruption flow field under the first height and the second height for being obtained according to emulation
Characterisitic parameter calculates the bright eruption flow field scale with height change to be evaluated, and according to the bright eruption flow field ruler with height change to be evaluated
The spectral absorptance of degree and each subregion obtains the infrared intensity function I changed with height H to be evaluatednew(H).The change
The operation that the execution of unit 502 is estimated in change is identical as step S102 in preceding method.
Intensity amending unit 503 is corrected to obtain bright eruption infrared intensity I for passing through following formula:
Wherein, Inew(H1) and Inew(H2) it is respectively by the first height H1With the second height H2Substitute into infrared intensity letter
Number Inew(H) infrared intensity obtained after.Step S103 in the operation of the intensity amending unit 503 execution and preceding method
It is identical.
Optionally, variation estimates unit 502 for performing the following operations to calculate the bright eruption stream with height change to be evaluated
Field scale:
Corresponding expansion ratio N is calculated according to height to be evaluatedPR, select the expansion ratio N under the first height and the second heightPR1With
NPR2In immediate height as basic template, the spray of foundation forms is calculated according to the bright eruption characteristic parameter of the flow field of foundation forms
Flame maximum equivalent radius RtemplateWith bright eruption equivalent length Ltemplate;
According to the bright eruption maximum equivalent radius R of the foundation formstemplateWith bright eruption equivalent length Ltemplate, by with
Lower formula calculates the bright eruption flow field scale changed with height H to be evaluated:
Wherein, RnewFor the bright eruption maximum equivalent radius of height H to be evaluated variations, LnewFor what is changed with height H to be evaluated
Bright eruption equivalent length, NPRAnd NPR_templateFor the expansion ratio for expanding when foundation forms of height to be evaluated, f (U∞template) and
g(U∞template) and for by the adjoint flow velocity degree U of foundation forms∞templateAs with flow velocity degree U∞High speed is substituted into stream shadow
It rings function f (U ∞) and g (U ∞) is obtained, kL and kRIt is fitting constant.
Optionally, the high speed is with stream influence function f (U∞) and g (U∞) be respectively:
Wherein, U∞For with flow velocity degree, U* is bright eruption limit speed of expansion, andγ
Exit is the specific heat ratio of engine export, and R is gas constant, TexitFor engine jet pipe outlet temperature.
Optionally, the variation the estimate unit 502 infrared spoke to change with height H to be evaluated for performing the following operations
Penetrate intensity function Inew(H):
According to the spectral absorption system of each bright eruption flow field subregion of bright eruption flow field dimension calculation changed with height H to be evaluated
Number;
The infrared intensity changed with height H to be evaluated is obtained according to the spectral absorptance of each bright eruption flow field subregion
Function Inew(U∞)。
It should be noted that above-mentioned each embodiment is provided estimates dress with height change low latitude bright eruption infrared signature
The contents such as information exchange, the implementation procedure between interior constituent parts are set, due to being based on same design with the method for the present invention embodiment,
Particular content can be found in the narration in the method for the present invention embodiment, and details are not described herein again.
It is further to note that provided in an embodiment of the present invention estimate with height change low latitude bright eruption infrared signature
Device can also be realized by software realization by way of hardware or software and hardware combining.For hardware view,
It is provided in an embodiment of the present invention with height change low latitude bright eruption infrared signature estimating device in addition to processor, memory, network
Except interface and nonvolatile memory, the equipment in embodiment where device usually can also include other hardware, such as negative
The forwarding chip etc. of duty processing message.For implemented in software, as shown in figure 5, as the device on a logical meaning, it is
Corresponding computer program instructions in nonvolatile memory are read to run in memory by the CPU of equipment where it and are formed
's.
In conclusion provided in an embodiment of the present invention with height change low latitude bright eruption infrared signature predictor method and dress
It sets, is based on liquid-propellant rocket engine low latitude bright eruption flow field characteristic, establish a kind of low latitude bright eruption infrared signature with height change
Fast and effectively predictor method, be mainly used in Infrared Targets scene it is real-time/quasi real time generate system.
Finally it should be noted that:The above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations;Although
Present invention has been described in detail with reference to the aforementioned embodiments, it will be understood by those of ordinary skill in the art that:It still may be used
With technical scheme described in the above embodiments is modified or equivalent replacement of some of the technical features;
And these modifications or replacements, various embodiments of the present invention technical solution that it does not separate the essence of the corresponding technical solution spirit and
Range.
Claims (10)
1. one kind is with height change low latitude bright eruption infrared signature predictor method, which is characterized in that including:
The bright eruption characteristic parameter of the flow field under the first height and the second height is obtained based on emulation, and calculates the first height and second high
The infrared intensity I of the lower bright eruption of degreetemplate(H1) and Itemplate(H2);
Bright eruption characteristic parameter of the flow field under the first height and the second height that are obtained according to emulation is calculated with height change to be evaluated
Bright eruption flow field scale, and obtained according to the bright eruption flow field scale of height change to be evaluated and the spectral absorptance of each subregion
To the infrared intensity function I changed with height H to be evaluatednew(H);
It corrects to obtain bright eruption infrared intensity I by following formula:
Wherein, Inew(H1) and Inew(H2) it is respectively by the first height H1With the second height H2Substitute into infrared intensity function Inew
(H) infrared intensity obtained after.
2. according to the method described in claim 1, it is characterized in that, first height obtained according to emulation and the second height
Under bright eruption characteristic parameter of the flow field calculate with height change to be evaluated bright eruption flow field scale, including:
Corresponding expansion ratio N is calculated according to height to be evaluatedPR, select the expansion ratio N under the first height and the second heightPR1And NPR2In
For immediate height as basic template, the bright eruption that foundation forms is calculated according to the bright eruption characteristic parameter of the flow field of foundation forms is maximum
Equivalent redius RtemplateWith bright eruption equivalent length Ltemplate;
According to the bright eruption maximum equivalent radius R of the foundation formstemplateWith bright eruption equivalent length Ltemplate, pass through following formula
Calculate the bright eruption flow field scale changed with height H to be evaluated:
Wherein, RnewFor the bright eruption maximum equivalent radius of height H to be evaluated variations, LnewFor the bright eruption changed with height H to be evaluated
Equivalent length, NPRAnd NPR_templateFor the expansion ratio for expanding when foundation forms of height to be evaluated, f (U∞template) and g
(U∞template) and for by the adjoint flow velocity degree U of foundation forms∞templateAs with flow velocity degree U∞Substitute into high speed influences with stream
Function f (U∞) and g (U∞) obtain, kLAnd kRIt is fitting constant.
3. according to the method described in claim 2, it is characterized in that, the high speed is with stream influence function f (U∞) and g (U∞) point
It is not:
Wherein, U∞For with flow velocity degree, U*For bright eruption limit speed of expansion, andγexitFor
The specific heat ratio of engine export, R are gas constant, TexitFor engine jet pipe outlet temperature.
4. method described in any one of claim 1 to 3, which is characterized in that the basis is with height change to be evaluated
Bright eruption flow field scale and the spectral absorptance of each subregion obtain the infrared intensity function changed with height H to be evaluated
Inew(H), including:
According to the spectral absorptance of each bright eruption flow field subregion of bright eruption flow field dimension calculation changed with height H to be evaluated;
The infrared intensity function changed with height H to be evaluated is obtained according to the spectral absorptance of each bright eruption flow field subregion
Inew(H)。
5. according to the method described in claim 4, it is characterized in that, the basis is with flow velocity degree U∞The bright eruption flow field of variation
The spectral absorptance of each bright eruption flow field subregion of dimension calculation, including:
The spectral absorptance in the area is calculated according to the jet pipe engine export parameter in flow field exits area;
The spectral absorptance in the area is calculated according to the flow field parameter of the incoming zone of influence;
According to the first height H1With the second height H2Flow field parameter calculate height H to be evaluated resume combustion area flow field parameter, and root
The spectral absorptance in the area is calculated according to the flow field parameter in resume combustion area.
6. according to the method described in claim 5, it is characterized in that, described according to the first height H1With the second height H2Flow field
The flow field parameter that parameter calculates the resume combustion area of height H to be evaluated includes:
(1) it is calculated by the following formula the temperature T in resume combustion areaafterburning:
Tafterburning=T1_max+(T2_max-T1_max)(H-H1)/(H2-H1);
Wherein T1_maxFor the first height H1Template in temperature maximum, T2_maxFor the second height H2Template in maximum temperature
Value;
(2) it is calculated by the following formula the pressure P in resume combustion areaafterburning:
Pafterburning=0.5* (P1(x″,0)+P2(x″′,0));
Wherein P1(x ", 0) is the first height H1Template in temperature maximum corresponding position pressure, P2(x " ', 0) is second high
Spend H2Template in temperature maximum corresponding position pressure;
(3) it is calculated by the following formula the density p in resume combustion areaafterburning:
ρafterburning=0.5* (ρ1(x″,0)+ρ2(x″′,0));
Wherein ρ1(x ", 0) is the first height H1Template in temperature maximum corresponding position density, ρ2(x " ', 0) is second high
Spend H2Template in temperature maximum corresponding position density;
(4) it is calculated by the following formula the speed U in resume combustion areaafterburning:
Uafterburning=0.5* (U1(x″,0)+U2(x″′,0));
Wherein U1(x ", 0) is the first height H1Template in temperature maximum corresponding position speed, U2(x " ', 0) is second high
Spend H2Template in temperature maximum corresponding position speed;
It is calculated by the following formula the CO in resume combustion area2Constituent mass concentration XCO2_afterburning:
XCO2_afterburning=(X1_co2(x″,0)-X1_co2(x′,0))*αafterburning/α1+X1_co2(x′,0)
Wherein, X1_co2(x ", 0) is the first height H1Template in CO at the position (x ", 0)2Constituent mass concentration, X1_co2(x′,0)
For the first height H1Template in CO at the position (x ', 0)2Constituent mass concentration, wherein (x ', 0) is the first height H1Template in
OH constituent mass concentration X1_OHPosition when (x, 0) >=0, αafterburningFor according to the speed U in resume combustion areaafterburningAnd resume combustion
The temperature T in areaafterburningThe volume of calculating inhales coefficient, α1According to the first height H1Template in temperature maximum corresponding position speed
Spend U1(x ", 0) and temperature T1_maxThe volume of calculating inhales coefficient;The CO constituent mass concentration in resume combustion area is calculated using same method
XCO_afterburningAnd H2O constituent mass concentration
7. one kind is with height change low latitude bright eruption infrared signature estimating device, which is characterized in that including:Simulation unit, change
Unit and intensity amending unit are estimated in change;
The simulation unit obtains the bright eruption characteristic parameter of the flow field under the first height and the second height based on emulation, and calculates the
The infrared intensity I of bright eruption under one height and the second heighttemplate(H1) and Itemplate(H2);
Unit is estimated in the variation, the bright eruption characteristic parameter of the flow field under the first height and the second height for being obtained according to emulation
Calculate the bright eruption flow field scale with height change to be evaluated, and according to bright eruption flow field scale with height change to be evaluated and each
The spectral absorptance of subregion obtains the infrared intensity function I changed with height H to be evaluatednew(H);
The intensity amending unit is corrected to obtain bright eruption infrared intensity I for passing through following formula:
Wherein, Inew(H1) and Inew(H2) it is respectively by the first height H1With the second height H2Substitute into infrared intensity function Inew
(H) infrared intensity obtained after.
8. device according to claim 7, which is characterized in that the variation estimates unit for performing the following operations in terms of
Calculate the bright eruption flow field scale with height change to be evaluated:
Corresponding expansion ratio N is calculated according to height to be evaluatedPR, select the expansion ratio N under the first height and the second heightPR1And NPR2In
For immediate height as basic template, the bright eruption that foundation forms is calculated according to the bright eruption characteristic parameter of the flow field of foundation forms is maximum
Equivalent redius RtemplateWith bright eruption equivalent length Ltemplate;
According to the bright eruption maximum equivalent radius R of the foundation formstemplateWith bright eruption equivalent length Ltemplate, pass through following formula
Calculate the bright eruption flow field scale changed with height H to be evaluated:
Wherein, RnewFor the bright eruption maximum equivalent radius of height H to be evaluated variations, LnewFor the bright eruption changed with height H to be evaluated
Equivalent length, NPRAnd NPR_templateFor the expansion ratio for expanding when foundation forms of height to be evaluated, f (U∞template) and g
(U∞template) and for by the adjoint flow velocity degree U of foundation forms∞templateAs with flow velocity degree U∞Substitute into high speed influences with stream
Function f (U∞) and g (U∞) obtain, kLAnd kRIt is fitting constant.
9. device according to claim 8, which is characterized in that the high speed is with stream influence function f (U∞) and g (U∞) point
It is not:
Wherein, U∞For with flow velocity degree, U*For bright eruption limit speed of expansion, andγexitFor
The specific heat ratio of engine export, R are gas constant, TexitFor engine jet pipe outlet temperature.
10. the device according to any one of claim 7~9, which is characterized in that the variation estimates unit for executing
The infrared intensity function I to change with height H to be evaluated is operated belownew(H):
According to the spectral absorptance of each bright eruption flow field subregion of bright eruption flow field dimension calculation changed with height H to be evaluated;
The infrared intensity function changed with height H to be evaluated is obtained according to the spectral absorptance of each bright eruption flow field subregion
Inew(U∞)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810413157.5A CN108647419B (en) | 2018-05-03 | 2018-05-03 | Low-altitude flame infrared radiation characteristic estimation method and device along with height change |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810413157.5A CN108647419B (en) | 2018-05-03 | 2018-05-03 | Low-altitude flame infrared radiation characteristic estimation method and device along with height change |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108647419A true CN108647419A (en) | 2018-10-12 |
CN108647419B CN108647419B (en) | 2022-01-28 |
Family
ID=63748613
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810413157.5A Active CN108647419B (en) | 2018-05-03 | 2018-05-03 | Low-altitude flame infrared radiation characteristic estimation method and device along with height change |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108647419B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109948284A (en) * | 2019-03-29 | 2019-06-28 | 北京环境特性研究所 | Tail bright eruption infrared spectral radiant intensity is with viewing directional angle the Fitting Calculation method |
CN113743033A (en) * | 2021-08-30 | 2021-12-03 | 北京航空航天大学 | Prediction method for height of supersonic jet Mach disk |
CN114896684A (en) * | 2022-04-02 | 2022-08-12 | 西安电子科技大学 | Correction function-based method for calculating infrared radiation characteristic of jet flame of high-speed aircraft |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1761588A (en) * | 2003-01-22 | 2006-04-19 | 瓦斯特能量系统有限公司 | Thermodynamic cycles using thermal diluent |
CN101976275A (en) * | 2010-09-21 | 2011-02-16 | 北京航空航天大学 | Airplane infrared radiation and atmospheric transmittance modeling method |
JP2015161444A (en) * | 2014-02-27 | 2015-09-07 | Jfeエンジニアリング株式会社 | waste incinerator |
CN105096323A (en) * | 2015-07-28 | 2015-11-25 | 中国石油天然气股份有限公司 | Pool fire flame height measurement method based on visible light image processing |
CN106570253A (en) * | 2016-10-26 | 2017-04-19 | 中国运载火箭技术研究院 | Real-time space-based infrared visual simulation method |
CN106570225A (en) * | 2016-10-19 | 2017-04-19 | 长春理工大学 | Analog simulation method for infrared jamming bomb |
CN107247823A (en) * | 2017-05-12 | 2017-10-13 | 北京环境特性研究所 | Bright eruption flow field predictor method based on accurate one-dimensional chemical dynamics process and self moulding |
-
2018
- 2018-05-03 CN CN201810413157.5A patent/CN108647419B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1761588A (en) * | 2003-01-22 | 2006-04-19 | 瓦斯特能量系统有限公司 | Thermodynamic cycles using thermal diluent |
CN101976275A (en) * | 2010-09-21 | 2011-02-16 | 北京航空航天大学 | Airplane infrared radiation and atmospheric transmittance modeling method |
JP2015161444A (en) * | 2014-02-27 | 2015-09-07 | Jfeエンジニアリング株式会社 | waste incinerator |
CN105096323A (en) * | 2015-07-28 | 2015-11-25 | 中国石油天然气股份有限公司 | Pool fire flame height measurement method based on visible light image processing |
CN106570225A (en) * | 2016-10-19 | 2017-04-19 | 长春理工大学 | Analog simulation method for infrared jamming bomb |
CN106570253A (en) * | 2016-10-26 | 2017-04-19 | 中国运载火箭技术研究院 | Real-time space-based infrared visual simulation method |
CN107247823A (en) * | 2017-05-12 | 2017-10-13 | 北京环境特性研究所 | Bright eruption flow field predictor method based on accurate one-dimensional chemical dynamics process and self moulding |
Non-Patent Citations (2)
Title |
---|
ZEUTHEN, ED等: ""Radiation Emissions from Turbulent Diffusion Flames Burning Vaporized Jet and Jet-like Fuels"", 《ENERGY & FUELS》 * |
尹雪梅 等: ""火箭尾喷焰辐射远程探测信号随飞行参数的变化规律"", 《郑州轻工业学院学报(自然科学版)》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109948284A (en) * | 2019-03-29 | 2019-06-28 | 北京环境特性研究所 | Tail bright eruption infrared spectral radiant intensity is with viewing directional angle the Fitting Calculation method |
CN109948284B (en) * | 2019-03-29 | 2023-01-20 | 北京环境特性研究所 | Fitting calculation method for infrared spectrum radiation intensity of tail jet flame along with viewing angle |
CN113743033A (en) * | 2021-08-30 | 2021-12-03 | 北京航空航天大学 | Prediction method for height of supersonic jet Mach disk |
CN113743033B (en) * | 2021-08-30 | 2023-12-12 | 北京航空航天大学 | Prediction method for supersonic jet Mach disk height |
CN114896684A (en) * | 2022-04-02 | 2022-08-12 | 西安电子科技大学 | Correction function-based method for calculating infrared radiation characteristic of jet flame of high-speed aircraft |
Also Published As
Publication number | Publication date |
---|---|
CN108647419B (en) | 2022-01-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Gier et al. | Interaction of shroud leakage flow and main flow in a three-stage LP turbine | |
Luo et al. | Tip leakage flow and aeroacoustics analysis of a low-speed axial fan | |
CN108647419A (en) | One kind is with height change low latitude bright eruption infrared signature predictor method and device | |
Klein et al. | A fully coupled approach for the integration of 3D-CFD component simulation in overall engine performance analysis | |
Gao et al. | Unsteady behavior of corner separation in a compressor cascade: Large eddy simulation and experimental study | |
Roumeliotis et al. | Development and integration of rain ingestion effects in engine performance simulations | |
Mehdi | Effect of swirl distortion on gas turbine operability | |
Tang et al. | High altitude low Reynolds number effect on the matching performance of a turbofan engine | |
Dufour et al. | Validation and flow structure analysis in a turbofan stage at windmill | |
Ferrer-Vidal et al. | Generating axial compressor maps to zero speed | |
Piccirillo et al. | Guidelines for the LTO Noise Assessment of Future Civil Supersonic Aircraft in Conceptual Design | |
Klein et al. | Integration of 3D-CFD component simulation into overall engine performance analysis for engine condition monitoring purposes | |
Mo et al. | Numerical investigations of the slot blowing technique on the hypersonic vehicle for drag reduction | |
Thollet et al. | Assessment of body force methodologies for the analysis of intake-fan aerodynamic interactions | |
Tweedt | Computational aerodynamic simulations of an 840 ft/sec tip speed advanced ducted propulsor fan system model for acoustic methods assessment and development | |
Righi et al. | Low-order models for the calculation of compressor sub-idle characteristics | |
Villarreal-Valderrama et al. | Applicability of correlational data-mining to small-scale turbojet performance map generation | |
CN108595867A (en) | One kind is with low latitude bright eruption infrared signature predictor method and device under the influence of stream | |
Lin et al. | Effective boundary conditions and numerical method for flow characteristics of aeroengine compressor at high Mach flight | |
Gurrola Arrieta et al. | Improved local scale generic cycle model for aerothermodynamic simulations of gas turbine engines for propulsion | |
Lee et al. | Part II: Validation and Verification of CFD Analysis for NASA SDT Transonic Fan Stage Coupled With Nacelle | |
Guo et al. | A mixed-fidelity computational model of aero engine for inlet distortion | |
Peters | Assessment of propfan propulsion systems for reduced environmental impact | |
Kächele et al. | Numerical simulations of an intake-compressor system | |
Tu et al. | Experimental and numerical investigation on total pressure distortion generated by a rectangular movable flat baffle |
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 |