CN115130323B - Interphase resistance analysis method and device suitable for rod bundle channel - Google Patents

Interphase resistance analysis method and device suitable for rod bundle channel Download PDF

Info

Publication number
CN115130323B
CN115130323B CN202210873680.2A CN202210873680A CN115130323B CN 115130323 B CN115130323 B CN 115130323B CN 202210873680 A CN202210873680 A CN 202210873680A CN 115130323 B CN115130323 B CN 115130323B
Authority
CN
China
Prior art keywords
phase
flow
resistance
inter
bundle channel
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202210873680.2A
Other languages
Chinese (zh)
Other versions
CN115130323A (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.)
Nuclear Power Institute of China
Original Assignee
Nuclear Power Institute of China
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 Nuclear Power Institute of China filed Critical Nuclear Power Institute of China
Priority to CN202210873680.2A priority Critical patent/CN115130323B/en
Publication of CN115130323A publication Critical patent/CN115130323A/en
Application granted granted Critical
Publication of CN115130323B publication Critical patent/CN115130323B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • 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
    • 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)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Algebra (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Operations Research (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Probability & Statistics with Applications (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Monitoring And Testing Of Nuclear Reactors (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The application discloses an interphase resistance analysis method and device suitable for a rod bundle channel, comprising the following steps: judging the flow pattern of the rod bundle channel according to the rod bundle channel, and obtaining a flow pattern judgment result; according to the flow pattern judgment result, calculating the phase-to-phase resistance by adopting different phase-to-phase resistance models aiming at different flow patterns; the flow pattern judgment result includes bubble flow, cap flow, and annular flow. According to the application, based on drag model development, the influence of strong stirring characteristics in the rod bundle channel is considered on drag coefficient and phase interface concentration, and according to the structural form of the rod bundle channel, only the interface concentration and drag coefficient under different flow patterns are corrected, so that the consistency of the inter-phase resistance models in the channels with different structures is ensured; meanwhile, the influence of strong stirring characteristics in the bundle channels is considered on the drag coefficient and the phase interface concentration, so that the inter-phase resistance calculation and analysis of the bundle channels in the nuclear reactor core fuel assembly are more accurate.

Description

Interphase resistance analysis method and device suitable for rod bundle channel
Technical Field
The application relates to the field of reactor thermal hydraulic design and safety analysis, in particular to an interphase resistance analysis method and device suitable for a rod bundle channel.
Background
When the gas phase and the vapor phase are mixed and flow together, the difference of the physical properties of the two phases enables the flow characteristics of the gas phase and the vapor phase to be quite different, so that the interaction force exists between the gas phase and the vapor phase, and the interaction force is interphase resistance. The interphase resistance has important influence on the deformation of a two-phase interface, the phase distribution, the flow pressure drop and the like, and further influences the thermodynamic and hydraulic characteristics of the reactor. Accurate simulation of interphase resistance characteristic-related phenomena requires an accurate interphase resistance model. As an important basis for basic equations and model selection in the two-phase flow and heat transfer calculation process, the inter-phase resistance model is a basic composition model of the current main stream thermodynamic and hydraulic programs (such as RELAP5, CATHARE, COBRA-TF and the like).
The nuclear reactor core fuel assembly is mainly of a complex rod bundle structure, and is used as a special open channel, wherein the rod bundle structure is different from a conventional closed channel (such as a circular tube and a rectangular flow channel), fluid has stronger mixing strength, a phase interface is easier to deform, and phase distribution is more complex. The main flow system analysis program mainly adopts a drift flow model and a drag force model to simulate the inter-phase resistance in two-phase flow, but most of the main flow system analysis program adopts a flow pattern division mode which is the same as a circular pipe, so that the stirring characteristics among all channels in a complex channel are difficult to embody.
Disclosure of Invention
The application aims to solve the technical problem that the inter-phase resistance of a rod bundle channel in the existing nuclear reactor core fuel assembly is analyzed without considering the influence of the rod bundle channel stirring effect on a phase interface, so that the inter-phase resistance is calculated inaccurately.
The application aims to provide an interphase resistance analysis method and an interphase resistance analysis device suitable for a rod bundle channel, and develops a set of interphase resistance calculation analysis method suitable for bubble flow, cap flow and annular flow in the rod bundle channel. The inter-phase resistance calculation of the rod bundle channel is identical to the inter-phase resistance calculation model of the circular tube in form, and the influence of the strong stirring characteristic in the rod bundle channel is considered on the drag coefficient and the phase interface concentration, which is the core of the patent. The application has wider application range and higher prediction precision of the interphase resistance model, so as to accurately predict the complex two-phase flow characteristic in the reactor core.
The application is realized by the following technical scheme:
in a first aspect, the present application provides a method of interphase resistance analysis for a bundle channel, the method comprising:
judging the flow pattern of the rod bundle channel according to the rod bundle channel, and obtaining a flow pattern judgment result;
according to the flow pattern judgment result, calculating the phase-to-phase resistance by adopting different phase-to-phase resistance models aiming at different flow patterns; wherein: the flow pattern judgment result includes bubble flow, cap flow and annular flow.
The working principle is as follows: the interphase resistance analysis based on the bundle channels in the existing nuclear reactor core fuel assembly does not consider the influence of the bundle channel stirring effect on the phase interface, so that the interphase resistance calculation is not accurate enough.
According to the application, a set of inter-phase resistance calculation analysis method suitable for bubble flow, cap flow and annular flow in the rod bundle channel is developed, on one hand, the judgment of the flow pattern of the rod bundle channel is considered, the judgment result of the flow pattern is obtained, the division judgment of the flow pattern is more accurate, and further, the follow-up inter-phase resistance calculation of the corresponding flow pattern is more accurate. On the other hand, the application is mainly developed based on a drag force model, and the drag force coefficient and the phase interface concentration in the model are mainly closed by a semi-empirical-semi-theoretical relation. The inter-phase resistance calculation of the rod bundle channel is identical to the inter-phase resistance calculation model of the circular tube in form, but the influence of the strong stirring property in the rod bundle channel is considered on the drag coefficient and the phase interface concentration, which is the core of the patent; according to the structure form of the rod bundle channel, only the interface concentration and the drag coefficient under different flow patterns are corrected, so that the consistency of inter-phase resistance models in different structure channels is ensured; meanwhile, the influence of strong stirring characteristics in the bundle channels is considered on the drag coefficient and the phase interface concentration, so that the inter-phase resistance calculation and analysis of the bundle channels in the nuclear reactor core fuel assembly are more accurate.
The application has wider application range and higher prediction precision of the interphase resistance model, so as to accurately predict the complex two-phase flow characteristic in the reactor core.
Further, judging the flow pattern of the rod bundle channel, and obtaining a flow pattern judgment result; the judgment basis is as follows:
judging the flow pattern according to the flow rate and the cavitation share of the rod bundle channel;
if the cavitation share of the rod bundle channel is smaller than 0.2, the flow pattern judgment result is bubble flow;
if the flow rate of the rod bundle channel is less than or equal to 3000 kg/(m) 2 S) and the cavitation fraction is between 0.2 and 0.85, or the flow rate of the rod bundle channel is greater than 3000 kg/(m) 2 S) and the cavitation proportion is between 0.4 and 0.85, and the judgment result of the flow pattern is cap-shaped flow;
if the cavitation proportion of the rod bundle channel is greater than 0.85, the flow pattern judgment result is annular flow.
Further, calculating the inter-phase resistance by adopting a first inter-phase resistance model for the bubble flow, wherein the first inter-phase resistance model is obtained by taking the concentration of a phase interface formed by the stirring characteristic of a rod bundle channel into consideration;
calculating inter-phase resistance by adopting a second inter-phase resistance model aiming at cap-shaped flow;
and calculating the inter-phase resistance by adopting a third inter-phase resistance model aiming at the annular flow.
Further, the first inter-phase resistance model and the second inter-phase resistance model are both related by the following formula:
wherein F is i Is inter-phase resistance; ρ c Is the density of continuous phase, namely the density rho of liquid phase f Or gas/vapor phase density ρ g ；v r For the relative velocity between the gas/vapor phase and the liquid phase, i.e. v r ＝v g -v f ，v g And v f Gas/vapor phase and liquid phase velocities, respectively; c (C) D For drag coefficient, ishii is used&The Chawlea (1979) relation is calculated; a, a gf Is the phase interface concentration.
Further, the phase interface concentration calculation formula of the first interphase resistance model is as follows:
wherein α is the cavitation fraction; sigma is the coefficient of viscosity of the adhesive,calculated by a Chexal-Lellouche relational expression; c (C) j1 Is the first key parameter obtained by considering the influence of the whipping of a rod bundle channel on the bubble characteristics, C j2 Is a second key parameter that is derived in view of the impact of the rod-channel whipping on the bubble characteristics.
Further, the second inter-phase resistance model corrects the bubble-shaped fluid phase interface concentration according to the cap-shaped fluid flow pattern interface concentration and the ratio of the cavitation share and the actual cavitation share at the transition of the flow pattern;
the phase interface concentration calculation formula of the second phase-to-phase resistance model is as follows:
wherein alpha is bc A cavitation fraction that is converted from bubble flow to cap flow;
the drag coefficient calculation formula of the second inter-phase resistance model is as follows:
wherein C is D For drag coefficient, alpha bc Cavitation fraction for hat flow to annular flow transition, Δρ is the liquid phase density ρ f Density ρ of gas/vapor phase g And (3) a difference.
Further, the third inter-phase resistance model comprises an air/steam core inter-phase resistance sub-model and a liquid drop inter-phase resistance sub-model; the expression of the third inter-phase resistance model is:
F i ＝F i_c +F i_d (5)
wherein F is i_c Is the inter-phase resistance of the gas/steam core, F i_d Is the inter-phase resistance of the liquid film, F i_c And F i_d The same relation calculation as the first interphase resistance model is adopted, and the corresponding drag coefficient and interface concentration are corrected.
Further, the correcting the corresponding drag coefficient and phase interface concentration includes:
for the gas/steam core, the drag coefficient C of the gas/steam core D_c Phase interface concentration alpha of gas/steam core ff The calculation formula is as follows:
1) When the vapor phase Reynolds number Re g When the number of the times is less than 500,
2) When the vapor phase Reynolds number Re g When the number of the carbon atoms is more than 1500,
C D_c ＝0.02{1+150[1-(1-α ff ) 0.5 ]} (7)
wherein alpha is f As liquid film/liquid phase fraction, C f Is a constant related to the structure of the bundle channels; alpha f As cavitation share, v g For vapor phase velocity, alpha ff Is the entrainment share;
3) Vapor phase Reynolds number Re g When the drag coefficient is between 500 and 1500, the calculation formula of the drag coefficient of the gas/steam core adopts the difference between the formula (6) and the formula (7);
concentration of phase interface a gf_c The calculation formula is as follows:
wherein C is j3 To take into account the constant of the influence of the rod-channel stirring on the phase interface.
In a second aspect, the present application further provides an interphase resistance analysis device suitable for a bundle channel, the device supporting the interphase resistance analysis method suitable for a bundle channel; the device comprises:
the flow pattern judging unit is used for judging the flow pattern of the bar bundle channel according to the bar bundle channel and obtaining a flow pattern judging result;
the inter-phase resistance calculation unit is used for calculating the inter-phase resistance by adopting different inter-phase resistance models according to the flow pattern judgment result; wherein: the flow pattern judgment result includes bubble flow, cap flow and annular flow.
Further, the flow pattern judging unit performs the following steps:
judging the flow pattern according to the flow rate and the cavitation share of the rod bundle channel;
if the cavitation share of the rod bundle channel is smaller than 0.2, the flow pattern judgment result is bubble flow;
if the flow rate of the rod bundle channel is less than or equal to 3000 kg/(m) 2 S) and the cavitation fraction is between 0.2 and 0.85, or the flow rate of the rod bundle channel is greater than 3000 kg/(m) 2 S) and the cavitation proportion is between 0.4 and 0.85, and the judgment result of the flow pattern is cap-shaped flow;
if the cavitation proportion of the rod bundle channel is greater than 0.85, the flow pattern judgment result is annular flow.
Compared with the prior art, the application has the following advantages and beneficial effects:
1. according to the application, on one hand, the judgment of the flow pattern of the rod bundle channel is considered, and the flow pattern judgment result is obtained, so that the division judgment of the flow pattern is more accurate, and further, the calculation of the phase-to-phase resistance of the subsequent corresponding flow pattern is more accurate. On the other hand, the application is mainly developed based on a drag force model, and the drag force coefficient and the phase interface concentration in the model are mainly closed by a semi-empirical-semi-theoretical relation. The inter-phase resistance calculation of the rod bundle channel is identical to the inter-phase resistance calculation model of the circular tube in form, but the influence of the strong stirring property in the rod bundle channel is considered on the drag coefficient and the phase interface concentration, which is the core of the patent; according to the structure form of the rod bundle channel, only the interface concentration and the drag coefficient under different flow patterns are corrected, so that the consistency of inter-phase resistance models in different structure channels is ensured; meanwhile, the influence of strong stirring characteristics in the bundle channels is considered on the drag coefficient and the phase interface concentration, so that the inter-phase resistance calculation and analysis of the bundle channels in the nuclear reactor core fuel assembly are more accurate.
2. According to the application, the interface concentration under the bubble flow is corrected according to the stirring effect of the rod bundle channel, the drag coefficient and the interface concentration under the bullet flow are corrected, and the interface concentration under the annular flow is corrected, so that an interphase resistance model applicable to the rod bundle channel is formed, and a basic model support can be provided for the research and development of a thermal hydraulic program based on two fluids-six equations.
3. Compared with the inter-phase resistance model in the circular tube channel, the application has relatively smaller difference of the bar bundle channel model, and is convenient for program integration. The model can improve the prediction accuracy of the phenomena of two-phase flow characteristics, cavitation share distribution and the like of thermal hydraulic characteristics.
Drawings
The accompanying drawings, which are included to provide a further understanding of embodiments of the application and are incorporated in and constitute a part of this specification, illustrate embodiments of the application and together with the description serve to explain the principles of the application. In the drawings:
FIG. 1 is a flow chart of a method for analyzing interphase resistance for a bundle channel according to the present application.
FIG. 2 is a detailed flow chart of a method of analyzing interphase resistance for a bundle channel according to the present application.
FIG. 3 is a schematic diagram of an interphase resistance analyzer suitable for a bundle channel according to the present application.
Detailed Description
For the purpose of making apparent the objects, technical solutions and advantages of the present application, the present application will be further described in detail with reference to the following examples and the accompanying drawings, wherein the exemplary embodiments of the present application and the descriptions thereof are for illustrating the present application only and are not to be construed as limiting the present application.
Example 1
The interphase resistance analysis based on the bundle channels in the existing nuclear reactor core fuel assembly does not consider the influence of the bundle channel stirring effect on the phase interface, so that the interphase resistance calculation is not accurate enough.
According to the application, a set of inter-phase resistance calculation analysis method suitable for bubble flow, cap flow and annular flow in the rod bundle channel is developed, on one hand, the judgment of the flow pattern of the rod bundle channel is considered, the judgment result of the flow pattern is obtained, the division judgment of the flow pattern is more accurate, and further, the follow-up inter-phase resistance calculation of the corresponding flow pattern is more accurate. On the other hand, the application is mainly developed based on a drag force model, and the drag force coefficient and the phase interface concentration in the model are mainly closed by a semi-empirical-semi-theoretical relation. The inter-phase resistance calculation of the rod bundle channel is identical to the inter-phase resistance calculation model of the circular tube in form, but the influence of the strong stirring property in the rod bundle channel is considered on the drag coefficient and the phase interface concentration, which is the core of the patent; according to the structure form of the rod bundle channel, only the interface concentration and the drag coefficient under different flow patterns are corrected, so that the consistency of inter-phase resistance models in different structure channels is ensured; meanwhile, the influence of strong stirring characteristics in the bundle channels is considered on the drag coefficient and the phase interface concentration, so that the inter-phase resistance calculation and analysis of the bundle channels in the nuclear reactor core fuel assembly are more accurate. The application has wider application range and higher prediction precision of the interphase resistance model, so as to accurately predict the complex two-phase flow characteristic in the reactor core.
As shown in fig. 1 and 2, the method for analyzing interphase resistance applied to a bundle channel according to the present application comprises the steps of:
judging the flow pattern of the rod bundle channel according to the rod bundle channel, and obtaining a flow pattern judgment result;
according to the flow pattern judgment result, calculating the phase-to-phase resistance by adopting different phase-to-phase resistance models aiming at different flow patterns; wherein: the flow pattern judgment result includes bubble flow, cap flow and annular flow.
As a further implementation, judging the flow pattern of the rod bundle channel, and obtaining a flow pattern judgment result; the judgment basis is as follows:
judging the flow pattern according to the flow rate and the cavitation share of the rod bundle channel;
if the cavitation share of the rod bundle channel is smaller than 0.2, the flow pattern judgment result is bubble flow;
if the flow rate of the rod bundle channel is less than or equal to 3000 kg/(m) 2 S) and the cavitation fraction is between 0.2 and 0.85, or the flow rate of the rod bundle channel is greater than 3000 kg/(m) 2 S) and the cavitation proportion is between 0.4 and 0.85, and the judgment result of the flow pattern is cap-shaped flow;
if the cavitation proportion of the rod bundle channel is greater than 0.85, the flow pattern judgment result is annular flow.
As a further implementation, the calculation of the inter-phase resistance is performed for the bubble flow using a first inter-phase resistance model obtained by using a phase interface concentration formed by taking into consideration the rod-bundle channel stirring characteristics;
calculating inter-phase resistance by adopting a second inter-phase resistance model aiming at cap-shaped flow;
and calculating the inter-phase resistance by adopting a third inter-phase resistance model aiming at the annular flow.
Specifically, the first inter-phase resistance model and the second inter-phase resistance model are both related by the following formula:
wherein F is i Is inter-phase resistance; ρ c Is the density of continuous phase, namely the density rho of liquid phase f Or gas/vapor phase density ρ g ；v r For the relative velocity between the gas/vapor phase and the liquid phase, i.e. v r ＝v g -v f ，v g And v f Gas/vapor phase and liquid phase velocities, respectively; c (C) D For drag coefficient, ishii is used&The Chawlea (1979) relation is calculated; a, a gf For the phase interface concentration, by taking into consideration the rod-bundle channel stirring characteristics, based on the bubble chord length distribution characteristics, the phase interface concentration is calculated by:
wherein α is the cavitation fraction; sigma is the coefficient of viscosity of the adhesive,calculated by a Chexal-Lellouche relational expression; c (C) j1 Is a first key parameter obtained by considering the influence of the whipping of the rod-bundle channel on the bubble characteristics (the first key parameter is a constant which is obtained by analyzing the relation between the bubble characteristic parameter and the rod-bundle channel size measured by the low flow rate test process), C j2 Is a second key parameter obtained by considering the influence of the whipping of the bundle channel on the bubble characteristics (the second key parameter is a constant which is obtained by analyzing the relation between the bubble characteristic parameter measured by the high flow rate test process and the bundle channel size).
Specifically, the cap flow in the bundle channel is actually a transition flow between bubble flow and bullet flow, and the transition flow comprises cap flow, stirring flow and the like in the general sense. It was found by research that the biggest difference from the bubble flow is that the interface concentration of the flow pattern gradually decreases with increasing cavitation fraction. Based on the above, the second inter-phase resistance model corrects the bubble-shaped fluid phase interface concentration according to the cap-shaped fluid flow type interface concentration closure and the ratio of the void fraction at the transition of the flow pattern to the actual void fraction;
the phase interface concentration calculation formula of the second phase-to-phase resistance model is as follows:
wherein alpha is bc A cavitation fraction that is converted from bubble flow to cap flow;
the drag coefficient calculation formula of the second inter-phase resistance model is as follows:
wherein C is D For drag coefficient, alpha bc Cavitation fraction for hat flow to annular flow transition, Δρ is the liquid phase density ρ f Density ρ of gas/vapor phase g And (3) a difference.
Specifically, the annular flow in the rod bundle channel mainly comprises two parts of an air/steam core and a liquid film. The interphase resistance under this flow pattern is mainly divided into two parts: interphase resistance of gas/steam core and liquid film;
the third inter-phase resistance model comprises an air/steam core inter-phase resistance sub-model and a liquid drop inter-phase resistance sub-model;
the expression of the third inter-phase resistance model is:
F i ＝F i_c +F i_d (5)
wherein F is i_c Is the inter-phase resistance of the gas/steam core, F i_d Is the inter-phase resistance of the liquid film, F i_c And F i_d All adopt the same relation meter as the first interphase resistance modelAnd calculating, and correcting the corresponding drag coefficient and the interface concentration.
The correction of the corresponding drag coefficient and phase interface concentration comprises the following steps:
first, for the gas/steam core, the drag coefficient C of the gas/steam core D_c Phase interface concentration alpha of gas/steam core ff The calculation formula is as follows:
1) When the vapor phase Reynolds number Re g When the number of the times is less than 500,
2) When the vapor phase Reynolds number Re g When the number of the carbon atoms is more than 1500,
C D_c ＝0.02{1+150[1-(1-α ff ) 0.5 ]} (7)
wherein alpha is f As liquid film/liquid phase fraction, C f Is a constant related to the structure of the bundle channels; alpha f As cavitation share, v g For vapor phase velocity, alpha ff Is the entrainment share;
3) Vapor phase Reynolds number Re g When the drag coefficient is between 500 and 1500, the calculation formula of the drag coefficient of the gas/steam core adopts the difference between the formula (6) and the formula (7);
concentration of phase interface a gf_c The calculation formula is as follows:
wherein C is j3 To take into account the constant of the influence of the rod-channel stirring on the phase interface.
Second oneDrag coefficient C for a droplet D_d Using Ishii&The Chawlea (1979) relation is calculated; interface concentration a gf_d The calculation relation is the same as that of the formula (2), and only alpha in the formula (2) is required to be modified into the liquid drop share alpha f 。
The method can be used for calculating the acting force between gas and gas liquid phases in the rod bundle channel, and can be implanted into a thermal hydraulic analysis program based on a two-fluid-six equation. Compared with an inter-phase resistance model in a mainstream system analysis program RELAP5, the application has the following characteristics:
1) The interphase resistance model provided by the application considers the influence of the stirring effect of the rod bundle channel on the phase interface.
2) Adopting a basic form of a drag force model, and considering the influence of a rod bundle channel stirring effect on the interface concentration and drag force coefficient; according to the stirring effect of the rod bundle channel, the interface concentration under the bubble flow is corrected, the drag coefficient and the interface concentration under the bullet flow are corrected, and the interface concentration under the annular flow is corrected, so that an interphase resistance model applicable to the rod bundle channel is formed, and a basic model support can be provided for the research and development of a thermodynamic hydraulic program based on two fluids-six equations.
3) Compared with the inter-phase resistance model in the circular tube channel, the difference of the bar bundle channel model is relatively smaller, and the program integration is facilitated. The model can improve the prediction accuracy of the phenomena of two-phase flow characteristics, cavitation share distribution and the like of thermal hydraulic characteristics.
Example 2
As shown in fig. 3, the difference between the present embodiment and embodiment 1 is that the present embodiment further provides an interphase resistance analysis device suitable for a bundle channel, which supports the interphase resistance analysis method suitable for a bundle channel described in embodiment 1; the device comprises:
the flow pattern judging unit is used for judging the flow pattern of the bar bundle channel according to the bar bundle channel and obtaining a flow pattern judging result;
the inter-phase resistance calculation unit is used for calculating the inter-phase resistance by adopting different inter-phase resistance models according to the flow pattern judgment result; wherein: the flow pattern judgment result includes bubble flow, cap flow and annular flow.
As a further implementation, the flow pattern judging unit performs the following steps:
judging the flow pattern according to the flow rate and the cavitation share of the rod bundle channel;
if the cavitation share of the rod bundle channel is smaller than 0.2, the flow pattern judgment result is bubble flow;
if the flow rate of the rod bundle channel is less than or equal to 3000 kg/(m) 2 S) and the cavitation fraction is between 0.2 and 0.85, or the flow rate of the rod bundle channel is greater than 3000 kg/(m) 2 S) and the cavitation proportion is between 0.4 and 0.85, and the judgment result of the flow pattern is cap-shaped flow;
if the cavitation proportion of the rod bundle channel is greater than 0.85, the flow pattern judgment result is annular flow.
The execution process of each unit is performed according to the steps of the phase-to-phase resistance analysis method applicable to the rectangular channel described in embodiment 1, and in this embodiment, details are not repeated.
It will be appreciated by those skilled in the art that embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The foregoing description of the embodiments has been provided for the purpose of illustrating the general principles of the application, and is not meant to limit the scope of the application, but to limit the application to the particular embodiments, and any modifications, equivalents, improvements, etc. that fall within the spirit and principles of the application are intended to be included within the scope of the application.

