A parallel implementation of the lattice solid model for the simulation of rock mechanics and earthquake dynamics

The Lattice Solid Model has been used successfully as a virtual laboratory to simulate fracturing of rocks, the dynamics of faults, earthquakes and gouge processes. However, results from those simulations show that in order to make the next step towards more realistic experiments it will be necessary to use models containing a significantly larger number of particles than current models. Thus, those simulations will require a greatly increased amount of computational resources. Whereas the computing power provided by single processors can be expected to increase according to Moore's law, i.e., to double every 18-24 months, parallel computers can provide significantly larger computing power today. In order to make this computing power available for the simulation of the microphysics of earthquakes, a parallel version of the Lattice Solid Model has been implemented. Benchmarks using large models with several millions of particles have shown that the parallel implementation of the Lattice Solid Model can achieve a high parallel-efficiency of about 80% for large numbers of processors on different computer architectures.

Information flow analysis for fail-secure devices

Information security devices must preserve security properties even in the presence of faults. This in turn requires a rigorous evaluation of the system behaviours resulting from component failures, especially how such failures affect information flow. We introduce a compositional method of static analysis for fail-secure behaviour. Our method uses reachability matrices to identify potentially undesirable information flows based on the fault modes of the system's components.

Dynamic rupture in a 3-D particle-based simulation of a rough planar fault

An appreciation of the physical mechanisms which cause observed seismicity complexity is fundamental to the understanding of the temporal behaviour of faults and single slip events. Numerical simulation of fault slip can provide insights into fault processes by allowing exploration of parameter spaces which influence microscopic and macroscopic physics of processes which may lead towards an answer to those questions. Particle-based models such as the Lattice Solid Model have been used previously for the simulation of stick-slip dynamics of faults, although mainly in two dimensions. Recent increases in the power of computers and the ability to use the power of parallel computer systems have made it possible to extend particle-based fault simulations to three dimensions. In this paper a particle-based numerical model of a rough planar fault embedded between two elastic blocks in three dimensions is presented. A very simple friction law without any rate dependency and no spatial heterogeneity in the intrinsic coefficient of friction is used in the model. To simulate earthquake dynamics the model is sheared in a direction parallel to the fault plane with a constant velocity at the driving edges. Spontaneous slip occurs on the fault when the shear stress is large enough to overcome the frictional forces on the fault. Slip events with a wide range of event sizes are observed. Investigation of the temporal evolution and spatial distribution of slip during each event shows a high degree of variability between the events. In some of the larger events highly complex slip patterns are observed.

In situ observation of bainite growth during isothermal holding

The mechanism of bainite growth has been investigated using in situ transmission electron microscopy observations. It was found that, in a number of alloys studied, a bainitic embryo is made of basic transformation units. These units are either a group of stacking faults or, in two dimensions, a series of parallelograms of different sizes. Thickening/widening of the bainite embryo takes place through shear along the stacking fault planes or twining planes. The bainite embryo is elongated by the formation of new transformation units at both tips of the bainite plate. The three-dimensional morphology of bainite is a convex tens-like lath. It is believed that the bainite embryo grows by shearing, which is controlled by the diffusion of solute atoms during the transformation. As the growth rate is much lower than that of martensite, it is therefore detectable. (c) 2006 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Multiagent realization of prediction-based diagnosis and loss prevention

A multiagent diagnostic system implemented in a Protege-JADE-JESS environment interfaced with a dynamic simulator and database services is described in this paper. The proposed system architecture enables the use of a combination of diagnostic methods from heterogeneous knowledge sources. The process ontology and the process agents are designed based on the structure of the process system, while the diagnostic agents implement the applied diagnostic methods. A specific completeness coordinator agent is implemented to coordinate the diagnostic agents based on different methods. The system is demonstrated on a case study for diagnosis of faults in a granulation process based on HAZOP and FMEA analysis.

Standardising bull breeding soundness evaluations and reporting in Australia

There is substantial variation in bull breeding soundness evaluation procedures and reports in Australia; the situation is compounded by difficulties in interpretation and the validity of many reports. In an effort to overcome this, the scientific literature was reviewed [Fordyce G. In: Fordyce G, editor. Bull fertility: selection and management in Australia. Eight Mile Plains, Australia: Australian Cattle Vets; 2002] and the needs of stakeholders were considered in preparing a manual, Evaluating and Reporting Bull Fertility [Entwistle KW, Fordyce G. Evaluating and reporting bull fertility. Eight Mile Plains, Australia: Australian Cattle Vets; 2003.] that outlined standards for assessing and reporting bull breeding soundness. A new recording and reporting system, called Bull Reporter, is based on standards from this manual and groups bull fertility traits into five summary categories: Scrotum, Physical, Crush-side Semen, Sperm Morphology, and Serving. The client will generally select which categories they wish to have included in the evaluation to suit their specific purposes. While there is adequate room for comments, the veterinarian is not required to make an overall judgment of whether the bull has normal capacity to sire calves under natural mating management, but ensures the standards for each selected category are met. Professional, standardised, easy-to-read reports are produced either electronically [Entwistle KW, Fordyce G. Evaluating and reporting bull fertility. Eight Mile Plains, Australia: Australian Cattle Vets; 2003.] or manually. A bull owner or their agent signs the certificate to affirm that bulls have not undergone procedures to rectify faults which may have otherwise caused them to fail the standards. An accreditation system for assessing sperm morphology was established because of its demonstrated relationship with pregnancy rates and because of the difficulties in achieving consistent and accurate assessments among laboratories. It is considered that Bull Reporter is applicable to beef and dairy bulls across all levels of management, genotypes and environments throughout Australia, with substantial potential for application elsewhere in the world. Crown Copyright (c) 2006 Published by Elsevier Inc. All rights reserved.

The effect of energy feedbacks on continental strength

The classical strength profile of continents(1,2) is derived from a quasi-static view of their rheological response to stress-one that does not consider dynamic interactions between brittle and ductile layers. Such interactions result in complexities of failure in the brittle-ductile transition and the need to couple energy to understand strain localization. Here we investigate continental deformation by solving the fully coupled energy, momentum and continuum equations. We show that this approach produces unexpected feedback processes, leading to a significantly weaker dynamic strength evolution. In our model, stress localization focused on the brittle-ductile transition leads to the spontaneous development of mid-crustal detachment faults immediately above the strongest crustal layer. We also find that an additional decoupling layer forms between the lower crust and mantle. Our results explain the development of decoupling layers that are observed to accommodate hundreds of kilometres of horizontal motions during continental deformation.

Zinc deposits and related mineralization of the Burketown mineral field, including the world-class Century Deposit, Northern Australia: Fluid inclusion and stable isotope evidence for basin fluid sources

The stratiform Century Zn-Pb deposit and the discordant Zn-Pb lode deposits of the Burketown mineral field, northern Australia, host ore and gangue minerals with primary fluid inclusions that have not been affected by the Isan orogeny, thus providing a unique opportunity to investigate the nature of the ore-forming brines. All of the deposits are hosted in shales and siltstones belonging to the Isa superbasin and comprise sphalerite, pyrite, carbonate, quartz, galena, minor chalcopyrite, and minor illite. According to Pb model ages, the main ore stage of mineralization at Century formed at I575 Ma, some 20 m.y. after deposition of the host shale sequence. Microthermometry on undeformed, primary fluid inclusions hosted in porous sphalerite shows that the Zn at Century was transported to the deposit by a homogeneous, Ca2+- and Na+-bearing brine with a salinity of 21.6 wt percent NaCl equiv. delta D-fluid of the fluid inclusion water ranges from -89 to -83 per mil, consistent with a basinal brine that evolved from meteoric water. Fluid inclusion homogenization temperatures range between 74 degrees and 125 degrees C, which are lower than the 120 degrees to 160 degrees C range calculated from vitrinite reflectance and illite crystallinity data from the deposit. This discrepancy indicates that mineralization likely formed at 50 to 85 Mpa, corresponding to a depth of 1,900 to 3,100 m. Transgressive galena-sphalerite veins that cut stratiform mineralization at Century and breccia-filled quartz-dolomite-sphalerite-galena veins in the discordant Zn-Pb lodes have Pb model ages between 1575 and 1485 Ma. Raman spectroscopy and microthermometry reveal that the primary fluid inclusions in these veins contain Ca2+, Na+. but they have lower salinities between 23 and 10 wt percent NaCl equiv and higher delta D-fluid values ranging from -89 to -61 per mil than fluid inclusions in porous sphalerite from Century. Fluid inclusion water from sphalerite in one of the lode deposits has delta O-18(fluid) values of 1.6 and 2.4 per mil, indistinguishable from delta O-18(fluid) values between -0.3 to +7.4 per mil calculated from the isotopic composition of coexisting quartz, dolomite, and illite. The trend toward lower salinities and higher delta D-fluid values relative to the earlier mineralizing fluids is attributed to mixing between the fluid that formed Century and a seawater-derived fluid from a different source. Based on seismic data from the Lawn Hill platform and paragenetic and geochemical results from the Leichhardt River fault trough to the south, diagenetic aquifers in the Underlying Calvert superbasin appear to have been the most likely sources for the fluids that formed Century and the discordant lode deposits. Paragenetically late sphalerite and calcite cut sphalerite, quartz, and dolomite in the lode deposits and contain Na+-dominated fluid inclusions with much lower salinities than their older counterparts. The isotopic composition of calcite also indicates delta O-18(fluid) from 3.3 to 10.7 per mil, which is larger than the range obtained from synmineralization minerals, supporting the idea that a unique fluid source was involved. The absolute timing of this event is unclear, but a plethora of Pb model, K-Ar, and Ar-40/Ar-39 ages between 1440 and 1300 Ma indicate that a significant volume of fluid was mobilized at this time. The deposition of the Roper superbasin from ca. 1492 +/- 4 Ma suggests that these late veins formed from fluids that may have been derived from aquifers in overlying sediments of the Roper superbasin. Clear, buck, and drusy quartz in veins unrelated to any form of Pb-Zn mineralization record the last major fluid event in the Burketown mineral field and form distinct outcrops and ridges in the district. Fluid inclusions in these veins indicate formation from a low-salinity, 300 degrees +/- 80 degrees C fluid. Temperatures approaching 300 degrees C recorded in organic matter adjacent to faults and at sequence boundaries correspond to K-Ar ages spanning 1300 to 1100 Ma, which coincides with regional hydrothermal activity in the northern Lawn Hill platform and the emplacement of the Lakeview Dolerite at the time of assemblage of the Rodinia supercontinent.

