Vortex dynamics fluid value solving method, device, equipment and medium

Abstract

The invention discloses a vortex dynamics fluid value solving method, device and equipment and a medium, and relates to the technical field of fluid mechanics, and the method comprises the steps: constructing a corresponding initial vortex dynamics equation; performing numerical discretization on the initial vortex kinetic equation by using the initial velocity vector, the initial vortex vector and the equation boundary condition according to a preset numerical discretization format and a convection splitting algorithm to obtain a target vortex kinetic equation; performing updating calculation on the target vortex kinetic equation to obtain an updated vortex vector; updating a flow function control equation associated with the target vortex kinetic equation by using the updated vortex vector to obtain a target velocity vector determined based on the updated flow function vector; and taking the target velocity vector and the updated vortex vector as an initial velocity vector and an initial vortex vector, skipping to execute the step of obtaining a target vortex kinetic equation for updating the vortex vector at the next moment, and stopping iteration and outputting corresponding fluid numerical calculation results at different moments until a preset iteration condition is met.

Description

Technical Field

[0001] The present invention relates to the technical field of fluid mechanics, and in particular to a method, device, equipment and medium for numerically solving vortex dynamics fluid. Background Art

[0002] The traditional incompressible fluid solver is based on the Navier-Stokes equations, which requires large-scale matrix iteration to solve the pressure Poisson equation, which brings risks to the overall numerical robustness of the solver. Furthermore, for the complex geometric shapes encountered in the aviation industry, the complex vortex structure is the key to affecting its dynamic and thermodynamic processes. Using the vortex dynamics equation as the control equation to develop numerical software can more clearly characterize the key processes in industrial fluids and provide more powerful tools for analyzing the mechanism of complex physical processes and further optimizing designs.

[0003] The traditional incompressible fluid solver is based on the Navier-Stokes equations. It uses the step-by-step solution of the nonlinear convection term, linear viscous dissipation term, and pressure Poisson term of the NS equations, and uses the pressure-velocity correction method to obtain the velocity information of the flow field. This is the basic principle and process of the numerical simulation of the fluid solver. However, the pressure Poisson equation needs to be solved through large-scale iterative matrix calculation. When the matrix is ​​singular, the iterative solution often cannot obtain satisfactory results, resulting in the iterative matrix calculation process to cycle back and forth, causing the entire fluid solver to deadlock. This is a problem often encountered in industrial software. The existing solution is to sacrifice the accuracy of the flow field and add non-physical artificial numerical viscosity to the control equation to make the iterative matrix converge as soon as possible. However, this makes the numerical simulation results given by the solver deviate from the real physics to some extent, which is also a compromise caused by the bottleneck of traditional technology. Summary of the invention

[0004] In view of this, the purpose of the present invention is to provide a method, device, equipment and medium for numerical solution of vortex dynamics fluid, which can directly and accurately solve the numerical value of vortex dynamics fluid and avoid adding pressure and velocity correction in the solution process. The specific scheme is as follows:

[0005] In a first aspect, the present application discloses a method for numerically solving vortex dynamics fluid, comprising:

[0006] Constructing the corresponding initial vortex dynamics equation based on the fluid numerical solution requirements of the industrial fluid application scenario; wherein the industrial fluid application scenario is an industrial scenario related to a vortex system whose structure meets the preset complexity conditions;

[0007] According to a preset numerical discretization format and a convection splitting algorithm and using an initial velocity vector, an initial vortex vector and equation boundary conditions, each component in the initial vortex dynamics equation is numerically discretized to obtain a target vortex dynamics equation for updating the vortex vector at a next moment;

[0008] Performing an update calculation on the target vortex dynamics equation to obtain an updated vortex vector;

[0009] Using the updated vortex vector, the stream function control equation associated with the target vortex dynamics equation is updated to obtain a target velocity vector determined based on the updated stream function vector;

[0010] The target velocity vector and the updated vortex vector are used as the initial velocity vector and the initial vortex vector respectively, and the steps of numerically discretizing the components in the initial vortex dynamics equation according to the preset numerical discretization format and the convection splitting algorithm and using the initial velocity vector, the initial vortex vector and the boundary conditions of the equation are jumped to obtain the target vortex dynamics equation for updating the vortex vector at the next moment, until the preset iteration condition is met, the iteration is stopped and the corresponding fluid numerical calculation results at different moments are output.

[0011] Optionally, the fluid numerical solution requirement based on the industrial fluid application scenario is used to construct the corresponding initial vortex dynamics equation, including:

[0012] Based on the initial velocity vector, initial vortex vector, equation boundary conditions, and fluid kinematic viscosity in the fluid numerical solution requirements of industrial fluid application scenarios, the corresponding initial vortex dynamics equation is constructed; the initial vortex dynamics equation is:

[0013] ;

[0014] in, represents the initial vortex vector, represents the initial velocity vector, represents the fluid kinematic viscosity, Indicates time, represents the displacement vector, represents the velocity vector constrained by the solid wall, and the initial vortex vector at the solid wall The curl of the initial velocity vector constraint, represents the curl operator, represents the directions of the three coordinate axes in the three-dimensional space represented by the Einstein tensor, and The values ​​are 1, 2, and 3 respectively.

