指定流量分配系数的多回流出口边界算法

doi:10.7527/S1000-6893.2022.26830

流体力学与飞行力学

本期目录 | 过刊浏览 | 高级检索

前一篇 | 后一篇

指定流量分配系数的多回流出口边界算法

许开龙¹^,²(), 刘再刚¹^,², 姜胜利¹^,², 王星¹^,², 张磐¹^,²

^1.中物院高性能数值模拟软件中心，北京 100088
^2.北京应用物理与计算数学研究所，北京 100088

收稿日期:2021-12-17 修回日期:2022-01-10 接受日期:2022-01-26 出版日期:2023-03-15 发布日期:2022-02-17
通讯作者: 许开龙 E-mail:xu_kailong@iapcm.ac.cn
基金资助:
国家自然科学基金(12102060)

Treatment of boundary condition at multiple outlets with recirculating flow and specified flow ratios

Kailong XU¹^,²(), Zaigang LIU¹^,², Shengli JIANG¹^,², Xing WANG¹^,², Pan ZHANG¹^,²

^1.CAEP Software Center for High Performance Numerical Simulation，Beijing 100088，China
^2.Institute of Applied Physics and Computational Mathematics，Beijing 100088，China

Received:2021-12-17 Revised:2022-01-10 Accepted:2022-01-26 Online:2023-03-15 Published:2022-02-17
Contact: Kailong XU E-mail:xu_kailong@iapcm.ac.cn
Supported by:
National Natural Science Foundation of China(12102060)

摘要/Abstract

摘要：

燃烧室出口处物理量的边界条件处理方法需满足工程设计对出口流量的要求，同时要保证出口回流时的数值稳定性，是工程级燃烧室模拟中重点关注的问题之一。基于流量分配系数发展了一种出口边界处理算法，支持不可压缩流动模拟中指定多个出口流量分配系数、出口存在回流时的稳定模拟。研究表明，在四面体/六面体网格算例中，该方法可实现出口流量分配系数在0~1.0范围内的稳定模拟，与设定值的最大相对偏差<0.001%。在保持出口流量分配系数不变时，改变计算域出口位置使得流动处于未充分发展、回流或旋流状态，计算结果不受明显影响。将该方法应用于TECFLAM旋流燃烧器冷态流场模拟中，计算域缩减至1/3后模拟仍能稳定收敛，预测值和实验数据保持较高的一致性。相比于使用常规出口边界的方法，该方法支持仿真人员根据需要设置出口流量分配系数和计算域大小，可显著减小数值计算开销。

关键词: 自由出口边界, 流量分配系数, 回流, 燃烧室模拟, 不可压缩流动

Abstract:

The treatment of boundary condition at multiple outlets is one of the key issues in the numerical simulation of aero-engine combustors， as the mass flow ratios need to be within the engineering design limits， and the numerical stability needs to be ensured with recirculating flow at the outlets. In the present paper， an incompressible flow outlet boundary condition algorithm is developed to treat the boundary condition at multiple outlets with recirculating flow and specified flow ratios. The results of a series of numerical experiments show that stable and convergent solutions can be obtained using the proposed method. Cases of tetrahedral/hexahedral meshes and the mass flow ratio ranging from 0 to 1.0 are simulated， and the obtained maximum deviation is <0.001%. When the mass flow ratios at the multiple outlets are fixed， the resulted velocity profiles are not obviously affected by the location of outlet boundaries when the flows at the boundaries are not fully-developed or even recirculating. The proposed method is applied to the simulation of TECFLAM swirl burner. The predicted velocity components agree well with the experimental data even when the computational domain is truncated to one third of the original size. Compared to the conventional treatment of outlet boundary conditions， the method proposed enables to specify the mass flow ratios at multiple outlets and to adjust computational domain sizes according to the engineering design limits， which contributes to significant reduction of the overall computational cost.

Key words: free outlet boundary, mass flow ratio, recirculating flow, combustor simulation, incompressible flow

中图分类号:

V231.3

许开龙, 刘再刚, 姜胜利, 王星, 张磐. 指定流量分配系数的多回流出口边界算法[J]. 航空学报, 2023, 44(5): 126830-126830.

Kailong XU, Zaigang LIU, Shengli JIANG, Xing WANG, Pan ZHANG. Treatment of boundary condition at multiple outlets with recirculating flow and specified flow ratios[J]. ACTA AERONAUTICAET ASTRONAUTICA SINICA, 2023, 44(5): 126830-126830.

图/表 16

图 1

图 2

图 3

表 1

90°T形分叉管算例流量统计结果

