CN109033740A - A kind of prediction technique and system of low temperature phase change cold storage materials thermophysics performance - Google Patents

A kind of prediction technique and system of low temperature phase change cold storage materials thermophysics performance Download PDF

Info

Publication number
CN109033740A
CN109033740A CN201811036081.5A CN201811036081A CN109033740A CN 109033740 A CN109033740 A CN 109033740A CN 201811036081 A CN201811036081 A CN 201811036081A CN 109033740 A CN109033740 A CN 109033740A
Authority
CN
China
Prior art keywords
eutectic
phase
temperature
proportion
enthalpy
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
Application number
CN201811036081.5A
Other languages
Chinese (zh)
Other versions
CN109033740B (en
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.)
South China University of Technology SCUT
Original Assignee
South China University of Technology SCUT
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 South China University of Technology SCUT filed Critical South China University of Technology SCUT
Priority to CN201811036081.5A priority Critical patent/CN109033740B/en
Publication of CN109033740A publication Critical patent/CN109033740A/en
Application granted granted Critical
Publication of CN109033740B publication Critical patent/CN109033740B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • Y02E60/14Thermal energy storage

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

The invention discloses a kind of prediction technique of low temperature phase change cold storage materials thermophysics performance, specific steps include: the proportion that (1) selectes all groups of composite phase-change cool storage material;(2) initial a certain group of proportion and initial temperature are chosen under normal pressure, calculate solid phase fugacity and liquid phase fugacity;(3) cycle calculations seek liquid-solid equilibria temperature, make freezing point curve according to the solid liquid phase fugacity under different materials proportion;(4) according to freezing point curve, the reference state easily obtained is chosen as initial state point, referring to the eutectic enthalpy of the calculation method forecast system of virtual route reference state method.The present invention is suitable for the calculating of the polynary low temperature phase change cold storage materials such as binary ternary, realize that the liquid-solid equilibria to low temperature phase change cold storage material carries out high-precision forecast, efficiently solve the problems, such as experimental method inefficiency, a large amount of resource and time cost can be saved, provides some theoretical reference foundations for forecasting research cool storage material performance.

Description

