A method, device and terminal device for flow field analysis of an afterburner

Abstract

This application is applicable to the field of numerical simulation technology, and provides a method, device and terminal device for flow field analysis of an afterburner. In the embodiments of this application, a calculation region of the afterburner is obtained, the above calculation region is meshed to obtain a plurality of target meshes; a first flow state of the fluid in the above calculation region is determined; based on preset boundary conditions, the above first flow state is discretized to obtain a second flow state of the above fluid in each of the above target meshes; the flow field distribution information of the above calculation region is determined according to a plurality of the above second flow states, thereby improving the efficiency of determining the flow field distribution of the afterburner.

Description

Technical Field

[0001] This application belongs to the technical field of numerical simulation, and particularly relates to a method, device, and terminal device for analyzing the flow field of an afterburner. Background Art

[0002] The combustion chamber is an important part of an aeroengine. Different from ordinary combustion chambers, the afterburner is mainly used to meet the demand for the engine to increase thrust in a short period, thereby greatly increasing the maneuverability of the aircraft. Generally, the design rationality of the afterburner is verified or the afterburner is optimized by determining the indoor flow field distribution of the afterburner.

[0003] However, with the improvement of engine performance, the working environment of the afterburner becomes more and more complex. The existing method of studying the flow field distribution through experiments will result in a large time cost, and thus the efficiency of determining the flow field distribution of the afterburner is relatively low. Summary of the Invention

[0004] The embodiments of this application provide a method, device, and terminal device for analyzing the flow field of an afterburner, which can solve the problem of relatively low efficiency in determining the flow field distribution of the afterburner.

[0005] In a first aspect, the embodiments of this application provide a method for analyzing the flow field of an afterburner, including:

[0006] Obtain the calculation area of the afterburner, perform grid division on the above calculation area to obtain a plurality of target grids;

[0007] Determine the first flow state of the fluid in the above calculation area;

[0008] Based on the preset boundary conditions, perform discretization processing on the above first flow state to obtain the second flow state of the fluid in each of the above target grids;

[0009] Determine the flow field distribution information of the above calculation area according to a plurality of the above second flow states.

[0010] In a second aspect, the embodiments of this application provide a device for analyzing the flow field of an afterburner, including:

[0011] A grid division module, configured to obtain the calculation area of the afterburner, perform grid division on the above calculation area to obtain a plurality of target grids;

[0012] A state determination module, configured to determine the first flow state of the fluid in the above calculation area;

[0013] A discrete module, configured to discretize the first flow state based on preset boundary conditions to obtain the second flow state of the fluid in each of the target grids;

[0014] A distribution information determination module, configured to determine the flow field distribution information of the calculation region according to a plurality of the second flow states.

[0015] In a third aspect, an embodiment of the present application provides a terminal device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the steps of any one of the flow field analysis methods of the afterburner are implemented.

[0016] In a fourth aspect, an embodiment of the present application provides a computer-readable storage medium. The computer-readable storage medium stores a computer program. When the computer program is executed by a processor, the steps of any one of the flow field analysis methods of the afterburner are implemented.

[0017] In a fifth aspect, an embodiment of the present application provides a computer program product. When the computer program product runs on a terminal device, the terminal device is enabled to execute any one of the flow field analysis methods in the first aspect.

[0018] In the embodiment of the present application, the calculation region of the afterburner, that is, the flow field region of the afterburner, is obtained, and then the calculation region is meshed to obtain a plurality of target grids, so as to simulate and calculate the flow fields corresponding to each grid respectively, and then determine the first flow state of the fluid in the calculation region, discretize the first flow state based on preset boundary conditions to obtain the second flow state of the fluid in each of the target grids, and finally determine the flow field distribution information of the calculation region according to a plurality of the second flow states, thereby improving the efficiency of determining the flow field distribution of the afterburner by simulating and calculating the flow field region of the afterburner. Description of the Drawings

[0019] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the following will briefly introduce the drawings required for use in the embodiments or the description of the prior art. Obviously, the following drawings are only some embodiments of the present application. For those of ordinary skill in the art, other drawings can be obtained based on these drawings without creative efforts.

[0020] Figure 1 It is a schematic flowchart of the flow field analysis method of the afterburner provided by the embodiment of the present application;

[0021] Figure 2It is a schematic diagram of the calculation area where struts exist in the afterburner provided by the embodiment of the present application;

[0022] Figure 3 It is a schematic diagram of the grid distribution where struts exist in the afterburner provided by the embodiment of the present application;

