An element-based displacement preconditioner for linear elasticity problems

Finite element analysis of problems in structural and geotechnical engineering results in linear systems where the unknowns are displacements and rotations at nodes. Although the solution of these systems can be carried out using either direct or iterative methods, in practice the matrices involved are usually very large and sparse (particularly for 3D problems) so an iterative approach is often advantageous in terms of both computational time and memory requirements. This memory saving can be further enhanced if the method used does not require assembly of the full coefficient matrix during the solution procedure. One disadvantage of iterative methods is the need to apply preconditioning to improve convergence. In this paper, we review a range of established element-based preconditioning methods for linear elastic problems and compare their performance with a new method based on preconditioning with element displacement components. This new method appears to offer a significant improvement in performance.     

CYLINDRICAL AND CONICAL FOLD GEOMETRIES IN THE CANTARELL STRUCTURE, SOUTHERN GULF OF MEXICO: IMPLICATIONS FOR HYDROCARBON EXPLORATION

The NW‐SE trending Cantarell structure in the Gulf of Campeche hosts the largest oilfield in Mexico. The oil occurs predominantly in latest Cretaceous – earliest Tertiary breccias with subsidiary reserves in Upper Jurassic (Oxfordian and Kimmeridgian) and Lower Cretaceous oolitic and partially dolomitized limestones, dolomites and shaly limestones. Cantarell has been interpreted both as a fold‐and‐thrust zone and as a dextral transpressional structure. Analysis of structure contours at 100m intervals, on the tops of the Tertiary breccia and the Kimmeridgian (Upper Jurassic) dolomite, indicates that the structure is an upright cylindrical fold with gently plunging conical terminations; there is also a conical portion in the central part of the structure. The axes of the central, NW and SE cones are subvertical. This geometry indicates that the two fold terminations and the central cone are aprons rather than points, with the NW and central cone axes intersecting the cylindrical fold axis at the point where the geometry switches from conical to cylindrical. The apical angle (i.e. the angle between the fold and cone axes) varies as follows: (i) in the NW cone, it is ～70° in the breccia and ～76° in the Kimmeridgian dolomite; (ii) in the central cone, it is ～77° in the breccia and ～73° in the Kimmeridgian dolomite; and (iii) in the SE cone, it is ～64° in the breccia and ～57° in the Kimmeridgian dolomite. This indicates that whereas the fold opens with depth in the NW cone, it tightens with depth in the central and SE cones. Assuming a parallel fold geometry, these apical angles indicate an increase in volume in the NW cone (i.e. larger hydrocarbon reservoirs), compared to the central and SE cones. Theoretical considerations indicate that the curvature increases dramatically towards the point of the cone. In the case of the Cantarell structure, the apices of the cones are located at the conical‐cylindrical fold junctions, where the highest curvature may have resulted in a higher degree of fracturing. The coincidence of maximum curvature and the intersection of the conical and cylindrical fold axes in the fold culminations with porous and permeable reservoir rocks may have made these locations favourable for the accumulation of hydrocarbons.     

Identification of the adverse effect of nitrate on the phosphate release rate and improvement of EBPR process models.

The adverse effect of nitrate on the phosphate release rate in the anaerobic phase was observed and was hardly explainable with conventional EBPR process models. Four possible mechanisms were proposed including substrate competition, reduced fermentation, parallel reaction and sequential reaction. Batch experiments were designed and conducted to identify the dominant mechanism. Results showed that the sequential reaction was the only possible mechanism where only denitrification occurred if any nitrate existed in the anaerobic phase. Then the phosphate release following after the nitrate was completely removed. Nitrate inhibition effect was added into the PHA storage rate to incorporate the sequential reaction in the conventional ASM3 plus EAWAG bio-P module (ASM3 + P). Nitrate inhibition coefficient, K(I,NO,PAO) was found to be as low as 0.05 mg/L. This correlated well with experimental observation where no also meant that the anaerobic compartment of a continuous flow reactor could be seriously affected by the residual nitrate contained in the sludge recycle flow. This phenomenon caused overestimation of the phosphate uptake rate and consequently underestimation of PO4(3-) -P concentration. This problem was resolved by incorporation of a nitrate inhibition term in the ASM3 + P for more accurate simulation of the EBPR process.     

Modeling and verification of the double frequency effect using a MEMS device

The derivation, verification, and implication of the nonlinear dynamic and frequency response of electrostatic actuator due to the double frequency effect (DFE) were reported in this study. In particular, an extra mode called half mode was observed and measured in various studies. However, a complete in-depth discussion of the effect was not reported in the past. In the present study, a second-order dynamical equation was adapted firstly to model the dynamic and frequency response of electrostatic actuator where typical harmonic input signal with a dc bias was used. Secondly, by solving the equation, complex waveform in dynamic response and an extra half mode in frequency response due to the double frequency effect can be observed and discussed. To verify the simulated result, an electrostatic driving device was fabricated using PolyMUMPS^© process. Note that in frequency response, when dc bias is equal to the amplitude of ac signal, simulated and experimental results indicated that the amplitude of half mode was one-fourth of first mode.     

Design of a retention tank: comparison of stormwater quality models with various levels of complexity.

Stormwater quality simulation models are useful tools for the design and management of sewer systems. Modelling results are highly sensitive to experimental data used for calibration. This sensitivity is examined for three modelling approaches of various complexities (site mean concentration approach, event mean concentration approach and build-up, washoff and transport modelling approach) applied to a typical case study (design of a dry detention tank), accounting for the variability of calibration data and their effect on simulation results. Calibrated models with different calibration data sets were used to simulate 3 years of rainfall with different retention tank specific volumes. Annual pollutant load interception efficiencies were determined. Simulations results revealed (i) that there is no advantage in using the EMC model compared to the SMC model and (ii) that the BWT model resulted in higher design ratios than those given by the SMC/hydraulic approach. For both EMC and BWT models, using an increasing number n of events for calibration leads to narrower confidence intervals for the design ratios. It is crucial for design ratios to account for successive storm events in chronological order and to account for the maximum allowable flow to be transferred to the downstream WWTP.     

An inventory analysis of rural pollution loads in China.

Rural pollution has attracted increasing attention over the past decade for its important consequences on surface and groundwater quality. This is particularly so in China due to the wide and effective control of industrial and municipal wastewater discharges over the past decade. Based on the discussion of rural activities in China, this paper develops a new method of inventory analysis through the definition of elementary unit (EU). This inventory analysis approach is then applied to estimate the total rural pollution loads of COD, TN and TP, as well as their spatial distributions in China. Further analysis is also conducted to investigate the correlations between rural pollution loads and agricultural production outputs, as well as population. It is found that China's agriculture is developing largely at the cost of environment quality, and a high population, as well as the increasing demand for agricultural products, is one of the dominant driving forces. The constructed EKC, which describes the relationship between rural pollution loads and agricultural income, also reveals that current agricultural development in China is mostly far from de-locked from the environmental problem.     

PETROLEUM MIGRATION, FAULTS AND OVERPRESSURE, PART I: CALIBRATING BASIN MODELLING USING PETROLEUM IN TRAPS — A REVIEW

This paper considers the principles of deciphering basin-scale hydrocarbon migration patterns using the geochemical information which is present in trapped petroleum. Petroleum accumulations in subsiding basins can be thought of as "data archives" within which stored information can help us to understand aspects of hydrocarbon formation and migration. This information can impart a time-resolved picture of hydrocarbon migration in a basin in response to processes associated with progressive burial, particularly in the context of the occurrence and periodic activity of faults.
This review, which includes a series of tentative models of migration-related processes in the extensional Halten Terrace area, offshore mid-Norway, illustrates how we can use information from the migrating mobile hydrocarbon phase to improve our knowledge of the static geological system. Of particular importance is the role of sub-seismic heterogeneities and faults in controlling migration processes. We focus on how the secondary migration process can be enhanced in a multi-source rock basin such as the Halten Terrace, thereby increasing prospectivity.     

The Applied Mathematics and Computer Science Schism

Isolating applied mathematics from computer science harms the profession and those who depend on it.     

The effect of grain size, strain rate, and temperature on the mechanical behavior of commercial purity aluminum

Commercial purity aluminum AA1050 was subjected to equal channel angular extrusion (ECAE) that resulted in an ultrafine-grained (UFG) microstructure with an as-received grain size of 0.35 μm. This UFG material was then annealed to obtain microstructures with grain sizes ranging from 0.47 to 20 μm. Specimens were compressed at quasi-static, intermediate, and dynamic strain rates at temperatures of 77 and 298 K. The mechanical properties were found to vary significantly with grain size, strain rate, and temperature. Yield stress was found to increase with decreasing grain size, decreasing temperature, and increasing strain rate. The work hardening rate was seen to increase with increasing grain size, decreasing temperature, and increasing strain rate. The influence of strain rate and temperature is most significant in the smallest grain size specimens. The rate of work hardening is also influenced by strain rate, temperature, and grain size with negative rates of work hardening observed at 298 K and quasi-static strain rates in the smallest grain sizes and increasing rates of work hardening with increasing loading rate and grain size. Work hardening behavior is correlated with the substructural evolution of these specimens.     

Sensor guided handling and assembly of active micro-systems

Positioning accuracies in the range of a few micrometers and below are necessary for the assembly of active micro-systems. In order to reach these accuracies, an assembly system for sensor guided micro-assembly is developed in the Collaborative Research Centre 516 “Design and manufacturing of active micro-systems”. The combination of a parallel robot with an integrated 3D vision sensor and a micro-gripper enables to reach relative positioning accuracies below 1 μm.