Generalized Viscoplastic Modeling of Debris Flow

Various concepts have been proposed or used in the development of rheological models for debris flow. The earliest model developed by Bagnold was based on the concept of the “dispersive” pressure generated by grain collisions. Bagnold’s concept appears to be theoretically sound, but his empirical model has been found to be inconsistent with most theoretical models developed from non-Newtonian fluid mechanics. Although the generality of Bagnold’s model is still at issue, debris-flow modelers in Japan have generally accepted Takahashi’s formulas derived from Bagnold’s model. Some efforts have recently been made by theoreticians in non-Newtonian fluid mechanics to modify or improve Bagnold’s concept or model. A viable rheological model should consist both of a rate-independent part and a rate-dependent part. A generalized viscoplastic fluid (GVF) model that has both parts as well as two major rheological properties (i.e., the normal stress effect and soil yield criterion) is shown to be sufficiently accurate, yet practical, for general use in debris-flow modeling. In fact, Bagnold’s model is found to be only a particular case of the GVF model. Analytical solutions for (steady) uniform debris flows in wide channels are obtained from the GVF model based on Bagnold’s simplified assumption of constant grain concentration.

Birth-death processes with mass annihilation and state-dependent immigration

A birth-death process is subject to mass annihilation at rate β with subsequent mass immigration occurring into state j at rateα j . This structure enables the process to jump from one sector of state space to another one (via state 0) with transition rate independent of population size. First, we highlight the difficulties encountered when using standard techniques to construct both time-dependent and equilibrium probabilities. Then we show how to overcome such analytic difficulties by means of a tool developed in Chen and Renshaw (1990, 1993b); this approach is applicable to many processes whose underlying generator on E\{0} has known probability structure. Here we demonstrate the technique through application to the linear birth-death generator on which is superimposed an annihilation/immigration process.

Influence of the driving mechanism on the response of systems with athermal dynamics: The example of the random-field Ising model

We investigate the influence of the driving mechanism on the hysteretic response of systems with athermal dynamics. In the framework of local mean-field theory at finite temperature (but neglecting thermally activated processes), we compare the rate-independent hysteresis loops obtained in the random field Ising model when controlling either the external magnetic field H or the extensive magnetization M. Two distinct behaviors are observed, depending on disorder strength. At large disorder, the H-driven and M-driven protocols yield identical hysteresis loops in the thermodynamic limit. At low disorder, when the H-driven magnetization curve is discontinuous (due to the presence of a macroscopic avalanche), the M-driven loop is reentrant while the induced field exhibits strong intermittent fluctuations and is only weakly self-averaging. The relevance of these results to the experimental observations in ferromagnetic materials, shape memory alloys, and other disordered systems is discussed.

Variation of yield loci in finite element analysis by considering texture evolution for AA5042 aluminum sheets

 Yield function has various material parameters that describe how materials respond plastically in given conditions. However, a significant number of mechanical tests are required to identify the many material parameters for yield function. In this study, an effective method using crystal plasticity through a virtual experiment is introduced to develop the anisotropic yield function for AA5042. The crystal plasticity approach was used to predict the anisotropic response of the material in order to consider a number of stress or strain modes that would not otherwise be evident through mechanical testing. A rate-independent crystal plasticity model based on a smooth single crystal yield surface, which removes the innate ambiguity problem within the rate-independent model and Taylor model for polycrystalline deformation behavior were employed to predict the material’s response in the balanced biaxial stress, pure shear, and plane strain states to identify the parameters for the anisotropic yield function of AA5042.

Precipitation reaction in alpha-Cu-Al-Ag alloys

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

A generalized finite element method with global-local enrichment functions for confined plasticity problems

The main feature of partition of unity methods such as the generalized or extended finite element method is their ability of utilizing a priori knowledge about the solution of a problem in the form of enrichment functions. However, analytical derivation of enrichment functions with good approximation properties is mostly limited to two-dimensional linear problems. This paper presents a procedure to numerically generate proper enrichment functions for three-dimensional problems with confined plasticity where plastic evolution is gradual. This procedure involves the solution of boundary value problems around local regions exhibiting nonlinear behavior and the enrichment of the global solution space with the local solutions through the partition of unity method framework. This approach can produce accurate nonlinear solutions with a reduced computational cost compared to standard finite element methods since computationally intensive nonlinear iterations can be performed on coarse global meshes after the creation of enrichment functions properly describing localized nonlinear behavior. Several three-dimensional nonlinear problems based on the rate-independent J (2) plasticity theory with isotropic hardening are solved using the proposed procedure to demonstrate its robustness, accuracy and computational efficiency.