Volcanic setting of the Bajo de la Alumbrera porphyry Cu-Au deposit, Farallon Negro Volcanics, Northwest Argentina

The late Miocene Farallon Negro volcanics, comprising basaltic to rhyodacitic volcano-sedimentary rocks, host the Bajo de la Alumbrera porphyry copper-gold deposit in northwest Argentina. Early studies of the geology of the district have underpinned the general model for porphyry ore deposits where hydrothermal alteration and mineralization develop in and around porphyritic intrusions emplaced at shallow depths (2.5-3.5 km) into stratovolcanic assemblages. The Farallon Negro succession is dominated by thick sequences of volcano-sedimentary breccias, with lavas forming a minor component volumetrically. These volcaniclastic rocks conformably overlie crystalline basement-derived sedimentary rocks deposited in a developing foreland basin southeast of the Puna-Altiplano plateau. Within the Farallon Negro volcanics, volcanogenic accumulations evolved from early mafic to intermediate and silicic compositions. The younger and more silicic rocks are demonstrably coeval and comagmatic with the earliest group of mineralized porphyritic intrusions at Bajo de la Alumbrera. Our analysis of the volcanic stratigraphy and facies architecture of the Farallon Negro volcanics indicates that volcanic eruptions evolved from effusive to mixed effusive and explosive styles, as magma compositions changed to more intermediate and silicic compositions. Air early phase of mafic to intermediate voleanism was characterized by small synsedimentary intrusions with peperitic contacts, and lesser lava units scattered widely throughout the district, and interbedded with thick and extensive successions of coarse-grained sedimentary breccias. These sedimentary breccias formed from numerous debris- and hyperconcentrated flow events. A later phase of silicic volcanism included both effusive eruptions, forming several areally restricted lavas, and explosive eruptions, producing more widely dispersed (up to 5 kin) tuff units, some tip to 30-m thickness in proximal sections. Four key features of the volcanic stratigraphy suggest that the Farallon Negro volcanics need not simply record the construction of a large steep-sided polygenetic stratovolcano: (1) sheetlike, laterally continuous debris-flow and other coarse-grained sedimentary deposits are dominant, particularly in the lower sections; (2) mafic-intermediate composition lavas are volumetrically minor; (3) peperites are present throughout the sequence; and (4) fine-grained lacustrine sandstone-siltstone sequences occur in areas previously thought to be proximal to the summit region of the stratovolcano. Instead, the nature, distribution, and geometry of volcanic and volcaniclastic facies suggest that volcanism occurred as a relatively low relief, multiple-vent volcanic complex at the eastern edge of a broad, > 200-km-wide late Miocene volcanic belt and oil ail active foreland sedimentary basin to the Puna-Altiplano. Volcanism that occurred synchronously with the earliest stages of porphyry-related mineralization at Bajo de la Alumbrera apparently developed in an alluvial to ring plain setting that was distal to larger volcanic edifices.

