Method, System and Equipment for Measuring OH Ultraviolet Spectral Radiation Characteristics in Unbalanced Flow

Abstract

The present invention belongs to the technical field of related calculations of molecular spectral radiation, and discloses a method, system and device for measuring the ultraviolet spectral radiation characteristics of OH in non-equilibrium flow. A simulation software is used to simulate the distributions of the temperature, pressure and component concentration of the gas flow field of a high-speed aircraft flying at high altitude; according to the flow field distribution, the flow field is stratified along the determined line-of-sight direction to obtain the temperature, pressure and concentration distributions of each layer; the chemical reaction rate in the OH collisional-radiative model is calculated; the flow field parameters are brought into the collisional-radiative model, and programming is used to calculate the number density of different vibrational energy levels of OH(A); according to the apparent integration method, programming is used to calculate the spectral radiance of OH in the tail flame or shock layer. The present invention takes into account the influence of the refined energy level structure of the OH molecule on the chemical reaction rate in the collisional-radiative model, thereby greatly improving the accuracy of related calculations and more accurately giving the calculation results of the OH spectral radiation characteristics in non-equilibrium flow.

Description

Technical Field

[0001] The present invention belongs to the technical field of related calculations of molecular spectral radiation, and particularly relates to a method, system and device for measuring the ultraviolet spectral radiation characteristics of OH in non-equilibrium flow. By associating the collision-radiation model with the radiation transfer equation, the spectral radiation characteristics of OH molecules in non-equilibrium flow are obtained. Background Art

[0002] Currently, with the continuous development of the aerospace industry, the requirements for the accuracy and rapidity of detecting high-speed aircraft are also constantly increasing. Studying the radiation characteristics of high-speed aircraft bodies has extremely important strategic and military significance. During the flight of a high-speed aircraft, a bow shock layer will form at the warhead part of the aircraft; in addition, the propellant burns in the nozzle and undergoes chemical reactions to produce a tail flame. There are vibrational-rotational transitions of molecules and atoms in the shock layer and the tail flame, making the flow field in a chemically non-equilibrium state. Therefore, obtaining the energy level distribution of molecular shock states in the shock layer or tail flame flow field of a high-speed aircraft has become the primary problem in studying the ultraviolet radiation of the aircraft body.

[0003] In the process of calculating the energy level distribution of molecular excited states in a non-equilibrium flow field based on the collision-radiation model, the refined energy level structure of the molecule itself is the main reason affecting the chemical reaction rate in the model. Therefore, studying the dependence of the chemical reaction rate on the refined energy level structure of the molecule is an important part of studying the energy level distribution of molecular excited states.

[0004] Currently, the collision-radiation model is mainly used to solve the energy level distribution of excited states. However, usually, the transition process of vibrational energy levels in substitution reactions in the model is simplified, and it is assumed that the predissociation reaction rate depends on the vibrational energy level. Such simplification and assumption usually have large errors, and the analysis of problems is relatively rough, and it cannot truly describe the energy level distribution of molecular excited states.

[0005] Through the above analysis, the problems and defects existing in the prior art are: the simplification and assumption in the prior art usually have large errors, the analysis of problems is relatively rough, and it cannot truly describe the energy level distribution of molecular excited states.

[0006] The difficulty in solving the above problems and defects is: due to the non-equilibrium characteristics of the flow field around a hypersonic aircraft, the accurate solution of the density of excited state energy levels of molecules and atoms is particularly crucial. The existing classical OH "collision-radiation" model only considers the influence of 2 N 2 and O

[0007] and ignores the coupling effect of substitution reactions and predissociation reactions in the model on the detailed energy level structure of OH molecules. And the accurate solution of the chemical reaction rate is a necessary prerequisite for accurately calculating the density of excited state energy levels. Therefore, it is difficult to accurately solve the reaction rates of various particles participating in various reactions in the flow field.The significance of solving the above problems and deficiencies is as follows: The method for solving the chemical reaction rate involved in the OH collision-radiation model: The refined energy level structure of the OH molecule is coupled with the chemical reaction rate to clarify the refined reaction rate values of each particle participating in the reaction. By expanding the vibrational energy level transition of OH(A) in the substitution reaction and precise dissociation reaction to the rotational energy level, the accurate solution of the chemical reaction rate is achieved, the calculation accuracy of the OH(A) energy level density is improved, thereby significantly enhancing the accuracy of calculating the OH ultraviolet spectral radiance of hypersonic vehicles, providing a theoretical basis for the precise detection of high-speed targets. Summary of the Invention