A kind of prediction technique and system of low temperature phase change cold storage materials thermophysics performance
Technical field
The present invention relates to the pre- of composite phase-change material field more particularly to a kind of low temperature phase change cold storage materials thermophysics performance Survey method and system.
Background technique
Low temperature phase change cold storage technology is that the environmental protection and energy saving technology of energy storage is carried out using the latent heat of phase change of phase-change material, benefit Extra cold energy is stored when load is had more than needed with cold-storage device, provides cooling capacity to device for cooling in underload, from And guarantee the technique even running of cold demand and device, it realizes " peak load shifting " of energy, is primarily adapted for use in duty ratio Compared with concentrating, fluctuate bigger occasion etc..But current Cool Storage Technology only has application in building air-conditioning field, and for 123 The low-temperature cold accumulation investigation of materials of~273K range is less, hinder in cold, deep cooling potential temperature Cool Storage Technology development.Cool storage material is The medium of energy storage is the important factor in order in Cool Storage Technology, thus the thermophysical property for studying cool storage material is exploitation Novel cold-storage material, Cool Storage Technology are able to the key point broken through.
Mainly there are experimental method and simulation calculating method for the research method of phase-change material for cold storage performance at present.Experimental method is mesh The more common method of the preceding domestic and international hot physical performance of research phase-change material, but there are defects below for experimental method: first, measurement is a large amount of When the hot physical performance of the composite material being made of different materials, component proportion is difficult to find, and time-consuming;Second, experimental implementation is complicated, It is influenced by feed purity and equipment precision;Third, experiment test used in sample it is few, the hot physical performance that sample measures is often There is difference with massive material in practice;If the phase-change material potential temperature of research is low, low temperature test condition is harsh to be not easy to realize, is needed Expend a large amount of manpower, material resources and financial resources.Compared with experimental method, simulation is less to phase-change material research, due to the mixing rule of component And interaction parameter, cause the simulation in the physical property method to calculate there are precision deficiency, there is the matching with practical application Not high problem is spent, while being also only used for the research of single phase-change material.
Summary of the invention
In view of the deficiencies of the prior art, the present invention provides a kind of prediction sides of low temperature phase change cold storage materials thermophysics performance Method.Method provided by the invention can fast and effeciently predict the thermophysical property of compound cool storage material, can both solve experimental method Time-consuming, susceptible, low temperature test condition are not easy to realize and test the problems such as big with practical application gap, is further promoted and is stored The Efficiency of cold material;It can also be used as and the fast speed ratio of phase-change material for cold storage is selected in practical engineering project, practical application meaning Justice is great.
The purpose of the present invention can be achieved through the following technical solutions:
A kind of prediction technique of low temperature phase change cold storage materials thermophysics performance, specific steps include:
(1) proportion of all groups of composite phase-change cool storage material is selected;
(2) initial a certain group of proportion and initial temperature are chosen under normal pressure, calculate solid phase fugacity and liquid phase fugacity;
(3) cycle calculations seek liquid-solid equilibria temperature, make solidification according to the solid liquid phase fugacity under different materials proportion Point curve;
(4) according to freezing point curve, the reference state easily obtained is chosen as initial state point, referring to virtual route reference state The eutectic enthalpy of the calculation method forecast system of method.
Preferably, in step (2), the solid phase fugacity of phase-change material for cold storage is calculated using virtual route reference state method;It adopts With the liquid phase fugacity of SRK-UNIFAC equation calculation phase-change material for cold storage.
Specifically, in step (2), in order to seek the pure component solid-state fugacity under specific low temperature degree pressure, by virtual The virtual reference state of each component solid-liquid two-phase establishes each component solid-state fugacity and virtual reference state fugacity by thermodynamic cycle Relationship, calculate and seek solid-state fugacity under virtual condition.
Further, the equation of virtual route reference state method is embodied as:
Wherein, fi sIndicate the solid phase fugacity of component i, fi lIndicate the liquid phase fugacity of component i, Δ HfIndicate the phase transformation of reference state Enthalpy, TfIndicate the phase transition temperature of reference state,Indicate that the liquid specific heat of component i holds,Indicate the solid-state specific heat capacity of component i, R Indicate ideal gas constant, T indicates selected temperature, and P indicates normal pressure.
The solid liquid phase fugacity of composite phase-change cool storage material is calculated separately by above-mentioned equation, is done to obtain eutectic status data Prepare.
Further, in the virtual route reference state method for calculating solid phase fugacity, selected reference state point includes But it is not limited to three phase point and standard conditions point.
Specifically, in step (3), the loop calculation are as follows: whether the solid liquid phase fugacity that judgment step (2) obtains It is equal, if solid phase fugacity and liquid phase fugacity are equal, enable initial temperature and pressure constant, repeat step (2), calculates composite phase-change Solid liquid phase fugacity in the case of next group of proportion of cool storage material;If solid phase fugacity and liquid phase fugacity are unequal, temperature is modified, Recalculate solid liquid phase fugacity of the composite phase-change cool storage material in group proportion.Circulation is until all proportion situations are whole Stop when completion.
Specifically, in step (4), the eutectic temperature and eutectic of composite phase-change cool storage material are chosen from freezing point curve Proportion, with reference to virtual reference state Path Method, calculates the eutectic enthalpy under the state, calculation formula are as follows:
Wherein, Δ HmIt indicates under eutectic proportion, eutectic enthalpy of the composite material under eutectic temperature;It indicates Under eutectic proportion, the solid-state enthalpy difference of eutectic temperature and reference state temperature;It indicates under eutectic proportion, reference state Enthalpy of phase change;It indicates under eutectic proportion, the liquid enthalpy difference of eutectic temperature and reference state temperature;It indicates The solid-state enthalpy of reference state point;Indicate the liquid enthalpy of reference state point;It indicates under eutectic temperature The heat of mixing, TmIndicate selected eutectic temperature, TfIndicate the temperature of selected reference state, n indicates proportion group number, xiExpression group Divide the molar fraction of i.
Further, in the calculation method of the composite phase-change cool storage material eutectic enthalpy, selected reference state point packet Include but be not limited to three phase point and standard conditions point.
Specifically, in the present invention, the phase transformation that the low temperature composite phase-change material is related to is mainly phase transformation volume change Small, the manageable solid-liquid phase change of system temperature.
Further, the low-temperature cold accumulation material is selected from alkane, alcohols etc. of the phase transition temperature between 123K~273K One or more of, can be used for shallow cold (223K~273K), in cold (173K~223K) and deep cooling (123K~173K) etc. no The storing cold energy in synthermal section and utilization.
Specifically, the alkane includes but is not limited to C4~C12, paraffin hydro carbons;The alcohols include but is not limited to ethyl alcohol, Butanol, enanthol, ethylene glycol.
Another object of the present invention is to provide a kind of forecasting systems of low temperature phase change cold storage materials thermophysics performance.This hair The system of bright offer proposes on the basis of former UNIFAC model, in virtual route reference state method replacement master mould Antoine equation establishes a kind of novel thermodynamical model combination comprising virtual route reference state method, compound efficiently to predict The thermophysical property of type low temperature phase change cold storage material.
Another object of the present invention can be achieved through the following technical solutions:
A kind of forecasting system of low temperature phase change cold storage materials thermophysics performance, including database module, eutectic temperature and match Than computing module, eutectic enthalpy computing module and interaction interface.
The database module, for supplying the reference state data transmission of collected cool storage material to interaction interface User's selection or correction, the reference state data after correction will be transported to eutectic temperature and proportion computing module and eutectic enthalpy Computing module;
The eutectic temperature and proportion computing module, are basic model with UNIFAC, establish with virtual route reference state method Combined with the novel thermodynamics of SRK-UNIFAC equation, for combine it is set by user include reference state data, initial temperature, match Than data such as group number and step-lengths, iterative calculation obtains the liquid-solid equilibria situation of low temperature phase change cold storage material, by freezing point curve It is output to interaction interface, the data such as eutectic temperature and eutectic proportion are output to eutectic enthalpy computing module;
The eutectic enthalpy computing module, for integrating eutectic temperature and matching the reference state data after data and correction, Eutectic enthalpy is calculated according to eutectic enthalpy calculation formula and exports calculated result to interaction interface;
The interaction interface, for initial temperature set by user, proportion group number and its iteration step length to be input to Eutectic temperature and proportion computing module, and the information for receiving the output of other modules is fed back into user.
The present invention compared to the prior art, have it is below the utility model has the advantages that
High-precision thermodynamical model proposed by the invention and calculation method, by the powerful physical data of simulation softward The thermophysical property of low temperature composite phase-change material is predicted in library, and the precision for obtaining result can compare favourably with experimental method result.It can It solves experimental method time-consuming, susceptible, low temperature test condition are not easy to realize, inefficiency and experiment and practical application gap The problems such as big.Its high efficiency can be to select the fast speed ratio of the phase-change material for cold storage of cold-storage device in engineering and design provides simultaneously Theory applies foundation.
Detailed description of the invention
Fig. 1 is a kind of specific steps schematic diagram of the prediction technique of low temperature phase change cold storage materials thermophysics performance;
Fig. 2 is the path profile for calculating the virtual thermodynamic calculation method of solid phase fugacity;
Fig. 3 is a kind of structural schematic diagram of low temperature phase change cold storage materials thermophysics performance prediction system;
Fig. 4 is the freezing point curve figure of n-heptanol-n-decane binary shallowly cold system;
Fig. 5 is the freezing point curve figure of cold system in n-butanol-n-hexane binary;
Fig. 6 is n-hexane-normal heptane binary deep cooling system freezing point curve figure.
Specific embodiment
Present invention will now be described in further detail with reference to the embodiments and the accompanying drawings, but embodiments of the present invention are unlimited In this.
Embodiment 1
Prediction for binary shallowly cold phase-change material for cold storage thermophysical property.
It is as shown in Figure 1 a kind of flow chart of the prediction technique of low temperature phase change cold storage materials thermophysics performance, comprising:
(1) proportion of all groups of composite phase-change cool storage material is selected;
(2) initial a certain group of proportion and initial temperature are chosen under normal pressure, calculate solid phase fugacity and liquid phase fugacity;
(3) cycle calculations seek liquid-solid equilibria temperature, make solidification according to the solid liquid phase fugacity under different materials proportion Point curve;
(4) according to freezing point curve, the reference state easily obtained is chosen as initial state point, referring to virtual route reference state The eutectic enthalpy of the calculation method forecast system of method.
In step (2), in order to seek the pure component solid-state fugacity under specific low temperature degree pressure, pass through virtual each component The virtual reference state of solid-liquid two-phase establishes the relationship of each component solid-state fugacity and virtual reference state fugacity by thermodynamic cycle, Calculate the solid-state fugacity sought under virtual condition.Calculate path profile such as Fig. 2 institute of the virtual thermodynamic calculation method of solid phase fugacity Show.
For realizing a kind of schematic diagram such as Fig. 3 institute of the system of the prediction technique of low temperature phase change cold storage materials thermophysics performance Show, the system comprises database module, eutectic temperature and proportion computing module, eutectic enthalpy computing module and information exchanges circle Face.
The database module, for supplying the reference state data transmission of collected cool storage material to interaction interface User's selection or correction, the reference state data after correction will be transported to eutectic temperature and proportion computing module and eutectic enthalpy Computing module;
The eutectic temperature and proportion computing module, are basic model with UNIFAC, establish with virtual route reference state method Combined with the novel thermodynamics of SRK-UNIFAC equation, for combine it is set by user include reference state data, initial temperature, match Than data such as group number and step-lengths, iterative calculation obtains the liquid-solid equilibria situation of low temperature phase change cold storage material, by freezing point curve It is output to interaction interface, the data such as eutectic temperature and eutectic proportion are output to eutectic enthalpy computing module;
The eutectic enthalpy computing module, for integrating eutectic temperature and matching the reference state data after data and correction, Eutectic enthalpy is calculated according to eutectic enthalpy calculation formula and exports calculated result to interaction interface;
The interaction interface, for initial temperature set by user, proportion group number and its iteration step length to be input to Eutectic temperature and proportion computing module, and the information for receiving the output of other modules is fed back into user.
According to the above method and system, in the present embodiment by taking n-heptanol-n-decane system as an example, with mass fraction 0.1 For step-length, totally 1001 groups of combinations calculate separately under different ratio proportion by powerful Data Base of Chemical Compound from 0:100 to 100:0 Solid liquid phase fugacity, capture solid liquid phase fugacity it is equal when temperature, obtain freezing point curve, as shown in Figure 4.
By freezing point curve it is found that eutectic proportion is n-heptanol: n-decane=76.7:23.3, eutectic point 234.55K. Using solid, liquid standard molar formation enthalpy as reference state point, two end number mixing heat is obtained with UNIFAC model database, calculates eutectic enthalpy, Obtaining eutectic enthalpy is 152.8kJ/kg.
The calculated result of some binary materials is as shown in table 1 in shallow cold-zone domain.
The hot physical performance computational chart of the shallow cold material of table 1
Cool storage material Congruent melting composition ratio Eutectic point (K) Fusion enthalpy (kJ/kg)
N-heptanol-n-decane 0.77:0.23 234.55 165.1
N-heptanol-n-undecane 0.78:0.22 238.25 152.8
N-decane-n-undecane 0.48:0.52 227.65 138.2
N-decane-n-dodecane 0.78:0.22 239.75 194.1
N-undecane-n-dodecane 0.76:0.24 241.75 155.9
Embodiment 2
The prediction of cold phase-change material for cold storage thermophysical property in binary.Method and system based in the present embodiment and reality It applies consistent in example 1.
Be step-length with mass fraction 0.1 in the present embodiment by taking n-butanol-n-hexane system as an example, proportion from 0:100 to 100:0 totally 1001 groups of combinations calculate separately the solid liquid phase fugacity under different ratio by powerful Data Base of Chemical Compound, capture solid Temperature when liquid phase fugacity is equal obtains freezing point curve, as shown in Figure 5.
By freezing point curve it is found that eutectic proportion is n-butanol: n-hexane=59.8:40.2, eutectic point 175.05K. Using solid, liquid standard molar formation enthalpy as reference state point, two end number mixing heat is obtained with UNIFAC model database, calculates eutectic enthalpy, Obtaining eutectic enthalpy is 132.1kJ/kg.
Some binary material calculated results are as shown in table 2 in middle cold-zone domain.
The hot physical performance computational chart of cold material in table 2
Embodiment 3
The prediction of binary deep cooling phase-change material for cold storage thermophysical property.Method and system based in the present embodiment and reality It applies consistent in example 1.
The present embodiment is step-length with mass fraction 0.1 by taking n-hexane-n-heptane system as an example, is matched from 0:100 to 100: 0 totally 1001 groups of combinations calculate separately the solid liquid phase fugacity under different ratio by powerful Data Base of Chemical Compound, capture solid liquid phase Temperature when fugacity is equal obtains freezing point curve, as shown in Figure 6.
By freezing point curve it is found that eutectic proportion is n-hexane: normal heptane=56.0:44.0, eutectic point 165.65K. Using solid, liquid standard molar formation enthalpy as reference state point, two end number mixing heat is obtained with UNIFAC model database, calculates eutectic enthalpy, Obtaining eutectic enthalpy is 141.5kJ/kg.
Some binary material calculated results are as shown in table 3 in deep cooling region.
The hot physical performance computational chart of 3 deep cooling material of table
Cool storage material Congruent melting composition ratio Eutectic point (K) Fusion enthalpy (kJ/kg)
N-hexane-normal heptane 0.56:0.44 165.65 141.5
Normal butane-pentane 0.62:0.38 122.75 82.7
Normal butane-n-hexane 0.91:0.09 132.65 85.0
Pentane-n-hexane 0.86:0.14 140.95 118.2
Pentane-normal heptane 0.89:0.11 141.75 116.1
Embodiment 4
Cool storage material is adapted to Journal of Sex Research with LNG cold energy air separation unit.Method and system based in the present embodiment with It is consistent in embodiment 1.
Liquefied natural gas (LNG) is the cryogenic liquid of 111K a kind of under normal pressure, needed before being supplied to downstream user by It gasifies and is heated to 273K or more, a large amount of cold energy can be released in gasification, 1.0tLNG about can release about when gasifying The cold energy of 200kWh gets up this part cold storage using to economizing on resources and improves energy efficiency with highly important Meaning.
But defeated load frequently changes with the dynamic with wave of downstream user outside LNG, defeated load is relatively large outside daytime, LNG gasification amount is smaller when night gas low ebb, or even stops gasification, causes to require continuous, stable operation LNG cold energy use dress It sets and is unable to stable operation because of cold energy supply discontinuity.Therefore, the present embodiment is carried out by cold-storage device to cooling supply and with device for cooling It is adapted to Journal of Sex Research, briefly explains the scope of application of the method for the present invention.
Cold energy air is the process for making air liquefaction in liquefaction process constituent of air to separate using cold energy.Certain cold energy The production scale of air separation unit is as shown in table 4.It provides the LNG of cooling capacity for air separation unit to provide for periphery receiving station, principal states Parameter is as follows: pressure 66.0barg, temperature 128.15K;LNG moles of composition are as follows: CH4: 96.30%;C2H6: 2.58%;C3H8: 0.49%;N-C4H10: 0.12%;I-C4H10: 0.10%;N-C5H12:0.003;I-C5H12:0.003;N2: 0.40%.Normal work Under condition, the air mass flow of air separation unit charging is 51950Nm3/ h needs the LNG of about 50.0t/h to provide gasification cooling capacity and meets The normal operation of cold energy air separation unit.
Air separation unit scale under 4 normal design operating condition of table
Under normal circumstances, the produce load adjustable range of air separation unit is 75%~105%.Under nominal situation, LNG is utilized Air separation unit LNG inlet temperature be 128K, when cooling capacity deficiency, air separation unit is with minimum 75% load operation, at this time most Low acceptable LNG inlet temperature is 153K.Therefore, for maintaining the potential temperature of the phase-change material of LNG cold energy space division to should be at 128K Between~153K.The thermal physical characteristic parameter of embodiment binary cold-storage phase-change material calculated before comprehensive is considering safety, warp Ji property etc. under the premise of, select cocrystallization point for 140.95K pentane-n-hexane system as the device cool storage material more Properly.
Outside gasifying defeated trough period typical day in addition for defeated load, since demand of the air separation unit to LNG is 50.0t/ H, comprehensively considers function of environment heat emission, cool storage heat exchanger heat transfer temperature difference and heat exchange efficiency etc., and the LNG cooling capacity which needs is 9115kW。
The same day lower than 50.0t/h gas supply underestimate the time be 5 hours, cold-storage device cooling supply is needed, because of the coolness storing material of use Material eutectic point is 140.95K, so the cooling capacity that the cooling capacity of storage discharges during being warming up to 140.95K by the LNG of 123K. The cold energy that pentane per ton-n-hexane cool storage material can store is 11.8 × 104KJ, so the demand of cool storage material is 278t/h, aggregate demand 1390t.LNG amount needed for cold-storage is 2477t.In the daytime gasification load is small up to 19 higher than 50.0t/h When, total amount of vaporization is 7614.8t, 2477t is much larger than, so the enough night gas supply of LNG gasification cooling capacity that receiving station has more than needed in the daytime The refrigeration requirement of low ebb.Cool storage material dosage can guarantee completely to liquefy the natural gas refrigerant of 50t/h when being 1390t, and will Cooling capacity brings air separation unit.
Defeated trough period cool storage material is adapted to table with LNG cold energy air separation unit outside table 5
In the present embodiment, it is using the energy consumption of the process production unit liquid yield of LNG cool storage material cold-storage 0.41kWh/m3, with traditional process about 1.00kWh/m3Specific energy consumption compared to there is apparent power savings advantages.
To sum up four embodiments, the novel thermodynamical model combination of this patent exploitation are right using powerful Data Base of Chemical Compound Preferable result can be obtained by calculating hydro carbons and alcohols thermophysical property.Stepping of the present invention can overcome the various of experimental method to ask Topic, additionally it is possible to provide theory using foundation for the selection and design of the cool storage material of cold-storage device.
The above embodiment is a preferred embodiment of the present invention, but embodiments of the present invention are not by above-described embodiment Limitation, other any changes, modifications, substitutions, combinations, simplifications made without departing from the spirit and principles of the present invention, It should be equivalent substitute mode, be included within the scope of the present invention.

