On the possibility of elastic strain localisation in a fault

The phenomenon of strain localisation is often observed in shear deformation of particulate materials, e.g., fault gouge. This phenomenon is usually attributed to special types of plastic behaviour of the material (e.g., strain softening or mismatch between dilatancy and pressure sensitivity or both). Observations of strain localisation in situ or in experiments are usually based on displacement measurements and subsequent computation of the displacement gradient. While in conventional continua the symmetric part of the displacement gradient is equal to the strain, it is no longer the case in the more realistic descriptions within the framework of generalised continua. In such models the rotations of the gouge particles are considered as independent degrees of freedom the values of which usually differ from the rotation of an infinitesimal volume element of the continuum, the latter being described for infinitesimal deformations by the non-symmetric part of the displacement gradient. As a model for gouge material we propose a continuum description for an assembly of spherical particles of equal radius in which the particle rotation is treated as an independent degree of freedom. Based on this model we consider simple shear deformations of the fault gouge. We show that there exist values of the model parameters for which the displacement gradient exhibits a pronounced localisation at the mid-layers of the fault, even in the absence of inelasticity. Inelastic effects are neglected in order to highlight the role of the independent rotations and the associated additional parameters. The localisation-like behaviour occurs if (a) the particle rotations on the boundary of the shear layer are constrained (this type of boundary condition does not exist in a standard continuum) and (b) the contact moment-or bending stiffness is much smaller than the product of the effective shear modulus of the granulate and the square of the width of the gouge layer. It should be noted however that the virtual work functional is positive definite over the range of physically meaningful parameters (here: contact stiffnesses, solid volume fraction and coordination number) so that strictly speaking we are not dealing with a material instability.

Capillary break-up rheometry of low-viscosity elastic fluids

We investigate the dynamics of the capillary thinning and break-up process for low viscosity elastic fluids such as dilute polymer solutions. Standard measurements of the evolution of the midpoint diameter of the necking fluid filament are augmented by high speed digital video images of the break up dynamics. We show that the successful operation of a capillary thinning device is governed by three important time scales (which characterize the relative importance of inertial, viscous and elastic processes), and also by two important length scales (which specify the initial sample size and the total stretch imposed on the sample). By optimizing the ranges of these geometric parameters, we are able to measure characteristic time scales for tensile stress growth as small as 1 millisecond for a number of model dilute and semi-dilute solutions of polyethylene oxide (PEO) in water and glycerol. If the final aspect ratio of the sample is too small, or the total axial stretch is too great, measurements are limited, respectively, by inertial oscillations of the liquid bridge or by the development of the well-known beads-on-a-string morphology which disrupt the formation of a uniform necking filament. By considering the magnitudes of the natural time scales associated with viscous flow, elastic stress growth and inertial oscillations it is possible to construct an operability diagram characterizing successful operation of a capillary break-up extensional rheometer. For Newtonian fluids, viscosities greater than approximately 70 mPas are required; however for dilute solutions of high molecular weight polymer, the minimum Viscosity is substantially lower due to the additional elastic stresses arising from molecular extension. For PEO of molecular weight 2.10(6) g/mol, it is possible to measure relaxation times of order 1 ms in dilute polymer solutions with zero-shear-rate viscosities on the order of 2-10 mPas.

Change in woolgrower profits from adopting staple strength-enhancing management strategies

Wool tenderness is a significant problem in Australia, especially in areas where sheep graze under highly seasonal conditions. In this study, a profit function model is specified, estimated and simulated to assess the economic impact of staple strength-enhancing research on the profits of Australian woolgrowers. The model is based on a number of fundamental characteristics of the Australian wool industry and the staple-strength enhancing technology being assessed. The model consists of a system of demand and supply equations that are specified in terms of effective, rather than actual, prices. The interrelationships between the inputs and outputs are allowed for in the model in a manner that is consistent with theoretical restrictions. The adoption of the new feed management strategy results in a 4.4% increase in the expected profits of Australian wool producers in the short-run, and a 2.2% increase in expected profits in the long-run.

Drop formation and breakup of low viscosity elastic fluids: Effects of molecular weight and concentration

