方腔内自然对流与混合对流的有限差分格子Boltzmann研究

基于有限差分法,建立了贴体坐标系下求解流体流动和传热的双分布格子Boltzmann模型.在密度分布函数和温度分布函数对应的离散速度方程中,时间项采用四阶Runge-Kutta法离散,空间离散采用二阶迎风和二阶中心差分的混合形式.采用此模型分别对瑞利数为10~3、10~4、10~5、10~6的方腔自然对流以及理查森数为0.1、1、10的方腔混合对流进行了数值模拟,获得了流体速度与温度分布的典型特征,得到的努塞尔数也与基准解高度吻合.计算结果表明了本文采用的数值方法和计算程序的有效性.     

周期结构的部分组集有限元法

本文分析了周期结构经有限元离散所形成系数矩阵的元素分布，提出部分组集的有限元法．算例表明部分组集大大降低了整体结构系数矩阵形成所需内外存空间与计算时间     

波浪场中物质输移离散系数的理论分析

海岸水域波浪引起的物质输移扩散存在着不同于水流引起的物质输移扩散的特征. 通过解析的方法对垂向扩散方程进行求解, 推导出波浪和潮流共同作用下的海岸水域物质输移离散系数的理论公式. 分析中将总的离散系数分为潮流、质量输移流、纯波浪波动作用以及潮流与质量输移流、潮流与波浪波动、质量输移流与波浪波动的相互作用6部分, 通过分析各组成部分特征, 得出海岸水域物质输移中波浪所产生的离散效应主要是质量输移流贡献, 波浪水质点运动及波流相互作用对时间平均离散系数的贡献较小的结论.仅波浪的作用(不包括潮流)的有关结果与他人数值模拟和实验结果的吻合良好, 证明了该理论推导的正确性. 给出了时间平均离散系数及离散系数波动幅值随波浪周期和波高的变化特征,同时将结果应用于垂向环流, 给出了垂向环流对应于质量输移流所产生的离散系数的特征.      

90°弯道明渠横向环流的大涡模拟研究

采用气液两相流大涡模型模拟90°弯道明渠水流的水力特性。使用有限体积法对控制方程进行离散;使用压力隐式算子分裂PISO(Pressure-Implicit with Splitting of Operators)算法求解速度与压力的耦合;采用了VOF(Volume of Fluid)法模拟自由水面。通过计算得到弯道处流速、压强系数、横向环流强度等水力参数的分布规律。分别比较弯道纵向流速与横向流速的模拟值与实验值,二者吻合较好。通过分析弯道环流的流态特征得出:因受重力作用与水流在流经弯道时发生的离心现象的影响,凸岸附近较凹岸存在较大的二次流强度的结果;凸岸与凹岸压强系数均沿程增大,当达到弯顶时,水面横比降达到最大值;弯道横向环流强度随弯道角度增大而逐渐增大,其最大值出现在明渠底部。     

稀薄流到连续流的气体运动论模型方程算法研究

通过引入碰撞松弛参数和当地平衡态分布函数对BGK模型方程进行修正，确定含流态控制参数可描述不同流域气体流动特性的气体分子速度分布函数的简化控制方程。发展和应用离散速度坐标法于气体分子速度空间，利用一套在物理空间和时间上连续而速度空间离散的分布函数来代替原分布函数对速度空间的连续依赖性。基于非定常时间分裂数值计算方法和无波动、无自由参数的NND耗散差分格式，建立直接求解气体分子速度分布函数的气体运动论有限差分数值方法。推广应用改进的Gauss-Hermite无穷积分法和华罗庚－王元提出的以单和逼近重积分的黄金分割数论积分方法等，对离散速度空间进行宏观取矩获取物理空间各点的气体流动参数，由此发展一套从稀薄流到连续流各流域统一的气体运动论数值算法。通过对不同Knudsen数下一维激波管问题、二维圆柱绕流和三维球体绕流的初步数值实验表明文中发展的数值算法是可行的。     

一维溃坝洪水波的高精度数值模拟

