Method, device, equipment and medium for simulating transport properties of mixed gases in high-temperature ionization region

Abstract

The application discloses a high-temperature ionization zone mixed gas transport property simulation method, device, equipment and medium, relates to the field of aerodynamics, and comprises the following steps: obtaining collision integral fitting data of a mixed gas in a high-temperature ionization zone, and calculating component viscosity coefficients based on the collision integral fitting data; the high-temperature ionization zone is an ionization zone with a temperature meeting a preset high-temperature judgment condition; the ion viscosity coefficient in the component viscosity coefficients is corrected to obtain a corrected component viscosity coefficient, and the corrected component viscosity coefficient is used to calculate component translational mode thermal conductivity and component internal mode thermal conductivity; the corrected component viscosity coefficient, the component translational mode thermal conductivity and the component internal mode thermal conductivity are used to determine the mixed gas viscosity coefficient and the mode thermal conductivities in the high-temperature ionization zone; and the mixed gas viscosity coefficient and the mode thermal conductivities are used to realize the transport property simulation of the mixed gas in the high-temperature ionization zone, so that the transport property simulation precision and efficiency of the mixed gas in the high-temperature ionization zone are improved.

Description

Technical Field

[0001] This invention relates to the field of aerodynamics, and in particular to a method, apparatus, equipment, and medium for simulating the transport properties of mixed gases in a high-temperature ionization region. Background Technology

[0002] Transport properties are one of the important physicochemical models in numerical simulations of nonequilibrium thermochemical flows of high-temperature gases, including component transport properties and mixed gas transport properties. For the transport properties of gases undergoing high-temperature thermochemical reactions, the traditional Wilke mixing law model, while simple in calculation and highly efficient, suffers from poor accuracy in simulating the transport properties of mixed gases in the high-temperature ionization region. When ionization reactions are strong, binary collisions are not simply collisions between ionized components, but also include collisions between ions and neutral components. The collision integrals of these two types differ in magnitude; therefore, using the collision integral of ionized components as a representative model will overestimate the overall collision integral level.

[0003] As can be seen from the above, how to improve the accuracy and efficiency of simulating the transport properties of mixed gases in the high-temperature ionization region and reduce the amount of computation is a problem that needs to be solved in this field. Summary of the Invention

[0004] In view of this, the purpose of this invention is to provide a method, apparatus, device, and medium for simulating the transport properties of mixed gases in a high-temperature ionization region, which can improve the accuracy and efficiency of simulating the transport properties of mixed gases in a high-temperature ionization region and reduce the computational load. The specific solution is as follows:

[0005] In a first aspect, this application discloses a method for simulating the transport properties of mixed gases in a high-temperature ionization region, including:

[0006] The collision integral fitting data of the mixed gas in the high-temperature ionization region is obtained, and the viscosity coefficient of the components is calculated based on the collision integral fitting data; the high-temperature ionization region is the ionization region whose temperature meets the preset high-temperature judgment condition.

[0007] The ionic viscosity coefficient in the viscosity coefficient of the component is corrected to obtain the corrected component viscosity coefficient. The translational mode thermal conductivity and the internal mode thermal conductivity of the component are calculated using the corrected component viscosity coefficient.

[0008] The viscosity coefficient and thermal conductivity of each mode of the mixed gas in the high-temperature ionization region are determined based on the corrected component viscosity coefficient, the translational mode thermal conductivity of the component, and the internal mode thermal conductivity of the component.

[0009] The transport properties of the mixed gas in the high-temperature ionization region are simulated by using the viscosity coefficient of the mixed gas and the thermal conductivity of each mode.

[0010] Optionally, calculating the component viscosity coefficient based on the collision integral fitting data includes:

[0011] The collision integral fitting data is calculated based on the collision integral database to obtain the single-component collision integral;

[0012] The component viscosity coefficient is obtained by calculating the collision integral of the single component using the formula for calculating the component viscosity coefficient.

[0013] The formula for calculating the viscosity coefficient of the component is as follows:

[0014] ;

[0015] in, Let i be the viscosity coefficient of component i. Pi It is a universal gas constant. and The temperature and molar mass of component i are respectively. Let Avogadro's constant be 1. This is the collision integral for a single component.

[0016] Optionally, before correcting the ionic viscosity coefficient in the component viscosity coefficient, the method further includes:

[0017] Determine the molar concentration of electrons in the mixed gas in the high-temperature ionization region;

[0018] The molar concentrations of the ionized component and the neutral component are calculated based on the electron molar concentration.

[0019] Optionally, the correction of the ionic viscosity coefficient in the viscosity coefficient of the component includes:

[0020] The viscosity coefficient of the nitrogen-neutral component in the viscosity coefficient of the components is used as the basis for weighted correction, and the ionic viscosity coefficient is corrected using the correction formula.

[0021] The corrected formula is:

[0022] ;

[0023] ;

[0024] ;

[0025] in, The molar concentration of the ionized component. Electron molar concentration, This refers to the molar concentration of the neutral component. It is the ionic viscosity coefficient. The viscosity coefficient of the nitrogen-neutral component. This is the corrected ionic viscosity coefficient.

[0026] Optionally, the step of calculating the translational modal thermal conductivity and the internal modal thermal conductivity of the component using the corrected component viscosity coefficient includes:

[0027] The viscosity coefficient of the corrected component is calculated using the component thermal conductivity calculation formula to obtain the component translational mode thermal conductivity and the component internal mode thermal conductivity. The component translational mode thermal conductivity includes the component translational mode thermal conductivity of heavy particles and the component translational mode thermal conductivity of free electrons. The component internal mode thermal conductivity includes the component rotational mode thermal conductivity of molecules, the component vibrational mode thermal conductivity of molecules, and the component bound electron excited mode thermal conductivity of heavy particles.

[0028] The formula for calculating the thermal conductivity of the component is:

[0029] ;

[0030] in, Let i be the translational mode thermal conductivity of component i. The viscosity coefficient of component i after correction. Let i be the specific heat of the translational mode of component i. It is a universal gas constant. Let i be the molar mass of component i. Let i be the internal modal thermal conductivity of component i. Let be the specific heat of the internal modes of component i.

[0031] Optionally, determining the viscosity coefficient and thermal conductivity of the mixed gas in the high-temperature ionization region based on the corrected component viscosity coefficient, the translational modal thermal conductivity of the component, and the internal modal thermal conductivity of the component includes:

[0032] The viscosity coefficient of the modified components was calculated using the formula for calculating the viscosity coefficient of a mixed gas, thus obtaining the viscosity coefficient of the mixed gas in the high-temperature ionization region.

