均匀、振荡及组合流绕平板流动的尾迹研究

本文利用离散涡模型及改进的新生涡产生机制对三种不同来流绕平板的近尾迹进行数值研究。计算讨论了定常流中平板绕流流动的总体特性和近尾迹流场;对于简谐振荡来流,相应于K_c=2.0、4.0 和10.0 分别得到两种不同的尾迹形态。给出了小 K_c 数平板尾迹涡配对、运动的新模式而相应的阻力、惯性力系数计算比以前涡模拟结果更接近于 U 型管实验结果。对于流向组合来流本文模拟了涡锁定及其动力特性并于实验相符,给出了流向扰动对平板绕流流动的影响。     

高雷诺数下多柱绕流特性研究

采用改进的延迟分离涡方法数值模拟了高雷诺数下的柱体绕流,包括单圆柱绕流、单方柱绕流、串列双圆柱绕流和串列双方柱绕流,研究了不同雷诺数下圆柱绕流与方柱绕流的水动力特性.计算结果与实验数据及其他文献的数值计算结果吻合良好,研究表明,单方柱绕流在2.0×10~3Re1.0×10~7范围内未出现类似于单圆柱绕流的阻力危机现象,其平均阻力系数C_d、升力系数均方根C'_1及斯特劳哈尔数S t维持在一定范围内波动.串列双圆柱绕流与串列双方柱绕流中,均选取L/D=2.0,2.5,3.0,3.5和4.0这五中间距比进行计算.串列双圆柱绕流中,当Re=2.2×10~4时,在3.0L/D3.5内存在一临界间距比(L_c/D)使得L_c/D前后上下游圆柱的升阻力系数发生跳跃性变化,且当L/DL_c/D时,下游圆柱的阻力系数为负数.而当Re=3.0×10~6时,则不存在临界间距比,且下游圆柱的阻力系数始终为正数.串列双方柱绕流在Re=1.6×10~4和Re=1.0×10~6两种工况下的临界间距比分别处于3.0L/D3.5和3.5L/D4.0区间内,且当L/DL_c/D时,两个雷诺数下的下游方柱阻力系数均为负数.     

棱柱绕流流动的数值模拟

利用数值方法对长宽比为1／3, 1和3的棱柱绕流在雷诺数为100的非稳态流动特性进行了分析和研究。采用有限体积法对棱柱绕流的二维流动N-S方程进行离散求解,分析和研究了非稳态的棱柱绕流流场,升力系数,阻力系数和涡动特性,数值模拟的结果与相关文献的数据比较吻合。通过上述研究能够为了解棱柱绕流的非稳态流动特性提供有力的帮助。而对棱柱三维流动的模拟分析和对雷诺数的变化对棱柱流动特性的影响进行研究,将为掌握棱柱绕流的工程特性打下基础。     

绕振荡水翼流动及其转捩特性的数值计算研究简

通过对比标准k-ω SST 湍流模型和基于标准k-ω SST 湍流模型修正的γ-Re_θ 转捩湍流模型对绕振荡NACA66 水翼流动的数值计算结果与实验结果,对水翼振荡过程的水动力特性和流场结构变化进行了分析研究. 结果表明：与标准k-ω SST 湍流模型的数值计算结果相比,基于标准k-ω SST 湍流模型修正的γ-Re_θ 转捩湍流模型能有效预测绕振荡翼型流场结构和水动力特性,捕捉流场边界层发生的流动分离和转捩现象;绕振荡水翼的流动过程可分为5 个特征阶段,当来流攻角较小时,在水翼前缘发生层流向湍流的转捩现象,水翼动力特征曲线出现变化拐点;随着来流攻角的增大,顺时针尾缘涡逐渐形成并向水翼前缘发展;当攻角较大时,前缘涡分离导致动力失速,水翼的动力特征曲线出现大幅波动;水翼处于顺时针向下旋转阶段,绕水翼的流动状态逐渐由湍流过渡为层流.     

绕振荡水翼流动及其转捩特性的数值计算研究

通过对比标准k-ω SST 湍流模型和基于标准k-ω SST 湍流模型修正的γ-Re_θ 转捩湍流模型对绕振荡NACA66 水翼流动的数值计算结果与实验结果,对水翼振荡过程的水动力特性和流场结构变化进行了分析研究. 结果表明:与标准k-ω SST 湍流模型的数值计算结果相比,基于标准k-ω SST 湍流模型修正的γ-Re_θ 转捩湍流模型能有效预测绕振荡翼型流场结构和水动力特性,捕捉流场边界层发生的流动分离和转捩现象;绕振荡水翼的流动过程可分为5 个特征阶段,当来流攻角较小时,在水翼前缘发生层流向湍流的转捩现象,水翼动力特征曲线出现变化拐点;随着来流攻角的增大,顺时针尾缘涡逐渐形成并向水翼前缘发展;当攻角较大时,前缘涡分离导致动力失速,水翼的动力特征曲线出现大幅波动;水翼处于顺时针向下旋转阶段,绕水翼的流动状态逐渐由湍流过渡为层流.      

环境流体力学 (上)

1.本文讨论了环境流体力学的一般领域和这个研究课题的应用。2.在理论、计算、物理模拟和现场研究方面,评述了大气和海洋中流体流动的一些问题。因为它们对于环境问题实际上很重要。例如在表面层、混合层、空气流过粗糙度不断变化的表面和流体流过沙波或海湾岩床等情况下出现的分层效应。3.也讨论了人类活动和环境的相互作用产生的一些流体力学问题。例如作用在建筑物上的定常流动、振荡流动和旋涡流动的流体力(简略涉及波能);障碍物的绕流和下游流动(包括高层建筑物附近风对人的影响);空气中(包括稠密气体的影响)和水道中(也括“热”污染)污染物的扩散。      

基于格子Boltzmann方法的低雷诺数四圆柱绕流流动模式与受力特性模拟

利用格子Boltzmann方法模拟了雷诺数为100时,均匀来流条件下的二维菱形排布的四柱绕流现象,得到了不同柱间距比下的绕流流动模式及阻力变化规律。结果表明:圆柱互扰效应与柱间距比有关,当L/D≤1.2时为单钝体模式,圆柱互扰效应以临近效应为主;当1.2相似文献   

表面活性剂溶液中弹性胶束断裂作用下圆柱绕流特征研究

表面活性剂分子溶于水自组装形成的蠕虫状胶束溶液由于其高表面活性和高黏弹性,具有独特的力学响应特性,呈现出与牛顿流体迥异的流动特征。采用考虑溶液中胶束可逆断裂和重构效应的双组分VCM本构方程,研究了溶液流变学特性,发现溶液拉伸应力随应变率呈现出拉伸硬化-稀化-硬化的强非线性特征。研究了蠕虫状胶束溶液绕圆柱流动的特性,在低雷诺数(Re～0.01)下,当维森贝格数Wi>1时,流场出现较强的不稳定性,绕流尾迹区流场呈现随时间波动的特征,通过功率谱分析发现流场波动具有拟周期特征,当胶束断裂引起流体拉伸应力对弹性响应的非线性增强时,流场波动的不规则性增强。圆柱后缘附近出现一对周期性变化的对称涡,表明流体弹性会促进流动分离的发生,流体拉伸应力对弹性变化的强非线性响应是造成低雷诺数绕流流动不稳定性的根本原因,揭示了蠕虫状胶束溶液在低雷诺数高黏弹性时绕流流动不稳定性的产生机理。     

具有水平流动的Rayleigh-Benard对流斑图

本文利用数值模拟,探讨了普朗特数 时有水平流动的Rayleigh-Benard对流结构．当水平流动强度 时,发现定常对流的多重稳定性．当 时,Rayleigh-Benard(RB)对流中存在三种对流斑图．它们的出现依赖于水平流动强度 和相对瑞利数 ．与 相比,在 时的行波具有不同的动力学特性．     

