Implementation of a coupled plastic damage distinct lattice spring model for dynamic crack propagation in geomaterials

This paper studies dynamic crack propagation by employing the distinct lattice spring model (DLSM) and 3‐dimensional (3D) printing technique. A damage‐plasticity model was developed and implemented in a 2D DLSM. Applicability of the damage‐plasticity DLSM was verified against analytical elastic solutions and experimental results for crack propagation. As a physical analogy, dynamic fracturing tests were conducted on 3D printed specimens using the split Hopkinson pressure bar. The dynamic stress intensity factors were recorded, and crack paths were captured by a high‐speed camera. A parametric study was conducted to find the influences of the parameters on cracking behaviors, including initial and peak fracture toughness, crack speed, and crack patterns. Finally, selection of parameters for the damage‐plasticity model was determined through the comparison of numerical predictions and the experimentally observed cracking features.     

Assessment of alternative adsorption models and global sensitivity analysis to characterize hexavalent chromium loss from soil to surface runoff

We investigate our ability to assess transfer of hexavalent chromium, Cr(VI), from the soil to surface runoff by considering the effect of coupling diverse adsorption models with a two‐layer solute transfer model. Our analyses are grounded on a set of two experiments associated with soils characterized by diverse particle size distributions. Our study is motivated by the observation that Cr(VI) is receiving much attention for the assessment of environmental risks due to its high solubility, mobility, and toxicological significance. Adsorption of Cr(VI) is considered to be at equilibrium in the mixing layer under our experimental conditions. Four adsorption models, that is, the Langmuir, Freundlich, Temkin, and linear models, constitute our set of alternative (competing) mathematical formulations. Experimental results reveal that the soil samples characterized by the finest grain sizes are associated with the highest release of Cr(VI) to runoff. We compare the relative abilities of the four models to interpret experimental results through maximum likelihood model calibration and four model identification criteria (i.e., the Akaike information criteria [AIC and AIC_C] and the Bayesian and Kashyap information criteria). Our study results enable us to rank the tested models on the basis of a set of posterior weights assigned to each of them. A classical variance‐based global sensitivity analysis is then performed to assess the relative importance of the uncertain parameters associated with each of the models considered, within subregions of the parameter space. In this context, the modelling strategy resulting from coupling the Langmuir isotherm with a two‐layer solute transfer model is then evaluated as the most skilful for the overall interpretation of both sets of experiments. Our results document that (a) the depth of the mixing layer is the most influential factor for all models tested, with the exception of the Freundlich isotherm, and (b) the total sensitivity of the adsorption parameters varies in time, with a trend to increase as time progresses for all of the models. These results suggest that adsorption has a significant effect on the uncertainty associated with the release of Cr(VI) from the soil to the surface runoff component.     

Area usage estimation and spatiotemporal variability in distribution patterns of southern right whales,Eubalaena australis,of southern Brazil

Southern right whales—Eubalaena australis (Desmoulins, 1822)—migrate seasonally from high‐latitude feeding grounds to coastal breeding and calving grounds at lower latitudes such as the southern coast of Brazil. Understanding how these whales are distributed along the coast is important for monitoring their postwhaling recovery and defining management strategies. In this study, we applied Kernel density estimators to aerial survey data to determine main occurrence and concentration areas of right whales in southern Brazil and investigate inter‐ and intra‐annual distribution patterns between 2003 and 2012. Our results show considerable variation in area usage within and among years, and changes in the general distribution pattern of right whales in the last years of the study. Intra‐annually, higher concentration area tended to expand from July to September and decrease in November. Some areas stood out as high‐density areas for right whales: Ribanceira/Ibiraquera, Itapirubá Sul/Sol, and from Arroio to Gaivota. Some evidences also suggest preferential areas for mother–calf pairs. The higher concentration area of right whales in southern Brazil was estimated at 52,541 km² and the occurrence area was 682.69 km², which is the whole study area. As right whale distribution in the region is likely expanding due to this population's current recovery, our study provides essential information for management plan of the Right Whale Environmental Protection Area.     

Small hypervelocity particles captured in aerogel collectors: Location,extraction, handling and storage

Abstract— It has now been about a decade since the first demonstrations that hypervelocity particles could be captured, partially intact, in aerogel collectors. But the initial promise of a bonanza of partially‐intact extraterrestrial particles, collected in space, has yet to materialize. One of the difficulties that investigators have encountered is that the location, extraction, handling and analysis of very small (10 μm and less) grains, which constitute the vast majority of the captured particles, is challenging and burdensome. Furthermore, current extraction techniques tend to be destructive over large areas of the collectors. Here we describe our efforts to alleviate some of these difficulties. We have learned how to rapidly and efficiently locate captured particles in aerogel collectors, using an automated microscopic scanning system originally developed for experimental nuclear astrophysics. We have learned how to precisely excavate small access tunnels and trenches using an automated micromanipulator and glass microneedles as tools. These excavations are only destructive to the collector in a very small area—this feature may be particularly important for excavations in the precious Stardust collectors. Using actuatable silicon microtweezers, we have learned how to extract and store “naked” particles—essentially free of aerogel—as small as 3 μm in size. We have also developed a technique for extracting particles, along with their terminal tracks, still embedded in small cubical aerogel blocks. We have developed a novel method for storing very small particles in etched nuclear tracks. We have applied these techniques to the extraction and storage of grains captured in aerogel collectors (Particle Impact Experiment, Orbital Debris Collector Experiment, Comet‐99) in low Earth orbit.     

Consolidation of double‐layered ground with vertical drains

