Benchmark Test of CFD Software Packages for Sunroof Buffeting in Hyundai Simplified Model

Transactions of the Korean Society for Noise and Vibration Engineering.
2014.
Mar,
24(3):
171-179

- Received : December 17, 2013
- Accepted : January 22, 2014
- Published : March 20, 2014

Download

PDF

e-PUB

PubReader

PPT

Export by style

Article

Metrics

Cited by

TagCloud

Sunroof buffeting is one of the most critical issues in the vehicle wind noise phenomena. The experimental approach to solve this issue typically requires a lot of time and resources. To reduce time and cost, the numerical approach could be taken, which can also privide more insights into physical phenomena involved in sunroof buffeting, only if the accuracy in its predictions are guranteed. The benchmark test of various numerical solvers is carried out for the buffeting behavior of a simplified vehicle body, the Hyundai simplified model(HSM). The results of each solver are compared to the experimental measurements in a Hyundai aeroacoustic wind tunnel(HAWT) at various wind speeds. In particular, acoustic response tests were performed and the results were provided prior to all simulations in order to consider the real world effects that could introduce discrepancies between the numerical and experimental approaches. Through this study, most solvers can demonstrate an acceptable accuracy level for actual commercial development and high precision experimental data and computational prediction priories can be shared in order to promote the numerical accuracy level of each numerical solver.
Real world effects
These effects should be taken into consideration to progress the accuracy and effectiveness of the numerical methodologies that are required in the road vehicle development process. Besides these effects, the accuracy issues in the actual development process can be caused by constraints of computational resources. The fully extensive simulation of turbulent flow around an entire car is nearly impossible, because it requires tremendous computational resources. However, it is not actually necessary to predict all properties at every single frequency range in the practical development process. For example, in a simulation of sunroof buffeting, only a peak frequency in which maximum sound pressure occurs is important. Outside the resonance area, it is not of interest and some errors can be acceptable in the practical development process. For the efficiency and accuracy of the numerical approach, computational prediction priority should be decided between S/W providers and commercial users.
In order to overcome these accuracy and efficiency issues, an open benchmark study about the buffeting phenomena was conducted with several CFD solvers. Seven pieces of worldwide commercial software containing major solvers and prototype codes under construction were in attendance for this study. The primary objective of the benchmark test was to predict buffeting frequencies and sound pressure levels at various wind speeds. Computational prediction priorities for the buffeting phenomena can be proposed to improve the efficiency of the simulation. It is assumed that only the damping effect is considered out of the real world effects. In order to apply the damping effect to the simulation, an acoustic response test (ART) was conducted. The results of ART, geometries of the simplified vehicle model and experimental conditions were provided as preliminary data for each numerical study.
Hyundai simplified model(HSM)
The experimental study consisted of three steps. First, the acoustic response test(ART) was carried out to consider the damping effect of the interior model volume. Second, in order to validate basic prediction of flow field, velocity profiles were measured at three points on the upper part of the model using by a pitot tube array at 60 km/h as shown in
Fig. 3
. Finally, the frequencies and sound pressure levels of buffeting were measured through a microphone positioned at the center of the interior volume at various wind speeds from 20 km/h to 100 km/h.
Measurement points of velocity profiles
All experiments were conducted in Hyundai aeroacoustic wind tunnel(HAWT). HAWT is a closed circuit type wind tunnel for aerodynamic and aeroacoustic measurement. Its maximun wind speed is 200 km/h with a 28 m
^{2}
nozzle area and the level of the background noise is 56 dBA at 100 km/h.
^{(1,2)}
.
As the acoustic excitation, a loud speaker was used to generate white noise of very high sound pressure level. According to the position of a response microphone, the acoustic response test was performed in two steps as shown in
Fig. 4
. The position of the reference microphone was fixed to the center point between the loud speaker and the upper surface of HSM. In the first step, the response microphone was located at the center of the interior volume of the model and the upper- side rectangular area was opened. When white noise was radiated, both the reference and response microphones received noise signals at the same time. In the second step, the upper opening area was closed, and the response microphone was located just above the upper surface of the model. In the same manner as in the first step, the sound pressure levels in the two microphones were measured simultaneously while under the radiation of white noise.
Acoustic response test
For each step, the complex transfer function (
ctf
) is defined as
where
g
is signal of a reference microphone,
f
is one of a response microphone,
G_{k}
and
F_{k}
are the Fourier transform of
g
(
t
) and
f
(
t
), the complex conjugate is denoted by ‘*’ and
psd
means power spectral density.
The total complex transfer function of ART can be derived from complex transfer functions (
ctf
) of each step as
From Eq. (2), the resonance frequency and Q-factor could be determined. “Q-factor” describes the characteristic of the damping effect of the whole system. Resonance frequency is the value of frequency when the phase angle is equal to -π/2 in the complex transfer function. From the magnitude of the complex transfer function, the Q-factor can be obtained using the half-power bandwidth method as
where the frequencies,
ω
_{0}
,
ω
_{1}
, and
ω
_{2}
could be determined as shown in
Fig. 5
. In this test, the resonance frequency was 26.7 Hz, and the Q-factor was 10.52 shown in
Fig. 6
.
Half-power bandwidth method
Complex transfer functions of the experimental ART
The buffeting phenomena can be assumed as a Helmholtz resonator, resonance frequency and Q-factor can be calculated by acoustic theories
^{(2)}
as
where
c
is the speed of sound,
S
is the neck area,
V
is the cavity volume,
L
is the neck length and
L
' is determined according to
L
' =
L
+ 0.96 •
S
^{1/2}
.
The theoretical values of the resonance frequency and Q-factor are 26.6 Hz and 48.6 respectively. The resonance frequencies from the theory and experimental data are quite similar. It can be shown that the maximum sound pressure of buffeting would be generated at the resonance frequency. The Q-factors are not as similar, the experimental value is significantly less than the theoretical one. This difference shows the damping effect caused by interior sound absorption materials. These values were provided to each solution maker prior to the numerical simulation in order to consider the real world effects of sound absorption materials.
Peak frequencies and peak sound pressure levels at each velocity from 20 to 100 km/h
Summary of CFD sSolvers
Comparison of experimental and simulation results of velocity profiles of 7 solvers from ‘A’ to ‘G’
Comparison of experimental and simulation results of peak sound pressure level(dB)
Solver A and C can predict the velocity of maximum SPL value, 50 km/h. Solver A, B, C and D can predict the tendencies of SPL. Solver B, C and E can make good agreement with experimental data under 50 km/h.
The capability to predict resonance frequencies was validated as shown in
Fig. 10
. Computational prediction priority can be proposed like this: First, the frequency at the maximum SPL value, second, the tendencies of frequencies according to the velocity and third, the quantities of resonance frequencies.
Comparison of experimental and simulation results of peak frequency(Hz)
Solver A and C can predict the frequency of the maximum SPL value. Solver A and F can predict the tendencies of peak frequencies, in particular, the local decrease at 30 km/h. Solver A, B, C, and E can make good agreement with experiments. Solver ‘G’ could not show enough accuracy because it was still in development.
Every numerical study was performed without any of the results of buffeting tests except the wind tunnel test conditions and ART results. The results of ART were provided before the simulation in order to take the real world effects into consideration.
After this benchmark study, all the results were posted to every participating maker. Computational prediction priorities could be shared among CFD solution suppliers and commercial users. Each solver could be calibrated by the accurate wind tunnel test results. Therefore, the know-how of a computational approach for commercial vehicles could be promoted and some of solvers could progress the accuracy level in the prediction of sunroof buffeting behaviors by modifying their simulation setup. As shown in
Fig. 11
, solver A can obtain much more accurate results from its second simulation than its first. It can be said that solver A can be well-calibrated with very high precision and show sufficient accuracy to apply real vehicle development with various real world effects. Solver B and C show the ability to predict the velocity at maximum value of sound pressure level(50 km/h) and the frequency at maximum SPL value with much better agreement with experimental results than their first numerical study.
Modified simulation results of solver A, B and C
c : Speed of sound ctf : Complex transfer function F : Fourier transform of f G : Fourier transform of g g : Signal of response microphone H : Total complex transfer function L : Neck length S : Neck area V : Cavity volume * : Complex conjugate
# A part of this paper was presented at the KSNVE 2012 Annual Autumn Conference
Munhwan Cho is a senior research engineer in R＆D division of Hyundai Motor Group, who is in charge of development of sound and vibration of a road vehicle.
Hyoung Gun Kim is a senior research engineer in R＆D division of Hyundai Motor Group, who is in charge of development of sound and vibration of a road vehicle. His recent research focus is to develop a virtual test system reproducing actual road driving conditions.
Chisung Oh is a senior research engineer in R＆D division of Hyundai Motor Group, who is in charge of development of sound and vibration of a road vehicle, especially the noise reduction system by active noise control.
Kang-duck Ih is a NVH fellow in R＆D division of Hyundai Motor Group, whose research interest is covering a variety of phenomena related with sound and vibration in a road vehicle.

