基于大涡模拟的移动粒子半隐式法研究

将大涡模拟(LES)和无网格的移动粒子半隐式法(MPS)相结合, 以求解湍流中的自由表面问题. 对N-S方程进行滤波计算可得到大涡模拟的控制方程, 大涡模拟的控制方程相对于以往的移动粒子半隐式法而言仅多出雷诺应力项, 通过亚粒子应力(sub-particle-scale,SPS)模型并引入Smagorinsky涡黏模型将雷诺应力模型化, 可实现移动粒子半隐式法的大涡模拟. 将MPS-LES应用至具有大变形自由表面的共振晃荡中, 其模拟结果同实验及其他数值模拟结果都相当接近.      

超薄膜磁头滑块气动力特性

采用有限差分法对广义润滑方程进行数值求解,计算出计算机磁头滑块压强场的分布情况。分析、研究了其稳态和动态气动力特性,并将计算结果分别与求解一阶、二阶修正雷诺方程所得到的结果进行了比较,得到如下三个结论:(1)当飞行高度很小,飞行速度较低时,必须采用广义润滑方程进行磁头滑块的气动力计算,(2)与广义润滑方程结果比较,求解一阶修正雷诺方程所得到的计算结果总是偏高,而求解二阶修正雷诺方程所得到的计算结果总是偏低。此外,还解决了大压缩数下数值失稳问题,使得压缩数可以计算到120万,足以适应任何实际工程的需要。     

分块隐式有限差分法计算弯管紊流

本文利用贴体坐标系中分块隐式有限差分法计算矩形截面90°弯管中不可压恒定紊流.在计算中,雷诺方程的数值离散采用混合差分格式,局部联立求解雷诺方程和连续方程而得到速度压力解.在全流场的迭求解过程中采用对称联立Gauss—Seidel法.利用标准K-ε紊流模型模拟紊流.计算结果与有关试验进行了对比.     

层状两相流动的基本方程式

本文将两相流动的雷诺方程沿深度方向进行平均,出求了平均后的运动方程、连续方程和能量方程式,在适当的物理假定下,简化得到可适用于考虑固-液流动及气-液流动两种情况下的两相层状流动基本方程.     

有限长圆柱轴承非稳态油膜力的级数解法

本文用分离变量法求解雷诺方程，在π油膜的假设下，求得雷诺方程应满足的特征值与用傅立叶级数表达的特征函数，进而求得有限长轴的非特急油膜力解析表达式。为分析轴承转子系统的非线性动力特性提供了帮助。     

基于S-A湍流模型的Runge-Kutta有限元算法

采用一方程S-A模型(Spalart-Allmaras模型)封闭雷诺时均N-S方程(RANS方程)进行湍流数值计算,可以减少方程求解数量,节约计算时间。本文对其进行了有限元数值算法研究,首先通过沿流线坐标变换,得到无对流项RANS方程,并引入三阶Runge-Kutta法对其进行时间离散;然后利用沿流线的Taylor展开解决坐标变换带来的网格更新的困难;最后采用Galerkin法进行空间离散,得到湍流模型的有限元算法。基于方柱绕流和覆冰输电线绕流模型,与试验结果进行对比,验证了该算法的有效性,与一阶数值算法相比,该算法在精度和收敛性方面更具优势。     

积分方程的加权残数配点法

对于用积分方程式表示的物理问题,本文应用加权残数的配点法进行求解,使得求解方法简便、可靠、正确。与常规的求解微分方程的加权残数配点法相比较,本方法不需假定任何试函数。文中给出较多的算例,论证了本文提出的积分方程的加权残数配点法的有效性,可靠性和精确性。     

磁头/盘界面超薄气膜挤压效应和动压效应研究

针对求解磁头/磁盘界面动态气膜力时出现的不易收敛和编程复杂的问题,本文提出了基于PDE工具求解气体润滑的瞬态和稳态雷诺方程的方法,计算了具有不同最小气膜厚度的Tri-pad正压型和Tri-pad负压型浮动块空气轴承在加载和卸载过程中的瞬态和稳态气膜压力分布,求解了作用面上的轴承力,并将2种状态下超薄气膜挤压效应和动压效应对轴承力的影响进行了对比.仿真结果表明:该方法具有足够的求解精度,且收敛速度快,为研究具有复杂磁头形貌特征的超薄气体润滑的动态特性提供了方便、准确的方法.     

NORMAL GRID RULE FOR HYPERSONIC HEAT FLUX NUMERICAL SIMULATION

