CN110083952A - Carbon dioxide train air-conditioning emulation mode - Google Patents
Carbon dioxide train air-conditioning emulation mode Download PDFInfo
- Publication number
- CN110083952A CN110083952A CN201910364215.4A CN201910364215A CN110083952A CN 110083952 A CN110083952 A CN 110083952A CN 201910364215 A CN201910364215 A CN 201910364215A CN 110083952 A CN110083952 A CN 110083952A
- Authority
- CN
- China
- Prior art keywords
- carbon dioxide
- compressor
- conditioning
- air
- heat transfer
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Automation & Control Theory (AREA)
- Air Conditioning Control Device (AREA)
Abstract
The present invention relates to a kind of emulation modes for carbon dioxide train air-conditioning, this method comprises: the status predication of carbon dioxide piston compressor, the information such as outlet temperature, flow, the wasted work of compressor can be calculated according to the inlet and outlet pressure of compressor and the dimension information of compressor itself；The performance of carbon dioxide fin-tube type heat exchanger calculates, and can calculate the heat exchange amount and air side outlet information of heat exchanger according to the refrigerant side access information and air side access information of carbon dioxide fin-tube type heat exchanger；The steady state predictions of carbon dioxide train air conditioning system can calculate state when carbon dioxide train air conditioning system is run to stable state in given external environment work information and system operational parameters.The present invention can assist optimization design carbon dioxide system and its components, predict before the experiments system running state, and parameter setting when can optimize the system operation, so that system realizes the optimization of performance under different operating conditions.
Description
Technical field
The present invention relates to refrigeration and numerical arts, and in particular, to it is a kind of can continuous iteration optimization titanium dioxide
Carbon train air-conditioning emulation mode.
Background technique
Due to using the system with high GWP (Global warming potential) value in traditional train air conditioning system
Cryogen, meeting increase from global warming, so using having the refrigerant of low GWP value imperative.Carbon dioxide is as a kind of nature
The advantages that working medium refrigerant has GWP value low (GWP=1), nontoxic non-combustible, at low cost, is current train air conditioner refrigerating
The main selection of agent substitution.However the operational mode of carbon dioxide air conditioning system and there are larger using the systems of conventional refrigerants
Difference: in most cases, carbon dioxide air conditioning system is run under Trans-critical cycle mode, and uses the system of conventional refrigerants logical
Often run under sub-critical mode.Due to this characteristic, there is no degree of supercooling, and its exhaust pressure in carbon dioxide air conditioning system
Power can also change in a larger range, and the emulation for carbon dioxide air conditioning system is caused to be different from using conventional refrigeration
The system of agent.
Since conventional refrigeration emulation mode and carbon dioxide system mismatch, it is therefore desirable to be directed to carbon dioxide train
The characteristic of air-conditioning system individually develops a set of emulation mode, for probing into the operation feelings inside carbon dioxide train air conditioning system
Condition, and designed for optimization system.
Summary of the invention
In view of the deficiencies of the prior art, the present invention proposes a kind of carbon dioxide train air-conditioning emulation modes, can assist excellent
Change design carbon dioxide system and its components, system running state is predicted before the experiments, and can optimization system fortune
Parameter setting when row, so that system realizes the optimization of performance under different operating conditions.
The present invention is achieved through the following technical solutions:
The first aspect of the present invention provides a kind of trend prediction method of carbon dioxide piston compressor, can be according to pressure
The inlet and outlet pressure information of contracting machine obtains the isentropic efficiency of compressor, mechanical efficiency and volumetric efficiency in conjunction with empirical equation,
And combination operating condition and sizecalculation go out the discharge state of compressor, energy consumption and flow information.
Preferably, the form that the formula of compressor efficiency is combined using multinomial and logarithm is calculated.
Preferably, compressor efficiency calculation formula is fitted using existing experimental data, and use variance as
Loss function in fit procedure.
Preferably, after calculating several efficiency of compressor, using flow -> outlet enthalpy -> outlet temperature -> compressor
Computation sequence as wasted work realizes compressor status predication.
The second aspect of the present invention provides a kind of Calculation Methods for Performance of carbon dioxide fin-tube type heat exchanger, by titanium dioxide
Carbon fin-tube type heat exchanger splits into lesser heat transfer unit, along refrigerant and air side flow sequence successively Numerical heat transfer list
The heat exchange amount and discharge state of member, finally obtain total heat exchange amount and discharge state.
Preferably, the heat transfer unit length of one section of carbon dioxide fin-tube type heat exchanger should be no more than 1cm, to guarantee to calculate essence
Degree.
Preferably, by the inlet and outlet state average value of heat transfer unit as the integrality of heat transfer unit.
Preferably, it when the heat exchange amount state of Numerical heat transfer unit, needs first to judge the state of heat transfer unit (main
It is whether to be overcritical, if be the judgement of two-phase).
Preferably, by postcritical carbon dioxide be approximately monophasic fluid go carry out pressure drop and conduct heat formula calculating.
It preferably, is approximately that dry air calculates by air in the calculating process of air cooler；And in the meter of evaporator
Needed during calculating first according to evaporating temperature and air humidity, distinguish heat transfer unit be in be dry state (no analysis is wet) or
Wet condition (has analysis wet), then is calculated.
Preferably, it in the iterative process of heat transfer unit, is optimized using Newton iteration method, calculates effect to improve
Rate.
The third aspect of the present invention provides a kind of steady state predictions method of carbon dioxide train air conditioning system, given outer
When boundary's environmental working condition and system operational parameters, the operating status of system entirety energy consumption and each components can be calculated.
Preferably, using the pressure at expulsion of system and the degree of superheat as the system operational parameters set.
Preferably, it is calculated the pressure of inspiration(Pi) of system as assumption value input system, and in the iterative process in later period
In constantly adjust.
Preferably, according to room temperature, initial value design is carried out to the pressure of inspiration(Pi) of system.
Preferably, using evaporator inlet enthalpy and air cooler outlet enthalpy deviation as calculate whether it is convergent judge according to
According to.
Preferably, it in the stable state calculating process of carbon dioxide train air conditioning system, is optimized using Newton iteration method,
To improve computational efficiency.
Compared with prior art, the present invention has the advantages that are as follows: the present invention can assist optimization design titanium dioxide
Carbon system and its components, before the experiments predict system running state, and parameter when can optimize the system operation is set
It is fixed, so that system realizes the optimization of performance under different operating conditions.
Detailed description of the invention
Fig. 1 is structure and functional schematic of the invention；
Fig. 2 is compressor state computation flow chart of the invention；
Fig. 3 is that the heat transfer unit of carbon dioxide fin-tube type heat exchanger of the invention iterates to calculate flow chart；
Fig. 4 is the overall calculation flow chart of carbon dioxide fin-tube type heat exchanger of the invention；
Fig. 5 is that the stable state of carbon dioxide train air conditioning system of the invention iterates to calculate flow chart.
Specific embodiment
The present invention is described in detail combined with specific embodiments below.Following embodiment will be helpful to the technology of this field
Personnel further understand the present invention, but the invention is not limited in any way.It should be pointed out that the ordinary skill of this field
For personnel, without departing from the inventive concept of the premise, various modifications and improvements can be made.These belong to the present invention
Protection scope.
Embodiment
Shown in specific embodiment Fig. 1 to Fig. 5.
Fig. 1 is structure and functional schematic of the invention.By external environment operating condition, system parts dimensional parameters and system
Input of the operating parameter as this emulation mode, by carbon-dioxide gas compressor state and carbon dioxide fin-tube type heat exchanger
The performance of (including air cooler and evaporator) calculates, and can get the operating status of carbon dioxide train air-conditioning steady-state system entirety.
By different size of heat exchanger Parameters Input Units, it can get performance of the heat exchanger under different dimensional parameters, selected most with this
Excellent dimensional parameters；By different system operational parameters input systems (including system discharge pressure and degree of superheat), can get
Performance of the system under different operating parameters, optimal system operational parameters are selected with this.The summation of the position Fig. 1 system function, figure
2- Fig. 5 is its concrete implementation method.
Fig. 2 is the calculation process of compressor state.It first can in conjunction with empirical equation according to the inlet and outlet pressure of compressor
Calculate the isentropic efficiency of compressor, mechanical efficiency and volumetric efficiency.Wherein, the empirical equation of compressor efficiency can be according to pressure
Contracting machine experimental data is fitted using multinomial and natural logrithm using minimum variance as target.Volumetric efficiency, which calculates, to be completed
Afterwards, compressor flowrate can be calculated in conjunction with compressor volume and revolving speed, for using when subsequent heat exchanger design；Isentropic efficiency calculates
After the completion, the enthalpy that can calculate compressor outlet can get compressor exit temperature in conjunction with compressor delivery pressure；Finally
The practical wasted work of compressor can be calculated according to mechanical efficiency, compressor flowrate and compressor outlet enthalpy.
Fig. 3 is the iterative calculation flow chart of one heat transfer unit of carbon dioxide fin-tube type heat exchanger.Given heat transfer unit
After inlet condition (refrigerant side and air side), it is first assumed that the outlet enthalpy of refrigerant side, and according to enthalpy and refrigerant
Pressure judges whether carbon dioxide is in two-phase section in this unit；According to flow and outlet enthalpy, calculating system that can be explicit
The theoretical heat exchange amount of cryogen side, and according to energy conservation principle, calculate the discharge state of air side.It is worth noting that, working as
When carbon dioxide is in a supercritical state in heat transfer unit, one-way fluid processing can be considered as.Based on being made in heat transfer unit
The state of cryogen and air is approximate constant it is assumed that use the mean value of entrance and exit state as the state of entire heat transfer unit,
Can be according to the heat transfer empirical equation of refrigerant side and air side, the explicit heat transfer coefficient for calculating air side and refrigerant side,
And the attainable practical heat exchange amount of heat exchange unit institute is calculated based on heat transfer coefficient.Pass through the inclined of theoretical heat exchange amount and practical heat exchange amount
Whether difference can determine whether to calculate and restrain.The outlet enthalpy that refrigerant side is constantly corrected using Newton iteration method can accelerate convergence speed
Degree.
Fig. 4 is the calculation flow chart of carbon dioxide fin-tube type heat exchanger entirety.One-row heat exchanger with four processes is
Example, each process are divided into five junior units.Refrigerant enters from the entrance of heat exchanger, after flowing through five heat transfer units, covers
First pass, and enter downstream；This process is repeated, until flowing out from heat exchanger exit.It is worth noting that, a list
Member is actually a bit of of a flat tube, and may include a more than flat tube in a process.A but because stream
Each flat tube has equivalence in journey, so need to only calculate the heat exchange amount of a heat transfer unit, can know that a process neutralizes
The heat exchange amount of other parallel all heat transfer units of this heat transfer unit.
Fig. 5 is that the stable state of carbon dioxide train air conditioning system iterates to calculate flow chart.Set system discharge pressure and mistake
After temperature, it is assumed that suction pressure of compressor can calculate compressor outlet state and compressor flowrate；Using flow information and
The calculation method of above-mentioned carbon dioxide fin-tube type heat exchanger can calculate separately the discharge state of air cooler and the entrance shape of evaporator
Whether state judges to calculate and restrain by the deviation that evaporator inlet enthalpy and air cooler export enthalpy.Not using Newton iteration method
Disconnected amendment suction pressure of compressor, can accelerate convergence rate.Can be obtained after convergence system entirety energy consumption and each components
Operation information.
Specific embodiments of the present invention are described above.It is to be appreciated that the invention is not limited to above-mentioned
Particular implementation, those skilled in the art can make various deformations or amendments within the scope of the claims, this not shadow
Ring substantive content of the invention.
Claims (11)
1. a kind of carbon dioxide train air-conditioning emulation mode, which comprises the following steps:
First, using external environment operating condition, system parts dimensional parameters and system operational parameters as the input of emulation mode, lead to
The state for crossing carbon dioxide piston compressor trend prediction method prediction carbon-dioxide gas compressor, passes through carbon dioxide fin-tube type
The Calculation Methods for Performance of heat exchanger calculates the performance of carbon dioxide fin-tube type heat exchanger；
Second, the trend prediction method of carbon dioxide piston compressor and the performance of carbon dioxide fin-tube type heat exchanger are calculated
Method is encapsulated as function, is calculated by calling them to complete the stable state of carbon dioxide train air conditioning system；
Different size of heat exchanger parameters and different system operational parameters input systems it is optimal to be selected heat exchanger by third
Dimensional parameters and optimal system operational parameters.
2. carbon dioxide train air-conditioning emulation mode according to claim 1, it is characterised in that the carbon dioxide piston
Formula compressor trend prediction method are as follows: according to the inlet and outlet pressure information of compressor, in conjunction with empirical equation, obtain compressor etc.
Entropic efficiency, mechanical efficiency and volumetric efficiency calculate discharge state, energy consumption and the flow information of compressor.
3. carbon dioxide train air-conditioning emulation mode according to claim 2, it is characterised in that the empirical equation, with big
Based on capacity compressor experimental data, using multinomial and natural logrithm, it is fitted by target of minimum variance.
4. carbon dioxide train air-conditioning emulation mode according to claim 1, it is characterised in that the carbon dioxide pipe wing
The Calculation Methods for Performance of formula heat exchanger are as follows: carbon dioxide fin-tube type heat exchanger is split into lesser heat transfer unit, along refrigeration
The heat exchange amount and discharge state of agent and air side flow sequence successively Numerical heat transfer unit, finally obtain total heat exchange amount and outlet
State.
5. carbon dioxide train air-conditioning emulation mode according to claim 4, it is characterised in that the carbon dioxide pipe wing
The heat transfer unit of formula heat exchanger, since heat transfer unit is smaller, internal state is considered approximately equal.
6. carbon dioxide train air-conditioning emulation mode according to claim 4, it is characterised in that the carbon dioxide pipe wing
The heat transfer unit of formula heat exchanger, is calculated using alternative manner: assuming that heat transfer unit discharge state, obtains theoretical heat exchange amount,
And practical heat exchange amount is obtained according to refrigerant side and the heat transfer area and heat transfer coefficient of air side, pass through theoretical heat exchange amount and reality
Whether the deviation of heat exchange amount restrains to judge to calculate, and heat transfer unit discharge state is adjusted if not restraining, iteration goes on.
7. carbon dioxide train air-conditioning emulation mode according to claim 6, it is characterised in that the heat transfer unit iteration
Calculation method accelerates convergence process using Newton iteration method.
8. carbon dioxide train air-conditioning emulation mode according to claim 6, it is characterised in that the refrigerant side heat transfer
Coefficient calculates, and according to the difference of carbon dioxide status, selects different carbon dioxide heat transfer coefficient calculation formula.
9. carbon dioxide train air-conditioning emulation mode according to claim 6, it is characterised in that the air side heat transfer
Coefficient calculates, and according to the dimensional parameters of fin and speed air flow, temperature, humidity information, completes to calculate using fitting formula.
10. carbon dioxide train air-conditioning emulation mode according to claim 1, which is characterized in that the carbon dioxide column
The stable state of vehicle air-conditioning system calculates, and is calculated using alternative manner: by assuming that compressor inlet pressure, calculates compressor and go out
Mouthful, air cooler outlet and evaporator inlet state are judged by the deviation that evaporator inlet enthalpy and air cooler export enthalpy
It calculates and whether restrains.
11. carbon dioxide train air-conditioning emulation mode according to claim 1, it is characterised in that the carbon dioxide train
The stable state of air-conditioning system calculates, and inputs mould using compressor delivery pressure and the suction port of compressor degree of superheat as system operational parameters
Type is not changed in an iterative process.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910364215.4A CN110083952A (en) | 2019-04-30 | 2019-04-30 | Carbon dioxide train air-conditioning emulation mode |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910364215.4A CN110083952A (en) | 2019-04-30 | 2019-04-30 | Carbon dioxide train air-conditioning emulation mode |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110083952A true CN110083952A (en) | 2019-08-02 |
Family
ID=67418278
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910364215.4A Pending CN110083952A (en) | 2019-04-30 | 2019-04-30 | Carbon dioxide train air-conditioning emulation mode |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110083952A (en) |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103158490A (en) * | 2011-12-19 | 2013-06-19 | 杭州三花研究院有限公司 | Car air-conditioning system |
| CN103542621A (en) * | 2013-09-26 | 2014-01-29 | 西安交通大学 | Design method of general combination pipe diameter air conditioner heat exchange equipment fluid passage |
| CN104572246A (en) * | 2014-12-08 | 2015-04-29 | 江苏大学 | Interactive optimization method for improving internal environment of building |
| US20150316302A1 (en) * | 2014-05-02 | 2015-11-05 | Gary Giubilo | Refrigeration apparatus and method |
| CN105678429A (en) * | 2016-02-19 | 2016-06-15 | 西南交通大学 | Method for achieving multi-discipline multi-target optimization of air-conditioning system of electric car |
| CN106169018A (en) * | 2016-06-20 | 2016-11-30 | 珠海格力电器股份有限公司 | Heat exchange quantity simulation method and device |
| CN106503388A (en) * | 2016-11-09 | 2017-03-15 | 东南大学 | air conditioning system characteristic recognition method |
| CN106871391A (en) * | 2017-04-26 | 2017-06-20 | 上海科凌能源科技有限公司 | Performance for Air Conditioning Systems online test method based on limited measure node |
| CN109323475A (en) * | 2018-09-11 | 2019-02-12 | 西安交通大学 | A kind of critical-cross carbon dioxide system and its optimization coupling process |
-
2019
- 2019-04-30 CN CN201910364215.4A patent/CN110083952A/en active Pending
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103158490A (en) * | 2011-12-19 | 2013-06-19 | 杭州三花研究院有限公司 | Car air-conditioning system |
| CN103542621A (en) * | 2013-09-26 | 2014-01-29 | 西安交通大学 | Design method of general combination pipe diameter air conditioner heat exchange equipment fluid passage |
| US20150316302A1 (en) * | 2014-05-02 | 2015-11-05 | Gary Giubilo | Refrigeration apparatus and method |
| CN104572246A (en) * | 2014-12-08 | 2015-04-29 | 江苏大学 | Interactive optimization method for improving internal environment of building |
| CN105678429A (en) * | 2016-02-19 | 2016-06-15 | 西南交通大学 | Method for achieving multi-discipline multi-target optimization of air-conditioning system of electric car |
| CN106169018A (en) * | 2016-06-20 | 2016-11-30 | 珠海格力电器股份有限公司 | Heat exchange quantity simulation method and device |
| CN106503388A (en) * | 2016-11-09 | 2017-03-15 | 东南大学 | air conditioning system characteristic recognition method |
| CN106871391A (en) * | 2017-04-26 | 2017-06-20 | 上海科凌能源科技有限公司 | Performance for Air Conditioning Systems online test method based on limited measure node |
| CN109323475A (en) * | 2018-09-11 | 2019-02-12 | 西安交通大学 | A kind of critical-cross carbon dioxide system and its optimization coupling process |
Non-Patent Citations (3)
| Title |
|---|
| Y.B.TAO；Y.L.HE；W.Q.TAO；Z.G.WU: ""Experimental study on the performance of CO2 residential air-conditioning system with an internal heat exchanger"", 《ENERGY CONVERSION AND MANAGEMENT》 * |
| 周成君: ""CO2冷库制冷系统设计与实验研究"", 《中国优秀硕士学位论文全文数据库 (工程科技Ⅱ辑)》 * |
| 杨涛，陈江平，陈芝久: ""二氧化碳汽车空调系统的动态仿真"", 《上海汽车》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wang et al. | Optimization of refrigeration system with gas-injected scroll compressor | |
| Elbel et al. | Experimental validation of a prototype ejector designed to reduce throttling losses encountered in transcritical R744 system operation | |
| CN114838531B (en) | Method and device for adjusting temperature of refrigerant in heat exchanger and air conditioner | |
| Kim et al. | Optimization of the injection-port geometries of a vapor injection scroll compressor based on SCOP under various climatic conditions | |
| CN102831302B (en) | Performance calculating method of finned tube evaporator under frosting working condition | |
| CN104634009B (en) | Control method of air conditioner circulating device | |
| CN107366621B (en) | Rolling rotor compressor with three-stage air supplement and air conditioning system | |
| CN101504177A (en) | Air conditioning apparatus and method for determining the amount of refrigerant of air-conditioning apparatus | |
| CN109543221B (en) | On-site fault diagnosis method for air conditioner water chilling unit | |
| Cho et al. | Optimization of injection holes in symmetric and asymmetric scroll compressors with vapor injection | |
| CN105627612B (en) | The refrigeration control method of outdoor unit refrigerant pipe-line system, air conditioner and air conditioner | |
| CN103229004B (en) | Aircondition | |
| CN104634008B (en) | Control method of air conditioner | |
| CN107461874A (en) | Air conditioner defrosting control method and air conditioner | |
| CN110243047A (en) | Method and device for calculating heat exchange capacity of heat exchange equipment, heat exchange equipment and system | |
| CN105910220B (en) | VRF divided gauging method | |
| CN103168204A (en) | Device for estimating flowrate of heating medium, heat source device, and method for estimating flowrate of heating medium | |
| Jeon et al. | Numerical study on the optimal design of injection-hole geometries of a twin rotary compressor in a liquid injection heat pump | |
| CN110175403A (en) | It is a kind of applied to refrigeration or air-conditioning equipment parameterisable dynamic simulation computing system | |
| Liu et al. | A modeling study on developing the condensing-frosting performance maps for a variable speed air source heat pump | |
| CN105588365B (en) | A kind of heat-flash type outdoor unit, heat pump system and its control method | |
| Cho et al. | Refrigerant charge reduction in R600a domestic refrigerator-freezer by optimizing hot-wall condenser geometry | |
| Niu et al. | Numerical simulation of CO2 two-stage compression refrigeration system with external intercooler | |
| CN110083952A (en) | Carbon dioxide train air-conditioning emulation mode | |
| CN110023684A (en) | Determination method, the determining device of tube diameter and the refrigerating plant of tube diameter |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190802 |
|
| WD01 | Invention patent application deemed withdrawn after publication |