[0033] The modified component viscosity coefficient, the component translational mode thermal conductivity, and the component internal mode thermal conductivity are calculated using the mixed gas modal thermal conductivity calculation formula to obtain the mixed gas modal thermal conductivity. The mixed gas modal thermal conductivity includes heavy particle translational mode thermal conductivity, rotational mode thermal conductivity, vibrational mode thermal conductivity, bound electron excited mode thermal conductivity, and free electron translational mode thermal conductivity.

[0034] The formula for calculating the viscosity coefficient of the mixed gas is as follows:

[0035] ;

[0036] ;

[0037] in, Where N is the viscosity coefficient of the gas mixture, and Ns is the number of components. and The molar concentrations of components i and j are respectively. and Let i and j be the corrected viscosity coefficients of components i and j, respectively. For weighted functions, and denoted as i and j, respectively, are the molar masses of components i and j.

[0038] Optionally, the formula for calculating the thermal conductivity of the mixed gas modes includes formulas for calculating the thermal conductivity of the translational mode of heavy particles, formulas for calculating the thermal conductivity of the rotational mode, formulas for calculating the thermal conductivity of the vibrational mode, formulas for calculating the thermal conductivity of the bound electron excited mode, and formulas for calculating the thermal conductivity of the free electron translational mode.

[0039] The formula for calculating the thermal conductivity of the translational mode of heavy particles is as follows:

[0040] ;

[0041] ;

[0042] in, Let H be the translational thermal conductivity of the heavy particles in the gas mixture, where H represents the heavy particles. and The molar concentrations of components i and j are respectively. It is the translational mode thermal conductivity of component i. For weighted functions, and Let be the corrected viscosity coefficients of components i and j, respectively. and Let i and j be the molar masses of components i and j, respectively.

[0043] The formula for calculating the thermal conductivity of the rotational mode is:

[0044] ;

[0045] ;

[0046] in, The rotational mode thermal conductivity of the gas mixture. It is the rotational mode thermal conductivity of component i;

[0047] The formula for calculating the thermal conductivity of the vibration mode is:

[0048] ;

[0049] ;

[0050] in, The vibrational mode thermal conductivity of the gas mixture. It is the vibrational mode thermal conductivity of component i;

[0051] The formula for calculating the thermal conductivity of the bound electron excited mode is as follows:

[0052] ;

[0053] ;

[0054] in, The bound electron excited mode thermal conductivity of the mixed gas. It is the bound electron excited mode thermal conductivity of component i;

[0055] The formula for calculating the thermal conductivity of the free electron translational mode is as follows:

[0056] ;

[0057] in, The free electron translational mode thermal conductivity of the gas mixture. It is the translational mode thermal conductivity of the free electron component. The molar concentration of the ionized component. This refers to the molar concentration of the neutral component. and These are the electronic viscosity coefficient and molar mass, respectively. and These are the viscosity coefficient and molar mass of the nitrogen-neutral component, respectively.

[0058] Secondly, this application discloses a device for simulating the transport properties of a mixed gas in a high-temperature ionization region, comprising:

[0059] The component viscosity coefficient calculation module is used to acquire collision integral fitting data of the mixed gas in the high-temperature ionization region, and calculate the component viscosity coefficient based on the collision integral fitting data; the high-temperature ionization region is an ionization region whose temperature meets the preset high-temperature judgment condition.

[0060] The ionic viscosity coefficient correction module is used to correct the ionic viscosity coefficient in the component viscosity coefficient to obtain the corrected component viscosity coefficient, and to calculate the translational mode thermal conductivity and internal mode thermal conductivity of the component using the corrected component viscosity coefficient.

[0061] The viscosity coefficient and thermal conductivity determination module is used to determine the viscosity coefficient and thermal conductivity of the mixed gas in the high-temperature ionization region based on the corrected component viscosity coefficient, the component translational mode thermal conductivity, and the component internal mode thermal conductivity.

[0062] The simulation module is used to simulate the transport properties of the mixed gas in the high-temperature ionization region by utilizing the viscosity coefficient of the mixed gas and the thermal conductivity of each mode.

[0063] Thirdly, this application discloses an electronic device, including:

[0064] Memory, used to store computer programs;

[0065] A processor is used to execute the computer program to implement the aforementioned method for simulating the transport properties of mixed gases in a high-temperature ionization region.

[0066] Fourthly, this application discloses a computer storage medium for storing a computer program; wherein, when the computer program is executed by a processor, it implements the steps of the aforementioned method for simulating the transport properties of a mixed gas in a high-temperature ionization region.

[0067] As can be seen, this application provides a method for simulating the transport properties of a mixed gas in a high-temperature ionization region, including acquiring collision integral fitting data of the mixed gas in the high-temperature ionization region, calculating the component viscosity coefficient based on the collision integral fitting data; the high-temperature ionization region is an ionization region whose temperature meets a preset high-temperature determination condition; correcting the ionic viscosity coefficient in the component viscosity coefficient to obtain a corrected component viscosity coefficient, calculating the component translational modal thermal conductivity and the component internal modal thermal conductivity using the corrected component viscosity coefficient; determining the mixed gas viscosity coefficient and each modal thermal conductivity in the high-temperature ionization region based on the corrected component viscosity coefficient, the component translational modal thermal conductivity, and the component internal modal thermal conductivity; and simulating the transport properties of the mixed gas in the high-temperature ionization region using the mixed gas viscosity coefficient and each modal thermal conductivity. This application calculates the component viscosity coefficient based on collision integral fitting data, corrects the ionic viscosity coefficient in the component viscosity coefficient to solve the problem of poor simulation accuracy of the transport properties of mixed gases in the high-temperature ionization region. Using the corrected component viscosity coefficient, the translational mode thermal conductivity and the internal mode thermal conductivity of the component are calculated. The original Wilke mixing law is used to determine the viscosity coefficient of the mixed gas and the thermal conductivity of each heavy particle mode in the high-temperature ionization region. The corrected free electron thermal conductivity mixing law is used to determine the translational mode thermal conductivity of free electrons. By using the viscosity coefficient of the mixed gas and the thermal conductivity of each mode, the transport properties of the mixed gas in the high-temperature ionization region can be simulated, improving the simulation accuracy and efficiency of the transport properties of the mixed gas in the high-temperature ionization region and reducing the amount of calculation. Attached Figure Description

[0068] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on the provided drawings without creative effort.

