Method and device for predicting flow characteristics of a turbine disc cavity

Abstract

The application provides a method and device for predicting flow characteristics of a turbine disc cavity, the method comprising: modeling the turbine disc cavity to obtain a disc cavity model in a target coordinate system; integrating a circumferential momentum equation of the disc cavity model in the target coordinate system in an axial direction to obtain a one-dimensional integral equation; processing the one-dimensional integral equation based on friction of a rotating disc and a stationary disc to obtain a rotating ratio radial control equation; in response to an input known parameter set, solving the rotating ratio radial control equation based on the known parameter set to obtain circumferential velocity distribution information of the rotating disc in a radius direction, wherein the known parameter set at least includes one of the following: geometric parameters of the disc, working condition parameters of the turbine disc cavity, physical property parameters of a fluid flowing in the turbine disc cavity, and boundary conditions, the boundary conditions including an inner radius rotating ratio of the rotating disc; and storing the circumferential velocity distribution information in a memory, wherein the flow characteristics of the turbine disc cavity include the circumferential velocity distribution information.

Description

Technical Field

[0001] This application relates to the field of compressor technology, and more specifically, to a method and device for predicting the flow characteristics of a turbine disk cavity. Background Technology

[0002] The turbine disk cavity is the enclosed space surrounding the turbine disk in aero engines and gas turbines, primarily used to house the rotating turbine disk and manage its operating environment. After the high-temperature, high-pressure gas expands and performs work in the turbine, some of the gas flow enters the disk cavity through the blade gaps, causing thermal load on the turbine disk and affecting structural safety and lifespan. To prevent overheating, cooling air (secondary gas flow) is typically drawn from the compressor to cool the disk cavity while simultaneously suppressing the intrusion of high-temperature gas. The flow within the disk cavity is extremely complex, involving the interaction of rotating and stationary boundary layers, directly affecting cooling efficiency and thermal stress distribution. Therefore, accurately understanding its flow characteristics is crucial for optimizing cooling design and improving engine efficiency and reliability.

[0003] The relevant technologies require a long time and have poor accuracy when predicting the flow characteristics of turbine disk cavities. Summary of the Invention

[0004] In view of this, this application provides a method and apparatus for predicting the flow characteristics of a turbine disk cavity.

[0005] One aspect of this application provides a method for predicting the flow characteristics of a turbine disk cavity, comprising: modeling the turbine disk cavity in response to a characteristic prediction command to obtain a disk cavity model in a target coordinate system, wherein the target coordinate system is a coordinate system established with the center of the disk in the disk cavity model as the origin and the perpendicular line, radius, and rotation direction of the disk, wherein the disk includes a rotating disk and / or a stationary disk; integrating the circumferential momentum equation of the disk cavity model in the target coordinate system in the axial direction to obtain a one-dimensional integral equation; processing the one-dimensional integral equation based on the friction of the rotating disk and the stationary disk to obtain a rotation ratio radial control equation, wherein the rotation... The radial control equation characterizes the radial variation of the rotation ratio of the aforementioned disk. In response to the input set of known parameters, the radial control equation for the rotation ratio is solved based on the known parameter set to obtain the circumferential velocity distribution information of the aforementioned disk along its radius. The known parameter set includes at least one of the following: the geometric parameters of the aforementioned disk, the operating parameters of the aforementioned turbine disk cavity, the physical property parameters of the fluid flowing within the aforementioned turbine disk cavity, and boundary conditions. The boundary conditions include the rotation ratio of the inner radius of the aforementioned disk. The circumferential velocity distribution information is stored in memory, wherein the flow characteristics of the aforementioned turbine disk cavity include the circumferential velocity distribution information.

[0006] According to an embodiment of this application, the circumferential velocity distribution information includes the circumferential velocity information processed at each position of the disk on the radius.

[0007] According to an embodiment of this application, the method further includes: integrating the radial dynamic equation based on the circumferential velocity information to obtain fluid pressure distribution information at different positions of the disk, and storing the fluid pressure distribution information in the memory, wherein the flow characteristics of the turbine disk cavity also include the fluid pressure distribution information.

[0008] According to an embodiment of this application, the circumferential momentum equation of the disk cavity model in the target coordinate system is integrated along the axial direction to obtain a one-dimensional integral equation, including: integrating the three-dimensional circumferential momentum equation along the axial direction to obtain a one-dimensional initial integral equation, wherein the initial integral equation includes a convection term integral and a viscous term integral, wherein the convection term integral represents the net angular momentum flux through the fluid space region of the turbine disk cavity, and the viscous term integral represents the net viscous torque acting on the boundary of the fluid space region; the initial integral equation is corrected based on the boundary conditions of the rotating disk and the stationary disk to obtain the one-dimensional integral equation.

[0009] According to an embodiment of this application, the initial integral equation is corrected based on the boundary conditions of the turntable and the stationary disk to obtain the one-dimensional integral equation, including: correcting the viscous term integral of the initial integral equation based on the no-slip condition of the turntable and the stationary disk to obtain the corrected viscous term integral; and generating the one-dimensional integral equation based on the convection term integral and the corrected viscous term integral.

[0010] According to an embodiment of this application, the one-dimensional integral equation is processed based on the friction of the turntable and the stationary disk to obtain the rotation ratio radial control equation, including: determining the turntable friction coefficient and the stationary disk friction coefficient; and processing the one-dimensional integral equation based on the turntable friction coefficient and the stationary disk friction coefficient to obtain the rotation ratio radial control equation.

[0011] According to an embodiment of this application, determining the friction coefficient of the turntable and the friction coefficient of the stationary disk includes: establishing a friction coefficient correlation formula for any target disk among the turntable and the stationary disk based on boundary layer theory; and decoupling the friction coefficient correlation formula to obtain the friction coefficient of the target disk.

[0012] According to an embodiment of this application, decoupling the above-mentioned friction coefficient correlation to obtain the friction coefficient of the target disk includes: decoupling the above-mentioned friction coefficient correlation to obtain a friction coefficient expression, wherein the above-mentioned friction coefficient expression includes an unknown rotation ratio and a known rotation Reynolds number; substituting the obtained rotation ratio parameter value into the above-mentioned friction coefficient expression to obtain the friction coefficient of the target disk.

[0013] According to an embodiment of this application, the aforementioned rotation ratio parameter value is generated as follows: the inner radius of the disk cavity at the i-th position is calculated based on the outer radius of the turntable and the i-th position of the turntable on the inner radius; the rotation ratio parameter value at the i-th position is calculated based on the inner radius of the disk cavity and the rotation ratio parameter values ​​at the previous i-1 positions, wherein the rotation ratio parameter value at the first position is calculated based on the initial value of the rotation ratio and the inner radius of the disk cavity at the first position.

[0014] According to an embodiment of this application, the rotation ratio radial control equation is solved based on the known parameter set to obtain the circumferential velocity distribution information of the turntable on the radius, including: calculating the rotation ratio parameter value at the j-th position on the turntable based on the known parameter set; solving the rotation ratio radial control equation based on the rotation ratio parameter value at the j-th position to obtain the circumferential velocity information at the j-th position; calculating the rotation ratio parameter value at the (j+1)-th position based on the rotation ratio parameter values ​​at the previous j positions, and calculating the circumferential velocity information at the (j+1)-th position based on the rotation ratio parameter value at the (j+1)-th position; and generating the circumferential velocity distribution information of the turntable on the radius based on the circumferential velocity information at multiple positions.