[0008] Aiming at the problems existing in the prior art, the present invention provides a method, system and device for measuring the OH ultraviolet spectral radiation characteristics in non-equilibrium flow.

[0009] The present invention is implemented as follows. A method for measuring the OH ultraviolet spectral radiation characteristics in non-equilibrium flow, the method for measuring the OH ultraviolet spectral radiation characteristics in non-equilibrium flow includes:

[0010] In the first step, a simulation software is used to simulate the distribution of the gas flow field temperature, pressure and component concentration of a high-speed aircraft flying at high altitude; the shock layer or wake flow field distribution of the hypersonic vehicle is obtained, providing a basis for the flow field stratification in the second step.

[0011] In the second step, according to the flow field distribution, the flow field is stratified along the determined line-of-sight direction to obtain the temperature, pressure and concentration distribution of each component in each layer; the temperature, pressure and concentration distribution of each layer are obtained, providing an accurate input source for the calculations in the fourth and fifth steps.

[0012] In the third step, the chemical reaction rate in the OH collision-radiation model is calculated; an accurate chemical reaction rate is obtained, providing an accurate input source for the fourth step.

[0013] In the fourth step, the flow field parameters are brought into the collision-radiation model, and programming is used to calculate the number density of different vibrational energy levels of OH(A); the number density of each vibrational energy level of OH(A) is obtained, providing an accurate input source for the fifth step, where A represents the first excited state, and the transition of the OH molecule A-X is the main radiation source in the ultraviolet band, simply referred to as excited state OH(A).

[0014] In the fifth step, according to the apparent integration method, programming is used to calculate the spectral radiance of OH in the wake or shock layer. An accurate target body radiation is obtained, providing a theoretical basis for the detection of high-speed targets.

[0015] Furthermore, in the third step, the chemical reaction rate in the OH collision-radiation model is calculated as follows:

[0016] (1) Based on the refined energy level transition theory of OH molecules, calculate the substitution reaction rate values corresponding to the transitions between different vibrational energy levels of the excited state OH(A), and use programming to solve the 144×144 dimensional substitution reaction rate matrix;

[0017] (2) Calculate the predissociation reaction rate values corresponding to different rotational energy levels;

[0018] (3) Calculate the chemical reaction rate values for reactions involving different molecular atoms.

[0019] Further includes:

[0020] 1) The substitution reaction rate strongly depends on the transition process between different vibrational energy levels of OH(A); among them, the transition processes from high energy levels to low energy levels mainly include ν′ = 1 → ν″ = 0, ν′ = 2 → ν″ = 0, ν′ = 2 → ν″ = 1, ν′ = 3 → ν″ = 2; the corresponding transition processes from low energy levels to high energy levels are ν′ = 0 → ν″ = 1, ν′ = 0 → ν″ = 2, ν′ = 1 → ν″ = 2, ν′ = 2 → ν″ = 3; in the process of programming to calculate the chemical reaction rate corresponding to different vibrational energy level transitions, expand the 144×1 dimensional rate matrix in the original model to 144×144 dimensions; use the formula Couple the vibrational energy level and rotational energy level transitions, and note that J″ > J′ should be satisfied in this formula; according to the formula Calculate the reverse rate of the substitution reaction;

[0021] 2) The predissociation reaction rate strongly depends on the rotational energy level. When the vibrational energy level ν′ = 0 and the rotational energy level J′ ≥ 22, when the vibrational energy level ν′ = 1 and the rotational energy level J′ ≥ 11, and for all rotational energy levels with the vibrational energy level ν′ ≥ 2, the predissociation rates have a strong dependence on the rotational quantum number; use the formula Calculate the predissociation reaction rate values corresponding to different rotational energy levels at the same vibrational energy level, where f(J′) = c 1 exp[c 2 J′(J′ + 1)];

[0022] 3) The reactions involving molecular atoms include de-excitation reactions and substitution reactions. The molecular atoms participating in the reactions mainly include O 2 , H 2 , H 2 O, CO, CO 2 , H, O and N 2 ; Obtain the de-excitation reaction rate values depending on the vibration ν′ = 0 - 3, and use the formula Calculate all de-excitation reaction rate values accurate to ν′ = 0 - 3, J′ = 0 - J′ max (ν′); The solution of the substitution reaction rate is calculated according to 1).

[0023] Further, the fifth step calculates the spectral radiance of OH in the tail flame or shock layer using programming based on the apparent integration method, including:

[0024] (1) Calculate the line strength at different temperatures based on the line-by-line method and the HITRAN high-temperature molecular database The absorption coefficient α of a single layer can be calculated using the energy level density obtained in the fourth step, α = S·N·Φ, where Φ is a linear function, and a more accurate absorption coefficient can be obtained using the Voigt line shape;

[0025] (2) Calculate the absorptance of each layer, 1 - e -α·l , where l represents the thickness of a single layer; the transmittance is e -α·l ;

[0026] (3) Select blackbody radiation as the initial radiance I of the target 0 , and calculate the spectral radiance after passing through the entire path based on the apparent integration method , where the optical thickness τ λ = κ sλ ·l; when calculating the radiance of the nth layer, considering the attenuation effect of passing through the paths of the previous n - 1 layers, the discrete form of the equation is

[0027] Another object of the present invention is to provide a computer device, which includes a memory and a processor. The memory stores a computer program. When the computer program is executed by the processor, the processor executes the steps of the method for measuring the OH ultraviolet spectral radiation characteristics in non-equilibrium flow.

[0028] Another object of the present invention is to provide a computer-readable storage medium storing a computer program. When the computer program is executed by a processor, the processor executes the steps of the method for measuring the OH ultraviolet spectral radiation characteristics in non-equilibrium flow.

[0029] Another object of the present invention is to provide an information data processing terminal for implementing the steps of the method for measuring the OH ultraviolet spectral radiation characteristics in non-equilibrium flow.

[0030] Another object of the present invention is to provide a system for measuring the OH ultraviolet spectral radiation characteristics in non-equilibrium flow for implementing the method for measuring the OH ultraviolet spectral radiation characteristics in non-equilibrium flow. The system for measuring the OH ultraviolet spectral radiation characteristics in non-equilibrium flow includes:

[0031] A gas parameter distribution module for simulating the distribution of the temperature, pressure, and component concentration of the gas flow field of a high-speed aircraft flying at high altitude using simulation software;

[0032] A flow field stratification processing module, which is used to stratify the flow field along a determined line-of-sight direction according to the flow field distribution, and obtain the temperature, pressure, and concentration distributions of each component in each layer;

[0033] A chemical reaction rate calculation module, which is used to calculate the chemical reaction rate in the OH collision-radiation model;

[0034] A number density calculation module, which is used to substitute the flow field parameters into the collision-radiation model and programmatically calculate the number density of different vibrational energy levels of OH(A);

[0035] A spectral radiance calculation module, which is used to calculate the spectral radiance of OH in the wake or shock layer by programming according to the apparent integration method.

[0036] Another object of the present invention is to provide a winged aircraft, and the winged aircraft is equipped with the OH ultraviolet spectral radiation characteristic determination system in the non-equilibrium flow.

[0037] Another object of the present invention is to provide a wingless aircraft, and the winged aircraft is equipped with the OH ultraviolet spectral radiation characteristic determination system in the non-equilibrium flow.

[0038] Combining all the above technical solutions, the advantages and positive effects of the present invention are as follows: The present invention associates the collision-radiation model with the radiation transfer equation to obtain the spectral radiation characteristics of OH molecules in the non-equilibrium flow, considers the influence of the refined energy level structure of OH molecules on the chemical reaction rate in the collision-radiation model, mainly expands the vibrational energy level transition of OH(A) in the substitution reaction, and expands the reaction rate matrix to a dimension of 144×144; in addition, the pre-dissociation reaction rate is also refined, and the reaction rate depending on the vibrational energy level in the original model is accurate to depend on the rotational energy level. Thereby greatly improving the accuracy of related calculations and more accurately giving the calculation results of the OH spectral radiation characteristics in the non-equilibrium flow. Description of the Drawings

[0039] Figure 1 is a flowchart of the method for determining the OH ultraviolet spectral radiation characteristics in the non-equilibrium flow provided by the embodiment of the present invention.

[0040] Figure 2 is a schematic structural diagram of the OH ultraviolet spectral radiation characteristic determination system provided by the embodiment of the present invention;

[0041] Figure 2 In the figure: 1. Gas parameter distribution module; 2. Flow field stratification processing module; 3. Chemical reaction rate calculation module; 4. Number density calculation module; 5. Spectral radiance calculation module.

[0042] Figure 3It is the implementation flowchart of the method for measuring the OH ultraviolet spectral radiation characteristics in non-equilibrium flow provided by the embodiments of the present invention.

[0043] Figure 4 It is the schematic diagram of the refined transition of the vibrational energy level in the substitution reaction in the OH collision-radiation model provided by the embodiments of the present invention.

[0044] Figure 5 It is the schematic diagram of the correspondence between the elements in the substitution reaction rate matrix and the vibrational energy level transition provided by the embodiments of the present invention.

[0045] Figure 6 It is the result diagram of the pre-dissociation reaction rate calculation depending on the rotational energy level provided by the embodiments of the present invention.

[0046] Figure 7 It is the verification diagram of the calculation result of the OH(A) energy level density in the afterburner and the Atlas III test provided by the embodiments of the present invention.

[0047] Figure 8 It is the verification diagram of the calculation result of the OH spectral radiance in the afterburner and the Atlas III test provided by the embodiments of the present invention.

[0048] Figure 9 It is the verification diagram of the calculation result of the normalized OH spectral radiance in the shock layer and the BSUV-2 test provided by the embodiments of the present invention. Detailed implementation manners

[0049] In order to make the objectives, technical solutions and advantages of the present invention clearer and more understandable, the present invention will be further described in detail below in conjunction with embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not used to limit the present invention.

[0050] In view of the problems existing in the prior art, the present invention provides a method, system and device for measuring the OH ultraviolet spectral radiation characteristics in non-equilibrium flow. The present invention will be described in detail below in conjunction with the accompanying drawings.

[0051] As Figure 1 shown, the method for measuring the OH ultraviolet spectral radiation characteristics in non-equilibrium flow provided by the present invention includes the following steps:

[0052] S101: Use simulation software to simulate the distribution of the gas flow field temperature, pressure and component concentration of a high-speed aircraft flying at high altitude;

[0053] S102: According to the flow field distribution, perform stratification processing on the flow field along the determined line of sight direction to obtain the temperature, pressure and concentration distribution of each component in each layer;

[0054] S103: Calculate the chemical reaction rate in the OH collision-radiation model;

[0055] S104: Substitute the flow field parameters into the collision-radiation model and programmatically calculate the number densities of different vibrational energy levels of OH(A);

[0056] S105: According to the apparent integration method, programmatically calculate the spectral radiance of OH in the wake or shock layer.

[0057] In step S103: (1) Based on the refined energy level transition theory of OH molecules, calculate the substitution reaction rate values corresponding to the transitions between different vibrational energy levels of OH(A), and programmatically solve a 144×144-dimensional substitution reaction rate matrix. (2) Calculate the predissociation reaction rate values corresponding to different rotational energy levels. (3) Calculate the chemical reaction rate values for the reactions involving different molecules and atoms.

[0058] Ordinary technicians in the industry can also implement the method for measuring the OH ultraviolet spectral radiation characteristics in non-equilibrium flows provided by the present invention using other steps. Figure 1 The method for measuring the OH ultraviolet spectral radiation characteristics in non-equilibrium flows provided by the present invention is only a specific embodiment.

[0059] As Figure 2 shown, the system for measuring the OH ultraviolet spectral radiation characteristics in non-equilibrium flows provided by the present invention includes:

[0060] Gas parameter distribution module 1, which is used to simulate the distributions of the temperature, pressure, and component concentration of the gas flow field of a high-speed aircraft flying at high altitude using simulation software;

[0061] Flow field stratification processing module 2, which is used to stratify the flow field along the determined line-of-sight direction according to the flow field distribution to obtain the temperature, pressure, and concentration distributions of each layer;

[0062] Chemical reaction rate calculation module 3, which is used to calculate the chemical reaction rates in the OH collision-radiation model;

[0063] Number density calculation module 4, which is used to substitute the flow field parameters into the collision-radiation model and programmatically calculate the number densities of different vibrational energy levels of OH(A);

[0064] Spectral radiance calculation module 5, which is used to programmatically calculate the spectral radiance of OH in the wake or shock layer according to the apparent integration method.

[0065] The technical solution of the present invention will be further described below with reference to the accompanying drawings.

[0066] As Figure 3 shown, the method for measuring the OH ultraviolet spectral radiation characteristics in non-equilibrium flows provided by the present invention includes the following steps:

[0067] (1) Select a software from numerous simulation softwares for fluid simulation. First, establish a model similar to a blunted body or a nozzle, and let it fly under high-altitude and hypersonic flight conditions to simulate the temperature, pressure, and concentration characteristics of each component in the shock layer or the afterbody flow field.

[0068] (2) Determine the line-of-sight direction, use programming software to perform stratification processing on the flow field, set an appropriate number of stratifications while ensuring calculation efficiency and accuracy, and obtain the temperature, pressure, and concentration values of each layer.

[0069] (3) Calculate the chemical reaction rate in the OH collision-radiation model;

[0070] (3a) The rate of the displacement reaction strongly depends on the transition process between different vibrational energy levels of OH(A). Among them, the transition processes from high energy levels to low energy levels mainly include ν′ = 1 → ν″ = 0, ν′ = 2 → ν″ = 0, ν′ = 2 → ν″ = 1, ν′ = 3 → ν″ = 2; the corresponding transition processes from low energy levels to high energy levels are ν′ = 0 → ν″ = 1, ν′ = 0 → ν″ = 2, ν′ = 1 → ν″ = 2, ν′ = 2 → ν″ = 3, as Figure 4 shown. During the programming calculation of the chemical reaction rates corresponding to different vibrational energy level transition processes, the original 144×1-dimensional rate matrix in the model is expanded to 144×144 dimensions, as Figure 5 shown. The formula can be used to couple the vibrational energy level and rotational energy level transitions, where h(J′) and g(J″ → J′) represent the rotational energy level J′-dependent function and the rotational energy level transition-dependent function respectively. It should be noted that J″ > J′ should be satisfied in this formula. According to the formula calculate the reverse rate of the displacement reaction, where θ represents the oscillator characteristic temperature, taking 4301.5K; represents the forward rate of the displacement reaction, represents the forward rate of the displacement reaction accurate to the vibrational energy level transition.

[0071] (3b) The rate of the predissociation reaction strongly depends on the rotational energy level. When the vibrational energy level ν′ = 0, the rotational energy level J′ ≥ 22, the vibrational energy level ν′ = 1, the rotational energy level J′ ≥ 11, and all rotational energy levels with the vibrational energy level ν′ ≥ 2 have a strong dependence of the predissociation rate on the rotational quantum number. The formula is used to calculate the predissociation reaction rate values corresponding to different rotational energy levels at the same vibrational energy level. The calculation results are as Figure 6 shown, where f(J′) = c 1 exp[c 2 J′(J′ + 1)], represents the predissociation rate dependent on the vibrational energy level, and c 1 represents the first radiation constant.

[0072] (3c) The reactions involving molecules and atoms include de-excitation reactions and displacement reactions. To make the model applicable to both the wake and the shock layer simultaneously, the main molecules and atoms participating in the reactions are O 2 , H 2 , H 2 O, CO, CO 2 , H, O, and N 2 . Taking O 2 as an example, the de-excitation reaction rate values depending on the vibration ν′ = 0 - 3 are obtained by referring to the literature. Using the formula , all the de-excitation reaction rate values accurate to ν′ = 0 - 3 and J′ = 0 - J′ max (ν′) can be calculated, where f Q，J′ represents the rotational energy level fitting function; the solution of the displacement reaction rate is calculated according to step (3a).

[0073] (4) The flow field parameters required in this model are the concentrations, temperatures, and pressures of the heavy particles participating in the reactions. When flying at low altitude (≤100 km), the particles colliding in the flow field are mainly neutral particles, so the temperature calculated by substituting into the collisional-radiative model is the translational temperature. Substitute the flow field parameters obtained in step (2) into the matrix equation of the collisional-radiative model, and the number density of different vibrational energy levels of OH(A) can be obtained by programming.

[0074] (5) Calculation of the spectral radiation characteristics of OH;

[0075] (5a) Calculate the line strength at different temperatures based on the line-by-line method and the HITRAN high-temperature molecular database where Q(T) represents the partition function at temperature T, E″ represents the energy of the lower state of the transition, v represents the wave number, c 2 represents the second radiation constant, and T 0 represents the reference temperature of 296 K. The absorption coefficient of a single layer α = S·N·Φ can be calculated using the number density of energy levels obtained in step (4), where Φ is a linear function, and a more accurate absorption coefficient can be obtained using the Voigt line shape.

[0076] (5b) Calculate the absorptivity of each layer 1 - e -α·l , where l represents the thickness of a single layer; the transmittance e -α·l .

[0077] (5c) Select the initial radiance I 0 of the blackbody radiation as the target, and calculate the spectral radiance after the entire path based on the apparent integration method , where the optical thickness τ λ = κ sλ ·l, and I bRepresents the blackbody radiation function. It should be noted that when calculating the radiance of the nth layer, the attenuation effect of the paths of the previous n - 1 layers will be passed through, that is, the discrete form of the equation is

[0078] When calculating the OH ultraviolet radiation characteristics in the non-equilibrium flow field of a high-speed aircraft, the present invention calculates the more refined vibrational-rotational energy level transition process of OH(A) by expanding the substitution reaction rate matrix in the collision-radiation model; in addition, the predissociation reaction rate that depends on the vibrational energy level is refined to depend on the rotational energy level, thereby improving the calculation result of the energy level density of the OH excited state and providing a more accurate data input source for subsequent spectral radiance calculations. In addition, the types of participating molecular atoms are increased, improving the applicability of the model.

[0079] The technical effects of the present invention will be described in detail below in combination with simulations.

[0080] 1. Simulation software and basic parameters used in the example

[0081] The simulation software used is Fluent software under Ansys, and ICEM is used for modeling and grid meshing. The flow field of the Atlas II aircraft tail flame is simulated to obtain the temperature field, pressure field, and the distribution fields of the concentrations of each component. Select the direction along the stagnation line, and import the flow field and grid into the tecplot software to perform hierarchical processing on the flow field. Based on the flow field parameters, the matrix equation of the OH collision-radiation model is calculated by programming with matlab software to obtain the population numbers of the OH(A) energy levels.

[0082] Taking the Atlas II aircraft tail flame as an example, the basic parameters of the model flight are as follows:

[0083] Flight altitude: 21 km, flight speed: 2.3 Ma, ambient air pressure: 4.808×10 3 Pa, ambient air temperature: 218 K.

[0084] 2. Data results

[0085] Figure 7 For the comparison of the calculated results of the OH(A) energy level density in the Atlas II aircraft tail flame with the experimental data, it can be seen that the density of the vibrational energy levels at the position 100 - 600 m from the nozzle exit is in good agreement with the experimental data, especially the agreement of the calculated results with ν′ = 2 is the highest. Figure 8 For the comparison of the calculated results of the OH spectral radiance with the experiment, it can be seen that the orders of magnitude of the two are the same. In addition, the standard error of the calculated results is 0.24582 relative to the experimental data in the wavelength range of 280 - 340 nm. As Figure 9As shown, in addition, when calculating the OH spectral radiance of the BSUV-2 shock layer, the standard error reaches 0.06768, indicating that the calculation results of the OH spectral brightness in the shock wave and the wake are quite accurate.

[0086] In summary, the present invention provides an accurate spectral line prediction model, which has important significance for the research of target early warning, detection and recognition. It supplements the fine vibration level transitions of the substitution reactions in the OH collision-radiation model and calculates the refined predissociation reaction rate coupled with the rotational energy levels. It describes the energy level density of the OH excited state more accurately. At the same time, during the programming calculation process, the reaction rate value coupled with the refined level transition is realized by expanding the substitution reaction rate matrix. When using programming to calculate the matrix equation of the collision-radiation model and subsequently solve the spectral radiance based on the apparent integration method, it has the significant advantages of short calculation time and accurate calculation results.

[0087] It should be noted that the embodiments of the present invention can be implemented by hardware, software, or a combination of software and hardware. The hardware part can be implemented using dedicated logic; the software part can be stored in a memory and executed by an appropriate instruction execution system, such as a microprocessor or dedicated design hardware. Those of ordinary skill in the art can understand that the above devices and methods can be implemented using computer-executable instructions and / or included in processor control code, for example, such code is provided on a carrier medium such as a disk, CD, or DVD-ROM, a programmable memory such as a read-only memory (firmware), or a data carrier such as an optical or electronic signal carrier. The devices and modules of the present invention can be implemented by hardware circuits of programmable hardware devices such as very large scale integrated circuits or gate arrays, semiconductors such as logic chips, transistors, etc., or field programmable gate arrays, programmable logic devices, etc., can also be implemented by software executed by various types of processors, or can be implemented by a combination of the above hardware circuits and software such as firmware.

[0088] The above is only the specific implementation manner of the present invention, but the protection scope of the present invention is not limited thereto. Any modification, equivalent replacement, and improvement made within the spirit and principle of the present invention by those skilled in the art within the technical scope disclosed by the present invention shall be covered by the protection scope of the present invention.

Claims

1. A method for measuring the OH ultraviolet spectral radiation characteristics in a non-equilibrium flow, characterized in that, the method for measuring the OH ultraviolet spectral radiation characteristics in the non-equilibrium flow includes: The first step is to use simulation software to simulate the distribution of the gas flow field temperature, pressure and component concentration of a high-speed aircraft flying at high altitude; The second step is to perform stratification processing on the flow field along the determined line-of-sight direction according to the flow field distribution to obtain the temperature, pressure and concentration distribution of each component in each layer; The third step is to calculate the chemical reaction rate in the OH collision-radiation model; The fourth step is to substitute the flow field parameters into the collision-radiation model and program to calculate the number density of different vibrational energy levels under OH(A); The fifth step is to use programming to calculate the spectral radiance of OH in the tail flame or shock layer according to the apparent integration method; The third step calculates the chemical reaction rate in the OH collision-radiation model: (1) According to the refined energy level transition theory of the OH molecule, calculate the substitution reaction rate values corresponding to the transitions between different vibrational energy levels of OH(A), and use programming to solve the 144×144-dimensional substitution reaction rate matrix; (2) Calculate the predissociation reaction rate values corresponding to different rotational energy levels; (3) Calculate the chemical rate values of different molecules and atoms participating in the reaction; Further includes: 1) The rate of the substitution reaction strongly depends on the transition process between different vibrational energy levels of OH(A); among them, the transition processes from high energy levels to low energy levels mainly include ν′ = 1 → ν″ = 0, ν′ = 2 → ν″ = 0, ν′ = 2 → ν″ = 1, ν′ = 3 → ν″ = 2; the corresponding transition processes from low energy levels to high energy levels are ν′ = 0 → ν″ = 1, ν′ = 0 → ν″ = 2, ν′ = 1 → ν″ = 2, ν′ = 2 → ν″ = 3; in the process of programming to calculate the chemical reaction rates corresponding to different vibrational energy level transition processes, the 144×1-dimensional rate matrix in the original model is expanded to 144×144 dimensions; using the formula To couple the vibrational energy level and rotational energy level transitions, it should be noted that J″ > J′ should be satisfied in this formula; according to the formula Calculate the reverse rate of the substitution reaction. 2) The pre-dissociation reaction rate strongly depends on the rotational energy level. When the vibrational energy level ν′ = 0 and the rotational energy level J′ ≥ 22, when the vibrational energy level ν′ = 1 and the rotational energy level J′ ≥ 11, and for all rotational energy levels with the vibrational energy level ν′ ≥ 2, the pre-dissociation rate has a strong dependence on the rotational quantum number; the formula is used to calculate the pre-dissociation reaction rate values corresponding to different rotational energy levels at the same vibrational energy level, where f(J′) = c 1 exp[c 2 J′(J′ + 1)]; 3) The reactions involving molecules and atoms include de-excitation reactions and substitution reactions. The main molecules and atoms participating in the reactions are O 2 , H 2 , H 2 O, CO, CO 2 , H, O, and N 2 ; Obtain the de-excitation reaction rate values depending on the vibration ν′ = 0 - 3, and use the formula to calculate all the de-excitation reaction rate values accurate to ν′ = 0 - 3, J′ = 0 - J′ max (ν′); The substitution reaction rate is solved according to 1) for calculation.

2. The method for measuring the OH ultraviolet spectral radiation characteristics in a non-equilibrium flow according to claim 1, characterized in that, the fifth step uses programming to calculate the spectral radiance of OH in the tail flame or shock layer according to the apparent integration method, including: (1) Calculate the line intensity at different temperatures based on the line-by-line method and the HITRAN high-temperature molecular database The absorption coefficient α of a single layer can be calculated using the energy level density obtained in the fourth step, where α = S·N·Φ, Φ is a linear function, and a more accurate absorption coefficient can be obtained using the Voigt line shape; (2) Calculate the absorption rate of each layer, 1 - e -α·l , where l represents the thickness of a single layer; the transmittance is e -α·l ; (3) Select blackbody radiation as the initial radiance I of the target 0 , and calculate the spectral radiance after the entire path based on the apparent integration method , where the optical thickness τ λ = κ sλ · l; when calculating the radiance of the nth layer, due to the attenuation effect of the path through the previous n - 1 layers, the discrete form of the equation is 3. A computer device, characterized in that, the computer device includes a memory and a processor, the memory stores a computer program, and when the computer program is executed by the processor, the processor executes the steps of the method for measuring the OH ultraviolet spectral radiation characteristics in the non-equilibrium flow according to any one of claims 1 to 2.

4. A computer-readable storage medium storing a computer program, and when the computer program is executed by a processor, the processor executes the steps of the method for measuring the OH ultraviolet spectral radiation characteristics in the non-equilibrium flow according to any one of claims 1 to 2.

5. An information data processing terminal, characterized in that, the information data processing terminal is used to implement the steps of the method for measuring the OH ultraviolet spectral radiation characteristics in the non-equilibrium flow according to any one of claims 1 to 2.

6. A system for measuring the OH ultraviolet spectral radiation characteristics in a non-equilibrium flow implementing the method for measuring the OH ultraviolet spectral radiation characteristics in the non-equilibrium flow according to any one of claims 1 to 2, characterized in that, the system for measuring the OH ultraviolet spectral radiation characteristics in the non-equilibrium flow includes: A gas parameter distribution module for using simulation software to simulate the distribution of the gas flow field temperature, pressure and component concentration of a high-speed aircraft flying at high altitude; A flow field stratification processing module for performing stratification processing on the flow field along the determined line-of-sight direction according to the flow field distribution to obtain the temperature, pressure and concentration distribution of each component in each layer; A chemical reaction rate calculation module for calculating the chemical reaction rate in the OH collision-radiation model; A number density calculation module for substituting the flow field parameters into the collision-radiation model and programming to calculate the number density of different vibrational energy levels of OH(A); A spectral radiance calculation module, which is used to calculate the spectral radiance of OH in the exhaust plume or shock layer by programming according to the apparent integration method.

7. A winged aircraft, characterized in that the winged aircraft is equipped with the OH ultraviolet spectral radiation characteristic determination system according to claim 6.

8. A wingless aircraft, characterized in that the wingless aircraft is equipped with the OH ultraviolet spectral radiation characteristic determination system according to claim 6.

Images