[0069] Figure 1 This is a diagram of the equilibrium air viscosity coefficient at standard atmospheric pressure calculated using a mixing law model with different gas transport properties disclosed in this application.

[0070] Figure 2 This is a thermal conductivity diagram at standard atmospheric pressure calculated using a mixing law model of different gas transport properties disclosed in this application.

[0071] Figure 3 This is a distribution diagram of the main components under standard atmospheric pressure as disclosed in this application;

[0072] Figure 4 This is a collision integral diagram between ionized components dominated by Coulomb potential under standard atmospheric pressure, as disclosed in this application.

[0073] Figure 5 This is a collision integral representation of the interaction of non-ionized Coulomb potentials under standard atmospheric pressure as disclosed in this application;

[0074] Figure 6 This is a flowchart of a method for simulating the transport properties of a mixed gas in a high-temperature ionization region, as disclosed in this application.

[0075] Figure 7 This is a comparison chart of the viscosity coefficient calculated using the Wilke correction method proposed in this application with other mixed-law models.

[0076] Figure 8 This is a comparison chart of thermal conductivity calculation results using the Wilke correction method disclosed in this application and calculation results using other hybrid law models.

[0077] Figure 9 This is a schematic diagram of the structure of a device for simulating the transport properties of a mixed gas in a high-temperature ionization region, as disclosed in this application.

[0078] Figure 10 This application provides a structural diagram of an electronic device. Detailed Implementation

[0079] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0080] Transport properties are one of the important physicochemical models in numerical simulations of nonequilibrium thermochemical flows of high-temperature gases, including component transport properties and mixed gas transport properties. For the transport properties of gases undergoing high-temperature thermochemical reactions, the traditional Wilke mixing law model, while simple in calculation and highly efficient, suffers from poor accuracy in simulating the transport properties of mixed gases in the high-temperature ionization region. When ionization reactions are strong, binary collisions are not simply collisions between ionized components; there are also collisions between ions and neutral components, etc. The collision integrals of these two types differ in magnitude. Therefore, using the collision integral of ionized components as a representative model will overestimate the overall collision integral level.

[0081] Based on theoretical derivations such as Hirschfelder (the theory of molecular kinetics and molecular interactions) and Chapman (microscopic transport properties of rarefied gases), the transport properties of components can be calculated based on binary collision integrals and Sonine expansions (series expansions with Sonine polynomials as basis functions).

[0082] For high-temperature gas mixtures, based on the theoretical derivations of Hirschfelder, Chapman, etc., it is necessary to solve the system of equations for the multi-component system using the binary collision integrals of all components and the corresponding Sonine expansion coefficients. The equations are in the following form:

[0083] (1) Viscosity coefficient: ,in, These are the coefficients of the Sonine expansion term. The viscosity coefficient of the mixed gas. Let Ni be the molar concentration of component i, and Ns be the number of components. This is the viscosity coefficient matrix.

[0084] (2) Translational thermal conductivity: ,in, These are the coefficients of the Sonine expansion term. This is the translational thermal conductivity coefficient matrix.

[0085] (3) The theoretical solution for internal thermal conductivity is relatively complex. A more accurate approximation method is to use the mixing law given by Brokaw (classical semi-empirical model) and the Eucken (classical semi-empirical) relation for calculation: , It is Boltzmann's constant. It is a universal gas constant. It is the mass fraction of component i. The internal specific heat of component i corresponds to the specific heat of rotational, vibrational, and bound electron excited degrees of freedom, respectively.

[0086] H ijv and L ij The calculations require the use of binary collision integrals for all components, including diffusion collision integrals. Viscous collision integral and collision integral ratio , These are usually provided in a dedicated database. For ease of calculation, the corresponding equivalent collision integral form is also commonly used:

[0087] ;

[0088] As can be seen, the Chapman theory method not only requires binary collision information for all components but also needs to solve for matrices, resulting in a high computational cost. Therefore, approximation methods were proposed. The most commonly used approximation methods are mainly divided into two categories: Yos's mixing law and Wilke's mixing law.

[0089] The Yos mixing law simplifies the coefficient matrix of a multi-component system equation system diagonally or expands it using a low-order approximation, preserving the binary collision integrals between different components without solving for the matrix. Specifically:

[0090] (1) Viscosity coefficient , It is Avogadro's constant;

[0091] (2) Translational thermal conductivity: ,

[0092] ;

[0093] (3) Internal thermal conductivity, the same as the Chapman theory method.

[0094] Wilke's mixing law, by assuming that different binary collision processes are identical and introducing the Eucken relation, further necessitates considering only the viscous collision integral between the same components. This reduces the amount of data and computation required. Specifically:

[0095] (1) Viscosity coefficient .

[0096] (2) Thermal conductivity: .

[0097] , ;

[0098] Here k i The thermal conductivity of component i in any energy mode, including translational, rotational, vibrational, and bound electron excitation modes, is calculated using the Eucken relation. .in Let i be the translational specific heat of component i; Let be the internal specific heat of component i, and be the specific heat corresponding to rotational, vibrational, and bound electron excited degrees of freedom, respectively.

[0099] The equilibrium air viscosity at standard atmospheric pressure calculated by the mixing law model of different gas transport properties is as follows: Figure 1 As shown, the thermal conductivity at standard atmospheric pressure calculated by the mixing law model for different gas transport properties is as follows: Figure 2 As shown, standard atmospheric pressure is 101325 Pa. From Figure 1 and Figure 2 As can be seen, Wilke's mixing law differs significantly from Chapman's theoretical method after 8000 K. Using Wilke's mixing law in numerical simulation of high-temperature gas flow fields will lead to significant deviations in the numerical simulation results.

[0100] Distribution of main components under standard atmospheric pressure as follows Figure 3 As shown, below 8000 K, dissociation reactions dominate in gases, primarily involving neutral molecules and atoms; above 8000 K, ionization reactions occur in gases, mainly involving atoms and ions resulting from atomic ionization. Therefore, Wilke's mixing law begins to differ significantly from Chapman's theoretical method in the high-temperature ionization reaction region.

[0101] The collision integral between ionized components dominated by the Coulomb potential at standard atmospheric pressure is as follows: Figure 4 As shown, the value of this collision integral is very large. The collision integral of the interaction of the non-ionized Coulomb potential under standard atmospheric pressure represents, for example,... Figure 5 As shown, although the magnitudes of the collision integrals differ between different types of collisions, they are basically on the same order of magnitude and are much smaller than the collision integrals between ionized components.