[0015] Optionally, the numerical discretization of each component in the initial vortex dynamics equation according to a preset numerical discretization format and a convection splitting algorithm and using an initial velocity vector, an initial vortex vector and equation boundary conditions to obtain a target vortex dynamics equation for updating the vortex vector at the next moment includes:

[0016] The first numerical discretization format of the Runge-Kutta third-order format coupled with the Crank-Nicolson second-order format is used to calculate the unsteady time term of the initial vortex dynamics equation based on the initial velocity vector, the initial vortex vector and the boundary conditions of the equation. Performing numerical discretization to obtain a first discrete representation;

[0017] The nonlinear convection term of the initial vortex dynamics equation is solved by using the second numerical discretization format of WENO5 coupled with the sixth-order compact difference and the convection splitting algorithm based on the initial velocity vector, the initial vortex vector and the boundary conditions of the equation. Performing numerical discretization to obtain a second discrete representation;

[0018] The dissipation term of the initial vortex dynamics equation is numerically discretized using the third numerical discretization format of the second-order central difference format and based on the initial velocity vector, the initial vortex vector and the boundary conditions of the equation. and vortex stretching term Performing numerical discretization to obtain a third discrete representation;

[0019] The curl term of the initial vortex dynamics equation is calculated based on the initial velocity vector, the initial vortex vector and the boundary conditions of the equation using the fourth numerical discretization format of the second-order central difference format. Performing numerical discretization to obtain a fourth discrete representation;

[0020] A target vortex dynamics equation for updating the vortex vector at the next moment is constructed based on the first discrete representation, the second discrete representation, the third discrete representation, and the fourth discrete representation.

[0021] Optionally, before the method of updating the stream function control equation associated with the target vortex dynamics equation using the updated vortex vector to obtain a target velocity vector determined based on the updated stream function vector, the method further includes:

[0022] The associated stream function control equation is constructed based on the target vortex dynamics equation, and the stream function control equation is:

[0023] ;

[0024] in, represents the stream function vector, and n represents the normal vector of the computational domain boundary.

[0025] Optionally, the updating of a stream function control equation associated with the target vortex dynamics equation using the updated vortex vector to obtain a target velocity vector determined based on the updated stream function vector includes:

[0026] The stream function control equation is solved based on the updated vortex vector using the fifth numerical discretization format of the second-order central difference format to obtain an updated stream function vector, so as to calculate the target velocity vector based on the updated stream function vector.

[0027] Optionally, the calculating the target velocity vector based on the updated stream function vector includes:

[0028] Substitute the updated stream function vector into , and obtain the target velocity vector.

[0029] Optionally, the output of the corresponding fluid numerical calculation results at different times includes:

[0030] The vortex vector and velocity vector at each moment calculated based on the input conditions of the initial velocity vector and the initial vortex vector at the current moment are output to obtain the flow field vortex vector information and flow field velocity vector information at each moment.

[0031] In a second aspect, the present application discloses a numerical solution device for vortex dynamics fluid, comprising:

[0032] An equation building module, used to build a corresponding initial vortex dynamics equation based on the fluid numerical solution requirements of an industrial fluid application scenario; wherein the industrial fluid application scenario is an industrial scenario related to a vortex system whose structure meets a preset complexity condition;

[0033] A numerical discretization module is used to perform numerical discretization on each component in the initial vortex dynamics equation according to a preset numerical discretization format and a convection splitting algorithm and using an initial velocity vector, an initial vortex vector and equation boundary conditions to obtain a target vortex dynamics equation for updating the vortex vector at the next moment;

[0034] An updating module, used for performing an updating calculation on the target vortex dynamics equation to obtain an updated vortex vector;

[0035] A secondary updating module, used for updating the stream function control equation associated with the target vortex dynamics equation using the updated vortex vector, so as to obtain a target velocity vector determined based on the updated stream function vector;

[0036] The result output module is used to use the target velocity vector and the updated vortex vector as the initial velocity vector and the initial vortex vector respectively, and jump to execute the steps of numerically discretizing the components in the initial vortex dynamics equation according to the preset numerical discretization format and the convection splitting algorithm and using the initial velocity vector, the initial vortex vector and the boundary conditions of the equation to obtain the target vortex dynamics equation for updating the vortex vector at the next moment, until the preset iteration condition is met, stop the iteration and output the corresponding fluid numerical calculation results at different moments.

[0037] In a third aspect, the present application discloses an electronic device, comprising:

[0038] Memory, used to store computer programs;

[0039] The processor is used to execute the computer program to implement the steps of the aforementioned disclosed method for numerical solution of vortex dynamics fluid.

[0040] In a fourth aspect, the present application discloses a computer-readable storage medium for storing a computer program; wherein, when the computer program is executed by a processor, the steps of the aforementioned disclosed method for numerical solution of vortex dynamics fluid are implemented.

