Computational methods for acoustic propagation based on linearized equations
-
摘要:
复杂流场如剪切层、旋涡等会改变气动噪声的传播特性,引起折射、反射和散射等现象,对声源识别、测量产生影响。基于线化方程的声传播计算方法是研究声波在复杂流场中传播的重要手段。本文针对非均匀流中声传播计算存在的问题,介绍了近年来课题组开展的研究工作:提出了改进梯度项抑制方法,以抑制基于线化方程模拟声波穿过剪切层时出现的数值不稳定波;发展了适用于线化方程的基于Boltzmann模型的有限体积法通量计算格式,以模拟声波在包含复杂外形流场中的传播问题;发展了简化的线化格子Boltzmann方法,改善了格子Boltzmann方法模拟声传播时内存占用大的问题;研究了剪切层对声源定位的影响规律,建立了定位误差随射流马赫数、斯特劳哈尔数变化的数学模型,为提高实验测量精度提供参考。
-
关键词:
- 线化方程 /
- 声传播 /
- 格子Boltzmann模型 /
- 计算气动声学
Abstract:Complex flows such as shear layers and vortices can change the propagation characteristics of aerodynamic noise, resulting in phenomena of refraction, reflection, and scattering, which have an impact on acoustic source measurement and identification. The linearized equation-based algorithm is an important tool for simulating the propagation of acoustic waves in complex flows. This paper briefly outlines the research carried out by our group in recent years to address the challenges associated with acoustic propagation in non-uniform flow. Specifically, we propose an improved gradient term suppression method to mitigate numerical instability when calculating acoustic wave propagating through shear layers using the linearized equations; develop a Boltzmann model-based flux solver under the finite volume method for the linearized equations to simulate acoustic wave propagation in flow fields containing complex geometries; establish a simplified linearized lattice Boltzmann method to reduce the computational memory requirements when using the lattice Boltzmann method; and evaluate the influences of shear layer characteristics on sound source localization, and establish a mathematical model of localization error with jet Mach number and Strouhal number to provide reference for improving the experimental measurement accuracy.
-
0 引言
随着人们消费水平的提升和消费观念的升级, 汽车NVH (Noise, Vibration, Harshness)性能中的车内噪声水平越来越受到消费者的关注。汽车在高速行驶时, 由外部气流引起的气动噪声是车内噪声的主要来源之一[1]。
在汽车的开发过程中, 由于风洞实验成本高, 数值模拟已成为研究车辆气动噪声的主要手段, 而准确的噪声分析是建立在精确的流场数据基础之上的。因此, 提高流场计算精度十分必要。
大涡模拟(LES)是计算非定常流场的常用方法。该方法利用滤波函数对流场中的涡进行尺度划分, 对大尺度的涡结构进行直接求解, 对小尺度的涡结构则采用亚格子模型进行建模, 而汽车气动噪声的产生与流场中小尺度涡旋的形成与运动有着重要的关系。滤波处理后的N-S方程将产生表征大、小尺度相互作用的亚格子应力项, 对亚格子应力项的求解通过构建亚格子模型进行。不同亚格子模型的选取会获得不同的流场结果, 进而导致不同的声场结果。
1 数值计算理论
1.1 流场数值计算方法
汽车运动的马赫数小于0.3, 其绕流可看作不可压缩流, 不可压N-S方程经过滤波后得到LES控制方程:
(1) (2) 式中, τij为亚格子应力项, 表示可求解的大尺度涡和不可求解的小尺度涡的动量交换:
(3) 为使方程封闭, 需构建亚格子应力模型, 即τij的表达式。采用Smagorinsky-Lilly、Dynamic Smagorinsky、WALE、WMLES、WMLES S-Omega以及Dynamic Kinetic Energy Subgrid-Scale Model 6种亚格子模型进行数值模拟。亚格子模型的τij表达式可写为:
(4) 式中,μt为湍流涡粘度,
为应变率张量:(5) 在Smagorinsky-Lilly模型[2]中, 湍流涡粘度:
(6) 其中,
(7) 式中, k为冯卡尔曼常数; d为近壁面的距离; Cs为Smagorinsky常数, 通常取0.1;Δ为亚网格尺寸, 由网格尺寸V确定, Δ=V1/3; Ls为亚格子尺度的混合长度。
动态Smagorinsky(Dynamic Smagorinsky)亚格子尺度模型[3]中, Cs基于运动尺度求解得出。该模型添加了网格滤波和检验滤波函数, 进一步得出湍流应力张量:
(8) 式中, Tij为滤波后的亚格子应力项。
通过Lilly最小二乘分析, 最终得到Cs表达式为:
(9) 式中:
(10) 局部涡粘度壁面自适应模型[4](Wall-Adapting Local Eddy-Viscosity Model, WALE)中湍流涡粘度为:
(11) 式中, 亚格子尺度的混合长度Ls和应变率张量Sijd的表达式为:
(12) (13) 式中, 常数Cw=0.325,
表达式为:(14) 代数形式的壁面函数模型[5](Algebraic Wall-Modeled LES Model, WMLES)在边界层内部使用雷诺平均方法进行模拟, 在边界层外部使用LES进行模拟。WMLES可应用于更高雷诺数的流场模拟, 其湍流涡粘度:
(15) 式中, dw为计算点与壁面的距离; S为应变率; 常数k=0.41, CSmag=0.2;Δ为亚网格尺度:
(16) 式中, hmax是六面体单元的最大边长, hwn是壁面法线方向网格间距, 常数Cw=0.15。
由于无法为恒定的剪切流提供零涡流粘度, WMLES模型在分离剪切层时会产生过大的涡粘性。在式(15)中, 若将S替换为|S-Ω|的差分形式, 则可以改进WMLES模型, 这种改进模型称为WMLES S-Omega模型, 其中Ω为旋涡强度。
动态动能亚格子模型[6](Dynamic Kinetic Energy Subgrid-Scale Model, KET)在传统Smagorinsky-Lilly模型上, 进一步考虑亚格子尺度湍动能的传输, 能更好地模拟亚格子尺度的湍流。亚格子湍流动能ksgs为:
(17) 湍流涡粘度μt通过亚格子湍流动能计算得出:
(18) 式中滤波尺寸Δf=V1/3, Ck根据流场状态获得。该模型的亚格子应力表达式为:
(19) 1.2 声学数值计算方法
采用莱特希尔声类比理论(Lighthill’s acoustic analogy)对气动声学进行求解。从流体力学的连续方程和动量方程推导得出莱特希尔方程:
(20) 其中:
(21) 式中, ρ'为流体密度波动量, c0为声音在均质介质中的速度, Tij为莱特希尔张量, ui为速度矢量沿i方向的分量, pi、ρ为流体的压力和密度, p0、ρ0分别为参考压力和密度,
为向量微分算子, δij为克罗内克函数。应用格林函数进行积分可对莱特希尔声类比方程求解, 但对于复杂边界, 格林函数往往难以获得。Oberai等人提出一种新的求解方法[7], 通过获得频域下的莱特希尔方程变分形式, 采用有限元法进行求解。
如图 1所示, 气流流经的固体区域为ΩS, 固体区域的边界为ΓS。流体区域为Ωf, 流体是非均匀的或者湍流的, 受N-S方程控制。在Ωa和Ω∞区域, 流体状态是均匀的。另外, 湍流区域(声源区)的外边界为Γf, 声学截断面Γa为无反射边界。n为区域Ωt=Ωf∪Ωa向外的法向向量。对于任意速度矢量u, 可用un= u · n =uini来表示其法向分量。
基于莱特希尔方程, 利用转换等式:
(22) 得到莱特希尔方程的频域形式:
(23) 式中, ψ=-b/iω, b为总熵, i(非下标)为虚数单位。
基于Oberai提出的变分公式, 可得:
(24) 式中, δψ为测试函数。利用格林函数进行变分转换, 同时结合动量方程, 最终得到弱解形式:
(25) 采用有限元法进行内声场求解。在一个由边界面Ω所包围、声源分布为q的流体域V中, 任意空间点(x, y, z)的声压p受二阶赫姆霍兹方程控制:
(26) 式中, k=ω/c=2πf/c为波数, c为声速, ρ为流体密度, ω为圆频率。
将式(26)转换为等效的加权残差积分方程:
(27) 式中,
为加权函数。根据发散定理, 由Ω面所包围的区域V内的发散向量场的积分等价于封闭面Ω上向量场φ的垂直分量的积分。由此得到式(26)的弱解为:(28) 式中nV为区域V向外的法向向量。
在有限元方法中, 连续体V离散成一系列有限小单元以及节点ne, 在每个单元内, 利用形函数Nie, 声压p可用一定数量np的
来近似表达。对于常用的四面体或六面体单元, 单元的形函数个数与单元节点的个数相等(np=ne), 且一一对应。同时, 构造一个基于整个区域V的全局形函数Ni。全局压力 可表示为:(29) 式中, nf为整个离散域的节点总数, N为全局形函数的(1×nf)向量,
为未知的节点声压值的(nf×1)向量。相应的压力梯度近似变为:(30) 式中, B为全局形函数梯度分量的(3×nf)矩阵。将式(29)和(30)代入到式(28)中, 得到:
(31) 式中, K为声学刚度矩阵, M为声学质量矩阵, Q i为声学激励向量, C为声学阻尼矩阵。未知节点的声压值
可由式(31)求解得到。1.3 外部气流引发车内噪声的机理
在密封良好情况下, 泄露噪声、风振噪声的影响可以忽略。此时车内气动噪声可归结为以下两种机理[8]:(1)空气中的旋涡收缩与变形及其相互作用, 以及气流与壁面相互作用会产生气动噪声源, 多种气动噪声源向四周传递噪声, 其中到达玻璃壁面的声压会引起壁面的振动从而向模型内部传递噪声; (2)当气流流经玻璃表面时, 引起了壁面振动, 从而向模型内部传递噪声。
虽然作用在玻璃表面上的湍流压力的压力级数要比声压的压力级数大, 但是在传递效率上, 湍流压力的传递损失却比声压的传递损失要大。所以, 在研究外部流场诱发汽车内部气动噪声时, 需要同时考虑上述2种作用机理。
2 计算模型与设置
2.1 计算模型
采用的计算模型为现代汽车公司建立的简易汽车模型Hyundai Simplified Model(HSM), 参考相关文献[9-10]建立HSM模型。如图 2所示, 模型长2 m, 宽1 m, 高1 m, 壁面从外到内分别为12 mm铝、40 mm吸声材料、4 mm隔膜以及20 mm吸声材料, 玻璃均为4 mm厚。材料属性如表 1和2所示。
表 1 材料属性Table 1 The material properties材料 密度/(kg·m-3) 杨氏模量/GPa 泊松比 铝 2700 69 0.33 隔膜 2000 0.04 0.45 车窗玻璃 2500 70 0.22 空气 1.2 - - 表 2 吸声材料属性Table 2 Properties of sound absorbing materials密度/(kg·m-3) 声速/(m·s-1) 孔隙率 孔曲率 827 343 0.879 3.31 实验中, 通过垫片、密封剂、螺栓等将车窗玻璃安装在铝制外壳上。为了验证车窗玻璃边界条件对噪声传播和数值模拟的影响, 文献[11]对车窗玻璃表面内外的频率响应函数(FRF)、玻璃传递函数和车内混响时间进行了测量, 结果表明, 玻璃内外的FRF结果几乎相同。可以认为, 玻璃安装到铝制外壳上是两种材料的简单附着, 仿真模型的简化对计算结果的影响可以忽略。
2.2 流场计算
参考现代汽车公司风洞实验段的外形[12]建立的流场计算域如图 3所示, 计算域长20 m, 宽14 m, 高10 m, 模型距离入口4 m。模型表面网格尺寸为5 mm, 计算域外表面网格尺寸为500 mm。利用密度盒对模型附近网格进行加密, 2个加密区的最大网格尺寸为20和100 mm, 其中, 第一个密度盒包含了湍流较为剧烈的声源区。在模型表面划分5层边界层网格, 第一层厚0.2 mm, 增长率为1.2。整个计算域网格数量约为1500万, 边界条件和参数设置见表 3。
表 3 边界条件设置Table 3 Boundary condition setting边界 类型 数值 Inlet Velocity inlet 110 km/h Outlet Pressure outlet 0 Pa Others Wall - 采用Fluent进行流场计算。利用稳态结果作为瞬态的初始值进行计算, 以减少瞬态计算达到收敛的时间。分别使用上述6种亚格子模型进行瞬态计算, 以获得流场的速度和压力结果。考虑到气动噪声集中在中低频段, 本次计算的最高频率初步定为1000 Hz。根据奈奎斯特采样定理, 确定采样时长为0.08 s, 计算时间步长为2×10-4 s。
2.3 声场计算
采用Actran进行声学计算, 使用积分插值的方法将流场信息映射至声学网格。外声场模型如图 4所示, 分为声源区、声传播区和无限元边界, 无限元边界为无反射边界。实验风洞为声学风洞, 壁面进行了吸声处理, 采用无反射边界和实验保持一致。根据分析, 频率的波长内至少包含6个网格单位的要求, 确定声场的网格尺寸为40 mm。采用四面体网格对外声场进行划分, 最终网格数约为140万。
图 5为内声场模型局部网格。玻璃和隔膜使用三棱柱网格划分, 其余部分使用四面体网格划分, 网格尺寸和外声场相同。
为了和实验数据[9, 11]对比, 在模型内部建立2个监测点, 监测点位置如图 6所示。
3 结果分析
3.1 流场结果
由流场决定的CFD截止频率fc为:
(32) 式中, ε为湍动能耗散率, Δx为网格尺度。
图 7为CFD截止频率分布图, 由图可知在湍流核心区域, CFD截止频率大部分在1000 Hz左右, 符合计算要求。
3.2 声场结果
图 8和9为1/3倍频程频率中心点的声压级。总体来看, 声压级Lp随频率f的上升而逐渐降低, 不同亚格子模型计算得到的结果的变化趋势和实验一致, 但在数值上有偏差。上述2个监测点的总声压级以及和实验值的误差如表 4和5所示。
表 4 监测点1总声压级及其误差Table 4 Overall sound pressure level and error of monitoring point 1SLM DSLM WALE WMLES WMLES S-Omega KET 实验值 总声压级/dB 90.1 90.8 80.4 77.2 80.2 85.9 71.6 误差/% 25.8 26.8 12.3 7.8 12.0 19.9 - 表 5 监测点2总声压级及其误差Table 5 Overall sound pressure level and error of monitoring point 2SLM DSLM WALE WMLES WMLES S-Omega KET 实验值 总声压级/dB 90.3 90.4 83.7 77.5 76.1 86.2 73.1 误差/% 23.5 23.7 14.5 6.0 4.1 17.9 - 3.3 结果分析
从总声压级看, 只有WMLES模型的结果误差在10%以内, 其余模型误差都较大。所以采用WMLES模型计算总声压级的精度更高。
由频谱曲线可知, 100~160 Hz, 6种亚格子模型得到的计算结果均大于实验值; 160~500 Hz, 计算结果均小于实验值; 500~1000 Hz, 计算结果和实验结果基本一致。在200 Hz时, 实验曲线出现峰值, 原因可能是该频率点为模型的固有频段, 共振使得声压级增大, WMLES和WMLES S-Omega模型准确地捕捉到了这一信息。
声场结果和流场结果有着直接的联系, 高频段结果和流场中小旋涡的变化及相互作用有关。出现低频段计算结果高于实验值、高频段计算结果低于实验值的原因可能是受计算资源的限制, 数值模拟没能精确捕捉到流场中的小旋涡运动, 导致部分高频段能量计入低频段。
在近壁面附近, 湍流尺度随壁面距离的减小而线性减小。所以, 壁面附近小尺度的旋涡占主要部分。由1.3节可知, 壁面附近的流场对车内噪声有着显著的影响。壁面附近气流模拟的精确性将直接影响车内噪声的计算结果。
LES在处理近壁面的流场时需要非常精细的网格尺度和非常小的时间步长(y+=1)。WMLES模型对LES进行了改进, 在处理近壁面流场时使用RANS进行模拟, 其余部分使用LES进行模拟。RANS在求解近壁面流场时, 不需要LES所要求的非常精细的网格尺度, 使其在求解高雷诺数下近壁面湍流更加有效。
在LES中, 亚格子模型表示了大、小旋涡的相互作用。几种亚格子模型的区别在于对湍流涡粘度μt的建模方式不同。湍流模拟中, μt对小尺度旋涡的形成与运动有很大的影响。在对μt建模过程中, 和网格尺寸有关的亚网格尺度Δ是一个重要的参数。其余亚格子模型采用网格体积的1/3次方定义Δ, WMLES模型改变了Δ的定义, 综合考虑了网格边长、网格中心点离壁面距离和网格法向尺寸, 此外建模时还考虑了y+值, 使得WMLES模型对壁面附近的旋涡捕捉更为精确, 而其余模型计算得到的湍流衰减过大, 这一点在文献[5]中已有体现。200 Hz时, WMLES及其改进模型WMLES S-Omega能够捕捉到峰值, 和其能够准确捕捉流场信息有关。
4 结论
基于不同亚格子模型, 对HSM模型进行了流场和声场的计算, 并与实验结果进行了对比分析, 结果表明:采用相同网格尺度时, 基于WMLES亚格子模型的LES数值模拟能较为准确地获得车外气动声源。
-
-
[1] 张树海, 李虎, 王益民. 含激波、旋涡和声波的复杂多尺度流动数值模拟研究[J]. 空气动力学学报, 2018, 36(3): 449–462. ZHANG S H, LI H, WANG Y M. Numerical study for complex multi-scale flow with shock, vortex and sound wave[J]. Acta Aerodynamica Sinica, 2018, 36(3): 449–462.
[2] COLONIUS T, LELE S K, MOIN P. The scattering of sound waves by a vortex: numerical simulations and analytical solutions[J]. Journal of Fluid Mechanics, 1994, 260: 271–298. doi: 10.1017/s0022112094003514
[3] HONG Z L, WANG X Y, JING X D, et al. Vortex sound interaction in acoustic resonance of a flow duct containing a plate[J]. Journal of Sound and Vibration, 2020: 115482. doi: 10.1016/j.jsv.2020.115482
[4] CLAIR V, GABARD G. Spectral broadening of acoustic waves by convected vortices[J]. Journal of Fluid Mechanics, 2018, 841: 50–80. doi: 10.1017/jfm.2018.94
[5] DU Y L, YU H W, LIU Y C, et al. Hydrodynamic sources of the vortex sound in a two-dimensional shear layer[J]. International Journal of Aeroacoustics, 2023, 22: 41–59. doi: 10.1177/1475472X221150177
[6] MA R X, ZHANG S H, LUO Y, et al. Numerical study on multiple acoustic scattering by a vortex array[J]. Journal of Sound and Vibration, 2022, 527: 116815. doi: 10.1016/j.jsv.2022.116815
[7] BAILLY C, JUVÉ D. Numerical solution of acoustic propagation problems using linearized Euler equations[J]. AIAA Journal, 2000, 38(1): 22–29. doi: 10.2514/2.949
[8] 李晓东, 江旻, 高军辉, 等. 计算气动声学进展与展望[J]. 中国科学: 物理学 力学 天文学, 2014, 44(3): 234-248. LI X D, JIANG M, GAO J H, et al. Progress and prospective of computational aeroacoustics[J]. Scientia Sinica (Physica, Mechanica & Astronomica), 2014, 44(3): 234-248.
[9] 吕宏强, 朱国祥, 宋江勇, 等. 线化欧拉方程的高阶间断有限元数值解法研究[J]. 力学学报, 2011, 43(3): 621–624. DOI: 10.6052/0459-1879-2011-3-lxxb2010-077 LYU H Q, ZHU G X, SONG J Y, et al. High-order discontinuous Galerkin solution of linearized Euler equations[J]. Chinese Journal of Theoretical and Applied Mechanics, 2011, 43(3): 621–624. doi: 10.6052/0459-1879-2011-3-lxxb2010-077
[10] 余培汛, 白俊强, 杨海, 等. 耦合RANS/LES模型与LEE方程的气动噪声混合预测方法[J]. 西北工业大学学报, 2017, 35(6): 990–997. DOI: 10.3969/j.issn.1000-2758.2017.06.009 YU P X, BAI J Q, YANG H, et al. Hybrid aero-acoustics prediction method coupled with RANS/LES model and LEE equation[J]. Journal of Northwestern Polytechnical University, 2017, 35(6): 990–997. doi: 10.3969/j.issn.1000-2758.2017.06.009
[11] XIN B, SUN D K, JING X D, et al. Numerical study of acoustic instability in a partly lined flow duct using the full linearized Navier–Stokes equations[J]. Journal of Sound and Vibration, 2016, 373: 132–146. doi: 10.1016/j.jsv.2016.02.042
[12] BOGEY C, BAILLY C, JUVÉ D. Computation of flow noise using source terms in linearized Euler’s equations[J]. AIAA Journal, 2002, 40(2): 235–243. doi: 10.2514/2.1665
[13] ZHANG X, CHEN X X, GILL J R. Gradient term filtering for stable sound propagation with linearized Euler equations[C]//Proc of the 20th AIAA/CEAS Aeroacoustics Conference. 2014: 3306. doi: 10.2514/6.2014-3306
[14] EWERT R, SCHRÖDER W. Acoustic perturbation equations based on flow decomposition via source filtering[J]. Journal of Computational Physics, 2003, 188(2): 365–398. doi: 10.1016/S0021-9991(03)00168-2
[15] BAUER M, DIERKE J, EWERT R. Application of a discontinuous Galerkin method to discretize acoustic perturbation equations[J]. AIAA Journal, 2011, 49(5): 898–908. doi: 10.2514/1.J050333
[16] EWERT R, SCHRÖDER W. On the simulation of trailing edge noise with a hybrid LES/APE method[J]. Journal of Sound and Vibration, 2004, 270(3): 509–524. doi: 10.1016/j.jsv.2003.09.047
[17] SEO J H, MOON Y J. Linearized perturbed compressible equations for low Mach number aeroacoustics[J]. Journal of Computational Physics, 2006, 218(2): 702–719. doi: 10.1016/j.jcp.2006.03.003
[18] HU F Q, LI X D, LI X Y, et al. Time domain wave packet method and suppression of instability waves in aeroacoustic computations[J]. Journal of Fluids Engineering, 2014, 136(6): 060905. doi: 10.1115/1.4025866
[19] 李旦望, 李晓东, 李小艳, 等. 管道后传声数值模拟的不稳定波抑制[J]. 工程热物理学报, 2012, 33(4): 587–590. LI D W, LI X D, LI X Y, et al. Suppression of instability waves in numerical simulations of sound propagation from aft ducts[J]. Journal of Engineering Thermophysics, 2012, 33(4): 587–590.
[20] PRAX C, GOLANSKI F, NADAL L. Control of the vorticity mode in the linearized Euler equations for hybrid aeroacoustic prediction[J]. Journal of Computational Physics, 2008, 227(12): 6044–6057. doi: 10.1016/j.jcp.2008.02.022
[21] TESTER B, GABARD G, ÖZYÖRÜK Y. Influence of mean flow gradients on fan exhaust noise predictions[C]//Proc of the 14th AIAA/CEAS Aeroacoustics Conference (29th AIAA Aeroacoustics Conference). 2008: 2825. doi: 10.2514/6.2008-2825
[22] DENG Y Y, ALOMAR A, DRAGNA D, et al. Characterization and suppression of the hydrodynamic instability in the time domain for acoustic propagation in a lined flow duct[J]. Journal of Sound and Vibration, 2021, 500: 115999. doi: 10.1016/j.jsv.2021.115999
[23] XU K. A gas-kinetic BGK scheme for the Navier-Stokes equations and its connection with artificial dissipation and Godunov method[J]. Journal of Computational Physics, 2001, 171(1): 289–335. doi: 10.1006/jcph.2001.6790
[24] Guo Z L, Wang L P, Qi Y M. Discrete unified gas kinetic scheme for continuum compressible flows[J]. Physical Review E, 2023, 107(2): 025304. doi: 10.1103/PhysRevE.107.025304
[25] Zhong M L, Zou S, Pan D X, et al. A simplified discrete unified gas-kinetic scheme for compressible flow[J]. Physics of Fluids, 2021, 33(3): 036103. doi: 10.1063/5.0033911
[26] SHU C, WANG Y, TEO C J, et al. Development of lattice boltzmann flux solver for simulation of incompressible flows[J]. Advances in Applied Mathematics and Mechanics, 2014, 6(4): 436–460. doi: 10.4208/aamm.2014.4.s2
[27] WANG Y, YANG L M, SHU C. From lattice boltzmann method to lattice boltzmann flux solver[J]. Entropy, 2015, 17(12): 7713–7735. doi: 10.3390/e17117713
[28] 杨鲤铭. 基于格子和连续Boltzmann模型的通量求解器研究及其应用[D]. 南京: 南京航空航天大学, 2016. YANG L M. Research and application of flux solver based on lattice and continuous Boltzmann model[D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2016.
[29] ZHU Y J, ZHONG C W, XU K. GKS and UGKS for high-speed flows[J]. Aerospace, 2021, 8(5): 141. doi: 10.3390/aerospace8050141
[30] PAN D X, ZHONG C W, ZHUO C S. An implicit gas-kinetic scheme for turbulent flow on unstructured hybrid mesh[J]. Computers & Mathematics with Applications, 2018, 75(11): 3825–3848. doi: 10.1016/j.camwa.2018.02.032
[31] LI J, ZHONG C W, WANG Y, et al. Implementation of dual time-stepping strategy of the gas-kinetic scheme for unsteady flow simulations[J]. Physical Review E, 2017, 95(5): 053307. doi: 10.1103/physreve.95.053307
[32] ZHANG C, LI Q B, FU S, et al. A third-order gas-kinetic CPR method for the Euler and Navier–Stokes equations on triangular meshes[J]. Journal of Computational Physics, 2018, 363: 329–353. doi: 10.1016/j.jcp.2018.02.040
[33] CHEN S Z, GUO Z L, XU K. Simplification of the unified gas kinetic scheme[J]. Physical Review E, 2016, 94(2): 023313. doi: 10.1103/physreve.94.023313
[34] YANG L M, SHU C, WU J. Extension of lattice Boltzmann flux solver for simulation of 3D viscous compressible flows[J]. Computers & Mathematics with Applications, 2016, 71(10): 2069–2081. doi: 10.1016/j.camwa.2016.03.027
[35] YANG Y Q, PAN L, XU K. High-order gas-kinetic scheme on three-dimensional unstructured meshes for compressible flows[J]. Physics of Fluids, 2021, 33(9): 09610. doi: 10.1063/5.0062368
[36] ZHANG C, LI Q B. A third-order subcell finite volume gas-kinetic scheme for the Euler and Navier-Stokes equations on triangular meshes[J]. Journal of Computational Physics, 2021, 436: 110245. doi: 10.1016/j.jcp.2021.110245
[37] ZHAN N Y, CHEN R Q, YOU Y C. Meshfree method based on discrete gas-kinetic scheme to simulate incompressible/compressible flows[J]. Physics of Fluids, 2021, 33(1): 017112. doi: 10.1063/5.0033770
[38] LIU Y Y, SHU C, ZHANG H W, et al. A high order least square-based finite difference-finite volume method with lattice Boltzmann flux solver for simulation of incompressible flows on unstructured grids[J]. Journal of Computational Physics, 2020, 401: 109019. doi: 10.1016/j.jcp.2019.109019
[39] MARIÉ S, RICOT D, SAGAUT P. Comparison between lattice Boltzmann method and Navier-Stokes high order schemes for computational aeroacoustics[J]. Journal of Computational Physics, 2009, 228(4): 1056–1070. doi: 10.1016/j.jcp.2008.10.021
[40] LATT J, MALASPINAS O, KONTAXAKIS D, et al. Palabos: parallel lattice boltzmann solver[J]. Computers & Mathematics with Applications, 2021, 81: 334–350. doi: 10.1016/j.camwa.2020.03.022
[41] CASALINO D, HAZIR A, MANN A. Turbofan broadband noise prediction using the lattice boltzmann method[J]. AIAA Journal, 2017, 56(2): 609–628. doi: 10.2514/1.J055674
[42] SCHÄFER R, BÖHLE M. Validation of the lattice boltzmann method for simulation of aerodynamics and aeroacoustics in a centrifugal fan[J]. Acoustics, 2020, 2(4): 735–752. doi: 10.3390/acoustics2040040
[43] 司海青, 朱卫军. 气动噪声计算方法及其应用[M]. 北京: 科学出版社, 2017. SI H Q, ZHU W J. Calculation method of aerodynamic noise and its application[M]. Beijing: Science Press, 2017.
[44] VERGNAULT E, MALASPINAS O, SAGAUT P. A lattice Boltzmann method for nonlinear disturbances around an arbitrary base flow[J]. Journal of Computational Physics, 2012, 231(24): 8070–8082. doi: 10.1016/j.jcp.2012.07.021
[45] PÉREZ J M, AGUILAR A, THEOFILIS V. Lattice Boltzmann methods for global linear instability analysis[J]. Theoretical and Computational Fluid Dynamics, 2017, 31(5): 643–664. doi: 10.1007/s00162-016-0416-7
[46] CHEN Z, SHU C, WANG Y, et al. A simplified lattice boltzmann method without evolution of distribution function[J]. Advances in Applied Mathematics and Mechanics, 2017, 9(1): 1–22. doi: 10.4208/aamm.oa-2016-0029
[47] VIGGEN E M. The lattice Boltzmann method with applications in acoustics[D]. Trondheim: Norwegian University of Science and Technology, 2009.
[48] 王李璨. 剪切流中声波传播问题的计算方法改进与影响规律研究[D]. 厦门: 厦门大学, 2021. WANG L C. Computational method improvement and influential assessment about sound propagation problems in A shear layer[D]. Xiamen: Xiamen University, 2021.
[49] WANG L C, CHEN R Q, YOU Y C, et al. Investigation of acoustic propagation and source localization in a hot jet flow[J]. Journal of Sound and Vibration, 2021, 492: 115801. doi: 10.1016/j.jsv.2020.115801
[50] ZHAN N Y, CHEN R Q, YOU Y C. Linear lattice Boltzmann flux solver for simulating acoustic propagation[J]. Computers & Mathematics with Applications, 2022, 114: 21–40. doi: 10.1016/j.camwa.2022.03.034
[51] ZHAN N Y, CHEN R Q, SONG Q C, et al. Linear discrete velocity model-based lattice Boltzmann flux solver for simulating acoustic propagation in fluids[J]. Physical Review E, 2022, 105(6): 065303. doi: 10.1103/physreve.105.065303
[52] QIAN Y H, D’HUMIÈRES D, LALLEMAND P. Lattice BGK models for Navier-Stokes equation[J]. Europhysics Letters (EPL), 1992, 17(6): 479–484. doi: 10.1209/0295-5075/17/6/001
[53] QU K, SHU C, CHEW Y T. Simulation of shock-wave propagation with finite volume lattice boltzmann method[J]. International Journal of Modern Physics C, 2007, 18(4): 447–454. doi: 10.1142/S012918310701067X
[54] CHEN Z, SHU C, TAN D. Highly accurate simplified lattice Boltzmann method[J]. Physics of Fluids, 2018, 30(10): 103605. doi: 10.1063/1.5050185
[55] CHEN Z, SHU C, TAN D. Immersed boundary-simplified lattice Boltzmann method for incompressible viscous flows[J]. Physics of Fluids, 2018, 30(5): 053601. doi: 10.1063/1.5028353
[56] YANG L M, SHU C, CHEN Z, et al. An improved multiphase lattice Boltzmann flux solver for the simulation of incompressible flow with large density ratio and complex interface[J]. Physics of Fluids, 2021, 33(3): 033306. doi: 10.1063/5.0038617
[57] LU J H, DAI C S, YU P. Analysis and reconstruction of the simplified thermal lattice Boltzmann method[J]. International Journal of Heat and Mass Transfer, 2022, 187: 122576. doi: 10.1016/j.ijheatmasstransfer.2022.122576
[58] SONG Q C, CHEN R Q, CAO S Q, et al. A simplified linearized lattice boltzmann method for acoustic propagation simulation[J]. Entropy, 2022, 24(11): 1622. doi: 10.3390/e24111622
[59] 张军, 王勋年, 张俊龙, 等. 开口风洞声阵列测量的剪切层修正方法[J]. 实验流体力学, 2018, 32(4): 39–46. DOI: 10.11729/syltlx20180013 ZHANG J, WANG X N, ZHANG J L, et al. Shear layer correction methods for open-jet wind tunnel phased array test[J]. Journal of Experiments in Fluid Mechanics, 2018, 32(4): 39–46. doi: 10.11729/syltlx20180013
[60] 张雪, 陈宝, 卢清华. Amiet剪切层理论的角度折射验证研究[J]. 应用声学, 2014, 33(5): 433–438. DOI: 10.11684/j.issn.1000-310X.2014.05.009 ZHANG X, CHEN B, LU Q H. Verification of angle refraction correction based on Amiet’s shear layer theory[J]. Journal of Applied Acoustics, 2014, 33(5): 433–438. doi: 10.11684/j.issn.1000-310X.2014.05.009
[61] 王李璨, 陈荣钱, 尤延铖, 等. 剪切层形态对声波传播与声源定位的影响研究[J]. 西北工业大学学报, 2019, 37(6): 1148–1157. DOI: 10.3969/j.issn.1000-2758.2019.06.008 WANG L C, CHEN R Q, YOU Y C, et al. Effects of shear layer characteristics on acoustic propagation and source localization[J]. Journal of Northwestern Polytechnical University, 2019, 37(6): 1148–1157. doi: 10.3969/j.issn.1000-2758.2019.06.008
-
期刊类型引用(3)
1. 许建林,公沛霖,张春博,贾永兴,李睿,杨永刚. 不同边界条件下单圆柱气动噪声大涡模拟研究. 空气动力学学报. 2024(12): 78-87 . 百度学术
2. 吴鹏,张华. 高速矩形水舌表面破碎的大涡模拟. 水利技术监督. 2022(04): 106-112 . 百度学术
3. 姚远,赵文彬,卢武,党朋,魏建华. 多风速段下低风阻导线抗风能力实验分析及导线截面结构参数优化. 电力系统保护与控制. 2021(21): 97-106 . 百度学术
其他类型引用(14)