二层密度分层流体重力流的数值模拟

从RANS方程和RNGk-ε湍流模型出发,采用流体体积法(VOF)来模拟密度分层流动,对盐水和淡水因密度差异导致的分层重力流动现象进行了数值模拟.文中报道了平底水槽重力流、狭孔交换流的数值模拟结果,分层重力流锋面运动速度的计算值与现有的半理论半经验公式一致.为了揭示地形变化对分层重力流的影响,对设有缓变潜堤的水槽内分层重力流动的形成过程进行了数值模拟,给出了重力头推进速度和局部流场的计算结果,并讨论了分层流界面、流量和锋面附近的流速分布特征.     

Resolving unsteady micro-scale atmospheric flows by nesting a CFD simulation into wide range numerical weather prediction models

The open boundary conditions for the CFD simulation of the micro- and meso-scale flow and temperature distribution around the Grimming mountain, Austria, are determined by an optimisation approach from interior observations. The numerical weather prediction model ALADIN–Austria provides wind speed and wind direction at those spatially arbitrarily defined observations. Furthermore, the plausibility of the resulting flow over the Grimming mountain is checked (a) by data of a ground station at the top of the Grimming mountain showing suitable correlation with the measurements. Besides, (b) the vertical turbulence profiles at the ground station are qualitatively compared to extensive studies covering the Askervein hill. In addition, (c) the wavelength of internal gravity waves obtained from the numerical CFD model is examined by analytic approximations at the Grimming mountain and by linear mountain wave theory at a mathematically idealised hill. Finally, the grid independence of the presented CFD model is shown.     

Analytical investigation on 3D non-Boussinesq mountain wave drag for wind profiles with vertical variations

A new analytical model was developed to predict the gravity wave drag (GWD) induced by an isolated 3-dimensional mountain, over which a stratified, non-rotating non-Boussinesq sheared flow is impinged. The model is confined to small amplitude motion and assumes the ambient velocity varying slowly with height. The modified Taylor-Goldstein equation with variable coefficients is solved with a Wentzel-Kramers-Brillouin (WKB) approximation, formally valid at high Richardson numbers. With this WKB solution, generic formulae of second order accuracy, for the GWD and surface pressure perturbation (both for hydrostatic and non-hydrostatic flow) are presented, enabling a rigorous treatment on the effects by vertical variations in wind profiles. In an ideal test to the circular bell-shaped mountain, it was found that when the wind is linearly sheared, that the GWD decreases as the Richardson number decreases. However, the GWD for a forward sheared wind (wind increases with height) decreases always faster than that for the backward sheared wind (wind deceases with height). This difference is evident whenever the model is hydrostatic or not.     

海堤越浪的数值模拟

基于RANS方程和两方程湍流模型,采用有限体积法,将人射波波场作为人工的分布源项加人动量方程,提出了适用于VOF方法的源造波一消波技术。通过对行波及驻波的计算,分别考察了数值波浪水槽前端及末端消波段的有效性。在本文建立的数值波浪水槽内对规则波在海堤上爬高和越浪过程进行了数值模拟,并将计算结果与现有实验结果进行了比较。验证计算结果表明,数值模拟结果较好地复演了海堤越浪过程。为了研究模型尺度对越浪量的影响,文中设计了一组满足重力相似但具有不同几何比尺的数值实验模型。系列数值实验结果表明,若按重力相似换算越浪量,计算结果与实验预报值间的偏差随模型比尺的增大和堤前波浪破碎强度的增强而增大,建议在进行越浪物理模型实验时需进一步考虑模型比尺对原型预报值的影响。     

Stratified flow under an ice keel—a numerical approach