[0041] It can be seen that the present application discloses a method for numerical solution of vortex dynamics fluid, including: constructing a corresponding initial vortex dynamics equation based on the fluid numerical solution requirements of an industrial fluid application scenario; wherein the industrial fluid application scenario is an industrial scenario related to a vortex system whose structure satisfies a preset complexity condition; numerically discretizing each component in the initial vortex dynamics equation according to a preset numerical discretization format and a convection splitting algorithm and using an initial velocity vector, an initial vortex vector and equation boundary conditions to obtain a target vortex dynamics equation for updating the vortex vector at the next moment; updating the target vortex dynamics equation to obtain an updated vortex vector; and using the updated vortex vector The stream function control equation associated with the target vortex dynamics equation is updated to obtain a target velocity vector determined based on the updated stream function vector; the target velocity vector and the updated vortex vector are used as the initial velocity vector and the initial vortex vector respectively, and the steps of numerically discretizing the components in the initial vortex dynamics equation according to the preset numerical discretization format and the convection splitting algorithm and using the initial velocity vector, the initial vortex vector and the boundary conditions of the equation are jumped to obtain the step of updating the target vortex dynamics equation for the vortex vector at the next moment, until the preset iteration condition is met, the iteration is stopped and the corresponding fluid numerical calculation results at different moments are output. It can be seen that the present invention uses the convection splitting algorithm and the hybrid format to numerically discretize the nonlinear convection terms, improves the accuracy and robustness, and greatly reduces the confusion error, so that the solver is very suitable for solving turbulent problems; on the other hand, there is no need to solve the pressure Poisson equation, there is no deadlock problem caused by the solution of the pressure equation, and there is no need to sacrifice accuracy in exchange for numerical robustness by adding non-physical artificial numerical viscosity. BRIEF DESCRIPTION OF THE DRAWINGS

[0042] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on the provided drawings without paying creative work.

[0043] Figure 1 A flow chart of a method for numerical solution of vortex dynamics fluid disclosed in this application;

[0044] Figure 2 This is a schematic diagram of the structure of a numerical solution device for vortex dynamics fluid disclosed in this application;

[0045] Figure 3 This is a structural diagram of an electronic device disclosed in this application. DETAILED DESCRIPTION

[0046] The following will be combined with the drawings in the embodiments of the present application to clearly and completely describe the technical solutions in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of the present invention.

[0047] The traditional incompressible fluid solver is based on the Navier-Stokes equations, which requires large-scale matrix iteration to solve the pressure Poisson equation, which brings risks to the overall numerical robustness of the solver. Furthermore, for the complex geometric shapes encountered in the aviation industry, the complex vortex structure is the key to affecting its dynamic and thermodynamic processes. Using the vortex dynamics equation as the control equation to develop numerical software can more clearly characterize the key processes in industrial fluids and provide more powerful tools for analyzing the mechanism of complex physical processes and further optimizing designs.

[0048] The traditional incompressible fluid solver is based on the Navier-Stokes equations. It uses the step-by-step solution of the nonlinear convection term, linear viscous dissipation term, and pressure Poisson term of the NS equations, and uses the pressure-velocity correction method to obtain the velocity information of the flow field. This is the basic principle and process of the numerical simulation of the fluid solver. However, the pressure Poisson equation needs to be solved through large-scale iterative matrix calculation. When the matrix is ​​singular, the iterative solution often cannot obtain satisfactory results, resulting in the iterative matrix calculation process to cycle back and forth, causing the entire fluid solver to deadlock. This is a problem often encountered in industrial software. The existing solution is to sacrifice the accuracy of the flow field and add non-physical artificial numerical viscosity to the control equation to make the iterative matrix converge as soon as possible. However, this makes the numerical simulation results given by the solver deviate from the real physics to some extent, which is also a compromise caused by the bottleneck of traditional technology.

[0049] To this end, the present invention provides a numerical solution for vortex dynamics fluid, which can directly and accurately solve the vortex dynamics fluid numerical value and avoid adding pressure and velocity corrections in the solution process.

[0050] Reference Figure 1 As shown, an embodiment of the present invention discloses a method for numerically solving vortex dynamics fluid, comprising:

[0051] Step S11: constructing a corresponding initial vortex dynamics equation based on the fluid numerical solution requirements of an industrial fluid application scenario; wherein the industrial fluid application scenario is an industrial scenario related to a vortex system whose structure meets a preset complexity condition.

[0052] In this embodiment, for compressor application scenarios, the top gap leakage vortex has an important influence on the compressor performance. It is necessary to solve the vortex distribution to analyze the formation and development of the top gap leakage vortex and its influence on the internal flow field of the compressor. Therefore, the corresponding initial vortex dynamics equation is the compressor vortex dynamics equation; for hydraulic machinery application scenarios, the complex spatiotemporal evolution of the vortex system in the key flow field has an important influence on the hydraulic characteristics of the fluid machinery. In order to study the spatiotemporal evolution process of the vortex system, analyze the relationship between the vortex distribution and the hydraulic characteristics, and create the corresponding initial vortex dynamics equation as the hydraulic machinery vortex dynamics equation; for aircraft application scenarios, the large-scale separation flow and complex vortex flow that appear in the large angle of attack turbulence have an important influence on the performance of the aircraft. In order to analyze the vortex distribution in these special flow field structures and provide a basis for the performance research and optimization design of the aircraft, the corresponding initial vortex dynamics equation is created as the aircraft vortex dynamics equation.

[0053] In this embodiment, the initial velocity vector, initial vortex vector, equation boundary conditions, and fluid kinematic viscosity in the fluid numerical solution requirements of the industrial fluid application scenario are used to construct the corresponding initial vortex dynamics equation; the initial vortex dynamics equation is:

[0054] ;

