CN111027155A - Simulation analysis method of airplane air circulation refrigeration system based on CFD technology - Google Patents

Simulation analysis method of airplane air circulation refrigeration system based on CFD technology Download PDF

Info

Publication number
CN111027155A
CN111027155A CN201911293882.4A CN201911293882A CN111027155A CN 111027155 A CN111027155 A CN 111027155A CN 201911293882 A CN201911293882 A CN 201911293882A CN 111027155 A CN111027155 A CN 111027155A
Authority
CN
China
Prior art keywords
refrigeration system
formula
equation
air circulation
gas
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201911293882.4A
Other languages
Chinese (zh)
Inventor
赵艺博
肖天雷
茅晓晨
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Institute of Measurement and Testing Technology
Original Assignee
Shanghai Institute of Measurement and Testing Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai Institute of Measurement and Testing Technology filed Critical Shanghai Institute of Measurement and Testing Technology
Priority to CN201911293882.4A priority Critical patent/CN111027155A/en
Publication of CN111027155A publication Critical patent/CN111027155A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a simulation analysis method of an aircraft air circulation refrigeration system based on a CFD technology, which comprises the following steps: establishing a geometric model of an air circulation refrigeration system of the airplane; dividing and setting a control area; setting a boundary condition; dividing grids; establishing a control equation; carrying out numerical solution; and analyzing the simulation result and predicting the refrigeration performance of the air circulation refrigeration system of the airplane. The two-dimensional calculation simulation of the air circulation refrigeration system of the airplane is realized through the commercial calculation fluid software Fluent, the refrigeration performance of the refrigeration system can be predicted in a short time, and the research and development cost is greatly reduced; and the operation mechanism and the fluid characteristics in the refrigeration system are deeply understood by analyzing the distribution of a temperature field, a pressure field and a flow field in the refrigeration system, and the method has important significance for the design optimization of the air circulation refrigeration system of the airplane.

Description

