k-ε双方程湍流模型对制退机内流场计算的适用性分析

为研究不同双方程湍流模型对制退机内复杂流场计算的适用性,以某火炮制退机为研究对象,建立了实际结构下的三维计算模型,利用动网格与滑移网格技术,实现了火炮实际后坐速度下的制退机内部三维运动流场的数值计算。分别采用标准k-ε模型、RNGk-ε模型和Realizable k-ε模型计算制退机内部各腔室压力,与实验曲线对比,结果表明,应用标准k-ε模型对后坐冲击过程的制退机内部压力计算的误差最小,与实验结果吻合最好。     

火炮制退机流场的数值模拟

利用Rodi ASM模型和SIMPLER计算程式对接近于实际结构的制退机流场做了数值模拟,获得了制退机湍流流场的结构,并对数值模拟结果进行了分析,数值模拟结果揭示了制退机流场中存在流动截面收缩现象。     

离轴受限圆射流的三维湍流流场分析

利用原始变量,交叉网格的方法,以及采用κ-ε湍流模型,计算了离轴受限射流三维湍流流场,计算结果与用LDV测量结果作了比较,结果是满意的。计算结果提供了这样一个复杂的三维流场的清晰图象,为研究具有离轴射流的预燃室的稳燃机理提供了很多有用的信息,对计算结果作了详细的分析和讨论后,指出这样一个流场的重要三维特征是在θ方向梳线有大角度的扭曲。     

非正交网格控制容积法分析复杂形状池内的流动

介绍一种基于非正交网格控制容积法的数学模型，及其在圆形沉沙池流动研究中的应用．该模型求解轴对称流动的连续方程及时均Ｎ－Ｓ方程，并采用标准ｋ－ε紊流模型，模拟圆形池内的流动．由于采用非正交网格，此计算模型可精确模拟几何形状较复杂的沉沙池内的流动，利用上述模型对某实际沉沙池进行了流场计算，计算所得流场与模型试验实测值符合良好．     

非正交网格控制容积法分析复杂形状池内的流动

介绍一种基于非正交网格控制容积法的数学模型，及其在圆形沉沙池流动研究中的应用．该模型求解轴对称流动的连续方程及时均Ｎ－Ｓ方程，并采用标准ｋ－ε紊流模型，模拟圆形池内的流动．由于采用非正交网格，此计算模型可精确模拟几何形状较复杂的沉沙池内的流动，利用上述模型对某实际沉沙池进行了流场计算，计算所得流场与模型试验实测值符合良好．     

基于分区拼接网格技术高升力装置流场数值模拟

针对高升力装置构型模型结构复杂、流场变化剧烈等特点,本文采用分区拼接网格的思想分别按照流场和结构拓扑特点对高升力装置进行了网格分区。在分区的基础上逐块生成网格,减小了增升装置网格的生成难度,提高了网格质量,减少了网格数目。首先,研究了高升力装置的分区策略及流场特点;接着以MD30P-30N多段翼型为研究算例研究了网格比例和插值方法对计算结果的影响。经过分析对比可知:外部区域与近壁区域之间的比例不宜小于1:5;内部域网格比例不宜超过1:1.8,最好保持在1:1左右;计算中应该采用高阶精度插值以保证计算精度。采用某四段翼型进行了验证;最后采用NASA标准高升力装置进行了三维高升力装置流场数值模拟并与相应风洞实验数据及对接网格计算结果进行了比较与分析,验证了拼接网格技术的高效性与可靠性。同时分析研究了绕三维增升装置的流动及其周围复杂的粘性流动现象。     

三维复杂流场内能激发松弛的DSMC数值模拟

本文研究了带内能激发松弛和不同来流条件下的复杂外形高超声速飞行器过渡区三维流场特性.采用一种对飞行器物面网格与DSMC计算域网格分别标识的方法,通过判断模拟分子与表面碰撞来完成飞行器物面网格与DSMC计算域网格间的信息传递和信息存贮,对于复杂外形飞行器精确描述的物面网格不需做进一步处理,直接应用于不依赖于飞行器外形的DSMC计算的通用子程序中.采用L-B碰撞模型,分析了内能激发松弛的流场参数特性,从分子碰撞、分子与飞行器表面碰撞的角度,分析了不同来流条件下的流场参数的变化情况.     

带制退器膛口流场的三维数值模拟

利用三维任意坐标 系下Euler方程组和LU分解的数值方法模拟了不带膛口装置的某步枪膛口流场,计算结果与 实验结果比较吻合. 对某口径火炮带有制退器的膛口流场进行了数值模拟,其结果反映了制 退器内外流场的主要特征和冲击波波系之间的相互作用,在数值模拟的基础上分析 了3排侧 孔对制退器性能的影响. 本文提供了有关膛口流场的详细数据,将对制退器总体性能研究有 重要的参考意义.     

