Modeling Method of Anisotropic Turbulence Model Based on Turbulent Fluctuation Velocity and Dissipation Rate

Abstract

The present invention discloses a method for modeling an anisotropic turbulence model based on turbulent pulsating velocity and dissipation rate, including: establishing a transport equation for pulsating velocity; solving the turbulent kinetic energy equation and the turbulent dissipation rate equation to obtain the numerical value of the turbulent dissipation rate in the flow field; solving the transport equation for the pulsating velocity to obtain the turbulent pulsating velocity; constructing the viscous coefficient of turbulence with the obtained turbulent dissipation rate and turbulent pulsating velocity; calculating each stress term of turbulence by using the constructed turbulent viscous coefficient; substituting the obtained turbulent stress terms into the mean flow equation of turbulence and the turbulent pulsating velocity equation to obtain the anisotropic turbulence model. The anisotropic turbulence model equation established by the present invention improves the numerical calculation accuracy of the turbulent boundary layer, and the constructed anisotropic turbulence model equation only contains a single artificial parameter, which can greatly reduce the experimental error caused by human factors, make the results more accurate, and the data 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] In this regard, 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 turbulent model equation, further improving the numerical calculation accuracy of the turbulent boundary layer.

[0007] To achieve the above object, on the one hand, the present invention provides an anisotropic turbulent model modeling method based on turbulent pulsation velocity and dissipation rate, including:

[0008] Establish the transport equation of the pulsation velocity ;

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

[0010] Solve the transport equation of the pulsation velocity to obtain the turbulent pulsation velocity

[0011] Using the obtained turbulent dissipation rate ω and turbulent pulsation velocity Construct the viscous coefficient μ of the turbulence T ;

[0012] Use the constructed turbulent viscous coefficient μ T to calculate each stress term of the turbulence;

[0013] Substitute the obtained turbulent stress terms into the mean flow equation of the turbulence and the turbulent pulsation velocity equation to obtain the anisotropic turbulent model.

[0014] A further preferred technical solution of the present invention is that the method for establishing the transport equation of the pulsation velocity is specifically as follows:

[0015] Divide the turbulent pulsation velocity in the flow field into two groups, and the average value of each group is not zero. Establish the transport equation of the pulsation velocity , and the expression is:

[0016]

[0017] where ρ represents the fluid density, u i , u j represent the mean flow velocity, represents the drift flow velocity, τ ij represents the mean flow turbulent viscous stress, represents the drift flow turbulent viscous stress, t is time, and x i is the spatial coordinate.

[0018] Preferably, the method for solving the turbulent kinetic energy K equation and the turbulent dissipation rate ω equation to obtain the numerical value of the turbulent dissipation rate ω in the flow field is specifically as follows:

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

[0020]

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

[0022]

[0023] where, 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 、β * 、γ represent model constants, ψ is the specific heat ratio, ψ = 0.41, σ ω1 = 0.5, β = 0.075, β * = 0.09, σ ω2 = 0.856; μ, μ i represent the laminar viscosity coefficient and the turbulent viscosity coefficient respectively, F1 and F2 are mixing functions, which are related to the distance from the point to the wall; D is the total differential symbol, t is the time, and x j is the spatial coordinate.

[0024] Preferably, the expression of the mixing function F1 is:

[0025]

[0026] The expression of the mixing function F2 is:

[0027]

[0028] where, CD kω represents the positive direction of the orthogonal diffusion term added on the basis of 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 object surface, and ε represents the turbulent energy dissipation rate.

[0029] Preferably, the turbulent viscosity coefficient μ constructed using the turbulent fluctuating velocity T and the turbulent dissipation rate ω has the following expression:

[0030]

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

[0032] Preferably, the constructed turbulent viscosity coefficient μ is used T to calculate the turbulent stress terms, including:

[0033]

[0034] where μ T represents the viscous turbulent coefficient; μ L represents the viscous laminar coefficient; S represents the deformation rate, and δ ij represents the Kronecker operator; represents the drift flow turbulent viscous stress, represents the drift flow turbulent strain rate, is the lateral mean pulsating velocity, represents the deformation rate of the drift flow velocity.

[0035] On the other hand, the present invention provides a non-transitory computer-readable storage medium, on which computer instructions are stored, and the computer instructions cause the computer to execute the above-mentioned anisotropic turbulence model modeling method based on turbulent pulsating velocity and dissipation rate.

