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.     

DNS study of passive plume interference emitting from two parallel line sources in a turbulent channel flow

Turbulent mixing of dual plumes emitting simultaneously from line sources in a turbulent channel flow has been studied using direct numerical simulation (DNS). Three test cases have been compared to investigate the effects of the source separation on turbulent mixing of the two instantaneous plumes. The dispersion and interference of dual plumes are investigated in both physical and spectral spaces, which include an analysis of statistical moments of the concentration field, cross-correlation between the two instantaneous plumes, pre-multiplied spectra of the velocity and concentration fields, and co-spectrum and coherency spectrum of the dual plumes. As the downstream distance from the line source increases, the plume development associated with a single source emission transitions from a turbulent convective stage to a turbulent diffusive stage. It is observed that a plume released from a ground-level source reaches the turbulent diffusive stage faster than that released from an elevated source. It is also observed that a smaller separation between the two line sources tends to facilitate a more rapid growth in the cross-correlation coefficient of two instantaneous plumes. In the near-source region, the maximum coherency spectrum is produced at lower frequencies indicating that dual-plume mixing is dominated by the external flapping effects of large-scale eddy motions. However, in the far downstream region of the sources, the coherency spectrum in the higher frequency range increases significantly, indicating that the spread of the total plume is larger than all scales of turbulent eddies, such that they all contribute to the in-plume mixing of the dual plumes.     

封闭和半封闭通风空间的环境流体问题

半封闭与全封闭空间的环境流体力学包括系统内的空气流动、温度和湿度分布以及污染物的凝积等问题, 其中涉及控制空间的空气质量、通风设计以及预防水汽凝结等多个方面的研究. 本文着重介绍国内外近年来相关的工作,研究如何从理论分析、实验研究以及数值模拟等方面入手, 解决通风问题,以达到获得高通风效率、低能耗和提高空气舒适度的目的.研究的对象包括房间, 建筑物, 城市等半封闭空间, 以及汽车、 地铁、飞机、潜艇、载人航天飞行器等不同程度的全封闭空间.现有的理论研究已经从解释单点源纯浮力对流的基础模型发展到能够描述具有体积流量和动量流量的扩散流的湍流喷泉模型.然而, 理论研究上的进步还不足以使它能够处理实际问题中多样的通风情况和复杂的空间结构,工程应用中还需更多依靠实验研究和数值模拟的结论. 对载人航天飞行器中的环境流体问题的研究虽然刚刚起步,但其研究手段和经验可以从相对成熟的半封闭与全封闭空间的相关研究中获得借鉴.      

Temperature spectra from a turbulent thermal plume by ultrasound scattering

This paper reports the results of a study on temperature inhomogeneities conducted on a thermal plume by using ultrasound scattering as a non-intrusive measurement technique. The plume rises from a metallic disk which can be heated up to 800 °C. The working fluid is air at atmospheric pressure. In the measurement technique, an incoming ultrasound wave is emitted towards the thermal plume. The incident wave is scattered because of non-linear couplings with the flow instabilities present in the measurement region. The scattered wave carries information about those flow instabilities. The technique allows for the retrieving of this information. The shape of the obtained spectrum of temperature fluctuations as a function of wave vector modulus is consistent with previous theoretical analysis. Three qualitatively different regions were identified: first, a production region characterized by a q² law; secondly, a region with behavior as per q⁻³ associated with a buoyancy region and; finally, a dissipation region associated with a q⁻⁷ law. These spectral regions characterize the energy transfers mechanisms among the length scales of flow investigated here. A coefficient of anisotropy γ was defined to analyze anisotropic features of the flow.     

Nonlinear fluid flow laws for orthotropic porous media

Nonlinear fluid flow laws for orthotropic porous media are written in invariant tensor form. As usual in the theory of fluid flow through porous media [1, 2], the equations contain the flow velocity up to the second power. Expressions that determine the nonlinear resistances to fluid flow are presented and it is shown that, on going over from linear to nonlinear flow laws, the asymmetry effect may manifest itself, that is, the fluid flow characteristics may differ along the same straight line in the positive and negative directions. It is shown that, as compared with the linear fluid flow law for orthotropic media when for three symmetry groups a single flow law is sufficient, in nonlinear laws the anisotropy manifestations are much more variable and each symmetry group must be described by specific equations. A system of laboratory measurements for finding the nonlinear flow characteristics for orthotropic porous media is considered.     

The Effect of Sudden Permeability Changes in Porous Media Filling Box Flows

