A method, device, equipment and medium for numerical solution of eddy dynamics fluid

Abstract

The present application discloses a method, apparatus, equipment and medium for numerically solving vortex dynamics fluid, which relates to the technical field of fluid mechanics, and includes: constructing a corresponding initial vortex dynamics equation; numerically discretizing 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; updating the target vortex dynamics equation and obtaining an updated vortex vector; using the updated vortex vector to update a flow function control equation associated with the target vortex dynamics equation and obtaining a target velocity vector determined based on the updated flow function vector; using the target velocity vector and the updated vortex vector as the initial velocity vector and the initial vortex vector, and jumping to execute the steps of obtaining the target vortex dynamics equation for updating the vortex vector at the next moment until a preset iteration condition is met, stopping the iteration and outputting the corresponding fluid numerical calculation results at different moments.

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] Traditional incompressible fluid solvers, based on the Navier-Stokes equations, require large-scale matrix iterations to solve the Poisson pressure equation, which poses risks to the solver's overall numerical robustness. Furthermore, for the complex geometries encountered in the aviation industry, the complex vortex structures are key to influencing their dynamic and thermodynamic processes. Developing numerical software using vortex dynamics equations as the governing equations can more clearly characterize key processes in industrial fluids, providing a more powerful tool for analyzing the mechanisms of complex physical processes and further optimizing designs.

[0003] Traditional incompressible fluid solvers are based on the Navier-Stokes equations. They employ a step-by-step approach to solve the nonlinear convection term, linear viscous dissipation term, and pressure Poisson term of the Navier-Stokes equations, and employ a pressure-velocity correction method to obtain velocity information for the flow field. This is the fundamental principle and process of numerical simulation using fluid solvers. However, the pressure Poisson equation requires large-scale iterative matrix calculations to solve. When the matrix is ​​singular, the iterative solution often fails to produce satisfactory results, causing the iterative matrix calculation process to loop repeatedly and deadlock the entire fluid solver. This is a common problem encountered in industrial software. Existing solutions sacrifice flow field accuracy by introducing non-physical, artificial numerical viscosity into the governing equations to accelerate convergence of the iterative matrix. However, this approach results in numerical simulation results that deviate somewhat from real physics, a compromise driven by bottlenecks in traditional technologies. Summary of the Invention

[0004] In view of this, the present invention aims to provide a method, apparatus, device, and medium for numerically solving vortex dynamics fluids, which can directly and accurately solve the numerical value of vortex dynamics fluids, avoiding the addition of pressure and velocity corrections during 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 corresponding initial vortex dynamics equations based on the fluid numerical solution requirements of industrial fluid application scenarios; wherein the industrial fluid application scenarios are industrial scenarios related to vortex systems whose structures meet preset complexity conditions;

[0007] numerically discretizing each component of the initial vortex dynamics equation according to a preset numerical discretization format and a convection splitting algorithm and utilizing an initial velocity vector, an initial vortex vector, and boundary conditions of the equation 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] updating 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;

