A second-order dynamic subgrid-scale stress model

A second-order dynamic model based on the general relation between the subgrid-scale stress and the velocity gradient tensors was proposed. A priori test of the second-order model was made using moderate resolution direct numerical simulation date at high Reynolds number ( Taylor microscale Reynolds number R-lambda = 102 similar to 216) for homogeneous, isotropic forced flow, decaying flow, and homogeneous rotating flow. Numerical testing shows that the second-order dynamic model significantly improves the correlation coefficient when compared to the first-order dynamic models.

Features of the supercritical transition in the wake of a circular cylinder

The features of the wake behind a uniform circular cylinder at Re = 200, which is just beyond the critical Reynolds number of 3-D transition, are investigated in detail by direct numerical simulations by solving 3-D incompressible Navier-Stokes equations using mixed spectral-spectral-element method. The high-order splitting algorithm based on the mixed stiffly stable scheme is employed in the time discretization. Due to the nonlinear evolution of the secondary instability of the wake, the spanwise modes with different wavelengths emerge. The spanwise characteristic length determines the transition features and global properties of the wake. The existence of the spanwise phase difference of the primary vortices shedding is confirmed by Fourier analysis of the time series of the spanwise vorticity and attributed. to the dominant spanwise mode. The spatial energy distributions of various modes and the velocity profiles in the near wake are obtained. The numerical results indicate that the near wake is in 3-D quasi-periodic laminar state with transitional behaviors at this supercritical Reynolds number.

Supersonic dusty-gas flows with Knudsen effect in interphase momentum exchange

On the basis of the two-continuum model of dilute gas-solid suspensions, the dynamic behavior of inertial particles in supersonic dusty-gas flows past a blunt body is studied for moderate Reynolds numbers, when the Knudsen effect in the interphase momentum exchange is significant. The limits of the inertial particle deposition regime in the space of governing parameters are found numerically under the assumption of the slip and free-molecule flow regimes around particles. As a model problem, the flow structure is obtained for a supersonic dusty-gas point-source flow colliding with a hypersonic flow of pure gas. The calculations performed using the full Lagrangian approach for the near-symmetry-axis region and the free-molecular flow regime around the particles reveal a multi-layer structure of the dispersed-phase density with a sharp accumulation of the particles in some thin regions between the bow and termination shock waves.

Numerical exploration of new features on three-dimensional transition of a cylinder wake

The three-dimensional transition of the wake flow behind a circular cylinder is studied in detail by direct numerical simulations using 3D incompressible N-S equations for Reynolds number ranging from 200 to 300. New features and vortex dynamics of the 3D transition of the wake are found and investigated. At Re = 200, the flow pattern is characterized by mode A instability. However, the spanwise characteristic length of the cylinder determines the transition features. Particularly for the specific spanwise characteristic length linear stable mode may dominate the wake in place of mode A and determine the spanwise phase difference of the primary vortices shedding. At Re = 250 and 300 it is found that the streamwise vortices evolve into a new type of mode - "dual vortex pair mode" downstream. The streamwise vortex structures switch among mode A, mode B and dual vortex pair mode from near wake to downstream wake. At Re = 250, an independent low frequency f(m) in addition to the vortex shedding frequency f(s) is identified. Frequency coupling between f(m) and f(s) occurs. These result in the irregularity of the temporal signals and become a key feature in the transition of the wake. Based on the formation analysis of the streamwise vorticity in the vicinity of cylinder, it is suggested that mode A is caused by the emergence of the spanwise velocity due to three dimensionality of the incoming flow past the cylinder. Energy distribution on various wave numbers and the frequency variation in the wake are also described.

Numerical estimate of the stability curve for the flow past a rotating cylinder

It is demonstrated that the primary instability of the wake of a two-dimensional circular cylinder rotating with constant angular velocity can be qualitatively well described by the Landau equation. The coefficients of the Landau equation are determined by means of numerical simulations for the Navier-Stokes equations. The critical Reynolds numbers, which depend on the angular velocity of the cylinder, are evaluated correctly by linear regression. (C) 2004 American Institute of Physics.

Numerical simulation of drop Marangoni migration under microgravity

Thermocapillary motion of a drop in a uniform temperature gradient is investigated numerically. The three-dimensional incompressible Navier-Stokes and energy equations are solved by the finite-element method. The front tracking technique is employed to describe the drop interface. To simplify the calculation, the drop shape is assumed to be a sphere. It has been verified that the assumption is reasonable under the microgravity environment. Some calculations have been performed to deal with the thermocapillary motion for the drops of different sizes. It has been verified that the calculated results are in good agreement with available experimental and numerical results. (C) 2003 Elsevier Ltd. All rights reserved.