[0023] Figure 4 It is a schematic diagram of the grid distribution where no struts exist in the afterburner provided by the embodiment of the present application;

[0024] Figure 5 It is a schematic diagram of the flow field velocity where struts exist in the afterburner provided by the embodiment of the present application;

[0025] Figure 6 It is a schematic diagram of the flow field velocity where no struts exist in the afterburner provided by the embodiment of the present application;

[0026] Figure 7 It is a vector diagram of the velocity field of the cavity cross-section where struts exist in the afterburner provided by the embodiment of the present application;

[0027] Figure 8 It is a vector diagram of the velocity field of the cavity cross-section where no struts exist in the afterburner provided by the embodiment of the present application;

[0028] Figure 9 It is a schematic diagram of the result comparison provided by the embodiment of the present application;

[0029] Figure 10 It is a schematic diagram of the structure of the flow field analysis device of the afterburner provided by the embodiment of the present application;

[0030] Figure 11 It is a schematic diagram of the structure of the terminal device provided by the embodiment of the present application. Detailed implementation manners

[0031] In the following description, specific details such as specific system structures and technologies are presented for the purpose of illustration rather than limitation, so as to thoroughly understand the embodiments of the present application. However, those skilled in the art should clearly understand that the present application can also be implemented in other embodiments without these specific details. In other cases, detailed descriptions of well-known systems, devices, circuits, and methods are omitted to avoid unnecessary details from interfering with the description of the present application.

[0032] It should be understood that when used in the specification and appended claims of the present application, the term "comprising" indicates the presence of the described features, wholes, steps, operations, elements, and / or components, but does not exclude the presence or addition of one or more other features, wholes, steps, operations, elements, components, and / or their combinations.

[0033] It should also be understood that the term "and / or" as used in the specification and appended claims of this application refers to any combination and all possible combinations of one or more of the associated listed items, and includes such combinations.

[0034] As used in the specification and appended claims of this application, the term "if" can be interpreted as "when" or "once" or "in response to determining" or "in response to detecting" depending on the context. Similarly, the phrase "if determined" or "if [the described condition or event] is detected" can be interpreted as meaning "once determined" or "in response to determining" or "once [the described condition or event] is detected" or "in response to detecting [the described condition or event]" depending on the context.

[0035] In addition, in the description of the specification and appended claims of this application, the terms "first", "second", "third", etc. are only used for distinguishing descriptions and cannot be understood as indicating or implying relative importance.

[0036] Currently, to determine the flow field distribution of the afterburner, it is necessary to determine it through experimental methods, resulting in a low efficiency in determining the flow field distribution of the afterburner. This application provides a method for analyzing the flow field of the afterburner, which can quickly determine the flow field distribution of the afterburner through numerical simulation, thereby improving the efficiency of determining the flow field distribution of the afterburner.

[0037] Figure 1 The following shows a schematic flowchart of a method for analyzing the flow field of an afterburner in an embodiment of this application. The execution subject of this method can be a terminal device, such as Figure 1 As shown, the above-mentioned method for analyzing the flow field of the afterburner may include the following steps:

[0038] Step S101: Obtain the calculation region of the afterburner, and perform mesh division on the above-mentioned calculation region to obtain a plurality of target meshes.

[0039] In this embodiment, the terminal device needs to obtain the geometric shape of the afterburner corresponding to the experiment, and this geometric shape exists in the form of a volume, so as to determine the calculation region of the afterburner according to its geometric shape. The calculation region can be a calculation region including all regions of the afterburner, or a calculation region only including part of the regions of the afterburner. The above-mentioned divided target meshes can be tetrahedral meshes, hexahedral meshes, etc., and can be specifically set according to requirements.

[0040] Specifically, the calculation region of the afterburner determined according to the geometric shape can be to further determine the calculation region of the afterburner from the geometric shape by determining the characteristics of the geometric shape. Generally, since the characteristics of the afterburner include that the geometric shape is cylindrical and there are many struts in the afterburner that can divide the afterburner along the circumference, the terminal device determines the target direction based on this characteristic, that is, the direction of the struts in the above-mentioned afterburner, so as to obtain the above-mentioned calculation region by cutting the above-mentioned geometric shape in the target direction. The calculation region can be a calculation region with struts, such as Figure 2 shown Figure 2 The direction indicated by the black arrow in is the direction of gas flow in the afterburner. The airflow in the afterburner can move backward through the inner and outer bypass ducts. Among them, the faster the gas velocity in the outer bypass duct, the more fully the air can flow out, Figure 2 There is also a cavity in. Since deepening the cavity can accelerate the inlet velocity, after obtaining the flow field distribution information of the entire calculation region, the cavity shape can be designed based on the flow field distribution information corresponding to the cavity position. In addition, the above-mentioned calculation region can also be a calculation region without struts, that is, the mainstream region of the afterburner.