基于多孔介质模型的养殖网箱周围流场特性研究

为探讨鱼类养殖网箱对周围海域流场的影响,基于雷诺时均方程和Realizable k-ε湍流模型,耦合多孔介质模型,建立了研究鱼类养殖网箱周围流场的三维水动力模型.根据典型养殖场的实际数据,分别对仅网衣作用、仅鱼类作用及网衣和鱼类共同作用下网箱周围流场进行数值模拟.结果表明综合考虑网衣和鱼类的作用时,网箱下游流速衰减最大...     

非结构混合网格高超声速绕流与磁场干扰数值模拟

对均匀磁场干扰下的二维钝头体无粘高超声速流场进行了基于非结构混合网格的数值模拟.受磁流体力学方程组高度非线性的影响及考虑到数值模拟格式的精度,目前在此类流场的数值模拟中大多使用结构网格及有限差分方法,因而在三维复杂外形及复杂流场方面的研究受到限制.本文主要探索使用非结构网格(含混合网格)技术时的数值模拟方法.控制方程为耦合了Maxwell方程及无粘流体力学方程的磁流体力学方程组,数值离散格式采用Jameson有限体积格心格式,5步Runge-Kutta显式时间推进.计算模型为二维钝头体,初始磁场均匀分布.对不同磁感应强度影响下的高超声速流场进行了数值模拟,并与有限的资料进行了对比,得到了较符合的结果.     

古树包滑坡滑带土蠕变经验模型

针对清江古树包滑坡滑带土的室内蠕变试验资料,结合经验模型理论,以双曲线函数描述蠕变应力-应变关系,并引入Morgan Mercer Florin函数,用来描述应变-时间关系,建立了一种新的蠕变经验模型。采用ε/ -ε关系曲线和遗传算法求解出了该模型的所有参数。经过计算,所建经验模型的计算结果与实际实验数据的相关系数在0.96以上,说明所建经验模型能有效地反映古树包滑坡滑带土的蠕变特性。     

Domain perturbation method and shape of a bubble in a uniform flow of an inviscid liquid

In many problems with a free boundary there is defined a small parameter, , for which the solution is sometimes known for a particular value, =0, and the general solution is obtained as a series in the parameter. To find this solution, the equations can be written on a reference configuration and solved in a fixed domain. The purpose of this study is to show that this method of domain perturbation is a good one. The range of validity of this method will be studied on the model example of the irrotational flow of a perfect fluid around a bubble. The radius of convergence of the series solution will be determined, as will the nature of the solution in the neighbourhood of the first real singularity.     

2:1:1 Resonance in the Quasi-Periodic Mathieu Equation

We present a small $ \ddot x + (\delta + \epsilon \cos t + \epsilon \cos \omega t) x=0 $ \ddot x + (\delta + \epsilon \cos t + \epsilon \cos \omega t) x=0      

CFD-simulation of a circulating fluidized bed riser

In the current work, a model of the fluid mechanics in the riser of a circulating fluidized bed （CFB） has been implemented using computational fluid dynamics （CFD）. The model developed shall be used in future as the basis of 3D-reactor model for the simulation of large scale CFB combustors. The two-fluid model （TFM） approach is used to represent the fluid mechanics involved in the flow. The computational implementation is accomplished by the commercial software FLUENT. Different closure formulations are tested on a simplified geometry. Two different turbulence formulations, namely the swirl modified RNG k-e model and the Realizable k-e model, are tested in combination with two different approaches to solid phase turbulence, namely the dispersion and per phase approach. One focus of the current work is put on the study of different drag correlations. Besides the drag correlations by Syamlal et al. [Syamlal, M., Rogers, W., ＆ O＇Brien, T. J. （1993）. MFIX documentation theory guide. Technical Report DOE/METC-9411004, U.S. Department of Energy （DOE）. Morgantown Energy Technology Center： Morgantown, WV] and Gidaspow [Gidaspow, D. （1994）. Multiphaseflow andfluidization. New York： Academic Press] the EMMS model has been used to determine the momentum exchange between the two phases. The resulting formulation is then used to simulate a 1-m × 0.3-m cold CFB setup and is validated by experimental results [Schlichtharle, P. （2000）. Fluid dynamics and mixing of solids and gas in the bottom zone of circulating fluidized beds. Unoublished doctoral dissertation, Technische Universitaet Hamburg-Harburg, Shaker Verlag： Aachen].     

Incompressible 3D Navier–Stokes Equations as a Limit of a Nonlinear Parabolic System

In this paper, we consider the Cauchy problem for a nonlinear parabolic system ${u^\epsilon_t - \Delta u^\epsilon + u^\epsilon \cdot \nabla u^\epsilon + \frac{1}{2}u^\epsilon\, {\rm div}\, u^\epsilon - \frac{1}{\epsilon}\nabla\, {\rm div}\, u^\epsilon = 0}$ in ${\mathbb {R}^3 \times (0,\infty)}$ with initial data in Lebesgue spaces ${L^2(\mathbb {R}^3)}$ or ${L^3(\mathbb {R}^3)}$ . We analyze the convergence of its solutions to a solution of the incompressible Navier?CStokes system as ${\epsilon \to 0}$ .     