Two-way coupling model for shock-induced laminar boundary-layer flows of a dusty gas

The present paper describes a numerical two-way coupling model for shock-induced laminar boundary-layer flows of a dust-laden gas and studies the transverse migration of fine particles under the action of Saffman lift force. The governing equations are formulated in the dilute two-phase continuum framework with consideration of the finiteness of the particle Reynolds and Knudsen numbers. The full Lagrangian method is explored for calculating the dispersed-phase flow fields (including the number density of particles) in the regions of intersecting particle trajectories. The computation results show a significant reaction of the particles on the two-phase boundary-layer structure when the mass loading ratio of particles takes finite values.

负压磁头气动力静态特性

该文用有限差分法离散Reynolds润滑方程,采用交替方向隐式格式(ADI)求解离散得到的代数方程组,计算出了负压滑块的静态气动力特性.

磁头润滑面形状对气动力特性的影响

该文采用微分离散差分方法求解滑流条件下的修正Reynolds方程.对于一阶导数间断的物理量-气膜厚度采用高精度保凸性磨光法进行处理.给出了气膜压强分布立体图以及几种因加工误差造成磁头滑块润滑面变形的误差型面的飞升曲线,中心线上压强分布,压力中心和负载随加工误差大小的变化情况,指出了磁头滑块润滑面加工时宁凸勿凹的倾向性意见.

内凹式负压磁头气体动力特性

本文给出了一种用于磁盘机读/写功能的非接触起停式负压磁头滑块模型。通过求解修正Reynolds方程,计算出其稳态气动力特性。给出了磁头滑块与磁盘之间润滑气膜压强分布的立体图,滑块工作曲线及压力载荷随磁盘线速度、气膜厚度、减压槽深度的变化规律曲线。

钝体分离旋涡流动的区域分解、杂交数值模拟——Ⅰ.理论方法及其应用

为克服涡旋法不能精确预计物体附近小尺度流动结构的理论缺陷,减少高Reynolds数流动N-S方程差分解的困难,本文提出一种区域分解、杂交耦合N-S方程有限差分解及涡旋法的新的数值模型和理论方法.将流场分解为内外两区,在靠近物体表面、范围为O(R)的内区进行N-S方程有限差分解,外区作Lagrange-Euler涡旋法解,建立了分区流动的联结、耦合条件,给出了杂交耦合求解的数值计算方法.用本方法作了Re=10~2,10~3的圆柱绕流计算,考察了区域交界面位置变化时解的稳定性.与全场N-S方程解及实验结果的比较表明本文方法能精确预计流动分离及近场流动的详细结构,并可有效地计算流动的总体特性,且比全场N-S方程解显著节省机时和计算量.

可瘪管定常流动的若干新现象

在较宽的参数范围内进行了可瘪管定常流动实验。当流速-波速比较高时,发现了若干新现象:(ⅰ)当管道开始被压瘪,或已压瘪的管道即将完全开启时,流态发生突变,从定常流变为脉动流(或相反);同时,流量、压力发生阶跃;(ⅱ)压差-流量曲线出现双峰;(ⅲ)在正阻尼区亦可能发生持续的振荡。讨论了实变发生的原因和激振机理。

非线性流动稳定性理论

<正> 流体的湍流运动普遍存在于大气、海洋、飞行器周围、推进装置和流体机械中,探讨湍流形成的条件和过程是流动稳定性理论的研究对象。 早在1843年,Stokes就预见到流体状态转捩的原因是失稳。整整100年前(1883年),Reynolds让液体流入不同口径的圆管,并在对称轴上注入一股纤细的颜色水,以便明显、敏感地观察到湍流的发生。他发现: (1)流速较低时,可以看到层次分明的层流,流速增加到一定程度后,就转变成高度无序的湍流状态;

短熔体柱的定常轴对称热毛细运动

本文处理了短熔体柱的定常轴对称热毛细运动。假定Marangoni数M1,Reynolds数Re■1以及毛细作用数■1。本文求得了级数形式的零级近似分析解,并对两种典型的外加温度分布的情况作出了数值计算。结果表明了端部效应对热毛细运动的重要影响。

Orr-Sommerfeld方程的Green函数和积分方程解法

本文提出了一种求解线性稳定性理论Orr-Sommerfeld方程的方法.我们首先定义了该方程的Green函数,并将它表达成矩阵形式;然后证明了Green函数的互易性;最后导出了等价于原方程的线性积分方程.该方法适用于两固壁间任意Reynolds数下各种主流速度分布的情况.