Simulation analysis method of airplane air circulation refrigeration system based on CFD technology
Technical Field
The invention relates to the technical field of refrigeration, in particular to a simulation analysis method of an aircraft air circulation refrigeration system based on a CFD (computational fluid dynamics) technology.
Background
The air conditioning system on an aircraft is divided into four subsystems: a refrigeration system, a temperature control system, an air distribution system, and a pressurization control system. The most important of them is the refrigeration system, and the core of the refrigeration system is the air circulation refrigerator. The two companies that are currently most competitive in aircraft manufacture and development in the world are the united states boeing company and the european airbus company, respectively. Table 1 summarizes the air cycle chiller versions employed by boeing and airbus aircraft. The air circulation refrigerator of the airplane which is most widely applied at present adopts a three-wheel or four-wheel boosting type high-pressure water removal circulation system, has stronger ground refrigerating capacity, and keeps the advantages of small air supply pressure and power saving. The latest Boeing 787 aircraft air circulation refrigerator adopts an electrically-driven four-wheel boosting high-pressure water removal refrigeration system, engine air entraining is omitted, and energy consumption of civil aircraft is greatly saved, so that fuel consumption of the engine is reduced.
TABLE 1 air circulation refrigerator form for the cockpit of a civil aircraft in the world
Figure BDA0002319948720000011
For China, the air circulation refrigerator can only adopt two-stage bleed air and three-wheel boosting high-pressure water removal systems of an engine according to the development level of the current environment control accessories in China. As the first civil aircraft ARJ21-700 with independent intellectual property rights in China, the air refrigerating machine of the civil aircraft adopts the circulating refrigerating system. For the main line civil aircraft C919 with independent intellectual property rights according to the latest international airworthiness standard in the first money of China, the air entraining system still adopts an engine air entraining mode, and a certain gap still exists compared with the most advanced technology in foreign countries.
Disclosure of Invention
The invention overcomes the defects of the prior art and provides a simulation analysis method and a simulation analysis system of an aircraft air circulation refrigeration system based on a CFD technology.
The method adopts commercial calculation fluid software Fluent to establish a two-dimensional calculation model of the air circulation refrigeration system of the airplane, performs numerical simulation and predicts the refrigeration performance of the refrigeration system. Furthermore, the distribution of a temperature field, a pressure field and a flow field in the refrigerating system is analyzed, and theoretical and practical application basis is provided for the design of a practical and efficient refrigerating system.
The invention relates to a simulation analysis method of an aircraft air circulation refrigeration system based on a CFD technology, which comprises the following steps:
the method comprises the following steps: establishing a geometric model of an air circulation refrigeration system of the airplane;
step two: dividing and setting a control area;
step three: setting a boundary condition;
step four: dividing grids;
step five: establishing an unsteady state control equation, and respectively establishing a control equation for a gas area and a porous medium area by taking working medium gas air in a refrigeration system as an object;
step six: carrying out numerical solution;
step seven: and analyzing the simulation result and predicting the refrigeration performance of the air circulation refrigeration system of the airplane.
In the invention, further, the distribution of the temperature field, the pressure field and the flow field in the aircraft air circulation refrigeration system is observed, so that the operation mechanism of the refrigeration system and the flow state of the fluid can be analyzed.
The simulation analysis method comprises the following specific steps:
the first step is as follows: a geometric model of the air circulation refrigeration system of the airplane is established, and the model is a two-dimensional CFD model and comprises a precooler 1, a regulating valve 2, a primary heat exchanger 3, a compressor 4, a secondary heat exchanger 5, an expansion machine 6, a water separator 7, a one-way valve 8 and a mixing chamber 9.
The second step is as follows: the control regions are divided and set, and the internal control region where the precooler 1, the primary heat exchanger 3, and the secondary heat exchanger 5 are located is set as a porous medium region, the internal control region where the water separator 7 is located is set as a gas-liquid mixing region, and the other internal control regions are set as gas regions.
The third step is that: setting boundary conditions, setting the left end boundary of the precooler 1 as a speed inlet boundary, and setting the outer boundaries of the regulating valve 2, the compressor 4, the expander 6, the water separator 7, the one-way valve 8 and the mixing chamber 9 as isothermal boundaries; the outer boundaries of the primary heat exchanger 3 and the secondary heat exchanger 5 are set as adiabatic boundaries.
The fourth step is that: and dividing grids, adopting structured grids for all control domains of the model, and simultaneously carrying out local grid refinement near the boundary.
The fifth step is as follows: and establishing an unsteady state control equation, and establishing a control equation for a gas area and a porous medium area respectively by taking working medium gas air in the refrigerating system as an object.
First, a governing equation is established for the gas region:
continuity equation:
Figure BDA0002319948720000021
in the formula (1), ρfAnd
Figure BDA0002319948720000031
respectively the density and velocity of the gas.
The momentum equation:
Figure BDA0002319948720000032
in the formula (2), p is a static pressure,
Figure BDA0002319948720000033
is the pressure tensor;
energy equation:
Figure BDA0002319948720000034
in the formula (3), kfIs gas thermal conductivity, T is temperature, EfThe total energy of the gas;
then, a control equation is established for the porous medium region:
continuity equation:
Figure BDA0002319948720000035
the momentum equation:
Figure BDA0002319948720000036
in the formula (4), ε represents the porosity of the porous medium, and the porous medium region is filled with the screen, so that the porosity ε is determined by the mesh number m and the diameter d of the screenwDetermining:
Figure BDA0002319948720000037
hydraulic diameter d of the wire meshhCan be calculated by equation (7):
Figure BDA0002319948720000038
in the formula (5), SiIs a source term in the momentum equation and consists of two parts of viscous loss term and inertial loss term, and for isotropic media, SiCan be expressed as:
Figure BDA0002319948720000039
in the formula (8), μ is the dynamic viscosity of the gas, α is the degree of permeability, C2Coefficient of inertial resistance, α and C2Can be calculated by the following method:
the axial pressure drop of a region of porous media can be expressed as:
Figure BDA00023199487200000310
in the formula (9), foscFor the friction factor, according to empirical formula:
Figure BDA00023199487200000311
in equation (10), Re is the reynolds number of the fluid in the porous medium region and can be expressed as:
Figure BDA00023199487200000312
by combining the formula (9), the formula (10) and the formula (11), the degree of penetration α and the coefficient of inertial resistance C can be determined2
Figure BDA0002319948720000041
Figure BDA0002319948720000042
Energy equation:
Figure BDA0002319948720000043
in the formula (14), psAnd EsDensity and total energy, k, respectively, of the solideffEffective thermal conductivity for porous media;
effective thermal conductivity k of porous media due to the effect of thermal contact resistanceeffMuch less than its average thermal conductivity, the following modifications can be made:
Figure BDA0002319948720000044
in formula (15), ksAnd kfThe thermal conductivity of the solid and gas respectively.
The sixth step: and (4) carrying out numerical solution, setting solver control parameters, a discrete format and a residual convergence standard according to the CFD model of the aircraft air circulation refrigeration system established in the steps, initializing a flow field, and starting to carry out numerical solution.
The seventh step is as follows: and analyzing the simulation result, and predicting the refrigeration performance of the air circulation refrigeration system of the airplane through analyzing the simulation result after the calculation is stable.
In addition, further, by observing the distribution conditions of the temperature field, the pressure field and the flow field inside the refrigeration system, the operation mechanism of the refrigeration system and the flow state of the fluid can be further analyzed.
The invention also provides a simulation analysis system of the aircraft air circulation refrigeration system based on the CFD technology, which comprises: the device comprises a temperature simulation analysis module, a pressure simulation analysis module and a speed simulation analysis module. The temperature simulation analysis module is used for simulating and analyzing a temperature field in the model, the pressure simulation analysis module is used for analyzing a pressure field in the model, and the speed simulation analysis module is used for simulating and analyzing a flow field in the model.
The model establishing module is used for establishing a geometric model of the aircraft air circulation refrigeration system;
the area setting module is used for dividing and setting a control area;
the boundary setting module is used for setting boundary conditions;
the grid division module is used for dividing grids;
the equation establishing module is used for establishing an unsteady state control equation, and respectively establishing control equations for a gas area and a porous medium area by taking working medium gas air in the refrigerating system as an object;
the numerical solving module is used for carrying out numerical solving;
and the result analysis module is used for analyzing the simulation result and predicting the refrigeration performance of the air circulation refrigeration system of the airplane.
Wherein the result analysis module comprises: the device comprises a temperature simulation analysis module, a pressure simulation analysis module and a speed simulation analysis module; the temperature simulation analysis module is used for simulating a temperature field in the analysis model, the pressure simulation analysis module is used for analyzing a pressure field in the model, and the speed simulation analysis module is used for simulating a flow field in the analysis model.
The simulation analysis system is used for realizing the simulation analysis method of the airplane air circulation refrigeration system based on the CFD technology.
The invention has the following beneficial effects and advantages: by the simulation analysis method and the simulation analysis system, the refrigeration performance of the air circulation refrigeration system of the airplane can be predicted in a short time, and a large amount of test time is saved; by the simulation analysis method and the system, a series of simulation results such as a temperature field, a pressure field, a flow field and the like in the air circulation refrigeration system of the airplane can be visually displayed, and the operation mechanism and the fluid characteristics in the refrigeration system can be deeply understood. The invention greatly reduces the research and development cost of the air circulation refrigeration system of the airplane and shortens the research and development period. The method and the device realize the rapid prediction of the refrigeration performance of the aircraft air circulation refrigeration system, deeply understand the operation mechanism and the fluid characteristic in the refrigeration system, and have important significance for the design optimization of the aircraft air circulation refrigeration system.
Drawings
Fig. 1 is a schematic view of a geometric model of an aircraft air cycle refrigeration system according to the present invention, wherein: 1-precooler, 2-regulating valve, 3-primary heat exchanger, 4-compressor, 5-secondary heat exchanger, 6-expander, 7-water separator, 8-one-way valve and 9-mixing chamber.
FIG. 2 is a schematic flow chart of the simulation analysis method of the present invention.
FIG. 3 is a schematic flow chart of a simulation analysis method according to the present invention.
FIG. 4 is a diagram of a result analysis module according to the present invention.
Fig. 5 is a temperature field profile inside a heat exchanger.
Fig. 6 is a pressure field profile inside a heat exchanger.
Fig. 7 is a flow field profile inside a heat exchanger.
FIG. 8 is a schematic diagram of a simulation analysis system according to the present invention.
Detailed Description
The present invention will be described in further detail with reference to the following specific examples and the accompanying drawings. The procedures, conditions, experimental methods and the like for carrying out the present invention are general knowledge and common general knowledge in the art except for the contents specifically mentioned below, and the present invention is not particularly limited. The following detailed description of the embodiments of the present invention is provided in connection with the accompanying drawings and the examples, but not limited thereto, and it is to be understood that the present invention is not limited to the modifications and equivalents, and that the modifications and equivalents may be made without departing from the spirit and scope of the invention.
In order to overcome the defects of the prior art, the invention provides a simulation analysis method and a simulation analysis system of an aircraft air circulation refrigeration system based on a CFD technology, wherein the method comprises the following steps: establishing a geometric model of an air circulation refrigeration system of the airplane; dividing and setting a control area; setting a boundary condition; dividing grids; establishing a control equation; carrying out numerical solution; and analyzing the simulation result and predicting the refrigeration performance of the air circulation refrigeration system of the airplane. The two-dimensional calculation simulation of the air circulation refrigeration system of the airplane is realized through the commercial calculation fluid software Fluent, the refrigeration performance of the refrigeration system can be predicted in a short time, and the research and development cost is greatly reduced; and the operation mechanism and the fluid characteristics in the refrigeration system are deeply understood by analyzing the distribution of a temperature field, a pressure field and a flow field in the refrigeration system, and the method has important significance for the design optimization of the air circulation refrigeration system of the airplane.
The method adopts commercial calculation fluid software Fluent to establish a two-dimensional calculation model of the air circulation refrigeration system of the airplane, performs numerical simulation, predicts the refrigeration performance of the refrigeration system, analyzes the distribution of a temperature field, a pressure field and a flow field in the refrigeration system, and provides a theoretical basis for the design of the practical and efficient refrigeration system.
As shown in fig. 3, the simulation analysis method of the present invention includes the following steps:
the method comprises the following steps: establishing a geometric model of an air circulation refrigeration system of the airplane;
step two: dividing and setting a control area;
step three: setting a boundary condition;
step four: dividing grids;
step five: establishing an unsteady state control equation, and respectively establishing a control equation for a gas area and a porous medium area by taking working medium gas air in a refrigeration system as an object;
step six: carrying out numerical solution;
step seven: analyzing the simulation result and predicting the refrigeration performance of the air circulation refrigeration system of the airplane
In this embodiment, the object of the simulation analysis is a three-wheel boost air cycle refrigeration system.
Specifically, as shown in fig. 2, the simulation analysis method of the present invention includes the following steps:
the method comprises the following steps: the geometric model of the three-wheel boosting type air circulation refrigeration system is established, and is a two-dimensional CFD model and comprises a precooler 1, a regulating valve 2, a primary heat exchanger 3, a compressor 4, a secondary heat exchanger 5, an expander 6, a water separator 7, a one-way valve 8 and a mixing chamber 9, and is shown in figure 1.
Step two: dividing and setting control areas, setting the internal control areas where the precooler 1, the primary heat exchanger 3 and the secondary heat exchanger 5 are located as porous medium areas, setting the internal control area where the water separator 7 is located as a gas-liquid mixing area, and setting the other internal control areas as gas areas;
step three: the boundary condition is set, and the left end boundary of the precooler 1 is set as the velocity inlet boundary, and in the present embodiment, the velocity is set to 30 m/s. The outer boundaries of the regulating valve 2, the compressor 4, the expander 6, the water separator 7, the one-way valve 8 and the mixing chamber 9 are set as isothermal boundaries; setting the outer boundaries of the primary heat exchanger 3 and the secondary heat exchanger 5 as heat insulation boundaries;
step four: dividing grids, adopting structured grids for all control domains of the model, and simultaneously carrying out local grid refinement near the boundary;
step five: establishing an unsteady state control equation, taking working medium gas air in a refrigerating system as an object, and respectively establishing control equations for a gas area and a porous medium area:
firstly, a control equation is established for a gas area:
continuity equation:
Figure BDA0002319948720000071
in the formula: rhofAnd
Figure BDA0002319948720000072
respectively the density and velocity of the gas.
The momentum equation:
Figure BDA0002319948720000073
in the formula: p is the static pressure of the molten steel,
Figure BDA0002319948720000074
is the pressure tensor;
energy equation:
Figure BDA0002319948720000075
in the formula: k is a radical offIs gas thermal conductivity, T is temperature, EfThe total energy of the gas;
and then establishing a control equation for the porous medium region:
continuity equation:
Figure BDA0002319948720000076
the momentum equation:
Figure BDA0002319948720000077
in the formula: epsilon is the porosity of the porous medium, and the porous medium area is filled with the silk screen, so that the porosity epsilon is formed by the mesh number m and the silk diameter d of the silk screenwDetermining:
Figure BDA0002319948720000078
hydraulic diameter d of the wire meshhCan be calculated by the following formula:
Figure BDA0002319948720000079
Siis a source term in the momentum equation and consists of two parts of viscous loss term and inertial loss term, and for isotropic media, SiCan be expressed as:
Figure BDA00023199487200000710
wherein μ is the dynamic viscosity of the gas, α is the degree of permeability, C2Coefficient of inertial resistance, α and C2Can be calculated by the following method:
the axial pressure drop of a region of porous media can be expressed as:
Figure BDA00023199487200000711
in the formula: f. ofoscFor the friction factor, according to empirical formula:
Figure BDA00023199487200000712
in the formula: re is the Reynolds number of the fluid in the region of porous media and can be expressed as:
Figure BDA00023199487200000713
by combining the formula (9), the formula (10) and the formula (11), the degree of penetration α and the coefficient of inertial resistance C can be determined2
Figure BDA0002319948720000081
Figure BDA0002319948720000082
Energy equation:
Figure BDA0002319948720000083
in the formula: rhosAnd EsDensity and total energy, k, respectively, of the solideffEffective thermal conductivity for porous media;
effective thermal conductivity k of porous media due to the effect of thermal contact resistanceeffMuch less than its average thermal conductivity, the following modifications can be made:
Figure BDA0002319948720000084
in this embodiment, the porous medium region inside the precooler 1 is filled with a 635-mesh stainless steel wire mesh, and the porous medium regions inside the primary heat exchanger 3 and the secondary heat exchanger 5 are filled with a 100-mesh copper wire mesh. Table 1 summarizes the various parameters of the web obtained by the above calculations.
TABLE 1 internal Screen parameters of porous Medium zone
Region(s) Material m dw(μm) ε dh(μm) α(m2) C2(m-1)
Precooler 1 SS304 635 20 0.607 31.4 1.5×10-11 81460
Primary heat exchanger 3 Copper (Cu) 100 100 0.691 222.6 7.7×10-10 12173
Secondary heat exchanger 5 Copper (Cu) 100 100 0.691 222.6 7.7×10-10 12173
Step six: performing numerical solution, setting solver control parameters, a discrete format and a residual convergence standard according to the CFD model of the aircraft air circulation refrigeration system established in the step, initializing a flow field, and starting to perform numerical solution;
in the embodiment, a laminar flow model is selected for calculation, a solver selects a PISO algorithm of pressure and speed coupling, and discrete format sampling is performedWith a second order windward format, the residual convergence standard energy term is 10-6All other parameters are 10-3The initial temperature was 300K.
Step seven: and analyzing the simulation result, and predicting the refrigeration performance of the air circulation refrigeration system of the airplane through analyzing the simulation result after the calculation is stable.
As shown in fig. 8, the present invention further provides a simulation analysis system of an aircraft air cycle refrigeration system based on CFD technology, where the system adopts a simulation analysis method of an aircraft air cycle refrigeration system based on CFD technology, and the method includes:
the model establishing module is used for establishing a geometric model of the aircraft air circulation refrigeration system;
the area setting module is used for dividing and setting a control area;
the boundary setting module is used for setting boundary conditions;
the grid division module is used for dividing grids;
the equation establishing module is used for establishing an unsteady state control equation, and respectively establishing control equations for a gas area and a porous medium area by taking working medium gas air in the refrigerating system as an object;
the numerical solving module is used for carrying out numerical solving;
the result analysis module is used for analyzing the simulation result and predicting the refrigeration performance of the air circulation refrigeration system of the airplane; the temperature simulation analysis module is used for simulating a temperature field in an analysis model, the pressure simulation analysis module is used for analyzing a pressure field in the model, and the speed simulation analysis module is used for simulating a flow field in the analysis model.
Furthermore, by observing the distribution of the temperature field, the pressure field and the flow field inside the refrigeration system, the operation mechanism of the refrigeration system and the flow state of the fluid can be further analyzed.
As shown in fig. 5, the temperature field distribution inside the primary heat exchanger was simulated.
As shown in FIG. 6, the pressure field distribution inside the primary heat exchanger was simulated
As shown in fig. 7, the flow field distribution inside the primary heat exchanger was simulated.
The protection of the present invention is not limited to the above embodiments. Variations and advantages that may occur to those skilled in the art may be incorporated into the invention without departing from the spirit and scope of the inventive concept, which is set forth in the following claims.