$α 1 ∶ α 2$	$0 ∶ 1.0$	$0.2 ∶ 0.8$	$0.4 ∶ 0.6$	$0.5 ∶ 0.5$	$0.6 ∶ 0.4$	$0.8 ∶ 0.2$	$1.0 ∶ 0$
入口流量/（kg∙s^-1）	0.667 805	0.667 805	0.667 805	0.667 805	0.667 805	0.667 805	0.667 804
出口1流量/（kg∙s^-1）	-0.000 005	0.133 561	0.267 119	0.333 895	0.400 683	0.534 243	0.667 805
出口2流量/（kg∙s^-1）	0.667 810	0.534 244	0.400 686	0.333 909	0.267 121	0.133 562	0.000 001
计算 $α 1$	-0.000 008	0.200 000	0.399 995	0.499 989	0.600 000	0.799 999	0.999 999
计算 $α 2$	1.000 008	0.800 000	0.600 005	0.500 010	0.399 999	0.200 001	0.000 001
流量误差/%	0.000 76	0.000 04	0.000 46	0.001 03	0.000 08	0.000 08	0.000 09

表 1

图 4

图 5

表 2

60°Y形分叉管算例流量统计结果

$α 1$ ∶ $α 2$	0∶1.0 图 5　60°Y形分叉管算例不同流量分配下速度分布 Fig. 5　Velocity distribution of duct case with 60° Y-branch at various mass flow ratios	0.1∶0.9 图 5　60°Y形分叉管算例不同流量分配下速度分布 Fig. 5　Velocity distribution of duct case with 60° Y-branch at various mass flow ratios	0.2∶0.8	0.3∶0.7	0.4∶0.6	0.5∶0.5
入口流量/（kg∙s^-1）	9.238 5×10^-4	9.238 5×10^-4	9.238 5×10^-4	9.238 5×10^-4	9.238 5×10^-4	9.238 5×10^-4
出口1流量/（kg∙s^-1）	-6.301 5×10^-9	9.238 4×10^-5	1.847 7×10^-4	2.771 6×10^-4	3.695 4×10^-4	4.619 3×10^-4
出口2流量/（kg∙s^-1）	9.238 6×10^-4	8.314 7×10^-4	7.390 8×10^-4	6.467 0×10^-4	5.543 1×10^-4	4.619 3×10^-4
计算 $α 1$	-0.000 007	0.099 999	0.200 000	0.300 001	0.400 000	0.500 000
计算 $α 2$	1.000 007	0.900 001	0.800 000	0.699 999	0.600 000	0.500 000
流量误差/%	0.000 7	0.000 1	0	-0.000 1	0	0

表 2

表 3

计算域缩减影响算例和取样位置设置

算例编号	W/m	$L i n / W$	$L m / W$	$L s / W$	出口边界条件	取样位置x/m
算例编号	W/m	$L i n / W$	$L m / W$	$L s / W$	出口边界条件	1.0	1.5	3.0	6.0
1	1	3	10	10	充分发展	√	√	√	√
2	1	3	10	10	本文方法	√	√	√	√
3	1	3	10	5	本文方法	√	√	√	√
4	1	3	10	3	本文方法	√	√	√
5	1	3	10	2	本文方法	√	√	√
6	1	3	10	1	本文方法	√	√
7	1	3	10	0.5	本文方法	√	√

表 3

图 6

图 7

图 8

图 9

图 10

图 11

图 12

图 13

参考文献 18

1	HASAN N， ANWER S F， SANGHI S. On the outflow boundary condition for external incompressible flows： A new approach［J］. Journal of Computational Physics， 2005， 206（2）： 661-683.
2	DONG S， KARNIADAKIS G E， CHRYSSOSTOMIDIS C. A robust and accurate outflow boundary condition for incompressible flow simulations on severely-truncated unbounded domains［J］. Journal of Computational Physics， 2014， 261： 83-105.
3	SANI R L， GRESHO P M. Résumé and remarks on the open boundary condition minisymposium［J］. International Journal for Numerical Methods in Fluids， 1994， 18（10）： 983-1008.
4	ORLANSKI I. A simple boundary condition for unbounded hyperbolic flows［J］. Journal of Computational Physics， 1976， 21（3）： 251-269.
5	SANI M， SAIDI M S. A lagged implicit segregated data reconstruction procedure to treat open boundaries［J］. Journal of Computational Physics， 2010， 229（14）： 5418-5431.
6	陶文铨. 数值传热学［M］. 2版. 西安：西安交通大学出版社， 2001： 232.
	TAO W Q. Numerical heat transfer［M］. 2nd ed. Xi’an： Xi’an Jiaotong University Press， 2001： 232 （in Chinese）.
