基于CFD的整車涉水性能仿真預(yù)測(cè)

胡穎 汪曉虎 董丹丹

摘 要：發(fā)動(dòng)機(jī)氣缸進(jìn)水將直接導(dǎo)致車輛熄火甚至發(fā)動(dòng)機(jī)報(bào)廢，因此避免發(fā)動(dòng)機(jī)進(jìn)氣口水入侵是整車涉水性能開(kāi)發(fā)的關(guān)注重點(diǎn)。本文采用流體力學(xué)的兩相流模型和動(dòng)網(wǎng)格模型建立整車涉水仿真模型。通過(guò)實(shí)時(shí)監(jiān)測(cè)發(fā)動(dòng)機(jī)艙水位高度、發(fā)動(dòng)機(jī)進(jìn)氣口的水體積分?jǐn)?shù)和進(jìn)水量來(lái)預(yù)測(cè)整車涉水性能。對(duì)標(biāo)發(fā)動(dòng)機(jī)艙內(nèi)水位高度的實(shí)驗(yàn)與仿真結(jié)果，驗(yàn)證了此方法可應(yīng)用于車輛涉水性能的前期預(yù)測(cè)。整車涉水仿真為發(fā)動(dòng)機(jī)進(jìn)氣口水流路徑優(yōu)化提供設(shè)計(jì)指導(dǎo)。

關(guān)鍵詞：計(jì)算流體力學(xué) 整車涉水仿真 發(fā)動(dòng)機(jī)進(jìn)氣口 兩相流 動(dòng)網(wǎng)格

Abstract：The water in the engine cylinder will directly cause the vehicle to stall or even the engine to be scrapped. Therefore， avoiding the water intrusion of the engine intake is the focus of the development of the vehicle's wading performance. In this paper， a two-phase flow model and a dynamic mesh model of fluid mechanics are used to establish a simulation model of the entire vehicle wading. Through real-time monitoring of the water level of the engine compartment， the water volume fraction of the engine air intake and the water inflow， the water performance of the vehicle can be predicted. The experimental and simulation results of the water level in the standard engine compartment verify that this method can be applied to the early prediction of the vehicle's wading performance. The vehicle wading simulation provides design guidance for the optimization of the water flow path of the engine air intake.

Key words：computational fluid dynamics， vehicle wading simulation， engine intake， two-phase flow， dynamic mesh

1 引言

車輛涉水性能是指車輛在暴雨或者一定深度積水路的惡劣環(huán)境下能夠正常行駛，車輛重要零部件不出現(xiàn)重大失效和功能性損壞的能力。其中，發(fā)動(dòng)機(jī)氣缸進(jìn)水將直接導(dǎo)致車輛熄火甚至發(fā)動(dòng)機(jī)報(bào)廢，因此避免發(fā)動(dòng)機(jī)進(jìn)氣口水侵入是車輛涉水性能開(kāi)發(fā)的關(guān)注重點(diǎn)。

目前車輛涉水性能開(kāi)發(fā)在項(xiàng)目前期主要依靠以往失效模式和經(jīng)驗(yàn)總結(jié)來(lái)預(yù)測(cè)風(fēng)險(xiǎn)并進(jìn)行設(shè)計(jì)規(guī)避，無(wú)法保證預(yù)測(cè)的精準(zhǔn)度，還會(huì)造成設(shè)計(jì)冗余;在項(xiàng)目后期，車輛涉水路駕駛試驗(yàn)是唯一的驗(yàn)證手段，存在驗(yàn)證困難、工程改動(dòng)成本高和更改周期長(zhǎng)等缺陷[1]。

本文采用流體力學(xué)方法模擬整車涉水過(guò)程中的發(fā)動(dòng)機(jī)艙水位高度、發(fā)動(dòng)機(jī)進(jìn)氣口的水體積分?jǐn)?shù)和進(jìn)水量來(lái)預(yù)測(cè)發(fā)動(dòng)機(jī)進(jìn)氣口水侵入風(fēng)險(xiǎn)，為整車涉水性能開(kāi)發(fā)提供有利手段[2-3]。

2 理論模型

2.1 兩相流模型

