Method for correcting differential pressure values, device and analyzer for analyzing concentration thereof

Abstract

The invention discloses a method for correcting differential pressure values of a fire extinguishing agent concentration analyzer in the engine compartment of an aircraft, the analyzer, and a readable storage media thereof. The method includes: correcting theoretical differential pressure values of the analyzer; measuring actual differential pressure values of the analyzer under different pressure; obtaining a series of corrected viscosity coefficients accordingly; obtaining a relation accordingly; setting a corrected model to replace the original model of the analyzer. It is solved for the problem that the noise influence caused by the pressure of flight environment cannot be overcome in the test of extinguishing agent concentration in aircraft engine compartment.

Claims

1. A method for correcting differential pressure values of a fire extinguishing agent concentration analyzer in an engine compartment of an aircraft, executed by a processor of a computer terminal, a differential pressure value being obtained by the analyzer using an original model related to a concentration of fire extinguishing agent and a differential pressure, the method comprising: step 1, obtaining a series of flow coefficients Cd according to mass flow-rates Q.sub.m of the fire extinguishing agent at an inlet of the analyzer under different pressures: $C_d=\frac{Q_m}{Q_{m0}}$ Q.sub.m0 meaning theoretical mass flow-rates of the analyzer; step 2, correcting theoretical differential pressure values P.sub.d0 according to C.sub.d, and obtaining corrected theoretical differential pressure values P.sub.d0 accordingly: $P_{d0}^{}=C_d.Math.P_{d0}$ $P_{d0}=\frac{32d_0^{2}l}{{\mathrm{nd}}^4}.Math.\sqrt{T_0}.Math.\frac{\sqrt{R_c/M}}{{\left(1+\frac{-1}{2}\right)}^{\frac{\left(+1\right)}{2\left(-1\right)}}}$ in the formula, meaning a theoretical dynamic viscosity of the fire extinguishing agent, meaning a theoretical gas specific heat ratio, M meaning a theoretical molecular mass; step 3, measuring actual differential pressure values P.sub.d of the analyzer under different pressures; step 4, obtaining a series of corrected viscosity coefficients C.sub. according to P.sub.d and P.sub.d0:
P.sub.d=C.sub.*P.sub.d0 step 5, obtaining a relation C.sub.(x,P) by fitting the corrected viscosity coefficients C.sub., C.sub.(x,P) being a function of C.sub. with regard to a pressure P and a volumetric concentration x of the fire extinguishing agent; step 6, setting a corrected model P.sub.d(x,P) with regard to P and x to replace the original model: ${P_d\left(x,P\right)}^{+}=C_d{C}_{}\left(x,P\right)\frac{32d_0^{2}l}{{\mathrm{nd}}^4}.Math.\sqrt{T_0}.Math.\frac{\left(x\right)\sqrt{\left(x\right)R_c/M\left(x\right)}}{{\left(1+\frac{\left(x\right)-1}{2}\right)}^{\frac{\left(\left(x\right)+1\right)}{2\left(\left(x\right)-1\right)}}}$ In the formula, d.sub.0 meaning an aperture of an internal current limiting element of the analyzer; d meaning an inner diameter of a capillary; l meaning a length of an internal capillary of the analyzer; n meaning a total number of the capillary; T.sub.0 meaning a stagnation temperature; (x) meaning a dynamic viscosity of the fire extinguishing agent about the volumetric concentration x; (x) meaning a theoretical gas specific heat ratio about the volumetric concentration x; M(x) meaning a theoretical molecular mass about the volumetric concentration x; R.sub.c meaning a universal gas constant, wherein in the step 5, the method for fitting C.sub.(x,P) linearly comprises: fitting C.sub. linearly by relation curves with regard to C.sub. and P of the fire extinguishing agent under different x:
C.sub.=a.Math.P+b in the formula, a meaning a slope of C.sub.; b meaning an intercept of C.sub.; fitting a(x) by relation curves with many a under different concentrations, and fitting b(x) by relation curves with many b under different concentrations; obtaining C.sub.(x,P) by different C under different a(x) and b(x):
C.sub.(x,P)=a(x).Math.P+b(x); and Step 7, correcting the differential pressure values using the corrected model P.sub.d(x,P).