将ENO(Essentially Non-Oscillatory)格式和Runge-Kutta时间离散的思想应用于一维Saint-Venant方程组的求解,数值模拟溃坝洪水,得出了水位和流速的沿程分布。经与理论解比较可见,数值解在间断波附近没有出现数值振荡,水位和流速大小均符合较好,表明ENO格式是一类新的高精度无振荡差分格式,采用ENO格式所建立的高分辨率模型能够很好地模拟溃坝波的演进过程。     

弯道水流二维浅水模型的修正研究

借鉴有关弯道水流流速分布的研究成果,计入深度平均流速与真实流流速分布差值引起的扩散效应,在正交曲线坐标系下建立了平面二维浅水模型.采用以标准κ-ε模型为基础的曲率效应修正紊流模型模拟紊动应力项,在一定程度上考虑了流线弯曲水流紊动应力的各向异性.应用控制体积法和交错网格法离散方程,并用SIMPLEC算法求解离散方程;同时采用修正后的模型对90°弯道水流进行了数值模拟,并与原模型的计算结果及实测资料进行了比较,结果表明该模型能够有效地模拟流线弯曲水流的水力特性.     

热线近壁测速的实验研究

本文分别以导热系数相差很大的材料做成实验平板，在常温与加热的情况下，用恒温热线风速仪和不同热线探头，测量了靠近固壁１ｍｍ以内的平板边界层时均流速分布。实验结果表明，有多种因素影响粘性次层流速分布的测量结果。本文展示这些因素，并进行了分析讨论。     

黏性阻尼系统时域响应灵敏度及其一致性研究

时域响应灵敏度分析是时域梯度优化算法的基础. 灵敏度分析通常只涉及对设计变量的微分运算, 但时域响应灵敏度问题还涉及时间域的离散化. 因此, 微分和离散的先后顺序可能对时域响应灵敏度结果产生影响. 针对黏性阻尼系统时域响应灵敏度求解问题, 基于改进精细积分方法, 分别推导了先微分后离散和先离散后微分两种伴随变量方法. 其中, 先微分后离散法首先对由伴随变量构造的增广函数微分, 再利用改进精细积分方法在各离散时间点求解时域响应灵敏度; 而先离散后微分方法则首先在各离散时间点引入残值方程构造增广函数, 再对各增广函数进行微分以求解时域响应灵敏度. 通过数值算例验证了所提出方法的有效性和准确性, 并与传统基于Newmark的方法进行比较. 结果表明, 积分方案、数值离散误差以及离散和微分的先后顺序共同影响灵敏度的一致性误差. 综合考虑精度、效率和一致性问题, 基于改进精细积分的先微分后离散伴随变量法表现更优, 最适合应用于黏性阻尼系统时域梯度优化算法.      

移动最小二乘浸没边界法中的权函数影响分析

考察了浸没边界法中运用移动最小二乘法构造的插值形函数的影响。通过在移动最小二乘法中采用不同的权函数,分析了相应的插值形函数的性质,并与传统的离散delta函数做了比较。以静止流体中的水平振荡圆柱为例,阐明了形函数对圆柱阻力系数的幅值和光滑性的影响,得到了较优的形函数分布,并将结果与文献对比验证了本文方法的可靠性。     

Dispersion of fine settling particles from an elevated source in an oscillatory turbulent flow

The present paper examines the stream-wise dispersion of suspended fine particles with settling velocities in an oscillatory turbulent shear flow with or without a non-zero mean over a rough-bed surface when the particles are being released from an elevated continuous source. A finite-difference implicit method is employed to solve the unsteady turbulent convective-diffusion equation. A combined scheme of central and four-point upwind differences is used to solve the steady state equation and the Alternating Direction Implicit (ADI) method is adopted for unsteady equation. It is shown how the mixing of settling particles is influenced by the tidal oscillatory current and the corresponding eddy diffusivity when the initial distribution of concentration regarded as a line-source. The vertical concentration profiles of suspended fine particles with settling velocities are presented for different downstream stations for various values of settling velocity and the frequency of the oscillation in tidal flow. For two-dimensional unsteady dispersion equation, the behaviour of iso-concentration lines for different values of settling velocity, frequency of the oscillation, dispersion time and releasing height is studied in terms of the relative importance of convection and eddy diffusion.     