[0041] Specifically, the terminal device can use the open-source software OpenFOAM to divide the grid of the above-mentioned calculation region, so as to obtain a schematic diagram of the grid distribution of the calculation region including each target grid, such as Figure 3 and Figure 4 shown. It can be understood that by using the open-source software OpenFOAM to calculate the cold-state numerical value of the afterburner, the solution has high portability and scalability, avoiding the situation that commercial software cannot be well implanted into a large-scale server cluster when using commercial software. The calculation method can also be modified according to requirements, resulting in lower overall calculation costs and higher parallelism, and can also achieve high-precision numerical calculation of the cold-state flow field of the afterburner.

[0042] Step S102: Determine the first flow state of the fluid in the above-mentioned calculation region.

[0043] In this embodiment, the terminal device needs to construct a fluid control equation, which describes the flow state of the fluid in the calculation region, that is, the above-mentioned first flow state, so as to perform calculations according to the fluid control equation to obtain the flow field distribution information of the entire calculation region. Among them, the above-mentioned fluid control equation can be set according to user requirements and can be the Navier-Stokes equations (N-S equations).

[0044] Specifically, since the computational domain is an incompressible Newtonian fluid, the governing equations can adopt the N-S equations that satisfy incompressible fluids, and the formulas are as follows:

[0045]

[0046] Among them, v above is the velocity field; p above is the fluid pressure; v0 above is the kinematic viscosity, v0 = μ / ρ, μ above is the dynamic viscosity of the fluid; ρ above is the fluid density. It can be understood that the first equation in the above equations is the momentum conservation equation; the second equation in the equations is the continuity equation, that is, the mass conservation equation, which describes that the amount of fluid flowing out per unit volume should be equal to the amount of fluid flowing in.

[0047] In one embodiment, in order to facilitate the solution and simplification of the fluid governing equations describing the first flow state above, so as to obtain the flow field distribution information more accurately and quickly, the above step S102 may include: the terminal performs time-averaging processing on the preset fluid governing equations according to the Reynolds time-averaging method, that is, performs time-averaging processing on each flow field physical quantity in the governing equations in the time domain. The governing equations can be the N-S equations, so as to obtain the first governing equation, that is, the time-averaged governing equation, and the formulas are as follows:

[0048]

[0049] Among them, U above is the average velocity; P above is the average pressure; τ above is the turbulent stress term, and τ above needs to satisfy u above is the velocity pulsation term, and v = U + u. It can be understood that turbulence is a transient three-dimensional flow. In this embodiment, only the overall characteristics of the flow need to be considered. Therefore, the N-S equations above need to be time-averaged within a certain time, so as to obtain the time-averaged governing equation.

[0050] It can be understood that the turbulent stress term in the above first governing equation is composed of viscous stress and Reynolds stress. In most of the turbulent main body, the Reynolds stress is larger than the viscous stress. Therefore, the viscous stress can be ignored in many solution scenarios, so that the turbulent stress term in the above first governing equation can be determined according to the eddy viscosity coefficient method, resulting in the parameter to be solved being changed to the turbulent viscosity, and the specific formula is as follows:

[0051]

[0052] Among them, v T is the turbulent viscosity or the effective viscosity.

[0053] Further, the terminal device processes the first control equation according to a preset turbulence model to obtain a second control equation, thereby introducing a turbulence model to facilitate the modeling of τ and promoting more accurate calculation of the eddy current in the overall flow field. The turbulence model can be the turbulence kinetic energy - specific dissipation rate model of RANS (the turbulence kinetic energy - the specific rate of dissipation, k - ω model). In the k - ω model, the turbulence kinetic energy term k needs to satisfy and u rms 2 is proportional to the above - mentioned turbulence kinetic energy term k; the above - mentioned turbulence specific dissipation rate ω needs to satisfy where, the above - mentioned β * is a constant; by processing, the above - mentioned v T = k / ω.

[0054] And the control equation introducing the k - ω model is:

[0055]

[0056] where, the above - mentioned γ can take a value of 5 / 9; the above - mentioned β can take a value of 3 / 40; the above - mentioned β * can take a value of 9 / 100; the above - mentioned σ can take a value of 1 / 2; the above - mentioned σ * can take a value of 1 / 2.