A numerical investigation of the flow of two immiscible stratified fluids under an isolated keel has been undertaken. The investigation utilized the two-dimensional Euler equations for incompressible flow, and the solution of these equations has been obtained by using the well-known finite volume marker and cell approach. Experimental drag-force measurements are also presented for a family of two-dimensional topographic models of fixed height with increasing surface slopes in a two-layer density system. The range of flow speeds explored covers the Froude number range from subcritical to fully supercritical. The drag force measurements are augmented by detailed observations of the interface distortion. The results clearly show large drag increases arising from the internal wave systems generated in the stratified flow. Very good agreement has been found between the experimental and numerical results for both the interface shape between the two fluids and the drag force on a variety of keels.     

Numerical study on heavy rigid particle motion in a plane wake flow by spectral element method

Experimental particle dispersion patterns in a plane wake flow at a high Reynolds number have been predicted numerically by discrete vortex method (Phys. Fluids A 1992; 4 :2244–2251; Int. J. Multiphase Flow 2000; 26 :1583–1607). To address the particle motion at a moderate Reynolds number, spectral element method is employed to provide an instantaneous wake flow field for particle dynamics equations, which are solved to make a detail classification of the patterns in relation to the Stokes and Froude numbers. It is found that particle motion features only depend on the Stokes number at a high Froude number and depend on both numbers at a low Froude number. A ratio of the Stokes number to squared Froude number is introduced and threshold values of this parameter are evaluated that delineate the different regions of particle behavior. The parameter describes approximately the gravitational settling velocity divided by the characteristic velocity of wake flow. In order to present effects of particle density but preserve rigid sphere, hollow sphere particle dynamics in the plane wake flow is investigated. The evolution of hollow particle motion patterns for the increase of equivalent particle density corresponds to that of solid particle motion patterns for the decrease of particle size. Although the thresholds change a little, the parameter can still make a good qualitative classification of particle motion patterns as the inner diameter changes. Copyright © 2008 John Wiley & Sons, Ltd.     

可侵蚀地表上方含尘大气运动的数值模型

在稀相气固两相流的双流体模型框架下，导出含尘大气运动的基本方程，其中的源项考虑相间双向耦合作用；通过对基本方程进行无量纲化，求得控制气相和弥散相动力学行为的相似准则，包括弗洛德数、罗斯贝数、颗粒雷诺数、颗粒惯性参数、颗粒剪切横移参数、颗粒旋转横移参数和沙尘质量载荷率等. 作为一个模型问题，研究可侵蚀地表上方充分发展湍流大气边界层流动对土壤颗粒的气动卷扬作用，其中计及气动阻力、萨弗曼升力和重力等因素的影响. 为了克服由于沙尘跃移运动轨迹交叉导致流动参数多值性所造成的困难，引进拉格朗日坐标下弥散相连续方程. 对两种不同风速和4种不同粒径条件下的沙尘运动进行数值模拟，给出相应情况中沙尘运动特性和密度分布剖面并讨论风速和粒径等参数的影响，还细致分析了含尘大气边界层中沙尘与气流之间的能量传递过程. 所得的方程、准则和方法可以为土壤风蚀和沙尘暴等自然灾害的预报提供理论基础.     

A numerical investigation of developing flow in an eccentric curved annulus in the presence of gravity

In this article developing incompressible viscous flow in an eccentric curved annulus in the presence of gravity is numerically studied using a second order finite difference method based on the projection algorithm to solve the governing equations including the continuity and full Navier–Stokes equations. The equations written in a bipolar–toroidal coordinate system are discretized in a three dimensional staggered grid. The effects of governing non-dimensional parameters including the eccentricity, non-dimensional curvature ratio, Dean number, Froude number, aspect ratio, and the Reynolds number on the flow field in the entrance and fully developed region are investigated. The numerical results indicate that at the small Froude numbers, the flow field distorts from the symmetrical condition due to the larger body force effect and the axial velocity formation mostly takes place at the lower half of the annulus. In addition, at the constant Froude number, by decreasing the curvature radius, the peak axial velocity and its sharp gradient appear on the outer curvature region due to the larger centrifugal forces and by increasing the eccentricity the flow rate intensifies at the wider region and weakens at the narrower region due to the larger flow resistance. Furthermore, the friction factor increases by decreasing the Froude number and increasing the Dean number.     

