Nghiên cứu mô phỏng hiện tượng ồn khí động do xe ô tô chuyển động gây ra

Nghiên cứu này giải quyết bài toán tương tác giữa dòng chảy bao vật thể và hiện tượng phát sinh, lan truyền sóng âm trong không khí. Các tính toán mô phỏng được thực hiện trên phần mềm Siemens-StarCCM+ cho mô hình mẫu xe MERCEDES BENZ S-CLASS 2013.

HỘI NGHỊ KH&CN CƠ KHÍ - ĐỘNG LỰC 2021 NGHIÊN CỨU MƠ PHỎNG HIỆN TƯỢNG ỒN KHÍ ĐỘNG DO XE Ô TÔ CHUYỂN ĐỘNG GÂY RA A NUMERICAL MODELING STUDY FOR AUTOMOTIVE AERODYNAMIC NOISE PHẠM LÊ BẮC, PHẠM VĂN SÁNG* Viện Cơ khí Động lực, Trường Đại học Bách khoa Hà Nội *Email liên hệ: sang.phamvan@hust.edu.vn Tóm tắt Tiếng ồn yếu tố quan trọng nhà sản xuất xe ô tô quan tâm ảnh hưởng tới thoải mái hành khách môi trường xung quanh Bài báo nghiên cứu tiếng ồn sinh trình chuyển động xe tơ Nghiên cứu giải tốn tương tác dòng chảy bao vật thể tượng phát sinh, lan truyền sóng âm khơng khí Các tính tốn mơ thực phần mềm Siemens-StarCCM+ cho mơ hình mẫu xe MERCEDES BENZ S-CLASS 2013 Kết mô cho phép đánh giá độ ồn xe vận tốc di chuyển lớn Kết trường vận tốc trường áp suất dùng để giải thích nguyên nhân gây tiếng ồn, từ đề xuất biện pháp để giảm tiếng ồn Từ khóa: Ồn khí động học, tính tốn động lực học dịng chảy Abstract Noise is one of the most important factors that is interested by car manufacturers due to it’s effect on the comfort of passengers and the surrounding environment This research focuses on the effects of noise that is produced during the movement of cars We examine the interaction problem between the flow through a object and the phenomenon of generating and propagating sound waves in the air Simulation calculations are performed by Siemens-StarCCM+ software for MERCEDES BENZ S-CLASS 2013 models The obtained simulation results allow to assess the noise level in the vehicle at high speeds The results of the velocity and pressure fields are used to explain for the noise source According to the results, we offer several solutions to reduce the aerodynamic noise Keywords: Aerodynamic noise, CFD 368 Giới thiệu Ngày nay, thoải mái hành khách tham gia giao thông ngày coi trọng ngành công nghiệp ô tô ngày phát triển nhu cầu sử dụng khách hàng ngày gia tăng Tiếng ồn khí động học vấn đề hãng xe trọng cải tiến để nâng cao chất lượng sản phẩm Tại vận tốc lớn, dịng chảy khơng khí qua xe tạo tiếng ồn khí động học khiến hành khách xe cảm thấy khơng thoải mái Tiếng ồn khí động học chia làm hai loại: Tiếng ồn kẽ hở tiếng ồn nhiễu động gió [5] Tiếng ồn kẽ hở xảy vỏ xe làm kín hay cửa sổ hai bên sổ trời mở làm cho dịng gió qua khe hở vào bên khoang cabin Tiếng ồn nhiễu động gió gây va đập dịng khơng khí với vỏ xe, đặc biệt va đập phận như: Gương chiếu hậu, giá để hàng, ăng-ten, trụ A, … Các va chạm làm cho dịng khơng khí qua xe bị nhiễu động mạnh gây tiếng ồn lớn [6] Với phát triển máy tính ngày cho phép khảo sát tiếng ồn khí động học mẫu xe từ đưa cải tiến phù hợp nhằm hạn chế ảnh hưởng tiếng ồn khí động học Hơn phương pháp mô giúp thực tiện lợi hơn, có kết nhanh tiết kiệm so với thử nghiệm thực tế Phương pháp mô số Để thực mô ảnh hưởng tiếng ồn khí động học lên hành khách xe, lựa chọn phần mềm mô Siemen-STAR CCM+ mẫu xe mô mẫu MERCEDES BENZ SCLASS 2013 Từ việc mô đưa liệu trường vận tốc, phân bố áp suất, mức cường độ âm điểm khảo sát, thơng qua đưa kết luận nguyên nhân, mức độ ảnh hưởng đề xuất số giải pháp cải thiện vấn đề tiếng ồn khí động học Việc mơ trường vận tốc, trường âm thanh, phân bố áp suất qua việc sử dụng phần mềm STAR CCM+ dựa phương pháp tính toán động học chất lỏng (CFD) Phương pháp bao gồm SỐ ĐẶC BIỆT (10-2021) HỘI NGHỊ KH&CN CƠ KHÍ - ĐỘNG LỰC 2021 phương pháp: mơ dịng xốy lớn (LES), mơ dịng xốy khơng tập trung (DES), Reynolds Averaged Navier Stokes (RANS) [7] Ở mô sử dụng phương pháp Reynolds Average Navier Stokes để tính tốn thơng số dịng chảy 2.1 Phương trình Reynold Average NavierStoke ( RANS) [1] Phương trình RANS phương trình biểu diễn thời gian trung bình chuyển động dịng chảy chất lỏng khơng gian, chúng biểu diễn phương trình sau: + + ̅ + + ̅ + + ̅ + + + + + = + = − = − + + ̅ − + + ∆ − (1) ̅ + ∆ ̅− − ̅ + Δ ̅ + ∆ − − (4) = (∆) +) ∆ = (∆ ∆ ∆ ) 2.3 Động lực học mơ hình SLM (DSLM) Vấn đề gặp phải sử dụng mơ hình lưới SLM cần số Cs mà số đo đạc phương pháp thực nghiệm dẫn đến kết mơ thiếu xác Khi sử dụng phương pháp DSLM Cs khơng số ↔ − (∆) = ↔ (| | (9) ) − (10) Trong đó: áp suất dịng xốy lưới sau lọc hai lần biểu diễn phương trình sau đây: ↔ = ↔↔ - (11) 2.4 Phương pháp mô âm Phương trình vận tốc âm lan truyền mơi trường chất lỏng: = Trong đó: (12) +) tỉ số nhiệt dung trạng thái đẳng áp đẳng tích mơi trường chất lỏng khảo sát; +) R số khí phổ quát; nhiệt độ tuyệt đối (Kelvin) Sau tham số vận tốc, tần số, áp suất giải từ phương trình trên, phần mềm tiếp tục chuyển tham số dạng âm Phương trình FW-H ứng dụng để tính tốn mức cường độ âm, trường âm dịng khơng khí gây bề mặt vỏ xe Phương trình FW-H biểu diễn sau: (7) Trong đó: ∆ kích thước lưới (độ dài thành phần lưới theo phương x,y,z ∆ , ∆ , ∆ ), hệ số thực SLM ↔ +) tenxơ ứng suất dịng xốy xác định phương trình sau: +) Mơ hình lưới SLM dạng mơ hình lưới biểu diễn mối quan hệ độ dài vận tốc theo phương trình sau: = ( Δ) | ̅| (5) | ̅| = (6) (8) Trong đó: (3) 2.2 Mơ hình lưới dạng mơ Smagorinsky - Lily (SLM) [1] [3] SỐ ĐẶC BIỆT (10-2021) =− (2) Từ phương trình ta có cơng thức chung khơng gian: = mà thay vào chúng xác định việc giải nghiệm phương trình chuyển động xốy Từ phương pháp động học mơ hình SLM cho kết cấu trúc dòng xốy xác việc sử dung mơ hình SLM đơn Hệ số Cs DSLM xác định sau: { −∇ + = {[ ( ( )=0 ℎ − + ( ] ( ) + ( , )>0 d(f ) = − ¶H ( f ) ¶f – ) ( )]} ( ) (13) (14) (15) rong đó: c vận tốc âm thanh, p áp suất âm thanh, t thời gian, mật độ điểm vơ cùng, 369 HỘI NGHỊ KH&CN CƠ KHÍ - ĐỘNG LỰC 2021 vận tốc dòng vận tốc nguồn âm bề mặt vỏ xe, vecto đơn vị, tenxơ ứng suất Phương pháp dùng để tổng hợp tham số tính để tính tốn phân bố lan truyền âm Sau phương trình FW-H giải, có kết mức áp suất âm thanh, mức cường độ âm dòng khơng khí qua vỏ xe gây 2.5 Mơ hình mơ Mơ hình sử dụng mơ mẫu MERCEDES BENZ S-CLASS 2013 2.6 Mơ hình lưới tính tốn điều kiện biên bề mặt Lưới sử dụng mơ hình mơ lưới Trimmed Cell Mesh (các lưới có dạng hình hộp chữ nhật) Phương trình biểu diễn mức cường độ âm dạng Curle Noise: ,( , )= ( ) ∮ ( ) ( , − ) ( ) (16) Trong đó: + +) , mật độ âm thanh; vận tốc truyền âm môi trường khảo sát; +) − a) Hình Mơ hình mẫu xe: (a) Mơ hình mẫu xe thực tế; (b) Mơ hình mẫu xe mơ STAR thời gian lan truyền; CCM+ +) p áp suất mặt; Bảng Thơng số xe khảo sát +) y vị trí điểm nguồn âm; +) r = |x-y| khoảng cách điểm khảo sát nguồn âm; +) n vecto pháp tuyến hướng Phương trình biểu diễn mức cường độ âm dạng Proudman Noise: Trong đó: (17) +) số tương quan với vận tốc theo phương j; +) u bình phương trung bình vận tốc thành phần; +) mật độ; +) vận tốc âm Tổng quan lý thuyết toán học sử dụng mơ phỏng: Đầu tiên dịng tính tốn phương pháp RANS Sau thực tổng hợp thông số phân bố áp suất, vận tốc dịng khí Cuối cùng, dựa vào thơng số có để giải phương trình FW-H Kết thu mô trường âm tiếng ồn khí động học bề mặt xe Chiều dài Chiều rộng Chiều cao Chiều dài trục sở Khoảng sáng gầm xe 4.940 1.960 1.773 2.933 192 Kích thước lưới bề mặt mơ hình xe 10mm khơng gian miền khảo sát 20mm Lưới làm mịn 50% so với kích thước vùng khơng gian gương chiếu hậu khoảng không gian vùng khoảng sáng gầm xe Số lượng lưới vùng khảo sát 7,21x10 lưới, khu vực gương chiếu hậu 5,53x105 lưới khu vực khoảng sáng gầm xe 1,15x10 lưới Điều kiện biên vùng xác định theo Bảng ảng Điề ụ ệ Biên Biến vật lý Loại điều kện biên áp dụng Bề mặt mơ hình xe Vật rắn Wall Biên dịng vào Khí Velocity inlet Biên dịng Khí Presure Outlet Vật rắn Wall Bề mặt xung quanh miền khảo sát 370 Kích thước (mm) Thơng số +) x vị trí khảo sát; = b) SỐ ĐẶC BIỆT (10-2021) HỘI NGHỊ KH&CN CƠ KHÍ - ĐỘNG LỰC 2021 Vận tốc dòng vào bề mặt biên dòng vào 33,33m/s tương đương với vận tốc 120km/h Thời gian mô (timestep) 0,1s Kết thảo luận a) Kết mô thu cho thấy tồn hai vùng xoáy với vận tốc thấp hình thành khu vực phía sau gương chiếu hậu Vịng xốy vận tốc thấp vùng khơng gian gương chiếu hậu bề mặt xe, dịng chảy qua bề mặt ngồi gương chiếu hậu tạo thành vùng xốy có vận tốc lớn Các dịng chảy sau tạo vùng xốy tiếp tục tăng tốc va đập bề mặt xe gây tiếng ồn khí động học đáng kể Độ lớn vận tốc vùng gương chiếu hậu bị biến đổi lớn Trong vận tốc dòng vào 33,33m/s (120km/h) qua bề mặt gương chiếu hậu tạo dao động độ lớn vận tốc từ khoảng 9,64m/s lớn 67,5m/s Áp suất ghi vùng dao động từ -1,8x10 Pa tới 1,15x103 Pa Các vùng có giá trị áp suất âm chủ yếu vùng xốy phía sau gương chiếu hậu Chính điều tạo thay đổi lớn cấu trúc dịng hình thành xốy vùng gương chiếu hậu a) b) b) c) Hình Biểu diễn lưới mô phỏng: (a) Lưới làm mịn khu vực gương chiếu hậu; (b) Lưới làm mịn khu vực khoảng sáng gầm xe; (c) Tên biên vật lí định dạng 3.1 Ảnh hưởng bề mặt xe tới trường vận tốc trường áp suất Kết trường vận tốc sử dụng phương pháp RANS qua mặt cắt vùng gương chiếu hậu biểu diễn Hình SỐ ĐẶC BIỆT (10-2021) c) Hình Trường vận tốc qua mặt cắt vùng gương chiếu hậu theo thời gian: (a) t=1,1s; (b) t=1,3s; (c) t=1,5s 371 HỘI NGHỊ KH&CN CƠ KHÍ - ĐỘNG LỰC 2021 a) b) Hình Biểu diễn dịng qua gương chiếu hậu: a, Vận tốc dòng; b, Áp suất dòng Khi dòng chảy qua trụ A, luồng khơng khí chảy bám vào cửa kính bên kính chắn gió phía a) b) 372 trước từ phân hai luồng tách kéo theo vùng xốy hình thành phần giao hai dòng trụ A Dao động vùng xốy có tần số dao động khơng lớn vùng xốy vùng gương chiếu hậu vùng xoáy rộng hơn, trải dọc theo trụ A nên chúng gây lượng tiếng ồn khí động học đáng kể Mặc dù hình dạng cấu trúc dòng bị biến thiên nhỏ so với vùng gương chiếu hậu độ lớn vận tốc biến thiên dao động khoảng 9,71m/s đến 68,6m/s vận tốc -1,35x10 3Pa tới 309Pa áp suất Các vùng có giá trị vận tốc, áp suất biến thiên lớn tập trung gần sát bề mặt trụ A Qua phần kết mô thấy phần lớn dịng chảy sau qua bề mặt xe giảm đáng kể vận tốc hình thành dao dộng vùng xe Các dịng chảy vùng xe khơng tạo xốy hai vùng trước chúng hình thành dao động với biên độ tần số tương đối lớn, từ phần dịng chảy hội tụ phía sau xe ảnh hưởng tương đối lớn tới mức độ ồn người ngồi hàng nghế sau Phần gầm xe phần có biến đổi dịng lớn độ lớn diện tích vùng bị biến đổi Từ kết mơ Hình ta thấy rõ thay đổi trường vận tốc áp suất vùng gầm xe a) b) Hình Biểu diễn dịng qua trụ A: Hình Biểu diễn dòng qua bánh xe: a, Vận tốc dòng; b, Áp suất dòng a, Vận tốc dòng; b, Áp suất dòng SỐ ĐẶC BIỆT (10-2021) HỘI NGHỊ KH&CN CƠ KHÍ - ĐỘNG LỰC 2021 Dịng chảy bị biến động va đập với gầm xe bánh xe gây nhiễu động dòng vùng diện tích lớn kéo dài từ phần đầu xe khoảng lớn phía sau xe Có thể thấy vận tốc áp suất biến thiên lớn nhiều so với vùng khảo sát trước Độ lớn dao động vận tốc vùng tương đối giống hai vùng trước khoảng dao động vận tốc vùng 9,91m/s 69,4m/s Trong khoảng dao động áp suất từ -1,75x103Pa tới 839Pa So sánh cấu trúc dòng chảy vùng khơng gian phía sau gương chiếu hậu vùng khơng gian xung quanh trụ A qua hai Hình Hình ta thấy cấu trúc xốy vùng phía sau gương chiếu hậu rõ ràng nhiều so với vùng không gian trụ A Hiện tượng cho ảnh hưởng hình dạng gương lên mức độ ồn khí động học lớn nhiều so với trụ A (2) điểm vùng cửa số bên phía trước - nơi gần với tai người lái, điểm điểm gần với hành khách ngồi phía sau, điểm điểm nằm vùng xốy phía sau gương chiếu hậu Các điểm khảo sát biểu diễn Hình 11 Hình Cấu trúc xốy vùng bánh xe phía trước Hình 10 Cấu trúc xốy vùng bánh xe phía sau Hình Cấu trúc dịng vùng trụ A Hình 11 Các điểm thực khảo sát mơ Hình Cấu trúc xốy vùng gương chiếu hậu mức độ tiếng ồn khí động học Như trình bày phần trước độ biến thiên vận tốc vùng gầm xe lớn độ lớn lẫn vùng ảnh hưởng Các kết cấu trúc dòng vùng biểu diễn Hình Hình 10 cho thấy dịng chảy vùng có cấu trúc xốy phức tạp nhiều so với vùng kính chiếu hậu vùng trụ A Kết phân tách luồng qua vị trí bề mặt xe mơ qua phương pháp RANS (1) chuyển thành liệu cường độ âm mức cường âm qua việc tính tốn phương trình Flows William-Hawking (13) thơng qua hai kiểu thị Curle Acoustic (16) Proudman Acoustic(17) 3.2 Kết thảo luận mức cường độ âm vị trí khác Các kết mức cường độ âm thu phù hợp với kết trường vận tốc, trường áp suất cấu trúc dòng Mức cường độ âm ồn khí động học lớn vùng có biến đổi dịng chảy lớn Khảo sát mức cường độ âm điểm (1) SỐ ĐẶC BIỆT (10-2021) 373 HỘI NGHỊ KH&CN CƠ KHÍ - ĐỘNG LỰC 2021 Sự biến thiên độ lớn vận tốc, áp suất, cấu trúc dòng xảy lớn khu vực gầm xe kết mức cường độ âm vùng lớn nhất, với mức cường độ âm lớn tính tốn vùng 83,3dB, vùng có tiếng ồn lớn phân bố rộng, gần bao trọn vùng gầm xe Đối với khu vực gương chiếu hậu, xoáy xuất sau khu vực gương chiếu hậu dẫn đến kết mức cường độ âm tính tốn vùng lớn nhất, đạt ngưỡng 87,4dB Tuy nhiên vùng gây tiếng ồn vùng diện tích khơng q lớn, kéo dài từ gương chiếu hậu tới hết phần trước xe Càng xa khu vực xoáy mức cường độ âm giảm dần mức gần 0dB xa vùng xoáy Đối với vùng trụ A- vùng có biến thiên dịng nhỏ so với hai khu vực lại, mức cường độ âm lớn ghi nhận 73,3dB (thấp so với hai vùng cịn lại) Vùng diện tích gây tiếng ồn chủ yếu tập trung sát phần trụ A kính chắn gió phía trước, khơng lan rộng hai vùng cịn lại Hình 13 Phân bố mức cường độ âm bề mặt kính So sánh phân bố tiếng ồn do gương chiếu hậu gây mô nghiên cứu mô thực Papoutsis-Kiachagias E M., Magoulas N., Mueller J [2] Hình 14 Hình 15 ta thấy rõ tương đồng phân bố âm ồn khí động học Các hiển thị phân bố mức cường độ âm lớn vị trí phía sau gương chiếu hậu bề mặt kính xe đồng a) a) b) b) Hình 14 Phân bố mức cường độ âm ồn khí động học khu vực gương chiếu hậu: a) Bài nghiên cứu c) Hình 12.Phân bố mức cường độ âm mặt cắt: b) Bài nghiên cứu Papoutsis-Kiachagias E M., Magoulas N., Mueller J [2] (a) mặt cắt trụ A; (b) mặt cắt gương chiếu hậu; (c) mặt cắt gầm xe 374 SỐ ĐẶC BIỆT (10-2021) HỘI NGHỊ KH&CN CƠ KHÍ - ĐỘNG LỰC 2021 a) Điểm điểm gần vị trí gương nhất- nơi xảy biến đổi lớn dòng chảy so với điểm khảo sát lại nên mức cường âm điểm tăng lên nhanh có giá trị lớn thời điểm ổn định (đạt ngưỡng khoảng 73dB) Điểm điểm hai vị trí gần tai người lái nên ảnh hưởng lớn tới thoải mái tai người lái Hai điểm có tốc độ tăng mức cường độ âm gần giống nhau, nhiên điểm tốc độ tăng mức cường độ âm nhanh điểm gần với vị trí vùng xốy Lúc ổn định mức cường độ âm điểm (khoảng 57dB) lớn so với điểm (khoảng 50dB) Hai điểm hai điểm gần với tai người ngồi phía sau Mức cường độ âm thời điểm ổn định hai điểm là: 59dB 70dB Thông thường, mức cường độ âm mà tai người chịu khoảng 125dB trở xuống, nhiên ngưỡng khoảng 60dB trở lên bắt đầu gây cảm giác khó chịu cho người nghe [4] Như vậy, âm điểm khảo sát truyền trực tiếp tới tai gây cảm giác khó chịu cho người lái lẫn hành khách ngồi phía sau b) Hình 15 Phân bố mức cường độ âm ồn khí động học khu vực kính cửa bên phía trước xe: a) Bài nghiên cứu b) Bài nghiên cứu Papoutsis-Kiachagias E M., Magoulas N., Mueller J [2] Kết luận Vấn đề tiếng ồn dòng chảy bao quanh ô tô nghiên cứu thông qua mơ số tính tốn dịng chảy bao quanh xe tiếng ồn dao động áp suất dịng chảy Kết mơ cho biết trường âm tạo trình xe di chuyển, nguồn phát âm khảo sát Mơ hình mơ số tiếng ồn khí động ứng dụng thiết kế hình dạng khí động học xe nhằm giảm tiếng ổn xe gây trình di chuyển, tạo thoải mái cho hành khách xe TÀI LIỆU THAM KHẢO [1] Xingjun Hu, Peng Guo, Zewei Wang, Jingyu Wang, Mo Wang, Jia Zhu and Dejiu Wu, Calculation of External Vehicle Aerodynamic Noise Based on LES Subgrid Model Hình 16 Sơ đồ biểu diễn thay đổi mức cường độ âm điểm khảo sát Hình 16 biểu diễn thay đổi mức cường độ âm theo thời gian điểm khảo sát (Hình 11) Qua sơ đồ ta thấy thời điểm ban đầu dịng khí bắt đầu làm cho mức cường độ âm ồn khí động học điểm khảo sát tăng dần lên Sau thời gian, dòng chảy ồn định, mức cường độ âm từ định Tuy nhiên, dịng chảy lúc tiếp tục dao động nhỏ lúc trước gây biến động nhỏ biên độ mức cường độ âm điểm SỐ ĐẶC BIỆT (10-2021) [2] Papoutsis-Kiachagias E M., Magoulas N., Mueller J., et al Noise reduction in automobile aerodynamics using a surrogate objective function and the continuous adjoint method with wall functions Computers and Fluids, Vol.122, pp.223232, 2015 [3] Ye-gang Chen , Xiang-hong Wang , Yao-ming Zhou, Numerical investigation on aerodynamic noises of the lateral window in vehicles Journal of Vibroengineering, Vol.19, Issue 8, pp.6502-6518, 2017 375 HỘI NGHỊ KH&CN CƠ KHÍ - ĐỘNG LỰC 2021 [4] Papoutsis-Kiachagias E M., Magoulas N., Mueller J., et al Noise reduction in automobile aerodynamics using a surrogate objective function and the continuous adjoint method with wall functions Computers and Fluids, Vol.122, pp.223232, 2015 [5] Yang, Z.D.; Gu, Z.Q.; Xie, C.; Zong, Y.Q.; Jiang, C.M Experimental research on noise reduction mechanism of a groove spoiler in vehicle sunroof J Hunan Univ (Nat Sci.), Vol.45, pp.26-34, 2018 376 [6] GEORGE A R Automobile aero-acoustics[G] AIAA Paper - 1067, 1989 [7] MuradN., Naser J., Alam F., et al Computational fluid dynamics study of vehicle A-pillar aeroacoustics Applied Acoustics, Vol.74, Issue 6, pp.882-896, 2013 Ngày nhận bài: Ngày nhận sửa: Ngày duyệt đăng: 30/6/2021 05/8/2021 19/8/2021 SỐ ĐẶC BIỆT (10-2021) ... tiếng ồn do gương chiếu hậu gây mô nghiên cứu mô thực Papoutsis-Kiachagias E M., Magoulas N., Mueller J [2] Hình 14 Hình 15 ta thấy rõ tương đồng phân bố âm ồn khí động học Các hiển thị phân bố... âm khảo sát Mơ hình mơ số tiếng ồn khí động ứng dụng thiết kế hình dạng khí động học xe nhằm giảm tiếng ổn xe gây trình di chuyển, tạo thoải mái cho hành khách xe TÀI LIỆU THAM KHẢO [1] Xingjun... tai gây cảm giác khó chịu cho người lái lẫn hành khách ngồi phía sau b) Hình 15 Phân bố mức cường độ âm ồn khí động học khu vực kính cửa bên phía trước xe: a) Bài nghiên cứu b) Bài nghiên cứu