[0015] Another aspect of this application provides a flow characteristic prediction device for a turbine disk cavity, comprising: a modeling module for modeling the turbine disk cavity in response to a characteristic prediction command, obtaining a disk cavity model in a target coordinate system, wherein the target coordinate system is a coordinate system established with the center of the disk in the disk cavity model as the origin, and the perpendicular line, radius, and rotation direction of the disk, wherein the disk includes a rotating disk and / or a stationary disk; an integration module for integrating the circumferential momentum equation of the disk cavity model in the target coordinate system in the axial direction, obtaining a one-dimensional integral equation; and a processing module for processing the one-dimensional integral equation based on the friction of the rotating disk and the stationary disk, obtaining a rotation ratio radial control equation, wherein... The aforementioned radial control equation for the rotation ratio characterizes the radial variation of the rotation ratio of the aforementioned disk. The solution module, in response to an input set of known parameters, solves the radial control equation for the rotation ratio based on the known parameter set to obtain the circumferential velocity distribution information of the aforementioned disk along its radius. The known parameter set includes at least one of the following: the geometric parameters of the aforementioned disk, the operating parameters of the aforementioned turbine disk cavity, the physical property parameters of the fluid flowing within the aforementioned turbine disk cavity, and boundary conditions. The boundary conditions include the rotation ratio of the inner radius of the aforementioned disk. The storage module stores the circumferential velocity distribution information in memory, wherein the flow characteristics of the aforementioned turbine disk cavity include the aforementioned circumferential velocity distribution information.

[0016] Another aspect of this application provides an electronic device comprising: one or more processors; and a memory for storing one or more programs, wherein, when the one or more programs are executed by the one or more processors, the one or more processors cause the one or more processors to perform the method as described above.

[0017] Another aspect of this application provides a computer-readable storage medium storing computer-executable instructions that, when executed, are used to implement the method described above.

[0018] Another aspect of this application provides a computer program product comprising computer-executable instructions which, when executed, are used to implement the method described above.

[0019] According to an embodiment of this application, a one-dimensional integral equation is obtained by integrating the circumferential momentum equation of the disk cavity model in the target coordinate system obtained from modeling. Based on the friction between the turntable and the stationary disk, the rotation ratio radial control equation is derived. The rotation ratio radial control equation is solved using a known parameter set to obtain the circumferential velocity distribution information of the turntable along the radius. Since integrating the circumferential momentum equation in the axial direction can reduce the three-dimensional problem to a one-dimensional problem, and solving the rotation ratio radial control equation can improve the accuracy of the solution speed and the prediction of flow characteristics. Attached Figure Description

[0020] The above and other objects, features and advantages of this application will become clearer from the following description of embodiments with reference to the accompanying drawings, in which:

[0021] Figure 1 An exemplary system architecture for predicting the flow characteristics of a turbine disk cavity, according to an embodiment of this application, is shown;

[0022] Figure 2 A flowchart of a method for predicting the flow characteristics of a turbine disk cavity according to an embodiment of this application is shown;

[0023] Figure 3 A schematic diagram of the structure of a turbine disk cavity according to an embodiment of this application is shown;

[0024] Figure 4 A block diagram of a flow characteristic prediction device for a turbine disk cavity according to an embodiment of this application is shown;

[0025] Figure 5 A block diagram of an electronic device suitable for implementing the methods described above, according to an embodiment of this application, is shown. Detailed Implementation

[0026] The embodiments of this application will now be described with reference to the accompanying drawings. However, it should be understood that these descriptions are exemplary only and are not intended to limit the scope of this application. In the following detailed description, numerous specific details are set forth to provide a thorough understanding of the embodiments of this application for ease of explanation. However, it will be apparent that one or more embodiments may be implemented without these specific details. Furthermore, descriptions of well-known structures and technologies are omitted in the following description to avoid unnecessarily obscuring the concepts of this application.

[0027] The terminology used herein is for the purpose of describing particular embodiments only and is not intended to limit the scope of this application. The terms “comprising,” “including,” etc., as used herein indicate the presence of the stated features, steps, operations, and / or components, but do not exclude the presence or addition of one or more other features, steps, operations, or components.

[0028] All terms used herein (including technical and scientific terms) have the meanings commonly understood by those skilled in the art, unless otherwise defined. It should be noted that the terms used herein are to be interpreted in a manner consistent with the context of this specification, and not in an idealized or overly rigid way.

[0029] When using expressions such as "at least one of A, B and C", they should generally be interpreted in accordance with the meaning that is commonly understood by those skilled in the art (e.g., "a system having at least one of A, B and C" should include, but is not limited to, a system having A alone, a system having B alone, a system having C alone, a system having A and B, a system having A and C, a system having B and C, and / or a system having A, B and C, etc.).

[0030] In the embodiments of this application, the collection, updating, analysis, processing, use, transmission, provision, disclosure, and storage of data (e.g., including but not limited to user personal information) comply with relevant laws and regulations, are used for legitimate purposes, and do not violate public order and good morals. In particular, necessary measures have been taken to prevent unauthorized access to user personal information data and to safeguard user personal information security and network security.

[0031] Currently, flow analysis of turbine disk cavities mainly relies on the following two methods:

[0032] 1) Three-dimensional CFD simulation method

[0033] Specifically, numerical methods such as Reynolds-Averaged Navier-Stokes (RANS) or Large Eddy Simulation (LES) are used to solve the flow inside the disk cavity in three dimensions.

[0034] Technical limitations: High computational resource consumption: Fine mesh generation leads to long computation times, making it difficult to use for rapid engineering analysis and optimization design; High modeling complexity: Factors such as rotation effects, turbulence models, and heat transfer coupling need to be considered, affecting computational convergence. Not suitable for early design stages: CFD methods are too inefficient in the conceptual design or parameter optimization stages.

[0035] 2) Empirical Model Method

[0036] The specific method involves approximate calculations based on simplified flow control equations (such as mass conservation, momentum conservation, and energy conservation) combined with empirical coefficients (such as friction factor and heat transfer coefficient).

[0037] Technical limitations: Limited accuracy: Ignoring the three-dimensional effects of flow (such as secondary flow and vortex structure), treating the rotating disk cavity as a 0-dimensional parameter system, it is impossible to predict the radial distribution characteristics of parameters, resulting in large prediction deviations of key parameters (such as pressure distribution and temperature field); Narrow applicability: Empirical coefficients depend on specific working conditions or geometric conditions, making it difficult to extend to different turbine disk cavity structures; Difficult to reflect complex flow phenomena: such as the Coriolis force effect caused by rotation, the interaction of sealing flow, etc.