[0102] Wilke's mixing law uses collision integrals of the same component to characterize all collision integrals. When ionization reactions are strong, neutral components undergo neutral-ion and neutral-neutral collisions, and their collision integrals are close, thus fulfilling the premise of Wilke's mixing law. However, ions do not only undergo collisions between ionized components but also ion-neutral component collisions, and these collision integrals differ in magnitude. Therefore, using the collision integral of ionized components as a representative would overestimate the overall collision integral level, resulting in a significant difference between Wilke's mixing law and Chapman's theoretical method.

[0103] As can be seen from the above, how to improve the accuracy and efficiency of simulating the transport properties of mixed gases in the high-temperature ionization region and reduce the amount of computation is a problem that needs to be solved in this field.

[0104] To improve the prediction accuracy of Wilke's mixing law while taking into account the calculation efficiency of transport properties, this application proposes to modify Wilke's mixing law to simulate the transport properties of mixed gases in the high-temperature ionization region.

[0105] See Figure 6 As shown in the figure, this invention discloses a method for simulating the transport properties of a mixed gas in a high-temperature ionization region, which may specifically include:

[0106] Step S11: Obtain the collision integral fitting data of the mixed gas in the high-temperature ionization region, and calculate the component viscosity coefficient based on the collision integral fitting data; the high-temperature ionization region is the ionization region whose temperature meets the preset high-temperature judgment condition.

[0107] In this embodiment, collision integral fitting data of the mixed gas in the high-temperature ionization region is obtained. The collision integral fitting data is then calculated based on the collision integral database to obtain the single-component collision integral. The single-component collision integral is then calculated using the component viscosity coefficient calculation formula to obtain the component viscosity coefficient. The component viscosity coefficient calculation formula is as follows:

[0108] ;

[0109] in, Let i be the viscosity coefficient of component i. Pi It is a universal gas constant. and The temperature and molar mass of component i are respectively. Let Avogadro's constant be 1. This is the collision integral for a single component.

[0110] In addition, the data acquisition standards for the collision integral of the mixed gas in the high-temperature ionization region can also refer to the data in the Gupta (embedded relational database) and Wright databases.

[0111] Step S12: Correct the ionic viscosity coefficient in the viscosity coefficient of the component to obtain the corrected viscosity coefficient of the component, and use the corrected viscosity coefficient of the component to calculate the translational mode thermal conductivity and the internal mode thermal conductivity of the component.

[0112] In this embodiment, the electron molar concentration of the mixed gas in the high-temperature ionization region is determined; the molar concentration of the ionized component and the molar concentration of the neutral component are calculated based on the electron molar concentration; the nitrogen atom neutral component viscosity coefficient in the component viscosity coefficient is used as the weighted correction basis, and the ionic viscosity coefficient is corrected using the correction formula;

[0113] The corrected formula is:

[0114] ;

[0115] ;

[0116] ;

[0117] in, The molar concentration of the ionized component. Electron molar concentration, This refers to the molar concentration of the neutral component. It is the ionic viscosity coefficient. The viscosity coefficient of the nitrogen-neutral component. This is the corrected ionic viscosity coefficient.

[0118] The corrected component viscosity coefficient includes the corrected ionic viscosity coefficient.

[0119] This application selects the viscosity coefficient of the neutral N (nitrogen) atom component in the gas mixture. As a basis for weighted correction, the ionic viscosity coefficient Make corrections to obtain the corrected ionic viscosity coefficient. The main reasons for weighted correction are as follows: Although there are some differences between the collision integrals of neutral particles and neutral particles and ions, they are relatively small, unlike the difference in magnitude between the collision integrals of ionized components. Therefore, any choice of the collision integrals of neutral particles and neutral particles and ions as the basis for weighted correction can achieve good correction results.

[0120] This application chose atom-to-atom collisions because they are representative and do not require the introduction of new data. Since the main dissociation atoms in high-temperature air are nitrogen atoms, [the following method was chosen]. As the basis for weighted correction.

[0121] Then, the viscosity coefficient of the corrected component is calculated using the component thermal conductivity calculation formula to obtain the translational modal thermal conductivity of the component. and the internal modal thermal conductivity of the components The component translational mode thermal conductivity includes the component translational mode thermal conductivity of heavy particles. Thermal conductivity of free electron composition in translational mode The internal modal thermal conductivity of the component includes the rotational modal thermal conductivity of the molecular structure. Thermal conductivity of molecular component vibration modes Thermal conductivity of bound electron excited modes of heavy particles Heavy particles include neutral particles and ions; molecules include neutral molecules and ionic molecules. Specific heat also has a corresponding classification.

[0122] This application calculates the component thermal conductivity based on the Eucken relation (a classical semi-empirical formula for the transport properties of correlated gases). The corresponding formula for calculating the component thermal conductivity is as follows:

[0123] ;

[0124] in, Let i be the translational mode thermal conductivity of component i. The viscosity coefficient of component i after correction. Let i be the specific heat of the translational mode of component i. It is a universal gas constant. Let i be the molar mass of component i. Let i be the internal modal thermal conductivity of component i. Let be the specific heat of the internal modes of component i.

[0125] Step S13: Determine the viscosity coefficient and thermal conductivity of the mixed gas in the high-temperature ionization region based on the corrected component viscosity coefficient, the translational mode thermal conductivity of the component, and the internal mode thermal conductivity of the component.

[0126] Based on the modified component viscosity coefficient, component translational mode thermal conductivity, and component internal mode thermal conductivity, this application uses the original Wilke mixing law (a classical model used to estimate the transport coefficient of gas mixtures) to calculate the modified mixed gas transport properties and corrects the free electron thermal conductivity.

[0127] In this embodiment, the viscosity coefficient of the modified component is calculated using the formula for calculating the viscosity coefficient of the mixed gas to obtain the viscosity coefficient of the mixed gas in the high-temperature ionization region; the viscosity coefficient of the modified component, the translational mode thermal conductivity of the component, and the internal mode thermal conductivity of the component are calculated using the formula for calculating the modal thermal conductivity of the mixed gas to obtain the modal thermal conductivity of the mixed gas; the modal thermal conductivity of the mixed gas includes the translational mode thermal conductivity of heavy particles, the rotational mode thermal conductivity, the vibrational mode thermal conductivity, the bound electron excited mode thermal conductivity, and the free electron translational mode thermal conductivity;

[0128] The formula for calculating the viscosity coefficient of the mixed gas is as follows:

[0129] ;

[0130] ;

[0131] in, Where N is the viscosity coefficient of the gas mixture, and Ns is the number of components. and The molar concentrations of components i and j are respectively. and Let i and j be the corrected viscosity coefficients of components i and j, respectively. For weighted functions, and denoted as i and j, respectively, are the molar masses of components i and j.

