凹坑地形风流结构对污染物散布的模拟研究

利用建立的三维非静力高分辨率高阶湍流闭合模式与随机游动扩散模式研究了一个深凹露天矿区污染物散布的规律,同时在风洞中进行了示踪实验。结果表明,由于凹坑内复环流结构的存在,使得坑内污染物浓度较大,且浓度最大值出现在源的上风侧。数值试验与风洞试验结果吻合较好。     

城市街渠大气湍流结构的小波分析

利用小波分析方法对城市街渠大气湍流场特征进行了初步研究，发现小波变换能够很好地揭示湍流的串级结构和多尺度结构。结合Fourier分析对湍流中不同尺度涡旋对能量的贡献进行了比较，发现噪声能量主要集中在高频处。     

滇西北羊拉地区金沙江古洋盆的穿时碰撞闭合：来自花岗岩年龄的制约

金沙江缝合带是特提斯东段重要的缝合带之一,羊拉地区的金沙江缝合带处在青藏高原东构造结东侧,是研究金沙江古洋盆的碰撞闭合过程等构造事件的重要窗口。在前人工作的基础上,通过野外系统观察取样,结合花岗岩侵入体的锆石SHRIMP U-Pb年龄测定结果,对金沙江古洋盆的碰撞闭合过程做初步研究。结果表明,从南部加仁岩体的通吉格花岗闪长岩(加仁岩体之一部分),经中部的路农—里农—江边岩体的花岗闪长岩,到北部的贝吾花岗闪长岩,岩体的侵位年龄分别为通吉格(246.1±3.5) Ma、路农花岗闪长岩(238.1±5.3) Ma、里农花岗闪长岩(239.0±5.7) Ma、江边岗闪长岩(227.9±5.1) Ma和贝吾花岗闪长岩(213.6±6.9) Ma,表明金沙江缝合带的碰撞关闭是从南向北逐渐闭合的,闭合的时限可达33 Ma,但闭合碰撞不是连续而是幕式的,碰撞闭合作用主要发生在约246 Ma、239 Ma、228 Ma、213 Ma 4个幕次。里农大沟还出现222 Ma的辉绿岩脉,与江边岩体时间上较为接近,似乎表明江边岩体侵入时期本区处于局部拉张的构造背景。     

Floc size variability under strong turbulence: Observations and artificial neural network modeling

The flocculation of cohesive sediment in the presence of waves is investigated using high-resolution field observations and a newly-developed flocculation model based on artificial neural networks. Vertical profiles of suspended sediment concentration and turbulent intensity are estimated using measurements of current profile and acoustic backscatter. The vertical distribution of floc size is estimated using an artificial neural network (ANN) that is trained and validated using floc size measurements at one vertical level. Data analysis suggests a linear correlation between suspended sediment concentration and turbulence intensity. Observations and numerical simulations show that floc size is inversely related to sediment concentration, turbulence intensity and water temperature. The numerical results indicate that floc growth is supported by low concentration and low turbulence. In the vertical direction, mean size of flocs decreases toward the bottom, suggesting floc breakage due to increasing turbulence intensity toward the bed. A significant decrease in turbulent shear could occur within the bottom few-cm, related to increased damping of turbulence by sediment induced density stratification. The results of the numerical simulations presented here are consistent with the concept of a cohesive sediment particle undergoing aggregation-fragmentation processes, and suggest that the ANN can be a precise tool to study flocculation processes.     

Investigation of drag crisis phenomenon using CFD methods

When fluid flow passes a cylinder, the drag crisis phenomenon occurs between the sub-critical and the super-critical Reynolds numbers. The focus of the present studies was on the numerical prediction of the drag crisis based on CFD methods. In this work, block structured meshes with refined grids near the cylinder surface and in the downstream were employed. Both 2D and 3D simulations were performed using various turbulence models, including the SST k − ω model, the k − ϵ model, the SST with LCTM, the DES model, and the LES model. In the convergence studies, the effects of the grid size, the time step, the first grid size and the aspect ratio (for 3D simulations) on the solutions were examined. The errors due to spatial and time discretizations were quantified according to a V&V procedure. Validation studies were carried out for various Reynolds numbers between Re = 6.31 × 10⁴ and 7.57 × 10⁵. The averaged drag force, the RMS of lift force and the Strouhal number were compared with experimental data. The studies indicated that standard 2D and 3D RANS methods were inadequate to capture the drag crisis phenomenon. The LES method however has the potential to address the problem.     

The statistics of a passive scalar in field-guided magnetohydrodynamic turbulence

A variety of studies of magnetised plasma turbulence invoke theories for the advection of a passive scalar by turbulent fluctuations. Examples include modelling the electron density fluctuations in the interstellar medium, understanding the chemical composition of galaxy clusters and the intergalactic medium, and testing the prevailing phenomenological theories of magnetohydrodynamic turbulence. While passive scalar turbulence has been extensively studied in the hydrodynamic case, its counterpart in MHD turbulence is significantly less well understood. Herein we conduct a series of high-resolution direct numerical simulations of incompressible, field-guided, MHD turbulence in order to establish the fundamental properties of passive scalar evolution. We study the scalar anisotropy, establish the scaling relation analogous to Yaglom’s law, and measure the intermittency of the passive scalar statistics. We also assess to what extent the pseudo Alfvén fluctuations in strong MHD turbulence can be modelled as a passive scalar. The results suggest that the dynamics of a passive scalar in MHD turbulence is considerably more complicated than in the hydrodynamic case.     

Fifty-five-year longevity for the largest member of the family Sparidae: the endemic red steenbras Petrus rupestris from South Africa

