CN114235891B - Heat exchange optimization method for high-pressure plate-fin heat exchanger of air separation device - Google Patents
Heat exchange optimization method for high-pressure plate-fin heat exchanger of air separation device Download PDFInfo
- Publication number
- CN114235891B CN114235891B CN202111554014.4A CN202111554014A CN114235891B CN 114235891 B CN114235891 B CN 114235891B CN 202111554014 A CN202111554014 A CN 202111554014A CN 114235891 B CN114235891 B CN 114235891B
- Authority
- CN
- China
- Prior art keywords
- heat exchanger
- plate
- fin heat
- heat exchange
- fin
- 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
- 238000000926 separation method Methods 0.000 title claims abstract description 47
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000005457 optimization Methods 0.000 title claims abstract description 34
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims abstract description 96
- 238000009826 distribution Methods 0.000 claims abstract description 66
- 229910052757 nitrogen Inorganic materials 0.000 claims abstract description 48
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims abstract description 37
- 239000001301 oxygen Substances 0.000 claims abstract description 37
- 229910052760 oxygen Inorganic materials 0.000 claims abstract description 37
- 238000004458 analytical method Methods 0.000 claims abstract description 21
- 238000012795 verification Methods 0.000 claims abstract description 21
- 238000004088 simulation Methods 0.000 claims description 13
- 238000004422 calculation algorithm Methods 0.000 claims description 11
- 241000254032 Acrididae Species 0.000 claims description 8
- 230000004044 response Effects 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 7
- 102100021807 ER degradation-enhancing alpha-mannosidase-like protein 1 Human genes 0.000 claims description 5
- 101000895701 Homo sapiens ER degradation-enhancing alpha-mannosidase-like protein 1 Proteins 0.000 claims description 5
- 230000009466 transformation Effects 0.000 claims description 5
- 230000002068 genetic effect Effects 0.000 claims description 3
- 239000002245 particle Substances 0.000 claims description 3
- 238000013461 design Methods 0.000 abstract description 6
- 230000009286 beneficial effect Effects 0.000 abstract description 4
- 230000001276 controlling effect Effects 0.000 description 36
- 230000000694 effects Effects 0.000 description 19
- 239000012530 fluid Substances 0.000 description 10
- 238000010992 reflux Methods 0.000 description 10
- 230000008569 process Effects 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 238000012546 transfer Methods 0.000 description 5
- MYMOFIZGZYHOMD-UHFFFAOYSA-N Dioxygen Chemical compound O=O MYMOFIZGZYHOMD-UHFFFAOYSA-N 0.000 description 4
- 239000007788 liquid Substances 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 238000012986 modification Methods 0.000 description 4
- 230000009467 reduction Effects 0.000 description 4
- 230000006872 improvement Effects 0.000 description 3
- 230000001105 regulatory effect Effects 0.000 description 3
- XKRFYHLGVUSROY-UHFFFAOYSA-N Argon Chemical compound [Ar] XKRFYHLGVUSROY-UHFFFAOYSA-N 0.000 description 2
- 238000001816 cooling Methods 0.000 description 2
- 239000003337 fertilizer Substances 0.000 description 2
- 229910052786 argon Inorganic materials 0.000 description 1
- 238000005219 brazing Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000013499 data model Methods 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 238000007710 freezing Methods 0.000 description 1
- 230000008014 freezing Effects 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 239000011261 inert gas Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000010926 purge Methods 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 238000010223 real-time analysis Methods 0.000 description 1
- 238000011897 real-time detection Methods 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N25/00—Investigating or analyzing materials by the use of thermal means
- G01N25/20—Investigating or analyzing materials by the use of thermal means by investigating the development of heat, i.e. calorimetry, e.g. by measuring specific heat, by measuring thermal conductivity
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Separation By Low-Temperature Treatments (AREA)
Abstract
The invention discloses a heat exchange optimization method for a high-pressure plate-fin heat exchanger of an air separation device, which comprises the following steps: s1, establishing a dynamic model of an internal flow channel structure of a plate-fin heat exchanger; s2, simulating and analyzing flow state parameters A and B when medium-pressure air, high-pressure air, pure nitrogen and high-pressure oxygen flow through the inside of the plate-fin heat exchanger to obtain the flow velocity distribution conditions of heat sources and cold sources in the plate-fin heat exchanger and the temperature distribution conditions of heat exchange surfaces of the heat sources and the cold sources; s3, substituting the flow state parameter A and the flow state parameter B into a dynamic model of an internal flow channel structure of the plate-fin heat exchanger, and calculating optimal solutions of V1, V2, V3 and V4 affecting heat exchange loss of the plate-fin heat exchanger; s4, substituting the optimal solutions of V1, V2, V3 and V4 into a dynamic model of the internal flow channel structure of the plate-fin heat exchanger for verification analysis; the invention has reasonable design, is beneficial to reducing the cold energy loss of the plate-fin heat exchanger, thereby improving the productivity of the air separation device and being suitable for mass popularization.
Description
Technical Field
The invention relates to the technical field of heat exchangers, in particular to a heat exchange optimization method for a high-pressure plate-fin heat exchanger of an air separation device.
Background
The air separation device is a device which takes air as a raw material, turns the air into liquid state by a compression circulation deep freezing method, and gradually separates and produces inert gases such as oxygen, nitrogen, argon and the like from the liquid air by rectification; the air separation plant produced in China has various forms and types, and is provided with a device for producing gaseous oxygen and nitrogen and a device for producing liquid oxygen and nitrogen; however, as for the basic flow, there are mainly four kinds of flow, namely, high pressure, medium pressure, high and low pressure and full low pressure flow; the air separation device is a large complex system and mainly comprises a power system, a purifying system, a refrigerating system, a heat exchange system, a rectifying system, a product conveying system, a liquid storage system, a control system and the like.
The plate-fin heat exchanger is used as one of high-efficiency heat exchange equipment, has very light weight and low cost due to high heat transfer efficiency and compact structure, and has high mechanical strength due to a full-brazing structure; in addition, the heat transfer area density at two sides of the plate-fin heat exchanger can be different by more than one order of magnitude, the heat transfer area density can adapt to the difference of heat transfer of media at two sides, the utilization rate of the heat transfer surface is improved, multi-strand fluid heat exchange can be organized, and the exchange between the whole heat exchanger and the surrounding environment is reduced to the greatest extent. With numerous advantages, plate-fin heat exchangers are increasingly being used.
However, the high pressure plate-fin heat exchanger in the prior art has the following problems in use: 1. the high-pressure plate-fin heat exchanger has large temperature difference, so that the cold loss is large, and the liquid yield is reduced year by year; 2. when the system cooling capacity loss exceeds the cooling capacity, the air separation device is forced to stop for heating and purging, and the long-period stable operation of the large fertilizer is affected.
Disclosure of Invention
Aiming at the technical problems, the invention provides a heat exchange optimization method for a high-pressure plate-fin heat exchanger of an air separation device.
The technical scheme of the invention is as follows: a heat exchange optimization method of a high-pressure plate-fin heat exchanger of an air separation device comprises the following steps:
s1, establishing a model
Modeling the internal flow passage structure of the plate-fin heat exchanger by adopting FLUENT software to obtain a dynamic model of the internal flow passage structure of the plate-fin heat exchanger;
s2, parameter acquisition
S2-1, according to the dynamic model of the internal flow channel structure of the plate-fin heat exchanger obtained in the step S1, adopting FLUENT software to simulate and analyze flow state parameters A of medium-pressure air and high-pressure air when the medium-pressure air and the high-pressure air flow pass through the inside of the plate-fin heat exchanger, and obtaining the flow velocity distribution condition of a heat source in the plate-fin heat exchanger and the temperature distribution condition of a heat exchange surface of the heat source;
s2-2, simulating and analyzing flow state parameters B of pure nitrogen and high-pressure oxygen when the pure nitrogen and the high-pressure oxygen flow through the inside of the plate-fin heat exchanger by adopting FLUENT software to obtain the flow velocity distribution condition of a cold source inside the plate-fin heat exchanger and the temperature distribution condition of a heat exchange surface of the cold source;
s3, heat exchange loss simulation
Substituting the flow state parameter A obtained in the step S2-1 and the flow state parameter B obtained in the step S2-2 into a dynamic model of the internal flow channel structure of the plate-fin heat exchanger obtained in the step S1, taking the opening V1 of a valve for controlling medium-pressure air flow, the opening V2 of a valve for controlling high-pressure air flow, the opening V3 of a valve for controlling pure nitrogen flow and the opening V4 of a valve for controlling high-pressure oxygen as parameter variables, and simulating the interrelationship of the flow state parameter A, the flow state parameter B, the V1, the V2, the V3 and the V4 by using the dynamic model of the internal flow channel structure of the plate-fin heat exchanger to calculate the corresponding optimal solution of the V1, the V2, the V3 and the V4 affecting the heat exchange loss of the plate-fin heat exchanger;
s4, analysis and verification
And substituting the optimal solutions corresponding to V1, V2, V3 and V4 obtained in the step S3 into the dynamic model of the internal flow channel structure of the plate-fin heat exchanger obtained in the step S1, and performing verification analysis on the distribution conditions of the flow velocity of the heat source and the cold source in the plate-fin heat exchanger and the distribution conditions of the temperature of the heat exchange surfaces of the heat source and the cold source.
Further, in the steps S2-1 and S2-2, the flow state parameters A and B comprise temperature variable parameters, pressure variable parameters and flow variable parameters, and the heat exchange effect of the plate-fin heat exchanger is designed in an intervention way by setting multiple parameters, so that the optimization effect of the invention is better outstanding.
Further, the specific operation of step S4 is:
s4-1, constructing a heat exchange function based on a sine and cosine function according to a dynamic model of an internal flow channel structure of the plate-fin heat exchanger; the expression of the heat exchange function is as follows:
XL＝x-bCos[(x/2)^2]
YL＝y-bCos[(y/2)^2]
ZL＝z-bCos[(z/2)^2]
p1＝ContourPlot3D[Cos[XL]Sin[YL]+Con[YL]Sin[ZL]+Cos[ZL]Sin[XL]＝0，{x，-Pi，Pi}，{y，-Pi，Pi}，{z，-Pi，Pi}]
wherein, p1 is a heat exchange function, contourPlot3D is a three-dimensional space transformation function, x is an x-axis coordinate value, y is a y-axis coordinate value, z is a z-axis coordinate value, and b is a variable parameter;
s4-2: performing three-dimensional modeling on the heat exchange function obtained in the step S4-1, and constructing corresponding unit bodies of V1, V2, V3 and V4;
s4-3, adjusting the unit bodies corresponding to V1, V2, V3 and V4 obtained in the step S4-2, obtaining the optimal unit bodies corresponding to V1, V2, V3 and V4 meeting preset conditions, and realizing verification analysis of the flow velocity distribution conditions of the heat source and the cold source inside the plate-fin heat exchanger and the temperature distribution conditions of the heat exchange surfaces of the heat source and the cold source; by obtaining the optimal unit bodies corresponding to V1, V2, V3 and V4, the resistance of the fluid in the plate-fin heat exchanger can be reduced, and the productivity of the air separation device is improved.
Further, leading the optimal unit bodies corresponding to V1, V2, V3 and V4 into a Grasshopper parametric modeling platform, and obtaining an aggregate unit consisting of the optimal unit bodies corresponding to V1, V2, V3 and V4 through a rapid iteration array; the V1, V2, V3 and V4 of the invention are regulated more rapidly by using a Grasshopper parameterized modeling platform.
Further, the specific operation of step S1 is:
s1-1, modeling a plate-fin heat exchanger by using Gambit software, selecting a grid type, setting grid boundary conditions, and deriving a grid file;
s1-2, importing the grid file obtained in the step S1-1 into an EDEM, and then filling the internal flow passage structural parameters of the plate-fin heat exchanger to obtain a plate-fin heat exchanger internal flow passage structural parameter model;
s1-3, introducing the parameter model of the internal flow passage structure of the plate-fin heat exchanger obtained in the step S1-2 into FLUENT, and obtaining a dynamic model of the internal flow passage structure of the plate-fin heat exchanger; by the modeling method, a reliable research environment can be provided for the optimal design of the heat exchange effect of the plate-fin heat exchanger.
Further, after the step S1 is completed, the dynamic model of the internal flow channel structure of the plate-fin heat exchanger is divided by adopting a structured grid; the structural grid is adopted to divide the dynamic model of the internal flow channel structure of the plate-fin heat exchanger, so that the optimization calculation time of each parameter is shortened, and the modification of complex variables in the model is facilitated.
Further, after step S3 is completed, according to the adjustable parameter ranges of the flow state parameter a and the flow state parameter B, performing up-down floating adjustment on the flow state parameter a and the flow state parameter B within the adjustable parameter ranges of the flow state parameter a and the flow state parameter B, and then verifying the response speeds of V1, V2, V3 and V4 by using the variable amounts of the flow state parameter a and the flow state parameter B; by verifying the response speeds of V1, V2, V3 and V4, the problem of capacity reduction of the air separation device caused by heat exchange loss can be avoided.
Further, after S4 is completed, performing reflux nitrogen temperature simulation on the dynamic model of the internal flow channel structure of the plate-fin heat exchanger by utilizing LUENT software, calculating the theoretical value of the reflux nitrogen temperature of the plate-fin heat exchanger according to the optimal solution corresponding to V1, V2, V3 and V4 in the step S4, and calculating the heat exchange coefficient of the plate-fin heat exchanger according to the theoretical value; the heat exchange coefficient of the plate-fin heat exchanger is calculated, so that the improvement of the plate-fin heat exchanger is facilitated, and the heat exchange effect is improved.
Further, in step S3, the algorithm used for calculation is one of a genetic algorithm, a particle swarm algorithm, and an ant colony algorithm.
Further, in step S3, the parameter variables further include a medium pressure air temperature T1, a high pressure air temperature T2, a pure nitrogen temperature T3, and a high pressure oxygen temperature T4; the influence of the flow and the temperature of the fluid entering the plate-fin heat exchanger on the heat exchange loss of the air separation device is simulated, so that the plate-fin heat exchanger can be subjected to multi-dimensional adjustment in the use process.
Compared with the prior art, the invention has the beneficial effects that:
firstly, the FLUENT software is utilized to simulate the internal flow channel structure of the plate-fin heat exchanger, so that the optimization method is more convenient to implement, manpower and material resources are greatly saved, and the requirements of energy conservation and emission reduction are met;
secondly, the invention can obtain the result through the operation of a computer, has simple and convenient operation and convenient model modification, and can adapt to complex models; meanwhile, the optimization method is convenient for real-time detection and analysis of the capacity of the air separation device, so that the parameters of the plate-fin heat exchanger can be timely adjusted, the loss of the plate-fin heat exchanger is reduced, and the capacity of the air separation device is improved.
Thirdly, the dynamic model of the internal flow channel structure of the plate-fin heat exchanger is convenient for transmitting real-time data of the air separation device production to a data model for optimization calculation, so that real-time optimization of the air separation device production process is realized; the optimization method disclosed by the invention has small environment dependence and can meet the complex requirements of large-scale air separation device production.
Drawings
FIG. 1 is a flow chart of a heat exchange optimization method of a plate fin heat exchanger of embodiment 1 of the present invention;
Detailed Description
Example 1
A heat exchange optimization method of a high-pressure plate-fin heat exchanger of an air separation device comprises the following steps:
s1, establishing a model
Modeling the internal flow passage structure of the plate-fin heat exchanger by adopting FLUENT software to obtain a dynamic model of the internal flow passage structure of the plate-fin heat exchanger;
s2, parameter acquisition
S2-1, according to the dynamic model of the internal flow channel structure of the plate-fin heat exchanger obtained in the step S1, adopting FLUENT software to simulate and analyze flow state parameters A of medium-pressure air and high-pressure air when the medium-pressure air and the high-pressure air flow pass through the inside of the plate-fin heat exchanger, and obtaining the flow velocity distribution condition of a heat source in the plate-fin heat exchanger and the temperature distribution condition of a heat exchange surface of the heat source;
s2-2, simulating and analyzing flow state parameters B of pure nitrogen and high-pressure oxygen when the pure nitrogen and the high-pressure oxygen flow through the inside of the plate-fin heat exchanger by adopting FLUENT software to obtain the flow velocity distribution condition of a cold source inside the plate-fin heat exchanger and the temperature distribution condition of a heat exchange surface of the cold source;
s3, heat exchange loss simulation
Substituting the flow state parameter A obtained in the step S2-1 and the flow state parameter B obtained in the step S2-2 into a dynamic model of the internal flow channel structure of the plate-fin heat exchanger obtained in the step S1, taking the opening V1 of a valve for controlling medium-pressure air flow, the opening V2 of a valve for controlling high-pressure air flow, the opening V3 of a valve for controlling pure nitrogen flow and the opening V4 of a valve for controlling high-pressure oxygen as parameter variables, and simulating the interrelationship of the flow state parameter A, the flow state parameter B, the V1, the V2, the V3 and the V4 by using the dynamic model of the internal flow channel structure of the plate-fin heat exchanger to calculate the corresponding optimal solution of the V1, the V2, the V3 and the V4 affecting the heat exchange loss of the plate-fin heat exchanger;
s4, analysis and verification
And substituting the optimal solutions corresponding to V1, V2, V3 and V4 obtained in the step S3 into the dynamic model of the internal flow channel structure of the plate-fin heat exchanger obtained in the step S1, and performing verification analysis on the distribution conditions of the flow velocity of the heat source and the cold source in the plate-fin heat exchanger and the distribution conditions of the temperature of the heat exchange surfaces of the heat source and the cold source.
Example 2
A heat exchange optimization method of a high-pressure plate-fin heat exchanger of an air separation device comprises the following steps:
s1, establishing a model
S1-1, modeling a plate-fin heat exchanger by using Gambit software, selecting a grid type, setting grid boundary conditions, and deriving a grid file;
s1-2, importing the grid file obtained in the step S1-1 into an EDEM, and then filling the internal flow passage structural parameters of the plate-fin heat exchanger to obtain a plate-fin heat exchanger internal flow passage structural parameter model;
s1-3, introducing the parameter model of the internal flow passage structure of the plate-fin heat exchanger obtained in the step S1-2 into FLUENT, and obtaining a dynamic model of the internal flow passage structure of the plate-fin heat exchanger; by the modeling method, a reliable research environment can be provided for the optimal design of the heat exchange effect of the plate-fin heat exchanger;
s2, parameter acquisition
S2-1, according to the dynamic model of the internal flow channel structure of the plate-fin heat exchanger obtained in the step S1, adopting FLUENT software to simulate and analyze flow state parameters A of medium-pressure air and high-pressure air when the medium-pressure air and the high-pressure air flow pass through the inside of the plate-fin heat exchanger, and obtaining the flow velocity distribution condition of a heat source in the plate-fin heat exchanger and the temperature distribution condition of a heat exchange surface of the heat source;
s2-2, simulating and analyzing flow state parameters B of pure nitrogen and high-pressure oxygen when the pure nitrogen and the high-pressure oxygen flow through the inside of the plate-fin heat exchanger by adopting FLUENT software to obtain the flow velocity distribution condition of a cold source inside the plate-fin heat exchanger and the temperature distribution condition of a heat exchange surface of the cold source;
s3, heat exchange loss simulation
Substituting the flow state parameter A obtained in the step S2-1 and the flow state parameter B obtained in the step S2-2 into a dynamic model of the internal flow channel structure of the plate-fin heat exchanger obtained in the step S1, taking the opening V1 of a valve for controlling medium-pressure air flow, the opening V2 of a valve for controlling high-pressure air flow, the opening V3 of a valve for controlling pure nitrogen flow and the opening V4 of a valve for controlling high-pressure oxygen as parameter variables, and simulating the interrelationship of the flow state parameter A, the flow state parameter B, the V1, the V2, the V3 and the V4 by using the dynamic model of the internal flow channel structure of the plate-fin heat exchanger to calculate the corresponding optimal solution of the V1, the V2, the V3 and the V4 affecting the heat exchange loss of the plate-fin heat exchanger;
s4, analysis and verification
And substituting the optimal solutions corresponding to V1, V2, V3 and V4 obtained in the step S3 into the dynamic model of the internal flow channel structure of the plate-fin heat exchanger obtained in the step S1, and performing verification analysis on the distribution conditions of the flow velocity of the heat source and the cold source in the plate-fin heat exchanger and the distribution conditions of the temperature of the heat exchange surfaces of the heat source and the cold source.
Example 3
A heat exchange optimization method of a high-pressure plate-fin heat exchanger of an air separation device comprises the following steps:
s1, establishing a model
S1-1, modeling a plate-fin heat exchanger by using Gambit software, selecting a grid type, setting grid boundary conditions, and deriving a grid file;
s1-2, importing the grid file obtained in the step S1-1 into an EDEM, and then filling the internal flow passage structural parameters of the plate-fin heat exchanger to obtain a plate-fin heat exchanger internal flow passage structural parameter model;
s1-3, introducing the parameter model of the internal flow passage structure of the plate-fin heat exchanger obtained in the step S1-2 into FLUENT, and obtaining a dynamic model of the internal flow passage structure of the plate-fin heat exchanger; dividing a dynamic model of an internal flow channel structure of the plate-fin heat exchanger by adopting a structured grid; the structural grid is adopted to divide the dynamic model of the internal flow channel structure of the plate-fin heat exchanger, so that the optimization calculation time of each parameter is shortened, and the modification of complex variables in the model is facilitated; by the modeling method, a reliable research environment can be provided for the optimal design of the heat exchange effect of the plate-fin heat exchanger;
s2, parameter acquisition
S2-1, according to the dynamic model of the internal flow channel structure of the plate-fin heat exchanger obtained in the step S1, adopting FLUENT software to simulate and analyze flow state parameters A of medium-pressure air and high-pressure air when the medium-pressure air and the high-pressure air flow pass through the inside of the plate-fin heat exchanger, and obtaining the flow velocity distribution condition of a heat source in the plate-fin heat exchanger and the temperature distribution condition of a heat exchange surface of the heat source;
s2-2, simulating and analyzing flow state parameters B of pure nitrogen and high-pressure oxygen when the pure nitrogen and the high-pressure oxygen flow through the inside of the plate-fin heat exchanger by adopting FLUENT software to obtain the flow velocity distribution condition of a cold source inside the plate-fin heat exchanger and the temperature distribution condition of a heat exchange surface of the cold source;
s3, heat exchange loss simulation
Substituting the flow state parameter A obtained in the step S2-1 and the flow state parameter B obtained in the step S2-2 into a dynamic model of the internal flow channel structure of the plate-fin heat exchanger obtained in the step S1, taking the opening V1 of a valve for controlling medium-pressure air flow, the opening V2 of a valve for controlling high-pressure air flow, the opening V3 of a valve for controlling pure nitrogen flow and the opening V4 of a valve for controlling high-pressure oxygen as parameter variables, and simulating the interrelationship of the flow state parameter A, the flow state parameter B, the V1, the V2, the V3 and the V4 by using the dynamic model of the internal flow channel structure of the plate-fin heat exchanger to calculate the corresponding optimal solution of the V1, the V2, the V3 and the V4 affecting the heat exchange loss of the plate-fin heat exchanger;
s4, analysis and verification
And substituting the optimal solutions corresponding to V1, V2, V3 and V4 obtained in the step S3 into the dynamic model of the internal flow channel structure of the plate-fin heat exchanger obtained in the step S1, and performing verification analysis on the distribution conditions of the flow velocity of the heat source and the cold source in the plate-fin heat exchanger and the distribution conditions of the temperature of the heat exchange surfaces of the heat source and the cold source.
Example 4
A heat exchange optimization method of a high-pressure plate-fin heat exchanger of an air separation device comprises the following steps:
s1, establishing a model
Modeling the internal flow passage structure of the plate-fin heat exchanger by adopting FLUENT software to obtain a dynamic model of the internal flow passage structure of the plate-fin heat exchanger;
s2, parameter acquisition
S2-1, according to the dynamic model of the internal flow channel structure of the plate-fin heat exchanger obtained in the step S1, adopting FLUENT software to simulate and analyze flow state parameters A of medium-pressure air and high-pressure air when the medium-pressure air and the high-pressure air flow pass through the inside of the plate-fin heat exchanger, and obtaining the flow velocity distribution condition of a heat source in the plate-fin heat exchanger and the temperature distribution condition of a heat exchange surface of the heat source;
s2-2, simulating and analyzing flow state parameters B of pure nitrogen and high-pressure oxygen when the pure nitrogen and the high-pressure oxygen flow through the inside of the plate-fin heat exchanger by adopting FLUENT software to obtain the flow velocity distribution condition of a cold source inside the plate-fin heat exchanger and the temperature distribution condition of a heat exchange surface of the cold source;
s3, heat exchange loss simulation
Substituting the flow state parameter A obtained in the step S2-1 and the flow state parameter B obtained in the step S2-2 into a dynamic model of the internal flow channel structure of the plate-fin heat exchanger obtained in the step S1, taking the opening V1 of a valve for controlling medium-pressure air flow, the opening V2 of a valve for controlling high-pressure air flow, the opening V3 of a valve for controlling pure nitrogen flow and the opening V4 of a valve for controlling high-pressure oxygen as parameter variables, and simulating the interrelationship of the flow state parameter A, the flow state parameter B, the V1, the V2, the V3 and the V4 by using the dynamic model of the internal flow channel structure of the plate-fin heat exchanger to calculate the optimal solution of the V1, the V2, the V3 and the V4 affecting the heat exchange loss of the plate-fin heat exchanger;
s4, analysis and verification
S4-1, constructing a heat exchange function based on a sine and cosine function according to a dynamic model of an internal flow channel structure of the plate-fin heat exchanger; the expression of the heat exchange function is as follows:
XL＝x-bCos[(x/2)^2]
YL＝y-bCos[(y/2)^2]
ZL＝z-bCos[(z/2)^2]
p1＝ContourPlot3D[Cos[XL]Sin[YL]+Con[YL]Sin[ZL]+Cos[ZL]Sin[XL]＝0，{x，-Pi，Pi}，{y，-Pi，Pi}，{z，-Pi，Pi}]
wherein, p1 is a heat exchange function, contourPlot3D is a three-dimensional space transformation function, x is an x-axis coordinate value, y is a y-axis coordinate value, z is a z-axis coordinate value, and b is a variable parameter;
s4-2: carrying out three-dimensional modeling on the heat exchange function obtained in the step S4-1, and constructing unit bodies corresponding to V1, V2, V3 and V4 by utilizing V1, V2, V3 and V4 obtained in the step S3;
s4-3, adjusting the unit bodies corresponding to V1, V2, V3 and V4 obtained in the step S4-2, and obtaining the optimal unit bodies corresponding to V1, V2, V3 and V4 meeting preset conditions, thereby realizing verification analysis of the flow velocity distribution conditions of the heat source and the cold source in the plate-fin heat exchanger and the temperature distribution conditions of the heat exchange surfaces of the heat source and the cold source; by obtaining the optimal unit bodies corresponding to V1, V2, V3 and V4, the resistance of the fluid in the plate-fin heat exchanger can be reduced, and the productivity of the air separation device is improved.
Example 5
A heat exchange optimization method of a high-pressure plate-fin heat exchanger of an air separation device comprises the following steps:
s1, establishing a model
Modeling the internal flow passage structure of the plate-fin heat exchanger by adopting FLUENT software to obtain a dynamic model of the internal flow passage structure of the plate-fin heat exchanger;
s2, parameter acquisition
S2-1, according to the dynamic model of the internal flow channel structure of the plate-fin heat exchanger obtained in the step S1, adopting FLUENT software to simulate and analyze flow state parameters A of medium-pressure air and high-pressure air when the medium-pressure air and the high-pressure air flow pass through the inside of the plate-fin heat exchanger, and obtaining the flow velocity distribution condition of a heat source in the plate-fin heat exchanger and the temperature distribution condition of a heat exchange surface of the heat source;
s2-2, simulating and analyzing flow state parameters B of pure nitrogen and high-pressure oxygen when the pure nitrogen and the high-pressure oxygen flow through the inside of the plate-fin heat exchanger by adopting FLUENT software to obtain the flow velocity distribution condition of a cold source inside the plate-fin heat exchanger and the temperature distribution condition of a heat exchange surface of the cold source; the flow state parameters A and the flow state parameters B comprise temperature variable parameters, pressure variable parameters and flow variable parameters, and the heat exchange effect of the plate-fin heat exchanger is designed in an intervention way through setting multiple parameters, so that the optimization effect of the invention is better outstanding;
s3, heat exchange loss simulation
Substituting the flow state parameter A obtained in the step S2-1 and the flow state parameter B obtained in the step S2-2 into a dynamic model of the internal flow channel structure of the plate-fin heat exchanger obtained in the step S1, taking the opening V1 of a valve for controlling medium-pressure air flow, the opening V2 of a valve for controlling high-pressure air flow, the opening V3 of a valve for controlling pure nitrogen flow and the opening V4 of a valve for controlling high-pressure oxygen as parameter variables, simulating the interrelationship of the flow state parameter A, the flow state parameter B, the V1, the V2, the V3 and the V4 by using the dynamic model of the internal flow channel structure of the plate-fin heat exchanger, and calculating the optimal solution corresponding to the V1, the V2, the V3 and the V4 affecting the heat exchange loss of the plate-fin heat exchanger by using an ant algorithm;
s4, analysis and verification
S4-1, constructing a heat exchange function based on a sine and cosine function according to a dynamic model of an internal flow channel structure of the plate-fin heat exchanger; the expression of the heat exchange function is as follows:
XL＝x-bCos[(x/2)^2]
YL＝y-bCos[(y/2)^2]
ZL＝z-bCos[(z/2)^2]
p1＝ContourPlot3D[Cos[XL]Sin[YL]+Con[YL]Sin[ZL]+Cos[ZL]Sin[XL]＝0，{x，-Pi，Pi}，{y，-Pi，Pi}，{z，-Pi，Pi}]
wherein, p1 is a heat exchange function, contourPlot3D is a three-dimensional space transformation function, x is an x-axis coordinate value, y is a y-axis coordinate value, z is a z-axis coordinate value, and b is a variable parameter;
s4-2: carrying out three-dimensional modeling on the heat exchange function obtained in the step S4-1, and constructing unit bodies corresponding to V1, V2, V3 and V4 by utilizing V1, V2, V3 and V4 obtained in the step S3;
s4-3, adjusting the unit bodies corresponding to V1, V2, V3 and V4 obtained in the step S4-2, and obtaining the optimal unit bodies corresponding to V1, V2, V3 and V4 meeting preset conditions, thereby realizing verification analysis of the flow velocity distribution conditions of the heat source and the cold source in the plate-fin heat exchanger and the temperature distribution conditions of the heat exchange surfaces of the heat source and the cold source; the optimal unit bodies corresponding to V1, V2, V3 and V4 are obtained, so that the resistance of fluid in the plate-fin heat exchanger can be reduced, the productivity of the air separation device is improved, and finally the optimal unit bodies corresponding to V1, V2, V3 and V4 are led into a Grasshopper parametric modeling platform, and an aggregate unit consisting of the optimal unit bodies corresponding to V1, V2, V3 and V4 is obtained through a rapid iteration array; the V1, V2, V3 and V4 of the invention are regulated more rapidly by using a Grasshopper parameterized modeling platform.
Example 6
A heat exchange optimization method of a high-pressure plate-fin heat exchanger of an air separation device comprises the following steps:
s1, establishing a model
Modeling the internal flow passage structure of the plate-fin heat exchanger by adopting FLUENT software to obtain a dynamic model of the internal flow passage structure of the plate-fin heat exchanger;
s2, parameter acquisition
S2-1, according to the dynamic model of the internal flow channel structure of the plate-fin heat exchanger obtained in the step S1, adopting FLUENT software to simulate and analyze flow state parameters A of medium-pressure air and high-pressure air when the medium-pressure air and the high-pressure air flow pass through the inside of the plate-fin heat exchanger, and obtaining the flow velocity distribution condition of a heat source in the plate-fin heat exchanger and the temperature distribution condition of a heat exchange surface of the heat source;
s2-2, simulating and analyzing flow state parameters B of pure nitrogen and high-pressure oxygen when the pure nitrogen and the high-pressure oxygen flow through the inside of the plate-fin heat exchanger by adopting FLUENT software to obtain the flow velocity distribution condition of a cold source inside the plate-fin heat exchanger and the temperature distribution condition of a heat exchange surface of the cold source;
s3, heat exchange loss simulation
Substituting the flow state parameter A obtained in the step S2-1 and the flow state parameter B obtained in the step S2-2 into a dynamic model of an internal flow channel structure of the plate-fin heat exchanger obtained in the step S1, and taking the opening V1 of a valve for controlling medium-pressure air flow, the opening V2 of a valve for controlling high-pressure air flow, the opening V3 of a valve for controlling pure nitrogen flow and the opening V4 of a valve for controlling high-pressure oxygen as parameter variables, wherein the parameter variables also comprise medium-pressure air temperature T1, high-pressure air temperature T2, pure nitrogen temperature T3 and high-pressure oxygen temperature T4; the influence of the flow and the temperature of the fluid entering the plate-fin heat exchanger on the heat exchange loss of the air separation device is simulated, so that the plate-fin heat exchanger can be subjected to multi-dimensional adjustment in the use process; simulating the interrelationship of a flow state parameter A, a flow state parameter B, V1, V2, V3, V4, T1, T2, T3 and T4 by using a dynamic model of an internal flow channel structure of the plate-fin heat exchanger, and calculating the optimal solution corresponding to V1, V2, V3, V4, T1, T2, T3 and T4 affecting the heat exchange loss of the plate-fin heat exchanger by using a genetic algorithm; finally, according to the adjustable parameter ranges of the flow state parameters A and B, carrying out up-down floating adjustment on the flow state parameters A and B in the adjustable parameter ranges of the flow state parameters A and B, and then verifying the response speeds of V1, V2, V3, V4, T1, T2, T3 and T4 by utilizing the variable quantities of the flow state parameters A and B; by verifying the response speeds of V1, V2, V3, V4, T1, T2, T3 and T4, the problem of capacity reduction of the air separation device caused by heat exchange loss can be avoided;
s4, analysis and verification
And (3) substituting the optimal solutions corresponding to V1, V2, V3, V4, T1, T2, T3 and T4 obtained in the step (S3) into the dynamic model of the internal flow channel structure of the plate-fin heat exchanger obtained in the step (S1), and performing verification analysis on the distribution condition of the flow velocity of the heat source and the cold source in the plate-fin heat exchanger and the distribution condition of the temperature of the heat exchange surfaces of the heat source and the cold source.
Example 7
A heat exchange optimization method of a high-pressure plate-fin heat exchanger of an air separation device comprises the following steps:
s1, establishing a model
Modeling the internal flow passage structure of the plate-fin heat exchanger by adopting FLUENT software to obtain a dynamic model of the internal flow passage structure of the plate-fin heat exchanger;
s2, parameter acquisition
S2-1, according to the dynamic model of the internal flow channel structure of the plate-fin heat exchanger obtained in the step S1, adopting FLUENT software to simulate and analyze flow state parameters A of medium-pressure air and high-pressure air when the medium-pressure air and the high-pressure air flow pass through the inside of the plate-fin heat exchanger, and obtaining the flow velocity distribution condition of a heat source in the plate-fin heat exchanger and the temperature distribution condition of a heat exchange surface of the heat source;
s2-2, simulating and analyzing flow state parameters B of pure nitrogen and high-pressure oxygen when the pure nitrogen and the high-pressure oxygen flow through the inside of the plate-fin heat exchanger by adopting FLUENT software to obtain the flow velocity distribution condition of a cold source inside the plate-fin heat exchanger and the temperature distribution condition of a heat exchange surface of the cold source;
s3, heat exchange loss simulation
Substituting the flow state parameter A obtained in the step S2-1 and the flow state parameter B obtained in the step S2-2 into a dynamic model of the internal flow channel structure of the plate-fin heat exchanger obtained in the step S1, taking the opening V1 of a valve for controlling medium-pressure air flow, the opening V2 of a valve for controlling high-pressure air flow, the opening V3 of a valve for controlling pure nitrogen flow and the opening V4 of a valve for controlling high-pressure oxygen as parameter variables, and simulating the interrelationship of the flow state parameter A, the flow state parameter B, the V1, the V2, the V3 and the V4 by using the dynamic model of the internal flow channel structure of the plate-fin heat exchanger to calculate the corresponding optimal solution of the V1, the V2, the V3 and the V4 affecting the heat exchange loss of the plate-fin heat exchanger;
s4, analysis and verification
Substituting the optimal solutions corresponding to V1, V2, V3 and V4 obtained in the step S3 into the dynamic model of the internal flow channel structure of the plate-fin heat exchanger obtained in the step S1, and performing verification analysis on the distribution conditions of the flow velocity of the heat source and the cold source in the plate-fin heat exchanger and the distribution conditions of the temperature of the heat exchange surfaces of the heat source and the cold source; finally, carrying out reflux nitrogen temperature simulation on an internal flow channel structure dynamic model of the plate-fin heat exchanger by utilizing LUENT software, calculating a theoretical value of the reflux nitrogen temperature of the plate-fin heat exchanger according to the optimal solution corresponding to V1, V2, V3 and V4, and calculating a heat exchange coefficient of the plate-fin heat exchanger according to the theoretical value; the heat exchange coefficient of the plate-fin heat exchanger is calculated, so that the improvement of the plate-fin heat exchanger is facilitated, and the heat exchange effect is improved.
Example 8
A heat exchange optimization method of a high-pressure plate-fin heat exchanger of an air separation device comprises the following steps:
s1, establishing a model
S1-1, modeling a plate-fin heat exchanger by using Gambit software, selecting a grid type, setting grid boundary conditions, and deriving a grid file;
s1-2, importing the grid file obtained in the step S1-1 into an EDEM, and then filling the internal flow passage structural parameters of the plate-fin heat exchanger to obtain a plate-fin heat exchanger internal flow passage structural parameter model;
s1-3, introducing the parameter model of the internal flow passage structure of the plate-fin heat exchanger obtained in the step S1-2 into FLUENT, and obtaining a dynamic model of the internal flow passage structure of the plate-fin heat exchanger; dividing a dynamic model of an internal flow channel structure of the plate-fin heat exchanger by adopting a structured grid; the structural grid is adopted to divide the dynamic model of the internal flow channel structure of the plate-fin heat exchanger, so that the optimization calculation time of each parameter is shortened, and the modification of complex variables in the model is facilitated; by the modeling method, a reliable research environment can be provided for the optimal design of the heat exchange effect of the plate-fin heat exchanger;
s2, parameter acquisition
S2-1, according to the dynamic model of the internal flow channel structure of the plate-fin heat exchanger obtained in the step S1, adopting FLUENT software to simulate and analyze flow state parameters A of medium-pressure air and high-pressure air when the medium-pressure air and the high-pressure air flow pass through the inside of the plate-fin heat exchanger, and obtaining the flow velocity distribution condition of a heat source in the plate-fin heat exchanger and the temperature distribution condition of a heat exchange surface of the heat source;
s2-2, simulating and analyzing flow state parameters B of pure nitrogen and high-pressure oxygen when the pure nitrogen and the high-pressure oxygen flow through the inside of the plate-fin heat exchanger by adopting FLUENT software to obtain the flow velocity distribution condition of a cold source inside the plate-fin heat exchanger and the temperature distribution condition of a heat exchange surface of the cold source; the flow state parameters A and the flow state parameters B comprise temperature variable parameters, pressure variable parameters and flow variable parameters, and the heat exchange effect of the plate-fin heat exchanger is designed in an intervention way through setting multiple parameters, so that the optimization effect of the invention is better outstanding;
s3, heat exchange loss simulation
Substituting the flow state parameter A obtained in the step S2-1 and the flow state parameter B obtained in the step S2-2 into a dynamic model of an internal flow channel structure of the plate-fin heat exchanger obtained in the step S1, and taking the opening V1 of a valve for controlling medium-pressure air flow, the opening V2 of a valve for controlling high-pressure air flow, the opening V3 of a valve for controlling pure nitrogen flow and the opening V4 of a valve for controlling high-pressure oxygen as parameter variables, wherein the parameter variables also comprise medium-pressure air temperature T1, high-pressure air temperature T2, pure nitrogen temperature T3 and high-pressure oxygen temperature T4; the influence of the flow and the temperature of the fluid entering the plate-fin heat exchanger on the heat exchange loss of the air separation device is simulated, so that the plate-fin heat exchanger can be subjected to multi-dimensional adjustment in the use process; simulating the interrelationship of a flow state parameter A, a flow state parameter B, V1, V2, V3, V4, T1, T2, T3 and T4 by using a dynamic model of an internal flow channel structure of the plate-fin heat exchanger, and calculating the optimal solution corresponding to V1, V2, V3, V4, T1, T2, T3 and T4 affecting the heat exchange loss of the plate-fin heat exchanger by using a particle swarm algorithm; finally, according to the adjustable parameter ranges of the flow state parameters A and B, carrying out up-down floating adjustment on the flow state parameters A and B in the adjustable parameter ranges of the flow state parameters A and B, and then verifying the response speeds of V1, V2, V3, V4, T1, T2, T3 and T4 by utilizing the variable quantities of the flow state parameters A and B; by verifying the response speeds of V1, V2, V3, V4, T1, T2, T3 and T4, the problem of capacity reduction of the air separation device caused by heat exchange loss can be avoided;
s4, analysis and verification
S4-1, constructing a heat exchange function based on a sine and cosine function according to a dynamic model of an internal flow channel structure of the plate-fin heat exchanger; the expression of the heat exchange function is as follows:
XL＝x-bCos[(x/2)^2]
YL＝y-bCos[(y/2)^2]
ZL＝z-bCos[(z/2)^2]
p1＝ContourPlot3D[Cos[XL]Sin[YL]+Con[YL]Sin[ZL]+Cos[ZL]Sin[XL]＝0，{x，-Pi，Pi}，{y，-Pi，Pi}，{z，-Pi，Pi}]
wherein, p1 is a heat exchange function, contourPlot3D is a three-dimensional space transformation function, x is an x-axis coordinate value, y is a y-axis coordinate value, z is a z-axis coordinate value, and b is a variable parameter;
s4-2: performing three-dimensional modeling on the heat exchange function obtained in the step S4-1, and constructing corresponding unit bodies of V1, V2, V3 and V4;
s4-3, adjusting the unit bodies corresponding to V1, V2, V3 and V4 obtained in the step S4-2 to obtain the optimal unit bodies corresponding to V1, V2, V3 and V4 meeting preset conditions, and finally verifying and analyzing the distribution condition of the flow velocity of the heat source and the cold source in the plate-fin heat exchanger and the distribution condition of the temperature of the heat exchange surfaces of the heat source and the cold source; the optimal unit bodies corresponding to V1, V2, V3, V4, T1, T2, T3 and T4 are obtained, so that the resistance of the fluid in the plate-fin heat exchanger can be reduced, and the productivity of the air separation device is improved; the optimal unit bodies corresponding to V1, V2, V3, V4, T1, T2, T3 and T4 are imported into a Grasshopper parametric modeling platform, and an aggregate unit consisting of the optimal unit bodies corresponding to V1, V2, V3, V4, T1, T2, T3 and T4 is obtained through a rapid iteration array; the Grasshopper parameterized modeling platform is utilized, so that V1, V2, V3, V4, T1, T2, T3 and T4 are regulated more rapidly; carrying out reflux nitrogen temperature simulation on an internal flow channel structure dynamic model of the plate-fin heat exchanger by utilizing LUENT software, and finally calculating a theoretical value of the reflux nitrogen temperature of the plate-fin heat exchanger according to the optimal solution corresponding to V1, V2, V3, V4, T1, T2, T3 and T4, and calculating a heat exchange coefficient of the plate-fin heat exchanger according to the theoretical value; the heat exchange coefficient of the plate-fin heat exchanger is calculated, so that the improvement of the plate-fin heat exchanger is facilitated, and the heat exchange effect is improved.
Application example
The heat exchange optimization is carried out on the plate-fin heat exchanger of a large chemical fertilizer air separation workshop in the south of China by using the method of the embodiments 1-8, after the optimization is finished, the temperature difference of the reflux nitrogen of the plate-fin heat exchanger and the liquid oxygen yield increase of the air separation device are measured, and the measurement results are shown in the table 1:
table 1 effects of various examples on plate-fin heat exchanger reflux nitrogen temperature differential and air separation plant liquid oxygen production increase
As can be seen from the data in table 1: compared with the distinguishing scheme of the embodiment 1, the embodiment 2 can provide a reliable research environment for the optimal design of the heat exchange effect of the plate-fin heat exchanger, thereby improving the heat exchange optimal effect of the plate-fin heat exchanger and the productivity of the air separation unit;
compared with the embodiment 1, the embodiment 3 is beneficial to shortening the optimization calculation time of each parameter and modifying complex variables in the model by dividing the dynamic model of the internal flow channel structure of the plate-fin heat exchanger by adopting the structured grid;
in example 4, compared with example 1, by obtaining the optimal unit bodies corresponding to V1, V2, V3, and V4, the resistance of the fluid in the plate-fin heat exchanger can be reduced, thereby improving the capacity of the air separation unit; the heat exchange effect of the plate-fin heat exchanger is interfered and designed by arranging multiple parameters, so that the optimization effect of the invention is better outstanding;
example 5 compared with example 1, by introducing the optimal unit bodies corresponding to V1, V2, V3, V4 into a Grasshopper parametric modeling platform, and by rapidly iterating the array, the aggregate unit consisting of the optimal unit bodies corresponding to V1, V2, V3, V4 is obtained, so that V1, V2, V3, V4 of the present invention is more rapid in adjustment;
in example 6, compared with example 1, by verifying the response speeds corresponding to V1, V2, V3, and V4, the problem of capacity drop of the air separation plant due to heat exchange loss can be avoided; the influence of the flow and the temperature of the fluid entering the plate-fin heat exchanger on the heat exchange loss of the air separation device is simulated, so that the plate-fin heat exchanger can be subjected to multi-dimensional adjustment in the use process;
compared with the embodiment 1, the embodiment 7 utilizes the LUENT software to simulate the reflux nitrogen temperature of the dynamic model of the internal flow channel structure of the plate-fin heat exchanger, calculates the heat exchange coefficient of the plate-fin heat exchanger, and is beneficial to improving the plate-fin heat exchanger and improving the heat exchange effect;
compared with the embodiment 1-7, in the heat exchange optimization process of the plate-fin heat exchanger, the reflux nitrogen temperature difference of the plate-fin heat exchanger is restrained by optimizing the dynamic model of the internal flow channel structure of the plate-fin heat exchanger, so that the liquid oxygen yield of the air separation device is effectively increased.
Claims (6)
1. The heat exchange optimization method for the high-pressure plate-fin heat exchanger of the air separation device is characterized by comprising the following steps of:
s1, establishing a model
Modeling the internal flow passage structure of the plate-fin heat exchanger by adopting FLUENT software to obtain a dynamic model of the internal flow passage structure of the plate-fin heat exchanger;
s2, parameter acquisition
S2-1, according to the dynamic model of the internal flow channel structure of the plate-fin heat exchanger obtained in the step S1, adopting FLUENT software to simulate and analyze flow state parameters A of medium-pressure air and high-pressure air when the medium-pressure air and the high-pressure air flow pass through the inside of the plate-fin heat exchanger, and obtaining the flow velocity distribution condition of a heat source in the plate-fin heat exchanger and the temperature distribution condition of a heat exchange surface of the heat source;
s2-2, simulating and analyzing flow state parameters B of pure nitrogen and high-pressure oxygen when the pure nitrogen and the high-pressure oxygen flow through the inside of the plate-fin heat exchanger by adopting FLUENT software to obtain the flow velocity distribution condition of a cold source inside the plate-fin heat exchanger and the temperature distribution condition of a heat exchange surface of the cold source;
s3, heat exchange loss simulation
Substituting the flow state parameter A obtained in the step S2-1 and the flow state parameter B obtained in the step S2-2 into a dynamic model of the internal flow channel structure of the plate-fin heat exchanger obtained in the step S1, taking the opening V1 of a valve for controlling medium-pressure air flow, the opening V2 of a valve for controlling high-pressure air flow, the opening V3 of a valve for controlling pure nitrogen flow and the opening V4 of a valve for controlling high-pressure oxygen as parameter variables, and simulating the interrelationship of the flow state parameter A, the flow state parameter B, the V1, the V2, the V3 and the V4 by using the dynamic model of the internal flow channel structure of the plate-fin heat exchanger to calculate the corresponding optimal solution of the V1, the V2, the V3 and the V4 affecting the heat exchange loss of the plate-fin heat exchanger;
s4, analysis and verification
Substituting the optimal solutions corresponding to V1, V2, V3 and V4 obtained in the step S3 into the dynamic model of the internal flow channel structure of the plate-fin heat exchanger obtained in the step S1, and performing verification analysis on the distribution conditions of the flow velocity of the heat source and the cold source in the plate-fin heat exchanger and the distribution conditions of the temperature of the heat exchange surfaces of the heat source and the cold source;
in the steps S2-1 and S2-2, the flow state parameters A and B comprise temperature variable parameters, pressure variable parameters and flow variable parameters;
the specific operation of step S4 is as follows:
s4-1, constructing a heat exchange function based on a sine and cosine function according to a dynamic model of an internal flow channel structure of the plate-fin heat exchanger; the expression of the heat exchange function is as follows:
XL＝x-bCos[(x/2)^2]
YL＝y-bCos[(y/2)^2]
ZL＝z-bCos[(z/2)^2]
p1＝ContourPlot3D[Cos[XL]Sin[YL]+Con[YL]Sin[ZL]+Cos[ZL]Sin[XL]＝0，{x，-Pi，Pi}，{y，-Pi，Pi}，{z，-Pi，Pi}]
wherein, p1 is a heat exchange function, contourPlot3D is a three-dimensional space transformation function, x is an x-axis coordinate value, y is a y-axis coordinate value, z is a z-axis coordinate value, and b is a variable parameter;
s4-2: performing three-dimensional modeling on the heat exchange function obtained in the step S4-1, and constructing corresponding unit bodies of V1, V2, V3 and V4;
s4-3, adjusting the unit bodies corresponding to V1, V2, V3 and V4 obtained in the step S4-2 to obtain the optimal unit bodies corresponding to V1, V2, V3 and V4 meeting preset conditions, and finally verifying and analyzing the distribution condition of the flow velocity of the heat source and the cold source in the plate-fin heat exchanger and the distribution condition of the temperature of the heat exchange surfaces of the heat source and the cold source;
the specific operation of step S1 is:
s1-1, modeling a plate-fin heat exchanger by using Gambit software, selecting a grid type, setting grid boundary conditions, and deriving a grid file;
s1-2, importing the grid file obtained in the step S1-1 into an EDEM, and then filling the internal flow passage structural parameters of the plate-fin heat exchanger to obtain a plate-fin heat exchanger internal flow passage structural parameter model;
s1-3, introducing the parameter model of the internal flow passage structure of the plate-fin heat exchanger obtained in the step S1-2 into FLUENT, and obtaining the dynamic model of the internal flow passage structure of the plate-fin heat exchanger.
2. The heat exchange optimization method for the high-pressure plate-fin heat exchanger of the air separation device according to claim 1, wherein the optimal unit bodies corresponding to V1, V2, V3 and V4 are led into a Grasshopper parameterized modeling platform, and an aggregate unit consisting of the optimal unit bodies corresponding to V1, V2, V3 and V4 is obtained through a rapid iteration array.
3. The method for optimizing heat exchange of high-pressure plate-fin heat exchanger of air separation device according to claim 1, wherein after step S1 is completed, the dynamic model of the internal flow channel structure of the plate-fin heat exchanger is divided by adopting a structured grid.
4. The heat exchange optimization method for the high-pressure plate-fin heat exchanger of the air separation device according to claim 1, wherein after the step S3 is completed, the flow state parameters a and B are adjusted in a floating manner up and down within the adjustable parameter ranges of the flow state parameters a and B, and then response speeds corresponding to V1, V2, V3 and V4 are verified by using the variable amounts of the flow state parameters a and B.
5. The method for optimizing heat exchange of high-pressure plate-fin heat exchanger of air separation unit according to claim 1, wherein in step S3, the algorithm used for calculation is one of genetic algorithm, particle swarm algorithm, and ant colony algorithm.
6. The method according to claim 1, wherein in step S3, the parameter variables further include a medium pressure air temperature T1, a high pressure air temperature T2, a pure nitrogen temperature T3, and a high pressure oxygen temperature T4.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111554014.4A CN114235891B (en) | 2021-12-17 | 2021-12-17 | Heat exchange optimization method for high-pressure plate-fin heat exchanger of air separation device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111554014.4A CN114235891B (en) | 2021-12-17 | 2021-12-17 | Heat exchange optimization method for high-pressure plate-fin heat exchanger of air separation device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114235891A CN114235891A (en) | 2022-03-25 |
| CN114235891B true CN114235891B (en) | 2024-03-19 |
Family
ID=80758306
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111554014.4A Active CN114235891B (en) | 2021-12-17 | 2021-12-17 | Heat exchange optimization method for high-pressure plate-fin heat exchanger of air separation device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114235891B (en) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103673695A (en) * | 2013-12-25 | 2014-03-26 | 山东大学 | Plate heat exchanger and method for calculating criterion numeral of plate heat exchanger |
| CN107391807A (en) * | 2017-06-28 | 2017-11-24 | 西安交通大学 | Plate-fin heat exchanger heat transfer flow performance value analogy method based on transient technology |
| CN107704705A (en) * | 2017-10-24 | 2018-02-16 | 广东唯金智能环境科技有限公司 | A kind of optimization method of the air source hot pump water heater condenser pipe based on kriging models |
| CN208075213U (en) * | 2018-04-10 | 2018-11-09 | 北京远大天益生态建筑设计院有限公司 | Air-source consumption reduction thermostatic equipment |
| CN109030039A (en) * | 2017-12-26 | 2018-12-18 | 上海齐耀动力技术有限公司 | A kind of regenerator performance detecting system and method |
| CN110993034A (en) * | 2019-11-26 | 2020-04-10 | 华南理工大学 | Simulation method of CFD-based cyclohexane non-catalytic oxidation reactor |
| CN111723472A (en) * | 2020-05-29 | 2020-09-29 | 同济大学 | Heat exchanger structure optimization method based on hot melt type gas-liquid two-phase heat exchange structure |
| CN112084591A (en) * | 2020-09-03 | 2020-12-15 | 西安电子科技大学 | Radiator cooling channel design method based on three-dimensional topological optimization |
-
2021
- 2021-12-17 CN CN202111554014.4A patent/CN114235891B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103673695A (en) * | 2013-12-25 | 2014-03-26 | 山东大学 | Plate heat exchanger and method for calculating criterion numeral of plate heat exchanger |
| CN107391807A (en) * | 2017-06-28 | 2017-11-24 | 西安交通大学 | Plate-fin heat exchanger heat transfer flow performance value analogy method based on transient technology |
| CN107704705A (en) * | 2017-10-24 | 2018-02-16 | 广东唯金智能环境科技有限公司 | A kind of optimization method of the air source hot pump water heater condenser pipe based on kriging models |
| CN109030039A (en) * | 2017-12-26 | 2018-12-18 | 上海齐耀动力技术有限公司 | A kind of regenerator performance detecting system and method |
| CN208075213U (en) * | 2018-04-10 | 2018-11-09 | 北京远大天益生态建筑设计院有限公司 | Air-source consumption reduction thermostatic equipment |
| CN110993034A (en) * | 2019-11-26 | 2020-04-10 | 华南理工大学 | Simulation method of CFD-based cyclohexane non-catalytic oxidation reactor |
| CN111723472A (en) * | 2020-05-29 | 2020-09-29 | 同济大学 | Heat exchanger structure optimization method based on hot melt type gas-liquid two-phase heat exchange structure |
| CN112084591A (en) * | 2020-09-03 | 2020-12-15 | 西安电子科技大学 | Radiator cooling channel design method based on three-dimensional topological optimization |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114235891A (en) | 2022-03-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111400958B (en) | ANSYS-based simulation method for oxygen content distribution of high-temperature carbonization furnace | |
| CN103699744B (en) | Wind power master control cabinet heat management analysis and optimization method based on finite element method | |
| CN111400934B (en) | Method for simulating oxygen content distribution of low-temperature carbonization furnace based on WORKBENCH | |
| CN112446177B (en) | Simulation method for heat insulation performance of external heat insulation material of high-temperature carbonization furnace | |
| CN112231945B (en) | Power battery system thermal diffusion joint simulation method based on star CCM + and Amesim | |
| CN111680432B (en) | Low-temperature carbonization furnace multi-coupling field stress distribution simulation method based on WORKBENCH | |
| CN112528569A (en) | Industrial heating furnace temperature field prediction method based on digital twinning | |
| CN114235891B (en) | Heat exchange optimization method for high-pressure plate-fin heat exchanger of air separation device | |
| WO2023087483A1 (en) | Greenhouse environment optimization design method considering crop production benefits and energy consumption reduction | |
| CN115034138A (en) | Boiler temperature field prediction method combining computational fluid dynamics and deep learning | |
| CN107506516A (en) | A kind of communications equipment room flow field model is established and analysis method and system | |
| WO2021207953A1 (en) | Workbench-based low-temperature carbonization furnace oxygen content distribution simulation method | |
| CN113158355B (en) | All-condition optimization design method for low-temperature liquid expander | |
| CN112270109B (en) | Method for simulating heating performance of graphite rod in high-temperature carbonization furnace | |
| CN111881492B (en) | CFD (computational fluid dynamics) method-based steam turbine valve steam distribution management function generation method | |
| CN110489912B (en) | Method for simulating numerical value of solar cross-season soil heat storage layered slice | |
| CN117113553A (en) | CFD-based gearbox heat dissipation simulation method | |
| CN115062560A (en) | Turbine multi-scale flow reduced-order coupling method based on machine learning | |
| CN111780491B (en) | Energy-saving air separation cooling system and method | |
| CN113420364B (en) | Electrophoresis process vehicle body structure deformation simulation method based on fluid-solid coupling | |
| CN105013837A (en) | Control method for improving high-speed wire rod air cooling uniformity | |
| CN110807280B (en) | Heat collecting device suitable for conduction, convection and radiation simultaneously | |
| CN107194018A (en) | The emulation mode of mould chilling temperature in thermoforming | |
| WO2022011724A1 (en) | Workbench-based method for simulating multi-coupling field stress distribution of low-temperature carbonization furnace | |
| Liu et al. | Numerical simulation on flow and heat transfer performance of air-cooler for a natural gas storage compressor unit |
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 |