Rheometry of granular materials in cylindrical Couette cells: Anomalous stress caused by gravity and shear

The cylindrical Couette device is commonly employed to study the rheology of fluids, but seldom used for dense granular materials. Plasticity theories used for granular flows predict a stress field that is independent of the shear rate, but otherwise similar to that in fluids. In this paper we report detailed measurements of the stress as a function of depth, and show that the stress profile differs fundamentally from that of fluids, from the predictions of plasticity theories, and from intuitive expectation. In the static state, a part of the weight of the material is transferred to the walls by a downward vertical shear stress, bringing about the well-known Janssen saturation of the stress in vertical columns. When the material is sheared, the vertical shear stress changes sign, and the magnitudes of all components of the stress rise rapidly with depth. These qualitative features are preserved over a range of the Couette gap and shear rate, for smooth and rough walls and two model granular materials. To explain the anomalous rheological response, we consider some hypotheses that seem plausibleapriori, but showthat none survive after careful analysis of the experimental observations. We argue that the anomalous stress is due to an anisotropic fabric caused by the combined actions of gravity, shear, and frictional walls, for which we present indirect evidence from our experiments. A general theoretical framework for anisotropic plasticity is then presented. The detailed mechanics of how an anisotropic fabric is brought about by the above-mentioned factors is not clear, and promises to be a challenging problem for future investigations. (C) 2013 AIP Publishing LLC.

Simple theoretical analysis of the effective electron mass in semiconductor nanowires

In this paper we study the effective electron mass (EEM) in Nano wires (NWs) of nonlinear optical materials on the basis of newly formulated electron dispersion relation by considering all types of anisotropies of the energy band constants within the framework of k . p formalism. The results for NWs of III-V, ternary and quaternary semiconductors form special cases of our generalized analysis. We have also investigated the EEM in NWs of Bi, IV-VI, stressed Kane type materials, Ge, GaSb and Bi2Te3 by formulating the appropriate 1D dispersion law in each case by considering the influence of energy band constants in the respective cases. It has been found that the 1D EEM in nonlinear optical materials depend on the size quantum numbers and Fermi energy due to the anisotropic spin orbit splitting constant and the crystal field splitting respectively. The 1D EEM is Bi, IV-VI, stressed Kane type semiconductors and Ge also depends on both the Fermi energy and the size quantum numbers which are the characteristic features of such NWs. The EEM increases with increase in concentration and decreasing film thickness and for ternary and quaternary compounds the EEM increases with increase in alloy composition. Under certain special conditions all the results for all the materials get simplified into the well known parabolic energy bands and thus confirming the compatibility test.

Vorticity moments in four numerical simulations of the 3D Navier-Stokes equations

The issue of intermittency in numerical solutions of the 3D Navier-Stokes equations on a periodic box 0, L](3) is addressed through four sets of numerical simulations that calculate a new set of variables defined by D-m(t) = (pi(-1)(0) Omega(m))(alpha m) for 1 <= m <= infinity where alpha(m) = 2m/(4m - 3) and Omega(m)(t)](2m) = L-3 integral(v) vertical bar omega vertical bar(2m) dV with pi(0) = vL(-2). All four simulations unexpectedly show that the D-m are ordered for m = 1,..., 9 such that Dm+1 < D-m. Moreover, the D-m squeeze together such that Dm+1/D-m NE arrow 1 as m increases. The values of D-1 lie far above the values of the rest of the D-m, giving rise to a suggestion that a depletion of nonlinearity is occurring which could be the cause of Navier-Stokes regularity. The first simulation is of very anisotropic decaying turbulence; the second and third are of decaying isotropic turbulence from random initial conditions and forced isotropic turbulence at fixed Grashof number respectively; the fourth is of very-high-Reynolds-number forced, stationary, isotropic turbulence at up to resolutions of 4096(3).

Taming of large diameter triaxial setup

Triaxial tests are essential to estimate the shear strength properties of the soil or rock. Normally triaxial tests are carried out on samples of 38 mm diameter and 76 mm height. Granular materials, predominantly used in base/sub-base construction of pavements or in railways have size range of 60-75 mm. Determination of shear strength parameters of those materials can be made possible only through triaxial tests on large diameter samples. This paper describes a large diameter cyclic triaxial testing facility set up in the Geotechnical Engineering lab of Indian Institute of Science. This setup consists of 100 kN capacity dynamic loading frame, which facilitates testing of samples of up to 300 mm diameter and 600 mm height. The loading ram can be actuated up to a maximum frequency of 10 Hz, with maximum amplitude of 100 mm. The setup is capable of carrying out static as well as dynamic triaxial tests under isotropic, anisotropic conditions with a maximum confining pressure of 1 MPa. Working with this setup is a difficult task because of the size of the sample. In this paper, a detailed discussion on the various problems encountered during the initial testing using the equipment, the ideas and solutions adopted to solve them are presented. Pilot experiments on granular sub-base material of 53 mm down size are also presented.