We report upon experimental and analytical investigations of filling box flows in non-uniform porous media characterized by a sudden change in permeability. The porous medium consists of two layers separated by a horizontal permeability jump and is initially filled with light ambient fluid. A line source located at the top of the upper layer supplies dense contaminated fluid that falls toward the bottom of the domain. Two configurations are studied, i.e., a low-permeability layer on top of a high-permeability layer and vice versa. In the former scenario, the flow dynamics are qualitatively similar to the case of a uniform porous medium. Thus, the analytical formulation of Sahu and Flynn (J Fluid Mech 782:455–478, 2015) can be adopted to compute the parameters of interest, e.g., the plume volume flux. In the latter scenario, the flow dynamics are significantly different from those of the uniform porous medium case; after reaching the permeability jump, some fraction of the dense plume propagates horizontally as a pair of oppositely directed interfacial gravity currents. Meanwhile, the remaining fraction of the plume flows downward into the lower layer where it accumulates along the bottom boundary in the form of a deepening layer of discharged plume fluid. Depending on the permeability ratio of the upper and lower layers and the source conditions, the gravity currents may become temporarily arrested after traveling some finite horizontal length. An analytical prediction for this so-called run-out length is derived, motivated, in part, by the immiscible analysis of Goda and Sato (J Fluid Mech 673:60–79, 2011). Finally, a prediction of the filling box time, consisting of the time required to fill the control volume up to the point of contaminated fluid overflow, is made. These predictions are compared with analog experimental measurements. Generally positive agreement is found when the higher-permeability layer is located below the lower-permeability layer. In the opposite circumstance, the agreement is conditional. If the run-out length of the gravity current is less than the horizontal dimensions of the control volume (or tank in case of the experiments), the agreement is good. By contrast, when the run-out length is large, comparatively poor agreement may be realized: In spite of the higher density of the contaminated fluid, it may occupy the entirety of the upper layer before filling the lower layer.     

Coherent motion of turbulent thermal plume in stably stratified fluid

The purpose of this study is to clarify the existence of an ordered and large scale coherent motion in a turbulent plane thermal plume in a thermally-stable stratified fluid inside a comparatively large enclosure. First, the upper part of the thermal plume was carefully observed by a flow visualization. Secondly, a wave form of plume temperature variation was measured. Thirdly, a spectrum analysis was carried out on time series data of the thermal plume. Finally, physical characteristics were investigated on vortices in the thermal plume based on results of the wave form and the spectrum analysis of the plume temperature. As a result, the main conclusions are obtained as follows. (1) An existence of vortices near the upper part of the thermal plume was firstly found by careful flow visualization. (2) From the wave form of temperature variation and the spectrum analysis of the thermal plume, it was clarified that the vortices are generated in the transition state and are transported to the turbulent state. (3) The vortices are ordered and they behave as a large scale coherent motion in the turbulent thermal plume.     

Modeling transients in the mechanical response of copper due to strain path changes

When copper is deformed to large strains its texture and microstructure change drastically, leading to plastic anisotropy and extended transients when it is reloaded along a different strain path. For predicting these transients, we develop a constitutive model for polycrystalline metals that incorporates texture and grain microstructure. The directional anisotropy in the single crystals is considered to be induced by variable latent hardening associated with cross-slip, cut-through of planar dislocation walls, and dislocation-based reversal mechanisms. These effects are introduced in a crystallographic hardening model which is, in turn, implemented into a polycrystal model. This approach successfully explains the flow response of OFHC Cu pre-loaded in tension (compression) and reloaded in tension (compression), and the response of OFHC Cu severely strained in shear by equal channel angular extrusion and subsequently compressed in each of the three orthogonal directions. This new theoretical framework applies to arbitrary strain path changes, and is fully anisotropic.     

Permeability anisotropy due to consolidation of compressible porous media

The characterisation of flow through porous media is important for all solid–liquid separation and fluid transport realms. The permeability of porous media can be anisotropic and furthermore, the extent of anisotropy can be increased as a result of an applied compressive force. However, the understanding of how anisotropy develops is incomplete. An overview of research on permeability anisotropy is given and an expression for predicting anisotropy as a function of void ratio is offered. The two underlying assumptions of the proposed model are: flow in different directions occurs within the same network of pores and deformation is primarily due to the compression of the particles in the direction of the applied force rather than due to particle rearrangement. The assumption of network connectivity allows permeability anisotropy to be described as a function of flow path tortuosity only. Results are presented for hydraulic anisotropy measured in lignite that has been upgraded by a compression dewatering method known as mechanical thermal expression. The lignite permeability is shown to be up to eight times greater in the direction perpendicular to compression, suggesting that the rate of dewatering could be significantly increased by choosing the drainage to also be perpendicular to the direction of the applied compressive force. It is illustrated that the proposed anisotropy model can be used to accurately predict the experimentally determined permeability anisotropy ratios for lignite, as well as for other materials including sand, clay and kaolin.     