Claims (10)

1. A simulation analysis method of an aircraft air circulation refrigeration system based on CFD technology is characterized by comprising the following steps:
the method comprises the following steps: establishing a geometric model of an air circulation refrigeration system of the airplane;
step two: dividing and setting a control area;
step three: setting a boundary condition;
step four: dividing grids;
step five: establishing an unsteady state control equation, and respectively establishing a control equation for a gas area and a porous medium area by taking working medium gas air in a refrigeration system as an object;
step six: carrying out numerical solution;
step seven: and analyzing the simulation result and predicting the refrigeration performance of the air circulation refrigeration system of the airplane.
2. The method as claimed in claim 1, wherein the method uses a commercial computational fluid software Fluent to build a two-dimensional computational model of an aircraft air cycle refrigeration system, to perform numerical simulations to predict the refrigeration performance of the refrigeration system.
3. The method of claim 1, further comprising analyzing an operating mechanism of the refrigeration system and a flow condition of the fluid in conjunction with observing a distribution of a temperature field, a pressure field, and a flow field within the aircraft air cycle refrigeration system.
4. The method of claim 1, wherein in the first step, the geometric model is a two-dimensional CFD model.
5. The method according to claim 1, wherein in the second step, the control area is divided and set, and comprises a porous medium area, a gas-liquid mixing area and a gas area; wherein the porous medium area is an internal control area where the precooler (1), the primary heat exchanger (3) and the secondary heat exchanger (5) are located; the gas-liquid mixing area is an internal control area where the water separator (7) is located; the gas zone is the other inner control zone.
6. The method of claim 1, wherein in step three, the boundary conditions comprise: velocity entrance boundaries, isothermal boundaries, adiabatic boundaries; wherein the speed inlet boundary is set as the left end boundary of the precooler (1); the isothermal boundary is set as the outer boundary of the regulating valve (2), the compressor (4), the expander (6), the water separator (7), the one-way valve (8) and the mixing chamber (9); the heat insulation boundary is set as the outer boundary of the primary heat exchanger (3) and the secondary heat exchanger (5).
7. The method of claim 1, wherein in the fourth step, the grid is divided, the structured grid is adopted for the control domain of the geometric model, and the local grid refinement is performed near the boundary.
8. The method of claim 1, wherein in step five,
establishing a governing equation for the gas region:
continuity equation:
Figure FDA0002319948710000011
in formula (1): rhofAnd
Figure FDA0002319948710000012
density and velocity of the gas, respectively, and t is time.
The momentum equation:
Figure FDA0002319948710000021
in formula (2): p is the static pressure of the molten steel,
Figure FDA0002319948710000022
is the pressure tensor.
Energy equation:
Figure FDA0002319948710000023
in formula (3): k is a radical offIs gas thermal conductivity, T is temperature, EfIs always the gas.
9. The method of claim 1, wherein in step five,
establishing a control equation for the porous medium region:
continuity equation:
Figure FDA0002319948710000024
in the formula (4), ε represents the porosity of the porous medium.
The momentum equation:
Figure FDA0002319948710000025
the porosity epsilon is determined by the mesh number m and the diameter d of the screenwDetermining:
Figure FDA0002319948710000026
hydraulic diameter d of the wire meshhCan be calculated by equation (7):
Figure FDA0002319948710000027
in the formula (5), SiIs a source term in the momentum equation and consists of two parts of viscous loss term and inertial loss term, and for isotropic media, SiCan be expressed as:
Figure FDA0002319948710000028
in the formula (8), μ is the dynamic viscosity of the gas, α is the degree of permeability, C2Coefficient of inertial resistance, α and C2Can be calculated by the following method:
the axial pressure drop of a region of porous media can be expressed as:
Figure FDA0002319948710000029
in the formula (9), foscFor the friction factor, according to empirical formula:
Figure FDA00023199487100000210
in equation (10), Re is the reynolds number of the fluid in the porous medium region and can be expressed as:
Figure FDA00023199487100000211
by combining the formula (9), the formula (10) and the formula (11), the degree of penetration α and the coefficient of inertial resistance C can be determined2
Figure FDA0002319948710000031
Figure FDA0002319948710000032
Energy equation:
Figure FDA0002319948710000033
in the formula (14), psAnd EsDensity and total energy, k, respectively, of the solideffEffective thermal conductivity for porous media;
effective thermal conductivity k of porous media due to the effect of thermal contact resistanceeffMuch less than its average thermal conductivity, the following modifications can be made:
Figure FDA0002319948710000034
in formula (15), ksAnd kfThe thermal conductivity of the solid and gas respectively.
10. The method as claimed in claim 1, wherein in the sixth step, the flow field is initialized and numerical solution is started according to the established CFD model of the aircraft air circulation refrigeration system, given solver control parameters, a discrete format and residual convergence criteria.
CN201911293882.4A 2019-12-13 2019-12-13 Simulation analysis method of airplane air circulation refrigeration system based on CFD technology Pending CN111027155A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911293882.4A CN111027155A (en) 2019-12-13 2019-12-13 Simulation analysis method of airplane air circulation refrigeration system based on CFD technology

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911293882.4A CN111027155A (en) 2019-12-13 2019-12-13 Simulation analysis method of airplane air circulation refrigeration system based on CFD technology

