陷落腔的流场数值计算研究

本文对有限分析解法的近似计算方法及其对计算精度的影响进行了数值计算和分析讨论；然后采用有限分析解法，对前、后壁壁缘高度不同的三种腔内流场运动发展过程及腔口剪切层运动情况作了数值计算，得到的腔口剪切层上腔口随边相互作用情况与对应腔的流场显示实验结果进行了比较，两者一致，计算结果表明，采用本文发展的近似方法计算可靠，可以较大地提高计算效率。     

集成光学陀螺的最佳腔长

作为谐振式光学陀螺的核心敏感器件,无源环形谐振腔(PRR)结构参数的选择,对陀螺灵敏度具有决定性作用.腔长正是这重要结构参数之一.基于多光束干涉理论的理论分析表明,谐振腔长的增加,一方面通过决定腔内干涉光程,优化PRR谐振特性,另一方面,又通过对腔光传输损耗的影响,决定着谐振条件的耦合比要求,从而成为腔内光传输总损耗重要决定因素,导致谐振特性劣化.根据谐振式光学陀螺的调频光谱测量原理,在陀螺灵敏度决定过程中,腔长的上述两方面影响所起的截然相反作用,从而决定了以陀螺最佳灵敏度为判据存在最佳腔长.仿真与实验结果与理论分析良好一致性,不仅证明了理论分析的正确性.也显示了其在集成光学陀螺结构设计与性能优化中的应用价值.     

有限长柱形药包土中爆腔特征尺寸的计算方法

为计算柱形药包土中爆腔尺寸，提出了一种有限长柱形药包在土中爆炸的特征尺寸近似计算方法，该方法利用球形药包爆腔膨胀准静态模型叠加的方式，给出了长径比较大情况下柱形药包爆腔特征尺寸及塑性区半径。与数值模拟对比表明，该方法的误差随长径比的增大而减小，当球形药包数量N=n、长径比在10及以上时，误差在12.2%以内，表明该方法能够较准确地预测有限长柱形药包爆腔的特征尺寸。     

平面装药模拟高超压长持续时间爆炸作用

平面装药长持续时间爆炸作用是一种有效的试验手段，对土中平面波传播和土介质与结构相互作用的研究具有十分重要的意义。从理论计算和大型野外试验两个方面对平面装药爆炸作用过程进行了深入研究。在理论计算上，建立了ＴＮＴ装药和真实空气双介质物理模型，利用一维不定常平面流动拉格朗日型偏微分方程组求解，给出了爆腔压力空间分布和地面压力随时间的变化结果。在野外试验中，克服了平面装药模拟高超压作用的许多技术困难，实现了高达４５ＭＰａ的压力加载。文中将理论计算和野外试验结果进行了综合比较，结果表明压力波形的振荡衰减特征完全一致，从而证明了平面装药模拟高超压长持续时间爆炸作用是完全可以实现的有效而可靠的手段。     

三维内埋式弹舱流动特性及形状影响数值分析

采用多块结构化网格生成方法,基于三维雷诺平均 N-S 方程和 SST k-ω湍流模型,模拟了三维内埋式弹舱的流动特性.通过改变弹舱各壁面的倾角及宽深比,分析其对弹舱流动特性的影响,并研究了前缘锯齿对气动特性的影响.计算结果表明,在适当范围内:增大弹舱后壁倾角、减小侧壁倾角、增大宽深比可以改善弹舱流动特性;舱门前缘增加锯齿,有明显的减阻效果.     

混合流体分离比对摆动行波的影响

利用Simple算法对流体力学基本方程组进行了数值模拟,研究了摆动行波的特性。结果表明:在摆动行波存在的范围内,当相对瑞利数r较小时,腔体内沿着空间的平均波数是随着时间周期变化的;当r较大时,腔体内沿着空间的平均波数随着时间增大保持为常数。摆动行波的摆动周期随r的增大而减小;分离比负值越大,变化越平缓,摆动行波出现的范围越大。随着长高比增加,较小分离比时,摆动周期的上、下限明显提高;分离比较大时,摆动行波存在区间对应的r上、下限明显增加。     

基于动网格技术的柔性后缘自适应机翼气动特性分析

研究了带柔性后缘的可变弯度自适应机翼在自适应变弯度过程中的气动特性.自适应变弯度过程中的气动力计算采用了基于弹簧理论的非结构动网格技术,求解NS方程时采用有限体积的二阶迎风格式离散,时间推进为隐格式双时间推进方法.通过计算柔性后缘机翼的升力特性、阻力特性及升阻比特性,并与带刚性后缘机翼的气动特性进行比较,发现柔性后缘机翼在后缘偏转时,其最大升阻比对应的迎角随着偏转角增大而降低.在中等迎角及接近失速迎角情况下,柔性后缘机翼升力系数明显优于刚性后缘机翼,并且其升力线变化较为平缓,有效迎角范围更大.     

