Method, device and equipment for solving isotropic uniform turbulence of vortex dynamics and medium
By using preset vortex dynamics equations and linear turbulence forced force models in isotropic uniform turbulence solution, the numerical rigidity problem caused by the Gibbs phenomenon in traditional methods is solved, and the robustness of the solver is improved and computational efficiency.
Patent Information
- Application Number
- CN202411768621.4
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2024-12-04
- Publication Date
- 2024-12-31
- Estimated Expiration
- Not applicable · inactive patent
AI Technical Summary
The traditional isotropic uniform turbulence solver will cause the Gibbs phenomenon when the Navier-Stokes equation is projected into the spectral space, resulting in the problem of weak numerical rigidity and robustness.
Using preset vortex dynamics equations based on the initial excitation coefficient, combined with the preset linear turbulence forced force model, the process is processed through numerical discrete format and convective splitting algorithm, avoiding projection into the spectral space, and directly updating the speed and vortex information until Stop iteration operation.
The Gibbs phenomenon is effectively avoided, the calculation efficiency and the robustness of the solver are improved, and the stability of numerical solutions is ensured.
Smart Images
Figure CN119227596A_ABST
Abstract
Description
Technical Field
[0001] The present invention relates to the technical field of turbulence solution, and in particular to a method, device, equipment and medium for solving isotropic uniform turbulence in vortex dynamics. Background Art
[0002] Uniform isotropic turbulence is an ideal simplified form of turbulence. Most of the existing results of turbulence statistical theory are based on the assumption of uniform isotropy. Therefore, it is of great significance to construct a solution method for uniform isotropic turbulence. Faced with industrial-scale numerical simulation scenarios, traditional isotropic uniform turbulence solvers are based on the Navier-Stokes equations on the one hand and the spectral method on the other. The disadvantage of the traditional method is that after the NS equations are projected into the spectral space and then the forced turbulent force source term is applied, it will lead to the Gibbs phenomenon, resulting in numerical rigidity and weak numerical robustness. Summary of the invention
[0003] In view of this, the purpose of the present invention is to provide a method, device, equipment and medium for solving isotropic uniform turbulence in vortex dynamics, which can avoid the Gibbs phenomenon caused by projecting the equation into the spectral space and applying the forced turbulence force source term, and improve the robustness of the numerical solution. The specific scheme is as follows:
[0004] In a first aspect, the present application discloses a method for solving isotropic uniform turbulence in vortex dynamics, comprising:
[0005] An initial excitation coefficient is set based on the isotropic uniform turbulence solution requirement, so as to create a preset vortex dynamics equation corresponding to the isotropic uniform turbulence solution requirement based on the initial excitation coefficient;
[0006] According to a preset numerical discretization format and a convection splitting algorithm and using the initial velocity information and the initial vortex information in the isotropic uniform turbulence solution requirement, the components of the preset vortex dynamics equation are numerically discretized to obtain updated velocity information and updated vortex information based on the discretization results, and the updated velocity information and the updated vortex information are used as the initial velocity information and the initial vortex information of the next time step, respectively;
[0007] The updated velocity information and the preset linear turbulence forcing force model are used to calculate and update the initial excitation coefficient of the isotropic uniform turbulence solution requirement, and the step of creating a preset vortex dynamics equation corresponding to the isotropic uniform turbulence solution requirement based on the initial excitation coefficient is jumped to execute until the iterative operation is stopped to obtain the updated velocity information and updated vortex information at each time step corresponding to the isotropic uniform turbulence solution requirement.
[0008] Optionally, the creating a preset vortex dynamics equation corresponding to the isotropic uniform turbulence solution requirement based on the initial excitation coefficient includes:
[0009] Create a preset vortex dynamics equation corresponding to the isotropic uniform turbulence solution requirement; wherein the preset vortex dynamics equation is:
[0010] ; Where A represents the initial excitation coefficient, Expressed as the vortex vector, represents the velocity vector, represents the fluid kinematic viscosity, Indicates time, represents the displacement vector; represents the velocity vector constrained by the solid wall, represents the curl, which is used to constrain the vortex vector at the solid wall. 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.
[0011] Optionally, the numerically discretizing each component of the preset vortex dynamics equation according to the preset numerical discretization format and the convection splitting algorithm and using the initial velocity information and the initial vortex information in the isotropic uniform turbulence solution requirement to obtain the updated velocity information and the updated vortex information based on the discretization result includes:
[0012] The first numerical discretization format of the Runge-Kutta third-order format coupled with the Crank-Nicolson second-order format is used to solve the unsteady time term of the preset vortex dynamics equation based on the initial velocity information and initial vortex information in the isotropic uniform turbulence solution requirement. Performing numerical discretization processing to obtain a first discrete result;
[0013] The nonlinear convection term of the preset eddy dynamics equation is solved by using the hybrid format of the sixth-order compact difference coupling WENO5 and the second numerical discretization format of the convection splitting algorithm and the initial velocity information and initial vortex information in the isotropic uniform turbulence solution requirement. Performing numerical discretization processing to obtain a second discrete result;
[0014] The third numerical discretization format of the second-order central difference format is used to solve the dissipation term of the preset vortex dynamics equation based on the initial velocity information and initial vortex information in the isotropic uniform turbulence solution requirement. , vortex stretch term , curl term Performing numerical discretization processing to obtain a third discrete result;
[0015] Updated velocity information and updated vortex information are acquired based on the first discrete result, the second discrete result, and the third discrete result.
[0016] Optionally, the acquiring updated velocity information and updated vortex information based on the first discrete result, the second discrete result, and the third discrete result includes:
[0017] Calculate and obtain process vortex information output by the preset vortex dynamics equation according to the first discrete result, the second discrete result, and the third discrete result, and use the process vortex information as updated vortex information;
[0018] numerically solving a stream function vector control equation associated with the preset vortex dynamics equation using the process vortex information to obtain an updated stream function vector;
[0019] The curl of the updated stream function vector is calculated to obtain updated velocity information.
[0020] Optionally, the performing curl calculation on the updated stream function vector to obtain updated velocity information includes:
[0021] use The updated stream function vector is subjected to curl calculation to obtain updated velocity information, wherein: Represents the updated stream function vector.
[0022] Optionally, before numerically solving a stream function vector control equation associated with the preset vortex dynamics equation using the process vortex information to obtain an updated stream function vector, the method further includes:
[0023] Build As the stream function vector governing equation; where, represents the stream function vector, and n represents the normal vector of the computational domain boundary.
[0024] Optionally, before using the updated velocity information and a preset linear turbulence forcing force model to calculate and update the initial excitation coefficient required for solving the isotropic uniform turbulence, the method further includes:
[0025] Build As the preset linear turbulence forcing model; where < > represents the statistical averaging of the computational domain, represents the statistical average dissipation rate, represents the statistically average turbulent kinetic energy, represents the target turbulent kinetic energy.
[0026] In a second aspect, the present application discloses an isotropic uniform turbulence solving device for vortex dynamics, comprising:
[0027] An equation creation module, used for setting an initial excitation coefficient based on an isotropic uniform turbulence solution requirement, so as to create a preset vortex dynamics equation corresponding to the isotropic uniform turbulence solution requirement based on the initial excitation coefficient;
[0028] A numerical discretization module, used for performing numerical discretization processing on each component of the preset vortex dynamics equation according to a preset numerical discretization format and a convection splitting algorithm and using the initial velocity information and initial vortex information in the isotropic uniform turbulence solution requirement, so as to obtain updated velocity information and updated vortex information based on the discretization result, and using the updated velocity information and the updated vortex information as the initial velocity information and initial vortex information of the next time step respectively;
[0029] A solution module is used to calculate and update the initial excitation coefficient of the isotropic uniform turbulence solution requirement using the updated velocity information and a preset linear turbulence forcing force model, and jump to execute the step of creating a preset vortex dynamics equation corresponding to the isotropic uniform turbulence solution requirement based on the initial excitation coefficient until the iterative operation is stopped to obtain the updated velocity information and updated vortex information at each time step corresponding to the isotropic uniform turbulence solution requirement.
[0030] In a third aspect, the present application discloses an electronic device, comprising:
[0031] Memory, used to store computer programs;
[0032] A processor is used to execute the computer program to implement the steps of the isotropic uniform turbulence solution method for vortex dynamics disclosed above.
[0033] 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 method for solving the isotropic uniform turbulence of vortex dynamics disclosed above are implemented.
[0034] It can be seen that the present application discloses a method for solving isotropic uniform turbulence in vortex dynamics, including: setting an initial excitation coefficient based on the isotropic uniform turbulence solution requirement, so as to create a preset vortex dynamics equation corresponding to the isotropic uniform turbulence solution requirement based on the initial excitation coefficient; performing numerical discretization processing on each component of the preset vortex dynamics equation according to a preset numerical discretization format and a convection splitting algorithm and using the initial velocity information and initial vortex information in the isotropic uniform turbulence solution requirement, so as to obtain updated velocity information and updated vortex information based on the discretization results, and The updated velocity information and the updated vortex information are used as the initial velocity information and initial vortex information of the next time step respectively; the updated velocity information and the preset linear turbulence forcing force model are used to calculate and update the initial excitation coefficient of the isotropic uniform turbulence solution requirement, and the step of creating the preset vortex dynamics equation corresponding to the isotropic uniform turbulence solution requirement based on the initial excitation coefficient is jumped to execute until the iterative operation is stopped to obtain the updated velocity information and updated vortex information at each time step corresponding to the isotropic uniform turbulence solution requirement. It can be seen that by adopting the preset vortex dynamics equation combined with the preset linear turbulence forcing force model, the method of projecting to the spectral space is abandoned, and the complex Fourier transform step is avoided, thereby improving the calculation efficiency, and completely avoiding the occurrence of the Gibbs phenomenon, greatly improving the robustness of the solver. BRIEF DESCRIPTION OF THE DRAWINGS
[0035] 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.
[0036] Figure 1 A flow chart of a method for solving isotropic uniform turbulence of vortex dynamics disclosed in the present application;
[0037] Figure 2 A turbulence solution calculation process of coupling a linear force calculation model with a vortex dynamics solver disclosed in the present application to become an isotropic uniform turbulence solver;
[0038] Figure 3 This is a schematic structural diagram of an isotropic uniform turbulence solving device for vortex dynamics disclosed in the present application;
[0039] Figure 4 This is a structural diagram of an electronic device disclosed in this application. DETAILED DESCRIPTION
[0040] 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.
[0041] In the face of industrial-scale numerical simulation scenarios, the traditional isotropic uniform turbulence solver is based on the Navier-Stokes equations on the one hand and the spectral method on the other. The disadvantage of the traditional method is that the NS equations are projected into the spectral space and then the forced turbulence source term is applied, which will lead to the Gibbs phenomenon, produce numerical rigidity, and have weak numerical robustness.
[0042] To this end, the present invention provides an isotropic uniform turbulence solution for vortex dynamics, which can avoid the Gibbs phenomenon caused by projecting the equation into the spectral space and applying a forced turbulent force source term, thereby improving the robustness of the numerical solution.
[0043] Reference Figure 1 As shown, an embodiment of the present invention discloses a method for solving isotropic uniform turbulence of vortex dynamics, comprising:
[0044] Step S11: setting an initial excitation coefficient based on the isotropic uniform turbulence solution requirement, so as to create a preset vortex dynamics equation corresponding to the isotropic uniform turbulence solution requirement based on the initial excitation coefficient.
[0045] It is understandable that, referring to Figure 2 As shown, according to the isotropic uniform turbulence solution requirement, the initial excitation coefficient is set to 0, and then a preset vortex dynamics equation corresponding to the isotropic uniform turbulence solution requirement is created; wherein the preset vortex dynamics equation is:
[0046] ; Where A represents the initial excitation coefficient, Expressed as the vortex vector, represents the velocity vector, represents the fluid kinematic viscosity, Indicates time, represents the displacement vector; represents the velocity vector constrained by the solid wall, represents the curl, which is used to constrain the vortex vector at the solid wall. 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.
[0047] In this way, the preset vortex dynamics equation is constructed based on the basic principles of vortex dynamics. Vortex dynamics describes the movement of fluids by studying the generation, transmission and dissipation of vorticity (vortex vector). This equation can accurately describe the relationship between vorticity and physical quantities such as velocity vector, fluid kinematic viscosity, time and displacement vector, and provides a theoretical framework for studying the turbulent characteristics of fluids. Therefore, the above preset vortex dynamics equation is constructed to solve the isotropic uniform turbulence problem.
[0048] Step S12: According to the preset numerical discretization format and convection splitting algorithm and using the initial velocity information and initial vortex information in the isotropic uniform turbulence solution requirement, the components of the preset vortex dynamics equation are numerically discretized to obtain updated velocity information and updated vortex information based on the discretization results, and the updated velocity information and the updated vortex information are used as the initial velocity information and initial vortex information of the next time step, respectively.
[0049] 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 solve the unsteady time term of the preset vortex dynamics equation based on the initial velocity information and initial vortex information in the isotropic uniform turbulence solution requirement. Performing numerical discretization to obtain a first discrete result; using a hybrid format of the sixth-order compact difference coupling WENO5, a second numerical discretization format of the convection splitting algorithm and based on the initial velocity information and initial vortex information in the isotropic uniform turbulence solution requirement, the nonlinear convection term of the preset vortex dynamics equation is calculated. Performing numerical discretization processing to obtain a second discrete result; using a third numerical discretization format of a second-order central difference format to solve the dissipation term of the preset vortex dynamics equation based on the initial velocity information and initial vortex information in the isotropic uniform turbulence solution requirement; , vortex stretch term , curl term Perform numerical discretization to obtain a third discrete result; obtain updated velocity information and updated vortex information based on the first discrete result, the second discrete result, and the third discrete result. It can be understood that in order to numerically simulate and solve isotropic uniform turbulence, the continuous physical equation (preset vortex dynamics equation) needs to be discretized. The construction form of the preset vortex dynamics equation is convenient for adopting various numerical discretization methods, such as using the Runge-Kutta 3rd order format coupled with the Crank-Nicolson 2nd order format to numerically discretize the unsteady time term, using the hybrid format of the 6th order compact difference coupling WENO5 and the convective splitting technology to numerically discretize the nonlinear convection term of the equation, using the second order central difference format to numerically discretize the dissipative term and vortex stretching term, and using the second order central difference to discretize the curl term. These discrete methods can convert the equation into a form suitable for computer calculation, thereby realizing the numerical solution of turbulence problems. There is no need to project the NS equation into the spectral space and then apply the forced turbulence source term, that is, it will not cause the Gibbs phenomenon, numerical rigidity, and weak numerical robustness.
[0050] In this embodiment, the process vortex information output by the preset vortex dynamics equation is calculated and obtained according to the first discrete result, the second discrete result, and the third discrete result, and the process vortex information is used as the updated vortex information; the stream function vector control equation associated with the preset vortex dynamics equation is numerically solved using the process vortex information to obtain the updated stream function vector; the curl of the updated stream function vector is calculated to obtain the updated velocity information. It can be understood that the process vortex information output by the discrete result calculation equation, that is, the updated vortex information , where in the process of obtaining the updated vortex information, the process velocity information of the equation output is often calculated at the same time. The vortex information is defined by the curl of the velocity information (the vortex vector is constrained by the curl of the velocity vector at the solid wall), and they jointly describe the vortex motion state of the fluid. From a physical essence, the generation, development and evolution of the vortex are closely related to the change of the velocity field. Therefore, when solving the preset vortex dynamics equation, it is necessary to update the vortex vector and the velocity vector at the same time to accurately reflect the real physical process of the vortex motion. Then, according to the updated vortex information, the associated stream function vector control equation (Poission control equation) is further updated. It should be noted that the updated vortex information obtained after the preset vortex dynamics equation is updated is further used as the right source term of the stream function vector control equation. Specifically, construct As the stream function vector governing equation; where, represents the stream function vector, and n represents the normal vector of the computational domain boundary. Therefore, the stream function vector in the above stream function vector control equation can be solved Get the updated stream function vector . Further, when obtaining the utilization The updated stream function vector is subjected to curl calculation to obtain updated velocity information, wherein: Represents the updated stream function vector. It should be noted that the preset vortex dynamics equation and the stream function vector control equation are both used to describe the flow state of the fluid. The preset vortex dynamics equation describes the vortex characteristics and flow of the fluid through the vortex vector and the velocity vector, while the stream function vector control equation describes the flow of the fluid through the stream function vector. The two describe the fluid flow from different angles and complement each other.
[0051] Step S13: Utilize the updated velocity information and the preset linear turbulence forcing model to calculate and update the initial excitation coefficient of the isotropic uniform turbulence solution requirement, and jump to execute the step of creating the preset vortex dynamics equation corresponding to the isotropic uniform turbulence solution requirement based on the initial excitation coefficient until the iterative operation is stopped to obtain the updated velocity information and updated vortex information at each time step corresponding to the isotropic uniform turbulence solution requirement.
[0052] In this embodiment, in order to solve the isotropic turbulence problem, in addition to setting the initial excitation coefficient A to 0 in the first calculation, the initial excitation coefficient at different time steps needs to be continuously updated with the linear turbulence forcing model and the updated velocity information. For details, refer to Figure 2 As shown in the figure, the expression of the adaptive excitation coefficient A of isotropic uniform turbulence calculated by the linear force model is shown. This coefficient A is very important for the simulation of isotropic uniform forced turbulence. In isotropic uniform forced turbulence, an external forcing force is required to maintain the turbulent state, and the coefficient A is related to the effect of this forcing force. It is determined based on the dissipation rate, turbulent kinetic energy, and the target turbulent kinetic energy given by the user, which are statistically averaged over the entire computational domain.
[0053] In this embodiment, construct As the preset linear turbulence forcing model; where < > represents the statistical averaging of the computational domain, represents the statistical average dissipation rate, represents the statistically average turbulent kinetic energy, Represents the target turbulent kinetic energy. It is understandable that after obtaining the updated velocity vector, only part of the description of the turbulent flow state is completed. Isotropic uniform forced turbulence is a special turbulent state, which needs to consider the influence of external forcing force. By presetting the linear turbulence forcing force model to determine the adaptive excitation coefficient A, the effect of the forcing force on turbulence is further considered, thereby more perfectly simulating the characteristics of isotropic uniform forced turbulence. For example, different values of the excitation coefficient A may lead to different turbulence intensities and structures. By reasonably determining the value of A, the simulation results can be made more consistent with the actual physical situation. Moreover, the determined adaptive excitation coefficient A is a key parameter in the entire solution process. It, together with the physical quantities in the previously preset vortex dynamics equation and the stream function vector control equation, constitutes a complete description of isotropic uniform forced turbulence. For example, in the preset vortex dynamics equation, the coefficient A will affect the update process of the vortex vector and the velocity vector (although it is not explicitly reflected, it is related in the actual physical and numerical relationship). Through this synergistic effect, the various equations cooperate with each other to gradually solve various physical quantities of isotropic uniform forced turbulence, including the turbulence excitation coefficient A, vortex vector, velocity vector, etc., to achieve accurate simulation of this special turbulent state.
[0054] It can be seen that the present application discloses a method for solving isotropic uniform turbulence in vortex dynamics, including: setting an initial excitation coefficient based on the isotropic uniform turbulence solution requirement, so as to create a preset vortex dynamics equation corresponding to the isotropic uniform turbulence solution requirement based on the initial excitation coefficient; performing numerical discretization processing on each component of the preset vortex dynamics equation according to a preset numerical discretization format and a convection splitting algorithm and using the initial velocity information and initial vortex information in the isotropic uniform turbulence solution requirement, so as to obtain updated velocity information and updated vortex information based on the discretization results, and The updated velocity information and the updated vortex information are used as the initial velocity information and initial vortex information of the next time step respectively; the updated velocity information and the preset linear turbulence forcing force model are used to calculate and update the initial excitation coefficient of the isotropic uniform turbulence solution requirement, and the step of creating the preset vortex dynamics equation corresponding to the isotropic uniform turbulence solution requirement based on the initial excitation coefficient is jumped to execute until the iterative operation is stopped to obtain the updated velocity information and the updated vortex information at each time step corresponding to the isotropic uniform turbulence solution requirement. It can be seen that by adopting the preset vortex dynamics equation combined with the preset linear turbulence forcing force model, the method of projecting to the spectral space is abandoned, and the complex Fourier transformation step is avoided, thereby improving the calculation efficiency, and completely avoiding the occurrence of the Gibbs phenomenon, greatly improving the robustness of the solver.
[0055] Reference Figure 3As shown, the present invention also discloses an isotropic uniform turbulence solving device for vortex dynamics, comprising:
[0056] An equation creation module 11 is used to set an initial excitation coefficient based on the isotropic uniform turbulence solution requirement, so as to create a preset vortex dynamics equation corresponding to the isotropic uniform turbulence solution requirement based on the initial excitation coefficient;
[0057] A numerical discretization module 12 is used to perform numerical discretization processing on each component of the preset vortex dynamics equation according to a preset numerical discretization format and a convection splitting algorithm and using the initial velocity information and initial vortex information in the isotropic uniform turbulence solution requirement, so as to obtain updated velocity information and updated vortex information based on the discretization result, and use the updated velocity information and the updated vortex information as the initial velocity information and initial vortex information of the next time step respectively;
[0058] The solution module 13 is used to calculate and update the initial excitation coefficient of the isotropic uniform turbulence solution requirement using the updated velocity information and the preset linear turbulence forcing force model, and jump to execute the step of creating a preset vortex dynamics equation corresponding to the isotropic uniform turbulence solution requirement based on the initial excitation coefficient until the iterative operation is stopped to obtain the updated velocity information and updated vortex information at each time step corresponding to the isotropic uniform turbulence solution requirement.
[0059] It can be seen that the present application discloses setting an initial excitation coefficient based on the isotropic uniform turbulence solution requirement, so as to create a preset vortex dynamics equation corresponding to the isotropic uniform turbulence solution requirement based on the initial excitation coefficient; numerically discretizing each component of the preset vortex dynamics equation according to a preset numerical discretization format and a convection splitting algorithm and using the initial velocity information and initial vortex information in the isotropic uniform turbulence solution requirement, so as to obtain updated velocity information and updated vortex information based on the discretization result, and use the updated velocity information and the updated vortex information as the initial velocity information and initial vortex information of the next time step respectively; using the updated velocity information and the preset linear turbulence forcing force model to calculate and update the initial excitation coefficient of the isotropic uniform turbulence solution requirement, and jumping to execute the step of creating a preset vortex dynamics equation corresponding to the isotropic uniform turbulence solution requirement based on the initial excitation coefficient until the iterative operation is stopped, so as to obtain the updated velocity information and updated vortex information at each time step corresponding to the isotropic uniform turbulence solution requirement. It can be seen that by adopting the preset eddy dynamics equation combined with the preset linear turbulence forcing model, the method of projection to the spectral space is abandoned, the complex Fourier transform steps are avoided, thereby improving the calculation efficiency, completely avoiding the occurrence of the Gibbs phenomenon, and greatly improving the robustness of the solver.
[0060] Furthermore, the present application also discloses an electronic device. Figure 4 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.
[0061] Figure 4 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 isotropic uniform turbulence solution method for vortex dynamics disclosed in any of the aforementioned embodiments. In addition, the electronic device 20 in this embodiment may specifically be an electronic computer.
[0062] 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.
[0063] 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.
[0064] 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.
[0065] 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 massive data 223 in the memory 22 by the processor 21, which can be Windows Server, Netware, Unix, Linux, etc. In addition to including a computer program that can be used to complete the isotropic uniform turbulence solution method of eddy dynamics 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 transmitted from an external device received by the electronic device, the data 223 can also include data collected by its own input and output interface 25, etc.
[0066] 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 isotropic uniform turbulence solution method for vortex dynamics disclosed above is implemented. For the specific steps of the method, reference may be made to the corresponding contents disclosed in the aforementioned embodiments, and no further description will be given here.
[0067] 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.
[0068] 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.
[0069] 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.
[0070] 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 solving isotropic uniform turbulence in vortex dynamics, characterized in that: include: An initial excitation coefficient is set based on the isotropic uniform turbulence solution requirement, so as to create a preset vortex dynamics equation corresponding to the isotropic uniform turbulence solution requirement based on the initial excitation coefficient; According to a preset numerical discretization format and a convection splitting algorithm and using the initial velocity information and the initial vortex information in the isotropic uniform turbulence solution requirement, the components of the preset vortex dynamics equation are numerically discretized to obtain updated velocity information and updated vortex information based on the discretization results, and the updated velocity information and the updated vortex information are used as the initial velocity information and the initial vortex information of the next time step, respectively; The updated velocity information and the preset linear turbulence forcing force model are used to calculate and update the initial excitation coefficient of the isotropic uniform turbulence solution requirement, and the step of creating a preset vortex dynamics equation corresponding to the isotropic uniform turbulence solution requirement based on the initial excitation coefficient is jumped to execute until the iterative operation is stopped to obtain the updated velocity information and updated vortex information at each time step corresponding to the isotropic uniform turbulence solution requirement.
2. The isotropic uniform turbulence solution method of vortex dynamics according to claim 1, characterized in that: The creating of a preset vortex dynamics equation corresponding to the isotropic uniform turbulence solution requirement based on the initial excitation coefficient comprises: Create a preset vortex dynamics equation corresponding to the isotropic uniform turbulence solution requirement; wherein the preset vortex dynamics equation is: ； Where A represents the initial excitation coefficient, Expressed as the vortex vector, represents the velocity vector, represents the fluid kinematic viscosity, Indicates time, represents the displacement vector; represents the velocity vector constrained by the solid wall, represents the curl, which is used to constrain the vortex vector at the solid wall. 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 isotropic uniform turbulence solution method of vortex dynamics according to claim 2 is characterized in that: The method of performing numerical discretization processing on each component of the preset vortex dynamics equation according to the preset numerical discretization format and the convection splitting algorithm and utilizing the initial velocity information and the initial vortex information in the isotropic uniform turbulence solution requirement to obtain updated velocity information and updated vortex information based on the discretization results includes: The first numerical discretization format of the Runge-Kutta third-order format coupled with the Crank-Nicolson second-order format is used to solve the unsteady time term of the preset vortex dynamics equation based on the initial velocity information and initial vortex information in the isotropic uniform turbulence solution requirement. Performing numerical discretization processing to obtain a first discrete result; The nonlinear convection term of the preset eddy dynamics equation is solved by using the hybrid format of the sixth-order compact difference coupling WENO5 and the second numerical discretization format of the convection splitting algorithm and the initial velocity information and initial vortex information in the isotropic uniform turbulence solution requirement. Performing numerical discretization processing to obtain a second discrete result; The third numerical discretization format of the second-order central difference format is used to solve the dissipation term of the preset eddy dynamics equation based on the initial velocity information and initial vortex information in the isotropic uniform turbulence solution requirement. , vortex stretch term , curl term Performing numerical discretization processing to obtain a third discrete result; Updated velocity information and updated vortex information are acquired based on the first discrete result, the second discrete result, and the third discrete result.
4. The isotropic uniform turbulence solution method of vortex dynamics according to claim 3 is characterized in that: The obtaining updated velocity information and updated vortex information based on the first discrete result, the second discrete result, and the third discrete result includes: Calculate and obtain process vortex information output by the preset vortex dynamics equation according to the first discrete result, the second discrete result, and the third discrete result, and use the process vortex information as updated vortex information; numerically solving a stream function vector control equation associated with the preset vortex dynamics equation using the process vortex information to obtain an updated stream function vector; The curl of the updated stream function vector is calculated to obtain updated velocity information.
5. The isotropic uniform turbulence solution method of vortex dynamics according to claim 4 is characterized in that: The performing curl calculation on the updated stream function vector to obtain updated velocity information includes: use The updated stream function vector is subjected to curl calculation to obtain updated velocity information, wherein: Represents the updated stream function vector.
6. The isotropic uniform turbulence solution method of vortex dynamics according to claim 4, characterized in that: Before numerically solving the stream function vector control equation associated with the preset vortex dynamics equation using the process vortex information to obtain the updated stream function vector, the method further includes: Build As the stream function vector governing equation; where, represents the stream function vector, and n represents the normal vector of the computational domain boundary.
7. The isotropic uniform turbulence solution method of vortex dynamics according to claim 2, characterized in that: Before calculating and updating the initial excitation coefficient of the isotropic uniform turbulence solution requirement using the updated velocity information and the preset linear turbulence forcing force model, the method further includes: Build As the preset linear turbulence forcing model; where < > represents the statistical averaging of the computational domain, represents the statistical average dissipation rate, represents the statistically average turbulent kinetic energy, represents the target turbulent kinetic energy.
8. An isotropic uniform turbulence solver for vortex dynamics, characterized in that: include: An equation creation module, used for setting an initial excitation coefficient based on an isotropic uniform turbulence solution requirement, so as to create a preset vortex dynamics equation corresponding to the isotropic uniform turbulence solution requirement based on the initial excitation coefficient; A numerical discretization module, used for performing numerical discretization processing on each component of the preset vortex dynamics equation according to a preset numerical discretization format and a convection splitting algorithm and using the initial velocity information and initial vortex information in the isotropic uniform turbulence solution requirement, so as to obtain updated velocity information and updated vortex information based on the discretization result, and using the updated velocity information and the updated vortex information as the initial velocity information and initial vortex information of the next time step respectively; A solution module is used to calculate and update the initial excitation coefficient of the isotropic uniform turbulence solution requirement using the updated velocity information and a preset linear turbulence forcing force model, and jump to execute the step of creating a preset vortex dynamics equation corresponding to the isotropic uniform turbulence solution requirement based on the initial excitation coefficient until the iterative operation is stopped to obtain the updated velocity information and updated vortex information at each time step corresponding to the isotropic uniform turbulence solution requirement.
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 isotropic uniform turbulence solution method for vortex dynamics as described in 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 isotropic uniform turbulence solution method for vortex dynamics as described in any one of claims 1 to 7 are implemented.