[0038] In view of this, embodiments of this application provide a method and apparatus for predicting the flow characteristics of a turbine disk cavity. The method includes modeling the turbine disk cavity to obtain a disk cavity model in a target coordinate system, wherein the target coordinate system is a coordinate system established with the center of the disk in the disk cavity model as the origin and the perpendicular line, radius, and rotation direction of the disk, and the disk includes a rotating disk and / or a stationary disk; integrating the circumferential momentum equation of the disk cavity model in the target coordinate system in the axial direction to obtain a one-dimensional integral equation; and processing the one-dimensional integral equation based on the friction of the rotating disk and the stationary disk to obtain the rotation ratio radial control equation. The rotation ratio radial control equation characterizes the radial variation of the disk's rotation ratio. In response to the input set of known parameters, the rotation ratio radial control equation is solved based on the known parameter set to obtain the circumferential velocity distribution information of the disk along its radius. The known parameter set includes at least one of the following: the disk's geometric parameters, the turbine disk cavity's operating parameters, the physical property parameters of the fluid flowing within the turbine disk cavity, and boundary conditions, including the inner radius rotation ratio of the disk. The circumferential velocity distribution information is stored in memory, where the flow characteristics of the turbine disk cavity include the circumferential velocity distribution information.

[0039] Figure 1 An exemplary system architecture for predicting the flow characteristics of a turbine disk cavity, according to embodiments of this application, is shown. It should be noted that... Figure 1The examples shown are merely examples of system architectures that can be applied to the embodiments of this application, in order to help those skilled in the art understand the technical content of this application, but do not mean that the embodiments of this application cannot be used in other devices, systems, environments or scenarios.

[0040] like Figure 1 As shown, the system architecture 100 according to this embodiment may include a first terminal device 101, a second terminal device 102, a third terminal device 103, a network 104, and a server 105. The network 104 serves as a medium for providing communication links between the first terminal device 101, the second terminal device 102, the third terminal device 103, and the server 105. The network 104 may include various connection types, such as wired and / or wireless communication links, etc.

[0041] Users can use the first terminal device 101, the second terminal device 102, and the third terminal device 103 to interact with the server 105 via the network 104 to receive or send messages, etc. Various communication client applications can be installed on the first terminal device 101, the second terminal device 102, and the third terminal device 103, such as shopping applications, web browser applications, search applications, instant messaging tools, email clients, and / or social media platform software, etc. (for example only).

[0042] The first terminal device 101, the second terminal device 102, and the third terminal device 103 can be various electronic devices with displays and support web browsing, including but not limited to smartphones, tablets, laptops, and desktop computers.

[0043] Server 105 can be a server that provides various services, such as a backend management server that supports websites browsed by users using the first terminal device 101, the second terminal device 102, and the third terminal device 103 (this is just an example). The backend management server can analyze and process data such as received user requests, and feed back the processing results (such as web pages, information, or data obtained or generated according to user requests) to the terminal devices.

[0044] It should be noted that the turbine disk cavity flow characteristic prediction method provided in this application embodiment can generally be executed by server 105. Correspondingly, the turbine disk cavity flow characteristic prediction device provided in this application embodiment can generally be located in server 105. The turbine disk cavity flow characteristic prediction method provided in this application embodiment can also be executed by a server or server cluster that is different from server 105 and capable of communicating with the first terminal device 101, the second terminal device 102, the third terminal device 103, and / or server 105. Correspondingly, the turbine disk cavity flow characteristic prediction device provided in this application embodiment can also be located in a server or server cluster that is different from server 105 and capable of communicating with the first terminal device 101, the second terminal device 102, the third terminal device 103, and / or server 105. Alternatively, the turbine disk cavity flow characteristic prediction method provided in this application embodiment can also be executed by the first terminal device 101, the second terminal device 102, or the third terminal device 103, or by other terminal devices different from the first terminal device 101, the second terminal device 102, or the third terminal device 103. Accordingly, the turbine disk cavity flow characteristic prediction device provided in this application embodiment can also be set in the first terminal device 101, the second terminal device 102 or the third terminal device 103, or in other terminal devices different from the first terminal device 101, the second terminal device 102 or the third terminal device 103.

[0045] It should be understood that Figure 1 The number of terminal devices, networks, and servers shown is merely illustrative. Depending on implementation needs, any number of terminal devices, networks, and servers can be included.

[0046] Before describing the flow characteristic prediction method of the turbine disk cavity in this application, the parameters involved in this application and their meanings are explained as shown in Table 1:

[0047] Table 1

[0048]

[0049]

[0050]

[0051]

[0052]

[0053] Figure 2 A flowchart of a method for predicting the flow characteristics of a turbine disk cavity according to an embodiment of this application is shown. Figure 3 A schematic diagram of the structure of a turbine disk cavity according to an embodiment of this application is shown.

[0054] like Figure 2 As shown, the method for predicting the flow characteristics of the turbine disk cavity includes operations S201~S205.

[0055] In operation S201, in response to the characteristic prediction command, the turbine disk cavity is modeled to obtain the disk cavity model in the target coordinate system. The target coordinate system is a coordinate system established with the center of the disk in the disk cavity model as the origin, and the perpendicular line, radius and rotation direction of the disk. The disk includes a turntable and / or a stationary disk.

[0056] In operation S202, the circumferential momentum equation of the disk cavity model in the target coordinate system is integrated along the axial direction to obtain a one-dimensional integral equation.

[0057] In operation S203, the one-dimensional integral equation is processed based on the friction between the turntable and the stationary disk to obtain the radial control equation for the rotation ratio. The radial control equation for the rotation ratio characterizes the variation law of the rotation ratio of the disk along the radial direction.

[0058] In operation S204, in response to the input known parameter set, the radial control equation of the rotation ratio is solved based on the known parameter set to obtain the circumferential velocity distribution information of the turntable on the radius. The known parameter set includes at least one of the following: the geometric parameters of the disk, the operating parameters of the turbine disk cavity, the physical property parameters of the fluid flowing in the turbine disk cavity, and the boundary conditions, including the inner radius rotation ratio of the turntable.

[0059] In operation S205, circumferential velocity distribution information is stored in memory, where the flow characteristics of the turbine disk cavity include circumferential velocity distribution information.

[0060] In one specific embodiment, a typical turbine disk cavity is used as an example for illustration, such as... Figure 3 As shown, the turbine disk cavity includes the following main components: Rotor Disc 1: a disc that rotates with the turbine rotor, with a radius of r_o and an angular velocity Ω; Stationary Disc 2: a disc fixed on the stator casing, with a radius typically the same as or slightly larger than the rotor disc; Disc Cavity 3: the axial clearance between the rotor and stationary discs, with a height of s (axial distance); Seal 4: located at the disk cavity outlet, used to prevent mainstream combustion gas from entering the disk cavity; Coolant Flow 5: drawn from the compressor, passing through the seal and entering the disk cavity, with a mass flow rate of ṁ.

[0061] According to the embodiments of this application, the turbine disk cavity is first modeled to obtain the disk cavity model, and a cylindrical coordinate system is first used. The coordinate system of the disk cavity model is defined as follows: r is the radial coordinate, pointing towards the outer edge of the disk; θ is the circumferential coordinate, along the rotation direction; z is the axial coordinate, perpendicular to the disk surface, with the origin being the center of the turntable.

[0062] According to embodiments of this application, in order to establish universally applicable governing equations, the following dimensionless parameters are introduced:

[0063] (1) Dimensionless radial coordinates , ,in, The outer radius of the turntable. The range of values ​​is This refers to the i-th position of the turntable on the inner radius, typically [0.3, 1]. That is, the inner radius of the disk cavity at position i.

[0064] (2) Rotation ratio : ,in: : The circumferential velocity of the fluid at radius r (m / s). : The linear velocity (m / s) of the turntable at radius r, β represents the degree of rotation of the fluid relative to the turntable, and its value ranges from [0, 1]: The fluid is completely stationary (relative to the inertial coordinate system). The fluid rotates at the same speed as the turntable (rigid body rotation).

[0065] (3) Dimensionless radial flow rate : ,in: Mass flow rate (kg / s) of radial flow within the disc cavity. : Fluid dynamic viscosity (Pa·s). It represents the magnitude of the sealing flow rate, with a typical value range of 10^4 - 10^6.

[0066] (4) Rotational Reynolds number : ,in: Fluid density (kg / m³). Angular velocity of the turntable (rad / s) It represents the ratio of rotational inertial force to viscous force, with a typical value range of 10^5 - 10^7.

[0067] According to an embodiment of this application, by adding velocity components to the cylindrical coordinate system, a disk cavity model in a rotating coordinate system (i.e., the target coordinate system) can be obtained. The governing equations for the flow within the turbine disk cavity are the Navier-Stokes equations in the cylindrical coordinate system. In the rotating coordinate system, the circumferential momentum equation is shown in the following formula:

[0068]

[0069] in: Velocity vector The fluid velocity at hi , The velocity component in the direction, τ → _θ: Circumferential viscous stress tensor, first term Transient term, second term : Convection term, third term Pressure gradient term, fourth term : Viscous diffusion term, fifth term Centrifugal force term. Among them... Let t represent the partial differential, and t be the time interval.

[0070] According to an embodiment of this application, the circumferential momentum equation of the disk cavity model in the target coordinate system is integrated in the axial direction (i.e., the z-direction in the cylindrical coordinate system) to obtain a one-dimensional integral equation. Based on the friction between the turntable and the stationary disk, the one-dimensional integral equation is derived to obtain the rotation ratio radial control equation.

[0071] According to an embodiment of this application, based on the known parameter set such as the geometric parameters of the input disk, the operating parameters of the turbine disk cavity, the physical property parameters of the fluid flowing in the turbine disk cavity, and the boundary conditions, the radial control equation of the rotation ratio is solved to obtain the circumferential velocity distribution information of the turntable on the radius, and the information is stored in memory.

[0072] According to an embodiment of this application, a one-dimensional integral equation is obtained by integrating the circumferential momentum equation of the disk cavity model in the target coordinate system obtained from modeling. Based on the friction between the turntable and the stationary disk, the rotation ratio radial control equation is derived. The rotation ratio radial control equation is solved using a known parameter set to obtain the circumferential velocity distribution information of the turntable along the radius. Since integrating the circumferential momentum equation in the axial direction can reduce the three-dimensional problem to a one-dimensional problem, and solving the rotation ratio radial control equation can improve the accuracy of the solution speed and the prediction of flow characteristics.

[0073] According to an embodiment of this application, the circumferential velocity distribution information includes circumferential velocity information processed at each position of the disk on the radius.

[0074] According to an embodiment of this application, the method further includes: integrating the radial dynamic equation based on the circumferential velocity information to obtain fluid pressure distribution information at different positions of the disk, and storing the fluid pressure distribution information in memory, wherein the flow characteristics of the turbine disk cavity also include the fluid pressure distribution information.

[0075] According to embodiments of this application, circumferential velocity distribution information The following formula can be used for calculation:

[0076]

[0077] Where β(r*) represents the radial distribution of the rotation ratio obtained from the solution. ω is the angular velocity of the turntable.

[0078] According to an embodiment of this application, the radial dynamic equation is integrated based on the circumferential velocity information to obtain the fluid pressure distribution information at different positions of the disk, as shown in the following formula:

[0079] .

[0080] According to an embodiment of this application, the circumferential momentum equation of the disk cavity model in the target coordinate system is integrated along the axial direction to obtain a one-dimensional integral equation. This includes: integrating the three-dimensional circumferential momentum equation along the axial direction to obtain a one-dimensional initial integral equation, wherein the initial integral equation includes convection term integration and viscous term integration. The convection term integration represents the net angular momentum flux through the fluid space region of the turbine disk cavity, and the viscous term integration represents the net viscous torque acting on the boundary of the fluid space region. The initial integral equation is corrected based on the boundary conditions of the rotating disk and the stationary disk to obtain the one-dimensional integral equation.

[0081] According to the embodiments of this application, based on the physical characteristics of the flow in the disk cavity, the following reasonable simplification assumptions are made: (1) Axisymmetry assumption: Since the disk cavity is an axisymmetric structure, it is assumed that all flow parameters do not change with the circumferential position θ: (2) Steady-state assumption: Consider steady-state conditions and ignore transient terms: (3) Radial flow dominance assumption: Radial flow is the main flow in the disk cavity. The axial velocity is much smaller than the radial velocity. The axial convection term can be ignored.

[0082] According to embodiments of this application, conventional methods for directly solving the three-dimensional circumferential momentum equation require knowledge of the three-dimensional velocity field distribution. This application integrates the circumferential momentum equation along the axial direction (z-direction), reducing the three-dimensional problem to a one-dimensional problem.

[0083] From the circumferential momentum equation (Stationary disk surface) to Integrating (on the surface of the turntable), we obtain the following formula:

[0084]

[0085] Where the left side of the equal sign is the integral of the convection term, and the right side is the integral of the viscous term.

[0086] The initial integral equation is corrected based on the boundary conditions of the rotating disk and the stationary disk to obtain the one-dimensional integral equation.

[0087] According to an embodiment of this application, the initial integral equation is corrected based on the boundary conditions of the rotating disk and the stationary disk to obtain a one-dimensional integral equation, including: correcting the viscous term integral of the initial integral equation based on the no-slip condition of the rotating disk and the stationary disk to obtain the corrected viscous term integral; and generating a one-dimensional integral equation based on the convection term integral and the corrected viscous term integral.

[0088] According to an embodiment of this application, the boundary conditions include the turntable surface (z=s): (No slip condition), on the stationary disk surface (z=0): (Without slip condition), using the boundary conditions, the viscous term on the right side can be expressed as the following formula:

[0089]

[0090] in: Circumferential shear stress on the surface of the turntable. : Circumferential shear stress on the surface of the stationary disc.

[0091] According to an embodiment of this application, the one-dimensional integral equation is processed based on the friction of the turntable and the stationary disk to obtain the rotation ratio radial control equation, including: determining the turntable friction coefficient and the stationary disk friction coefficient; and processing the one-dimensional integral equation based on the turntable friction coefficient and the stationary disk friction coefficient to obtain the rotation ratio radial control equation.

[0092] According to embodiments of this application, a definition of friction coefficient is introduced: (Coefficient of friction of turntable) (Stationary disk friction coefficient), at this point, the specific values ​​of the turntable friction coefficient and the stationary disk friction coefficient are determined. The one-dimensional integral equation is then processed using the turntable friction coefficient and the stationary disk friction coefficient to obtain the dimensionless form of the rotation ratio radial control equation, as shown in the following formula:

[0093]

[0094] This equation describes the variation of the rotation ratio β along the radial coordinate r*, reflecting the following physical mechanism:

[0095] (1) The first term on the right side of the equal sign: :

[0096] Physical meaning: The effect of wall friction on the rotation ratio; Item: The friction of the stationary disk slows down the fluid (the stationary disk does not rotate, thus applying a resistance effect to the fluid); ² Item: The friction of the turntable accelerates the fluid (the turntable rotates, and the fluid rotates through friction); Relative magnitude: The relative magnitude of the two determines whether the net friction effect is acceleration or deceleration.

[0097] (2) The second term on the right side of the equal sign: :

[0098] Physical meaning: Angular momentum conservation effect; Physical explanation: When a fluid flows radially outward, in order to maintain angular momentum... Conservation, circumferential velocity It must decrease as the radius r increases, resulting in Decrease; Symbol: a negative sign indicates that the term always makes the decrease. Decrease.

[0099] According to an embodiment of this application, determining the friction coefficient of the turntable and the friction coefficient of the stationary disk includes: establishing a friction coefficient correlation formula for any target disk among the turntable and the stationary disk based on boundary layer theory, wherein the target disk includes the turntable and the stationary disk; and decoupling the friction coefficient correlation formula to obtain the friction coefficient of the target disk.

[0100] According to the embodiments of this application, the friction coefficients of rotating and stationary disks are decoupled by boundary layer theory, and the friction coefficient correlation is determined by CFD data fitting. The friction coefficient correlation is then decoupled to obtain the friction coefficient of the target disk, thereby accurately predicting the radial distribution characteristics of the disk cavity flow parameters in a one-dimensional framework.

[0101] According to an embodiment of this application, the friction coefficient correlation is decoupled to obtain the friction coefficient of the target disk, including: decoupling the friction coefficient correlation to obtain the friction coefficient expression, wherein the friction coefficient expression includes an unknown rotation ratio and a known rotation Reynolds number; substituting the obtained rotation ratio parameter value into the friction coefficient expression to obtain the friction coefficient of the target disk.

[0102] In one specific embodiment, near the surfaces of the rotating and stationary disks, the flow is influenced by wall viscosity, forming a boundary layer. According to boundary layer theory, the velocity profile and wall shear stress within the boundary layer are closely related to the boundary layer thickness.

[0103] Turntable boundary layer thickness :

[0104]

[0105] Static disk boundary layer thickness :

[0106]

[0107] in: : Dimensionless boundary layer thickness coefficient, which depends on the flow state; Standard scale law for turbulent boundary layers.

[0108] According to an embodiment of this application, the turntable operates at a speed Rotation, while the fluid moves at a speed Rotation, there is a relative velocity between the two. The frictional force between the turntable and the fluid is proportional to this relative velocity. Based on the above principle, the correlation formula for the turntable friction coefficient can be obtained, as shown in the following formula:

[0109]

[0110] in, : Rotary friction coefficient constant, dimensionless, with a value range of 0.15-0.25, and a typical value of 0.20; The Reynolds number dependence reflects that the viscous effect weakens as Re increases (a characteristic of turbulent boundary layers). The rotation ratio depends on the relative velocity; when β→1 (the fluid and the turntable move at the same speed). →0 (no relative motion, no friction), when β→0 (the fluid does not rotate). It has reached its maximum value.

[0111] The power (1 / 4) of the above correlation comes from the classical power-law distribution of the velocity profile within the turbulent boundary layer. Combined with the relationship between shear stress and velocity gradient, the relationship between friction coefficient and the 1 / 4 power of relative velocity can be derived.

[0112] According to an embodiment of this application, the stationary disk remains stationary, and the fluid flows at a velocity... The relative velocity between the two is the rotation. The frictional force exerted by the stationary disk on the fluid is proportional to this relative velocity. The correlation for the stationary disk friction coefficient is shown in the following formula:

[0113]

[0114] in, : Static disc friction coefficient constant, dimensionless, with a value range of 0.15-0.25, and a typical value of 0.20; : Reynolds number dependency, with the same physical meaning as λ_R; The rotation ratio depends on the relative velocity; when β→0 (the fluid does not rotate), λ_S→0 (no relative motion, no friction); when β→1 (the fluid and the turntable move at the same speed), λ_S reaches its maximum value.

[0115] According to an embodiment of this application, the constant of the turntable friction coefficient is... and static friction coefficient constant It can be calculated in the following way:

[0116] Method 1: Based on boundary layer thickness correlation

[0117] According to boundary layer theory, the relationship between the friction coefficient and the boundary layer thickness is shown in the following formula:

[0118]

[0119]

[0120] The coefficient 0.18 is derived from the empirical constant of the turbulent boundary layer.

[0121] Boundary layer thickness coefficient and Values ​​for: For a smooth turntable: ≈ 0.08-0.12, typical value 0.10; for smooth stationary disks: ≈ 0.08-0.12, typical value 0.10; for rough surfaces or disks with bosses: γ value increases by 10%-30%.

[0122] Method 2: Data fitting based on Computational Fluid Dynamics (CFD) (data-driven method)

[0123] For a specific disk cavity geometry, the following steps can be used to determine it. and :

