A discretization and multigrid solver for a Darcy–Stokes system of three dimensional vuggy porous media

We develop a finite element discretization and multigrid solver for a Darcy–Stokes system of three-dimensional vuggy porous media, i.e., porous media with cavities. The finite element method uses low-order mixed finite elements in the Darcy and Stokes domains and special transition elements near the Darcy–Stokes interface to allow for tangential discontinuities implied by the Beavers–Joseph boundary condition. We design a multigrid method to solve the resulting saddle point linear system. The intertwining of the Darcy and Stokes subdomains makes the resulting matrix highly ill-conditioned. The velocity field is very irregular, and its discontinuous tangential component at the Darcy–Stokes interface makes it difficult to define intergrid transfer operators. Our definition is based on mass conservation and the analysis of the orders of magnitude of the solution. The coarser grid equations are defined using the Galerkin method. A new smoother of Uzawa type is developed based on taking an optimal step in a good search direction. Our algorithm has a measured convergence factor independent of the size of the system, at least when there are no disconnected vugs. We study the macroscopic effective permeability of a vuggy medium, showing that the influence of vug orientation; shape; and, most importantly, interconnectivity determine the macroscopic flow properties of the medium. This work was supported by the U.S. National Science Foundation under grants DMS-0074310 and DMS-0417431.     

Large Aperture Scintillometer Intercomparison Study

Two field studies with six large aperture scintillometers (LASs) were performed using horizontal and slant paths. The accuracy of this novel and increasingly popular technique for measuring sensible heat fluxes was quantified by comparing measurements from different instruments over nearly identical transects. Random errors in LAS measurements were small, since correlation coefficients between adjacent measurements were greater than 0.995. However, for an ideal set-up differences in linear regression slopes of up to 21% were observed with typical inter-instrument differences of 6%. Differences of 10% are typical in more realistic measurement scenarios over homogeneous natural vegetation and different transect heights and locations. Inaccuracies in the optics, which affect the effective aperture diameter, are the most likely explanation for the observed differences.     

Fluid–rock interaction at the roots of an Eopaleozoic epithermal system, Sierra Norte de Córdoba, Argentina: Paragenesis, thermometry and fluid sources

A deep epithermal vein system hosted in Late Proterozoic to Cambrian granodiorite has been identified in the Sierra Norte de Córdoba, the easternmost range of the Sierras Pampeanas Orientales of Argentina. The vein swarm extends over an area of 3 km² parallel to a mylonitic belt and formed in fractured granodiorite. Thicknesses of veins are less than 0.5 m and their visible strike length is less than 100 m. Veins are either barren or weakly mineralized in base-metal sulfides. Most veins have mineral associations dominated by calcite and quartz with lesser amounts of chlorite, sericite, pyrite, and minor illite. In other less exposed albite-rich, adularia-bearing veins, chalcopyrite, bornite, galena, sphalerite, chalcocite and covellite may occur. The widespread occurrence of bladed calcite without any petrographic or microthermometric evidence of boiling implies that this particular habit of calcite may also develop under sub-near boiling fluid conditions. Thermometric calculations based on fluid inclusion data, chlorite composition and oxygen isotopes in the quartz–calcite pair, constrain the formation of the system between 300 and 350 °C, at pressures between 42 and 64 MPa (1.5–2.3 km). Stable isotope data suggest that W/R interaction might have been the most probable mechanism of alteration, involving the participation of meteoric fluids; nevertheless, the metallic signature of some weakly mineralized veins as well as intermediate fluid inclusion salinities favor a magmatic input and a mixed origin for the fluids. Textures and mineral associations, as well as the absence of evidence of boiling in fluid inclusions, all suggest that the silica–carbonate vein system formed deeper than typically shallow Au and Ag-bearing boiling solutions. A 485 (±25) Ma lamprophyre dike crosscuts some of these veins locally producing metasomatic reactions and skarn formation, which constrains the age of the hydrothermal system to the Cambrian-Early Ordovician time span.     

On-line recognition of supernova neutrino bursts in the LVD

In this paper we show the capabilities of the Large Volume Detector (INFN Gran Sasso National Laboratory) to identify a neutrino burst associated with a supernova explosion, in the absence of an “external trigger”, e.g., an optical observation. We describe how the detector trigger and event selection have been optimized for this purpose, and we detail the algorithm used for the on-line burst recognition. The on-line sensitivity of the detector is defined and discussed in terms of supernova distance and intensity at the source.     

Explicit Fracture Network Modelling: From Multiple Point Statistics to Dynamic Simulation

This paper addresses the problem of explicit fractured media modelling in an operational case. On one side, realistic fracture models are mainly used for research purposes in order to investigate better the flow behaviour impacted by the complex multi-scale fracture network. Often, a very fine grid and hence an increased computation time are needed. On the other hand, an operational fractured reservoir is still generally modelled using an implicit fracture media representation. The upscaled petrophysical properties and dual media are defined on a coarse grid to limit the computational time of dynamic simulation. The challenge of this work is to demonstrate that an explicit fracture modelling is not reserved only for the research domain, but can be applied to an operational case study. The static model is constructed using a multiple point statistics approach in order to represent complex interaction patterns of fractures and faults observed at the analogue outcrop. The dynamic behaviour is simulated based on this spatial fracture network representation.     

The effects of acoustic-gravity waves on the K i 7699 line

We examine the effects of acoustic-gravity waves with long and short periods on the solar profile of the K i7699 line using a dynamic model of line formation.First we studied the kinetic equilibrium of the K i atoms in a static atmosphere confirming, with up-to-date atomic data and atmospheric model, that a good fit of the resonance line 7699 is possible only when non-LTE effects are accounted for.Then the static non-LTE line source function and lower-level population are used as input data for calculating the line formation in the presence of waves.The time behaviour of the synthetic profiles corresponding to 300 s and 30 s waves is extensively discussed. The characteristic redshift induced by the 30 s wave is explained within the framework of the S-S line formation model. Long-period waves yield an anticorrelation between the asymmetry at different residual intensities and the line core shift, as observed. The short-period waves with velocity amplitude of about 100 m s^–1 (at the base of the photosphere) produce a mean bisector whose lower part has a slope in agreement with the observed one. The efficiency of waves to produce macro and microturbulence is also discussed.Currently NAS/NRC Senior Research Associate at Sacramento Peak Obs., New Mexico, U.S.A.     

Review: Hydrogeology of Northern Mendoza (Argentina), from the Andes to the eastern plains,in the context of climate change

Hydrogeology Journal - In the drylands of Northern Mendoza, Argentina, water supply depends on rivers and groundwater. Climate change makes this region vulnerable due to the snow-glacial-melt...     

Geotechnical zonation and soil–structure interaction at Puerto Vallarta,México

Natural Hazards - We performed a seismic vulnerability assessment that involves geotechnical and building structure analysis for Puerto Vallarta, Mexico, a city located along the pacific coast....