Evolution of a wall-attached buoyant plume in confined boxes: Direct numerical simulations,entrainment coefficient and an integral model

The evolution of a wall-attached plume in a confined box is studied here with the aid of three dimensional direct numerical simulations (DNS). The plume originates from a local line heat source of length, L, placed at the bottom left corner of the box. The Reynolds number of the wall plume, based on box height and buoyant velocity scale, is ${\mathit{Re}}_H=14530$ and boxes of two different aspect ratios (ratio of box width to height) for a particular value of L are simulated. We observe that the plume develops along the vertical sidewall while remaining attached to it before spreading across the top wall to form a buoyant fluid layer and eventually moving downwards and filling the whole box. The original filling box model of Baines and Turner (1969) is modified to incorporate the wall shear stress, and the results from the DNS are compared against the new model. In modelling plumes, we find that the entrainment coefficient ($\alpha$) for wall-attached plumes is reduced to approximately half of that in the free plume, and the main reason is a diminished contribution of turbulence production to $\alpha$ resulting from a restricted ability of the large-scale eddies to transport momentum. Also, unlike the free plume where away from the source inertial forces balances buoyancy forces, here in our simulations of wall-attached plumes this balance is marginally off, likely due to wall friction. A reasonable agreement is observed between our model and DNS data for the volume and momentum fluxes in the quiescent uniform environment and also for the time-dependent buoyancy profile calculated far away from the plume.     

π平面上原状土的空心扭剪试验

为了通过空心扭剪试验获得原状土在π 平面上的屈服特性,在分析空心扭剪试验仪加载过程和加载特点的基础上,研究了空心圆柱试样可实现的应力状态,并设计了相应的加载路径．与重塑试样不同,原状土样的初始应力状态位于K0 固结线上而不是等倾线上,因此,其加载路径应从原点开始沿K0 固结线行至原始应力状态,然后再在π 平面上进行以初始应力为起点的等p 试验．在考虑原状土的初始应力状态和各项异性的基础上,针对具有不同主应力方向角的平行试样,给出了等p 条件下保持主应力方向角不变的加载路径实现方法,采用该方法可得到任意平均主应力时π 平面上六分之一范围内的屈服曲线,该范围对应于大主应力方向角从-45o~45o．本文设计的加载方案可完成真三轴试验的部分功能,但经济成本相对较低,操作方法简单易行．     

The Linear Vortex Instability of the Near-vertical Line Source Plume in Porous Media

Free convection plumes usually rise vertically, but do not do so when in an asymmetrical environment. In such cases they are susceptible to a thermoconvective instability because warmer fluid lies below cooler fluid in the upper half of the plume. We analyse the behaviour of streamwise vortex disturbances in plumes that are close to being vertical. The linearised equations subject to the boundary layer approximation are parabolic and are solved using a marching method. Our computations indicate that disturbances tend to be centred in the upper half of the plume. A neutral curve is determined and an asymptotic theory is developed to describe the right hand branch of this curve. The left hand branch is not amenable to an asymptotic analysis, and it is found that the onset of convection for small wavenumbers is very sensitively dependent on both the profile of the initiating disturbance and where it is introduced.     

Mixing and Entrainment of Transitional Non-Circular Buoyant Reactive Plumes

Three-dimensional spatial direct numerical simulation is used to investigate the evolution of reactive plumes established on non-circular sources. Simulations are performed for three cases: a rectangular plume with an aspect ratio of 2:1, a square plume, and the square plume in a corner configuration. Buoyancy-induced large scale vortical structures evolve spatially in the flow field. A stronger tendency of transition to turbulence is observed for the free rectangular plume than the free square case due to the aspect ratio effect. Dynamics of the corner square plume differs significantly from the corresponding free case due to the enhanced mixing by the side-wall effects. A turbulent inertial subrange has been observed for the free rectangular and corner square plumes. Mean flow properties are also calculated. The study shows significant effects of source geometry and side-wall boundary on the flow transition and entrainment of reactive plumes. This revised version was published online in July 2006 with corrections to the Cover Date.     

A note on mixed convection of laminar plane water plumes at temperatures around the density extremum

