CN113486482A - Method for calculating speed and temperature boundary layer of sweepforward spiral tube bundle of liquid lead bismuth - Google Patents
Method for calculating speed and temperature boundary layer of sweepforward spiral tube bundle of liquid lead bismuth Download PDFInfo
- Publication number
- CN113486482A CN113486482A CN202110784158.2A CN202110784158A CN113486482A CN 113486482 A CN113486482 A CN 113486482A CN 202110784158 A CN202110784158 A CN 202110784158A CN 113486482 A CN113486482 A CN 113486482A
- Authority
- CN
- China
- Prior art keywords
- liquid lead
- lead bismuth
- temperature
- boundary layer
- spiral tube
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- 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
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- 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/08—Fluids
-
- 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
-
- 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/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Computer Networks & Wireless Communication (AREA)
- Computational Mathematics (AREA)
- Complex Calculations (AREA)
Abstract
The invention discloses a method for calculating a liquid lead bismuth sweepforward spiral tube bundle speed and temperature boundary layer, which comprises the following steps: 1. establishing a simplified geometric model of a fluid domain of the sweepforward spiral tube bundle of the liquid lead and bismuth; 2. establishing a liquid lead bismuth sweepforward spiral tube beam body area grid model; 3. performing numerical simulation calculation in the flowing heat exchange process of the liquid lead and bismuth; 4. extracting the fluid section speed and temperature field distribution information of the liquid lead-bismuth swept-out spiral tube bundle, and calculating the speed and temperature gradient distribution condition; 5. determining the position of a speed boundary layer and extracting the thickness of the boundary layer; 6. the temperature boundary layer position is determined and the temperature boundary layer thickness is extracted. The method can be applied to the calculation of the speed and temperature boundary layer of the sweepforward spiral tube bundle of the liquid lead bismuth to obtain the distribution of the speed and temperature boundary layer of the inner wall surface of the fluid domain.
Description
Technical Field
The invention belongs to the technical field of methods and particularly relates to a method for calculating a liquid lead bismuth sweepforward spiral tube bundle speed temperature boundary layer.
Background
Compared with common fluids such as water and air, the flowing heat exchange of the liquid metal is greatly different. The reason for the difference is mainly the stronger thermal conductivity of the liquid metal in comparison, which also makes the prandtl number of the liquid metal in the order of 0.01-0.1 (conventional fluid ≈ 1). The molecular prandtl number is commonly used to characterize the ratio of fluid flow to diffusion capacity, i.e., the ratio of the fluid flow boundary layer to the temperature boundary layer. For a conventional fluid, the exploration of the change law of the boundary layer of the liquid metal can help to deeply analyze the change law of the flow heat exchange.
In a common turbulence model, a reynolds-like method is commonly used to simulate the energy exchange in turbulence. The accuracy of the turbulent plandt model will greatly affect the accuracy of the numerical simulation. In the currently common RANS model, the default set turbulence Plantt number is typically 0.85, which has proven its high accuracy in a number of studies for conventional fluids. For liquid metals, this arrangement is no longer suitable, and therefore, many scholars have proposed a corresponding model of the planckian number of turbulence for liquid metal flow.
Therefore, when experimental research is carried out on the flowing heat exchange of the liquid metal, particularly when fluid domains with complex structures such as an sweepforward spiral tube bundle and the like are aimed at, the temperature position of the main flow fluid needs to be determined. The analysis method provided by the invention is helpful for exploring the speed and temperature boundary layer of the liquid metal sweepback spiral tube bundle by using a numerical simulation method, finely describes the real physical process in the speed and temperature boundary layer of the liquid metal, and has important significance for exploring the flow heat transfer rule of the liquid metal sweepback spiral tube bundle.
Disclosure of Invention
The invention provides a method for calculating a liquid lead bismuth sweepback spiral tube bundle speed and temperature boundary layer.
In order to achieve the purpose, the invention adopts the following technical scheme:
a liquid lead bismuth sweepback spiral tube bundle speed temperature boundary layer calculation method aims at the liquid lead bismuth sweepback spiral tube bundle flow and realizes the calculation of a speed temperature boundary layer by a fluid dynamics calculation program;
the method comprises the following steps:
step 1: establishing software SOLIDWORKS by using a geometric model, establishing a simplified model of the calculation domain of the liquid lead bismuth swept-out spiral tube bundle to obtain a fluid domain geometric model of the liquid lead bismuth swept-out spiral tube bundle, selecting a representative calculation unit in consideration of model symmetry, and simplifying the calculation domain of the liquid lead bismuth swept-out spiral tube bundle into a circumferential 1/12 cylindrical sleeve of the liquid lead bismuth swept-out spiral tube bundle;
step 2: performing mesh division on the liquid lead bismuth swept-out spiral tube beam body area on the basis of the liquid lead bismuth swept-out spiral tube bundle fluid area geometric model obtained in the step 1 by using mesh division software ANSYS-ICEM to obtain a liquid lead bismuth swept-out spiral tube beam body area mesh model, wherein the liquid lead bismuth swept-out spiral tube beam body area mesh model comprises boundary layer area meshes and integral meshes;
and step 3: performing numerical simulation calculation on the flowing heat exchange process of the liquid lead bismuth in the liquid lead bismuth sweepforward spiral tube beam body area grid model obtained in the step 2, and specifically comprising the following steps:
step 3-1: setting the boundary of the upper surface of the beam body area grid model of the sweepforward spiral tube of the liquid lead bismuth as an inlet boundary of the liquid lead bismuth, and setting the speed, the temperature and the pressure of the liquid lead bismuth at the inlet position according to actual conditions;
step 3-2: setting the boundary of the lower surface of the beam body area grid model of the sweepforward spiral tube of the liquid lead bismuth as the outlet boundary of the liquid lead bismuth, and setting the pressure of the liquid lead bismuth at the outlet position according to actual conditions;
step 3-3: setting the boundary of the outer surface of the tube bundle of the liquid lead bismuth sweepback spiral tube beam body area grid model as a constant wall temperature boundary, simulating the process of cooling the liquid lead bismuth by the tube bundle, setting the wall surfaces of the inner sleeve and the outer sleeve as heat insulation boundaries, and not exchanging heat between the wall surfaces of the inner sleeve and the outer sleeve and the liquid lead bismuth;
step 3-4: the physical parameters of the liquid lead bismuth obtained by calculation comprise the density of the liquid lead bismuth, the dynamic viscosity of the liquid lead bismuth alloy, the specific heat capacity of the liquid lead bismuth alloy and the dynamic thermal conductivity of the liquid lead bismuth alloy, and the density of the liquid lead bismuth is calculated as follows:
ρLBE＝11096-1.3236T (1)
t-temperature of liquid lead bismuth, K
The calculation formula of the specific heat capacity of the liquid lead-bismuth alloy is as follows:
cp,LBE＝159-2.72×10-2T+7.12×10-6T2 (2)
the dynamic viscosity of the liquid lead-bismuth alloy is calculated as follows:
the calculation formula of the dynamic thermal conductivity of the liquid lead-bismuth alloy is as follows:
λLBE＝3.61+1.517×10-2T-1.741×10-6T2 (4)
step 3-5: selecting a turbulent flow Plantt digital model suitable for liquid lead bismuth flow heat exchange simulation, loading the turbulent flow Plantt digital model into computational fluid dynamics software, and correcting the turbulent flow Plantt digital model in the computational fluid dynamics software FLUENT, wherein the selected model is as follows:
Prt＝4.12Pe≤1000 (6)
pe-liquid lead bismuth flow peclet number;
Prt-turbulent prandtl number;
step 3-6: solving the mass, momentum, energy and component transport equation of the liquid lead bismuth to obtain a liquid lead bismuth sweepforward spiral tube bundle calculation domain velocity field and a liquid lead bismuth temperature field;
and 4, step 4: establishing a cross section superposed with the spiral tube bundle by utilizing post-processing software TECPLOT from the calculation result of the speed field and the temperature field of the liquid lead bismuth swept-out spiral tube bundle calculation domain obtained in the step 3, extracting the fluid cross section speed and the temperature field distribution information of the liquid lead bismuth swept-out spiral tube bundle, calculating the fluid cross section speed and the temperature gradient distribution condition of the liquid lead bismuth swept-out spiral tube bundle, and solving the equation of the speed gradient and the temperature gradient;
u, V, W-x, y, z direction component velocity, m/s;
t-liquid lead bismuth temperature, K;
a-velocity gradient, s-1；
b-temperature gradient, K/m;
and 5: calculating the distribution information of the domain speed field according to the liquid lead bismuth sweepback spiral tube bundle obtained in the step 3, calculating the gradient distribution condition of the domain speed according to the liquid lead bismuth sweepback spiral tube bundle, determining the position of a speed boundary layer, and extracting the thickness of the speed boundary layer, wherein the specific steps are as follows:
step 5-1: calculating the distribution information of the domain velocity field according to the liquid lead bismuth sweepback spiral tube bundle obtained in the step 3, determining a velocity boundary layer region, extracting the position information of the velocity boundary layer, and determining the region as the velocity boundary layer region when the local flow velocity of the liquid lead bismuth is less than 99% of the main flow velocity according to the judgment basis:
U＝0.99U∞ (10)
U∞-the main flow velocity of liquid lead bismuth, m/s;
step 5-2: according to the velocity gradient distribution condition of the fluid section of the liquid lead bismuth sweepback spiral tube bundle obtained in the step 4, when the velocity gradient of the liquid lead bismuth is sharply reduced, determining the region as a velocity boundary layer region, comparing the velocity boundary layer region with the velocity boundary layer region determined in the step 5-1, and supplementing the region which is not covered in the step 5-1:
step 5-3: extracting the coordinates of the velocity boundary layer area determined in the step 5-2, and calculating to obtain the thickness of the velocity boundary layer, wherein the calculation formula is as follows:
δv-velocity boundary layer thickness, m;
x, y, z-velocity boundary layer boundary coordinates;
step 6: according to the distribution information of the temperature field of the liquid lead bismuth sweepback spiral tube bundle calculation domain obtained in the step 3, the temperature gradient distribution condition of the liquid lead bismuth sweepback spiral tube bundle calculation domain is determined, the position of a temperature boundary layer is determined, and the thickness of the temperature boundary layer is extracted, and the specific steps are as follows:
step 6-1: calculating the distribution information of the domain temperature field according to the liquid lead bismuth sweepforward spiral tube bundle obtained in the step 3, determining the temperature boundary layer area, extracting the position information of the temperature boundary layer, and judging that the fluid is positioned in the temperature boundary layer when the temperature difference between the liquid lead bismuth temperature and the wall surface temperature reaches below 99 percent of the temperature difference between the main flow temperature and the wall surface temperature:
T-Tw＝0.99(T∞-Tw) (12)
T∞-liquid lead bismuth main stream temperature, K;
Tw-cooling the wall temperature, K;
step 6-2: according to the temperature gradient distribution condition of the fluid section of the liquid lead bismuth sweepback spiral tube bundle obtained in the step 4, when the temperature gradient of the liquid lead bismuth is sharply reduced, determining that the area is a temperature boundary layer area, comparing the temperature boundary layer area with the temperature boundary layer area determined in the step 6-1, and supplementing the area which is not covered in the step 6-1:
step 6-3: extracting the coordinates of the temperature boundary layer area determined in the step 6-2, and calculating to obtain the thickness of the temperature boundary layer, wherein the calculation formula is as follows:
δt-temperature boundary layer thickness, m;
x, y, z-velocity boundary layer boundary coordinates;
the invention has the following advantages and beneficial effects:
1. the method provides a liquid metal sweepforward spiral tube bundle speed temperature boundary layer calculation method, the model is independent, the method has strong universality, and the method is suitable for most of the existing fluid dynamics calculation programs;
2. the method can carry out numerical simulation calculation aiming at the liquid metal flow heat transfer boundary layer under the complex structure, and provides method reference for selection of the main flow position in the related experiment.
3. The method is described in detail aiming at liquid lead and bismuth, and can be applied to the exploration of speed and temperature boundary layers of other kinds of liquid metals.
Drawings
FIG. 1 is a flow chart of the method of the present invention.
FIG. 2 is a schematic drawing of a selected cross-section of a boundary layer study.
FIG. 3 is a simplified mesh model of a pressure shell region.
FIG. 4 is a schematic drawing of a selected cross-section of a boundary layer study.
Detailed Description
The invention is described in further detail below with reference to the following figures and detailed description:
the invention provides a method for calculating a speed and temperature boundary layer of a liquid lead bismuth swept-out spiral tube bundle as shown in figure 1, which comprises the following specific steps:
step 1: establishing a simplified model of a liquid lead bismuth swept-out spiral tube bundle calculation domain by using a geometric model establishing software SOLIDWORKS, obtaining a fluid domain geometric model of the liquid lead bismuth swept-out spiral tube bundle, simplifying the calculation region by considering model symmetry and calculation cost, selecting a representative calculation unit, and simplifying the liquid lead bismuth swept-out spiral tube bundle calculation domain into a circumferential 1/12 liquid lead bismuth swept-out spiral tube bundle cylindrical sleeve as shown in FIG. 2;
step 2: applying mesh division software ANSYS-ICEM, on the basis of the liquid lead bismuth sweepout spiral tube bundle fluid domain geometric model obtained in the step 1, carrying out mesh division on the liquid lead bismuth sweepout spiral tube bundle fluid domain to obtain a liquid lead bismuth sweepout spiral tube bundle fluid domain mesh model, wherein the liquid lead bismuth sweepout spiral tube bundle fluid domain mesh model comprises boundary layer area meshes and integral meshes, and the boundary layer meshes are added on the outer side of the spiral tube bundle and the inner wall surface and the outer wall surface of the cylindrical sleeve, as shown in fig. 3;
and step 3: and (3) introducing the liquid lead bismuth sweepout spiral tube beam body area grid model obtained in the step (2) into computational fluid dynamics software ANSYS-FLUENT to perform numerical simulation calculation in the process of liquid lead bismuth flowing heat exchange, and specifically comprising the following steps:
step 3-1: introducing the liquid lead bismuth sweepback spiral tube beam flow body area grid model into computational fluid mechanics software ANSYS-FLUENT, setting the surface boundary of the upper part of the model as the inlet boundary of the liquid lead bismuth, and setting the speed, the temperature and the pressure of the liquid lead bismuth at the inlet position according to actual conditions;
step 3-2: leading the liquid lead bismuth sweepback spiral tube beam flow body area grid model into computational fluid mechanics software ANSYS-FLUENT, setting the lower surface boundary as a liquid lead bismuth outlet boundary, and setting the liquid lead bismuth pressure at the outlet position according to actual conditions;
step 3-3: setting the boundary of the outer surface of the liquid lead bismuth sweepout spiral tube beam flow body area grid model tube bundle as a constant wall temperature boundary in computational fluid dynamics software ANSYS-FLUENT, simulating the process of cooling the liquid lead bismuth by the tube bundle, setting the wall surfaces of the inner sleeve and the outer sleeve as heat insulation boundaries, and not exchanging heat between the wall surface of the inner sleeve and the wall surface of the outer sleeve;
step 3-4: the physical parameters of the liquid lead bismuth obtained by calculation comprise the density of the liquid lead bismuth, the dynamic viscosity of the liquid lead bismuth alloy, the specific heat capacity of the liquid lead bismuth alloy and the dynamic thermal conductivity of the liquid lead bismuth alloy, and the density of the liquid lead bismuth is calculated as follows:
ρLBE＝11096-1.3236T (1)
t-temperature of liquid lead bismuth, K
The calculation formula of the specific heat capacity of the liquid lead-bismuth alloy is as follows:
cp,LBE＝159-2.72×10-2T+7.12×10-6T2 (2)
the dynamic viscosity of the liquid lead-bismuth alloy is calculated as follows:
the calculation formula of the dynamic thermal conductivity of the liquid lead-bismuth alloy is as follows:
λLBE＝3.61+1.517×10-2T-1.741×10-6T2 (4)
step 3-5: selecting a turbulent flow Plantt digital model suitable for liquid lead bismuth flow heat exchange simulation, loading the turbulent flow Plantt digital model into computational fluid dynamics software, and correcting the turbulent flow Plantt digital model in the computational fluid dynamics software FLUENT, wherein the selected model is as follows:
Prt＝4.12Pe≤1000 (6)
pe-liquid lead bismuth flow peclet number;
Prt-turbulent prandtl number;
step 3-6: solving the mass, momentum, energy and component transport equation of the liquid lead bismuth to obtain a liquid lead bismuth sweepforward spiral tube bundle calculation domain velocity field and a liquid lead bismuth temperature field;
pe-liquid lead bismuth flow peclet number;
and 4, step 4: establishing a cross section superposed with the spiral tube bundle by utilizing post-processing software TECPLOT from the calculation result of the speed field and the temperature field of the liquid lead bismuth swept-out spiral tube bundle calculation domain obtained in the step 3, as shown in FIG. 4, extracting the information of the fluid cross section speed and the temperature field distribution of the liquid lead bismuth swept-out spiral tube bundle, further calculating the temperature field speed field distribution to obtain the fluid cross section speed and the temperature gradient distribution condition of the liquid lead bismuth swept-out spiral tube bundle, and solving the equation of the speed gradient and the temperature gradient into the equation;
u, V, W-x, y, z direction component velocity, m/s;
t-liquid lead bismuth temperature, K;
a-velocity gradient, s-1；
b-temperature gradient, K/m;
and 5: calculating domain speed field distribution information according to the liquid lead bismuth swept-out spiral tube bundle obtained in the step 3, calculating domain speed gradient distribution conditions according to the liquid lead bismuth swept-out spiral tube bundle, determining the position of a speed boundary layer, acquiring boundary coordinates by using TECPLOT software, and calculating to obtain the thickness of the speed boundary layer, wherein the specific steps are as follows:
step 5-1: calculating the distribution information of the domain velocity field according to the liquid lead bismuth sweepback spiral tube bundle obtained in the step 3, determining a velocity boundary layer region, extracting the position information of the velocity boundary layer, and determining the region as the velocity boundary layer region when the local flow velocity of the liquid lead bismuth is less than 99% of the main flow velocity according to the judgment basis:
U＝0.99U∞ (10)
U∞-the main flow velocity of liquid lead bismuth, m/s;
step 5-2: according to the velocity gradient distribution condition of the fluid section of the liquid lead bismuth sweepback spiral tube bundle obtained in the step 4, when the velocity gradient of the liquid lead bismuth is sharply reduced, determining the region as a velocity boundary layer region, comparing the velocity boundary layer region with the velocity boundary layer region determined in the step 5-1, and supplementing the region which is not covered in the step 5-1:
step 5-3: acquiring the coordinate of the speed boundary layer region determined in the step 5-2 by using TECPLOT software, and calculating to obtain the thickness of the speed boundary layer, wherein the calculation formula is as follows:
δv-velocity boundary layer thickness, m;
x, y, z-velocity boundary layer boundary coordinates;
step 6: calculating domain temperature field distribution information according to the liquid lead bismuth outside-swept spiral tube bundle obtained in the step 3, calculating domain temperature gradient distribution conditions according to the liquid lead bismuth outside-swept spiral tube bundle, determining the position of a temperature boundary layer, acquiring boundary coordinates by using TECPLOT software, and extracting the thickness of the temperature boundary layer, wherein the method specifically comprises the following steps:
step 6-1: calculating the distribution information of the domain temperature field according to the liquid lead bismuth sweepforward spiral tube bundle obtained in the step 3, determining the temperature boundary layer area, extracting the position information of the temperature boundary layer, and judging that the fluid is positioned in the temperature boundary layer when the temperature difference between the liquid lead bismuth temperature and the wall surface temperature reaches below 99 percent of the temperature difference between the main flow temperature and the wall surface temperature:
T-Tw＝0.99(T∞-Tw) (12)
T∞-liquid lead bismuth main stream temperature, K;
Tw-cooling the wall temperature, K;
step 6-2: according to the temperature gradient distribution condition of the fluid section of the liquid lead bismuth sweepback spiral tube bundle obtained in the step 4, when the temperature gradient of the liquid lead bismuth is sharply reduced, determining that the area is a temperature boundary layer area, comparing the temperature boundary layer area with the temperature boundary layer area determined in the step 6-1, and supplementing the area which is not covered in the step 6-1:
step 6-3: acquiring the coordinates of the temperature boundary layer determined in the step 6-2 by using TECPLOT software, and calculating to obtain the thickness of the temperature boundary layer, wherein the calculation formula is as follows:
δt-temperature boundary layer thickness, m;
x, y, z-velocity boundary layer boundary coordinates.
Claims (1)
1. A method for calculating the speed and temperature boundary layer of a sweepforward spiral tube bundle of liquid lead and bismuth is characterized by comprising the following steps of: aiming at the flow of the liquid lead bismuth sweepforward spiral tube bundle, the calculation of a speed temperature boundary layer is realized by a fluid dynamics calculation program;
the method comprises the following steps:
step 1: establishing a simplified model of the calculation domain of the liquid lead bismuth swept-out spiral tube bundle by using a geometric model establishing software SOLIDWORKS to obtain a fluid domain geometric model of the liquid lead bismuth swept-out spiral tube bundle, selecting a representative calculation unit in consideration of the model symmetry, and simplifying the calculation domain of the liquid lead bismuth swept-out spiral tube bundle into a circumferential 1/12 cylindrical sleeve of the liquid lead bismuth swept-out spiral tube bundle;
step 2: performing mesh division on the fluid domain of the liquid lead bismuth swept-out spiral tube bundle by using mesh division software ANSYS-ICEM on the basis of the geometric model of the fluid domain of the liquid lead bismuth swept-out spiral tube bundle obtained in the step 1 to obtain a beam fluid body domain mesh model of the liquid lead bismuth swept-out spiral tube bundle, wherein the beam fluid body domain mesh model of the liquid lead bismuth swept-out spiral tube bundle comprises boundary layer area meshes and integral meshes;
and step 3: performing numerical simulation calculation on the flowing heat exchange process of the liquid lead bismuth in the liquid lead bismuth sweepforward spiral tube beam body area grid model obtained in the step 2, and specifically comprising the following steps:
step 3-1: setting the boundary of the upper surface of the beam body area grid model of the sweepforward spiral tube of the liquid lead bismuth as the inlet boundary of the liquid lead bismuth, and setting the speed, the temperature and the pressure of the liquid lead bismuth at the inlet position according to actual conditions;
step 3-2: setting the boundary of the lower surface of the beam body area grid model of the sweepforward spiral tube of the liquid lead bismuth as the outlet boundary of the liquid lead bismuth, and setting the pressure of the liquid lead bismuth at the outlet position according to actual conditions;
step 3-3: setting the boundary of the outer surface of the tube bundle of the liquid lead bismuth sweepback spiral tube beam body area grid model as a constant wall temperature boundary, simulating the process of cooling the liquid lead bismuth by the tube bundle, setting the wall surfaces of the inner sleeve and the outer sleeve as heat insulation boundaries, and not exchanging heat between the wall surfaces of the inner sleeve and the outer sleeve and the liquid lead bismuth;
step 3-4: the physical parameters of the liquid lead bismuth obtained by calculation comprise the density of the liquid lead bismuth, the dynamic viscosity of the liquid lead bismuth alloy, the specific heat capacity of the liquid lead bismuth alloy and the dynamic thermal conductivity of the liquid lead bismuth alloy, and the density of the liquid lead bismuth is calculated as follows:
ρLBE＝11096-1.3236T (1)
t-temperature of liquid lead bismuth, K
The calculation formula of the specific heat capacity of the liquid lead-bismuth alloy is as follows:
cp,LBE＝159-2.72×10-2T+7.12×10-6T2 (2)
the dynamic viscosity of the liquid lead-bismuth alloy is calculated as follows:
the calculation formula of the dynamic thermal conductivity of the liquid lead-bismuth alloy is as follows:
λLBE＝3.61+1.517×10-2T-1.741×10-6T2 (4)
step 3-5: selecting a turbulent flow Plantt digital model suitable for liquid lead bismuth flow heat exchange simulation, loading the turbulent flow Plantt digital model into computational fluid dynamics software, and correcting the turbulent flow Plantt digital model in the computational fluid dynamics software FLUENT, wherein the selected model is as follows:
Prt＝4.12 Pe≤1000 (6)
pe-liquid lead bismuth flow peclet number;
Prt-turbulent prandtl number;
step 3-6: solving the mass, momentum, energy and component transport equation of the liquid lead bismuth to obtain a liquid lead bismuth sweepforward spiral tube bundle calculation domain velocity field and a liquid lead bismuth temperature field;
and 4, step 4: establishing a cross section superposed with the spiral tube bundle by utilizing post-processing software TECPLOT from the calculation result of the speed field and the temperature field of the liquid lead bismuth swept-out spiral tube bundle calculation domain obtained in the step 3, extracting the fluid cross section speed and the temperature field distribution information of the liquid lead bismuth swept-out spiral tube bundle, calculating the fluid cross section speed and the temperature gradient distribution condition of the liquid lead bismuth swept-out spiral tube bundle, and solving the equation of the speed gradient and the temperature gradient;
u, V, W-x, y, z direction component velocity, m/s;
t-liquid lead bismuth temperature, K;
a-velocity gradient, s-1；
b-temperature gradient, K/m;
and 5: calculating the domain velocity field distribution information according to the liquid lead bismuth sweepback spiral tube bundle obtained in the step 3, calculating the domain velocity gradient distribution condition according to the liquid lead bismuth sweepback spiral tube bundle, determining the position of a velocity boundary layer, and extracting the thickness of the velocity boundary layer, wherein the specific steps are as follows:
step 5-1: calculating the distribution information of the domain velocity field according to the liquid lead bismuth sweepback spiral tube bundle obtained in the step 3, determining a velocity boundary layer area, extracting the position information of the velocity boundary layer, and determining the area as the velocity boundary layer area when the local flow velocity of the liquid lead bismuth is less than 99% of the main flow velocity according to the judgment basis:
U＝0.99U∞ (10)
U∞-the main flow velocity of liquid lead bismuth, m/s;
step 5-2: according to the velocity gradient distribution condition of the fluid section of the liquid lead bismuth sweepback spiral tube bundle obtained in the step 4, when the velocity gradient of the liquid lead bismuth is sharply reduced, determining the region as a velocity boundary layer region, comparing the velocity boundary layer region with the velocity boundary layer region determined in the step 5-1, and supplementing the region which is not covered in the step 5-1:
step 5-3: extracting the coordinates of the speed boundary layer area determined in the step 5-2, and calculating to obtain the thickness of the speed boundary layer, wherein the calculation formula is as follows:
δv-velocity boundary layer thickness, m;
x, y, z-velocity boundary layer boundary coordinates;
step 6: according to the temperature field distribution information of the calculation domain of the liquid lead bismuth swept-out spiral tube bundle obtained in the step 3, the temperature gradient distribution condition of the calculation domain of the liquid lead bismuth swept-out spiral tube bundle is determined, the position of a temperature boundary layer is determined, and the thickness of the temperature boundary layer is extracted, and the method specifically comprises the following steps:
step 6-1: calculating the distribution information of the domain temperature field according to the liquid lead bismuth sweepforward spiral tube bundle obtained in the step 3, determining the temperature boundary layer area, extracting the position information of the temperature boundary layer, and judging that the fluid is positioned in the temperature boundary layer when the temperature difference between the liquid lead bismuth temperature and the wall surface temperature reaches below 99 percent of the temperature difference between the main flow temperature and the wall surface temperature:
T-Tw＝0.99(T∞-Tw) (12)
T∞-liquid lead bismuth main stream temperature, K;
Tw-cooling the wall temperature, K;
step 6-2: according to the temperature gradient distribution condition of the fluid section of the liquid lead bismuth sweepback spiral tube bundle obtained in the step 4, when the temperature gradient of the liquid lead bismuth is sharply reduced, determining the area as a temperature boundary layer area, comparing the temperature boundary layer area with the temperature boundary layer area determined in the step 6-1, and supplementing the area which is not covered in the step 6-1:
step 6-3: extracting the coordinates of the temperature boundary layer area determined in the step 6-2, and calculating to obtain the thickness of the temperature boundary layer, wherein the calculation formula is as follows:
δt-temperature boundary layer thickness, m;
x, y, z-velocity boundary layer boundary coordinates.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110784158.2A CN113486482B (en) | 2021-07-12 | 2021-07-12 | Liquid lead bismuth sweepforward spiral tube bundle speed and temperature boundary layer calculation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110784158.2A CN113486482B (en) | 2021-07-12 | 2021-07-12 | Liquid lead bismuth sweepforward spiral tube bundle speed and temperature boundary layer calculation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113486482A true CN113486482A (en) | 2021-10-08 |
| CN113486482B CN113486482B (en) | 2022-10-28 |
Family
ID=77938122
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110784158.2A Active CN113486482B (en) | 2021-07-12 | 2021-07-12 | Liquid lead bismuth sweepforward spiral tube bundle speed and temperature boundary layer calculation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113486482B (en) |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2460160C1 (en) * | 2011-07-11 | 2012-08-27 | Федеральное государственное унитарное предприятие "Научно-исследовательский технологический институт имени А.П. Александрова" | Cleaning and deactivation method of reactor plant equipment with liquid-metal lead-bismuth heat carrier |
| CN102654513A (en) * | 2012-03-30 | 2012-09-05 | 中国科学院合肥物质科学研究院 | Method for measuring speed distribution of liquid metal boundary layer under high-temperature strong magnetic field |
| US20120245903A1 (en) * | 2011-03-23 | 2012-09-27 | Desktop Aeronautics, Inc. | Generating inviscid and viscous fluid flow simulations over a surface using a quasi-simultaneous technique |
| CN105627964A (en) * | 2016-01-20 | 2016-06-01 | 四川大学 | Sound field enhanced air flow drying and boundary layer measurement integration experiment system |
| CN108444900A (en) * | 2018-05-31 | 2018-08-24 | 中国科学院近代物理研究所 | Control the corrosion of oxygen liquid lead bismuth and ion irradiation cooperating research experimental provision and method |
| CN111261232A (en) * | 2020-01-14 | 2020-06-09 | 西安热工研究院有限公司 | Reactor primary loop coolant flow field, temperature field and stress field coupling calculation method |
| CN111832201A (en) * | 2020-06-23 | 2020-10-27 | 西安交通大学 | Method for calculating cooling and evaporation coupling of two sides of liquid metal spiral pipe steam generator |
| CN111881631A (en) * | 2020-07-22 | 2020-11-03 | 西安交通大学 | Oxygen concentration diffusion numerical simulation method in lead-cooled fast reactor oxygen control process |
| US20210124861A1 (en) * | 2019-10-28 | 2021-04-29 | Vorcat, Inc. | System and method for simulating turbulence |
| CN113095004A (en) * | 2021-03-29 | 2021-07-09 | 西安交通大学 | Liquid metal flow heat transfer calculation method |
-
2021
- 2021-07-12 CN CN202110784158.2A patent/CN113486482B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120245903A1 (en) * | 2011-03-23 | 2012-09-27 | Desktop Aeronautics, Inc. | Generating inviscid and viscous fluid flow simulations over a surface using a quasi-simultaneous technique |
| RU2460160C1 (en) * | 2011-07-11 | 2012-08-27 | Федеральное государственное унитарное предприятие "Научно-исследовательский технологический институт имени А.П. Александрова" | Cleaning and deactivation method of reactor plant equipment with liquid-metal lead-bismuth heat carrier |
| CN102654513A (en) * | 2012-03-30 | 2012-09-05 | 中国科学院合肥物质科学研究院 | Method for measuring speed distribution of liquid metal boundary layer under high-temperature strong magnetic field |
| CN105627964A (en) * | 2016-01-20 | 2016-06-01 | 四川大学 | Sound field enhanced air flow drying and boundary layer measurement integration experiment system |
| CN108444900A (en) * | 2018-05-31 | 2018-08-24 | 中国科学院近代物理研究所 | Control the corrosion of oxygen liquid lead bismuth and ion irradiation cooperating research experimental provision and method |
| US20210124861A1 (en) * | 2019-10-28 | 2021-04-29 | Vorcat, Inc. | System and method for simulating turbulence |
| CN111261232A (en) * | 2020-01-14 | 2020-06-09 | 西安热工研究院有限公司 | Reactor primary loop coolant flow field, temperature field and stress field coupling calculation method |
| CN111832201A (en) * | 2020-06-23 | 2020-10-27 | 西安交通大学 | Method for calculating cooling and evaporation coupling of two sides of liquid metal spiral pipe steam generator |
| CN111881631A (en) * | 2020-07-22 | 2020-11-03 | 西安交通大学 | Oxygen concentration diffusion numerical simulation method in lead-cooled fast reactor oxygen control process |
| CN113095004A (en) * | 2021-03-29 | 2021-07-09 | 西安交通大学 | Liquid metal flow heat transfer calculation method |
Non-Patent Citations (5)
| Title |
|---|
| NIRMAL S. RAJPUT等: "Liquid-gas heat exchanger for low pressure refrigerant application", 《2017 INTERNATIONAL CONFERENCE ON ELECTRONICS PACKAGING (ICEP)》 * |
| 孙畅等: "基于CFD方法的铅铋冷却燃料棒束的热工水力特性分析", 《原子能科学技术》 * |
| 席文宣等: "圆管内液态铅铋合金强制对流换热特性研究", 《工程热物理学报》 * |
| 杨宇鹏等: "液态金属螺旋管式直流蒸汽发生器数值模拟研究", 《原子能科学技术--网络首发版》 * |
| 王淑香等: "超临界压力下CO_2在螺旋管中沿程传热的实验研究", 《核动力工程》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113486482B (en) | 2022-10-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Cong et al. | Study on secondary side flow of steam generator with coupled heat transfer from primary to secondary side | |
| CN105067661B (en) | Gas liquid exchanger heat transfer coefficient determining device | |
| Arat et al. | A comprehensive numerical investigation of unsteady-state two-phase flow in gravity assisted heat pipe enclosure | |
| Qu et al. | Experimental and numerical investigations on the air–steam mixture bubble condensation characteristics in stagnant cool water | |
| CN113609795B (en) | Three-dimensional calculation method for thermal hydraulic characteristics of high-temperature gas cooled reactor spiral tube steam generator | |
| CN111141400B (en) | Method for measuring temperature of pipe wall of thermal fatigue sensitive area of bent pipe of nuclear power station | |
| CN104655549B (en) | Low-temperature acid dew point corrosion experimental device for power station boiler flue gas | |
| CN112016214B (en) | Heat transfer calculation method and system for buried pipe heat exchanger based on resistance-capacitance model | |
| CN113486482B (en) | Liquid lead bismuth sweepforward spiral tube bundle speed and temperature boundary layer calculation method | |
| Zhang et al. | Unsteady-state thermal stress and thermal deformation analysis for a pressurizer surge line subjected to thermal stratification based on a coupled CFD-FEM method | |
| Liu et al. | Numerical study on the heat transfer characteristics of a liquid lead–bismuth eutectic in a D-type channel | |
| Zhang et al. | Analytical modeling for vapor condensation in the presence of noncondensable gas and experimental validation | |
| CN109063370B (en) | Method for predicting water depth average water flow direction of curved riverway of flood plain | |
| Wang et al. | Influence of buoyancy on turbulent mixed convection of LBE in a uniform cooled inclined tube | |
| Lahey Jr | The prediction of phase distribution and separation phenomena using two-fluid models | |
| Al-Rifai et al. | Heat and mass transfer correlations for staggered nanoporous membrane tubes in flue gas crossflow | |
| Gorpinyak et al. | Vapor–gas mixture condensation in tubes | |
| Lubis et al. | Design of a heat exchanger of three concentric tube layer on contrary flow | |
| Tao et al. | Heat transfer evaluation of saturated pool boiling outside slender tube bundle in passive residual heat removal system of molten salt reactor | |
| Zhao et al. | Experimental study and analysis on the interfacial drag of two-phase flow in porous media | |
| Du et al. | Numerical study of swirling effect on enhancing the bath dynamics in the top-submerged-lance furnace | |
| CN115455847A (en) | Method for judging generation point of supercooling boiling clean steam | |
| Wongwises et al. | Heat-mass transfer and flow characteristics of two-phase countercurrent annular flow in a vertical pipe | |
| Habib et al. | On the development of deadleg criterion | |
| Chu et al. | Numerical study on supercritical co 2 in printed circuit heat exchanger with asymmetrical airfoil fins |
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 |