BOS测量密度场的泊松方程源项快速计算方法

  • 张建 ,
  • 张征宇 ,
  • 杨洋 ,
  • 钱丰学 ,
  • 李小霞 ,
  • 王材钢 ,
  • 罗家杰
展开
  • 1. 西南科技大学
    2. 中国空气动力研究与发展中心

收稿日期: 2024-12-25

  修回日期: 2025-02-23

  网络出版日期: 2025-02-25

基金资助

国家自然科学基金

A fast calculation method of poisson equation source term for BOS measuring density field

  • ZHANG Jian ,
  • ZHANG Zheng-Yu ,
  • YANG Yang ,
  • QIAN Feng-Xue ,
  • LI Xiao-Xia ,
  • WANG Cai-Gang ,
  • LUO Jia-Jie
Expand

Received date: 2024-12-25

  Revised date: 2025-02-23

  Online published: 2025-02-25

Supported by

National Natural Science Foundation of China

摘要

摘 要:背景纹影是非接触测量流场密度的重要手段,通常采用有限差分法求解泊松方程计算投影密度场,但现有方法 计算泊松方程源项(即网格节点上光线偏折角的一阶偏导数)耗时长、突变处精度差。为此,本文基于BOS各测量点的光 线偏折角数据,创建测量点坐标与光线偏折角场间的哈希函数,旨在快速查找以给定点为中心的局部区域BOS测量点集 合及其光线偏折角;建立光线偏折角场中偏折角突变测量点捕捉方法,按突变测量点划分区域,推导并构造区域内偏折 角插值型求导公式,基于哈希表和插值型求导公式,分别计算各区域均匀网格节点上光线偏折角的偏导数。仿真实验表 明:较现有泊松方程源项计算方法,本方法的残差绝对值减小了56.66%、峰值误差减小了75.8%,在密度变化的峰值处 重建效果更好。风洞实验表明:较现有方法,本方法解得的空腔模型密度场更精细,加速比高达3582.72;2米超风洞 “7°(半锥角)的锥柱体模型和某音爆模型”头部微弱激波的测量结果与理论值吻合,激波两侧密度比测量的最大相对误 差为3.9%,流场密度解算正确。因此,本方法提高了密度场测量的速率、精度与微弱激波捕捉能力,工程应用价值大。

本文引用格式

张建 , 张征宇 , 杨洋 , 钱丰学 , 李小霞 , 王材钢 , 罗家杰 . BOS测量密度场的泊松方程源项快速计算方法[J]. 航空学报, 0 : 1 -0 . DOI: 10.7527/S1000-6893.2025.31713

Abstract

Abstract: Background-oriented schlieren (BOS) is an important non-contact method for measuring flow density field. The den sity field is typically calculated by solving the Poisson equation using finite difference methods. However, the current methods for calculating the Poisson equation source term (i.e., the first-order derivatives of ray deflection angles at grid nodes) are time-con suming and low accuracy where deflection angles change hardly. To address these issues, a hash function is created to establish a mapping between point coordinates and the deflection angles by this paper based on ray deflection angles measured by BOS, to find rapidly the BOS measurement point set and their deflection angles within a local region whose center is the given point. A method is developed to capture measurement points where deflection angles change hardly, used to divide the field into regions accordingly, and the interpolation-based derivative formula is derived to calculate the deflection angles’ derivatives within these regions. The simulation experiments demonstrate the proposed method reduces the residuals by 56.66% and peak error by 75.8%, compared to existing methods for calculating the source term of the Poisson equation, have the better results at peaks of density changing. The wind tunnel experiments further demonstrate this method achieves a better refinement density field of a cavity model with significant speed-up computing ratio up to 3582.72, compared to existing approaches. The measuring results of the weak shock wave at the head of the " 7°(semi-apex angle)cone-cylinder model and a sonic boom model" in the 2-meter supersonic wind tunnel are consistent with the theoretical values, and the maximum relative error of the density ratio on both sides of this weak shock wave is 3.9%, and the flow field density reconstruction is correct. Therefore, this method improves the rate, accuracy and weak shock wave capture ability for BOS measuring density field, and has great engineering application value.