Claims (10)

1. a kind of prediction technique of low temperature phase change cold storage materials thermophysics performance, which is characterized in that specific steps include:
(1) proportion of all groups of composite phase-change cool storage material is selected;
(2) initial a certain group of proportion and initial temperature are chosen under normal pressure, calculate solid phase fugacity and liquid phase fugacity;
(3) cycle calculations seek liquid-solid equilibria temperature according to the solid liquid phase fugacity under different materials proportion, make freezing point song Line;
(4) according to freezing point curve, the reference state easily obtained is chosen as initial state point, referring to virtual route reference state method The eutectic enthalpy of calculation method forecast system.
2. a kind of prediction technique of low temperature phase change cold storage materials thermophysics performance according to claim 1, which is characterized in that In step (2), the solid phase fugacity of phase-change material for cold storage is calculated using virtual route reference state method;Using SRK-UNIFAC equation Calculate the liquid phase fugacity of phase-change material for cold storage.
3. a kind of prediction technique of low temperature phase change cold storage materials thermophysics performance according to claim 2, which is characterized in that The equation of virtual route reference state method is embodied as:
Wherein, fi sIndicate the solid phase fugacity of component i, fi lIndicate the liquid phase fugacity of component i, Δ HfIndicate the enthalpy of phase change of reference state, TfIndicate the phase transition temperature of reference state,Indicate that the liquid specific heat of component i holds,Indicate the solid-state specific heat capacity of component i, R table Show ideal gas constant, T indicates selected temperature, and P indicates normal pressure.
4. a kind of prediction technique of low temperature phase change cold storage materials thermophysics performance according to claim 1, which is characterized in that In the virtual route reference state method of the calculating solid phase fugacity of step (2), selected reference state point includes but is not limited to three phase point With standard conditions point.
5. a kind of prediction technique of low temperature phase change cold storage materials thermophysics performance according to claim 1, which is characterized in that In step (3), the loop calculation are as follows: whether the solid liquid phase fugacity that judgment step (2) obtains is equal, if solid phase fugacity It is equal with liquid phase fugacity, then it enables initial temperature constant, repeats step (2), calculate next group of proportion of composite phase-change cool storage material In the case of solid liquid phase fugacity;If solid phase fugacity and liquid phase fugacity are unequal, temperature is modified, recalculates composite phase-change cold-storage Solid liquid phase fugacity of the material when the group matches.Circulation is until stopping when all proportion situations are fully completed.
6. a kind of prediction technique of low temperature phase change cold storage materials thermophysics performance according to claim 1, which is characterized in that In step (4), the eutectic temperature and eutectic proportion of composite phase-change cool storage material are chosen from freezing point curve, according to virtual ginseng State path is examined, the eutectic enthalpy under the state, calculation formula are calculated are as follows:
Wherein, Δ HmIt indicates under eutectic proportion, eutectic enthalpy of the composite material under eutectic temperature;It indicates in eutectic Under proportion, the solid-state enthalpy difference of eutectic temperature and reference state temperature;It indicates under eutectic proportion, the phase transformation of reference state Enthalpy;It indicates under eutectic proportion, the liquid enthalpy difference of eutectic temperature and reference state temperature;Indicate reference The solid-state enthalpy of state point;Indicate the liquid enthalpy of reference state point;It indicates under eutectic temperature The heat of mixing, TmIndicate selected eutectic temperature, TfIndicate the temperature of selected reference state,Indicate that the liquid specific heat of component i holds,Indicate the solid-state specific heat capacity of component i, n indicates proportion group number, xiIndicate the molar fraction of component i.
7. a kind of prediction technique of low temperature phase change cold storage materials thermophysics performance according to claim 1, which is characterized in that In the calculation method of the composite phase-change cool storage material eutectic enthalpy of step (4), selected reference state point includes but is not limited to three-phase Point and standard conditions point.
8. a kind of prediction technique of low temperature phase change cold storage materials thermophysics performance according to claim 1, which is characterized in that The low-temperature cold accumulation material is selected from one or more of alkane, the alcohols of phase transition temperature between 123K~273K etc., can use Yu Qianleng (223K~273K), in the different temperatures section such as cold (173K~223K) and deep cooling (123K~173K) storing cold energy And utilization.
9. a kind of prediction technique of low temperature phase change cold storage materials thermophysics performance according to claim 8, which is characterized in that The alkane includes but is not limited to C4~C12, paraffin hydro carbons;The alcohols includes but is not limited to ethyl alcohol, butanol, enanthol, second two Alcohol.
10. a kind of forecasting system of the low temperature phase change cold storage materials thermophysics performance for realizing claim 1-9, including data Library module, eutectic temperature and proportion computing module, eutectic enthalpy computing module and interaction interface;It is characterized in that,
The database module, for by the reference state data transmission of collected cool storage material to interaction interface for user Selection or correction, the reference state data after correction will be transported to eutectic temperature and proportion computing module and eutectic enthalpy calculates Module;
The eutectic temperature and proportion computing module, are basic model with UNIFAC, establish with virtual route reference state method and The novel thermodynamics of SRK-UNIFAC equation combines, and includes reference state data, initial temperature, proportion for combining set by user The group data such as number and step-length, iterative calculation obtains the liquid-solid equilibria situation of low temperature phase change cold storage material, and freezing point curve is defeated Interaction interface is arrived out, and the data such as eutectic temperature and eutectic proportion are output to eutectic enthalpy computing module;
The eutectic enthalpy computing module, for integrating eutectic temperature and matching the reference state data after data and correction, according to Eutectic enthalpy calculation formula calculates eutectic enthalpy and exports calculated result to interaction interface;
The interaction interface, for initial temperature set by user, proportion group number and its iteration step length to be input to eutectic Temperature and proportion computing module, and the information for receiving the output of other modules is fed back into user.
CN201811036081.5A 2018-09-06 2018-09-06 Method and system for predicting thermal physical properties of low-temperature phase-change cold storage material Active CN109033740B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811036081.5A CN109033740B (en) 2018-09-06 2018-09-06 Method and system for predicting thermal physical properties of low-temperature phase-change cold storage material

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811036081.5A CN109033740B (en) 2018-09-06 2018-09-06 Method and system for predicting thermal physical properties of low-temperature phase-change cold storage material

