CN112287579A - Binary non-azeotropic mixed working medium condensation evaporation simulation method considering component distribution - Google Patents
Binary non-azeotropic mixed working medium condensation evaporation simulation method considering component distribution Download PDFInfo
- Publication number
- CN112287579A CN112287579A CN202011153247.9A CN202011153247A CN112287579A CN 112287579 A CN112287579 A CN 112287579A CN 202011153247 A CN202011153247 A CN 202011153247A CN 112287579 A CN112287579 A CN 112287579A
- Authority
- CN
- China
- Prior art keywords
- working medium
- component
- mixed working
- azeotropic mixed
- binary non
- 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.)
- Granted
Links
Images
Classifications
-
- 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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
-
- 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/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Engine Equipment That Uses Special Cycles (AREA)
- Organic Low-Molecular-Weight Compounds And Preparation Thereof (AREA)
Abstract
A binary non-azeotropic mixed working medium condensation and evaporation simulation method considering component distribution is used for researching and developing the refrigeration performance of a binary non-azeotropic mixed working medium in the condensation or evaporation process, and according to the difference ratio between the ideal equilibrium state and the actual state of the component content of the two components of the binary non-azeotropic mixed working medium in the condensation or evaporation process, the ratio of the condensation or evaporation rate of the two components is used as the ratio of the condensation or evaporation rate of the two components, the mass and energy transfer relational expression between the whole binary non-azeotropic mixed working medium and the gas phase and the liquid phase in the condensation or evaporation process of the two components is set, and finally the heat exchange coefficients of the two components of the binary non-azeotropic mixed working medium and the component distribution diagram in equipment are obtained; the method can be used for two-dimensional model calculation and three-dimensional model calculation, and can save the research and development cost and period of the binary non-azeotropic mixed working medium.
Description
Technical Field
The invention relates to the technical field of thermodynamics and heating ventilation air conditioners, in particular to a binary non-azeotropic mixed working medium condensation evaporation simulation method considering component distribution.
Technical Field
The refrigerating working medium is blood of a refrigerating air-conditioning system, and is mainly used in the fields of heating ventilation air-conditioning such as household air-conditioning, automobile air-conditioning, water chilling unit systems and the like. The refrigeration working media mainly comprise two categories of pure working media and mixed working media, and the refrigeration working media influencing global warming gradually face to be eliminated along with the rise of environmental protection requirements in recent years; however, it is also not wise to simply pursue the effect of the refrigeration duty on global warming and ignore its refrigeration performance. Compared with pure working medium, the binary mixed working medium can not only reduce the effect of the working medium on global warming, but also maintain the refrigeration performance, so that the binary mixed working medium is selected as an excellent alternative method.
The refrigeration performance of the binary mixed working medium is obviously influenced by the mixing proportion of the binary mixed working medium, however, as the types and the mixing proportion of the mixed working medium are numerous, the refrigeration performance of the binary mixed working medium is estimated by simply adopting an experimental measurement method, so that a large amount of time and economic cost are spent, and the significance of the prediction method for researching the refrigeration performance of the binary mixed working medium is great. The binary mixed working medium comprises an azeotropic working medium and a non-azeotropic working medium, the condensation temperature and the evaporation temperature of the azeotropic working medium are the same, and when the condensation or evaporation phenomenon occurs, the content ratio of the two components in the gas phase and the liquid phase does not change and is always equal to the content ratio in the initial mixed state; the condensation temperature and the evaporation temperature of the non-azeotropic working medium are different, and the content ratio of two components in a gas phase and a liquid phase can be changed continuously when condensation or evaporation occurs, so that the calculation of the condensation or evaporation rate of the two components of the binary non-azeotropic mixed working medium has certain difficulty.
In the existing scientific research and production process, a binary non-azeotropic mixed working medium condensation and evaporation simulation method considering component distribution is lacked, and refrigeration performance parameters such as a component distribution diagram, a heat exchange coefficient and the like of the binary non-azeotropic mixed working medium in equipment cannot be obtained.
Disclosure of Invention
In order to overcome the defects of the prior art, the invention aims to provide a binary non-azeotropic mixed working medium condensation and evaporation simulation method considering component distribution, which can calculate the rate ratio of two components in the condensation or evaporation process of a binary non-azeotropic mixed working medium mixed in any proportion, and further combines a finite element simulation technology to obtain the refrigeration performance parameters such as a component distribution diagram, a heat exchange coefficient and the like of the binary non-azeotropic mixed working medium in equipment.
In order to achieve the purpose, the invention adopts the technical scheme that:
the binary non-azeotropic mixed working medium condensation evaporation simulation method considering component distribution comprises the following steps:
1) according to the type, initial ratio and pressure of the binary non-azeotropic mixed working medium A-B, consulting and calculating software NIST to draw a 'saturation temperature-mass fraction' diagram of the binary non-azeotropic mixed working medium A-B under a certain fixed pressure;
2) when the temperature of the binary non-azeotropic mixed working medium A-B is higher than the dew point temperature or lower than the bubble point temperature, setting the condensation or evaporation rate ratio of the binary non-azeotropic mixed working medium A-B as the mass fraction ratio of the mixed working medium in the initial filling state, wherein the condensation or evaporation amount of the component A and the component B of the binary non-azeotropic mixed working medium A-B occupies the ratio alpha of the total condensation or evaporation amount of the binary non-azeotropic mixed working medium A-BAAnd alphaBRespectively as follows:
αArepresents the ratio of the condensation or evaporation capacity of the component A to the total condensation or evaporation capacity of the binary non-azeotropic mixed working medium A-B, alphaBRepresents the ratio of the condensation or evaporation capacity of the component B to the total condensation or evaporation capacity of the binary non-azeotropic mixed working medium A-B, w0The component B occupies the mass ratio of the total binary non-azeotropic mixed working medium A to B in the initial filling state;
at the moment, the saturation temperature of the binary non-azeotropic mixed working medium A-B is set as follows:
Tsatthe saturated temperature of the binary non-azeotropic mixed working medium A-B is represented, and when the condensation process is calculated, the bubble point temperature corresponding to the mass fraction of the component A and the component B in the initial filling state is a fixed value; when the evaporation process is calculated, the dew point temperature corresponding to the mass fractions of the component A and the component B in the initial mixing state is a fixed value; p represents the pressure of the binary non-azeotropic mixed working medium A-B and is set as a constant value in the whole calculation process;
when the temperature of the binary non-azeotropic mixed working medium A-B is lower than the dew point temperature and higher than the bubble point temperature, setting the condensation or evaporation rate ratio of the binary non-azeotropic mixed working medium A-B as DeltamA/ΔmB,ΔmAAnd Δ mBIs based on the following formula:
ΔmArepresents the difference between the actual mass of component A in the gas phase and the ideal equilibrium mass, Δ mBRepresents the difference between the mass of component B in the gas phase in the actual state and the mass in the ideal equilibrium state,mis the total mass of binary non-azeotropic mixed working medium A-B, x3Is the dryness, x, of an ideal equilibrium state 3 Is the dryness of the actual state, w3”Is the mass fraction of component B in the gas phase at the ideal equilibrium state, w 3” Is the mass fraction of component B in the gas phase in the actual state;
the condensation or evaporation capacity of the component A and the component B of the binary non-azeotropic mixed working medium A-B occupies the ratio alpha of the total condensation or evaporation capacity of the binary non-azeotropic mixed working medium A-BAAnd alphaBRespectively as follows:
at the moment, the saturation temperature of the binary non-azeotropic mixed working medium A-B is set as follows:
Tsatthe saturation temperature of the binary non-azeotropic mixed working medium A-B is represented, and when the condensation process is calculated, the bubble point temperature corresponding to the mass fraction of the component A and the component B in the actual calculation process is a variable; when the evaporation process is calculated, the dew point temperature corresponding to the mass fractions of the component A and the component B in the actual calculation process is a variable; w is the ratio of the component B to the total binary non-azeotropic mixed working medium A-B in the actual calculation process;
3) obtaining alphaAAnd alphaBAnd then, combining finite element simulation software FLUENT, selecting a component transport model, a VOF two-phase flow model and an SST turbulence model, setting the binary non-azeotropic mixed working medium A-B and the mass and energy transfer between the component A and the component B thereof according to the following formula in a source item mode, and carrying out calculation to obtain the heat exchange coefficient of the binary non-azeotropic mixed working medium A-B and the component distribution diagram of the component A and the component B in the equipment.
If the temperature of the binary non-azeotropic mixed working medium A-B is higher than the saturation temperature, the mass transfer amounts of the component A, the component B and the binary non-azeotropic mixed working medium A-B in gas phase and liquid phase are as follows:
in the formula (I), the compound is shown in the specification,represents the mass transfer of the gas phase of component a,represents the mass transfer amount of the liquid phase of component a,represents the mass transfer of the gas phase of component B,represents the mass transfer of the liquid phase of component B,rthe condensing evaporation coefficient is 0.01-3000000, betalRepresents the volume fraction, rho, of the liquid phase of the binary non-azeotropic mixed working medium A-BlExpressing the density of the liquid phase of the binary non-azeotropic mixed working medium A-B, T expressing the temperature of the binary non-azeotropic mixed working medium A-B, TsatRepresenting the saturation temperature of the binary non-azeotropic mixed working medium A-B.
If the temperature of the binary non-azeotropic mixed working medium A-B is lower than the saturation temperature, the mass transfer amounts of the component A, the component B and the binary non-azeotropic mixed working medium A-B in gas phase and liquid phase are as follows:
in the formula, betaνRepresenting the volume fraction, rho, of the binary non-azeotropic mixed working medium A-B gas phaseνThe density of the gas phase of the binary non-azeotropic mixed working medium A-B is shown.
The energy transfer amounts of the component A, the component B and the binary non-azeotropic mixed working medium A-B in gas phase and liquid phase are as follows:
in the formula, EARepresents the amount of energy transfer of component a,denotes the mass transfer of the liquid phase of component A, γADenotes the latent heat of component A, EBRepresents the amount of energy transfer of component B,denotes the mass transfer of the liquid phase of component B, γBThe latent heat of the component B is shown, and E represents the total energy transfer of the binary non-azeotropic mixed working medium A-B.
The two components of the binary non-azeotropic mixed working medium A-B are refrigeration working media or low-temperature working media, the refrigeration working media are R22, R32, R123, R125, R134a, R152a, R161, R1234yf and R1234ze (E), and the low-temperature working media are methane, ethane, propane and butane.
The ideal equilibrium state in the step 2) refers to the state of the binary non-azeotropic mixed working medium A-B at each temperature in the 'saturation temperature-mass fraction' diagram under a certain fixed pressure in the step 1), the bubble point temperature is selected in the calculation of the condensation process, and the dew point temperature is selected in the calculation of the evaporation process.
The binary non-azeotropic mixed working medium condensation evaporation simulation method considering component distribution is used for two-dimensional model calculation and can also be used for three-dimensional model calculation.
The invention has the beneficial effects that: the invention has wide application range, can obtain the refrigeration performance such as component distribution diagram, heat exchange coefficient and the like of the binary non-azeotropic mixed working medium in the equipment, provides reference for researching the performance of the novel binary non-azeotropic mixed working medium, and saves the research and development cost and period of the novel binary non-azeotropic mixed working medium.
Drawings
FIG. 1 is a flow chart of the present invention.
FIG. 2 is a diagram of "saturation temperature-mass fraction" of binary zeotropic mixture working fluids A-B in the examples.
Detailed Description
The present invention will be described in detail below with reference to the accompanying drawings by taking the condensation and evaporation of a mixture of binary non-azeotropic mixture R1234ze (E) and R32 as an example.
As shown in FIG. 1, the binary non-azeotropic mixed working medium condensation and evaporation simulation method considering component distribution comprises the following steps:
1) according to the type, initial ratio and pressure reference calculation software NIST of the binary non-azeotropic mixed working medium A-B, drawing a 'saturation temperature-mass fraction' diagram of the binary non-azeotropic mixed working medium A-B under a certain fixed pressure, as shown in figure 2, a component A is R1234ze (E), a component B is R32, the left side of an abscissa represents the mass fraction of the component R1234ze (E), the right side represents the mass fraction of the component R32, the ordinate represents the temperature of the binary non-azeotropic mixed working medium, a solid line represents a bubble point temperature, and a dotted line represents a dew point temperature; the initial mixture ratio is 54 percent of R1234ze (E) and 46 percent of R32, and the pressure is 1.74 MPa;
in FIG. 2 point 1 is superheated vapor, point 2 or 2 "is saturated vapor, point 2 'is saturated liquid, point 4 or 4' is saturated liquid, point 4" is saturated gas, and point 5 is subcooled liquid; a point 3 between the points 2 and 4 is in a gas-liquid coexisting state;
2) when the temperature of the binary non-azeotropic mixed working medium is higher than the dew point temperature, the binary non-azeotropic mixed working medium is in a superheated steam state, namely a state of point 1, and the condensation or evaporation rate ratio of A-B of the binary non-azeotropic mixed working medium is set to be the mass fraction ratio (0.54/0.46) of the mixed working medium in the initial filling state; when the temperature of the binary non-azeotropic mixed refrigerant is lower than the bubble point temperature, the binary non-azeotropic mixed refrigerant is in a supercooled liquid state, namely a state of point 5, the condensation or evaporation rate ratio of A-B of the binary non-azeotropic mixed refrigerant is set to be the mass fraction ratio of the mixed refrigerant in the initial filling state, and then the ratio alpha of the condensation amount of the component A and the condensation amount of the component B of the binary non-azeotropic mixed refrigerant to the total condensation amount of the binary non-azeotropic mixed refrigerant is determinedAAnd alphaBRespectively as follows:
αArepresents the ratio of the condensation or evaporation capacity of the component A to the total condensation or evaporation capacity of the binary non-azeotropic mixed working medium A-B, alphaBThe ratio of the condensation or evaporation capacity of the component B to the total condensation or evaporation capacity of the binary non-azeotropic mixed working medium A-B is expressed;
in addition, the saturation temperature of the binary non-azeotropic mixed working medium is set as follows:
Tsatthe saturated temperature of binary non-azeotropic mixed working medium A-B is represented, and when the condensation process is calculated, the component is in the initial filling stateThe bubble point temperature corresponding to the mass fraction of A and B is a fixed value; when the evaporation process is calculated, the dew point temperature corresponding to the mass fractions of the component A and the component B in the initial filling state is a fixed value;
when the temperature of the binary non-azeotropic mixed working medium is lower than the dew point temperature and higher than the bubble point temperature, the binary non-azeotropic mixed working medium is in a gas-liquid coexisting state, namely a state of point 3, and the mass fraction ratio of the component A to the component B is changed all the time; point 3 is an ideal equilibrium state, the mass fraction of point 3 in the liquid phase is the value corresponding to point 3', and the mass fraction in the gas phase is the value corresponding to point 3 ″; however, in the actual iterative calculation process, its mass fraction in the liquid phase is a value corresponding to point 3' and its mass fraction in the gas phase is a point3”The corresponding value is not in an equilibrium state; thus setting the actual state (point) of its two components A and B3’And point3”) The ratio of the difference from the ideal equilibrium state (point 3 'and point 3') is taken as the ratio of the condensation or evaporation rates Δ m of component A and component BA/ΔmB,ΔmAAnd Δ mBIs based on the following formula:
mis the total mass, x, of the binary non-azeotropic mixed working medium3Is the dryness, x, of an ideal equilibrium state 3 Is the dryness of the actual state, w3”Is the mass fraction of component B in the gas phase at the ideal equilibrium state, w 3” Is the mass fraction of component B in the gas phase in the actual state; the ideal equilibrium state refers to the state of the binary non-azeotropic mixed working medium at each temperature in the saturation temperature-mass fraction diagram at the working pressure in the step 1), the bubble point temperature is selected by calculation in the condensation process, and the dew point temperature is selected by calculation in the evaporation process;
the condensation or evaporation capacity of the components A and B of the binary non-azeotropic mixed working medium A-B occupies the ratio alpha of the total condensation or evaporation capacity of the binary non-azeotropic mixed working medium A-BAAnd alphaBRespectively as follows:
ΔmArepresents the difference between the actual mass of component A in the gas phase and the ideal equilibrium mass, Δ mBRepresents the difference between the actual state mass of component B in the gas phase and the ideal equilibrium state mass;
at this time, the saturation temperature of the binary non-azeotropic mixed working medium A-B is set to be
TsatThe saturated temperature of the binary non-azeotropic mixed working medium A-B is represented, and when the condensation process is calculated, the bubble point temperature corresponding to the mass fraction of the components A and B in the actual calculation process is a variable; when calculating the evaporation process, the dew point temperature corresponding to the mass fraction of the components A and B in the actual calculation process is a variable; w is the ratio of the component B to the total binary non-azeotropic mixed working medium A-B in the actual calculation process;
3) obtaining alphaAAnd alphaBThen, combining with finite element simulation software FLUENT, selecting a component transportation model, a VOF two-phase flow model and an SST turbulence model, and setting the binary non-azeotropic mixed working medium A-B and the mass and energy transfer between the component A and the component B thereof according to the following formula in the form of source terms:
if the temperature of the binary non-azeotropic mixed working medium is higher than the saturation temperature, the mass transfer amounts of the component A, B and the binary non-azeotropic mixed working medium gas phase and liquid phase are as follows:
in the formula (I), the compound is shown in the specification,represents the mass transfer of the gas phase of component a,represents the mass transfer amount of the liquid phase of component a,represents the mass transfer of the gas phase of component B,the mass transfer amount of the liquid phase of component B is represented, and r represents the condensation evaporation coefficient, and is usually 0.01-3000000, betalRepresents the volume fraction, rho, of the liquid phase of the binary non-azeotropic mixed working medium A-BlExpressing the density of the liquid phase of the binary non-azeotropic mixed working medium A-B, T expressing the temperature of the binary non-azeotropic mixed working medium A-B, TsatRepresenting the saturation temperature of the binary non-azeotropic mixed working medium A-B;
if the temperature of the binary non-azeotropic mixed working medium is lower than the saturation temperature, the mass transfer amounts of the component A, the component B and the binary non-azeotropic mixed working medium A-B gas phase and liquid phase are as follows:
in the formula, betaνRepresenting the volume fraction, rho, of the binary non-azeotropic mixed working medium A-B gas phaseνRepresenting the density of the gas phase of the binary non-azeotropic mixed working medium A-B;
the energy transfer capacity of the component A, B and the binary non-azeotropic mixed working medium gas phase and liquid phase is as follows:
in the formula, EARepresents the amount of energy transfer of component a,denotes the mass transfer of the liquid phase of component A, γADenotes the latent heat of component A, EBRepresents the amount of energy transfer of component B,denotes the mass transfer of the liquid phase of component B, γBThe latent heat of the component B is expressed, and the E represents the total energy transfer amount of the binary non-azeotropic mixed working medium;
and carrying out calculation to obtain the heat exchange coefficient of the binary non-azeotropic mixed working medium and the component distribution diagram of the component A and the component B in the equipment.
Claims (7)
1. The binary non-azeotropic mixed working medium condensation evaporation simulation method considering component distribution is characterized by comprising the following steps of:
1) according to the type, initial ratio and pressure of the binary non-azeotropic mixed working medium A-B, consulting and calculating software NIST to draw a 'saturation temperature-mass fraction' diagram of the binary non-azeotropic mixed working medium A-B under a certain fixed pressure;
2) when the temperature of the binary non-azeotropic mixed working medium A-B is higher than the dew point temperature or lower than the bubble point temperature, setting the condensation or evaporation rate ratio of the binary non-azeotropic mixed working medium A-B as the mass fraction ratio of the mixed working medium in the initial filling state, wherein the condensation or evaporation amount of the component A and the component B of the binary non-azeotropic mixed working medium A-B occupies the ratio alpha of the total condensation or evaporation amount of the binary non-azeotropic mixed working medium A-BAAnd alphaBRespectively as follows:
αArepresents the ratio of the condensation or evaporation capacity of the component A to the total condensation or evaporation capacity of the binary non-azeotropic mixed working medium A-B, alphaBRepresents the ratio of the condensation or evaporation capacity of the component B to the total condensation or evaporation capacity of the binary non-azeotropic mixed working medium A-B, w0The component B occupies the mass ratio of the total binary non-azeotropic mixed working medium A to B in the initial filling state;
at the moment, the saturation temperature of the binary non-azeotropic mixed working medium A-B is set as follows:
Tsatthe saturated temperature of the binary non-azeotropic mixed working medium A-B is represented, and when the condensation process is calculated, the bubble point temperature corresponding to the mass fraction of the component A and the component B in the initial filling state is a fixed value; when the evaporation process is calculated, the dew point temperature corresponding to the mass fractions of the component A and the component B in the initial filling state is a fixed value; p represents the pressure of the binary non-azeotropic mixed working medium A-B and is set as a constant value in the whole calculation process;
when the temperature of the binary non-azeotropic mixed working medium A-B is lower than the dew point temperature and higher than the bubble point temperature, setting the condensation or evaporation rate ratio of the binary non-azeotropic mixed working medium A-B as DeltamA/ΔmB,ΔmAAnd Δ mBIs based on the following formula:
ΔmArepresents the difference between the actual mass of component A in the gas phase and the ideal equilibrium mass, Δ mBRepresents the difference between the actual mass of component B in the gas phase and the ideal equilibrium mass. m is the total mass of binary non-azeotropic mixture A-B, x3Is the dryness, x, of an ideal equilibrium state 3 Is the dryness of the actual state, w3”Is the mass fraction of component B in the gas phase at the ideal equilibrium state, w 3” Is the mass fraction of component B in the gas phase in the actual state;
the condensation or evaporation capacity of the component A and the component B of the binary non-azeotropic mixed working medium A-B occupies the ratio alpha of the total condensation or evaporation capacity of the binary non-azeotropic mixed working medium A-BAAnd alphaBRespectively as follows:
at the moment, the saturation temperature of the binary non-azeotropic mixed working medium A-B is set as follows:
Tsatthe saturation temperature of the binary non-azeotropic mixed working medium A-B is represented, and when the condensation process is calculated, the bubble point temperature corresponding to the mass fraction of the component A and the component B in the actual calculation process is a variable; when the evaporation process is calculated, the dew point temperature corresponding to the mass fractions of the component A and the component B in the actual calculation process is a variable; w is the ratio of the component B to the total binary non-azeotropic mixed working medium A-B in the actual calculation process;
3) obtaining alphaAAnd alphaBAnd then, combining finite element simulation software FLUENT, selecting a component transport model, a VOF two-phase flow model and an SST turbulence model, setting the binary non-azeotropic mixed working medium A-B and the mass and energy transfer between the component A and the component B thereof according to the following formula in a source item mode, and carrying out calculation to obtain the heat exchange coefficient of the binary non-azeotropic mixed working medium A-B and the component distribution diagram of the component A and the component B in the equipment.
2. The binary non-azeotropic mixed working medium condensation and evaporation simulation method considering component distribution according to claim 1, wherein: if the temperature of the binary non-azeotropic mixed working medium A-B is higher than the saturation temperature, the mass transfer amounts of the component A, the component B and the binary non-azeotropic mixed working medium A-B in gas phase and liquid phase are as follows:
in the formula (I), the compound is shown in the specification,represents the mass transfer of the gas phase of component a,represents the mass transfer amount of the liquid phase of component a,represents the mass transfer of the gas phase of component B,represents the mass transfer of the liquid phase of component B,rthe condensing evaporation coefficient is 0.01-3000000, betalRepresents the volume fraction, rho, of the liquid phase of the binary non-azeotropic mixed working medium A-BlExpressing the density of the liquid phase of the binary non-azeotropic mixed working medium A-B, T expressing the temperature of the binary non-azeotropic mixed working medium A-B, TsatRepresenting the saturation temperature of the binary non-azeotropic mixed working medium A-B.
3. The binary non-azeotropic mixed working medium condensation and evaporation simulation method considering component distribution according to claim 1, wherein: if the temperature of the binary non-azeotropic mixed working medium A-B is lower than the saturation temperature, the mass transfer amounts of the component A, the component B and the binary non-azeotropic mixed working medium A-B in gas phase and liquid phase are as follows:
in the formula, betaνRepresenting the volume fraction, rho, of the binary non-azeotropic mixed working medium A-B gas phaseνThe density of the gas phase of the binary non-azeotropic mixed working medium A-B is shown.
4. The binary non-azeotropic mixed working medium condensation and evaporation simulation method considering component distribution according to claim 1, wherein: the energy transfer amounts of the component A, the component B and the binary non-azeotropic mixed working medium A-B in gas phase and liquid phase are as follows:
in the formula, EARepresents the amount of energy transfer of component a,denotes the mass transfer of the liquid phase of component A, γADenotes the latent heat of component A, EBRepresents the amount of energy transfer of component B,denotes the mass transfer of the liquid phase of component B, γBThe latent heat of the component B is shown, and E represents the total energy transfer of the binary non-azeotropic mixed working medium A-B.
5. The binary non-azeotropic mixed working medium condensation and evaporation simulation method considering component distribution according to claim 1, wherein: the two components of the binary non-azeotropic mixed working medium A-B are refrigeration working media or low-temperature working media, the refrigeration working media are R22, R32, R123, R125, R134a, R152a, R161, R1234yf and R1234ze (E), and the low-temperature working media are methane, ethane, propane and butane.
6. The binary non-azeotropic mixed working medium condensation and evaporation simulation method considering component distribution according to claim 1, wherein: the ideal equilibrium state in the step 2) refers to the state of the binary non-azeotropic mixed working medium A-B at each temperature in the 'saturation temperature-mass fraction' diagram under a certain fixed pressure in the step 1), the bubble point temperature is selected in the calculation of the condensation process, and the dew point temperature is selected in the calculation of the evaporation process.
7. The binary non-azeotropic mixed working medium condensation and evaporation simulation method considering component distribution according to claim 1, wherein: the binary non-azeotropic mixed working medium condensation evaporation simulation method considering component distribution is used for two-dimensional model calculation and can also be used for three-dimensional model calculation.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011153247.9A CN112287579B (en) | 2020-10-26 | 2020-10-26 | Binary non-azeotropic mixed working medium condensation evaporation simulation method considering component distribution |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011153247.9A CN112287579B (en) | 2020-10-26 | 2020-10-26 | Binary non-azeotropic mixed working medium condensation evaporation simulation method considering component distribution |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112287579A true CN112287579A (en) | 2021-01-29 |
CN112287579B CN112287579B (en) | 2022-12-09 |
Family
ID=74424933
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011153247.9A Active CN112287579B (en) | 2020-10-26 | 2020-10-26 | Binary non-azeotropic mixed working medium condensation evaporation simulation method considering component distribution |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112287579B (en) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104946206A (en) * | 2014-03-27 | 2015-09-30 | 北京工业大学 | Binary non-azeotropic mixed refrigerant containing difluoromethane and dichlorotrifluoroethane |
CN108694274A (en) * | 2018-04-18 | 2018-10-23 | 青岛海聚仿真软件技术有限公司 | Condense the Numeral Emulation System and its method of phase transformation secondary development analysis |
-
2020
- 2020-10-26 CN CN202011153247.9A patent/CN112287579B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104946206A (en) * | 2014-03-27 | 2015-09-30 | 北京工业大学 | Binary non-azeotropic mixed refrigerant containing difluoromethane and dichlorotrifluoroethane |
CN108694274A (en) * | 2018-04-18 | 2018-10-23 | 青岛海聚仿真软件技术有限公司 | Condense the Numeral Emulation System and its method of phase transformation secondary development analysis |
Non-Patent Citations (1)
Title |
---|
张永欣等: "非共沸工质R134a/R245fa的管内冷凝换热特性研究", 《工程热物理学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN112287579B (en) | 2022-12-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Lee et al. | Experimental studies of isobutane (R600a) as the refrigerant in domestic refrigeration system | |
Chen et al. | Performance of a new refrigeration cycle using refrigerant mixture R32/R134a for residential air-conditioner applications | |
Corberán et al. | Influence of the source and sink temperatures on the optimal refrigerant charge of a water-to-water heat pump | |
Chen et al. | Experimental study on R-22, R-427A, R-161 and R-290 in air-source heat pump for space heating at low ambient temperatures | |
Deymi-Dashtebayaz et al. | Investigation of the effect of using various HFC refrigerants in geothermal heat pump with residential heating applications | |
Thu et al. | Drop-in experiments and exergy assessment of HFC-32/HFO-1234yf/R744 mixture with GWP below 150 for domestic heat pumps | |
Aprea et al. | Change in energy performance as a result of a R422D retrofit: An experimental analysis for a vapor compression refrigeration plant for a walk-in cooler | |
JPH1163747A (en) | Freezing air conditioner and method for attaining composition of refrigerant in freezing air conditioner | |
Llopis et al. | Experimental evaluation of HCFC-22 replacement by the drop-in fluids HFC-422A and HFC-417B for low temperature refrigeration applications | |
Song et al. | Experimental investigation on a capillary tube based transcritical CO2 heat pump system | |
Fuadi et al. | Chiller performance study with refrigerant R290 | |
Xie et al. | Numerical investigation on thermodynamic performance of CO2-based mixed refrigerants applied in transcritical system | |
Liu et al. | Influence of secondary fluid on the performance of indirect refrigeration system | |
CN112287579B (en) | Binary non-azeotropic mixed working medium condensation evaporation simulation method considering component distribution | |
Li | Optimal analysis of gas cooler and intercooler for two-stage CO2 trans-critical refrigeration system | |
Yıldız et al. | Experimental investigation of exergy, exergoenvironmental and exergoenviroeconomic analysis of the heat pump system | |
Arcaklıoğlu et al. | An algorithmic approach towards finding better refrigerant substitutes of CFCs in terms of the second law of thermodynamics | |
Liu et al. | The occurrence of pinch point and its effects on the performance of high temperature heat pump | |
Haberschill et al. | Performance prediction of a refrigerating machine using R‐407C: the effect of the circulating composition on system performance | |
Pons et al. | Energy analysis of two-phase secondary refrigeration in steady-state operation, Part 1: global optimization and leading parameter | |
Chen et al. | Concentration shift simulation for the mixed refrigerants R-404A, R-32/134a, and R-407C in an air conditioning system | |
Meewisse et al. | Freezing point depression of various ice slurries | |
EP4065660A1 (en) | Refrigerant compositions and use thereof | |
Cui et al. | Pinch point variation in the gas cooler of the CO2 heat pump water heater: An experimental and simulation comparative study | |
Hwang et al. | Comparison of hydrocarbon R-290 and two HFC blends R-404A and R-410A for medium temperature refrigeration applications |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |