Transverse jet mixing characteristics prediction method and device

Abstract

The application provides a transverse jet flow mixing characteristic prediction method and device, wherein the method comprises the following steps: obtaining physical parameters of a main flow, initial physical parameters and initial geometric parameters of a transverse jet flow, and a position of a target prediction plane; obtaining positions of penetration centers of the transverse jet flow at each intermediate prediction plane based on the physical parameters of the main flow, the initial physical parameters and the initial geometric parameters of the transverse jet flow; and obtaining a trajectory of the transverse jet flow and coordinates of a penetration center of the transverse jet flow at the target prediction plane based on the positions of the penetration centers of the transverse jet flow at the intermediate prediction planes. The transverse jet flow mixing characteristic prediction method and device provided by the application can realize more efficient, rapid and relatively accurate prediction of the mixing characteristic of a small-flow transverse jet flow by regarding the trajectory flow of the transverse jet flow as a series of steady flows composed of mutually non-interfering jet flow elements, performing force analysis on any jet flow element, and performing recursion based on the Lagrange recursion algorithm.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of aero-engine, and particularly relates to a transverse jet mixing characteristic prediction method and device. BACKGROUND

[0002] The development of high-performance aero-engine technology makes the aero-engine develop in the direction of high pressure ratio and high turbine inlet temperature (i.e. combustor outlet temperature). This development trend affects the reliability and durability of turbine blades. In order to ensure the safe and reliable operation of turbine blades, the hot spot of the outlet temperature field of the main combustor cannot be too high and must meet certain distribution characteristics, which requires the ability of accurate prediction, careful design and precise adjustment of the combustor outlet temperature distribution. Therefore, the regulation of the combustor outlet temperature distribution has become one of the key technologies in the design of the combustor.

[0003] The flame tube of the combustor is divided into a combustion zone and a mixing zone, and air enters each zone through each group of holes and grooves. Among them, the mixing gas is mixed with the high-temperature gas generated in the combustion zone to adjust the outlet temperature distribution. The mixing effect of this transverse jet plays an important role in the regulation of the combustor outlet temperature distribution. However, with the development of military and civilian high-performance engine technology, there is a demand for high-temperature rise and low-emission combustion technology. The amount of mixing gas of these advanced combustor technologies has decreased from the traditional 20% to less than 10%. Generally, the combustor with a mixing gas ratio less than 10% is defined as a low-mixing air combustor.

[0004] In the low-mixing air combustor, due to the decrease of the mixing gas ratio to less than 10%, the interaction between the main flow and the transverse jet is dominated by the main flow swirl, which cannot be ignored. At this time, the flow mixing of the non-uniform swirl inflow and the transverse jet and the strong momentum exchange between them significantly affect the flow state downstream of the transverse jet, and further change the temperature distribution at the combustor outlet relative to the case when the inflow is uniform. In this case, the regulation of the combustor outlet temperature distribution depends on the combined action of the non-uniform swirl inflow and the small flow transverse jet mixing. Therefore, for the low-mixing air combustor, it is very important to explore the mixing characteristics of the small flow transverse jet under the condition of swirl inflow.

[0005] The current gas turbine combustor is provided with a plurality of groups of hole grooves for air to enter into different functional areas from the head to the outlet, wherein the mixing jet flow is affected by the transverse airflow from the flame tube, the jet flow trajectory is curved, and a series of vortex structures are generated downstream of the jet flow under the action of the transverse airflow. For the mixing section of the combustor, the transverse jet flow characteristics of the mixing hole affect the flow field structure of the mixing section of the combustor, and further have an indiscernible connection with the outlet temperature distribution of the combustor. The important characteristics of the transverse jet flow mixing include the flow trajectory of the transverse jet flow and the jet flow penetration depth.

[0006] For the estimation of the transverse jet flow mixing characteristics, the conventional research method is to derive the dimensionless control equation set from the conservation equation set, and to obtain the accurate solution of the jet flow trajectory under specific conditions. Under the condition of uniform flow, since the main flow is simple, the problem can be regarded as a two-dimensional problem, and the partial differential form of the basic control equation set can be simplified into a common differential form, so that the analytic solution is obtained by integration. However, for the small flow transverse jet flow in the low-mixing air combustor, the momentum change in three-dimensional space needs to be considered, so it is inevitable to solve the analytic solution of the partial differential equation set in three-dimensional problems, and such method is no longer feasible.

[0007] Therefore, how to estimate the mixing characteristics of the small flow transverse jet flow in the low-mixing air combustor becomes a topic to be solved in the field. SUMMARY

[0008] The present application provides a transverse jet flow mixing characteristic estimation method and device, which solves the defect that the characteristics of the small flow transverse jet flow in the low-mixing air combustor cannot be estimated in the prior art, and realizes the estimation of the characteristics of the small flow transverse jet flow in the low-mixing air combustor.

[0009] The present application provides a transverse jet flow mixing characteristic estimation method, comprising:

[0010] Obtaining the physical parameters of the main flow, the initial physical parameters and the initial geometric parameters of the transverse jet flow, and the position of the target estimation plane;

[0011] Based on the physical parameters of the main flow, the initial physical parameters and the initial geometric parameters of the transverse jet flow, the position of the penetration center of the transverse jet flow at each intermediate estimation plane is obtained;

[0012] Based on the position of the penetration center of the transverse jet flow at each intermediate estimation plane, the trajectory of the transverse jet flow and the coordinates of the penetration center of the transverse jet flow at the target estimation plane are obtained.

[0013] The method for predicting the mixing characteristics of a transverse jet flow provided by the application comprises the following steps: obtaining the position of the penetration center of the transverse jet flow on each intermediate prediction plane based on the physical parameters of the main flow, the initial physical parameters and the initial geometric parameters of the transverse jet flow, wherein the method comprises the following steps:

[0014] In the case that the position of the penetration center of the transverse jet flow on the kth intermediate prediction plane does not reach the position of the target prediction plane, the position of the penetration center of the transverse jet flow on the (k+1)th intermediate prediction plane is obtained based on the physical parameters of the main flow, the position of the penetration center of the transverse jet flow on the kth intermediate prediction plane, the physical parameters and the geometric parameters of the kth jet flow element included in the transverse jet flow.

[0015] The method for predicting the mixing characteristics of a transverse jet flow provided by the application comprises the following steps: obtaining the position of the penetration center of the transverse jet flow on each intermediate prediction plane based on the physical parameters of the main flow, the initial physical parameters and the initial geometric parameters of the transverse jet flow, wherein the method comprises the following steps:

[0016] The physical parameters and the geometric parameters of the (k+1)th jet flow element included in the transverse jet flow are obtained based on the physical parameters of the main flow, the physical parameters and the geometric parameters of the kth jet flow element.

[0017] The position of the penetration center of the transverse jet flow on the (k+1)th intermediate prediction plane is obtained based on the position of the penetration center of the transverse jet flow on the kth intermediate prediction plane and the physical parameters of the (k+1)th jet flow element.

[0018] The method for predicting the mixing characteristics of a transverse jet flow provided by the application comprises the following steps: obtaining the position of the penetration center of the transverse jet flow on each intermediate prediction plane based on the physical parameters of the main flow, the initial physical parameters and the initial geometric parameters of the transverse jet flow, wherein the method comprises the following steps:

[0019] The mass of the (k+1)th jet flow element and the total momentum increment of the kth jet flow element are obtained based on the physical parameters of the main flow, the physical parameters and the geometric parameters of the kth jet flow element.

[0020] The velocity of the (k+1)th jet flow element is obtained based on the mass of the (k+1)th jet flow element, the momentum increment of the kth jet flow element and the velocity of the kth jet flow element.

[0021] The transverse jet mixing characteristic prediction method provided by the application comprises the following steps:

[0022] The first momentum increment, the first mass increment, the second momentum increment and the second mass increment of the kth jet element are obtained based on the physical parameters of the main flow, the physical parameters and the geometric parameters of the kth jet element.

[0023] The mass of the (k+1)th jet element is obtained based on the mass of the kth jet element, the first mass increment and the second mass increment, and the total momentum increment of the kth jet element is obtained based on the first momentum increment and the second momentum increment.

[0024] The first momentum increment and the first mass increment are caused by the shear entrainment of the transverse jet, and the second momentum increment and the second mass increment are caused by the forced entrainment of the transverse jet.

[0025] The transverse jet mixing characteristic prediction method provided by the application comprises the following steps:

[0026] The position of the penetration center of the transverse jet in the (k+1)th intermediate prediction plane is obtained based on the position of the penetration center of the transverse jet in the kth intermediate prediction plane and the physical parameters of the (k+1)th jet element.

[0027] The transverse jet mixing characteristic prediction device provided by the application comprises:

[0028] The acquisition module is configured to acquire the physical parameters of the main flow, the initial physical parameters and the initial geometric parameters of the transverse jet, and the position of the target prediction plane.

[0029] The recursion module is configured to acquire the position of the penetration center of the transverse jet in each intermediate prediction plane based on the physical parameters of the main flow, the initial physical parameters and the initial geometric parameters of the transverse jet.

[0030] The prediction module is configured to acquire the trajectory of the transverse jet and the coordinates of the penetration center of the transverse jet in the target prediction plane based on the positions of the penetration centers of the transverse jet in the intermediate prediction planes.

[0031] The application further provides an electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the transverse jet flow mixing characteristic prediction method according to any one of the above when executing the program.

[0032] The application further provides a non-transitory computer readable storage medium, which stores a computer program, wherein the computer program is executed by a processor to implement the transverse jet flow mixing characteristic prediction method according to any one of the above.

[0033] The application further provides a computer program product, comprising a computer program, wherein the computer program is executed by a processor to implement the transverse jet flow mixing characteristic prediction method according to any one of the above.

[0034] The transverse jet flow mixing characteristic prediction method and device provided by the application can realize more efficient, fast and relatively accurate prediction of the mixing characteristics of small flow transverse jet flow, and can provide technical support for the design of the mixing area of an advanced low-mixing air combustion chamber. BRIEF DESCRIPTION OF DRAWINGS

[0035] In order to more clearly illustrate the technical solutions of the present application or prior art, the following will briefly introduce the drawings needed to be used in the embodiments or prior art description. Obviously, the drawings in the following description are some embodiments of the present application, and for those skilled in the art, other drawings can also be obtained without creative labor.

[0036] Figure 1 is one of the flow schematic diagrams of the transverse jet flow mixing characteristic prediction method provided by the application;

[0037] Figure 2 is a calculation area schematic diagram of the transverse jet flow mixing characteristic prediction method provided by the application;

[0038] Figure 3 is a flow schematic diagram of the transverse jet flow mixing characteristic prediction method provided by the application;

[0039] Figure 4 is a schematic diagram of the transverse jet flow mixing characteristic prediction method provided by the application, wherein the transverse jet flow is composed of a finite number of micro-element segments;

[0040] Figure 5 is another flow schematic diagram of the transverse jet flow mixing characteristic prediction method provided by the application;

[0041] Figure 6 is a schematic diagram of the trajectory of the transverse jet in the three-dimensional space in the transverse jet mixing characteristic prediction method provided by the present application;

[0042] Figure 7 is a structural schematic diagram of the transverse jet mixing characteristic prediction device provided by the present application.

[0043] Figure 8 is a structural schematic diagram of the electronic device provided by the present application. DETAILED DESCRIPTION

[0044] To make the objectives, technical solutions and advantages of the present application clearer, the technical solutions in the present application will be described below in conjunction with the accompanying drawings in the present application. Obviously, the described embodiments are only some of the embodiments of the present application, rather than all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those of ordinary skill in the art without creative work fall within the protection scope of the present application.

[0045] In the description of the embodiments of the present application, the terms “first”, “second”, “third” are only for the purpose of description, and cannot be understood as indicating or implying relative importance, and do not involve sequence.

[0046] The transverse jet mixing characteristic prediction method and device provided by the present application will be described below in conjunction with Figures 1 to 8

[0047] Figure 1 is one of the flow schematic diagrams of the transverse jet mixing characteristic prediction method provided by the present application. As shown in Figure 1 the execution subject of the transverse jet mixing characteristic prediction method provided by the embodiment of the present application can be a transverse jet mixing characteristic prediction device, and the method comprises steps 101, 102 and 103.

[0048] For the application scene of the advanced gas turbine low-mixing air combustion chamber, important features are extracted in combination with the engineering actual situation to obtain a simplified physical model. The model represents the flow structure of the porous jet and the mixing of cold and hot air, and the main flow inflow condition is swirl.

[0049] In the embodiment of the present application, the main flow inflow condition is a non-uniform temperature field and a swirl flow field.

[0050] In the embodiment of the present application, the transverse jet (which can be simply referred to as “jet”) is a small-flow transverse jet. The small-flow transverse jet refers to the transverse jet as a mixing gas, and the mixing gas ratio is less than a ratio threshold. Optionally, the ratio threshold can be 5% to 20%.

[0051] ​Step 101, obtaining the physical parameters of the main flow, the initial physical parameters and the initial geometric parameters of the lateral jet flow, and the position of the target estimated plane.

[0052] Specifically, as shown in Figure 2 and Figure 3 , the physical parameters of the main flow can include temperature, velocity and density. The temperature field and velocity field of the main flow are both non-uniform fields. Among them, the velocity fields of the main flow in three directions in the rectangular coordinate system are U x , U y , U z , and can be given by a fitting function; the temperature field of the main flow is T m ; and the density of the main flow is p m . The three directions in the rectangular coordinate system are the axial direction (X-axis direction), the radial direction (Y-axis direction) and the spanwise direction (Z-axis direction) respectively.

[0053] The physical model of the advanced gas turbine main combustion chamber is simplified to obtain the calculation region under the condition of the main flow rotational flow inflow as shown in Figure 2 . The model contains the following assumptions: ① assuming that the control body wall is an adiabatic wall; ② assuming that the control body inlet plane is a square, i.e. the channel width W is equal to the channel height H; ③ assuming that the rotational flow center of the main flow is the geometric center of the inlet plane; ④ assuming that the origin of the coordinate system is the rotational flow center of the inlet plane, i.e. the coordinates of the point are (0, 0, 0); ⑤ assuming that the coordinates of the center of the circular jet hole are (x j , y j , z j ), and the position of the jet hole of the initial research model is taken as z j = 0, i.e. the center of the jet hole has no eccentricity relative to the rotational flow center in the spanwise direction; ⑥ assuming that the flow condition is: steady, no combustion, viscous, rotational flow inflow.

[0054] Determination of the non-uniform field of the main flow, the temperature field is a non-uniform data matrix T m , and the velocity type in three directions can be expressed by a fitting function as follows:

[0055] Tangential velocity: where the angular velocity Ω = a1x + Ω inlet , Ω inlet is the rotational flow angular velocity at the inlet, R is the vortex core range, and a1 is a constant;

[0056] Axial velocity: where U1 = a2x + U 1,inlet , U 1,inlet is the characteristic axial velocity at the inlet representing the non-uniformity of the velocity, U0 is the characteristic axial velocity representing r→∞, and a2 is a constant;

[0057] Radial velocity: vr (r) = Cr, where C is a known constant.

[0058] like Figure 4 As shown, in this embodiment of the invention, the trajectory shape of the three-dimensional transverse jet is simplified to a curved cone, and the trajectory flow is regarded as a steady flow composed of a series of non-interfering jet elements. That is, the transverse jet may include multi-step jet elements. The properties of the jet elements change over a discrete time step Δt.

[0059] A transverse jet is incident perpendicularly into the main stream from a jet orifice. The initial physical parameters of the transverse jet may include initial velocity, initial temperature, and initial density, which are u j T j and ρ j The initial geometric parameters of the transverse jet may include the radius b0 and thickness h0 of the jet element in step 0. The diameter of the jet orifice is D.

[0060] The initial conditions for the transverse jet are:

[0061] (v x ,v y ,v z )0=(0,u j ,0)

[0062] (b,h)0=(0.5D,0.5D)

[0063] (T,ρ)0=(T j ,ρ j )

[0064]

[0065] The position of the target prediction plane can be a preset position. Optionally, the position of the target prediction plane can be represented by the axial position Xout / H of the target prediction plane relative to the center of the jet hole.

[0066] Step 102: Based on the mainstream physical parameters, the initial physical parameters and initial geometric parameters of the transverse jet, obtain the position of the penetration center of the transverse jet in each intermediate estimated plane.

[0067] Specifically, for any jet element at any step, the characteristic values ​​of the jet element at that step, such as position, velocity, and temperature, can all be obtained recursively from the initial physical and geometric parameters of the transverse jet.

[0068] For the k-th jet element, the position of the penetration center is (x k ,y k ,z k ), speed is (v x,k ,v y,kv z,k , the velocity is The temperature is T k , the density is p k The angle between the jet axis and the horizontal plane is The angle between the projection of the jet axis on the horizontal plane and the X axis is k The radius of the k-th step jet element is b k The thickness is h k The mass is

[0069] Optionally, based on the Lagrange recursion method, the recursive relationship of the flow characteristic values of the adjacent two-step jet elements can be analyzed from the perspective of mass conservation, momentum conservation, energy conservation and state equation, and the recursive relationship of the geometric characteristic values of the adjacent two-step jet elements can be analyzed from the perspective of the geometric characteristics of the model, so as to obtain the characteristic values of the (k+1) -th step jet element.

[0070] It should be noted that the jet element corresponds to the intermediate estimation plane one by one, and the position of the penetration center of the k-th step jet element (x k , y k , z k ) is the position of the penetration center of the transverse jet in the k-th intermediate estimation plane.

[0071] In the case that the position of the penetration center of the transverse jet in a certain intermediate estimation plane reaches the position of the target estimation plane, the recursion can be stopped.

[0072] The position of the penetration center of the transverse jet in the k-th intermediate estimation plane can be represented by X k / H.

[0073] Step 103, based on the position of the penetration center of the transverse jet in each intermediate estimation plane, the trajectory of the transverse jet and the coordinate of the penetration center of the transverse jet in the target estimation plane are obtained.

[0074] Specifically, after obtaining the position of the penetration center of the transverse jet in each intermediate estimation plane, the mixing characteristics of the transverse jet can be integrated, that is, the obtained position of the penetration center of the transverse jet in each intermediate estimation plane is re-integrated to obtain the mixing characteristics of the transverse jet, including the trajectory of the transverse jet and the coordinate of the penetration center in the target estimation plane.

[0075] Optionally, the positions of the penetration centers of the transverse jets in the intermediate estimation planes can be integrated into a matrix.

[0076] Optionally, any data fitting method can be used to fit the positions of the penetration centers of the transverse jets in the intermediate estimation planes to obtain the trajectory of the transverse jet.

[0077] Optionally, based on the trajectory of the lateral jet and the position of the target estimation plane, the coordinate of the penetration center of the lateral jet in the target estimation plane can be obtained.

[0078] Optionally, after step 103, the mixing characteristics of the lateral jet can be visualized, the trajectory of the lateral jet is visualized and output based on the mixing characteristics of the lateral jet obtained in step 103, and the coordinate of the penetration center of the lateral jet in the target estimation plane is displayed.

[0079] Optionally, after step 103, the coordinate of the penetration center of the lateral jet in the target estimation plane can also be obtained, and the jet penetration depth can be obtained. Based on the coordinate of the penetration center of the lateral jet in the target estimation plane, the jet penetration depth can include: obtaining the vertical distance (i.e. the distance in the Y-axis direction) between the penetration center of the lateral jet in the target estimation plane and the center of the jet hole as the jet penetration depth.

[0080] The embodiment of the present application regards the trajectory flow of the lateral jet as a series of steady flow composed of mutually non-interfering jet elements, analyzes the force of any jet element, recursively obtains the trajectory of the lateral jet and the coordinate of the penetration center of the lateral jet in the target estimation plane based on the Lagrange recursive algorithm, can more efficiently, quickly and accurately predict the mixing characteristics of small-flow lateral jet, and can provide technical support for the design of the mixing area of the advanced low-mixing air combustion chamber.

[0081] Based on the content of any of the above embodiments, based on the physical parameters of the main flow, the initial physical parameters and the initial geometric parameters of the lateral jet, the position of the penetration center of the lateral jet in each intermediate estimation plane is obtained, including: in the case that the position of the penetration center of the lateral jet in the kth intermediate estimation plane does not reach the position of the target estimation plane, based on the physical parameters of the main flow, the position of the penetration center of the lateral jet in the kth intermediate estimation plane, the physical parameters and the geometric parameters of the kth step jet element included in the lateral jet, the position of the penetration center of the lateral jet in the (k+1)th intermediate estimation plane is obtained.

[0082] Specifically, after obtaining the position of the penetration center of the kth intermediate estimation plane, it can be determined whether the position of the penetration center of the kth intermediate estimation plane reaches the position of the target estimation plane.

[0083] If (X k / H)≥(Xout / H), it has reached; if (X k / H)<(Xout / H), it has not reached.

[0084] If not, the position of the penetration center of the lateral jet in the (k+1)th intermediate estimation plane can be recursively calculated based on the position of the penetration center of the lateral jet in the kth intermediate estimation plane.

[0085] The position of the penetration center of the lateral jet in the (k+1)th intermediate estimation plane can be recursively calculated based on the position of the penetration center of the lateral jet in the kth intermediate estimation plane according to the physical law followed by the kth step jet element and the (k+1)th step jet element, based on the physical parameters of the main flow and the physical parameters and geometric parameters of the kth step jet element included in the lateral jet.

[0086] Optionally, a recursive model can be established according to the physical law followed by the kth step jet element and the (k+1)th step jet element. The recursive model is used to map the relationship between the position of the penetration center of the lateral jet in the kth intermediate estimation plane and the position of the penetration center of the lateral jet in the (k+1)th intermediate estimation plane.

[0087] The position of the penetration center of the lateral jet in the (k+1)th intermediate estimation plane can be obtained by inputting the position of the penetration center of the lateral jet in the kth intermediate estimation plane and the physical parameters and geometric parameters of the kth step jet element into the recursive model.

[0088] The physical parameters and geometric parameters of the kth step jet element together constitute the characteristic value of the aforementioned kth step jet element. The aforementioned physical parameters can include velocity, temperature and density. The aforementioned geometric parameters can include radius b k and thickness h k .

[0089] The thickness is proportional to the magnitude of the local jet velocity, and h k = V k Δt.

[0090] Optionally, according to the physical law followed by the kth step jet element and the (k+1)th step jet element, based on the physical parameters of the main flow, the position of the penetration center of the lateral jet in the kth intermediate estimation plane, the physical parameters and geometric parameters of the kth step jet element included in the lateral jet, the position of the penetration center of the lateral jet in the (k+1)th intermediate estimation plane is calculated.

[0091] The embodiment of the application regards the trajectory flow of the transverse jet as a series of steady flows composed of mutually non-interfering jet element groups, performs force analysis on the kth step jet element, performs recursion based on the Lagrange recursion algorithm, and obtains the position of the penetration center of the transverse jet in the (k+1)th intermediate estimation plane, so as to obtain the trajectory of the transverse jet and the coordinates of the penetration center of the transverse jet in the target estimation plane, and the mixing characteristics of the small-flow transverse jet can be more efficiently, quickly and accurately predicted, and technical support can be provided for the design of the mixing area of the advanced low-mixing air combustion chamber.

[0092] Based on the content of any of the above embodiments, based on the physical parameters of the main flow, the position of the penetration center of the transverse jet in the kth intermediate estimation plane, the physical parameters and geometric parameters of the kth step jet element included in the transverse jet, the position of the penetration center of the transverse jet in the (k+1)th intermediate estimation plane is obtained, including: based on the physical parameters of the main flow, the physical parameters and geometric parameters of the kth step jet element, obtaining the physical parameters and geometric parameters of the (k+1)th step jet element included in the transverse jet.

[0093] Specifically, the physical parameters and geometric parameters of the (k+1)th step jet element included in the transverse jet can be obtained based on the physical parameters of the main flow, the physical parameters and geometric parameters of the kth step jet element, according to the physical laws followed by the kth step jet element and the (k+1)th step jet element.

[0094] The temperature T of the (k+1)th step jet element k+1 satisfies the following energy conservation equation:

[0095]

[0096] The density of the (k+1)th step jet element satisfies the following state equation (assuming constant pressure):

[0097]

[0098] The thickness h of the (k+1)th step jet element k+1 and the radius b k+1 respectively satisfy the following recursive equations:

[0099]

[0100]

[0101] wherein V k+1 represents the resultant velocity of the (k+1)th step jet element.

[0102] The velocities v x,k+1 , v y,k+1 , vz,k+1 The following momentum conservation equation is satisfied:

[0103]

[0104]

[0105]

[0106] where ΔP x,k , ΔP y,k , and ΔP z,k represent the total momentum increment of the k-th step jet element in three directions, respectively; M k+1 represents the mass of the (k+1)-th step jet element. The above total momentum increment is increased due to the momentum introduced by turbulent entrainment.

[0107] The mass M k+1 of the (k+1)-th step jet element satisfies the following mass conservation equation:

[0108] M k+1 = M k + ΔM k

[0109]

[0110] where ΔM k represents the mass increment of the k-th step jet element caused by turbulent entrainment.

[0111] Based on the position of the penetration center of the lateral jet in the k-th intermediate estimation plane and the physical parameters of the (k+1)-th step jet element, the position of the penetration center of the lateral jet in the (k+1)-th intermediate estimation plane is obtained.

[0112] Specifically, the position (x k+1 , y k+1 , z k+1 ) of the penetration center of the lateral jet in the (k+1)-th intermediate estimation plane can be obtained by the following formula:

[0113] x k+1 = x k + v x,k+1 · Δt

[0114] y k+1 = y k + v y,k+1 · Δt

[0115] z k+1 = z k + v z,k+1 · Δt

[0116] Δs k+1= V k+1 Delta t

[0117] Wherein, Delta s k+1 Indicates the displacement of the (k+1)th jet element; Delta t can be based on the thickness of the kth jet element and y k Determination.

[0118] For the kth jet element in the strong turbulent high-temperature gas, the recursive relationship of the physical parameters of the adjacent two jet elements is analyzed from the aspects of mass conservation, momentum conservation, energy conservation and state equation, and the recursive relationship of the geometric parameters of the adjacent two jet elements is analyzed from the aspect of the geometric characteristics of the model.

[0119] Mass conservation is used to determine the mass increment Delta M k Of the kth jet element; momentum conservation is used to determine the momentum increment Delta P k Of the kth jet element, which comes from the momentum increase caused by turbulent entrainment, that is, the increment obtained by the momentum transport of the main flow; energy conservation is used to determine the energy increment Delta H k Of the kth jet element, which comes from the energy increase caused by turbulent entrainment, that is, the increment obtained by the energy transport of the main flow; state equation is used to establish the proportional relationship between the density value and the temperature value of the kth and (k+1)th jet elements; the recursive method of the geometric parameters includes the recursive method of the thickness determined by the proportional relationship between the velocity and the thickness of any jet element, the recursive method of the angle determined by the three direction velocity components and the resultant velocity of the (k+1)th jet element, and the recursive method of the center position determined by the displacement of the kth to (k+1)th jet elements.

[0120] The embodiment of the application regards the trajectory flow of the transverse jet as a series of jet element groups that do not interfere with each other, analyzes the force of any jet element, recursively calculates based on the Lagrange recursive algorithm, mathematically describes the physical process by listing the basic control equation group and obtaining the recursive formula group, establishes the small-flow transverse jet mixing model, and obtains the trajectory of the transverse jet and the coordinates of the penetration center of the transverse jet in the target estimation plane based on the model, so that the mixing characteristics of the small-flow transverse jet can be more efficiently, quickly and accurately predicted, and technical support can be provided for the design of the mixing area of the advanced low-mixing air combustion chamber.

[0121] Based on the content of any of the above embodiments, the physical parameters and geometric parameters of the (k+1)th jet element included in the transverse jet are obtained based on the physical parameters of the main flow, the physical parameters and geometric parameters of the kth jet element, including: based on the physical parameters of the main flow, the physical parameters and geometric parameters of the kth jet element, the mass of the (k+1)th jet element and the total momentum increment of the kth jet element are obtained.

[0122] Specifically, after obtaining the physical parameters of the main flow, the physical parameters and the geometric parameters of the k-th jet element, the mass of the (k+1)-th jet element can be obtained according to the aforementioned mass conservation equation.

[0123] The mass of the (k+1)-th jet element and the total momentum increment of the k-th jet element can be obtained based on the physical law according to the physical parameters of the main flow and the physical parameters of the k-th jet element.

[0124] The velocity of the (k+1)-th jet element is obtained based on the mass of the (k+1)-th jet element, the momentum increment of the k-th jet element and the velocity of the k-th jet element.

[0125] Specifically, after obtaining the mass of the (k+1)-th jet element, the momentum increment of the k-th jet element and the velocity of the k-th jet element, the velocity of the (k+1)-th jet element can be obtained according to the aforementioned momentum conservation equation.

[0126] The velocity of the k-th jet element and the velocity of the (k+1)-th jet element can include the aforementioned velocities in the three directions.

[0127] The embodiment of the present application can realize more efficient, fast and relatively accurate prediction of the mixing characteristics of the small-flow transverse jet by regarding the trajectory flow of the transverse jet as a series of steady flows composed of mutually non-interfering jet elements, performing force analysis on any jet element, performing recursion based on the Lagrange recursion algorithm, listing the basic control equation set and deriving the recursion formula set to mathematize the physical process, establishing a small-flow transverse jet mixing model, and obtaining the trajectory of the transverse jet and the coordinates of the penetration center of the transverse jet in the target estimation plane based on the model.

[0128] Based on the content of any of the above embodiments, the mass of the (k+1)-th jet element and the total momentum increment of the k-th jet element are obtained based on the physical parameters of the main flow, the physical parameters and the geometric parameters of the k-th jet element, including: based on the physical parameters of the main flow, the physical parameters and the geometric parameters of the k-th jet element, obtaining the first momentum increment, the first mass increment, the second momentum increment and the second mass increment of the k-th jet element; wherein the first momentum increment and the first mass increment are generated by the shear entrainment of the transverse jet; and the second momentum increment and the second mass increment are generated by the forced entrainment of the transverse jet.

[0129] Specifically, the mass increment of the jet element comes from the first mass increment ΔM s and the second mass increment ΔM f of the transverse jet. The first mass increment ΔM s, which is the mass of the main flow entrained into the jet element in the time range of Δt due to the velocity difference between the main flow and the jet element in the kth step, and its calculation contains an entrainment coefficient determined by an empirical formula; and a second mass increment ΔM f , which is the mass of all the main flow fluid on the windward surface of the jet element in the kth step entrained into the jet, and its calculation contains the projection value of the main flow velocity in the normal direction of the jet element surface obtained by vector calculation.

[0130] Shear entrainment, which is the shear effect due to the velocity difference between the main flow and the jet element in the kth step.

[0131] Forced entrainment, which is the mass of all the main flow fluid on the windward surface of the jet element entrained into the jet.

[0132] Alternatively, the shear entrainment amount (i.e., the first mass increment) and the forced entrainment amount (i.e., the second mass increment) are calculated respectively for the two types of entrainment, and the shear entrainment amount ΔM s and the forced entrainment amounts ΔM f,x , ΔM f,y , and ΔM f,z in three directions are determined according to the physical meaning of the two types of entrainment and geometric analysis, so as to obtain the momentum increments of the forced entrainment in three directions as ΔP x,k , ΔP y,k , and ΔP z,k . On this basis, the entrainment amount ΔM k of the jet element in the kth step is determined by the maximum assumption method.

[0133]

[0134]

[0135]

[0136] where Δρ k = ρ m - ρ k .

[0137]

[0138]

[0139]

[0140]

[0141]

[0142]

[0143]

[0144]

[0145]

[0146] wherein β x , β y , β z respectively represent the included angle between the main flow and the windward surface of the transverse jet in three directions; i, j, k respectively represent the direction vector of the velocity of the main flow; A represents the surface area of the transverse jet.

[0147] wherein a1, a2, a3, a4 are all coefficient values, and the calculation formula is:

[0148]

[0149]

[0150] a3 = R cos ψ

[0151] a4 = R sin ψ

[0152]

[0153]

[0154]

[0155] The calculation formula of the jet direction includes:

[0156]

[0157]

[0158]

[0159]

[0160]

[0161] Based on the mass of the kth step jet element, the first mass increment and the second mass increment, the mass of the (k+1)th step jet element is obtained, and based on the first momentum increment and the second momentum increment, the total momentum increment of the kth step jet element is obtained.

[0162] Specifically, the mass increment ΔM k may be determined by a weight selection method based on the size of the first mass increment ΔM s and the second mass increment ΔM f .

[0163] Optionally, the total momentum increment of the kth step jet flow element can be the sum of the first momentum increment and the second momentum increment.

[0164] The embodiment of the present application regards the trajectory flow of the transverse jet flow as a series of steady flows composed of mutually non-interfering jet flow elements, analyzes the force acting on any jet flow element, recursively calculates based on the Lagrange recursive algorithm, lists the basic control equation set and obtains the recursive formula set to mathematize the physical process, establishes the small-flow transverse jet flow mixing model, and obtains the trajectory of the transverse jet flow and the coordinates of the penetration center of the transverse jet flow in the target estimation plane based on the model, so that the mixing characteristics of the small-flow transverse jet flow can be more efficiently, quickly and accurately predicted, and technical support can be provided for the design of the mixing area of the advanced low-mixing air combustion chamber.

[0165] Based on the content of any of the above embodiments, the position of the penetration center of the transverse jet flow in the (k+1)th intermediate estimation plane is obtained based on the position of the penetration center of the transverse jet flow in the kth intermediate estimation plane and the physical parameters of the (k+1)th step jet flow element, including: obtaining the position of the penetration center of the transverse jet flow in the (k+1)th intermediate estimation plane based on the position of the penetration center of the transverse jet flow in the kth intermediate estimation plane and the velocity of the (k+1)th step jet flow element.

[0166] Specifically, the position of the penetration center of the transverse jet flow in the (k+1)th intermediate estimation plane (x k+1 ,y k+1 ,z k+1 ) can be obtained by the following formula:

[0167] x k+1 =x k +v x,k+1 ·Δt

[0168] y k+1 =y k +v y,k+1 ·Δt

[0169] z k+1 =z k +v z,k+1 ·Δt

[0170] Δs k+1 =V k+1 Δt

[0171] Wherein, Δs k+1 represents the displacement of the (k+1)th step jet flow element.

[0172] The embodiment of the application regards the trajectory flow of the transverse jet flow as a series of steady flows composed of mutually non-interfering jet flow elements, analyzes the force of any jet flow element, recursively calculates based on the Lagrange recursive algorithm, lists the basic control equation set and obtains the recursive formula set to mathematize the physical process, establishes the small-flow transverse jet flow mixing model, and obtains the trajectory of the transverse jet flow and the coordinates of the penetration center of the transverse jet flow in the target estimation plane based on the model, so that the mixing characteristics of the small-flow transverse jet flow can be more efficiently, quickly and accurately predicted, and technical support can be provided for the design of the mixing area of the advanced low-mixing air combustion chamber.

[0173] In order to facilitate the understanding of the above-mentioned embodiments of the application, the implementation process of the transverse jet flow mixing characteristic estimation method can be as shown in the flowchart of Fig. 1. Figure 5 Figure 6 The trajectory of the transverse jet flow obtained based on any of the above-mentioned embodiments is shown in Fig. 2.

[0174] In view of the small-flow transverse jet flow environment in the mixing area of the advanced low-mixing air combustion chamber, the non-uniform main flow inflow is abstracted as a non-uniform temperature field matrix and a swirl velocity field in three directions which can be mathematically expressed as a fitting function; the transverse jet flow is abstracted as a three-dimensional trajectory composed of a limited number of jet flow micro-element segments based on the Lagrange recursive algorithm analysis of the transverse jet flow, the force of any jet flow element is analyzed, and the recursive formula is written according to the conservation equation, so that the characteristic values of all jet flow elements in the entire jet flow trajectory can be recursively obtained from the initial jet flow element. The physical process is mathematically expressed, and the small-flow transverse jet flow mixing estimation model is realized through programming, so that the mixing situation of the transverse jet flow can be predicted in the application scenario of the low-mixing air combustion chamber, including two characteristic values of the transverse jet flow trajectory and the jet penetration depth.

[0175] The advantages of the above-mentioned embodiments of the application are as follows:

[0176] (1) The main flow characteristics in the engineering application scenario are extracted, and the velocity type of the swirl velocity field in three directions, i.e. the fitting function, is determined, effectively solving the problem that the main flow non-uniformity is not considered in the previous model;

[0177] (2) In view of the small-flow jet flow condition in the mixing area of the advanced low-mixing air combustion chamber, the momentum influence of the main flow swirl on the transverse jet flow in three directions is considered, the jet flow trajectory is regarded as a flow composed of a series of mutually non-interfering jet flow elements, the conservation equation is established through the force analysis of the jet flow elements, and a set of jet flow trajectory solving formula system under the main flow swirl condition is obtained, which makes up for the vacancy that the traditional method can only predict the jet flow trajectory under the uniform inflow;

[0178] ​(3) The low-flow transverse jet flow mixing prediction model can quickly and efficiently and relatively accurately predict the jet flow mixing characteristics, and the operation efficiency is improved by about two orders of magnitude compared with the numerical simulation calculation method, thereby providing support for the rapid prediction and design of the mixing zone flow in the low-air mixing combustion chamber.

[0179] The transverse jet flow mixing characteristic prediction device provided by the present application is described below, and the transverse jet flow mixing characteristic prediction device described below can be referred to in correspondence with the transverse jet flow mixing characteristic prediction method described above.

[0180] Figure 7 is a structural schematic diagram of the transverse jet flow mixing characteristic prediction device provided by the present application. Based on the content of any of the above embodiments, as shown in Figure 7 The device includes an acquisition module 701, a recursion module 702, and a prediction module 703.

[0181] The acquisition module 701 is configured to acquire the physical parameters of the main flow, the initial physical parameters and the initial geometric parameters of the transverse jet flow, and the position of the target prediction plane.

[0182] The recursion module 702 is configured to acquire the position of the penetration center of the transverse jet flow at each intermediate prediction plane based on the physical parameters of the main flow, the initial physical parameters and the initial geometric parameters of the transverse jet flow.

[0183] The prediction module 703 is configured to acquire the trajectory of the transverse jet flow and the coordinates of the penetration center of the transverse jet flow at the target prediction plane based on the positions of the penetration centers of the transverse jet flow at the intermediate prediction planes.

[0184] Specifically, the acquisition module 701, the recursion module 702, and the prediction module 703 can be sequentially electrically connected.

[0185] Optionally, the recursion module 702 can be specifically configured to, in a case where the position of the penetration center of the transverse jet flow at the kth intermediate prediction plane does not reach the position of the target prediction plane, acquire the position of the penetration center of the transverse jet flow at the (k+1)th intermediate prediction plane based on the physical parameters of the main flow, the position of the penetration center of the transverse jet flow at the kth intermediate prediction plane, the physical parameters and the geometric parameters of the kth jet element included in the transverse jet flow.

[0186] Optionally, the recursion module 702 can include:

[0187] The first acquisition unit is configured to acquire the physical parameters and the geometric parameters of the (k+1)th jet element included in the transverse jet flow based on the physical parameters of the main flow, the physical parameters and the geometric parameters of the kth jet element.

[0188] The second acquisition unit is configured to acquire the position of the penetration center of the transverse jet flow at the (k+1)th intermediate estimation plane based on the position of the penetration center of the transverse jet flow at the kth intermediate estimation plane and the physical parameter of the (k+1)th jet flow element.

[0189] Optionally, the first acquisition unit can include:

[0190] The first acquisition subunit is configured to acquire the mass of the (k+1)th jet flow element and the total momentum increment of the kth jet flow element based on the physical parameter of the main flow, the physical parameter and the geometric parameter of the kth jet flow element.

[0191] The second acquisition subunit is configured to acquire the velocity of the (k+1)th jet flow element based on the mass of the (k+1)th jet flow element, the momentum increment of the kth jet flow element and the velocity of the kth jet flow element.

[0192] Optionally, the first acquisition subunit can be specifically configured to:

[0193] acquire the first momentum increment, the first mass increment, the second momentum increment and the second mass increment of the kth jet flow element based on the physical parameter of the main flow, the physical parameter and the geometric parameter of the kth jet flow element;

[0194] acquire the mass of the (k+1)th jet flow element based on the mass, the first mass increment and the second mass increment of the kth jet flow element, and acquire the total momentum increment of the kth jet flow element based on the first momentum increment and the second momentum increment;

[0195] The first momentum increment and the first mass increment are generated by the shear entrainment of the transverse jet flow, and the second momentum increment and the second mass increment are generated by the forced entrainment of the transverse jet flow.

[0196] Optionally, the second acquisition unit can be specifically configured to acquire the position of the penetration center of the transverse jet flow at the (k+1)th intermediate estimation plane based on the position of the penetration center of the transverse jet flow at the kth intermediate estimation plane and the velocity of the (k+1)th jet flow element.

[0197] The transverse jet flow mixing characteristic estimation device provided by the embodiments of the present application is used for executing the transverse jet flow mixing characteristic estimation method provided by the present application, and has the same beneficial effects as the transverse jet flow mixing characteristic estimation method provided by the present application, and thus will not be described here.

[0198] The transverse jet flow mixing characteristic estimation device is used for the transverse jet flow mixing characteristic estimation method of the above-mentioned embodiments. Therefore, the description and definition in the transverse jet flow mixing characteristic estimation method in the above-mentioned embodiments can be used for understanding of each execution module in the embodiments of the present application.

[0199] The embodiment of the present application can realize more efficient, fast and relatively accurate prediction of the mixing characteristics of small flow transverse jet by regarding the trajectory flow of the transverse jet as a series of steady flows composed of mutually non-interfering jet element, performing force analysis on any jet element, and performing recursion based on Lagrange recursion algorithm to obtain the trajectory of the transverse jet and the coordinates of the penetration center of the transverse jet in the target estimation plane, and can provide technical support for the design of the mixing area of an advanced low-mixing air combustion chamber.

[0200] Figure 8 is a structural schematic diagram of an electronic device provided by the present application, as Figure 8 indicated, the electronic device can include a processor 810, a communications interface 820, a memory 830 and a communications bus 840, wherein the processor 810, the communications interface 820 and the memory 830 complete mutual communication through the communications bus 840. The processor 810 can invoke the logical instructions in the memory 830 to execute the transverse jet mixing characteristic estimation method, which includes: obtaining the physical parameters of the main flow, the initial physical parameters and the initial geometric parameters of the transverse jet, and the position of the target estimation plane; obtaining the position of the penetration center of the transverse jet in each intermediate estimation plane based on the physical parameters of the main flow, the initial physical parameters and the initial geometric parameters of the transverse jet; and obtaining the trajectory of the transverse jet and the coordinates of the penetration center of the transverse jet in the target estimation plane based on the position of the penetration center of the transverse jet in each intermediate estimation plane.

[0201] In addition, the logical instructions in the memory 830 described above can be implemented in the form of a software functional unit and sold or used as an independent product, which can be stored in a computer readable storage medium. Based on such understanding, the technical solutions of the present application essentially or the part that contributes to the prior art or part of the technical solutions can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including a plurality of instructions to make a computer device (which can be a personal computer, a server, or a network device, etc.) execute all or part of the steps of the embodiments of the present application. The foregoing storage medium includes: a U disk, a mobile hard disk, a read-only memory (ROM, Read-Only Memory), a random access memory (RAM, Random Access Memory), a magnetic disk or an optical disk, and various media that can store program codes.

[0202] The processor 810 in the electronic device provided in the embodiments of the present application can invoke the logical instructions in the memory 830, the implementation of which is consistent with the implementation of the method for estimating the mixing characteristics of the lateral jet provided in the present application, and the same beneficial effects can be achieved, and thus will not be described herein again.

[0203] In another aspect, the present application also provides a computer program product, which comprises a computer program stored on a non-transitory computer-readable storage medium, and the computer program comprises program instructions, when the program instructions are executed by a computer, the computer can execute the method for estimating the mixing characteristics of the lateral jet provided in the above-mentioned methods, and the method comprises: obtaining the physical parameters of the main flow, the initial physical parameters and the initial geometric parameters of the lateral jet, and the position of the target estimation plane; obtaining the position of the penetration center of the lateral jet at each intermediate estimation plane based on the physical parameters of the main flow, the initial physical parameters and the initial geometric parameters of the lateral jet; and obtaining the trajectory of the lateral jet and the coordinates of the penetration center of the lateral jet at the target estimation plane based on the position of the penetration center of the lateral jet at each intermediate estimation plane.

[0204] The computer program product provided in the embodiments of the present application is executed to implement the method for estimating the mixing characteristics of the lateral jet, the specific implementation is consistent with the implementation described in the embodiments of the aforementioned method, and the same beneficial effects can be achieved, and thus will not be described herein again.

[0205] In another aspect, the present application also provides a non-transitory computer-readable storage medium, which stores a computer program, and the computer program is executed by a processor to implement the method for estimating the mixing characteristics of the lateral jet provided in the above-mentioned methods, and the method comprises: obtaining the physical parameters of the main flow, the initial physical parameters and the initial geometric parameters of the lateral jet, and the position of the target estimation plane; obtaining the position of the penetration center of the lateral jet at each intermediate estimation plane based on the physical parameters of the main flow, the initial physical parameters and the initial geometric parameters of the lateral jet; and obtaining the trajectory of the lateral jet and the coordinates of the penetration center of the lateral jet at the target estimation plane based on the position of the penetration center of the lateral jet at each intermediate estimation plane.

[0206] The computer program stored on the non-transitory computer-readable storage medium provided in the embodiments of the present application is executed to implement the method for estimating the mixing characteristics of the lateral jet, the specific implementation is consistent with the implementation described in the embodiments of the aforementioned method, and the same beneficial effects can be achieved, and thus will not be described herein again.

[0207] The device embodiments described above are merely illustrative, wherein the units illustrated as separate components can or can not be physically separated, and the components illustrated as units can or can not be physical units, i.e., can be located in one place, or can be distributed to multiple network units. Part or all of the modules can be selected to achieve the purposes of the embodiments according to actual needs. Those skilled in the art can understand and implement without creative labor.

[0208] Through the description of the above embodiments, those skilled in the art can clearly understand that the embodiments can be realized by means of software and the necessary general hardware platform, and of course can also be realized by hardware. Based on such understanding, the above technical solutions can be embodied in the form of software products, and the computer software products can be stored in a computer readable storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., and include a plurality of instructions to make a computer device (which can be a personal computer, a server, or a network device, etc.) execute the methods of the embodiments or some parts of the embodiments.

[0209] Finally, it should be noted that: the above embodiments are only used to illustrate the technical solutions of the present application, and not to limit them; although the present application has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that: it can still modify the technical solutions recorded in the foregoing embodiments, or make equivalent replacement for part of the technical features; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present application.

Claims

1. A method for predicting the characteristics of a cross-flow mixing, characterized by, The method comprises the following steps: acquiring physical parameters of a main flow, initial physical parameters and initial geometric parameters of a lateral jet flow, and a position of a target estimation plane; wherein the physical parameters of the main flow comprise temperature, velocity and density; the initial physical parameters of the lateral jet flow comprise initial velocity, initial temperature and initial density; the initial geometric parameters comprise radius and thickness of a 0th step jet flow element; and the position of the target estimation plane is an axial position of the target estimation plane relative to a center of a jet flow hole; based on the physical parameters of the main flow, the initial physical parameters and the initial geometric parameters of the lateral jet flow, acquiring a position of a penetration center of the lateral jet flow at each intermediate estimation plane; wherein the position of the penetration center is a position of a penetration center of the lateral jet flow at each intermediate estimation plane; based on the positions of the penetration centers of the lateral jet flow at the intermediate estimation planes, acquiring a trajectory of the lateral jet flow and a coordinate of a penetration center of the lateral jet flow at the target estimation plane.

2. The cross-flow mixing characteristic estimation method according to claim 1, characterized by, The step of acquiring the position of the penetration center of the lateral jet flow at each intermediate estimation plane based on the physical parameters of the main flow, the initial physical parameters and the initial geometric parameters of the lateral jet flow comprises the following steps: in a case where the position of the penetration center of the lateral jet flow at the kth intermediate estimation plane does not reach the position of the target estimation plane, acquiring the position of the penetration center of the lateral jet flow at the (k+1)th intermediate estimation plane based on the physical parameters of the main flow, the position of the penetration center of the lateral jet flow at the kth intermediate estimation plane, physical parameters and geometric parameters of a kth step jet flow element included in the lateral jet flow.

3. The cross-flow mixing characteristic estimation method according to claim 2, characterized by, The step of acquiring the position of the penetration center of the lateral jet flow at the (k+1)th intermediate estimation plane based on the physical parameters of the main flow, the position of the penetration center of the lateral jet flow at the kth intermediate estimation plane, the physical parameters and the geometric parameters of the kth step jet flow element included in the lateral jet flow comprises the following steps: acquiring physical parameters and geometric parameters of a (k+1)th step jet flow element included in the lateral jet flow based on the physical parameters of the main flow, the physical parameters and the geometric parameters of the kth step jet flow element; acquiring the position of the penetration center of the lateral jet flow at the (k+1)th intermediate estimation plane based on the position of the penetration center of the lateral jet flow at the kth intermediate estimation plane and the physical parameters of the (k+1)th step jet flow element.

4. The cross-flow mixing characteristic estimation method according to claim 3, characterized by, The step of acquiring the physical parameters and the geometric parameters of the (k+1)th step jet flow element included in the lateral jet flow based on the physical parameters of the main flow, the physical parameters and the geometric parameters of the kth step jet flow element comprises the following steps: acquiring mass of the (k+1)th step jet flow element and total momentum increment of the kth step jet flow element based on the physical parameters of the main flow, the physical parameters and the geometric parameters of the kth step jet flow element; acquiring velocity of the (k+1)th step jet flow element based on the mass of the (k+1)th step jet flow element, the momentum increment of the kth step jet flow element and the velocity of the kth step jet flow element.

5. The cross-flow mixing characteristic estimation method according to claim 4, characterized by, The process of obtaining the mass of the (k+1)th step jet element and the total momentum increment of the kth step jet element based on the mainstream physical parameters, the physical parameters and geometric parameters of the kth step jet element includes: Based on the mainstream physical parameters, the physical parameters and geometric parameters of the k-th step jet element, the first momentum increment, the first mass increment, the second momentum increment and the second mass increment of the k-th step jet element are obtained. Based on the mass of the jet element at step k, the first mass increment, and the second mass increment, the mass of the jet element at step (k+1) is obtained, and based on the first momentum increment and the second momentum increment, the total momentum increment of the jet element at step k is obtained. Wherein, the first momentum increment and the first mass increment are generated by the shear entrainment of the transverse jet; the second momentum increment and the second mass increment are generated by the forced entrainment of the transverse jet.

6. The cross-flow mixing characteristic estimation method according to claim 4 or 5, characterized by, The method of obtaining the position of the penetration center of the transverse jet in the (k+1)th intermediate estimated plane based on the position of the penetration center of the transverse jet in the kth intermediate estimated plane and the physical parameters of the jet element in the (k+1)th step includes: Based on the position of the penetration center of the transverse jet in the kth intermediate estimated plane and the velocity of the jet element in the (k+1)th step, the position of the penetration center of the transverse jet in the (k+1)th intermediate estimated plane is obtained.

7. A lateral jet mixing characteristic prediction device characterized by comprising: include: The acquisition module is used to acquire the physical parameters of the main stream, the initial physical parameters and initial geometric parameters of the transverse jet, and the position of the target prediction plane; wherein, the physical parameters of the main stream include temperature, velocity, and density; the initial physical parameters of the transverse jet include initial velocity, initial temperature, and initial density; the initial geometric parameters include the radius and thickness of the jet element in step 0; and the position of the target prediction plane is the axial position of the target prediction plane relative to the center of the jet orifice. The recursive module is used to obtain the position of the penetration center of the transverse jet in each intermediate estimated plane based on the physical parameters of the main stream, the initial physical parameters and the initial geometric parameters of the transverse jet; wherein, the position of the penetration center is the position of the penetration center of the transverse jet in each intermediate estimated plane. The estimation module is used to obtain the trajectory of the transverse jet and the coordinates of the transverse jet's penetration center on the target estimation plane based on the position of the transverse jet's penetration center on each of the intermediate estimation planes.

8. An electronic device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized in that, When the processor executes the program, it implements the method for predicting the mixing characteristics of transverse jets as described in any one of claims 1 to 6. 9.A non-transitory computer-readable storage medium having stored thereon a computer program, characterized in that, When the computer program is executed by the processor, it implements the method for predicting the mixing characteristics of transverse jets as described in any one of claims 1 to 6.

10. A computer program product comprising a computer program, characterized in that, When the computer program is executed by the processor, it implements the method for predicting the mixing characteristics of transverse jets as described in any one of claims 1 to 6.

Images