Numerical Calculation of Equivalent Cell Permeability Tensors for General Quadrilateral Control Volumes

A new method for upscaling fine scale permeability fields to general quadrilateral-shaped coarse cells is presented. The procedure, referred to as the conforming scale up method, applies a triangle-based finite element technique, capable of accurately resolving both the coarse cell geometry and the subgrid heterogeneity, to the solution of the local fine scale problem. An appropriate averaging of this solution provides the equivalent permeability tensor for the coarse scale quadrilateral cell. The general level of accuracy of the technique is demonstrated through application to a number of flow problems. The real strength of the conforming scale up method is demonstrated when the method is applied in conjunction with a flow-based gridding technique. In this case, the approach is shown to provide results that are significantly more accurate than those obtained using standard techniques.     

Phase relationships in Buchan facies series pelitic assemblages: calculations with application to andalusite-staurolite parageneses in the Mount Lofty Ranges,South Australia

Low-pressure, medium- to high-temperature (Buchan-type) regional metamorphism of pelitic rocks in the Mount Lofty Ranges, South Australia, is defined by the development of biotite, staurolite-andalusite, fibrolite, prismatic sillimanite and migmatite zones. K-feldspar makes its first appearance in the prismatic sillimanite zone and here we restrict our discussion to lower grade assemblages containing prograde muscovite, concentrating particularly on well-developed andalusitestaurolite parageneses. In general, the spatial distribution and mineral chemical variation of these assemblages accord with the predictions of petrogenetic grids and P-T and T-X _Fe pseudo-sections constructed from the internally consistent data set of Holland and Powell (1990) in the system KFMASH, assuming a(H₂O) 1, although analysed white mica compositions are systematically more aluminous than predicted. Importantly, the stability ranges of most critical assemblages predicted by these grids and pseudo-sections coincide closely with P-T estimates calculated using the data set of Holland and Powell (1990) from the Mount Lofty Ranges assemblages. With the exception of Mn in garnet and Zn in one staurolite-cordierite-muscovite assemblage non-KFMASH components do not significantly appear to have affected the stability ranges of the observed assemblages. An apparent local reversal in isograd zonation in which andalusite first appears down-grade of staurolite suggests a metamorphic field gradient concave towards the temperature axis and, together with evidence for essentially isobaric heating of individual rocks, is consistent with the exposures representing an oblique profile through a terrain in which heat was dissipated from intrusive bodies at discrete structural levels.Mineral abbreviations used in figures als Al₂SiO₅ phase - bi biotite - chl chlorite - ky kyanite - ph phengite - sill sillimanite - and andalusite - cd cordieritc - gt garnet - mu muscovite - q quartz - st staurolite     

The Kaidun meteorite: Mineralogy of an unusual CM1 lithology

Abstract Kaidun is a breccia of disparate enstatite and carbonaceous chondrite clasts that continues to provide real surprises. Many Kaidun clasts have been intensely altered by aqueous fluids, as evidenced by the widespread occurrence of ferromagnesian phyllosilicates and by the presence of carbonate- and phyllo-silicate-filled veins. In this report, we describe an unusual CM lithology containing many mineralogical features not previously reported from any meteorite, including pyrrhotite, with exclusive needlelike morphologies and thick mantles of phyllosilicate, and complex aggregates of phyllosilicate, melanite garnet, crosscut by pentlandite veins. The latter features appear to be due in large part to extensive hydrothermal alteration at temperatures on the order of 450 °C, which is significantly higher than that attained during secondary processing from other known CM material.     

Petrophysics and mineral exploration: a workflow for data analysis and a new interpretation framework

