CN114068051A - Method for calculating temperature and flow of coolant of main pipeline of nuclear reactor based on ultrasonic array - Google Patents

Method for calculating temperature and flow of coolant of main pipeline of nuclear reactor based on ultrasonic array Download PDF

Info

Publication number
CN114068051A
CN114068051A CN202111251821.9A CN202111251821A CN114068051A CN 114068051 A CN114068051 A CN 114068051A CN 202111251821 A CN202111251821 A CN 202111251821A CN 114068051 A CN114068051 A CN 114068051A
Authority
CN
China
Prior art keywords
temperature
ultrasonic
coolant
flow
model
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
CN202111251821.9A
Other languages
Chinese (zh)
Other versions
CN114068051B (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.)
Sichuan University
Original Assignee
Sichuan 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 Sichuan University filed Critical Sichuan University
Priority to CN202111251821.9A priority Critical patent/CN114068051B/en
Publication of CN114068051A publication Critical patent/CN114068051A/en
Application granted granted Critical
Publication of CN114068051B publication Critical patent/CN114068051B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G21NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21CNUCLEAR REACTORS
    • G21C17/00Monitoring; Testing ; Maintaining
    • G21C17/017Inspection or maintenance of pipe-lines or tubes in nuclear installations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/27Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G21NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21CNUCLEAR REACTORS
    • G21C17/00Monitoring; Testing ; Maintaining
    • G21C17/02Devices or arrangements for monitoring coolant or moderator
    • G21C17/032Reactor-coolant flow measuring or monitoring
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00Energy generation of nuclear origin
    • Y02E30/30Nuclear fission reactors

Abstract

The invention relates to the field of ultrasonic measurement, in particular to a method for calculating the temperature and the flow of a coolant of a nuclear reactor main pipeline based on an ultrasonic array, which solves the problem of low measurement precision in the process of measuring the flow by ultrasonic in the prior art. The invention sets ultrasonic transducer group on the outer side of the pipeline, and connects with the remote server through the network module, one through CFD to build model to calculate the model field of temperature and flow speed under the input parameter; secondly, establishing a relation description model; measuring and calculating a three-reconstruction field; and fourthly, performing feedback setting on the topological structure of the ultrasonic transducer group through errors of the reconstruction field and the model field. According to the invention, the measurement precision of the transducer group is improved by eliminating the influence factors of ultrasonic waves in a non-aqueous medium; the temperature and the flow velocity under the input parameters are calculated through the CFD and the neural network model, variance calculation is carried out on the temperature and the flow velocity measured by the transducer group, and the measurement precision is further improved by adjusting the topological structure of the transducer group. The remote server is connected with the transducer group to realize real-time monitoring of measurement.

Description

