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.     

PROPAGATION OF 30—60 DAY LOW FREQUENCY OSCILLATIONS AND THEIR INFLUENCE UPON THE SUBTROPICAL WESTERLIES JET STREAM DURING NORTHERN HEMISPHERE WINTER

Based on daily ECMWF gridpoint data of two winters during 1981—1983 including an ENSOyear,propagation of low frequency oscillations(LFO)during Northern Hemisphere winters andtheir influences upon 30—60 day oscillations of the subtropical jet stream are studied with the sta-tistical methods as complex empirical orthogonal function(CEOF)and so on.Results show that inthe winter of a normal year(1981—1982),30—60 day oscillations in the subtropical zone aremainly in the northern and southern flanks of exit region of jet stream.In the ENSO year(1982—1983),they are mainly in the vicinity of entrance and exit regions of jet stream.Intraseasonalchanges of subtropical jet stream manifested themselves as latitudinal fluctuation or longitudinalprogression or regression of about 40 day period.There are marked differences between propagat-ing passages of low frequency modes responsible for changes of subtropical jet stream in the normalyear(1981—1982)and in the ENSO year(1982—1983).Changes of oscillation amplitude showobvious phases.In general,the one in late winter is stronger than that in early winter,strongestone occurs in February.     

Evolution of 2:1 layered silicates in low-grade metamorphosed Liassic shales of Central Switzerland

White mica from the Liassic black shales and slates in Central Switzerland was analysed by transmission electron microscopy (TEM) and electron microprobe to determine its textural and compositional evolution during very low-grade prograde metamorphism. Samples were studied from the diagenetic zone, anchizone and epizone (T ≈100°–450 °C). Phyllosilicate minerals analysed include illite/smectite (I/S), phengite, muscovite, brammallite, paragonite, margarite and glauconite. Textural evolution primarily is towards larger, more defect-free grains with compositions that approach those of their respective end-members. The smectite-to-illite transformation reduced the amounts of the exchange components SiK_?1Al_?1, MgSiAl_?2, and Fe³⁺Al_?1. These trends continue to a lesser degree in the anchizone and epizone. Correlations between the proportion of smectite in I/S and the composition of I/S indicate that smectite layers may contain a high layer charge. Illite in I/S bears a compositional resemblance to macrocrystalline phengite in some samples, but is different in others. Paragonite first appears in the upper diagenetic zone or lower anchizone as an interlayer-deficient brammallite, and it may be mixed with muscovite on the nanometre scale. Owing to the small calculated structure factor for paragonite-muscovite superstructures, conventional X-ray powder diffraction cannot distinguish between mixed-layer structures and a homogeneous compositionally intermediate solid solutions. However, indirect TEM evidence shows that irregularly shaped domains of Na- and K-rich mica exist below 10 nm. Subsequent coarsening of domains at higher grades produced discrete paragonite grains at the margins of muscovite crystals or in laths parallel to the basal plane of the host muscovite. Margarite appears in the epizone and follows a textural evolution similar to paragonite in that mixtures of margarite, paragonite, and muscovite may initially occur on the nanometre scale. However, no evidence of interlayer-poor margarite has been found.     

Optimal distribution of linear control intensity over the frequency range

In this paper, a control strategy for structural systems is proposed and developed in the frequency domain. The algorithm is substantially based on a linear derivative feedback and a convolution of the control parameter, whose distribution in the frequency field is chosen in such a manner as to comply with the requirements of an ad hoc formulated constrained optimum problem, with the response data monitored until the instant of control action application. Some numerical testing is carried out by referring to given recorded accelerograms, showing a good performance of the proposed approach. Copyright © 2002 John Wiley & Sons, Ltd.     

塑膜温室油桃休眠期的温度变化

本文详尽分析了冰岛低压和大西洋高压的长期变化规律，发现两者均具有明显的阶段性变化，而且冰岛低压和大西洋高压的强度具有明显的反相关。同时，还分析了大西洋高压和冰岛低压与我国气候的关系，发现大西洋高压与我国冬季气温基本呈正相关，而冰岛低压与我国气温大致呈反相关，但均没有显著影响，所以两者都不是影响我国气候变化的重要天气系统。     

环形构造及其地震预报意义

根据遥感图象的解译结果，介绍了环形构造的基本特征及其对中强地震的控震作用，指出了环形构造的地震构造模式和环形构造对分析地震构造环境和地震预报具有一定的意义。     

一次辽宁秋季暴雨天气的诊断分析

