CN103970992B - The visual analysis method of flow concussion after a kind of lead bismuth fluid gas injection - Google Patents

Abstract

The present invention discloses the visual analysis method of flow concussion after a kind of lead bismuth fluid gas injection, the method includes: inject lead bismuth alloy in (1) Double-tube cyclic device, respectively to lead bismuth alloy injects different amounts of helium or argon, then calculate the thermal physical property parameter under different temperatures according to measurement apparatus；(2) Double-tube cyclic device is carried out stress and strain model, set the Structure matrix of CFD computation model, and set convergence control criterion；(3) according to the Structure matrix set, solving equation is listed, according to the free convection characteristic of lead bismuth alloy, the selected isoparametric computation model of the coefficient of heat transfer, by the way of iterative cycles, the calculating that it is carried out load data and transmission；(4) utilize CFD software that data message is converted into image information, it is shown intuitively.

Description

The visual analysis method of flow concussion after a kind of lead bismuth fluid gas injection

Technical field

The invention belongs to nuclear energy power generation and reactor safety analysis field, be specifically related to neck Territory is that lead bismuth fluid injects the hottest physical property and flow concussion after gas with various Analysis.

Background technology

This special material of lead bismuth alloy, its physicochemical properties aspect exists many The difference of common heap-type coolant, carried out the thermal-hydraulic of lead bismuth alloy relatively in the past Analyzing and researching the most necessary, the physicochemical properties aspect of lead bismuth alloy has many special Different part, it is high that it has boiling point, low-melting feature, the most also have good in Son learns the advantages such as characteristic, Flouride-resistani acid phesphatase and good heat-transfer.The one of lead bismuth Natural Circulation Form be inject gas, gas be filled with primarily to improve the ability of Natural Circulation, It is referred to as SAPMAC method.As important coolant, its hot-working hydraulic characteristic is for institute Experimental stand or the reactor core built have particularly significant, and this can ensure that its heat Normal derivation.And on the one hand the flow concussion of coolant can affect the reliability of heat transfer, On the other hand may make bringing-up section that heat exhaustion occurs, thus affect its service life, Unfavorable to system.

The country that countries in the world carry out studying to this now is not the most a lot, the most individually Several countries are strong builds complete test-bed and studies.Domestic starting ratio Later, the research of these mechanisms is mainly carried out by the experimental stand having built up The experimentation of lead bismuth Natural Circulation, by experimentation natural circulation flow, heat transfer Etc. the feature of aspect.

One-dimensional lead bismuth Natural Circulation is tested by P.agostini in 2002 etc.； G.Coccoluto in 2011 etc. have studied lead bismuth to NAICE loop injection argon and naturally follow Circulation；A.V.Besnosov in 2011 etc. are to the lead under high temperature and the speed of lead bismuth fluid Degree field and temperature field have carried out experimentation.Current Domestic ADS has been in step Section, the research to lead bismuth heat exchange is less, North China Electric Power University's nuclear heat work safe level standard Fluid interchange after lead bismuth gas injection is studied by the Zou Wenchong changing institute；Fructus Perillae Prestige have studied the thermal stratification of lead bismuth；Chinese Academy of Sciences's Hefei nuclear energy safe practice is ground Study carefully the corrosion problems of institute's primary study lead bismuth alloy material.At present, relate to both at home and abroad And the research of the lead bismuth natural circulation flow concussion of thermal-hydraulic safety is considerably less, the most still Do not form a kind of analysis method.After the most at different temperatures to lead bismuth fluid gas injection, The visual analysis method studying the concussion of its natural circulation flow is the most necessary.

In existing periodical file, in electric power journal volume 28 the 4th interim disclosure in 2013 The periodical of article numbered 1005-6548 (2013) 04-0336-05, this journal title It is referred to as liquid lead bismuth alloy hot physical property program development research, this periodical discloses physical property Calculation procedure, but change procedural knowledge and simply describe the calculating process of hot physical property, There is no the calculating of the hot physical property of in-depth study, this periodical does not has systematic to the conjunction of lead bismuth yet The flow reforming phenomena produced after gold gas injection is studied.

Summary of the invention

In order to solve the problems referred to above, present inventor has performed and study with keen determination, found that: After being completed lead bismuth fluid gas injection by measurement apparatus and Double-tube cyclic device by use The calculating of flow concussion and analysis so that the calculating of flow concussion after lead bismuth fluid gas injection Convenient with the process of analysis, procedure.

It is an object of the invention to provide following aspect:

The visual analyzing side of flow concussion after first aspect, a kind of lead bismuth fluid gas injection Method, it is characterised in that: the method is the completeest by measurement apparatus and Double-tube cyclic device Becoming, described measurement apparatus includes:

Initialization module, the temperature value of its reception Double-tube cyclic device, and by temperature value It is delivered to Numerical Simulation Module,

Numerical Simulation Module, the temperature value of its reception initialization module transmission, and according to Temperature value calculates lead bismuth alloy thermal physical property parameter at such a temperature, and by hot physical property Parameter is delivered to output module,

Output module, it receives the thermal physical property parameter of Numerical Simulation Module transmission and by it Output CFD computation model the method includes:

(1) Double-tube cyclic device is carried out stress and strain model, set CFD computation model Structure matrix, and set convergence control criterion；

(2) Double-tube cyclic device injects lead bismuth alloy, note in lead bismuth alloy respectively Enter different amounts of helium or argon, then calculate under different temperatures according to measurement apparatus Thermal physical property parameter；

(3) according to the Structure matrix set, solving equation is listed, according to lead bismuth alloy Free convection characteristic, the selected isoparametric computation model of the coefficient of heat transfer, pass through iteration The mode of circulation, the calculating that it is carried out load data and transmission；

(4) utilize CFD software that data message is converted into image information, it is directly perceived Show.

Flow after second aspect, a kind of lead bismuth fluid gas injection as described in above-mentioned first aspect The visual analysis method of amount concussion, it is characterised in that: the heat solved in step (1) Physical parameter includes density, specific heat capacity, thermal conductivity, dynamic viscosity, kinematic viscosity.

Flow after the third aspect, a kind of lead bismuth fluid gas injection as described in above-mentioned first aspect The visual analysis method of amount concussion, it is characterised in that: Double-tube cyclic device includes:

Precooling zone 1, increase in pipeline 2 is passed through with double in the pipeline one end of the thermal medium that wherein circulates Pipe bringing-up section 3 top is connected, and the body other end flows in cooling section 5 by declining pipeline 4 The body top of body heat medium is connected；

Cooling section 5, the body top of the thermal medium that wherein circulates is connected with declining pipeline 4, Its other end is connected with preheating section 7 by the first pipeline 6；

Preheating section 7, its one end is connected with cooling section 5 by the first pipeline 6, its another Hold by being connected bottom the second two-tube bringing-up section of pipeline 83,

Two-tube bringing-up section 3, it includes electric heating tube 3c, bottom electric heating tube 3c and the second line of a couplet The top of case 3b is connected, and electric heating tube 3c top is connected bottom upper header 3a；The first line of a couplet The top of case 3a is connected away from one end of precooling zone 1 with increase in pipeline 2, lower collecting box 3b's Bottom is connected away from one end of preheating section 7 with the second pipeline 8；

Wherein, the second pipeline 8 arranging gas injection port 9, lead bismuth alloy flows successively through preheating It flow back into preheating section after section, two-tube bringing-up section, precooling zone, cooling section.

Flow after fourth aspect, a kind of lead bismuth fluid gas injection as described in the above-mentioned third aspect The visual analysis method of amount concussion, it is characterised in that:

Outside increase in pipeline 2, decline pipeline the 4, first pipeline 6 and the first pipeline 8 all use Footpath 32mm, the steel pipe of internal diameter 25mm.

Flow after 5th aspect, a kind of lead bismuth fluid gas injection as described in the above-mentioned third aspect The visual analysis method of amount concussion, it is characterised in that:

