CN112287515B - Superconducting energy pipeline overall quench prediction method based on multi-field coupling model - Google Patents
Superconducting energy pipeline overall quench prediction method based on multi-field coupling model Download PDFInfo
- Publication number
- CN112287515B CN112287515B CN202010986813.8A CN202010986813A CN112287515B CN 112287515 B CN112287515 B CN 112287515B CN 202010986813 A CN202010986813 A CN 202010986813A CN 112287515 B CN112287515 B CN 112287515B
- Authority
- CN
- China
- Prior art keywords
- temperature
- current
- energy pipeline
- superconducting
- superconducting energy
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/14—Pipes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
A superconducting energy pipeline overall quench prediction method based on a multi-field coupling model comprises the steps of firstly measuring critical currents of high-temperature superconducting tapes at different temperatures, and fitting to obtain a critical current I c Data curve I as a function of temperature c (T); then establishing an electromagnetic field model and calculating the current distribution I of the copper framework former (x) (ii) a Then establishing a flow field analysis model according to the design working condition of the low-temperature fuel, and calculating to obtain the axial distribution h (x) of the heat exchange coefficient of the low-temperature fuel; then according to the current distribution I former (x) Establishing a thermal analysis model, and calculating to obtain axial temperature distribution T (x) of the superconducting energy pipeline; updating model parameters according to the axial temperature distribution T (x), and repeating until the set time T is reached total When the temperature is higher than the preset temperature, the calculation is finished, and a curve T (T) of the change of the highest temperature along with the time is extracted; finally, judging whether the maximum temperature change T (T) of the superconducting energy pipeline meets the overall quench judgment criterion of the energy pipeline; the invention improves the running safety of the superconducting energy pipeline.
Description
Technical Field
The invention relates to the field of natural gas and energy transportation, in particular to a superconducting energy pipeline overall quench prediction method based on a multi-field coupling model.
Background
Due to uneven distribution of energy resources, fossil energy and electric energy are required to be transported in a long distance in order to meet energy use requirements of various regions. The high-temperature superconducting cable has the advantages of zero impedance, high current-carrying density, small space occupation ratio, environmental friendliness and the like, and can realize long-distance transmission with zero loss. On the other hand, long distance transmission of fossil energy such as hydrogen, natural gas, shale gas, etc. usually takes a liquid form to achieve low transmission loss and high transmission capacity; the two are jointly conveyed, the liquid fuel is utilized to provide a low-temperature operation environment for the superconducting cable, a superconducting energy pipeline is assembled, the integrated conveying of electric power and fossil energy is realized, and the energy conveying efficiency can be greatly improved.
The superconducting energy pipeline is widely researched and discussed in recent years, the concept of an LNG/power mixed transmission superconducting energy pipeline is proposed in Chinese patent CN201210118316.1, a high-temperature superconducting cable is cooled by using Liquefied Natural Gas (LNG), a cable and natural gas combined transmission pipeline system is designed, and the energy loss of the combined transmission system is greatly reduced compared with that of the two independent transmission systems. Chinese patent CN201810804587.X provides a superconducting energy pipeline for low-temperature fuel conduction cooling, which can realize safe and efficient combined delivery of low-temperature liquid fuel and electric energy. However, none of the above solutions consider the stability safety issue of the superconducting energy pipeline in case of short circuit fault or thermal disturbance. Under a certain condition, the short-circuit fault can cause the quench of the superconducting cable, the quench is divided into local quench and integral quench, and when the fault disturbance impact exceeds a certain range, the superconducting cable can generate the integral quench. For a superconducting energy pipeline, the phenomenon of quench involves a complex coupling relationship among an electric field, a magnetic field, a flow field and a temperature field: the operating state of the cryogenic fuel affects the operating state of the cable, which also affects the flow state of the cryogenic fuel.
At present, various simulation models aiming at multi-physical fields have been proposed in academia, for example, chinese patent CN201910997413.4, etc., and the main method is to utilize commercial software to calculate the electromagnetic field, flow field and temperature field distribution separately in steps, so that the strong coupling relationship among the multi-physical fields is not embodied, and the calculation accuracy of these methods applied to the superconducting energy pipeline will be greatly reduced. In addition, the methods for detecting the quench of the superconducting cable mainly include a voltage method, a current method, wen Shengfa, a pressure method, a flow velocity method and an ultrasonic wave method, but these methods can only detect a local quench after the quench occurs, and cannot predict the entire quench.
Disclosure of Invention
In order to overcome the defects of the prior art, the invention aims to provide a superconducting energy pipeline overall quench prediction method based on a multi-field coupling model, wherein an electromagnetic heat flow multi-field coupling model is established to reflect the mutual influence relation among an electromagnetic field, a flow field and a temperature field, the superconducting energy pipeline is subjected to overall quench prediction, and the prediction is carried out before the overall quench of the energy pipeline occurs, so that the operation configuration is reduced in advance, the quench is prevented, and the safety and reliability of system operation are improved.
In order to achieve the purpose, the invention adopts the technical scheme that:
a superconducting energy pipeline overall quench prediction method based on a multi-field coupling model comprises the following steps:
(1) Measuring the critical current of the high-temperature superconducting tape at different temperatures in a laboratory, and fitting to obtain the critical current I c Data curve I as a function of temperature c (T); designing structural parameters and operation conditions of the superconducting energy pipeline; according to the operating temperature T f Calculating the physical properties of the low-temperature fuel and the cable material by calling NIST;
real-time acquisition of cable running current I in superconducting energy pipeline op ；
(2) According to the structural parameters and the critical current I of the superconducting energy pipeline in the step (1) c (T) establishing an electromagnetic field model, and calculating to obtain the current distribution I of the copper skeleton in the superconducting energy pipeline former (x)；
(3) Firstly, according to the current distribution I obtained in the step (2) former (x) Generating an electromagnetic heat exchange component h EM Thereby realizing the coupling effect between the electromagnetic field and the flow field; establishing a flow field analysis model according to the operation condition and the physical properties of the low-temperature fuel in the step (1), and calculating to obtain the axial distribution h (x) of the heat exchange coefficient of the low-temperature fuel;
(4) According to the current distribution I obtained in the step (2) former (x) And (4) establishing a thermal analysis model according to the heat exchange coefficient distribution h (x) obtained in the step (3), and calculating to obtain axial temperature distribution T (x) of the superconducting energy pipeline;
(5) Updating model parameters according to the axial temperature distribution T (x) obtained in the step (4) so as to realize the coupling effect between the temperature field and other physical fields, calculating the next time step, and repeating the steps (2) to (5) until the set time T is reached total When the temperature is higher than the preset temperature, the calculation is finished, and a curve T (T) of the change of the highest temperature along with the time is extracted;
(6) And (4) judging whether the change curve T (T) of the highest temperature obtained in the step (5) along with time meets the judgment criterion of the overall quench of the superconducting energy pipeline, so as to predict whether the superconducting energy pipeline is quenched integrally.
In the step (2), the electromagnetic field model is a conservation model based on kirchhoff's law, and the equation is as follows:
I op ＝I f +∑I i
wherein U is terminal voltage, I i Is a flowing current, R f Is the resistance of the copper skeleton, R 1 -R m Is the resistance of the respective superconducting tapes, I f Is the current on the copper skeleton, I 1 -I m Is the current of each superconducting tape; i is op Is the running current, L f Is self-inductance of the copper skeleton, L 1 -L m Is the self-inductance of the respective superconducting layer, M i,j Is the mutual inductance between superconducting tapes i and j;
the current-voltage characteristics of the superconducting strip material follow an n-exponential power relation as follows:
wherein E is 0 Is the quench voltage, l is the strip length, I c Is the critical current and n is the power exponent.
The flow field analysis model in the step (3) is established by adopting an empirical correlation method, and specifically comprises the following steps:
judging whether the mass gas fraction a of the low-temperature fuel at the current calculation node is equal to 0 or 1, if so, calculating the heat exchange coefficient of the low-temperature fuel at the node by adopting a Dittus-Boelter correlation formula, wherein the calculation formula is h = h fc +h EM Where h is the heat transfer coefficient at the current node, h fc Is single-phase forced convection heat exchangeComponent, h EM Is the electromagnetic heat exchange component;
if 0<a<1, calculating the nucleate boiling initiation temperature T by using an empirical correlation ONB Modal boiling onset temperature T MFB And maximum heat flux density of nucleate boiling q CHF ；
If the current node temperature T w Greater than T MFB If the heat exchange coefficient of the low-temperature fuel at the node is calculated by adopting a Darr correlation formula, the calculation formula is h = h b +h tc +h dw +h EM Wherein h is b 、h tc 、h dw 、h EM Respectively a modal boiling component, a convective heat transfer component, a heat flow component and an electromagnetic heat transfer component;
if the current node temperature T w Less than T ONB If so, calculating the heat exchange coefficient of the low-temperature fuel at the node by adopting a Dittus-Boelter correlation formula;
if the current node temperature T w Between T and MFB and T ONB Then comparing the current node heat flux density q w And maximum heat flux q of nucleate boiling CHF The size relationship of (1):
if q is w Less than q CHF If the heat exchange coefficient of the low-temperature fuel at the node is calculated by adopting a Chen correlation formula, the calculation formula is h = Eh fc +Sh p +h EM Wherein E is a forced convection heat transfer factor, S is a nucleate boiling heat transfer factor, h p Is the nucleate boiling heat transfer component, h EM Is an electromagnetic heat exchange component;
if q is w Greater than q CHF And calculating the heat exchange coefficient of the low-temperature fuel at the node by adopting an interpolation method.
The thermal analysis model in the step (4) is a one-dimensional thermal diffusion model, and the equation is as follows:
wherein, A c And P c Is the cross-sectional area and perimeter of the copper skeleton, T is the temperature value at the current node of the superconducting energy pipeline, and rho is the density of the copper skeletonDegree c p Is the specific heat capacity, k is the thermal conductivity, p R Is the resistivity, I former Is the copper skeleton current, h is the heat transfer coefficient at the current node, T f Is the operating temperature.
The model parameters to be updated in step (5) include: density rho and specific heat capacity c of copper skeleton p Thermal conductivity k, resistivity ρ R Density of cryogenic fuel ρ l Specific heat capacity c pl And coefficient of thermal conductivity k l 。
Setting time t in the step (5) total The prediction result is more accurate when the time is more than 80 seconds.
The overall quench determination criterion in the step (6) is as follows: extracting a slope value at the last point of the curve T (T) of the change curve of the highest temperature along with time, and if the slope value is greater than zero, determining that the superconducting energy pipeline has integral quench; if the slope value is less than or equal to zero, the superconducting energy pipeline is considered not to be subjected to integral quench, and the pipeline system is safe to operate.
Compared with the prior art, the invention has at least the following beneficial effects:
the method adopts different heat exchange correlation formulas to calculate different flowing conditions of the low-temperature fuel, avoids the defect of inaccurate prediction of nucleate boiling by traditional numerical simulation, and combines an electromagnetic field model and flow and heat transfer analysis, so that the method can reflect complex multi-physical field effects of mutual influence of an electric field, a magnetic field, a flow field and a temperature field. The invention provides a judgment criterion of the whole quench of the energy pipeline, is applied to the prediction of the whole quench of the superconducting energy pipeline, has strong logic, is scientific and reliable, and can obviously improve the calculation accuracy, thereby improving the running safety of the superconducting energy pipeline.
Drawings
FIG. 1 is a flow chart of the present invention.
FIG. 2 is a diagram of a multi-field coupling model of the present invention.
FIG. 3 is a graph showing the time-dependent change T (T) of the maximum temperature in example 1 of the present invention.
FIG. 4 is a graph showing the time-dependent variation T (T) of the maximum temperature in example 2 of the present invention.
Detailed Description
The invention is further described with reference to the following figures and examples, which are not to be construed as limiting the invention in any way, and any limited number of modifications which can be made within the scope of the claims of the invention are still within the scope of the claims of the invention.
As shown in fig. 1, a superconducting energy pipeline overall quench prediction method based on a multi-field coupling model includes the following steps:
(1) Measuring the critical current of the high-temperature superconducting tape at different temperatures in a laboratory, and fitting to obtain the critical current I c Data curve I as a function of temperature c (T) =518-5.63T; designing the structural parameters and the operation conditions of the superconducting energy pipeline, which comprises the following steps: rated current I of superconducting cable e =3000A, operating temperature T of superconducting energy pipeline f = 80K-90K, superconducting cable winding coefficient a i ＝1,a j =1 superconducting cable winding pitch L p =362mm; according to the operating temperature T f And (3) calling NIST to calculate the physical properties of the low-temperature fuel and the cable material, wherein the method comprises the following steps: density of copper skeleton ρ =9000kg/m 3 Specific heat capacity c p = 380J/(kg · K), thermal conductivity K = 400W/(m · K), resistivity ρ R ＝2x10 -9 Omega.m; density of cryogenic fuel ρ l ＝486.55kg/m 3 Specific heat capacity c pl = 3077J/(kg. K), thermal conductivity K l ＝0.234W/(m·K)；
Real-time acquisition of running current I of cable in superconducting energy pipeline op ；
As shown in fig. 2, the multi-field coupling model consists of three submodels: an electromagnetic field model, a flow field analysis model and a thermal analysis model are respectively established by the subsequent steps;
(2) According to the structural parameters and the critical current I of the superconducting energy pipeline in the step (1) c (T) establishing an electromagnetic field model, and calculating to obtain the current distribution I of the copper skeleton in the superconducting energy pipeline former (x)；
The electromagnetic field model adopts a conservation model based on kirchhoff's law, and the equation is as follows:
I op ＝I f +∑I i
wherein U is terminal voltage, I i Is a flowing current, R f Is the resistance of the copper skeleton, R 1 -R m Is the resistance of the respective superconducting tapes, I f Is the current on the copper skeleton, I 1 -I m Is the current of each superconducting tape; I.C. A op Is the running current, L f Is self-inductance of the copper skeleton, L 1 -L m Is self-inductance of the respective superconducting layer, M i,j Is the mutual inductance between superconducting tapes i and j;
the current-voltage characteristics of the superconducting strip follow the relation of n exponential power, the power exponent n is 32, the quench voltage is E 0 =0.0001V/m, expression:
wherein E is 0 Is the quench voltage, l is the strip length, I c Is the critical current;
the critical current of the superconducting tape is obtained in the step (1), and the expression is as follows:
(3) Firstly, according to the current distribution I obtained in the step (2) former (x) Generating an electromagnetic heat exchange component h EM Thereby realizing the coupling effect between the electromagnetic field and the flow field; according to the operation condition and the physical property of the low-temperature fuel in the step (1), establishing a flow field analysis model, and calculating to obtain the axial distribution h (x) of the heat exchange coefficient of the low-temperature fuel;
the flow field analysis model is established by adopting an empirical correlation method, and specifically comprises the following steps:
judging whether the mass air content a of the low-temperature fuel at the current calculation node is equal to 0 or 1, if so, judging the heat exchange coefficient of the low-temperature fuel at the nodeCalculated by a Dittus-Boelter correlation formula, and the calculation formula is h = h fc +h EM Where h is the heat transfer coefficient at the current node, h fc Is a single-phase forced convection heat transfer component, h EM Is the electromagnetic heat exchange component;
if 0<a<1, calculating the nucleate boiling initiation temperature T by using an empirical correlation ONB Modal boiling onset temperature T MFB And maximum heat flux q of nucleate boiling CHF ；
If the current node temperature T w Greater than T MFB If the heat exchange coefficient of the low-temperature fuel at the node is calculated by adopting a Darr correlation formula, the calculation formula is h = h b +h tc +h dw +h EM Wherein h is b 、h tc 、h dw 、h EM Respectively a modal boiling component, a convective heat transfer component, a heat flow component and an electromagnetic heat transfer component;
if the current node temperature T w Less than T ONB If so, calculating the heat exchange coefficient of the low-temperature fuel at the node by adopting a Dittus-Boelter correlation formula;
if the current node temperature T w Between T MFB And T ONB Then comparing the current node heat flux density q w And maximum heat flux q of nucleate boiling CHF The size relationship of (1):
if q is w Less than q CHF If the heat exchange coefficient of the low-temperature fuel at the node is calculated by adopting a Chen correlation formula, the calculation formula is h = Eh fc +Sh p +h EM Wherein E is a forced convection heat transfer factor, S is a nucleate boiling heat transfer factor, h p Is the heat transfer component of nucleate boiling, h EM Is an electromagnetic heat exchange component;
if q is w Greater than q CHF Calculating the heat exchange coefficient of the low-temperature fuel at the node by adopting an interpolation method;
(4) According to the current distribution I obtained in the step (2) former (x) And (4) establishing a thermal analysis model according to the heat exchange coefficient distribution h (x) obtained in the step (3), and calculating to obtain axial temperature distribution T (x) of the superconducting energy pipeline;
the thermal analysis model adopts a one-dimensional thermal diffusion model, and the equation is as follows:
wherein A is c And P c Is the cross-sectional area and perimeter of the copper skeleton, T is the temperature value at the current node of the superconducting energy pipeline, ρ is the density of the copper skeleton, c p Is the specific heat capacity, k is the thermal conductivity, p R Is the resistivity, I former Is the copper skeleton current, h is the heat transfer coefficient at the current node, T f Is the operating temperature;
(5) Updating model parameters according to the axial temperature distribution T (x) obtained in the step (4) so as to realize the coupling effect between the temperature field and other physical fields, calculating the next time step, and repeating the steps (2) to (5) until the set time T is reached total When the temperature is higher than the preset temperature, the calculation is finished, and a curve T (T) of the change of the highest temperature along with the time is extracted;
wherein, t total Taking 100 seconds; the model parameters that need to be updated include: density rho and specific heat capacity c of copper skeleton p Thermal conductivity k, resistivity ρ R Density of cryogenic fuel ρ l Specific heat capacity c pl And coefficient of thermal conductivity k l ；
6) Judging whether the maximum temperature change T (T) of the superconducting energy pipeline obtained in the step (5) meets the overall quench judgment criterion of the energy pipeline, so as to predict whether the superconducting energy pipeline is quenched integrally;
the overall quench determination criterion is as follows: extracting a slope value at the last point of the curve T (T) of the change curve of the highest temperature along with time, and if the slope value is greater than zero, determining that the superconducting energy pipeline has integral quench; if the slope value is less than or equal to zero, the superconducting energy pipeline is considered not to be subjected to integral quench, and the pipeline system is safe to operate.
As shown in fig. 3, in example 1, the superconducting tape is Yttrium Barium Copper Oxide (YBCO) tape, the cryogenic fuel is liquefied natural gas, and the rated operating current of the cable is 3000A; in the graph of the change curve T (T) of the maximum temperature along with time, which is calculated by the method of the present invention, and the instantaneous working current of the cable acquired in the embodiment 1 is 20000A, the slope value at the last point is-0.012, which is less than zero, which indicates that the heat flow generated by the acquired fault current can be taken away by the low-temperature fuel, and the superconducting energy pipeline does not generate the integral quench.
As shown in fig. 4, in example 2, the superconducting tape is Yttrium Barium Copper Oxide (YBCO) tape, the cryogenic fuel is liquefied natural gas, and the rated operating current of the cable is 3000A; the instantaneous working current of the cable collected in the embodiment 2 is 26000A, and in a curve T (T) of the change of the highest temperature with time calculated by the method of the present invention, the slope value at the last point is 0.036 which is greater than zero, which indicates that the superconducting energy pipeline is subjected to the overall quench, and at this time, the working current should be immediately reduced to prevent the overall quench from occurring.
Claims (7)
1. A superconducting energy pipeline overall quench prediction method based on a multi-field coupling model is characterized by comprising the following steps:
(1) Measuring the critical current of the high-temperature superconducting tape at different temperatures in a laboratory, and fitting to obtain the critical current I c Data curve I as a function of temperature c (T); designing the structural parameters and the operation conditions of the superconducting energy pipeline; according to the operating temperature T f Calculating the physical properties of the low-temperature fuel and the cable material by calling NIST;
real-time acquisition of cable running current I in superconducting energy pipeline op ；
(2) According to the structural parameters and the critical current I of the superconducting energy pipeline in the step (1) c (T) establishing an electromagnetic field model, and calculating to obtain the current distribution I of the copper skeleton in the superconducting energy pipeline former (x)；
(3) Firstly, according to the current distribution I obtained in the step (2) former (x) Generating an electromagnetic heat exchange component h EM Thereby realizing the coupling effect between the electromagnetic field and the flow field; establishing a flow field analysis model according to the operation condition and the physical properties of the low-temperature fuel in the step (1), and calculating to obtain the axial distribution h (x) of the heat exchange coefficient of the low-temperature fuel;
(4) According to the current distribution I obtained in the step (2) former (x) And step (3)) Establishing a thermal analysis model according to the obtained heat exchange coefficient distribution h (x), and calculating to obtain axial temperature distribution T (x) of the superconducting energy pipeline;
(5) Updating model parameters according to the axial temperature distribution T (x) obtained in the step (4) so as to realize the coupling effect between the temperature field and other physical fields, calculating the next time step, and repeating the steps (2) to (5) until the set time T is reached total When the temperature is higher than the preset temperature, the calculation is finished, and a curve T (T) of the change of the highest temperature along with the time is extracted;
(6) And (4) judging whether the change curve T (T) of the highest temperature obtained in the step (5) along with time meets the judgment criterion of the integral quench of the superconducting energy pipeline or not so as to predict whether the superconducting energy pipeline is totally quenched or not.
2. The superconducting energy pipeline overall quench prediction method based on the multi-field coupling model according to claim 1, characterized in that: in the step (2), the electromagnetic field model is a conservation model based on kirchhoff's law, and the equation is as follows:
I op ＝I f +∑I i
wherein U is terminal voltage, I i Is a flowing current, R f Is the resistance of the copper skeleton, R 1 -R m Is the resistance of the respective superconducting tape, I f Is the current on the copper skeleton, I 1 -I m Is the current of each superconducting tape; I.C. A op Is the running current, L f Is self-inductance of the copper skeleton, L 1 -L m Is the self-inductance of the respective superconducting layer, M i,j Is the mutual inductance between superconducting tapes i and j;
the current-voltage characteristics of the superconducting tape follow an n-exponential power relation as follows:
wherein, E 0 Is the quench voltage, l is the strip length, I c Is the critical current, and n is the power exponent.
3. The superconducting energy pipeline overall quench prediction method based on the multi-field coupling model according to claim 1, characterized in that: the flow field analysis model in the step (3) is established by adopting an empirical correlation method, and specifically comprises the following steps:
judging whether the mass gas content a of the low-temperature fuel at the current calculation node is equal to 0 or 1, if so, calculating the heat exchange coefficient of the low-temperature fuel at the node by adopting a Dittus-Boelter correlation formula, wherein the calculation formula is h = h fc +h EM Where h is the heat transfer coefficient at the current node, h fc Is a single-phase forced convection heat transfer component, h EM Is the electromagnetic heat exchange component;
if 0<a<1, calculating the nucleate boiling initiation temperature T by using an empirical correlation ONB Modal boiling onset temperature T MFB And maximum heat flux density of nucleate boiling q CHF ；
If the current node temperature T w Greater than T MFB If the heat exchange coefficient of the low-temperature fuel at the node is calculated by adopting a Darr correlation formula, the calculation formula is h = h b +h tc +h dw +h EM Wherein h is b 、h tc 、h dw 、h EM Respectively a modal boiling component, a convective heat transfer component, a heat flow component and an electromagnetic heat transfer component;
if the current node temperature T w Less than T ONB If so, calculating the heat exchange coefficient of the low-temperature fuel at the node by adopting a Dittus-Boelter correlation formula;
if the current node temperature T w Between T MFB And T ONB Then comparing the current node heat flux density q w And maximum heat flux q of nucleate boiling CHF The size relationship of (1):
if q is w Less than q CHF If the heat exchange coefficient of the low-temperature fuel at the node is calculated by adopting a Chen correlation formula, the calculation formula is h = Eh fc +Sh p +h EM Wherein E is strongPreparing convection heat transfer factor, S is nucleate boiling heat transfer factor, h p Is the nucleate boiling heat transfer component, h EM Is an electromagnetic heat exchange component;
if q is w Greater than q CHF And calculating the heat exchange coefficient of the low-temperature fuel at the node by adopting an interpolation method.
4. The superconducting energy pipeline overall quench prediction method based on the multi-field coupling model according to claim 1, characterized in that: the thermal analysis model in the step (4) is a one-dimensional thermal diffusion model, and the equation is as follows:
wherein A is c And P c Is the cross-sectional area and perimeter of the copper skeleton, T is the temperature value at the current node of the superconducting energy pipeline, ρ is the density of the copper skeleton, c p Is the specific heat capacity, k is the thermal conductivity, p R Is the resistivity, I former Is the copper skeleton current, h is the heat transfer coefficient at the current node, T f Is the operating temperature.
5. The superconducting energy pipeline overall quench prediction method based on the multi-field coupling model according to claim 1, characterized in that: the model parameters to be updated in the step (5) comprise: density rho and specific heat capacity c of copper skeleton p Thermal conductivity k, resistivity ρ R Density of cryogenic fuel ρ l Specific heat capacity c pl And coefficient of thermal conductivity k l 。
6. The superconducting energy pipeline overall quench prediction method based on the multi-field coupling model according to claim 1, characterized in that: setting time t in the step (5) total The prediction result is more accurate when the time is more than 80 seconds.
7. The superconducting energy pipeline overall quench prediction method based on the multi-field coupling model according to claim 1, characterized in that: the overall quench determination criterion in the step (6) is as follows: extracting a slope value at the last point of the curve T (T) of the change curve of the highest temperature along with time, and if the slope value is greater than zero, determining that the superconducting energy pipeline has integral quench; if the slope value is less than or equal to zero, the superconducting energy pipeline is considered not to be subjected to integral quench, and the pipeline system is safe to operate.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010986813.8A CN112287515B (en) | 2020-09-18 | 2020-09-18 | Superconducting energy pipeline overall quench prediction method based on multi-field coupling model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010986813.8A CN112287515B (en) | 2020-09-18 | 2020-09-18 | Superconducting energy pipeline overall quench prediction method based on multi-field coupling model |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112287515A CN112287515A (en) | 2021-01-29 |
| CN112287515B true CN112287515B (en) | 2022-12-09 |
Family
ID=74420049
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010986813.8A Active CN112287515B (en) | 2020-09-18 | 2020-09-18 | Superconducting energy pipeline overall quench prediction method based on multi-field coupling model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112287515B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112906268B (en) * | 2021-02-07 | 2024-03-29 | 中国电力科学研究院有限公司 | Method and system for calculating quench resistivity of YBCO high-temperature superconducting unit |
| CN113064108A (en) * | 2021-03-15 | 2021-07-02 | 西安交通大学 | Device for measuring quenching and recovery characteristics of superconducting tape under forced convection cooling |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107895095A (en) * | 2017-12-15 | 2018-04-10 | 西安交通大学 | The coupling analytical method of cluster subchannel and critical heat flux density mechanism model |
| CN110688782A (en) * | 2019-09-12 | 2020-01-14 | 清华大学 | Long-distance superconducting energy pipeline simulation method considering axial and radial field distribution |
-
2020
- 2020-09-18 CN CN202010986813.8A patent/CN112287515B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107895095A (en) * | 2017-12-15 | 2018-04-10 | 西安交通大学 | The coupling analytical method of cluster subchannel and critical heat flux density mechanism model |
| CN110688782A (en) * | 2019-09-12 | 2020-01-14 | 清华大学 | Long-distance superconducting energy pipeline simulation method considering axial and radial field distribution |
Non-Patent Citations (1)
| Title |
|---|
| 不同载流运行下冷绝缘高温超导电缆的温度场数值分析;栗会峰等;《低温与超导》;20120924(第09期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112287515A (en) | 2021-01-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112287515B (en) | Superconducting energy pipeline overall quench prediction method based on multi-field coupling model | |
| Chen et al. | Energy-saving superconducting power delivery from renewable energy source to a 100-MW-class data center | |
| Thadela et al. | Computational investigation on thermohydraulic characteristics of High-Temperature Superconducting (HTS) power cables | |
| Nie et al. | An improved natural frequency based transmission line fault location method with full utilization of frequency spectrum information | |
| Lee et al. | Thermo-hydraulic analysis on long three-phase coaxial HTS power cable of several kilometers | |
| Liu et al. | Robust optimization in HTS cable based on design for six sigma | |
| Sato et al. | Temperature and Pressure Simulation of a High-Temperature Superconducting Cable Cooled by Subcooled $\mbox {LN} _2 $ With Fault Current | |
| CN113297750A (en) | Cable model selection method and system for offshore wind power plant | |
| Zhu et al. | Transient modeling of quench and recovery of LNG-HTS hybrid energy transmission system based on multi-field coupled analysis | |
| Miyagi et al. | Suitable cable structure of HTS triaxial cable cooled by counter flow cooling method for long-distance power transmission | |
| Yasui et al. | Temperature and pressure distribution simulations of 3-km-long high-temperature superconducting power cable system with fault current for 66-kV-class transmission lines | |
| Kalsia et al. | Statistical correlations for thermophysical properties of Supercritical Argon (SCAR) used in cooling of futuristic High Temperature Superconducting (HTS) cables | |
| Pi et al. | Numerical study of current distribution and stability of LTS/HTS hybrid superconductor | |
| Bang et al. | Critical current, critical temperature and magnetic field based EMTDC model component for HTS power cable | |
| Yasui et al. | Temperature and Pressure Simulation of a 1.5-km HTS Power Cable Cooled by Subcooled $\text {LN} _ {2} $ With a Fault Current | |
| Miyagi et al. | Thermal analysis of co-axial multi-layered BSCCO HTS power cable | |
| Li et al. | Comparison of liquid nitrogen flow resistance in corrugated pipe with smooth pipe for HTS cable | |
| Adachi et al. | Design of 22-kV 10-kA HTS triaxial superconducting bus | |
| Chen et al. | Thermal analysis of tri-axial HTS cable | |
| Savoldi et al. | Thermal-hydraulic models for the cooling of HTS power-transmission cables: status and needs | |
| Bian et al. | Temperature estimation of water cooling stator windings considering twisted structure and resistance distribution | |
| Liu et al. | Study of quench behavior of no-insulation REBCO pancake considering complex critical current density distribution | |
| Sitko et al. | 3-D numerical analysis of the 68 kA heat exchanger of the ITER HTS current leads | |
| Zhu et al. | Magneto-thermal coupling modeling for the stability analysis of a 110 kV/3 kA high temperature superconducting cable under fault current conditions | |
| Ge et al. | Hybrid current lead design of HTS SMES |
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 |