[0055] in, represents the initial vortex vector, represents the initial velocity vector, represents the fluid kinematic viscosity, Indicates time, represents the displacement vector, represents the velocity vector constrained by the solid wall, and the initial vortex vector at the solid wall The curl of the initial velocity vector constraint, represents the curl operator, represents the directions of the three coordinate axes in the three-dimensional space represented by the Einstein tensor, and The values ​​are 1, 2, and 3 respectively.

[0056] Initial vortex vector and the initial velocity vector pass Determined. It is understandable that all vectors in the initial vortex dynamics equation are expressed in the form of Einstein indicators. Since vorticity is a physical quantity that describes the rotation characteristics of the fluid, it reflects the degree of rotation of the fluid micro-group. When studying fluid flow, vorticity is a key physical quantity because it is closely related to the dynamic characteristics of the fluid. Through the study of vorticity, we can gain a deeper understanding of the processes of fluid rotation, mixing, and energy transfer.

[0057] In traditional technology, the solution of the NS equations for incompressible fluid flow is usually based on the original variables , it is necessary to solve the pressure Poisson equation through large-scale matrix iteration. When the matrix is ​​singular, the iterative solution often cannot obtain satisfactory results, which causes the iterative matrix calculation process to cycle back and forth, thereby causing the entire fluid solver to deadlock. This is a problem often encountered in industrial software. The existing solution is to sacrifice the accuracy of the flow field and add non-physical artificial numerical viscosity to the control equation to make the iterative matrix converge as soon as possible. However, this makes the numerical simulation results given by the solver deviate from the real physics to some extent. Therefore, the present invention further designs and creates an initial vortex dynamics equation based on the NS equation, which converts the original NS equation into a vortex equation. The equation form with φ as the core not only does not need to solve the pressure Poisson equation and abandons the pressure-velocity correction process, but is also more conducive to studying the vortex characteristics of the fluid.

[0058] It should be noted that the constraints at the solid wall are considered in the initial vortex dynamics equation. The velocity vector represents the solid wall constraint, and the vortex vector at the solid wall is constrained by the curl of the velocity vector, which is the curl operator. This consideration of the solid wall constraint is based on the actual physical situation, because in the actual fluid flow, the fluid close to the wall is restricted by the wall, and the distribution of its velocity and vorticity will be affected. By considering the solid wall constraint, the equation can more accurately describe the actual fluid flow.

[0059] Step S12: numerically discretize each component in the initial vortex dynamics equation according to a preset numerical discretization format and a convection splitting algorithm and using an initial velocity vector, an initial vortex vector and equation boundary conditions to obtain a target vortex dynamics equation for updating the vortex vector at the next moment.

[0060] In this embodiment, the first numerical discretization format of the Runge-Kutta third-order format coupled with the Crank-Nicolson second-order format is used to calculate the unsteady time term of the initial vortex dynamics equation based on the initial velocity vector, the initial vortex vector and the boundary conditions of the equation. Numerical discretization is performed to obtain a first discrete representation; a second numerical discretization format of WENO5 coupled with a sixth-order compact difference and a convection splitting algorithm are used to solve the nonlinear convection term of the initial vortex dynamics equation based on the initial velocity vector, the initial vortex vector and the boundary conditions of the equation. numerically discretize to obtain a second discrete representation; numerically discretize the dissipation term of the initial vortex dynamics equation based on the initial velocity vector, the initial vortex vector and the boundary conditions of the equation using a third numerical discretization format of the second-order central difference format; and vortex stretching term numerically discretize to obtain a third discrete representation; and use a fourth numerical discretization format of a second-order central difference format and based on an initial velocity vector, an initial vortex vector, and boundary conditions of the equation to calculate the curl term of the initial vortex dynamics equation. Numerical discretization is performed to obtain a fourth discrete representation; based on the first discrete representation, the second discrete representation, the third discrete representation, and the fourth discrete representation, a target vortex dynamics equation for updating the vortex vector at the next moment is constructed. It can be understood that in order to improve the numerical robustness and accuracy of the solver, a series of specific numerical discretization methods are used for the initial vortex dynamics equation. The selection of these discretization methods is based on the understanding of the physical properties of different terms and the expectation of the numerical solution effect. Unsteady time term: reflects the characteristics of the fluid changing with time, and is an important part of describing the time dimension in the fluid dynamics process. Nonlinear convection term: reflects the nonlinear characteristics of fluid convection motion, and has an important influence on the mixing, transmission and other processes of the fluid. Dissipation term: describes the process of fluid energy dissipation, which is an important manifestation of energy conservation in fluid dynamics. Vortex stretching term: related to the stretching change of vorticity, plays an important role in the development and evolution of vortex. Curl term: directly related to vorticity, is a key part of calculating vortex curl, and is crucial for accurately describing the vortex characteristics of the fluid. These terms together constitute the initial vortex dynamics equation, each describing the physical characteristics and dynamics of the fluid from different angles.