As mineral exploration seeks deeper targets, there will be a greater reliance on geophysical data and a better understanding of the geological meaning of the responses will be required, and this must be achieved with less geological control from drilling. Also, exploring based on the mineral system concept requires particular understanding of geophysical responses associated with altered rocks. Where petrophysical datasets of adequate sample size and measurement quality are available, physical properties show complex variations, reflecting the combined effects of various geological processes. Large datasets, analysed as populations, are required to understand the variations. We recommend the display of petrophysical data as frequency histograms because the nature of the data distribution is easily seen with this form of display. A petrophysical dataset commonly contains a combination of overlapping sub-populations, influenced by different geological factors. To understand the geological controls on physical properties in hard rock environments, it is necessary to analyse the petrophysical data not only in terms of the properties of different rock types. It is also necessary to consider the effects of processes such as alteration, weathering, metamorphism and strain, and variables such as porosity and stratigraphy. To address this complexity requires that much more supporting geological information be acquired than in current practice. The widespread availability of field portable instruments means quantitative geochemical and mineralogical data can now be readily acquired, making it unnecessary to rely primarily on categorical rock classification schemes. The petrophysical data can be combined with geochemical, petrological and mineralogical data to derive explanations for observed physical property variations based not only on rigorous rock classification methods, but also in combination with quantitative estimates of alteration and weathering. To understand how geological processes will affect different physical properties, it is useful to define three end-member forms of behaviour. Bulk behaviour depends on the physical properties of the dominant mineral components. Density and, to a lesser extent, seismic velocity show such behaviour. Grain and texture behaviour occur when minor components of the rock are the dominate controls on its physical properties. Grain size and shape control grain properties, and for texture properties the relative positions of these grains are also important. Magnetic and electrical properties behave in this fashion. Thinking in terms of how geological processes change the key characteristics of the major and minor mineralogical components allows the resulting changes in physical properties to be understood and anticipated.     

A three-dimensional,iterative mapping procedure for the implementation of an ionosphere-magnetosphere anisotropic Ohm’s law boundary condition in global magnetohydrodynamic simulations

The mathematical formulation of an iterative procedure for the numerical implementation of an ionosphere-magnetosphere (IM) anisotropic Ohm’s law boundary condition is presented. The procedure may be used in global magnetohydrodynamic (MHD) simulations of the magnetosphere. The basic form of the boundary condition is well known, but a well-defined, simple, explicit method for implementing it in an MHD code has not been presented previously. The boundary condition relates the ionospheric electric field to the magnetic field-aligned current density driven through the ionosphere by the magnetospheric convection electric field, which is orthogonal to the magnetic field B, and maps down into the ionosphere along equipotential magnetic field lines. The source of this electric field is the flow of the solar wind orthogonal to B. The electric field and current density in the ionosphere are connected through an anisotropic conductivity tensor which involves the Hall, Pedersen, and parallel conductivities. Only the height-integrated Hall and Pedersen conductivities (conductances) appear in the final form of the boundary condition, and are assumed to be known functions of position on the spherical surface R=R₁ representing the boundary between the ionosphere and magnetosphere. The implementation presented consists of an iterative mapping of the electrostatic potential , the gradient of which gives the electric field, and the field-aligned current density between the IM boundary at R=R₁ and the inner boundary of an MHD code which is taken to be at R₂>R₁. Given the field-aligned current density on R=R₂, as computed by the MHD simulation, it is mapped down to R=R₁ where it is used to compute by solving the equation that is the IM Ohm’s law boundary condition. Then is mapped out to R=R₂, where it is used to update the electric field and the component of velocity perpendicular to B. The updated electric field and perpendicular velocity serve as new boundary conditions for the MHD simulation which is then used to compute a new field-aligned current density. This process is iterated at each time step. The required Hall and Pedersen conductances may be determined by any method of choice, and may be specified anew at each time step. In this sense the coupling between the ionosphere and magnetosphere may be taken into account in a self-consistent manner.     

An orthotropic interface damage model for simulating drying processes in soils

The study of drying process in soils has received an increased attention in the last few years. This is very complex phenomenon that generally leads to the formation and propagation of desiccation cracks in the soil mass. In recent engineering applications, high aspect ratio elements have proved to be well suited to tackle this type of problem using finite elements. However, the modeling of interfaces between materials with orthotropic properties that generally exist in this type of problem using standard (isotropic) constitutive model is very complex and challenging in terms of the mesh generation, leading to very fine meshes that are intensive CPU demanding. A novel orthotropic interface mechanical model based on damage mechanics and capable of dealing with interfaces between materials in which the strength depends on the direction of analysis is proposed in this paper. The complete mathematical formulation is presented together with the algorithm suggested for its numerical implementation. Some simple yet challenging synthetic benchmarks are analyzed to explore the model capabilities. Laboratory tests using different textures at the contact surface between materials were conducted to evaluate the strengths of the interface in different directions. These experiments were then used to validate the proposed model. Finally, the approach is applied to simulate an actual desiccation test involving an orthotropic contact surface. In all the application cases the performance of the model was very satisfactory.