Gravity Affected Lateral Dispersion and Diffusion in a Stationary Horizontal Porous Medium Shear Flow

Transport of dissolved species by a carrier fluid in a porous medium comprises advection and diffusion/dispersion processes. Hydrodynamic dispersion is commonly characterized by an empirical relationship, in which the dispersion mechanism is described by contributions of molecular diffusion and mechanical dispersion expressed as a function of the molecular Peclét number. Mathematically these two phenomena are modeled by a constant diffusion coefficient and by velocity dependent dispersion coefficients, respectively. Here, the commonly utilized Bear--Scheidegger dispersion model of linear proportionality between mechanical dispersion and velocity, and the more complicated Bear--Bachmat model derived on a streamtube array model porous medium and better describing observed dispersion coefficients in the moderate molecular Peclét number range, will be considered. Analyzing the mixing flow of two parallelly flowing confluent fluids with different concentrations of a dissolved species within the frames of boundary layer theory one has to deal with transverse mixing only. With the Boussinesq approximation being adopted approximate analytical solutions of the corresponding boundary layer system of equations show that there is no effect of density coupling on concentration distributions across the mixing layer in the pure molecular diffusion regime case. With the Peclét number of the oncoming flow growing beyond unity, density coupling has an increasing influence on the mixing zone. When the Peclét number grows further this influence is successively reduced until its disappearance in the pure mechanical dispersion regime.     

Approximate solutions for the problem of liquid film flow over an unsteady stretching sheet with thermal radiation and magnetic field

The proposed method is based on replacement of the unknown function by a truncated series of the shifted Legendre polynomial expansion. An approximate formula of the integer derivative is introduced. Special attention is given to study the convergence analysis and derive an upper bound of the error for the presented approximate formula. The introduced method converts the proposed equation by means of collocation points to a system of algebraic equations with shifted Legendre coefficients. Thus, after solving this system of equations, the shifted Legendre coefficients are obtained. This efficient numerical method is used to solve the system of ordinary differential equations which describe the thin film flow and heat transfer with the effects of the thermal radiation, magnetic field, and slip velocity.     

Impact of local diffusion on macroscopic dispersion in three-dimensional porous media

While macroscopic longitudinal and transverse dispersion in three-dimensional porous media has been simulated previously mostly under purely advective conditions, the impact of diffusion on macroscopic dispersion in 3D remains an open question. Furthermore, both in 2D and 3D, recurring difficulties have been encountered due to computer limitation or analytical approximation. In this work, we use the Lagrangian velocity covariance function and the temporal derivative of second-order moments to study the influence of diffusion on dispersion in highly heterogeneous 2D and 3D porous media. The first approach characterizes the correlation between the values of Eulerian velocity components sampled by particles undergoing diffusion at two times. The second approach allows the estimation of dispersion coefficients and the analysis of their behaviours as functions of diffusion. These two approaches allowed us to reach new results. The influence of diffusion on dispersion seems to be globally similar between highly heterogeneous 2D and 3D porous media. Diffusion induces a decrease in the dispersion in the direction parallel to the flow direction and an increase in the dispersion in the direction perpendicular to the flow direction. However, the amplification of these two effects with the permeability variance is clearly different between 2D and 3D. For the direction parallel to the flow direction, the amplification is more important in 3D than in 2D. It is reversed in the direction perpendicular to the flow direction.     

Simulation of turbulent dispersion using Weierstrass modes

A simple fractal model is proposed for the dispersion of passive scalars in an incompressible homogeneous turbulent flow field. The dispersion process is based on a three-dimensional velocity field which is assumed to be a linear superposition of Taylor-Green vortices with wave numbers and amplitudes as those in a Weierstrass function. The chosen velocity field satisfies the continuity equation as well as Kolmogorov's inertial range power law, and has a fractal dimension D between upper and lower length-scale bounds. A cloud of tracer particles is relased into the flow field and dispersed by the fluid motion, the trajectories are numerically integrated from the velocity field. The puffs which result from this process are used to examine some aspects of turbulent dispersion, through comparisons with integrated concentration wind-tunnel measurements. The agreement between numerical and experimental results indicates the significance of the proposed simulation model.     