采用数值求解雷诺平均Navier-Stokes 方程方法对高超声速钝双锥模型进行了计算,比较了不同法向网格尺度对热流计算的影响,并在此基础上提出了一种适于热流计算的网格准则. 研究表明:在实际外形的热流计算中,使用过密或过粗的法向网格均不利于得到准确的热流数据;采用该文提出的基于当地温度梯度的法向网格尺度准则可以在保证热流计算精度的同时实现较高的收敛速度.     

强非线性动力系统的频率增量法

提出一类强非线性动力系统的暧时频率增量法,将描述动力系统的二阶常微分方程,化为以相位为自变量、瞬廛频率为未知函数的积分方程;用谐波平衡原理,将求解瞬时频率的积分问题,归结为求解以频率增量的Fourier系数为独立变量的线性代数方程组;给出了若干例子。     

圆形可倾瓦推力轴承入油边界条件的确定

采用有限元法对Reynolds方程、能量方程及弹性方程进行了联立求解,通过计算瓦边法向平均速度,较精确地确定了瓦的入油边界条件,从而更精确地计算出圆形瓦的压力分布及温度分布,通过轴承试验台试验发现,计算结果和实验结果较吻合。     

Generalized Reynolds and integrated energy equations for fluid-film lubrication

A generalized form of Reynolds equation governing pressure distribution and an integrated form of energy equation in the fluid-film between porous bounding walls has been developed from fundamental equations of thermohydrodynamics with a minimum number of restricting assumptions and slip conditions. In particular cases these reduce to various forms which have been found to date. It is suggested that the derived equations will find applications in discussing the various models for synovial joints and porous bearings.     

A flexible rod in annular leakage flow: Influence of turbulence and equilibrium offset,and analysis of instability mechanisms

The linear stability of a flexible, cylindrical rod subjected to annular leakage flow is studied. The mathematical models developed by Li, Kaneko, and Hayama in 2002 and Fujita and Shintani in 2001 are bridged and extended, to account for a flexible rod with equilibrium offset (eccentricity) in laminar or turbulent leakage flow. Stability characteristics are analyzed numerically for a variety of configurations. It is found that simply supported rods may become unstable at a certain critical flow speed by either divergence or flutter, depending on dimensions and fluid/solid properties. It is furthermore found that the critical flow speed is quite insensitive to use of a laminar friction model at high Reynolds numbers in cases of divergence, but sensitive to it in cases of flutter. These findings are verified analytically though analysis of an energy equation. This equation shows that (i) divergence instability is independent of fluid friction; (ii) flutter instability is caused solely by fluid friction. It also suggests a possible explanation to the question of why a ‘wrong’ fluid friction assumption gives a too large critical flow speed in cases of flutter instability at a high Reynolds number.     

Unsteady laminar separated flow through constricted channel

In the present paper unsteady Navier-Stokes equations have been solved numerically by finite-difference technique in staggered grid distribution for a flow through a channel with locally symmetric and asymmetric constrictions. A coordinate stretching has been made to map the infinite irregular geometry into a finite regular computational domain. Pressure and pressure-velocity corrections scheme have been developed. Convergence criteria (in terms of continuity equation) has been achieved after few time iterations. The critical Reynolds number for asymmetric flow through a symmetric constriction has been found. Critical values depend on the area reduction and the length of the constriction. The increment of Reynolds number grows the asymmetry of the flow. The root mean square (r.m.s.) centreline vertical velocity for asymmetric flow through a symmetric constriction has been drawn at different Reynolds numbers. For flow through symmetric constriction the centreline vertical velocity shows finite oscillation behind the constriction at high Reynolds number.     

非牛顿体通用模型线接触弹流润滑的数值分析

基于润滑剂在弹流润滑状态下表现为非牛顿体特性,根据弹流润滑理论,采用一种新的非牛顿体流变模型,建立了适用于非牛顿体的修正Ｒｅｙｎｏｌｄｓ方程,进行了等温弹流润滑的数值计算,并在等温解的基础地温度场分析。数值分析结果表明,由于滑滚比和模型参数对剪应力影响较大,因而在滑滚比和模型参数较大时应进行热弹流计算。通过温度场分析可证明：非牛顿体通过模型可用于等温弹流润滑和热弹流润滑计算。     

Development of a turbulent near-wall temperature model and its application to channel flow