[0061] For the unsteady time terms, the Runge-Kutta 3rd order format coupled with the Crank-Nicolson 2nd order format is used for numerical discretization. This combination can better handle time-related changes and accurately capture the state of the fluid at different times. For the nonlinear convection terms, the hybrid format of the 6th order compact difference coupled WENO5 and the convection splitting technology are used for numerical discretization. This helps to improve the calculation accuracy of the convection terms, reduce aliasing errors, and make the solver more suitable for solving multi-scale turbulence problems. For the dissipation terms and vortex stretching terms, the second order central difference format is used for numerical discretization. This format can reasonably describe the dissipation and vortex stretching characteristics of the fluid. For the curl terms, the second order central difference discretization is used. This helps to accurately calculate the curl of the vorticity and maintain the physical accuracy of the equation.

[0062] Step S13: performing update calculation on the target vortex dynamics equation to obtain an updated vortex vector.

[0063] In this embodiment, after discretizing each term of the initial vortex dynamics equation using the above-mentioned discrete method, the target vortex dynamics equation can be obtained, and then the velocity vector update calculation is performed on the target vortex dynamics equation to realize the update of the vortex vector based on the update calculation of the velocity vector to obtain the updated vortex vector.

[0064] Step S14: using the updated vortex vector to update the stream function control equation associated with the target vortex dynamics equation, so as to obtain a target velocity vector determined based on the updated stream function vector.

[0065] In this embodiment, before the stream function control equation associated with the target vortex dynamics equation is updated using the updated vortex vector to obtain the target velocity vector determined based on the updated stream function vector, the method further includes: constructing an associated stream function control equation based on the target vortex dynamics equation, wherein the stream function control equation is:

[0066] ;

[0067] in, represents the stream function vector, and n represents the normal vector of the computational domain boundary.

[0068] Then, the stream function control equation is solved based on the updated vortex vector using the fifth numerical discretization format of the second-order central difference format to obtain the updated stream function vector, so as to calculate the target velocity vector based on the updated stream function vector. It can be understood that the right-hand source term The updated vortex vector is provided by updating the target vortex dynamics equation. The second-order derivative and the first-order derivative of the stream function control equation are both discretized using the second-order central difference. The stream function vector can be updated by numerically solving the stream function control equation.

[0069] In this embodiment, the updated stream function vector is substituted into , and the target velocity vector is obtained. It can be understood that the above formula gives the expression for updating the velocity vector. By calculating the curl of the updated stream function vector, the updated velocity vector, that is, the target velocity vector, is obtained. Specifically, the second-order central difference discretization curl term is used , in order to find the curl.

[0070] In this way, through the above three steps, namely, processing based on the initial vortex dynamics equation, solving the stream function control equation and applying the expression for updating the velocity vector, the updating of the vortex vector and the velocity vector is achieved.

[0071] Step S15: Use the target velocity vector and the updated vortex vector as the initial velocity vector and the initial vortex vector respectively, and jump to execute the steps of numerically discretizing the components in the initial vortex dynamics equation according to the preset numerical discretization format and convection splitting algorithm and using the initial velocity vector, the initial vortex vector and the boundary conditions of the equation to obtain the target vortex dynamics equation for updating the vortex vector at the next moment, until the preset iteration condition is met, stop the iteration and output the corresponding fluid numerical calculation results at different moments.

[0072] In this embodiment, since the fluid numerical solution requirement involves the flow of the fluid, and the flow of the fluid is a dynamic process, its vortex vector and velocity vector will change continuously over time. Therefore, it is necessary to gradually solve this information at each moment (iterative update) to accurately capture the flow state of the fluid at different time points and reflect the unsteady characteristics of the fluid. For example, in the aircraft application scenario, the flow field around the aircraft will undergo complex changes with the change of flight time, including the formation, development and dissipation of the vortex. Only by gradually solving the vortex vector and velocity vector information at each moment can we fully understand these dynamic changes. Therefore, after meeting the iterative update stop condition, the vortex vector and velocity vector at each moment calculated under the input conditions of the initial velocity vector and the initial vortex vector at the current moment are output to obtain the flow field vortex vector information and flow field velocity vector information at each moment.

[0073] In this way, the NS equations based on vortex dynamics do not need to solve the Poisson pressure problem, thus eliminating the numerical stiffness problem caused by the difficulty in determining the pressure boundary conditions and the pressure-velocity correction; in addition, the study found that in the study of boundary layer flow, since the vorticity can be ignored in the area far away from the wall, the vortex transport equation only needs to be solved in vorticity-related areas such as the boundary layer and wake in the numerical simulation.

[0074] It can be seen that the present application discloses a method for numerical solution of vortex dynamics fluid, including: constructing a corresponding initial vortex dynamics equation based on the fluid numerical solution requirements of an industrial fluid application scenario; wherein the industrial fluid application scenario is an industrial scenario related to a vortex system whose structure satisfies a preset complexity condition; numerically discretizing each component in the initial vortex dynamics equation according to a preset numerical discretization format and a convection splitting algorithm and using an initial velocity vector, an initial vortex vector and equation boundary conditions to obtain a target vortex dynamics equation for updating the vortex vector at the next moment; updating the target vortex dynamics equation to obtain an updated vortex vector; and using the updated vortex vector The stream function control equation associated with the target vortex dynamics equation is updated to obtain a target velocity vector determined based on the updated stream function vector; the target velocity vector and the updated vortex vector are used as the initial velocity vector and the initial vortex vector respectively, and the steps of numerically discretizing the components in the initial vortex dynamics equation according to the preset numerical discretization format and the convection splitting algorithm and using the initial velocity vector, the initial vortex vector and the boundary conditions of the equation are jumped to obtain the step of updating the target vortex dynamics equation for the vortex vector at the next moment, until the preset iteration condition is met, the iteration is stopped and the corresponding fluid numerical calculation results at different moments are output. It can be seen that the present invention uses the convection splitting algorithm and the hybrid format to numerically discretize the nonlinear convection terms, improves the accuracy and robustness, and greatly reduces the confusion error, so that the solver is very suitable for solving turbulent problems; on the other hand, there is no need to solve the pressure Poisson equation, there is no deadlock problem caused by the solution of the pressure equation, and there is no need to sacrifice accuracy in exchange for numerical robustness by adding non-physical artificial numerical viscosity.