The dynamics of drop formation and pinch-off have been investigated for a series of low viscosity elastic fluids possessing similar shear viscosities, but differing substantially in elastic properties. On initial approach to the pinch region, the viscoelastic fluids all exhibit the same global necking behavior that is observed for a Newtonian fluid of equivalent shear viscosity. For these low viscosity dilute polymer solutions, inertial and capillary forces form the dominant balance in this potential flow regime, with the viscous force being negligible. The approach to the pinch point, which corresponds to the point of rupture for a Newtonian fluid, is extremely rapid in such solutions, with the sudden increase in curvature producing very large extension rates at this location. In this region the polymer molecules are significantly extended, causing a localized increase in the elastic stresses, which grow to balance the capillary pressure. This prevents the necked fluid from breaking off, as would occur in the equivalent Newtonian fluid. Alternatively, a cylindrical filament forms in which elastic stresses and capillary pressure balance, and the radius decreases exponentially with time. A (0+1)-dimensional finitely extensible nonlinear elastic dumbbell theory incorporating inertial, capillary, and elastic stresses is able to capture the basic features of the experimental observations. Before the critical "pinch time" t(p), an inertial-capillary balance leads to the expected 2/3-power scaling of the minimum radius with time: R-min similar to(t(p)-t)(2/3). However, the diverging deformation rate results in large molecular deformations and rapid crossover to an elastocapillary balance for times t>t(p). In this region, the filament radius decreases exponentially with time R-min similar to exp[(t(p)-t)/lambda(1)], where lambda(1) is the characteristic time constant of the polymer molecules. Measurements of the relaxation times of polyethylene oxide solutions of varying concentrations and molecular weights obtained from high speed imaging of the rate of change of filament radius are significantly higher than the relaxation times estimated from Rouse-Zimm theory, even though the solutions are within the dilute concentration region as determined using intrinsic viscosity measurements. The effective relaxation times exhibit the expected scaling with molecular weight but with an additional dependence on the concentration of the polymer in solution. This is consistent with the expectation that the polymer molecules are in fact highly extended during the approach to the pinch region (i.e., prior to the elastocapillary filament thinning regime) and subsequently as the filament is formed they are further extended by filament stretching at a constant rate until full extension of the polymer coil is achieved. In this highly extended state, intermolecular interactions become significant, producing relaxation times far above theoretical predictions for dilute polymer solutions under equilibrium conditions. (C) 2006 American Institute of Physics

Differential Priors for Elastic Nets

The elastic net and related algorithms, such as generative topographic mapping, are key methods for discretized dimension-reduction problems. At their heart are priors that specify the expected topological and geometric properties of the maps. However, up to now, only a very small subset of possible priors has been considered. Here we study a much more general family originating from discrete, high-order derivative operators. We show theoretically that the form of the discrete approximation to the derivative used has a crucial influence on the resulting map. Using a new and more powerful iterative elastic net algorithm, we confirm these results empirically, and illustrate how different priors affect the form of simulated ocular dominance columns.

Bifurcation in growth patterns for arrays of parallel Griffith, edge and sliding cracks

This paper presents the recent finding by Muhlhaus et al [1] that bifurcation of crack growth patterns exists for arrays of two-dimensional cracks. This bifurcation is a result of the nonlinear effect due to crack interaction, which is, in the present analysis, approximated by the dipole asymptotic or pseudo-traction method. The nonlinear parameter for the problem is the crack length/ spacing ratio lambda = a/h. For parallel and edge crack arrays under far field tension, uniform crack growth patterns (all cracks having same size) yield to nonuniform crack growth patterns (i.e. bifurcation) if lambda is larger than a critical value lambda(cr) (note that such bifurcation is not found for collinear crack arrays). For parallel and edge crack arrays respectively, the value of lambda(cr) decreases monotonically from (2/9)(1/2) and (2/15.096)(1/2) for arrays of 2 cracks, to (2/3)(1/2)/pi and (2/5.032)(1/2)/pi for infinite arrays of cracks. The critical parameter lambda(cr) is calculated numerically for arrays of up to 100 cracks, whilst discrete Fourier transform is used to obtain the exact solution of lambda(cr) for infinite crack arrays. For geomaterials, bifurcation can also occurs when array of sliding cracks are under compression.

Changes in three-dimensional architecture of microfilaments in cultured vascular smooth muscle cells during phenotypic modulation