2. The method for correcting differential pressure values according to claim 1, wherein in the step 6, following formulas are established if a pressure influence is ignored: $\left(x\right)=\frac{{\mathrm{xc}}_{p1}+\left(1-x\right)c_{p2}}{{\mathrm{xc}}_{p1}+\left(1-x\right)c_{p2}-R}$ $M\left(x\right)={\mathrm{xM}}_1+\left(1-x\right)M_2$ $\left(x\right)=\frac{x{}_1}{x+\left(1-x\right){}_{12}}+\frac{\left(1-x\right){}_2}{\left(1-x\right)+x{}_{21}}$ ${}_{12}=\frac{{\left[1+{\left(\frac{{}_1}{{}_2}\right)}^{\frac{1}{2}}{\left(\frac{M_2}{M_1}\right)}^{\frac{1}{4}}\right]}^2}{{\left[8\left(1+\frac{M_1}{M_2}\right)\right]}^{\frac{1}{2}}}$ ${}_{21}=\frac{{}_2}{{}_1}\frac{M_1}{M_2}{}_{12}$ in the formula, c.sub.p1 meaning a specific heat at constant pressure of the fire extinguishing agent with pure component; c.sub.p2 meaning a specific heat at constant pressure of a carrier gas; R meaning a gas constant, and R=R.sub.c/M; M.sub.1 meaning a molar mass of the fire extinguishing agent with pure component; M.sub.2 meaning a molar mass of the carrier gas; .sub.1 meaning a dynamic viscosity of the fire extinguishing agent with pure component; .sub.2 meaning a dynamic viscosity of the carrier gas.

3. The method for correcting differential pressure values according to claim 1, wherein Q.sub.m0 is shown as follow: $Q_{m0}=\frac{d_0^2}{4}\frac{P_0}{\sqrt{T_0}}\sqrt{\frac{M}{R_c}}\frac{1}{{\left(1+\frac{-1}{2}\right)}^{\frac{\left(+1\right)}{2\left(-1\right)}}}$ in the formula, P.sub.0 meaning a stagnation pressure.

4. The method for correcting differential pressure values according to claim 3, wherein P.sub.0 means a pressure of the fire extinguishing agent at the entrance of the analyzer, T.sub.0 means a gas flow temperature.

5. The method for correcting differential pressure values according to claim 1, wherein the analyzer uses a standard gas prepared by a gravimetric method, the standard gas is a binary-component gas mixed by a carrier gas and the fire extinguishing agent with pure component, the carrier gas is air or nitrogen.

6. A device for correcting differential pressure values of a fire extinguishing agent concentration analyzer in the engine compartment of an aircraft by using the method of claim 1, comprising: a first obtaining unit, electrically connected with the processor, and used for obtaining a series of flow coefficients C.sub.d according to the mass flow-rates Q.sub.m of the fire extinguishing agent at the inlet of the analyzer under different pressures: $C_d=\frac{Q_m}{Q_{m0}}$ Q.sub.m0 meaning the theoretical mass flow-rates of the analyzer; a correcting unit, electrically connected with the processor, and used for correcting theoretical differential pressure values P.sub.d0 according to C.sub.d, and obtaining corrected theoretical differential pressure values P.sub.d0 accordingly: $P_{d0}^{}=C_d.Math.P_{d0}$ $P_{d0}=\frac{32d_0^{2}l}{{\mathrm{nd}}^4}.Math.\sqrt{T_0}.Math.\frac{\sqrt{R_c/M}}{{\left(1+\frac{-1}{2}\right)}^{\frac{\left(+1\right)}{2\left(-1\right)}}}$ in the formula, meaning the theoretical dynamic viscosity of the fire extinguishing agent, meaning the theoretical gas specific heat ratio, M meaning the theoretical molecular mass; a measuring unit, electrically connected with the processor, and used for measuring actual differential pressure values P.sub.d by the analyzer under different pressure pressures; a second obtaining unit, electrically connected with the processor, and used for obtaining a series of corrected viscosity coefficients C.sub. according to P.sub.d and P.sub.d0:
P.sub.d=C.sub..Math.P.sub.d0 a third obtaining unit, electrically connected with the processor, and used for obtaining a relation C.sub.(x,P) by fitting the corrected viscosity coefficients C.sub., C.sub.(x,P) being a function of C.sub. with regard to the pressure P and the volumetric concentration x of the fire extinguishing agent; a setting unit, electrically connected with the processor, and used for setting a corrected model P.sub.d(x,P) with regard to P and x to replace the original model: $P_d\left(x,P\right)=C_d{C}_{}\left(x,P\right)\frac{32d_0^{2}l}{{\mathrm{nd}}^4}.Math.\sqrt{T_0}.Math.\frac{\left(x\right)\sqrt{\left(x\right)R_c/M\left(x\right)}}{{\left(1+\frac{\left(x\right)-1}{2}\right)}^{\frac{\left(\left(x\right)+1\right)}{2\left(\left(x\right)-1\right)}}}$ in the formula, d.sub.0 meaning the aperture diameter of the internal current limiting element of the analyzer; l meaning the length of the internal capillary of the analyzer; n meaning the total number of the capillary; T.sub.0 meaning the stagnation temperature; (x) meaning the dynamic viscosity of the fire extinguishing agent about the volumetric concentration x; (x) meaning the theoretical gas specific heat ratio about the volumetric concentration x; M(x) meaning the theoretical molecular mass about the volumetric concentration x; R.sub.c meaning a universal gas constant.

