Fast flow allocation computation method, system and computer readable medium thereof

Abstract

This application provides a fast flow allocation calculation method and its applicable system and computer-readable medium, relating to the field of flow allocation technology. The fast flow allocation calculation method includes the following steps: acquiring the three-dimensional shape and size data corresponding to each pipe; dividing multiple pipes into multiple pipe groups based on the three-dimensional shape and size data; acquiring the actual total inlet mass flow rate of the pipe assembly; constructing an inlet mass flow rate-pressure difference database for each pipe group through three-dimensional CFD simulation modeling and based on the actual total inlet mass flow rate, the inlet mass flow rate-pressure difference database including multiple different single-pipe inlet mass flow rates and corresponding pressure differences; and obtaining the single-pipe inlet mass flow rate of each pipe in the pipe assembly through one-dimensional iterative calculation based on the multiple inlet mass flow rate-pressure difference databases and the actual total inlet mass flow rate, and using the predicted pressure difference corresponding to the single-pipe inlet mass flow rate as the inlet and outlet pressure difference of the pipe assembly.

Description

Technical Field

[0001] This application relates primarily to the field of traffic allocation technology, and in particular to a fast traffic allocation calculation method and its applicable system and computer-readable medium. Background Technology

[0002] In the thermal-hydraulic design of nuclear reactors, core flow distribution is a critical safety and performance design element. Core flow distribution requires flow calculations for a large number of parallel pipes of varying shapes. Because different pipe geometries result in different flow resistance coefficients, this directly affects the flow distribution results.

[0003] Currently, simplified calculation methods can only obtain the flow resistance coefficients of circular pipes and a limited number of rectangular cross-section straight pipes with fixed aspect ratios. They cannot obtain the flow resistance coefficients of other pipes with special cross-sectional shapes or bends. In addition, modeling and calculating a large number of pipes requires processing massive amounts of data, resulting in excessive computational load and time. This makes it unsuitable for the rapid iteration of core power design and fails to meet the engineering application requirements of core design and model development.

[0004] Therefore, there is an urgent need for an efficient and accurate method for calculating fast traffic allocation. Summary of the Invention

[0005] The technical problem to be solved by this application is to provide a fast traffic allocation calculation method and its applicable system and computer-readable medium, which can achieve efficient and accurate traffic allocation calculation.

[0006] To address the aforementioned technical problems, this application provides a rapid flow allocation calculation method applicable to pipe assemblies. The pipe assembly comprises multiple pipes connected in parallel, with at least two pipes having different three-dimensional shape and size data. The rapid flow allocation calculation method includes the following steps: obtaining the three-dimensional shape and size data corresponding to each pipe; dividing the multiple pipes into multiple pipe groups based on the three-dimensional shape and size data, wherein each pipe in each pipe group has the same three-dimensional shape and size data; obtaining the actual total inlet mass flow rate of the pipe assembly; constructing an inlet mass flow rate-pressure difference database for each pipe group through three-dimensional CFD simulation modeling and based on the actual total inlet mass flow rate, the inlet mass flow rate-pressure difference database including multiple different single-pipe inlet mass flow rates and corresponding pressure differences; and obtaining the single-pipe inlet mass flow rate of each pipe in the pipe assembly through one-dimensional iterative calculation based on the multiple inlet mass flow rate-pressure difference databases and the actual total inlet mass flow rate, and using the predicted pressure difference corresponding to the single-pipe inlet mass flow rate as the inlet and outlet pressure difference of the pipe assembly.

[0007] Optionally, the step of constructing an inlet mass flow rate-pressure difference database for each pipeline group through 3D CFD simulation modeling and based on the actual total inlet mass flow rate further includes: for each pipeline group, determining the lower limit inlet mass flow rate and upper limit inlet mass flow rate corresponding to the pipeline group based on the actual total inlet mass flow rate; for each pipeline group, constructing a corresponding simulation pipeline model based on the corresponding 3D shape and size data; acquiring fluid data of the pipeline components; and for each pipeline group, calculating multiple different inlet mass flow rates and corresponding pressure differences between the lower limit inlet mass flow rate and the upper limit inlet mass flow rate based on the fluid data and the corresponding simulation pipeline model.

[0008] Optionally, for each pipeline group, the step of determining the lower limit inlet mass flow rate and the upper limit inlet mass flow rate corresponding to the pipeline group based on the actual total inlet mass flow rate further includes: for each pipeline group, determining the maximum inlet mass flow rate corresponding to the pipeline group based on the actual total inlet mass flow rate; and for each pipeline group, determining the corresponding lower limit inlet mass flow rate and the upper limit inlet mass flow rate based on the corresponding maximum inlet mass flow rate.

[0009] Optionally, the first The maximum inlet mass flow rate corresponding to each pipeline group The calculation expression is: In the formula This represents the actual total import mass flow rate. For the first The total number of pipes in a pipe group.

[0010] Optionally, the fluid data includes the fluid density and viscosity data of the fluid flowing through the pipe assembly. The viscosity data includes the effective viscosity and the bulk viscosity coefficient. The calculation expression of the Navier-Stokes equations corresponding to the simulated pipe model is as follows: , ,in, , In the formula For fluid density, For fluid velocity, For stress tensor, For pressure, It is the acceleration due to gravity. Effective adhesion, The bulk viscosity coefficient, It is a unit tensor.

[0011] Optionally, the steps of obtaining the single-pipe inlet mass flow rate of each pipe in the pipeline assembly through one-dimensional iterative calculation based on multiple inlet mass flow rate-pressure difference databases and the actual total inlet mass flow rate, and using the predicted pressure difference corresponding to the single-pipe inlet mass flow rate as the inlet and outlet pressure difference of the pipeline assembly, further include: step S51, determining the maximum and minimum calculated pressure difference based on multiple inlet mass flow rate-pressure difference databases; step S52, calculating the actual calculated pressure difference based on the maximum and minimum calculated pressure difference; step S53, determining whether the absolute value of the flow difference between the sum of the inlet mass flow rates of each pipe under the actual calculated pressure difference and the actual total inlet mass flow rate is less than a preset error value. If the determination result is yes, then the actual calculated pressure difference is used as the predicted pressure difference; if the determination result is no, then the maximum or minimum calculated pressure difference is updated based on the flow difference, and step S52 is continued.

[0012] Optionally, the step of updating the maximum or minimum calculated pressure difference based on the flow difference further includes: determining whether the flow difference is greater than zero; if the determination result is yes, then the actual calculated pressure difference is used as the maximum calculated pressure difference; if the determination result is no, then the actual calculated pressure difference is used as the minimum calculated pressure difference.

[0013] Optionally, the actual calculated pressure difference The calculation expression is: In the formula To maximize the calculated pressure difference, This is the minimum calculated pressure difference.

[0014] To address the aforementioned technical problems, this application provides a fast traffic allocation calculation system, comprising: a memory for storing instructions executable by a processor; and a processor for executing the instructions to implement the aforementioned fast traffic allocation calculation method.

[0015] To address the aforementioned technical problems, this application provides a computer-readable medium storing computer program code, which, when executed by a processor, implements the aforementioned fast flow allocation calculation method.

[0016] Compared with existing technologies, this application has the following advantages: Based on the three-dimensional shape and size data of each pipe in the pipeline assembly, all pipes are divided into multiple pipe groups, and an inlet mass flow rate-pressure difference database corresponding to each pipe group is constructed. This avoids the need for overall simulation modeling of the pipeline assembly containing a large number of pipes, while simultaneously obtaining multiple accurate inlet mass flow rates and pressure differences corresponding to different pipes within the pipeline assembly. Furthermore, a one-dimensional iterative approach is used to quickly obtain the predicted pressure difference corresponding to the actual total inlet mass flow rate and the predicted flow rate corresponding to each pipe using data from all inlet mass flow rate-pressure difference databases. This achieves efficient calculation of flow allocation and ensures the accuracy of the calculation results. Attached Figure Description

[0017] The accompanying drawings are included to provide a further understanding of this application; they are incorporated into and constitute a part of this application. The drawings illustrate embodiments of this application and, together with this specification, serve to explain the principles of this application. In the drawings: Figure 1 This is a flowchart illustrating a fast traffic allocation calculation method according to an embodiment of this application; Figure 2 This is a schematic diagram of the structure of a pipe assembly according to an embodiment of this application; Figure 3 yes Figure 1 A flowchart illustrating the sub-steps of step S4 in the middle section; Figure 4 yes Figure 3 A flowchart illustrating the sub-steps of step S41. Figure 5 This is a schematic diagram of a simulated pipeline model according to an embodiment of this application; Figure 6 yes Figure 1 A flowchart illustrating the sub-steps of step S5; and Figure 7 This is a schematic diagram of the structure of a fast traffic allocation calculation system according to an embodiment of this application. Detailed Implementation