To investigate changes in the three-dimensional microfilament architecture of vascular smooth muscle cells (SMC) during the process of phenotypic modulation, rabbit aortic SMCs cultured under different conditions and at different time points were either labelled with fluorescein-conjugated probes to cytoskeletal and contractile proteins for observation by confocal laser scanning microscopy, or extracted with Triton X-100 for scanning electron microscopy. Densely seeded SMCs in primary culture, which maintain a contractile phenotype, display prominent linear myofilament bundles (stress fibres) that are present throughout the cytoplasm with alpha-actin filaments predominant in the central part and beta-actin filaments in the periphery of the cell. Intermediate filaments form a meshed network interconnecting the stress fibres and linking directly to the nucleus. Moderately and sparsely seeded SMCs, which modulate toward the synthetic phenotype during the first 5 days of culture, undergo a gradual redistribution of intermediate filaments from the perinuclear region toward the peripheral cytoplasm and a partial disassembly of stress fibres in the central part of the upper cortex of the cytoplasm, with an obvious decrease in alpha-actin and myosin staining. These changes are reversed in moderately seeded SMCs by day 8 of culture when they have reached confluence. The results reveal two changes in microfilament architecture in SMCs as they undergo a change in phenotype: the redistribution of intermediate filaments probably due to an increase in synthetic organelles in the perinuclear area, and the partial disassembly of stress fibres which may reflect a degradation of contractile components.

Effects of beta-escin and saponin on the transverse-tubular system and sarcoplasmic reticulum membranes of rat and toad skeletal muscle

Mechanically skinned skeletal muscle fibres from rat and toad were exposed to the permeabilizing agents beta-escin and saponin. The effects of these agents on the sealed transverse tubular system (t-system) and sarcoplasmic reticulum (SR) were examined by looking at changes in the magnitude of the force responses to t-system depolarization, the time course of the fluorescence of fura-2 trapped in the sealed t-system, and changes in the magnitude of caffeine-induced contractures following SR loading with Ca2+ under defined conditions. In the presence of 2 mu g ml(-1) beta-escin and saponin, the response to t-system depolarization was not completely abolished, decreasing to a plateau, and a large proportion of fura-2 remained in the sealed t-system. At 10 mu g ml(-1), both agents abolished the ability of both rat and toad preparations to respond to t-system depolarization after 3 min of exposure, but a significant amount of fura-2 remained in sealed t-tubules even after exposure to 100 mu g ml(-1) beta-escin and saponin for 10 min. beta-Escin took longer than saponin to reduce the t-system depolarizations and fura-2 content of the sealed t-system to a similar level. The ability of the SR to load Ca2+ was reduced to a lower level after treatment with beta-escin than saponin. This direct effect on the SR occurred at much lower concentrations for rat (2 mu g ml(-1) beta-escin and 10 mu g ml(-1) saponin) than toad (10 mu g ml(-1) beta-escin and 150 mu g ml(-1) saponin). The reverse order in sensitivities to beta-escin and saponin of t-system and SR membranes indicates that the mechanisms of action of beta-escin and saponin are different in the two types of membrane. In conclusion, this study shows that: (1) beta-escin has a milder action on the surface membrane than saponin; (2) beta-escin is a more potent modifier of SR function; (3) simple permeabilization of membranes is not sufficient to explain the effects of beta-escin and saponin on muscle membranes; and (4) the t-system network within muscle fibres is not a homogeneous compartment.

Quantum open-systems approach to current noise in resonant tunneling junctions

A quantum Markovian master equation is derived to describe the current noise in resonant tunneling devices. This equation includes both incoherent and coherent quantum tunneling processes. We show how to obtain the population master equation by adiabatic elimination of quantum coherences in the presence of elastic scattering. We calculate the noise spectrum for a double well device and predict subshot noise statistics for strong tunneling between the wells. The method is an alternative to Green's function methods and population master equations for very small coherently coupled quantum dots.

Modelling the large deformations in stratified media - the Cosserat continuum approach

Methods employing continuum approximation in describing the deformation of layered materials possess a clear advantage over explicit models, However, the conventional implicit models based on the theory of anisotropic continua suffers from certain difficulties associated with interface slip and internal instabilities. These difficulties can be remedied by considering the bending stiffness of the layers. This implies the introduction of moment (couple) stresses and internal rotations, which leads to a Cosserat-type theory. In the present model, the behaviour of the layered material is assumed to be linearly elastic; the interfaces are assumed to be elastic perfectly plastic. Conditions of slip or no slip at the interfaces are detected by a Coulomb criterion with tension cut off at zero normal stress. The theory is valid for large deformation analysis. The model is incorporated into the finite element program AFENA and validated against analytical solutions of elementary buckling problems in layered medium. A problem associated with buckling of the roof and the floor of a rectangular excavation in jointed rock mass under high horizontal in situ stresses is considered as the main application of the theory. Copyright (C) 1999 John Wiley & Sons, Ltd.