大气边界层内羽流扩散研究

主要研究了高架点源的污染,即``羽流扩散'. 由于大气边界层内的湍流运动,引 起污染物扩散的主要因素是湍流扩散, 研究湍流扩散有两种基本方法：统计理论和梯度理论. 采用拉格朗日粒子扩散模型 (Lagrangian particle dispersion model, LPDM)来评价羽流扩散的平均浓度,分别计算了表征扩散程度的3个参数：沿烟轴方向的 着地浓度、水平方向的羽流扩散宽度和垂直方向的羽流扩散宽度. 从计算结果和风洞实验的 数据对比来看,用该模式可以较好地模拟大气边界层内复杂气候条件下的羽流扩散.     

Pollution dispersion analysis using the puff model with numerical flow field data

The computations of the flowfield and pollutant dispersion over a flat plate and the Russian hills of various slopes are described. The Gaussian plume and the puff model have been used to calculate concentration of pollutant. The Reynolds-averaged unsteady incompressible Navier–Stokes equation with low Reynolds k– model has been used to calculate the flowfield. The flow data of a flat plate and the Russian hills from Navier–Stokes equation solutions has been used as the input data for the puff model. The computational results of flowfield agree well with experimental results of both a flat plate and Russian hills. The concentration prediction by the Gaussian plume model and the Gaussian puff model also agrees fairly well with experiments.     

A comparison between a plume path model and a virtual point source model for a stack plume

Some years ago we developed a theoretical model for the calculation of the rise of a stack plume in the atmosphere. Recently we have extended this model in such a way that the ground level concentration of gases emitted by a stack can also be calculated. The cross-section of the plume, which was originally assumed to be a circle, is now assumed to be an ellipse. The height-to-width ratio of the ellipse is chosen in accordance with full-scale data. The entrainment of air into the plume due to atmospheric turbulence can now be calculated not only for neutral atmospheric conditions, but also for stable and unstable conditions.Predictions of the ground level concentration made with this model have been compared with predictions made with the often used virtual point source model. For neutral and unstable atmospheric conditions the agreement is rather good. For stable conditions the differences can be rather large. The sensitivity of the ground level concentration with respect to the atmospheric stability condition, the stack height, the plume exit momentum flux and the plume exit buoyancy flux has also been investigated with the aid of the model. The influence of the plume exit momentum flux is negligible; however, the influence of the other parameters can be large.     

A fast non‐Fickian particle‐tracking diffusion simulator and the effect of shear on the pollutant diffusion process

This paper presents a fast method for the generation of non‐Fickian particle paths within a particle‐tracking pollutant diffusion model based on a Fourier spectral representation of fractional Brownian motion (fBm), a generalization of ordinary Brownian motion. Correlated diffusive components in a particle‐tracking algorithm are modelled using fBm increments that have long‐range correlations over numerous spatial and/or temporal scales; hence producing non‐Fickian diffusion. A fast algorithm to generate fBm and its increment by using its power spectral density S(f) in a fast Fourier transform algorithm is given. A general equation for the scaling of fBm within a velocity flow field with simple linear shear is presented. An initial numerical study of the nature of fBm shear dispersion has been conducted by incorporating fBm increments into a non‐Fickian particle‐tracking algorithm. It is shown that the effect of simple (i.e. linear) shear on the diffusion process is to produce enhanced diffusive phenomena with the longitudinal spreading of the plume scaling with exponent ∼1+H, where H is the Hurst exponent used to describe fBm. Finally, a more complex shear zone at the entrance of a coastal bay model is investigated using both a traditional particle‐tracking method and the fBm‐based method. Copyright © 2000 John Wiley & Sons, Ltd.     

Large eddy simulations of a turbulent thermal plume