Method for calculating temperature and flow of coolant of main pipeline of nuclear reactor based on ultrasonic array
Technical Field
The invention relates to the field of ultrasonic wave application, in particular to a method for calculating the temperature and the flow of a coolant of a main pipeline of a nuclear reactor based on an ultrasonic array.
Background
The nuclear reactor is a complex and efficient nuclear heat conversion device, coolant in a hot section of a main pipeline has the characteristics of high temperature, high pressure, high radiation and high flow rate, and meanwhile, because the enrichment degrees of all fuel assemblies of a reactor core are different, the temperature of the coolant flowing into the hot section of the main pipeline from different fuel channels is different, and the obvious temperature stratification phenomenon exists in the hot section of the main pipeline. The temperature and the flow of the coolant of the main pipeline of the nuclear reactor directly reflect the nuclear power and the heat conducting capacity of the reactor core, and are one of the very important thermal parameters in the safety protection and the operation control of the reactor. The accurate measurement and calculation of the temperature and the flow of the coolant of the main pipeline have important significance for the safety and the economical efficiency of the reactor.
The problems of insufficient reactor output, frequent control rod action and even unplanned shutdown caused by inaccurate and untimely measurement of the temperature of the coolant of the main reactor pipeline become important bottlenecks which seriously restrict further improving the safety and the economical efficiency of the reactor. Research data at home and abroad show that the thermal stratification temperature difference of the outlet of the pressurized water reactor core can reach 15 ℃, the thermal stratification temperature difference of the hot section of the main pipeline 1.5 meters away from the outlet of the core still reaches 10 ℃, and the stratification state continuously changes along with the operation of the reactor, so that the accurate and timely measurement of a representative coolant thermal physical field is very difficult.
The method for measuring the temperature of the coolant of the reactor main pipeline in the prior art is generally a thermal resistance method and an ultrasonic temperature measuring method: the thermal resistance point type temperature measurement technology has the following problems: 1. because the coolant in the hot section has a dynamic change layering phenomenon and is limited by the installation space of a main pipeline of the reactor, the problems of insufficient representativeness and large measurement error exist, the typical measurement error reaches 1.95 ℃ (equivalent to 5.5% of full power), and the operation capacity of the reactor is severely limited; 2. because the thermometer has the problem of larger thermal inertia, the response time of a temperature measurement channel exceeds 10s, and the timeliness of safety protection and control is seriously restricted; 3. the intrinsic safety of the reactor is reduced due to the need to open holes in the pressure boundary for installing the thermometer. Ultrasonic temperature measurement is used as a non-contact temperature measurement method, has the characteristics of timely response and high safety, can form surface temperature and even body temperature by measured linear temperature through a reconstruction algorithm, can greatly reduce calculation errors caused by temperature stratification compared with traditional point temperature measurement, and improves the measurement precision of average temperature. However, because the temperature, the pressure and the concentration of the reactor coolant have dynamic change characteristics, a composite relation between the ultrasonic propagation speed and multiple elements of the coolant is not established at present; the problem of large error caused by unreasonable division of sub-temperature zones also exists.
The traditional pipeline flow measurement method is widely applied to calculating the pipeline flow by measuring the pipeline flow speed based on ultrasonic waves, however, the propagation speed of the ultrasonic waves in a medium is influenced by various factors such as temperature, pressure, medium density and the like. When the influence factors change greatly, the calculation accuracy of the conventional flowmeter is greatly reduced, and the problems of low robustness, complex structure and the like of the ultrasonic flowmeter exist;
in the ultrasonic pipeline flow meter designed by the time difference principle, the measurement accuracy of the average flow velocity is determined by the measurement accuracy of the propagation time of the ultrasonic waves in the pipeline, and the measurement of the pipeline flow is influenced finally. However, in the actual measurement of the propagation time, the propagation time in the solid medium such as the pipe wall and the transducer is also included in the measurement of the propagation time, and the measurement accuracy is difficult to guarantee because the propagation speed of the ultrasonic wave in the pipe is interfered by factors such as temperature, pressure, liquid density and the like.
A new method for calculating temperature and flow rate that can solve the above problems is desired.
Disclosure of Invention
The invention provides a method for calculating the temperature and the flow of a coolant of a main pipeline of a nuclear reactor based on an ultrasonic array, which solves the problem of low energy calculation precision of the nuclear reactor in the prior art, and can realize the simultaneous online high-precision measurement of average temperature and flow.
The technical scheme of the invention is realized as follows: the method for calculating the temperature and the flow of the coolant of the main pipeline of the nuclear reactor based on the ultrasonic array comprises the following steps that an ultrasonic transducer group is arranged on the outer side of the pipeline and is connected with a remote server through a network module, and a model field is established: establishing a complete and continuous main pipeline coolant air-time domain temperature layered diffusion model and a flow velocity distribution model under input parameters through a CFD (computational fluid dynamics) construction model; secondly, establishing a relation description model of three elements of the ultrasonic propagation speed and the coolant; three-reconstruction field measurement and calculation: measuring the flight time of each effective sound wave propagation path; the network module transmits the signals of the transducer group to a remote server through the network module; reconstructing the real-time temperature distribution and the flow velocity distribution of the coolant in the pipeline; calculating the average temperature and flow; and fourthly, performing feedback setting on the topological structure of the ultrasonic transducer group through errors of the reconstruction field and the model field.
The third step comprises ultrasonic transmission time calibration, which comprises simulating an experiment platform and acquiring experiment data; the simulation experiment platform comprises a pipeline with controllable temperature and an ultrasonic transducer arranged on the outer side of the pipeline.
The ultrasonic transmission time calibration specifically comprises the following steps:
1 in the empty tube state and at a temperature T0Then, τ at this temperature is obtained by changing the propagation path of the ultrasonic wave0The result of (1);
Figure RE-GDA0003428298860000031
wherein L is1And L2The lengths of the two propagation paths, t1And t2The propagation time of the ultrasonic waves on the two propagation paths is respectively;
2 changing the temperature, and repeatedly calculating the temperature tau at different temperatures by experiments0The result of (1);
3, fitting the experimental calculation results to obtain a description model tau between the time delay and the temperature of the ultrasonic wave in the non-aqueous medium0=F(T)。
The second step is specifically as follows: (1) establishing a relational expression of three elements of ultrasonic sound velocity and coolant; (2) simulating an experiment platform and acquiring experiment data; (3) establishing a relation description model between the ultrasonic sound velocity and the temperature, pressure and concentration of the coolant by adopting a machine learning method; the simulation experiment platform also comprises a thermostatic bath with a water inlet, a water outlet and a concentration control device, and a metal pipeline for controlling pressure; the ultrasonic transducer is connected with the micro-control device.
The fourth step is specifically sub-area division: adjusting the ultrasonic transducer topology: adjusting the topological structure of the ultrasonic transducer array based on the error distribution of the reconstruction field and the model field; dynamic adjustment based on feedback mechanism: the radius R of the inner ring of the subarea by taking the reconstruction error as a target variableirAnd the radial angle theta is used as a control variable to determine the optimal inner ring radius RirAnd a radial angle theta.
The ultrasonic transducer array topological structure is in a broadcast type, and specifically, the ultrasonic transducers are uniformly arranged on the cross section of the pipeline to be measured.
The first step is specifically as follows: a, calculating the temperature and flow rate data of the coolant of the typical section of the pipeline under set parameters by calculating the fluid mechanics (CFD); b, training a CFD simulation result by using a neural network model: training the simulation result in the step a through a neural network model; c, prediction result: and c, predicting to obtain continuous and complete temperature and flow velocity distribution of the coolant in the pipeline by using the neural network model in the step b.
The fitting of the step 3 is a machine learning method combining more than two ultrasonic transmission characteristic individual learner learning and individual learner strategies.
The invention discloses a method for calculating the temperature and the flow of a coolant of a main pipeline of a nuclear reactor based on an ultrasonic array, which improves the measurement precision of a transducer group by eliminating the influence factors of ultrasonic waves in a non-aqueous medium; the temperature and the flow velocity under the input parameters are calculated through a CFD and a neural network model, variance calculation is carried out on the temperature and the flow velocity measured by the transducer group, a main pipeline hot section coolant space-time domain temperature layered diffusion model is reconstructed through ultrasonic measurement, a relation description model of three elements of an ultrasonic sound velocity coolant is obtained through a simulation experiment platform, and the measurement precision is further improved through adjusting the topological structure of the transducer group. The remote server is connected with the transducer group to realize real-time monitoring of measurement.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to these drawings without creative efforts.
FIG. 1: simulating an experiment platform;
FIG. 2: a pipeline subregion division schematic diagram;
FIG. 3: establishing a flow chart of a space-time domain temperature hierarchical diffusion description model of the coolant in the hot section of the main pipeline;
FIG. 4: an ultrasonic transducer array diagram;
wherein: 1. an ultrasonic transducer; 2. a metal pipe; 3. a micro-control device; 4. a thermostatic bath; 5. a pump; 6. A boric acid concentration control device; 7. a water inlet; 8. a water outlet; 9. a pressure controller.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The invention discloses a method for calculating the temperature and the flow of a coolant of a main pipeline of a nuclear reactor based on an ultrasonic array, wherein an ultrasonic transducer 1 group is arranged outside the pipeline, the transducer group is connected with a remote server through a network module, and a model field is established: establishing a complete and continuous main pipeline coolant air-time domain temperature layered diffusion model and a flow velocity distribution model under input parameters through a CFD (computational fluid dynamics) construction model; secondly, establishing a relation description model of three elements of the ultrasonic propagation speed and the coolant; three-reconstruction field measurement and calculation: measuring the flight time of each effective sound wave propagation path; the network module transmits the signals of the transducer group to a remote server through the network module; reconstructing the real-time temperature distribution and the flow velocity distribution of the coolant in the pipeline; calculating the average temperature and flow; specifically, three-dimensional flow velocity field reconstruction and three-dimensional temperature field reconstruction; fourthly, feedback setting of the topological structure of the ultrasonic transducer 1 group is carried out through errors of the reconstruction field and the model field.
The third step comprises ultrasonic transmission time calibration, which comprises simulating an experiment platform and acquiring experiment data; the simulation experiment platform comprises a pipeline with controllable temperature and an ultrasonic transducer 1 arranged on the outer side of the pipeline.
The ultrasonic transmission time calibration specifically comprises the following steps:
1 in the empty tube state and at a temperature T0Then, τ at this temperature is obtained by changing the propagation path of the ultrasonic wave0The result of (1);
Figure RE-GDA0003428298860000061
wherein L is1And L2The lengths of the two propagation paths, t1And t2The propagation time of the ultrasonic waves on the two propagation paths is respectively;
2 changing the temperature, and repeatedly calculating the temperature tau at different temperatures by experiments0The result of (1);
3, fitting the experimental calculation results to obtain a description model tau between the time delay and the temperature of the ultrasonic wave in the non-aqueous medium0=F(T)。
The second step is specifically as follows: (1) establishing a relational expression of three elements of ultrasonic sound velocity and coolant; (2) simulating an experiment platform and acquiring experiment data; (3) establishing a relation description model between the ultrasonic sound velocity and the temperature, pressure and concentration of the coolant by adopting a machine learning method; the simulation experiment platform also comprises a thermostatic bath 4 with a water inlet 7, a water outlet 8 and a concentration control device, and a metal pipeline 2 for controlling pressure; the ultrasonic transducer 1 is connected with a micro-control device 3.
The first step is specifically as follows: a, calculating the temperature and flow rate data of the coolant of the typical section of the pipeline under set parameters by calculating the fluid mechanics (CFD); b, training a CFD simulation result by using a neural network model: training the simulation result in the step a through a neural network model; c, prediction result: and c, predicting to obtain continuous and complete temperature and flow velocity distribution of the coolant in the pipeline by using the neural network model in the step b.
The fitting of the step 3 is a machine learning method combining more than two ultrasonic transmission characteristic individual learner learning and individual learner strategies.
The temperature data prediction using ELM was performed in the following manner: the cross-section data obtained by CFD analysis is used as a training sample with the format of { locmn,TmnWhere M is 1, …, M is the total number of temperature measurement cross sections calculated by the CFD software, N is 1, …, N is the total number of grids divided on the temperature measurement cross sections, and a grid diagram of the main pipe and the temperature measurement cross sections is shown in fig. 3. locmn={[(xmn,ymn),dm,tm]In which d ismDistance from m cross section to core outlet, tmFor the m-th cross-sectional coolant relative movement time, (x)mn,ymn) Is the coordinate of the center point of the nth region on the mth section, TmnCorresponding to the nth region of the m-th cross-sectionA temperature value. The ELM regression model for an activation function g (x) containing R hidden layer nodes is represented as:
Figure RE-GDA0003428298860000071
in the formula, ωr=[ωr1r2,L,ωrm]TRepresenting a weight between the input node and the r-th hidden layer node; beta is ar=[βr1r2,L,βrm]TRepresenting a weight between the output node and the r-th hidden layer node; brIs the ith hidden layer threshold; r is the number of hidden layer nodes.
To improve the model prediction capability, g (x) can be selected to train the samples in a zero-error approximation manner, i.e.
Figure RE-GDA0003428298860000072
Then there is
Figure RE-GDA0003428298860000073
The rewriting of the above formula into a matrix form can be expressed as
Gβ=T
Figure RE-GDA0003428298860000081
Where G is the hidden layer output matrix, the input weight ωrAnd hidden layer threshold brAll the values are randomly assigned, so that the output matrix G of the hidden layer is known, and the training is completed after the output weight matrix beta is obtained. Setting the number of the predicted sections as H, and predicting the coordinate information loc of the H predicted sectionshnInputting into the trained ELM model, and outputting the predicted temperature value ThnAnd the establishment of a space-time domain temperature hierarchical diffusion description model of the coolant in the hot section of the main pipeline is realized. The whole process flow of the method is shown in FIG. 3;
establishing a relation description model of three elements of ultrasonic propagation speed and coolant;
(1) and establishing a relation between the ultrasonic sound velocity and the coolant. According to the wave equation, the propagation velocity of ultrasonic waves in a liquid medium is related to the adiabatic coefficient of compression K and the density ρ of the liquid, which are expressed as follows:
Figure RE-GDA0003428298860000082
the adiabatic compression coefficient K and the density ρ are related to temperature and pressure. The coolant of the hot section of the reactor main pipeline is boric acid solution, and the temperature, the pressure and the boron concentration of the boric acid solution can influence the propagation speed of the ultrasonic wave. Therefore, the relationship between the propagation speed of the ultrasonic wave and the coolant is expressed as follows:
w=f(T,P,S) (X)
where w is the ultrasonic sound velocity, T is the coolant temperature, P is the coolant pressure, and S is the coolant boron concentration. And (3) giving an ultrasonic sound velocity expression under the multi-element composite influence in the coolant by referring to the ultrasonic transmission characteristics in seawater:
w=w0+wT+wP+wS+wTPS (XI)
in the formula w0Is a constant value, wTThe relationship between the speed of sound of the ultrasonic wave and the temperature of the coolant, wPThe relationship between the speed of sound of the ultrasonic wave and the pressure of the coolant, wSThe relationship between the ultrasonic sound velocity and the boron concentration of the coolant, wTPSIs the relationship between the ultrasonic sound velocity and the coolant temperature, pressure and boron concentration.
(2) Simulating an experiment platform and acquiring experiment data; and establishing a relation description model between the ultrasonic sound velocity and the temperature, pressure and concentration of the coolant by adopting a machine learning method. A real environment is simulated, a multi-element ultrasonic transmission characteristic experiment platform is built, the incidence relation between the ultrasonic sound velocity and the coolant is researched, and the schematic diagram of the experiment platform is shown in figure 1. The experimental platform takes boric acid solution (coolant) as an ultrasonic wave propagation medium, and the device can control and adjust the temperature, the pressure and the concentration of the coolant. In the experiment, the influence of each element of the coolant on the propagation speed of the ultrasonic wave is respectively researched by adopting a controlled variable method, and a large amount of experimental data is obtained.
Further, the simulation experiment platform comprises a thermostatic bath 4 with a water inlet 7, a water outlet 8 and a concentration control device, and a metal pipeline 2 for controlling pressure, wherein the metal pipeline 2 is also provided with an ultrasonic transducer 1; the ultrasonic transducer 1 is connected with a micro-control device 3.
Preferably, the method of machine learning employs more than two ultrasound transmission characteristics individual learner learning and individual learner combination strategies.
(1) Determining the topological structure and the effective sound wave path of the ultrasonic transducer 1; a radial axial array of several ultrasonic transducers is shown in fig. 4(a), taking 20 ultrasonic transducers 1 as an example, the transducers are uniformly mounted around the pipe, as indicated by the black dots in fig. 4 (a). These transducers are considered transceivers, which are controlled to transmit and detect ultrasonic signals at different times. The ultrasonic signal propagates from one transducer to another, resulting in an effective acoustic path as shown in fig. 4 (a). In theory, an ultrasonic path exists between every two transducers. However, since the ultrasonic path on the edge or surface does not contribute much to the reconstruction of the internal temperature, they are not necessary.
(2) Reconstructing a three-dimensional temperature field;
and after the flight time on each effective sound wave path is obtained, the three-dimensional cylindrical temperature field is inverted by using a radial basis temperature field reconstruction algorithm according to the mapping relation between the ultrasonic wave propagation speed and the three elements of the coolant. The several ultrasonic transducers 1 generate between each other I effective ultrasonic paths, which are divided into J sub-temperature zones, as shown in fig. 2. The obtained sonic wave flight time on each ultrasonic wave propagation path can be expressed as:
Figure RE-GDA0003428298860000101
where a (x, y, z) represents the reciprocal of the sound velocity of the ultrasonic wave, and when radial basis functions are used, a (x, y, z) is expressed as a linear combination of J radial basis functions:
Figure RE-GDA0003428298860000102
Figure RE-GDA0003428298860000103
in the formula: epsilonjIs the undetermined coefficient; (x)j,yj,zj) The coordinates of the central point of the jth sub-temperature zone; alpha is a radial basis function phijThe shape parameters of (x, y, z) are determined as appropriate during the course of a particular experiment.
Order to
Figure RE-GDA0003428298860000104
Then can obtain
Figure RE-GDA0003428298860000105
Rewriting formula (XVI) to matrix form
t=F·E (XVII)
In the formula: t ═ t1,t2,...,tI]T,F=(fkj)k=1,2,...I;j=1,2,...J,E=[ε12,...,εJ]T
Using singular value decomposition of the matrix F and Tikhonov regularization technique, the regularization solution of equation (XVII) is:
Figure RE-GDA0003428298860000106
in the formula: sigma is a nonzero singular value of the coefficient matrix F, and J is the total number of the nonzero singular values; u. ofj、vjThe left and right singular vectors of F are respectively; μ is a regularization parameter. When the ultrasonic wave is receivedAfter the array position is determined, the shape parameter alpha of the radial basis function is given, and the coefficient matrix F and the singular value thereof can be obtained. After the time-of-flight matrix t on each acoustic path is obtained, the parameter vector epsilon can be determined according to formula (XVIII). After obtaining the reciprocal of the sound velocity a (x, y, z) of the ultrasonic wave, the reciprocal is substituted into the acquired experimental data to establish the mapping relation between the propagation velocity of the ultrasonic wave and the three elements of the coolant:
T=f-1(w,P,S) (XIX)
the three-dimensional temperature distribution of the coolant can be reconstructed from the formula (XIX).
And calculating the error between the radial direction and the axial direction of the reconstructed field and the model field according to the established pipeline space-time domain temperature distribution model. And substituting the central coordinates of the divided sub-temperature areas to obtain the temperature value of the central coordinate of each sub-temperature area, solving the temperature value of the central coordinate point of each sub-temperature area of the reconstructed temperature field and calculating the maximum absolute error, the minimum absolute error, the average relative error and the root mean square error of the temperature value by using a model. The expression is as follows:
Emax=max|T(xi,yi,zi)-Tm(xi,yi,zi)| (XX)
Emin=min|T(xi,yi,zi)-Tm(xi,yi,zi)| (XXI)
Figure RE-GDA0003428298860000111
Figure RE-GDA0003428298860000112
in the formulas (XX) - (XXIII), n is the number of all central coordinate temperature values in the region to be measured; t isaIs the average temperature of the simulated field; t (x)i,yi,zi) For model field in coordinates (x)i,yi,zi) The temperature value of (a); t ism(xi,yi,zi) For reconstructing the field in coordinates (x)i,yi,zi) The temperature value of (2).
(3) Reconstructing a three-dimensional flow velocity field;
aiming at the complex flow state of different positions of a reactor pipeline, the existence of turbulent vortex in a pipeline flow field is considered, and compared with a two-dimensional pipeline section average flow velocity calculation scheme, the flow calculation method based on three-dimensional cylindrical flow field distribution is more accurate.
And (3) installing an ultrasonic detection array on a reactor pipeline, selecting a proper ultrasonic signal emission period, measuring the flight time of a sound wave path, calculating the flow velocity according to the difference of the countercurrent flight time and the downstream flight time after obtaining the flight time, and inverting the three-dimensional flow velocity field.
Time t of propagation of ultrasonic wave from downstream to upstreamuAnd upstream to downstream propagation time tdCan be respectively expressed as:
Figure RE-GDA0003428298860000121
the relationship between the ultrasonic flight time and the flow velocity can be obtained by the formula (1):
Figure RE-GDA0003428298860000122
therefore, after the flight time on each effective sound wave path is obtained, the three-dimensional flow velocity field is inverted by using a reconstruction algorithm based on a radial basis according to the relation between the ultrasonic flight time and the flow velocity.
Several ultrasound transducers 1 generate between each other I active ultrasound paths, dividing them into J sub-regions, as shown in fig. 2. The sub-regions being divided by an inner radius RirAnd the divided sub-regions are changed along with two key parameters of the radial angle theta, the radius of the inner ring is different, and the radial angle is different. The relationship between the flight time of the ultrasonic wave along a specific propagation path in the region to be measured and the flow velocity can be expressed as a line integral as follows:
Figure RE-GDA0003428298860000123
wherein a (x, y, z) represents the inverse of the flow rate, lkIs an ultrasonic wave propagation path. When radial basis functions are used, a (x, y, z) is expressed as a linear combination of J radial basis functions:
Figure RE-GDA0003428298860000124
Figure RE-GDA0003428298860000125
in the formula: epsilonjIs the undetermined coefficient; (x)j,yj,zj) The coordinates of the center point of the jth sub-region; alpha is a radial basis function phijThe shape parameters of (x, y, z) are determined as appropriate during the course of a particular experiment.
Order to
Figure RE-GDA0003428298860000131
Then can obtain
Figure RE-GDA0003428298860000132
Rewriting the formula (6) to a matrix form
Δt=F·E (8)
In the formula: Δ t ═ Δ t1,Δt2,...,ΔtI]T,F=(fkj)k=1,2,...I;j=1,2,...J,E=[ε12,...,εJ]T
Using the singular value decomposition of the matrix F and the Tikhonov regularization technique, the regularization solution of equation (27) is:
Figure RE-GDA0003428298860000133
in the formula: sigma is a non-zero singular value of the matrix F, and J is the total number of the non-zero singular values; u. ofj、vjThe left and right singular vectors of F are respectively; μ is a regularization parameter. After the position of the ultrasonic receiving and transmitting array is determined, the shape parameter alpha of the radial basis function is given, and then the matrix F and the singular value thereof can be obtained. The delta t can also be determined by actually measuring the sound wave flight time and the sound channel angle theta, the parameter vector epsilon can be determined according to the formula (28), and the parameter vector epsilon is substituted into the formula (23) to obtain the distribution function of the flow velocity reciprocal, so that the distribution function of the flow velocity can be obtained through the following formula, and the reconstruction of the flow velocity field is realized.
Figure RE-GDA0003428298860000134
Calculating the error between a reconstructed field and a model field according to the established pipeline space-time domain flow velocity distribution model, introducing the central coordinates in the divided sub-regions to obtain the flow velocity value of the central coordinate of each sub-region, and solving the flow velocity value of the central coordinate point of each sub-region of the reconstructed flow velocity field and the model to calculate the maximum absolute error, the minimum absolute error, the average relative error and the root mean square error of the flow velocity value. The expression is as follows:
Emax=max|v(xi,yi,zi)-vm(xi,yi,zi)| (11)
Emin=min|v(xi,yi,zi)-vm(xi,yi,zi)| (12)
Figure RE-GDA0003428298860000141
Figure RE-GDA0003428298860000142
in formulae (11) to (14):
n: the number of all central coordinate flow velocity values in the region to be measured
va: mean flow velocity of the simulated field
v(xi,yi,zi): model field in coordinates (x)i,yi,zi) Flow rate value of
vm(xi,yi,zi): flow velocity field in coordinate (x)i,yi,zi) Flow rate value of
Further, the third step further includes: a, calculating the average temperature in the measurement area according to the reconstructed temperature field; and b, calculating the average flow in the measurement area according to the reconstructed flow velocity field.
As shown in the modeling flowchart of the three-dimensional cylindrical pipe of the reactor in fig. 3, specifically: the reactor cylindrical pipeline is modeled by a space-time domain temperature and flow velocity distribution model, and the fluid can be modeled by single-phase turbulence because the fluid does not generate phase change in the pipeline. The following time-averaged monophasic continuity equations, momentum equations and energy equations may thus be established for describing the fluid flow process in the pipe.
Continuity equation:
Figure RE-GDA0003428298860000143
the momentum equation:
Figure RE-GDA0003428298860000144
wherein mueff=μ+μturb
Energy equation:
Figure RE-GDA0003428298860000145
wherein k iseff=k+Prtt/ρα)。
Since turbulence affects both mass and energy transfer, it is necessary to build a suitable turbulence model to quantify the turbulence viscosity. The method can be used for assisting analysis based on a standard k-epsilon turbulence model, the k-epsilon model solves the equation of turbulence momentum k and dissipation rate epsilon thereof, and the expression is as follows:
Figure RE-GDA0003428298860000151
Figure RE-GDA0003428298860000152
wherein the turbulent viscosity is defined as mut=ρCμ(k2/. epsilon.). The turbulence prandtl number is also introduced in the above equations (4) and (5) to quantify the turbulence energy dissipation coefficient of the turbulence viscosity.
As shown in a grid division diagram of the reactor cylindrical pipeline modeling in fig. 4, because the calculation for the above equation is converted into discrete equation calculation, the modeling area is subjected to grid division by using CFD software to realize discretization, and the temperature and flow velocity distribution of the pipeline flow measurement section is obtained.
In order to establish an overall space-time domain temperature distribution model and a flow velocity distribution model of a flow measuring section of a three-dimensional cylindrical pipeline of a reactor, data prediction needs to be carried out on the basis of existing data. And predicting the flow speed data by adopting a radial basis function RBF neural network, and predicting the temperature data by adopting an extreme learning machine ELM.
The flow rate data prediction by using the RBF neural network can adopt the following modes: and using the flow velocity data obtained by the CFD analysis as a training sample, wherein the format of the training sample is { X, V }, and X is [ X ═ X [ ]1,x2,x3,x4]T,x1Is the abscissa, x, of the position of the model mesh2Is the ordinate, x, of the position of the model mesh3Being the vertical coordinate, x, of the position of the model mesh4The relative movement time t of water in a reactor pipeline is taken as the time t; v is a one-dimensional output speed obtained by CFD simulation; RBF networkThe output is Y ═ Y](ii) a p is the number of hidden layer nodes (p)>4) N is the total number of input samples; the iteration is terminated with precision. The output of the RBF neural network can be expressed as:
Figure RE-GDA0003428298860000153
the structure of the radial basis neuron is shown in fig. 4, and the formula is as follows:
Figure RE-GDA0003428298860000161
the purpose of model training is to obtain three parameters W in the formulaj、CjAnd DjAnd obtaining a space-time domain flow velocity distribution model. The training of the RBF neural network is mainly divided into two steps, firstly, the parameter C between an input layer and a hidden layer needs to be determined through unsupervised learningjAnd DjThen, supervised learning is carried out to determine the weight W between the hidden layer and the output layerj
First, the connection weight W from the hidden layer to the output layer is [ W ]1,w2,…,wp]TAnd center parameter C of each node of the hidden layerj=[cj1,cj2,cj3,cj4]TWidth vector D of each node of hidden layerj=[dj1,dj2,dj3,dj4]T(j ═ 1,2, …, p) is initialized, with the initialization equations:
Figure RE-GDA0003428298860000162
in the formula (22), p is the number of hidden layer nodes (p > 4), min (V)k) Is the minimum value of all expected outputs in the output neurons in the training set, max (V)k) Is the maximum of all expected outputs in the output neurons in the training set.
Figure RE-GDA0003428298860000163
In formula (23), min (x)i) Is the minimum value, max (x), of all input information in the ith input feature of the training seti) Is the maximum value of all input information in the ith input feature of the training set.
Figure RE-GDA0003428298860000164
In the formula (24), dfThe width adjustment coefficient is generally less than 1, so that each hidden layer neuron can easily sense local information, and the RBF neural network is helped to improve the local response capability.
Then, a gradient descent method is used for carrying out self-adaptive adjustment on the central parameters, the width vectors and the adjustment weights, and the iterative calculation formulas are respectively as follows:
Figure RE-GDA0003428298860000165
Figure RE-GDA0003428298860000166
Figure RE-GDA0003428298860000171
wherein, the E is a learning factor, and the E is an RBF network evaluation function:
Figure RE-GDA0003428298860000172
and finally, judging whether the condition of quitting the iterative operation is met, namely the root mean square error RMS is less than or equal to epsilon, if so, finishing the training, and if not, continuing to perform the weight iterative calculation.
After training is finished, values of relevant parameters can be obtained, namely, the establishment of the space-time domain flow velocity distribution model of the reactor three-dimensional cylindrical pipeline is realized.
The fourth step is specifically sub-area division: adjusting the ultrasound transducer 1 topology: adjusting the topological structure of the ultrasonic transducer 1 array based on the error distribution of the reconstructed field and the model field; dynamic adjustment based on feedback mechanism: the radius R of the inner ring of the subarea by taking the reconstruction error as a target variableirAnd the radial angle theta is used as a control variable to determine the optimal inner ring radius RirAnd a radial angle theta.
The ultrasonic transducer 1 array topological structure is in a broadcast type, and specifically, the ultrasonic transducers 1 are uniformly arranged on the section of the pipeline to be measured.
Further, the fitting of step 3 is a machine learning method using more than two ultrasound transmission characteristics individual learner learning and individual learner strategy combination. The method divides experimental data for ultrasonic transmission time calibration into a training set, a verification set and a test set, wherein the training set is used for model training and optimizing parameters such as weight, bias and the like; the verification set is used for optimizing hyper-parameters, learning rate, regularization coefficients and the like; the test set is used to evaluate the generalization ability of the model. And determining the optimal data classification method, the individual learner model and the integration strategy through experiments to finally obtain the mapping relation description model.
The method adopts an average method for integration, and comprises a simple average method, a weighted average method and the like; when the data volume is very large, the data can be integrated by adopting a learning method, namely, a plurality of individual learners are combined by a new learner so as to improve the learning accuracy. And finally, the model with the minimum root mean square error obtained by using the test set is the optimal model.
The invention discloses a method for calculating the temperature and the flow of a coolant of a main pipeline of a nuclear reactor based on an ultrasonic array, which improves the measurement precision of a transducer group by eliminating the influence factors of ultrasonic waves in a non-aqueous medium; the temperature and the flow velocity under the input parameters are calculated through the CFD and the neural network model, variance calculation is carried out on the temperature and the flow velocity measured by the transducer group, and the measurement precision is further improved by adjusting the topological structure of the transducer group. The remote server is connected with the transducer group to realize real-time monitoring of measurement.
It is understood that various changes and modifications may be made by those skilled in the art without departing from the spirit and scope of the invention, and it is intended to cover in the appended claims all such changes and modifications.