1. Introduction

The buffeting phenomenon is considered one of the most important issues of wind noise problems. It is very difficult to apply the numerical approach because of the strong coupling between aerodynamic behavior and aeroacoustic propagation. The real world effects of actual vehicle driving conditions can particularly highlight the differences between the numerical method and experimental results. These are challenging obstacles when the numerical approach is applied to the commercial vehicle development process.
One of the real world effects is leakage flow from the vehicle body as shown in
Fig. 1(a)
. Additional flow is made by pressure differences between interior and exterior parts of the driving vehicle. This auxiliary flow can affect the shear layer at the trailing edge of the sunroof opening. Moreover, sound absorption material inside a car can make a damping effect to change the characteristics of buffeting as depicted in
Fig. 1(b)
. It is very important to consider the damping effect, which can affect the shear layers on the sunroof opening area, since a ground vehicle has a lot of sound absorption material such as seats, trim covers and floor carpets. In addition, elastic deformations of vehicle body panels through strong resonance, which is referred to as the compliance effect is shown in
Fig. 1(c)
.
PPT Slide

Lager Image

2. Experimental Results

- 2.1 Experimental Setup

The upper part of a conventional road vehicle body can be simplified to a wedge-like shape for numerical and experimental studies. A Hyundai simple model(HSM) is deviced to reproduce complex vortex structures at important parts of a vehicle including the wind shield glass, roof, A-pillars and the rear wake area. For the buffeting phenomenon, HSM has a hollow interior volume with sound absorption pads on its interior walls and a rectangular opening area on the upper surface. The panels(12 mm) and frames of the body are made of aluminium and the absorption pad is 20 mm thick, as
Fig. 2
.
PPT Slide