Large eddy simulations of a three-dimensional turbulent thermal plume in an open environment have been carried out using a self-developed parallel computational fluid dynamics code SMAFS (smoke movement and flame spread) to study the thermal plume’s dynamics including its puffing, self-preserving and air entrainment. In the simulation, the sub-grid stress was modeled using both the standard Smagorinsky and the buoyancy modified Smagorinsky models, which were compared. The sub-grid scale (SGS) scalar flux in the filtered enthalpy transport equation was modeled based on a simple gradient transport hypothesis with constant SGS Prandtl number. The effect of the Smagorinsky model constant and the SGS Prandtl number were examined. The computation results were compared with experimental measurements, thermal plume theory and empirical correlations, showing good agreement. It is found that both the buoyancy modification and the SGS turbulent Prandtl number have little influence on simulation. However, the SGS model constant C _s has a significant effect on the prediction of plume spreading, although it does not affect much the prediction of puffing.     

大涡模拟建筑物对近源大气污染物扩散的影响

采用大涡模拟(LES)方法对一个由上风低源排放的烟流在孤立建筑物周围扩散的现象进行了数值模拟。模拟中,离散网格的尺度控制在100倍Kolmogorov湍流尺度的量级以内,亚格子应力由Smagorinsky模型加以模拟,并采用了随机数方法生成符合指定湍流特征的入流边界条件。数值模拟结果表明:建筑物尾流区平均浓度场的分布和大小不仅取决于局地流场的影响,在很大程度上也受到因建筑物阻挡出现的烟流抬升和马蹄涡对烟流的捕获的影响。通过对比风洞试验数据显示:大涡模拟方法能较好地再现气流绕建筑物的流动形态和建筑物周围的浓度分布,因而具有广泛的应用前景。     

饱和非饱和土壤中溶质运移的数值模拟

本文结合拉格朗日观点的运动坐标与欧拉观点的固定坐标,分别求解对流—弥敌方程的对流效应和弥散效应,并在从固定坐标向运动坐标投影时利用了三次样条插值方法,使模型更为精确而简单。通过与特定情况下的理论解比较,说明模型是准确可靠的。通过与实验室室内试验结果的比较,模型的可靠性得到了进一步的验证。利用模型对淡水洗盐结合浅群井抽排咸水治理盐碱地过程的模拟分析,认为结合浅群井抽排可较大地提高淡水洗盐的效果。     

A Fickian diffusion model for the spreading of liquid plumes infiltrating in heterogeneous media

Infiltration of water and non-aqueous phase liquids (NAPLs) in the vadose zone gives rise to complex two- and three-phase immiscible displacement processes. Physical and numerical experiments have shown that ever-present small-scale heterogeneities will cause a lateral broadening of the descending liquid plumes. This behavior of liquid plumes infiltrating in the vadose zone may be similar to the familiar transversal dispersion of solute plumes in single-phase flow. Noting this analogy we introduce a mathematical model for ‘phase dispersion’ in multiphase flow as a Fickian diffusion process. It is shown that the driving force for phase dispersion is the gradient of relative permeability, and that addition of a phase-dispersive term to the governing equations for multiphase flow is equivalent to an effective capillary pressure which is proportional to the logarithm of the relative permeability of the infiltrating liquid phase. The relationship between heterogeneity-induced phase dispersion and capillary and numerical dispersion effects is established. High-resolution numerical simulation experiments in heterogeneous media show that plume spreading tends to be diffusive, supporting the proposed convection-dispersion model. Finite difference discretization of the phase-dispersive flux is discussed, and an illustrative application to NAPL infiltration from a localized source is presented. It is found that a small amount of phase dispersion can completely alter the behavior of an infiltrating NAPL plume, and that neglect of phase-dispersive processes may lead to unrealistic predictions of NAPL behavior in the vadose zone.     

Investigation of use of reach-back characteristics method for 2D dispersion equation

