Anisotropic turbulence model modeling method based on turbulence pulsation speed and dissipation rate

Abstract

The invention discloses an anisotropic turbulence model modeling method based on turbulence pulsation speed and dissipation rate. The method comprises the following steps: establishing a pulsation speed transport equation; the turbulent flow kinetic energy equation and the turbulent flow dissipation rate equation are solved, and the numerical value of the turbulent flow dissipation rate in the flow field is obtained; solving the transportation equation of the pulsation speed to obtain the turbulence pulsation speed; constructing a viscosity coefficient of the turbulent flow according to the obtained turbulent flow dissipation rate and the turbulent flow pulsation speed; calculating each stress item of the turbulent flow by using the constructed turbulent flow viscosity coefficient; and substituting the obtained turbulence stress item into an average flow equation and a turbulence pulsation speed equation of turbulence to obtain the anisotropic turbulence model. According to the anisotropic turbulence model equation established by the invention, the numerical calculation precision of the turbulence boundary layer is improved, and the established anisotropic turbulence model equation only contains a unique artificial parameter, so that the experimental error caused by human factors can be greatly reduced, the result is more accurate, and the data is more persuasive.

Description

Technical Field

[0001] The present invention relates to the field of numerical model construction, and in particular to a method for modeling an anisotropic turbulence model based on turbulent pulsation velocity and dissipation rate. Background Art

[0002] The turbulent boundary layer contains a large number of vortex structures with random properties. The irregular movement of these vortex structures in the flow field induces the pulsation of various basic physical quantities in the flow field. Its complex multi-scale characteristics make it very difficult to predict turbulence. The turbulent boundary layer has a significant impact on the aerodynamic shape design, structural optimization, flow and combustion organization inside the engine, etc. In the development of high-speed aircraft and advanced engines, the shock wave / turbulent boundary layer interaction (SWTBLI) problem has become one of the difficulties that need to be solved urgently.

[0003] Considering the timeliness and economy of engineering calculations, the turbulence model method based on the Reynolds-averaged Navier-Stokes (NS) equations still plays an important role. The turbulence model is one of the key factors affecting the accurate calculation of shock wave / turbulent boundary layer interference. Therefore, it is of great significance to develop effective turbulence models in a targeted manner.

[0004] At present, the turbulence model developed under the traditional turbulence theory and mathematical averaging method is widely used in fluid numerical calculation. Its development is mainly based on the consideration of simple turbulence problems. Therefore, there are obvious difficulties in simulating more complex turbulence such as separated flow, shock wave, and complex vortex system. The main reason for the lack of simulation ability of these complex flows is that the turbulence model is not accurate enough. At present, the derivation of Reynolds stress turbulence model is relatively rigorous. In theory, it can be applied to anisotropic turbulence calculation. However, due to its shortcomings such as excessive complexity, numerous empirical coefficients, high computational cost and poor robustness, and the calculation accuracy of various types of flows is not completely superior to the eddy viscosity model, the linear eddy viscosity model is still widely used in the actual numerical calculation of aviation compressors. Since the eddy viscosity turbulence model is based on the eddy viscosity assumption, this type of model has many shortcomings such as "isotropy" and "instantaneous equilibrium", which is powerless to calculate anisotropic flows caused by corner secondary flow, end wall leakage flow, rotation effect, shock wave boundary layer interference, etc. in the compressor. In addition, the eddy viscosity turbulence model that is currently widely used also has the disadvantages of having many artificial empirical coefficients and poor versatility. Summary of the invention

[0005] The present invention aims to solve one of the technical problems existing in the related art at least to a certain extent.

[0006] To this end, the present invention provides a method for constructing a turbulent viscosity coefficient using turbulent pulsation velocity and turbulent dissipation rate parameters and establishing an anisotropic turbulence model equation, thereby further improving the numerical calculation accuracy of the turbulent boundary layer.

