基于高斯烟羽模型的船舶尾气扩散研究

为有效对港区大气污染进行治理、分析船舶尾气,本文详细介绍了一种基于高斯烟羽模型,通过MATLAB模拟仿真模型,其包括实验仿真过程、技术原理及理论模型对船舶尾气扩散进行的研究。该模型是在传统的高斯烟羽模型的基础上,通过对实源像源进行加权选择输入参数;通过矢量合成确定了气体扩散的方向,利用合成后的"风速"进行计算仿真,有效模拟了船舶尾气在港区或者海洋环境中的气体扩散模型。其模型简单且可以有效模拟船舶尾气扩散。并且进一步对后续模型的精确优化进行分析。     

Ambipolar diffusion in the solar atmosphere

The process of non-linear ambipolar diffusion in the region overlying the solar surface can be an effective mechanism for producing sharp magnetic structures and current sheets. These may be the sites responsible for the occurrence of connectivity of magnetic field lines, and the subsequent explosive input of energy for heating of some of the features in the atmosphere of the Sun..     

MONITORING THE DIFFUSION OF SUSPENDED SEDIMENTS AND STABILITY OF TIDAL RADIAL SAND RIDGES AREA USING MULTI-SATELLITE REMOTE SENSING DATA

Nine Landsat TM tapes and images and MSS images, 10 NOAA tapes and images and 1 SAR image from 1973 to 1997 were used to analyse the diffusion of suspended sediments and the change of tidal radial sand ridges in the northern part of the Changjiang River delta, the South Yellow Sea. The results showed that the diffusion of suspended sediments was controlled by the tide, net current, and submarine topography in this area. The distribution of suspended sediments had close relationship with thesubmarine topography. The old Huanghe River delta and the Changjiang River comprise the main sediment supply for the formation of radial sand ridges, whose evolution can be divided into three stages since the Huangbe River changed its course and flowed northward into the Bohai Sea.     

Sediment transport in analogue flume models compared with real‐world sedimentary systems: a new look at scaling evolution of sedimentary systems in a flume

Evolution of sedimentary systems at large temporal and spatial scales cannot be scaled down to laboratory dimensions by conventional hydraulic Froude scaling. Therefore, many researchers question the validity of experiments aiming to simulate this evolution. Yet, it has been shown that laboratory experiments yield stratigraphic responses to allocyclic forcing that are remarkably similar to those in real‐world prototypes, hinting at scale independency with strong dependence on boundary conditions but weak dependence on the actual sediment transport dynamics. This paper addresses the dilemma by contrasting sediment transport rules that apply in the laboratory with those that apply in real‐world geological systems. It is demonstrated that the generation of two‐dimensional stratigraphy in a flume can be simulated numerically by the non‐linear diffusion equation. Sediment transport theory is used to demonstrate that only suspension‐dominated meandering rivers should be simulated with linear diffusion. With increasing grain‐size (coarse sand to gravel) and shallowness of river systems, the prediction of long‐term transport must be simulated by non‐linear, slope‐dependent diffusion to allow for increasing transport rates and thus change in stratigraphic style. To point out these differences in stratigraphic style, three stages in infill of accommodation have been defined here: (i) a start‐up stage, when the system is prograding to base level (e.g. the shelf edge) with no sediment flux beyond the base‐level point; (ii) a fill‐up stage, when the system is further aggrading while progressively more sediment is bypassing base level with the progression of the infill; and (iii) a keep‐up stage, when more than 90% of the input is bypassing the base level and less than 10% is used for filling the accommodation. By plotting the rate of change in flux for various degrees of non‐linearity (varying the exponent in the diffusion equation) it was found that the error between model and real‐world prototype is largest for the suspension‐dominated prototypes, although never more than 30% and only at the beginning of the fill‐up stage. The error reduces to only 10% for the non‐linear sandy‐gravelly and gravelly systems. These results are very encouraging and open up ways to calibrate numerical models of sedimentary system evolution by such experiments.     

