Fine reconstruction method of airfoil surface pressure based on multi-source data fusion
-
摘要: 在风洞试验模型表面布置测压孔是获得表面压力分布的重要手段,但受限于空间位置和试验成本,通常难以在复杂模型表面布置足量的测压孔获得完整的表面压力分布信息,直接积分获得的升力和力矩精度不足,因此提出了一种融合稀疏的风洞试验数据和数值计算(CFD)数据的方法,通过较少的风洞测压试验数据获得高精度的压力分布。首先通过本征正交分解技术提取数值计算数据的压力分布低维特征(POD基函数),然后利用稀疏的试验测压数据,通过压缩感知算法获得基函数的坐标,最后将坐标转化到物理空间重构出压力分布。利用定常固定翼型变状态以及变几何变来流状态算例验证该方法的精度,重构结果均能精确匹配试验结果。该重构方法可在一定程度上解决空间受限稀疏观测条件下的分布载荷精细化重构难题。Abstract: Laying out pressure taps on the surface of the wind tunnel test model is an important means to obtain the surface pressure distribution. However, due to the limited space location and experimental cost, it is usually difficult to arrange enough pressure taps on the complex model surface to obtain complete surface pressure distribution information. Hence, the accuracy of the lift and moment calculated by the direct integration may fail to meet expectations. In this paper, to obtain the high-precision pressure distribution through less pressure test data, a method of combining the sparse wind tunnel test data and the numerical simulation data is proposed. Firstly, the proper orthogonal decomposition (POD) technique is used to extract the low-dimensional feature of the pressure distribution of the numerical simulation data, which is called the POD basis. Then, by applying the compressed sensing algorithm, the coordinates of the basis function are obtained with the sparse wind tunnel pressure measurement data, and finally transformed into the physical space to reconstruct the pressure distribution. The accuracy of the method is verified by the steady fixed airfoil variable state or with variable geometry in conjunction with variable flow state examples, and the reconstructed results can accurately match the experimental results. The developed reconstruction method largely solves the problem of fine reconstruction of distributed load under the condition of limited space and sparse observation.
-
表 1 重构误差
Table 1. Reconstruction error
${C_{{l} } }$ ${C_{{m} } }$ 实验 0.3910 0.0030 重构 0.3862 0.0017 Absolute error 0.0048 0.0013 Relative error 0.0122 0.4333 表 2 误差区间
Table 2. Error interval
上限 下限 $C_{l\_{\rm{error}}}$ 0.0050 0.0010 $C_{m\_{\rm{error}}}$ 0.0035 0.0010 表 3 重构误差
Table 3. Reconstruction error
(a) (b) ${C_{l,{\rm{CFD}}} }$ 0.7650 0.6730 ${C_{l,{\rm{exp}}} }$ 0.7040 0.6130 ${C_{l,{\rm{r}}} }$ 0.7058 0.6143 $C_{l\_{\rm{error} } }$ 0.0018 0.0013 Relative error(Cl) 0.0026 0.0021 RMSE 0.0705 0.0655 表 4 重构误差
Table 4. Reconstruction error
(a) (b) (c) ${C_{l,{\rm{exp}}} }$ 0.5330 0.4080 0.1270 ${C_{l,{\rm{r} } } }$ 0.5292 0.4047 0.1269 $C_{l\_{\rm{error} } }$ 0.0038 0.0033 0.0001 Relative error($ {C_{{l}}} $) 0.0071 0.0069 0.0001 ${C_{m,{\rm{exp}}} }$ 0.0015 0.0043 0.0036 ${C_{m,{\rm{r}}} }$ 0.0074 0.0089 0.0050 $C_{m\_{\rm{error} } }$ 0.0059 0.0046 0.0014 Relative error($ {C_m} $) 3.9333 1.0698 0.3888 RMSE 0.0578 0.0380 0.0318 表 5 误差区间
Table 5. Error interval
上限 下限 $C_{l\_{\rm{error} } }$ 0.005 0.001 $C_{m\_{\rm{error} } }$ 0.006 0.001 表 6 重构误差
Table 6. Reconstruction error
(a) (b) (c) ${C_{l,{\rm{exp}}} }$ 1.0240 0.2580 0.5660 ${C_{l,{\rm{r}}} }$ 1.0242 0.2628 0.5627 $C_{l\_{\rm{error} } }$ 0.0002 0.0048 0.0033 Relative error(${C_{{l} } }$) 0.0002 0.0186 0.0058 ${C_{m,{\rm{exp}}} }$ –0.0405 0.0055 –0.082 ${C_{m,{\rm{r}}} }$ –0.0400 0.0035 –0.0795 $C_{m\_{\rm{error} } }$ 0.0005 0.002 0.0025 Relative error(${C_{{m} } }$) 0.0123 0.3636 0.0305 表 7 误差区间
Table 7. Error interval
上限 下限 $C_{l\_{\rm{error} } }$ 0.005 0.001 $C_{m\_{\rm{error} } }$ 0.005 0.001 -
[1] 郑亚青. WDPSS缩比模型的低速风洞测力试验[J]. 华侨大学学报(自然科学版),2009,30(2):119-122.ZHENG Y Q. Force-measuring experiment for the scale model of WDPSS in low-speed wind tunnel[J]. Journal of Huaqiao University (Natural Science),2009,30(2):119-122. [2] 李平,谢艳,杨奇磷. 2.4 m风洞大规模测压试验技术及应用[J]. 实验流体力学,2002,16(2):92-96. doi: 10.3969/j.issn.1672-9897.2002.02.018LI P,XIE Y,YANG Q L. Test technique and application of large-scale pressure measurement in the 2.4 m × 2.4 m transonic wind tunnel[J]. Experiments and Measurements in Fluid Mechanics,2002,16(2):92-96. doi: 10.3969/j.issn.1672-9897.2002.02.018 [3] BELYAEV M,BURNAEV E,KAPUSHEV E,et al. Building data fusion surrogate models for spacecraft aerodynamic problems with incomplete factorial design of experiments[J]. Advanced Materials Research,2014,1016:405-412. doi: 10.4028/www.scientific.net/amr.1016.405 [4] GHOREYSHI M,BADCOCK K J,WOODGATE M A. Accelerating the numerical generation of aerodynamic models for flight simulation[J]. Journal of Aircraft,2009,46(3):972-980. doi: 10.2514/1.39626 [5] 王文正,桂业伟,何开锋,等. 基于数学模型的气动力数据融合研究[J]. 空气动力学学报,2009,27(5):524-528. doi: 10.3969/j.issn.0258-1825.2009.05.004WANG W Z,GUI Y W,HE K F,et al. Aerodynamic data fusion technique exploration[J]. Acta Aerodynamica Sinica,2009,27(5):524-528. doi: 10.3969/j.issn.0258-1825.2009.05.004 [6] WANG X,KOU J Q,ZHANG W W. Multi-fidelity surrogate reduced-order modeling of steady flow estimation[J]. International Journal for Numerical Methods in Fluids,2020,92(12):1826-1844. doi: 10.1002/fld.4850 [7] KOU J Q,ZHANG W W. Multi-fidelity modeling framework for nonlinear unsteady aerodynamics of airfoils[J]. Applied Mathematical Modelling,2019,76:832-855. doi: 10.1016/j.apm.2019.06.034 [8] HE L, ZHOU Y, QIAN W, et al. Aerodynamic data fusion with a multi-fidelity surrogate modeling method[C]//Proc of the 7th European Conference for Aeronautics and Space Sciences. 2017. [9] MIFSUD M,VENDL A,HANSEN L U,et al. Fusing wind-tunnel measurements and CFD data using constrained gappy proper orthogonal decomposition[J]. Aerospace Science and Technology,2019,86:312-326. doi: 10.1016/j.ast.2018.12.036 [10] PERRON C. Multi-fidelity reduced-order modeling applied to fields with inconsistent representations[D]. Atlanta : Georgia Institute of Technology, 2020. [11] RENGANATHAN S A,HARADA K,MAVRIS D N. Aerodynamic data fusion toward the digital twin paradigm[J]. AIAA Journal,2020,58(9):3902-3918. doi: 10.2514/1.J059203 [12] SUN S X,LIU S,LIU J,et al. Wind field reconstruction using inverse process with optimal sensor placement[J]. IEEE Transactions on Sustainable Energy,2019,10(3):1290-1299. doi: 10.1109/TSTE.2018.2865512 [13] ZHAO X,DU L,PENG X H,et al. Research on refined reconstruction method of airfoil pressure based on compressed sensing[J]. Theoretical and Applied Mechanics Letters,2021,11(2):100223. doi: 10.1016/j.taml.2021.100223 [14] LI K,KOU J Q,ZHANG W W. Deep learning for multifidelity aerodynamic distribution modeling from experimental and simulation data[J]. AIAA Journal,2022:1-15. doi: 10.2514/1.J061330 [15] DONOHO D L. Compressed sensing[J]. IEEE Transactions on Information Theory,2006,52(4):1289-1306. doi: 10.1109/TIT.2006.871582 [16] 寇家庆,张伟伟,高传强. 基于POD和DMD方法的跨声速抖振模态分析[J]. 航空学报,2016,37(9):2679-2689.KOU J Q,ZHANG W W,GAO C Q. Modal analysis of transonic buffet based on POD and DMD method[J]. Acta Aeronautica et Astronautica Sinica,2016,37(9):2679-2689. [17] 罗杰,段焰辉,蔡晋生. 基于本征正交分解的流场快速预测方法研究[J]. 航空工程进展,2014,5(3):350-357. doi: 10.3969/j.issn.1674-8190.2014.03.014LUO J,DUAN Y H,CAI J S. A quick method of flow field prediction based on proper orthogonal decomposition[J]. Advances in Aeronautical Science and Engineering,2014,5(3):350-357. doi: 10.3969/j.issn.1674-8190.2014.03.014 [18] 余路,曲建岭,高峰,等. 基于过完备字典的缺失振动数据压缩感知重构算法[J]. 系统工程与电子技术,2017,39(8):1871-1877. doi: 10.3969/j.issn.1001-506X.2017.08.29YU L,QU J L,GAO F,et al. Missing vibration data reconstruction using compressed sensing based on over-complete dictionary[J]. Systems Engineering and Electronics,2017,39(8):1871-1877. doi: 10.3969/j.issn.1001-506X.2017.08.29 [19] 田引黎,杨林华,张鹏嵩,等. 基于半张量积压缩感知的形变数据重构在航天器结构健康监测中的应用[J]. 航天器环境工程,2019,36(2):134-138. doi: 10.12126/see.2019.02.005TIAN Y L,YANG L H,ZHANG P S,et al. Deformation data reconstruction based on semi-tensor compressed sensing in structural health monitoring of spacecraft[J]. Spacecraft Environment Engineering,2019,36(2):134-138. doi: 10.12126/see.2019.02.005 [20] 吴超,王勇,田洪伟,等. 基于盲压缩感知模型的图像重构方法[J]. 系统工程与电子技术,2014,36(6):1050-1056. doi: 10.3969/j.issn.1001-506X.2014.06.06WU C,WANG Y,TIAN H W,et al. Image reconstruction method based on blind compressed sensing model[J]. Systems Engineering and Electronics,2014,36(6):1050-1056. doi: 10.3969/j.issn.1001-506X.2014.06.06 [21] Langley Research Center. 2DN00: 2D NACA 0012 airfoil validation case[EB/OL]. (2021-11-12) [2021-12-02]. https://turbmodels.larc.nasa.gov/naca0012_val.html. [22] Thibert J J, Grandjacques M. Experimental Data Base for Computer Program Assessment[R]. AGARD-AR-138, 1979. [23] HARRIS C D. Two-dimensional aerodynamic character-istics of the NACA 0012 airfoil in the Langley 8 foot transonic pressure tunnel[R]. NASA-TM-81927, 1981. [24] HELTON J C,DAVIS F J. Latin hypercube sampling and the propagation of uncertainty in analyses of complex systems[J]. Reliability Engineering & System Safety,2003,81(1):23-69. doi: 10.1016/S0951-8320(03)00058-9