Stochastic modelling of groundwater flow and solute transport in aquifers

This paper presents an introductory overview of recently developed stochastic theories for tackling spatial variability problems in predicting groundwater flow and solute transport. Advantages and limitations of the theories are discussed. Lastly, strategies based on the stochastic approaches to predict solute transport in aquifers are recommended.     

Three dimensional simulation of fluid flow in X‐ray CT images of porous media

A numerical scheme is developed in order to simulate fluid flow in three dimensional (3‐D) microstructures. The governing equations for steady incompressible flow are solved using the semi‐implicit method for pressure‐linked equations (SIMPLE) finite difference scheme within a non‐staggered grid system that represents the 3‐D microstructure. This system allows solving the governing equations using only one computational cell. The numerical scheme is verified through simulating fluid flow in idealized 3‐D microstructures with known closed form solutions for permeability. The numerical factors affecting the solution in terms of convergence and accuracy are also discussed. These factors include the resolution of the analysed microstructure and the truncation criterion. Fluid flow in 2‐D X‐ray computed tomography (CT) images of real porous media microstructure is also simulated using this numerical model. These real microstructures include field cores of asphalt mixes, laboratory linear kneading compactor (LKC) specimens, and laboratory Superpave gyratory compactor (SGC) specimens. The numerical results for the permeability of the real microstructures are compared with the results from closed form solutions. Copyright © 2004 John Wiley & Sons, Ltd.     

Reliability assessment of granular filters in embankment dams

Empirical criteria have been used successfully to design filters of most embankment large dam projects throughout the world. However, these empirical rules are only applicable to a particular range of soils tested in laboratory and do not take into account the variability of the base material and filter particle sizes. In addition, it is widely accepted that the safety of fill dams is mainly dependent on the reliability of their filter performance. The work herein presented consists in a new general method for assessing the probability of fulfilling any empirical filter design criteria accounting for base and filter heterogeneity by means of first‐order reliability methods (FORM), so that reliability indexes and probabilities of fulfilling any particular criteria are obtained. This method will allow engineers to estimate the safety of existing filters in terms of probability of fulfilling their design criteria and might also be used as a decision tool on sampling needs and material size tolerances during construction. In addition, sensitivity analysis makes possible to analyse how reliabilities are influenced by different sources of input data. Finally, in case of a portfolio risk assessment, this method will allow engineers to compare the safety of several existing dams in order to prioritize safety investments and it is expected to be a very useful tool to evaluate probabilities of failure due to internal erosion. Copyright © 2006 John Wiley & Sons, Ltd.     

Fluid infiltration during contact metamorphism of interbedded marble and calc-silicate hornfels, Twin Lakes area, central Sierra Nevada, California

Abstract In the Twin Lakes area, central Sierra Nevada, California, most contact metamorphosed marbles contain calcite + dolomite + forsterite ± diopside ± phlogopite ± tremolite, and most calc-silicate hornfelses contain calcite + diopside + wollastonite + quartz ± anorthite ± K-feldspar ± grossular ± titanite. Mineral-fluid equilibria involving calcite + dolomite + tremolite + diopside + forsterite in two marble samples and wollastonite + anorthite + quartz + grossular in three hornfels samples record P± 3 kbar and T± 630° C. Various isobaric univariant assemblages record CO₂-H₂O fluid compositions of χCO₂= 0.61–0.74 in the marbles and χCO₂= 0.11 in the hornfelses. Assuming a siliceous dolomitic limestone protolith consisting of dolomite + quartz ° Calcite ± K-feldspar ± muscovite ± rutile, all plausible prograde reaction pathways were deduced for marble and hornfels on isobaric T-XCO₂ diagrams in the model system K₂O-CaO-MgO-Al₂O₃-SiO₂-H₂O-CO₂. Progress of the prograde reactions was estimated from measured modes and mass-balance calculations. Time-integrated fluxes of reactive fluid which infiltrated samples were computed for a temperature gradient of 150 °C/km along the fluid flow path, calculated fluid compositions, and estimated reaction progress using the mass-continuity equation. Marbles and hornfelses record values in the range 0.1–3.6 × 10⁴ cm³/cm² and 4.8–12.9 × 10⁴ cm³/cm², respectively. For an estimated duration of metamorphism of 10⁵ years, average in situ metamorphic rock permeabilities, calculated from Darcy's Law, are 0.1–8 × 10^?6 D in the marbles and 10–27 × 10^?6 D in the hornfelses. Reactive metamorphic fluids flowed up-temperature, and were preferentially channellized in hornfelses relative to the marbles. These results appear to give a general characterization of hydrothermal activity during contact metamorphism of small pendants and screens (dimensions ± 1 km or less) associated with emplacement of the Sierra Nevada batholith.     

Hysteresis‐based analysis of overland metal transport

Introducing a concept of equivalent mass depth of flow, this study describes the phenomenon of non‐point source pollutant (metal) transport for pavement (or overland) flow in analogy with wave propagation in wide open channels. Hysteretic and normal mass rating curves are developed for runoff rate and mass of 12 dissolved and particulate‐bound metal elements (pollutants) using the rainfall‐runoff and water quality data of the 15 × 20 m² instrumented pavement in Cincinnati, USA. Normal mass rating curves developed for easy computation of pollutant load are found to be of a form similar to Manning's, which is valid for open channel flows. Based on the hysteresis analysis, wave types for dissolution and mixing of particulate‐bound metals are identified. The analysis finds that the second‐order partial‐differential equation normally used for metal transport does not have the efficacy to describe fully the strong non‐linear phenomena such as is described for various metal elements by dynamic waves. In addition, the proportionality concept of the popular SCS‐CN concept is extended for determining the potential maximum metal mass M_p of all the 12 elements transported by a rain storm and related to the antecedent dry period (ADP). For the primary metal zinc element, M_p is found to increase with the ADP and vice versa. Copyright © 2003 John Wiley & Sons, Ltd.     

全球海洋与大气界面的熵通量估算

给出并讨论了海洋系统中由于热量交换、物质交换以及外力做功引起的熵通量的表达式;同时,根据新近调整的全球海 气界面的气候平均资料,估计了全球海洋与大气界面气候年平均的熵通量(熵流),这一熵通量主要取决于海 气界面的热量和物质交换的空间不均匀性。计算结果显示,由热量交换分量引起的熵通量对海洋系统的总熵通量起到主要的贡献作用,其值约为-555.6mW·(m2·K)-1;由海表风应力做功引起的熵通量相对较小,约为-0.09mW·(m2·K)-1;由物质交换引起的熵通量最小,仅为-0.02mW·(m2·K)-1。总的来说,海洋系统从外界获得的总熵通量为-555.7mW·(m2·K)-1,这也就意味着在气候平均定态下,海洋系统内部的熵产生在量值上等于系统的熵交换,即为555.7mW·(m2·K)-1。海洋系统的负熵流与其内部的各种不可逆过程引起的熵产生取得平衡,确保了全球海洋系统处于非平衡热力学定态,并维持着海洋系统中各种尺度的时 空有序现象的消长过程。     

Freshening of the Intermediate Waters in the Northern South China Sea over the PastSix Decades

The properties of salinity in the South China Sea (SCS), a significant marginal sea connecting the Pacific andIndian Oceans, are greatly influenced by the transport of fresh water flux between the two oceans. However, the long-termchanges in the intermediate water in the SCS have not been thoroughly studied due to limited data, particularly in relationto its thermodynamic variations. This study utilized reanalysis data products to identify a 60-year trend of freshening in theintermediate waters of the northern South China Sea (NSCS), accompanied by an expansion of low-salinity water. Thestudy also constructed salinity budget terms, including advection and entrainment processes, and conducted an analysis ofthe salinity budget to understand the impacts of external and internal dynamic processes on the freshening trend of theintermediate water in the NSCS. The analysis revealed that the freshening in the northwest Pacific Ocean and theintensification of intrusion through the Luzon Strait at intermediate levels are the primary drivers of the salinity changes inthe NSCS. Additionally, a weakened trend in the intensity of vertical entrainment also contributes to the freshening in theNSCS. This study offers new insights into the understanding of regional deep sea changes in response to variations in boththermodynamics and oceanic dynamic processes.     

Transport and Housing Security in the Klang Valley, Malaysia

This paper explores links between transport and housing security issues for the urban poor using the example of the Klang Valley in Malaysia. The interface between these issues is identified as a gap in the literature, including policy debates, on both housing and transport. A number of linkages are shown to be important and likely to be relevant in many cities of the South, especially those with rapid motorisation and large numbers of "squatters". A simple framework for understanding these linkages is presented. Key examples include displacement to make way for transport infrastructure and the impact on transport problems for the poor of policies affecting the location of urban poor housing, including relocation sites and transit accommodation. The case study of the Klang Valley is used to illustrate and test the relevance of a focus on this issue and the utility of the conceptual framework. Some policy implications of the investigation and case study are suggested.     

图形系统中微机制图整饰的应用研究

本文着重阐述了微机图形系统中图面动态整饰的应用研究，其中包括对图面内容、图名、图例，以及内、外图廓的整饰。图形系统的图面动态整饰，主要是针对计算机图形学的研究对象、其整饰的内容以及现在多数微机图形系统普遍存在的问题，采用不同的数学算法及模型，研究图面整饰的技巧和方法。图面内容的动态整饰是全文的重点，它主要包含：多要素地图的依比例分级动态显示、汉字注记的变体、注记的配置整饰和注记的沿线性要素动态分布显示等。因此，图形系统中微机制图整饰对其传统工艺进行革新，使人们从繁锁复杂的工艺劳动中解脱出来，极大地提高了工作效率，同时改进了图形系统的制图工艺。     

基于复模态的有限元模型修正算法

针对地下结构地震响应分析中无限地基辐射阻尼问题，引入复模态情况下的具有非简化的堆积阻尼矩阵的阻尼模型，并针对具有集中质量阵的阻尼模型提出了合并与质量有关的阻尼和堆积阻尼的思想，并据此提出了一种修正此类有限元模型的两步法，首先从复模态参数中提取实模态参数，采用基于模态残余力的识别算法修正刚度矩阵，然后根据复模态参数和已得的刚度矩阵来识别阻尼模型中的刚度参与系数和质量阻尼堆积阻尼联合矩阵。