Publications (1)

Publication Number Publication Date
CN111027155A true CN111027155A (en) 2020-04-17

Family

ID=70209512

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911293882.4A Pending CN111027155A (en) 2019-12-13 2019-12-13 Simulation analysis method of airplane air circulation refrigeration system based on CFD technology

Country Status (1)

Country Link
CN (1) CN111027155A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112214937A (en) * 2020-09-09 2021-01-12 黄冈格罗夫氢能汽车有限公司 Optimal design method and system for hydrogen energy automobile cooling pipeline
CN112765850A (en) * 2021-01-19 2021-05-07 潍柴(扬州)亚星新能源商用车有限公司 Design method for air conditioning system in large VAN vehicle cabin
CN114048697A (en) * 2021-12-31 2022-02-15 中国飞机强度研究所 Airplane test airflow organization design method

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130018641A1 (en) * 2011-07-12 2013-01-17 Ingrain, Inc. Method For Simulating Fractional Multi-Phase/Multi-Component Flow Through Porous Media
US20140257771A1 (en) * 2011-11-30 2014-09-11 Ming Lu Numerical simulation method for aircrasft flight-icing
CN105526729A (en) * 2015-12-01 2016-04-27 中国科学院上海技术物理研究所 Simulation analysis method of light-weight miniature coaxial pulse tube refrigerator based on CFD technology
CN106682346A (en) * 2017-01-05 2017-05-17 中南大学 Method for optimizing complicated member gas-quenching system based on CFD software
WO2017084106A1 (en) * 2015-11-20 2017-05-26 田川 System and method for numerical simulation of aircraft flow field
CN110020465A (en) * 2019-03-17 2019-07-16 浙江大学 A kind of analog analysing method of the liquid tinned food heat-sterilization process based on CFD

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130018641A1 (en) * 2011-07-12 2013-01-17 Ingrain, Inc. Method For Simulating Fractional Multi-Phase/Multi-Component Flow Through Porous Media
US20140257771A1 (en) * 2011-11-30 2014-09-11 Ming Lu Numerical simulation method for aircrasft flight-icing
WO2017084106A1 (en) * 2015-11-20 2017-05-26 田川 System and method for numerical simulation of aircraft flow field
CN105526729A (en) * 2015-12-01 2016-04-27 中国科学院上海技术物理研究所 Simulation analysis method of light-weight miniature coaxial pulse tube refrigerator based on CFD technology
CN106682346A (en) * 2017-01-05 2017-05-17 中南大学 Method for optimizing complicated member gas-quenching system based on CFD software
CN110020465A (en) * 2019-03-17 2019-07-16 浙江大学 A kind of analog analysing method of the liquid tinned food heat-sterilization process based on CFD

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
周晨阳等: "基于CFD的通信机房流场模型建立与分析", 《测控技术》 *
杨磊等: "微型保鲜库预冷气流组织的仿真分析", 《计算机仿真》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112214937A (en) * 2020-09-09 2021-01-12 黄冈格罗夫氢能汽车有限公司 Optimal design method and system for hydrogen energy automobile cooling pipeline
CN112765850A (en) * 2021-01-19 2021-05-07 潍柴(扬州)亚星新能源商用车有限公司 Design method for air conditioning system in large VAN vehicle cabin
CN114048697A (en) * 2021-12-31 2022-02-15 中国飞机强度研究所 Airplane test airflow organization design method
CN114048697B (en) * 2021-12-31 2022-04-15 中国飞机强度研究所 Airplane test airflow organization design method

