Method, apparatus, device, and medium for foreign object fluid force prediction in a tube array environment

Abstract

The application provides a method, device, equipment and medium for foreign object fluid force prediction in a tube array environment, and relates to the technical field of foreign object fluid force prediction in a tube array environment. The method comprises the following steps: obtaining information parameters of a foreign object to be predicted, the information parameters comprising geometric parameters and working condition parameters; determining the category of the foreign object to be predicted according to the geometric parameters of the foreign object to be predicted; matching the conservative lift envelope line and the conservative drag envelope line corresponding to the category of the foreign object to be predicted; performing dimension reduction on the conservative lift envelope line and the conservative drag envelope line according to the information parameters of the foreign object to be predicted, to obtain a dimensional lift power spectral density and a dimensional drag power spectral density; generating random phase spectrums based on the dimensional lift power spectral density and the dimensional drag power spectral density respectively, and obtaining fluctuating lift time history signals and fluctuating drag time history signals in the time domain respectively through inverse Fourier transform. The method improves the defects of low calculation efficiency, incomplete working condition coverage and poor universality of traditional technologies.

Description

Technical Field

[0001] This application relates to the field of heterogeneous material force prediction technology in tube array environments, and in particular to a method, apparatus, equipment and medium for heterogeneous material force prediction in tube array environments. Background Technology

[0002] The presence of foreign objects (such as metal debris and welding residue) in tube array equipment such as steam generators in nuclear power plants and heat exchangers in petrochemical plants can lead to leakage risks. Current mainstream methods for predicting the forces of such foreign objects include computational fluid dynamics (CFD) simulation and direct experimental measurement. However, simulating the flow field by numerically solving the Navier-Stokes equations requires individual modeling of each foreign object in the tube array and high-resolution mesh simulation of vortex shedding between tubes. Due to the strong nonlinearity of vortex shedding between tubes, the time step must be extremely short (t < 10⁻⁴ s), with a single calculation taking hours to days, making it difficult to support real-time engineering decisions. Direct experimental measurement can measure the lift of a single foreign object, i.e., measuring the lift time history signal FL(t) of the foreign object in an experimental setup. However, due to high-frequency turbulent pulsation noise interference, repeated tests are required to cover extreme conditions, resulting in high costs and poor data generalization. Current simulation and experimental methods cannot simultaneously achieve efficiency, accuracy, and universality, exhibiting poor adaptability and leading to inaccurate predictions of foreign object forces, thus affecting the safety of nuclear power plants. Summary of the Invention

[0003] This application provides a method, apparatus, device, and medium for predicting the physical strength of heterogeneous materials in a tube array environment, in order to improve the technical problems of low efficiency and poor versatility of existing methods for predicting the physical strength of heterogeneous materials in tube array equipment.

[0004] This application provides a method for predicting the physical properties of heterogeneous materials in a tube array environment, the method comprising: Acquire information parameters of the foreign object to be predicted, including geometric parameters and operating condition parameters; The category of the foreign object to be predicted is determined based on its geometric parameters. Match the conservative lift envelope and conservative drag envelope corresponding to the foreign object category to be predicted; Based on the information parameters of the foreign object to be predicted, the dimensions of the conservative lift envelope and the conservative drag envelope are restored to obtain the dimensional lift power spectral density and the dimensional drag power spectral density. Random phase spectra are generated based on the dimensional lift power spectral density and the dimensional drag power spectral density, respectively, and the time-domain pulsating lift time history signal and pulsating drag time history signal are obtained by inverse Fourier transform. The steady-state resistance of the foreign object to be predicted is obtained, and the total resistance of the foreign object to be predicted is obtained based on the steady-state resistance of the foreign object to be predicted and the time history signal of the pulsating resistance.

[0005] In an exemplary embodiment of this application, obtaining the conservative lift envelope includes: Acquire the lift time history signal of the foreign object; The lift time history signal of the same type of foreign object is subjected to Fourier transform to generate a dimensional power spectral density. The dimensional power spectral density is made dimensionless to obtain the power spectral density under different operating parameters. By mapping the power spectral density to a unified dimensionless coordinate system, a conservative lift envelope can be obtained.

[0006] In an exemplary embodiment of this application, the dimensional power spectral density is made dimensionless to achieve power spectral density under different operating parameters. The formula is as follows:

[0007] Where Φ(f) is the dimensionless PSD lift spectrum; Up is the inter-tube flow velocity; ρ is the fluid density; d is the equivalent diameter of the foreign object; l is the length of the foreign object; and D is the diameter of the heat exchange tube.

[0008] In an exemplary embodiment of this application, obtaining the conservative lift envelope includes: Map the power spectral density to a unified dimensionless coordinate system; The maximum power spectral density of each segment is taken as a node, and spline interpolation is used to connect the nodes to generate a conservative lift envelope.

[0009] In an exemplary embodiment of this application, the category of the foreign object to be predicted includes cylindrical foreign objects, plate-shaped foreign objects, spiral gasket-shaped foreign objects, and three-dimensional complex foreign objects.

[0010] In an exemplary embodiment of this application, obtaining the steady-state resistance includes: The parameterized drag coefficient C is established based on the category of the foreign object to be predicted and the Reynolds number Re. D Associative, ; Calculate the steady-state resistance based on the correlation formula. ; Where A is the windward area, Up is the flow velocity between pipes, ρ is the fluid density, d is the equivalent diameter of the foreign object, and l is the length of the foreign object.

[0011] In an exemplary embodiment of this application, obtaining the conservative resistance envelope includes: Obtain the time history signal of the resistance of the foreign object; The resistance time history signal of the same type of foreign object is subjected to Fourier transform to generate a dimensional power spectral density. The dimensional power spectral density is made dimensionless to obtain the power spectral density under different operating parameters. By mapping the power spectral density to a unified dimensionless coordinate system, a conservative drag envelope can be obtained.

[0012] A second aspect of this application provides an apparatus for predicting the physical properties of foreign objects in a tube array environment. The apparatus includes: an information parameter acquisition module for acquiring information parameters of the foreign object to be predicted, the information parameters including geometric parameters and operating condition parameters; a foreign object category determination module for determining the category of the foreign object to be predicted based on its geometric parameters; an envelope acquisition module for acquiring a conservative lift envelope and a conservative drag envelope corresponding to the category of the foreign object to be predicted; a power spectral density acquisition module for restoring the dimensions of the conservative lift envelope and the conservative drag envelope based on the information parameters of the foreign object to be predicted, to obtain dimensional lift power spectral density and dimensional drag power spectral density; a time history signal calculation module for generating random phase spectra based on the dimensional lift power spectral density and the dimensional drag power spectral density, and obtaining time-domain pulsating lift time history signals and pulsating drag time history signals respectively through inverse Fourier transform; and a total drag calculation module for obtaining the steady-state drag of the foreign object to be predicted, and obtaining the total drag of the foreign object to be predicted based on the steady-state drag of the foreign object to be predicted and the pulsating drag time history signals.

[0013] A third aspect of this application provides an electronic device comprising: one or more processors; and a storage device for storing one or more programs, which, when executed by the one or more processors, cause the electronic device to implement the method for predicting the physical properties of heterogeneous materials in a tube array environment as described in any of the preceding claims.

[0014] A fourth aspect of this application provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a computer processor, causes the computer to perform the method for predicting the forces of heterogeneous materials in a tube array environment as described above.

[0015] In combination with existing technologies, the beneficial effects of this application are as follows: Existing methods for predicting the dynamics of foreign matter in tube array environments fail to balance efficiency, accuracy, and versatility, resulting in poor adaptability. This application obtains conservative lift and drag envelopes through experiments or fluid dynamics simulations. By substituting the parameters of the foreign matter to be predicted into these envelopes, and then using dimensional reduction and Fourier transform, the pulsating lift and drag time histories are obtained. This application effectively overcomes the shortcomings of traditional methods, such as ignoring lift pulsations and failing to adapt to tube array layouts. It also addresses the limitations of traditional techniques, including low computational efficiency, incomplete operating condition coverage, and poor versatility, thus significantly improving the efficiency and versatility of fluid dynamics prediction. Attached Figure Description

[0016] The accompanying drawings, which are incorporated in and form part of this specification, illustrate embodiments consistent with this application and, together with the description, serve to explain the principles of this application. It is obvious that the drawings described below are merely some embodiments of this application, and those skilled in the art can obtain other drawings based on these drawings without any inventive effort.

[0017] In the attached diagram: Figure 1 A schematic diagram illustrating a method for predicting the physical forces of heterogeneous materials in a tube array environment according to an embodiment of this application; Figure 2 This is a schematic diagram of the process for obtaining the conservative lift envelope provided in one embodiment of this application; Figure 3 This is a schematic diagram of the process for obtaining the conservative resistance envelope provided in one embodiment of this application; Figure 4 This is an exemplary conservative lift envelope provided in one embodiment of this application; Figure 5 This is an exemplary conservative resistance envelope provided in one embodiment of this application; Figure 6 A schematic diagram of the structure of a computer system suitable for implementing the electronic device of the present application is shown. Detailed Implementation

[0018] The following specific examples illustrate the implementation of this application. Those skilled in the art can easily understand other advantages and effects of this application from the content disclosed in this specification. This application can also be implemented or applied through other different specific embodiments. Various details in this specification can also be modified or changed based on different viewpoints and applications without departing from the spirit of this application. In the absence of conflict, the following embodiments and features in the embodiments can be combined with each other.

[0019] It should be noted that the illustrations provided in the following embodiments are only schematic representations of the basic concept of this application. The drawings only show the components related to this application and are not drawn according to the actual number, shape and size of the components in the actual implementation. In the actual implementation, the form, quantity and proportion of each component can be arbitrarily changed, and the layout of the components may also be more complex.

[0020] In the following description, numerous details are explored to provide a more thorough explanation of embodiments of the present application. However, it will be apparent to those skilled in the art that embodiments of the present application may be practiced without these specific details. In other embodiments, well-known structures and devices are shown in block diagram form rather than in detail to avoid obscuring embodiments of the present application.

[0021] It should be noted that the terms such as "upper," "lower," "left," "right," "middle," and "one" used in this specification are merely for clarity of description and are not intended to limit the scope of this application. Any changes or adjustments to their relative relationships, without substantially altering the technical content, shall also be considered as part of the scope of this application.

[0022] This application provides a method for predicting the physical properties of heterogeneous materials in a tube array environment, in order to improve the problems of low efficiency, incomplete coverage of operating conditions, and poor versatility of traditional prediction methods.

[0023] Please see Figure 1 , Figure 1 A schematic diagram of a method for predicting the physical properties of heterogeneous materials in an exemplary tube array environment is shown. The method specifically includes the following steps: Step S110: Obtain information parameters of the foreign object to be predicted, including geometric parameters and operating condition parameters; Based on the geometric parameters of the foreign objects, they are classified into cylindrical foreign objects, plate-shaped foreign objects, spiral gasket-shaped foreign objects, and three-dimensional complex foreign objects. By classifying the foreign objects, conservative lift envelopes and conservative drag envelopes are obtained for each category, making the prediction of pulsating lift and pulsating drag of the foreign object to be predicted more accurate.

[0024] Taking a cylindrical foreign object as an example, the geometric parameters of the cylindrical foreign object include the equivalent diameter, and the formula for calculating the equivalent diameter d is:

[0025] Where Ap is the projected area of ​​the object in the direction of resistance, and l is the length of the object.

[0026] Geometric parameters also include the aspect ratio γ of the foreign object, which is an equivalent parameter. The formula for calculating the aspect ratio γ is:

[0027] The geometric parameters of a foreign object may also include the area ratio, etc.

[0028] Fluid operating parameters include fluid density ρ, inter-pipe gap velocity Up, etc. The inter-pipe gap velocity Up can be converted from the incoming flow velocity U and the pipe diameter ratio using the formula Up=U×P / (PD); By classifying foreign objects, an accurate envelope can be selected to obtain more accurate hydrodynamic predictions.

[0029] Step S120: Determine the category of the foreign object to be predicted based on its geometric parameters.

[0030] The category of the foreign object to be predicted is determined based on geometric parameters. The categories include one of the following: cylindrical foreign object, plate-shaped foreign object, spiral gasket-shaped foreign object, and three-dimensional complex foreign object.

[0031] Step S130: Match the conservative lift envelope and conservative drag envelope corresponding to the foreign object category to be predicted.

[0032] Different types of foreign objects correspond to different conservative lift envelopes and conservative drag envelopes. Select the appropriate conservative lift envelope and conservative drag envelope based on the type of foreign object to be predicted.

[0033] Please see Figure 2 , Figure 2 A schematic diagram illustrating an exemplary process for obtaining the conservative lift envelope is shown. The steps for obtaining the conservative lift envelope are as follows: Step S210: Obtain the lift time history signal of the foreign object.

[0034] The lift time history signals of four typical foreign objects were obtained through experimental measurement or numerical simulation.

[0035] Step S220: Perform Fourier transform on the lift time history signal of the same type of foreign object to generate a dimensional power spectral density.

[0036] Step S230: The dimensional power spectral density is transformed into a dimensionless form to obtain the power spectral density under different operating parameters.

[0037] By eliminating the influence of geometric and flow parameters through dimensionless transformation, power spectral density (PSD) data under different flow velocities and pipe diameters are unified, making PSDs under different conditions comparable.

[0038] In some embodiments, the dimensional power spectral density is made dimensionless to achieve power spectral density under different operating parameters. The formula is as follows:

[0039] Where Φ(f) is the dimensionless PSD lift spectrum; Up is the inter-tube flow velocity; ρ is the fluid density; d is the equivalent diameter of the foreign object; l is the length of the foreign object; and D is the diameter of the heat exchange tube.

[0040] The original PSD's dimension ΦL(f) represents the square of the force per hertz (dimension: N). 2 / Hz), through the square of the characteristic force scale ( )2 (Dimension: N) 2 ) and the reciprocal of the characteristic time scale Introducing a characteristic time scale (dimension: Hz), the final overall Φ L Becoming a dimensionless quantity .

[0041] Step S240: Map the power spectral density to a unified dimensionless coordinate system to obtain the conservative lift envelope.

[0042] Mapping the normalized PSD generated from the time-domain lift data of the same shaped foreign object under the same tube array arrangement to a unified fD / U p Dimensionless coordinate system. Obtain its conservative upper bound envelope. , which represents the maximum pulsation intensity of the pulsating lift.

[0043] Obtaining the conservative lift envelope involves mapping the power spectral density to a uniform dimensionless coordinate system; taking the maximum power spectral density value of each segment as a node, and connecting the nodes using spline interpolation to generate the conservative lift envelope.

[0044] Please see Figure 4 , Figure 4 The normalized lift-drag (PSD) spectrum of a cylindrical foreign object is shown. The envelope formula (represented by the dashed line) obtained by connecting the nodes using spline interpolation is shown below:

[0045] This envelope represents a conservative upper limit for the pulsating lift of all cylindrical foreign objects and can be used to predict the turbulence-induced lift experienced by foreign objects due to wear.

[0046] Please see Figure 3 , Figure 3 The diagram illustrates the process of obtaining the conservative resistance envelope, which includes: Step S310: Obtain the resistance time history signal of the foreign object.

[0047] Resistance time history signals of four typical foreign objects were obtained through experimental measurement or numerical simulation.

[0048] Step S320: Perform Fourier transform on the resistance time history signal of the same type of foreign object to generate a dimensional power spectral density.

[0049] Step S330: The dimensional power spectral density is transformed into a dimensionless form to obtain the power spectral density under different operating parameters.

[0050] Step S340: Map the power spectral density to a unified dimensionless coordinate system to obtain the conservative drag envelope.

[0051] Please see Figure 5 , Figure 5 The normalized drag PSD spectrum of the cylindrical foreign object is shown. The envelope formula (represented by the dashed line) obtained by connecting the nodes using spline interpolation is shown below:

[0052] The conservative resistance envelope is a conservative upper limit of the turbulent resistance of all cylindrical foreign objects. It can be used to predict the turbulent-induced resistance experienced by foreign objects during wear, and the envelope can cover the pulsating peak characteristics of all flow velocities under this shape and size.

[0053] Please continue reading. Figure 1 Step S140: Based on the information parameters of the foreign object to be predicted, the dimensions of the conservative lift envelope and the conservative drag envelope are restored to obtain the dimensional lift power spectral density and the dimensional drag power spectral density.

[0054] Step S150: Generate random phase spectra based on the dimensional lift power spectral density and the dimensional drag power spectral density, and obtain the time-domain pulsating lift time history signal and pulsating drag time history signal respectively through inverse Fourier transform.

[0055] First, the dimensions of the normalized drag PSD spectrum are restored, and a random phase spectrum is generated based on the dimensional PSD. Then, the time-domain turbulent force signal is reconstructed through inverse Fourier transform to obtain the fluctuating lift. Time history signal, pulsation resistance Time history signal;

[0056]

[0057] The actual lift of the foreign object is the pulsating lift. .

[0058] Step S160: Obtain the steady-state resistance of the foreign object to be predicted, and obtain the total resistance of the foreign object to be predicted based on the steady-state resistance of the foreign object to be predicted and the pulsating resistance time history signal.

[0059] In some embodiments, obtaining the steady-state resistance includes: The parameterized drag coefficient C is established based on the category of the foreign object to be predicted and the Reynolds number Re. D Associative, ; Calculate the steady-state resistance based on the correlation formula. ; Where A is the windward area, Up is the flow velocity between pipes, ρ is the fluid density, d is the equivalent diameter of the foreign object, and l is the length of the foreign object.

[0060] The total resistance of the foreign object to be predicted is the sum of the steady-state resistance and the pulsating resistance. .

[0061] A second aspect of this application provides an apparatus for predicting the physical properties of foreign objects in a tube array environment. The apparatus includes: an information parameter acquisition module for acquiring information parameters of the foreign object to be predicted, the information parameters including geometric parameters and operating condition parameters; a foreign object category determination module for determining the category of the foreign object to be predicted based on its geometric parameters; an envelope acquisition module for acquiring a conservative lift envelope and a conservative drag envelope corresponding to the category of the foreign object to be predicted; a power spectral density acquisition module for restoring the dimensions of the conservative lift envelope and the conservative drag envelope based on the information parameters of the foreign object to be predicted, to obtain dimensional lift power spectral density and dimensional drag power spectral density; a time history signal calculation module for generating random phase spectra based on the dimensional lift power spectral density and the dimensional drag power spectral density, and obtaining time-domain pulsating lift time history signals and pulsating drag time history signals respectively through inverse Fourier transform; and a total drag calculation module for obtaining the steady-state drag of the foreign object to be predicted, and obtaining the total drag of the foreign object to be predicted based on the steady-state drag of the foreign object to be predicted and the pulsating drag time history signals.

[0062] It should be noted that the heterogeneous material force prediction device in the tube array environment provided in the above embodiments and the heterogeneous material force prediction method in the tube array environment provided in the above embodiments belong to the same concept. The specific operation methods of each module and unit have been described in detail in the method embodiments and will not be repeated here. In practical applications, the heterogeneous material force prediction device in the tube array environment provided in the above embodiments can be assigned to different functional modules as needed, that is, the internal structure of the device can be divided into different functional modules to complete all or part of the functions described above. This is not a limitation here.

[0063] Embodiments of this application also provide an electronic device, including: one or more processors; and a storage device for storing one or more programs, which, when executed by the one or more processors, cause the electronic device to implement the heterogeneous material force prediction method in the tube array environment provided in the above embodiments.

[0064] Figure 6 A schematic diagram of a computer system suitable for implementing the embodiments of this application is shown. It should be noted that... Figure 6 The computer system 600 of the electronic device shown is merely an example and should not impose any limitation on the functionality and scope of use of the embodiments of this application.

[0065] like Figure 6As shown, the computer system 600 includes a Central Processing Unit (CPU) 601, which can perform various appropriate actions and processes based on programs stored in Read-Only Memory (ROM) 602 or programs loaded from Storage Unit 608 into Random Access Memory (RAM) 603, such as performing the methods described in the above embodiments. The RAM 603 also stores various programs and data required for system operation. The CPU 601, ROM 602, and RAM 603 are interconnected via a bus 604. An Input / Output (I / O) interface 605 is also connected to the bus 604.

[0066] The following components are connected to I / O interface 605: an input section 606 including a keyboard, mouse, etc.; an output section 607 including a cathode ray tube (CRT), liquid crystal display (LCD), etc., and speakers, etc.; a storage section 608 including a hard disk, etc.; and a communication section 609 including a network interface card such as a LAN (Local Area Network) card, modem, etc. The communication section 609 performs communication processing via a network such as the Internet. A drive 610 is also connected to I / O interface 605 as needed. A removable medium 611, such as a disk, optical disk, magneto-optical disk, semiconductor memory, etc., is installed on drive 610 as needed so that computer programs read from it can be installed into storage section 608 as needed.

[0067] Specifically, according to embodiments of this application, the processes described above with reference to the flowcharts can be implemented as computer software programs. For example, embodiments of this application include a computer program product comprising a computer program carried on a computer-readable medium, the computer program including a computer program for performing the methods shown in the flowcharts. In such embodiments, the computer program can be downloaded and installed from a network via communication section 609, and / or installed from removable medium 611. When the computer program is executed by central processing unit (CPU) 601, it performs various functions defined in the system of this application.

