Turbulence simulation verification method and device, computer equipment and medium

Abstract

The invention relates to the technical field of simulation, and discloses a turbulence simulation verification method and device, computer equipment and a medium, which are used for improving turbulence simulation precision. The method comprises the steps of obtaining turbulence characteristic parameters of a target object in an actual environment, wherein the turbulence characteristic parameters comprise a turbulence speed and a yaw angle; taking the turbulence characteristic parameter as a boundary condition of a simulation calculation domain entrance, and substituting the turbulence characteristic parameter into a numerical simulation geometric model of the target object for simulation analysis to obtain a turbulence simulation numerical value; comparing the turbulence simulation numerical value with an actual turbulence numerical value corresponding to the numerical value type; and when the difference value between the turbulence simulation numerical value and the actual turbulence numerical value is smaller than a preset value, it is confirmed that the simulation precision of the numerical simulation geometric model of the target object meets the requirement.

Description

Technical Field

[0001] The present invention relates to the field of simulation technology, and in particular to a turbulence simulation verification method, device, computer equipment and medium. Background Art

[0002] CFD (Computational Fluid Dynamics) simulation is one of the important means of aerodynamic research. It mainly evaluates the aerodynamic performance of the target object through numerical simulation methods, especially in the early stage of automobile development when there is no real car, so as to optimize the automobile's exterior styling, improve aerodynamic performance, and achieve the purpose of reducing drag and noise, reducing fuel consumption or increasing cruising range.

[0003] The algorithms in the existing simulation software can simulate some turbulent wind characteristics measured on actual roads or in wind tunnels, such as wind speed pulsation, turbulence intensity, and turbulence integral scale. However, in traditional simulation schemes, due to the incompleteness of the algorithms themselves and some inherent defects of the algorithms, the pulsating wind speed spectrum will drop sharply, which is significantly different from the turbulent wind characteristics measured on actual roads or in wind tunnels, and the accuracy of the simulation results is relatively low. Summary of the invention

[0004] The present invention provides a turbulence simulation verification method, device, computer equipment and medium to solve the technical problem that the simulation results of turbulent wind characteristics are of relatively low accuracy in traditional solutions.

[0005] In a first aspect, a turbulence simulation verification method is provided, comprising:

[0006] Acquire turbulence characteristic parameters of the target object in the actual environment, wherein the turbulence characteristic parameters include turbulence velocity and yaw angle;

[0007] Substituting the turbulence characteristic parameters as boundary conditions at the entrance of the simulation calculation domain into the numerical simulation geometric model of the target object for simulation analysis to obtain turbulence simulation values;

[0008] Comparing the turbulence simulation value with the actual turbulence value corresponding to the value type;

[0009] When the difference between the turbulence simulation value and the actual turbulence value is less than a preset value, it is confirmed that the simulation accuracy of the numerical simulation geometric model of the target object meets the requirements.

[0010] Furthermore, the yaw angle is obtained by:

[0011] Obtaining the average speed and pulsating wind speed of the target object in the actual environment;

[0012] Acquire the pulsating wind speed in the x direction and the pulsating wind speed in the y direction according to the pulsating wind speed;

[0013] The yaw angle is obtained according to an average speed of the target object in an actual environment, the pulsating wind speed in the x direction, and the pulsating wind speed in the y direction.

[0014] Further, the obtaining of the yaw angle according to the average speed of the target object in the actual environment, the fluctuating wind speed in the x direction, and the fluctuating wind speed in the y direction includes:

[0015] The yaw angle is calculated by the following method:

[0016]

[0017] Wherein, β represents the yaw angle, u′ y represents the fluctuating wind speed in the y direction, u′ x represents the fluctuating wind speed in the x direction, Represents the average speed of the target object in the actual environment.

[0018] Furthermore, the turbulent velocity is obtained by:

[0019] Obtaining the average speed and pulsating wind speed of the target object in the actual environment;

[0020] The average velocity and the fluctuating wind speed are added together to obtain the turbulent velocity.

[0021] Furthermore, the pulsating wind speed is obtained by:

[0022] Acquire a first wind speed spectrum synthesis parameter and a second wind speed spectrum synthesis parameter of the target object in the actual environment, wherein the first wind speed spectrum synthesis parameter and the second wind speed spectrum synthesis parameter are respectively determined based on a simulated wind speed spectrum of the target object in the actual environment and a corresponding weight coefficient;

[0023] The pulsating wind speed is calculated according to the first wind speed spectrum synthesis parameter and the second wind speed spectrum synthesis parameter.

[0024] Further, the calculating the pulsating wind speed according to the first wind speed spectrum synthesis parameter and the second wind speed spectrum synthesis parameter includes:

[0025] Preset a first random number and a second random number;

[0026] Calculate a first speed parameter according to the first random number and a first wind speed spectrum synthesis parameter, and calculate a second speed parameter according to the second random number and a second wind speed spectrum synthesis parameter;

[0027] The pulsating wind speed is calculated according to the first speed parameter and the second speed parameter.

[0028] Further, the calculating the pulsating wind speed according to the first speed parameter and the second speed parameter includes:

[0029] The pulsating wind speed is calculated as follows:

[0030]

[0031]

[0032]