The heating power of preheating section is 6KW, and or

The heating power of two-tube bringing-up section is 7KW-15KW.

There is advantages that

The first, by the present invention in that to enter with by measurement apparatus physical property hot to lead bismuth alloy Row calculates, and improves the efficiency of calculating, saves the calculating time；

The second, the present invention uses Double-tube cyclic device so that the research of hot physical property is more pasted Near actual, meet the situation of reality, not only include that data are calculated, also include true Calculating；

3rd, image is converted to by new for data number by calling the built-in relative program of CFD Signal, directly perceived, concrete shows on screen by result, after having calculated, protects Deposit and carry out post processing, in conjunction with hydromechanical knowledge, it being completed oscillation of flows Analyze so that flow concussion after lead bismuth fluid gas injection can be displayed intuitively.

Accompanying drawing explanation

After Fig. 1 illustrates a kind of according to the preferred embodiment of the present invention lead bismuth fluid gas injection The flow chart of the visual analysis method of flow concussion；

After Fig. 2 illustrates a kind of according to the preferred embodiment of the present invention lead bismuth fluid gas injection The structural representation of Double-tube cyclic device in the visual analysis method of flow concussion；

After Fig. 3 illustrates a kind of according to the preferred embodiment of the present invention lead bismuth fluid gas injection The workflow diagram of measurement apparatus in the visual analysis method of flow concussion.

Drawing reference numeral illustrates:

1-precooling zone

2-increase in pipeline

The two-tube bringing-up section of 3-

3a-upper header

3b-lower collecting box

3c-electric heating tube

4-declines pipeline

5-cooling section

6-the first pipeline

7-preheating section

8-the second pipeline

9-gas injection port

Detailed description of the invention

Below by the present invention is described in detail, the features and advantages of the invention will Become more apparent from, clearly along with these explanations.

The most special word " exemplary " means " as example, embodiment or say Bright property ".The here as any embodiment illustrated by " exemplary " should not necessarily be construed as excellent In or be better than other embodiments.Although the various aspects of embodiment shown in the drawings, But unless otherwise indicated, it is not necessary to accompanying drawing drawn to scale.

According in a preferred embodiment of the invention, as illustrated in fig. 1-3, it is provided that The visual analysis method of flow concussion after a kind of lead bismuth fluid gas injection, the method is passed through Measurement apparatus and Double-tube cyclic device complete the meter of flow concussion after lead bismuth fluid gas injection Calculate and analyze.

In one preferred embodiment, as shown in Figure 3, measurement apparatus is by defeated The temperature value entering the outlet of two-tube hot arc calculates the hot physical property ginseng that this temperature value is corresponding Number, and is delivered to CFD software by this thermal physical property parameter, thermal physical property parameter include density, Specific heat capacity, thermal conductivity, dynamic viscosity, kinematic viscosity.

Wherein, the Numerical Simulation Module in measurement apparatus is according to lead bismuth alloy temperature value meter Calculate the thermal physical property parameter of lead bismuth alloy, wherein,

(1) for the calculating of lead bismuth alloy density according to equation below

ρ=11096-1.3236 T

The range of formula is 400K-1273K.

Wherein, ρ represents that density, unit are kg/m³；T represents that temperature, unit are K.

(2) for the calculating of lead bismuth alloy specific heat capacity according to equation below

C=159-2.72 × 10^-2·T+7.12×10^-6·T²

The range of formula is 400K-1273K.

Wherein, C is specific heat capacity, and unit is J/ (kg K)；T represents that temperature, unit are K。

(3) for the calculating of lead bismuth alloy thermal conductivity according to equation below

λ=3.61+1.517 × 10^-2·T-1.741×10^-6·T²

The range of formula is 400K-1273K.

Wherein, λ is thermal conductivity, and unit is W/ (m K)；T represents temperature, unit It is K.

(4) for the calculating of lead bismuth alloy dynamic viscosity according to equation below

$\mathrm{\η}=4.94\×{10}^{-4}\·\exp \left(\frac{754.1}{T}\right)$