The use of the Holly-Preissmann two-point scheme has been very popular for the calculation of the dispersion equation. The key to this scheme is to use the characteristics method incorporating the Hermite cubic interpolation technique to approximate the trajectory foot of the characteristics. This method can avoid the excessive numerical damping and oscillation associated with most finite difference schemes for advection computation. On the basis of the fundamental idea of the Holly-Preissmann two-point scheme, a new technique is introduced herein for the computation of the two-dimensional dispersion equation. This new scheme allows the characteristics projecting back several time steps to fall on the spatial or temporal axis, while the characteristics foot is still solved by the Holly-Preissmann two-point method. The diffusion portion of the dispersion equation is solved by the commonly used Crank-Nicholson method. The calculation for these two processes consisting of advection and diffusion is carried out separately but consecutively in one time step, a method known as the split operator algorithm. A hypothetical model was constructed to demonstrate the applicability of this new technique for the calculation of the pure advection and dispersion equation in two dimensions.     

Large eddy simulation of aircraft wake vortex with self-adaptive grid method

A self-adaptive-grid method is applied to numerical simulation of the evolution of aircraft wake vortex with the large eddy simulation(LES). The Idaho Falls(IDF)measurement of run 9 case is simulated numerically and compared with that of the field experimental data. The comparison shows that the method is reliable in the complex atmospheric environment with crosswind and ground effect. In addition, six cases with different ambient atmospheric turbulences and Brunt V¨ais¨al¨a(BV) frequencies are computed with the LES. The main characteristics of vortex are appropriately simulated by the current method. The onset time of rapid decay and the descending of vortices are in agreement with the previous measurements and the numerical prediction. Also, secondary structures such as baroclinic vorticity and helical structures are also simulated.Only approximately 6 million grid points are needed in computation with the present method, while the number can be as large as 34 million when using a uniform mesh with the same core resolution. The self-adaptive-grid method is proved to be practical in the numerical research of aircraft wake vortex.     

A method to compute ship exhaust plumes with waves and wind

Thermal and concentration transport models are implemented in CFDShip‐Iowa version 4.5, a semi‐coupled solver for air/water free surface flow (Int. J. Numer. Meth. Fluids 2008; 58 (6):591–624), to investigate the exhaust plume around ship superstructures. An incompressible, variable density approximation is implemented where the density can change in all governing equations due to temperature variations only. The thermal and concentration models are tested for the cases of steady and unsteady flow with thermal and solution transport in a 2D square cavity, and for a 3D thermal plume in an open environment, showing good agreement between computational results and experimental data. To test the method in an extreme motions condition, the exhaust plume of the ONR Tumblehome model DTMB 5613 is studied, showing complicated vortical structures in air including a pair of counter‐rotating vortices downstream of the stack for cross‐flow, and bended bird‐plume shape in the symmetry plane and varying arc‐shape in axial sections both for temperature and NO_x concentration fields. Effects of smoke exhaust speed and wind speed on the temperature and concentration distributions are studied. Finally, a smoke downwash computation is performed for a ship free to move in 6 degrees of freedom in a sea state 8 condition. Copyright © 2010 John Wiley & Sons, Ltd.     

AN EFFICIENT SOLUTION FOR 2D RAYLEIGH-BÉNARD CONVECTION USING FFT

研究提高二维方腔瑞利-贝纳德对流 直接数值模拟求解方法的计算效率问题.对于非定常湍流热对流, 压力泊松方程的求解是影响整个计算效率的关键. 利用快速傅里叶变换(fast Fourier transform,FFT)解耦并结合追赶法, 可实现压力泊松方程的直接求解.通过与跳点超松弛迭代法在求解精度和计算速度对比, 可以看到, 利用FFT压力泊松方程直接方法计算热对流问题是高效的.还给出了典型状态的热对流初始羽流和大尺度环流温度场, 以及系列瑞利数(Ra)计算结果的宏观传热努塞数(Nu)变化.     

The numerical analysis of water and solute movement in scale heterogeneous profiles

