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.     

P‐wave transmission across fractures with nonlinear deformational behaviour

Stress wave attenuation across fractured rock masses is a great concern of underground structure safety. When the wave amplitude is large, fractures experience nonlinear deformation during the wave propagation. This paper presents a study on normal transmission of P‐wave across parallel fractures with nonlinear deformational behaviour (static Barton–Bandis model). The results show that the magnitude of transmission coefficient is a function of incident wave amplitude, nondimensional fracture spacing and number of fractures. Two important indices of nondimensional fracture spacing are identified, and they divide the area of nondimensional fracture spacing into three parts (individual fracture area, transition area and small spacing area). In the different areas, the magnitude of transmission coefficient has different trends with nondimensional fracture spacing and number of fractures. In addition, the study reveals that under some circumstances, the magnitude of transmission coefficient increases with increasing number of fractures, and is larger than 1. Copyright © 2006 John Wiley & Sons, Ltd.     

Porphyroblast inclusion-trail orientation data: eppure non son girate*!

Extensive examination of large numbers of spatially orientated thin sections of orientated samples from orogens of all ages around the world has demonstrated that porphyroblasts do not rotate relative to geographical coordinates during highly non-coaxial ductile deformation of the matrix subsequent to their growth. This has been demonstrated for all tectonic environments so far investigated. The work also has provided new insights and data on metamorphic, structural and tectonic processes including: (1) the intimate control of deformation partitioning on metamorphic reactions; (2) solutions to the lack of correlation between lineations that indicate the direction of movement within thrusts and shear zones, and relative plate motion; and (3) a possible technique for determining the direction of relative plate motion that caused orogenesis in ancient orogens.     

Reliability assessment of granular filters in embankment dams

Empirical criteria have been used successfully to design filters of most embankment large dam projects throughout the world. However, these empirical rules are only applicable to a particular range of soils tested in laboratory and do not take into account the variability of the base material and filter particle sizes. In addition, it is widely accepted that the safety of fill dams is mainly dependent on the reliability of their filter performance. The work herein presented consists in a new general method for assessing the probability of fulfilling any empirical filter design criteria accounting for base and filter heterogeneity by means of first‐order reliability methods (FORM), so that reliability indexes and probabilities of fulfilling any particular criteria are obtained. This method will allow engineers to estimate the safety of existing filters in terms of probability of fulfilling their design criteria and might also be used as a decision tool on sampling needs and material size tolerances during construction. In addition, sensitivity analysis makes possible to analyse how reliabilities are influenced by different sources of input data. Finally, in case of a portfolio risk assessment, this method will allow engineers to compare the safety of several existing dams in order to prioritize safety investments and it is expected to be a very useful tool to evaluate probabilities of failure due to internal erosion. Copyright © 2006 John Wiley & Sons, Ltd.     

Chemical alteration of tephra in the depositional environment: theoretical stability modelling

The study of the chemical stability of vitreous material in aqueous media is well‐established. There has to date been little consideration of the implications of variations in the chemical durability of tephra in Quaternary tephrochronology. Chemical alteration can take the form of cationic leaching from the matrix, or complete destruction of the silica network, either of which could constrain the ability to chemically identify distal tephra. Here we apply established models of vitreous durability to the published chemical analyses of a large number of Icelandic tephras in order to predict their relative durabilities under equivalent conditions. This suggests that some important tephras have relatively poor chemical stability, and that rhyolitic tephras are, in general, more stable than basaltic. We conclude that tephras should be expected to show predictable differential chemical stability in the post‐depositional environment. Copyright © 2003 John Wiley & Sons, Ltd.     

Observations on the ecological importance of salt marshes in the Cumberland Basin,a macrotidal estuary in the Bay of Fundy

The Cumberland Basin, a 118 km² estuary at the head of the Bay of Fundy which has an average tidal range of about 11m, contains large tracts of salt marsh (15% of the area below highest high water). Low marsh (below about 0·9 m above mean high water) is composed almost exclusively of Spartina alterniflora while the vegetation on high marsh is more diverse but dominated by Spartina patens. Because of its higher elevation, high marsh is flooded infrequently for short periods by only extreme high tides. Low marsh is inundated much more frequently by water as much as 4m deep for periods as long as 4 h per tide. Temporal variability in the occurrence of extreme tides influences the flooding frequency of high marsh for any given month and year. Using a modification of Smalley's method, the mean annual net aerial primary production (NAPP) of low and high marsh is estimated to be 272 and 172 g C m^?2, respectively. Vegetation turnover times average 1·0 and 2·0 y for low and high marsh, respectively. Because of abundant tidal energy, much of the low marsh production appears to be exported and distributed widely about the estuary. Since high levels of turbidity suppress phytoplankton production, salt marshes produce approximately half of the carbon fixed photosynthetically in the Cumberland Basin. It is concluded that salt marshes play a major ecological role in the Cumberland Basin.     