Parallel 3D Simulation of a Fault Gouge using the Lattice Solid Model

Despite the insight gained from 2-D particle models, and given that the dynamics of crustal faults occur in 3-D space, the question remains, how do the 3-D fault gouge dynamics differ from those in 2-D? Traditionally, 2-D modeling has been preferred over 3-D simulations because of the computational cost of solving 3-D problems. However, modern high performance computing architectures, combined with a parallel implementation of the Lattice Solid Model (LSM), provide the opportunity to explore 3-D fault micro-mechanics and to advance understanding of effective constitutive relations of fault gouge layers. In this paper, macroscopic friction values from 2-D and 3-D LSM simulations, performed on an SGI Altix 3700 super-cluster, are compared. Two rectangular elastic blocks of bonded particles, with a rough fault plane and separated by a region of randomly sized non-bonded gouge particles, are sheared in opposite directions by normally-loaded driving plates. The results demonstrate that the gouge particles in the 3-D models undergo significant out-of-plane motion during shear. The 3-D models also exhibit a higher mean macroscopic friction than the 2-D models for varying values of interparticle friction. 2-D LSM gouge models have previously been shown to exhibit accelerating energy release in simulated earthquake cycles, supporting the Critical Point hypothesis. The 3-D models are shown to also display accelerating energy release, and good fits of power law time-to-failure functions to the cumulative energy release are obtained.

An automated failure mode and effect analysis based on high-level design specificication with behavior trees

Formal methods have significant benefits for developing safety critical systems, in that they allow for correctness proofs, model checking safety and liveness properties, deadlock checking, etc. However, formal methods do not scale very well and demand specialist skills, when developing real-world systems. For these reasons, development and analysis of large-scale safety critical systems will require effective integration of formal and informal methods. In this paper, we use such an integrative approach to automate Failure Modes and Effects Analysis (FMEA), a widely used system safety analysis technique, using a high-level graphical modelling notation (Behavior Trees) and model checking. We inject component failure modes into the Behavior Trees and translate the resulting Behavior Trees to SAL code. This enables us to model check if the system in the presence of these faults satisfies its safety properties, specified by temporal logic formulas. The benefit of this process is tool support that automates the tedious and error-prone aspects of FMEA.

Implementation of a triple modular redundant FPGA based safety critical system for reliable software execution

This paper describes the implementation of a TMR (Triple Modular Redundant) microprocessor system on a FPGA. The system exhibits true redundancy in that three instances of the same processor system (both software and hardware) are executed in parallel. The described system uses software to control external peripherals and a voter is used to output correct results. An error indication is asserted whenever two of the three outputs match or all three outputs disagree. The software has been implemented to conform to a particular safety critical coding guideline/standard which is popular in industry. The system was verified by injecting various faults into it.

Test strategies on functionally partitioned module-based programmable architecture for base-band processing

A specialised reconfigurable architecture for telecommunication base-band processing is augmented with testing resources. The routing network is linked via virtual wire hardware modules to reduce the area occupied by connecting buses. The number of switches within the routing matrices is also minimised, which increases throughput without sacrificing flexibility. The testing algorithm was developed to systematically search for faults in the processing modules and the flexible high-speed routing network within the architecture. The testing algorithm starts by scanning the externally addressable memory space and testing the master controller. The controller then tests every switch in the route-through switch matrix by making loops from the shared memory to each of the switches. The local switch matrix is also tested in the same way. Next the local memory is scanned. Finally, pre-defined test vectors are loaded into local memory to check the processing modules. This algorithm scans all possible paths within the interconnection network exhaustively and reports all faults. Strategies can be inserted to bypass minor faults

Stratigrafia fisica ed analisi di facies dei depositi torbiditici della Formazione delle Arenarie del Monte Cervarola tra la Val Secchia e la Val Scoltenna (Aquitaniano-Burdigaliano, Appennino Settentrionale)