[0007] In order to achieve the above-mentioned purpose, the present invention provides, on the one hand, a method for modeling an anisotropic turbulence model based on turbulent pulsation velocity and dissipation rate, comprising:

[0008] Establishing the pulsation speed The transport equation of

[0009] Solve the turbulent kinetic energy K equation and the turbulent dissipation rate ω equation to obtain the value of the turbulent dissipation rate ω in the flow field;

[0010] Solve for the pulsation speed The transport equation of turbulent flow is used to obtain the turbulent pulsation velocity

[0011] To obtain the turbulent dissipation rate ω and turbulent pulsation velocity Constructing the viscosity coefficient μ of turbulence T ;

[0012] Using the constructed turbulent viscosity coefficient μ T Calculate various turbulent stress terms;

[0013] The obtained turbulent stress term is substituted into the mean flow equation of turbulence and the turbulent fluctuation velocity equation to obtain the anisotropic turbulence model.

[0014] A further preferred technical solution of the present invention is that the pulsating speed is established The transport equation of , the specific method is:

[0015] The turbulent pulsating velocity in the flow field is divided into two groups, the average value of each group is not zero, and the pulsating velocity is established. The transport equation is expressed as:

[0016]

[0017] Where ρ represents the fluid density, u i 、u j represents the average flow velocity, represents the drift flow velocity, τ ij represents the mean flow turbulent viscous stress, represents the turbulent viscous stress of the drift flow, t is time, x i is the space coordinate.

[0018] Preferably, the turbulent kinetic energy K equation and the turbulent dissipation rate ω equation are solved to obtain the value of the turbulent dissipation rate ω in the flow field; the specific method is:

[0019] The turbulent kinetic energy K equation is established and the expression is:

[0020]

[0021] The turbulent dissipation rate ω equation is established, and the expression is:

[0022]

[0023] Among them, v t represents the eddy kinematic viscosity coefficient, Ω is the vorticity, a 1 =0.31; ρ represents fluid density, K represents turbulent kinetic energy, ω represents turbulent dissipation rate, β, σ ω , σ ω2 , β * , γ represents the model constant, ψ is the specific heat ratio, ψ = 0.41, σ ω1 =0.5,β=0.075,β * =0.09,σ ω2 =0.856; μ, μ i Indicates laminar viscosity coefficient and turbulent viscosity coefficient, F 1 and F 2 is a mixed function, which is related to the distance from the point to the wall; D is the total differential symbol, t is time, x j is the space coordinate.

[0024] Preferably, the mixing function F 1 The expression is:

[0025]

[0026] Mixing function F 2 The expression is:

[0027]

[0028] Among them, CD kω It represents the positive direction of the orthogonal diffusion term added to the standard k-ω turbulence model and the standard k-ε turbulence model. y represents the minimum distance from a point in the flow field to the surface of an object, and ε represents the turbulent energy dissipation rate.

[0029] As a preference, use a turbulent pulsating velocity and turbulence dissipation rate ω, the viscosity coefficient of constructed turbulence μ T , the expression is:

[0030]

[0031] Where ρ represents the density of the selected fluid, and the coefficient c is an artificial parameter, which is initially set to c=1.0.

[0032] As a preference, the constructed turbulent viscosity coefficient μ T The calculated turbulent stress terms include:

[0033]

[0034] Among them, μ T represents the viscous turbulence coefficient; μ L represents the viscous laminar coefficient; S represents the deformation rate, δ ij stands for Kronecker operator; represents the turbulent viscous stress of the drift flow, represents the turbulent strain rate of the drift flow, is the average pulsating velocity of the lateral deviation, Represents the rate of change of drift flow velocity.

[0035] Another aspect of the present invention provides a non-transitory computer-readable storage medium having computer instructions stored thereon, which enable a computer to execute the above-mentioned anisotropic turbulence modeling method based on turbulent pulsation velocity and dissipation rate.