冲击作用下刚塑性拱的大变形响应分析

本文对受集中冲击作用的深圆拱的刚塑性动力响应进行了理论分析和数值计算，用瞬时构形法得到了问题的全程解，提出发生反向弯曲的必要条件和反向弯曲变形的近似分析方法，确定了反向弯曲出现的临界冲击速度范围，并讨论质量比，能量比和支承条件对结构的响应时间，塑性形区域和最变形的影响。本文理论分析结果与实验数据吻合。     

煤(岩)体埋深及倾角对压应力型冲击地压的影响研究

为了研究不同埋深和煤岩倾角对煤岩体开采引起的压应力型冲击地压过程,利用颗粒流方法进行冲击地压模拟。煤岩体与顶板倾角工况设置为7种,分别为0°,5°,10°,15°,20°,30°和40°;煤岩层埋深深度设置为8种,分别为-120m～-820m,间隔100m。在上述56种工况下进行压应力型冲击地压过程模拟,计算至平衡时统计飞石颗粒数量和变形颗粒数量,研究这两个数量与埋深和倾角之间的关系。结果表明,飞石数量和变形颗粒数量与埋深均为幂函数关系,前者随倾角增加幂次数增加且均大于1,后者随倾角增加幂次数增加但均小于1;在深度不变时两者数量增长与倾角为线性关系,且随着深度增加,斜率和截距均增加。     

NUMERICAL INVESTIGATION ON FLOW STABILITY OF RAYLEIGH-B?NARD CONVECTION OF COLD WATER IN A RECTANGULAR CAVITY COOLED AND HEATED SYMMETRICALLY RELATIVE TO THE TEMPERATURE OF DENSITY MAXIMUM

为了解具有密度极值流体瑞利-贝纳德对流特有现象和规律,利用有限容积法对长方体腔内关于密度极值温度对称加热-冷却时冷水瑞利-贝纳德对流的分岔特性进行了三维数值模拟,得到了不同条件下的对流结构型态及其分岔序列,分析了密度极值特性、瑞利数、热边界条件以及宽深比对瑞利-贝纳德对流的影响. 结果表明:具有密度极值冷水瑞利-贝纳德对流系统较常规流体更加稳定,且流动型态及其分岔序列更加复杂;相同瑞利数下多种流型可以稳定共存,各流型在相互转变中存在滞后现象;随着宽深比的增加,流动更易失稳,对流传热能力增强;系统在导热侧壁时比绝热侧壁更加稳定,对流传热能力有所减弱;基于计算结果,采用线性回归方法,得到了热壁传热关联式.     

Investigation of three-dimensional effects on a cavitating Venturi flow

A numerical investigation of the behaviour of a cavitation pocket developing along a Venturi geometry has been performed using a compressible one-fluid hybrid RANS/LES solver. The interplay between turbulence and cavitation regarding the unsteadiness and structure of the flow is complex and not well understood. This constitutes a determinant point to accurately simulate the dynamic of sheet cavities. Various turbulent approaches are tested: a new Scale-Adaptive model and the Detached Eddy Simulation. 2D and 3D simulations are compared with the experimental data. An oblique mode of the sheet is put in evidence.     

Mixing Characteristics in a Supersonic Combustor with Gaseous Fuel Injection Upstream of a Cavity Flameholder

Non-reacting experiments and numerical simulations have been performed to investigate the mixing characteristics in a supersonic combustor with gaseous fuel injection upstream of a flameholding cavity in a supersonic vitiated air flow with stream Mach number 1.7. Using helium as simulated fuel, the acetone vapor is adulterated into the fuel jet. The fuel distribution in spanwise and streamwise direction is imaged by the planar laser-induced fluorescence (PLIF) measurement. According to the similarity of experimental observations with different cavities, the typical L/D = 7 cavity with aft wall angle 45° is chosen and the flowfield with the injection is calculated by Large Eddy Simulation. Experimental and numerical results have shown that most of the fuel flow away upon the open cavity with the lifting counter-rotating vortex structures induced by the transverse jet. Only a small portion of the fuel is convected into the cavity shear layer by the vortex interaction of the jet with cavity shear layer, and then transported into the cavity due to the cavity shear layer motion and the interaction of the shear layer with the cavity trailing edge.     