[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 step of performing numerical discretization on 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 is 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 construction of the corresponding initial vortex dynamics equation based on the fluid numerical solution requirements of the industrial fluid application scenario includes:

[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, The velocity vector representing the solid wall constraint, 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 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 based on the initial velocity vector, initial vortex vector and equation boundary conditions by using the second numerical discretization format of WENO5 coupled with the 6th order compact difference and the convection splitting algorithm. 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 Perform numerical discretization to obtain a third discrete representation;

[0019] The curl term of the initial vortex dynamics equation is obtained by using the fourth 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. 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 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 a fifth numerical discretization format of a second-order central difference format to obtain an updated stream function vector, so as to calculate a target velocity vector based on the updated stream function vector.

[0027] Optionally, 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 outputting of corresponding fluid numerical calculation results at different moments 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 construction module for constructing corresponding initial vortex dynamics equations 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 for numerically discretizing each component of the initial vortex dynamics equation according to a preset numerical discretization format and a convection splitting algorithm and using the initial velocity vector, the initial vortex vector, and the boundary conditions of the equation to obtain a target vortex dynamics equation for updating the vortex vector at the next moment;

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

[0035] a secondary updating module, configured to update a 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 convection splitting algorithm and using the initial velocity vector, the initial vortex vector and the equation boundary conditions 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 meets a preset complexity condition; according to a preset numerical discrete format and a convection splitting algorithm and using an initial velocity vector, an initial vortex vector and equation boundary conditions, the components in the initial vortex dynamics equation are numerically discretized 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; 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 step of performing numerical discretization on the components of 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 is 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 numerical calculation results of the fluid at the corresponding 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 aliasing error, making the solver very suitable for solving turbulent 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. 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 following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are merely embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on the provided drawings without paying any creative work.

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

[0044] Figure 2 This is a schematic structural diagram 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 embodiments described 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 efforts are within the scope of protection of the present invention.

[0047] Traditional incompressible fluid solvers, based on the Navier-Stokes equations, require large-scale matrix iterations to solve the Poisson pressure equation, which poses risks to the solver's overall numerical robustness. Furthermore, for the complex geometries encountered in the aviation industry, the complex vortex structures are key to influencing their dynamic and thermodynamic processes. Developing numerical software using vortex dynamics equations as the governing equations can more clearly characterize key processes in industrial fluids, providing a more powerful tool for analyzing the mechanisms of complex physical processes and further optimizing designs.

[0048] Traditional incompressible fluid solvers are based on the Navier-Stokes equations. They employ a step-by-step approach to solve the nonlinear convection term, linear viscous dissipation term, and pressure Poisson term of the Navier-Stokes equations, and employ a pressure-velocity correction method to obtain velocity information for the flow field. This is the fundamental principle and process of numerical simulation using fluid solvers. However, the pressure Poisson equation requires large-scale iterative matrix calculations to solve. When the matrix is ​​singular, the iterative solution often fails to produce satisfactory results, causing the iterative matrix calculation process to loop repeatedly and deadlock the entire fluid solver. This is a common problem encountered in industrial software. Existing solutions sacrifice flow field accuracy by introducing non-physical, artificial numerical viscosity into the governing equations to accelerate convergence of the iterative matrix. However, this approach results in numerical simulation results that deviate somewhat from real physics, a compromise driven by bottlenecks in traditional technologies.

[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 performance of the compressor. 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 domain has an important influence on the hydraulic characteristics of the fluid machinery. In order to study the spatiotemporal evolution process of the vortex system and analyze the relationship between the vortex distribution and the hydraulic characteristics, the corresponding initial vortex dynamics equation is created 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 corresponding initial vortex dynamics equation is constructed based on 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; 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, The velocity vector representing the solid wall constraint, 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 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's understandable that all vectors in the initial vortex dynamics equations are expressed using Einstein indices. Vorticity is a physical quantity that describes the rotational properties of a fluid and reflects the degree of rotation of fluid clusters. Vorticity is a key physical quantity in the study of fluid flow because it is closely related to the fluid's dynamic properties. Studying vorticity provides insights into processes such as 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 iterative solution. 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's important to note that the initial vortex dynamics equations take into account the constraints at the solid wall. The velocity vector, denoted by the solid wall, and the vortex vector at the solid wall, constrained by the curl of the velocity vector, are the curl operators. This consideration of the solid wall constraint is based on practical physics: in real fluid flow, fluid near a wall is confined by the wall, affecting its velocity and vorticity distribution. By accounting for the solid wall constraint, the equations more accurately describe real-world 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 the initial velocity vector, the initial vortex vector and the boundary conditions of the equation 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 6th-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. Numerical discretization is performed to obtain a second discrete representation; a third numerical discretization format of the second-order central difference format is used to 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. and vortex stretching term Numerical discretization is performed to obtain a third discrete representation; a curl term of the initial vortex dynamics equation is obtained by using a fourth numerical discretization format of a second-order central difference format and based on the initial velocity vector, the initial vortex vector and the boundary conditions of the equation. Numerical discretization is performed to obtain a fourth discrete representation; and based on the first, second, third, and fourth discrete representations, the target vortex dynamics equation for updating the vortex vector at the next moment is constructed. It will be appreciated that to improve the numerical robustness and accuracy of the solver, a series of specific numerical discretization methods are employed for the initial vortex dynamics equation. These discretization methods are selected based on an understanding of the physical properties of different terms and the desired numerical solution performance. The unsteady time term reflects the time-varying characteristics of the fluid and is an important component in describing the time dimension of fluid dynamics. The nonlinear convection term reflects the nonlinear characteristics of fluid convection and has a significant impact on processes such as mixing and transport. The dissipation term describes the process of fluid energy dissipation and is a key manifestation of energy conservation in fluid dynamics. The vortex stretch term is related to the stretching change of vorticity and plays an important role in the development and evolution of vortices. The curl term is directly related to vorticity and is a key component in calculating vortex curl, which is crucial for accurately describing the vortex characteristics of the fluid. These terms together constitute the initial vortex dynamics equation, each describing the physical properties and dynamics of the fluid from different perspectives.

[0061] For the unsteady time terms, a Runge-Kutta third-order scheme coupled with a Crank-Nicolson second-order scheme is used for numerical discretization. This combination can better handle time-dependent changes and accurately capture the state of the fluid at different times. For the nonlinear convection terms, a hybrid scheme of sixth-order compact differences coupled with WENO5 and convection splitting technology is used for numerical discretization. This helps improve the computational accuracy of the convection terms, reduces aliasing errors, and makes the solver more suitable for solving multiscale turbulence problems. For the dissipation terms and vortex stretching terms, a second-order central difference scheme is used for numerical discretization. This scheme can reasonably describe the dissipation and vortex stretching characteristics of the fluid. For the curl terms, a second-order central difference discretization is used. This helps accurately calculate the curl of the vorticity and maintain the physical accuracy of the equations.

[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 of the target vortex dynamics equation is performed to realize the update of the vortex vector based on the updated 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 to obtain a target velocity vector determined based on the updated stream function vector.

[0065] In this embodiment, before using the updated vortex vector to update the stream function control equation associated with the target vortex dynamics equation 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, and 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 source term on the right is The updated vortex vector is provided by updating the target vortex dynamics equation. The second-order and first-order derivatives of the stream function governing equation are both discretized using second-order central differences. The stream function vector can be updated by numerically solving the stream function governing equation.

[0069] In this embodiment, the updated stream function vector is substituted into , get the target velocity vector. It can be understood that the above formula gives the expression for velocity vector update. 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. , to find the rotation.

[0070] In this way, the update of the vortex vector and the velocity vector is achieved through the above three steps, namely, the processing of the initial vortex dynamics equation, the solution of the stream function control equation and the application of the expression for the velocity vector update.

[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 equation boundary conditions 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 requires 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 an aircraft application scenario, the flow field around the aircraft will undergo complex changes as the flight time changes, including processes such as the formation, development, and dissipation of vortices. Only by gradually solving the vortex vector and velocity vector information at each moment can we fully understand these dynamic changes. Therefore, after the iterative update stop condition is met, the vortex vector and velocity vector at each moment calculated based on the input conditions of the initial velocity vector and 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 the vortex-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 meets a preset complexity condition; according to a preset numerical discrete format and a convection splitting algorithm and using an initial velocity vector, an initial vortex vector and equation boundary conditions, the components in the initial vortex dynamics equation are numerically discretized 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; 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 step of performing numerical discretization on the components of 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 is 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 numerical calculation results of the fluid at the corresponding 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 aliasing error, making the solver very suitable for solving turbulent 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. 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, comprising:

[0076] An equation construction module 11 is used to construct corresponding initial vortex dynamics equations 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 for numerically discretizing each component of the initial vortex dynamics equation according to a preset numerical discretization format and a convection splitting algorithm and using the initial velocity vector, the initial vortex vector, and the boundary conditions of the equation 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 configured to update 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;

[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 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.

[0081] It can be seen that the present application discloses a method for 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, the initial velocity vector, the initial vortex vector and the boundary conditions of the equation are used to numerically discretize the components in the initial vortex dynamics equation 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 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 a target vortex dynamics equation for updating the vortex vector at the next moment is jumped to execute 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 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 aliasing 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. There is no need to sacrifice accuracy in exchange for numerical robustness by adding non-physical artificial numerical viscosity.

[0082] Furthermore, the embodiment of 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 should not be considered as any limitation to the scope of application of the present application.

[0083] Figure 3 This is 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 operating 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. The communication protocol it follows is any communication protocol that can be applied to the technical solution of this 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. 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 resource storage, can be a read-only memory, random access memory, disk or CD, 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 various hardware devices and computer programs 222 on the electronic device 20 to enable the processor 21 to calculate and process the massive data 223 in the memory 22. It 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 computer programs that can be used to complete other specific tasks. In addition to including 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.

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

[0089] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on its differences from the other embodiments. Reference can be made to the descriptions of the identical or similar parts between the various embodiments. For the devices disclosed in the embodiments, since they correspond to the methods disclosed in the embodiments, the descriptions are relatively simple, and the relevant parts can be referred to the descriptions of the methods.

[0090] Professionals may further appreciate that the units and algorithmic steps of each example described in conjunction with the embodiments disclosed herein can be implemented in 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. Professionals and technicians can 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 using 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 (Electric Erasable Programmable Read Only Memory), registers, hard disk, removable disk, CD-ROM (Compact Disc-Read Only Memory), or any other form of storage medium known in the technical field.

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

[0092] The above is a detailed introduction to the solution provided by the present invention. Specific examples are used herein to illustrate the principles and implementation methods 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 ideas. At the same time, for those skilled in the art, according to the ideas of the present invention, there may be changes in the specific implementation methods and application scopes. In summary, the contents of this specification should not be understood as limiting the present invention.

Claims

1. A numerical solution method for eddy dynamics fluid, characterized in that: include: Constructing corresponding initial vortex dynamics equations based on the fluid numerical solution requirements of industrial fluid application scenarios; wherein the industrial fluid application scenarios are industrial scenarios related to vortex systems whose structures meet preset complexity conditions; numerically discretizing each component of the initial vortex dynamics equation according to a preset numerical discretization format and a convection splitting algorithm and utilizing an initial velocity vector, an initial vortex vector, and boundary conditions of the equation 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; updating 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; 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 performing numerical discretization on 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 is 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 numerically solving eddy-dynamic fluid according to claim 1, characterized in that: The fluid numerical solution requirements based on industrial fluid application scenarios are used to construct the corresponding initial vortex dynamics equations, 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, The velocity vector representing the solid wall constraint, 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 three-dimensional space represented by the Einstein tensor, and The values ​​are 1, 2, and 3 respectively.

3. The method for numerically solving eddy-dynamic 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 utilizing an initial velocity vector, an initial vortex vector, and boundary conditions of the equation 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 based on the initial velocity vector, initial vortex vector and equation boundary conditions by using the second numerical discretization format of WENO5 coupled with the 6th order compact difference and the convection splitting algorithm. 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 Perform numerical discretization to obtain a third discrete representation; The curl term of the initial vortex dynamics equation is obtained by using the fourth 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. 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 eddy-dynamic fluid according to claim 3, characterized in that: Before 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: 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 numerically solving eddy-dynamic 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 a fifth numerical discretization format of a second-order central difference format to obtain an updated stream function vector, so as to calculate a target velocity vector based on the updated stream function vector.

6. The method for numerical solution of eddy-dynamic fluid according to claim 5, characterized in that: The calculating the target velocity vector based on the updated stream function vector includes: 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 of 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 eddy dynamics fluid, characterized in that: include: An equation construction module for constructing corresponding initial vortex dynamics equations 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 for numerically discretizing each component of the initial vortex dynamics equation according to a preset numerical discretization format and a convection splitting algorithm and using the initial velocity vector, the initial vortex vector, and the boundary conditions of the equation to obtain a target vortex dynamics equation for updating the vortex vector at the next moment; An updating module, configured to perform an updating calculation on the target vortex dynamics equation to obtain an updated vortex vector; a secondary updating module, configured to update a 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 convection splitting algorithm and using the initial velocity vector, the initial vortex vector and the equation boundary conditions 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 a computer program; wherein, when the computer program is executed by a processor, the steps of the numerical solution method for vortex dynamics fluid are implemented as described in any one of claims 1 to 7.

Images