[0036] Yet another aspect of the present invention provides an electronic device, comprising: a processor, a communication interface, a memory and a communication bus, wherein the processor, the communication interface and the memory communicate with each other via the communication bus, and the processor calls the logic instructions in the memory to execute the above-mentioned anisotropic turbulence modeling method based on turbulent pulsation velocity and dissipation rate.

[0037] On the other hand, the present invention provides a computer program product, which includes a computer program stored on a non-transitory computer-readable storage medium. When the computer program is executed by a processor, the computer executes the above-mentioned anisotropic turbulence model modeling method based on turbulent pulsation velocity and dissipation rate.

[0038] Beneficial effects: The present invention adopts the turbulent viscosity coefficient constructed by using the turbulent pulsation velocity and turbulent dissipation rate parameters, and establishes the anisotropic turbulence model equation, so as to improve the numerical calculation accuracy of the turbulent boundary layer. The constructed anisotropic turbulence model equation only contains a unique artificial parameter, which can greatly reduce the experimental error caused by human factors, making the results more accurate and the data more convincing. BRIEF DESCRIPTION OF THE DRAWINGS

[0039] Figure 1 A flowchart of the construction and application of an anisotropic turbulence model in the embodiment;

[0040] Figure 2A three-dimensional diagram of a zero pressure gradient plate, a geometric model established in the embodiment;

[0041] Figure 3 A mesh division diagram of a geometric model - a zero pressure gradient plate in an embodiment;

[0042] Figure 4 A code form diagram of the pulsating velocity equation in the new turbulence model established in the embodiment;

[0043] Figure 5 is a code form diagram of the turbulent viscosity coefficient in the new turbulence model established in the embodiment;

[0044] Figure 6 It is the residual convergence diagram of the velocity, pulsating velocity, etc. of the new turbulence model established in the embodiment;

[0045] Figure 7 is a comparison diagram of the dimensionless average velocity profile of the new turbulence model established in the embodiment;

[0046] Figure 8 It is a comparison diagram of the dimensionless average velocity profile corrected by the new turbulence model coefficient C established in the embodiment. DETAILED DESCRIPTION

[0047] In order to make the purpose, technical solution and advantages of the present invention clearer, the technical solution of the present invention will be clearly and completely described below in conjunction with the drawings in the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments, and they should not be understood as limitations on the present invention. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention. In the description of the present invention, it should be understood that the terms used are only for descriptive purposes and cannot be understood as indicating or implying relative importance.

[0048] Combine the following Figure 1-Figure 8 The present invention describes the anisotropic turbulence modeling method based on turbulent pulsation velocity and dissipation rate.

[0049] Embodiment 1: This embodiment first provides a method for modeling an anisotropic turbulence model based on turbulent pulsation velocity and dissipation rate to establish a new turbulence model. Then, the established new turbulence model is applied in practice and numerical simulation calculations are performed.

[0050] Anisotropic turbulence modeling methods based on turbulent pulsation velocity and dissipation rate, such as Figure 1 As shown, the steps include:

[0051] S1. Establish pulsation speed The transport equation of ; the specific method is:

[0052] The turbulent pulsating velocity in the flow field is divided into two groups, the average value of each group is not zero, and the pulsating velocity is established. The transport equation is expressed as:

[0053]

[0054] Where ρ represents the fluid density, u i 、u j represents the average flow velocity, represents the drift flow velocity, τ ij represents the mean flow turbulent viscous stress, represents the turbulent viscous stress of the drift flow, t is time, x i is the space coordinate.

[0055] S2. Solve the turbulent kinetic energy K equation and the turbulent dissipation rate ω equation to obtain the value of the turbulent dissipation rate ω in the flow field; the specific method is:

[0056] The turbulent kinetic energy K equation is established and the expression is:

[0057]

[0058] The turbulent dissipation rate ω equation is established, and the expression is:

[0059]