[0068] It should be noted that the computer-readable medium shown in the embodiments of this application can be a computer-readable signal medium or a computer-readable storage medium, or any combination of the two. A computer-readable storage medium can be, for example, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination thereof. More specific examples of a computer-readable storage medium may include, but are not limited to: an electrical connection having one or more wires, a portable computer disk, a hard disk, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM), flash memory, optical fiber, portable compact disc read-only memory (CD-ROM), optical storage device, magnetic storage device, or any suitable combination thereof. In this application, a computer-readable signal medium may include a data signal propagated in baseband or as part of a carrier wave, carrying a computer-readable computer program. Such propagated data signals can take various forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination thereof. Computer-readable signal media can also be any computer-readable medium other than computer-readable storage media, which can send, propagate, or transmit a program for use by or in connection with an instruction execution system, apparatus, or device. The computer program contained on the computer-readable medium can be transmitted using any suitable medium, including but not limited to wireless, wired, etc., or any suitable combination thereof.

[0069] The flowcharts and block diagrams in the accompanying drawings illustrate the architecture, functionality, and operation of possible implementations of systems, methods, and computer program products according to various embodiments of this application. Each block in a flowchart or block diagram may represent a module, segment, or portion of code, which contains one or more executable instructions for implementing a specified logical function. It should also be noted that in some alternative implementations, the functions indicated in the blocks may occur in a different order than those indicated in the drawings. For example, two consecutively indicated blocks may actually be executed substantially in parallel, and they may sometimes be executed in reverse order, depending on the functions involved. It should also be noted that each block in a block diagram or flowchart, and combinations of blocks in a block diagram or flowchart, can be implemented using a dedicated hardware-based system that performs the specified function or operation, or using a combination of dedicated hardware and computer instructions.

[0070] The units described in the embodiments of this application can be implemented in software or hardware, and the described units can also be located in a processor. The names of these units do not necessarily limit the specific unit itself.

[0071] Another aspect of this application provides a computer-readable storage medium storing a computer program that, when executed by a computer's processor, causes the computer to perform the heterogeneous material force prediction method in a tube array environment as described above. This computer-readable storage medium may be included in the electronic device described in the above embodiments, or it may exist independently and not assembled into the electronic device.

[0072] This application separates the lift components of different foreign objects through projection, extracts the conservative envelope of the power spectral density (PSD) of turbulent lift, generates a random phase spectrum based on the envelope, and reconstructs the time-domain turbulent lift signal through inverse Fourier transform, finally synthesizing a high-precision lift time-domain signal FL(t). This solves the shortcomings of traditional methods, such as ignoring lift pulsations and not adapting to tube array layouts, and addresses the limitations of low computational efficiency, incomplete operating condition coverage, and poor versatility of traditional techniques. It provides an efficient and reliable core input for assessing the wear life of foreign objects in high-risk scenarios such as nuclear power plant steam generators. This application is applicable to the accurate prediction method and system for the lift and drag of foreign objects such as metal scrap, welding slag, and gaskets in the tube bundle gap flow field of key heat exchange equipment such as nuclear power plant steam generators, thermal power plant condensers, and chemical heat exchangers. It is suitable for rapid calculation and conservative load assessment of the mass forces of foreign objects under various sizes, shapes, and flow velocity conditions. Therefore, this application effectively overcomes some practical problems in the prior art, thus having high utilization value and practical significance.

[0073] The above embodiments are merely illustrative of the principles and effects of this application and are not intended to limit this application. Any person skilled in the art can modify or alter the above embodiments without departing from the spirit and scope of this application. Therefore, all equivalent modifications or alterations made by those skilled in the art without departing from the spirit and technical concept disclosed in this application should still be covered by the claims of this application.

Claims

1. A method for predicting the physical properties of heterogeneous materials in a tube array environment, characterized in that, The method includes: Acquire information parameters of the foreign object to be predicted, including geometric parameters and operating condition parameters; The category of the foreign object to be predicted is determined based on its geometric parameters. Match the conservative lift envelope and conservative drag envelope corresponding to the foreign object category to be predicted; Based on the information parameters of the foreign object to be predicted, the dimensions of the conservative lift envelope and the conservative drag envelope are restored to obtain the dimensional lift power spectral density and the dimensional drag power spectral density. Random phase spectra are generated based on the dimensional lift power spectral density and the dimensional drag power spectral density, respectively, and the time-domain pulsating lift time history signal and pulsating drag time history signal are obtained by inverse Fourier transform. The steady-state resistance of the foreign object to be predicted is obtained, and the total resistance of the foreign object to be predicted is obtained based on the steady-state resistance of the foreign object to be predicted and the time history signal of the pulsating resistance.

