A grid-adaptive turbulence simulation method based on Vreman dynamic coefficient coupled k-ω series model

Abstract

The present invention discloses a grid adaptive turbulence simulation method based on Vreman dynamic coefficient coupled k-ω series model, and the main process of its implementation includes: step 1, judging whether to apply shielding function; step 2, calculating Vreman dynamic coefficient; step 3, determining local dynamic grid length scale; step 4, constructing scale-related adjustment function; step 5, reconstructing turbulent viscosity of k-ω series model; step 6, using reconstructed turbulent viscosity to simulate turbulence. The present invention constructs Vreman dynamic coefficient based on local velocity gradient tensor, identifies local grid size, and determines local grid length scale Δ * , and then reconstruct the turbulent viscosity by constructing scale-related functions through turbulent energy spectrum integration to achieve grid adaptive simulation, effectively overcoming the problem of high empirical dependence of existing turbulent mixing models on grids, further broadening the grid adaptability, effectively improving the calculation accuracy, greatly reducing the calculation cost, and significantly accelerating the turbulence simulation process.

Description

Technical Field

[0001] The invention relates to the field of complex flow calculation of aircraft engines and gas turbines, and in particular to a grid adaptive turbulence simulation method based on Vreman dynamic coefficient coupled k-ω series models. Background Art

[0002] Aircraft engines and gas turbines are a major national demand and a key development area for my country's high-end equipment manufacturing. Turbulence is ubiquitous in nature and engineering, and accurate prediction of turbulence is a major challenge in studying complex flow problems. Turbulent flows in aircraft engines and gas turbines are extremely complex, with multi-scale, nonlinear, and unsteady characteristics. These complex turbulences have a great impact on the performance of aircraft engines and gas turbines. Therefore, in engineering design, it is urgent to develop turbulence simulation methods with high prediction accuracy and high computational efficiency, so that they can accurately simulate turbulent flows in complex engineering problems such as the internal flow of aircraft engines and gas turbines.

[0003] The turbulence simulation method commonly used in existing engineering projects is mainly based on solving the Reynolds-averaged NS equations (RANS) method. Although the amount of calculation is small, it is based on simple basic flows and has poor prediction accuracy for complex turbulent flows such as multi-scale, unsteady, and large separation in engineering problems. Therefore, it is difficult to achieve a detailed study of the complex flow mechanism in engineering problems such as the internal flow of aircraft engines and gas turbines, which seriously restricts the improvement of the design level in the fields of aircraft engines and gas turbines. As a high-precision numerical simulation method, the large eddy simulation (LES) method has a high requirement for the number of grids, which makes its computational cost increase exponentially compared to the RANS method, which is much higher than the computational cost that can be tolerated in engineering applications. For engineering flow problems such as the internal flow of aircraft engines and gas turbines, the Reynolds number is often high. Under the existing computing power, the LES method is still difficult to apply to the flow prediction in complex engineering fields at the engineering design level.

[0004] In the past two decades, the RANS-LES hybrid simulation method was born. By using the RANS method in the near-wall area and the LES method in the mainstream area, the calculation accuracy and efficiency are balanced, which provides a good solution to the problem of high-precision simulation and high computational cost of complex flows. However, the current classic RANS-LES hybrid model has relatively strict requirements on the grid, requiring users to have rich experience in high-precision numerical simulation. The existing RANS-LES hybrid simulation method is highly dependent on the experience-based grid distribution for the prediction accuracy of complex engineering flows such as the internal flow of aircraft engines and gas turbines, and the grid demand is large, making it difficult to apply to daily engineering design. Therefore, reducing the empirical dependence of the RANS-LES hybrid model on the grid is of great significance to better realize the accurate, efficient and high-precision simulation of multi-scale, nonlinear, unsteady and other turbulent flow phenomena in complex engineering flow problems.

[0005] The present invention further improves the prediction accuracy and grid adaptability of the grid adaptive turbulence simulation method by introducing the Vreman dynamic coefficient. Compared with the disclosed invention patent CN 115186607 A, the prediction accuracy of the grid adaptive turbulence simulation method can be further improved and the calculation cost can be greatly reduced. Summary of the invention