[0060] Among them, v t represents the eddy kinematic viscosity coefficient, Ω is the vorticity, a 1 =0.31; ρ represents fluid density, K represents turbulent kinetic energy, ω represents turbulent dissipation rate, β, σ ω , σ ω2 , β * , γ represents the model constant, ψ is the specific heat ratio, ψ = 0.41, σ ω1 =0.5,β=0.075,β * =0.09,σ ω2 =0.856; μ, μ i Indicates laminar viscosity coefficient and turbulent viscosity coefficient, F 1 and F 2 is a mixed function, which is related to the distance from the point to the wall; D is the total differential symbol, t is time, x j is the space coordinate.

[0061] Mixing function F 1 The expression is:

[0062]

[0063] Mixing function F 2 The expression is:

[0064]

[0065] Among them, CD kω It represents the positive direction of the orthogonal diffusion term added to the standard k-ω turbulence model and the standard k-ε turbulence model. y represents the minimum distance from a point in the flow field to the surface of an object, and ε represents the turbulent energy dissipation rate.

[0066] S3. Solving the pulsation speed The transport equation of turbulent flow is used to obtain the turbulent pulsation velocity

[0067] S4, to obtain the turbulent dissipation rate ω and turbulent pulsation velocity Constructing the viscosity coefficient μ of turbulence T ; The viscosity coefficient of the constructed turbulence μ T , the expression is:

[0068]

[0069] Where ρ represents the density of the selected fluid, and the coefficient c is an artificial parameter, which is initially set to c=1.0.

[0070] S5. Using the constructed turbulent viscosity coefficient μ T Calculate various turbulent stress terms, including:

[0071]

[0072]

[0073] Among them, μ T represents the viscous turbulence coefficient; μ L represents the viscous laminar coefficient; S represents the deformation rate, δ ij stands for Kronecker operator; represents the turbulent viscous stress of the drift flow, represents the turbulent strain rate of the drift flow, is the average pulsating velocity of the lateral deviation, Represents the rate of change of drift flow velocity.

[0074] S6. Substituting the obtained turbulent stress term into the mean flow equation and the turbulent pulsating velocity equation of the turbulence to obtain the anisotropic turbulence model.

[0075] The above method is based on the NS equation, uses the turbulent viscosity coefficient constructed by the turbulent pulsation velocity and turbulent dissipation rate parameters, and establishes an anisotropic turbulence model. The following is a specific implementation using a flat plate example. The implementation process is as follows: Figure 1 As shown, the following steps are included:

[0076] (1) Establishing geometric model - zero pressure gradient plate

[0077] The 3D model of the flat panel case is as follows Figure 2 As shown, the flat plate model has M = 0.3, Re = 6 × 10 6 The grid size is 65×97; the length of the viscous insulation wall is 1 unit; a region with a length of x=0.333 is set at the front of the plate and set as a symmetric boundary condition; the grid domain height is about 1.0, and the minimum normal grid spacing of the wall is 1×10 6 , the grid stretching ratio is set to 1.18; the horizontal grid spacing is constant, and the spacing is set to ΔX = 0.0208333; the grid division of the flat plate case is as follows Figure 3 shown.

[0078] (2) Setting boundary conditions for flat plate boundary layer flow

[0079] Inflow boundary settings: Total temperature is The total pressure is Other flow variables are determined using one-dimensional characteristic boundary conditions.

[0080] Outflow boundary setting: linear extrapolation method is used.

[0081] Setting at the top boundary: adopt Riemann invariant one-dimensional characteristic boundary condition.

[0082] (3) Import the turbulence pulsation formula of the new turbulence model into Open FOAM software

[0083] Since Open FOAM is an application written in C++, the pulsation speed needs to be The transport equation (1), the turbulent dissipation rate ω equation (3) and the viscosity coefficient μ T Formula (6) is imported into OpenFOAM software in the form of code, and a separate folder needs to be created to store the formula code. Figure 4 The code form of the turbulence pulsation formula in Open FOAM is shown. Figure 5 The code form of the turbulent viscosity coefficient in Open FOAM is shown.