[0033] in, represents the pulsating wind speed, x represents the coordinates of the simulation grid loading point, L represents the turbulence integral length, t represents time, k m represents the number of segments, U1 represents the first random number, U2 represents the second random number, N represents the number of samples between the first random number and the second random number, P m,n represents the first wind speed spectrum synthesis parameter, q m,n represents the first wind speed spectrum synthesis parameter, ω m,n represents the angular frequency of the nth wave number in the mth segment of the simulated wind speed spectrum, represents the random phase that changes with time t and coordinates;

[0034] represents the first speed parameter, Represents the second speed parameter.

[0035] Furthermore, the first wind speed spectrum synthesis parameter and the second wind speed spectrum synthesis parameter are determined by:

[0036]

[0037] Among them, P m,n represents the first wind speed spectrum synthesis parameter, q m,n represents the first wind speed spectrum synthesis parameter, a represents the weight coefficient, E(k m ) represents the von Karman wind speed spectrum, and N represents the number of samples between the first random number and the second random number.

[0038] Furthermore, the turbulence characteristic parameters also include turbulence intensity and turbulence integral size.

[0039] In a second aspect, a turbulence simulation verification device is provided, comprising:

[0040] An acquisition module, used to acquire turbulence characteristic parameters of the target object in the actual environment, wherein the turbulence characteristic parameters include turbulence velocity and yaw angle;

[0041] A simulation module, used to substitute the turbulence characteristic parameters as boundary conditions of the entrance of the simulation calculation domain into the numerical simulation geometric model of the target object for simulation analysis to obtain turbulence simulation values;

[0042] A comparison module, used for comparing the turbulence simulation value with an actual turbulence value corresponding to the value type;

[0043] The determination module is used to confirm that the simulation accuracy of the numerical simulation geometric model of the target object meets the requirements when the difference between the turbulence simulation value and the actual turbulence value is less than a preset value.

[0044] In a third aspect, a computer device is provided, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the steps of the turbulence simulation verification method as described in any one of the preceding items when executing the computer program.

[0045] In a fourth aspect, a computer-readable storage medium is provided, wherein the computer-readable storage medium stores a computer program, wherein when the computer program is executed by a processor, the steps of turbulence simulation verification as described in any one of the above items are implemented.

[0046] One of the solutions provided above provides a turbulence simulation verification method, in which the influence of the yaw angle on the consistency of turbulence characteristics such as fluctuating wind speed is taken into account, and it is used as a turbulence characteristic parameter to substitute into the numerical simulation geometric model for simulation verification. Considering the yaw angle in the actual environment can better imitate the actual turbulent wind characteristics, thereby reducing turbulent wind differences, solving the problem of inconsistency between the pulsating wind speed spectrum of the simulation algorithm and the pulsating wind speed spectrum of the actual environment, and improving the simulation precision and accuracy. BRIEF DESCRIPTION OF THE DRAWINGS

[0047] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the accompanying drawings required for use in the description of the embodiments of the present invention will be briefly introduced below. Obviously, the accompanying drawings in the following description are only some embodiments of the present invention. For ordinary technicians in this field, other accompanying drawings can be obtained based on these accompanying drawings without paying creative labor.

[0048] Figure 1 It is a schematic diagram of a process of a turbulence simulation verification method in one embodiment of the present invention;

[0049] Figure 2 is a schematic diagram of a coordinate system of a turbulence simulation verification method in one embodiment of the present invention;

[0050] Figure 3It is a numerical diagram of actual fluctuating wind speed obtained in road tests and wind tunnel tests;

[0051] Figure 4 It is a numerical diagram of a pulsating wind speed simulation obtained by a turbulence simulation verification method in one embodiment of the present invention;

[0052] Figure 5 It is a schematic diagram of the actual yaw angle obtained in road tests and wind tunnel tests;

[0053] Figure 6 It is a numerical diagram of a yaw angle simulation obtained by a turbulence simulation verification method in one embodiment of the present invention;

[0054] Figure 7 It is a schematic diagram of an actual pulsating wind speed spectrum obtained under road tests and wind tunnel tests;

[0055] Figure 8 is a schematic diagram of a pulsating wind speed simulation spectrum obtained by a turbulence simulation verification method in one embodiment of the present invention;

[0056] Fig. 9 It is a structural schematic diagram of a turbulence simulation verification device in one embodiment of the present invention;

[0057] Fig.10 It is a structural schematic diagram of a computer device in one embodiment of the present invention. DETAILED DESCRIPTION

[0058] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0059] It is understood that the target object on the actual road will usually have turbulent wind due to the presence of other objects such as surrounding trees and buildings. Therefore, performance evaluation under such turbulent wind environment conditions is essential. The target object can be a mobile object such as a car or a robot, or an object that moves relative to other reference objects, such as a fixed building. Taking the target object as a car as an example, when the car is driving on the actual road, the driving environment is relatively complicated due to the presence of other objects such as surrounding trees, buildings and vehicles. Therefore, if the car wants to show good aerodynamic performance when driving on the actual road, performance evaluation under the turbulent wind environment conditions of the actual road is also essential.