[0057] Optionally, the above - mentioned turbulence model can also be the Reynolds - averaged Navier - Stokes (RANS) and Large Eddy Simulation (LES).

[0058] When the terminal device obtains the first control equation, the above - mentioned turbulence stress term, and the second control equation, it determines a target equation set according to the first control equation, the above - mentioned turbulence stress term, and the second control equation, that is, the equation set includes the first control equation, the representation equation corresponding to the above - mentioned turbulence stress term, and the second control equation, and the target equation set is used to describe the first flow state.

[0059] Step S103: Discretize the first flow state based on preset boundary conditions to obtain the second flow state of the fluid in each of the above - mentioned target grids.

[0060] In this embodiment, since the equations in the target equation set describing the first flow state are partial differential equations, in order to obtain the analytical solution or approximate analytical solution corresponding to each target grid, it is necessary to discretize the above target equation set based on preset boundary conditions, that is, to perform an integration on the partial differential equations in the above target equation set, and convert the continuous partial differential equation set and its definite solution conditions into an algebraic equation set on the discrete grid of the calculation region according to a preset method following specific rules, so as to obtain the discrete numerical approximation solution corresponding to each continuous target grid in the calculation region, and convert this discrete numerical approximation solution to obtain the algebraic equation set Ax = b corresponding to each target grid respectively, that is, to obtain the discrete equations corresponding to each target grid respectively for describing the second flow state. Moreover, since many physical quantities in the above equation set are in the form of derivatives with respect to x, y, and z, and the unit volume V = x * y * z, a triple integral of the volume is required when discretizing the target equation set.

[0061] Specifically, the terminal device can perform discretization processing through the finite volume method, that is, integrate the partial differential equations in the target equation set within the unit volume. In order to calculate the integral of the control volume, it is necessary to assume the variation law of the variable values between the grid points. The terminal device sets discrete formats such as second-order central difference to improve the stability and accuracy of the calculation during the discretization process. The above discretization processing of the above target equation set may include: discretizing the velocity gradient term in the above target equation set through a preset second-order central difference format, where, in the N-S equation, the above velocity gradient term is Discretizing the velocity convection term in the above target equation set through a preset second-order upwind format, where, in the N-S equation, the above velocity convection term is Discretizing the turbulent kinetic energy term in the above target equation set through a preset first-order upwind format, where, in the first control equation, the above turbulent kinetic energy term is τ; discretizing the Laplace term in the above target equation set through a preset second-order central format, and a restricted second-order central format can be further used for discretization, where, in the N-S equation, the above Laplace term is In addition, the interpolation method during the discretization process can adopt a second-order central discretization format.

[0062] Specifically, the terminal device determines the unknown quantities in the target equation set by presetting boundary conditions to achieve the discretization of the target equation set. The boundary conditions are the initial conditions. The inlet boundary condition in the above boundary conditions is:

[0063]

[0064] Among them, the above U inlet is the inlet velocity.

[0065] The outlet boundary condition in the above boundary conditions is as follows:

[0066]

[0067] Among them, the above U outlet is the outlet velocity.

[0068] Step S104: Determine the flow field distribution information of the above calculation region according to multiple second flow states.

[0069] In this embodiment, after the terminal device obtains multiple discrete equations for describing the above second flow state, the flow field distribution information of the calculation region corresponding to the afterburner can be obtained by calculating the multiple discrete equations. The flow field distribution information includes but is not limited to the velocity or pressure in the flow field, and can be displayed in the form of a flow field diagram, so as to design the shape of the combustion chamber or judge the rationality of the combustion chamber according to the flow field diagram. For example, Figure 5 and Figure 6 as shown, Figure 5 and Figure 6 show the velocity changes at each position in the flow field of the calculation region at a certain moment. Figure 5 and Figure 6 show the X-section of the calculation region, that is, the velocity changes corresponding to the Y-axis and Z-axis directions. Figure 5 and Figure 6 The values in each represent the velocity values at the corresponding positions.

[0070] In one embodiment, the terminal device can further display the flow field diagram corresponding to the local position in the calculation region according to user requirements, such as Figure 7 and Figure 8 , Figure 7 and Figure 8 which show the velocity schematic diagram of the cavity section, demonstrating the air flow velocity direction at the cavity section.