[0084] (4) Use the k-ωSST turbulence model provided by Open FOAM software to obtain the initial field conditions of ω

[0085] Set the initial conditions: the inlet velocity in the 0 file is set to U = 33.8 m / s, the temperature T is set to 298 K, the inlet pressure P = 0, and the ω value is given. The k-ωSST turbulence model is set in the Constant file, the time step is set to 5000 steps, and the time interval is 500, that is, the data is output and saved every 500 times.

[0086] (5) Apply the initial field conditions of ω obtained from the k-ωSST turbulence model to the initial conditions of the new turbulence model;

[0087] Since the new model constructs the ω equation, the initial parameters of the ω equation can be obtained through the existing k-ωSST turbulence model, so that the new turbulence model is calculated in a physical field close to the data, which will speed up the calculation and make it convenient.

[0088] (6) Use Open FOAM software to perform numerical simulation of the new turbulence model

[0089] Set initial conditions: The 0 file sets the boundary condition parameters, such as initial velocity U, drift velocity U2, temperature T, pressure P, viscosity coefficient mut, etc.; the constant file sets the turbulence model and thermal physical properties, such as Prandtl number Pr = 0.71, Cp = 1005, and the value of coefficient C. In the system file, set the time step to 5000, the time interval to 500, and set the output files of residuals, speed, etc. During the operation, the Open FOAM command can be used to obtain the data residual graph of parameters such as velocity and drift velocity during the calculation process, such as Figure 6 As shown, we can observe whether the model converges in the flow field and whether the results are reasonable.

[0090] (7) Output the velocity through the Open FOAM command and make an average velocity profile

[0091] The Open FOAM command was used to extract the calculated velocity, and the Y+ and U+ data results were obtained using the solution formula. The average velocity profile was drawn using Origin software. Figure 7 It is mainly a data comparison chart of the new turbulence model, the SpalartAllmaras model, the k-ωSST model and the theoretical value. From the figure, it can be seen that the data results of the new turbulence model are more consistent with the theoretical value. However, the gap between the results and the theoretical value is a little bit large, so it can be corrected by modifying the only artificial parameter - coefficient C. Figure 8 As can be seen in the figure, for the flat plate case, it is found that when the coefficient C is in the range of 0.5 to 6.0, the fitting effect between the result and the theoretical value is relatively perfect.

[0092] The calculations for the flat plate case yielded the following conclusions:

[0093] Through numerical simulation calculations of the flat plate case, the calculation results of the anisotropic turbulence model equation constructed with only a single artificial parameter were compared with the SA turbulence model, k-ωSST turbulence model and theoretical values. It was found that the calculation results were close to the theoretical values, and the calculation process was simple and the calculation time was short.

[0094] By adjusting the only artificial parameter - coefficient C, the numerical results can be better fitted with the theoretical values, and the coefficient C can be adjusted within a certain range, and the data will be more accurate. The new turbulence model has fewer artificial parameters, which can greatly reduce the experimental errors caused by human factors, making the results more accurate and the data more convincing.

[0095] Embodiment 2: This embodiment provides a non-transitory computer-readable storage medium having computer instructions stored thereon, the computer instructions causing the computer to execute an anisotropic turbulence modeling method based on turbulent pulsation velocity and dissipation rate, the method comprising the following steps:

[0096] Establishing the pulsation speed The transport equation of

[0097] Solve the turbulent kinetic energy K equation and the turbulent dissipation rate ω equation to obtain the value of the turbulent dissipation rate ω in the flow field;

[0098] Solve for the pulsation speed The transport equation of turbulent flow is used to obtain the turbulent pulsation velocity

[0099] To obtain the turbulent dissipation rate ω and turbulent pulsation velocity Constructing the viscosity coefficient μ of turbulence T ;

[0100] Using the constructed turbulent viscosity coefficient μ T Calculate various turbulent stress terms;