[0132] Specifically, the calculation formulas for the thermal conductivity of mixed gas modes include the calculation formulas for the thermal conductivity of heavy particle translational modes, rotational modes, vibrational modes, bound electron excited modes, and free electron translational modes.

[0133] The formula for calculating the thermal conductivity of the translational mode of heavy particles is as follows:

[0134] ;

[0135] ;

[0136] in, Let H be the translational thermal conductivity of the heavy particles in the gas mixture, where H represents the heavy particles. and The molar concentrations of components i and j are respectively. It is the translational mode thermal conductivity of component i. For weighted functions, and Let be the corrected viscosity coefficients of components i and j, respectively. and Let i and j be the molar masses of components i and j, respectively.

[0137] The formula for calculating the thermal conductivity of the rotational mode is:

[0138] ;

[0139] ;

[0140] in, The rotational mode thermal conductivity of the gas mixture. It is the rotational mode thermal conductivity of component i;

[0141] The formula for calculating the thermal conductivity of the vibration mode is:

[0142] ;

[0143] ;

[0144] in, The vibrational mode thermal conductivity of the gas mixture. It is the vibrational mode thermal conductivity of component i;

[0145] The formula for calculating the thermal conductivity of the bound electron excited mode is as follows:

[0146] ;

[0147] ;

[0148] in, The bound electron excited mode thermal conductivity of the mixed gas. It is the bound electron excited mode thermal conductivity of component i;

[0149] The formula for calculating the thermal conductivity of the free electron translational mode is as follows:

[0150] ;

[0151] in, The free electron translational mode thermal conductivity of the gas mixture. It is the translational mode thermal conductivity of the free electron component. The molar concentration of the ionized component. This refers to the molar concentration of the neutral component. and These are the electronic viscosity coefficient and molar mass, respectively. and These are the viscosity coefficient and molar mass of the nitrogen-neutral component, respectively.

[0152] Step S14: Simulate the transport properties of the mixed gas in the high-temperature ionization region using the viscosity coefficient of the mixed gas and the thermal conductivity of each mode.

[0153] Under standard atmospheric pressure, the viscosity coefficient calculated using the Wilke correction method proposed in this application is compared with the results of other mixing law models, for example... Figure 7 As shown, the thermal conductivity calculation results using the Wilke correction method in this application are compared with those calculated using other hybrid law models, for example... Figure 8 As shown, the Wilke correction method proposed in this application achieves better agreement with the Chapman theoretical method in the high-temperature ionization region, thus improving the accuracy of the original Wilke method.

[0154] This application proposes a method using a weighted average of the viscosity coefficients of the ionized component and the neutral component. This weighted averaged viscosity coefficient is then used to replace the original ionic viscosity coefficient, which is finally substituted into the original Wilke formula to obtain the corrected transport properties of the mixed gas. Furthermore, the mixing law of thermal conductivity in the free electron translational mode is corrected. Compared to the original Wilke formula, this application adds a viscosity coefficient weighting step, with almost no increase in computational load, but improves the simulation accuracy of transport coefficients in the high-temperature ionization reaction region.

[0155] In this embodiment, collision integral fitting data of the mixed gas in the high-temperature ionization region is obtained, and the component viscosity coefficient is calculated based on the collision integral fitting data; the high-temperature ionization region is an ionization region whose temperature meets the preset high-temperature judgment condition; the ionic viscosity coefficient in the component viscosity coefficient is corrected to obtain the corrected component viscosity coefficient, and the component translational mode thermal conductivity and component internal mode thermal conductivity are calculated using the corrected component viscosity coefficient; the viscosity coefficient of the mixed gas in the high-temperature ionization region and the thermal conductivity of each mode are determined based on the corrected component viscosity coefficient, the component translational mode thermal conductivity, and the component internal mode thermal conductivity; the transport properties of the mixed gas in the high-temperature ionization region are simulated using the mixed gas viscosity coefficient and the thermal conductivity of each mode. This application calculates the component viscosity coefficient based on collision integral fitting data, corrects the ionic viscosity coefficient in the component viscosity coefficient to solve the problem of poor simulation accuracy of the transport properties of mixed gases in the high-temperature ionization region. Using the corrected component viscosity coefficient, the translational mode thermal conductivity and the internal mode thermal conductivity of the component are calculated. The original Wilke mixing law is used to determine the viscosity coefficient of the mixed gas and the thermal conductivity of each heavy particle mode in the high-temperature ionization region. The corrected free electron thermal conductivity mixing law is used to determine the translational mode thermal conductivity of free electrons. By using the viscosity coefficient of the mixed gas and the thermal conductivity of each mode, the transport properties of the mixed gas in the high-temperature ionization region can be simulated, improving the simulation accuracy and efficiency of the transport properties of the mixed gas in the high-temperature ionization region and reducing the amount of calculation.

[0156] See Figure 9 As shown in the figure, this invention discloses a device for simulating the transport properties of a mixed gas in a high-temperature ionization region, which may specifically include:

[0157] The component viscosity coefficient calculation module 11 is used to acquire collision integral fitting data of the mixed gas in the high-temperature ionization region and calculate the component viscosity coefficient based on the collision integral fitting data; the high-temperature ionization region is an ionization region whose temperature meets the preset high-temperature judgment condition.

[0158] The ionic viscosity coefficient correction module 12 is used to correct the ionic viscosity coefficient in the component viscosity coefficient to obtain the corrected component viscosity coefficient, and to calculate the translational mode thermal conductivity and the internal mode thermal conductivity of the component using the corrected component viscosity coefficient.

[0159] The viscosity coefficient and thermal conductivity determination module 13 is used to determine the viscosity coefficient and thermal conductivity of the mixed gas in the high-temperature ionization region based on the corrected component viscosity coefficient, the component translational mode thermal conductivity, and the component internal mode thermal conductivity.

[0160] The simulation module 14 is used to simulate the transport properties of the mixed gas in the high-temperature ionization region using the viscosity coefficient and the thermal conductivity of each mode.