A computer based numerical method is presented for the analysis of water and solute movement in unsaturated heterogeneous porous materials. Such a method is necessary since, for those field studies where solute movement is of concern, the soil profiles under consideration are invariably heterogeneous. The numerical analysis is based on a general one-dimensional finite difference soil water flow model which includes a numerical technique combining the concepts of scale heterogeneity with an interpolative soil water hysteresis model. An explicit finite difference solute movement subroutine is incorporated into the unsaturated flow model to describe the transport of nonreactive solutes. A velocity dependent longitudinal dispersion coefficient is used in the solution of the hydrodynamic dispersion equation. The resulting hysteretic scale heterogeneous solute movement model permits the study of solute dynamics during infiltrating and redistribution in realistically complex spatially varying soil profiles. Results are presented for the leaching of both coarse grading to fine and fine grading to coarse sand profiles. Both vertical and horizontal profiles are studied using either a constant flux or a constant concentration input boundary condition. The four cases studied demonstrate the versatility of the numerical method and emphasise the substantial differences in transport behavior that can arise between heterogeneous and homogeneous profiles.Now with BHP Petroleum Pty. Ltd., GPO Box 1911R, Melbourne, Vic. 3001, Australia.     

Dispersion error reduction for acoustic problems using the smoothed finite element method (SFEM)

The smoothed finite element method (SFEM), which was recently introduced for solving the mechanics and acoustic problems, uses the gradient smoothing technique to operate over the cell‐based smoothing domains. On the basis of the previous work, this paper reports a detailed analysis on the numerical dispersion error in solving two‐dimensional acoustic problems governed by the Helmholtz equation using the SFEM, in comparison with the standard finite element method. Owing to the proper softening effects provided naturally by the cell‐based gradient smoothing operations, the SFEM model behaves much softer than the standard finite element method model. Therefore, the SFEM can significantly reduce the dispersion error in the numerical solution. Results of both theoretical and numerical experiments will support these important findings. It is shown clearly that the SFEM suits ideally well for solving acoustic problems, because of the crucial effectiveness in reducing the dispersion error. Copyright © 2015 John Wiley & Sons, Ltd.     

复杂大气入流下海上风机力学特性研究

随着风能技术的不断进步, 风机叶片逐渐向大型化发展, 这使得真实复杂大气入流对风机运行性能的影响愈发显著. 为研究真实复杂大气入流下海上风机的力学特性响应, 利用基于大涡模拟的域前模拟方法生成复杂大气入流, 并结合致动线模型模拟风机叶片, 对中性复杂大气入流下海上固定式风机进行数值模拟, 重点分析风机的气动性能及转子和叶片根部的力学特性, 并与均匀入流计算工况进行对比. 计算结果表明, 中性复杂大气入流中的大尺度低速气流团使得风机气动功率输出值在较长一段时间处于较低水平, 此外, 中性复杂大气入流的高湍流强度特征使得风机气动功率的变化幅值和标准差较均匀入流工况大幅增加; 风机轴向推力的标准差值增加到均匀入流的53倍, 中性复杂大气入流的来流流场扰动引起偏航力矩的最大值、均方根和标准差分别增加到均匀入流的10、4.4和4.3倍; 速度垂向分布的不均匀性以及轮毂高度附近的大尺度低速羽流结构导致摆振剪力和弯矩的标准差响应值分别为均匀入流的2倍和4.6倍.      

Finite element solution to passive scalar transport behind line sources under neutral and unstable stratification

This study employed a direct numerical simulation (DNS) technique to contrast the plume behaviours and mixing of passive scalar emitted from line sources (aligned with the spanwise direction) in neutrally and unstably stratified open‐channel flows. The DNS model was developed using the Galerkin finite element method (FEM) employing trilinear brick elements with equal‐order interpolating polynomials that solved the momentum and continuity equations, together with conservation of energy and mass equations in incompressible flow. The second‐order accurate fractional‐step method was used to handle the implicit velocity–pressure coupling in incompressible flow. It also segregated the solution to the advection and diffusion terms, which were then integrated in time, respectively, by the explicit third‐order accurate Runge–Kutta method and the implicit second‐order accurate Crank–Nicolson method. The buoyancy term under unstable stratification was integrated in time explicitly by the first‐order accurate Euler method. The DNS FEM model calculated the scalar‐plume development and the mean plume path. In particular, it calculated the plume meandering in the wall‐normal direction under unstable stratification that agreed well with the laboratory and field measurements, as well as previous modelling results available in literature. Copyright © 2005 John Wiley & Sons, Ltd.     

