高温真实气体底部流动的NS方程数值计算

本文数值模拟了高超音速飞行时钝锥的底部流动。采用轴对称NS方程并考虑真实气体效应。湍流模型采用修正的Baldwin-Lomax涡粘性代数模型,数值方法空间离散对流项采用显式NND格式,粘性项采用中心差。时间离散采用三阶的龙格-库塔法。真实气体模型采用考虑七种组分四种反应的汉森模型。给出了底部流场的压力和温度分布及各组分的浓度分布。可以看出在近底部区域高速流-绕过拐角就产生一回流旋涡区。由于温度变化很大,气体的热力学特性受气体离解、复合和振动能激发的影响。所以整个流动过程变得十分复杂。     

风压作用下房屋自然通风的PIV实验研究

为认识屋内流动结构和改善房屋自然通风效果,本文介绍了风压作用下大开口的房屋自然通风模型的风洞实验结果.首先通过表面压力分布测量得到了建筑物表面风压系数的分布规律,以此为依据选取合理的窗口位置.而后采用PIV粒子图像测速技术,对房屋模型内自然通风的流场进行定量测量,研究不同窗口大小和窗口位置对室内流场速度分布的影响;通过...     

平屋盖风压分布的数值模拟

基于Reynolds时均N-S方程和RSM模型对平屋面的风压分布进行了数值模拟,在此基础上系统研究了风向角、跨高比、地面粗糙度、风速等因素对屋面风压分布的影响,探讨了结构周围流场的绕流特性,最后根据屋面的结构形式及风压分布特点将屋面进行分区,给出了屋面在不同风向角下的分区风载体型系数以供工程设计参考。     

Equation of state calculations for hot,dense matter at arbitrary densities and temperatures

Within the approximations of spherical lattice cell, central-field, and relativistic Fermi statistics, an algorithm with average atom model is presented to calculate the electronic energy levels and equation of state for hot and dense matter at arbitrary densities and temperatures. Choosing Zink's analytical potential as initial potential, we have solved the Dirac-Slater equation which satisfies the Weigner-Seitz boundary condition. The electronic energy bands are not taken into account. Taking energy level degeneracy as a continuous function of density, we have considered the pressure ionization effects for highly dense matter. Results for¹³Al atom are shown.     

基于雷诺应力模型的平衡大气边界层模拟

基于雷诺应力湍流模型（简称RSM模型）,研究了平衡大气边界层风场数值模拟问题．假设流体不可压,且不计雷诺应力输运方程中的对流项、浮力产生项、系统旋转产生项和扩散项,在准各向同性的条件下,推导出RSM模型湍动能k的表达式是标准k-ε模型k常数表达式的0.893倍．考虑k沿高度变化的修正,根据在标准k-ε模型中满足水平均匀性的湍流来流边界条件,提出在RSM模型中产生平衡大气边界层的湍流来流边界条件．基于空风洞的数值模拟结果表明,与工程上常用的湍流来流边界条件相比,基于本文提出的湍流来流边界条件得到的风场水平均匀性更优,且在整个流域内,得到的雷诺应力剖面更合适．从而验证了该湍流来流边界条件的适用性．     

低层房屋风压分布特性的模拟与分析

基于Reynolds时均N—S方程和标准κ-ε双方程封闭模型,对一类带挑檐的低层双坡屋面房屋在不同风向角下的表面风压进行了敷值模拟,模拟结果与风洞试验结果作了比较。在此基础上,根据模拟结果对房屋各表面及挑檐处的风压分布特性和规律进行了分析和计算,获得了不同风向角下各表面的风载体型系数,为该类房屋的工程抗风设计提供了依据。     

石英光纤探针在非金属材料冲击实验中的应用

使用石英光纤探针成功测量了Teflon样品中的冲击波速度和电炮加载的厚度为75 m的Mylar膜飞片的平均速度和到达时间一致性。实验证明,石英光纤探针技术为非金属材料的冲击实验提供了一种直接的、不受电磁干扰的、精确的测量手段。     

各向异性本构关系在板料成形数值模拟中的应用

对几种能表达面内各向异性的屈服准则Hill、Barlat-Lian、Barlat进行了比较。以弹性变形服从各向同性广义虎克定律的情况下，给出了基于张量算法推导的弹塑性本构关系的一般表达式，并由此导出了相应屈服准则的弹塑性本构关系的显式表达。借助ABAQUS软件本构模块用户子程序接口，分别实现了这些屈服准则在ABAQUS的嵌入。以模拟方形盒的拉延过程为例，分析了不同的屈服准则在板料成形过程数值模拟中的应用。模拟结果表明，基于弹塑性本构关系一般表达所列出的相应屈服准则的显式表达式是正确的；在采用壳元来模拟板料成形时，采用Barlat准则的模拟结果和采用Barlat-Lian准则的结果差别不大。     

Analysis of velocity equation of steady flow of a viscous incompressible fluid in channel with porous walls

Steady flow of a viscous incompressible fluid in a channel, driven by suction or injection of the fluid through the channel walls, is investigated. The velocity equation of this problem is reduced to nonlinear ordinary differential equation with two boundary conditions by appropriate transformation and convert the two‐point boundary‐value problem for the similarity function into an initial‐value problem in which the position of the upper channel. Then obtained differential equation is solved analytically using differential transformation method and compare with He's variational iteration method and numerical solution. These methods can be easily extended to other linear and nonlinear equations and so can be found widely applicable in engineering and sciences. Copyright © 2009 John Wiley & Sons, Ltd.     

A boundary integral equation method for Navier-Stokes equations—application to flow in annulus of eccentric cylinders

A novel Navier-Stokes solver based on the boundary integral equation method is presented. The solver can be used to obtain flow solutions in arbitrary 2D geometries with modest computational effort. The vorticity transport equation is modelled as a modified Helmholtz equation with the wave number dependent on the flow Reynolds number. The non-linear inertial terms partly manifest themselves as volume vorticity sources which are computed iteratively by tracking flow trajectories. The integral equation representations of the Helmholtz equation for vorticity and Poisson equation for streamfunction are solved directly for the unknown vorticity boundary conditions. Rapid computation of the flow and vorticity field in the volume at each iteration level is achieved by precomputing the influence coefficient matrices. The pressure field can be extracted from the converged streamfunction and vorticity fields. The solver is validated by considering flow in a converging channel (Hamel flow). The solver is then applied to flow in the annulus of eccentric cylinders. Results are presented for various Reynolds numbers and compared with the literature.