[0075] Reference Figure 2 As shown, the present invention also discloses a numerical solution device for vortex dynamics fluid, including:

[0076] An equation construction module 11 is used to construct a corresponding initial vortex dynamics equation based on the fluid numerical solution requirements of an industrial fluid application scenario; wherein the industrial fluid application scenario is an industrial scenario related to a vortex system whose structure meets a preset complexity condition;

[0077] A numerical discretization module 12 is used to perform numerical discretization on each component in the initial vortex dynamics equation according to a preset numerical discretization format and a convection splitting algorithm and using an initial velocity vector, an initial vortex vector and equation boundary conditions to obtain a target vortex dynamics equation for updating the vortex vector at the next moment;

[0078] An updating module 13 is used to perform an update calculation on the target vortex dynamics equation to obtain an updated vortex vector;

[0079] A secondary updating module 14 is used to update the stream function control equation associated with the target vortex dynamics equation using the updated vortex vector to obtain a target velocity vector determined based on the updated stream function vector;

[0080] The result output module 15 is used to use the target velocity vector and the updated vortex vector as the initial velocity vector and the initial vortex vector respectively, and jump to execute the steps of numerically discretizing each component in the initial vortex dynamics equation according to the preset numerical discretization format and convection splitting algorithm and using the initial velocity vector, the initial vortex vector and the boundary conditions of the equation to obtain the target vortex dynamics equation for updating the vortex vector at the next moment, until the preset iteration condition is met, stop the iteration and output the corresponding fluid numerical calculation results at different moments.

[0081] It can be seen that the present application discloses constructing a corresponding initial vortex dynamics equation based on the fluid numerical solution requirements of an industrial fluid application scenario; wherein the industrial fluid application scenario is an industrial scenario related to a vortex system whose structure meets a preset complexity condition; according to a preset numerical discretization format and a convection splitting algorithm, each component in the initial vortex dynamics equation is numerically discretized using an initial velocity vector, an initial vortex vector, and an equation boundary condition to obtain a target vortex dynamics equation for updating the vortex vector at the next moment; the target vortex dynamics equation is updated and calculated to obtain an updated vortex vector; the stream function control equation associated with the target vortex dynamics equation is updated using the updated vortex vector to obtain a target velocity vector determined based on the updated stream function vector; the target velocity vector and the updated vortex vector are used as the initial velocity vector and the initial vortex vector, respectively, and the step of numerically discretizing each component in the initial vortex dynamics equation according to a preset numerical discretization format and a convection splitting algorithm using an initial velocity vector, an initial vortex vector, and an equation boundary condition to obtain a target vortex dynamics equation for updating the vortex vector at the next moment is jumped to execute until a preset iteration condition is met, the iteration is stopped, and the corresponding fluid numerical calculation results at different moments are output. It can be seen that the present invention adopts the convection splitting algorithm and the hybrid format to numerically discretize the nonlinear convection terms, which improves the accuracy and robustness while greatly reducing the confusion error, making the solver very suitable for solving turbulence problems; on the other hand, there is no need to solve the pressure Poisson equation, and there is no deadlock problem caused by solving the pressure equation, and there is no need to sacrifice accuracy in exchange for numerical robustness by adding non-physical artificial numerical viscosity.

[0082] Furthermore, the present application also discloses an electronic device. Figure 3 This is a structural diagram of an electronic device 20 according to an exemplary embodiment. The content in the diagram cannot be regarded as any limitation on the scope of use of the present application.

[0083] Figure 3 A schematic diagram of the structure of an electronic device 20 provided in an embodiment of the present application. The electronic device 20 may specifically include: at least one processor 21, at least one memory 22, a power supply 23, a communication interface 24, an input / output interface 25, and a communication bus 26. The memory 22 is used to store a computer program, which is loaded and executed by the processor 21 to implement the relevant steps in the numerical solution method for vortex dynamics fluid disclosed in any of the aforementioned embodiments. In addition, the electronic device 20 in this embodiment may specifically be an electronic computer.

[0084] In this embodiment, the power supply 23 is used to provide working voltage for each hardware device on the electronic device 20; the communication interface 24 can create a data transmission channel between the electronic device 20 and the external device, and the communication protocol it follows is any communication protocol that can be applied to the technical solution of the present application, and is not specifically limited here; the input and output interface 25 is used to obtain external input data or output data to the outside world, and its specific interface type can be selected according to specific application needs and is not specifically limited here.