[0060] In view of the above-mentioned turbulence phenomenon, it is necessary to evaluate the aerodynamic performance of the vehicle under turbulent environment conditions. Based on this, the embodiment of the present invention mainly provides a turbulence simulation verification method, which is applied to the turbulence verification analysis of various target objects. The target object can be a mobile object such as a car or a robot, or an object that moves relative to other reference objects, such as a fixed building, without specific limitation, so that the aerodynamic performance of the target object can be evaluated. The following is a description through a specific embodiment.

[0061] In one embodiment, if Figure 1 As shown, a turbulence simulation verification method is provided, comprising the following steps:

[0062] S10: Acquire turbulence characteristic parameters of the target object in the actual environment, wherein the turbulence characteristic parameters include turbulence velocity and yaw angle.

[0063] In this embodiment, it is necessary to first obtain the turbulence condition data of the target object in the actual environment, and extract the turbulence characteristic parameters from these turbulence condition data. The turbulence characteristic parameters characterize the turbulent wind characteristics when the target object moves or moves relative to each other in the actual environment. In this embodiment, the turbulence characteristic parameters include turbulence velocity and yaw angle. Preferably, the turbulence characteristic parameters also include turbulence intensity and turbulence integral size. Among them, the turbulence integral size is also called turbulence integral scale, which is collectively referred to as turbulence integral size here. The turbulence velocity characterizes the speed under this condition, which can be understood as the speed or wind speed of the target object, or it can be the speed determined separately under this turbulence condition, which is not limited here. The yaw angle is a key parameter during the movement or relative movement of the target object, which is used to describe the motion state of the target object. It should be understood that during the movement or relative movement of the target object, due to the influence of turbulence, etc., the angle of the target object body will change. Therefore, the yaw angle can be used to characterize the balance of the target object.

[0064] Exemplarily, the target object may be a mobile object such as a car or a robot, or an object that moves relative to other reference objects, such as a fixed building. For ease of explanation, the target object is taken as an example of a car. Exemplarily, the turbulence characteristic parameters of the car measured on an actual road or in a wind tunnel may be collected. Specifically, the turbulence characteristic parameters may be collected under driving environments such as urban conditions, national highway conditions, or high-speed conditions, without specific limitation. It is also worth noting that in actual applications, the road turbulent wind environment is more complex due to the influence of surrounding buildings, trees, vehicles, etc., and the collection and extraction of data must have certain repeatability and representativeness. Therefore, in this process, in order to extract valid data, the collected data may be classified, analyzed, and extracted to extract valid turbulence characteristic parameters.

[0065] S20: Substituting the turbulence characteristic parameters as boundary conditions at the entrance of the simulation calculation domain into the numerical simulation geometric model of the target object for simulation analysis to obtain turbulence simulation values.

[0066] In this embodiment, the numerical simulation geometric model of the target object can be pre-constructed according to the type of the target object to be analyzed and verified. Exemplarily, taking the target object as a car, a numerical simulation geometric model corresponding to the real car will be pre-constructed, and then the turbulence characteristic parameters of the car in the actual environment will be used as the boundary conditions of the simulation calculation domain entrance of the numerical simulation geometric model, and substituted into the numerical simulation geometric model of the real car for simulation analysis, thereby obtaining turbulence simulation values. Taking the target object as a robot, a numerical simulation geometric model corresponding to the physical robot will be pre-constructed, and then the turbulence characteristic parameters of the physical robot in the actual environment will be used as the boundary conditions of the simulation calculation domain entrance of the numerical simulation geometric model, and substituted into the numerical simulation geometric model of the physical robot for simulation analysis, thereby obtaining turbulence simulation values.

[0067] It should be noted that the numerical simulation geometric model corresponding to the target object can be constructed based on the target object and simulation software, among which CFD (Computational Fluid Dynamics) simulation is one of the important means of aerodynamic research. It mainly evaluates the aerodynamic performance of the vehicle through numerical simulation methods, especially in the early stage of automobile development when there is no real car stage, so as to optimize the automobile exterior styling scheme, improve aerodynamic performance, and achieve the purpose of reducing drag and noise, reducing fuel consumption or increasing cruising range. Therefore, the numerical simulation geometric model of the target object can be constructed based on CFD simulation. The construction process includes meshing, parameter selection and solution setting, which are not limited to the specifics and will not be explained in detail.

[0068] It should be noted that, as an example, during the construction process, the following conditions must be met: the numerical simulation geometry model must be as consistent as possible with the target object, including retaining the exterior details as much as possible; the simulation domain size must meet the blocking ratio requirements, and the key components of the target object must be meshed during meshing. For example, taking a car as an example, the key components may include the grille, A-pillars, rearview mirrors, tires, and rear spoilers of the car. It should be noted that the above conditions are met only to further improve the accuracy of the numerical simulation geometry model, and there is no specific limitation.

[0069] Of course, the numerical simulation geometric model of a car or other target object can also be constructed by other simulation methods, and the construction process is not described in detail or limited here.

[0070] S30: Compare the turbulence simulation value with the actual turbulence value corresponding to the value type.

[0071] S40: When the difference between the turbulence simulation value and the actual turbulence value is less than a preset value, it is confirmed that the simulation accuracy of the numerical simulation geometric model of the target object meets the requirements.