7	DONG S C， SHEN J. A pressure correction scheme for generalized form of energy-stable open boundary conditions for incompressible flows［J］. Journal of Computational Physics， 2015， 291： 254-278.
8	STEGGEL N， ROCKLIFF N. Simulation of the effects of body shape on lock-in characteristics in pulsating flow by the discrete vortex method［J］. Journal of Wind Engineering and Industrial Aerodynamics， 1997， 69-71： 317-329.
9	PAPANASTASIOU T C， MALAMATARIS N， ELLWOOD K. A new outflow boundary condition［J］. International Journal for Numerical Methods in Fluids， 1992， 14（5）： 587-608.
10	何吉欢. 简论开口边界条件［J］. 上海大学学报（自然科学版）， 1998， 4（2）： 213-217.
	HE J H. On open boundary condition—A variational approach［J］. Journal of Shanghai University （Natural Science Edition）， 1998， 4（2）： 213-217 （in Chinese）.
11	LI P W， TAO W Q. Effects of outflow boundary condition on convective heat transfer with strong recirculating flow［J］. Wärme - und Stoffübertragung， 1994， 29（8）： 463-470.
12	XUE S C， W BARTON G. Incompressible fluid flow simulations with flow rate as the sole information at synthetic inflow and outflow boundaries［J］. International Journal for Numerical Methods in Fluids， 2015， 78（12）： 739-760.
13	FORMAGGIA L， VENEZIANI A， VERGARA C. Flow rate boundary problems for an incompressible fluid in deformable domains： Formulations and solution methods［J］. Computer Methods in Applied Mechanics and Engineering， 2010， 199（9-12）： 677-688.
14	YANG Z G， DONG S C. Multiphase flows of N immiscible incompressible fluids： An outflow/open boundary condition and algorithm［J］. Journal of Computational Physics， 2018， 366： 33-70.
15	LIU Q K， MO Z Y， ZHANG A Q， et al. JAUMIN： A programming framework for large-scale numerical simulation on unstructured meshes［J］. CCF Transactions on High Performance Computing， 2019， 1： 35-48.
16	HAYES R E， NANDAKUMAR K， NASR-EL-DIN H. Steady laminar flow in a 90 degree planar branch［J］. Computers & Fluids， 1989， 17（4）： 537-553.
17	SCHNEIDER C， DREIZLER A， JANICKA J. Fluid dynamical analysis of atmospheric reacting and isothermal swirling flows［J］. Flow， Turbulence and Combustion， 2005， 74（1）： 103-127.
18	FREITAG M， JANICKA J. Investigation of a strongly swirled unconfined premixed flame using LES［J］. Proceedings of the Combustion Institute， 2007， 31（1）： 1477-1485.

E-mail：hkxb@buaa.edu.cn

关于我们

期刊社服务

专业学科

封面文章

友情链接

主管单位：中国科学技术协会主办单位：中国航空学会北京航空航天大学

[1]	李伟平, 杨龙金. 某回流燃烧室火焰筒冷却性能分析[J]. 航空学报, 2023, 44(9): 127326-127326.
[2]	杨占刚, 隋政, 张起浩, 刘建英. 复合材料飞机接地回流网络网内压降分析[J]. 航空学报, 2022, 43(1): 324859-324859.
[3]	李超群, 李易, 张晨曦, 唐硕. 沿流向微结构沟槽流场直接数值模拟[J]. 航空学报, 2020, 41(11): 123628-123628.
[4]	王彤, 孙亮亮, 邵昱昌. 惯性粒子分离器清除流道PIV测量与分析[J]. 航空学报, 2016, 37(10): 2961-2969.
[5]	吴子牛, 白晨媛, 徐珊姝, 李娟, 林景, 陈梓钧, 姚瑶. 突然启动流动问题:从不可压到高超声速流动[J]. 航空学报, 2015, 36(8): 2588-2600.
[6]	王建勇, 谢旅荣, 赵昊, 滕瑜琳. 一种改善高超声速进气道自起动能力的流场控制研究[J]. 航空学报, 2015, 36(5): 1401-1410.
[7]	葛立新;李椿萱. 旋涡破裂引起的流动非定常特性[J]. 航空学报, 1998, 19(5): 553-555.
[8]	李概奇;王家骅. 低压下V形火焰稳定器回流区流动特性研究[J]. 航空学报, 1989, 10(8): 447-449.
[9]	沈熊;于和生. 二维激光多普勒测速系统用于测量湍流分离流动[J]. 航空学报, 1984, 5(4): 393-400.
[10]	林其勋;丁毅生;杜琴芳;唐明;肖宁芳. 用双色激光多卜勒风速仪研究燃烧室的复杂流场[J]. 航空学报, 1984, 5(2): 233-240.

指定流量分配系数的多回流出口边界算法

Treatment of boundary condition at multiple outlets with recirculating flow and specified flow ratios

RichHTML

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

图/表 16

参考文献 18

相关文章 10

编辑推荐

Metrics

本文评价