[0018] To more clearly illustrate the technical solutions of the embodiments of this application, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the drawings described below are merely some examples or embodiments of this application. For those skilled in the art, these drawings can be applied to other similar scenarios without creative effort. Unless obvious from the context or otherwise specified, the same reference numerals in the drawings represent the same structures or operations.

[0019] As indicated in this application and claims, unless the context clearly indicates otherwise, the words "a," "an," "an," and / or "the" are not specifically singular and may include plural forms. Generally speaking, the terms "comprising" and "including" only indicate the inclusion of explicitly identified steps and elements, which do not constitute an exclusive list, and the method or apparatus may also include other steps or elements.

[0020] Unless otherwise specifically stated, the relative arrangement, numerical expressions, and values ​​of the components and steps described in these embodiments do not limit the scope of this application. It should also be understood that, for ease of description, the dimensions of the various parts shown in the drawings are not drawn to actual scale. Techniques, methods, and devices known to those skilled in the art may not be discussed in detail, but where appropriate, such techniques, methods, and devices should be considered part of the specification. In all examples shown and discussed herein, any specific values ​​should be interpreted as merely exemplary and not as limitations. Therefore, other examples of exemplary embodiments may have different values. It should be noted that similar reference numerals and letters in the following drawings denote similar items; therefore, once an item is defined in one drawing, it need not be further discussed in subsequent drawings.

[0021] In the description of this application, it should be understood that the orientation or positional relationship indicated by directional terms such as "front, back, up, down, left, right", "horizontal, vertical, horizontal" and "top, bottom" is usually based on the orientation or positional relationship shown in the accompanying drawings, and is only for the convenience of describing this application and simplifying the description. Unless otherwise stated, these directional terms do not indicate or imply that the device or element referred to must have a specific orientation or be constructed and operated in a specific orientation, and therefore should not be construed as a limitation on the scope of protection of this application; the directional terms "inner" and "outer" refer to the inner and outer contours relative to the outline of each component itself.

[0022] For ease of description, spatial relative terms such as "above," "on top of," "on the upper surface of," "above," etc., are used herein to describe the spatial positional relationship of a device or feature as shown in the figures to other devices or features. It should be understood that spatial relative terms are intended to encompass different orientations in use or operation beyond the orientation of the device as described in the figures. For example, if the device in the figures were inverted, a device described as "above" or "on top of" other devices or structures would subsequently be positioned as "below" or "under" other devices or structures. Thus, the exemplary term "above" can include both "above" and "below." The device may also be positioned in other different ways (rotated 90 degrees or in other orientations), and the spatial relative descriptions used herein will be interpreted accordingly.

[0023] Furthermore, it should be noted that the use of terms such as "first" and "second" to define components is merely for the purpose of distinguishing the corresponding components. Unless otherwise stated, these terms have no special meaning and therefore should not be construed as limiting the scope of protection of this application. In addition, although the terminology used in this application is selected from commonly known and used terms, some terms mentioned in this application's specification may have been chosen by the applicant according to his or her judgment, and their detailed meanings are explained in the relevant sections of this description. Moreover, this application should be understood not only through the actual terms used, but also through the meaning implied by each term.

[0024] It should be understood that when a component is referred to as "on another component," "connected to another component," "coupled to another component," or "in contact with another component," it can be directly on, connected to, coupled to, or in contact with that other component, or there may be an intervening component. In contrast, when a component is referred to as "directly on another component," "directly connected to," "directly coupled to," or "directly in contact with" another component, there is no intervening component. Similarly, when a first component is referred to as "electrically contacting" or "electrically coupled to" a second component, there is an electrical path between the first and second components that allows current to flow. This electrical path may include capacitors, coupled inductors, and / or other components that allow current to flow, even if there is no direct contact between the conductive components.

[0025] Flowcharts are used in this application to illustrate the operations performed by the methods according to embodiments of this application. It should be understood that the preceding or following operations are not necessarily performed in exact order. Instead, various steps can be processed in reverse order or simultaneously. Furthermore, other operations may be added to these processes, or one or more steps may be removed from these processes.

