Strong shear turbulence dynamic grid adaptive simulation method coupled with k-epsilon series model

Abstract

The invention discloses a strong shear turbulence dynamic grid adaptive simulation method coupled with a k-epsilon series model. The method mainly comprises the following steps of: 1, judging whether a shielding function is applied or not; step 2, calculating a Vman dynamic coefficient; step 3, identifying the self-adaptive length scale of the shear layer; step 4, constructing a scale-related adjustment function based on the turbulence energy spectrum integral; 5, reconstructing the turbulence viscosity of the k-epsilon series model by using an adjustment function; and 6, performing turbulence simulation by using the reconstructed turbulence viscosity. According to the method, the size of a local grid is identified by constructing a Vman dynamic coefficient, the length scale of the local grid is jointly determined in combination with the length scale of the shear layer adaptive sub-grid, and then the scale correlation function is constructed through the turbulence energy spectrum integral to reconstruct the turbulence viscosity, so that grid adaptive simulation is realized; an effective numerical simulation method is provided for solving low-consumption and high-precision simulation of strong shear turbulence in a complex engineering flow problem.

Description

Technical Field

[0001] The present invention relates to the field of complex fluid mechanics calculation of aircraft engines and gas turbines, and in particular to a dynamic grid adaptive simulation method for strong shear turbulence coupled with a k-ε series model. Background Art

[0002] Strong shear flows in turbulent phenomena pose a challenge to engineering turbulence prediction due to their large velocity gradients, multi-scale turbulent interactions, and non-equilibrium transport. The current mainstream Reynolds-averaged NS (RANS) method, while computationally efficient, lacks sufficient accuracy for predicting multi-scale unsteady flows and anisotropic turbulence. Large eddy simulation (LES), while highly accurate, has an excessively high computational cost, far exceeding the requirements of engineering applications. Over the past two decades, the RANS-LES hybrid simulation method has found a balance between accuracy and computational efficiency by employing the RANS method in near-wall regions and the LES method in mainstream regions. This approach resolves the contradiction between simulation accuracy and computational cost for complex flows and offers certain advantages. However, the traditional RANS-LES hybrid model places high demands on the mesh and requires the user to have extensive experience in high-precision numerical simulations, thus limiting its ability to predict complex engineering flows.

[0003] Although the emerging grid adaptive method has improved the ability to simulate complex flows, with the further increase in engineering demand for low-cost and high-precision flow prediction, typical flows such as jets and separation and reattachment have a large number of anisotropic grids under sparse grids, and the flow structure has multi-scale nonlinear characteristics during its development and evolution. Existing technologies require a large computational cost to analyze such complex flow scenarios coupled with strong shear flows, and still face challenges in predicting the strong velocity gradients and complex vortex structures within the sparse grids. It is urgent to develop more intelligent adaptive algorithms to improve the ability to simulate complex turbulence in lower-cost scenarios.

[0004] In order to improve the accuracy of RANS-LES hybrid simulation of strong shear turbulence, the present invention proposes an adaptive sub-grid length scale for the shear layer and simultaneously introduces an improved scheme for the Vreman dynamic coefficient. By introducing adaptive dynamic coefficients and improved grid scale definitions, the ability of sparse anisotropic grids to capture complex flows can be better improved. By coupling improvements to the grid scale and dynamic coefficients, the versatility of the grid adaptation method is further enhanced, and the cost of high-fidelity prediction of strong shear turbulence and complex multi-scale vortex structures is reduced. Compared with the disclosed invention patents CN 115034162 A and CN118133701 A, it can further optimize flow simulation without significantly increasing computational costs, providing higher prediction accuracy and better computational efficiency. Summary of the Invention

[0005] (1) Technical issues to be resolved