Claims (8)

1. A method for calculating the temperature and the flow of a coolant of a main pipeline of a nuclear reactor based on an ultrasonic array is characterized in that an ultrasonic transducer group is arranged on the outer side of the pipeline and is connected with a remote server through a network module, and the method comprises the following steps:
a model field is established: establishing a complete and continuous main pipeline coolant air-time domain temperature layered diffusion model and a flow velocity distribution model under input parameters through a CFD (computational fluid dynamics) construction model;
secondly, establishing a relation description model of three elements of the ultrasonic propagation speed and the coolant;
three-reconstruction field measurement and calculation: measuring the flight time of each effective sound wave propagation path; the network module transmits the signals of the transducer group to a remote server through the network module; reconstructing the real-time temperature distribution and the flow velocity distribution of the coolant in the pipeline; calculating the average temperature and flow;
and fourthly, performing feedback setting on the topological structure of the ultrasonic transducer group through errors of the reconstruction field and the model field.
2. The method for calculating the temperature and the flow of the coolant of the main pipeline of the nuclear reactor based on the ultrasonic array as claimed in claim 1, wherein the method comprises the following steps: the third step comprises ultrasonic transmission time calibration, which comprises simulating an experiment platform and acquiring experiment data; the simulation experiment platform comprises a pipeline with controllable temperature and an ultrasonic transducer arranged on the outer side of the pipeline.
3. The method for calculating the temperature and the flow of the coolant of the main pipeline of the nuclear reactor based on the ultrasonic array as claimed in claim 2, wherein the method comprises the following steps: the ultrasonic transmission time calibration specifically comprises the following steps:
1 in the empty tube state and at a temperature T0Then, τ at this temperature is obtained by changing the propagation path of the ultrasonic wave0The result of (1);
Figure FDA0003320528950000011
wherein L is1And L2The lengths of the two propagation paths, t1And t2The propagation time of the ultrasonic waves on the two propagation paths is respectively;
2 changing the temperature, and repeatedly calculating the temperature tau at different temperatures by experiments0The result of (1);
3, fitting the experimental calculation results to obtain a description model tau between the time delay and the temperature of the ultrasonic wave in the non-aqueous medium0=F(T)。
4. The method for calculating the temperature and the flow of the coolant of the main pipeline of the nuclear reactor based on the ultrasonic array as claimed in claim 3, wherein the method comprises the following steps: the second step is specifically as follows: (1) establishing a relational expression of three elements of ultrasonic sound velocity and coolant; (2) simulating an experiment platform and acquiring experiment data; (3) establishing a relation description model between the ultrasonic sound velocity and the temperature, pressure and concentration of the coolant by adopting a machine learning method; the simulation experiment platform also comprises a thermostatic bath with a water inlet, a water outlet and a concentration control device, and a metal pipeline for controlling pressure; the ultrasonic transducer is connected with the micro-control device.
5. The method for calculating the temperature and the flow of the coolant of the main pipeline of the nuclear reactor based on the ultrasonic array as claimed in claim 3, wherein the method comprises the following steps: the fourth step is specifically sub-area division:
adjusting the ultrasonic transducer topology: adjusting the topological structure of the ultrasonic transducer array based on the error distribution of the reconstruction field and the model field;
dynamic adjustment based on feedback mechanism: the radius R of the inner ring of the subarea by taking the reconstruction error as a target variableirAnd the radial angle theta is used as a control variable to determine the optimal inner ring radius RirAnd a radial angle theta.
6. The method for calculating the temperature and the flow of the coolant of the main pipeline of the nuclear reactor based on the ultrasonic array as claimed in claim 5, wherein the method comprises the following steps: the ultrasonic transducer array topological structure is in a broadcast type, and specifically, the ultrasonic transducers are uniformly arranged on the cross section of the pipeline to be measured.
7. The method for calculating the temperature and the flow of the coolant of the main pipeline of the nuclear reactor based on the ultrasonic array as claimed in claim 6, wherein the method comprises the following steps: the first step is specifically as follows:
a, calculating the temperature and flow rate data of the coolant of the typical section of the pipeline under set parameters by calculating the fluid mechanics (CFD);
b, training a CFD simulation result by using a neural network model: training the simulation result in the step a through a neural network model;
c, prediction result: and c, predicting to obtain continuous and complete temperature and flow velocity distribution of the coolant in the pipeline by using the neural network model in the step b.
8. The method for calculating the temperature and the flow of the coolant of the main pipeline of the nuclear reactor based on the ultrasonic array as claimed in claim 3, wherein the method comprises the following steps: the fitting of the step 3 is a machine learning method combining more than two ultrasonic transmission characteristic individual learner learning and individual learner strategies.
CN202111251821.9A 2021-10-26 2021-10-26 Method for calculating temperature flow of nuclear reactor main pipeline coolant based on ultrasonic array Active CN114068051B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111251821.9A CN114068051B (en) 2021-10-26 2021-10-26 Method for calculating temperature flow of nuclear reactor main pipeline coolant based on ultrasonic array

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111251821.9A CN114068051B (en) 2021-10-26 2021-10-26 Method for calculating temperature flow of nuclear reactor main pipeline coolant based on ultrasonic array