Claims (5)

1. A method of interphase drag analysis for a bundle channel, the method comprising:
judging the flow pattern of the rod bundle channel according to the rod bundle channel, and obtaining a flow pattern judgment result;
according to the flow pattern judgment result, calculating the phase-to-phase resistance by adopting different phase-to-phase resistance models aiming at different flow patterns; wherein: the flow pattern judgment result comprises bubble flow, cap flow and annular flow;
the flow pattern of the bar bundle channel is judged, and a flow pattern judgment result is obtained; the judgment basis is as follows:
judging the flow pattern according to the flow rate and the cavitation share of the rod bundle channel;
if the cavitation share of the rod bundle channel is smaller than 0.2, the flow pattern judgment result is bubble flow;
if the flow rate of the rod bundle channel is less than or equal to 3000 kg/(m) 2 S) and the cavitation fraction is between 0.2 and 0.85, or the flow rate of the rod bundle channel is greater than 3000 kg/(m) 2 S) and the cavitation proportion is between 0.4 and 0.85, and the judgment result of the flow pattern is cap-shaped flow;
if the cavitation share of the rod bundle channel is greater than 0.85, the flow pattern judgment result is annular flow;
calculating phase-to-phase resistance by adopting a first phase-to-phase resistance model aiming at bubble flow, wherein the first phase-to-phase resistance model is obtained by taking phase interface concentration formed by taking the stirring characteristic of a rod bundle channel into consideration;
calculating inter-phase resistance by adopting a second inter-phase resistance model aiming at cap-shaped flow;
calculating inter-phase resistance by adopting a third inter-phase resistance model aiming at the annular flow;
the first inter-phase resistance model and the second inter-phase resistance model are both in relation:
wherein F is i Is inter-phase resistance; ρ c Is the density of continuous phase, namely the density rho of liquid phase f Or gas/vapor phase density ρ g ；v r For the relative velocity between the gas/vapor phase and the liquid phase, i.e. v r ＝v g -v f ，v g And v f Gas/vapor phase and liquid phase velocities, respectively; c (C) D Is the drag coefficient; a, a gf Is the concentration of phase interface;
the phase interface concentration calculation formula of the first interphase resistance model is as follows:
wherein α is the cavitation fraction; sigma is the coefficient of viscosity of the adhesive,calculated by a Chexal-Lellouche relational expression; c (C) j1 Is the first key parameter obtained by considering the influence of the whipping of a rod bundle channel on the bubble characteristics, C j2 Is a second key parameter that is derived in view of the impact of the rod-channel whipping on the bubble characteristics.
2. The method for analyzing interphase resistance applicable to a bundle channel according to claim 1, wherein the second interphase resistance model is used for correcting the bubble-like fluid phase interface concentration according to the cap-like fluid flow pattern interface concentration closure by the ratio of the void fraction at the transition of the flow pattern to the actual void fraction;
the phase interface concentration calculation formula of the second phase-to-phase resistance model is as follows:
wherein alpha is bc A cavitation fraction that is converted from bubble flow to cap flow;
the drag coefficient calculation formula of the second inter-phase resistance model is as follows:
wherein C is D For drag coefficient, alpha ca Cavitation fraction for hat flow to annular flow transition, Δρ is the liquid phase density ρ f Density ρ of gas/vapor phase g And (3) a difference.
3. The method for analyzing interphase resistance applicable to a bundle channel according to claim 1, wherein the third interphase resistance model comprises an air/steam core interphase resistance submodel and a liquid film interphase resistance submodel;
the expression of the third inter-phase resistance model is:
F i ＝F i_c +F i_d
wherein F is i_c Is the inter-phase resistance of the gas/steam core, F i_d Is the inter-phase resistance of the liquid film, F i_c And F i_d The same relation calculation as the first interphase resistance model is adopted, and the corresponding drag coefficient and interface concentration are corrected.
4. A method of analyzing interphase resistance for a bundle tunnel according to claim 3, wherein said modifying the corresponding drag coefficient and phase interface concentration comprises:
for the gas/steam core, the drag coefficient C of the gas/steam core D_c Phase interface concentration a of gas/steam core gf_c The calculation formula is as follows:
1) When the vapor phase Reynolds number Re g When the number of the times is less than 500,
2) When the vapor phase Reynolds number Re g When the number of the carbon atoms is more than 1500,
C D_c ＝0.02{1+150[1-(1-α ff ) 0.5 ]} (7)
wherein alpha is f As liquid film/liquid phase fraction, C f Is a constant related to the structure of the bundle channels; v g For gas/vapor phase velocity, alpha ff Is the entrainment share;
3) Vapor phase Reynolds number Re g When the drag coefficient is between 500 and 1500, the calculation formula of the drag coefficient of the gas/steam core adopts formula (6) and formula #7) A difference between them;
phase interface concentration a of gas/steam core gf_c The calculation formula is as follows:
wherein C is j3 To take into account the constant of the influence of the rod-channel stirring on the phase interface.
5. An interphase resistance analysis device suitable for a bundle channel, characterized in that the device supports an interphase resistance analysis method suitable for a bundle channel according to any one of claims 1 to 4; the device comprises:
the flow pattern judging unit is used for judging the flow pattern of the bar bundle channel according to the bar bundle channel and obtaining a flow pattern judging result;
the inter-phase resistance calculation unit is used for calculating the inter-phase resistance by adopting different inter-phase resistance models according to the flow pattern judgment result; wherein: the flow pattern judgment result includes bubble flow, cap flow and annular flow.
CN202210873680.2A 2022-07-21 2022-07-21 Interphase resistance analysis method and device suitable for rod bundle channel Active CN115130323B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210873680.2A CN115130323B (en) 2022-07-21 2022-07-21 Interphase resistance analysis method and device suitable for rod bundle channel

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210873680.2A CN115130323B (en) 2022-07-21 2022-07-21 Interphase resistance analysis method and device suitable for rod bundle channel

