Pore scale modeling of immiscible and miscible fluid flows using smoothed particle hydrodynamics

A numerical model based on smoothed particle hydrodynamics (SPH) was developed and used to simulate immiscible and miscible fluid flows in porous media and to study effects of pore scale heterogeneity and anisotropy on such flows.     

A bed‐load transport model for rough turbulent open‐channel flows on plane beds

Data from flume studies are used to develop a model for predicting bed‐load transport rates in rough turbulent two‐dimensional open‐channel flows moving well sorted non‐cohesive sediments over plane mobile beds. The object is not to predict transport rates in natural channel flows but rather to provide a standard against which measured bed‐load transport rates influenced by factors such as bed forms, bed armouring, or limited sediment availability may be compared in order to assess the impact of these factors on bed‐load transport rates. The model is based on a revised version of Bagnold's basic energy equation i_bs_b = e_bω, where i_b is the immersed bed‐load transport rate, ω is flow power per unit area, e_b is the efficiency coefficient, and s_b is the stress coefficient defined as the ratio of the tangential bed shear stress caused by grain collisions and fluid drag to the immersed weight of the bed load. Expressions are developed for s_b and e_b in terms of G, a normalized measure of sediment transport stage, and these expressions are substituted into the revised energy equation to obtain the bed‐load transport equation i_b = ω G ^3·4. This equation applies regardless of the mode of bed‐load transport (i.e. saltation or sheet flow) and reduces to i_b = ω where G approaches 1 in the sheet‐flow regime. That i_b = ω does not mean that all the available power is dissipated in transporting the bed load. Rather, it reflects the fact that i_b is a transport rate that must be multiplied by s_b to become a work rate before it can be compared with ω. It follows that the proportion of ω that is dissipated in the transport of bed load is i_bs_b/ω, which is approximately 0·6 when i_b = ω. It is suggested that this remarkably high transport efficiency is achieved in sheet flow (1) because the ratio of grain‐to‐grain to grain‐to‐bed collisions increases with bed shear stress, and (2) because on average much more momentum is lost in a grain‐to‐bed collision than in a grain‐to‐grain one. Copyright © 2006 John Wiley & Sons, Ltd.     

裂隙岩质边坡地下水渗流数值模拟

本文从反倾向岩质高边坡的渗流特征和渗流普遍规律出发,根据渗流力学知识,建立数学模型,通过渗流数值模拟计算分析了渗透压力、水头变化规律以及对岩质边坡稳定性的影响.     

鄂尔多斯风沙滩地区土壤凝结水试验研究

鄂尔多斯风沙滩地区,常年干旱少雨,蒸发强烈,地下水资源消耗主要源于表土蒸发.而土壤凝结水的形成与转化又与表土蒸发之间存在着十分密切的关系.作者在内蒙古自治区乌审旗气象局院内,采用自制的笼屉式测渗仪进行土壤凝结水的观测试验研究,测量方法为称重法.并同步观测地表以上50cm处的空气温湿度以及地表和地表以下5cm、10cm、15cm、20cm、30cm、40cm、50cm处的土壤温度.通过试验,初步探明了该区土壤凝结水的形成与转化的基本过程,以及土壤凝结水的转化与气象因素之间的相关关系.     

基于GIS的滑坡、泥石流灾害危险性区划关键问题研究

随着GIS技术的引入,滑坡、泥石流灾害危险性区划的效率和准确性得以大大提高。依据工程地质类比原则,在灾害学理论指导下,结合专家打分、层次分析、人工神经网络、信息量、Logistic回归、统计量等模型方法,以MAPGIS软件为平台,利用C++语言开发了滑坡、泥石流灾害危险性区划评价分析系统;并重点探讨了GIS支持下的滑坡、泥石流灾害危险性区划过程中的因子分析、模型选取、模型复合、单元划分、系统集成、结果评价等关键问题,建立了一整套基于GIS的滑坡、泥石流灾害评价方法体系。应用该系统对长江三峡库区和辽宁省鞍山市分别开展了滑坡、泥石流灾害危险性区划研究,取得了较好的效果。     

Mantle heat flow and thermal structure of the northern block of Southern Granulite Terrain, India