[0101] The obtained turbulent stress term is substituted into the mean flow equation of turbulence and the turbulent fluctuation velocity equation to obtain the anisotropic turbulence model.

[0102] Embodiment 3: This embodiment provides an electronic device, which may include: a processor, a communications interface, a memory, and a communications bus, wherein the processor, the communications interface, and the memory communicate with each other via the communications bus. The processor may call logic instructions in the memory to execute an anisotropic turbulence modeling method based on turbulent pulsation velocity and dissipation rate, the method comprising the following steps:

[0103] Establishing the pulsation speed The transport equation of

[0104] Solve the turbulent kinetic energy K equation and the turbulent dissipation rate ω equation to obtain the value of the turbulent dissipation rate ω in the flow field;

[0105] Solve for the pulsation speed The transport equation of turbulent flow is used to obtain the turbulent pulsation velocity

[0106] To obtain the turbulent dissipation rate ω and turbulent pulsation velocity Constructing the viscosity coefficient μ of turbulence T ;

[0107] Using the constructed turbulent viscosity coefficient μ T Calculate various turbulent stress terms;

[0108] The obtained turbulent stress term is substituted into the mean flow equation of turbulence and the turbulent fluctuation velocity equation to obtain the anisotropic turbulence model.

[0109] In addition, the logic instructions in the above-mentioned memory can be implemented in the form of software functional units and can be stored in a computer-readable storage medium when sold or used as an independent product. Based on such an understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art or the part of the technical solution, can be embodied in the form of a software product, which is stored in a storage medium and includes several instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes: various media that can store program codes, such as a USB flash drive, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk or an optical disk.

[0110] Embodiment 4: This embodiment provides a computer program product, which includes a computer program. The computer program can be stored in a non-transitory computer-readable storage medium. When the computer program is executed by a processor, the computer can execute an anisotropic turbulence modeling method based on turbulent pulsation velocity and dissipation rate. The method includes the following steps:

[0111] Establishing the pulsation speed The transport equation of

[0112] Solve the turbulent kinetic energy K equation and the turbulent dissipation rate ω equation to obtain the value of the turbulent dissipation rate ω in the flow field;

[0113] Solve for the pulsation speed The transport equation of turbulent flow is used to obtain the turbulent pulsation velocity

[0114] To obtain the turbulent dissipation rate ω and turbulent pulsation velocity Constructing the viscosity coefficient μ of turbulence T ;

[0115] Using the constructed turbulent viscosity coefficient μ T Calculate various turbulent stress terms;

[0116] The obtained turbulent stress term is substituted into the mean flow equation of turbulence and the turbulent fluctuation velocity equation to obtain the anisotropic turbulence model.

[0117] The device embodiments described above are merely illustrative, wherein the units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the scheme of this embodiment. Ordinary technicians in this field can understand and implement it without paying creative labor.

[0118] Through the description of the above implementation methods, those skilled in the art can clearly understand that each implementation method can be implemented by means of software plus a necessary general hardware platform, and of course, it can also be implemented by hardware. Based on this understanding, the above technical solution is essentially or the part that contributes to the prior art can be embodied in the form of a software product, and the computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, a disk, an optical disk, etc., including a number of instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to execute the methods described in each embodiment or some parts of the embodiments.

[0119] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for modeling anisotropic turbulence model based on turbulent pulsation velocity and dissipation rate, characterized in that: include: Establishing the pulsation speed The transport equation of Solve the turbulent kinetic energy K equation and the turbulent dissipation rate ω equation to obtain the value of the turbulent dissipation rate ω in the flow field; Solve for the pulsation speed The transport equation of turbulent flow is used to obtain the turbulent pulsation velocity To obtain the turbulent dissipation rate ω and turbulent pulsation velocity Constructing the viscosity coefficient μ of turbulence T ; Using the constructed turbulent viscosity coefficient μ T Calculate various turbulent stress terms; The obtained turbulent stress term is substituted into the mean flow equation of turbulence and the turbulent fluctuation velocity equation to obtain the anisotropic turbulence model.