7. An analyzer for analyzing the fire extinguishing agent concentration, in the engine compartment of an aircraft, and using the method of claim 1 for correcting its differential pressure values.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 is a flow diagram of a method for correcting differential pressure values of a fire extinguishing agent concentration analyzer, according to the first embodiment of the present invention.

(2) FIG. 2 is a graph of relation with flow coefficients C.sub.d under different concentrations and corresponding pressures, according to the first embodiment of the present invention.

(3) FIG. 3 is a flow diagram of a method for setting the model with related to the concentration of fire extinguishing agent and the differential pressure, according to the second embodiment of the present invention.

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS

(4) The present invention will now be described more specifically with reference to the following embodiments. It is to be noted that the following descriptions of preferred embodiments of this invention are presented herein for illustration and description only. It is not intended to be exhaustive or to be limited to the precise disclosed form.

(5) It should be noted that when a component is said to be installed in another component, it can be directly on the other component or it can also exist in a resident component. When a component is considered to be set on another component, it may be set directly on the other component or there may also be a centered component. When a component is considered to be fixed to another component, it may be fixed directly to another component or there may also be a centered component.

(6) Unless otherwise defined, all technical and scientific terms used herein have the same meaning as would normally be understood by a person skilled in the technical field belonging to the present invention. The terms used herein in the specification of the invention are for the purpose of describing specific embodiments only and are not intended to limit the invention. The term or/and as used herein includes any and all combinations of one or more related listed items.

(7) The first embodiment of the present invention is as hereunder mentioned.

(8) The fire extinguishing agent concentration analyzer of the present invention is mainly applied in the effectiveness verification test of the fire extinguishing system in the engine compartment of an aircraft, and is used to measure whether the fire extinguishing agent in the cabin can reach the fire extinguishing concentration and maintain it for a period of time, so as to extinguish the fire and prevent the fire from reignition. The analyzer adopts a modified model with related to the concentration of fire extinguishing agent and the differential pressure to replace the original model correspondingly. Then it is obtained for the measured fire extinguishing agent concentration according to the measured differential pressure and fire extinguishing agent pressure. It is solved for the traditional problem that the noise influence caused by the flight environment pressure cannot be overcome.

(9) The method for correcting differential pressure values of the present invention uses a model to get differential pressure values. The model is with related to the concentration of fire extinguishing agent and the differential pressure. The present invention is generally directed to the form of software applied in the computer terminal. The software can be designed as an independent APP, or embedded software and so on. The computer terminal can be the controller of the fire extinguishing agent concentration analyzer of the invention, namely the CPU. The computer terminal can includes a memory, a processor, and a computer program. The computer program is stored on the memory and can be run on the processor. The computer terminal also can be smart phones, tablets, laptops, etc. capable of executing programs. In some embodiments the processor may be a Central Processing Unit (CPU), a controller, a microcontroller, a microprocessor, or other data processing chip. This processor is typically used to control the overall operation of computer equipment. In the present embodiment, the processors are used to run program code stored in memory or to process data. When the program is executed by the processor, the function of the method can be realized.

(10) In application, the method of the present invention can also be designed as a readable storage medium, such as the U-shield. The computer program instructions are stored in the readable storage medium. When the computer program instructions are read and run by a processor, the method of the analyzer is executed. In the form of U-shield, it can be plugged into the traditional extinguishing agent concentration analyzer by electronic plugging. The analyzer in the aircraft engine compartment can accurately measure the extinguishing agent concentration when the aircraft is flying by reading and executing the computer program instructions in the U-shield, and solve the traditional technical problem that the noise influence caused by the flight environment pressure cannot be overcome. Therefore, the upgrade and transformation of the traditional aircraft engine compartment fire extinguishing agent concentration analyzer is realized by means of U-shield, which is conducive to the promotion and application of the invention.

(11) Whether it is non-embedded or embedded, it can be summarized as a correction device for the differential pressure value of the corresponding extinguishing agent concentration analyzer. The device includes a first obtaining unit, a correcting unit, a measuring unit, a second obtaining unit, a third obtaining unit and a setting unit.

(12) The model with related to the concentration of fire extinguishing agent and the differential pressure is obtained by a method for correcting differential pressure values. Referring to FIG. 1, The method includes following steps: step 1 to step 6.

(13) Step 1, obtain a series of flow coefficients C.sub.d according to the mass flow-rates Q.sub.m of the fire extinguishing agent at the inlet of the analyzer under different pressures.

(14) $C_d=\frac{Q_m}{Q_{m0}}$

(15) Q.sub.m0 means the theoretical mass flow-rates of the analyzer.

(16) The step 1 can be executed by the first obtaining unit. A graph of relation with flow coefficients C.sub.d under different concentrations and corresponding pressures is shown in FIG. 2.

(17) The theoretical mass flow-rates Q.sub.m0 can be calculated by a theoretical model combined with the internal geometry of the extinguishing agent concentration analyzer, gas concentration, air pressure at the inlet, air temperature and related gas property parameters. The theoretical model is with related to the concentration of fire extinguishing agent and the differential pressure. The there are many ways to obtain Q.sub.m0. In the present embodiment, Q.sub.m0 is shown as following formula.

(18) 0$Q_{m0}=\frac{d_0^2}{4}\frac{P_0}{\sqrt{T_0}}\sqrt{\frac{M}{R_C}}\frac{1}{{\left(1+\frac{-1}{2}\right)}^{\frac{\left(+1\right)}{2\left(-1\right)}}}$

(19) In the formula, P.sub.0 means the stagnation pressure, approximate to the pressure of the fire extinguishing agent at the entrance of the analyzer.

(20) T.sub.0 means the stagnation temperature, approximate to the gas flow temperature.

(21) d.sub.0 means the aperture of the internal current limiting element of the analyzer.

(22) means the theoretical gas specific heat ratio.

(23) M means the theoretical molecular mass.

(24) R.sub.c means an universal gas constant.

(25) Step 2, correct theoretical differential pressure values P.sub.d0 according to C.sub.d, and obtaining corrected theoretical differential pressure values P.sub.d0 accordingly.
P.sub.d0=C.sub.d.Math.P.sub.d0

(26) The step 2 can be executed by the correcting unit.

(27) In the present embodiment, P.sub.d0 is shown as following formula.

(28) $P_{d0}=\frac{32d_0^{2}l}{n{d}^4}.Math.\sqrt{T_0}.Math.\frac{\sqrt{R_C/M}}{{\left(1+\frac{-1}{2}\right)}^{\frac{\left(+1\right)}{2\left(-1\right)}}}$

(29) In the formula, l means the length of the internal capillary of the analyzer.

(30) n means the total number of the capillary.

(31) d means the inner diameter of the capillary.

(32) means the theoretical dynamic viscosity of the fire extinguishing agent.

(33) Step 3, measure actual differential pressure values P.sub.d of the analyzer under different pressure. The step 3 can be executed by the measuring unit. It is only necessary to receive the corresponding measurement results from the extinguishing agent concentration analyzer at the time of execution.

(34) Step 4, obtain a series of corrected viscosity coefficients C.sub. according to P.sub.d and P.sub.d0.
P.sub.d=C.sub..Math.P.sub.d0

(35) The step 4 can be executed by the second obtaining unit. In the present embodiment, the second obtaining unit can be a multiplier.

(36) Step 5, obtain a relation C.sub.(x,P) by fitting different C.sub. under different pressures and concentrations. C.sub.(x,P) of C.sub., is with regard to the pressure P and the volumetric concentration x of the fire extinguishing agent.

(37) The step 5 can be executed by the third obtaining unit.

(38) There are many ways for fitting C.sub.(x,P). Then various forms of curve function can be adopted. In the present embodiment, linear fitting is adopted. The method for fitting C.sub.(x,P) linearly includes following steps.

(39) (1) Fitting C.sub. linearly by relation curves with regard to C.sub. and P of the fire extinguishing agent under different x.
C.sub.=a.Math.P+b

(40) In the formula, a means the slope of C.sub. and b means the intercept of C.sub..

(41) (2) fit a(x) by relation curves with many a under different concentrations, and fit b(x) by relation curves with many b under different concentrations.

(42) (3) obtain C.sub.(x,P) by different C.sub. under different a(x) and b(x).
C.sub.(x,P)=a(x).Math.P+b(x)

(43) Step 6, set a corrected model P.sub.d(x,P) with regard to P and x to replace the original model.

(44) $P_d\left(x,P\right)=C_d{C}_{}\left(x,P\right)\frac{32d_0^{2}l}{{\mathrm{nd}}^4}.Math.\sqrt{T_0}.Math.\frac{\left(x\right)\sqrt{\left(x\right)R_C/M\left(x\right)}}{{\left(1+\frac{\left(x\right)-1}{2}\right)}^{\frac{\left(\left(x\right)+1\right)}{2\left(\left(x\right)-1\right)}}}$

(45) In the formula, (x) means the dynamic viscosity of the fire extinguishing agent about the volumetric concentration x.

(46) (x) means the theoretical gas specific heat ratio about the volumetric concentration x.

(47) M(x) means the theoretical molecular mass about the volumetric concentration x.

(48) The step 6 can be executed by the setting unit.

(49) Following formulas are established if the pressure influence is ignored.

(50) $\left(x\right)=\frac{{\mathrm{xc}}_{p1}+\left(1-x\right)c_{p2}}{{\mathrm{xc}}_{p1}+\left(1-x\right)c_{p2}-R}$

(51) $M\left(x\right)={\mathrm{xM}}_1+\left(1-x\right)M_2\left(x\right)=\frac{x{}_1}{x+\left(1-x\right){}_{12}}+\frac{\left(1-x\right){}_2}{\left(1-x\right)+x{}_{21}}{}_{12}=\frac{{\left[1+{\left(\frac{{}_1}{{}_2}\right)}^{\frac{1}{2}}{\left(\frac{M_2}{M_1}\right)}^{\frac{1}{4}}\right]}^2}{{\left[8\left(1+\frac{M_2}{M_1}\right)\right]}^{\frac{1}{2}}}{}_{21}=\frac{{}_1}{{}_2}\frac{M_2}{M_1}{}_{12}$

(52) In the formula, c.sub.p1 means the specific heat at constant pressure of the fire extinguishing agent with pure component. c.sub.p2 means the specific heat at constant pressure of a carrier gas. R means the gas constant, R.sub.c/M. M.sub.1 means the molar mass of the fire extinguishing agent with pure component. M.sub.2 means the molar mass of the carrier gas. .sub.1 means the dynamic viscosity of the fire extinguishing agent with pure component. .sub.2 means the dynamic viscosity of the carrier gas.

(53) Aiming at the problem that the analyzer based on the pressure difference principle is limited to normal pressure calibration and testing. The present invention discloses a method for correcting the differential pressure value of the analyzer, which mainly includes preparing a series of measured fire extinguishing agents with the same temperature and different pressures, and feeding them into the analyzer respectively. The mass flow rate of the air at the inlet was measured against the differential pressure generated by the laminar flow structure of the fire extinguishing agent concentration analyzer, and compared with the theoretical mass flow rate and theoretical differential pressure value obtained based on the theoretical model of fire extinguishing agent concentration-differential pressure. The flow coefficient and laminar viscosity correction coefficient were introduced to explain the differences between the theory and the experiment, and finally The corrected model of the coupled pressure was established. Fixed the effect of low pressure environment on concentration test.

(54) The second embodiment of the present invention is as hereunder mentioned.

(55) The second embodiment is basically the same as the first embodiment, except that all the steps before and after the corrected model establishment are introduced in detail in the second embodiment. Please referring to FIG. 3, the method for establishing the corrected model includes following steps.

(56) Step S1, prepare a series of fire extinguishing agents having different volume concentration under the same temperature and different pressures P.sub.i (i=1, 2, 3, . . . ). The analyzer uses a standard gas prepared by the gravimetric method. The standard gas is a binary-component gas mixed by the carrier gas and the fire extinguishing agent with pure component. The carrier gas is air or nitrogen.

(57) Step S2, fed the fire extinguishing agents into the analyzer respectively under the same temperature.

(58) Step S3, measure the pressure and the temperature of the fire extinguishing agent at the entrance of the analyzer.

(59) Step S4, obtain Q.sub.m0 and P.sub.d0 combined with the internal geometry of the extinguishing agent concentration analyzer, gas concentration, air pressure at the inlet, air temperature and related gas property parameters, according to the theoretical model.

(60) $Q_{m0}=\frac{d_0^2}{4}\frac{P_0}{\sqrt{T_0}}\sqrt{\frac{M}{R_C}}\frac{1}{{\left(1+\frac{-1}{2}\right)}^{\frac{\left(+1\right)}{2\left(-1\right)}}}{P}_{d0}=\frac{32d_0^{2}l}{n{d}^4}.Math.\sqrt{T_0}.Math.\frac{\sqrt{R_C/M}}{{\left(1+\frac{-1}{2}\right)}^{\frac{\left(+1\right)}{2\left(-1\right)}}}$

(61) In the formula, d.sub.0 means the aperture of the internal current limiting element of the analyzer. P.sub.0 means the stagnation pressure. In the present embodiment, P.sub.0 means the pressure of the fire extinguishing agent at the entrance of the analyzer. T.sub.0 means the stagnation temperature. In the present embodiment, T.sub.0 means the gas flow temperature. means the theoretical gas specific heat ratio. M means the theoretical molecular mass. R.sub.c means an universal gas constant. l means the length of the internal capillary of the analyzer. n means the total number of the capillary. means the dynamic viscosity of the fire extinguishing agent.

(62) The parameters of gas properties are related to gas concentration, and the influence of pressure can be ignored. Then the parameters includes (x), M(x), (x).

(63) $\left(x\right)=\frac{{\mathrm{xc}}_{p1}+\left(1-x\right)c_{p2}}{{\mathrm{xc}}_{p1}+\left(1-x\right)c_{p2}-R}$

