CN108134156B - Parameter calculation method, system, medium, terminal and battery pack for flowing heat exchange of refrigerant in pipe - Google Patents
Parameter calculation method, system, medium, terminal and battery pack for flowing heat exchange of refrigerant in pipe Download PDFInfo
- Publication number
- CN108134156B CN108134156B CN201711216003.9A CN201711216003A CN108134156B CN 108134156 B CN108134156 B CN 108134156B CN 201711216003 A CN201711216003 A CN 201711216003A CN 108134156 B CN108134156 B CN 108134156B
- Authority
- CN
- China
- Prior art keywords
- vector
- node
- calculating
- heat exchange
- refrigerant
- 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.)
- Active
Links
- 239000003507 refrigerant Substances 0.000 title claims abstract description 67
- 238000004364 calculation method Methods 0.000 title claims abstract description 55
- 239000012071 phase Substances 0.000 claims abstract description 50
- 239000011555 saturated liquid Substances 0.000 claims abstract description 39
- 238000000034 method Methods 0.000 claims description 23
- 238000012546 transfer Methods 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 9
- 229920006395 saturated elastomer Polymers 0.000 claims description 6
- 239000007791 liquid phase Substances 0.000 claims description 4
- 238000012545 processing Methods 0.000 description 12
- 238000001816 cooling Methods 0.000 description 11
- 230000008859 change Effects 0.000 description 7
- 230000008569 process Effects 0.000 description 7
- 230000000704 physical effect Effects 0.000 description 6
- 239000007788 liquid Substances 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000005514 two-phase flow Effects 0.000 description 3
- 238000004378 air conditioning Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000005494 condensation Effects 0.000 description 1
- 238000009833 condensation Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 235000019800 disodium phosphate Nutrition 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000008020 evaporation Effects 0.000 description 1
- 238000001704 evaporation Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- BWHOZHOGCMHOBV-BQYQJAHWSA-N trans-benzylideneacetone Chemical compound CC(=O)\C=C\C1=CC=CC=C1 BWHOZHOGCMHOBV-BQYQJAHWSA-N 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01M—PROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
- H01M10/00—Secondary cells; Manufacture thereof
- H01M10/60—Heating or cooling; Temperature control
- H01M10/61—Types of temperature control
- H01M10/613—Cooling or keeping cold
-
- 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
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01M—PROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
- H01M10/00—Secondary cells; Manufacture thereof
- H01M10/60—Heating or cooling; Temperature control
- H01M10/62—Heating or cooling; Temperature control specially adapted for specific applications
- H01M10/625—Vehicles
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01M—PROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
- H01M10/00—Secondary cells; Manufacture thereof
- H01M10/60—Heating or cooling; Temperature control
- H01M10/63—Control systems
- H01M10/633—Control systems characterised by algorithms, flow charts, software details or the like
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01M—PROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
- H01M10/00—Secondary cells; Manufacture thereof
- H01M10/60—Heating or cooling; Temperature control
- H01M10/65—Means for temperature control structurally associated with the cells
- H01M10/655—Solid structures for heat exchange or heat conduction
- H01M10/6554—Rods or plates
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01M—PROCESSES OR MEANS, e.g. BATTERIES, FOR THE DIRECT CONVERSION OF CHEMICAL ENERGY INTO ELECTRICAL ENERGY
- H01M10/00—Secondary cells; Manufacture thereof
- H01M10/60—Heating or cooling; Temperature control
- H01M10/65—Means for temperature control structurally associated with the cells
- H01M10/656—Means for temperature control structurally associated with the cells characterised by the type of heat-exchange fluid
- H01M10/6567—Liquids
-
- 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
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
- Y02E60/10—Energy storage using batteries
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Chemical & Material Sciences (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Electrochemistry (AREA)
- General Chemical & Material Sciences (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Secondary Cells (AREA)
Abstract
The invention provides a parameter calculation method, a system, a medium, a terminal and a battery pack for in-tube refrigerant flow heat exchange, which comprises the steps of setting pressure in a state equation of a thermodynamic system as an independent variable; calculating a saturated liquid phase prandtl number vector; calculating a heat exchange coefficient vector of each node according to the saturated liquid phase Prandtl number vector; based on the energy relation, calculating a first relation curve of a coordinate vector and a saturation pressure vector of a certain node according to a heat exchange coefficient vector of the node and heat flow loaded on the node; calculating a second relation curve of the coordinate vector and the saturation pressure vector of the node based on the pressure relation; determining the coordinate parameter of the node according to the intersection point of the first relation curve and the second relation curve; and calculating the heat exchange coefficient of each node according to the coordinate parameters of the node. The parameter calculation method, the system, the medium, the terminal and the battery pack for the flowing heat exchange of the refrigerant in the pipe can quickly and accurately calculate the parameter of the flowing heat exchange of the refrigerant in the pipe.
Description
Technical Field
The invention relates to the technical field of data processing, in particular to a parameter calculation method, a system, a medium, a terminal and a battery pack for flow heat exchange of a refrigerant in a pipe.
Background
The electric automobile is a vehicle which takes a vehicle-mounted power supply as power and drives wheels to run by using a motor, and meets various requirements of road traffic and safety regulations. Because the influence on the environment is smaller than that of the traditional automobile, the prospect is widely seen. Under the influence of the driving habits of the traditional automobiles, the requirement on the long driving range of the electric automobile gradually becomes the mainstream trend, for example, the highest pure electric driving range of a Tesla Model-S long driving range version is nearly 500 km, and people also claim that the electric automobile exceeding 400km is proposed in the future.
The electric vehicle provides kinetic energy through a Battery Pack (Battery Pack) arranged on the electric vehicle. With the increasing energy density required for battery packs, some manufacturers, particularly european automotive manufacturers, have used refrigerants as cooling media to directly cool power battery systems, such as BMW i3 and BENZ S400. This cooling is known as direct cooling of the battery pack. Compared with liquid cooling, direct cooling needs a cooling working medium with smaller mass to realize the cooling of the battery system, and meanwhile, the requirement of high energy density of the battery pack can be met. With the mature technology of the electric compressor for the automobile, the cost advantage of the direct cooling mode becomes more and more obvious.
Therefore, for the power battery of the electric vehicle, the two-phase flow phase change and cooling technology becomes the current hot research topic. The core processes involved in the heat exchange and flow of two-phase flow are the condensation and evaporation processes. The calculation of the heat exchange of the two-phase flow is complicated because the conversion process and the convection process of the gas-liquid two states are carried out simultaneously. In the prior art, the formulas for calculating heat exchange and pressure change are all one-dimensional empirical formulas, and all the formulas relate to a plurality of intermediate parameters and physical properties, such as dryness, cavitation coefficient, slip ratio, Martin's number, Lewy number, Reynolds number, buoyancy lift coefficient, Prandtl number, friction coefficient and the like. These intermediate parameters and properties are strongly coupled to the heat transfer coefficient, pressure and temperature to be solved.
At present, in the process of numerical simulation of two-phase regions, the physical properties (such as refrigerant temperature) of a plurality of refrigerants in two-phase states are assumed to be unchanged, such as numerical simulation and analysis of air-conditioning condensers of Sushuyu, Liuwei, Xujing and the like, academic newspaper of China science and technology (Nature science edition), 2007,35(9): 42-44; further examples are Garcia-Cascales and Josef Ramon "A model for the analysis of compact heat exchanges", Purdue University of USA: proceedings of the 12th reference and Air Conditioning Conference, 2008. The above processing manner is sufficient for calculating the air conditioner. However, this approach is insufficient for the cold plate of the battery pack because the refrigerant flowing in the cold plate of the battery is likely to have a temperature difference of more than 2degC due to the change in pressure drop.
Disclosure of Invention
In view of the above-mentioned drawbacks of the prior art, an object of the present invention is to provide a method, a system, a medium, a terminal, and a battery pack for calculating a parameter of heat exchange of refrigerant flowing in a tube, which can calculate a parameter of heat exchange of refrigerant flowing in a tube quickly and accurately without assuming a constant for any change in the physical properties and intermediate parameters of the refrigerant along the flowing direction.
In order to achieve the above and other related objects, the present invention provides a method for calculating parameters of refrigerant flowing heat exchange in a tube, comprising the steps of: setting the pressure in a state equation of the thermodynamic system as an independent variable; calculating a saturated liquid phase prandtl number vector; dividing the tube into N equal parts according to equal enthalpy difference by adopting N-1 nodes, and calculating the heat exchange coefficient vector of each node according to the saturated liquid phase Prandtl number vector; based on the energy relation, calculating a first relation curve of a coordinate vector and a saturation pressure vector of a certain node according to a heat exchange coefficient vector of the node and heat flow loaded on the node; calculating a second relation curve of the coordinate vector and the saturation pressure vector of the node based on the pressure relation; determining the coordinate parameter of the node according to the intersection point of the first relation curve and the second relation curve; and calculating the heat exchange coefficient of each node according to the coordinate parameters of the node.
In an embodiment of the present invention, the equation of state of the thermal system is f (p, v, T) ═ 0, where p is pressure, v is volume, and T is temperature.
In one embodiment of the present invention, the saturated liquidus prandtl number vector is calculated according to Pry Ndy/Dry; wherein, Pry is a saturated liquid phase Prandtl number vector, Cpy is a saturated liquid phase constant pressure specific heat capacity vector, Dry is a saturated liquid phase heat conductivity coefficient vector, and Ndy is a saturated liquid phase dynamic viscosity vector.
In one embodiment of the present invention, the heat transfer coefficient vector of the node is
Wherein Pry is the Prandtl number vector of the saturated liquid phase, Dry is the heat conductivity coefficient vector of the saturated liquid phase, and Reynolds number vector ReEQ＝GEQ.Dr./NdlEquivalent mass flow vector GEQ＝G((1-X)+X.*(Mdy./MdgΛ (0.5)), dry-weight vector X ═ hi./Qr-Hl./Qr；hiIs the specific enthalpy value of the ith node, HlIs a specific enthalpy vector, Qr is a flow vector, G is a mass flow vector, MdyIs the refrigerant liquid phase density vector MdgIs the refrigerant gas phase density vector, DrIs hydraulic diameter, NdlIs the viscosity vector.
In an embodiment of the present invention, based on the energy relationship, calculating a first relationship curve between a coordinate vector of a node and a saturation pressure vector according to a heat transfer coefficient vector of the node and a heat flow loaded on the node includes the following steps:
calculating the length vector from the previous node to a certain node according to the heat exchange coefficient vector of the node and the heat flow loaded on the node;
calculating the coordinate parameter of the node according to the coordinate parameter of the previous node and the length vector;
and calculating a first relation curve of the coordinate vector of the node and the saturation pressure vector based on the energy relation.
In an embodiment of the present invention, when calculating the second relationship curve between the coordinate vector of the node and the saturation pressure vector based on the pressure relationship, the second relationship curve is obtained according toCalculating a pressure vector at the node, wherein DrIs hydraulic diameter, /)i-1Is the coordinate parameter, rho, of the i-1 th nodei-1The refrigerant density of the i-1 th node, Frp the average friction coefficient vector, Mdp the average density vector, L "the next node coordinate vector, and Md the current node density vector; g is the mass flow vector, DrIs the hydraulic diameter.
Correspondingly, the invention provides a parameter calculation system for in-tube refrigerant flowing heat exchange, which comprises a setting module, a first calculation module, a second calculation module, a third calculation module, a fourth calculation module, a determination module and a fifth calculation module, wherein the setting module is used for setting the parameter of the in-tube refrigerant flowing heat exchange;
the setting module is used for setting the pressure in a state equation of the thermodynamic system as an independent variable;
the first calculation module is used for calculating a saturated liquid phase Prandtl number vector;
the second calculation module is used for dividing the tube into N equal parts according to the equal enthalpy difference by adopting N-1 nodes, and calculating the heat exchange coefficient vector of each node according to the saturated liquid phase Plantt number vector;
the third calculation module is used for calculating a first relation curve of a coordinate vector and a saturation pressure vector of a certain node according to the heat exchange coefficient vector of the node and the heat flow loaded on the node based on the energy relation;
the fourth calculation module is used for calculating a second relation curve of the coordinate vector of the node and the saturated pressure vector based on the pressure relation;
the determining module is used for determining the coordinate parameter of the node according to the intersection point of the first relation curve and the second relation curve;
and the fifth calculation module is used for calculating the heat exchange coefficient at each node according to the coordinate parameters of the node.
The invention provides a storage medium, which stores a computer program, and is characterized in that the program is executed by a processor to realize the parameter calculation method for the flowing heat exchange of the refrigerant in the pipe.
The invention provides a terminal, which comprises a processor and a memory;
the memory is used for storing a computer program;
the processor is used for executing the computer program stored in the memory so as to enable the terminal to execute the parameter calculation method for the flowing heat exchange of the refrigerant in the pipe.
Finally, the invention provides a battery pack which comprises a temperature-equalizing cold plate, wherein a refrigerant in a pipe is arranged in the temperature-equalizing cold plate, and the refrigerant in the pipe sets flowing heat exchange parameters according to the parameter calculation method for flowing heat exchange of the refrigerant in the pipe.
As mentioned above, the parameter calculation method, system, medium, terminal and battery pack for the flowing heat exchange of the refrigerant in the pipe have the following beneficial effects:
(1) the constant assumption is not carried out on the change of the physical property and the intermediate parameter of any related refrigerant along the flowing direction, and the flowing heat exchange parameter of the refrigerant in the pipe can be quickly and accurately calculated;
(2) iteration and interpolation in the numerical calculation process are avoided or reduced, and the calculation time and the calculation precision are greatly improved;
(3) the normal realization of the flowing heat exchange of the refrigerant in the refrigerating pipe is ensured, so that the liquid cooling plate component has simple, light and compact structure, excellent heat transfer performance, lower flow resistance and low cost;
(4) the flowing heat exchange requirement of the battery pack is met.
Drawings
FIG. 1 is a flow chart illustrating a method of calculating heat transfer parameters for refrigerant flow in a tube according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a parameter calculation system for heat exchange of refrigerant flowing in a tube according to an embodiment of the present invention;
fig. 3 is a schematic structural diagram of a terminal according to an embodiment of the invention;
fig. 4 is a schematic structural diagram of a battery pack according to an embodiment of the invention.
Description of the element reference numerals
21 setting module
22 first computing Module
23 second calculation Module
24 third computing Module
25 fourth calculation Module
26 determining module
27 fifth calculation Module
31 processor
32 memory
41 warm and cold plate
411 intraductal refrigerant
Detailed Description
The embodiments of the present invention are described below with reference to specific embodiments, and other advantages and effects of the present invention will be easily understood by those skilled in the art from the disclosure of the present specification. The invention is capable of other and different embodiments and of being practiced or of being carried out in various ways, and its several details are capable of modification in various respects, all without departing from the spirit and scope of the present invention.
It should be noted that the drawings provided in the present embodiment are only for illustrating the basic idea of the present invention, and the components related to the present invention are only shown in the drawings rather than drawn according to the number, shape and size of the components in actual implementation, and the type, quantity and proportion of the components in actual implementation may be changed freely, and the layout of the components may be more complicated.
The parameter calculation method, the system, the medium, the terminal and the battery pack for the flowing heat exchange of the refrigerant in the pipe do not perform constant assumption on the change of any related refrigerant physical property and intermediate parameter along the flowing direction, can quickly and accurately calculate the flowing heat exchange parameter of the refrigerant in the pipe, and completely meet the flowing heat exchange requirement of the battery pack.
As shown in fig. 1, in an embodiment, the method for calculating the parameter of the refrigerant flowing heat exchange in the tube of the present invention includes the following steps:
step S1, setting the pressure in the state equation of the thermodynamic system as an independent variable.
In describing the strength parameters in a constant-mass thermodynamic system, only a certain number of parameters are independent state parameters, the number of independent state parameters is called the degree of freedom of the system, in a thermodynamic equilibrium state, the number N of independent state parameters describing the system is α - β +2, wherein α is the component and β is the phase state number, therefore, for a single-component constant-mass thermodynamic system, there is only one independent variable in two phase states (α is 1 and β is 2).
Any one of three parameters p, v and T of a thermodynamic system with a state equation f (p, v and T) equal to 0 serves as an independent variable of the system, and all other physical parameters can be expressed as functions of the independent variable. Where p is pressure, v is volume, and T is temperature. In the present invention, p is chosen as an independent variable.
And step S2, calculating a saturated liquid phase Prandtl number vector.
In one embodiment of the present invention, the saturated liquidus prandtl number vector is calculated according to Pry Ndy/Dry; wherein, Pry is a saturated liquid phase Prandtl number vector, Cpy is a saturated liquid phase constant pressure specific heat capacity vector, Dry is a saturated liquid phase heat conductivity coefficient vector, and Ndy is a saturated liquid phase dynamic viscosity vector.
And step S3, dividing the tube into N equal parts according to equal enthalpy difference by adopting N-1 nodes, and calculating the heat exchange coefficient vector of each node according to the saturated liquid phase Plantt number vector.
In one embodiment of the present invention, the heat transfer coefficient vector of the node is
Wherein Pry is the Prandtl number vector of the saturated liquid phase, Dry is the heat conductivity coefficient vector of the saturated liquid phase, and Reynolds number vector ReEQ＝GEQ.Dr./NdlEquivalent mass flow vector GEQ＝G((1-X)+X.*(Mdy./MdgΛ (0.5)), dry-weight vector X ═ hi./Qr-Hl./Qr；hiIs the specific enthalpy value of the ith node, HlIs a specific enthalpy vector, Qr is a flow vector, G is a mass flow vector, MdyIs the refrigerant liquid phase density vector MdgIs the refrigerant gas phase density vector, DrIs hydraulic diameter, NdlIs the viscosity vector.
And step S4, based on the energy relationship, calculating a first relationship curve of the coordinate vector and the saturation pressure vector of a certain node according to the heat exchange coefficient vector of the node and the heat flow loaded on the node.
In an embodiment of the present invention, based on the energy relationship, calculating a first relationship curve between a coordinate vector of a node and a saturation pressure vector according to a heat transfer coefficient vector of the node and a heat flow loaded on the node includes the following steps:
41) calculating the length vector from the previous node to a certain node according to the heat exchange coefficient vector of the node and the heat flow loaded on the node;
42) calculating the coordinate parameter of the node according to the coordinate parameter of the previous node and the length vector;
43) and calculating a first relation curve of the coordinate vector of the node and the saturation pressure vector based on the energy relation.
And step S5, calculating a second relation curve of the coordinate vector of the node and the saturated pressure vector based on the pressure relation.
In an embodiment of the present invention, when calculating the second relationship curve between the coordinate vector of the node and the saturation pressure vector based on the pressure relationship, the second relationship curve is obtained according to
Calculating a pressure vector at the node, wherein DrIs hydraulic diameter, /)i-1Is the coordinate parameter, rho, of the i-1 th nodei-1The refrigerant density of the i-1 th node, Frp the average friction coefficient vector, Mdp the average density vector, L "the next node coordinate vector, and Md the current node density vector; g is the mass flow vector, DrIs the hydraulic diameter.
And step S6, determining the coordinate parameter of the node according to the intersection point of the first relation curve and the second relation curve.
And step S7, calculating the heat exchange coefficient at each node according to the coordinate parameters of the node.
Specifically, the heat exchange coefficient calculation mode has various empirical formulas in the prior art, and the method is suitable for any empirical formula.
As shown in fig. 2, in an embodiment, the parameter calculation system for heat exchange of refrigerant flowing in a tube according to the present invention includes a setting module 21, a first calculation module 22, a second calculation module 23, a third calculation module 24, a fourth calculation module 25, a determination module 26, and a fifth calculation module 27.
The setting module 21 is used for setting the pressure in the state equation of the thermodynamic system as an independent variable.
In describing the strength parameters in a constant-mass thermodynamic system, only a certain number of parameters are independent state parameters, the number of independent state parameters is called the degree of freedom of the system, in a thermodynamic equilibrium state, the number N of independent state parameters describing the system is α - β +2, wherein α is the component and β is the phase state number, therefore, for a single-component constant-mass thermodynamic system, there is only one independent variable in two phase states (α is 1 and β is 2).
Any one of three parameters p, v and T of a thermodynamic system with a state equation f (p, v and T) equal to 0 serves as an independent variable of the system, and all other physical parameters can be expressed as functions of the independent variable. Where p is pressure, v is volume, and T is temperature. In the present invention, p is chosen as an independent variable.
The first calculating module 22 is connected to the setting module 21 and configured to calculate a saturated liquid phase prandtl number vector.
In one embodiment of the present invention, the saturated liquidus prandtl number vector is calculated according to Pry Ndy/Dry; wherein, Pry is a saturated liquid phase Prandtl number vector, Cpy is a saturated liquid phase constant pressure specific heat capacity vector, Dry is a saturated liquid phase heat conductivity coefficient vector, and Ndy is a saturated liquid phase dynamic viscosity vector.
The second calculation module 23 is connected to the first calculation module 22, and is configured to divide the tube into N equal parts according to the equal enthalpy difference by using N-1 nodes, and calculate the heat exchange coefficient vector of each node according to the saturated liquid phase prandtl number vector.
In one embodiment of the present invention, the heat transfer coefficient vector of the node is
Wherein Pry is the Prandtl number vector of the saturated liquid phase, Dry is the heat conductivity coefficient vector of the saturated liquid phase, and Reynolds number vector ReEQ＝GEQ.Dr./NdlEquivalent mass flow vector GEQ＝G((1-X)+X.*(Mdy./MdgΛ (0.5)), dry-weight vector X ═ hi./Qr-Hl./Qr；hiIs the specific enthalpy value of the ith node, HlIs a specific enthalpy vector, Qr is a flow vector, G is a mass flow vector, MdyIs the refrigerant liquid phase density vector MdgTo makeGas phase density vector of refrigerant, DrIs hydraulic diameter, NdlIs the viscosity vector.
The third calculating module 24 is connected to the second calculating module 23, and is configured to calculate, based on an energy relationship, a first relationship curve between a coordinate vector of a node and a saturation pressure vector according to a heat exchange coefficient vector of the node and a heat flow loaded on the node.
In an embodiment of the present invention, based on the energy relationship, calculating a first relationship curve between a coordinate vector of a node and a saturation pressure vector according to a heat transfer coefficient vector of the node and a heat flow loaded on the node includes the following steps:
41) calculating the length vector from the previous node to a certain node according to the heat exchange coefficient vector of the node and the heat flow loaded on the node;
42) calculating the coordinate parameter of the node according to the coordinate parameter of the previous node and the length vector;
43) and calculating a first relation curve of the coordinate vector of the node and the saturation pressure vector based on the energy relation.
The fourth calculating module 25 is configured to calculate a second relation curve of the coordinate vector of the node and the saturation pressure vector based on the pressure relation.
In an embodiment of the present invention, when calculating the second relationship curve between the coordinate vector of the node and the saturation pressure vector based on the pressure relationship, the second relationship curve is obtained according to
Calculating a pressure vector at the node, wherein DrIs hydraulic diameter, /)i-1Is the coordinate parameter, rho, of the i-1 th nodei-1The refrigerant density of the i-1 th node, Frp the average friction coefficient vector, Mdp the average density vector, L "the next node coordinate vector, and Md the current node density vector; g is the mass flow vector, DrIs the hydraulic diameter.
The determining module 26 is connected to the third calculating module 24 and the fourth calculating module 25, and is configured to determine the coordinate parameter of the node according to the intersection of the first relation curve and the second relation curve.
The fifth calculating module 27 is connected to the determining module 26, and is configured to calculate the heat exchange coefficient at each node according to the coordinate parameter of the node.
Specifically, the heat exchange coefficient calculation mode has various empirical formulas in the prior art, and the method is suitable for any empirical formula.
It should be noted that the division of the modules of the above system is only a logical division, and the actual implementation may be wholly or partially integrated into one physical entity, or may be physically separated. And these modules can be realized in the form of software called by processing element; or may be implemented entirely in hardware; and part of the modules can be realized in the form of calling software by the processing element, and part of the modules can be realized in the form of hardware. For example, the x module may be a processing element that is set up separately, or may be implemented by being integrated in a chip of the apparatus, or may be stored in a memory of the apparatus in the form of program code, and the function of the x module may be called and executed by a processing element of the apparatus. Other modules are implemented similarly. In addition, all or part of the modules can be integrated together or can be independently realized. The processing element described herein may be an integrated circuit having signal processing capabilities. In implementation, each step of the above method or each module above may be implemented by an integrated logic circuit of hardware in a processor element or an instruction in the form of software.
For example, the above modules may be one or more integrated circuits configured to implement the above methods, such as: one or more Application Specific Integrated Circuits (ASICs), or one or more microprocessors (DSPs), or one or more Field Programmable Gate Arrays (FPGAs), etc. For another example, when one of the above modules is implemented in the form of a processing element scheduler code, the processing element may be a general-purpose processor, such as a Central Processing Unit (CPU) or other processor capable of calling program code. For another example, these modules may be integrated together and implemented in the form of a system-on-a-chip (SOC).
The storage medium of the present invention stores thereon a computer program which, when executed by a processor, implements the above-described method for calculating a parameter for flowing heat exchange of a refrigerant in a tube. Preferably, the storage medium includes: various media that can store program codes, such as ROM, RAM, magnetic or optical disks.
As shown in fig. 3, in one embodiment, the terminal of the present invention includes a processor 31 and a memory 32.
The memory 32 is used for storing computer programs.
Preferably, the memory 32 comprises: various media that can store program codes, such as ROM, RAM, magnetic or optical disks.
The processor 31 is connected to the memory 32 and is configured to execute a computer program stored in the memory 32, so that the terminal performs the above-mentioned parameter calculation method for refrigerant flow heat exchange in the pipe.
Preferably, the processor 31 may be a general-purpose processor, including a Central Processing Unit (CPU), a Network Processor (NP), and the like; the integrated circuit may also be a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a Field Programmable Gate Array (FPGA) or other programmable logic device, discrete gate or transistor logic device, or discrete hardware components.
As shown in fig. 4, in an embodiment, the battery pack of the present invention includes a temperature-equalizing cold plate 41, wherein a refrigerant 411 in a tube is disposed in the temperature-equalizing cold plate 41, and the refrigerant in the tube sets flowing heat exchange parameters according to the parameter calculation method for flowing heat exchange of the refrigerant in the tube.
Specifically, the flowing heat exchange parameters are set by the parameter calculation method for flowing heat exchange of the refrigerant in the tube, so that the requirement for flowing heat exchange of the battery pack can be fully met.
In summary, the parameter calculation method, system, medium, terminal and battery pack for the flowing heat exchange of the refrigerant in the pipe do not perform constant assumption on the physical properties and the intermediate parameter of any related refrigerant along the flowing direction, and can quickly and accurately calculate the flowing heat exchange parameter of the refrigerant in the pipe; iteration and interpolation in the numerical calculation process are avoided or reduced, and the calculation time and the calculation precision are greatly improved; the normal realization of the flowing heat exchange of the refrigerant in the refrigerating pipe is ensured, so that the liquid cooling plate component has simple, light and compact structure, excellent heat transfer performance, lower flow resistance and low cost; the flowing heat exchange requirement of the battery pack is met. Therefore, the invention effectively overcomes various defects in the prior art and has high industrial utilization value.
The foregoing embodiments are merely illustrative of the principles and utilities of the present invention and are not intended to limit the invention. Any person skilled in the art can modify or change the above-mentioned embodiments without departing from the spirit and scope of the present invention. Accordingly, it is intended that all equivalent modifications or changes which can be made by those skilled in the art without departing from the spirit and technical spirit of the present invention be covered by the claims of the present invention.
Claims (10)
1. A parameter calculation method for flow heat exchange of a refrigerant in a pipe is characterized by comprising the following steps: the method comprises the following steps:
setting the pressure in a state equation of the thermodynamic system as an independent variable;
calculating a saturated liquid phase prandtl number vector;
dividing the tube into N equal parts according to equal enthalpy difference by adopting N-1 nodes, and calculating the heat exchange coefficient vector of each node according to the saturated liquid phase Prandtl number vector;
based on the energy relation, calculating a first relation curve of a coordinate vector and a saturation pressure vector of a certain node according to a heat exchange coefficient vector of the node and heat flow loaded on the node;
calculating a second relation curve of the coordinate vector and the saturation pressure vector of the node based on the pressure relation;
determining the coordinate parameter of the node according to the intersection point of the first relation curve and the second relation curve;
and calculating the heat exchange coefficient of each node according to the coordinate parameters of the node.
2. The method of calculating the parameter of the flowing heat exchange of the refrigerant in the pipe according to claim 1, wherein: the state equation of the thermodynamic system is f (p, v, T) ═ 0, wherein p is pressure, v is volume, and T is temperature.
3. The method of calculating the parameter of the flowing heat exchange of the refrigerant in the pipe according to claim 1, wherein: calculating a saturated liquid phase prandtl number vector according to Pry-Cpy Ndy/Dry; wherein, Pry is a saturated liquid phase Prandtl number vector, Cpy is a saturated liquid phase constant pressure specific heat capacity vector, Dry is a saturated liquid phase heat conductivity coefficient vector, and Ndy is a saturated liquid phase dynamic viscosity vector.
4. The method of calculating the parameter of the flowing heat exchange of the refrigerant in the pipe according to claim 1, wherein: the heat transfer coefficient vector of the node is
Wherein Pry is the Prandtl number vector of the saturated liquid phase, Dry is the heat conductivity coefficient vector of the saturated liquid phase, and Reynolds number vector ReEQ＝GEQDr/NdlEquivalent mass flow vector GEQ＝G((1-X)+X*(Mdy/Mdg) Λ (0.5)), and the dryness vector X ═ hi/Qr-Hl/Qr；hiIs the specific enthalpy value of the ith node, HlIs a specific enthalpy vector, Qr is a flow vector, G is a mass flow vector, MdyIs the refrigerant liquid phase density vector MdgIs the refrigerant gas phase density vector, DrIs hydraulic diameter, NdlIs the viscosity vector.
5. The method of calculating the parameter of the flowing heat exchange of the refrigerant in the pipe according to claim 1, wherein: based on the energy relation, calculating a first relation curve of a coordinate vector and a saturation pressure vector of a certain node according to a heat exchange coefficient vector of the node and heat flow loaded on the node comprises the following steps:
calculating the length vector from the previous node to a certain node according to the heat exchange coefficient vector of the node and the heat flow loaded on the node;
calculating the coordinate parameter of the node according to the coordinate parameter of the previous node and the length vector;
and calculating a first relation curve of the coordinate vector of the node and the saturation pressure vector based on the energy relation.
6. The method of calculating the parameter of the flowing heat exchange of the refrigerant in the pipe according to claim 1, wherein: based on the pressure relation, calculating a second relation curve between the coordinate vector of the node and the saturation pressure vector
Calculating a pressure vector at the node, wherein DrIs hydraulic diameter, /)i-1Is the coordinate parameter, rho, of the i-1 th nodei-1The refrigerant density of the i-1 th node, Frp the average friction coefficient vector, Mdp the average density vector, L "the next node coordinate vector, and Md the current node density vector; g is the mass flow vector, DrIs the hydraulic diameter.
7. A parameter calculation system for in-tube refrigerant flow heat exchange is characterized in that: the device comprises a setting module, a first calculating module, a second calculating module, a third calculating module, a fourth calculating module, a determining module and a fifth calculating module;
the setting module is used for setting the pressure in a state equation of the thermodynamic system as an independent variable;
the first calculation module is used for calculating a saturated liquid phase Prandtl number vector;
the second calculation module is used for dividing the tube into N equal parts according to the equal enthalpy difference by adopting N-1 nodes, and calculating the heat exchange coefficient vector of each node according to the saturated liquid phase Plantt number vector;
the third calculation module is used for calculating a first relation curve of a coordinate vector and a saturation pressure vector of a certain node according to the heat exchange coefficient vector of the node and the heat flow loaded on the node based on the energy relation;
the fourth calculation module is used for calculating a second relation curve of the coordinate vector of the node and the saturated pressure vector based on the pressure relation;
the determining module is used for determining the coordinate parameter of the node according to the intersection point of the first relation curve and the second relation curve;
and the fifth calculation module is used for calculating the heat exchange coefficient at each node according to the coordinate parameters of the node.
8. A storage medium having stored thereon a computer program, characterized in that the program when executed by a processor implements the parameter calculating method for refrigerant flow heat exchange in tubes of any one of claims 1 to 6.
9. A terminal comprising a processor and a memory;
the memory is used for storing a computer program;
the processor is used for executing the computer program stored in the memory to enable the terminal to execute the parameter calculation method for the heat exchange of the refrigerant flowing in the pipe as claimed in any one of claims 1 to 6.
10. A battery pack is characterized by comprising a temperature equalizing cold plate, wherein a refrigerant in a pipe is arranged in the temperature equalizing cold plate, and the refrigerant in the pipe sets flowing heat exchange parameters according to the parameter calculation method for flowing heat exchange of the refrigerant in the pipe of any one of claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711216003.9A CN108134156B (en) | 2017-11-28 | 2017-11-28 | Parameter calculation method, system, medium, terminal and battery pack for flowing heat exchange of refrigerant in pipe |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711216003.9A CN108134156B (en) | 2017-11-28 | 2017-11-28 | Parameter calculation method, system, medium, terminal and battery pack for flowing heat exchange of refrigerant in pipe |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108134156A CN108134156A (en) | 2018-06-08 |
| CN108134156B true CN108134156B (en) | 2020-02-14 |
Family
ID=62388747
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711216003.9A Active CN108134156B (en) | 2017-11-28 | 2017-11-28 | Parameter calculation method, system, medium, terminal and battery pack for flowing heat exchange of refrigerant in pipe |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108134156B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109783904B (en) * | 2018-12-28 | 2020-10-27 | 西安交通大学 | Wide parameter range carbon dioxide physical property solving method |
| CN113468463B (en) * | 2021-09-02 | 2022-02-22 | 中兴通讯股份有限公司 | Parameter calculation method for two-phase cold plate |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003068944A (en) * | 2001-08-28 | 2003-03-07 | Nec Corp | Method and device for computing optimum fin-interval of forced convection cooling type heat sink and heat sink |
| US7908872B2 (en) * | 2006-08-03 | 2011-03-22 | Machflow Energy Inc. | Rare-gas-based bernoulli heat pump and method |
| CN103542621B (en) * | 2013-09-26 | 2016-01-20 | 西安交通大学 | A kind of method for designing of general combination pipe diameter air conditioner heat exchange equipment fluid passage |
| DE102015115761A1 (en) * | 2015-09-18 | 2017-03-23 | Endress + Hauser Flowtec Ag | Method for on-site calibration of a thermal flow measuring device, method for carrying out a temperature-compensated flow measurement and thermal flow meter |
-
2017
- 2017-11-28 CN CN201711216003.9A patent/CN108134156B/en active Active
Non-Patent Citations (2)
| Title |
|---|
| Heat transfer and hydraulic resistance of supercritical-pressure coolants. Part II:;V.A. Kurganov等;《International Journal of Heat and Mass Transfer》;20121208;第58卷(第1-2期);第152-167页 * |
| 多元混合工质相平衡及传热模拟与系统实验研究;邹晓超;《中国优秀硕士学位论文全文数据库工程科技II辑》;20100515;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108134156A (en) | 2018-06-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Shen et al. | System simulation on refrigerant-based battery thermal management technology for electric vehicles | |
| Kuang et al. | Research on control strategy for a battery thermal management system for electric vehicles based on secondary loop cooling | |
| Oliet et al. | Parametric studies on automotive radiators | |
| Junqi et al. | Experimental study of R744 heat pump system for electric vehicle application | |
| Sanaye et al. | Modeling and multi-objective optimization of parallel flow condenser using evolutionary algorithm | |
| Li et al. | Numerical simulation and analysis of transport air conditioning system integrated with passenger compartment | |
| CN108134156B (en) | Parameter calculation method, system, medium, terminal and battery pack for flowing heat exchange of refrigerant in pipe | |
| Krüger et al. | Energy consumption of battery cooling in hybrid electric vehicles | |
| Wang et al. | Parametric investigation on the performance of a direct evaporation cooling battery thermal management system | |
| CN105387650A (en) | Two-pass evaporator | |
| CN210778890U (en) | Battery cooling integrated system and electric vehicle | |
| Moria et al. | Exergoeconomic analysis of a Peltier effect air cooler using experimental data | |
| Hamut et al. | Experimental and theoretical efficiency investigation of hybrid electric vehicle battery thermal management systems | |
| Bagheri et al. | Theoretical and experimental investigation into anti-idling A/C system for trucks | |
| Shen et al. | Simulation and analysis of dual-evaporator refrigeration system for electric vehicles | |
| Wang et al. | Thermal performance analysis and burning questions of refrigerant direct cooling for electric vehicle battery | |
| Wang et al. | Performance investigation of integrated thermal management system based on a pumped two-phase cooling system for electric vehicles | |
| Awad et al. | Thermodynamic optimization | |
| CN115081161A (en) | Temperature analysis method, device, equipment and medium for power battery system | |
| Sevilgen et al. | The investigation of the innovative hybrid heat pump system designed and prototyped for heating process of electric vehicles | |
| Song et al. | Air conditioning system sizing for pure electric vehicle | |
| Ismael et al. | Radiator heat dissipation performance | |
| CN114778142A (en) | Electric automobile heat management waterway system testing device | |
| Bayram et al. | Thermal analysis of different heat pump systems for heating process of electric vehicles | |
| Fan et al. | Modeling and Simulation Analysis of Electric Vehicle Thermal Management System Based on Distributed Parameter Method |
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 | ||
| CB02 | Change of applicant information |
Address after: 334000 Xingyuan West Avenue, Shangrao economic and Technological Development Zone, Shangrao, Jiangxi, China Applicant after: AIWAYS AUTOMOBILE Co.,Ltd. Address before: 334000 Xingyuan West Avenue, Shangrao economic and Technological Development Zone, Shangrao, Jiangxi, China Applicant before: JIANGXI AICHI YIWEI INDUSTRIAL CO.,LTD. |
|
| CB02 | Change of applicant information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| PP01 | Preservation of patent right |
Effective date of registration: 20230131 Granted publication date: 20200214 |
|
| PP01 | Preservation of patent right | ||
| PD01 | Discharge of preservation of patent |
Date of cancellation: 20240108 Granted publication date: 20200214 |
|
| PD01 | Discharge of preservation of patent | ||
| PP01 | Preservation of patent right |
Effective date of registration: 20240227 Granted publication date: 20200214 |
|
| PP01 | Preservation of patent right |