[0036] On another aspect, the present invention provides an electronic device, including: a processor, a communication interface, a memory, and a communication bus. Among them, the processor, the communication interface, and the memory complete mutual communication through the communication bus, and the processor calls the logical instructions in the memory to execute the above-mentioned anisotropic turbulence model modeling method based on turbulent pulsating velocity and dissipation rate.

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

[0038] Advantageous effects: The present invention constructs a turbulent viscosity coefficient using turbulent pulsating velocity and turbulent dissipation rate parameters, and establishes an anisotropic turbulence model equation, which improves the numerical calculation accuracy of the turbulent boundary layer. Moreover, the constructed anisotropic turbulence model equation only contains a single artificial parameter, which can greatly reduce the experimental error caused by human factors, make the results more accurate, and the data more persuasive. Description of the Drawings

[0039] Figure 1 is the flow chart of the construction and application of the anisotropic turbulence model in the embodiment;

[0040] Figure 2 is the three-dimensional view of the geometric model - zero pressure gradient flat plate established in the embodiment;

[0041] Figure 3Mesh generation diagram of the geometric model - zero pressure gradient flat plate in the embodiment;

[0042] Figure 4 Code form diagram of the pulsating velocity equation in the newly established turbulence model in the embodiment;

[0043] Figure 5 Code form diagram of the turbulent viscosity coefficient in the newly established turbulence model in the embodiment;

[0044] Figure 6 Residual convergence diagram of velocity, pulsating velocity, etc. of the newly established turbulence model in the embodiment;

[0045] Figure 7 Dimensionless mean velocity profile comparison diagram of the newly established turbulence model in the embodiment;

[0046] Figure 8 Dimensionless mean velocity profile comparison diagram of the coefficient C correction of the newly established turbulence model in the embodiment. Detailed implementation method

[0047] To make the objectives, technical solutions and advantages of the present invention clearer, the technical solutions in the present invention will be clearly and completely described below in conjunction with the accompanying drawings in the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, rather than all of them, and they should not be construed as limiting the present invention. Based on the embodiments in the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts fall 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 the purpose of description and cannot be construed as indicating or implying relative importance.

[0048] The following combines Figures 1-8 Describe the anisotropic turbulence model modeling method based on turbulent pulsating velocity and dissipation rate provided by the present invention.

[0049] Embodiment 1: This embodiment first provides an anisotropic turbulence model modeling method based on turbulent pulsating velocity and dissipation rate to establish a new turbulence model. Then, the newly established turbulence model is applied to actual applications and numerical simulation calculations are carried out.

[0050] The anisotropic turbulence model modeling method based on turbulent pulsating velocity and dissipation rate, as Figure 1 shown, includes the steps:

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

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

[0053]

[0054] where ρ represents the fluid density, u i and u j represent the mean flow velocity, represents the drift flow velocity, τ ij represents the mean flow turbulent viscosity stress, represents the drift flow turbulent viscosity stress, t is time, and x i is the spatial coordinate.

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

[0056] Establish the turbulent kinetic energy K equation, which is expressed as:

[0057]

[0058] Establish the turbulent dissipation rate ω equation, which is expressed as:

[0059]

[0060] where 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, β, σ ω and σ ω2 and β * and γ represent model constants, ψ is the specific heat ratio, ψ = 0.41, σ ω1 = 0.5, β = 0.075, β * = 0.09, σ ω2 = 0.856; μ, μ i represent the laminar viscosity coefficient and the turbulent viscosity coefficient, F1 and F2 are mixing functions, which are related to the distance from the point to the wall; D is the total differential symbol, t is time, and x j is the spatial coordinate.

[0061] The expression of the mixing function F1 is:

[0062]

[0063] The expression of the mixing function F2 is:

[0064]

[0065] where CD kωIndicates the positive direction of the orthogonal diffusion term added on the basis of 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 object surface, and ε represents the turbulent energy dissipation rate.

[0066] S3. Solve the transport equation of the pulsating velocity to obtain the turbulent pulsating velocity

[0067] S4. Construct the turbulent viscosity coefficient μ of the turbulence with the obtained turbulent dissipation rate ω and turbulent pulsating velocity ; the constructed turbulent viscosity coefficient μ T ; the expression is: T where ρ represents the density of the selected fluid, and the coefficient c is an artificial parameter, initially set to c = 1.0.

[0068]