Highly Resolved LES and URANS of Turbulent Buoyancy-Driven Flow Within Inclined Differentially-Heated Enclosures

In the present study, both LES and unsteady RANS computations are presented, for turbulent natural convection of air inside differentially-heated rectangular tilted cavities using a finite volume code (Code_Saturne), for an aspect ratio of H/L?=?28.6 and Rayleigh number of 0.86×10⁶. Attention is focused on two angles of inclination: 15° to the horizontal with hot lower and cold upper wall, the 15° unstable case, and the mirror image of this case where the angle is the same but with a hot upper and cold lower wall, the 15° stable case. In accordance with recent experimental data, the LES computations for both the stable and unstable tilted cavities returned three-dimensional time-averaged flow fields. In the case of the unstably stratified enclosure, the flow is highly unsteady with coherent turbulent structures in the core of the enclosure. Time-averaged temperature, velocity and resolved turbulence intensities resulting from LES computations show close agreement to measured data. Subsequent comparisons of different URANS schemes with LES are used in order to explore to what extent these models are able to reproduce the large-scale unsteady flow structures. All URANS schemes have been found to be able to reproduce the 3-D unsteady flow features present in the 15° unstable cavity. However, the low-Reynolds-number model tested, as well as requiring a high resolution near-wall grid, also needed a finer grid in the core region than the high-Reynolds-number models, thus making it computationally very expensive. Flow within the 15° stable cavity also shows some 3-D features, although it is significantly less unsteady, and the URANS models tested here have been less successful in reproducing this flow pattern. The overall heat transfer is presented here for both differentially heated enclosures.     

Characteristics of Oscillations in Supersonic Open Cavity Flows

Characteristics of oscillations in supersonic open cavity flows are investigated numerically using hybrid RANS/LES (Reynolds-Averaged Navier-Stokes/Large Eddy Simulation) method. The oscillation regimes and feedback mechanisms for the supersonic cavity flows are identified and analyzed. The calculation captures a mixed shear-layer/wake oscillation mode in the flow of Ma = 1.75, where these two modes occur alternately. The shear-layer mode and wake mode are driven by vortex convection-acoustic feedback and absolute instability, respectively. In particular, the results indicate that the feedback-acoustic-wave in the shear-layer mode is probably generated by the reflection of the downstream-traveling pressure wave, associated with the shed vortex in the shear layer, on the aft wall. The cavity flow of Ma = 2.52 is then simulated to see the influence of Mach number. It is found that the increase of Mach number may decrease the amplitude of the fluctuations in the shear layer, inhibiting the transition to wake mode. Furthermore, the influence of upstream injection is also studied, where the results show that the injection only weakens the oscillations and faintly shifts the resonant frequencies.     

用基于M-SST模型的DES数值模拟喷流流场

脱体涡数值模拟方法(dettached eddy simulation,DES)是把雷诺平均Navier-Stokes方程(RANS)方法及大涡模拟方法(LES)结合起来模拟有脱体涡的湍流流场的数值模拟方法,其主要思想是在物面附近解雷诺平均Navier-Stokes方程、在其他区域采用Smagorinski大涡模拟方法。本文在剪切应力传输(SST)湍流模型的基础上用DES及混合非结构网格数值模拟具有横向喷流的湍流流场,算法采用Osher逆风格式,利用该套程序(包括网格生成及算法),对导弹在不同马赫数下的喷流流场进行了数值模拟,并与同时开展的实验研究的结果进行了对比,结果表明用该方法处理这类问题是较准确的。     

Prediction of separation flows around a 6：1 prolate spheroid using RANS/LES hybrid approaches

This paper presents hybrid Reynolds-averaged Navier–Stokes (RANS) and large-eddy-simulation (LES) methods for the separated flows at high angles of attack around a 6:1 prolate spheroid. The RANS/LES hybrid methods studied in this work include the detached eddy simulation (DES) based on Spalart–Allmaras (S–A), Menter’s k–ω shear-stress-transport (SST) and k–ω with weakly nonlinear eddy viscosity formulation (Wilcox–Durbin+, WD+) models and the zonal-RANS/LES methods based on the SST and WD+ models. The switch from RANS near the wall to LES in the core flow region is smooth through the implementation of a flow-dependent blending function for the zonal hybrid method. All the hybrid methods are designed to have a RANS mode for the attached flows and have a LES behavior for the separated flows. The main objective of this paper is to apply the hybrid methods for the high Reynolds number separated flows around prolate spheroid at high-incidences. A fourth-order central scheme with fourth-order artificial viscosity is applied for spatial differencing. The fully implicit lower–upper symmetric-Gauss–Seidel with pseudo time sub-iteration is taken as the temporal differentiation. Comparisons with available measurements are carried out for pressure distribution, skin friction, and profiles of velocity, etc. Reasonable agreement with the experiments, accounting for the effect on grids and fundamental turbulence models, is obtained for the separation flows. The project supported by the National Natural Science Foundation of China (10502030 and 90505005).     