Lager Image

PPT Slide

Lager Image

- 2.2 Acoustic Response Test(ART)

In order to consider the damping effect by the interior sound absorption material, one of the important real world effects, complex transfer functions, can be examined through the ART
PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

- 2.3 Buffeting Test

In the buffeting test, peak frequencies and peak sound pressure levels were measured for a range of wind speeds from 20 km/h to 100 km/h, as shown in
Fig. 7
. The peak sound pressure level had maximum value at 50 km/h and the buffeting phenomena became weak at wind speeds over 80 km/h.
PPT Slide

Lager Image

3. Benchmark Study Results

- 3.1 Overview of CFD solvers

Seven commercial solvers attended the benchmark test, which included a variety of software packages covered from world-wide CFD software providers to prototype code under development.
Each solver applied a variety of numerical methods and turbulence models that covered finite volume, element, difference method, and the lattice Boltzmann method as shown in
Table 1
.
Summary of CFD sSolvers

PPT Slide

Lager Image

- 3.2 Validation of Velocity Profiles

When the sunroof area on the upper surface was closed, the velocity profiles were measured at three points as shown in
Fig. 3
. The comparison of experimental and numerical results of each solver is shown in
Fig. 8
. Most solvers can predict the thickness of the boundary layer with enough accuracy, but the shapes of profiles are not matched well in any of the solvers. These discrepancies can introduce some errors into prediction of the buffeting behavior.
PPT Slide

Lager Image

- 3.3 Benchmark Results of Buffeting

A comparison of measurements and numerical studies is shown in
Fig. 9
and
10
. At each velocity range from 20 km/h to 100 km/h, the accuracy of each solver is validated.
The ability to predict peak sound pressure level can be examined by experimental data as shown in
Fig. 9
. In order to apply practical vehicle development process, computational prediction priority can be proposed like this: First, the velocity at the maximum level of sound pressure, second, the tendencies of sound pressure levels according to the velocities, and third, the quantities of peak sound pressure levels.
PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

4. Conclusion

The digital development of commercial vehicles can be crucial to the reduction of development time and costs. But accuracy issues should be overcome through the application of highly efficient simulations. To solve the issues, real world effects and prediction priorities should be considered. Therefore, this open benchmark study was conducted to minimize the issues.
As the main theme of this study, the prediction of the sunroof buffeting behaviors were presented, which has strong coupling of aerodynamics and aeroacoustics. To consider the real world effects, an acoustic response test(ART) was carried out and the results were attained before each solver carried out its simulation.
Seven solvers tried to simulate the acoustic behavior and some of them showed enough accuracy to apply to the process of commercial vehicle development. All of the results of each simulation were validated by high precision experimental measurements in a Hyundai aeroacoustic wind tunnel(HAWT). Almost all solvers could modify the errors of the first trial and get much more accurate results that should be described in the following papers. In this benchmark test, a variety of the CFD solvers could be calibrated to promote their accuracy level to the standard required in the actual process of road vehicle development.
Nomenclature

Acknowledgements

The authors would like to thank ANSYS India ＆ Korea, CD-adapco, Cedic Co., ESI Korea, EXA, Flow&Noise Inc. and Next Limit Technologies(listed in alphabetical order) for their supports of the numerical simulation.

BIO

Crouse B.
,
Senthooran S.
,
Balasubramanian G.
,
Freed D.
,
Nolting S.
,
Mongeau L.
,
Hong J. S.
2005
Sunroof Buffeting of a Simplified Car Model: Simulations of the Acoustic and Flow-induced Responses, SAE Paper No. 2005-01-2498

Kinsler L.
,
Frey A.
,
Coppens A.
,
Sanders J.
2000
Fundamentals of Acoustics
4th ed.
John Wiley & Sons Inc

Citing 'Benchmark Test of CFD Software Packages for Sunroof Buffeting in Hyundai Simplified Model
'

@article{ SOJDCM_2014_v24n3_171}
,title={Benchmark Test of CFD Software Packages for Sunroof Buffeting in Hyundai Simplified Model}
,volume={3}
, url={http://dx.doi.org/10.5050/KSNVE.2014.24.3.171}, DOI={10.5050/KSNVE.2014.24.3.171}
, number= {3}
, journal={Transactions of the Korean Society for Noise and Vibration Engineering}
, publisher={The Korean Society for Noise and Vibration Engineering}
, author={조, 문 환
and
오, 치 성
and
김, 형 건
and
이, 강 덕}
, year={2014}
, month={Mar}