In all studies concerning mixed convection in plane laminar plumes a linear relationship between fluid density and temperature has been used. However, it is known that the water density-temperature relationship is non-linear at low temperatures with a density maximum at 3.98°C for pure water. In this note the problem of plane laminar water plume in a coflowing vertical free stream has been investigated taking into account the non-linearity between density and temperature. This is the first work in the literature which treats plane mixed convection plumes with nonlinear relation between density and temperature. Both rising and descending plumes have been investigated. It was found that the ambient water temperature plays an important role on the results. When the ambient temperature is greater than maximum density temperature (T _a > T _m), the water plume behavior is similar to that of the classical plume with linear density-temperature relationship. However, when the ambient temperature is equal or lower than the maximum density temperature the water plume behavior is completely different from the classical plume with linear density-temperature relationship. The centerline velocity shows a series of maxima and minima which are produced by the combination of the nonlinear density-temperature relation and the free stream.     

Dynamics of a spiral convective plume in a high-Prandtl-number fluid

The results of an experimental investigation of a developed convective plume proceeding from a laser-radiation-generated point heat source in a fluid with a high Prandtl number Pr = 2×10³ in the presence of background cellular convective flow are presented. It is found that when the plume growth velocity is similar in value with the characteristic velocity of the cellular convective flow the plume can take the shape of a vertical plane spiral.     

Effects of a transverse flows with variable viscosity in micropolar fluids

A regular perturbation analysis is presented for three laminar natural convection flows in micropolar fluids in liquids with temperature dependent viscosity: a freely-rising plane plume, the flow above a horizontal line source on an adiabatic surface (a plane wall plume) and the flow adjacent to a vertical uniform flux surface. While these flows have well-known power-low similarity solutions when the fluid viscosity is taken to be constant, they are non-similar when the viscosity is considered to a function of temperature. A single similar flow, that adjacent to a vertical isothermal surface, is also analysed for comparison in order to estimate the extent of validity of perturbation analysis. The formulation used here provides a unified treatment of variable viscosity effects on those four flows. Computed first-order perturbation quantities are presented for all four flows. Numerical results for velocity, angular velocity and thermal functions has been shown graphically or tabulated for different values of micropolar parameters. Received on 20 October 1997     

Laboratory observations of interactions of forced plumes with stratified shear layers

Plumes of fluid are often observed in nature to interact with stratified shear layers. Examples of this include chimney plumes hitting inversion-layer ceilings; sewage plumes impinging on unmixed fresh/saltwater interfaces; descending plumes of cold water formed at ice-leads interacting with the oceanic thermocline; and volcano plumes interacting with atmospheric interfaces. Controlled laboratory studies of these phenomena have not previously been described in the literature, and as a result there is a lack of understanding regarding their morphology and dynamics. Thus, a novel set of experiments is described here in which the behaviour of a turbulent plume is observed in the presence of a two-layer ambient. The lower layer, into which the plume initially emerges, is quiescent and at a relatively high density. The upper layer is forced to flow uniformly across the top of the lower layer, and has a lower density. The flow of the resulting plume is characterised by (a) its vertical and lateral spreading in the lower layer; (b) the nature of its extension upstream and downstream at the interface; and (c) the extent to which it penetrates into the upper layer. The behaviour is found to be governed by three non-dimensional parameters: the initial gradient Richardson number of the interface Ri_G, the ratio of the upper layer crossflow speed to the speed of the plume when it first impinges on the interface U_F/U_PI, and the ratio of the plume Monin–Obukhov lengthscale to the lower layer depth L_MO/H_L. Regime diagrams are presented showing the effects of changing these parameters on the plume flow, quantitative relationships are determined, and practical applications of the results are considered.     

Investigation of CO<Subscript>2</Subscript> Plume Behavior for a Large-Scale Pilot Test of Geologic Carbon Storage in a Saline Formation

The hydrodynamic behavior of carbon dioxide (CO₂) injected into a deep saline formation is investigated, focusing on trapping mechanisms that lead to CO₂ plume stabilization. A numerical model of the subsurface at a proposed power plant with CO₂ capture is developed to simulate a planned pilot test, in which 1,000,000 metric tons of CO₂ is injected over a 4-year period, and the subsequent evolution of the CO₂ plume for hundreds of years. Key measures are plume migration distance and the time evolution of the partitioning of CO₂ between dissolved, immobile free-phase, and mobile free-phase forms. Model results indicate that the injected CO₂ plume is effectively immobilized at 25 years. At that time, 38% of the CO₂ is in dissolved form, 59% is immobile free phase, and 3% is mobile free phase. The plume footprint is roughly elliptical, and extends much farther up-dip of the injection well than down-dip. The pressure increase extends far beyond the plume footprint, but the pressure response decreases rapidly with distance from the injection well, and decays rapidly in time once injection ceases. Sensitivity studies that were carried out to investigate the effect of poorly constrained model parameters permeability, permeability anisotropy, and residual CO₂ saturation indicate that small changes in properties can have a large impact on plume evolution, causing significant trade-offs between different trapping mechanisms.