Asymmetrically zoned reaction rims: assessment of grain boundary diffusivities and growth rates related to natural diffusion-controlled mineral reactions

This study explores garnet coronas around hedenbergite, which were formed by the reaction plagioclase + hedenbergite→garnet + quartz, to derive information about diffusion paths that allowed for material redistribution during reaction progress. Whereas quartz forms disconnected single grains along the garnet/hedenbergite boundaries, garnet forms ～20‐μm‐wide continuous polycrystalline rims along former plagioclase/hedenbergite phase boundaries. Individual garnet crystals are separated by low‐angle grain boundaries, which commonly form a direct link between the reaction interfaces of the plagioclase|garnet|hedenbergite succession. Compositional variations in garnet involve: (i) an overall asymmetric compositional zoning in Ca, Fe²⁺, Fe³⁺ and Al across the garnet layer; and (ii) micron‐scale compositional variations in the near‐grain boundary regions and along plagioclase/garnet phase boundaries. These compositional variations formed during garnet rim growth. Thereby, transfer of the chemical components occurred by a combination of fast‐path diffusion along grain boundaries within the garnet rim, slow diffusion through the interior of the garnet grains, and by fast diffusion along the garnet/plagioclase and the garnet/hedenbergite phase boundaries. Numerical simulation indicates that diffusion of Ca, Al and Fe²⁺ occurred about three to four, four and six to seven orders of magnitude faster along the grain boundaries than through the interior of the garnet grains. Fast‐path diffusion along grain boundaries contributed substantially to the bulk material transfer across the growing garnet rim. Despite the contribution of fast‐path diffusion, bulk diffusion through the garnet rim was too slow to allow for chemical equilibration of the phases involved in garnet rim formation even on a micrometre scale. Based on published garnet volume diffusion data the growth interval of a 20‐μm‐wide garnet rim is estimated at ～10³–10⁴ years at the inferred reaction conditions of 760 ± 50 °C at 7.6 kbar. Using the same parameterization of the growth law, 100‐μm‐ and 1‐mm‐thick garnet rims would grow within 10⁵–10⁶ and 10⁶–10⁷ years respectively.     

THERIA_G: a software program to numerically model prograde garnet growth

We present the software program THERIA_G, which allows for numerical simulation of garnet growth in a given volume of rock along any pressure–temperature–time (P–T–t) path. THERIA_G assumes thermodynamic equilibrium between the garnet rim and the rock matrix during growth and accounts for component fractionation associated with garnet formation as well as for intracrystalline diffusion within garnet. In addition, THERIA_G keeps track of changes in the equilibrium phase relations, which occur during garnet growth along the specified P–T–t trajectory. This is accomplished by the combination of two major modules: a Gibbs free energy minimization routine is used to calculate equilibrium phase relations including the volume and composition of successive garnet growth increments as P and T and the effective bulk rock composition change. With the second module intragranular multi-component diffusion is modelled for spherical garnet geometry. THERIA_G allows to simulate the formation of an entire garnet population, the nucleation and growth history of which is specified via the garnet crystal size frequency distribution. Garnet growth simulations with THERIA_G produce compositional profiles for the garnet porphyroblasts of each size class of a population and full information on equilibrium phase assemblages for any point along the specified P–T–t trajectory. The results of garnet growth simulation can be used to infer the P–T–t path of metamorphism from the chemical zoning of garnet porphyroblasts. With a hypothetical example of garnet growth in a pelitic rock we demonstrate that it is essential for the interpretation of the chemical zoning of garnet to account for the combined effects of the thermodynamic conditions of garnet growth, the nucleation history and intracrystalline diffusion. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.

Possible Density Change of the Origin of North Pacific Intermediate Water Due to Double Diffusion in the Mixed Water Region