Publications (2)

Publication Number Publication Date
CN109033740A true CN109033740A (en) 2018-12-18
CN109033740B CN109033740B (en) 2023-09-29

Family

ID=64624327

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811036081.5A Active CN109033740B (en) 2018-09-06 2018-09-06 Method and system for predicting thermal physical properties of low-temperature phase-change cold storage material

Country Status (1)

Country Link
CN (1) CN109033740B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113420428A (en) * 2021-06-09 2021-09-21 中国科学院理化技术研究所 Solid phase cold accumulation medium screening method and design method of liquid air energy storage power station
CN115558469A (en) * 2022-10-10 2023-01-03 青岛中科未来健康研究院有限公司 Mixed phase change cold accumulation medium

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102495968A (en) * 2011-12-14 2012-06-13 中国科学院工程热物理研究所 Numerical simulation method for solid/liquid phase change in high-temperature heat storage container in gravity condition
CN102779217A (en) * 2012-08-06 2012-11-14 大连三洋压缩机有限公司 Computer simulation performance computation method of refrigeration system under frosting working condition
CN103063320A (en) * 2012-12-11 2013-04-24 西南石油大学 Determination method of throttle temperature and vapor-liquid ratio of high-temperature and high-pressure complex natural gas

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102495968A (en) * 2011-12-14 2012-06-13 中国科学院工程热物理研究所 Numerical simulation method for solid/liquid phase change in high-temperature heat storage container in gravity condition
CN102779217A (en) * 2012-08-06 2012-11-14 大连三洋压缩机有限公司 Computer simulation performance computation method of refrigeration system under frosting working condition
CN103063320A (en) * 2012-12-11 2013-04-24 西南石油大学 Determination method of throttle temperature and vapor-liquid ratio of high-temperature and high-pressure complex natural gas

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113420428A (en) * 2021-06-09 2021-09-21 中国科学院理化技术研究所 Solid phase cold accumulation medium screening method and design method of liquid air energy storage power station
CN113420428B (en) * 2021-06-09 2022-04-22 中国科学院理化技术研究所 Solid phase cold accumulation medium screening method and design method of liquid air energy storage power station
CN115558469A (en) * 2022-10-10 2023-01-03 青岛中科未来健康研究院有限公司 Mixed phase change cold accumulation medium