(64) $M\left(x\right)={\mathrm{xM}}_1+\left(1-x\right)M_2\left(x\right)=\frac{x{}_1}{x+\left(1-x\right){}_{12}}+\frac{\left(1-x\right){}_2}{\left(1-x\right)+x{}_{21}}{}_{12}=\frac{{\left[1+{\left(\frac{{}_1}{{}_2}\right)}^{\frac{1}{2}}{\left(\frac{M_2}{M_1}\right)}^{\frac{1}{4}}\right]}^2}{{\left[8\left(1+\frac{M_2}{M_1}\right)\right]}^{\frac{1}{2}}}{}_{21}=\frac{{}_1}{{}_2}\frac{M_2}{M_1}{}_{12}$

(65) In the formula, R means the gas constant, R.sub.c/M. x means the volumetric concentration of the fire extinguishing agent. c.sub.p1 means the specific heat at constant pressure of the fire extinguishing agent with pure component. c.sub.p2 means the specific heat at constant pressure of a carrier gas. M.sub.1 means the molar mass of the fire extinguishing agent with pure component. M.sub.2 means the molar mass of the carrier gas. .sub.1 means the dynamic viscosity of the fire extinguishing agent with pure component. .sub.2 means the dynamic viscosity of the carrier gas. The related gas properties of the dynamic viscosity and the carrier gas, such as specific heat capacity at constant pressure, molar mass, dynamic viscosity, etc., can be checked by some physical property parameter softwares.

(66) Step S5, obtain the mass flow-rates Q.sub.m of the fire extinguishing agent at the inlet of the analyzer under different pressures. Compared the mass flow-rates Q.sub.m with the theoretical mass flow-rates Q.sub.m0 to get a series of flow coefficients C.sub.d.

(67) $C_d=\frac{Q_m}{Q_{m0}}$

(68) Step S6, The pressure and gas concentration have little effect on the flow coefficient C.sub.d, which is mainly related to the internal current limiting structure of the fire extinguishing agent concentration analyzer, so the average value of the series of flow coefficients C.sub.d is taken. A corrected theoretical differential pressure values P.sub.d0 is obtained according to the theoretical differential pressure values P.sub.d0.
P.sub.d0=C.sub.d.Math.P.sub.d0

(69) Step S7, measure actual differential pressure values P.sub.d of the analyzer under different pressure.

(70) Step S8, obtain a series of corrected viscosity coefficients C.sub. according to P.sub.d and P.sub.d0.

(71) $P_d=C_{}.Math.P_{d0}^{}=C_d{C}_{}\frac{32d_0^{2}l}{{\mathrm{nd}}^4}.Math.\sqrt{T_0}.Math.\frac{\sqrt{R_C/M}}{{\left(1+\frac{-1}{2}\right)}^{\frac{\left(+1\right)}{2\left(-1\right)}}}$

(72) Step S9, obtain a relation C.sub.(x,P) by fitting different C.sub. under different pressures and concentrations, C.sub.(x,P) of C.sub. with regard to the pressure P and the volumetric concentration x of the fire extinguishing agent.

(73) The method for fitting C.sub.(x,P) linearly includes following steps.

(74) Step 9.1, fit C.sub. linearly by relation curves with regard to C.sub. and P of the fire extinguishing agent under different x.
C.sub.=a.Math.P+b

(75) In the formula, C.sub. means corrected viscosity coefficients. a means the slope of C.sub.. b means the intercept of C.sub..

(76) Step 9.2, fit a(x) by relation curves with many a under different concentrations, fit b(x) by relation curves with many b under different concentrations.

(77) Step 9.3, obtain C.sub.(x,P) by different C.sub. under different a(x) and b(x) accordingly.
C.sub.(x,P)=a(x).Math.P+b(x)

(78) Step 10, set a corrected model P.sub.d(x,P) with regard to P and x to replace the original model.

(79) 0$P_d\left(x,P\right)=C_d{C}_{}\left(x,P\right)\frac{32d_0^{2}l}{{\mathrm{nd}}^4}.Math.\sqrt{T_0}.Math.\frac{\left(x\right)\sqrt{\left(x\right)R_C/M\left(x\right)}}{{\left(1+\frac{\left(x\right)-1}{2}\right)}^{\frac{\left(\left(x\right)+1\right)}{2\left(\left(x\right)-1\right)}}}$

(80) Compared with the prior art, the present invention has the following technical effects. The invention takes the pressure as the hidden factor variable of separation to participate in the establishment of the concentration-differential pressure model, and can rely on the adaptability of the model to the pressure to realize the concentration test under different working pressures, so that the established model has better versatility. At the same time, the modified coefficient is introduced to describe the complex effects of kinetic energy, gas expansion, gas non-ideality, wall slip and other factors to simplify the model.

(81) The third embodiment of the present invention is as hereunder mentioned.

(82) The third embodiment is an introduction to a specific practical operation of the second embodiment.

(83) The entire flow chart of implementation steps is basically implemented in accordance with FIG. 1 of the second embodiment.