[0072] In this embodiment, by taking the turbulence characteristic parameters as the boundary conditions of the entrance of the simulation calculation domain and substituting them into the numerical simulation geometric model of the target object for simulation analysis, the turbulence simulation value can be obtained. This is equivalent to using the typical data of the actual environment, that is, the turbulence characteristic parameters actually obtained by the real vehicle as the boundary conditions of the model, and simulating the numerical simulation geometric model of the target object to imitate the actual turbulence conditions, thereby obtaining the turbulence simulation value. Finally, the turbulence simulation value is compared with the actual turbulence value corresponding to the numerical type. When the difference between the turbulence simulation value and the actual turbulence value is less than the preset value, it means that the simulation accuracy of the numerical simulation geometric model of the target object meets the requirements.

[0073] Exemplarily, as some examples, the turbulence simulation values ​​may include fluctuating wind speed simulation values, yaw angle simulation values, fluctuating wind speed simulation spectrum, and may also include turbulence intensity simulation values ​​and turbulence integral scale simulation values. In the verification process, for example, the fluctuating wind speed simulation values, and / or the fluctuating yaw angle simulation values, and the corresponding fluctuating wind speed simulation spectrum are compared with the actual fluctuating wind speed, the actual yaw angle simulation values, and the actual fluctuating wind speed simulation spectrum, respectively, so as to confirm whether the simulation accuracy of the numerical simulation geometric model of the target object meets the requirements.

[0074] It can be seen that when the simulation accuracy of the numerical simulation geometric model of the target object is high, the numerical simulation geometric model can be used to evaluate the dynamic performance of the target object. Taking a car as an example, in the early stage of model development, without the condition of a real car, the corresponding numerical simulation geometric model of the car can be used to simulate the real driving wind environment of the car, thereby improving the effectiveness of the aerodynamic scheme of drag reduction and noise reduction, improving design work efficiency and vehicle comfort, reducing the fuel consumption of the vehicle in actual use or increasing the cruising range, and improving the practicality of the real car.

[0075] It can be seen that in this embodiment, a turbulence simulation verification method is provided. In this turbulence simulation verification method, the influence of the yaw angle on the consistency of turbulence characteristics such as the fluctuating wind speed is taken into account, and the yaw angle is used as a turbulence characteristic parameter to be substituted into the numerical simulation geometric model for simulation verification. Considering the yaw angle in the actual environment can better imitate the actual turbulent wind characteristics, thereby reducing the turbulent wind difference, solving the problem of the inconsistency between the fluctuating wind speed spectrum of the simulation algorithm and the fluctuating wind speed spectrum of the actual environment, and improving the simulation precision and accuracy.

[0076] In one embodiment, the turbulent velocity u i Obtained through the following methods:

[0077] Get the average speed of the target object in the actual environment and fluctuating wind speed u′ i ;

[0078] The average wind speed and fluctuating wind speed u′ i The sum is used as the turbulent velocity.

[0079] It can be expressed as follows:

[0080]

[0081] Among them, u i represents the turbulent velocity in the x, y, and z directions, and the corresponding u′ i Indicates the pulsating wind speed in the x, y, and z directions. It represents the average speed in the x, y, and z directions. When calculating, add the average speed and the fluctuating wind speed in the same direction to obtain the turbulent speed in the corresponding direction.

[0082] like Figure 2 As shown, taking a car as an example, the position relationship of the three directions is as follows Figure 2 As shown, it can be seen that the center of gravity of the car is taken as the origin, the length direction of the vehicle is taken as the x direction, the length direction is taken as the y direction, and the height direction is taken as the z direction. Of course, for other target objects, the above directions are determined by analogy, and the specific details are not explained in detail here.

[0083] In this embodiment, in order to facilitate simulation calculation and improve simulation accuracy, the average speed of the target object in the actual environment is and fluctuating wind speed u′ i After conversion, the average speed It refers to the moving speed of the target object in the x direction. The speed parameter after adding the two is used as the boundary condition to set the parameters of the simulation calculation domain entrance, which can also reduce the setting of parameters during simulation. It should be noted that in actual applications, it is also possible to obtain the average wind speed and pulsating wind speed u′ of the target object in the actual environment. i ; The average wind speed and the pulsating wind speed u′ i The turbulent velocity is added together, and is not specifically limited. It should be understood that in the actual movement process, there is only a slight difference between the average wind speed and the movement speed of the target object, but the converted turbulent velocity reflects the velocity characteristics of turbulence very well.

[0084] Among them, in calculating the fluctuating wind speed u′ i , the embodiment of the present application provides two calculation methods, which are described below respectively.

[0085] First method

[0086] In this method, the Spectral Synthesizer simulation method is used. The pulsating velocity component in this method is calculated by Fourier harmonic synthesis of the non-divergent velocity vector field. The pulsating velocity u′ i The calculation formula is as follows:

[0087]

[0088]

[0089]

[0090] in, represents the coordinates of the simulation grid loading point, t represents time, l,τ represent the turbulence length scale and turbulence time scale respectively, ε ijk Refers to the permutation tensor during tensor operations. and ω n represent n wave number vectors and frequencies respectively.

[0091] The second calculation method

[0092] It should be noted that in order to further improve the simulation accuracy of the numerical simulation geometric model and simplify the calculation process, the embodiment of the present application also calculates the pulsating wind speed u′ i Specifically, the pulsating wind speed u′ can also be obtained by the following method: i :

[0093] S101: Acquire a first wind speed spectrum synthesis parameter and a second wind speed spectrum synthesis parameter of the target object in the actual environment, wherein the first wind speed spectrum synthesis parameter and the second wind speed spectrum synthesis parameter are respectively determined based on a simulated wind speed spectrum of the target object in the actual environment and a corresponding weight coefficient;

[0094] S102: Calculate the pulsating wind speed according to the first wind speed spectrum synthesis parameter and the second wind speed spectrum synthesis parameter.

[0095] In this embodiment, a weight coefficient for calculating the wind speed spectrum is set, and then based on the weight coefficient and the simulated wind speed spectrum corresponding to the target object in the actual environment, two wind speed spectrum synthesis parameters are first converted, namely the first wind speed spectrum synthesis parameter and the second wind speed spectrum synthesis parameter. Exemplarily, the simulated wind speed spectrum can be a von Karman wind speed spectrum or other similar wind speed spectrum expression forms, which are not specifically limited.

[0096] Taking the von Karman wind speed spectrum as an example, the von Karman wind speed spectrum is calculated as follows:

[0097]

[0098] Among them, E(k m ) represents the von Karman wind speed spectrum, I represents the turbulence intensity, which is a constant value in the actual environment, L represents the turbulence integral length, which is a constant value in the actual environment, u′ i represents the pulsating wind speed, and f represents the sampling frequency of the speed, which can be the wind speed or the moving speed of the target object, and is not specifically limited.

[0099] In this embodiment, based on the von Karman wind speed spectrum and the weight coefficient, the pulsating wind speed spectrum is decomposed into a first wind speed spectrum synthesis parameter and a second wind speed spectrum synthesis parameter. Finally, the pulsating wind speed is calculated according to the first wind speed spectrum synthesis parameter and the second wind speed spectrum synthesis parameter, which can improve the consistency of the pulsating wind speed spectrum.

[0100] Further, in step S102, that is, calculating the fluctuating wind speed according to the first wind speed spectrum synthesis parameter and the second wind speed spectrum synthesis parameter, includes:

[0101] Preset a first random number and a second random number;

[0102] Calculate a first speed parameter according to the first random number and a first wind speed spectrum synthesis parameter, and calculate a second speed parameter according to the second random number and a second wind speed spectrum synthesis parameter;

[0103] The pulsating wind speed is calculated according to the first speed parameter and the second speed parameter.

[0104] Specifically, calculating the pulsating wind speed according to the first speed parameter and the second speed parameter includes:

[0105] The pulsating wind speed is calculated as follows:

[0106]

[0107]

[0108]

[0109] in, represents the pulsating wind speed, x represents the coordinates of the simulation grid loading point, L represents the turbulence integral length, t represents time, k m represents the number of segments, U1 represents the first random number, U2 represents the second random number, P m,n represents the first wind speed spectrum synthesis parameter, q m,n represents the first wind speed spectrum synthesis parameter, ω m,n represents the angular frequency of the nth wave number in the mth segment of the simulated wind speed spectrum, Indicates the random phase that changes with time t and coordinates. For example, U1 and U2 respectively represent random numbers from -1 to 1, and their values ​​conform to the normal distribution. Random cosine function, Random sine function.

[0110] represents the first speed parameter, Represents the second speed parameter.

[0111] In this embodiment, two random numbers are set and two speed parameters are calculated respectively in combination with the corresponding first wind speed spectrum synthesis parameter and the second wind speed spectrum synthesis parameter, which are used to simplify the calculation of the two-dimensional speed parameters. The parameters considered are also different. After comparison, the calculated pulsating wind speed spectrum simulation results can be more in line with the actual situation and have higher accuracy.

[0112] Furthermore, the first wind speed spectrum synthesis parameter and the second wind speed spectrum synthesis parameter are determined by:

[0113]

[0114] Among them, P m,n represents the first wind speed spectrum synthesis parameter, q m,n represents the first wind speed spectrum synthesis parameter, a represents the weight coefficient, E(k m ) represents the von Karman wind speed spectrum, N represents the number of samples between the first random number and the second random number, and illustratively, N is generally taken as 100.

[0115] In this embodiment, it is described how to perform calculations based on the von Karman wind speed spectrum and the corresponding weight coefficient a, and the calculation method of the first wind speed spectrum synthesis parameter and the second wind speed spectrum synthesis parameter is clarified to ensure the feasibility of the solution. It is also worth noting that, on the one hand, in practical applications, the value of the weight coefficient can be configured, and on the other hand, the weight coefficient is set in this embodiment, and after comparison and verification, it can improve the consistency between the calculated pulsating wind speed and the actual one. It should also be noted that if other simulated wind speed spectra are used, the calculation process can be based on this corresponding experimental transformation, and there is no specific limitation.

[0116] In one embodiment, the yaw angle is obtained by:

[0117] Obtaining the average speed and pulsating wind speed of the target object in the actual environment;

[0118] Acquire the pulsating wind speed in the x direction and the pulsating wind speed in the y direction according to the pulsating wind speed;

[0119] The yaw angle is obtained according to an average speed of the target object in an actual environment, the pulsating wind speed in the x direction, and the pulsating wind speed in the y direction.

[0120] Specifically, obtaining the yaw angle according to the average speed of the target object in the actual environment, the fluctuating wind speed in the x direction, and the fluctuating wind speed in the y direction includes:

[0121] The yaw angle is calculated by the following method:

[0122] Specifically, the yaw angle is obtained in the following manner:

[0123]

[0124] Wherein, β represents the yaw angle, u′ y represents the fluctuating wind speed in the y direction, u′ x represents the fluctuating wind speed in the x direction, Represents the average speed of the target object in the actual environment.

[0125] In this embodiment, the calculation method of the yaw angle is clarified to ensure the feasibility of the solution, wherein the yaw angle is calculated using the average speed of the target object in the actual environment to ensure the validity of the yaw angle.

[0126] The following is the test comparison data:

[0127] See also Figure 3 and Figure 4 As shown, Figure 3 Schematic diagram of the pulsating wind speed of a car in a road test (natural wind) and a wind tunnel test (reproduced wind). Figure 4 This is a schematic diagram of the pulsating wind speed output by the numerical simulation geometric model of the present application, wherein the horizontal axis represents time (time), the unit is second (sec), and the vertical axis represents the pulsating wind speed Δu, the unit is U. By comparison, it can be found that the pulsating wind speed spectrum simulated and output by the turbulence simulation verification method of the present application is consistent with the pulsating wind speed measured on the actual road and in the wind tunnel.

[0128] See also Figure 5 and Figure 6 As shown, Figure 5 Schematic diagram of the yaw angle of the car in the road test (natural wind) and wind tunnel test (reproduced wind). Figure 6 This is a schematic diagram of the yaw angle output by the numerical simulation geometric model of the present application, wherein the horizontal axis represents time (time), the unit is second (sec), and the vertical axis represents the pulsating wind speed Δψ, the unit is deg. By comparison, it can be found that the yaw angle simulated and output by the turbulence simulation verification method of the present application is also consistent with the yaw angle measured on the actual road and in the wind tunnel.

[0129] See also Figure 7 and Figure 8 As shown, Figure 7 The following is a schematic diagram of the pulsating wind speed spectrum of a car in a road test (natural wind) and a wind tunnel test (reproduced wind). Figure 7 There are three curves in it. The middle curve is the von Karman wind speed spectrum curve. The one above the Karman wind speed spectrum curve is the pulsating wind speed spectrum under the road test. The one above the Karman wind speed spectrum curve is the pulsating wind speed spectrum under the wind test. Figure 8 Schematic diagram of the fluctuating wind speed spectrum output by the numerical simulation geometric model of the present application, wherein the abscissa represents the frequency (f) in Hertz (Hz), and the ordinate represents Δu 2 , unit is (m / s) 2 / Hz. By comparison, it can be found that the pulsating wind speed spectrum simulated and output by the turbulence simulation verification method of the present application is also consistent with the pulsating wind speed spectrum measured on the actual road and in the wind tunnel.

[0130] Please see Table 1 below:

[0131] Road test Wind tunnel testing Results of this Example Turbulence intensity 10％ 8％ 8％ Turbulence integral scale 1.8m 1.8m 1.8m

[0132] Table 1

[0133] It can also be seen from the table that, by comparison, it can be found that the turbulence intensity and turbulence integral scale simulated and outputted by the turbulence simulation verification method of the present application are also relatively consistent with the turbulence intensity and turbulence integral scale measured on actual roads and in wind tunnels.

[0134] In summary, it can be seen that in the embodiments of the present application, a turbulence simulation verification method is provided, which mainly solves the problem that the simulation data simulated by the existing simulation algorithm does not match the actual data by considering the yaw angle, improves the degree of consistency, and thus improves the simulation accuracy. It is mainly used in the field of automotive aerodynamics, specifically for automotive wind resistance and wind noise performance evaluation, and aims to simulate the turbulent wind characteristics measured on actual roads or in wind tunnels through numerical simulation methods. In the early stages of vehicle development, when there is no actual vehicle, this algorithm is used to simulate the actual driving wind environment, thereby improving the effectiveness of the aerodynamic solution for drag reduction and noise reduction, improving design work efficiency and vehicle comfort, and reducing the vehicle's fuel consumption in actual use or increasing the cruising range.

[0135] It should be understood that the order of execution of the steps in the above embodiment does not necessarily mean the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiment of the present invention.

[0136] In one embodiment, a turbulence simulation verification device is provided, and the turbulence simulation verification device corresponds one-to-one to the turbulence simulation verification method in the above embodiment. Fig. 9 As shown, the turbulence simulation verification device includes an acquisition module 101, a simulation module 102, a comparison module 103 and a determination module 104. The detailed description of each functional module is as follows:

[0137] An acquisition module 101 is used to acquire turbulence characteristic parameters of a target object in an actual environment, wherein the turbulence characteristic parameters include turbulence velocity and yaw angle;

[0138] A simulation module 102 is used to substitute the turbulence characteristic parameters as boundary conditions of the entrance of the simulation calculation domain into the numerical simulation geometric model of the target object for simulation analysis to obtain turbulence simulation values;

[0139] A comparison module 103, used for comparing the turbulence simulation value with the actual turbulence value corresponding to the value type;

[0140] The determination module 104 is used to confirm that the simulation accuracy of the numerical simulation geometric model of the target object meets the requirements when the difference between the turbulence simulation value and the actual turbulence value is less than a preset value.