Publications (2)

Publication Number Publication Date
CN114068051A true CN114068051A (en) 2022-02-18
CN114068051B CN114068051B (en) 2023-06-16

Family

ID=80235653

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111251821.9A Active CN114068051B (en) 2021-10-26 2021-10-26 Method for calculating temperature flow of nuclear reactor main pipeline coolant based on ultrasonic array

Country Status (1)

Country Link
CN (1) CN114068051B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113990539A (en) * 2021-10-26 2022-01-28 四川大学 Ultrasonic measurement reconstruction method for temperature distribution of coolant in hot section of nuclear reactor main pipeline
CN114877954A (en) * 2022-07-12 2022-08-09 杭州春来科技有限公司 Fixed pollution source measuring method and system
CN117057212A (en) * 2023-06-05 2023-11-14 四川大学 Acoustic reconstruction method for dynamic temperature field of nuclear power device
CN117473873A (en) * 2023-11-13 2024-01-30 上海交通大学 Nuclear thermal coupling realization method based on deep M & Mnet neural network
CN117473873B (en) * 2023-11-13 2024-04-26 上海交通大学 Nuclear thermal coupling realization method based on DeepM & Mnet neural network

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4318776A (en) * 1978-07-21 1982-03-09 Kraftwerk Union Aktiengesellschaft Measuring lance for boiling-water nuclear reactors
GB8529446D0 (en) * 1985-11-29 1986-01-08 Univ Aberdeen Divergent ultrasound arrays
US20100286956A1 (en) * 2009-05-08 2010-11-11 American Power Conversion Corporation System and method for predicting cooling performance of arrangements of equipment in a data center
EP2657658A1 (en) * 2012-04-25 2013-10-30 General Electric Company Ultrasonic flow measurement system
WO2014166847A1 (en) * 2013-04-10 2014-10-16 Areva Np Methods for simulating the flow of a fluid in a vessel of a nuclear reactor and for calculating the mechanical deformation of assemblies of a nuclear reactor core, and associated computer program products
KR101776567B1 (en) * 2017-01-11 2017-09-11 대한민국 Method for improving temperature stratification in order to efficient heating and cooling and energy saving in buildings and heating control system using thereof
US20170263341A1 (en) * 2016-03-10 2017-09-14 Westinghouse Electric Company Llc Reactor coolant system piping temperature distribution measurement system
CN207132929U (en) * 2017-08-11 2018-03-23 北京清环智慧水务科技有限公司 A kind of flow probe
CN108154672A (en) * 2018-02-07 2018-06-12 天津中德应用技术大学 System and method is configured in a kind of ultrasonic liquid flow meter Automatic parameter
CN208061413U (en) * 2018-02-07 2018-11-06 天津中德应用技术大学 A kind of ultrasonic liquid flow meter Automatic parameter configuration system
CN110268478A (en) * 2016-12-15 2019-09-20 斯特凡Tto有限公司 The method and process of subject's specificity computation model of decision support and diagnosis for cardiovascular disease are provided
CN110618197A (en) * 2019-02-01 2019-12-27 中国石油化工股份有限公司 Long-distance oil and gas pipeline full-defect combined ultrasonic array nondestructive scanning method
CN113001070A (en) * 2015-03-26 2021-06-22 克里凯文斯管线国际有限公司 Pipeline cooling system
CN113375828A (en) * 2021-05-10 2021-09-10 电子科技大学 Temperature field reconstruction method based on ultrasonic waves

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4318776A (en) * 1978-07-21 1982-03-09 Kraftwerk Union Aktiengesellschaft Measuring lance for boiling-water nuclear reactors
GB8529446D0 (en) * 1985-11-29 1986-01-08 Univ Aberdeen Divergent ultrasound arrays
US20100286956A1 (en) * 2009-05-08 2010-11-11 American Power Conversion Corporation System and method for predicting cooling performance of arrangements of equipment in a data center
EP2657658A1 (en) * 2012-04-25 2013-10-30 General Electric Company Ultrasonic flow measurement system
WO2014166847A1 (en) * 2013-04-10 2014-10-16 Areva Np Methods for simulating the flow of a fluid in a vessel of a nuclear reactor and for calculating the mechanical deformation of assemblies of a nuclear reactor core, and associated computer program products
CN113001070A (en) * 2015-03-26 2021-06-22 克里凯文斯管线国际有限公司 Pipeline cooling system
US20170263341A1 (en) * 2016-03-10 2017-09-14 Westinghouse Electric Company Llc Reactor coolant system piping temperature distribution measurement system
CN110268478A (en) * 2016-12-15 2019-09-20 斯特凡Tto有限公司 The method and process of subject's specificity computation model of decision support and diagnosis for cardiovascular disease are provided
KR101776567B1 (en) * 2017-01-11 2017-09-11 대한민국 Method for improving temperature stratification in order to efficient heating and cooling and energy saving in buildings and heating control system using thereof
CN207132929U (en) * 2017-08-11 2018-03-23 北京清环智慧水务科技有限公司 A kind of flow probe
CN108154672A (en) * 2018-02-07 2018-06-12 天津中德应用技术大学 System and method is configured in a kind of ultrasonic liquid flow meter Automatic parameter
CN208061413U (en) * 2018-02-07 2018-11-06 天津中德应用技术大学 A kind of ultrasonic liquid flow meter Automatic parameter configuration system
CN110618197A (en) * 2019-02-01 2019-12-27 中国石油化工股份有限公司 Long-distance oil and gas pipeline full-defect combined ultrasonic array nondestructive scanning method
CN113375828A (en) * 2021-05-10 2021-09-10 电子科技大学 Temperature field reconstruction method based on ultrasonic waves

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
XINZHI ZHOU: "Ultrasonic Thermometry Algorithm Based on Inverse Quadratic Function", 《IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL》, vol. 68, no. 5, XP011852017, DOI: 10.1109/TUFFC.2020.3036116 *
崔岢: "基于TDC7200的超声测温系统设计", 《仪表技术与传感器》, vol. 439, no. 8 *
石文泽: "高温铸锻件内部温度场电磁超声重构方法研究", 《仪器仪表学报》, vol. 41, no. 8 *
秦洁: "反应堆冷却剂系统不可隔离滞流分支管热疲劳瞬态的CFD研究", 《核动力工程》, vol. 32, no. 1 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113990539A (en) * 2021-10-26 2022-01-28 四川大学 Ultrasonic measurement reconstruction method for temperature distribution of coolant in hot section of nuclear reactor main pipeline
CN113990539B (en) * 2021-10-26 2023-05-23 四川大学 Ultrasonic measurement reconstruction method for temperature distribution of coolant in hot section of main pipeline of nuclear reactor
CN114877954A (en) * 2022-07-12 2022-08-09 杭州春来科技有限公司 Fixed pollution source measuring method and system
CN114877954B (en) * 2022-07-12 2022-09-23 杭州春来科技有限公司 Fixed pollution source measuring method and system
CN117057212A (en) * 2023-06-05 2023-11-14 四川大学 Acoustic reconstruction method for dynamic temperature field of nuclear power device
CN117473873A (en) * 2023-11-13 2024-01-30 上海交通大学 Nuclear thermal coupling realization method based on deep M & Mnet neural network
CN117473873B (en) * 2023-11-13 2024-04-26 上海交通大学 Nuclear thermal coupling realization method based on DeepM & Mnet neural network