Late Holocene great earthquakes and relative sea‐level change at Kenai,southern Alaska

Kenai, located on the west coast of the Kenai Peninsula, Alaska, subsided during the great earthquake of AD 1964. Regional land subsidence is recorded within the estuarine stratigraphy as peat overlain by tidal silt and clay. Reconstructions using quantitative diatom transfer functions estimate co‐seismic subsidence (relative sea‐level rise) between 0.28±0.28 m and 0.70±0.28 m followed by rapid post‐seismic recovery. Stratigraphy records an earlier co‐seismic event as a second peat‐silt couplet, dated to ～1500–1400 cal. yr BP with 1.14±0.28 m subsidence. Two decimetre‐scale relative sea‐level rises are more likely the result of glacio‐isostatic responses to late Holocene and Little Ice Age glacier expansions rather than to co‐seismic subsidence during great earthquakes. Comparison with other sites around Cook Inlet, at Girdwood and Ocean View, helps in constructing regional patterns of land‐level change associated with three great earthquakes, AD 1964, ～950–850 cal. yr BP and ～1500–1400 cal. yr BP. Each earthquake has a different spatial pattern of co‐seismic subsidence which indicates that assessment of seismic hazard in southern Alaska requires an understanding of multiple great earthquakes, not only the most recent. All three earthquakes show a pre‐seismic phase of gradual land subsidence that marked the end of relative land uplift caused by inter‐seismic strain accumulation. Copyright © 2005 John Wiley & Sons, Ltd.     

Survival,growth and shell deformities of copper- and cadmium-exposed mussels (Mytilus edulis L.) in brackish water

Mussels, Mytilus edulis L., were exposed to elevated concentrations of copper or cadmium in the laboratory, then placed in cages in the sea (salinity 7‰). One year later maximum lengths of the mussels were measured and shells screened for deformities. Growth was 0·6 cm year^?1 in the control cages and retarded in cages of exposed mussels. A total of 63% of cadmium-exposed and 46% of copper-exposed mussels had shell deformities. In the control cages 26% were deformed while in a natural population only 3% were deformed. The proportion of deformities to growth was inversely related. Low concentrations showed their injuriousness in this long-term test, probably due to the closing of the valves during exposure to high concentrations.     

Numerical prediction of blast‐induced stress wave from large‐scale underground explosion

This paper presents a numerical model for predicting the dynamic response of rock mass subjected to large‐scale underground explosion. The model is calibrated against data obtained from large‐scale field tests. The Hugoniot equation of state for rock mass is adopted to calculate the pressure as a function of mass density. A piecewise linear Drucker–Prager strength criterion including the strain rate effect is employed to model the rock mass behaviour subjected to blast loading. A double scalar damage model accounting for both the compression and tension damage is introduced to simulate the damage zone around the charge chamber caused by blast loading. The model is incorporated into Autodyn3D through its user subroutines. The numerical model is then used to predict the dynamic response of rock mass, in terms of the peak particle velocity (PPV) and peak particle acceleration (PPA) attenuation laws, the damage zone, the particle velocity time histories and their frequency contents for large‐scale underground explosion tests. The computed results are found in good agreement with the field measured data; hence, the proposed model is proven to be adequate for simulating the dynamic response of rock mass subjected to large‐scale underground explosion. Extended numerical analyses indicate that, apart from the charge loading density, the stress wave intensity is also affected, but to a lesser extent, by the charge weight and the charge chamber geometry for large‐scale underground explosions. Copyright © 2004 John Wiley & Sons, Ltd.     

Discrimination of sedimentary facies by association analysis

The ability of association analysis to discriminate sedimentary facies was tested on Purdy's modal analyses of modern sediments of the Great Bahama Bank. Purdy's data set has served in the past as a standard reference for evaluating various multivariate classification algorithms. In order to adapt Purdy's data to association analysis, the percent abundance of the 12 constituents was converted to binary form by dichotomizing each variable on its mean value. The results obtained by association analysis are virtually identical to those obtained by Purdy and other authors. The same four main sedimentary facies were discriminated; 86% of the samples were identically classified (97% when misclassified borderline cases are counted as matches); the total partition variance of the classification is only negligibly greater (4%); and the grouping of the variables yielded the same four groups. The rank order of the three division-attributes responsible for the sample classification is fines, oolites, and corals. Association analysis has been employed by other authors to differentiate meaningful facies groups in studies of ancient reef carbonates, modern reef sediments, and heavy minerals in stream sediments. In all these studies, the results were found to be compatible with those obtained by using the continuous quantitative measurements, indicating that qualitative binary data may often be sufficient for the purpose of facies discrimination in many branches of geology and that association analysis is an effective method for this purpose.