CN112287515A - 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 PDF

Info

Publication number
CN112287515A
CN112287515A CN202010986813.8A CN202010986813A CN112287515A CN 112287515 A CN112287515 A CN 112287515A CN 202010986813 A CN202010986813 A CN 202010986813A CN 112287515 A CN112287515 A CN 112287515A
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.)
Granted
Application number
CN202010986813.8A
Other languages
Chinese (zh)
Other versions
CN112287515B (en
Inventor
厉彦忠
祝乘风
谭宏博
史佳敏
聂杨
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN202010986813.8A priority Critical patent/CN112287515B/en
Publication of CN112287515A publication Critical patent/CN112287515A/en
Application granted granted Critical
Publication of CN112287515B publication Critical patent/CN112287515B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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 IcData curve I as a function of temperaturec(T); then establishing an electromagnetic field model and calculating the current distribution I of the copper frameworkformer(x) (ii) a 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 Iformer(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 axial temperature distribution T (x)Repeating the steps until the set time t is reachedtotalWhen the calculation is finished, extracting a change curve T (t) of the highest temperature along with the time; finally, judging whether the maximum temperature change T (t) of the superconducting energy pipeline accords with the overall quench judgment criterion of the energy pipeline; the invention improves the running safety of the superconducting energy pipeline.

Description

Superconducting energy pipeline overall quench prediction method based on multi-field coupling model
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. generally adopts a liquid form to achieve low transmission loss and high transmission capacity; the two are jointly transported, a low-temperature operation environment is provided for the superconducting cable by utilizing the liquid fuel, and the superconducting cable is assembled into a superconducting energy pipeline, so that the integrated transportation of electric power and fossil energy is realized, and the energy transportation efficiency can be greatly improved.
The superconducting energy pipeline is widely researched and discussed in recent years, the concept of an LNG/power mixed transportation 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 transportation pipeline system is designed, and the energy loss of the combined transportation system is greatly reduced compared with that of the two independent transportation 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 quench phenomenon relates to a complex coupling relation 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, and the operating state of the cable also affects the flowing state of the cryogenic fuel.
At present, various simulation models for multiple physical fields have been proposed in academia, for example, chinese patent CN201910997413.4, etc., and the main method is to separately calculate the electromagnetic field, flow field and temperature field distribution step by using commercial software, so that the strong coupling relationship between multiple 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 quench of a superconducting cable mainly include a voltage method, a current method, a temperature rise method, a pressure method, a flow velocity method, and an ultrasonic method, but these methods can only detect a local quench after the quench occurs, and cannot predict a total 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 IcData curve I as a function of temperaturec(T); designing structural parameters and operation conditions of the superconducting energy pipeline; according to the operating temperature TfCalculating 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 pipelineop；
(2) According to the structural parameters and the critical current I of the superconducting energy pipeline in the step (1)c(T) modeling an electromagnetic fieldAnd calculating to obtain the current distribution I of the copper framework in the superconducting energy pipelineformer(x)；
(3) Firstly, according to the current distribution I obtained in the step (2)former(x) Generating an electromagnetic heat exchange component hEMThereby 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 property 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) Establishing a thermal analysis model, 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 reachedtotalWhen the calculation is finished, extracting a change curve T (t) of the highest temperature along with the time;
(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:
Figure BDA0002689543690000021
Iop＝If+∑Ii
wherein U is terminal voltage, IiIs a flowing current, RfIs the resistance of the copper skeleton, R1-RmIs the resistance of the respective superconducting tapes, IfIs the current on the copper skeleton, I1-ImIs the current of each superconducting tape; i isopIs the running current, LfIs self-inductance of the copper skeleton, L1-LmIs the self-inductance of the respective superconducting layer, Mi,jIs a superMutual inductance between conduction bands i and j;
the current-voltage characteristics of the superconducting strip material follow an n-exponential power relation as follows:
Figure BDA0002689543690000022
wherein E is0Is the quench voltage, l is the strip length, IcIs 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 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 ═ hfc+hEMWhere h is the heat transfer coefficient at the current node, hfcIs a single-phase forced convection heat transfer component, hEMIs the electromagnetic heat exchange component;
if 0<a<1, calculating the nucleate boiling initiation temperature T by using an empirical correlationONBModal boiling onset temperature TMFBAnd maximum heat flux q of nucleate boilingCHF；
If the current node temperature TwGreater than TMFBIf 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 ═ hb+htc+hdw+hEMWherein h isb、htc、hdw、hEMRespectively a modal boiling component, a convective heat transfer component, a heat flow component and an electromagnetic heat transfer component;
if the current node temperature TwLess than TONBIf 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 TwBetween TMFBAnd TONBThen comparing the current node heat flux density qwAnd maximum heat flux q of nucleate boilingCHFThe size relationship of (1):
if q iswLess than qCHFIf the heat exchange coefficient of the low-temperature fuel at the node is calculated by using Chen correlation formula, the calculation formula is h ═ Ehfc+Shp+hEMWherein E is a forced convection heat transfer factor, S is a nucleate boiling heat transfer factor, hpIs the nucleate boiling heat transfer component, hEMIs an electromagnetic heat exchange component;
if q iswGreater than qCHFAnd 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:
Figure BDA0002689543690000031
wherein A iscAnd PcIs 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, cpIs the specific heat capacity, k is the thermal conductivity, pRIs the resistivity, IformerIs the copper skeleton current, h is the heat transfer coefficient at the current node, TfIs the operating temperature.
The model parameters to be updated in the step (5) comprise: density rho and specific heat capacity c of copper skeletonpThermal conductivity k, resistivity ρRDensity of cryogenic fuel ρlSpecific heat capacity cplAnd coefficient of thermal conductivity kl。
Setting time t in the step (5)totalThe 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 a curve T (t) of the change 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 of the variation T (t) of the maximum temperature with time in example 1 of the present invention.
FIG. 4 is a graph of the variation T (t) of the maximum temperature with time 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 IcData curve I as a function of temperaturec(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 cablee3000A, operating temperature T of superconducting energy pipelinef80-90K, superconducting cable winding coefficient ai＝1,aj1, superconducting cable winding pitch Lp362 mm; according to the operating temperature TfAnd (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 rho of copper skeleton 9000kg/m3Specific heat capacity cp380J/(kg. K), thermal conductivityK400W/(m.K), resistivity ρR＝2x10-9Omega.m; density of cryogenic fuel ρl＝486.55kg/m3Specific heat capacity cpl3077J/(kg. K), thermal conductivity Kl＝0.234W/(m·K)；
Real-time acquisition of running current I of cable in superconducting energy pipelineop；
As shown in fig. 2, the multi-field coupling model consists of three sub-models: 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 pipelineformer(x)；
The electromagnetic field model adopts a conservation model based on kirchhoff's law, and the equation is as follows:
Figure BDA0002689543690000041
Iop＝If+∑Ii
wherein U is terminal voltage, IiIs a flowing current, RfIs the resistance of the copper skeleton, R1-RmIs the resistance of the respective superconducting tapes, IfIs the current on the copper skeleton, I1-ImIs the current of each superconducting tape; i isopIs the running current, LfIs self-inductance of the copper skeleton, L1-LmIs the self-inductance of the respective superconducting layer, Mi,jIs 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 E00.0001V/m, the expression is:
Figure BDA0002689543690000051
wherein E is0Is the quench voltage, l is the strip length, IcIs the critical current;
the critical current of the superconducting tape is obtained by the step (1), and the expression is as follows:
Figure BDA0002689543690000052
(3) firstly, according to the current distribution I obtained in the step (2)former(x) Generating an electromagnetic heat exchange component hEMThereby 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 property 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;
the flow field analysis model 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 ═ hfc+hEMWhere h is the heat transfer coefficient at the current node, hfcIs a single-phase forced convection heat transfer component, hEMIs the electromagnetic heat exchange component;
if 0<a<1, calculating the nucleate boiling initiation temperature T by using an empirical correlationONBModal boiling onset temperature TMFBAnd maximum heat flux q of nucleate boilingCHF；
If the current node temperature TwGreater than TMFBIf 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 ═ hb+htc+hdw+hEMWherein h isb、htc、hdw、hEMRespectively a modal boiling component, a convective heat transfer component, a heat flow component and an electromagnetic heat transfer component;
if the current node temperature TwLess than TONBIf 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 TwBetween TMFBAnd TONBBetweenThen comparing the current node heat flux density qwAnd maximum heat flux q of nucleate boilingCHFThe size relationship of (1):
if q iswLess than qCHFIf the heat exchange coefficient of the low-temperature fuel at the node is calculated by using Chen correlation formula, the calculation formula is h ═ Ehfc+Shp+hEMWherein E is a forced convection heat transfer factor, S is a nucleate boiling heat transfer factor, hpIs the nucleate boiling heat transfer component, hEMIs an electromagnetic heat exchange component;
if q iswGreater than qCHFCalculating 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) Establishing a thermal analysis model, 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:
Figure BDA0002689543690000053
wherein A iscAnd PcIs 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, cpIs the specific heat capacity, k is the thermal conductivity, pRIs the resistivity, IformerIs the copper skeleton current, h is the heat transfer coefficient at the current node, TfIs the working 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 reachedtotalWhen the calculation is finished, extracting a change curve T (t) of the highest temperature along with the time;
wherein, ttotalTaking 100 seconds; the model parameters that need to be updated include: density rho and specific heat capacity c of copper skeletonpThermal conductivity k, resistivity ρRLower, lowerDensity of warm fuel ρlSpecific heat capacity cplAnd coefficient of thermal conductivity kl；
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 a curve T (t) of the change 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 with time, which is calculated by the method of the present invention and is obtained in the example 1 when the instantaneous working current of the cable is 20000A, the slope value at the last point is-0.012 and is less than zero, which indicates that the heat flow generated by the collected fault current can be taken away by the low-temperature fuel, and the superconducting energy pipeline can 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 IcData curve I as a function of temperaturec(T); designing structural parameters and operation conditions of the superconducting energy pipeline; according to the operating temperature TfCalculating 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 pipelineop；
(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 pipelineformer(x)；
(3) Firstly, according to the current distribution I obtained in the step (2)former(x) Generating an electromagnetic heat exchange component hEMThereby 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 property 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) Establishing a thermal analysis model, 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 reachedtotalWhen the calculation is finished, extracting a change curve T (t) of the highest temperature along with the time;
(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.
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:
Figure FDA0002689543680000011
Iop＝If+∑Ii
wherein U is terminal voltage, IiIs a flowing current, RfIs the resistance of the copper skeleton, R1-RmIs the resistance of the respective superconducting tapes, IfIs the current on the copper skeleton, I1-ImIs the current of each superconducting tape; i isopIs the running current, LfIs self-inductance of the copper skeleton, L1-LmIs the self-inductance of the respective superconducting layer, Mi,jIs 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:
Figure FDA0002689543680000012
wherein E is0Is the quench voltage, l is the strip length, IcIs 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 ═ hfc+hEMWhere h is the heat transfer coefficient at the current node, hfcIs a single-phase forced convection heat transfer component, hEMIs the electromagnetic heat exchange component;
if 0<a<1, calculating the nucleate boiling initiation temperature T by using an empirical correlationONBModal boiling onset temperature TMFBAnd maximum heat flux q of nucleate boilingCHF；
If the current node temperature TwGreater than TMFBThen, thenThe heat exchange coefficient of the low-temperature fuel at the node is calculated by adopting a Darr correlation formula, wherein the calculation formula is h ═ hb+htc+hdw+hEMWherein h isb、htc、hdw、hEMRespectively a modal boiling component, a convective heat transfer component, a heat flow component and an electromagnetic heat transfer component;
if the current node temperature TwLess than TONBIf 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 TwBetween TMFBAnd TONBThen comparing the current node heat flux density qwAnd maximum heat flux q of nucleate boilingCHFThe size relationship of (1):
if q iswLess than qCHFIf the heat exchange coefficient of the low-temperature fuel at the node is calculated by using Chen correlation formula, the calculation formula is h ═ Ehfc+Shp+hEMWherein E is a forced convection heat transfer factor, S is a nucleate boiling heat transfer factor, hpIs the nucleate boiling heat transfer component, hEMIs an electromagnetic heat exchange component;
if q iswGreater than qCHFAnd 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:
Figure FDA0002689543680000021
wherein A iscAnd PcIs 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, cpIs the specific heat capacity, k is the thermal conductivity, pRIs the resistivity, IformerIs the copper skeleton current, h is the heat transfer coefficient at the current node, TfIs to workAnd (3) 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 skeletonpThermal conductivity k, resistivity ρRDensity of cryogenic fuel ρlSpecific heat capacity cplAnd coefficient of thermal conductivity kl。
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)totalThe 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 a curve T (t) of the change 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.
CN202010986813.8A 2020-09-18 2020-09-18 Superconducting energy pipeline overall quench prediction method based on multi-field coupling model Active CN112287515B (en)

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 true CN112287515A (en) 2021-01-29
CN112287515B 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)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112906268A (en) * 2021-02-07 2021-06-04 中国电力科学研究院有限公司 Calculation method and system for 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)

* Cited by examiner, † Cited by third party
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

Patent Citations (2)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
栗会峰等: "不同载流运行下冷绝缘高温超导电缆的温度场数值分析", 《低温与超导》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112906268A (en) * 2021-02-07 2021-06-04 中国电力科学研究院有限公司 Calculation method and system for quench resistivity of YBCO high-temperature superconducting unit
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

Also Published As

Publication number Publication date
CN112287515B (en) 2022-12-09

Similar Documents

Publication Publication Date Title
CN112287515B (en) Superconducting energy pipeline overall quench prediction method based on multi-field coupling model
CN107451381B (en) Transformer thermal stability calculation method based on magnetic-thermal-current coupling
CN103616578A (en) Method for determining multi-state reliability of high temperature superconducting cable system in power system
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
CN107944079B (en) Method for rapidly acquiring temperature rise of pipe-arranging laying three-phase cable group based on transfer matrix
Ballarino HTS current leads: Performance overview in different operating modes
CN116662721B (en) COMSOL-based heat pipe cooling reactor nuclear thermoelectric internal coupling numerical calculation method
Zhu et al. Transient modeling of quench and recovery of LNG-HTS hybrid energy transmission system based on multi-field coupled analysis
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
Pi et al. Numerical study of current distribution and stability of LTS/HTS hybrid superconductor
Kalsia et al. Statistical correlations for thermophysical properties of Supercritical Argon (SCAR) used in cooling of futuristic High Temperature Superconducting (HTS) cables
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
Zhu et al. Multi-field coupled effect of thermal disturbance on quench and recovery characteristic along the hybrid energy pipe
Li et al. Comparison of liquid nitrogen flow resistance in corrugated pipe with smooth pipe for HTS cable
Savoldi et al. Thermal-hydraulic models for the cooling of HTS power-transmission cables: status and needs
Zanino et al. Pressure drop analysis in the CS insert coil
Bian et al. Temperature estimation of water cooling stator windings considering twisted structure and resistance distribution
Kan et al. Numerically thermal stability analysis on a geometrically symmetrical strand fabricated by 2G wires in liquid helium temperature
Ge et al. Hybrid current lead design of HTS SMES
Yokoo et al. Safety verification of a 275-kV HTS cable system under short-circuit current accidents
Takeda et al. Temperature and pressure simulations of 66 kV 40 m HTS cable system in short-circuit current accidents compared with experimental results
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

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