A numerical study of fluid flow and heat transfer in a two-dimensional channel under fully developed turbulent conditions is reported. A computer program which is capable of treating both forced and natural convection problems under turbulent conditions has been developed. The code uses the high-Reynolds-number form of the two equation turbulent model(k-?) in which a turbulent kinetic energy near-wall model is incorporated in order to accurately represent the behavior of the flow near the wall, particularly in the viscous sublayer where the turbulent Reynolds number is small. A near-wall temperature model has been developed and incorporated into the energy equation to allow accurate prediction of the temperature distribution near the wall and, therefore, accurate calculation of heat transfer coefficients. The sensitivity of the prediction of flow and heat transfer to variations in the coefficients used in the turbulence model is investigated. The predictions of the model are compared to available experimental and theoretical results; good agreement is obtained. The inclusion of the near-wall temperature model has further improved the predictions of the temperature profile and heat transfer coefficient. The results indicate that the turbulent kinetic energy Prandtl number should be a function of Reynolds number.     

Numerical methods for theoretical analysis of foil bearing dynamics

The bearing under consideration is a self-acting, gas bearing with a flexibly supported foil. In order to conduct the theoretical analysis, three systems were identified: a rotor, a gas film and a flexible structure. Each of these elements has been mathematically analysed and analytical equations have been formulated. It was found, that the gas flow in the bearing can be described by the Reynolds equation, whereas a spring-damper model was selected for the structural analysis. The Reynolds equation is a differential equation the exact solution to which is unknown. The work describes the finite difference method in detail, where the partial derivatives in the Reynolds equation are replaced by a system of algebraic equations. In order to solve the resulting system, the Alternating Direction Implicit method (ADI) was used. Based on that and on the analysis of interactions between the three systems the equations have been computed. The accuracy of methods has been verified by means of series of numerical tests. Chosen results are described in this paper.     

Modification of an implicit approach based on nonstandard rules for structural dynamics analysis

The present study aims to modify a recently suggested implicit approach consisted of the approximate Euler method and closed-form exponential mapping (herein referred to as the Liu scheme) for the dynamic analysis of structures. Such modification has been developed based upon nonstandard rules. The equation of motion is formulated in the augmented dynamic space to apply the exponential mapping as a group preserving scheme. The formulation of the proposed method involves the hyperbolic sine and cosine functions. The method is therefore prone to divergence due to the behavior of the hyperbolic functions in structures with a high ratio of stiffness to mass. In the present study, to consider the properties of the structural equation into the formulation of the time step size and thereby avoid the divergence, a parameter, known as stability parameter, is thus derived from the exact solution of the equation of motion based on nonstandard rules. Embedding this parameter into the proposed method improves its stability. Afterward, for evaluating the performance of the proposed method, it is applied to several structures with different loading patterns while implemented in programing environment of the Matlab software. The results are compared to those of several commonly used numerical methods in structural applications. It is found that the proposed method has acceptable convergence and accuracy, and low time consumption compared to several commonly used methods. Furthermore, its stability is guaranteed by embedding the stability parameter into the proposed method.     

Flow Regimes in Packed Beds of Spheres from Pre-Darcy to Turbulent

Characteristics of flow regimes in porous media, along the processes of energy dissipation in each regime, are critical for applications of such media. The current work presents new experimental data for water flow in packed steel spheres of 1- and 3-mm diameters. The porosity of the porous media was about 35 % for both cases. The extensive dataset covered a broad range of flow Reynolds number such that several important flow regimes were encountered, including the elusive pre-Darcy regime, which is rarely or never seen in porous-media literature and turbulent regime. When compared to previous information, the results of this study are seen to add to the divergence of available data on pressure drop in packed beds of spheres. The divergence was also present in the coefficients of Ergun equation and in the Kozeny–Carman constant. The porous media of the current work were seen to exhibit different values of permeability and Forchheimer coefficient in each flow regime. The current data correlated well using the friction factor based on the permeability (measured in the Darcy regime) and the Reynolds number based on the same length scale. An attempt was made to apply recent theoretical results regarding the applicability of the quadratic and cubic Forchheimer corrections in the strong and weak inertia regime.     

Numerical prediction of turbulent flow in a conical diffuser usingk-ε model

Fully developed incompressible turbulent flow in a conical diffuser having a total divergence angle of 8° and an area ratio of 4∶1 has been simulated by ak-ε turbulence model with high Reynolds number and adverse pressure gradient. The research has been done for pipe entry Reynolds numbers of 1.16×10⁵ and 2.93×10⁵. The mean flow velocity and turbulence energy are predicted successfully and the advantage of Boundary Fit Coordinates approach is discussed. Furthermore, thek-ε turbulence model is applied to a flow in a conical diffuser having a total divergence angle of 30° with a perforated screen. A simplified mathematical model, where only the pressure drop is considered, has been used for describing the effect of the perforated screen. The optimum combination of the resistance coefficient and the location of the perforated screen is predicted for high diffuser efficiency or the uniform velocity distribution.