[0006] 1. Technical issues to be resolved

[0007] The purpose of the present invention is to propose a grid-adaptive turbulence simulation method based on the Vreman dynamic coefficient coupled k-ω series model, which further improves the grid adaptability by introducing the Vreman dynamic coefficient, effectively overcomes the problem of the existing RANS-LES hybrid model's high empirical dependence on the grid, and while improving the calculation accuracy, greatly reduces the calculation cost, significantly speeds up the turbulence simulation process, and more quickly stimulates the unsteady characteristics of the flow, providing an efficient numerical simulation method for solving the fast and high-precision simulation of multi-scale, nonlinear, unsteady and other turbulent flows in complex engineering flow problems.

[0008] (II) Technical solution

[0009] In order to solve the above technical problems, the present invention provides a grid adaptive turbulence simulation method based on Vreman dynamic coefficient coupled k-ω series model, comprising the following steps:

[0010] Step 1, determine whether to apply the shielding function;

[0011] Step 2, calculate the Vreman dynamic coefficient;

[0012] Step 3, determine the local dynamic grid length scale;

[0013] Step 4: construct a scale-dependent adjustment function;

[0014] Step 5, reconstruct the turbulent viscosity of the k-ω series model;

[0015] Step 6, use the reconstructed turbulent viscosity to perform turbulence simulation;

[0016] ① The determination of whether to apply the shielding function includes:

[0017] Combined with the type of flow state being simulated, determine whether to use the shielding function F GAS Specifically, when the flow state type is free shear flow, the shielding function is not used. In this case, the shielding function F GAS = 0; when the flow state type is near-wall flow, a shielding function is used, and the shielding function F GAS Select one of the F1 shielding function and F2 shielding function from the DDES-SST model and the F from the DDES-SA model. d Masking function;

[0018] ② The calculation of the Vreman dynamic coefficient includes:

[0019] Based on the local velocity gradient tensor α ij , construct the Vreman dynamic coefficient C DGAS ;

[0020] The Vreman dynamic coefficient C DGAS It is obtained from the following formula:

[0021]

[0022] β ij =α mi α mj

[0023]

[0024] Where S is the strain rate, C GAS The empirical coefficient is taken as 0.6;

[0025] ③ The determination of the local dynamic grid length scale includes:

[0026] Combined with the shielding function F described in step 1 GAS , and the Vreman dynamic coefficient C described in step 2 DGAS , determine the local dynamic grid length scale Δ * ;

[0027] The local dynamic grid length scale Δ * Given by:

[0028] Δ * =C DGAS [(1-F GAS )Δ vol +F GAS Δ max ]

[0029] Δ max =max(Δ x ,Δ y ,Δ z )

[0030]

[0031] Among them, Δ x is the length of the local hexahedral grid, Δ y is the width of the local hexahedral grid, Δ z is the height of the local hexahedral grid;

[0032] ④ The scale-related adjustment function constructed based on the turbulence energy spectrum integral coupling k-ω series model includes:

[0033] According to the modeling method of turbulent kinetic energy in the k-ω series model, the original modeled turbulent kinetic energy k is obtained. m , according to the local dynamic grid length scale Δ * Based on the turbulence energy spectrum, the actual simulated turbulent kinetic energy k is obtained by integration u ;

[0034] The actual simulated turbulent kinetic energy k u It is obtained from the following formula:

[0035]

[0036] Among them, C k is the Kolmogorov constant coefficient, which is taken as 1.5; ε is the actual turbulence dissipation rate; κ c To solve the turbulent cutoff wave number, the local dynamic grid length scale Δ * Determine as follows:

[0037]

[0038] Among them, π is the ratio of pi, which is 3.14;

