期刊文献+
共找到671篇文章
< 1 2 34 >
每页显示 20 50 100
Numerical analysis of turbulent mixed convection air flow in inclined plane channel with k-εtype turbulence model 被引量:1
1
作者 XIE Zhengrui YANG Yanhua GU Hanyang CHENG Xu 《Nuclear Science and Techniques》 SCIE CAS CSCD 2008年第2期121-128,共8页
Numerical study on turbulent mixed convection in inclined plane channels,from 15° to 90° (vertical),was carried out to examine the effect of inclination on fluid flow and heat transfer distributions.The turb... Numerical study on turbulent mixed convection in inclined plane channels,from 15° to 90° (vertical),was carried out to examine the effect of inclination on fluid flow and heat transfer distributions.The turbulent air flows upward or downward into the duct with one wall heated from bottom.Calculation results with several kinds of k-εtype turbulence models were used to compare the experimental data with those in literatures to determine suitable model.The dependents of Nusselt number on the inclination angle of both the buoyancy-aided and buoyancy-opposed flow are discussed. 展开更多
关键词 混合对流 k-ε模型 湍流模型 主测线渠道
下载PDF
Numerical prediction of inner turbulent flow in conical diffuser by using a new five-point scheme and DLR k-ε turbulence model 被引量:1
2
作者 蒋光彪 何永森 +1 位作者 舒适 肖映雄 《Journal of Central South University》 SCIE EI CAS 2008年第S1期181-186,共6页
The internal turbulent flow in conical diffuser is a very complicated adverse pressure gradient flow.DLR k-ε turbulence model was adopted to study it.The every terms of the Laplace operator in DLR k-ε turbulence mod... The internal turbulent flow in conical diffuser is a very complicated adverse pressure gradient flow.DLR k-ε turbulence model was adopted to study it.The every terms of the Laplace operator in DLR k-ε turbulence model and pressure Poisson equation were discretized by upwind difference scheme.A new full implicit difference scheme of 5-point was constructed by using finite volume method and finite difference method.A large sparse matrix with five diagonals was formed and was stored by three arrays of one dimension in a compressed mode.General iterative methods do not work wel1 with large sparse matrix.With algebraic multigrid method(AMG),linear algebraic system of equations was solved and the precision was set at 10-6.The computation results were compared with the experimental results.The results show that the computation results have a good agreement with the experiment data.The precision of computational results and numerical simulation efficiency are greatly improved. 展开更多
关键词 conical DIFFUSER turbulENT flow DLR k turbulence model 5-point scheme ALGEBRAIC MULTIGRID method(AMG)
下载PDF
REVISED k-ε TURBULENCE MODEL IN ELECTROMAGNETIC CONTINUOUS CASTING OF MELT 被引量:1
3
作者 H.F. Huo and B.K. LiDepartment of Thermal Engineering, School of Materials and Metallurgy, Northeastern University, Shenyang 110004, China. 《Acta Metallurgica Sinica(English Letters)》 SCIE EI CAS CSCD 2003年第5期395-400,共6页
The research is motivated by the ongoing the electromagnetic continuous casting of molten metal. The revised k-ε model considering the effect of magnetic field application was derived. The specific model equations fo... The research is motivated by the ongoing the electromagnetic continuous casting of molten metal. The revised k-ε model considering the effect of magnetic field application was derived. The specific model equations for the electromagnetic braking were used to calculate the velocity distribution in the continuous casting mold of steel. The results show that the revised k-ε model considering the effect of magnetic field application tends to suppress the production of turbulence and difference between the conventional and revised k-e model is small. 展开更多
关键词 revised k turbulence model electromagnetic continuous casting magnetic field
下载PDF
A K-εTWO-EQUATION TURBULENCE MODEL FOR THE SOLID-LIQUID TWO-PHASE FLOWS  被引量:1
4
作者 刘小兵 程良骏 《Applied Mathematics and Mechanics(English Edition)》 SCIE EI 1996年第6期523-531,共9页
A two-equation turbulence model has been dereloped for predicting two-phase flow the two equations describe the conserration of turbulence kinetic energy and dissipation rate of that energy for the incompressible carr... A two-equation turbulence model has been dereloped for predicting two-phase flow the two equations describe the conserration of turbulence kinetic energy and dissipation rate of that energy for the incompressible carrier fluid in a two-phase flow The continuity, the momentum, K and εequations are modeled. In this model,the solid-liquid slip veloeites, the particle-particte interactions and the interactions between two phases are considered,The sandy water pipe turbulent flows are sueeessfuly predicted by this turbulince model. 展开更多
关键词 solid-liqtlid two-phase. k-εtwo-equation turbulence model
下载PDF
Assessment of k–ε models using tetrahedral grids to describe the turbulent flow field of a PBT impeller and validation through the PIV technique 被引量:3
5
作者 Victor X. Mendoza-Escamilla Alejandro Alonzo-Garcia +3 位作者 Helvio R. Mollinedo Israel Gonzalez-Neria J. Antonio Yanez-Varela Sergio A. Martinez-Delgadillo 《Chinese Journal of Chemical Engineering》 SCIE EI CAS CSCD 2018年第5期942-956,共15页
In turbulence modeling, the RNG and Realizable models have important improvements in the turbulent production and dissipation terms in comparison to the Standard. The selection of the appropriate turbulence model has ... In turbulence modeling, the RNG and Realizable models have important improvements in the turbulent production and dissipation terms in comparison to the Standard. The selection of the appropriate turbulence model has an impact on the convergence and solution in STRs, and they are used in mixing, multiphase modeling or as starting solution of transient models as DES and LES. Although there are several studies with the pitched blade turbine(PBT) impeller, most of them used the Standard model as representative of all k–ε models, using structured hexahedral grids composed of low number of cells, and in some cases under axial symmetry assumptions.Accordingly, in this work the assessment of the Standard, RNG and Realizable models to describe the turbulent flow field of this impeller, using the Multiple Reference Frame(MRF) and Sliding Mesh(SM) approaches with tetrahedral domains in dense grids, is presented. This kind of cell elements is especially suitable to reproduce complex geometries. Flow velocities and turbulent parameters were verified experimentally by PIV and torque measurements. The three models were capable of predicting fairly the pumping number, the power number based on torque, and velocities. Although the RNG improved the predictions of the turbulent kinetic energy and dissipation rate, the Realizable model presented better performance for both approaches. All models failed in the prediction of the total dissipation rate, and a dependence of its value on the number of cells for the MRF was found. 展开更多
关键词 流动速度 PIV 格子 评价 模特儿 PBT 标准模型 技术
下载PDF
PARAMETER ESTIMATION OF ENGINEERING TURBULENCE MODEL 被引量:1
6
作者 钱炜棋 蔡金狮 《Acta Mechanica Sinica》 SCIE EI CAS CSCD 2001年第4期302-309,共8页
A parameter estimation algorithm is introduced and used to determine the parameters in the standard k-epsilon two equation turbulence model (SKE). It can be found from the estimation results that although the paramete... A parameter estimation algorithm is introduced and used to determine the parameters in the standard k-epsilon two equation turbulence model (SKE). It can be found from the estimation results that although the parameter estimation method is an effective method to determine model parameters, it is. difficult to obtain a set of parameters for SKE to suit all kinds of separated flow and a modification of the turbulence model structure should be considered. So, a new nonlinear k-e two-equation model (NNKE) is put forward in this paper and the corresponding parameter estimation technique is applied to determine the model parameters. By implementing the NNKE to solve some engineering turbulent flows, it is shown that NNKE is more accurate and versatile than SKE. Thus, the success of NNKE implies that the parameter estimation technique may have a bright prospect in engineering turbulence model research. 展开更多
关键词 turbulence model parameter estimation k-epsilon two-equation model
下载PDF
Heat Transfer and Flow Characteristics Predictions with a Refined k-ε-f_u Turbulent Model in Impinging Jet 被引量:1
7
作者 Qinglin Niu Biao Chen +2 位作者 Zhihong He Jianfei Tong Shikui Dong 《Journal of Harbin Institute of Technology(New Series)》 EI CAS 2019年第4期9-17,共9页
Local heat transfer and flow characteristics in a round turbulent impinging jet for Re≈23 000 is predicted numerically with the RANS approach and a k-ε-fu turbulence model. The heat transfer predictions and turbulen... Local heat transfer and flow characteristics in a round turbulent impinging jet for Re≈23 000 is predicted numerically with the RANS approach and a k-ε-fu turbulence model. The heat transfer predictions and turbulence parameters are verified against the axis-symmetric free jet impingement measurements and compared with previous other turbulence models, and results show the k-ε-fu model has a good performance in predictions of the local wall heat transfer coefficient, and in agreement with measurements in mean velocity profiles at different radial positions as well. The numerical model is further used to examine the effect of the fully confined impingement jet on the local Nusselt number. Local Nusselt profiles in x and y-centerlines for the target plate over three separation distances are predicted. Compared with the experimental data, the numerical results are accurate in the central domain around the stagnation region and present a consistent structure distribution. 展开更多
关键词 HEAT TRANSFER IMPINGEMENT flow k-ε-fu turbulence model Nusselt NUMBER
下载PDF
CURVAURE SENSITIVE NONLINEAR TURBULENCE MODEL
8
作者 钱炜祺 符松 《Acta Mechanica Sinica》 SCIE EI CAS CSCD 2002年第1期31-41,共11页
By modifying the Rodi assumption to take account of the influence of flow curvature, a new curvature modified algebraic stress model(CMASM) is de- veloped from the second moment closure in the generalized curvilinear ... By modifying the Rodi assumption to take account of the influence of flow curvature, a new curvature modified algebraic stress model(CMASM) is de- veloped from the second moment closure in the generalized curvilinear coordinate system. And the explicit form of this ASM, a new curvature modified nonlinear k-ε model (CMNKE), is derived in the orthogonal curvilinear coordinate system. This new nonlinear k-ε model is further validated by a numerical simulation of a two- dimensional U-type turnaround duct flow. The results show that the CMNKE can effectively capture the main characteristic of this curvature flow and simulate the damping effect of the shear stress by a convex curvature and the enhancing effect by a concave curvature. So, this model is a rational and effective simplification to the second moment closure. 展开更多
关键词 turbulence curvature effect algebraic stress model nonlinear k model
下载PDF
Solution of Self-similar Equations of the k-ε Model in the Shear Turbulent Mixing Problem and Its Numerical Simulation
9
作者 Vyacheslav P. Statsenko Yulia V. Tret'yachenko Yury V. Yanilkin 《Journal of Physical Science and Application》 2015年第6期377-395,共19页
下载PDF
考虑旋转和曲率影响的SST k-ω湍流模型改进 被引量:22
10
作者 任芸 刘厚林 +2 位作者 舒敏骅 吴贤芳 吴登昊 《农业机械学报》 EI CAS CSCD 北大核心 2012年第11期123-128,共6页
分析了离心泵中常用湍流模型的优缺点,采用考虑旋转和曲率影响的RNG k-ε模型和Realizable k-ε模型分别对SST k-ω湍流模型进行了改进,生成了两种新的湍流模型。以开源代码库OpenFOAM作为平台对改进的湍流模型实现了程序化,应用新的湍... 分析了离心泵中常用湍流模型的优缺点,采用考虑旋转和曲率影响的RNG k-ε模型和Realizable k-ε模型分别对SST k-ω湍流模型进行了改进,生成了两种新的湍流模型。以开源代码库OpenFOAM作为平台对改进的湍流模型实现了程序化,应用新的湍流模型分别对一低、中、高比转数的离心泵进行数值计算,并将修改前后的数值计算结果与试验结果从能量性能的角度进行了详细对比。分析结果表明,改进后的湍流模型能成功用于离心泵内流计算中,且计算精度较修改前更接近于试验值。 展开更多
关键词 离心泵 湍流模型 sst k-ω模型 数值模拟
下载PDF
k-ω SST两方程湍流模型中参数影响的初步分析 被引量:87
11
作者 周宇 钱炜祺 +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-ω SST模型对翼型失速特性的模拟能力 被引量:15
12
作者 文晓庆 柳阳威 +1 位作者 方乐 陆利蓬 《北京航空航天大学学报》 EI CAS CSCD 北大核心 2013年第8期1127-1132,共6页
采用CFD软件fluent,k-ωSST模型对S825翼型进行了二维数值模拟.针对k-ωSST模型对翼型分离失速特性模拟不准的情况,对分离区域内湍流的强非平衡输运特性进行了分析.研究表明:原始模型由于不能准确模拟分离区内湍流的强非平衡输运特性,... 采用CFD软件fluent,k-ωSST模型对S825翼型进行了二维数值模拟.针对k-ωSST模型对翼型分离失速特性模拟不准的情况,对分离区域内湍流的强非平衡输运特性进行了分析.研究表明:原始模型由于不能准确模拟分离区内湍流的强非平衡输运特性,导致其对翼型失速特性模拟失效;提出了通过修正模型系数a1及β*∞,提高k-ωSST模型预测非平衡湍流输运特性的方法,从而提高对翼型失速特性的模拟精度. 展开更多
关键词 S825翼型 ksst模型 失速特性 湍流输运特性
下载PDF
应用SSTk-ω湍流模型计算Rushton搅拌釜流场 被引量:12
13
作者 韩路长 刘跃进 +2 位作者 李毅斌 李保华 罗和安 《化学工程》 CAS CSCD 北大核心 2006年第3期28-32,共5页
应用SSTkω-湍流模型对Rushton搅拌釜流场进行了计算,并与文献报道实验数据以及标准k-ε模型、RNGkω-模型预测结果进行了比较。结果表明SSTkω-模型所计算的排除流量准数与实验值相差较小,桨叶区速度场与实验数据吻合较好,略优于标准k... 应用SSTkω-湍流模型对Rushton搅拌釜流场进行了计算,并与文献报道实验数据以及标准k-ε模型、RNGkω-模型预测结果进行了比较。结果表明SSTkω-模型所计算的排除流量准数与实验值相差较小,桨叶区速度场与实验数据吻合较好,略优于标准k-ε模型,同时也给出了近壁区的计算结果。在湍动动能的预报上,SSTk-ω模型值与标准k-ε模型比较接近,优于RNGkε-模型,但与实验值都还不同程度地存在着一些差异。 展开更多
关键词 sst k-ω湍流模型 Rushton搅拌釜 流场计算
下载PDF
基于SST k-ω模型的汽车空调CFD仿真分析 被引量:9
14
作者 范平清 赵林林 邢彦锋 《机械强度》 CAS CSCD 北大核心 2014年第4期587-591,共5页
以汽车空调为研究目标,利用SST k-ω两方程模型,在除霜模式下分别对带HVAC及不带HVAC的空调系统进行仿真计算,分析了HVAC系统对仿真计算结果产生的影响。此外,分析结果显示位于驾驶员侧的前挡玻璃表面速度满足要求,位于副驾驶侧的前挡... 以汽车空调为研究目标,利用SST k-ω两方程模型,在除霜模式下分别对带HVAC及不带HVAC的空调系统进行仿真计算,分析了HVAC系统对仿真计算结果产生的影响。此外,分析结果显示位于驾驶员侧的前挡玻璃表面速度满足要求,位于副驾驶侧的前挡风玻璃表面速度大于1.5 m/s的区域覆盖率仅为75%左右,不满足国标要求。侧窗玻璃表面速度均值不到1 m/s,也不满足国标,仿真结果与试验结果对比相一致。风道存在较大问题,导致除霜风道气体集中从中间风道吹出,两侧风道的气流极小。根据结果提出进一步的改进建议,即在除霜风道增加必要的隔板装置,为汽车空调的改进设计提供依据。 展开更多
关键词 sst k-ω湍流 模型汽车空调 除霜风道 HVAC
下载PDF
基于改进k-ω SST模型的风力机尾流数值模拟 被引量:9
15
作者 杨祥生 赵宁 田琳琳 《太阳能学报》 EI CAS CSCD 北大核心 2017年第4期920-927,共8页
基于双方程k-ω SST模型提出两种改进湍流模型,采用丹麦科技大学的自编求解器EllipSys3D结合致动盘方法,对单台风力机尾流流场进行数值模拟。该文将SST-sust模型引入到风力机尾流的数值研究中,为了更进一步增大预测到的湍流强度和提高... 基于双方程k-ω SST模型提出两种改进湍流模型,采用丹麦科技大学的自编求解器EllipSys3D结合致动盘方法,对单台风力机尾流流场进行数值模拟。该文将SST-sust模型引入到风力机尾流的数值研究中,为了更进一步增大预测到的湍流强度和提高尾流恢复速度,提出SST-Csust模型和SST-Dsust模型,它们是在SST-sust模型的基础上分别调整湍流模型的封闭常系数、修正ω方程的耗散项。通过与理论、实验以及LES模型进行对比分析,结果表明该文提出的改进模型在尾流区速度分布和湍流强度分布的预测方面可取得较好的应用效果。 展开更多
关键词 k sst湍流模型 风力机尾流 EllipSys3D求解器 致动盘(AD)模型
下载PDF
基于SST k-ω湍流模型的平衡大气边界层模拟 被引量:14
16
作者 胡朋 李永乐 廖海黎 《空气动力学学报》 EI CSCD 北大核心 2012年第6期737-743,共7页
为实现SST k-ω湍流模型下空流场大气边界层的平衡,提出通过在湍动能输送方程与比耗散率输送方程中添加源项的方法使来流边界条件与湍流模型相协调,并通过相应的CFD数值算例分别考察了不考虑源项与考虑源项时沿流场不同位置处的平均速... 为实现SST k-ω湍流模型下空流场大气边界层的平衡,提出通过在湍动能输送方程与比耗散率输送方程中添加源项的方法使来流边界条件与湍流模型相协调,并通过相应的CFD数值算例分别考察了不考虑源项与考虑源项时沿流场不同位置处的平均速度、湍动能及比耗散率剖面的一致性。计算结果表明,当考虑源项时沿流场不同位置处的平均速度剖面、湍动能剖面及比耗散率剖面均有较高程度的一致性,较好地实现了大气边界层的平衡。研究结论为平衡大气边界层的研究提供了新的思路和研究方法,对结构风工程中CFD数值模拟有一定的理论及应用价值。 展开更多
关键词 平衡大气边界层 sst k-ω湍流模型 源项 计算流体动力学
下载PDF
SST k-ω-k_p两相湍流模型及其在湿蒸汽凝结流动数值模拟中的应用 被引量:3
17
作者 吴晓明 李国君 +1 位作者 丰镇平 李亮 《西安交通大学学报》 EI CAS CSCD 北大核心 2007年第5期526-529,584,共5页
首先,建立了湿蒸汽凝结流动的双流体模型,考虑了湿蒸汽两相流动中相间的速度滑移、耦合以及湍流扩散作用的影响.然后,针对蒸汽透平叶栅中流动的湍流特性,在单相湍流计算中数值模拟精度相对良好的两方程SST k-ω湍流模型基础上,参照颗... 首先,建立了湿蒸汽凝结流动的双流体模型,考虑了湿蒸汽两相流动中相间的速度滑移、耦合以及湍流扩散作用的影响.然后,针对蒸汽透平叶栅中流动的湍流特性,在单相湍流计算中数值模拟精度相对良好的两方程SST k-ω湍流模型基础上,参照颗粒湍能输运方程理论,推导建立了SST k-ω-k_p湿蒸汽两相流动湍流模型,模型中引入了液相黏性、导热及扩散系数等拟流体概念.对二维叶栅中存在自发凝结的湿蒸汽流动进行的数值模拟表明:建立的数值模型表现良好,明显优于单流体模型,特别是吸力面的压力分布模拟精度有很大提高,具有较高的精确性和可靠性.由于实验数据的匮乏,数值模拟中发现的汽液两相的动力特性以及液相的参数分布有待于进一步的实验来验证. 展开更多
关键词 sst k-ω-kp 两相湍流模型 湿蒸汽 凝结流动 数值模拟
下载PDF
基于SST k-ω湍流模型的高层建筑体型系数数值研究 被引量:4
18
作者 王强 陈超 胡荣华 《沈阳建筑大学学报(自然科学版)》 CAS 北大核心 2012年第3期417-422,共6页
目的研究建筑物高度和宽度的变化对其表面绕流流场的影响,分析建筑物风荷载体型系数的变化,为结构设计人员在结构抗风设计方面提供参考.方法利用计算流体动力学软件STAR-CCM+,选用基于雷诺平均(RANS)的SST k-ω湍流模型,对不同高度和宽... 目的研究建筑物高度和宽度的变化对其表面绕流流场的影响,分析建筑物风荷载体型系数的变化,为结构设计人员在结构抗风设计方面提供参考.方法利用计算流体动力学软件STAR-CCM+,选用基于雷诺平均(RANS)的SST k-ω湍流模型,对不同高度和宽度的建筑物进行数值模拟,计算建筑物风荷载体型系数.结果在宽度相同的情况下,随着高度的增加,正面、侧面和背面的体型系数绝对值增加;随着宽度的增加,侧面和背面的体型系数绝对值增加,而正面影响不大.结论高层建筑的体型和体量对风荷载体型系数有一定的影响,在结构抗风设计中应该适当予以考虑. 展开更多
关键词 高层建筑 体型系数 sst k-ω湍流模型 数值模拟
下载PDF
基于SST k-ω湍流模型的二维圆柱涡激振动数值仿真计算 被引量:6
19
作者 李骏 李威 《舰船科学技术》 北大核心 2015年第2期30-34,共5页
随着海洋工程逐渐向深海发展,广泛使用的柔性立管因为高频振动很容易受到严重的损伤。因此对于这类细长柔性结构的涡激振动(VIV)研究是当今的一个热点,同时随着计算机的快速发展,CFD(计算流体力学)技术成为研究涡激振动问题不可或缺的... 随着海洋工程逐渐向深海发展,广泛使用的柔性立管因为高频振动很容易受到严重的损伤。因此对于这类细长柔性结构的涡激振动(VIV)研究是当今的一个热点,同时随着计算机的快速发展,CFD(计算流体力学)技术成为研究涡激振动问题不可或缺的一种方法。本文采用雷诺平均纳维尔-斯托克斯(RANS)方程,并结合SST k-ω湍流模型,研究低质量比弹性支撑刚性圆柱体的涡激振动问题。从振幅响应、频率响应、3个响应分支的水动力性能、尾涡模式等方面和Williamson相关实验作对比。结果表明,SST k-ω湍流模型能够有效准确地模拟圆柱绕流的涡激振动。本文丰富了海洋工程的理论研究,为柔性立管的实际应用提供了一定的理论指导。 展开更多
关键词 涡激振动 sst k-ω湍流模型 立管
下载PDF
基于SST k-ω湍流模型与遗传算法的高速电梯井道风口优化设计 被引量:7
20
作者 邱望标 郭天水 +1 位作者 薛玉 付静 《液压与气动》 北大核心 2019年第5期51-57,共7页
为了改善和提高电梯的气动特性和运行稳定性,基于SST k-ω湍流模型对电梯进行空气流体仿真,计算气动阻力、侧翻力矩、俯仰力矩,分析实验数据变化趋势。通过最小二乘拟合,建立电梯气动特性的一元数值模型。根据正交实验结果与一元数值模... 为了改善和提高电梯的气动特性和运行稳定性,基于SST k-ω湍流模型对电梯进行空气流体仿真,计算气动阻力、侧翻力矩、俯仰力矩,分析实验数据变化趋势。通过最小二乘拟合,建立电梯气动特性的一元数值模型。根据正交实验结果与一元数值模型,基于MATLAB回归分析,建立电梯气动特性的多元数值模型,并以此作为目标函数,以气动阻力、侧翻力矩和俯仰力矩最小作为优化目标,基于遗传算法进行电梯气动特性的多目标优化,仿真验证优化后电梯气动特性得到较大改善,为电梯气动特性的改善提供了一定的参考价值和意义。 展开更多
关键词 高速电梯 遗传算法 sstk-ω湍流模型 回归分析 风口优化
下载PDF
上一页 1 2 34 下一页 到第
使用帮助 返回顶部