[0026] Reference Figure 1 One embodiment of this application proposes a fast flow allocation calculation method 100, which is applicable to pipe assemblies. (See also...) Figure 2 The pipe assembly 200 includes multiple pipes 21 connected in parallel, and at least two pipes 21 have different three-dimensional shape and size data. Understandably, the three-dimensional shape and size data includes outer surface data and inner surface data, and this data can be used to construct a three-dimensional model corresponding to the pipe. It should be noted that... Figure 2 The shape of pipe 21 is shown only as an example. In some embodiments, pipe 21 can be a special shape, such as a corrugated pipe, a bend, a serpentine coil, etc.

[0027] Continue to refer to Figure 1 and Figure 2The rapid flow allocation calculation method 100 includes the following steps: Step S1 is to obtain the three-dimensional shape and size data corresponding to each pipe 21. Step S2 is to divide the multiple pipes 21 into multiple pipe groups based on the three-dimensional shape and size data. Each pipe 21 in each pipe group has the same three-dimensional shape and size data. It can be understood that pipes 21 with the same three-dimensional shape and size data have the same flow rate change corresponding to the same pressure difference in the pipe assembly 200, that is, each pipe 21 in the same pipe group has the same impact on the pipe assembly 200. Therefore, through the above steps S1 and S2, a large number of pipes 21 can be grouped according to the consistency of their three-dimensional shape and size data, thereby avoiding the need for subsequent processing of each pipe 21, thus reducing the amount of data processing and improving processing efficiency.

[0028] Continue to refer to Figure 1 and Figure 2 Step S3 involves obtaining the actual total inlet mass flow rate of the pipe assembly 200. Specifically, in some examples, this is achieved by installing a flow sensor at the total inlet of the pipe assembly 200. Step S4 involves constructing an inlet mass flow rate-pressure differential database for each pipe group using 3D CFD simulation modeling and based on the actual total inlet mass flow rate. This database includes multiple different inlet mass flow rates and their corresponding pressure differentials. For example, some inlet mass flow rates and pressure differentials in the inlet mass flow rate-pressure differential database for a pipe group are shown in Table 1 below.

[0029] Table 1:

[0030] In Table 1 above, the first column represents the inlet mass flow rate in kg / s, and the second column represents the pressure difference corresponding to the inlet mass flow rate in Pa. For example, the cell in the second row and second column indicates that the pressure difference for any pipe 21 in this pipe group is 465.10 Pa when the inlet mass flow rate is 0.017045 kg / s. It is understandable that, because this embodiment utilizes three-dimensional shape and size data through high-fidelity 3D CFD simulation modeling, it is possible to obtain inlet mass flow rates and pressure differences accurate to multiple decimal places. This provides high-precision basic calculation data for subsequent calculations, ensuring the accuracy of the final result.

[0031] The following explains the specific methods for obtaining the inlet mass flow rate and pressure difference from the aforementioned inlet mass flow rate-pressure difference database. (Continue referring to...) Figure 1 and Figure 3 Step S4 includes the following sub-steps. Step S41 involves determining the lower and upper limits of the inlet mass flow rate for each pipeline group based on the actual total inlet mass flow rate. Further details can be found by referring to... Figure 4Step S41 includes the following sub-steps. Step S411 is to determine the maximum inlet mass flow rate corresponding to each pipe group based on the actual total inlet mass flow rate. Wherein, the... The maximum inlet mass flow rate corresponding to each pipeline group The calculation expression is: , In the formula This represents the actual total import mass flow rate. For the first The total number of pipes in the pipeline group. It should be noted that, using the above formula for calculating the maximum inlet mass flow rate, the number of pipes in the extreme case can be calculated. The upper limit of the maximum inlet mass flow rate of each pipe 21 in a pipe group when it accounts for the maximum flow rate of the entire pipe assembly 200. It is understandable that the above extreme case refers to the case excluding the first... In the case where there is almost no flow in the other 21 pipes outside the pipe group, at this time the... The sum of the flow rates of each pipe 21 in the pipeline group is almost equal to the actual total inlet mass flow rate.

