CN115221812A - Interphase resistance analysis method and device suitable for rectangular channel - Google Patents
Interphase resistance analysis method and device suitable for rectangular channel Download PDFInfo
- Publication number
- CN115221812A CN115221812A CN202210859528.9A CN202210859528A CN115221812A CN 115221812 A CN115221812 A CN 115221812A CN 202210859528 A CN202210859528 A CN 202210859528A CN 115221812 A CN115221812 A CN 115221812A
- Authority
- CN
- China
- Prior art keywords
- interphase resistance
- rectangular channel
- resistance
- interphase
- flow
- 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/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
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses an interphase resistance analysis method and device suitable for a rectangular channel, which comprises the following steps: judging the flow pattern of the narrow rectangular channel according to the narrow rectangular channel to obtain a flow pattern judgment result; according to the flow pattern judgment result, calculating interphase resistance by adopting different interphase resistance models aiming at different flow patterns; if the flow type judgment result is the bubble flow, calculating the interphase resistance by adopting a first interphase resistance model; if the flow pattern judgment result is elastic flow, calculating the interphase resistance by adopting a second interphase resistance model; and if the flow pattern judgment result is annular flow, calculating the interphase resistance by adopting a third interphase resistance model. The method is suitable for calculating and analyzing the interphase resistance of the bubbly flow, the elastic flow and the annular flow in the narrow rectangular channel, adopts a basic form of a drag force model, and considers the structural influence of the rectangular channel on the phase interface concentration and the drag force coefficient; the prediction accuracy of phenomena such as the two-phase flow characteristic of the thermal hydraulic characteristic and the void fraction distribution can be improved.
Description
Technical Field
The invention relates to the field of reactor thermohydraulic design and safety analysis, in particular to an interphase resistance analysis method and device suitable for a rectangular channel.
Background
When gas/vapor-liquid two phases are mixed and flow together, the flow characteristics of the gas/vapor-liquid two phases are completely different due to the difference of the physical properties of the two phases, so that the gas/vapor-liquid two phases have interaction force, and the force is interphase resistance. The interphase resistance has important influences on the deformation of a two-phase interface, the phase state distribution, the flow pressure drop and the like, and further influences the thermal hydraulic characteristics of the reactor. Accurate simulation of interphase resistance characteristic-related phenomena requires an accurate interphase resistance model. As an important basis for selecting a basic equation and a model in the two-phase flow and heat transfer calculation process, the interphase resistance model is a basic composition model of the current mainstream thermotechnical hydraulic program (such as RELAP5, CATARE, COBRA-TF and the like).
At present, an interphase resistance analysis method adopted by a mainstream system analysis program is mainly developed based on round tube or rod bundle experimental data and is not suitable for a narrow rectangular channel. The reason is as follows: for a narrow rectangular channel, the growth of bubbles is limited by the width of the rectangular channel, and circular bubbles which can be formed in a circular tube or a rod bundle can be extruded and deformed in the rectangular channel, so that the gas/vapor-liquid interaction area is changed, and the method for analyzing the interphase resistance in the circular tube or the rod bundle channel is not suitable for the rectangular channel.
Disclosure of Invention
The invention aims to solve the technical problem that the existing interphase resistance analysis method is mainly developed based on round tube or rod bundle experimental data and is not suitable for narrow rectangular channels. This is because, for a narrow rectangular channel, the growth of bubbles is limited by the width of the rectangular channel, and circular bubbles that can be formed in a round tube or a rod bundle are extruded and deformed in the rectangular channel, resulting in a change in the gas/vapor-liquid interaction area, and thus causing the method for analyzing the interphase resistance in the round tube or the rod bundle channel to be unsuitable for the rectangular channel.
The invention aims to provide an interphase resistance analysis method and an interphase resistance analysis device suitable for a rectangular channel, wherein the interphase resistance calculation and analysis method suitable for bubble flow, elastic flow and annular flow in a narrow rectangular channel is designed, a basic form of a drag force model is adopted, and the influence of the rectangular channel structure is considered in the aspects of interface concentration and drag force coefficient; the prediction accuracy of phenomena such as the two-phase flow characteristic of the thermal hydraulic characteristic and the void fraction distribution can be improved.
The invention is realized by the following technical scheme:
in a first aspect, the present invention provides an interphase resistance analysis method suitable for a rectangular channel, including:
judging the flow pattern of the narrow rectangular channel according to the narrow rectangular channel to obtain a flow pattern judgment result;
according to the flow pattern judgment result, calculating interphase resistance by adopting different interphase resistance models aiming at different flow patterns; wherein:
if the flow pattern judgment result is a bubble flow, calculating the interphase resistance by adopting a first interphase resistance model to obtain a first interphase resistance result; if the flow pattern judgment result is elastic flow, calculating the interphase resistance by adopting a second interphase resistance model to obtain a second interphase resistance result; and if the flow pattern judgment result is annular flow, calculating the inter-phase resistance by adopting a third inter-phase resistance model to obtain a third inter-phase resistance result.
The working principle is as follows: the existing interphase resistance analysis method is mainly developed based on round tube or rod bundle experimental data and is not suitable for narrow rectangular channels. The reason is that for the narrow rectangular channel, the growth of bubbles is limited by the width of the rectangular channel, and the circular bubbles which can be formed in the circular tube or the rod bundle can be extruded and deformed in the rectangular channel, so that the gas/vapor-liquid interaction area is changed, and the interphase resistance analysis method in the circular tube or the rod bundle channel is not suitable for the rectangular channel.
On one hand, the method considers that no formed interphase resistance analysis for the annular flow exists in the narrow rectangular channel in the prior art, if the annular flow is encountered in the narrow rectangular channel, the interphase resistance analysis for the annular flow in the narrow rectangular channel is also analyzed by adopting an interphase resistance analysis method for a circular tube or a rod bundle, and the method is not suitable for the annular flow in fact, so that the interphase resistance calculation and analysis for the annular flow are not accurate; the invention designs formed interphase resistance analysis aiming at annular flow in a narrow rectangular channel. On the other hand, the method is mainly developed based on a drag model, and drag coefficients and phase interface concentrations in the model are mainly closed by a relation of semi-experience-semi-theory; the interphase resistance analysis of the narrow rectangular channel is the same as the interphase resistance analysis in the circular tube and rod bundle channels in form, and the key point is that the structural influence of the geometric channel is considered on the drag coefficient and the phase interface concentration, so that the interphase resistance calculation and analysis of the rectangular channel without flow patterns are more accurate, and the core of the invention is also the interphase resistance analysis of the rectangular channel without flow patterns.
Further, the flow pattern of the narrow rectangular channel is judged according to the following criteria: and judging the flow pattern according to the flow and the void fraction of the narrow rectangular channel.
Further, the first inter-phase resistance model is an inter-phase resistance model obtained by considering a phase interface concentration of a rectangular channel characteristic;
the relation of the first interphase resistance model is as follows:
wherein, F i Is interphase resistance; rho c Is continuous phase density, i.e. liquid phase density rho f Or gas/vapor phase density ρ g ；v r Is the relative velocity between the gas/vapour 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 D For drag coefficient, ishii is used&Calculating a relation of Chawla (1979); a is a gf Is the phase interface concentration.
Further, the calculation formula of the phase interface concentration is as follows:
wherein alpha is the void fraction in the narrow rectangular channel, sigma is the viscosity coefficient,
the vector is obtained by the calculation of a Chexal-Lellouche relational expression; c j Is a key parameter that takes into account the effect of the rectangular channel structure size on the bubble characteristics.
Further, the second interphase resistance model comprises an interphase resistance submodel of the gas/steam bullet and the mainstream liquid and an interphase resistance submodel of the wake flow liquid block and the mainstream liquid;
the second interphase resistance model has the relation:
F i ＝F i_t +F i_s (3)
wherein, F i_t Is the interphase resistance of the gas/vapour bomb and the main flow of liquid, F i_s Is the interphase resistance of the wake flow liquid block and the main flow liquid; f i_t And F i_s Calculating by adopting the same relational expression as the first interphase resistance model;
meanwhile, for wake flow liquid block, the liquid block part can be approximately regarded as bubble flow, and the same calculation relation F as the bubble flow is adopted i_s Only the void fraction needs to be corrected into the liquid block part void fraction;
for gas/steam bullet, interphase resistance F of gas/steam bullet and main flow liquid i_t The drag coefficient and the phase interface concentration are calculated by the length of the narrow rectangular channel and the width of the narrow rectangular channel;
the calculation formula of the drag coefficient is as follows:
wherein, C D_t Is the drag coefficient of the gas/steam bomb, a is the length of the narrow rectangular channel, and b is the width of the narrow rectangular channel;
the calculation formula of the phase interface concentration is as follows:
wherein, a gf_t Is the phase interface concentration, alpha, of the gas/vapour bomb b Is the void fraction of a single bubble.
Further, the third interphase resistance model comprises a gas/vapor core interphase resistance submodel and a liquid drop interphase resistance submodel;
the expression of the third inter-phase resistance model is as follows:
F i ＝F i_c +F i_d (6)
wherein, F i_c Is the interphase resistance of the gas/steam core, F i_d Is the interphase resistance of the droplets, F i_c And F i_d And calculating by using the same relational expression as the first interphase resistance model, and correcting the corresponding drag coefficient and phase interface concentration.
Further, the correcting the corresponding drag coefficient and the phase interface concentration includes:
for gas/steam core, drag coefficient C of gas/steam core D_c Phase interface concentration alpha of gas/steam core ff The calculation formula is as follows:
reynolds number Re when vapor phase g When the average particle size is less than 500 a,
reynolds number Re when vapor phase g When greater than 1500, C D_c ＝0.02{1+150[1-(1-α ff ) 0.5 ]} (8)
Wherein alpha is f In terms of droplet/liquid phase fraction, C f Is a constant related to the size of the rectangular channel structure; alpha (alpha) ("alpha") f Is a void fraction, v g Is the vapor phase velocity, α ff Is the entrainment fraction;
reynolds number Re of vapor phase g In the range of 500 to 1500, the calculation formula of the drag coefficient of the gas/steam core adopts the difference between the formula (7) and the formula (8).
In a second aspect, the invention further provides an interphase resistance analysis device suitable for a rectangular channel, which supports the interphase resistance analysis method suitable for the rectangular channel; the device includes:
the flow pattern judging unit is used for judging the flow pattern of the narrow rectangular channel according to the narrow rectangular channel to obtain a flow pattern judging result;
the interphase resistance calculation unit is used for calculating interphase resistance by adopting different interphase resistance models aiming at different flow patterns according to the flow pattern judgment result; the interphase resistance calculating unit comprises a first interphase resistance calculating subunit, a second interphase resistance calculating subunit and a third interphase resistance calculating subunit;
the first interphase resistance calculating subunit is used for calculating the interphase resistance of the bubble-shaped flow by adopting a first interphase resistance model to obtain a first interphase resistance result;
the second interphase resistance calculation subunit is used for calculating the interphase resistance of the elastohydrodynamic flow by adopting a second interphase resistance model to obtain a second interphase resistance result;
and the third inter-phase resistance calculation subunit is used for calculating inter-phase resistance of the annular flow by adopting a third inter-phase resistance model to obtain a third inter-phase resistance result.
Further, the flow pattern of the narrow rectangular channel is judged according to the following conditions: and judging the flow pattern according to the flow and the void fraction of the narrow rectangular channel.
Further, the first interphase resistance model of the first interphase resistance calculating subunit is an interphase resistance model obtained by considering a phase interface concentration of a rectangular channel characteristic;
the relation of the first interphase resistance model is as follows:
wherein, F i Is interphase resistance; rho c As density of continuous phase, i.e. density of liquid phase, p f Or gas/vapor phase density ρ g ；v r Is the relative velocity between the gas/vapour 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 D For drag coefficient, ishii is used&Calculating a relation formula of Chawla (1979); a is gf Is the phase interface concentration;
the calculation formula of the phase interface concentration is as follows:
wherein alpha is the void fraction in the narrow rectangular channel, sigma is the viscosity coefficient,
the vector is obtained by the calculation of a Chexal-Lellouche relational expression; c j Is a key parameter that takes into account the effect of the rectangular channel structure size on the bubble characteristics.
Compared with the prior art, the invention has the following advantages and beneficial effects:
1. the invention designs a set of interphase resistance calculation and analysis method suitable for bubble flow, elastic flow and annular flow in a narrow rectangular channel, adopts a basic form of a drag force model, and considers the structural influence of the rectangular channel on interface concentration and drag force coefficients; the calculation of the interphase resistance under the bubble flow, the bullet flow and the annular flow in the narrow rectangular channel is more accurate.
2. The method adopts a basic form of a drag force model, and considers the structural influence of a rectangular channel on the interface concentration and drag force coefficient; according to the structural form of the rectangular channel, the phase interface concentration under the bubble flow is corrected, the drag coefficient and the phase interface concentration under the elastic flow are corrected, the phase interface concentration under the annular flow is corrected, an interphase resistance analysis method suitable for the rectangular channel is formed, and basic model support can be provided for the research and development of a thermal hydraulic program based on a two-fluid-six equation.
3. Compared with the calculation and analysis of the inter-phase resistance in the circular tube and the rod bundle channel, the difference of the calculation and analysis of the inter-phase resistance of the rectangular channel is relatively small, and the program integration is facilitated. The analysis method can improve the prediction accuracy of phenomena such as the two-phase flow characteristic of the thermal hydraulic characteristic and the void fraction distribution.
Drawings
The accompanying drawings, which are included to provide a further understanding of the embodiments of the invention and are incorporated in and constitute a part of this application, illustrate embodiment(s) of the invention and together with the description serve to explain the principles of the invention. In the drawings:
FIG. 1 is a schematic diagram of the invention in which the bullet flow in the narrow rectangular channel mainly comprises two parts of an air/steam bullet and a wake flow liquid block.
Fig. 2 is a flowchart of an interphase resistance analysis method suitable for a rectangular channel according to an embodiment of the present invention.
Fig. 3 is a detailed flowchart of an interphase resistance analysis method suitable for a rectangular channel according to an embodiment of the present invention.
Fig. 4 is a schematic structural diagram of an interphase resistance analysis device suitable for a rectangular channel according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail below with reference to examples and accompanying drawings, and the exemplary embodiments and descriptions thereof are only used for explaining the present invention and are not meant to limit the present invention.
Example 1
The existing interphase resistance analysis method is mainly developed based on round tube or rod bundle experimental data and is not suitable for narrow rectangular channels. This is because, for a narrow rectangular channel, the growth of bubbles is limited by the width of the rectangular channel, and circular bubbles that can be formed in a round tube or a rod bundle are extruded and deformed in the rectangular channel, resulting in a change in the gas/vapor-liquid interaction area, and thus causing the method for analyzing the interphase resistance in the round tube or the rod bundle channel to be unsuitable for the rectangular channel.
On one hand, the method considers that no formed interphase resistance analysis for the annular flow exists in the narrow rectangular channel in the prior art, if the annular flow is encountered in the narrow rectangular channel, the interphase resistance analysis for the annular flow in the narrow rectangular channel is also analyzed by adopting an interphase resistance analysis method for a circular tube or a rod bundle, and the method is not suitable for the annular flow in fact, so that the interphase resistance calculation and analysis for the annular flow are not accurate; the invention designs a formed interphase resistance analysis for annular flow in a narrow rectangular channel. On the other hand, the method is mainly developed based on a drag force model, and drag force coefficients and phase interface concentration in the model are mainly closed by a relation of semi-experience-semi-theory; the interphase resistance analysis of the narrow rectangular channel is the same as the interphase resistance analysis in the circular tube and rod bundle channels in form, and the important point is that the structural influence of the geometric channel is considered on the drag coefficient and the phase interface concentration, so that the interphase resistance calculation and analysis of the rectangular channel without the flow pattern is more accurate, and the core of the invention is also the interphase resistance analysis of the rectangular channel without the flow pattern.
As shown in fig. 2 and 3, the present invention provides an interphase resistance analysis method suitable for a rectangular channel, including:
judging the flow pattern of the narrow rectangular channel according to the narrow rectangular channel to obtain a flow pattern judgment result;
according to the flow pattern judgment result, calculating interphase resistance by adopting different interphase resistance models aiming at different flow patterns; wherein:
if the flow pattern judgment result is a bubble flow, calculating the interphase resistance by adopting a first interphase resistance model to obtain a first interphase resistance result; if the flow pattern judgment result is elastic flow, calculating the interphase resistance by adopting a second interphase resistance model to obtain a second interphase resistance result; and if the flow pattern judgment result is annular flow, calculating the inter-phase resistance by adopting a third inter-phase resistance model to obtain a third inter-phase resistance result.
By further implementation, the flow pattern of the narrow rectangular channel is judged according to the following steps: and judging the flow pattern according to the flow and the void fraction of the narrow rectangular channel.
As a further implementation, the first inter-phase resistance model is an inter-phase resistance model obtained by taking into account a phase interface concentration of a rectangular channel characteristic;
the relation of the first interphase resistance model is as follows:
wherein, F i Is interphase resistance; rho c Is continuous phase density, i.e. liquid phase density rho f Or gas/vapor phase density ρ g ；v r Is a phase between gas/vapor phase and liquid phaseTo velocity, i.e. v r ＝v g -v f ，v g And v f Gas/vapor phase and liquid phase velocities, respectively; c D For drag coefficient, ishii is used&Calculating a relation of Chawla (1979); a is gf The phase interface concentration is calculated by taking the rectangular channel characteristics into consideration and based on the bubble chord length distribution characteristics through the following calculation formula:
wherein alpha is the void fraction in the narrow rectangular channel, sigma is the viscosity coefficient,
the method is obtained by the calculation of a Chexal-Lellouche relational expression; c j Is a key parameter that takes into account the effect of the rectangular channel structure size on the bubble characteristics.
As a further implementation, the bullet flow in the narrow rectangular channel mainly comprises two parts of an air/steam bullet and a wake flow liquid block, as shown in FIG. 1. The interphase resistance under the flow pattern is mainly divided into two parts: the interphase resistance of the gas/steam bomb and the main flow liquid and the interphase resistance of the wake flow liquid block and the main flow liquid;
specifically, the second interphase resistance model comprises an interphase resistance submodel of the gas/steam bullet and the mainstream liquid and an interphase resistance submodel of the wake flow liquid block and the mainstream liquid;
the second interphase resistance model has the relation:
F i ＝F i_t +F i_s (3)
wherein, F i_t Is the interphase resistance of the gas/steam bomb and the main flow liquid, F i_s The inter-phase resistance of the wake flow liquid block and the main flow liquid; f i_t And F i_s All calculated by using the same relational expression (formula (1)) as the first interphase resistance model;
meanwhile, for wake flow liquid block, the liquid block part can be approximately regarded as bubble flow, and the same calculation relation F as the bubble flow is adopted i_s Only the vacuole portion is requiredCorrecting the blank liquid to be a vacuole portion of the liquid block; i.e. drag coefficient C of wake flow liquid mass D_s Using Ishii&Calculating the phase interface concentration a of the wake fluid block by using a relation of Chawla (1979) gf_s Calculated by the formula (2).
For the gas/steam bullet, the rectangular channel influences the shape of the gas/steam bullet. Drag coefficient C D_t No longer using Ishii&Chawla (1979) relational computation; but the interphase resistance F of the gas/steam bullet and the main flow liquid i_t The drag coefficient and the phase interface concentration are calculated by the length of the narrow rectangular channel and the width of the narrow rectangular channel;
the calculation formula of the drag coefficient is as follows:
wherein, C D_t The drag coefficient of the gas/steam bullet, wherein a is the length of the narrow rectangular channel, and b is the width of the narrow rectangular channel;
the calculation formula of the phase interface concentration is as follows:
wherein, a gf_t The phase interface concentration of the gas/steam bomb; alpha (alpha) ("alpha") b Is the void fraction of a single bubble, and is a basic parameter.
As a further implementation, the annular flow in the rectangular channel mainly comprises two parts of a gas/vapor core and liquid drops. The interphase resistance under the flow pattern is mainly divided into two parts: interphase resistance of gas/vapor core and liquid drop.
The third interphase resistance model comprises a gas/steam core interphase resistance submodel and a liquid drop interphase resistance submodel;
the expression of the third inter-phase resistance model is as follows:
F i ＝F i_c +F i_d (6)
wherein, F i_c Is the interphase resistance of the gas/steam core, F i_d Is the interphase resistance of the droplets, F i_c And F i_d And calculating by using the same relational expression as the first interphase resistance model, and correcting the corresponding drag coefficient and phase interface concentration.
Further, the correcting the corresponding drag coefficient and the phase interface concentration includes:
for gas/steam core, drag coefficient C of gas/steam core D_c Phase interface concentration alpha of gas/vapor core ff The calculation formula is as follows:
when vapor phase Reynolds number Re g When the average particle size is less than 500 a,
reynolds number Re when vapor phase g When greater than 1500, C D_c ＝0.02{1+150[1-(1-α ff ) 0.5 ]} (8)
Wherein alpha is f In droplet/liquid phase fraction, C f Is a constant related to the size of the rectangular channel structure; alpha (alpha) ("alpha") f Is a void fraction, v g Is the vapor phase velocity, α ff Is the entrainment fraction;
reynolds number Re of vapor phase 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 (7) and the formula (8).
For droplets, drag coefficient C D_d Using Ishii&Calculating a relation of Chawla (1979); concentration of phase interface a gf_d The calculation relation is the same as the formula (2), and only alpha in the formula (2) needs to be modified into the liquid drop share alpha f 。
The method can be used for calculating the gas/vapor-liquid interphase acting force in the rectangular channel and can be implanted into a thermal hydraulic analysis program based on two-fluid-six equation. Compared with the interphase resistance analysis in a mainstream system analysis program RELAP5, the method provided by the invention has the following characteristics:
1) The interphase resistance analysis method provided by the invention is suitable for the rectangular narrow channel, and provides the formed interphase resistance analysis for the annular flow in the narrow rectangular channel.
2) The method adopts a basic form of a drag force model, and considers the influence of a rectangular channel structure on the phase interface concentration and the drag force coefficient; according to the structural form of the rectangular channel, the phase interface concentration under the bubble flow is corrected, the drag coefficient and the phase interface concentration under the elastic flow are corrected, the phase interface concentration under the annular flow is corrected, an interphase resistance analysis method suitable for the rectangular channel is formed, and basic model support can be provided for the research and development of a thermal hydraulic program based on a two-fluid-six equation.
3) Compared with the calculation and analysis of the inter-phase resistance in the circular tube and the rod bundle channel, the difference of the calculation and analysis of the inter-phase resistance of the rectangular channel is relatively small, and the program integration is facilitated. The analysis method can improve the prediction accuracy of phenomena such as two-phase flow characteristics of thermal hydraulic characteristics, vacuole share distribution and the like.
Example 2
As shown in fig. 4, this embodiment is different from embodiment 1 in that this embodiment provides an interphase resistance analyzing apparatus suitable for a rectangular channel, which supports an interphase resistance analyzing method suitable for a rectangular channel described in embodiment 1; the device includes:
the flow pattern judging unit is used for judging the flow pattern of the narrow rectangular channel according to the narrow rectangular channel to obtain a flow pattern judging result;
the interphase resistance calculation unit is used for calculating interphase resistance by adopting different interphase resistance models according to different flow patterns according to the flow pattern judgment result; the interphase resistance calculating unit comprises a first interphase resistance calculating subunit, a second interphase resistance calculating subunit and a third interphase resistance calculating subunit;
the first interphase resistance calculating subunit is used for calculating the interphase resistance of the bubble flow by adopting a first interphase resistance model to obtain a first interphase resistance result;
the second interphase resistance calculation subunit is used for calculating the interphase resistance of the elastohydrodynamic flow by adopting a second interphase resistance model to obtain a second interphase resistance result;
and the third inter-phase resistance calculation subunit is used for calculating inter-phase resistance of the annular flow by adopting a third inter-phase resistance model to obtain a third inter-phase resistance result.
In this embodiment, the basis for judging the flow pattern of the narrow rectangular channel is as follows: and judging the flow pattern according to the flow and the void fraction of the narrow rectangular channel.
In the present embodiment, the first interphase resistance model of the first interphase resistance calculating subunit is an interphase resistance model obtained by considering a phase interface concentration of a rectangular channel characteristic;
the relation of the first interphase resistance model is as follows:
wherein, F i Is interphase resistance; rho c Is continuous phase density, i.e. liquid phase density rho f Or gas/vapor phase density ρ g ；v r Is the relative velocity between the gas/vapour 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 D For drag coefficient, ishii is used&Calculating a relation formula of Chawla (1979); a is a gf Is the phase interface concentration;
the calculation formula of the phase interface concentration is as follows:
wherein alpha is the void fraction in the narrow rectangular channel, sigma is the viscosity coefficient,
the method is obtained by the calculation of a Chexal-Lellouche relational expression; c j Is a key parameter that takes into account the effect of the rectangular channel structure size on the bubble characteristics.
The execution process of each unit is executed according to the flow steps of the interphase resistance analysis method applicable to the rectangular channel described in embodiment 1, and details are not repeated in this embodiment.
As will be appreciated by one skilled in the art, 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 so forth) 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 flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams 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 above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are merely exemplary embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (10)
1. An interphase resistance analysis method suitable for a rectangular channel, the method comprising:
acquiring a flow pattern judgment result of the narrow rectangular channel, and calculating interphase resistance by adopting different interphase resistance models for different flow patterns according to the flow pattern judgment result; wherein:
if the flow pattern judgment result is a bubble flow, calculating the interphase resistance by adopting a first interphase resistance model to obtain a first interphase resistance result; if the flow pattern judgment result is elastic flow, calculating the interphase resistance by adopting a second interphase resistance model to obtain a second interphase resistance result; and if the flow pattern judgment result is annular flow, calculating the inter-phase resistance by adopting a third inter-phase resistance model to obtain a third inter-phase resistance result.
2. The method for analyzing the interphase resistance in the rectangular channel according to claim 1, further comprising: judging the flow pattern of the narrow rectangular channel according to the narrow rectangular channel to obtain a flow pattern judgment result of the narrow rectangular channel;
the flow pattern of the narrow rectangular channel is judged according to the following steps: and judging the flow pattern according to the flow and the void fraction of the narrow rectangular channel.
3. The interphase resistance analysis method according to claim 1, wherein the first interphase resistance model is an interphase resistance model obtained by considering a phase interface concentration of a rectangular channel characteristic;
the relation of the first interphase resistance model is as follows:
wherein, F i Is interphase resistance; rho c Is continuous phase density, i.e. liquid phase density rho f Or gas/vapor phase density ρ g ；v r Is the relative velocity between the gas/vapour 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 D Is the drag coefficient; a is gf Is the phase interface concentration.
4. The method for analyzing interphase resistance according to claim 3, wherein the formula for the phase interface concentration is:
wherein alpha is the void fraction in the narrow rectangular channel, sigma is the viscosity coefficient,
the vector is obtained by the calculation of a Chexal-Lellouche relational expression; c j Is a key parameter that takes into account the effect of the rectangular channel structure size on the bubble characteristics.
5. The method for analyzing the interphase resistance applicable to the rectangular channel according to claim 3, wherein the second interphase resistance model comprises an interphase resistance submodel of a gas/steam bullet and a mainstream liquid and an interphase resistance submodel of a wake liquid block and the mainstream liquid;
the second interphase resistance model has the relation:
F i ＝F i_t +F i_s
wherein, F i_t Is the interphase resistance of the gas/steam bomb and the main flow liquid, F i_s Is the interphase resistance of the wake flow liquid block and the main flow liquid; f i_t And F i_s Calculating by adopting the same relational expression as the first interphase resistance model;
meanwhile, for wake flow liquid block, the same calculation relation formula F is adopted for the liquid block part as the bubble flow i_s Only the void fraction needs to be corrected into the liquid block part void fraction;
for gas/steam bullet, interphase resistance F of gas/steam bullet and main flow liquid i_t The drag coefficient and the phase interface concentration are calculated by the length of the narrow rectangular channel and the width of the narrow rectangular channel;
the calculation formula of the drag coefficient is as follows:
wherein, C D_t Is the drag coefficient of the gas/steam bomb, a is the length of the narrow rectangular channel, and b is the width of the narrow rectangular channel;
the calculation formula of the phase interface concentration is as follows:
wherein, a gf_t Is the phase interface concentration of gas/vapour bomb, alpha b Is the void fraction of a single bubble.
6. The method for analyzing the interphase resistance applicable to the rectangular channel according to claim 3, wherein the third interphase resistance model includes a gas/vapor core interphase resistance submodel and a liquid droplet interphase resistance submodel;
the expression of the third inter-phase resistance model is as follows:
F i ＝F i_c +F i_d
wherein, F i_c Is the interphase resistance of the gas/steam core, F i_d Is the interphase resistance of the droplets, F i_c And F i_d And calculating by using the same relational expression as the first interphase resistance model, and correcting the corresponding drag coefficient and phase interface concentration.
7. The method for analyzing the interphase resistance in the rectangular channel according to claim 6, wherein the correcting the corresponding drag coefficient and the phase interface concentration includes:
for gas/steam core, drag coefficient C of gas/steam core D_c Phase interface concentration alpha of gas/steam core ff The calculation formula is as follows:
reynolds number Re when vapor phase g When the average particle size is less than 500 a,
reynolds number Re when vapor phase g When greater than 1500, C D_c ＝0.02{1+150[1-(1-α ff ) 0.5 ]} (8)
Wherein alpha is f In droplet/liquid phase fraction, C f Is a constant related to the size of the rectangular channel structure; alpha is alpha f Is a void fraction, v g Is the vapor phase velocity, α ff Is the entrainment fraction;
reynolds number Re of vapor phase g In the range of 500 to 1500, the calculation formula of the drag coefficient of the gas/steam core adopts the difference between the formula (7) and the formula (8).
8. An interphase resistance analyzing apparatus adapted to a rectangular channel, characterized by supporting an interphase resistance analyzing method adapted to a rectangular channel according to any one of claims 1 to 7; the device comprises:
the flow pattern judging unit is used for judging the flow pattern of the narrow rectangular channel according to the narrow rectangular channel to obtain a flow pattern judging result;
the interphase resistance calculation unit is used for calculating interphase resistance by adopting different interphase resistance models aiming at different flow patterns according to the flow pattern judgment result; the interphase resistance calculating unit comprises a first interphase resistance calculating subunit, a second interphase resistance calculating subunit and a third interphase resistance calculating subunit;
the first interphase resistance calculating subunit is used for calculating the interphase resistance of the bubble flow by adopting a first interphase resistance model to obtain a first interphase resistance result;
the second interphase resistance calculating subunit is used for calculating the interphase resistance of the elastic flow by adopting a second interphase resistance model to obtain a second interphase resistance result;
and the third inter-phase resistance calculation subunit is used for calculating inter-phase resistance of the annular flow by adopting a third inter-phase resistance model to obtain a third inter-phase resistance result.
9. The interphase resistance analyzing apparatus according to claim 8, wherein the flow pattern of the narrow rectangular channel is determined by: and judging the flow pattern according to the flow and the void fraction of the narrow rectangular channel.
10. The interphase resistance analyzing apparatus according to claim 8, wherein the first interphase resistance model of the first interphase resistance calculating subunit is an interphase resistance model obtained by considering a phase interface concentration of a rectangular channel characteristic;
the relation of the first interphase resistance model is as follows:
wherein, F i Is interphase resistance; rho c As density of continuous phase, i.e. density of liquid phase, p f Or gas/vapor phase density ρ g ；v r Is the relative velocity between the gas/vapour 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 D Is the drag coefficient; a is gf Is the phase interface concentration;
the calculation formula of the phase interface concentration is as follows:
wherein alpha is the void fraction in the narrow rectangular channel, sigma is the viscosity coefficient,
the method is obtained by the calculation of a Chexal-Lellouche relational expression; c j Is a key parameter that takes into account the effect of the rectangular channel structure size on the bubble characteristics.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210859528.9A CN115221812B (en) | 2022-07-21 | 2022-07-21 | Interphase resistance analysis method and device suitable for rectangular channel |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210859528.9A CN115221812B (en) | 2022-07-21 | 2022-07-21 | Interphase resistance analysis method and device suitable for rectangular channel |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115221812A true CN115221812A (en) | 2022-10-21 |
| CN115221812B CN115221812B (en) | 2023-09-19 |
Family
ID=83614311
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210859528.9A Active CN115221812B (en) | 2022-07-21 | 2022-07-21 | Interphase resistance analysis method and device suitable for rectangular channel |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115221812B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060142401A1 (en) * | 2004-11-03 | 2006-06-29 | Tonkovich Anna L | Partial boiling in mini and micro-channels |
| 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 |
-
2022
- 2022-07-21 CN CN202210859528.9A patent/CN115221812B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060142401A1 (en) * | 2004-11-03 | 2006-06-29 | Tonkovich Anna L | Partial boiling in mini and micro-channels |
| 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 |
Non-Patent Citations (2)
| Title |
|---|
| JIAN DENG等: "Experimental research and model development on interfacial drag in rectangle channel bubbly and slug flow", 《EXPERIMENTAL THERMAL AND FLUID SCIENCE 130 (2022)》, pages 1 - 13 * |
| 张卢腾等: "基于竖直圆管空气-水两相流实验的相间曳力模型研究", 《原子能科学技术（网络版）》, pages 1 - 7 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115221812B (en) | 2023-09-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Luo et al. | A mass conserving level set method for detailed numerical simulation of liquid atomization | |
| Bugg et al. | A numerical model of Taylor bubbles rising through stagnant liquids in vertical tubes | |
| CN105512363B (en) | The analogy method of water column separation in pressure conduit based on Godunov formats | |
| CN109858105B (en) | Rough microchannel internal flow field simulation analysis method based on Fluent software | |
| Aboulhasanzadeh et al. | A validation of an embedded analytical description approach for the computations of high Schmidt number mass transfer from bubbles in liquids | |
| CN113128096B (en) | Method for acquiring direct navigation additional mass of underwater vehicle | |
| Liu et al. | Natural frequency analysis of a cantilevered piping system conveying gas–liquid two-phase slug flow | |
| CN110705189A (en) | Method for establishing sedimentation air flotation tank air flotation zone hydrodynamics model | |
| WO2014069421A1 (en) | Simulation device, simulation method, and program | |
| Li et al. | Spinning disk atomization: Theory of the ligament regime | |
| Ryu et al. | A comparative study of lattice Boltzmann and volume of fluid method for two-dimensional multiphase flows | |
| CN115221812A (en) | Interphase resistance analysis method and device suitable for rectangular channel | |
| CN113688580B (en) | Gas-liquid two-phase flow phase interface density calculation method, device, equipment and storage medium | |
| CN106777770B (en) | The analogy method of hole stream in aqueduct based on finite volume method | |
| Fang | Development of Advanced Analysis Toolkit for Turbulent Bubbly Flow Simulations. | |
| CN115130323B (en) | Interphase resistance analysis method and device suitable for rod bundle channel | |
| Tong et al. | Relaxation dynamics of a free elongated liquid ligament | |
| Antritter et al. | Suppressing artificial equilibrium states caused by spurious currents in droplet spreading simulations with dynamic contact angle model | |
| Pertant et al. | A finite-volume method for simulating contact lines on unstructured meshes in a conservative level-set framework | |
| JP6728994B2 (en) | Fluid analysis device, fluid simulation method, and fluid simulation program | |
| Laín et al. | Study of flow resistance coefficients acting on regular non-spherical particles in simple shear flow at moderate Reynolds numbers | |
| CN114169203A (en) | Two-phase flow total-hidden numerical method for transient safety analysis of nuclear power plant | |
| Narsipur | Low-order modeling of dynamic stall on airfoils in incompressible flow | |
| Berdowski et al. | 3D Lagrangian VPM: simulations of the near-wake of an actuator disc and horizontal axis wind turbine | |
| Ueda et al. | Numerical Simulation of Gas-Liquid Two-Phase Flow in a Horizontally Placed Hydrophobic Rectangular Channel (Part 1, Influence of Abrupt Expansion) |
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 |