[0085] Among them, the processor 21 may include one or more processing cores, such as a 4-core processor, an 8-core processor, etc. The processor 21 can be implemented in at least one hardware form of DSP (Digital Signal Processing), FPGA (Field-Programmable Gate Array), and PLA (Programmable Logic Array). The processor 21 may also include a main processor and a coprocessor. The main processor is a processor for processing data in the awake state, also known as a CPU (Central Processing Unit); the coprocessor is a low-power processor for processing data in the standby state. In some embodiments, the processor 21 may be integrated with a GPU (Graphics Processing Unit), which is responsible for rendering and drawing the content to be displayed on the display screen. In some embodiments, the processor 21 may also include an AI (Artificial Intelligence) processor, which is used to process computing operations related to machine learning.

[0086] In addition, the memory 22, as a carrier for storing resources, can be a read-only memory, a random access memory, a disk or an optical disk, etc. The resources stored thereon can include an operating system 221, a computer program 222, etc., and the storage method can be temporary storage or permanent storage.

[0087] Among them, the operating system 221 is used to manage and control the hardware devices and computer programs 222 on the electronic device 20, so as to realize the operation and processing of the processor 21 on the massive data 223 in the memory 22, which can be Windows Server, Netware, Unix, Linux, etc. In addition to including a computer program that can be used to complete the vortex dynamics fluid numerical solution method performed by the electronic device 20 disclosed in any of the aforementioned embodiments, the computer program 222 can further include a computer program that can be used to complete other specific tasks. In addition to data received by the electronic device and transmitted from an external device, the data 223 can also include data collected by its own input and output interface 25, etc.

[0088] Furthermore, the present application also discloses a computer-readable storage medium for storing a computer program; wherein, when the computer program is executed by a processor, the aforementioned disclosed method for numerically solving vortex dynamics fluid is implemented. The specific steps of the method can be referred to the corresponding contents disclosed in the aforementioned embodiments, and will not be repeated here.

[0089] In this specification, each embodiment is described in a progressive manner, and each embodiment focuses on the differences from other embodiments. The same or similar parts between the embodiments can be referred to each other. For the device disclosed in the embodiment, since it corresponds to the method disclosed in the embodiment, the description is relatively simple, and the relevant parts can be referred to the method part.

[0090] Professionals may further appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented with electronic hardware, computer software, or a combination of the two. In order to clearly illustrate the interchangeability of hardware and software, the composition and steps of each example have been generally described in terms of function in the above description. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel may use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of this application. The steps of the method or algorithm described in conjunction with the embodiments disclosed herein can be implemented directly with hardware, a software module executed by a processor, or a combination of the two. The software module can be placed in a random access memory RAM (Random Access Memory), memory, read-only memory ROM (Read Only Memory), electrically programmable EPROM (Electrically Programmable Read Only Memory), electrically erasable programmable EEPROM (ElectricErasable Programmable Read Only Memory), register, hard disk, removable disk, CD-ROM (CompactDisc-Read Only Memory), or any other form of storage medium known in the technical field.

[0091] Finally, it should be noted that, in this article, relational terms such as first and second, etc. are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, the elements defined by the sentence "comprise a ..." do not exclude the presence of other identical elements in the process, method, article or device including the elements.

[0092] The scheme provided by the present invention is introduced in detail above. Specific examples are used in this article to illustrate the principle and implementation mode of the present invention. The description of the above embodiments is only used to help understand the method of the present invention and its core idea. At the same time, for those skilled in the art, according to the idea of ​​the present invention, there will be changes in the specific implementation mode and application scope. In summary, the content of this specification should not be understood as limiting the present invention.

Claims

1. A method for numerical solution of vortex dynamics fluid, characterized in that: include: Constructing the corresponding initial vortex dynamics equation based on the fluid numerical solution requirements of the industrial fluid application scenario; wherein the industrial fluid application scenario is an industrial scenario related to a vortex system whose structure meets the preset complexity conditions; According to a preset numerical discretization format and a convection splitting algorithm and using an initial velocity vector, an initial vortex vector and equation boundary conditions, each component in the initial vortex dynamics equation is numerically discretized to obtain a target vortex dynamics equation for updating the vortex vector at a next moment; Performing an update calculation on the target vortex dynamics equation to obtain an updated vortex vector; Using the updated vortex vector, the stream function control equation associated with the target vortex dynamics equation is updated to obtain a target velocity vector determined based on the updated stream function vector; The target velocity vector and the updated vortex vector are used as the initial velocity vector and the initial vortex vector respectively, and the steps of numerically discretizing the components in the initial vortex dynamics equation according to the preset numerical discretization format and the convection splitting algorithm and using the initial velocity vector, the initial vortex vector and the boundary conditions of the equation are jumped to obtain the target vortex dynamics equation for updating the vortex vector at the next moment, until the preset iteration condition is met, the iteration is stopped and the corresponding fluid numerical calculation results at different moments are output.

2. The method for numerical solution of vortex dynamics fluid according to claim 1, characterized in that: The fluid numerical solution requirement based on the industrial fluid application scenario is used to construct the corresponding initial vortex dynamics equation, including: Based on the initial velocity vector, initial vortex vector, equation boundary conditions, and fluid kinematic viscosity in the fluid numerical solution requirements of industrial fluid application scenarios, the corresponding initial vortex dynamics equation is constructed; the initial vortex dynamics equation is: ； in, represents the initial vortex vector, represents the initial velocity vector, represents the fluid kinematic viscosity, Indicates time, represents the displacement vector, represents the velocity vector constrained by the solid wall, and the initial vortex vector at the solid wall The curl of the initial velocity vector constraint, represents the curl operator, represents the directions of the three coordinate axes in the three-dimensional space represented by the Einstein tensor, and The values ​​are 1, 2, and 3 respectively.