The range of formula is 400K-1273K.

Wherein, η is dynamic viscosity, and unit is Pa；T represents that temperature, unit are K.

(5) computing formula for kinematic viscosity is as follows:

$\mathrm{\ν}=\frac{\mathrm{\η}}{\mathrm{\ρ}}$

Wherein, η is dynamic viscosity, and unit is Pa；ρ is density, and unit is kg/m³； ν is kinematic viscosity, and unit is m²/s。

(6) expression formula for coefficient of heat transfer h is as follows:

$h=\mathrm{Nu}\frac{\mathrm{\λ}}{d}$

In formula, h is the coefficient of heat transfer, and unit is W/ (m²·K)；λ is thermal conductivity, single Position is W/ (m K)；D is pipeline equivalent diameter, and unit is m；Nu is nusselt number.

(7) for the heat transfer formula of lead bismuth alloy.

Nu=A+0.014Pe^0.8

$A=\left\{\begin{array}{cc}4.5& \mathrm{Pe}\≤1000\\ 5.4-9{10}^{-4}\mathrm{Pe}& 1000\≤\mathrm{Pe}\≤2000\\ 3.6& \mathrm{Pe}\≥2000\end{array}\right.$

Wherein, Nu is nusselt number；Pr is Prandtl number.

Draw Pr_tRelational expression as follows:

${\mathrm{Pr}}_t=\left\{\begin{array}{cc}4.12& \mathrm{Pe}\≤1000\\ \frac{0.01\mathrm{Pe}}{{[{0.018}^{0.8}-(7-A)]}^{1.25}}& 1000\≤\mathrm{Pe}\≤6000\end{array}\right.$

Wherein, Pr_tIt it is turbulent prandtl number；Pe is Berkeley number.

The most above-mentioned (1)-(7) are according to Liu Mengying. lead bismuth alloy natural circulation heat transfer research [D]. and Beijing: North China Electric Power University, 2013.

In one preferred embodiment, as in figure 2 it is shown, Double-tube cyclic device is used It is simulated, by Double-tube cyclic device in completing lead bismuth alloy is circulated state Lead bismuth alloy flow concussion feelings at different temperatures can be simulated more really Condition, described Double-tube cyclic device includes preheating section, two-tube bringing-up section, precooling zone, cold But section；

Wherein, preheating section is for preheating the lead bismuth alloy in device, in order to Making lead bismuth alloy reach the minimum temperature that Natural Circulation rises, described preheating section there is no Particular restriction, can with heating furnace or heating net,

Wherein, two-tube bringing-up section is for heating lead bismuth alloy so that lead bismuth closes Gold completes the ascent stage of Natural Circulation, and two-tube bringing-up section there is no particular restriction, permissible For being made up of the thermocouple arranged in bringing-up section and bringing-up section outer wall,

Wherein, described precooling zone and cooling section are the completeest for carrying out cooling down for lead bismuth alloy Becoming the Natural Circulation of lead bismuth alloy, described precooling zone and cooling section there is no particular restriction, It can be conventional cooler.

In one preferred embodiment, the present invention Double-tube cyclic device built with Lead bismuth alloy, measures lead bismuth alloy gas injection in Double-tube cyclic device by CFD software After temperature value；

Temperature value according to receiving utilizes self-editing program to calculate the hot thing of lead bismuth alloy Property parameter；

Utilize CFD software to divide the grid of two-tube model, set corresponding perimeter strip Part, inputs physical parameter on request, then calculates.Last CFD software will stream Amount concussion situation is shown.

There is advantages that

The first, by the present invention in that to enter with by measurement apparatus physical property hot to lead bismuth alloy Row calculates, and improves the efficiency of calculating, saves the calculating time；

The second, the present invention uses Double-tube cyclic device so that the research of hot physical property is more pasted Near actual, meet the situation of reality, be no longer data calculations；

3rd, image is converted to by new for data number by calling the built-in relative program of CFD Signal, directly perceived, concrete shows on screen by result, after having calculated, protects Deposit and carry out post processing, in conjunction with hydromechanical knowledge, it being completed oscillation of flows Analyze so that flow concussion after lead bismuth fluid gas injection can be displayed intuitively.

Describe the present invention above in association with specific embodiments.But, these Embodiment is merely illustrative, and protection scope of the present invention is not intended that any by it Limit.It will be appreciated by those skilled in the art that without departing from or deviation scope In the case of, technical solutions and their implementation methods of the present invention have multiple modification, improvement Or equivalent, these all should fall within the scope of protection of the present invention.

Claims (4)

1. a visual analysis method for flow concussion after the gas injection of lead bismuth fluid, it is special Levy and be: the method is completed jointly by measurement apparatus and Double-tube cyclic device, described Measurement apparatus includes:

Initialization module, the temperature value of its reception Double-tube cyclic device, and by temperature value It is delivered to Numerical Simulation Module,

Numerical Simulation Module, the temperature value of its reception initialization module transmission, and according to Temperature value calculates lead bismuth alloy thermal physical property parameter at such a temperature, and by hot physical property Parameter is delivered to output module, and

Output module, it receives the thermal physical property parameter of Numerical Simulation Module transmission and by it Export in CFD computation model；

The method includes:

(1) Double-tube cyclic device is carried out stress and strain model, set CFD computation model Structure matrix, and set convergence control criterion；

(2) Double-tube cyclic device injects lead bismuth alloy, note in lead bismuth alloy respectively Enter different amounts of helium or argon, then calculate under different temperatures according to measurement apparatus Thermal physical property parameter；

(3) according to the Structure matrix set, solving equation is listed, according to lead bismuth alloy Free convection characteristic, the computation model of selected coefficient of heat transfer parameter, followed by iteration The mode of ring, the calculating that it is carried out load data and transmission；

(4) utilize CFD software that data message is converted into image information, it is directly perceived Show,

Double-tube cyclic device includes:

Precooling zone (1), increase in pipeline (2) is passed through in the pipeline one end of the thermal medium that wherein circulates Being connected with two-tube bringing-up section (3) top, the body other end is by declining pipeline (4) It is connected with the body top of fluid thermal medium in cooling section (5)；

Cooling section (5), the body top of the thermal medium that wherein circulates and decline pipeline (4) Being connected, its other end is connected with preheating section (7) by the first pipeline (6)；

Preheating section (7), the first pipeline (6) and cooling section (5) phase are passed through in its one end Connection, its other end is connected by the second pipeline (8) two-tube bringing-up section (3) bottom,

Two-tube bringing-up section (3), it includes electric heating tube (3c), electric heating tube (3c) Bottom is connected with the top of lower collecting box (3b), electric heating tube (3c) top and the first line of a couplet Case (3a) bottom is connected；The top of upper header (3a) and increase in pipeline (2) away from One end of precooling zone (1) is connected, the bottom of lower collecting box (3b) and the second pipeline (8) One end away from preheating section (7) is connected；

Wherein, the second pipeline (8) arranging gas injection port (9), lead bismuth alloy flows successively It flow back into preheating section after crossing preheating section, two-tube bringing-up section, precooling zone, cooling section.

After a kind of lead bismuth fluid gas injection the most as claimed in claim 1, flow concussion can Depending on fractional analysis method, it is characterised in that: the thermal physical property parameter bag solved in step (1) Include density, specific heat capacity, thermal conductivity, dynamic viscosity, kinematic viscosity.

After a kind of lead bismuth fluid gas injection the most as claimed in claim 1, flow concussion can Depending on fractional analysis method, it is characterised in that:

Increase in pipeline (2), decline pipeline (4), the first pipeline (6) and the first pipeline (8) external diameter 32mm, the steel pipe of internal diameter 25mm are all used.

After a kind of lead bismuth fluid gas injection the most as claimed in claim 1, flow concussion can Depending on fractional analysis method, it is characterised in that:

The heating power of preheating section is 6KW, and or

The heating power of two-tube bringing-up section is 7KW-15KW.