Publications (2)

Publication Number Publication Date
CN115130323A CN115130323A (en) 2022-09-30
CN115130323B true CN115130323B (en) 2023-10-13

Family

ID=83384609

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210873680.2A Active CN115130323B (en) 2022-07-21 2022-07-21 Interphase resistance analysis method and device suitable for rod bundle channel

Country Status (1)

Country Link
CN (1) CN115130323B (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106897478A (en) * 2016-12-19 2017-06-27 杨保文 The subchannel analysis method of effect is mixed comprising screen work based on resistance distribution
CN107895095A (en) * 2017-12-15 2018-04-10 西安交通大学 The coupling analytical method of cluster subchannel and critical heat flux density mechanism model
CN113688580A (en) * 2021-07-21 2021-11-23 上海交通大学 Gas-liquid two-phase flow interface density calculation method, device, equipment and storage medium
CN114492244A (en) * 2022-02-10 2022-05-13 上海交通大学 Method for rapidly calculating distribution of void fraction of gas-liquid two-phase flow in rod bundle sub-channel

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106897478A (en) * 2016-12-19 2017-06-27 杨保文 The subchannel analysis method of effect is mixed comprising screen work based on resistance distribution
CN107895095A (en) * 2017-12-15 2018-04-10 西安交通大学 The coupling analytical method of cluster subchannel and critical heat flux density mechanism model
CN113688580A (en) * 2021-07-21 2021-11-23 上海交通大学 Gas-liquid two-phase flow interface density calculation method, device, equipment and storage medium
CN114492244A (en) * 2022-02-10 2022-05-13 上海交通大学 Method for rapidly calculating distribution of void fraction of gas-liquid two-phase flow in rod bundle sub-channel

Also Published As

Publication number Publication date
CN115130323A (en) 2022-09-30

Similar Documents

Publication Publication Date Title
Xu et al. Numerical investigation on thermal hydraulic characteristics of steam jet condensation in subcooled water flow in pipes
CN106897478A (en) The subchannel analysis method of effect is mixed comprising screen work based on resistance distribution
CN113128096B (en) Method for acquiring direct navigation additional mass of underwater vehicle
CN111985166A (en) Pipeline hydraulic transient simulation method and storage medium with implicit consideration of dynamic friction resistance
CN103488862A (en) Method for simulating high-viscosity fluid
Reddick et al. Effect of a mixture of carbon dioxide and steam on ejector performance: An experimental parametric investigation
CN115130323B (en) Interphase resistance analysis method and device suitable for rod bundle channel
CN111125972B (en) Hydraulic load analysis method for water loss accident of break of nuclear power plant
Benim et al. RANS predictions of turbulent flow past a circular cylinder over the critical regime
CN113688580A (en) Gas-liquid two-phase flow interface density calculation method, device, equipment and storage medium
CN112613247A (en) Flow field analysis method and system for filling high-porosity porous medium open channel
Vallabh et al. Investigation of nozzle contours in the CSIR supersonic wind tunnel
CN115221812B (en) Interphase resistance analysis method and device suitable for rectangular channel
Xu et al. Wet gas overreading characteristics of a long-throat Venturi at high pressure based on CFD
CN104573365B (en) A kind of flow conditioner performance estimating method
CN114943054B (en) Method for determining morphology and physical field of multiphase sintering by flow field coupling phase field
CN114169203B (en) Two-phase flow total hidden numerical method for transient safety analysis of nuclear power plant
McDonald et al. A survey of recent work on interacted boundary-layer theory for flow with separation
CN106017582A (en) A pitot tube flow measuring method based on the tube diameter dichotomy principle
JP2019164610A (en) Simulation method, simulation device, and program
Yuan et al. Numerical analyses of aerodynamic characteristics of integrated L-shaped high-rise building
Venter et al. Benchmarking the Performance of the Actuator-Disk Method for Low-Pressure Axial Flow Fan Simulation
CN113505544A (en) Bicycle motion virtual numerical value wind tunnel system based on finite volume method
Narsipur Low-order modeling of dynamic stall on airfoils in incompressible flow
CN113627099B (en) Method and device for predicting cavitation flow field

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