[0006] The purpose of this paper is to propose a dynamic mesh adaptive simulation method for strong shear turbulence coupled to the k-ε series model. By incorporating the Vreman dynamic coefficient, the method dynamically adjusts the mesh size under strong shear flow and high Reynolds number conditions to optimize accuracy in strong shear regions, enhance adaptability to complex flow structures, balance high precision with high efficiency, significantly reduce computational overhead, and significantly accelerate the turbulence simulation process. This method aims to provide a low-cost, high-precision, and effective numerical simulation method for solving the complex strong shear flow problems encountered in engineering flows.

[0007] (2) Technical solution

[0008] In order to solve the above technical problems, the present invention provides a dynamic grid adaptive simulation method for strong shear turbulence coupled with a k-ε series model, comprising the following steps:

[0009] Step 1: determine whether to apply the shielding function;

[0010] Step 2: Calculate the Vreman dynamic coefficient;

[0011] Step 3, identifying the shear layer adaptive length scale;

[0012] Step 4: construct a scale-dependent adjustment function based on the turbulence energy spectrum integral;

[0013] Step 5: Use the adjustment function to reconstruct the turbulent viscosity of the k-ε series model;

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

[0015] ① The step 1 of determining whether to apply the shielding function includes:

[0016] 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 Choose to use the F1 shielding function, F2 shielding function from the DDES-SST model and F from the DDES-SA model. d Shield function;

[0017] ② The calculation of the Vreman dynamic coefficient in step 2 includes:

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

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

[0020]