[0032] Continue to refer to Figure 4 Step S412 involves determining the lower limit and upper limit of the inlet mass flow rate for each pipeline group based on the corresponding maximum inlet mass flow rate. Specifically, in this embodiment, a positive number not greater than 1 and with a leading zero not less than one digit is selected as the lower limit of the inlet mass flow rate, for example, 0.017045 kg / s. This ensures that even when pipeline 21 in the pipeline group is under extreme conditions, the corresponding inlet mass flow rate-pressure difference database can still be used to achieve the effectiveness of subsequent calculations. Furthermore, in this embodiment, 1.05 times the maximum inlet mass flow rate is used as the upper limit of the inlet mass flow rate for the corresponding pipeline group. This avoids the actual inlet mass flow rate of the pipeline in the pipeline group exceeding the above-mentioned maximum inlet mass flow rate under special circumstances, thereby providing effective support for subsequent calculations under special circumstances. It should be noted that this application does not limit the specific values ​​of the lower limit and upper limit of the inlet mass flow rate. In some embodiments, in order to make the data in the inlet mass flow rate-pressure difference database usable for flow distribution calculations of pipeline components 200 under more operating conditions, the lower limit of the inlet mass flow rate can be further reduced and / or the upper limit of the inlet mass flow rate can be further increased, thereby improving the versatility of the data in the inlet mass flow rate-pressure difference database.

[0033] Continue to refer to Figure 3Step S42 involves constructing a corresponding simulation pipeline model for each pipeline group based on its corresponding three-dimensional shape and size data. In this embodiment, high-fidelity computational fluid dynamics (CFD) three-dimensional simulation software is used to construct the simulation pipeline model. (Refer to...) Figure 5 In some examples, the simulated pipeline model 300 is a single-phase flow CFD model, and the simulated pipeline model 300 includes a lower bottom tank 31, a simulated pipeline 32, and an upper top tank 33 connected in sequence. The fluid flows through the lower bottom tank 31, the simulated pipeline 32, and the upper top tank 33 in sequence.

[0034] Continue to refer to Figure 3 Step S43 involves acquiring fluid data for the pipe assembly 200. This fluid data includes the fluid density and viscosity data of the fluid flowing through the pipe assembly, with the viscosity data including effective viscosity and bulk viscosity coefficient. Step S44, for each pipe group, calculates multiple different inlet mass flow rates and corresponding pressure differentials between the lower and upper inlet mass flow rates, based on the fluid data and the corresponding simulated pipe model.

[0035] In this embodiment, the pressure difference in the simulated pipeline model The calculation expression is: , In the formula The average pressure at the pipe inlet section of the simulated pipe model. This represents the average pressure at the pipe outlet section of the simulated pipe model.

[0036] Furthermore, in the simulated pipeline model 300, the first... The maximum inlet velocity of the simulated equivalent pipe 32 corresponding to each pipe group The calculation expression is: , In the formula The density of the fluid at the inlet of the simulated pipe 32. Let 32 ​​be the cross-sectional area of ​​the simulated pipe.

[0037] Accordingly, the expression for calculating the inlet velocity of the bottom bucket 31 is:

[0038] In the formula For the first The inlet flow velocity of simulated pipe 32 corresponding to each pipe group. Let 31 be the cross-sectional area of ​​the bottom bucket. The inlet flow rate of the bottom barrel 31 is denoted as .

[0039] It should also be noted that, in these examples, the expression for the Navier-Stokes equations used to describe fluid flow for the simulated pipe model 300 is as follows: , ,in,

[0040] In the formula For fluid density, For fluid velocity, For stress tensor, For pressure, It is the acceleration due to gravity. For effective adhesion, = The bulk viscosity coefficient, It is a unit tensor. = lam + sit , = In the formula lam Due to the viscosity of fluid molecules, sit It has variable viscosity. It should be noted that the fluid density... The viscosity of fluid molecules can vary with pressure. And temperature changes. Furthermore, in laminar flow, the effective viscosity coefficient only includes molecular viscosity, i.e. When considering transitional and turbulent flows, it is necessary to include shear stress-induced turbulence. Therefore, the turbulence model corresponding to the simulated pipeline model 300 mentioned above can be implemented using large eddy simulation or... Turbulence model.

[0041] Continue to refer to Figure 1 and Figure 2 Step S5 involves using multiple inlet mass flow rate-pressure difference databases and the actual total inlet mass flow rate to calculate the individual inlet mass flow rate of each pipe in the pipeline assembly through one-dimensional iterative calculation. The predicted pressure difference corresponding to the individual inlet mass flow rate is then used as the inlet and outlet pressure difference of the pipeline assembly. Further reference... Figure 6Step S5 includes the following sub-steps. Step S51 determines the maximum and minimum calculated pressure difference based on multiple inlet mass flow rate-pressure difference databases. Specifically, in this embodiment, in step S51, the minimum pressure difference is selected as the minimum calculated pressure difference from the pressure difference data in all inlet mass flow rate-pressure difference databases, and the maximum pressure difference is selected as the maximum calculated pressure difference, thereby avoiding omission of pressure difference values ​​corresponding to the actual total inlet mass flow rate. Step S52 calculates the actual calculated pressure difference based on the maximum and minimum calculated pressure differences. Specifically, in this embodiment, the actual calculated pressure difference... The calculation expression is: , In the formula To maximize the calculated pressure difference, This is the minimum calculated pressure difference.