[0141] Furthermore, the turbulent velocity is obtained by:

[0142] Obtaining the average speed and pulsating wind speed of the target object in the actual environment;

[0143] The average velocity and the fluctuating wind speed are added together to obtain the turbulent velocity.

[0144] Furthermore, the pulsating wind speed is obtained by:

[0145] Acquire a first wind speed spectrum synthesis parameter and a second wind speed spectrum synthesis parameter of the target object in the actual environment, wherein the first wind speed spectrum synthesis parameter and the second wind speed spectrum synthesis parameter are respectively determined based on a simulated wind speed spectrum of the target object in the actual environment and a corresponding weight coefficient;

[0146] The pulsating wind speed is calculated according to the first wind speed spectrum synthesis parameter and the second wind speed spectrum synthesis parameter.

[0147] Further, the calculating the pulsating wind speed according to the first wind speed spectrum synthesis parameter and the second wind speed spectrum synthesis parameter includes:

[0148] Preset a first random number and a second random number;

[0149] Calculate a first speed parameter according to the first random number and a first wind speed spectrum synthesis parameter, and calculate a second speed parameter according to the second random number and a second wind speed spectrum synthesis parameter;

[0150] The pulsating wind speed is calculated according to the first speed parameter and the second speed parameter.

[0151] Further, the calculating the pulsating wind speed according to the first speed parameter and the second speed parameter includes:

[0152] The pulsating wind speed is calculated as follows:

[0153]

[0154]

[0155]

[0156] in, represents the pulsating wind speed, x represents the coordinates of the simulation grid loading point, L represents the turbulence integral length, t represents time, k m represents the number of segments, U1 represents the first random number, U2 represents the second random number, N represents the number of samples between the first random number and the second random number, P m,n represents the first wind speed spectrum synthesis parameter, q m,n represents the first wind speed spectrum synthesis parameter, ω m,n represents the angular frequency of the nth wave number in the mth segment of the simulated wind speed spectrum, represents the random phase that changes with time t and coordinates;

[0157] represents the first speed parameter, Represents the second speed parameter.

[0158] Furthermore, the first wind speed spectrum synthesis parameter and the second wind speed spectrum synthesis parameter are determined by:

[0159]

[0160] Among them, P m,n represents the first wind speed spectrum synthesis parameter, q m,n represents the first wind speed spectrum synthesis parameter, a represents the weight coefficient, E(k m ) represents the von Karman wind speed spectrum, and N represents the number of samples between the first random number and the second random number.

[0161] Furthermore, the yaw angle is obtained by:

[0162]

[0163] Wherein, β represents the yaw angle, u′ y represents the fluctuating wind speed in the y direction, u′ x represents the fluctuating wind speed in the x direction, Represents the average speed of the target object in the actual environment.

[0164] For the specific definition of the turbulence simulation verification device, please refer to the definition of the turbulence simulation verification method above, which will not be repeated here. Each module in the above-mentioned turbulence simulation verification device can be implemented in whole or in part by software, hardware and a combination thereof. The above-mentioned modules can be embedded in or independent of the processor in the computer device in the form of hardware, or can be stored in the memory of the computer device in the form of software, so that the processor can call and execute the operations corresponding to the above modules.

[0165] In one embodiment, a computer device is provided, whose internal structure diagram can be as follows: Fig.10 As shown. The vehicle-mounted gateway includes a processor, a memory, and a network interface connected through a system bus. Among them, the processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. When the computer program is executed by the processor, a turbulence simulation verification method is implemented.

[0166] In one embodiment, a computer device is provided, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the following steps when executing the computer program:

[0167] Acquire turbulence characteristic parameters of the target object in the actual environment, wherein the turbulence characteristic parameters include turbulence velocity and yaw angle;

[0168] Substituting the turbulence characteristic parameters as boundary conditions at the entrance of the simulation calculation domain into the numerical simulation geometric model of the target object for simulation analysis to obtain turbulence simulation values;

[0169] Comparing the turbulence simulation value with the actual turbulence value corresponding to the value type;

[0170] When the difference between the turbulence simulation value and the actual turbulence value is less than a preset value, it is confirmed that the simulation accuracy of the numerical simulation geometric model of the target object meets the requirements.

[0171] In one embodiment, a computer-readable storage medium is provided, on which a computer program is stored. When the computer program is executed by a processor, the turbulence simulation verification method mentioned in any of the above embodiments is implemented, which will not be repeated here.

[0172] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program, and the computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to memory, storage, database or other media used in the embodiments provided in this application can include non-volatile and / or volatile memory. Non-volatile memory can include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM) or flash memory. Volatile memory can include random access memory (RAM) or external cache memory. As an illustration and not limitation, RAM is available in many forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM).

[0173] Those skilled in the art can clearly understand that for the convenience and simplicity of description, only the division of the above-mentioned functional units and modules is used as an example. In actual applications, the above-mentioned functions can be distributed and completed by different functional units and modules as needed, that is, the internal structure of the device can be divided into different functional units or modules to complete all or part of the functions described above.

[0174] The embodiments described above are only used to illustrate the technical solutions of the present invention, rather than to limit the same. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that the technical solutions described in the aforementioned embodiments may still be modified, or some of the technical features may be replaced by equivalents. Such modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present invention, and should all be included in the protection scope of the present invention.

