饱和-非饱和土壤中污染物运移过程的数值模拟

本文提出了一个模拟饱和 非饱和土壤中溶和污染物运移过程的数值模型．模拟的控制污染物运移的物理 化学现象包括：对流，机械逸散，分子弥散，吸附，蜕变，不动水效应．发展了一个修正的特征线Ｇａｌｅｒｋｉｎ方法以离散污染物运移过程的控制方程并导出了一个用于有限元方程求解的显式算法．数值例题结果表明所提出模型和算法的功能     

An optical flow MTV based technique for measuring microfluidic flow in the presence of diffusion and Taylor dispersion

A novel technique is presented for accurately measuring flow fields in microfluidic flows from molecular tagging velocimetry (MTV). Limited optical access is frequently encountered in microfluidic systems. Therefore, in this contribution we analyze the special case of tagging a line across the thin dimension of a microchannel and subsequent imaging along this line. This represents a set-up that is applicable to a wide range of microfluidic applications. A volume illumination has to be used, resulting in an integration of the visualized dye across the flow profile. This leads to the well-known effect of Taylor dispersion. Our novel technique consists of measuring motion from digital image sequences in a gradient-based approach. A motion model is developed which explicitly deals with brightness changes due to Taylor dispersion and additional molecular diffusion of dyes. The presented approach is specific to the case of a parabolic velocity profile. In the presence of these effects, an accurate computation of motion is only possible by applying this novel motion model. Our technique is tested on simulated sequences corrupted with varying levels of noise and on actual measurements. Measurements were conducted in a microfluidic mixer of precisely known flow properties. In the same mixer, a comparative study of our MTV technique to μPIV was performed. Also, the results were compared to bulk measurements of the fluid flow velocity. The novel algorithm compared favorably and also, measurements were conducted on inhomogeneous flow configurations.     

A method for computing compressible,highly stratified flows in astrophysics based on operator splitting

A new numerical approach based on consistent operator splitting is presented for computing compressible, highly stratified flows in astrophysics. The algorithm is particularly designed to search for steady or almost steady solutions for the time-dependent Navier–Stokes equations, describing viscous flow under the influence of a strong gravitational field. The algorithm proposed is multidimensional and works in Cartesian, cylindrical or spherical co-ordinates. It uses a second-order finite volume scheme with third-order upwinding and a second-order time discretization. An adaptive time step control and monotonic multilevel grid distribution has been incorporated to speed up convergence. This method has been incorporated into a hydrodynamical code by which, for the first time, for two-dimensional models the dynamics of the boundary layer in the accretion disk around a compact star could be computed over the whole viscous time scale. © 1998 John Wiley & Sons, Ltd.     

非定常流动的快速拉格朗日涡方法数值模拟

采用快速拉格朗日涡方法数值模拟有复杂旋涡运动的非定常流动. 利用离散涡元模拟旋涡的产生、聚集和输送过程. 拉格朗日描述法用来计算离散涡元的移动,而移动速度则利用广义毕奥-萨伐尔公式结合快速多极子展开法计算,修正的涡半径扩散模型用来模拟离散涡元的黏性扩散. 突然起动圆柱和大攻角下突然起动翼型的非定常有涡流动的数值模拟,及其与试验结果的对比验证了方法的有效性. 另外,大攻角下突然起动翼型的计算结果给出了翼型起动后吸力面旋涡的产生、发展,周期性非定常流动的形成,以及尾流旋涡结构等一些重要的流动特征.[关键词] 非定常流有涡流动快速涡方法      

A multiblock Navier–Stokes algorithm using equal-order quadratic finite elements

A new multiblock pressure-based finite element algorithm has been developed. This methodology implements quadratic interpolation for both the elemental velocity and pressure fields. A direct streamline upwinding scheme previously developedby the authors is used to model the non-linear inertia effects. Details of the algorithm and its multiblock foundation are provided along with validating test cases. The results presented clearly demonstrate the accuracy of this new approach and the differences in the pressure field for an element using quadratic versus the traditional bi linear approximation of the pressure field.     

Turbulent transport dissimilarity with modulated turbulence structure in channel flow of viscoelastic fluid

This experimental study investigated the turbulent transport dissimilarity with a modulated turbulence structure in a channel flow of a viscoelastic fluid using simultaneous particle image velocimetry and planar laser-induced fluorescence measurements. An instantaneous dye concentration field with fluctuating velocity vectors showed that mass was transferred by hierarchically large-scale wavy motions with inclination. A co-spectral analysis showed that the spatial phase modulation of the streamwise velocity and dye concentration fluctuations for the wall-normal velocity fluctuation corresponded to the relaxation time. The occurrence of intense dye concentration fluctuation and small streamwise velocity fluctuation in a thin boundary layer caused dissimilar turbulent transport because of the non-zero negative correlation of the streamwise velocity and dye concentration fluctuations for the wall-normal velocity fluctuation only on large scales. This explains the turbulent transport dissimilarity which leads to the zero averaged Reynolds shear stress and non-zero wall-normal turbulent mass flux.     