The Cervarola Sandstones Formation (CSF), Aquitanian-Burdigalian in age, was deposited in an elongate, NW-stretched foredeep basin formed in front of the growing Northern Apennines orogenic wedge. The stratigraphic succession of the CSF, in the same way of other Apennine foredeep deposits, records the progressive closure of the basin due to the propagation of thrust fronts toward north-east, i.e. toward the outer and shallower foreland ramp. This process produce a complex foredeep characterized by synsedimentary structural highs and depocenters that can strongly influence the lateral and vertical turbidite facies distribution. Of consequence the main aim of this work is to describe and discuss this influence on the basis of a new high-resolution stratigraphic framework performed by measuring ten stratigraphic logs, for a total thickness of about 2000m, between the Secchia and Scoltenna Valleys (30km apart). In particular, the relationship between the turbidite sedimentation and the ongoing tectonic activity during the foredeep evolution has been describe through various stratigraphic cross sections oriented parallel and perpendicular to the main tectonic structures. On the basis of the high resolution physical stratigraphy of the studied succession, we propose a facies tract and an evolutionary model for the Cervarola Sandstones in the studied area. Thanks to these results and the analogies with others foredeep deposits of the northern Apennines, such as the Marnoso-arenacea Formation, the Cervarola basin has been interpreted as a highly confined foredeep controlled by an intense synsedimentary tectonic activity. The most important evidences supporting this hypothesis are: 1) the upward increase, in the studied stratigraphic succession (about 1000m thick), of sandstone/mudstone ratio, grain sizes and Ophiomorpha-type trace fossils testifying the high degree of flow deceleration related to the progressive closure and uplift of the foredeep. 2) the occurrence in the upper part of the stratigraphic succession of coarse-grained massive sandstones overlain by tractive structures such as megaripples and traction carpets passing downcurrent into fine-grained laminated contained-reflected beds. This facies tract is interpreted as related to deceleration and decoupling of bipartite flows with the deposition of the basal dense flows and bypass of the upper turbulent flows. 3) the widespread occurrence of contained reflected beds related to morphological obstacles created by tectonic structures parallel and perpendicular to the basin axis (see for example the Pievepelago line). 4) occurrence of intra-formational slumps, constituted by highly deformed portion of fine-grained succession, indicating a syn-sedimentary tectonic activity of the tectonic structures able to destabilize the margins of the basin. These types of deposits increase towards the upper part of the stratigraphic succession (see points 1 and 2) 5) the impressive lateral facies changes between intrabasinal topographic highs characterized by fine-grained and thin sandstone beds and marlstones and depocenters characterized by thick to very thick coarse-grained massive sandstones. 6) the common occurrence of amalgamation surfaces, flow impact structures and mud-draped scours related to sudden deceleration of the turbidite flows induced by the structurally-controlled confinement and morphological irregularities. In conclusion, the CSF has many analogies with the facies associations occurring in other tectonically-controlled foredeeps such as those of Marnoso-arenacea Formation (northern Italy) and Annot Sandstones (southern France) showing how thrust fronts and transversal structures moving towards the foreland, were able to produce a segmented foredeep that can strongly influence the turbidity current deposition.

The design of robust protocols for distributed real-time systems

Modern distributed control systems comprise of a set of processors which are interconnected using a suitable communication network. For use in real-time control environments, such systems must be deterministic and generate specified responses within critical timing constraints. Also, they should be sufficiently robust to survive predictable events such as communication or processor faults. This thesis considers the problem of coordinating and synchronizing a distributed real-time control system under normal and abnormal conditions. Distributed control systems need to periodically coordinate the actions of several autonomous sites. Often the type of coordination required is the all or nothing property of an atomic action. Atomic commit protocols have been used to achieve this atomicity in distributed database systems which are not subject to deadlines. This thesis addresses the problem of applying time constraints to atomic commit protocols so that decisions can be made within a deadline. A modified protocol is proposed which is suitable for real-time applications. The thesis also addresses the problem of ensuring that atomicity is provided even if processor or communication failures occur. Previous work has considered the design of atomic commit protocols for use in non time critical distributed database systems. However, in a distributed real-time control system a fault must not allow stringent timing constraints to be violated. This thesis proposes commit protocols using synchronous communications which can be made resilient to a single processor or communication failure and still satisfy deadlines. Previous formal models used to design commit protocols have had adequate state coverability but have omitted timing properties. They also assumed that sites communicated asynchronously and omitted the communications from the model. Timed Petri nets are used in this thesis to specify and design the proposed protocols which are analysed for consistency and timeliness. Also the communication system is mcxielled within the Petri net specifications so that communication failures can be included in the analysis. Analysis of the Timed Petri net and the associated reachability tree is used to show the proposed protocols always terminate consistently and satisfy timing constraints. Finally the applications of this work are described. Two different types of applications are considered, real-time databases and real-time control systems. It is shown that it may be advantageous to use synchronous communications in distributed database systems, especially if predictable response times are required. Emphasis is given to the application of the developed commit protocols to real-time control systems. Using the same analysis techniques as those used for the design of the protocols it can be shown that the overall system performs as expected both functionally and temporally.