[0039] According to the actual simulated turbulent kinetic energy k u , the original modeled turbulent kinetic energy k m and the shielding function F described in step 1 GAS , construct the dynamic scale-dependent regulation function D f ; Define the scale ratio as the actual simulated turbulent kinetic energy ku and the original modeled turbulent kinetic energy k m The dynamic scale-dependent regulation function D f is a function related to the ratio of the scales, which is obtained by the following formula:

[0040]

[0041] l GAS =(1-F GAS ) u +F GAS l m

[0042]

[0043] Among them, l u is the grid-related scale, l m is the turbulence length scale given by the k-ω series model, which is obtained from the following equation:

[0044]

[0045] Among them, ω m is the specific dissipation rate of the original model obtained by the k-ω series model; β * is the constant coefficient in the k-ω series model, which is taken as 0.09;

[0046] ⑤The turbulent viscosity of the reconstructed k-ω series model includes:

[0047] Use the dynamic scale-dependent adjustment function D described in step 4 f For the turbulent viscosity ν in the k-ω series model t By adjusting the turbulent viscosity ν, we can obtain the reconstructed turbulent viscosity ν. sfs , obtained from the following formula:

[0048] ν sfs =D f ·ν t

[0049] ⑥ The turbulence simulation using the reconstructed turbulent viscosity includes:

[0050] Using the reconstructed turbulent viscosity ν described in step 5 sfs , calculate the Reynolds stress, and update the transport equation of the k-ω series model, replacing the turbulent viscosity ν in the original transport equation of the k-ω series model t , combined with the k-ω series model, the grid adaptive turbulence simulation method based on Vreman dynamic coefficient coupling k-ω series model is obtained.

[0051] (III) Beneficial effects

[0052] The grid adaptive turbulence simulation method based on Vreman dynamic coefficient coupled k-ω series model provided by the present invention has the following beneficial effects: constructing Vreman dynamic coefficient based on local velocity gradient tensor, determining local dynamic grid length scale Δ by identifying local grid size * , and then reconstruct the turbulent viscosity by constructing the scale-related function through the turbulent energy spectrum integration to achieve grid adaptive simulation, which effectively overcomes the problem of high empirical dependence of the existing RANS-LES hybrid model on the grid.

[0053] Compared with the prior art CN 115186607 A, the grid adaptability of the grid adaptive turbulence simulation method is further broadened by introducing the Vreman dynamic coefficient. While improving the calculation accuracy, the calculation cost is greatly reduced, and the turbulence simulation process is significantly accelerated. This provides an efficient numerical simulation method for fast and high-precision simulation of multi-scale, nonlinear, unsteady and other turbulent flows in complex engineering flow problems.

[0054] The method of the present invention is concise in form and highly portable, can be well combined with the k-ω series model, is easy to be implanted into the existing CFD code for application and expansion, and has broad academic and engineering application prospects. BRIEF DESCRIPTION OF THE DRAWINGS

[0055] Figure 1 It is a flow chart of a grid adaptive turbulence simulation method based on Vreman dynamic coefficient coupled k-ω series model of the present invention;

[0056] Figure 2 A dynamic scale-dependent adjustment function D is constructed for a grid adaptive turbulence simulation method based on Vreman dynamic coefficient coupled k-ω series model of the present invention. f Schematic diagram of

[0057] Figure 3 Schematic diagram of the turbulent eddy structure of the back-step example calculated by the SST k-ω model using the k-ω series model;

[0058] Figure 4 The figure is a schematic diagram of the turbulent vortex structure of a back-step example calculated by a grid adaptive turbulence simulation method based on the Vreman dynamic coefficient coupled k-ω series model of the present invention. DETAILED DESCRIPTION

[0059] In the following, in conjunction with the accompanying drawings, the SST k-ω turbulence model in the k-ω series model is taken as an example to further describe the specific implementation of the present invention in detail. The following examples are used to illustrate the present invention but are not used to limit the scope of the present invention.

[0060] The back step example is selected as the calculation example, and the high-quality hexahedral mesh is used to spatially discretize the calculation domain. The total number of meshes is about 410,000.