[0161] In this embodiment, collision integral fitting data of the mixed gas in the high-temperature ionization region is obtained, and the component viscosity coefficient is calculated based on the collision integral fitting data; the high-temperature ionization region is an ionization region whose temperature meets the preset high-temperature judgment condition; the ionic viscosity coefficient in the component viscosity coefficient is corrected to obtain the corrected component viscosity coefficient, and the component translational mode thermal conductivity and component internal mode thermal conductivity are calculated using the corrected component viscosity coefficient; the viscosity coefficient of the mixed gas in the high-temperature ionization region and the thermal conductivity of each mode are determined based on the corrected component viscosity coefficient, the component translational mode thermal conductivity, and the component internal mode thermal conductivity; the transport properties of the mixed gas in the high-temperature ionization region are simulated using the mixed gas viscosity coefficient and the thermal conductivity of each mode. This application calculates the component viscosity coefficient based on collision integral fitting data, corrects the ionic viscosity coefficient in the component viscosity coefficient to solve the problem of poor simulation accuracy of the transport properties of mixed gases in the high-temperature ionization region. Using the corrected component viscosity coefficient, the translational mode thermal conductivity and the internal mode thermal conductivity of the component are calculated. The original Wilke mixing law is used to determine the viscosity coefficient of the mixed gas and the thermal conductivity of each heavy particle mode in the high-temperature ionization region. The corrected free electron thermal conductivity mixing law is used to determine the translational mode thermal conductivity of free electrons. By using the viscosity coefficient of the mixed gas and the thermal conductivity of each mode, the transport properties of the mixed gas in the high-temperature ionization region can be simulated, improving the simulation accuracy and efficiency of the transport properties of the mixed gas in the high-temperature ionization region and reducing the amount of calculation.

[0162] In some specific embodiments, the component viscosity coefficient calculation module 11 may specifically include:

[0163] A single-component collision integral calculation module is used to calculate the collision integral fitting data based on the collision integral database to obtain the single-component collision integral.

[0164] The component viscosity coefficient calculation module is used to calculate the single-component collision integral using the component viscosity coefficient calculation formula to obtain the component viscosity coefficient.

[0165] The formula for calculating the viscosity coefficient of the component is as follows:

[0166] ;

[0167] in, Let i be the viscosity coefficient of component i. Pi It is a universal gas constant. and The temperature and molar mass of component i are respectively. Let Avogadro's constant be 1. This is the collision integral for a single component.

[0168] In some specific embodiments, the ionic viscosity coefficient correction module 12 may specifically include:

[0169] The electron molar concentration determination module is used to determine the electron molar concentration of the mixed gas in the high-temperature ionization region;

[0170] The module for calculating the molar concentration of the ionized component and the neutral component is used to calculate the molar concentration of the ionized component and the molar concentration of the neutral component based on the electron molar concentration.

[0171] In some specific embodiments, the ionic viscosity coefficient correction module 12 may specifically include:

[0172] The weighted correction module is used to take the viscosity coefficient of the nitrogen-neutral component in the viscosity coefficient of the components as the basis for weighted correction, and to correct the ionic viscosity coefficient using the correction formula;

[0173] The corrected formula is:

[0174] ;

[0175] ;

[0176] ;

[0177] in, The molar concentration of the ionized component. Electron molar concentration, This refers to the molar concentration of the neutral component. It is the ionic viscosity coefficient. The viscosity coefficient of the nitrogen-neutral component. This is the corrected ionic viscosity coefficient.

[0178] In some specific embodiments, the viscosity coefficient and thermal conductivity determination module 13 may specifically include:

[0179] The thermal conductivity calculation module is used to calculate the corrected component viscosity coefficient using the component thermal conductivity calculation formula to obtain the component translational mode thermal conductivity and the component internal mode thermal conductivity. The component translational mode thermal conductivity includes the component translational mode thermal conductivity of heavy particles and the component translational mode thermal conductivity of free electrons. The component internal mode thermal conductivity includes the component rotational mode thermal conductivity of molecules, the component vibrational mode thermal conductivity of molecules, and the component bound electron excited mode thermal conductivity of heavy particles.

[0180] The formula for calculating the thermal conductivity of the component is:

[0181] ;

[0182] in, Let i be the translational mode thermal conductivity of component i. The viscosity coefficient of component i after correction. Let i be the specific heat of the translational mode of component i. It is a universal gas constant. Let i be the molar mass of component i. Let i be the internal modal thermal conductivity of component i. Let be the specific heat of the internal modes of component i.

[0183] In some specific embodiments, the viscosity coefficient and thermal conductivity determination module 13 may specifically include:

[0184] The mixed gas viscosity coefficient calculation module is used to calculate the viscosity coefficient of the corrected components using the mixed gas viscosity coefficient calculation formula, so as to obtain the viscosity coefficient of the mixed gas in the high temperature ionization region.

[0185] The module for calculating the thermal conductivity of each mode is used to calculate the corrected component viscosity coefficient, the component translational mode thermal conductivity, and the component internal mode thermal conductivity using the mixed gas mode thermal conductivity calculation formula, thereby obtaining the thermal conductivity of each mode of the mixed gas; the mixed gas mode thermal conductivity includes heavy particle translational mode thermal conductivity, rotational mode thermal conductivity, vibrational mode thermal conductivity, bound electron excited mode thermal conductivity, and free electron translational mode thermal conductivity;

[0186] The formula for calculating the viscosity coefficient of the mixed gas is as follows:

[0187] ;

[0188] ;

[0189] in, Where N is the viscosity coefficient of the gas mixture, and Ns is the number of components. and The molar concentrations of components i and j are respectively. and Let i and j be the corrected viscosity coefficients of components i and j, respectively. For weighted functions, and denoted as i and j, respectively, are the molar masses of components i and j.

[0190] In some specific embodiments, the formula for calculating the thermal conductivity of the mixed gas modes includes formulas for calculating the thermal conductivity of the translational mode of heavy particles, formulas for calculating the thermal conductivity of the rotational mode, formulas for calculating the thermal conductivity of the vibrational mode, formulas for calculating the thermal conductivity of the bound electron excited mode, and formulas for calculating the thermal conductivity of the free electron translational mode.

[0191] The formula for calculating the thermal conductivity of the translational mode of heavy particles is as follows:

[0192] ;

[0193] ;

[0194] in, Let H be the translational thermal conductivity of the heavy particles in the gas mixture, where H represents the heavy particles. and The molar concentrations of components i and j are respectively. It is the translational mode thermal conductivity of component i. For weighted functions, and Let be the corrected viscosity coefficients of components i and j, respectively. and Let i and j be the molar masses of components i and j, respectively.

[0195] The formula for calculating the thermal conductivity of the rotational mode is:

[0196] ;

[0197] ;

[0198] in, The rotational mode thermal conductivity of the gas mixture. It is the rotational mode thermal conductivity of component i;

[0199] The formula for calculating the thermal conductivity of the vibration mode is:

[0200] ;

[0201] ;

[0202] in, The vibrational mode thermal conductivity of the gas mixture. It is the vibrational mode thermal conductivity of component i;

[0203] The formula for calculating the thermal conductivity of the bound electron excited mode is as follows:

[0204] ;

[0205] ;

[0206] in, The bound electron excited mode thermal conductivity of the mixed gas. It is the bound electron excited mode thermal conductivity of component i;

[0207] The formula for calculating the thermal conductivity of the free electron translational mode is as follows:

[0208] ;

[0209] in, The free electron translational mode thermal conductivity of the gas mixture. It is the translational mode thermal conductivity of the free electron component. The molar concentration of the ionized component. This refers to the molar concentration of the neutral component. and These are the electronic viscosity coefficient and molar mass, respectively. and These are the viscosity coefficient and molar mass of the nitrogen-neutral component, respectively.

[0210] Figure 10 This is a schematic diagram of an electronic device provided in an embodiment of this application. The electronic device 20 may specifically include: at least one processor 21, at least one memory 22, a power supply 23, a communication interface 24, an input / output interface 25, and a communication bus 26. The memory 22 stores a computer program, which is loaded and executed by the processor 21 to implement the relevant steps in the simulation method for the transport properties of mixed gases in a high-temperature ionization region, performed by the electronic device as disclosed in any of the foregoing embodiments.

[0211] In this embodiment, the power supply 23 is used to provide operating voltage for each hardware device on the electronic device 20; the communication interface 24 can create a data transmission channel between the electronic device 20 and external devices, and the communication protocol it follows can be any communication protocol applicable to the technical solution of this application, and is not specifically limited here; the input / output interface 25 is used to acquire external input data or output data to the outside world, and its specific interface type can be selected according to specific application needs, and is not specifically limited here.

[0212] In addition, the memory 22, as a carrier for resource storage, can be a read-only memory, random access memory, disk or optical disk, etc. The resources stored on it include operating system 221, computer program 222 and data 223, etc., and the storage method can be temporary storage or permanent storage.

[0213] The operating system 221 manages and controls the various hardware devices on the electronic device 20 and the computer program 222 to enable the processor 21 to perform calculations and processing on the data 223 in the memory 22. It can be Windows, Unix, Linux, etc. The computer program 222, in addition to including a computer program capable of performing the simulation method for the transport properties of a mixed gas in a high-temperature ionization region executed by the electronic device 20 as disclosed in any of the foregoing embodiments, may further include computer programs capable of performing other specific tasks. The data 223 may include data received by the simulation device for the transport properties of a mixed gas in a high-temperature ionization region from external devices, and may also include data collected by its own input / output interface 25.

[0214] The steps of the methods or algorithms described in conjunction with the embodiments disclosed herein can be implemented directly by hardware, a software module executed by a processor, or a combination of both. The software module can be located in random access memory (RAM), main memory, read-only memory (ROM), electrically programmable ROM, electrically erasable programmable ROM, registers, hard disk, removable disk, CD-ROM, or any other form of storage medium known in the art.

[0215] Furthermore, this application also discloses a computer-readable storage medium storing a computer program. When the computer program is loaded and executed by a processor, it implements the steps of the method for simulating the transport properties of mixed gases in a high-temperature ionization region disclosed in any of the foregoing embodiments.

[0216] Finally, it should be noted that in this document, relational terms such as "first" and "second" are used only to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus. Without further limitations, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or apparatus that includes said element.

[0217] The foregoing has provided a detailed description of the method, apparatus, equipment, and storage medium for simulating the transport properties of a mixed gas in a high-temperature ionization region provided by the present invention. Specific examples have been used to illustrate the principles and implementation methods of the present invention. The descriptions of the above embodiments are only for the purpose of helping to understand the method and core ideas of the present invention. At the same time, for those skilled in the art, there will be changes in the specific implementation methods and application scope based on the ideas of the present invention. Therefore, the content of this specification should not be construed as a limitation of the present invention.

Claims

1. A method of simulating transport properties of a mixture of gases in a high temperature ionization region, characterized by, The method comprises the following steps: obtaining collision integral fitting data of a mixed gas in a high-temperature ionization region, calculating component viscosity coefficients based on the collision integral fitting data; the high-temperature ionization region is an ionization region with a temperature satisfying a preset high-temperature judgment condition; correcting ion viscosity coefficients in the component viscosity coefficients to obtain corrected component viscosity coefficients, and calculating component translational mode thermal conductivities and component internal mode thermal conductivities by using the corrected component viscosity coefficients; determining a mixed gas viscosity coefficient and each mode thermal conductivity in the high-temperature ionization region based on the corrected component viscosity coefficients, the component translational mode thermal conductivities and the component internal mode thermal conductivities; realizing simulation of transport properties of the mixed gas in the high-temperature ionization region by using the mixed gas viscosity coefficient and each mode thermal conductivity; wherein the correction of the ion viscosity coefficients in the component viscosity coefficients comprises taking a nitrogen atom neutral component viscosity coefficient in the component viscosity coefficients as a weighted correction basis, and correcting the ion viscosity coefficients by using a correction formula; the correction formula is: ； ； ； wherein, is the molar concentration of the ionized component, is the molar concentration of the electrons, is the molar concentration of the neutral component, is the ionic viscosity coefficient, is the viscosity coefficient of the neutral component of the nitrogen atom, is the corrected ionic viscosity coefficient; the calculation of the component translational mode thermal conductivities and the component internal mode thermal conductivities by using the corrected component viscosity coefficients comprises calculating the corrected component viscosity coefficients by using a component thermal conductivity calculation formula to obtain the component translational mode thermal conductivities and the component internal mode thermal conductivities; the component translational mode thermal conductivities comprise component translational mode thermal conductivities of heavy particles and component translational mode thermal conductivities of free electrons; the component internal mode thermal conductivities comprise component rotational mode thermal conductivities of molecules, component vibration mode thermal conductivities of molecules and component bound electron excitation mode thermal conductivities of heavy particles.

2. The method for simulating transport properties of a gas mixture in a high temperature ionization region according to claim 1, wherein the calculation of the component viscosity coefficients based on the collision integral fitting data comprises: calculating the collision integral fitting data based on a collision integral database to obtain single-component collision integrals; calculating the single-component collision integrals by using a component viscosity coefficient calculation formula to obtain the component viscosity coefficients; the component viscosity coefficient calculation formula is: ； wherein is the viscosity coefficient of component i, is the circle constant, is the universal gas constant, and are the temperature and the molar mass of component i, respectively, is the Avogadro constant, is the single component collision integral.