[0071] Furthermore, in order to verify whether the flow field distribution information obtained by the numerical simulation of the afterburner is accurate, a set of experimental results can be obtained by conducting experiments on the flow field and compared with the calculation results obtained by simulation. For example, Figure 9 as shown, Figure 9 what is compared is the cavity position with struts. Figure 9 The horizontal axis in is the velocity, and the vertical axis is the relative position with respect to the vortex center of the cavity, that is, the position of each position on the vertical axis with respect to the vortex center with the vortex center as the center point. Figure 9In this case, PIV represents the experimental results and CFD represents the calculation results. By comparing the calculation results with the experimental results, the error between the two can be obtained. Through the numerical simulation method in this embodiment, the error between the calculation results and the experimental results is relatively low. For example, the velocity error is only 7.8% and the distance error is only 0.5%. Therefore, the accuracy of the calculation results obtained through the numerical simulation method in this embodiment is relatively high.

[0072] In one embodiment, to improve the calculation speed, step S104 described above may include: the terminal device determines the processing cores corresponding to each of the above-mentioned target grids, and thus calculates the discrete equations of the multiple above-mentioned target grids through the processing cores corresponding to the multiple above-mentioned target grids respectively, that is, disperses the discrete equations of the multiple target grids to different processing cores, promotes the parallel processing of the multiple discrete equations, and obtains the corresponding calculation results simultaneously, so as to determine the flow field distribution information of the above-mentioned calculation area according to the calculation results of the multiple above-mentioned target grids.

[0073] In one embodiment, to improve the calculation efficiency, determining the processing cores corresponding to each of the above-mentioned target grids may include: obtaining the spatial position corresponding to the above-mentioned target grid or the calculation amount of the discrete equation, and thus determining the processing core corresponding to the above-mentioned target grid according to the above-mentioned spatial position or the above-mentioned calculation amount, that is, relatively more processing cores are used for calculation in the place where the grid is dense, and relatively fewer processing cores are used for calculation in the place where the grid is sparse; relatively more processing cores are used for calculating the discrete equation with a large calculation amount, and relatively fewer processing cores are used for calculating the discrete equation with a small calculation amount.

[0074] In one embodiment, when the processing core calculates the discrete equation, the multiple grid method can be used for solution, and when the preset convergence condition is reached, the calculation result corresponding to the target grid is determined.

[0075] In the embodiment of the present application, the calculation area of the afterburner is obtained, that is, the flow field area of the afterburner, and then the above-mentioned calculation area is meshed to obtain multiple target grids, so as to simulate and calculate the flow field corresponding to each grid respectively, and then determine the first flow state of the fluid in the above-mentioned calculation area, and discretize the above-mentioned first flow state based on the preset boundary conditions to obtain the second flow state of the above-mentioned fluid in each of the above-mentioned target grids respectively. Finally, the flow field distribution information of the above-mentioned calculation area is determined according to the multiple above-mentioned second flow states, so as to improve the efficiency of determining the flow field distribution of the afterburner by simulating and calculating the flow field area of the afterburner.

[0076] It should be understood that the magnitudes of the sequence numbers of the steps in the above embodiments do not mean the order of execution. The execution order of each process should be determined according to its function and internal logic, and should not constitute any limitation to the implementation process of the embodiments of the present application.

[0077] Corresponding to the flow field analysis method of an afterburner described above, Figure 10 The following is a schematic structural diagram of a flow field analysis device for an afterburner in an embodiment of the present application. As Figure 10 shown, the flow field analysis device for the above-mentioned afterburner may include:

[0078] A mesh generation module 1001, configured to obtain a calculation region of the afterburner, perform mesh generation on the above-mentioned calculation region, and obtain a plurality of target meshes.

[0079] A state determination module 1002, configured to determine a first flow state of the fluid in the above-mentioned calculation region.

[0080] A discretization module 1003, configured to perform discretization processing on the above-mentioned first flow state based on a preset boundary condition, and obtain a second flow state of the above-mentioned fluid in each of the above-mentioned target meshes.

[0081] A distribution information determination module 1004, configured to determine flow field distribution information of the above-mentioned calculation region according to a plurality of the above-mentioned second flow states.

[0082] In one embodiment, the above-mentioned mesh generation module 1001 may include:

[0083] A shape acquisition sub-module, configured to acquire the geometric shape of the above-mentioned afterburner.

[0084] A cutting sub-module, configured to cut the above-mentioned geometric shape in a target direction to obtain the above-mentioned calculation region, where the above-mentioned target direction is the direction of the strut in the above-mentioned afterburner.

[0085] In one embodiment, the above-mentioned state determination module 1002 may include:

[0086] A time-averaging processing sub-module, configured to perform time-averaging processing on a preset fluid control equation according to the Reynolds time-averaging method to obtain a first control equation.

