汽车空气阻力和散热性能的仿真及优化
汽车空气阻力和散热性能的仿真及优化/章林凤
设计-研究
doi:10.3969/j.issn.1005—2550,2013.01.009
汽车空气阻力和散热性能的仿真及优化
章林凤
北汽福田汽车股份有限公司集团产品规划本部,北京102206
摘要:基于先进成熟设计理论,以降低汽车空气阻力和提高散热性能为目标,应用PowerFLOW软件.采用基于数字
仿真分析结果与相关基本实验测试结果相结合的方法建立了仿真模型。结果表明,优化方法可以有效减少卡车空气
阻力及提高散热性能。
关键词:仿真及优化;发动机散热性能;燃油消耗;空气阻力.几何模型:冷却模块
中图分类号:U464.138".2
文献标志码:A
文章编号:1005—2550(2013)01—0031—06
PerformanceSimulationandOptimizationofTruckAir
ResistanceandThermalCapability
ZHANGLin-feng
(GroupProductPlanningDeptofBeiqiFotonMotorCo.Ltd.,Beijing102206,China)
Abstract:Inordertoreducevehicleairresistanceandimprovethethermalperformanceand
mature
as
thetargetbasedofmethods
Oil
advanced
Off
design
theory,itapplicatedPowerFLOWsoftwareand
a
used
a
combinationbased
digital
simulationresultswiththeexperimentaltestresultstobuild
can
simulationmodel.Theresultshowsthatoptimizationmethod
effectivelyreducethetruckairresistanceandimprovethethermalpefformance.
and
optimization;engine
cooling
Keywords:simulationmodel:coolingmodule
performance;fuel
consumption;air
resistance;the
geometric
在当今中国汽车市场,消费需求不断升级,国家的排放法规政策也越来越严格,对于整车在价格、性麓、可靠性、舒适性、燃油经济性、加速性、售后服务的及时性、低使用成本等方面的要求越来越高。空气阻力…是影响燃油消耗主要因素之一,所以最大限度地减小整车空气阻力是降低油耗的有效方法,降低油耗的同时也能减少排放并降低使用成本。但这会降低散热系统冷却空气流速.即影响整车散热性能,整车将无法满足在高温环境下的作业。
本文介绍的国内某自主品牌卡车开发项目,在产品实车性能测试阶段,发现整车油耗及散热系统水温偏高,难以形成产品竞争力。最后通过仿真分析,对整车驾驶室顶部导流罩和发动机舱的形状进行优化设计,减少了空气阻力,同时大大降低整车散热系统的水温,提高了产品竞争力。本文中采用的PowerFLOW分析软件是多年以来汽车行业中空气
案列已被证实[”】。
1仿真数学模型
流体仿真使用格子玻尔兹曼方程(LBE)算法进行。这种方法在文献E8,9]中有大量详细描述,这里作简要的概述。
从动力学理论出发,玻尔兹曼方程描述了一种代表性的颗粒以及它们如何演变为流体。玻尔兹曼方程可以写成下列形式,见式(1):
知矿.芒:0
~
(1)
aV
式中:,为速度的概率分布函数;日为碰撞算子。
在一个格子的形式中,这样的方程可以表示为代数方程组的每个状态的概率分布‘:
动力学的标准工具∞],拥有全面的、成熟的热分析
能力,其发动机散热性能的仿真精度在以往的成功
收稿日期:2012—12—20
fi(t+AtlX+ejAt)《(£傅)相(£l髫)
一个格子BGK碰撞项:
(2)
嘶二)=一等m;)骨(£。三)
(3)
・31・
设计・研究
式中:r为弛豫时间和为格子平衡态分布函数,用速度表示为[10l:
肿知棚+争+业2c)(警也…4)
通过定义速度分布函数,液体动力学变量例如密度和速度可以写为:
p(t。;)=∑S(t,;)
(5)
i
孤;)=—l_∑z(矗);。
(6)
P(t,戈)i
通过所有i状态碰撞项求和等于零,LBE进行质量和动量守恒:
∑鼠(矗)=0
(7)
oi(t,;);._0
(8)
使用修改后的基于原始的RNG公式K一占模型建立湍流效应模型…,121:
p瓦Dk=上tr乓Ox『(M)萼ax]+砷£
(9)
这基于湍流模型的LBE,传递出流动历史和上游信息,并包含考虑非线性雷诺应力的高阶项m]。通过解决下面的偏微分方程,评估在温度中的变化:
pcP瓦2_1、,‘pc,瓦DT=之O[(等+龇)婴1+Q(1。)x
lIu’
xOPptr
J’+l+VOx通过修改LAX—Wendroff格式的时间有限差分等来求解湍流和温度方程。
2仿真几何模型
开发facetized格式的卡车几何模型。此几何模型是来自车辆的CAD模型,用于空气动力学以及热学仿真。为了达到最大精度的目的,图1为卡车几何模型。
a1详细车身底板几何模型b)冷却模块系统几何模型
图l卡车几何模型
・32・
汽车科技第1期2013年1月
3仿真的边界条件
3.1
仿真优化工况
空气动力性能和散热性能仿真优化工况如下:a)发动机散热系统性能极限工况:以风扇转速
为3200r/min、车速13.3km/h进行:
b1空气动力学性能工况:以车速100km/h进行仿真计算。
仿真是基于车辆三维模型,仿真分析结果采用设计优化后理论分析数据和试验数据进行比较判定。3.2数字风洞设置
将facetized格式的几何模型放入数字风洞中,
数字风洞带有边界条件和更多工况参数,如指定的工作条件和热交换性能等。数字风洞是用来重现公
路行驶条件,即非常低的障碍和完全移动的地板。因此,卡车的车轮设置为一个旋转的面体边界条件。3.3冷却模块安装
冷却模块安装由中冷器及散热器两个热交换器
组成,这两个热交换器用多孔介质建模。对于热仿真,来自冷却液流动的热量和交换到空气中的热量通过“PowerCOOL”加入耦合。零部件供应商提供的冷却器压降特性以及热的特点,分别如图2和图3所示。
图2
中冷器和水箱的冷却气体压力降曲线
图3中冷器和水箱在不同冷却液质量流量速率下的传热系数
汽车空气阻力和散热性能的仿真及优化/章林凤
散热性能仿真使用的风扇模型为完整的几何形状,风扇是带有风扇转速的模块化模型。在空气动力学仿真中,风扇几何模型已被删除,且除去风扇的影响,通过对风扇模型的内部线删除来实现。
表1
CFD模型的工况
操作环境常温运行
高温运行
环境空气温度,℃2538大气压力,kPa101325lOl325空气流速,I【Ph10013.2风扇转速/r.rain。3005
347l散热器散热,kW52.9散热器的散热率/kg.s。3.593中冷器人口温度,℃176中冷器流量/I【g.s‘1
O.204
4整车空气阻力的仿真及优化
整车的总体空气阻力是影响空气动力性能的主要参数。然而,除了用总体空气阻力数据外,也可以充分利用流场数据,研究环绕整个卡车的空气流动行为,进行阻力的细化、量化描述,来指导设计优化。图4显示:优化前,导流罩的初始设计没有足够的停滞区来转移空气流量。达到避免汽车货厢前方空气存在的目的;除了大型高压区,在驾驶室顶上方的流动显示了一个大的分隔区,导致显著的总压力损失;此外,还有部分高压空气流量流入驾驶室后端和车厢之间的间隙。
一矽
图4环绕卡车几何基准线的流动区域(优化前)
在综合分析结果的基础上,需要重新设计一概念车顶导流罩。基于审美造型因素,重新设计的车项导流罩高度不得超过原车顶导流罩,最好解决办法是设计一种更加陡峭的导流罩角度。
图5中的概念分析表明,优化后重新设计的车顶导流罩以一种更好的方式转移空气流动。因此卡车货厢前方的空气停滞、车项上方及高压空气流量流入间隙大部分得以避免。表面压力显示,关于车厢
设计一砥究
前方的角落部分,仍然有一些需要特别优化的潜力这需要微调的车顶导流罩形状。
图5环绕卡车几何基准线的流动区域(优化后)
相对于原车,通过优化后设计变更获得减少阻力约14%。
由图6可见:驾驶室本身的阻力实际上是增加的,这是由于车顶导流罩角度较陡,造成更高的压力。显然,车厢前方压力较低,其中基线显示风阻进
一步增加,而优化修改后的设计风阻显著减少。
图6卡车一标准化风阻基线和概念卡车
(车顶导流罩优化概念设计)
基于这些结果,采用两种车项导流罩装车,进行道路试验燃油消耗测试数据显示,优化后的设计带来7%左右燃油效率的提高,占风阻降低效率的一半。这与卡车燃油效率研究的结果非常一致,即动力的50%~60%用于克服高速行驶时的空气阻力…。
表2为同一底盘不同高度货厢车型外流场分析。可以看到,相对于原车,无论哪一种车顶导流罩设计,车厢前方上都有较大的空气流体冲击。流体离开导流罩更多地冲向车厢,导致车厢前部的压力更高。尽管如此,导流板优化设计依然贡献了9.1%的改善。
表2标准化风阻对比
%
阻力(占原基本车型百分比)
基本车
加高货厢车
基础车型
1.133
概念车型
0:8611.03
・33・
设计-研究
5整车散热性能的仿真及优化
这部分工作涉及的主要参数是散热器进口冷却液温度、中冷器出口温度以及冷却空气流速。所有的数值分析,在相同的车辆运行工况下(见表1)进行。基准(原车)的仿真结果显示冷却液温度明显高于设计目标,如表3所示。在高温恶劣的工况下运行,整车散热性能测试结果数据显示无法满足性能目标要求。
表3基准仿真结果和冷却性能目标对比
散热性能要素基准仿真分析
目标
环境空气温度/℃
3838进口温度,℃
115.1loo.0
散热器
空气流量/kg・s’1
1.74
进气温度,℃
90.5
散热量,kW
19.722.1
空气流量/kg・s“1.5l中冷器
环境空气温度,℃81.669.O
进口温度,℃
63
表4测量数据和CFD预测数据散热性能要素模拟设计迭代
测试
环境空气温度/℃
3838进水温度/oC
101.9101.0
散热器
空气流量,kg・8-11.8l进气温度,℃
76.7散热量,kW
20.619.5
空气流量/kg・s。
1.56中冷器
环境空气温度,℃77.281.0
进口温度,℃
5l_0
图7发动机窗流场顶视图
・34・
汽车科技第1期2013年1月
从发动机窗流场的详细分析可以看出散热性能
不足的几个原因。如图7所示,冷却模块左、右侧有显著的热流量再循环以及散热模块内逆向流动。此外,风扇罩内的气流量再循环,阻塞了通过散热器的空气流动的途径,因此导致散热性能较差的主要原因是发动机窗、风扇护风罩的形状设计。
从分析看出:在冷却模块的侧面、顶部和底部,增加防止空气流再循环的隔板,可以改善冷却模块的散热性能。比较表3和表4,这个简单的优化设计增加隔板方案确实显著降低进入冷却模块空气的温度,实现了散热器顶部水室温度降低效果,散热器不
会过热并保持车况稳定。从表4可以看出,测得的数据与计算结果吻合良好。
因为在道路行驶测试工况下没有非常精确的控
制测试条件,以及考虑散热器性能的不稳定性,保留
2℃设计余量。如果能在精确定义工况下进行仿真分
析,可以减少这种设计余量。增加隔板方案实施后,进行仿真。图8和图9表明在冷却模块中仍然存在反向流动和漩流问题。因此,进行更多仿真分析用来解决这些问题,寻求全面优化的可行性方案。
图8增加隔板优化,冷却模块仍然存在少量逆向流动(散热器)
图9增加隔板优化,风扇导流罩内仍然存在部分漩流
设计・研究
nalFlowSimulations
to
汽车科技第1期2013年1月
Aerodynamics,Aeroacousticsand
a
[10]D.d’Humieres,P.Lallemand
BGKmodelsforNavier—Stokes
and
Y.H.Quian,“Lattice
ThermalManagementof01—0100.
Pickup
Truck”[C].SAE
2007—
equations”[J].Europhysics
Letters,17(6):479-484,1992.[11]V.Yakhot,and
lysis
S.A.I
[5]A.Alajbegovic,B.Xu,A.KonstantinovJ.AmodeoandW.
Jansen.“SimulationofcoolingairflOWunderdi舶rentdriv—
ing
Orszag,“RenormalizationGroupAna—
ofTurbulence.I.Basic
Theory”[J].Sci.Comput.,1
C.
conditions”[C].SAE
2007一l一0766.《2),3-51,1986.
[6]B.Bhatnagar,D.Schlesinger,A.Alajbegovic,J.BeedyIK.Ho-
rrigan,F.SarrazinandB.Xu,“Predictionandoptimizationof
a
[12]V.Yakhot,V.,S.A.Orszag,S.Thangam,T.Gatski,and
Speziale,“Developmentflowsby
a
ofturbulencemodelsforshear
class8
truckcoolingperformance”[C].SAE
2007-01—
doubleexpansiontechnique”[J].Phys.Fluids
4111.
A,4(7),1510—1520,1992.
and
a
[7]C.L.Lafferty,A.Alajbegovic
hood
Thermal
Simulation
of
K.Horrigan,“Under—
8
Tmck”『C].SAE
[13]H.ChenIS.A.Orszag,I.StaroselskyIand
edAnalogyand
between
Bohzmann
S.Succi,“Expand-
Theory
of
Fluid
Class
Kinetic
2007—01-4280.
Turbulence”[J].J.Fluid
Cooling
Mech.,519—307—314,2004.AirflowofHeavyTrucks—a2008-01-1171.
[8]S.Chen
fluid
andG.D.Doolen,“LatticeBoltzmannmethodfor
Reviewof
Fluid
Mechanics30:
[14]T.Hallqvist,“The
Parametric
flows”[J],Annual
Study”[C].SAE
329-364,1998.
[15]PowerFLOW[16]PowerCOOL
User’SGuide[C].Release4.1,ExaCorpora—
[9]S.Succi,The
icsand
LatticeBohzmannEquationforFluidDynam—
Numerical
Mathematics
and
tion,Boston,Massachusetts,2008.
User’sGuide,Release4.1,ExaCorporation,
Beyond[M],Series
ScientificComputation,ClarendonPress,0xford,2001.
Boston,Massachusetts,2008.
(上接第27页)
FFr运算是对亚抽样的结果进行运算,运算过程由SPE(信号处理引擎)完成。SPE属于内核(e2002335一COl'e)的一部分,具有强大的数学运算能
1.8检测结果的判断
通过程序命令将能量集成修正后的结果与单片机内库存数据进行比较,即可确定是否有爆震发生。
力,能够轻松完成FFF运算过程中最主要环节——
乘加运算(MAC)。SPE配备有32个64位的通用寄存器(GPR),支持单指令多数据操作(SIMD),在每个系统时钟内能够完成两个32位数的算术运算,包括乘加运算(MAC)和无符号乘法运算(MUL),通过
新一代的32位MCUSPC563M64中集成的eQADC模块可以完成增益放大、滤波、A/D转换、抽样滤波等对爆震信号的逐步处理,最后通过SPE的运算集成及相应的程序命令可以完成对爆震的检测。同时,在eTPU的实时控制及eDMA的数据直接传送,可以完成对爆震的定时检测,减少主控制器的干涉,降低系统的负荷。
2结论
编程,足以完成FW运算,满足实际的需要。
爆震信号频率集中在7kHz左右,故此处设计有一个带通滤波器来提取爆震信号。通过程序将该带通滤波器的中心频率设定为6.5~7kHz之中,通带宽度1
kHz
f6.15~7.15kHz),增益为l。
1.7能量集成
能量集成的目的是计算震动的强度,其计算方
法采取将通过F丌计算得到的各特征点的模值(取
平方值)在相应的爆震窗口内进行累加后取其均值,并对结果进行简单的修正,见式(1)。
(1)
参考文献:
[1]程佩清.数字信号处理教程[M].北京:清华大学出版社,
2011.
[2]张欣.LSI数字信号处理[M:.北京:科学出版社,2003.[3]武得钰,傅茂林,李建权,等.火花点火发动机爆震强度
评价指标的研究[J].内燃机学报,1997,15(1):62—69.
[4]高磊,康晓敦.高效实用的片上爆震检测解决方案fEB].
http://www.ed—china.COB2009.11.30.
式中值。
筋i为该点的模