Similar Documents

Publication Publication Date Title
CN111027155A (en) Simulation analysis method of airplane air circulation refrigeration system based on CFD technology
Chen et al. Pressure analysis on two-step high pressure reducing system for hydrogen fuel cell electric vehicle
CN105480433B (en) The test device and method of aircraft fuel system under a kind of simulated altitude environment
CN105205281A (en) Heatproof integrated design method for reentry trajectory
CN105526729A (en) Simulation analysis method of light-weight miniature coaxial pulse tube refrigerator based on CFD technology
Ye et al. Transient flow performance and heat transfer characteristic in the cylinder of hydraulic driving piston hydrogen compressor during compression stroke
Zhang et al. Heat transfer in a long, thin tube section of an air compressor: an empirical correlation from CFD and a thermodynamic modeling
Zhou et al. Aerodynamic characteristics of the evacuated tube maglev train considering the suspension gap
CN103105017B (en) Environment test device refrigeration system and energy-conservation operating mode fast determination method thereof
Li et al. Dynamic simulation model for three-wheel air-cycle refrigeration systems in civil aircrafts
CN111143939A (en) Simulation analysis system of airplane air circulation refrigeration system based on CFD technology
Fernández-Yáñez et al. Impact of relative position vehicle-wind blower in a roller test bench under climatic chamber
Zhang et al. Reliability analysis of bleed air anti-icing system based on subset simulation method
Hairong et al. Comparison of turbulence models for multiphase-flow oscillating heat transfer enhancement
Xue et al. Mixing characteristics in a ventilated room with non-isothermal ceiling air supply
Li et al. Effect of pipes in heat pump system on electric vehicle energy saving
Liu et al. Comparative study of numerical simulation of indoor thermal environment in the pattern of personalized ventilation and stratum ventilation
Mokheimer et al. Comparative analysis of different configuration domestic refrigerators: a computational fluid dynamics approach
CN202032792U (en) Fin type integrated evaporator
CN103334985A (en) Hydraulic oil cooling oil tank
Chen et al. Research on 1D-2D co-simulation for cabin air environment accident
CN112710058A (en) Excel-based air conditioner model selection method
Pang et al. Performance of an integrated thermal management system for helicopter
Aldas et al. The Investigation of Gas Ejector Performance using CFD Modelling
Chiappini et al. Sizing and Optimization of a Vortex Tube for Electric Vehicle HVAC Purposes

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination