Fast iterative solution of large undrained soil‐structure interaction problems

In view of rapid developments in iterative solvers, it is timely to re‐examine the merits of using mixed formulation for incompressible problems. This paper presents extensive numerical studies to compare the accuracy of undrained solutions resulting from the standard displacement formulation with a penalty term and the two‐field mixed formulation. The standard displacement and two‐field mixed formulations are solved using both direct and iterative approaches to assess if it is cost‐effective to achieve more accurate solutions. Numerical studies of a simple footing problem show that the mixed formulation is able to solve the incompressible problem ‘exactly’, does not create pressure and stress instabilities, and obviate the need for an ad hoc penalty number. In addition, for large‐scale problems where it is not possible to perform direct solutions entirely within available random access memory, it turns out that the larger system of equations from mixed formulation also can be solved much more efficiently than the smaller system of equations arising from standard formulation by using the symmetric quasi‐minimal residual (SQMR) method with the generalized Jacobi (GJ) preconditioner. Iterative solution by SQMR with GJ preconditioning also is more elegant, faster, and more accurate than the popular Uzawa method. Copyright © 2003 John Wiley & Sons, Ltd.     

Screening of polymers on selected Hawaii soils for erosion reduction and particle settling

In recent years, high‐molecular‐weight anionic polyacrylamides (PAMs) have been tested on a variety of soils, primarily in temperate climates. However, little information is available regarding the effectiveness of PAM for preventing soil loss through runoff in tropical settings. Screening tests were performed using three negatively charged PAMs and one positively charged PAM on five Hawaii soils (two Oxisols, one Vertisol, and two Aridisols) to determine erosion loss, sediment settling, and aggregate stability. A laboratory‐scale rainfall simulator was used to apply erosive rainfall at intensities from 5 to 8·5 cm h^?1 at various PAM doses applied in both dry and solution forms. Soil detachment due to splash and runoff, as well as the runoff and percolate water volumes, were measured for initial and successive storms. The impact of PAM on particle settling and aggregate stability was also evaluated for selected soil‐treatment combinations. Among the PAMs, Superfloc A‐836 was most effective, and significantly reduced runoff and splash sediment loss for the Wahiawa Oxisol and Pakini Andisol at rates varying between 10 and 50 kg ha^?1. Reduced runoff and splash sediment loss were also noted for PAM Aerotil‐D when applied in solution form to the Wahiawa Oxisol. Significant reductions in soil loss were not noted for either the Lualualei Vertisol or the Holomua Oxisol. It is believed that the high montmorillonite content of the Lualualei Vertisol and the low cation‐exchange capacity of the Holomua Oxisol diminished the effectiveness of the various PAMs tested. The polymers were also found to enhance sediment settling of all soils and helped improve their aggregate stability. This screening study shows the potential use of PAM for tropical soils for applications such as infiltration enhancement, runoff reduction, and enhanced sedimentation of detention ponds. Copyright © 2005 John Wiley & Sons, Ltd.     

Analytical solutions of one‐dimensional infiltration before and after ponding

This paper is a continuation of previous research, which obtained a convenient solution for arbitrary surface fluxes before ponding. By means of Fourier Transformation this has been extended to derive analytical solutions of a linearized Richards' equation for arbitrary input fluxes after surface saturation. Copyright © 2003 John Wiley & Sons, Ltd.     

Geochemical correlations between effusive and explosive silicic volcanics in the Saraykent region (Yozgat), central Anatolia,Turkey

Two main volcanic events are distinguished between Saraykent and Akçakışla in the Yozgat province of central Anatolia: (1) early Late Cretaceous–Palaeocene effusive activity, that produced a sequence of intermediate to felsic ‘basal lavas’; and (2) marginally later Palaeocene explosive activity that formed a series of covering ignimbrite flows. Due to their close temporal and spatial relation, geochemical comparisons were made between the silicic members of the lavas and ignimbrites, to identify chemical groups and their relative petrogenesis. The basal lavas range from calc‐alkaline basaltic andesites to dominant rhyolites. Based on trace element correlations three main geochemical groups were identified: the Akçakışla rhyolites (present as domes); Akçakışla rhyodacites‐dacites (lava flows); and Ozan‐Saraykent rhyolites (lava flows). Large‐ion lithophile elements have been mobile in all the groups, but mainly in the Akçakışla rhyolites. Rare earth element (REE) patterns show marked similarity between the Ozan and Saraykent basal lavas. The Akçakışla dome rhyolites are more fractionated with lower La_N/Yb_N ratios (c.10), whereas the Akçakışla basal lavas have much higher La_N/Yb_N ratios (c.30). The chemical coherence and petrographic similarities between the Saraykent and Ozan lavas suggest a single suite related via fractionation. Three geochemical groups were also established for the ignimbrites: Saraykent ignimbrite; Bağlıca ignimbrite‐Toklu‐Kızıldağ crystal tuffs; and Keklikpınar ignimbrite. The ignimbrites, like the basal lavas, display a pronounced depletion in Ba on ORG‐normalized plots. Relative to the basal lavas, chondrite‐normalized patterns for the ignimbrites are different in displaying negative Eu anomalies that indicate feldspar fractionation. The lack of geochemical overlap or coherence between any of the lava and ignimbrite groups suggests that they represent distinct eruptive events and are not related in any simple volcanic development and cogenetic sense. Two geochemical features are common to all the volcanic rock groups: (1) the presence of a Nb‐Ta anomaly, which is generally accepted as a crustal signature; and (2) the relatively low Y abundances which appear characteristic for the region as a whole. These fundamental features of the local silicic volcanism largely reflect source composition and effects. Copyright © 2002 John Wiley & Sons, Ltd.     