On Uniform Approximate Solutions in Bending of Symmetric Laminated Plates

A layer-wise theory with the analysis of face ply independent of lamination is used in the bending of symmetric laminates with anisotropic plies. More realistic and practical edge conditions as in Kirchhoff's theory are considered. An iterative procedure based on point-wise equilibrium equations is adapted. The necessity of a solution of an auxiliary problem in the interior plies is explained and used in the generation of proper sequence of two dimensional problems. Displacements are expanded in terms of polynomials in thickness coordinate such that continuity of transverse stresses across interfaces is assured. Solution of a fourth order system of a supplementary problem in the face ply is necessary to ensure the continuity of in-plane displacements across interfaces and to rectify inadequacies of these polynomial expansions in the interior distribution of approximate solutions. Vertical deflection does not play any role in obtaining all six stress components and two in-plane displacements. In overcoming lacuna in Kirchhoff's theory, widely used first order shear deformation theory and other sixth and higher order theories based on energy principles at laminate level in smeared laminate theories and at ply level in layer-wise theories are not useful in the generation of a proper sequence of 2-D problems converging to 3-D problems. Relevance of present analysis is demonstrated through solutions in a simple text book problem of simply supported square plate under doubly sinusoidal load.

Influence of Spatial Variation of Hydraulic Conductivity of Municipal Solid Waste on Performance of Bioreactor Landfill

The current study analyzes the leachate distribution in the Orchard Hills Landfill, Davis Junction, Illinois, using a two-phase flow model to assess the influence of variability in hydraulic conductivity on the effectiveness of the existing leachate recirculation system and its operations through reliability analysis. Numerical modeling, using finite-difference code, is performed with due consideration to the spatial variation of hydraulic conductivity of the municipal solid waste (MSW). The inhomogeneous and anisotropic waste condition is assumed because it is a more realistic representation of the MSW. For the reliability analysis, the landfill is divided into 10 MSW layers with different mean values of vertical and horizontal hydraulic conductivities (decreasing from top to bottom), and the parametric study is performed by taking the coefficients of variation (COVs) as 50, 100, 150, and 200%. Monte Carlo simulations are performed to obtain statistical information (mean and COV) of output parameters of the (1) wetted area of the MSW, (2) maximum induced pore pressure, and (3) leachate outflow. The results of the reliability analysis are used to determine the influence of hydraulic conductivity on the effectiveness of the leachate recirculation and are discussed in the light of a deterministic approach. The study is useful in understanding the efficiency of the leachate recirculation system. (C) 2013 American Society of Civil Engineers.

Floquet generation of Majorana end modes and topological invariants

We show how Majorana end modes can be generated in a one-dimensional system by varying some of the parameters in the Hamiltonian periodically in time. The specific model we consider is a chain containing spinless electrons with a nearest-neighbor hopping amplitude, a p-wave superconducting term, and a chemical potential; this is equivalent to a spin-1/2 chain with anisotropic XY couplings between nearest neighbors and a magnetic field applied in the (z) over cap direction. We show that varying the chemical potential (or magnetic field) periodically in time can produce Majorana modes at the ends of a long chain. We discuss two kinds of periodic driving, periodic delta-function kicks, and a simple harmonic variation with time. We discuss some distinctive features of the end modes such as the inverse participation ratio of their wave functions and their Floquet eigenvalues which are always equal to +/- 1 for time-reversal-symmetric systems. For the case of periodic delta-function kicks, we use the effective Hamiltonian of a system with periodic boundary conditions to define two topological invariants. The first invariant is a well-known winding number, while the second invariant has not appeared in the literature before. The second invariant is more powerful in that it always correctly predicts the numbers of end modes with Floquet eigenvalues equal to + 1 and -1, while the first invariant does not. We find that the number of end modes can become very large as the driving frequency decreases. We show that periodic delta-function kicks in the hopping and superconducting terms can also produce end modes. Finally, we study the effect of electron-phonon interactions (which are relevant at finite temperatures) and a random noise in the chemical potential on the Majorana modes.

Forming limit investigation of friction stir welded sheets: influence of shoulder diameter and plunge depth