2. The anisotropic turbulence modeling method based on turbulent pulsation velocity and dissipation rate according to claim 1 is characterized in that: The established pulsation speed The transport equation of , the specific method is: The turbulent pulsating velocity in the flow field is divided into two groups, the average value of each group is not zero, and the pulsating velocity is established. The transport equation is expressed as: Where ρ represents the fluid density, u i 、u j represents the average flow velocity, represents the drift flow velocity, τ ij represents the mean flow turbulent viscous stress, represents the turbulent viscous stress of the drift flow, t is time, x i is the space coordinate.

3. The anisotropic turbulence modeling method based on turbulent pulsation velocity and dissipation rate according to claim 1 is characterized in that: The turbulent kinetic energy K equation and the turbulent dissipation rate ω equation are solved to obtain the value of the turbulent dissipation rate ω in the flow field; the specific method is: The turbulent kinetic energy K equation is established and the expression is: The turbulent dissipation rate ω equation is established, and the expression is: Among them, v t represents the eddy kinematic viscosity coefficient, Ω is the vorticity, a1=0.31; ρ represents the fluid density, K represents the turbulent kinetic energy, ω represents the turbulent dissipation rate, β, σ ω , σ ω2 , β * , γ represents the model constant, ψ is the specific heat ratio, ψ = 0.41, σ ω1 =0.5,β=0.075,β * =0.09,σ ω2 =0.856; μ, μ i represents the laminar viscosity coefficient and turbulent viscosity coefficient, F1 and F2 are mixed functions, which are related to the distance from the point to the wall; D is the total differential symbol, t is time, x j is the space coordinate.

4. The anisotropic turbulence modeling method based on turbulent pulsation velocity and dissipation rate according to claim 3 is characterized in that: The expression of the mixing function F1 is: The expression of the mixed function F2 is: Among them, CD kω It represents the positive direction of the orthogonal diffusion term added to the standard k-ω turbulence model and the standard k-ε turbulence model. y represents the minimum distance from a point in the flow field to the surface of an object, and ε represents the turbulent energy dissipation rate.

5. The anisotropic turbulence modeling method based on turbulent pulsation velocity and dissipation rate according to claim 1 is characterized in that: Use turbulent pulse speed and turbulence dissipation rate ω, the viscosity coefficient of constructed turbulence μ T , the expression is: Where ρ represents the density of the selected fluid, and the coefficient c is an artificial parameter, which is initially set to c=1.

0.

6. The anisotropic turbulence modeling method based on turbulent pulsation velocity and dissipation rate according to claim 1 is characterized in that: Using the constructed turbulent viscosity coefficient μ T The calculated turbulent stress terms include: Among them, μ T represents the viscous turbulence coefficient; μ L represents the viscous laminar coefficient; S represents the deformation rate, δ ij stands for Kronecker operator; represents the turbulent viscous stress of the drift flow, represents the turbulent strain rate of the drift flow, is the average pulsating velocity of the lateral deviation, Represents the rate of change of drift flow velocity.

7. A non-transitory computer-readable storage medium having computer instructions stored thereon, the computer instructions enabling a computer to execute the anisotropic turbulence modeling method based on turbulent pulsation velocity and dissipation rate as described in any one of claims 1 to 6.

8. An electronic device comprising: A processor, a communication interface, a memory and a communication bus, wherein the processor, the communication interface and the memory communicate with each other through the communication bus, and the processor calls the logic instructions in the memory to execute the anisotropic turbulence model modeling method based on turbulent pulsation velocity and dissipation rate as described in any one of claims 1-6.

9. A computer program product, comprising a computer program stored on a non-transitory computer-readable storage medium, wherein when the computer program is executed by a processor, the computer executes the anisotropic turbulence model modeling method based on turbulent pulsation velocity and dissipation rate as described in any one of claims 1 to 6.

Images