参考文献

[1] Chihang Y, Zhang H, Weida F. Pattern control for large-scale spacecraft swarms in elliptie orbits via den-sityfields[J]. Chinese Journal of Aeronautics, 2022, 35(3): 367-379.
[2] Zhu W, Bons J, Gregory J. Boundary layer and near wake measurements of a two-dimensional airfoil with background-oriented schlieren method[J]. Experiments in Fluids, 2023, 64(1): 7.
[3] 吴军, 祝玉恒, 王豪爽, 等. 采用纹影法解耦流速与密度场的高速气流场压力分布重构方法[J]. 光学学报:2023,43(11):71-80.
WU J, ZHU Y H, WANG H S. et al. Reconstruetion of Pressure Distribution in High-Velocity Airflow Fields by Schlieren Decoupling of Velocity and Density Fields[J]. Acta Optica Sinica, 2023, 43(11): 1112004(in Chinese).
[4] ?umnu A, Güzelbey ?. The effects of different wing configurations on missile aerodynamics[J]. Journal of Thermal Engineering, 2023, 9(5): 1260-1271.
[5] 孙佳濛, 左光, 徐艺哲, 等. 临近空间高超声速飞行器武器投放方案数值模拟[J]. 航空学报, 2023,44(13):127808.
SUN J M, ZUO G, XU Y Z, et al. Numerical simulation of weapon delivery schemes for hypersonie vehicles in near space[J]. Acta Aeronautica et Astronautica Sinica, 2023, 44(13): 127808(in Chinese).
[6] Fang C, Hong L I U. Particle image velocimetry for com-bustion measurements: Applications and develop-ments[J]. Chinese Journal of Aeronauties, 2018, 31(7): 1407-1427.
[7] Guo G M, Liu H. Density and temperature reconstruc-tion of a flame-induced distorted flow feld based on back-ground-oriented schlieren (BOS) technique[J]. Chinese Physics B, 2017, 26(6): 064701.
[8] Rückert F U, Lehser-Pfeffermann D, Hübner D, et al. Qualitative visualization of butane and helium tracers in air using a digital background-oriented schlieren method[J]. Fuel, 2023, 342: 127818.
[9] Tziotziou K, Scullion E, et al. Vortex motions in the solar atmosphere: Definitions, theory, observations, and modelling[J]. Space Science Reviews, 2023, 219(1): 1.
[10] 许亮, 王鹿洋, 万自明, 等, 2023. 不同高度对0°~15°攻角气动光学成像偏移的影响[J].红外与激光工程, 52(5): 347-355.
XU L, WANG L Y, WAN Z M, et al. Influence of dif-ferent altitudes on deviation of aero-optics imaging of 0°-15° angle of attack [J]. Infrared and Laser Engineer-ing, 52(5), 347-355. (in Chinese).
[11] Bathel B, Weisberger J, Jones S, et al. Density-Based Optical Diagnostic Techniques at Nasa Langley Re-search Center[C]//AlAA Tennessee Section Monthly Seminar. 2023(20220017668).
[12] WANG Y, et al. Optical signal characteristics analysis of atmospheric disturbance density fields generated by high-speed aircraft, Chin J Aeronaut (2024), https://doi.org/10.1016/j.cja.2024.10.002
[13] 熊渊. 背景纹影测量技术研究与应用进展[J]. 实验流体力学, 2022, 36(02): 30-48.
XIONG Y. Recent advances in background oriented Schlieren and its applications[J]. Journal of Experiments in Fluid Mechanics, 2022, 36(02): 30-48(in Chinese).
[14] Weisberger J M, Bathel B F. Projection background -oriented schlieren[J]. Applied Optics, 2022, 61(20): 6006-6015.
[15] Akatsuka J, Nagai S. Flow visualization by a simplified BOS technique[C]//29th AIAA Applied Aerodynam-icsConference. 2011: 3653.
[16] 张俊,吴运刚,严来军,等. 基于BOS的超声速流场瞬态密度场的可视化[J].气体物理,2021,6(1):62-68.
ZHANG J, WU Y G, YAN L J, et al. Visualization of in-stantaneous density distribution based on BOS for su-personic flow[J]. Physis of Gases, 2021, 6(1): 62-68(in Chinese).
[17] Bron J A, Baars W J, Schrijer F F J. Density Field Reconstruction of an Overexpanded Supersonic Jet us-ing Tomographic Background-Oriented Schlieren[J]. arXiv preprint arXiv:2311.10332, 2023.
[18] Shoaib Amjad, Julio Soria, Callum Atkinson. Three-dimensional density measurements of a heated jet using laser-speckle tomographic background-oriented schlie-ren[J] Experimental Thermal and Fluid Science, 2023, 142: 110819.
[19] Jun Z, Di X, Long Z. Research on density measure-ment based on background oriented schlieren meth-od[J]. Journal of Experiments in Fluid Mechanics, 2015 (1): 77-82.
[20] Fisher T B, Quinn M K, Smith K L . An experimental sensitivity comparison of the schlieren and background oriented schlieren techniques applied to hypersonic flow[J]. Measurement science & technology, 2019, 30(6): 065202.
[21] Meier A H, Roesgen T. Improved background orient ed schlieren imaging using laser speckle illumination[J]. Experiments in fluids, 2013.54: 1-6.
[22] Iwamoto M, Miki Y, Kinefuchi K. Background-oriented schlieren and laser Rayleigh scattering com-plementary method for accurate density field visualiza-tion[J]. Experiments in Fluids, 2024, 65(6): 87.
[23] Akamine M, Teramoto S, Okamoto K. Formulation and demonstrations of three-dimensional background-oriented schlieren using a mirror for near-wall density measurements[J]. Experiments in Fluids, 2023, 64(7): 134.
[24] Gomez M, Grauer S J, Ludwigsen J, et al. Megahertz-rate background-oriented schlieren tomography in postdetonation blasts[J]. Applied optics, 2022, 61(10): 2444-2458.
[25] 张征宇, 王显圣, 黄叙辉, 等. 高速复杂流动结构的视频测量[J]. 航空学报, 2017,38(8):120989.
ZHANG Z Y, WANG X S, HUANG X H, et al.Videogrammetry measurement for high-speed com-plex flow structures[J]. Acta Aeronautica et Astronautica Sinica, 2017, 38(8): 120989(in Chinese).
[26] Bheeroo V, Mandel T L. Comparison of schlieren-based
techniques for measurements of a turbulent and wavy free
surface[J]. Experiments in Fluids, 2023, 64(6): 1-21.
[27] Amjad S, Karami S, Soria J, et al. Assessment of three-di mensional density measurements from tomographic background-oriented schlieren(BOS)[J]. Measurement Science and Technology,2020,31(11):114002.
[28] 周润, 唐亮, 李平, 等. 基于光线偏折的复杂流动度重建方法[J]. 光学学报, 2022, 42(09): 0901001.
ZHOU R, TANG L, LI P, et al. Density Reconstruc-tionMethod for Complex Flow Based on Light Deflec-tion[J]. Acta Optica Sinica, 2022, 42(9): 0901001(in Chinese).
[29] 李庆扬. 数值分析[M]. 第五版. 北京: 清华大学大学出版社, 2008: 25-29.
Li Q Y. Numerical Analysis[M]. 5th ed. Beijing: Tsinghua University Press, 2008: 25-29. (in Chinese).
[30] Thompson, P. A. Compressible-Fluid Dynamics[M]. 1st ed. New York: McGraw-Hill, 1972: 276-303.
[31] Anderson, J. D. Modern Compressible Flow: With His-torical Perspective[M]. 2nd ed. New York: McGraw-Hill, 1990: 294-303.
文章导航

/