Effect of hysteretic smoothness on inelastic response spectra with constant‐ductility

This paper presents an efficient methodology for computing constant‐ductility inelastic response spectra. The computation of constant‐ductility spectra involves numerical root‐finding algorithms to find the strongest structure providing a desired ductility response. Smooth inelastic structural behavior is modeled using a first‐order nonlinear differential equation and the transient structural response is solved using an implicit algorithm requiring Newton iterations at each time step. For structural models with smooth hysteretic behavior (not piece‐wise linear), a simple root‐finding method involving a combination of hyperbolic fits, linear interpolation, and Newton's method converges upon the highest strength (conservative) solution with a small number of iterations. The effect of the hysteretic smoothness on the occurrence of multiple roots is examined for two near‐fault and two far‐fault earthquake records, and for two measures of ductility and for normalized hysteretic energy. The results indicate how the smoothness of the hysteretic behavior affects ductility demand and constant‐ductility response spectra. Copyright © 2010 John Wiley & Sons, Ltd.     

Finite element analysis of the effect of corrosion on the behavior of reinforced earth walls

A simple model for the corrosion‐induced loss of stiffness and strength of the steel strips of earth‐reinforced walls was introduced in a finite element simulation of the long‐term behavior of the wall, in which the backfill‐strips interactions are taken into account by means of a generalized homogenization procedure (called a multiphase model). The results show an initial phase of slow displacements induced by the loss of stiffness, followed after a few decades by a steep acceleration of the displacements, leading to wall failure. The influences of the parameter controlling corrosion, the backfill cohesion and the heterogeneity of the corrosion process are discussed. Results are used to discuss a strategy for reinforced earth wall surveillance. Copyright © 2011 John Wiley & Sons, Ltd.     

A mass ratio limit for primordial black holes

We augment our scenario for the formation of astronomical objects from macroscopic superstrings by the assumption that the central matter keeps its identity in the fragmentation. From the condition that the angular momentum per mass squared of this matter should be less than the Kerr limit G/c, we obtain upper limits for the ratio of the mass of central black holes M(BH) to the mass M of the host object. This limit is M(BH)/M ≈ 0.001, and, expressed in observed quantities, approximately M(BH)/M ≈ σ²/(v · c) where σ is the r.m.s. velocity, v the rotational velocity and c the velocity of light. The valuesM(BH) agree with the observed behaviour both in order of magnitude and in the variation with velocity dispersion. (© 2003 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Molecular study of insoluble organic matter in Kainsaz CO3 carbonaceous chondrite: Comparison with CI and CM IOM

Abstract— Kainsaz CO3 insoluble organic matter (IOM) was studied using Curie point pyrolysis, electronic paramagnetic resonance (EPR), and high‐resolution transmission electron microscopy (HRTEM) to determine the effect of thermal metamorphism on molecular chondritic fingerprints. Pyrolysis released a very low amount of products that consist of one‐ and two‐ring aromatic units with methyl, dimethyl, and ethyl substituents. Moreover, Kainsaz IOM contains two orders of magnitude fewer radicals than Orgueil, Murchison, and Tagish Lake IOM. In addition, no diradicaloids were found in Kainsaz, although they are thought to constitute a specific signature for weakly organized extraterrestrial organic compounds in contrast to terrestrial ones. HRTEM reveals a very heterogeneous structure, with microporous disordered carbon, mesoporous graphitic carbons and graphite. Graphitization likely occurs and explains the differences between Kainsaz and CI or CM IOM. Heating stress experienced by Kainsaz IOM, on the parent body and/or prior its accretion, is likely responsible for the differences in molecular and structural organizations compared with those of CI and CM IOM.     

New Reference Materials and Assessment of Matrix Effects for SIMS Measurements of Oxygen Isotopes in Garnet

Accurate ion microprobe analysis of oxygen isotope ratios in garnet requires appropriate reference materials to correct for instrumental mass fractionation that partly depends on the garnet chemistry (matrix effect). The matrix effect correlated with grossular, spessartine and andradite components was characterised for the Cameca IMS 1280HR at the SwissSIMS laboratory based on sixteen reference garnet samples. The correlations fit a second‐degree polynomial with maximum bias of ca. 4‰, 2‰ and 8‰, respectively. While the grossular composition range 0–25% is adequately covered by available reference materials, there is a paucity of them for intermediate compositions. We characterise three new garnet reference materials GRS2, GRS‐JH2 and CAP02 with a grossular content of 88.3 ± 1.2% (2s), 83.3 ± 0.8% and 32.5 ± 3.0%, respectively. Their micro scale homogeneity in oxygen isotope composition was evaluated by multiple SIMS sessions. The reference δ¹⁸O value was determined by CO₂ laser fluorination (δ¹⁸O_LF). GRS2 has δ¹⁸O_LF = 8.01 ± 0.10‰ (2s) and repeatability within each SIMS session of 0.30–0.60‰ (2s), GRS‐JH2 has δ¹⁸O_LF = 18.70 ± 0.08‰ and repeatability of 0.24–0.42‰ and CAP02 has δ¹⁸O_LF = 4.64 ± 0.16‰ and repeatability of 0.40–0.46‰.