Dynamics of large solid particles in bioconvective sedimentation

Settling of one or two large solid particles in a bioconvection flow induced by gyrotactic motile microorganisms is investigated using a 2D numerical model. The results of varying the initial positions of large particles on the bioconvection flow pattern are investigated. The Chimera method is utilized to generate subgrids around the moving particles. It is demonstrated that the introduction of a single large particle displaces bioconvection plume and changes its shape. The introduction of two particles on the same side of the bioconvection plume further displaces the plume while the introduction of two particles on opposite sides reduces this displacement. The influence of the bioconvection plume on the particles' settling paths and particles' settling velocities is investigated. Copyright © 2006 John Wiley & Sons, Ltd.     

A Lagrangian vortex method for unbounded flows

A technique is presented for velocity calculations on the highly distorted node distributions typical of those found in Lagrangian vortex methods. The method solves the partial differential equation for streamfunction directly on the nodes, via a sparse, symmetric system of equations that can be solved using standard iterative solvers. When implemented in a triangulated vortex method, the technique gives computation times which scale as N^1.23, where N is the number of nodes. The computation scheme is derived for two‐dimensional problems and applied to the prediction of the evolution of perturbed multipolar vortices. Due to the numerical performance of the method, it has been possible to examine such evolution at higher and lower Reynolds numbers than have been considered in published numerical studies. Copyright © 2008 John Wiley & Sons, Ltd.     

Numerical analysis and modeling of plume meandering in passive scalar dispersion downstream of a wall-mounted cube

A DNS database is employed to examine the onset of plume meandering downstream of a wall-mounted cube and to address the impact of large-scale unsteadiness in modeling dispersion using the RANS equations. The cube is immersed in a uniform stream where the thin boundary-layer developing over the flat plate is responsible for the onset of vortex-shedding in the wake of the bluff-body. Spectra of velocity and concentration fluctuations exhibit a prominent peak in the energy content at the same frequency, showing that the plume meandering is established by the action of the vortex-shedding. The vortex-shedding and plume meandering display a low-frequency modulation where coherent fluctuations are suppressed at times with a quasi-regular period. The onset of the low-frequency modulation is indicated by a secondary peak in the energy spectrum and confirmed by the autocorrelation of velocity and scalar fluctuations. Unsteady RANS simulations performed with the v² − f model are able to detect the onset of the plume meandering and show remarkable improvement of the predicted decay rate and rate of spread of the scalar plume when compared to steady RANS solutions. By computing explicitly the periodic component of velocity and scalar fluctuations, the unsteady v² − f model is able to provide a representation of scalar flux components consistent with DNS statistics, where the counter-gradient transport mechanism that takes place in the streamwise component is also captured by URANS results. Nonetheless, the agreement with DNS statistics for the mean concentration and the plume width is limited by the onset of the low-frequency modulation in the vortex-shedding and plume meandering, giving a challenging modeling issue in the simulation of dispersion using the RANS equations.     

A FINITE DIFFERENCE METHOD FOR SIMULATING TRANSVERSE VIBRATIONS OF AN AXIALLY MOVING VISCOELASTIC STRING

A finite difference method is presented to simulate transverse vibrations of an axially moving string.By discretizing the governing equation and the equation of stress- strain relation at different frictional knots,two linear sparse finite difference equation systems are obtained.The two explicit difference schemes can be calculated alternatively, which make the computation much more efficient.The numerical method makes the nonlinear model easier to deal with and of truncation errors,O(Δt~2 Δx~2).It also shows quite good stability for small initial values.Numerical examples are presented to demonstrate the efficiency and the stability of the algorithm,and dynamic analysis of a viscoelastic string is given by using the numerical results.