Also Published As

Publication number Publication date
CN109033740B (en) 2023-09-29

Similar Documents

Publication Publication Date Title
She et al. Enhancement of round trip efficiency of liquid air energy storage through effective utilization of heat of compression
Zhou et al. A review of thermal energy storage in compressed air energy storage system
Noël et al. Phase change materials
CN105899766B (en) For storing/producing the therrmodynamic system of electric energy
Fu et al. Thermodynamic of a novel solar heat storage compressed carbon dioxide energy storage system
Kaşka et al. The performance assessment of a combined organic Rankine-vapor compression refrigeration cycle aided hydrogen liquefaction
Zhang et al. Technical and economic analysis of Brayton-cycle-based pumped thermal electricity storage systems with direct and indirect thermal energy storage
Dincer et al. Progress in Sustainable Energy Technologies Vol II
Walsh et al. Free-cooling thermal energy storage using phase change materials in an evaporative cooling system
Fan et al. Thermo-economic analysis of the integrated system of thermal power plant and liquid air energy storage
Ryu et al. A novel liquid air energy storage system using a combination of sensible and latent heat storage
Shen et al. Comprehensive review on cascaded latent heat storage technology: Recent advances and challenges
Lu et al. Design and thermodynamic analysis of an advanced liquid air energy storage system coupled with LNG cold energy, ORCs and natural resources
Chen et al. Construction and optimization of the cold storage process based on phase change materials used for liquid air energy storage system
Ouyang et al. Multi-energy flow cooperative dispatch for supply-demand balance of distributed power grid with liquid air energy storage system
Ouyang et al. Flexible dispatch strategy of purchasing-selling electricity for coal-fired power plant based on compressed air energy storage
CN109033740A (en) A kind of prediction technique and system of low temperature phase change cold storage materials thermophysics performance
She et al. Preliminary study of Liquid Air Energy Storage integrated with LNG cold recovery
CN204923571U (en) Thermoelectric cold trigeminy supplies system based on heat pump technique and compressed air electric power storage technique
Zhang et al. Heat-storage performance optimization for packed bed using cascaded PCMs capsules
Zhou et al. Performance study on a new solar aided liquid air energy storage system integrated with organic Rankine cycle and thermoelectric generator
Deng et al. Controlling effect of phase change material based heat exchanger on supercritical CO2 Brayton cycle
Lin et al. Thermodynamic analysis and operation investigation of a cross-border integrated energy system based on steam Carnot battery
Wang et al. Review on operation control of cold thermal energy storage in cooling systems
Bi et al. Exergy analysis of a gas-hydrate cool storage system

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