3. The method for simulating transport properties of a gas mixture in a high temperature ionization region according to claim 1, wherein before the correction of the ion viscosity coefficients in the component viscosity coefficients, the method further comprises: determining an electron molar concentration of the mixed gas in the high-temperature ionization region; calculating ion component molar concentrations and neutral component molar concentrations based on the electron molar concentration.

4. The method for simulating transport properties of a gas mixture in a high temperature ionization region according to claim 1, wherein the component thermal conductivity calculation formula is: ； wherein, is the translational modal thermal conductivity of component i, is the viscosity coefficient of component i after correction, is the specific heat of the translational modal of component i, is the universal gas constant, is the molar mass of component i, is the internal modal thermal conductivity of component i, is the specific heat of the internal modal of component i.

5. The method for simulating transport properties of a gas mixture in a high-temperature ionization region according to any one of claims 1 to 4, characterized in that, the determination of the mixed gas viscosity coefficient and each mode thermal conductivity in the high-temperature ionization region based on the corrected component viscosity coefficients, the component translational mode thermal conductivities and the component internal mode thermal conductivities comprises: calculating the corrected component viscosity coefficients by using a mixed gas viscosity coefficient calculation formula to obtain the mixed gas viscosity coefficient in the high-temperature ionization region; calculating the corrected component viscosity coefficients, the component translational mode thermal conductivities and the component internal mode thermal conductivities by using a mixed gas mode thermal conductivity calculation formula to obtain each mode thermal conductivity of the mixed gas; the mixed gas mode thermal conductivities comprise heavy particle translational mode thermal conductivities, rotational mode thermal conductivities, vibration mode thermal conductivities, bound electron excitation mode thermal conductivities and free electron translational mode thermal conductivities; the mixed gas viscosity coefficient calculation formula is: ； ； wherein, is the viscosity coefficient of the mixture, Ns is the number of components, and are the molar concentrations of components i and j, respectively, and are the modified viscosity coefficients of components i and j, respectively, is a weighting function, and are the molar masses of components i and j, respectively.

6. The method for simulating transport properties of a gas mixture in a high-temperature ionization region according to claim 5, wherein The mixed gas modal thermal conductivity calculation formula comprises a heavy particle translational modal thermal conductivity calculation formula, a rotational modal thermal conductivity calculation formula, a vibrational modal thermal conductivity calculation formula, a bound electron excitation modal thermal conductivity calculation formula and a free electron translational modal thermal conductivity calculation formula. The heavy particle translational modal thermal conductivity calculation formula is: ； ； wherein, H is the heavy particle, and are the molar concentrations of components i and j, respectively, is the translational modal thermal conductivity of component i, is a weighting function, and are the modified viscosity coefficients of components i and j, respectively, and are the molar masses of components i and j, respectively. The rotational modal thermal conductivity calculation formula is: ； ； wherein, is the rotational modal thermal conductivity of the mixture gas, is the rotational modal thermal conductivity of component i; The vibrational modal thermal conductivity calculation formula is: ； ； wherein is the vibrational modal thermal conductivity of the mixture gas, is the vibrational modal thermal conductivity of component i; The bound electron excitation modal thermal conductivity calculation formula is: ； ； wherein, is the bound electron excitation mode thermal conductivity of the mixture gas, is the bound electron excitation mode thermal conductivity of component i; The free electron translational modal thermal conductivity calculation formula is: ； wherein, is the free electron translational mode thermal conductivity of the mixture gas, is the constituent translational mode thermal conductivity of the free electrons, is the ionized constituent molar concentration, is the neutral constituent molar concentration, and are the electronic viscosity coefficient and molar mass, respectively, and are the nitrogen atom neutral constituent viscosity coefficient and molar mass, respectively.

7. An apparatus for simulating transport properties of a mixture of gases in a high temperature ionization region, characterized by, Comprise: The component viscosity coefficient calculation module is used for obtaining collision integral fitting data of the mixed gas in the high-temperature ionization region, and calculating the component viscosity coefficient based on the collision integral fitting data; the high-temperature ionization region is an ionization region with a temperature satisfying a preset high-temperature determination condition; The ion viscosity coefficient correction module is used for correcting the ion viscosity coefficient in the component viscosity coefficient to obtain a corrected component viscosity coefficient, and calculating the component translational modal thermal conductivity and the component internal modal thermal conductivity by using the corrected component viscosity coefficient; The viscosity coefficient and thermal conductivity determination module is used for determining the viscosity coefficient and the modal thermal conductivity of the mixed gas in the high-temperature ionization region based on the corrected component viscosity coefficient, the component translational modal thermal conductivity and the component internal modal thermal conductivity; The simulation module is used for simulating the transport property of the mixed gas in the high-temperature ionization region by using the viscosity coefficient and the modal thermal conductivity of the mixed gas; The correction of the ion viscosity coefficient in the component viscosity coefficient comprises: taking the nitrogen atom neutral component viscosity coefficient in the component viscosity coefficient as a weighted correction basis, and correcting the ion viscosity coefficient by using a correction formula. The correction formula is: ； ； ； wherein, is the molar concentration of the ionized component, is the molar concentration of the electrons, is the molar concentration of the neutral component, is the ionic viscosity coefficient, is the neutral component viscosity coefficient of the nitrogen atom, is the corrected ionic viscosity coefficient; The calculation of the component translational modal thermal conductivity and the component internal modal thermal conductivity by using the corrected component viscosity coefficient comprises: calculating the corrected component viscosity coefficient by using a component thermal conductivity calculation formula to obtain the component translational modal thermal conductivity and the component internal modal thermal conductivity; the component translational modal thermal conductivity comprises the component translational modal thermal conductivity of the heavy particle and the component translational modal thermal conductivity of the free electron; and the component internal modal thermal conductivity comprises the component rotational modal thermal conductivity of the molecule, the component vibrational modal thermal conductivity of the molecule and the component bound electron excitation modal thermal conductivity of the heavy particle.

8. An electronic device, comprising: Comprise: The memory is used for saving the computer program; The processor is used for executing the computer program to realize the transport property simulation method of the mixed gas in the high-temperature ionization region.

9. A computer-readable storage medium, characterized in that, The memory is used for saving the computer program; wherein the computer program is executed by the processor to realize the transport property simulation method of the mixed gas in the high-temperature ionization region. The memory is used for saving the computer program; wherein the computer program is executed by the processor to realize the transport property simulation method of the mixed gas in the high-temperature ionization region.

Images