[0124] (1) CFD baseline calculation: Select representative operating conditions (such as 3-5 different Re and (Combined) to perform full three-dimensional CFD calculations to obtain high-precision flow field data; extract the shear stress distribution on the surface of the rotating and stationary disks, and calculate the friction coefficient. (r) and (r).

[0125] (2) Parameter fitting: The parameters obtained from CFD and Substitute the following relation:

[0126]

[0127]

[0128] The optimal fit is determined using the least squares method or other optimization algorithms. and value.

[0129] (3) Verification and correction: using the fitted results and Solve the one-dimensional model; compare the prediction results with the CFD results; if the error is large, adjust the values ​​and iterate until the set accuracy is reached.

[0130] According to the embodiments of this application, the obtained rotation ratio parameter value is... Substituting into the friction coefficient expression, we obtain the friction coefficient of the target disk.

[0131] According to embodiments of this application, computational fluid dynamics calculations are performed. and It fully utilizes the high precision of CFD, requiring only a small number of CFD calculations (e.g., 3-5 operating conditions), and the fitted coefficients are specific and applicable to particular geometries; furthermore, this method is effective in determining... and After that, all subsequent operating conditions can be quickly calculated using a one-dimensional model.

[0132] According to an embodiment of this application, the rotation ratio parameter value is generated as follows: the inner radius of the disk cavity at the i-th position is calculated based on the outer radius of the turntable and the i-th position of the turntable on the inner radius; the rotation ratio parameter value at the i-th position is calculated based on the inner radius of the disk cavity and the rotation ratio parameter values ​​at the previous i-1 positions, wherein the rotation ratio parameter value at the first position is calculated based on the initial value of the rotation ratio and the inner radius of the disk cavity at the first position.

[0133] According to an embodiment of this application, within the radius of the disk cavity... place, Let be the inner radius of the disk cavity at position i. That is, the i-th position of the turntable on the inner radius. That is, the outer radius of the turntable, the sealing flow rate entering the turntable cavity, and the boundary conditions are shown in the following formula:

[0134]

[0135] in The initial value of the rotation ratio at the sealing position is usually determined based on the sealing structure: For simple radial clearance seals: For maze sealing or brush sealing: .

[0136] Based on the above formula, the rotation ratio parameter value at each position can be calculated. .

[0137] According to embodiments of this application, the rotation ratio control equation based on the aforementioned rotation ratio parameter values ​​simplifies the complex three-dimensional flow problem into a one-dimensional ordinary differential equation problem, laying the foundation for efficient solution. Furthermore, this equation is based on rigorous physical derivation, ensuring the accuracy of the predictions.

[0138] According to an embodiment of this application, the rotation ratio radial control equation is solved based on a known parameter set to obtain the circumferential velocity distribution information of the turntable on the radius, including: for the j-th position on the turntable, calculating the rotation ratio parameter value at the j-th position based on the known parameter set; solving the rotation ratio radial control equation based on the rotation ratio parameter value at the j-th position to obtain the circumferential velocity information at the j-th position; calculating the rotation ratio parameter value at the (j+1)-th position based on the rotation ratio parameter values ​​at the previous j positions, and calculating the circumferential velocity information at the (j+1)-th position based on the rotation ratio parameter value at the (j+1)-th position; and generating the circumferential velocity distribution information of the turntable on the radius based on the circumferential velocity information at multiple positions.

[0139] In one specific embodiment, substituting the friction coefficient correlation into the rotation ratio control equation yields a closed first-order ordinary differential equation, as shown in the following formula:

[0140]

[0141] Where f on the right side is a known function, namely the above rotation ratio control equation.

[0142] This ordinary differential equation can be solved using various numerical methods, such as the fourth-order Runge-Kutta method and the shooting method.

[0143] In practice, first input the following parameters: Geometric parameters: r_i (inner radius), r_o (outer radius), s (disc spacing); Operating parameters: (Speed) (Sealing flow rate); Fluid properties: (density), (Viscosity); Boundary conditions: (Inner radius rotation ratio). Then calculate the dimensionless parameter: ; ; .

[0144] Initialize the following parameters: radial grid number N (e.g., N = 100-1000); radial step size. Initial conditions .

[0145] Based on the above, the solution is solved step by step: from Beginning, gradually moving towards Advancement: For To N: Current = Calculate the coefficient of friction: ; .

[0146] The derivative is further calculated using the following formula:

[0147]

[0148] β can be updated based on the calculation results, for example, using the Runge-Kutta method, as shown in the following formula:

[0149]

[0150]

[0151]

[0152]

[0153]

[0154] Next step i

[0155] Finally, the radial distribution of the rotation ratio β(r*) is output. This allows for the calculation of the circumferential velocity distribution. Information on fluid pressure distribution.

[0156] In the above embodiments, when the number of grids N = 100, the computation time is approximately 0.1 seconds (single-core CPU); when the number of grids N = 1000, the computation time is approximately 1 second, representing an efficiency improvement of approximately 20,000 to 80,000 times compared to conventional CFD computation (6-24 hours). Regarding memory usage, this method only requires storing a one-dimensional array, with a memory usage of <1MB; compared to CFD grids (several GB), it saves more than 10^4 times the memory.

[0157] According to the embodiments of this application, compared with the prior art, the one-dimensional flow characteristic prediction model for turbine disk cavity proposed in this invention has the following advantages and positive effects:

[0158] 1. Significantly improves computational efficiency while maintaining good accuracy.

[0159] A one-dimensional flow model is constructed by decoupling the rotation ratio governing equation from the friction coefficient. This simplifies the complex three-dimensional flow problem into a one-dimensional rotation ratio differential equation and introduces a physical mechanism-based friction coefficient correlation, avoiding the need for full three-dimensional mesh generation and solving while retaining the ability to describe core physical phenomena. Compared to traditional CFD methods that require hours or even days of computation, this method completes the calculation in just seconds.

[0160] 2. Accurately capture radial parameter distribution characteristics

[0161] Based on the physical mechanism, the one-dimensional radial differential equation is combined with the boundary layer theory. The governing equation in this invention retains the radial distribution information and can predict the radial variation trend of key parameters such as rotation ratio β and pressure coefficient Cp, rather than just giving the average value. This overcomes the limitation of traditional 0-dimensional / empirical models that cannot reflect the radial distribution of parameters and can accurately predict the flow field structure characteristics inside the disk cavity.

[0162] 3. Wide range of applications and good parameter versatility.

[0163] The correlation between friction coefficient and dimensionless sealing flow rate is established based on physical dimensional analysis and data fitting. By introducing dimensionless parameters and combining theoretical derivation with CFD data fitting, a widely applicable correlation for friction coefficient is developed. This correlation is applicable to disc cavity structures with different geometric dimensions, rotational speeds, and sealing flow rates, thus eliminating the dependence of traditional empirical models on specific operating conditions.

[0164] Figure 4 A block diagram of a flow characteristic prediction device for a turbine disk cavity according to an embodiment of this application is shown.

[0165] like Figure 4 As shown, the flow characteristic prediction device 400 for turbine disk cavity includes a modeling module 410, an integration module 420, a processing module 430, a solution module 440, and a storage module 450.

[0166] The modeling module 410 is used to model the turbine disk cavity in response to the characteristic prediction command, and obtain the disk cavity model in the target coordinate system. The target coordinate system is a coordinate system established with the center of the disk in the disk cavity model as the origin, and the perpendicular line, radius and rotation direction of the disk. The disk includes a turntable and / or a stationary disk.

[0167] The integration module 420 is used to integrate the circumferential momentum equation of the disk cavity model in the target coordinate system in the axial direction to obtain a one-dimensional integral equation.

[0168] The processing module 430 is used to process the one-dimensional integral equation based on the friction of the turntable and the stationary disk to obtain the radial control equation of the rotation ratio, wherein the radial control equation of the rotation ratio characterizes the variation law of the rotation ratio of the disk along the radial direction.

[0169] The solver module 440 is used to respond to the input known parameter set and solve the radial control equation of the rotation ratio based on the known parameter set to obtain the circumferential velocity distribution information of the turntable on the radius. The known parameter set includes at least one of the following: the geometric parameters of the disk, the operating parameters of the turbine disk cavity, the physical property parameters of the fluid flowing in the turbine disk cavity, and the boundary conditions, including the inner radius rotation ratio of the turntable.

[0170] The storage module 450 is used to store circumferential velocity distribution information in memory, wherein the flow characteristics of the turbine disk cavity include circumferential velocity distribution information.

[0171] According to an embodiment of this application, a one-dimensional integral equation is obtained by integrating the circumferential momentum equation of the disk cavity model in the target coordinate system obtained from modeling. Based on the friction between the turntable and the stationary disk, the rotation ratio radial control equation is derived. The rotation ratio radial control equation is solved using a known parameter set to obtain the circumferential velocity distribution information of the turntable along the radius. Since integrating the circumferential momentum equation in the axial direction can reduce the three-dimensional problem to a one-dimensional problem, and solving the rotation ratio radial control equation can improve the accuracy of the solution speed and the prediction of flow characteristics.

[0172] Any one or more of the modules, submodules, units, and subunits according to the embodiments of this application, or at least part of the functions of any one or more of them, can be implemented in one module. Any one or more of the modules, submodules, units, and subunits according to the embodiments of this application can be implemented by dividing them into multiple modules. Any one or more of the modules, submodules, units, and subunits according to the embodiments of this application can be at least partially implemented as hardware circuits, such as field-programmable gate arrays (FPGAs), programmable logic arrays (PLAs), systems-on-a-chip, systems-on-a-substrate, systems-on-package, application-specific integrated circuits (ASICs), or implemented by hardware or firmware in any other reasonable manner by integrating or packaging circuits, or implemented in any one of software, hardware, and firmware, or in a suitable combination of any of these. Alternatively, one or more of the modules, submodules, units, and subunits according to the embodiments of this application can be at least partially implemented as computer program modules, which, when run, can perform corresponding functions.

[0173] For example, any multiple of the modeling module 410, integration module 420, processing module 430, solving module 440, and storage module 450 can be combined into one module / unit / subunit, or any one of these modules / units / subunits can be split into multiple modules / units / subunits. Alternatively, at least part of the functionality of one or more of these modules / units / subunits can be combined with at least part of the functionality of other modules / units / subunits and implemented in one module / unit / subunit. According to embodiments of this application, at least one of the modeling module 410, integration module 420, processing module 430, solving module 440, and storage module 450 can be at least partially implemented as hardware circuitry, such as a field-programmable gate array (FPGA), a programmable logic array (PLA), a system-on-a-chip, a system-on-a-substrate, a system-on-package, an application-specific integrated circuit (ASIC), or any other reasonable means of integrating or packaging the circuitry, or implemented in software, hardware, or firmware, or in any suitable combination of any of these three implementation methods. Alternatively, at least one of the modeling module 410, integration module 420, processing module 430, solving module 440, and storage module 450 may be implemented at least partially as a computer program module, which can perform corresponding functions when the computer program module is run.

[0174] It should be noted that the flow characteristic prediction device part of the turbine disk cavity in the embodiments of this application corresponds to the flow characteristic prediction method part of the turbine disk cavity in the embodiments of this application. The description of the flow characteristic prediction device part of the turbine disk cavity is specifically referred to the flow characteristic prediction method part of the turbine disk cavity, and will not be repeated here.

[0175] Figure 5 A block diagram of an electronic device suitable for implementing the methods described above, according to an embodiment of this application, is shown. Figure 5 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.

[0176] like Figure 5As shown, an electronic device 500 according to an embodiment of this application includes a processor 501, which can perform various appropriate actions and processes according to a program stored in a read-only memory (ROM) 502 or a program loaded from a storage portion 508 into a random access memory (RAM) 503. The processor 501 may include, for example, a general-purpose microprocessor (e.g., a CPU), an instruction set processor and / or an associated chipset and / or a special-purpose microprocessor (e.g., an application-specific integrated circuit (ASIC)), etc. The processor 501 may also include onboard memory for caching purposes. The processor 501 may include a single processing unit or multiple processing units for performing different actions of the method flow according to an embodiment of this application.

[0177] RAM 503 stores various programs and data required for the operation of electronic device 500. Processor 501, ROM 502, and RAM 503 are interconnected via bus 504. Processor 501 executes various operations of the method flow according to embodiments of this application by executing programs in ROM 502 and / or RAM 503. It should be noted that the programs may also be stored in one or more memories other than ROM 502 and RAM 503. Processor 501 may also execute various operations of the method flow according to embodiments of this application by executing programs stored in said one or more memories.

[0178] According to embodiments of this application, the electronic device 500 may further include an input / output (I / O) interface 505, which is also connected to a bus 504. The electronic device 500 may also include one or more of the following components connected to the input / output (I / O) interface 505: an input section 506 including a keyboard, mouse, etc.; an output section 507 including a cathode ray tube (CRT), liquid crystal display (LCD), etc., and a speaker, etc.; a storage section 508 including a hard disk, etc.; and a communication section 509 including a network interface card such as a LAN card, modem, etc. The communication section 509 performs communication processing via a network such as the Internet. A drive 510 is also connected to the input / output (I / O) interface 505 as needed. A removable medium 511, such as a disk, optical disk, magneto-optical disk, semiconductor memory, etc., is installed on the drive 510 as needed so that computer programs read from it can be installed into the storage section 508 as needed.

[0179] According to embodiments of this application, the method flow according to embodiments of this application can be implemented as a computer software program. For example, embodiments of this application include a computer program product comprising a computer program carried on a computer-readable storage medium, the computer program containing program code for performing the methods shown in the flowchart. In such embodiments, the computer program can be downloaded and installed from a network via communication section 509, and / or installed from removable medium 511. When the computer program is executed by processor 501, it performs the functions defined in the system of embodiments of this application. According to embodiments of this application, the systems, devices, apparatuses, modules, units, etc., described above can be implemented by computer program modules.

[0180] This application also provides a computer-readable storage medium, which may be included in the device / apparatus / system described in the above embodiments; or it may exist independently and not assembled into the device / apparatus / system. The computer-readable storage medium carries one or more programs, which, when executed, implement the method according to the embodiments of this application.

[0181] According to embodiments of this application, the computer-readable storage medium can be a non-volatile computer-readable storage medium. Examples include, but are not limited to: portable computer disks, hard disks, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), portable compact disk read-only memory (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination thereof. In this application, the computer-readable storage medium can be any tangible medium containing or storing a program that can be used by or in conjunction with an instruction execution system, apparatus, or device.

[0182] For example, according to embodiments of this application, a computer-readable storage medium may include the ROM 502 and / or RAM 503 described above and / or one or more memories other than ROM 502 and RAM 503.

[0183] Embodiments of this application also include a computer program product comprising a computer program containing program code for performing the methods provided in the embodiments of this application. When the computer program product is run on an electronic device, the program code is used to enable the electronic device to implement the methods provided in the embodiments of this application.

[0184] When the computer program is executed by the processor 501, it performs the functions defined in the system / apparatus of this application embodiment. According to the embodiments of this application, the systems, apparatuses, modules, units, etc., described above can be implemented by computer program modules.

[0185] In one embodiment, the computer program may rely on a tangible storage medium such as an optical storage device or a magnetic storage device. In another embodiment, the computer program may also be transmitted and distributed in the form of signals over a network medium, and may be downloaded and installed via the communication section 509, and / or installed from a removable medium 511. The program code contained in the computer program can be transmitted using any suitable network medium, including but not limited to: wireless, wired, etc., or any suitable combination thereof.

[0186] According to embodiments of this application, program code for executing the computer programs provided in the embodiments of this application can be written in any combination of one or more programming languages. Specifically, these computational programs can be implemented using high-level procedural and / or object-oriented programming languages, and / or assembly / machine languages. Programming languages ​​include, but are not limited to, languages ​​such as Java, C++, Python, "C", or similar programming languages. The program code can be executed entirely on the user's computing device, partially on the user's device, partially on a remote computing device, or entirely on a remote computing device or server. In cases involving remote computing devices, the remote computing device can be connected to the user's computing device via any type of network, including a local area network (LAN) or a wide area network (WAN), or it can be connected to an external computing device (e.g., via the Internet using an Internet service provider).

[0187] 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. In this regard, each block in a flowchart or block diagram may represent a module, segment, or portion of code containing 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, may 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. Those skilled in the art will understand that the features described in the various embodiments of this application can be combined and / or combined in various ways, even if such combinations are not explicitly described in this application. In particular, without departing from the spirit and teachings of this application, the features described in the various embodiments of this application can be combined and / or combined in various ways. All such combinations and / or combinations fall within the scope of this application.

[0188] The embodiments of this application have been described above. However, these embodiments are merely illustrative and not intended to limit the scope of this application. Although various embodiments have been described above, this does not mean that the measures in the various embodiments cannot be used advantageously in combination. This application does not depart from its scope, and those skilled in the art can make various substitutions and modifications, all of which should fall within the scope of this application.

Claims

1. A method for predicting the flow characteristics of a turbine disk cavity, characterized in that, include: In response to the characteristic prediction command, the turbine disk cavity is modeled to obtain a disk cavity model in the target coordinate system, wherein the target coordinate system is a coordinate system established with the center of the disk in the disk cavity model as the origin, and the perpendicular line, radius and rotation direction of the disk, and the disk includes a turntable and / or a stationary disk. By integrating the circumferential momentum equation of the disk cavity model in the target coordinate system along the axial direction, a one-dimensional integral equation is obtained. The one-dimensional integral equation is processed based on the friction between the turntable and the stationary disk to obtain the radial control equation for the rotation ratio, wherein the radial control equation for the rotation ratio characterizes the radial variation of the rotation ratio of the disk. In response to the input known parameter set, the radial control equation of the rotation ratio is solved based on the known parameter set to obtain the circumferential velocity distribution information of the turntable on the radius. The known parameter set includes at least one of the following: the geometric parameters of the turntable, the operating parameters of the turbine disk cavity, the physical property parameters of the fluid flowing in the turbine disk cavity, and the boundary conditions, wherein the boundary conditions include the inner radius rotation ratio of the turntable. The circumferential velocity distribution information is stored in memory, wherein the flow characteristics of the turbine disk cavity include the circumferential velocity distribution information.

2. The method according to claim 1, characterized in that, The circumferential velocity distribution information includes the circumferential velocity information processed at each position of the disk on the radius; The method further includes: Integrating the radial dynamic equation based on the circumferential velocity information, the fluid pressure distribution information at different positions of the disk is obtained, and the fluid pressure distribution information is stored in the memory. The flow characteristics of the turbine disk cavity also include the fluid pressure distribution information.

3. The method according to claim 1, characterized in that, Integrating the circumferential momentum equation of the disk cavity model in the target coordinate system along the axial direction yields a one-dimensional integral equation, including: Integrating the three-dimensional circumferential momentum equation along the axial direction yields a one-dimensional initial integral equation, which includes a convection term integral and a viscous term integral. The convection term integral represents the net angular momentum flux through the fluid space region of the turbine disk cavity, and the viscous term integral represents the net viscous torque acting on the boundary of the fluid space region. The initial integral equation is corrected based on the boundary conditions of the turntable and the stationary disk to obtain the one-dimensional integral equation.

4. The method according to claim 3, characterized in that, The initial integral equation is corrected based on the boundary conditions of the rotating disk and the stationary disk to obtain the one-dimensional integral equation, which includes: The viscous integral of the initial integral equation is corrected based on the no-slip condition of the turntable and the stationary disk to obtain the corrected viscous integral. The one-dimensional integral equation is generated based on the integral of the convection term and the integral of the corrected viscous term.

5. The method according to claim 1, characterized in that, The one-dimensional integral equation is processed based on the friction between the turntable and the stationary disk to obtain the radial control equation for the rotation ratio, which includes: Determine the coefficient of friction of the turntable and the coefficient of friction of the stationary disk; The one-dimensional integral equation is processed based on the rotational friction coefficient and the stationary friction coefficient to obtain the rotational ratio radial control equation.

6. The method according to claim 1, characterized in that, Determining the rotational friction coefficient of the turntable and the stationary friction coefficient of the stationary disk includes: For any target disk among the turntable and stationary disk, a correlation formula for the friction coefficient of the target disk is established based on boundary layer theory; The friction coefficient of the target disk is obtained by decoupling the correlation of the friction coefficient.

7. The method according to claim 6, characterized in that, Decoupling the friction coefficient correlation yields the friction coefficient of the target disk, including: The friction coefficient correlation is decoupled to obtain the friction coefficient expression, wherein the friction coefficient expression includes the unknown rotation ratio and the known rotation Reynolds number; Substituting the obtained rotation ratio parameter value into the friction coefficient expression, the friction coefficient of the target disk is obtained.

8. The method according to claim 7, characterized in that, The rotation ratio parameter value is generated in the following manner: Calculate the inner radius of the disk cavity at the i-th position based on the outer radius of the turntable and the i-th position of the turntable on the inner radius; The rotation ratio parameter value at position i is calculated based on the inner radius of the disk cavity and the rotation ratio parameter values ​​at the first i-1 positions. The rotation ratio parameter value at position i is calculated based on the initial value of the rotation ratio and the inner radius of the disk cavity at position i.

9. The method according to claim 1, characterized in that, Solving the radial control equation for the rotation ratio based on the known parameter set yields the circumferential velocity distribution information of the turntable along its radius, including: For the j-th position on the turntable, calculate the rotation ratio parameter value for the j-th position based on the known parameter set; The rotation ratio radial control equation is solved based on the rotation ratio parameter value at the j-th position to obtain the circumferential velocity information at the j-th position. Calculate the rotation ratio parameter value at position (j+1) based on the rotation ratio parameter values ​​at the first j positions, and then calculate the circumferential velocity information at position (j+1) based on the rotation ratio parameter value at position (j+1). The circumferential velocity distribution information of the turntable on the radius is generated based on the circumferential velocity information at multiple locations.

10. A flow characteristic prediction device for a turbine disk cavity, characterized in that, include: The modeling module is used to model the turbine disk cavity in response to the characteristic prediction command, and obtain the disk cavity model in the target coordinate system. The target coordinate system is a coordinate system established with the center of the disk in the disk cavity model as the origin, and the perpendicular line, radius and rotation direction of the disk. The disk includes a turntable and / or a stationary disk. The integration module is used to integrate the circumferential momentum equation of the disk cavity model in the target coordinate system in the axial direction to obtain a one-dimensional integral equation. The processing module is used to process the one-dimensional integral equation based on the friction between the turntable and the stationary disk to obtain the rotation ratio radial control equation, wherein the rotation ratio radial control equation characterizes the variation law of the rotation ratio of the disk along the radial direction. The solution module is used to solve the radial control equation of the rotation ratio based on the input known parameter set in response to the known parameter set, so as to obtain the circumferential velocity distribution information of the turntable on the radius. The known parameter set includes at least one of the following: the geometric parameters of the disk, the operating parameters of the turbine disk cavity, the physical property parameters of the fluid flowing in the turbine disk cavity, and the boundary conditions, wherein the boundary conditions include the inner radius rotation ratio of the turntable. A storage module is used to store the circumferential velocity distribution information in memory, wherein the flow characteristics of the turbine disk cavity include the circumferential velocity distribution information.

Images