[0087] A stress term determination sub-module, configured to determine a turbulent stress term in the above-mentioned first control equation according to the eddy viscosity coefficient method.

[0088] A processing sub-module, configured to process the above-mentioned first control equation according to a preset turbulence model to obtain a second control equation.

[0089] An equation set determination sub-module, configured to determine a target equation set according to the above-mentioned first control equation, the above-mentioned turbulent stress term, and the above-mentioned second control equation, where the above-mentioned target equation set is used to describe the above-mentioned first flow state.

[0090] In one embodiment, the above-mentioned discretization module 1003 may include:

[0091] A discrete sub-module for discretizing the above-mentioned target equations based on preset boundary conditions to obtain discrete equations corresponding to each of the above-mentioned target grids, where the discrete equations are used to describe the above-mentioned second flow state.

[0092] In one embodiment, the above-mentioned discrete sub-module may include:

[0093] A first discrete unit for discretizing the velocity gradient term in the above-mentioned target equations through a preset second-order central difference scheme.

[0094] A second discrete unit for discretizing the velocity convection term in the above-mentioned target equations through a preset second-order upwind scheme.

[0095] A third discrete unit for discretizing the turbulent kinetic energy term in the above-mentioned target equations through a preset first-order upwind scheme.

[0096] A fourth discrete unit for discretizing the Laplacian term in the above-mentioned target equations through a preset second-order central scheme.

[0097] In one embodiment, the above-mentioned distribution information determination module 1004 may include:

[0098] A processing core determination sub-module for determining the processing core corresponding to each of the above-mentioned target grids.

[0099] An equation calculation sub-module for calculating the discrete equations of multiple above-mentioned target grids respectively through the processing cores corresponding to the multiple above-mentioned target grids.

[0100] A distribution information determination sub-module for determining the flow field distribution information of the above-mentioned calculation region according to the calculation results of the multiple above-mentioned target grids.

[0101] In one embodiment, the above-mentioned processing core determination sub-module may include:

[0102] An information acquisition unit for acquiring the spatial position corresponding to the above-mentioned target grid or the computational amount of the discrete equation.

[0103] A processing core determination unit for determining the processing core corresponding to the above-mentioned target grid according to the above-mentioned spatial position or the above-mentioned computational amount.

[0104] In the embodiment of the present application, the calculation area of the afterburner, that is, the flow field area of the afterburner, is obtained, and then the above calculation area is meshed to obtain a plurality of target meshes, so as to facilitate the simulation calculation of the flow field corresponding to each mesh respectively. Then, the first flow state of the fluid in the above calculation area is determined, and the above first flow state is discretized based on the preset boundary conditions to obtain the second flow state of the above fluid in each of the above target meshes. Finally, the flow field distribution information of the above calculation area is determined according to a plurality of the above second flow states, thereby improving the efficiency of determining the flow field distribution of the afterburner by simulating the flow field area of the afterburner.

[0105] Those skilled in the art can clearly understand that for the convenience and simplicity of description, the specific working processes of the above-described devices and modules can refer to the corresponding processes in the foregoing system embodiments and method embodiments, and will not be described herein again.

[0106] Figure 11 It is a schematic structural diagram of the terminal device provided by the embodiment of the present application. For the convenience of description, only the parts related to the embodiment of the present application are shown.

[0107] As Figure 11 shown, the terminal device 11 of this embodiment includes: at least one processor 1100 ( Figure 11 only one is shown in the figure), a memory 1101 connected to the above processor 1100, and a computer program 1102 stored in the above memory 1101 and operable on the above at least one processor 1100, such as a flow field analysis program of the afterburner. When the above processor 1100 executes the above computer program 1102, the steps in the above method embodiments of the flow field analysis of each afterburner are implemented, such as Figure 1 the steps S101 to S104 shown. Or, when the above processor 1100 executes the above computer program 1102, the functions of each module in the above device embodiments are implemented, such as Figure 10 the functions of the modules 1001 to 1004 shown.

[0108] Exemplarily, the above computer program 1102 can be divided into one or more modules. The above one or more modules are stored in the above memory 1101 and executed by the above processor 1100 to complete the present application. The above one or more modules can be a series of computer program instruction segments capable of performing specific functions, and the instruction segments are used to describe the execution process of the above computer program 1102 in the above terminal device 11. For example, the above computer program 1102 can be divided into a mesh generation module 1001, a state determination module 1002, a discretization module 1003, and a distribution information determination module 1004. The specific functions of each module are as follows:

[0109] The mesh generation module 1001 is configured to obtain the calculation region of the afterburner, generate a mesh for the above-mentioned calculation region, and obtain a plurality of target meshes;

[0110] The state determination module 1002 is configured to determine the first flow state of the fluid in the above-mentioned calculation region;

[0111] The discretization module 1003 is configured to perform a discretization process on the above-mentioned first flow state based on a preset boundary condition, and obtain the second flow state of the above-mentioned fluid in each of the above-mentioned target meshes;

[0112] The distribution information determination module 1004 is configured to determine the flow field distribution information of the above-mentioned calculation region according to a plurality of the above-mentioned second flow states.

[0113] The above-mentioned terminal device 11 may include, but is not limited to, a processor 1100 and a memory 1101. Those skilled in the art can understand that Figure 11 This is only an example of the terminal device 11, and does not constitute a limitation on the terminal device 11. It may include more or fewer components than shown in the figure, or combine some components, or different components. For example, it may also include input / output devices, network access devices, buses, etc.

[0114] The so-called processor 1100 may be a central processing unit (CPU), and the processor 1100 may also be other general-purpose processors, digital signal processors (DSPs), application specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. The general-purpose processor may be a microprocessor or the processor may also be any conventional processor, etc.

[0115] In some embodiments, the above-mentioned memory 1101 may be an internal storage unit of the above-mentioned terminal device 11, such as the hard disk or memory of the terminal device 11. In some other embodiments, the above-mentioned memory 1101 may also be an external storage device of the above-mentioned terminal device 11, such as a plug-in hard disk equipped on the above-mentioned terminal device 11, a Smart Media Card (SMC), a Secure Digital (SD) card, a Flash Card, etc. Further, the above-mentioned memory 1101 may also include both the internal storage unit of the above-mentioned terminal device 11 and the external storage device. The above-mentioned memory 1101 is used to store an operating system, application programs, a Boot Loader, data, and other programs, such as the program code of the above-mentioned computer program, etc. The above-mentioned memory 1101 may also be used to temporarily store data that has been output or will be output.

[0116] Those skilled in the art can clearly understand that, for the convenience and simplicity of description, only the division of the above-mentioned functional units and modules is used as an example for illustration. In practical applications, the above-mentioned functions can be allocated to different functional units and modules according to needs, that is, the internal structure of the above-mentioned device can be divided into different functional units or modules to complete all or part of the functions described above. Each functional unit and module in the embodiments can be integrated into a processing unit, or each unit can exist physically alone, or two or more units can be integrated into one unit. The above-mentioned integrated unit can be implemented in the form of hardware or in the form of a software functional unit. In addition, the specific names of each functional unit and module are only for the convenience of mutual distinction and do not limit the protection scope of this application. The specific working processes of the units and modules in the above-mentioned system can refer to the corresponding processes in the foregoing method embodiments and will not be elaborated here. In the above-mentioned embodiments, each embodiment is described with emphasis. For the parts not detailed or recorded in a certain embodiment, reference can be made to the relevant descriptions of other embodiments.

[0117] Those of ordinary skill in the art can realize that the units and algorithm steps of each example described in combination with the embodiments disclosed herein can be implemented by electronic hardware, or by a combination of computer software and electronic hardware. Whether these functions are executed in hardware or software depends on the specific application and design constraints of the technical solution. Professional technicians can use different methods to implement the described functions for each specific application, but such implementation should not be considered to exceed the scope of this application.

[0118] In the embodiments provided in the present application, it should be understood that the disclosed device / terminal device and method can be implemented in other ways. For example, the device / terminal device embodiments described above are merely illustrative. For example, the above-mentioned division of modules or units is only a logical function division. In actual implementation, there may be other division methods. For example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the displayed or discussed coupling or direct coupling or communication connection between each other can be through some interfaces. The indirect coupling or communication connection of the device or unit can be in electrical, mechanical or other forms.

[0119] The units described above as separate components may or may not be physically separated. The components displayed as units may or may not be physical units, that is, they can be located in one place, or they can be distributed to multiple network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0120] If the above-mentioned integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on such an understanding, to implement all or part of the processes in the above-mentioned embodiment methods of the present application, a computer program can be used to instruct relevant hardware to complete. The above-mentioned computer program can be stored in a computer-readable storage medium. When the computer program is executed by a processor, the steps of the above-mentioned various method embodiments can be implemented. Among them, the above-mentioned computer program includes computer program code, and the above-mentioned computer program code can be in source code form, object code form, executable file or some intermediate form, etc. The above-mentioned computer-readable medium can at least include: any entity or device that can carry the computer program code to the photographing device / terminal device, recording medium, computer memory, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), electrical carrier signal, telecommunication signal, and software distribution medium. For example, a USB flash drive, a mobile hard disk, a magnetic disk or an optical disc, etc. In some jurisdictions, according to legislation and patent practice, the computer-readable medium cannot be an electrical carrier signal and a telecommunication signal.