A previous study that explored the age and growth of red steenbras Petrus rupestris (Valenciennes, 1830), a large sparid (family Sparidae, seabreams or porgies) endemic to South Africa and reported to approach 2 m in length, provided estimates as a moderately slow-growing species with a maximum age near 30 years. The stock is considered collapsed and a fishing moratorium was imposed in 2012, resulting in this species being assessed as Endangered by the IUCN. One consideration in addressing population status is validation of life-history traits. In this study, estimates of age for red steenbras from thin-sectioned otoliths were reassessed visually in terms of both the original ages and revised estimates using a different age-reading pattern. The revised ages exceeded the original ages by up to 30 years, with a maximum estimated age of 55 years from a well-defined otolith section that provided a basis for the revised age-reading procedure. Bomb radiocarbon (¹⁴C) dating revealed there was an offset in the timing of the ¹⁴C rise on the Agulhas Bank that, when coupled with considerations for regional oceanography, provided support for the revised age-reading pattern and an estimated longevity that exceeds 50 years. These findings were further supported by the fortuitous recapture of a tagged red steenbras that was at liberty for 22 years.     

考虑接触面积影响的粗糙节理渗流分析

节理中接触面积的存在使得渗流曲折效应更加明显,将对渗流产生显著影响。为了研究节理接触面积对节理渗流的影响,建立节理渗流概念模型,推导得到了节理渗流的接触面积影响系数;然后在节理渗流计算公式中引入该影响系数,得到了考虑接触面积影响的渗流计算公式。同时,由于接触面积的存在使得渗流曲折因子的计算更为复杂,在获得节理隙宽分布的基础上提出了计算节理渗流曲折因子的五点比较法。然后以人工大理岩节理试件为研究对象,在获得其三维表面形貌的基础上,一方面对其进行不同接触状态下5级法向应力作用下的渗流试验;另一方面,在根据节理试件的三维表面形貌数据计算节理的渗流平均曲折因子和接触面积比的基础上,分别采用考虑接触面积影响的渗流计算公式和Zimmerman计算公式对节理试件在不同法向应力作用下的渗流体积流量进行计算,并将计算结果与实测值进行比较。结果表明：考虑接触面积影响的渗流计算公式计算结果与实测值吻合较好,从而验证了考虑接触面积影响的渗流计算公式的正确性,而Zimmerman计算公式高估了节理的渗流流量。     

荒漠河岸林涡度相关通量测定中平均时间的确定及其对通量计算的影响

极端干旱区在全球生态系统中具有重要的地位而越来越被人们所重视。涡度相关技术是研究极端干旱区生态系统水热和CO2交换的有力工具之一。而涡度相关数据在实际应用中需要根据研究区的实际情况选择适宜于该地区的采样和计算参数（如采样频率和平均时间等）。根据2011年9月22日～10月8日（17个白天）原始数据（采样频率为10 Hz）,采用不同平均时间（1～120 min）对塔河下游柽柳河岸林生态系统的潜热通量[（LE）]、感热通量[（H）]和CO2通量[（Fc）]进行了重新计算,比较分析各通量不同平均时间的计算值与30 min通量值,以揭示平均时间对该地区通量计算结果的影响。结果发现：（1）平均时间在120 min以内,17个白天湍流通量及能量平衡比率[（EBR）]均值随平均时间的延长而增大,当平均时间>120 min后,湍流通量（及[EBR]）均值随平均时间的延长大幅减小,其中15～60 min的[EBR]增幅较小,仅在2%左右。（2）平均时间>60 min后,[LE]、[H]和[FC]各通量计算值发生了不同程度的变异。通过进一步对上午[EBR]较为接近的5个样本日进行对比分析,发现当平均时间取15～60 min之间时,样本日的上午[EBR]值变化趋势一致,样本日间[EBR]差异较小;而当平均时间取>60 min或<15 min时,样本日间[EBR]差异明显增大。结合Ogive函数计算分析,我们的结论是,对于通量的长期观测研究而言,该地区适宜平均时间为60 min;而对通量的日变化研究而言,该地区适宜平均时间为15 min。通过分析15、30和60 min平均时间对小时通量的影响,发现当通量为增加趋势时,平均时间延长能够进一步增大通量绝对值,而当通量为减小趋势时,平均时间延长能够进一步减小通量绝对值。     

Geostatistical modeling of the geological uncertainty in an iron ore deposit

This paper addresses the problem of quantifying the joint uncertainty in the grades of elements of interest (iron, silica, manganese, phosphorus and alumina), loss on ignition, granulometry and rock types in an iron ore deposit. Sampling information is available from a set of exploration drill holes. The methodology considers the construction of multiple rock type outcomes by plurigaussian simulation, then outcomes of the quantitative variables (grades, loss on ignition and granulometry) are constructed by multigaussian joint simulation, accounting for geological domains specific to each quantitative variable as well as for a stoichiometric closure formula linking these variables. The outcomes are validated by checking the reproduction of the data distributions and of the data values at the drill hole locations, and their ability to measure the uncertainty at unsampled locations is assessed by leave-one-out cross validation.Both the plurigaussian and multigaussian models offer much flexibility to the practitioner to face up to the complexity of the variables being modeled, in particular: (1) the contact relationships between rock types, (2) the geological controls exerted by the rock types over the quantitative variables, and (3) the cross-correlations and stoichiometric closure linking the quantitative variables. In addition to this flexibility, the use of efficient simulation algorithms turns out to be essential for a successful application, due to the high number of variables, data and locations targeted for simulation.