[0069]

[0070] S5. Use the constructed turbulent viscosity coefficient μ T to calculate each stress term of the turbulence, including:

[0071]

[0072]

[0073] where μ T represents the viscous turbulent coefficient; μ L represents the viscous laminar coefficient; S represents the deformation rate, and δ ij represents the Kronecker operator; represents the drift flow turbulent viscous stress, represents the drift flow turbulent strain rate, is the lateral deviation mean pulsating velocity, represents the deformation rate of the drift flow velocity.

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

[0075] Figure 1 The above method is based on the N-S equation, constructs the turbulent viscosity coefficient by using the turbulent pulsating velocity and the turbulent dissipation rate parameters, and establishes the anisotropic turbulence model. The following takes a flat plate example for specific implementation, and the implementation process is as shown, including the following steps:

[0076] (1). Establish a geometric model - zero pressure gradient flat plate (1). Establish a geometric model - zero pressure gradient flat plate

[0077] The three-dimensional model of the flat plate case is as Figure 2As 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, enabling the new turbulence model to be calculated in a physical field close to the data. The calculation speed will be fast and the calculation is also convenient.

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

[0089] Set the initial conditions: Set the boundary condition parameters in the 0 file, such as the initial velocity U, drift velocity U2, temperature T, pressure P, viscosity coefficient mut, etc.; Set the turbulence model and thermophysical properties in the constant file, such as the Prandtl number Pr = 0.71, Cp = 1005, and the magnitude of the coefficient C value. Set the number of time steps to 5000, the time interval to 500, and set the output files for residuals, velocity, etc. in the system file. During the running process, data residual plots of parameters such as velocity and drift velocity during the calculation can be obtained through Open FOAM commands, as Figure 6 shown, and it can be observed whether the model converges in the flow field and whether the results are reasonable.

[0090] (7), Output the velocity through Open FOAM commands and create an average velocity profile

[0091] Use the Open FOAM command to extract the calculated velocity, and use the solution formula to obtain the data results of Y+ and U+. Use the Origin software to plot the average velocity profile. Figure 7 Mainly it is a data comparison graph of the new turbulence model with the SpalartAllmaras model, k-ω SST model, and theoretical values. It can be seen from the graph that the data results of the new turbulence model are relatively consistent with the theoretical values. However, the gap between the results and the theoretical values is slightly large, so it can be corrected by modifying the only artificial parameter - the coefficient C. From Figure 8 it can be seen that for the flat plate case, it is found that when the value of the coefficient C is in the range of 0.5 - 6.0, the fitting effect of the results with the theoretical values is relatively perfect.

[0092] Calculations were carried out on the flat plate case, and the following conclusions can be obtained:

[0093] Through the numerical simulation calculations of the flat plate case, comparing the results calculated by the anisotropic turbulence model equation containing only the unique artificial parameter with the SA turbulence model, k-ω SST turbulence model, and theoretical values, it is found that the calculated results are close to the theoretical values, and the calculation process is simple and the calculation time is short.

[0094] By adjusting the only artificial parameter - coefficient C, the fitting effect between the numerical result and the theoretical value can be better, and the coefficient C can be adjusted within a certain range, making the data more accurate. The new turbulence model has fewer artificial parameters, which can greatly reduce the experimental error caused by human factors, making the results more accurate and the data more persuasive.

[0095] Example 2: This example provides a non - transient computer - readable storage medium, on which computer instructions are stored. These computer instructions enable a computer to execute an anisotropic turbulence model modeling method based on turbulent pulsation velocity and dissipation rate. The method includes the following steps:

[0096] Establish the transport equation of the pulsation velocity ;

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

[0098] Solve the transport equation of the pulsation velocity to obtain the turbulent pulsation velocity

[0099] With the obtained turbulent dissipation rate ω and turbulent pulsation velocity Construct the turbulent viscosity coefficient μ of the turbulence T ;

[0100] Use the constructed turbulent viscosity coefficient μ T to calculate each stress term of the turbulence;

[0101] Substitute the obtained turbulent stress terms into the mean flow equation of the turbulence and the turbulent pulsation velocity equation to obtain the anisotropic turbulence model.

[0102] Example 3: This example provides an electronic device, which may include: a processor, a communication interface, a memory, and a communication bus. Among them, the processor, the communication interface, and the memory complete mutual communication through the communication bus. The processor can call the logical instructions in the memory to execute an anisotropic turbulence model modeling method based on turbulent pulsation velocity and dissipation rate. The method includes the following steps:

[0103] Establish the transport equation of the pulsation velocity ;

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

[0105] Solve the transport equation of the pulsation velocity to obtain the turbulent pulsation velocity

[0106] to obtain the turbulent dissipation rate ω and the turbulent fluctuating velocity Construct the turbulent viscosity coefficient μ of the turbulence T ;

[0107] Utilize the constructed turbulent viscosity coefficient μ T Calculate each stress term of the turbulence;

[0108] Substitute the obtained turbulent stress terms into the mean flow equation of the turbulence and the turbulent fluctuating velocity equation to obtain the anisotropic turbulence model.

[0109] In addition, when the logical instructions in the above-mentioned memory are implemented in the form of software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on such an understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art, or a part of this technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods described in various embodiments of the present invention. The foregoing storage medium includes: various media such as USB flash drives, mobile hard disks, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), magnetic disks, or optical discs that can store program codes.

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

[0111] Establish the transport equation of the fluctuating velocity ;

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

[0113] Solve the transport equation of the fluctuating velocity to obtain the turbulent fluctuating velocity

[0114] to obtain the turbulent dissipation rate ω and the turbulent fluctuating velocity Construct the turbulent viscosity coefficient μ of the turbulence T ;

[0115] Utilize the constructed turbulent viscosity coefficient μ TCalculate each stress term of the turbulence;

[0116] Substitute the obtained turbulence stress terms into the mean flow equation of the turbulence and the turbulence pulsation velocity equation to obtain the anisotropic turbulence model.

[0117] The device embodiments described above are merely illustrative. The units described as separation components may or may not be physically separated. The components shown as units may or may not be physical units, that is, they may be located in one place or distributed to multiple network units. Some or all of the modules can be selected according to actual needs to achieve the purpose of the solution of this embodiment. Those of ordinary skill in the art can understand and implement it without creative efforts.

[0118] Through the description of the above embodiments, those skilled in the art can clearly understand that each embodiment can be implemented by means of software plus a necessary general hardware platform, and of course also by hardware. Based on this understanding, the above technical solutions, in essence, or the part that contributes to the prior art can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., and includes several instructions for causing 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 and are not intended to limit them. Although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that they can still modify the technical solutions described in the foregoing embodiments or perform equivalent replacements for some of the technical features. These modifications or replacements do not cause the essence of the corresponding technical solutions to deviate 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 ; specifically: 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: （1）； in, represents the fluid density, , represents the average flow velocity, , represents the drift flow velocity, represents the mean flow turbulent viscous stress, represents the turbulent viscous stress of the drift flow, t is the time, is the spatial coordinate; Solving for turbulent kinetic energy Equations and turbulent dissipation rates Equation, to obtain the turbulent dissipation rate in the flow field The value of; specifically: Building turbulent kinetic energy The equation is expressed as: （2）； Establishing the turbulent dissipation rate The equation is expressed as: （3）； in, represents the eddy kinematic viscosity coefficient, , is the vorticity, ; represents the fluid density, represents the turbulent kinetic energy, represents the turbulent dissipation rate, , , , , represents the model constant, , is the specific heat ratio, , , , , ; , represents the laminar viscosity coefficient and turbulent viscosity coefficient, and is a mixed function, which is related to the distance from the point to the wall; is the total differential symbol For time, is the spatial coordinate; 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 ; Using the constructed turbulent viscosity coefficient 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: Mixing functions The expression is: （4）； Mixing functions The expression is: （5）； in, It represents the positive direction of the orthogonal diffusion term added to the standard k-ω turbulence model and the standard k-ε turbulence model. , It represents the minimum distance from a point in the flow field to the surface of an object, and ε represents the turbulent energy dissipation rate.

3. 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 turbulent dissipation rate , the viscosity coefficient of the constructed turbulence , the expression is: （6）； in, Represents the density of the selected fluid, the coefficient It is an artificial parameter, and the initial setting is =1.

0.

4. 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 The calculated turbulent stress terms include: （7）； （8）； （9）； in, represents the viscous turbulence coefficient; represents the viscous laminar flow coefficient; represents the deformation rate, 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.

5. A non-transitory computer-readable storage medium having computer instructions stored thereon, wherein the computer instructions enable 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 4.

6. 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-4.

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

Images