Continental shield regions are normally characterized by low-to-moderate mantle heat flow. Archaean Dharwar craton of the Indian continental shield also follows the similar global pattern. However, some recent studies have inferred significantly higher mantle heat flow for the Proterozoic northern block of Southern Granulite Terrain (SGT) in the immediate vicinity of the Dharwar craton by assuming that the radiogenic elements depleted exposed granulites constitute the 45-km-thick crust. In this study, we use four-layered model of the crustal structure revealed by integrated geophysical studies along a geo-transect in this region to estimate the mantle heat flow. The results indicate that: (i) the mantle heat flow of the northern block of SGT is 17 ± 2 mW/m², supporting the global pattern, and (ii) the lateral variability of 10–12 mW/m² in the surface heat flow within the block is of crustal origin. In terms of temperature, the Moho beneath the eastern Salem–Namakkal region appears to be at 80–100 °C higher temperature than that beneath the western Avinashi region.     

Velocity distribution in a gradually accelerating free surface flow

This study investigates the interaction of the vertical velocity v and the streamwise velocity u in a gradually accelerating flow. The analytical result shows that the momentum of uv driven by the mean velocities in a non-uniform flow is not negligible. This additional momentum directly results in the concave profiles of Reynolds shear stress in gradually accelerating flows, a departure from the expected linear profile. Consequently, this momentum causes the maximum velocity to be located below the free surface, i.e., the dip-phenomenon. This paper investigated the interactions of the Reynolds shear stress, non-zero vertical velocity and dip-phenomenon, it is found that the non-zero vertical velocity causes the dip-phenomenon. The approach is tested using the experimental data of Song and others, and good agreements between the predicted and measured velocity profiles have been achieved.     

Mantle process beneath Philippine Sea back-arc spreading ridges: A synthesis of peridotite petrology and tectonics

Abstract In order to obtain a general view of the mantle process beneath a back‐arc basin spreading ridge, the diversity of peridotite petrology and tectonic occurrences in two back‐arc basin spreading ridges from the Philippine Sea were examined: the Parece Vela Rift and the Mariana Trough. The Parece Vela Basin spreading ridge (Parece Vela Rift) was a physically fast/intermediate‐spreading ridge, although many tectono‐magmatic features resemble those of slow‐ to ultraslow‐spreading ridges. Two unusual features of the Parece Vela Rift further demonstrate the uniqueness of the ridge: full‐axial development of oceanic core complexes and exposure of mantle peridotite at segment midpoints. The Parece Vela Rift yields a lithological assemblage of residual but still fertile lherzolite/harzburgite, plagioclase‐bearing harzburgite and dunite; similar assemblages are reported from the equatorial Mid‐Atlantic Ridge at the Romanche Fracture Zone and the ultraslow‐spreading ridges from the Indian and Arctic Oceans. The tectono‐magmatic characteristics of the Parece Vela Rift suggest that diffuse porous melt flow and pervasive melt–mantle interaction were the important mantle processes there. Globally, this ‘porous melt flow‐type’ mantle process is likely to occur beneath a segment midpoint of the ridge having a thick lithosphere, typically an ultraslow‐spreading ridge. In contrast, the Mariana Trough is a typical slow‐spreading ridge, exposing mantle peridotite at segment ends. The Mariana Trough yields a lithological assemblage of residual harzburgite and veined harzburgite, a common assemblage among the global abyssal peridotite suite. The tectono‐magmatic characteristics of the Mariana Trough suggest that channeled melt/fluid flow and limited melt–mantle interaction are the important mantle processes there, because of the colder wall‐rock peridotite in the segment end. This ‘channeled melt flow‐type’ mantle process is likely to occur in the shallow lithospheric mantle at the segment ends of any spreading ridges.     

一种海面风场反演的快速风矢量搜索算法

本文以SeaWinds散射计为例，从最大似然估计的海面风场反演方法出发，根据其目标函数的一般分布特征，提出并实现了一种快速风矢量搜索算法。然后采用美国SeaWinds散射计的部分LA实测数据对该搜索算法进行了验证，并与常规算法对同一实验数据的反演偏差及运行时间进行了对比分析，结果表明该算法在能够保证反演精度的前提下，在较大程度上提高了运行效率。