Numerical modelling of 3D stratified free surface flows: a case study of sediment dumping

A three‐dimensional numerical model has been developed to simulate stratified flows with free surfaces. The model is based on the Reynolds‐averaged Navier–Stokes (RANS) equations with variable fluid density. The equations are solved in a transformed σ‐coordinate system with the use of operator‐splitting method (Int. J. Numer. Meth. Fluids 2002; 38 :1045–1068). The numerical model is validated against the one‐dimensional diffusion problem and the two‐dimensional density‐gradient flow. Excellent agreements are obtained between numerical results and analytical solutions. The model is then used to study transport phenomena of dumped sediments into a water body, which has been modelled as a strongly stratified flow. For the two‐dimensional problem, the numerical results compare well with experimental data in terms of mean particle falling velocity and spreading rate of the sediment cloud for both coarse and medium‐size sediments. The model is also employed to study the dumping of sediments in a three‐dimensional environment with the presence of free surface. It is found that during the descending process an annulus‐like cloud is formed for fine sediments whereas a plate‐like cloud for medium‐size sediments. The model is proven to be a good tool to simulate strongly stratified free surface flows. Copyright © 2005 John Wiley & Sons, Ltd.     

Numerical viscosity reduction in the resolution of the shallow water equations with turbulent term

A first‐order finite volume model for the resolution of the 2D shallow water equations with turbulent term is presented. An upwind discretization of the equations that include the turbulent term is carried out. A method to reduce the excess of numerical viscosity (or diffusion) produced by the upwinding of the flux term is proposed. Two different discretizations of the turbulent term are compared, and results for uniform distributions of the viscosity are presented. Finally, two discretizations of the time derivative which are more efficient than Euler's are proposed and compared. Copyright © 2008 John Wiley & Sons, Ltd.     

A semi‐implicit scheme for large Eddy simulation of piston engine flow and combustion

A semi‐implicit scheme is presented for large eddy simulation of turbulent reactive flow and combustion in reciprocating piston engines. First, the governing equations in a deforming coordinate system are formulated to accommodate the moving piston. The numerical scheme is made up of a fourth‐order central difference for the diffusion terms in the transport equations and a fifth‐order weighted essentially nonoscillatory (WENO) scheme for the convective terms. A second‐ order Adams–Bashforth scheme is used for time integration. For higher density ratios, it is combined with a predictor–corrector scheme. The numerical scheme is explicit for time integration of the transport equations, except for the continuity equation which is used together with the momentum equation to determine the pressure field and velocity field by using a Poisson equation for the pressure correction field. The scheme is aimed at the simulation of low Mach number flows typically found in piston engines. An efficient multigrid method that can handle high grid aspect ratio is presented for solving the pressure correction equation. The numerical scheme is evaluated on two test engines, a laboratory four‐stroke engine with rectangular‐shaped engine geometry where detailed velocity measurements are available, and a modified truck engine with practical cylinder geometry where lean ethanol/air mixture is combusted under a homogeneous charge compression ignition (HCCI) condition. Copyright © 2012 John Wiley & Sons, Ltd.     

A finite difference self-adaptive mesh solution of flow in a sedimentation tank

A new numerical model has been developed to evaluate the removal efficiency of primary sedimentation clarifiers operating at neutral density condition. The velocity and concentration fields as well as the development in time and space of the settled particle bed thickness are simulated. The main difficulties in simulation of velocity and concentration fields are related to (1) numerical instabilities produced by the prevalence of convective terms in the unknown variable high-gradient regions and (2) turbulence effects on the suspension of solid particles from the settled bed. The need to overcome the numerical instabilities without the upwind difference approximation, which introduces high numerical viscosity, suggests the use of non-uniform grids of calculation. The velocity field is obtained by solving the motion equations in the vorticity and streamfunction formulation by means of a new numerical method based upon a dynamically self-adjusting calculation grid. These grids allow for a finer mesh following the evolution of the unknown quantities. A k–? model is used to simulate turbulence phenomena. The sedimentation field is found by solving the diffusion and transport equation of the solid particle concentration. Boundary conditions on the bottom line are imposed relating the amount of turbulence flux and sedimentation flux to the actual concentration and the reference concentration. Such an approach makes it possible to represent the solid particle suspension from the bottom, taking into account its dependence on (1) the characteristics and the evolution in time of the settled bed, (2) the velocity component parallel to the bottom line and (3) the turbulence structure.     