Claims

1. A turbulence simulation verification method, characterized in that: include: Acquire turbulence characteristic parameters of the target object in the actual environment, wherein the turbulence characteristic parameters include turbulence velocity and yaw angle; Substituting the turbulence characteristic parameters as boundary conditions at the entrance of the simulation calculation domain into the numerical simulation geometric model of the target object for simulation analysis to obtain turbulence simulation values; Comparing the turbulence simulation value with the actual turbulence value corresponding to the value type; When the difference between the turbulence simulation value and the actual turbulence value is less than a preset value, it is confirmed that the simulation accuracy of the numerical simulation geometric model of the target object meets the requirements.

2. The turbulence simulation verification method according to claim 1, characterized in that: The yaw angle is obtained by: Obtaining the average speed and pulsating wind speed of the target object in the actual environment; Acquire the pulsating wind speed in the x direction and the pulsating wind speed in the y direction according to the pulsating wind speed; The yaw angle is obtained according to an average speed of the target object in an actual environment, the pulsating wind speed in the x direction, and the pulsating wind speed in the y direction.

3. The turbulence simulation verification method according to claim 2, characterized in that: The obtaining the yaw angle according to the average speed of the target object in the actual environment, the fluctuating wind speed in the x direction, and the fluctuating wind speed in the y direction comprises: The yaw angle is calculated by the following method: Wherein, β represents the yaw angle, u′ y represents the fluctuating wind speed in the y direction, u′ x represents the fluctuating wind speed in the x direction, Represents the average speed of the target object in the actual environment.

4. The turbulence simulation verification method according to claim 1, characterized in that: The turbulent velocity is obtained as follows: Obtaining the average speed and pulsating wind speed of the target object in the actual environment; The average velocity and the fluctuating wind speed are added together to obtain the turbulent velocity.

5. The turbulence simulation verification method according to any one of claims 2 to 4, characterized in that: The pulsating wind speed is obtained by: Acquire a first wind speed spectrum synthesis parameter and a second wind speed spectrum synthesis parameter of the target object in the actual environment, wherein the first wind speed spectrum synthesis parameter and the second wind speed spectrum synthesis parameter are respectively determined based on a simulated wind speed spectrum of the target object in the actual environment and a corresponding weight coefficient; The pulsating wind speed is calculated according to the first wind speed spectrum synthesis parameter and the second wind speed spectrum synthesis parameter.

6. The turbulence simulation verification method according to claim 5, characterized in that: The calculating the pulsating wind speed according to the first wind speed spectrum synthesis parameter and the second wind speed spectrum synthesis parameter comprises: Preset a first random number and a second random number; Calculate a first speed parameter according to the first random number and a first wind speed spectrum synthesis parameter, and calculate a second speed parameter according to the second random number and a second wind speed spectrum synthesis parameter; The pulsating wind speed is calculated according to the first speed parameter and the second speed parameter.

7. The turbulence simulation verification method according to claim 6, characterized in that: The step of calculating the pulsating wind speed according to the first speed parameter and the second speed parameter comprises: The pulsating wind speed is calculated as follows: in, represents the pulsating wind speed, x represents the coordinates of the simulation grid loading point, L represents the turbulence integral length, t represents time, k m represents the number of segments, U1 represents the first random number, U2 represents the second random number, N represents the number of samples between the first random number and the second random number, P m,n represents the first wind speed spectrum synthesis parameter, q m,n represents the first wind speed spectrum synthesis parameter, ω m,n represents the angular frequency of the nth wave number in the mth segment of the simulated wind speed spectrum, represents the random phase that changes with time t and coordinates; represents the first speed parameter, Represents the second speed parameter.

8. The turbulence simulation verification method according to claim 5, characterized in that: The first wind speed spectrum synthesis parameter and the second wind speed spectrum synthesis parameter are determined by: Among them, P m,n represents the first wind speed spectrum synthesis parameter, q m,n represents the first wind speed spectrum synthesis parameter, a represents the weight coefficient, E(k m ) represents the von Karman wind speed spectrum, and N represents the number of samples between the first random number and the second random number.

9. The turbulence simulation verification method according to claim 1, characterized in that: The turbulence characteristic parameters also include turbulence intensity and turbulence integral size.

10. A turbulence simulation verification device, characterized in that: include: An acquisition module, used to acquire turbulence characteristic parameters of the target object in the actual environment, wherein the turbulence characteristic parameters include turbulence velocity and yaw angle; A simulation module, used to substitute the turbulence characteristic parameters as boundary conditions of the entrance of the simulation calculation domain into the numerical simulation geometric model of the target object for simulation analysis to obtain turbulence simulation values; A comparison module, used for comparing the turbulence simulation value with an actual turbulence value corresponding to the value type; The determination module is used to confirm that the simulation accuracy of the numerical simulation geometric model of the target object meets the requirements when the difference between the turbulence simulation value and the actual turbulence value is less than a preset value.

11. A computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the computer program, the steps of the turbulence simulation verification method according to any one of claims 1 to 9 are implemented.

12. A computer-readable storage medium storing a computer program, characterized in that: When the computer program is executed by a processor, the steps of the turbulence simulation verification method according to any one of claims 1 to 9 are implemented.

Images