Filtration model of longitudinal flow in a finned cylindrical channel

The problem of laminar flow of a viscous incompressible fluid in a finned circular tube is considered. A solution is obtained in the form of series in eigenfunctions of the Laplace operator; the coefficients in the series are found numerically. For the same problem, a simpler filtration approximation is proposed in which the system of fins is modeled by a radially inhomogeneous porous layer, and fluid flow in it is described by the Brinkman equation. A formula for the effective permeability of the porous medium is obtained by varying the number and height of fins. The formula provides an accurate evaluation of the mean flow velocity and viscous drag coefficient in finned channels.     

Equivalent Diffusion in a Thin Film Flow with a Non-One-Dimensional Velocity Field and Anisotropic Diffusion Tensor

For describing the mass transfer processes in channels, Taylor's dispersion theory is widely used. This theory makes it possible, with asymptotic rigor, to replace the complete diffusion (heat conduction) equation with a convective term that depends on the coordinate transverse to the flow by an effective diffusion (dispersion) equation with constant coefficients, averaged over the channel cross-section. In numerous subsequent studies, Taylor's theory was generalized to include more complex situations, and novel algorithms for constructing the dispersion equations were proposed. For thin film flows a theory similar to Taylor's leads to a matrix of dispersion coefficients.In this study, Taylor's theory is extended to film flows with a non-one-dimensional velocity field and anisotropic diffusion tensor. These characteristics also depend to a considerable extent on the spatial coordinates and time. The dispersion equations obtained can be simplified in regions in which the effective diffusion coefficient tensor changes sharply.     

Probabilistic approach for transport of contaminants through porous media

The channels formed between individual particles in porous media have variable dimensions and orientations. The porosity, permeability and its anisotropy exhibit random spatial distributions. The probabilistic approach can effectively describe the transport of contaminants through porous media and is analysed in this paper. Numerical results are obtained by considering (I) random dispersion coefficients without and with spatial structure, (II) random time distribution of concentration at the inlet boundary, (III) random velocity distribution in the flow field without and (IV) with variable dispersion coefficient, (V) non-linearity of the governing equation and (VI) anisotropy of the dispersion coefficient. Two methods are used for probabilistic predictions: (1) Gaussian field approach in conjunction with Monte Carlo method and (2) random walk method. The input random parameters are assumed to have normal and log-normal distributions according to available experimental data. The probability distribution functions of the contaminant concentration at different locations within the flow domain are calculated and compared with the input distributions as a function of the mean and fluctuation Peclet numbers. The one-dimensional case is analysed in detail and the illustrative numerical predictions are compared with analytical and experimental results. The extension to a two-dimensional domain is discussed in the last part of this paper.     

On model equations for particle dispersion in inhomogeneous turbulence

Comparisons are made between the Advection–Diffusion Equation (ADE) approach for particle transport and the two-fluid model approach based on the PDF method. In principle, the ADE approach offers a much simpler way of calculating the inertial deposition of particles in a turbulent boundary layer than that based on the PDF approach. However the ADE equations that have recently been used are only strictly valid for a simple Gaussian process when particle inertia is small. Using a prescribed, but in general non-Gaussian random particle velocity field, it is shown that the net particle mass flux contains a drift term in addition to that from the mean velocity of the particle velocity field, associated with the compressibility of the velocity field. Furthermore the diffusive flux in general depends not only upon the gradient of the mean concentration (true only for a Gaussian random flow field) but also upon higher order derivatives whose relative contribution depends on diffusion coefficients D_ijk… etc. These coefficients depend upon the statistical moments associated with random displacements and compressibility of the particle flow field along particle trajectories which in turn depend upon particle inertia. In contrast the PDF approach offers the advantage of using a simple gradient (Gaussian) approximation in particle phase space which can lead to a non-Gaussian spatial dispersion process when particle inertia is important. Conditions based on the particle mean free path are derived for which a simple ADE is appropriate. Some of the features of particle transport in an inhomogeneous turbulent flow are illustrated by examining particle dispersion in a random flow field composed of pairs of counter rotating vortices which has an rms velocity which increase linearly from a stagnation point.