[发明专利]过冷大水滴条件下结冰风洞试验中模型参数的计算方法

摘要

本发明公开了一种过冷大水滴条件下结冰风洞试验中模型参数的计算方法，目的在于解决现有过冷小水滴的结冰试验参数选取方法用于过冷大水滴结冰实验的效果不好，通过飞机的缩比模型上的结冰无法反映真实结冰的情况的问题。本发明在确定模型与真实飞机的缩比关系后，根据理论推导与数值计算的方法，对应于真实的结冰飞行条件，依次给出风洞试验段的空气速度、水滴平均粒径、空气压强、液态水含量、结冰时间、过冷水低温度这几个实验参数。根据换算得到的试验参数，进行飞行器结冰风洞实验，试验结果与真实结果能够满足结冰过程的主要特征的相似要求。可见，采用本发明的计算方法能够得到用于结冰风洞试验中的模型参数，相应参数用于结冰风洞实验能够反映真实结冰的情况。

主权项

过冷大水滴条件下结冰风洞试验中模型参数的计算方法，包括如下步骤：(一)根据被模拟物的要求，确定被模拟物的真实外形和真实气象条件，所述真实外形和真实气象条件包含以下七个参数：真实外形特征长度Lf、飞行速度Vf、空气中的过冷水滴直径df、大气压强Pf、空气中的液态水含量LWCf、结冰总时间tf、来流温度Tf；(二)将被模拟物的外形按照比例缩小，制作一个与被模拟物几何相似的物理模型；(三)根据步骤二制备的物理模型，确定物理模型的模型特征长度Lm，根据模型特征长度Lm和真实外形特征长度Lf，得到实验的缩比数，缩比数的值为(四)计算撞击因子Ki，计算公式如下：该式中，ρd为水密度，σd为水表面张力系数，μd为水粘性系数，Dd为水滴直径，V为来流速度，根据真实气象条件，计算得到Ki的值；(五)根据步骤四得到的Ki值，采用如下公式得到判断因子ε的值：ϵ=0.6Ki>57.7ϵ=0.5Ki<=57.7;]]>(六)根据步骤一中的Vf、步骤三得到的步骤五得到的ε的值，采用如下公式计算得到实验需要的速度Vm：该式中，γ＝1.4；(七)根据步骤一中的Vf、步骤三得到的步骤五得到的ε的值，采用如下公式计算得到实验需要的过冷水滴直径dm：dmdf=(VfVm)1ϵ;]]>(八)由于在开氏温标下，实验需要的温度和真实温度相差相近，令的初始值为1；(九)采用如下公式计算实验需要的压强Pm：PmPf=(VfVm)2(TmTf);]]>(十)采用如下公式计算，得到实验需要的液态水含量LWCm：LWCmLWCf=PmPf/((LmLf)0.5(VmVf)0.5);]]>(十一)采用如下公式计算，得到实验需要的结冰总时间tm：tmtf=(LmLf)1.5/((VmVf)0.5(PmPf)0.5);]]>(十二)采用如下公式计算，得到实验需要的来流温度Tm′：Tm′=Tf+θm-θfb+Vf2-Vm22*Cp,w,]]>该式中，Cp,w是水的比热，Cp,w＝4.184×103J/Kg.K‑1，θ是空气能量传输参数，b是相对温度因子，θm=(273.15-Tm)-rVm22Cp,]]>θf=(273.15-Tf)-rVf22Cp,]]>b=LWCf·Vf·Cp,whc,]]>其中，驻点处的热对流系数k是空气的热传导率；ρ是水的密度；r是恢复因子，r＝Pr1/2，其中Pr＝0.72，Cp是空气的比热，μ是空气动力粘性系数；(十三)采用如下公式计算温度的相对误差Tδ：Tδ＝|T′m‑Tm|/Tm，若Tδ＞0.000001,则令Tm＝T′m，返回到步骤九；如果Tδ≤0.000001,则Tm＝T′m，结束循环。