3. The method for numerical solution of vortex dynamics fluid according to claim 2, characterized in that: The method of numerically discretizing each component in the initial vortex dynamics equation according to a preset numerical discretization format and a convection splitting algorithm and using an initial velocity vector, an initial vortex vector and equation boundary conditions to obtain a target vortex dynamics equation for updating the vortex vector at the next moment includes: The first numerical discretization format of the Runge-Kutta third-order format coupled with the Crank-Nicolson second-order format is used to calculate the unsteady time term of the initial vortex dynamics equation based on the initial velocity vector, the initial vortex vector and the boundary conditions of the equation. Performing numerical discretization to obtain a first discrete representation; The nonlinear convection term of the initial vortex dynamics equation is solved by using the second numerical discretization format of WENO5 coupled with the sixth-order compact difference and the convection splitting algorithm based on the initial velocity vector, the initial vortex vector and the boundary conditions of the equation. Performing numerical discretization to obtain a second discrete representation; The dissipation term of the initial vortex dynamics equation is numerically discretized using the third numerical discretization format of the second-order central difference format and based on the initial velocity vector, the initial vortex vector and the boundary conditions of the equation. and vortex stretching term Performing numerical discretization to obtain a third discrete representation; The curl term of the initial vortex dynamics equation is calculated based on the initial velocity vector, the initial vortex vector and the boundary conditions of the equation using the fourth numerical discretization format of the second-order central difference format. Performing numerical discretization to obtain a fourth discrete representation; A target vortex dynamics equation for updating the vortex vector at the next moment is constructed based on the first discrete representation, the second discrete representation, the third discrete representation, and the fourth discrete representation.

4. The method for numerically solving vortex dynamics fluid according to claim 3, characterized in that: Before the stream function control equation associated with the target vortex dynamics equation is updated using the updated vortex vector to obtain the target velocity vector determined based on the updated stream function vector, the method further includes: The associated stream function control equation is constructed based on the target vortex dynamics equation, and the stream function control equation is: ； in, represents the stream function vector, and n represents the normal vector of the computational domain boundary.

5. The method for numerical solution of vortex dynamics fluid according to claim 4, characterized in that: The updating of the stream function control equation associated with the target vortex dynamics equation using the updated vortex vector to obtain a target velocity vector determined based on the updated stream function vector includes: The stream function control equation is solved based on the updated vortex vector using the fifth numerical discretization format of the second-order central difference format to obtain an updated stream function vector, so as to calculate the target velocity vector based on the updated stream function vector.

6. The method for numerical solution of vortex dynamics fluid according to claim 5, characterized in that: The calculating of the target velocity vector based on the updated stream function vector comprises: Substitute the updated stream function vector into , and obtain the target velocity vector.

7. The method for numerically solving vortex dynamics fluid according to any one of claims 1 to 6, characterized in that: The output corresponding to the fluid numerical calculation results at different times includes: The vortex vector and velocity vector at each moment calculated based on the input conditions of the initial velocity vector and the initial vortex vector at the current moment are output to obtain the flow field vortex vector information and flow field velocity vector information at each moment.

8. A numerical solution device for vortex dynamics fluid, characterized in that: include: An equation building module, used to build a corresponding initial vortex dynamics equation based on the fluid numerical solution requirements of an industrial fluid application scenario; wherein the industrial fluid application scenario is an industrial scenario related to a vortex system whose structure meets a preset complexity condition; A numerical discretization module is used to perform numerical discretization on each component in the initial vortex dynamics equation according to a preset numerical discretization format and a convection splitting algorithm and using an initial velocity vector, an initial vortex vector and equation boundary conditions to obtain a target vortex dynamics equation for updating the vortex vector at the next moment; An updating module, used for performing an updating calculation on the target vortex dynamics equation to obtain an updated vortex vector; A secondary updating module, used for updating the stream function control equation associated with the target vortex dynamics equation using the updated vortex vector, so as to obtain a target velocity vector determined based on the updated stream function vector; The result output module is used to use the target velocity vector and the updated vortex vector as the initial velocity vector and the initial vortex vector respectively, and jump to execute the steps of numerically discretizing the components in the initial vortex dynamics equation according to the preset numerical discretization format and the convection splitting algorithm and using the initial velocity vector, the initial vortex vector and the boundary conditions of the equation to obtain the target vortex dynamics equation for updating the vortex vector at the next moment, until the preset iteration condition is met, stop the iteration and output the corresponding fluid numerical calculation results at different moments.

9. An electronic device, characterized in that: include: Memory, used to store computer programs; A processor is used to execute the computer program to implement the steps of the method for numerical solution of vortex dynamics fluid according to any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that: Used to store computer programs; wherein, when the computer program is executed by a processor, the steps of the method for numerical solution of vortex dynamics fluid as described in any one of claims 1 to 7 are implemented.

Images