(84) Step S1, prepare a series of Pentafluoroethane gas having different volume concentration (0%, 2%, 4.02%, 9.02%, 13%, 17.99%, 25.04%, 30.03%) under the same temperature (25 C.) and different pressures (60 kPa, 70 kPa, 80 kPa, 90 kPa, 101 kPa). The gravimetric method is used to prepare the series of Pentafluoroethane gas. The carrier gas is air or nitrogen.

(85) Step S2, fed the fire extinguishing agents prepared under different pressure (60 kPa, 70 kPa, 80 kPa, 90 kPa, 101 kPa) into the analyzer respectively.

(86) Step S3, measure the pressure and the temperature of the fire extinguishing agent at the entrance of the analyzer.

(87) Step S4, obtain Q.sub.m0 and P.sub.d0 combined with the internal geometry of the extinguishing agent concentration analyzer, gas concentration, air pressure at the inlet, air temperature and related gas property parameters, according to the theoretical model.

(88) $Q_{m0}=\frac{d_0^2}{4}\frac{P_0}{\sqrt{T_0}}\sqrt{\frac{M}{R_c}}\frac{1}{{\left(1+\frac{-1}{2}\right)}^{\frac{\left(+1\right)}{2\left(-1\right)}}}$$P_{d0}=\frac{32d_0^{2}l}{{\mathrm{nd}}^4}.Math.\sqrt{T_0}.Math.\frac{\sqrt{R_c/M}}{{\left(1+\frac{-1}{2}\right)}^{\frac{\left(+1\right)}{2\left(-1\right)}}}$

(89) In the formula, d.sub.0 means the aperture of the internal current limiting element of the analyzer. P.sub.0 means the stagnation pressure. In the present embodiment, P.sub.0 means the pressure of the fire extinguishing agent at the entrance of the analyzer. T.sub.0 means the stagnation temperature. In the present embodiment, To means the gas flow temperature. means the theoretical gas specific heat ratio. M means the theoretical molecular mass. R.sub.c means an universal gas constant. l means the length of the internal capillary of the analyzer. n means the total number of the capillary. d means the inner diameter of the capillary. means the dynamic viscosity of the fire extinguishing agent.

(90) The parameters of gas properties are related to gas concentration, and the influence of pressure can be ignored. Then the parameters includes (x), M(x), (x).

(91) $\left(x\right)=\frac{{\mathrm{xc}}_{p1}+\left(1-x\right)c_{p2}}{{\mathrm{xc}}_{p1}+\left(1-x\right)c_{p2}-R}$

(92) $M\left(x\right)={\mathrm{xM}}_1+\left(1-x\right)M_2$$\left(x\right)=\frac{x{}_1}{x+\left(1-x\right){}_{12}}+\frac{\left(1-x\right){}_2}{\left(1-x\right)+x{}_{21}}$${}_{12}=\frac{{\left[1+{\left(\frac{{}_1}{{}_2}\right)}^{\frac{1}{2}}{\left(\frac{M_2}{M_1}\right)}^{\frac{1}{4}}\right]}^2}{{\left[8\left(1+\frac{M_1}{M_2}\right)\right]}^{\frac{1}{2}}}$${}_{21}=\frac{{}_2}{{}_1}\frac{M_1}{M_2}{}_{12}$

(93) In the formula, R means the gas constant, R.sub.c/M. x means the volumetric concentration of the fire extinguishing agent. c.sub.p1 means the specific heat at constant pressure of the fire extinguishing agent with pure component. c.sub.p2 means the specific heat at constant pressure of a carrier gas. M.sub.1 means the molar mass of the fire extinguishing agent with pure component. M.sub.2 means the molar mass of the carrier gas. .sub.1 means the dynamic viscosity of the fire extinguishing agent with pure component. .sub.2 means the dynamic viscosity of the carrier gas.

(94) Specifically, the specific heat capacity at constant pressure, molar mass and dynamic viscosity of pure component fire extinguishing agent and carrier gas at 25 C. were obtained by physical property parameter software, as shown in Table 1.