[0042] Continue to refer to Figure 2 and Figure 6 Step S53 involves determining whether the absolute value of the difference between the sum of the inlet mass flow rates of each pipeline under the actual calculated pressure difference and the actual total inlet mass flow rate is less than a preset error value. If the determination result is yes, the actual calculated pressure difference is used as the predicted pressure difference. If the determination result is no, the maximum or minimum calculated pressure difference is updated based on the difference, and step S52 continues. Specifically, the step of updating the maximum or minimum calculated pressure difference based on the difference is as follows: determining whether the flow difference is greater than zero. If the determination result is yes, the actual calculated pressure difference is used as the maximum calculated pressure difference. If the determination result is no, the actual calculated pressure difference is used as the minimum calculated pressure difference.

[0043] It should be noted that, for actual engineering flow scenarios where the flow is single-phase and does not reach congested flow, the inlet mass flow rate and pressure difference of pipeline 21 have a monotonically increasing relationship; that is, the larger the inlet mass flow rate, the larger the corresponding pressure difference. Therefore, when the difference between the sum of the inlet mass flow rates corresponding to the actual calculated pressure difference and the actual total inlet mass flow rate is greater than zero, it indicates that the actual calculated pressure difference is greater than the true pressure difference corresponding to the actual total inlet mass flow rate. Thus, this actual calculated pressure difference can be used as the new maximum calculated pressure difference, thereby selecting a new actual calculated pressure difference using the new pressure difference range. Conversely, when the difference between the sum of the inlet mass flow rates corresponding to the actual calculated pressure difference and the actual total inlet mass flow rate is less than zero, it indicates that the actual calculated pressure difference is less than the true pressure difference corresponding to the actual total inlet mass flow rate. Therefore, this actual calculated pressure difference can be used as the new minimum calculated pressure difference. It should also be noted that when the inlet mass flow rate corresponding to the actual calculated pressure difference does not exist in the inlet mass flow rate-pressure difference database, this embodiment obtains the inlet mass flow rate corresponding to the actual calculated pressure difference in the database through interpolation. Understandably, through the above step S5, accurate results can be obtained without performing extensive simulations.

[0044] An embodiment of this application also proposes a method such as Figure 7 The fast traffic allocation calculation system 400 shown is illustrated. According to... Figure 7 The fast traffic allocation computing system 400 may include an internal communication bus 41, a processor 42, a read-only memory (ROM) 43, a random access memory (RAM) 44, and a communication port 45. When applied to a personal computer, the fast traffic allocation computing system 400 may also include a hard disk 46.

[0045] The internal communication bus 41 enables data communication between components of the fast traffic allocation computing system 400. The processor 42 can perform judgments and issue prompts. In some embodiments, the processor 42 may consist of one or more processors. The communication port 45 enables data communication between the fast traffic allocation computing system 400 and external systems. In some embodiments, the fast traffic allocation computing system 400 can send and receive information and data from a network through the communication port 45.

[0046] The fast traffic allocation computing system 400 may also include different types of program storage units and data storage units, such as a hard disk 46, read-only memory (ROM) 43, and random access memory (RAM) 44, capable of storing various data files used for computer processing and / or communication, as well as possible program instructions executed by the processor 42. The processor executes these instructions to implement the main part of the method. The results of the processor processing are transmitted to the user equipment through a communication port and displayed on the user interface.

[0047] In addition, this application also proposes a computer-readable medium storing computer program code that implements the above-described fast flow allocation calculation method 100 when executed by a processor.

[0048] The basic concepts have been described above. Obviously, for those skilled in the art, the above disclosure is merely illustrative and does not constitute a limitation of this application. Although not explicitly stated herein, those skilled in the art may make various modifications, improvements, and corrections to this application. Such modifications, improvements, and corrections are suggested in this application, and therefore remain within the spirit and scope of the exemplary embodiments of this application.

[0049] Furthermore, this application uses specific terms to describe embodiments of the application. For example, "an embodiment," "one embodiment," and / or "some embodiments" refer to a particular feature, structure, or characteristic related to at least one embodiment of the application. Therefore, it should be emphasized and noted that "an embodiment," "one embodiment," or "an alternative embodiment" mentioned twice or more in different locations in this specification do not necessarily refer to the same embodiment. In addition, certain features, structures, or characteristics in one or more embodiments of the application can be appropriately combined.

[0050] Similarly, it should be noted that, in order to simplify the description of the present application and thus aid in the understanding of one or more embodiments, the foregoing description of the embodiments of the present application sometimes combines multiple features into a single embodiment, drawing, or description thereof. However, this disclosure method does not imply that the subject matter of the present application requires more features than those mentioned in the claims. In fact, the embodiments contain fewer features than all the features of the single embodiments disclosed above.

[0051] In some embodiments, numbers describing the quantity of components and attributes are used. It should be understood that such numbers used in the description of embodiments are modified in some examples with the terms "approximately," "approximately," or "generally." Unless otherwise stated, "approximately," "approximately," or "generally" indicates that the numbers are allowed to vary by ±20%. Accordingly, in some embodiments, the numerical parameters used in the specification and claims are approximate values, which may be changed depending on the characteristics required by individual embodiments. In some embodiments, numerical parameters should take into account specified significant digits and employ a general method of digit reservation. Although the numerical ranges and parameters used to confirm their breadth of scope in some embodiments of this application are approximate values, in specific embodiments, such values ​​are set as precisely as feasible.

[0052] Some aspects of this application can be executed entirely by hardware, entirely by software (including firmware, resident software, microcode, etc.), or by a combination of hardware and software. The aforementioned hardware or software may be referred to as a "data block," "module," "engine," "unit," "component," or "system." The processor may be one or more application-specific integrated circuits (ASICs), digital signal processors (DSPs), digital signal processing devices (DAPDs), programmable logic devices (PLDs), field-programmable gate arrays (FPGAs), processors, controllers, microcontrollers, microprocessors, or combinations thereof. Furthermore, aspects of this application may manifest as computer products residing in one or more computer-readable media, including computer-readable program code. For example, computer-readable media may include, but are not limited to, magnetic storage devices (e.g., hard disks, floppy disks, magnetic tapes, etc.), optical discs (e.g., compressed CDs, digital multifunction DVDs, etc.), smart cards, and flash memory devices (e.g., cards, sticks, key drives, etc.).

[0053] A computer-readable medium may contain a propagated data signal containing computer program code, for example, on baseband or as part of a carrier wave. This propagated signal may take various forms, including electromagnetic, optical, and so on, or suitable combinations thereof. A computer-readable medium can be any computer-readable medium other than a computer-readable storage medium, which can be connected to an instruction execution system, apparatus, or device to enable communication, propagation, or transmission of a program for use. The program code located on the computer-readable medium can be propagated through any suitable medium, including radio, cable, fiber optic cable, radio frequency signals, or similar media, or any combination of the above media.

[0054] Although this application has been described with reference to specific embodiments, those skilled in the art should recognize that the above embodiments are only used to illustrate this application, and various equivalent changes or substitutions can be made without departing from the spirit of this application. Therefore, any changes or modifications to the above embodiments within the essential spirit of this application will fall within the scope of the claims of this application.

Claims

1. A fast flow allocation calculation method, characterized in that, Applicable to a pipe assembly comprising multiple pipes connected in parallel, at least two of the pipes having different three-dimensional shape and size data, the fast flow distribution calculation method includes the following steps: Obtain the three-dimensional shape and size data corresponding to each of the pipes; The plurality of pipes are divided into a plurality of pipe groups based on the three-dimensional shape and size data, wherein each pipe in each pipe group has the same three-dimensional shape and size data; Obtain the actual total inlet mass flow rate of the pipeline assembly; By using 3D CFD simulation modeling and based on the actual total inlet mass flow rate, an inlet mass flow rate-pressure differential database is constructed for each of the pipe groups. The inlet mass flow rate-pressure differential database includes multiple different single-pipe inlet mass flow rates and corresponding pressure differentials. Based on the multiple inlet mass flow rate-pressure difference databases and the actual total inlet mass flow rate, the single-pipe inlet mass flow rate of each pipe in the pipeline assembly is obtained through one-dimensional iterative calculation, and the predicted pressure difference corresponding to the single-pipe inlet mass flow rate is used as the inlet and outlet pressure difference of the pipeline assembly.