In this study we test Talley's hypothesis that Oyashio winter mixed-layer water (26.5–26.6σ _θ) increases its density to produce the North Pacific Intermediate Water (NPIW) salinity minimum (26.7– 26.8σ_θ) in the Mixed Water Region, assuming a combination of cabbeling and double diffusion. The possible density change of Oyashio winter mixed-layer water is discussed using an instantaneous ratio of the change of temperature and salinity along any particular intrusion (R _l ). We estimate the range of R _l ^DD required to convert Oyashio winter mixed-layer water to the NPIW salinity minimum due to double diffusion, and then assume double-diffusive intrusions as this conversion mechanism. A double-diffusive intrusion model is used to estimate R _l ^DD in a situation where salt fingering dominates vertical mixing, as well as to determine whether Oyashio winter mixed-layer water can become the NPIW salinity minimum. Possible density changes are estimated from the model R _l ^DD by assuming the amount of density change due to cabbeling. From these results, we conclude that Oyashio winter mixed-layer water contributes to a freshening of the lighter layer of the NPIW salinity minimum (around 26.70σ_θ) in the MWR.     

Three-dimensional Layer-integrated Modelling of Estuarine Flows with Flooding and Drying

Details are given of the refinement and application of a thee-dimensional (3-D) layer-integrated numerical model of tidal circulation, with the aim of simulating severe tidal conditions for practical engineering applications. The mode splitting strategy has been used in the model. A set of depth-integrated 2-D equations are first solved to give the pressure gradient, and the layer-integrated 3-D equations are then solved to obtain the vertical distributions of the flow velocities. Attention has been given to maintaining consistency of the physical quantities derived from the 2-D and 3-D equations. A TWO=layer mixing length turbulence model for the vertical shear stress distribution has been included in the model. Emphasis has been focused on applying the model to a real estuary, which is geometrically complicated and has large tidal ranges giving rise to extensive flooding and drying. The model has been applied to three examples, including: wind-driven flow in a rectangular lake, tidal circulation in a model rectangular harbour, and tidal circulation in a large estuary. Favourable results have been obtained for both the simple and complex flow beds.     

西江河口重金属迁移扩散数值模拟的研究

本文以水体中重金属(包括溶解态和颗粒态)的物理迁移过程为基础,建立重金属迁移扩散的数学模式。模式中考虑悬浮物沉降和再悬浮对水体中重金属浓度的影响,并考虑咸潮活动引起悬浮物的絮凝沉淀作用。计算结果表明,这种潮汐河口区重金属迁移扩散模式基本上是成功的。并以镉为例,说明重金属在西江河口区迁移扩散的特点。     

近岸溢油漂移扩散预测方法研究——以胶州湾溢油事件为例

近些年,由于中国经济的高速发展,原油需求逐年增长,同时海上溢油事故发生的风险也在加剧。海上溢油会严重破坏海洋环境,危害我国经济发展,加强溢油漂移扩散预报研究能够为海上溢油应急响应提供技术支撑。因此,国家海洋信息中心研发了海上突发事件应急一体化预测预警系统。该系统基于GIS平台研发,能够在二维电子海图基础上叠加相关的海洋环境动力要素信息,模拟溢油扩散和漂移态势,同时计算剩余油量,估算溢油面积以及岸线吸附程度等。2013年11月22日,山东黄岛发生输油管道爆炸事故,造成大量原油溢入胶州湾。本文利用该系统结合高分辨率大气模型WRF的三重嵌套技术和海流模型SELFE的非结构化网格加密技术,对事故溢油漂移扩散开展了预测工作。预测结果显示,油污受胶州湾内往复流的影响极大,72 h后,胶州湾内外大部分海域都将存在油污,主要集中在黄岛沿岸海域、团岛及浮山湾附近,预测扫海面积为70.32 km2。经海事部门高分辨率卫星观测结果印证,系统模拟的溢油分布较为合理。该系统的业务化应用为相关海上溢油污染应急工作提供了一定理论依据和参考。