The main aim of the present work is to analyze the influence of shoulder diameter and plunge depth on the formability of friction stir welded sheets. The base material used for welding and forming was AA6061-T6. Formability evaluation was performed through limiting dome height tests. The forming limit curve, FLC (only in the stretching region), thickness distribution, and strain hardening exponent of the weld region were monitored during formability studies. It is found from the work that the forming limit of friction stir welded sheets is better than unwelded sheets. In general, with an increase in shoulder diameter and plunge depth, the forming limit is found to improve considerably. With a decrease in thickness gradient severity and an increase in strain hardening exponent (n) of the weld region, the forming limit is found to increase. The increase in n value of the weld region is believed to occur because of the reduction in dislocation density. The maximum thickness difference is higher in the retreating side, rather than in the advancing side, of the weld. This is due to the differential straining and hardness levels attained by both sides during friction stir welding.

A systematic investigation of fracture mechanisms in Al-Si based eutectic alloy-Effect of Si modification

The fracture characteristics of Al-Si based eutectic alloy are investigated in the unmodified and modified conditions under compression. The investigations are carried out at different strain rates and temperatures. Fracture of the alloy starts with eutectic Si particle fracture and modification plays an important role in particle fracture. The fraction of fractured particles is found to be always lesser in the modified condition than in the unmodified condition. Particle fracture increases with increase in strain. It is found that the Si particle fracture shows an increase with increase in strain rate and decreases with increase in temperature at 10% strain. Large and elongated particles show a greater tendency for fracture in the unmodified and modified conditions. Particle orientation plays an important role on fracture and the cracks are found to occur almost in a direction normal to the tensile strain imposed upon the particles by the deforming matrix in the unmodified alloy. The modified alloy shows a random distribution of fractured particles and crack orientation. The criteria of fracture based on dislocation pile-up mechanism and fiber loading explain the observed difference in particle fracture characteristics due to modification. The particle fracture for the modified alloy is also discussed in terms of Weibull statistics and the existing models of dispersion hardening. Particle/matrix interface decohesion is observed at higher strain rates and temperatures in the modified alloy. Dendritic rotation of 10 degrees is also observed at higher strain rates, which can increase the amount of particle fracture. (C) 2013 Elsevier B.V. All rights reserved.

Solid-solid collapse transition in a two dimensional model molecular system

Solid-solid collapse transition in open framework structures is ubiquitous in nature. The real difficulty in understanding detailed microscopic aspects of such transitions in molecular systems arises from the interplay between different energy and length scales involved in molecular systems, often mediated through a solvent. In this work we employ Monte-Carlo simulation to study the collapse transition in a model molecular system interacting via both isotropic as well as anisotropic interactions having different length and energy scales. The model we use is known as Mercedes-Benz (MB), which, for a specific set of parameters, sustains two solid phases: honeycomb and oblique. In order to study the temperature induced collapse transition, we start with a metastable honeycomb solid and induce transition by increasing temperature. High density oblique solid so formed has two characteristic length scales corresponding to isotropic and anisotropic parts of interaction potential. Contrary to the common belief and classical nucleation theory, interestingly, we find linear strip-like nucleating clusters having significantly different order and average coordination number than the bulk stable phase. In the early stage of growth, the cluster grows as a linear strip, followed by branched and ring-like strips. The geometry of growing cluster is a consequence of the delicate balance between two types of interactions, which enables the dominance of stabilizing energy over destabilizing surface energy. The nucleus of stable oblique phase is wetted by intermediate order particles, which minimizes the surface free energy. In the case of pressure induced transition at low temperature the collapsed state is a disordered solid. The disordered solid phase has diverse local quasi-stable structures along with oblique-solid like domains. (C) 2013 AIP Publishing LLC.

WEAR AND SLIDING-FRICTION CHARACTERISTICS OF STAINLESS STEEL SS316LN WITHOUT A COATING AND WITH A CrN COATING AT ELEVATED TEMPERATURES AND IN A VACUUM

Stainless steel of type AISI 316LN - one of the structural materials of fast neutron reactors - must have a long service life under conditions that subject it to different types of wear (galling, adhesion, fretting, and abrasion). Cobalt-based hard facings are generally avoided due to induced radioactivity. Nickel-based hard facings are strongly preferred instead. One alternative to both types of coatings is a hard-alloy coating of CrN. This article examines wear and friction characteristics during the sliding of uncoated steel SS316LN and the same steel with a CrN coating. In addition, a specially designed pin-on-disk tribometer is used to perform tests in a vacuum at temperatures of up to 1000 degrees C in order to study the effect of oxygen on the wear of these materials. The morphology of the wear surface and the structure of the subsurface were studied by scanning electron microscopy. The formation of an adhesion layer and the self-welding of mating parts are seen to take place in the microstructure at temperatures above 500 degrees C. It is also found that steel SS316LN undergoes shear strain during sliding wear. The friction coefficient depends on the oxygen content, load, and temperature, while the wear rate depends on the strain-hardening of the surface of the material being tested.