[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 step 3 of identifying the shear layer adaptive length scale includes:

[0026] Use the shielding function F described in step 1 GAS And the Vreman dynamic coefficient C described in step 2 DGAS , determine the local grid length scale Δ*; the local grid length scale Δ* is given by the following formula:

[0027] Δ * =C GAS [(1-F GAS )Δ SLA +F GAS Δ max ],

[0028]

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

[0030] Where m,n=1,2,...,8;

[0031] where i = 1, 2, ..., 8;

[0032] 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; is the distance from the cell center to the i-th vertex, is the unit normal vector of the vorticity vector, is the unit vertex position vector, is the unit center position vector;

[0033] According to the following formula:

[0034]

[0035] in, is the strain rate tensor, is the vorticity vector, ω 2 is the square of the vorticity vector; θ is taken as 0.99, a1=0.15,a2=0.3;ν is the viscosity coefficient of the fluid, v t is the turbulent viscosity given by the k-ε series model, v t∞ The far-field turbulent viscosity given by the k-ε series model;

[0036] <vtm>is the average value of the local VTM function value and the VTM function value of the adjacent grid, F KH Based on the <vtm>A two-stage linear function of the value;

[0037] ④ The scale-related adjustment function constructed based on the turbulence energy spectrum integral in step 4 includes:

[0038] 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 grid length scale Δ* described in step 3, the actual simulated turbulent kinetic energy k is obtained by integration based on the turbulent energy spectrum. u The actual modeled turbulent kinetic energy k u It is obtained from the following formula:

[0039]

[0040] Where E(κ) is the Kolmogorov turbulence energy spectrum, κ is the turbulence wave number, C k is the Kolmogorov constant coefficient, which is 1.5, ε is the actual turbulence dissipation rate, κ c The cutoff wave number for resolvable turbulence is determined by the local grid length scale Δ* described in step 3:

[0041]

[0042] Where π is the ratio of circumference to circle, which is 3.14;

[0043] According to the actual modeled 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-related adjustment function D f ; Define the scale ratio as the actual modeled turbulent kinetic energy k u and the original modeled turbulent kinetic energy k m The dynamic scale-related adjustment function D f is a function related to the ratio of the scales, which is obtained from the following formula:

[0044]

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

[0046]

[0047] Among them, l u is the grid-related scale, l m is the turbulence length scale given by the k-ε series model, which is given by:

[0048]

[0049] Among them, ε m is the dissipation rate of the original model given by the k-ε series model;

[0050] ⑤ The step 5 described in which the adjustment function is used to reconstruct the turbulent viscosity of the k-ε series model includes:

[0051] Adopt the dynamic scale-related adjustment function D described in step 4 f For the turbulent viscosity ν in the k-ε series model t By regulating, we can obtain the reconstructed turbulent viscosity ν sfs , which is obtained from the following formula:

[0052] ν sfs =D f ·ν t

[0053] ⑥ The turbulence simulation using the reconstructed turbulent viscosity described in step 6 includes:

[0054] Using the reconstructed turbulent viscosity ν described in step 5 sfs Calculate the Reynolds stress and update the transport equations of the k-ε series model; use the reconstructed turbulent viscosity ν described in step 5 sfs Substituting the turbulent viscosity ν in the k-ε series of models t , combined with the k-ε series model, the strong shear turbulence dynamic grid adaptive simulation method of the coupled k-ε series model is obtained.

[0055] (3) Beneficial effects

[0056] The present invention provides a method for dynamic grid adaptive simulation of strong shear turbulence of coupled k-ε series models, which has the following beneficial effects: constructing Vreman dynamic coefficients based on the local velocity gradient tensor, determining the local grid length scale by identifying the local grid size and combining it with the shear layer adaptive sub-grid length scale, and then reconstructing the turbulent viscosity by constructing a scale-related function through turbulent energy spectrum integration to achieve grid adaptive simulation, effectively improving the existing RANS-LES hybrid simulation method for the problems of high computational cost and poor accuracy in shear layer simulation, and significantly accelerating the turbulence simulation process.

[0057] Compared with the existing technologies CN 115034162 A and CN 118133701 A, the method of the present invention significantly reduces the influence of the grid on the development of the shear layer by introducing the SLA length scale and Vreman dynamic coefficient, and can dynamically adapt to strong shear turbulent flow, achieving more accurate simulation under high Reynolds number and strong shear flow conditions, further improving the prediction accuracy of the grid adaptive turbulence simulation method, and effectively reducing the computational cost, providing an important means for solving the strong shear multi-scale complex flow problem in aircraft engines.

[0058] This method boasts a concise form and strong portability. It integrates well with various grid-adaptive turbulence simulations and can be easily integrated into existing CFD codes for application and expansion. It has broad academic and engineering application prospects, providing an effective, low-cost, high-precision numerical simulation method for solving complex engineering flow problems involving strong shear. It holds significant application value in the prediction of complex flows involving strong shear, such as those in aircraft engines and gas turbines. BRIEF DESCRIPTION OF THE DRAWINGS

[0059] Figure 1 This is a flow chart of a dynamic grid adaptive simulation method for strong shear turbulence coupled with a k-ε series model of the present invention;

[0060] Figure 2 The invention is a method for dynamic grid adaptive simulation of strong shear turbulence coupled with k-ε series model. The method is applied to the specific embodiment of circular pipe jet flow based on the D of shear layer adaptive sub-grid length scale. f Distribution cloud map;

[0061] Figure 3 The turbulent vortex structure diagram of the circular pipe jet example calculated using the original standard k-ε model;

[0062] Figure 4 The figure shows the turbulent vortex structure of a circular pipe jet example calculated by using a strong shear turbulence dynamic grid adaptive simulation method coupled with a k-ε series model of the present invention. DETAILED DESCRIPTION

[0063] The following, in conjunction with the accompanying drawings and examples, uses the standard k-ε model in the k-ε series of models as an example and a circular pipe jet flow as a calculation example to further describe the specific embodiments of the present invention. The following examples are intended only to illustrate the present invention and are not intended to limit the scope of the present invention.

[0064] The present invention provides a dynamic grid adaptive simulation method for strong shear turbulence coupled with a k-ε series model, comprising the following steps:

[0065] Step 1: determine whether to apply the shielding function;

[0066] In this step, the type of flow state being simulated is used to 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 Choose to use the F1 shielding function, F2 shielding function from the DDES-SST model and F from the DDES-SA model. d Masking function.

[0067] In the circular pipe jet calculation example of the embodiment of the present invention, a high-quality hexahedral grid is used to spatially discretize the computational domain, and the total number of grids is approximately 6.05 million.

[0068] Using the F from the DDES-SA model d function, then:

[0069] F GAS =F d

[0070]

[0071] Among them, U ij is the velocity gradient tensor, ν t is the turbulent viscosity given by the standard k-ε model, ν is the viscosity coefficient of the fluid, d is the distance from the local grid to the wall, κ is the Karman constant, which is 0.41; C d1 、C d2 are empirical coefficients, which are 20 and 3 respectively.

[0072] Step 2: Calculate the Vreman dynamic coefficient

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

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

[0075]

[0076] β ij =α mi α mj

[0077]

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

[0079] Step 3, identifying the shear layer adaptive length scale;

[0080] In this step, the shielding function F in step 1 is used GAS And the Vreman dynamic coefficient C in step 2 DGAS , determine the local grid length scale Δ*; the local grid length scale Δ* is given by the following formula:

[0081] Δ * =C GAS [(1-F GAS )Δ SLA +F GAS Δ max ],

[0082]

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

[0084] Where m,n=1,2,...,8;

[0085] where i = 1, 2, ..., 8;

[0086] 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; is the distance from the cell center to the i-th vertex, is the unit normal vector of the vorticity vector, is the unit vertex position vector, is the unit center position vector;

[0087] According to the following formula:

[0088]

[0089] in, is the strain rate tensor, is the vorticity vector, ω 2 is the square of the vorticity vector; θ is taken as 0.99, a1=0.15,a2=0.3;ν is the viscosity coefficient of the fluid, v t is the turbulent viscosity given by the k-ε series model, and in this embodiment is the turbulent viscosity given by the standard k-ε model; v t∞ is the far-field turbulent viscosity given by the k-ε series model, and in this embodiment is the far-field turbulent viscosity given by the standard k-ε model.

[0090] The above VTM function ranges from 0 to 1. In the quasi-two-dimensional region at the initial stage of separation, the vorticity vector Also the strain rate tensor The characteristic vector of VTM function tends to 0. In fully developed three-dimensional turbulence, the vorticity vector Direction and strain rate tensors The eigenvector correlation is weak, and the VTM function value tends to 1.

[0091] <vtm>is the average value of the local VTM function value and the VTM function value of the adjacent grid, F KH Based on the <vtm>A two-stage linear function of the value. <vtm>When the value is less than a certain threshold, the VTM function value is small. <vtm>As the value gradually increases, the VTM function value rapidly increases to 1, thereby reducing the grid size in the early stage of separation. Through the sub-grid stress model, it can be seen that when the VTM function value increases to 1, the turbulent viscosity coefficient also obtains the expected reduction.

[0092] Step 4: construct a scale-dependent adjustment function based on the turbulence energy spectrum integral;

[0093] In this step, according to the modeling method of turbulent kinetic energy in the k-ε series model, the original modeled turbulent kinetic energy k is obtained. m In this embodiment, according to the modeling method of turbulent kinetic energy in the standard k-ε model, the original modeled turbulent kinetic energy k is obtained. m .

[0094] According to the local grid length scale Δ* in step 3, the actual modeled turbulent kinetic energy k is obtained by integration based on the turbulent energy spectrum. u The actual modeled turbulent kinetic energy k u It is obtained from the following formula:

[0095]

[0096] Where E(κ) is the Kolmogorov turbulence energy spectrum, κ is the turbulence wave number, C k is the Kolmogorov constant coefficient, which is 1.5, ε is the actual turbulence dissipation rate, κ c The cutoff wave number for resolvable turbulence is determined by the local grid length scale Δ* described in step 3:

[0097]

[0098] Where π is the ratio of circumference to circle, which is 3.14;

[0099] According to the actual 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-related adjustment function D f ,like Figure 2 shown.

[0100] The scale ratio is defined as the actual modeled turbulent kinetic energy k u and the original modeled turbulent kinetic energy k m The dynamic scale-related adjustment function D f is a function related to the scale ratio, which is obtained from the following formula:

[0101]

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

[0103]

[0104] Among them, l u is the grid-related scale, l m is the turbulence length scale given by the k-ε series model, which is given by:

[0105]

[0106] Among them, ε m is the dissipation rate of the original model given by the k-ε series model.

[0107] Step 5: Use the adjustment function to reconstruct the turbulent viscosity of the k-ε series model;

[0108] In this step, the dynamic scale-related adjustment function D described in step 4 is used. f For the turbulent viscosity ν in the k-ε series model t By regulating, we can obtain the reconstructed turbulent viscosity ν sfs , which is obtained from the following formula:

[0109] ν sfs =D f ·ν t

[0110] In this embodiment, the standard k-ε turbulence model in the k-ε series of models is taken as an example, and the turbulent viscosity ν in the standard k-ε turbulence model is calculated. t By regulating, we can obtain the reconstructed turbulent viscosity ν sfs As follows:

[0111]

[0112] Among them, C μ is the constant coefficient in the standard k-ε model and is taken as 0.09.

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

[0114] 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 model. In this embodiment, the standard k-ε model in the k-ε series model is taken as an example, and the turbulent viscosity ν reconstructed in step 5 is used. sfs Replacement of the turbulent viscosity ν in the standard k-ε model t , the new transport equation is shown as follows:

[0115]

[0116] Among them, the model coefficient C 1ε is 1.44, C 2ε It is 1.92.

[0117] The obtained new transport equation is combined with the standard k-ε model to obtain a dynamic grid adaptive simulation method for strong shear turbulence of the coupled k-ε series model, and is used for the numerical simulation of a circular pipe jet example.

[0118] Transient calculations were performed using a fully implicit coupled solution technique, with a time step that met the CFL condition in engineering computational fluid dynamics. A standard k-ε model from the k-ε series of models was also selected to perform numerical simulations of a circular pipe jet flow example from an embodiment of the present invention. The results were compared with the numerical simulation results of the dynamic mesh adaptive simulation method for strong shear turbulence coupled with the k-ε series of models of the present invention.

[0119] Figure 3 This is a diagram of the turbulent vortex structure of a circular pipe jet flow example calculated using the original standard k-ε model in the k-ε series of models, colored by the turbulent viscosity ratio.

[0120] Figure 4 This is a turbulent vortex structure diagram of a circular pipe jet example calculated using a strong shear turbulence dynamic grid adaptive simulation method coupled with a k-ε series model of the present invention, and is colored using the turbulent viscosity ratio.

[0121] Figure 3 、 Figure 4 Comparative analysis shows that the turbulent vortex structure analysis capability calculated by the dynamic grid adaptive simulation method of strong shear turbulence coupled with the k-ε series model proposed in the present invention is stronger than the vortex structure calculated by the standard k-ε model. With the same number of grids, it can capture richer turbulent structures and provide more accurate flow field details.

[0122] 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 scope of protection of the present invention.

[0123] In summary, the proposed method for dynamic grid-adaptive simulation of strong shear turbulence coupled to the k-ε series model introduces the Vreman dynamic coefficient to identify the local shear layer adaptive length scale, improving the identification of the local grid length scale Δ*. Furthermore, by integrating the turbulent energy spectrum to construct a scale-dependent function and reconstruct the turbulent viscosity, more efficient grid-adaptive simulation is achieved. This significantly reduces computational cost and significantly accelerates the turbulence simulation process. Furthermore, by dynamically adjusting the grid scale, it is more suitable for the numerical simulation of strong shear turbulence, providing a high-precision, low-cost numerical simulation method for solving high-precision simulations of strong shear turbulence in complex engineering flow problems.< / vtm> < / vtm> < / vtm> < / vtm> < / vtm> < / vtm>

Claims

1. A dynamic grid adaptive simulation method for strong shear turbulence coupled with k-ε series models, characterized by: The steps include: Step 1: determine whether to apply the shielding function; Step 2: Calculate the Vreman dynamic coefficient; Step 3, identifying the shear layer adaptive length scale; Step 4: construct a scale-dependent adjustment function based on the turbulence energy spectrum integral; Step 5: Use the adjustment function to reconstruct the turbulent viscosity of the k-ε series model; Step 6: Use the reconstructed turbulent viscosity to perform turbulence simulation; ① The step 1 of determining 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 Choose to use the F1 shielding function, F2 shielding function from the DDES-SST model and F from the DDES-SA model. d Shield function; ② The calculation of the Vreman dynamic coefficient in step 2 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 step 3 of identifying the shear layer adaptive length scale includes: Use the shielding function F described in step 1 GAS And the Vreman dynamic coefficient C described in step 2 DGAS , determine the local grid length scale Δ*; the local grid length scale Δ* is given by the following formula: D * =C GAS [(1-F GAS )D SLA +F GAS D max ], D max =max(Δ x ,D y ,D z ), Where m,n=1,2,...,8; where i = 1, 2, ..., 8; 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; is the distance from the cell center to the i-th vertex, is the unit normal vector of the vorticity vector, is the unit vertex position vector, is the unit center position vector; According to the following formula: in, is the strain rate tensor, is the vorticity vector, ω 2 is the square of the vorticity vector; θ is taken as 0.99, a1=0.15,a2=0.3;ν is the viscosity coefficient of the fluid, v t is the turbulent viscosity given by the k-ε series model, v t∞ The far-field turbulent viscosity given by the k-ε series model; <VTM> is the average value of the local VTM function value and the VTM function value of the adjacent grid, F KH is a two-stage linear function based on the <VTM> value; ④ The scale-related adjustment function constructed based on the turbulence energy spectrum integral in step 4 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 grid length scale Δ* described in step 3, the actual simulated turbulent kinetic energy k is obtained by integration based on the turbulent energy spectrum. u The actual modeled turbulent kinetic energy k u It is obtained from the following formula: Where E(κ) is the Kolmogorov turbulence energy spectrum, κ is the turbulence wave number, C k is the Kolmogorov constant coefficient, which is 1.5, ε is the actual turbulence dissipation rate, κ c The cutoff wave number for resolvable turbulence is determined by the local grid length scale Δ* described in step 3: Where π is the ratio of circumference to circle, which is 3.14; According to the actual modeled 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-related adjustment function D f ; Define the scale ratio as the actual modeled turbulent kinetic energy k u and the original modeled turbulent kinetic energy k m The dynamic scale-related adjustment function D f is a function related to the ratio of the scales, which is obtained from 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 given by: Among them, ε m is the dissipation rate of the original model given by the k-ε series model; ⑤ The step 5 described in which the adjustment function is used to reconstruct the turbulent viscosity of the k-ε series model includes: Adopt the dynamic scale-related adjustment function D described in step 4 f For the turbulent viscosity ν in the k-ε series model t By regulating, we can obtain the reconstructed turbulent viscosity ν sfs , which is obtained from the following formula: n sfs =D f ·n t ⑥ The turbulence simulation using the reconstructed turbulent viscosity described in step 6 includes: Using the reconstructed turbulent viscosity ν described in step 5 sfs Calculate the Reynolds stress and update the transport equations of the k-ε series model; use the reconstructed turbulent viscosity ν described in step 5 sfs Substituting the turbulent viscosity ν in the k-ε series of models t , combined with the k-ε series model, the strong shear turbulence dynamic grid adaptive simulation method of the coupled k-ε series model is obtained.

Images