本文模擬車輛涉水過(guò)程中水和空氣交界面變化，水和空氣屬于兩種連續(xù)介質(zhì)且不相融，故采用VOF兩相流模型。

VOF模型通過(guò)求解連續(xù)性方程（1）獲得不同時(shí)刻第q相流體的體積分?jǐn)?shù)αq來(lái)構(gòu)造和追蹤交界面[4]。

2.2 動(dòng)網(wǎng)格模型

本文采用動(dòng)網(wǎng)格模型模擬車輛在涉水路駕駛試驗(yàn)中的實(shí)際運(yùn)動(dòng)狀態(tài)。動(dòng)網(wǎng)格模型模擬由于流場(chǎng)域邊界剛體運(yùn)動(dòng)或者邊界變形引起的流體域形狀隨時(shí)間變化的流動(dòng)問(wèn)題[5]。

對(duì)任意一個(gè)邊界運(yùn)動(dòng)的控制體V，廣義標(biāo)量φ的守恒方程如下：

式中，ρ為流體密度;u為流體速度矢量;ug為動(dòng)網(wǎng)格的網(wǎng)格速度;為耗散系數(shù);為φ的源項(xiàng)。

3 仿真分析

3.1 網(wǎng)格劃分

整車面網(wǎng)格包含整車外表面、車外底盤(pán)零件和發(fā)動(dòng)機(jī)艙內(nèi)部零件，乘客艙形成一個(gè)封閉體。流體計(jì)算域分為4個(gè)部分：車前區(qū)域、車后區(qū)域、車周圍的核心區(qū)域和過(guò)渡區(qū)域，如圖1所示。核心區(qū)域采用貼體性較好的四面體網(wǎng)格離散并對(duì)發(fā)動(dòng)機(jī)艙區(qū)域加密[6]，其余三個(gè)區(qū)域采用結(jié)構(gòu)化網(wǎng)格離散并在涉水路的水面高度處加密。

3.2 仿真模型及邊界條件

本文采用Fluent軟件進(jìn)行整車涉水仿真。本文采用瞬態(tài)分析，湍流模型選擇k-omega模型，多相流模型選擇VOF模型。本文采用動(dòng)網(wǎng)格模型中的層鋪法，根據(jù)涉水路駕駛試驗(yàn)中的車速指定網(wǎng)格運(yùn)動(dòng)。仿真模型如圖2所示，整車表面、涉水路地面和兩側(cè)壁面都采用無(wú)滑移壁面邊界條件。涉水路上方采用壓力出口邊界條件，相對(duì)壓力為0。發(fā)動(dòng)機(jī)進(jìn)氣口采用壓力出口邊界條件，根據(jù)進(jìn)氣系統(tǒng)壓降定義負(fù)壓。

3.3 流場(chǎng)關(guān)鍵參數(shù)獲取

本文采用多孔介質(zhì)模型來(lái)模擬中冷器、冷凝器和散熱器。水流遇到冷卻模塊受阻，水位會(huì)沿著冷卻模塊前端上升，而發(fā)動(dòng)機(jī)進(jìn)氣口位置一般布置在冷卻模塊之上，是重要水流路徑。

本通過(guò)實(shí)驗(yàn)設(shè)計(jì)，如圖3所示，測(cè)得通過(guò)中冷器、冷凝器和散熱器的水流流速與壓降關(guān)系曲線。

多孔介質(zhì)壓降關(guān)系為：

式中，一次項(xiàng)系數(shù)b和二次項(xiàng)系數(shù)a與慣性阻力系數(shù)和粘性阻力系數(shù)的關(guān)系為：

根據(jù)公式（4）和（5）求得散熱器、冷凝器和中冷器對(duì)水流的慣性阻力系數(shù)和粘性阻力系數(shù)（沿流動(dòng)方向），如表1所示。

4 對(duì)標(biāo)驗(yàn)證

本文對(duì)某車型在涉水路駕駛過(guò)程中的發(fā)動(dòng)機(jī)艙內(nèi)水位高度進(jìn)行實(shí)驗(yàn)與仿真對(duì)標(biāo)。