Effect of dislocation structure on fracture toughness of strained BCC-metals

Mechanical properties of materials depend on their structure. Examined are the effects of dislocation structure on fracture toughness and mechanisms of fracture of BCC-metals. Fracture toughness was determined by depending specimens with cracks introduced into the plane perpendicular to the plane of rolling. Fracture toughness increases with decreasing yield stress (for =15–25%). This is due to instability of slightly misoriented cell structure under repeated loading. The peak of fracture toughness at the temperature 77 K was not observed. The increase of fracture toughness for high strained metals (>60% for Mo, and >85% for Cr) corresponds to cell size reduction and the change of fracture mechanisms.     

On the Darcy–Brinkman Flow Through a Channel with Slightly Perturbed Boundary

The goal of this paper is to study the effects of a slightly perturbed boundary on the Darcy–Brinkman flow through a porous channel. We start from a rectangular domain and then perturb the upper part of its boundary by the product of the small parameter \(\epsilon \) and arbitrary smooth function h. Using asymptotic analysis with respect to \(\epsilon \), the effective model has been formally derived. Being in the form of the explicit formulae for the velocity and pressure, the asymptotic approximation clearly shows the nonlocal effects of the small boundary perturbation. The error analysis is also conducted providing the order of accuracy of the asymptotic solution.     

马兰黄土的应力-应变本构方程

依据马兰黄土的室内CU试验,运用土力学的相关理论,建议性地提出了马兰黄土的应力-应变本构方程。通过对不同围压下的试验结果(应变小于15%)与理论模型的拟合发现,二者吻合较好。另外,与以往的模型相比,本文的模型参数a和b只需通过一次常规三轴试验即可全部求得,具有参数少、容易确定的优点。     

亚微米级型CL-20的制备、表征与性能

采用超临界气体抗溶剂（gas anti-solvent, GAS）技术制备了平均粒径为721.9 nm的亚微米级型六硝基六氮杂异伍兹烷（CL-20）,应用傅立叶变换红外光谱（Fourier translation infrared spectrum, FT-IR）,对亚微米级CL-20的晶型进行了鉴别,并进行了小隔板试验、撞击感度和爆发点测试。研究结果表明,细化后的亚微米级型CL-20冲击波感度的隔板厚度降低了58.6%,撞击感度的特性落高提高了84.1%,而爆发点保持不变。     

Progress in subgrid-scale transport modelling for continuous hybrid non-zonal RANS/LES simulations

The partially integrated transport modelling (PITM) method can be viewed as a continuous approach for hybrid RANS/LES modelling allowing seamless coupling between the RANS and the LES regions. The subgrid turbulence quantities are thus calculated from spectral equations depending on the varying spectral cutoff location [Schiestel, R., Dejoan, A., 2005. Towards a new partially integrated transport model for coarse grid and unsteady turbulent flow simulations. Theoretical and Computational Fluid Dynamics 18, 443–468; Chaouat, B., Schiestel, R., 2005. A new partially integrated transport model for subgrid-scale stresses and dissipation rate for turbulent developing flows. Physics of Fluids, 17 (6)] The PITM method can be applied to almost all statistical models to derive its hybrid LES counterpart. In the present work, the PITM version based on the transport equations for the turbulent Reynolds stresses together with the dissipation transport rate equation is now developed in a general formulation based on a new accurate energy spectrum function E(κ) valid in both large and small eddy ranges that allows to calibrate more precisely the c_sgs2 function involved in the subgrid dissipation rate _sgs transport equation. The model is also proposed here in an extended form which remains valid in low Reynolds number turbulent flows. This is achieved by considering a characteristic turbulence length-scale based on the total turbulent energy and the total dissipation rate taking into account the subgrid and resolved parts of the dissipation rate. These improvements allow to consider a large range of flows including various free flows as well as bounded flows. The present model is first tested on the decay of homogeneous isotropic turbulence by referring to the well known experiment of Comte-Bellot and Corrsin. Then, initial perturbed spectra E(κ) with a peak or a defect of energy are considered for analysing the model capabilities in strong non-equilibrium flow situations. The second test case is the classical fully turbulent channel flow that allows to assess the performance of the model in non-homogeneous flows characterised by important anisotropy effects. Different simulations are performed on coarse and refined meshes for checking the grid independence of solutions as well as the consistency of the subgrid-scale model when the filter width is changed. A special attention is devoted to the sharing out of the energy between the subgrid-scales and the resolved scales. Both the mean velocity and the turbulent stress computations are compared with data from direct numerical simulations.