Exhumation rates in the Archean from pressure-ime paths: Examplefrom the Skjoldungen Orogen (SE Greenland)

The discussions on the orogenic evolution during Earth's history converge to the question of a different thermal structure in the Archean compared to the Phanerozoic and the applicability of the plate tectonic paradigm. However, geothermal structures are transient in orogens and are difficult to translate into large-scale tectonics and exhumation rates. Therefore, we propose depth–time data in the Archean Skjoldungen Orogen (SE Greenland, North Atlantic Craton) that allow for reconstruction of an exhumation rate independent of geothermal gradients. The resulting exhumation rate of ca. 0.4 km/Ma is similar to exhumation rates during erosion-controlled processes in modern orogens. These exhumation rates can only be established by erosion time constants similar to modern orogens. The occurrence of erosion-controlled exhumation is best explained by a stiff foreland promoting localized deformation in the orogen. Therefore, a switch from magmatic-dominated processes to localized deformation is proposed in the Skjoldungen Orogen area. This is supported by a change in magma composition and volume, from widespread granodiorite to localized alkaline intrusions. In addition, the involved metasedimentary rocks include detrital zircons of the only 50 Ma older foreland, which also correspond to erosion and tectonics as in modern orogens, i.e. flysh-type sediments. Relatively fast exhumation rates and the structural-magmatic evolution of the Neoarchean Skjoldungen Orogen thus indicate modern-style tectonic processes where stiff Mesoarchean continental crust forms a foreland to a collisional orogen instead of typical accretionary tectonics of weak island arc-like terranes in granite-greenstone terranes.

Dynamic friction and the origin of the complexity of earthquake sources.

We study a simple antiplane fault of finite length embedded in a homogeneous isotropic elastic solid to understand the origin of seismic source heterogeneity in the presence of nonlinear rate- and state-dependent friction. All the mechanical properties of the medium and friction are assumed homogeneous. Friction includes a characteristic length that is longer than the grid size so that our models have a well-defined continuum limit. Starting from a heterogeneous initial stress distribution, we apply a slowly increasing uniform stress load far from the fault and we simulate the seismicity for a few 1000 events. The style of seismicity produced by this model is determined by a control parameter associated with the degree of rate dependence of friction. For classical friction models with rate-independent friction, no complexity appears and seismicity is perfectly periodic. For weakly rate-dependent friction, large ruptures are still periodic, but small seismicity becomes increasingly nonstationary. When friction is highly rate-dependent, seismicity becomes nonperiodic and ruptures of all sizes occur inside the fault. Highly rate-dependent friction destabilizes the healing process producing premature healing of slip and partial stress drop. Partial stress drop produces large variations in the state of stress that in turn produce earthquakes of different sizes. Similar results have been found by other authors using the Burridge and Knopoff model. We conjecture that all models in which static stress drop is only a fraction of the dynamic stress drop produce stress heterogeneity.

Size independent fracture energy from fracture energy release rate in plain concrete beams

Size independent fracture energy and size effect on fracture energy are the key concerns for characterization of concrete fracture. Although there have been inconsistencies in results, a consensual fact is that the fracture energy from a large specimen is size independent. The fracture energy is proportional to the size of the fracture process zone (FPZ). FPZ size increases with size of the specimen, but the rate of increase of FPZ size decreases with increase in specimen size 1] implying that rate of increase of fracture energy decreases with increase in specimen size, more appropriately with increase in un-cracked ligament length. The ratio of fracture energy to the un-cracked ligament length almost becomes a constant at larger un-cracked ligament lengths. In the present study an attempt is made to obtain size independent fracture energy from fracture energy release rate. (C) 2012 Elsevier Ltd. All rights reserved.

A method for estimating the turbulent kinetic energy dissipation rate from a vertically pointing Doppler lidar, and independent evaluation from balloon-borne in situ measurements

A method of estimating dissipation rates from a vertically pointing Doppler lidar with high temporal and spatial resolution has been evaluated by comparison with independent measurements derived from a balloon-borne sonic anemometer. This method utilizes the variance of the mean Doppler velocity from a number of sequential samples and requires an estimate of the horizontal wind speed. The noise contribution to the variance can be estimated from the observed signal-to-noise ratio and removed where appropriate. The relative size of the noise variance to the observed variance provides a measure of the confidence in the retrieval. Comparison with in situ dissipation rates derived from the balloon-borne sonic anemometer reveal that this particular Doppler lidar is capable of retrieving dissipation rates over a range of at least three orders of magnitude. This method is most suitable for retrieval of dissipation rates within the convective well-mixed boundary layer where the scales of motion that the Doppler lidar probes remain well within the inertial subrange. Caution must be applied when estimating dissipation rates in more quiescent conditions. For the particular Doppler lidar described here, the selection of suitably short integration times will permit this method to be applicable in such situations but at the expense of accuracy in the Doppler velocity estimates. The two case studies presented here suggest that, with profiles every 4 s, reliable estimates of ϵ can be derived to within at least an order of magnitude throughout almost all of the lowest 2 km and, in the convective boundary layer, to within 50%. Increasing the integration time for individual profiles to 30 s can improve the accuracy substantially but potentially confines retrievals to within the convective boundary layer. Therefore, optimization of certain instrument parameters may be required for specific implementations.

