期刊文献+
共找到25篇文章
< 1 2 >
每页显示 20 50 100
k-ω SST两方程湍流模型中参数影响的初步分析 被引量:93
1
作者 周宇 钱炜祺 +1 位作者 邓有奇 马明生 《空气动力学学报》 EI CSCD 北大核心 2010年第2期213-217,共5页
Menter的k-ωSST两方程湍流模型在流体力学计算中有良好表现。利用均匀试验设计方法对k-ωSST模型中八个参数在小攻角亚声速、跨声速和大攻角三种典型的翼型绕流流场计算中的影响规律进行了分析。分析结果表明在附着流中参数a1是对流场... Menter的k-ωSST两方程湍流模型在流体力学计算中有良好表现。利用均匀试验设计方法对k-ωSST模型中八个参数在小攻角亚声速、跨声速和大攻角三种典型的翼型绕流流场计算中的影响规律进行了分析。分析结果表明在附着流中参数a1是对流场的影响最大的参数;而在分离流中,参数a1对流场的影响与湍流产生与耗散等项的作用相当,模型中的四个参数:σω2、β2、β*和a1的取值都会对计算结果有较大的影响。 展开更多
关键词 k-ω SST方程湍流模型 翼型 均匀试验设计
在线阅读 下载PDF
基于新的描述湍流耗散方程的k-ζ两方程湍流模型的数值算法研究 被引量:6
2
作者 张强 杨永 李喜乐 《西北工业大学学报》 EI CAS CSCD 北大核心 2009年第4期466-470,共5页
对Navier-Stokes方程进行雷诺平均后出现的各关联项建模,通过新的描述湍流脉动耗散的变量构造耗散方程,建立k-ζ两方程湍流模型,研究了k-ζ两方程湍流模型的数值求解方法。通过求解有限体积法离散的RANS流动控制方程,数值模拟了平板,翼... 对Navier-Stokes方程进行雷诺平均后出现的各关联项建模,通过新的描述湍流脉动耗散的变量构造耗散方程,建立k-ζ两方程湍流模型,研究了k-ζ两方程湍流模型的数值求解方法。通过求解有限体积法离散的RANS流动控制方程,数值模拟了平板,翼型,机翼等不同湍流流场,并与理论解、实验值及SSTk-ω模型进行比较,全面考察了k-ζ两方程湍流模型在湍流流场计算中的准确性及适用性。数值计算表明,通过建立新的耗散方程研究湍流的方法是可行的,目前的k-ζ两方程湍流模型具有良好的数值稳定性,并且计算结果要优于或者至少与传统的两方程模型精度相当。 展开更多
关键词 NAVIER-STOKES方程 雷诺平均 k-ζ方程湍流模型
在线阅读 下载PDF
一种新的非线性K-ε两方程湍流模型 被引量:1
3
作者 钱炜祺 蔡金狮 《空气动力学学报》 CSCD 北大核心 1999年第4期391-397,共7页
本文基于理论分析和对标准K-ε两方程湍流模型数值模拟结果的分析,在标准K-ε两方程模型的雷诺应力表达式中加入高阶非线性项;并在模型的涡粘性系数μt 中考虑流动各向异性的影响,将涡粘性系数中参数Cμ取为表征流动各向异性... 本文基于理论分析和对标准K-ε两方程湍流模型数值模拟结果的分析,在标准K-ε两方程模型的雷诺应力表达式中加入高阶非线性项;并在模型的涡粘性系数μt 中考虑流动各向异性的影响,将涡粘性系数中参数Cμ取为表征流动各向异性的参数A 的一个线性函数,从而得到了一个新的非线性K-ε两方程湍流模型。此外,本文利用两种不同扩比下的后向台阶流动的实验结果来对新模型中引入的待定参数进行了参数辨识,并用另一种扩比下的后向台阶流动、绕方块的分离流动、U 型管道内流动、以及接近失速攻角情况下绕GAW- 1 翼型的低速流动等四种流动来对该非线性K-ε两方程模型作了进一步的验证,得到了令人满意的结果。 展开更多
关键词 非线性 k-ε方程模型 雷诺应力表达式 湍流
在线阅读 下载PDF
用非线性 K-ε 两方程湍流模型计算后向台阶分离流动
4
作者 钱炜祺 蔡金狮 《空气动力学学报》 CSCD 北大核心 1998年第4期426-432,共7页
本文利用TBGatski和CGSpeziale根据张量不变性理论推出的二维流动情况下的非线性K-ε两方程模型来对两种不同扩比下的后向台阶分离湍流进行了数值模拟。从计算结果可以看出:用非线性K-ε两方程模型计算得到的回... 本文利用TBGatski和CGSpeziale根据张量不变性理论推出的二维流动情况下的非线性K-ε两方程模型来对两种不同扩比下的后向台阶分离湍流进行了数值模拟。从计算结果可以看出:用非线性K-ε两方程模型计算得到的回流区长度、壁面压力分布和摩阻分布等物理量的精度比用标准K-ε两方程模型计算得到的结果有了明显改进。而且,该模型还具有标准K-ε两方程模型的结构简单、使用方便的优点,易于工程应用。 展开更多
关键词 湍流 非线性 k-ε方程模型 分离流动
在线阅读 下载PDF
高浓度固-液两相流紊流的动理学模型 被引量:13
5
作者 唐学林 徐宇 吴玉林 《力学学报》 EI CSCD 北大核心 2002年第6期956-962,共7页
采用分子动理学方法,基于固-液两相流液相分子或颗粒相颗粒的Boltzmann方程,对Boltzmann方程分别取零矩和一次矩,则得到高浓度固-液两相流紊流的连续方程和动量方程,再和较成熟的低浓度两相流连续方程和动量方程比较,取低浓度两相流控... 采用分子动理学方法,基于固-液两相流液相分子或颗粒相颗粒的Boltzmann方程,对Boltzmann方程分别取零矩和一次矩,则得到高浓度固-液两相流紊流的连续方程和动量方程,再和较成熟的低浓度两相流连续方程和动量方程比较,取低浓度两相流控制方程中较成熟合理的有关项和高浓度时由动理学方法推导出的颗粒间碰撞项,则得到高浓度固-液两相流紊流的最终控制方程:连续方程和动量方程. 展开更多
关键词 固-液相流 紊流 动理学模型 速度分布函数 BOLTZMANN方程
在线阅读 下载PDF
丁坝绕流的k-ε紊流模型数值解 被引量:11
6
作者 程年生 李昌华 《水利水运工程学报》 CSCD 1989年第3期11-23,共13页
本文采用雷诺方程组描述丁坝绕流运动,并用紊流k-ε双方程模型闭合雷诺方程组。根据控制体积法导出沿水深平均的二维偏微分方程的离散格式,然后用SIMPLER程式进行求解。通过计算,获得丁坝绕流流场内流速、水化、紊动能和紊动粘性系数的... 本文采用雷诺方程组描述丁坝绕流运动,并用紊流k-ε双方程模型闭合雷诺方程组。根据控制体积法导出沿水深平均的二维偏微分方程的离散格式,然后用SIMPLER程式进行求解。通过计算,获得丁坝绕流流场内流速、水化、紊动能和紊动粘性系数的断面分布规律,其中流速和水位的计算值与实测资料较为一致。 展开更多
关键词 k-ε紊流模型 数值解 离散格式 体积法 粘性系数 紊动 雷诺方程 方程模型 斯托克斯方程 雷诺应力
在线阅读 下载PDF
水沙两相流ASM模型在浑水异重流计算中的应用及模型试验研究 被引量:8
7
作者 李永 李嘉 安瑞冬 《四川大学学报(工程科学版)》 EI CAS CSCD 北大核心 2009年第4期102-108,共7页
挟沙水流与清水相遇时,由于前者的密度比后者大,在条件合适的时候,挟沙水流就会潜入清水底部继续向前流动,形成浑水异重流,是泥沙运动的一种特殊形式,特别是在水库的泥沙设计及运行排沙方面需要进行深入的研究。基于水沙两相流理论,将... 挟沙水流与清水相遇时,由于前者的密度比后者大,在条件合适的时候,挟沙水流就会潜入清水底部继续向前流动,形成浑水异重流,是泥沙运动的一种特殊形式,特别是在水库的泥沙设计及运行排沙方面需要进行深入的研究。基于水沙两相流理论,将含沙浑水作为一相均匀流体处理,引入相间滑移速度的概念,建立混合相连续方程、动量方程、第二相的体积分数方程和相间相对速度方程的ASM(Algebraic Slip Mixture)两相流模型,并耦合求解水动力学计算的标准k-ε方程组。还介绍了针对异重流研究设计的多相流试验水槽及模型试验方法,对ASM模型计算的结果进行试验验证并分析模型适用条件。结果表明:该模型可以对异重流潜入,前锋形态及推进速度进行有效模拟;对Re数大于2000及含沙量高于1 kg/m3的异重流条件可以进行较精确的模拟。研究为深入了解水库、湖泊及河流入海口水沙异重流问题提供了新的途径。 展开更多
关键词 异重流 水沙相流 相间滑移 ASM模型 标准k-ε方程
在线阅读 下载PDF
贯流式水轮机活动导叶内部三维固液两相紊流研究 被引量:4
8
作者 刘文俊 吴玉林 曹树良 《水利水电技术》 CSCD 北大核心 1999年第5期75-77,共3页
针对贯流式水轮机活动导叶内部固液两相亲流的三维数值模拟问题,基于N-S方程和两相紊流的。。一个模型,采用贴体坐标和交错网格系统,用SIMPLEC算法对水轮机活动导叶内部的固液两相亲流进行数值模拟,提出设计工况时活动导叶内部的... 针对贯流式水轮机活动导叶内部固液两相亲流的三维数值模拟问题,基于N-S方程和两相紊流的。。一个模型,采用贴体坐标和交错网格系统,用SIMPLEC算法对水轮机活动导叶内部的固液两相亲流进行数值模拟,提出设计工况时活动导叶内部的两相紊流特性.得出了导叶内部流场和压力场的分布.将计算结果用Tecplot7进行处理,用图形表示其内部的流动状态和压力分布. 展开更多
关键词 活动导叶 紊流 k-ε-Ap模型 磨损 水轮机
在线阅读 下载PDF
基于混合网格Navier-Stokes方程的并行隐式计算方法研究 被引量:11
9
作者 郑冠男 邓守春 +1 位作者 韩同来 杨国伟 《应用力学学报》 CAS CSCD 北大核心 2011年第3期211-218,323,共8页
针对结构网格很难处理复杂外形和非结构网格无法计算具有边界层的粘性流动的缺点,发展了基于混合网格格点的隐式算法,成功地解决了在工程应用中难于处理的复杂外形粘性流场计算和效率问题。同时针对大规模的工程问题,发展了基于MPI通信... 针对结构网格很难处理复杂外形和非结构网格无法计算具有边界层的粘性流动的缺点,发展了基于混合网格格点的隐式算法,成功地解决了在工程应用中难于处理的复杂外形粘性流场计算和效率问题。同时针对大规模的工程问题,发展了基于MPI通信技术的染色分层通讯并行计算方法。其中空间离散采用基于Roe格式发展的三阶迎风HLLEW(Harten-Lax-Van Leer-Einfeldt-Wada)或AUSM格式,湍流模型采用k??两方程湍流模型,时间推进考虑到LU-SGS并行等效较困难则采用基于DP-LUR(Data-Parallel Lower-Upper Relaxation)格式的隐式算法,计算CFL数可取到105量级,从2个到128个CPU的并行加速效率都保持在90%以上,大大提高了计算效率。算例对标模M6机翼模型流场进行计算,验证了方法的可靠性;然后对标模DLR-F6翼身组合体进行混合网格粘性与无粘计算结果进行比较,进一步验证混合网格方法;最后计算了DLR-WBNP外挂发动机翼身组合体模型,准确模拟了外挂和超临界机翼的相互干扰流动问题,采用4 CPU 16 CORE到24 CPU 96 CORE,2000步计算时间都不超过3小时。为民机跨声速气动弹性分析的计算效率提升提供了基本的数值模拟工具。 展开更多
关键词 混合网格 k-ω方程湍流模型 DP-LUR隐式方法 MPI 染色分层通讯
在线阅读 下载PDF
泄水建筑物体型优化与消能工模型试验研究 被引量:1
10
作者 张铎 党宏祖 《农业科技与信息》 2015年第12期106-109,共4页
本文利用引入VOF方法的RNG k-ε双方程紊流数学模型,对溢流表孔和冲砂底孔进行二维数值模拟。通过模拟,得出计算域的流速场和压强场水力特性,并结合物理模型试验对数值计算结果进行验证。作者分析了不同工况下水力特性的变化规律,论证... 本文利用引入VOF方法的RNG k-ε双方程紊流数学模型,对溢流表孔和冲砂底孔进行二维数值模拟。通过模拟,得出计算域的流速场和压强场水力特性,并结合物理模型试验对数值计算结果进行验证。作者分析了不同工况下水力特性的变化规律,论证泄水体型的合理性。另外,经系列物理模型试验对消能工进行优化。结果表明,物理模型试验与数值计算吻合较好。冲砂底孔出口段,采用合理的收缩比将水流纵向拉开,起到分散水流的作用,消能效果较好,较原设计方案有明显改善。同时,为类似工程泄水建筑物体型优化的研究提供参考依据。 展开更多
关键词 RNG k-ε两方程紊流模型 VOF方法 模型试验 体型优化 消能
在线阅读 下载PDF
对称开口空间紊流混合对流的数值模拟及分析
11
作者 黄永红 何杰 《长沙水电师院学报(自然科学版)》 1998年第4期378-383,共6页
对工程实际中广为应用的紊流混合对流换热进行了数值研究,具有计算对象为具有对称边界条件的多进出口空调房间的传热.针对此对象编制了一个应用程序,并运用此程序在较宽的参数范围内对不同的配风方式及边界条件进行了数值模拟,得到... 对工程实际中广为应用的紊流混合对流换热进行了数值研究,具有计算对象为具有对称边界条件的多进出口空调房间的传热.针对此对象编制了一个应用程序,并运用此程序在较宽的参数范围内对不同的配风方式及边界条件进行了数值模拟,得到了相应的速度场、温度场及其它物理量场的分析,进而总结出紊流混合对流换热的一些规律.通过与相关文献的对比,说明了本文数值模拟的有效性. 展开更多
关键词 混合对流换热 紊流 k-ε方程模型
在线阅读 下载PDF
考虑气温变化影响的引水渠道水内冰演变数值模拟 被引量:13
12
作者 王晓玲 周正印 +2 位作者 蒋志勇 周莎莎 张自强 《天津大学学报》 EI CAS CSCD 北大核心 2010年第6期515-522,共8页
高寒地区冰害问题是关系到引水式水电站能否正常运行的关键问题.为此,建立了三维非稳态Euler-Euler两相流k-ε紊流模型,在动量方程中考虑了相间曳力、升力、虚拟质量力,在热量传递方程中考虑了冰水之间的热量传递以及太阳辐射、有效辐... 高寒地区冰害问题是关系到引水式水电站能否正常运行的关键问题.为此,建立了三维非稳态Euler-Euler两相流k-ε紊流模型,在动量方程中考虑了相间曳力、升力、虚拟质量力,在热量传递方程中考虑了冰水之间的热量传递以及太阳辐射、有效辐射、蒸发热损失和水面对流热损失等.以新疆某水电站引水渠道为例,模拟分析了气温变化条件下流速、水温、冰温及冰体积分数的沿程分布.结果表明:水温沿程先降低后增大,沿程出现最低点,冰温沿程逐渐增大;冰体积分数沿程存在突变点,且突变点在水温最低点附近;随着气温降低,最终的产冰量和冰体积分数逐渐增大.讨论了太阳辐射、有效辐射和蒸发热损失对水内冰分布的影响,并通过与吴剑疆等的黄河上游河道中水内冰演变计算结果及经验公式计算产冰量的对比分析,验证了模型的可靠性. 展开更多
关键词 气温 水内冰演变 非稳态Euler-Euler相流 k-ε紊流模型 引水渠道
在线阅读 下载PDF
基于VOF法的长距离无压引水隧洞充水过程模拟 被引量:8
13
作者 于跃 王晨晨 +2 位作者 安娟 李涛 段琦琦 《水利水电技术》 CSCD 北大核心 2009年第4期36-40,共5页
结合长距离无压引水隧洞充水过程,建立了水气两相流耦合VOF(Volume of Fluid)法的k-ε紊流模型,近壁面采用考虑糙率影响的壁面函数,结合实际工程对长距离无压引水隧洞水气两相流进行了模拟。模拟结果表明:隧洞内水位随充水时间的增加逐... 结合长距离无压引水隧洞充水过程,建立了水气两相流耦合VOF(Volume of Fluid)法的k-ε紊流模型,近壁面采用考虑糙率影响的壁面函数,结合实际工程对长距离无压引水隧洞水气两相流进行了模拟。模拟结果表明:隧洞内水位随充水时间的增加逐渐上升并趋于稳定;对隧洞充水过程中沿程典型断面的水气两相流进行了动态分析,结果表明进口流量为65 m3/s时水流能顺利通过隧洞;对某水库放水兼放空洞水深实测数据进行了验证对比,模拟结果与实验结果基本吻合。 展开更多
关键词 长距离无压引水隧洞 水气相流 VOF法 k-ε紊流模型 充水过程
在线阅读 下载PDF
三棱柱形网格下自由表面流的三维数值模拟 被引量:6
14
作者 王昆 金生 +1 位作者 马志强 高述峰 《水科学进展》 EI CAS CSCD 北大核心 2009年第4期490-494,共5页
为了模拟自由表面流动问题,建立了基于三棱柱形网格的三维Navier-Stokes方程的半隐式动压离散模型。出流边界采用在动量方程右端附加线性衰减项的海绵层处理方式来消减波能。引入k-ε双方程紊流模型求解涡黏性系数,使方程达到封闭。用... 为了模拟自由表面流动问题,建立了基于三棱柱形网格的三维Navier-Stokes方程的半隐式动压离散模型。出流边界采用在动量方程右端附加线性衰减项的海绵层处理方式来消减波能。引入k-ε双方程紊流模型求解涡黏性系数,使方程达到封闭。用规则周期波通过淹没障碍物的传播变形和复合明渠流动算例,验证了当遇到短波高频问题时,所述模型模拟结果与测量值吻合良好,体现了模型的精确性和有效性。 展开更多
关键词 三棱柱形网格 自由表面流 NAVIER-STOKES方程 k-ε紊流模型
在线阅读 下载PDF
复杂无压隧洞输水过程三维数值模拟与方案优选 被引量:4
15
作者 盛代林 孙月峰 +2 位作者 安娟 段琦琦 刘雪朋 《天津大学学报》 EI CAS CSCD 北大核心 2009年第12期1048-1054,共7页
运用基于流体体积(VOF)法的三维k-ε双方程紊流模型模拟了无压引水隧洞水气两相流紊流流场,近壁面采用考虑糙率影响的壁函数法进行处理.以某复杂长距离无压引水隧洞的输水过程为例,模拟了3种方案的断面平均过水面积、流量、流速和气压... 运用基于流体体积(VOF)法的三维k-ε双方程紊流模型模拟了无压引水隧洞水气两相流紊流流场,近壁面采用考虑糙率影响的壁函数法进行处理.以某复杂长距离无压引水隧洞的输水过程为例,模拟了3种方案的断面平均过水面积、流量、流速和气压的变化,实现了多方案的比选.结果表明3,种方案水流均能顺利通过隧洞;方案二水流到达隧洞出口的时间最短;方案二的断面平均过水面积、流量、流速和气压均优于方案一;方案三与方案二差别较小,综合考虑过流能力及成本因素,方案二为最优方案.通过与某水库放水兼放空洞水深实测数据及采用传统水力学方法计算的经验值的对比分析,平均相对误差分别为5.01%和0.96%,验证了模型的可靠性. 展开更多
关键词 无压隧洞 输水过程 k-ε紊流模型 水气相流 方案优选
在线阅读 下载PDF
CFD法设计轴流式油气混输泵初探 被引量:15
16
作者 马希金 郭俊杰 吴蓓 《流体机械》 CSCD 2004年第7期15-18,共4页
选用标准的K -ε紊流模型 ,利用先进的CFD软件FLUENT及其预处理软件GAMBIT ,模拟了轴流式油气混输泵样机的叶轮及混流器内的三维紊流场 ,依据分析结果对叶轮及混流器相关的几何参数进行了优化 ,并对模型的选择、三维建模、网格划分、计... 选用标准的K -ε紊流模型 ,利用先进的CFD软件FLUENT及其预处理软件GAMBIT ,模拟了轴流式油气混输泵样机的叶轮及混流器内的三维紊流场 ,依据分析结果对叶轮及混流器相关的几何参数进行了优化 ,并对模型的选择、三维建模、网格划分、计算分析和试验结果作了介绍。 展开更多
关键词 油气混输泵 相流 流场分析 设计 CFD k-ε紊流模型 三维建模 网格划分 油田
在线阅读 下载PDF
绕翼型低雷诺数流动的数值分析研究 被引量:14
17
作者 张强 杨永 《空气动力学学报》 EI CSCD 北大核心 2006年第4期482-486,共5页
分离泡的产生是绕翼型低雷诺数流动的一个重要特征,分离泡通常都是非定常的,会对整个流场产生极大的影响,由于分离附面层的不稳定性,很快诱发转捩,并产生湍流再附。文中通过求解雷诺平均N-S方程,数值模拟了绕Eppler387翼型的低雷诺数非... 分离泡的产生是绕翼型低雷诺数流动的一个重要特征,分离泡通常都是非定常的,会对整个流场产生极大的影响,由于分离附面层的不稳定性,很快诱发转捩,并产生湍流再附。文中通过求解雷诺平均N-S方程,数值模拟了绕Eppler387翼型的低雷诺数非定常流动,并对两方程SSTk-ω湍流模型,代数B-L模型和层流的计算结果作了比较。N-S方程和两方程湍流模型的控制方程非耦合求解,时间推进均采用近似因子方法(AF);空间离散无粘项采用ROE格式,粘性项用中心差分方法计算;计算非定常流场时,采用伪时间子迭代(τt-s)方法保证二阶时间精度,湍流模型计算时都是在固定点转捩。最后分析了转捩对低雷诺数流动的影响、分离泡的存在引起的流动不稳定性和周期性的脱出涡,数值计算给出的时间平均升力系数、阻力系数和压力分布与实验结果比较吻合。 展开更多
关键词 分离泡 低雷诺数 转捩 方程SST k-ω湍流模型 伪时间子迭代(τ-ts)
在线阅读 下载PDF
泵站出水流道的数模分析 被引量:4
18
作者 朱荣生 李维斌 +1 位作者 黄道见 曹卫东 《农机化研究》 北大核心 2003年第3期79-80,84,共3页
将CFD技术应用于泵站出水流道的分析与设计,基于Navier-Stokes方程和标准k-ε紊流模型,求解泵站出水流道内湍流流场,数值模拟的结果对泵站出水流道水力优化设计具有指导意义。
关键词 泵站 出水流道 数模分析 Navier—Stokes方程 标准k-ε紊流模型 水力优化设计 计算流体力学 CFD技术
在线阅读 下载PDF
氧碘化学激光器组分超音速混合反应数值模拟 被引量:5
19
作者 张黎 叶正寅 王刚 《光电工程》 EI CAS CSCD 北大核心 2007年第6期20-24,29,共6页
采用数值方法求解二维可压缩非定常N-S方程及组分连续方程,计算氧碘化学激光器中的组分超音速混合反应流场。空间离散格式为AUSM+-up格式,用四步龙格-库塔方法作显式时间推进,湍流粘性系数使用k-ε两方程湍流模型进行求解。假设混合气... 采用数值方法求解二维可压缩非定常N-S方程及组分连续方程,计算氧碘化学激光器中的组分超音速混合反应流场。空间离散格式为AUSM+-up格式,用四步龙格-库塔方法作显式时间推进,湍流粘性系数使用k-ε两方程湍流模型进行求解。假设混合气体为热力学完全而热值非完全气体,化学反应模型采用有限速率反应模型。使用一种松弛迭代的方法来处理化学源项的刚性问题。给出了出口处碘分子的质量分数分布,数值模拟结果与文献结果符合良好。 展开更多
关键词 氧碘化学激光器 有限速率反应模型 AUSM+-up格式 k-ε方程湍流模型
在线阅读 下载PDF
舰船普通排气系统的红外特性 被引量:5
20
作者 林维德 邢平平 洪家琳 《红外与毫米波学报》 SCIE EI CAS CSCD 北大核心 2000年第2期125-128,共4页
采用 k-ε两方程模型和壁面函数法计算舰船普通排气系统缩比模型的湍流场和温度场分布 .经实验校模后利用自模化原理将数值模拟方法推广到实物装置上去 .在建立了观察点与系统有关区域间看视关系后 ,利用
关键词 k-ε方程模型 红外特性模型 舰船排气系统
在线阅读 下载PDF
上一页 1 2 下一页 到第
使用帮助 返回顶部