城市和工业区大气环境评估与流动实验--风环境模拟研究之二

借用实验流体力学分层流实验技术，采用均匀的水模拟中性大气，表示大气随高度没有变化，再用盐水进行密度分层来模拟大气温度随高度变化的稳定层结大气，丰富了大气环境的实验研究手段，本介绍用水来模拟大气的相似理论和实验研究方法，针对北京市地形对城市大气环境影响变化，来解释和说明首都在解决城市环境污染治理上应该关注的问题。同时结合福建省某核电厂厂址地区受地形影响，分析在厂运行期间排出的放射性气载释放物的大气输运和扩散影响，以及提供突发性污染事故应急预警方案，采用现场实测、水槽模拟和数值计算等3种手段，开展研究并取得成功、     

含颗粒污染物的油液均质流的数值分析

含颗粒污染物的油液是浓度很稀、粒径极小的伪均质流,为了掌握颗粒污染物在输送过程中的浓度分布,利用一维扩散方程构建了污染油液的数学模型;通过特征线法数值求解,获得了污染油液中各相的动态特征。结果表明,油液压力和速度沿管长呈脉动规律运动,且随着时间的延长逐渐衰减;颗粒污染物对油液速度具有极好的跟随特性;颗粒污染物的浓度分布也随着油液流速的变化而呈现规律性的变化;在不同运行时间内油液压力沿管长的衰减趋势不同,油液速度沿管长的变化趋势与压力的趋势相反;颗粒污染物速度和浓度分布沿管长与油液速度具有紧密的联系。     

Mie-Grneisen状态方程可压缩多流体流动的PPM方法

采用流体体积分数的混合型多流体数值模型，将piecewise parabolic method (PPM)方法应用于可压缩多流体流动的数值模拟，拓展了以前提出的模型和数值方法，使它能够处理一般的Mie-Grneisen状态方程。采用双波近似和两层迭代算法求解一般状态方程的Riemann问题；并根据多流体接触界面无振荡原则设计高精度计算格式，对典型的纯界面平移问题可以从理论上证明本算法在接触间断附近压力和速度没有振荡，而且数值模拟结果表明界面数值耗散也被控制在2~3个网格之内。模拟了多种复杂的可压缩多流体流动，算例结果表明本文方法可以有效地处理接触间断、激波等物理问题，且具有耗散小精度高的特点。     

A coupled multigrid-domain-splitting technique for simulating incompressible flows in geometrically complex domains

This paper describes a domain decomposition numerical procedure for solving the Navier-Stokes equations in regions with complex geometries. The numerical method includes a modified version of QUICK (quadratic upstream interpolation convective kinematics) for the formulation of convective terms and a central difference scheme for the diffusion terms. A second-order-accurate predictor-corrector scheme is employed for the explicit time stepping. Although the momentum equations are solved independently on each subdomain, the pressure field is computed simultaneously on the entire flow field. A multigrid technique coupled with a Schwarz-like iteration method is devised to solve the pressure equation over the composite domains. The success of this strategy depends crucially on appropriate methods for specifying intergrid pressure boundary conditions on subdomains. A proper method for exchanging information among subdomains during the Schwarz sweep is equally important to the success of the multigrid solution for the overall pressure field. These methods are described and subsequently applied to two forced convection flow problems involving complex geometries to demonstrate the power and versatility of the technique. The resulting pressure and velocity fields exhibit excellent global consistency. The ability to simulate complex flow fields with this method provides a powerful tool for analysis and prediction of mixing and transport phenomenon.     

Finite element technique to solve the elastic strain for Leonov fluid flow

The finite element method is used to find the elastic strain (and thus the stress) for given velocity fields of the Leonov model fluid. With a simple linearization technique and the Galerkin formulation, the quasi-linear coupled first-order hyperbolic differential equations together with a non-linear equality constraint are solved over the entire domain based on a weighted residual scheme. The proposed numerical scheme has yielded efficient and accurate convective integrations for both the planar channel and the diverging radial flows for the Leonov model fluid. Only the strain in the inflow plane is required to be prescribed as the boundary conditions. In application, it can be conveniently incorporated in an existing finite element algorithm to simulate the Leonov viscoelastic fluid flow with more complex geometry in which the velocity field is not known a priori and an iterative procedure is needed.     

A physically based flux limiter for QUICK calculations of advective scalar transport

Transient, advective transport of a contaminant into a clean domain will exhibit a moving sharp front that separates contaminated and clean regions. Due to ‘numerical diffusion’—the combined effects of ‘cross‐wind diffusion’ and ‘artificial dispersion’—a numerical solution based on a first‐order (upwind) treatment will smear out the sharp front. The use of higher‐order schemes, e.g. QUICK (quadratic upwinding) reduces the smearing but can introduce non‐physical oscillations in the solution. A common approach to reduce numerical diffusion without oscillations is to use a scheme that blends low‐order and high‐order approximations of the advective transport. Typically, the blending is based on a parameter that measures the local monotonicity in the predicted scalar field. In this paper, an alternative approach is proposed for use in scalar transport problems where physical bounds C_Low?C?C_High on the scalar are known a priori. For this class of problems, the proposed scheme switches from a QUICK approximation to an upwind approximation whenever the predicted upwind nodal value falls outside of the physical range [C_Low, C_High]. On two‐dimensional steady‐state and one‐dimensional transient test problems predictions obtained with the proposed scheme are essentially indistinguishable from those obtained with monotonic flux‐limiter schemes. An analysis of the modified equation explains the observed performance of first‐ and second‐order time‐stepping schemes in predicting the advective transport of a step. In application to the transient two‐dimensional problem of contaminate transport into a streambed, predictions obtained with the proposed flux‐limiter scheme agree with those obtained with a scheme from the literature. Copyright © 2007 John Wiley & Sons, Ltd.     

Optically sliced measurement of velocity and pH distribution in microchannel

A simultaneous measurement technique for the velocity and pH distribution was developed by using a confocal microscope and a 3CCD color camera for investigations of a chemical reacting flow field in a microchannel. Micron-resolution particle image velocimetry and laser induced fluorescence were utilized for the velocity and pH measurement, respectively. The present study employed fluorescent particles with 1 μm diameter and Fluorescein sodium salt whose fluorescent intensity increases with an increase in pH value over the range of pH 5.0–9.0. The advantages of the present system are to separate the fluorescence of particles from that of dye by using the 3CCD color camera and to provide the depth resolution of 5.0 μm by the confocal microscope. The measurement uncertainties of the velocity and pH measurements were estimated to be 5.5 μm/s and pH 0.23, respectively. Two aqueous solutions at different pH values were introduced into a T-shaped microchannel. The mixing process in the junction area was investigated by the present technique, and the effect of the chemical reaction on the pH gradient was discussed by a comparison between the proton concentration profiles obtained from the experimental pH distribution and those calculated from the measured velocity data. For the chemical reacting flow with the buffering action, the profiles from the numerical simulation showed smaller gradients compared with those from the experiments, because the production or extinction of protons was yielded by the chemical reaction. Furthermore, the convection of protons was evaluated from the velocity and pH distribution and compared with the diffusion. It is found that the ratio between the diffusion and convection is an important factor to investigate the mixing process in the microfluidic device with chemical reactions.     

稠密气固两相流各向异性颗粒相矩方法

基于气体分子动力学和颗粒动理学方法,考虑颗粒速度脉动各向异性,建立颗粒相二阶矩模型.应用初等输运理论,对三阶关联项进行模化和封闭.考虑颗粒与壁面之间的能量传递和交换,建立颗粒相边界条件模型.采用Koch等计算方法模拟气固脉动速度关联矩.考虑气体-颗粒间相互作用,建立稠密气体-颗粒流动模型.数值模拟提升管内气固两相流动特性,模拟结果表明提升管内颗粒相湍流脉动具有明显的各向异性.预测颗粒速度、浓度和颗粒脉动速度二阶矩与Tartan等实测结果相吻合.模拟结果表明轴向颗粒速度脉动强度约为平均颗粒相脉动强度的1.5倍,轴向颗粒脉动能大约是径向颗粒脉动能3.0倍.     

Transport with a Very Low Density Contrast in Hele–Shaw Cell and Porous Medium: Evolution of the Mixing Zone

The optimal concentration of a blue dye solution with 'tracer' properties, enabling a pollutant to be marked was determined by the use of numerical, theoretical and experimental approaches. Experimental investigations were performed on a transparent Hele–Shaw cell and the concentration distribution was analyzed using an optical technique based on dye light absorption properties. The injected optimal concentration was established thanks to a theoretical and experimental study carried out on the output signal dynamics. Using the same experimental conditions, numerical simulations were performed. The very good agreement between the data (experimental and numerical) clarified that: (i) the choice of the blue dye optimal concentration was valid and (ii) the concentration-dependent density should not be neglected in flow and transport equations even if it concerns a so-called 'tracer'. Following this remark, a theoretical aspect was developed in order to determine the analogous conditions between a Hele–Shaw cell and a porous medium for the variable density transport phenomenon. The structure of the concentration-dependent dispersion tensor used in the numerical code was obtained by homogenizing the Stokes flow of a bi-component mixture. The numerical results show that, as long as the tracer density does not exceed a certain value, it is not necessary to take into account a density contrast in terms of the dispersion tensor. The classical form of the Taylor dispersion tensor can be used successfully.