[0121] The embodiments described above are only used to illustrate the technical solutions of the present application, rather than to limit it; although the present application has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that they can still modify the technical solutions recorded in the foregoing embodiments, or perform equivalent replacements on some 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, and should all be included in the protection scope of the present application.

Claims

1. A method for analyzing the flow field of an afterburner, characterized in that Including: Obtain the calculation region of the afterburner, perform mesh division on the calculation region to obtain a plurality of target meshes; Determine the first flow state of the fluid within the calculation region; Based on preset boundary conditions, perform discretization processing on the first flow state to obtain the second flow state of the fluid within each of the target meshes; Determine the flow field distribution information of the calculation region according to the plurality of second flow states; The determining the first flow state of the fluid within the calculation region includes: Perform time-averaging processing on a preset fluid control equation according to the Reynolds time-averaging method to obtain a first control equation; Determine the turbulent stress term in the first control equation according to the eddy viscosity coefficient method; Process the first control equation according to a preset turbulence model to obtain a second control equation; Determine a target equation set according to the first control equation, the turbulent stress term, and the second control equation, and the target equation set is used to describe the first flow state.

2. The flow field analysis method of the afterburner according to claim 1, characterized in that The obtaining the calculation region of the afterburner includes: Obtain the geometric shape of the afterburner; Cut the geometric shape in a target direction to obtain the calculation region, where the target direction is the direction of the strut in the afterburner.

3. The flow field analysis method of the afterburner according to claim 1, characterized in that, The performing discretization processing on the first flow state based on preset boundary conditions to obtain the second flow state of the fluid within each of the target meshes includes: Perform discretization processing on the target equation set based on preset boundary conditions to obtain discrete equations corresponding to each of the target meshes, and the discrete equations are used to describe the second flow state.

4. The flow field analysis method of the afterburner according to claim 3, characterized in that The performing discretization processing on the target equation set includes: Perform discretization processing on the velocity gradient term in the target equation set through a preset second-order central difference format; Perform discretization processing on the velocity convection term in the target equation set through a preset second-order upwind format; Perform discretization processing on the turbulent kinetic energy term in the target equation set through a preset first-order upwind format; Perform discretization processing on the Laplace term in the target equation set through a preset second-order central format.

5. The flow field analysis method of the afterburner according to claim 3, characterized in that The determining the flow field distribution information of the calculation region according to the plurality of second flow states includes: Determine the processing cores corresponding to each of the target meshes; Calculate the discrete equations of the plurality of target meshes respectively through the processing cores corresponding to the plurality of target meshes; Determine the flow field distribution information of the calculation region according to the calculation results of the plurality of target meshes.

6. The flow field analysis method of the afterburner according to claim 5, characterized in that The determining the processing cores corresponding to each of the target meshes includes: Obtain the spatial position corresponding to the target mesh or the calculation amount of the discrete equation; Determine the processing core corresponding to the target mesh according to the spatial position or the calculation amount.

7. A flow field analysis device for an afterburner, characterized in that, Including: A mesh division module, configured to obtain the calculation region of the afterburner, perform mesh division on the calculation region to obtain a plurality of target meshes; A state determination module, configured to determine the first flow state of the fluid within the calculation region; A discrete module for discretizing the first flow state based on preset boundary conditions to obtain the second flow state of the fluid in each of the target grids; A distribution information determination module for determining the flow field distribution information of the calculation region according to the plurality of second flow states; The state determination module is further configured to: Perform time averaging on a preset fluid control equation according to the Reynolds time-averaging method to obtain a first control equation; Determine the turbulent stress term in the first control equation according to the eddy viscosity coefficient method; Process the first control equation according to a preset turbulence model to obtain a second control equation; Determine a target equation set according to the first control equation, the turbulent stress term, and the second control equation, and the target equation set is used to describe the first flow state.

8. A terminal device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the computer program, the steps of a method for analyzing the flow field of an afterburner as described in any one of claims 1 to 6 are implemented.

9. A computer-readable storage medium storing a computer program, characterized in that, When the computer program is executed by a processor, the steps of a method for analyzing the flow field of an afterburner as described in any one of claims 1 to 6 are implemented.

Images