Coarsening of three-dimensional structured and unstructured grids for subsurface flow

We present a generic, semi-automated algorithm for generating non-uniform coarse grids for modeling subsurface flow. The method is applicable to arbitrary grids and does not impose smoothness constraints on the coarse grid. One therefore avoids conventional smoothing procedures that are commonly used to ensure that the grids obtained with standard coarsening procedures are not too rough. The coarsening algorithm is very simple and essentially involves only two parameters that specify the level of coarsening. Consequently the algorithm allows the user to specify the simulation grid dynamically to fit available computer resources, and, e.g., use the original geomodel as input for flow simulations. This is of great importance since coarse grid-generation is normally the most time-consuming part of an upscaling phase, and therefore the main obstacle that has prevented simulation workflows with user-defined resolution. We apply the coarsening algorithm to a series of two-phase flow problems on both structured (Cartesian) and unstructured grids. The numerical results demonstrate that one consistently obtains significantly more accurate results using the proposed non-uniform coarsening strategy than with corresponding uniform coarse grids with roughly the same number of cells.     

Multi-scale modeling of three-phase–three-component processes in heterogeneous porous media

Flow and transport processes in porous media occur on different spatial and temporal scales and may also be locally different. Additionally, the structure of the porous medium itself generally shows a high dependence on the spatial scale.     

Flow and deposition of pyroclastic granular flows: A type example from the 1975 Ngauruhoe eruption,New Zealand

Small-volume pyroclastic density currents (PDCs) are generated frequently during explosive eruptions with little warning. Assessing their hazard requires a physical understanding of their transport and sedimentation processes which is best achieved by the testing of experimental and numerical models of geophysical mass flows against natural flows and/or deposits. To this end we report on one of the most detailed sedimentological studies ever carried out on a series of pristine small-volume PDC deposits from the 1975 eruption of Ngauruhoe volcano, whose emplacement were also witnessed during eruption. Using high-resolution GPS surveys, a series of lateral excavations across the deposits, and bulk sedimentological analysis we constrained the geomorphology, internal structure and texture of the deposits with respect to laterally varying modes of deposition.     

Magmatic structures in the Krkonoše–Jizera Plutonic Complex,Bohemian Massif: evidence for localized multiphase flow and small-scale thermal–mechanical instabilities in a granitic magma chamber

The present paper examines magmatic structures in the Jizera and Liberec granites of the Krkonoše–Jizera Plutonic Complex, Bohemian Massif. The magmatic structures are here interpreted to preserve direct field evidence for highly localized magma flow and other processes in crystal-rich mushes, and to capture the evolution of physical processes in an ancient granitic magma chamber. We propose that after chamber-wide mixing and hybridization, as suggested by recent petrological studies, laminar magma flow became highly localized to weaker channel-like domains within the higher-strength crystal framework. Mafic schlieren formed at flow rims, and their formation presumably involved gravitational settling and velocity gradient flow sorting coupled with interstitial melt escape. Local thermal or compositional convection may have resulted in the formation of vertical schlieren tubes and ladder dikes whereas subhorizontal tubes or channels formed during flow driven by lateral gradients in magma pressure. After the cessation or deceleration of channel flow, gravity-driven processes (settling of crystals and enclaves, gravitational differentiation, development of downward dripping instabilities), accompanied by compaction, filter pressing and melt segregation, dominated in the crystal mush within the flow channels. Subsequently, magmatic folds developed in schlieren layers and the magma chamber recorded complex, late magmatic strains at high magma crystallinities. Late-stage magma pulsing into localized submagmatic cracks represents the latest events of magmatic history of the chamber prior to its final crystallization. We emphasize that the most favorable environments for the formation and preservation of magmatic structures, such as those hosted in the Jizera and Liberec granites, are slowly cooling crystal-rich mushes. Therefore, where preserved in plutons, these structures may lend strong support for a “mush model” of magmatic systems.     

Magnetic petrology of the 1991–1995 dacite lava of Unzen volcano,Japan: Degree of oxidation and implications for the growth of lava domes

To understand the oxidation state and process of oxidation of lava domes, we carried out magnetic petrological analyses of lava samples obtained from domes and block-and-ash-flow deposits associated with the 1991–1995 eruption of Unzen volcano, Japan. As a result, we recognize three different types of magnetic petrology, each related to deuteric high-temperature oxidation during initial cooling. Type A oxides are characterized by homogenous titanomagnetite and titanohematite, indicating a low oxidation state and high titanomagnetite concentrations. Type B oxides are weakly exsolved and contain titanohematite laths and rutile lenses, indicating a higher oxidation state. Type C oxides, which represent the highest oxidation state, are completely exsolved and composed of Ti-poor titanomagnetite, titanohematite, rutile, and pseudobrookite, indicating high hematite concentrations. Some grains in Types A and B show indications of reduction, which was related to interaction with volcanic gases subsequent to high-temperature oxidation. In terms of geological occurrence, the oxidation processes probably differed for endogenous and exogenous domes. Endogenous dome lavas are oxidized concentrically and are classified into the three types according to their location within the dome: samples from the surface are strongly oxidized and classified as Type C, while the inner part is unoxidized and classified as Type A. Exogenous dome lavas are unoxidized and assigned to Type A. Some samples show signs of reduction, which may have occurred around fumaroles. We propose that location within the dome and the process of dome growth are the factors that control oxidation.     

Toward micro-scale spatial modeling of gentrification

 A simple preliminary model of gentrification is presented. The model is based on an irregular cellular automaton architecture drawing on the concept of proximal space, which is well suited to the spatial externalities present in housing markets at the local scale. The rent gap hypothesis on which the model's cell transition rules are based is discussed. The model's transition rules are described in detail. Practical difficulties in configuring and initializing the model are described and its typical behavior reported. Prospects for further development of the model are discussed. The current model structure, while inadequate, is well suited to further elaboration and the incorporation of other interesting and relevant effects. Received: 1 November 2001 / Accepted: 12 April 2002 Town centres data were made available by kind permission of the UK Department of the Environment, Transport and the Regions.     

南京城区夏季低层逆温对人体舒适状况影响初探

利用2000年夏季南京志区地面气象要素和低层探空资料，分析了该夏季南京城区早、晚低层逆温的频率、厚度和强度变化特征，着重从统计角度研究了逆温的厚度、强度与城区地面温度、相对湿度变化的关系及对人体舒适状况的影响。结果表明：低层逆温对城区地面温湿变化及其幅度影响较大，当07时有逆温层时，08时至14时的温度增加幅度变大，相对温度减小幅度也变大，人体不适指数增加幅度变大；当19时有逆温层时，20时至23时温度减小幅度变小，相对湿度增加幅度也变小，人体不适指数减小幅度变小，并通过相关分析建立了人体舒适度指数与逆温强度之间的回归方程。     

Design of the Urban Land Grading Information System

1　ＩｎｔｒｏｄｕｃｔｉｏｎＧｒａｄｉｎｇｕｒｂａｎｌａｎｄｃｏｒｒｅｃｔｌｙａｎｄｄｅｔｅｒｍｉｎｉｎｇｔｈｅｌａｎｄ ｐｒｉｃｅａｒｅｉｍｐｏｒｔａｎｔｂａｓｉｃｗｏｒｋｔｏｃａｒｒｙｏｕｔｔｈｅｃｏｍｐｅｎｓａｔｉｏｎａｌｅｘ ｐｌｏｉｔａｔｉｏｎａｎｄｄｅｖｅｌｏｐｔｈｅｒｅａｌｅｓｔａｔｅｍａｒｋｅｔ.Ｔｈｅｌａｎｄｇｒａｄｉｎｇｍｕｓｔｂｅｂｕｉｌｔｏｎｏｂｊｅｃｔｉｖｅｅｖａｌｕａｔｉｏｎａｎｄｓｃｉｅｎｔｉｆｉｃｅｖａｌｕａｔｉｏｎａｎｄｃａｌｃｕｌａｔｉｏｎｏｆｌａｎｄ ｑｕａｌｉ…     

非定常不稳定液固两相流动中旋涡对颗粒运动影响的数值研究

构建起双向耦合的液固两相流动旋涡动力学模型与数值方法;应用离散涡方法,计算非定常不稳定水流场;采用Lagrange方法模拟颗粒运动,颗粒对流体的反作用通过修正涡泡运动速度来实现。利用所建模型,计算了两种St数的泥沙粒子在圆柱绕流场中的运动。结果证明了液固两相流动中颗粒运动与旋涡存在着明确的相关结构:(1)当水沙混合物中的泥沙颗粒碰上旋涡时,泥沙颗粒被卷入旋涡中,被卷入旋涡中的泥沙颗粒在运动过程中始终分布于旋涡区;(2)均匀水沙混合物绕圆柱流动,由于流体流过圆柱时产生剧烈分离流动,使得在尾迹流内中等St数 (St～o (1))的泥沙颗粒从均匀水沙混合物中分离出来而往旋涡区聚集。     

长江口九段沙下段冲淤演变水动力机制分析

概述了长江口深水航道整治一期工程及工程对九段沙下段地形冲淤的影响.根据已建立的σ坐标系下三维非线性水流数学模型,用新测水文资料对模型进行验证.利用所建立的模型,结合由现场测量得到一期工程前后九段沙下段近期地形演变,数值模拟了北槽区域流通量、底层欧拉余流和北槽中下段平面水流特征,以综合分析北槽水域水动力变化对九段沙下段地形冲淤变化的影响.