Also Published As

Publication number Publication date
CN114068051B (en) 2023-06-16

Similar Documents

Publication Publication Date Title
CN114068051B (en) Method for calculating temperature flow of nuclear reactor main pipeline coolant based on ultrasonic array
Chamoli ANN and RSM approach for modeling and optimization of designing parameters for a V down perforated baffle roughened rectangular channel
CN113990539B (en) Ultrasonic measurement reconstruction method for temperature distribution of coolant in hot section of main pipeline of nuclear reactor
Heng et al. Transient thermal prediction methodology for parabolic trough solar collector tube using artificial neural network
Peng et al. Neural networks analysis of thermal characteristics on plate-fin heat exchangers with limited experimental data
Beigzadeh et al. Prediction of thermal and fluid flow characteristics in helically coiled tubes using ANFIS and GA based correlations
Wang et al. Data-driven CFD modeling of turbulent flows through complex structures
CN105444923A (en) Mechanical temperature instrument error prediction method based on genetic-algorithm optimized least square support vector machine
CN113947045A (en) Method for calculating pipeline flow based on setting of topological structure of neural network and energy converter
CN113551744A (en) Ultrasonic flowmeter performance online monitoring method and system
Sridharan Applications of artificial intelligence techniques in heat exchanger systems
Uluer et al. Using the artificial neural network model for modeling the performance of the counter flow vortex tube
CN103617563B (en) A kind of water supply network theoretical based on geo-statistic spatial analysis is without monitoring node pressure determination statement
CN103729534A (en) Particulate fouling experimental device, prediction method and prediction system for arc-tube heat exchanger
Sridharan et al. Application of generalized regression neural network in predicting the thermal performance of solar flat plate collector systems
Xi et al. Cooling performance analysis and structural parameter optimization of X-type truss array channel based on neural networks and genetic algorithm
Jing et al. Investigation of open channel flow with unsubmerged rigid vegetation by the lattice Boltzmann method
Esfandyari et al. Optimization of ultrasonic-excited double-pipe heat exchanger with machine learning and PSO
CN108021755A (en) The definite method of the convection transfer rate of porous media
Al-Salaymeh Optimization of hot-wire thermal flow sensor based on a neural net model
Liu et al. Surrogate-based optimization and experiment validation of a fan-shaped film cooling hole with a large lateral space
CN114216582A (en) High-precision rapid temperature field reconstruction method, system, equipment and terminal
CN113887142A (en) Average flow velocity calculation method based on channel cross section
CN116595845A (en) Structure optimization method of spiral fin tube type heat exchanger
Thummar et al. Machine learning for vortex flowmeter design

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