使用1.0°×1.0°NCEP再分析资料,对2006年10月21—22日深秋暴雨在天气形势分析的基础上,进行物理量诊断。结果表明:在有利的环境背景形势下,高位涡从对流层高层向低层伸展并形成湿位涡柱,引起气旋性环流与低涡环流叠加。对流层低层的湿斜压性增强,引起低层的锋区加强及垂直涡度发展,高空入侵干冷空气锲入底层,低层暖湿空气强迫抬升,使地面发展为气旋;高低空急流耦合产生上升气流,同时较强的补偿下沉运动激发上升运动加强,使次级环流加强,触发不稳定能量的释放;低空急流和超低空急流向辽宁输送暖湿空气及能量,对流层中低层形成湿柱并积聚高不稳定能量;中尺度气旋、高低空急流、湿位涡柱、次级环流上升支、地面高水汽含量湿区、高假相当位温出现的时间、强度、位置和结构决定了暴雨的时间和落区。     

下扬子地区若干上地幔橄榄岩的显微构造研究

本文对下扬子区内２个玄武岩产地（江苏六合瓜埠山、安徽当涂釜山）的幔源包体进行研究,结果表明本区上地幔尖晶石二辉橄榄岩具原生粒状结构、碎斑结构、板状镶嵌等粒结构及部分熔融结构四种类型。在光学岩组图上显示出从原生粒状结构→碎斑结构→板状镶嵌等粒结构,其橄榄石光性方位优选性逐渐增强。变形橄榄石滑移轴为［１００］,滑移面以｛０ｋ１｝为主,（０１０）次之。高温高压人工变形实验证实,这种滑移系形成于高温低应变速率条件下。橄榄石位错类型多样,（１００）、（００１）位错壁分布广泛,位错壁间距显示对数正态分布特征．估算的本区上地幔高温蠕交流动差异应力为（１０．９～１１．５ＭＰａ）。上地幔流变速率、流动应力、等效粘滞度与深度关系指示本区新生代大地构造性质为大陆拉张带或大陆裂谷带。     

Meso‐Cenozoic Tectonothermal History of Permian Strata,Southwestern Weibei Uplift: Insights from Thermochronology and Geothermometry

This study provides an integrated interpretation for the Mesozoic-Cenozoic tectonothermal evolutionary history of the Permian strata in the Qishan area of the southwestern Weibei Uplift, Ordos Basin. Apatite fission-track and apatite/zircon(U-Th)/He thermochronometry, bitumen reflectance, thermal conductivity of rocks, paleotemperature recovery, and basin modeling were used to restore the Meso-Cenozoic tectonothermal history of the Permian Strata. The Triassic AFT data have a pooled age of ～180±7 Ma with one age peak and P(χ2)=86%. The average value of corrected apatite(U-Th)/He age of two Permian sandstones is ～168±4 Ma and a zircon(U-Th)/He age from the Cambrian strata is ～231±14 Ma. Bitumen reflectance and maximum paleotemperature of two Ordovician mudstones are 1.81%, 1.57% and ～210°C, ～196°C respectively. After undergoing a rapid subsidence and increasing temperature in Triassic influenced by intrusive rocks in some areas, the Permian strata experienced four cooling-uplift stages after the time when the maximum paleotemperature reached in late Jurassic:(1) A cooling stage(～163 Ma to ～140 Ma) with temperatures ranging from ～132°C to ～53°C and a cooling rate of ～3°C/Ma, an erosion thickness of ～1900 m and an uplift rate of ～82 m/Ma;(2) A cooling stage(～140 Ma to ～52 Ma) with temperatures ranging from ～53°C to ～47°C and a cooling rate less than ～0.1°C/Ma, an erosion thickness of ～300 m and an uplift rate of ～3 m/Ma;(3)(～52 Ma to ～8 Ma) with ～47°C to ～43°C and ～0.1°C/Ma, an erosion thickness of ～500 m and an uplift rate of ～11 m/Ma;(3)(～8 Ma to present) with ～43°C to ～20°C and ～3°C/Ma, an erosion thickness of ～650 m and an uplift rate of ～81 m/Ma. The tectonothermal evolutionary history of the Qishan area in Triassic was influenced by the interaction of the Qinling Orogeny and the Weibei Uplift, and the south Qishan area had the earliest uplift-cooling time compared to other parts within the Weibei Uplift. The early Eocene at ～52 Ma and the late Miocene at ～8 Ma, as two significant turning points after which both the rate of uplift and the rate of temperature changed rapidly, were two key time for the uplift-cooling history of the Permian strata in the Qishan area of the southwestern Weibei Uplift, Ordos Basin.