[0061] The present invention provides a grid adaptive turbulence simulation method based on Vreman dynamic coefficient coupled k-ω series model, comprising the following steps:

[0062] Step 1, determine whether to apply the shielding function;

[0063] In this step, the type of flow state being simulated is considered to determine whether to use the screening function F GAS Specifically, when the flow state type is free shear flow, the shielding function is not used. In this case, the shielding function F GAS = 0; when the flow state type is near-wall flow, a shielding function is used, and the shielding function F GAS Select one of the F1 shielding function and F2 shielding function from the DDES-SST model and the F from the DDES-SA model. d Masking function.

[0064] In this embodiment, the F2 function from the DDES-SST model is used as an example, and then:

[0065] F GAS =F2

[0066] F2=tanh(ψ 2 )

[0067]

[0068] Among them, k m is the turbulent kinetic energy of the original model of the k-ω series model, ω m is the specific dissipation rate of the original model obtained by the k-ω series model, ν is the viscosity coefficient of the fluid, d is the distance from the local grid to the wall, β * is a constant coefficient in the SST k-ω turbulence model and is taken as 0.09.

[0069] Step 2, calculate the Vreman dynamic coefficient;

[0070] In this step, based on the local velocity gradient tensor α ij , construct the Vreman dynamic coefficient C DGAS ;

[0071] The Vreman dynamic coefficient C DGAS It is obtained from the following formula:

[0072]

[0073]

[0074] β ij =α mi α mj

[0075]

[0076] Where S is the strain rate, C GAS The empirical coefficient is taken as 0.6.

[0077] Step 3, determine the local dynamic grid length scale;

[0078] In this step, combined with the shielding function F in step 1 GAS , and the Vreman dynamic coefficient C in step 2 DGAS , we can get the local dynamic grid length scale Δ * ;

[0079] Local dynamic grid length scale Δ * Given by:

[0080] Δ * =C DGAS [(1-F GAS )Δ vol +F GAS Δ max ]

[0081] Δ max =max(Δ x ,Δ y ,Δ z )

[0082]

[0083] Where Δ x is the length of the local hexahedral grid, Δ y is the width of the local hexahedral grid, Δ z is the height of the local hexahedral grid.

[0084] Step 4: construct a scale-dependent adjustment function;

[0085] In this step, according to the modeling method of turbulent kinetic energy in the k-ω series turbulence model, the original modeled turbulent kinetic energy k is obtained. m In this embodiment, taking the SST k-ω turbulence model in the k-ω series model as an example, according to the modeling method of the turbulent kinetic energy in the SST k-ω turbulence model, the original modeled turbulent kinetic energy k is obtained. m ;

[0086] According to the local grid length scale, the actual simulated turbulent kinetic energy k is obtained by integration based on the turbulence energy spectrum. u, is given by:

[0087]

[0088] Among them, C k is the Kolmogorov constant coefficient, which is taken as 1.5; ε is the actual turbulence dissipation rate; κ c To solve the turbulent cutoff wave number, the local dynamic grid length scale Δ in step 3 is * Determine as follows:

[0089]

[0090] Among them, π is the ratio of pi, which is 3.14;

[0091] According to the actual simulated turbulent kinetic energy k u , the original modeled turbulent kinetic energy k m and the shielding function F in step 1 GAS , construct the dynamic scale-dependent regulation function D f ,like Figure 2 As shown; the scale ratio is defined as the actual simulated turbulent kinetic energy k u and the original modeled turbulent kinetic energy k m The dynamic scale-dependent regulation function D f is a function related to the scale ratio, which is obtained by the following formula:

[0092]

[0093] l GAS =(1-F GAS ) u +F GAS l m

[0094]

[0095] Among them, l u is the grid-related scale, l m is the turbulence length scale given by the SST k-ω turbulence model and is given by:

[0096]

[0097] Among them, ω m is the specific dissipation rate of the original model obtained by the SST k-ω turbulence model; β * is a constant coefficient in the SST k-ω turbulence model and is taken as 0.09.