2. The method for predicting the dynamics of heterogeneous materials in a tube array environment according to claim 1, characterized in that, The acquisition of the conservative lift envelope includes: Acquire the lift time history signal of the foreign object; The lift time history signal of the same type of foreign object is subjected to Fourier transform to generate a dimensional power spectral density. The dimensional power spectral density is made dimensionless to obtain the power spectral density under different operating parameters. By mapping the power spectral density to a unified dimensionless coordinate system, a conservative lift envelope can be obtained.

3. The method for predicting the dynamics of heterogeneous materials in a tube array environment according to claim 2, characterized in that, The dimensional power spectral density is made dimensionless to obtain the power spectral density under different operating parameters. The formula is as follows: Where Φ(f) is the dimensionless PSD lift spectrum; Up is the inter-tube flow velocity; ρ is the fluid density; d is the equivalent diameter of the foreign object; l is the length of the foreign object; and D is the diameter of the heat exchange tube.

4. The method for predicting the dynamics of heterogeneous materials in a tube array environment according to claim 2, characterized in that, Obtaining the conservative lift envelope includes: Map the power spectral density to a unified dimensionless coordinate system; The maximum power spectral density of each segment is taken as a node, and spline interpolation is used to connect the nodes to generate a conservative lift envelope.

5. The method for predicting the dynamics of heterogeneous materials in a tube array environment according to claim 1, characterized in that, The categories of foreign objects to be predicted include cylindrical foreign objects, plate-shaped foreign objects, spiral gasket-shaped foreign objects, and three-dimensional complex foreign objects.

6. The method for predicting the dynamics of heterogeneous materials in a tube array environment according to claim 1, characterized in that, The steady-state resistance is obtained by: The parameterized drag coefficient C is established based on the category of the foreign object to be predicted and the Reynolds number Re. D Associative, ； Calculate the steady-state resistance based on the correlation formula. ； Where A is the windward area, Up is the flow velocity between pipes, ρ is the fluid density, d is the equivalent diameter of the foreign object, and l is the length of the foreign object.

7. The method for predicting the dynamics of heterogeneous materials in a tube array environment according to claim 1, characterized in that, The acquisition of the conservative resistance envelope includes: Obtain the time history signal of the resistance of the foreign object; The resistance time history signal of the same type of foreign object is subjected to Fourier transform to generate a dimensional power spectral density. The dimensional power spectral density is made dimensionless to obtain the power spectral density under different operating parameters. By mapping the power spectral density to a unified dimensionless coordinate system, a conservative drag envelope can be obtained.

8. A device for predicting the physical properties of heterogeneous materials in a tube array environment, characterized in that, The device includes: The information parameter acquisition module acquires information parameters of the foreign object to be predicted, including geometric parameters and operating condition parameters. The foreign object category determination module determines the category of the foreign object to be predicted based on its geometric parameters. The envelope acquisition module acquires the conservative lift envelope and conservative drag envelope corresponding to the foreign object category to be predicted; The power spectral density acquisition module restores the dimensions of the conservative lift envelope and the conservative drag envelope based on the information parameters of the foreign object to be predicted, and obtains the dimensional lift power spectral density and the dimensional drag power spectral density. The time history signal calculation module generates random phase spectra based on the dimensional lift power spectral density and the dimensional drag power spectral density, and obtains the time-domain pulsating lift time history signal and pulsating drag time history signal respectively through inverse Fourier transform. The total resistance calculation module obtains the steady-state resistance of the foreign object to be predicted, and obtains the total resistance of the foreign object to be predicted based on the steady-state resistance of the foreign object to be predicted and the pulsating resistance time history signal.

9. An electronic device, characterized in that, The electronic device includes: One or more processors; A storage device for storing one or more programs, which, when executed by one or more processors, cause the electronic device to implement the method for predicting the forces of heterogeneous materials in a tube array environment as claimed in any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that, It stores a computer program that, when executed by the computer's processor, causes the computer to perform the method for predicting the forces of heterogeneous materials in a tube array environment as described in any one of claims 1 to 7.

Images