2. The fast traffic allocation calculation method as described in claim 1, characterized in that, The step of constructing an inlet mass flow rate-pressure difference database for each pipeline group through three-dimensional CFD simulation modeling and based on the actual total inlet mass flow rate further includes: For each of the pipeline groups, the lower limit inlet mass flow rate and the upper limit inlet mass flow rate corresponding to the pipeline group are determined based on the actual total inlet mass flow rate; For each of the pipe groups, a corresponding simulation pipe model is constructed based on the corresponding three-dimensional shape and size data; Obtain fluid data for the pipeline assembly; For each of the pipeline groups, based on the fluid data and the corresponding simulated pipeline model, multiple different inlet mass flow rates and corresponding pressure differences between the lower limit inlet mass flow rate and the upper limit inlet mass flow rate are calculated.

3. The fast flow allocation calculation method as described in claim 2, characterized in that, For each of the pipeline groups, the step of determining the lower limit inlet mass flow rate and the upper limit inlet mass flow rate corresponding to the pipeline group based on the actual total inlet mass flow rate further includes: For each of the pipeline groups, the maximum inlet mass flow rate corresponding to the pipeline group is determined based on the actual total inlet mass flow rate; For each of the pipeline groups, the corresponding lower limit inlet mass flow rate and the upper limit inlet mass flow rate are determined based on the corresponding maximum inlet mass flow rate.

4. The fast flow allocation calculation method as described in claim 3, characterized in that, No. The maximum inlet mass flow rate corresponding to each pipeline group The calculation expression is: ， In the formula The actual total inlet mass flow rate is given. For the first The total number of pipes described in each pipe group.

5. The fast flow allocation calculation method as described in claim 2, characterized in that, The fluid data includes the fluid density and viscosity data of the fluid flowing through the pipe assembly. The viscosity data includes the effective viscosity and the bulk viscosity coefficient. The calculation expression of the NS equation corresponding to the simulated pipe model is as follows: ， ,in, ， In the formula The fluid density is... For fluid velocity, For stress tensor, For pressure, It is the acceleration due to gravity. The effective viscosity, The viscosity coefficient of the volume is [missing information]. It is a unit tensor.

6. The fast traffic allocation calculation method as described in claim 1, characterized in that, The steps of obtaining the single-pipe inlet mass flow rate of each pipe in the pipeline assembly through one-dimensional iterative calculation based on multiple inlet mass flow rate-pressure difference databases and the actual total inlet mass flow rate, and using the predicted pressure difference corresponding to the single-pipe inlet mass flow rate as the inlet and outlet pressure difference of the pipeline assembly, further include: Step S51: Determine the maximum and minimum calculated pressure difference based on the multiple inlet mass flow rate-pressure difference databases; Step S52: Calculate the actual calculated pressure difference based on the maximum and minimum calculated pressure differences; Step S53: Determine whether the absolute value of the difference between the sum of the inlet mass flow rates of each of the pipelines under the actual calculated pressure difference and the actual total inlet mass flow rate is less than a preset error value. If the determination result is yes, then the actual calculated pressure difference is used as the predicted pressure difference. If the determination result is no, then the maximum calculated pressure difference or the minimum calculated pressure difference is updated according to the flow difference, and step S52 is continued.

7. The fast flow allocation calculation method as described in claim 6, characterized in that, The step of updating the maximum calculated pressure difference or the minimum calculated pressure difference based on the flow difference further includes: Determine whether the flow difference is greater than zero. If the determination result is yes, then the actual calculated pressure difference is taken as the maximum calculated pressure difference. If the determination result is no, then the actual calculated pressure difference is taken as the minimum calculated pressure difference.

8. The fast flow allocation calculation method as described in claim 6 or 7, characterized in that, The actual calculated pressure difference The calculation expression is: ， In the formula For the maximum calculated pressure difference, This is the minimum calculated pressure difference.

9. A fast traffic allocation calculation system, comprising: Memory is used to store instructions that can be executed by the processor; as well as A processor for executing the instructions to implement the fast traffic allocation calculation method as described in any one of claims 1-8.

10. A computer-readable medium storing computer program code that, when executed by a processor, implements the fast traffic allocation calculation method as described in any one of claims 1-8.

Images