Cortisol, blood pressure, and heart rate responses to food intake were independent of physical fitness levels in women

This research tested the hypothesis that women who had higher levels of physical fitness will have lower hypothalamo-pituitary-adrenal axis (cortisol) and sympatho-adrenal medullary system (blood pressure and heart rate) responses to food intake compared with women who had low levels of physical fitness. Lower fitness (n = 22; maximal oxygen consumption = 27.4 ± 1.0 mL∙kg(-1)·min(-1)) and higher fitness (n = 22; maximal oxygen consumption = 41.9 ± 1.6 mL∙kg(-1)·min(-1)) women (aged 30-50 years; in the follicular phase of the menstrual cycle) who participated in levels of physical activity that met (lower fitness = 2.7 ± 0.5 h/week) or considerably exceeded (higher fitness = 7.1 ± 1.4 h/week) physical activity guidelines made their own lunch using standardised ingredients at 1200 h. Concentrations of cortisol were measured in blood samples collected every 15 min from 1145-1400 h. Blood pressures and heart rate were also measured every 15 min between 1145 h and 1400 h. The meal consumed by the participants consisted of 20% protein, 61% carbohydrates, and 19% fat. There was a significant overall response to lunch in all of the parameters measured (time effect for all, p < 0.01). The cortisol response to lunch was not significantly different between the groups (time × treatment, p = 0.882). Overall, both groups showed the same pattern of cortisol secretion (treatment p = 0.839). Systolic blood pressure, diastolic blood pressure, mean arterial pressure, or heart rate responses (time × treatment, p = 0.726, 0.898, 0.713, and 0.620, respectively) were also similar between higher and lower fitness women. Results suggest that the physiological response to food intake in women is quite resistant to modification by elevated physical fitness levels.

Autonomic control of heart rate is virtually independent of temperature but seems related to the neuroanatomy of the efferent vagal supply to the heart in the bullfrog, Lithobathes catesbeianus

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Fishing directly selects on growth rate via behaviour: implications of growth-selection that is independent of size

Size-selective harvest of fish and crustacean populations has reduced stock numbers, and led to reduced growth rates and earlier maturation. In contrast to the focus on size-selective effects of harvest, here, we test the hypothesis that fishing may select on life-history traits (here, growth rate) via behaviour, even in the absence of size selection. If true, then traditional size-limits used to protect segments of a population cannot fully protect fast growers, because at any given size, fast-growers will be more vulnerable owing to bolder behaviour. We repeatedly measured individual behaviour and growth of 86 crayfish and found that fast-growing individuals were consistently bold and voracious over time, and were subsequently more likely to be harvested in single- and group-trapping trials. In addition, there was some indication that sex had independent effects on behaviour and trappability, whereby females tended to be less active, shyer, slower-growing and less likely to be harvested, but not all these effects were significant. This study represents, to our knowledge, the first across-individual support for this hypothesis, and suggests that behaviour is an important mechanism for fishing selectivity that could potentially lead to evolution of reduced intrinsic growth rates.

Raman microscopy of formamide-intercalated kaolinites treated by controlled-rate thermal analysis technology

Raman spectroscopy of formamide-intercalated kaolinites treated using controlled-rate thermal analysis technology (CRTA), allowing the separation of adsorbed formamide from intercalated formamide in formamide-intercalated kaolinites, is reported. The Raman spectra of the CRTA-treated formamide-intercalated kaolinites are significantly different from those of the intercalated kaolinites, which display a combination of both intercalated and adsorbed formamide. An intense band is observed at 3629 cm-1, attributed to the inner surface hydroxyls hydrogen bonded to the formamide. Broad bands are observed at 3600 and 3639 cm-1, assigned to the inner surface hydroxyls, which are hydrogen bonded to the adsorbed water molecules. The hydroxyl-stretching band of the inner hydroxyl is observed at 3621 cm-1 in the Raman spectra of the CRTA-treated formamide-intercalated kaolinites. The results of thermal analysis show that the amount of intercalated formamide between the kaolinite layers is independent of the presence of water. Significant differences are observed in the CO stretching region between the adsorbed and intercalated formamide.