Large eddy simulation of fuel injection and mixing process in a diesel engine

The large eddy simulation(LES) approach implemented in the KIVA-3V code and based on one-equation sub-grid turbulent kinetic energy model are employed for numerical computation of diesel sprays in a constant volume vessel and in a Caterpillar 3400 series diesel engine.Computational results are compared with those obtained by an RANS(RNG k-ε) model as well as with experimental data.The sensitivity of the LES results to mesh resolution is also discussed.The results show that LES generally provides flow and spray characteristics in better agreement with experimental data than RANS;and that small-scale random vortical structures of the in-cylinder turbulent spray field can be captured by LES.Furthermore,the penetrations of fuel droplets and vapors calculated by LES are larger than the RANS result,and the sub-grid turbulent kinetic energy and sub-grid turbulent viscosity provided by the LES model are evidently less than those calculated by the RANS model.Finally,it is found that the initial swirl significantly affects the spray penetration and the distribution of fuel vapor within the combustion chamber.     

Log-layer mismatch and commutation error in hybrid RANS/LES simulation of channel flow

Hybrid approach combining large eddy simulation (LES) with the Reynolds-averaged Navier–Stokes equation (RANS) is expected to accurately simulate wall-bounded turbulent flows at high Reynolds numbers. As an important issue in developing hybrid methods, it is known that the log layers in the RANS and LES regions are not lined up in hybrid RANS/LES simulations of channel flow. Although several methods including additional filtering near the RANS/LES interface have been proposed to eliminate the log-layer mismatch, there is no obvious physical justification for the methods and some ad hoc tuning is necessary. In this work, the commutation error terms in the filtered velocity equations are investigated to justify the method of additional filtering. It is shown that the additional filtering can be considered as a finite difference approximation to extra terms due to the non-commutivity between the hybrid filter and the spatial derivative. Moreover, an expression determining the filter width and its location for the additional filtering is obtained. To validate the expression, a hybrid simulation of channel flow is carried out. The additional filtering with the filter width derived is shown to be effective in eliminating the log-layer mismatch and improving the mean velocity profile.     

Hybrid LES–RANS technique based on a one-equation near-wall model

In order to reduce the high computational effort of wall-resolved large-eddy simulations (LES), the present paper suggests a hybrid LES–RANS approach which splits up the simulation into a near-wall RANS part and an outer LES part. Generally, RANS is adequate for attached boundary layers requiring reasonable CPU-time and memory, where LES can also be applied but demands extremely large resources. Contrarily, RANS often fails in flows with massive separation or large-scale vortical structures. Here, LES is without a doubt the best choice. The basic concept of hybrid methods is to combine the advantages of both approaches yielding a prediction method, which, on the one hand, assures reliable results for complex turbulent flows, including large-scale flow phenomena and massive separation, but, on the other hand, consumes much fewer resources than LES, especially for high Reynolds number flows encountered in technical applications. In the present study, a non-zonal hybrid technique is considered (according to the signification retained by the authors concerning the terms zonal and non-zonal), which leads to an approach where the suitable simulation technique is chosen more or less automatically. For this purpose the hybrid approach proposed relies on a unique modeling concept. In the LES mode a subgrid-scale model based on a one-equation model for the subgrid-scale turbulent kinetic energy is applied, where the length scale is defined by the filter width. For the viscosity-affected near-wall RANS mode the one-equation model proposed by Rodi et al. (J Fluids Eng 115:196–205, 1993) is used, which is based on the wall-normal velocity fluctuations as the velocity scale and algebraic relations for the length scales. Although the idea of combined LES–RANS methods is not new, a variety of open questions still has to be answered. This includes, in particular, the demand for appropriate coupling techniques between LES and RANS, adaptive control mechanisms, and proper subgrid-scale and RANS models. Here, in addition to the study on the behavior of the suggested hybrid LES–RANS approach, special emphasis is put on the investigation of suitable interface criteria and the adjustment of the RANS model. To investigate these issues, two different test cases are considered. Besides the standard plane channel flow test case, the flow over a periodic arrangement of hills is studied in detail. This test case includes a pressure-induced flow separation and subsequent reattachment. In comparison with a wall-resolved LES prediction encouraging results are achieved.