Open-source MATLAB implementation of consistent discretisations on complex grids   总被引：1，自引：0，他引：1  

Knut–Andreas Lie  Stein Krogstad  Ingeborg Skjelkv?le Ligaarden  Jostein Roald Natvig  Halvor M?ll Nilsen  B?rd Skaflestad 《Computational Geosciences》2012,16(2):297-322

Accurate geological modelling of features such as faults, fractures or erosion requires grids that are flexible with respect to geometry. Such grids generally contain polyhedral cells and complex grid-cell connectivities. The grid representation for polyhedral grids in turn affects the efficient implementation of numerical methods for subsurface flow simulations. It is well known that conventional two-point flux-approximation methods are only consistent for K-orthogonal grids and will, therefore, not converge in the general case. In recent years, there has been significant research into consistent and convergent methods, including mixed, multipoint and mimetic discretisation methods. Likewise, the so-called multiscale methods based upon hierarchically coarsened grids have received a lot of attention. The paper does not propose novel mathematical methods but instead presents an open-source Matlab? toolkit that can be used as an efficient test platform for (new) discretisation and solution methods in reservoir simulation. The aim of the toolkit is to support reproducible research and simplify the development, verification and validation and testing and comparison of new discretisation and solution methods on general unstructured grids, including in particular corner point and 2.5D PEBI grids. The toolkit consists of a set of data structures and routines for creating, manipulating and visualising petrophysical data, fluid models and (unstructured) grids, including support for industry standard input formats, as well as routines for computing single and multiphase (incompressible) flow. We review key features of the toolkit and discuss a generic mimetic formulation that includes many known discretisation methods, including both the standard two-point method as well as consistent and convergent multipoint and mimetic methods. Apart from the core routines and data structures, the toolkit contains add-on modules that implement more advanced solvers and functionality. Herein, we show examples of multiscale methods and adjoint methods for use in optimisation of rates and placement of wells.  相似文献   

The L3–6 chondritic regolith breccia Northwest Africa (NWA) 869: (I) Petrology,chemistry, oxygen isotopes,and Ar‐Ar age determinations     

Knut METZLER  Addi BISCHOFF  Richard C. GREENWOOD  Herbert PALME  Marko GELLISSEN  Jens HOPP  Ian A. FRANCHI  Mario TRIELOFF 《Meteoritics & planetary science》2011,46(5):652-680

Abstract– Northwest Africa (NWA) 869 consists of thousands of individual stones with an estimated total weight of about 7 metric tons. It is an L3–6 chondrite and probably represents the largest sample of the rare regolith breccias from the L–chondrite asteroid. It contains unequilibrated and equilibrated chondrite clasts, some of which display shock‐darkening. Impact melt rocks (IMRs), both clast‐free and clast‐poor, are strongly depleted in Fe,Ni metal, and sulfides. An unequilibrated microbreccia, two different light inclusions and two different SiO₂‐bearing objects were found. Although the matrix of this breccia appears partly clastic, it is not a simple mixture of fine‐grained debris formed from the above lithologies, but mainly represents an additional specific lithology of low petrologic type. We speculate that this material stems from a region of the parent body that was only weakly consolidated. One IMR clast and one SiO₂‐bearing object show Δ¹⁷O values similar to bulk NWA 869, suggesting that both are related to the host rock. In contrast, one light inclusion and one IMR clast appear to be unrelated to NWA 869, suggesting that the IMR clast is contaminated with impactor material. ⁴⁰Ar‐³⁹Ar analyses of a type 4 chondrite clast yield a plateau age of 4402 ± 7 Ma, which is interpreted to be the result of impact heating. Other impact events are recorded by an IMR clast at 1790 ± 36 Ma and a shock‐darkened clast at 2216 ± 40 Ma, demonstrating that NWA 869 escaped major reset in the course of the event at approximately 470 Ma that affected many L–chondrites.  相似文献   

Thermal history of Northwest Africa 5073––A coarse‐grained Stannern‐trend eucrite containing cm‐sized pyroxenes and large zircon grains     

Julia ROSZJAR  Knut METZLER  Addi BISCHOFF  Jean‐Alix BARRAT  Thorsten GEISLER  Richard C. GREENWOOD  Ian A. FRANCHI  Stephan KLEMME 《Meteoritics & planetary science》2011,46(11):1754-1773

Abstract– We report on the bulk chemical composition, petrology, oxygen isotopic composition, trace element composition of silicates, and degree of self‐irradiation damage on zircon grains of the eucrite Northwest Africa (NWA) 5073 to constrain its formation and postcrystallization thermal history, and to discuss their implications for the geologic history of its parent body. This unequilibrated and unbrecciated meteorite is a new member of the rare Stannern‐trend eucrites. It is mainly composed of elongated, zoned pyroxene phenocrysts up to 1.2 cm, plagioclase laths up to 0.3 cm in length, and is rich in mesostasis. The latter contains zircon grains up to 30 μm in diameter, metal, sulfide, tridymite, and Ca‐phosphates. Textural observations and silicate compositions, coupled with the occurrence of extraordinary Fe‐rich olivine veins that are restricted to large pyroxene laths, indicate that NWA 5073 underwent a complex thermal history. This is also supported by the annealed state of zircon grains inferred from μ‐Raman spectroscopic measurements along with U and Th data obtained by electron probe microanalyses.  相似文献   

Chemistry and oxygen isotopic composition of cluster chondrite clasts and their components in LL3 chondrites          下载免费PDF全文

Knut Metzler  Andreas Pack 《Meteoritics & planetary science》2016,51(2):276-302

Cluster chondrites are characterized by close‐fit textures of deformed and indented chondrules, taken as evidence for hot chondrule accretion (Metzler 2012 ). We investigated seven cluster chondrite clasts from six brecciated LL3 chondrites and measured their bulk oxygen isotopic and chemical composition, including REE, Zr, and Hf. The same parameters were measured in situ on 93 chondrules and 4 interchondrule matrix areas. The CI‐normalized REE patterns of the clasts are flat, showing LL‐chondritic concentrations. The mean chemical compositions of chondrules in clasts and other LL chondrites are indistinguishable and we conclude that cluster chondrite chondrules are representative of the normal LL chondrule population. Type II chondrules are depleted in MgO, Al₂O₃ and refractory lithophiles (REE, Zr, Hf) by factors between 0.65 and 0.79 compared to type I chondrules. The chondrule REE patterns are basically flat with slight LREE < HREE fractionations. Many chondrules exhibit negative Eu anomalies while matrix shows a complementary pattern. Chondrules scatter along a correlation line with a slope of 0.63 in the oxygen 3‐isotope diagram, interpreted as the result of O‐isotope exchange between chondrule melts and ¹⁸O‐rich nebular components. In one clast, a distinct anticorrelation between chondrule size and δ¹⁸O is found, which may indicate a more intense oxygen isotope exchange by smaller chondrules. In some clasts the δ¹⁸O values of type I chondrules are correlated with concentrations of SiO₂ and MnO and anticorrelated with MgO, possibly due to the admixture of a SiO₂‐ and MnO‐rich component to chondrule melts during oxygen isotope exchange. Two chondrules with negative anomalies in Sm, Eu, and Yb were found and may relate their precursors to refractory material known from group III CAIs. Furthermore, three chondrules with strong LREE > HREE and Zr/Hf fractionations were detected, whose formation history remains to be explained.  相似文献