Internal wave generation by a fountain in a stratified fluid

The purpose of the study is the direct numerical and theoretical modeling of fountain dynamics in a fluid with density stratification in the form of a pycnocline. The fountain is formed as a vertical jet penetrates through the pycnocline. In numerical simulation the jet flow is initiated by means of preassigning a boundary condition in the form of an upward-directed laminar flow of a neutral-buoyancy fluid with an axisymmetric Gaussian velocity profile. The calculations show that at a Froude number Fr greater than a certain critical value the flow becomes unstable and the fountain executes self-oscillations accompanied by internal wave generation in the pycnocline. Depending on Fr, two self-oscillation modes can be distinguished. At fairly low Fr the fountain executes circular motion in the horizontal plane, in the vicinity of the center of jet, its shape remaining almost invariant. In this case, internal waves in the form of unwinding spirals are radiated. At fairly high Fr another mode predominates, when the fountain top chaotically “strays” in the vicinity of the center of the jet and, periodically breaking down, generates wave packets propagating toward the periphery of the computation domain. In both cases, the main peak in the frequency spectrum of the internal waves coincides with the fountain top oscillation frequency which monotonically decreases with increase in Fr. In numerical simulation the Fr-dependence of the fountain top oscillation amplitude is in good agreement with that predicted by the theoretical model of the concurrence of the interacting modes in the soft self-excitation regime.     

Shallow flow simulation on dynamically adaptive cut cell quadtree grids

A computationally efficient, high‐resolution numerical model of shallow flow hydrodynamics is described, based on dynamically adaptive quadtree grids. The numerical model solves the two‐dimensional non‐linear shallow water equations by means of an explicit second‐order MUSCL‐Hancock Godunov‐type finite volume scheme. Interface fluxes are evaluated using an HLLC approximate Riemann solver. Cartesian cut cells are used to improve the fit to curved boundaries. A ghost‐cell immersed boundary method is used to update flow information in the smallest cut cells and overcome the time step restriction that would otherwise apply. The numerical model is validated through simulations of reflection of a surge wave at a wall, a low Froude number potential flow past a circular cylinder, and the shock‐like interaction between a bore and a circular cylinder. The computational efficiency is shown to be greatly improved compared with solutions on a uniform structured grid implemented with cut cells. Copyright © 2006 John Wiley & Sons, Ltd.     

Flutter and resonances of a flag near a free surface

We investigate the effects of a nearby free surface on the stability of a flexible plate in axial flow. Confinement by rigid boundaries is known to affect flag flutter thresholds and fluttering dynamics significantly, and this work considers the effects of a more general confinement involving a deformable free surface. To this end, a local linear stability is proposed for a flag in axial uniform flow and parallel to a free surface, using one-dimensional beam and potential flow models to revisit this classical fluid–structure interaction problem. The physical behaviour of the confining free surface is characterized by the Froude number, corresponding to the ratio of the incoming flow velocity to that of the gravity waves. After presenting the simplified limit of infinite span (i.e. two-dimensional problem), the results are generalized to include finite-span and lateral confinement effects. In both cases, three unstable regimes are identified for varying Froude number. Rigidly-confined flutter is observed for low Froude number, i.e. when the free surface behaves as a rigid wall, and is equivalent to the classical problem of the confined flag. When the flow and wave velocities are comparable, a new instability is observed before the onset of flutter (i.e. at lower reduced flow speed) and results from the resonance of a structural bending wave and one of the fundamental modes of surface gravity waves. Finally, for large Froude number (low effect of gravity), flutter is observed with significant but passive deformation of the free surface in response of the flag’s displacement.