[0098] Step 5, reconstruct the turbulent viscosity of the SST k-ω turbulence model;

[0099] In this step, the dynamic scale-dependent adjustment function D in step 4 is used f The turbulent viscosity ν in the k-ω series turbulence models t By adjusting the turbulent viscosity ν, we can obtain the reconstructed turbulent viscosity ν. sfs In this embodiment, the SST k-ω turbulence model in the k-ω series model is taken as an example, and the turbulent viscosity ν in the SST k-ω turbulence model is t By adjusting the turbulent viscosity ν, we can obtain the reconstructed turbulent viscosity ν. sfs , obtained from the following formula:

[0100]

[0101] Where a1 is a constant coefficient in the SST k-ω turbulence model and is taken as 0.31.

[0102] Step 6: Use the reconstructed turbulent viscosity to perform turbulence simulation:

[0103] In this step, the turbulent viscosity ν reconstructed in step 5 is used sfs , calculate the Reynolds stress, and update the transport equation of the k-ω series turbulence model, replacing the turbulent viscosity ν in the original transport equation of the k-ω series turbulence model t , combined with the k-ω series model, the grid adaptive turbulence simulation method based on Vreman dynamic coefficient coupling k-ω series model is obtained. In this embodiment, the SST k-ω turbulence model in the k-ω series model is taken as an example, and the turbulent viscosity ν reconstructed in step 5 is used. sfs , calculate the Reynolds stress, and update the transport equation of the SST k-ω turbulence model, replacing the turbulent viscosity ν in the original transport equation of the SST k-ω turbulence model t The new transport equation obtained is shown as follows:

[0104]

[0105] The new transport equation is combined with the SST k-ω turbulence model to obtain a grid adaptive turbulence simulation method based on Vreman dynamic coefficient coupled k-ω series models, which is used for the numerical simulation of the back-step example.

[0106] The fully implicit coupling solution technology is used for transient calculations, and the time steps all meet the CFL conditions in engineering computational fluid dynamics. In addition, in order to better demonstrate the advantages of the embodiments of the present invention, the SST k-ω turbulence model in the k-ω series model is selected to perform numerical simulations on the back step example, and the numerical simulation results are compared with those of the grid adaptive turbulence simulation method based on the Vreman dynamic coefficient coupled k-ω series model of the present invention.

[0107] Figure 3 Schematic diagram of the turbulent eddy structure of the back-step example calculated using the SST k-ω model in the k-ω series of models, colored by the turbulent viscosity ratio.

[0108] Figure 4 It is a schematic diagram of the turbulent eddy structure of a back-step example calculated by a grid adaptive turbulence simulation method based on the Vreman dynamic coefficient coupled k-ω series model of the present invention, and is colored by the turbulent viscosity ratio.

[0109] Figure 3 , Figure 4 Comparative analysis shows that the turbulent eddy structure analysis capability obtained by calculating the grid adaptive turbulence simulation method based on the Vreman dynamic coefficient coupled k-ω series model proposed in the present invention is stronger than the result calculated by the SST k-ω model. It can capture richer turbulent structures with the same number of grids and provide more accurate flow field details.

[0110] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the protection scope of the present invention.

[0111] In summary, the grid adaptive turbulence simulation method based on Vreman dynamic coefficient coupled k-ω series model proposed in the present invention constructs Vreman dynamic coefficient based on local velocity gradient tensor, identifies local grid size, and determines local dynamic grid length scale Δ * , and then reconstruct the turbulent viscosity by constructing scale-related functions through turbulent energy spectrum integration to achieve grid adaptive simulation, effectively overcoming the problem of high empirical dependence of existing turbulent mixing models on grids. By introducing the Vreman dynamic coefficient, the grid adaptability of the grid adaptive turbulence simulation method is further broadened, which can greatly reduce the computational cost while improving the calculation accuracy, significantly speed up the turbulence simulation process, and provide a new method for flow prediction in complex engineering fields.

Claims