圖4為仿真獲得的不同時(shí)刻車外水面變化情況。車輛行駛的前4秒，水面受到車輛沖擊，車頭水位高過(guò)格柵頂端，車輛受到水的阻力，車速下降并穩(wěn)定后，車外水面穩(wěn)定在上格柵下方。

圖5為不同時(shí)刻發(fā)動(dòng)機(jī)艙內(nèi)（選取冷卻模塊前端面的水體積分?jǐn)?shù)云圖，上方為冷凝器，下方為中冷器）水位高度。車輛行駛的前4秒，發(fā)動(dòng)機(jī)艙水位很高，水流沖刷了整個(gè)冷卻模塊前端;車速穩(wěn)定后，發(fā)動(dòng)機(jī)艙內(nèi)水位在冷凝器中間高度附近，遠(yuǎn)低于發(fā)動(dòng)機(jī)進(jìn)氣口位置，此時(shí)不存在發(fā)動(dòng)機(jī)進(jìn)氣口水侵入風(fēng)險(xiǎn)。

圖6為車速穩(wěn)定后發(fā)動(dòng)機(jī)進(jìn)氣口下方水位高度的實(shí)驗(yàn)與仿真對(duì)比，實(shí)驗(yàn)與仿真結(jié)果吻合良好，驗(yàn)證了本文仿真模型的可行性和可靠性。

5 應(yīng)用

某車型的發(fā)動(dòng)機(jī)進(jìn)氣口朝下扣在CRFM 導(dǎo)流板的開(kāi)孔上，如圖7所示，通過(guò)仿真分析發(fā)現(xiàn)該布置形式會(huì)導(dǎo)致大量的水持續(xù)涌入進(jìn)氣口。優(yōu)化發(fā)動(dòng)機(jī)進(jìn)氣口水流路徑，在CRFM導(dǎo)流板開(kāi)孔中間增加格柵并在開(kāi)孔下方增加擋板，如圖8所示，優(yōu)化后的發(fā)動(dòng)機(jī)進(jìn)氣口的進(jìn)水量明顯減少，如圖9所示。

6 結(jié)論

1）在整車開(kāi)發(fā)前期，可通過(guò)整車涉水仿真，評(píng)估整車涉水性能，預(yù)測(cè)發(fā)動(dòng)機(jī)進(jìn)氣口水侵入風(fēng)險(xiǎn);

2）通過(guò)車輛涉水過(guò)程中發(fā)動(dòng)機(jī)艙內(nèi)水位高度的實(shí)驗(yàn)與仿真對(duì)標(biāo)，驗(yàn)證了本文仿真模型的可行性和可靠性;

3）在車型開(kāi)發(fā)前期，整車涉水仿真分析對(duì)發(fā)動(dòng)機(jī)進(jìn)氣口及其水流路徑設(shè)計(jì)具有指導(dǎo)作用，降低工程更改成本，縮短整車開(kāi)發(fā)周期。

參考文獻(xiàn)：

[1]王東，陳書(shū)帆，余王培.計(jì)算流體動(dòng)力學(xué)在乘用車涉水性能開(kāi)發(fā)的應(yīng)用及發(fā)展[J].上海汽車，2020，（4）：9-15.

[2]ZHENG X， QIAO X， KONG F. Vehicle wading simulation with STAR-CCM+[A]. Proceedings of the FISITA 2012 World Automotive Congress[C]. Berlin： Springer， 2013. 157-165.

[3]KHAPANE P， GANESHWADE U. Wading simulation-challenges and solutions[C]. SAE Technical Paper 2014-01-0936， 2014.

[4]童亮，余罡，彭政，等.基于VOF模型與動(dòng)網(wǎng)格技術(shù)的兩相流耦合模擬[J].武漢理工大學(xué)學(xué)報(bào)，2008，30（4）： 525-528.

[5]隋洪濤，李鵬飛，馬世虎，等. 精通CFD動(dòng)網(wǎng)格工程仿真與案例實(shí)戰(zhàn)[M]. 北京：人民郵電出版社，2013.

[6]黃勁，余志毅，劉朝勛，等. 基于CFD的兩棲車輛水動(dòng)力導(dǎo)數(shù)的計(jì)算方法[J]. 車輛與動(dòng)力技術(shù)，2014，30（2）：39-43.