(95) TABLE-US-00001 TABLE 1 Molar specific heat at dynamic mass M constant pressure viscosity Gas type (g/mol) c.sub.p (J/(kg .Math. K) (Pa .Math. s) pentafluoroethane 120.0214 796.51 1.296 10.sup.5 air 28.96 1006.6 1.849 10.sup.5

(96) Step S5, obtain the mass flow-rates Q.sub.m of the fire extinguishing agent at the inlet of the analyzer under different pressures. Compared the mass flow-rates Q.sub.m with the theoretical mass flow-rates Q.sub.m0 to get a series of flow coefficients C.sub.d.

(97) $C_d=\frac{Q_m}{Q_{m0}}$

(98) Step S6, The pressure and gas concentration have little effect on the flow coefficient C.sub.d, which is mainly related to the internal current limiting structure of the fire extinguishing agent concentration analyzer, so the average value 0.86 of the series of flow coefficients C.sub.d is taken. A corrected theoretical differential pressure values P.sub.d0 is obtained according to the theoretical differential pressure values P.sub.d0.
P.sub.d0=C.sub.d.Math.P.sub.d0

(99) Step S7, measure actual differential pressure values P.sub.d of the analyzer under different pressure.

(100) Step S8, referring to FIG. 2 obtain a series of corrected viscosity coefficients C.sub. according to P.sub.d and P.sub.d0.

(101) $P_d=C_{}.Math.P_{d0}^{}=C_d{C}_{}\frac{32d_0^{2}l}{{\mathrm{nd}}^4}.Math.\sqrt{T_0}.Math.\frac{\sqrt{R_c/M}}{{\left(1+\frac{-1}{2}\right)}^{\frac{\left(+1\right)}{2\left(-1\right)}}}$

(102) Step S9, obtain a relation C.sub.(x,P) by fitting different C.sub. under different pressures and concentrations, C.sub.(x,P) of C.sub. with regard to the pressure P and the volumetric concentration x of the fire extinguishing agent.

(103) The method for fitting C.sub.(x,P) linearly includes following steps.

(104) Step 9.1, fit C.sub. linearly by relation curves with regard to C.sub. and P of the fire extinguishing agent under different x.
C.sub.=a.Math.P+b

(105) In the formula, C.sub. means corrected viscosity coefficients. a means the slope of C.sub.. b means the intercept of C.sub..

(106) Referring to table 2, slope, intercept values and fit R.sup.2 under different pressures and different concentrations are shown.

(107) TABLE-US-00002 TABLE 2 concentration slope a.sub.i intercept b.sub.i R.sup.2 0% 4.03 10.sup.4 0.6672 0.99673 2.00% 3.90 10.sup.4 0.67339 0.99515 4.02% 4.04 10.sup.4 0.68028 0.99740 9.02% 4.42 10.sup.4 0.6834 0.99910 13.00% 4.05 10.sup.4 0.69214 0.99936 17.99% 4.63 10.sup.4 0.69823 0.99921 25.04% 4.46 10.sup.4 0.70757 0.99989 30.03% 4.66 10.sup.4 0.71286 0.99944

(108) Through the above analysis, the concentration a has little effect on the slope, so the slope a is averaged as 4.2710.sup.4. Then at a fixed slope, linear fitting was performed on the relationship between the above viscosity correction coefficient and pressure again, and the fitting parameters obtained were shown in Table 3.

(109) TABLE-US-00003 TABLE 3 concentration slope a.sub.i intercept b.sub.i R.sup.2 0% 4.27 10.sup.4 0.66525 0.99310 2.00% 0.67044 0.98631 4.02% 0.67841 0.99405 9.02% 0.6846 0.99794 13.00% 0.69037 0.99638 17.99% 0.70114 0.99309 25.04% 0.70907 0.99813 30.03% 0.716 0.99240

(110) Step 9.2, fit a(x) by relation curves with many a under different concentrations, fit b(x) by relation curves with many b under different concentrations.
b(x)=0.66664+0.0022 x1.91342 x.sup.2

(111) Step 9.3, obtain C.sub.(x,P) by different C.sub. under different a(x) and b(x) accordingly.
C.sub.(x,P)=4.2710.sup.4.Math.P+0.66664+0,0022x1.91342 x.sup.2

(112) Step 10, set a corrected model P.sub.d(x,P) with regard to P and x to replace the original model.

(113) $P_d\left(x,P\right)=C_d{C}_{}\left(x,P\right)\frac{32d_0^{2}l}{{\mathrm{nd}}^4}.Math.\sqrt{T_0}.Math.\frac{\left(x\right)\sqrt{\left(x\right)R_c/M\left(x\right)}}{{\left(1+\frac{\left(x\right)-1}{2}\right)}^{\frac{\left(\left(x\right)+1\right)}{2\left(\left(x\right)-1\right)}}}$

(114) As the current calibration and testing of fire extinguishing agent concentration analyzer are limited to atmospheric pressure environment, and the impact of flight low pressure environment on concentration testing is unknown. The present invention adopts a relatively simplified method and introduces a correction factor to uniformly describe the complex effects brought by kinetic energy, gas expansion, gas non-ideal, wall slip and other factors. The concentration-differential pressure model of coupled pressure was constructed for the fire extinguishing agent concentration analyzer based on the principle of pressure difference. The influence of pressure on concentration test was corrected. The adaptability of the analyzer to the flight pressure variable environment was improved. The performance of the analyzer was optimized.

(115) While the invention has been described in terms of what is presently considered to be the most practical and preferred embodiments, it is to be understood that the invention needs not to be limited to the disclosed embodiment. On the contrary, it is intended to cover various modifications and similar arrangements included within the spirit and scope of the appended claims which are to be accorded with the broadest interpretation so as to encompass all such modifications and similar structures.