Vertical drains are usually installed in subsoil consisting of several layers. Due to the complex nature of the problem, over the past decades, the consolidation properties of multi‐layered ground with vertical drains have been analysed mainly by numerical methods. An analytical solution for consolidation of double‐layered ground with vertical drains under quasi‐equal strain condition is presented in this paper. The main steps for the computation procedure are listed. The convergence of the series solution is discussed. The comparisons between the results obtained by the present analytical method and the existing numerical solutions are described by figures. The orthogonal relation for the system of double‐layered ground with vertical drains is proven. Finally, some consolidation properties of double‐layered ground with vertical drains are analysed. Copyright © 2001 John Wiley & Sons, Ltd.     

Non‐equilibrium concepts lead to a unified explanation of the formation of chondrules and chondrites

Abstract— Calculations of the formation of seven types of chondrules in Semarkona from a gas of solar composition were performed with the FACT computer program to predict the chemistries of oxides (including silicates), developed by the authors and their colleagues. The constrained equilibrium theory was used in the calculations with two nucleation constraints suggested by nucleation theory. The first constraint was the blocking of Fe and other metal gaseous atoms from condensing to form solids or liquids because of very high surface free energies and high surface tensions of the solid and liquid metals, respectively. The second constraint was the blocking of the condensation of solids and the formation of metastable liquid oxides (including silicates) well below their liquidus temperatures. Our laboratory experiments suggested subcooling of type IIA chondrule compositions of 400 degrees or more below the liquidus temperature. The blocking of iron leads to a supersaturation of Fe atoms, so that the partial pressure of Fe (p_Fe) is larger than the partial pressure at equilibrium (p_Fe(eq)). The supersaturation ratio S = p_Fe/p_Fe(eq) becomes larger than 1 and increases rapidly with a decrease in temperature. This drives the reaction Fe + H₂O ? H₂ + FeO to the right. With S = 100, the activity of FeO in the liquid droplet is 100 times as large as the value at equilibrium. The FeO activities are a function of temperature and provide relative average temperatures of the crystallization of chondrules. Our calculations for the LL3.0 chondrite Semarkona and our study of some non‐equilibrium effects lead to accurate representations of the compositions of chondrules of types IA, IAB, IB, IIA, IIAB, IIB, and CC. Our concepts readily explain both the variety of FeO concentrations in the different chondrule types and the entire process of chondrule formation. Our theory is unified and could possibly explain the formation of chondrules in all chondritic meteorites as well as provide a simple explanation for the complex chemistries of chondrites, and especially for type 3 chondrites.     

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.     

格网DEM地形模拟的形态保真度研究

分析了现有格网DEM地形模拟的失真现象,研究DEM地形模拟失真的根源.提出了DEM地形形态保真度的概念,探讨了建设高保真DEM必须解决的问题.     

Observational evidence for solar coronal decimetric radio pulsations with very short periods

Using the decimetric (700–1500 MHz) radio spectrometer and the synchronous observational system with high temporal resolution at four frequencies (1420, 2130, 2840 and 4260 MHz) of Yunnan Observatory, two rare events were observed on 2001 June 24 and 1990 July 30. The former was a small radio burst exhibiting pulsations with short periods (about 29, 40 and 100 ms) in the impulsive phase. The latter was a large radio burst, which at 2840 MHz produced radio pulsations with period of about 30 ms. This paper focuses on pulsations with very short periods in the range of 29–40 ms. The mechanism of generation of such pulsations may be modulation of radio radiation by the periodic trains of whistler packets originating in unstable regions of the corona. Alternatively, these pulsations can be attributed to wave-wave non-linear interactions of electrostatic upper hybrid waves driven by beams of precipitating electrons in flaring loops.     

On the orbits and masses of the satellites of the Pluto-Charon system

Two small satellites of Pluto, S/2005 P1 (hereafter P1) and S/2005 P2 (hereafter P2), have recently been discovered outside the orbit of Charon, and their orbits are nearly circular and nearly coplanar with that of Charon. Because the mass ratio of Charon-Pluto is ∼0.1, the orbits of P2 and P1 are significantly non-Keplerian even if P2 and P1 have negligible masses. We present an analytic theory, with P2 and P1 treated as test particles, which shows that the motion can be represented by the superposition of the circular motion of a guiding center, the forced oscillations due to the non-axisymmetric components of the potential rotating at the mean motion of Pluto-Charon, the epicyclic motion, and the vertical motion. The analytic theory shows that the azimuthal periods of P2 and P1 are shorter than the Keplerian orbital periods, and this deviation from Kepler's third law is already detected in the unperturbed Keplerian fit of Buie and coworkers. In this analytic theory, the periapse and ascending node of each of the small satellites precess at nearly equal rates in opposite directions. From direct numerical orbit integrations, we show the increasing influence of the proximity of P2 and P1 to the 3:2 mean-motion commensurability on their orbital motion as their masses increase within the ranges allowed by the albedo uncertainties. If the geometric albedos of P2 and P1 are high and of order of that of Charon, the masses of P2 and P1 are sufficiently low that their orbits are well described by the analytic theory. The variation in the orbital radius of P2 due to the forced oscillations is comparable in magnitude to that due to the best-fit Keplerian eccentricity, and there is at present no evidence that P2 has any significant epicyclic eccentricity. However, the orbit of P1 has a significant epicyclic eccentricity, and the prograde precession of its longitude of periapse with a period of 5300 days should be easily detectable. If the albedos of P2 and P1 are as low as that of comets, the large inferred masses induce significant short-term variations in the epicyclic eccentricities and/or periapse longitudes on the 400-500-day timescales due to the proximity to the 3:2 commensurability. In fact, for the maximum inferred masses, P2 and P1 may be in the 3:2 mean-motion resonance, with the resonance variable involving the periapse longitude of P1 librating. Observations that sample the orbits of P2 and P1 well on the 400-500-day timescales should provide strong constraints on the masses of P2 and P1 in the near future.