1. A grid adaptive turbulence simulation method based on Vreman dynamic coefficient coupled k-ω series model, characterized in that: The steps include: Step 1, determine whether to apply the shielding function; Step 2, calculate the Vreman dynamic coefficient; Step 3, determine the local dynamic grid length scale; Step 4: construct a scale-dependent adjustment function; Step 5, reconstruct the turbulent viscosity of the k-ω series model; Step 6, use the reconstructed turbulent viscosity to perform turbulence simulation; ① The determination of whether to apply the shielding function includes: Combined with the type of flow state being simulated, determine whether to use the shielding function F GAS Specifically, when the flow state type is free shear flow, the shielding function is not used. In this case, the shielding function F GAS = 0; when the flow state type is near-wall flow, a shielding function is used, and the shielding function F GAS Select one of the F1 shielding function and F2 shielding function from the DDES-SST model and the F from the DDES-SA model. d Masking function; ② The calculation of the Vreman dynamic coefficient includes: Based on the local velocity gradient tensor α ij , construct the Vreman dynamic coefficient C DGAS ; The Vreman dynamic coefficient C DGAS It is obtained from the following formula: β ij ＝α mi α mj Where S is the strain rate, C GAS The empirical coefficient is taken as 0.6; ③ The determination of the local dynamic grid length scale includes: Combined with the shielding function F described in step 1 GAS , and the Vreman dynamic coefficient C described in step 2 DGAS , determine the local dynamic grid length scale Δ * ; The local dynamic grid length scale Δ * Given by: D * =C DGAS [(1-F GAS )D vol +F GAS D max ] D max =max(Δ x ,D y ,D z ) Among them, Δ x is the length of the local hexahedral grid, Δ y is the width of the local hexahedral grid, Δ z is the height of the local hexahedral grid; ④ The scale-related adjustment function constructed based on the turbulence energy spectrum integral coupling k-ω series model includes: According to the modeling method of turbulent kinetic energy in the k-ω series model, the original modeled turbulent kinetic energy k is obtained. m , according to the local dynamic grid length scale Δ * Based on the turbulence energy spectrum, the actual simulated turbulent kinetic energy k is obtained by integration u ; The actual simulated turbulent kinetic energy k u It is obtained from the following formula: Among them, C k is the Kolmogorov constant coefficient, which is taken as 1.5; ε is the actual turbulence dissipation rate; κ c To solve the turbulent cutoff wave number, the local dynamic grid length scale Δ * Determine as follows: Among them, π is the ratio of pi, which is 3.14; According to the actual simulated turbulent kinetic energy k u , the original modeled turbulent kinetic energy k m and the shielding function F described in step 1 GAS , construct the dynamic scale-dependent regulation function D f ; Define the scale ratio as the actual simulated turbulent kinetic energy k u and the original modeled turbulent kinetic energy k m The dynamic scale-dependent regulation function D f is a function related to the ratio of the scales, which is obtained by the following formula: the GAS =(1-F GAS )l u +F GAS the m Among them, l u is the grid-related scale, l m is the turbulence length scale given by the k-ω series model, which is obtained from the following equation: Among them, ω m is the specific dissipation rate of the original model obtained by the k-ω series model; β * is the constant coefficient in the k-ω series model, which is taken as 0.09; ⑤The turbulent viscosity of the reconstructed k-ω series model includes: Use the dynamic scale-dependent adjustment function D described in step 4 f For the turbulent viscosity ν in the k-ω series model t By adjusting the turbulent viscosity ν, we can obtain the reconstructed turbulent viscosity ν. sfs , obtained from the following formula: n sfs =D f ·n t ⑥ The turbulence simulation using the reconstructed turbulent viscosity includes: Using the reconstructed turbulent viscosity ν described in step 5 sfs , calculate the Reynolds stress, and update the transport equation of the k-ω series model, replacing the turbulent viscosity ν in the original transport equation of the k-ω series model t , combined with the k-ω series model, the grid adaptive turbulence simulation method based on Vreman dynamic coefficient coupling k-ω series model is obtained.

Images