ENDOR studies of hydrogen bonding of the AsO44- centre in x-ray irradiated KH2AsO4

We here report what we believe to be an important method for studying hydrogen bonding in systems containing a paramagnetic centre. The technique of electron-nuclear double resonance ( ENDOR) has been applied to study the hydrogen-bond network around the AsO44-. centre in X-ray irradiated KH2AsO4. ENDOR transitions from several sets of hydrogen nuclei surrounding the centre were observed at 4.2 degrees K and the spectra for two sets of neighbouring nuclei are identified. The angular dependences for these spectra are fitted with a spin-Hamiltonian to obtain the isotropic and anisotropic magnetic hyperfine constants. The results are discussed in terms of the available spectroscopic and crystallographic data on KH2AsO4 and the order-disorder model of ferroelectrictricity in this class of crystals.

Non-Destructive Evaluation of Defects and Damage in Composite Materials and Structures

The mechanical behaviour of composite materials differs from that of conventional structural materials owing to their heterogeneous and anisotropic nature. Different types of defects and anomalies get induced in these materials during the fabrication process. Further, during their service life, the components made of composite materials develop different types of damage. The performance and life of such components is governed by the combined effect of all these defects and damage. While porosity, voids, inclusions etc., are some defects those can get induced during the fabrication of composites, matrix cracks, interface debonds, delaminations and fiber breakage are major types of service induced damage which are of concern. During the service life of components made of composites, one type of damage can grow and initiate another type of damage. For example, matrix cracks can gradually grow to the interface and initiate debonds. Interface debonds in a particular plane can lead to delaminations. Consequently, the combined effect of different types of distributed damage causes the failure of the component. A set of non-destructive evaluation (NDE) methods is well established for testing conventional metallic materials. Some of them can also be utilized for composite materials as they are, and in some cases with a little different approach or modification. Ultrasonics, Radiography, Thermography, Fiber Optics, Acoustic Emision Techniques etc., to name a few. Detection, evaluation and characterization of different types of defects and damage encountered in composite materials and structures using different NDE tools is discussed briefly in this paper.

Micropillar and macropillar compression responses of magnesium single crystals oriented for single slip or extension twinning

Uniaxial compression experiments were conducted on two magnesium (Mg) single crystals whose crystallographic orientations facilitate the deformation either by basal slip or by extension twinning. Specimen size effects were examined by conducting experiments on mu m- and mm-sized samples. A marked specimen size effect was noticed, with micropillars exhibiting significantly higher flow stress than bulk samples. Further, it is observed that the twin nucleation stress exerts strong size dependence, with micropillars requiring substantially higher stress than the bulk samples. The flow curves obtained on the bulk samples are smooth whereas those obtained from micropillars exhibit intermittent and precipitous stress drops. Electron backscattered diffraction and microstructural analyses of the deformed samples reveal that the plastic deformation in basal slip oriented crystals occurs only by slip while twin oriented crystals deform by both slip and twinning modes. The twin oriented crystals exhibit a higher strain hardening during plastic deformation when compared to the single slip oriented crystals. The strain hardening rate, theta, of twin oriented crystals is considerably greater in micropillars compared to the bulk single crystals, suggesting the prevalence of different work hardening mechanisms at these different sample sizes. (C) 2013 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

A finite element method for the Saint-Venant torsion and bending problems for prismatic beams

In this work, we present a finite element formulation for the Saint-Venant torsion and bending problems for prismatic beams. The torsion problem formulation is based on the warping function, and can handle multiply-connected regions (including thin-walled structures), compound and anisotropic bars. Similarly, the bending formulation, which is based on linearized elasticity theory, can handle multiply-connected domains including thin-walled sections. The torsional rigidity and shear centers can be found as special cases of these formulations. Numerical results are presented to show the good coarse-mesh accuracy of both the formulations for both the displacement and stress fields. The stiffness matrices and load vectors (which are similar to those for a variable body force in a conventional structural mechanics problem) in both formulations involve only domain integrals, which makes them simple to implement and computationally efficient. (C) 2014 Elsevier Ltd. All rights reserved.