On numerical implementation of an isotropic elastic-viscoplastic constitutive model for bulk metallic glasses

An efficient algorithm within the finite deformation framework is developed for finite element implementation of a recently proposed isotropic, Mohr-Coulomb type material model, which captures the elastic-viscoplastic, pressure sensitive and plastically dilatant response of bulk metallic glasses. The constitutive equations are first reformulated and implemented using an implicit numerical integration procedure based on the backward Euler method. The resulting system of nonlinear algebraic equations is solved by the Newton-Raphson procedure. This is achieved by developing the principal space return mapping technique for the present model which involves simultaneous shearing and dilatation on multiple potential slip systems. The complete stress update algorithm is presented and the expressions for viscoplastic consistent tangent moduli are derived. The stress update scheme and the viscoplastic consistent tangent are implemented in the commercial finite element code ABAQUS/Standard. The accuracy and performance of the numerical implementation are verified by considering several benchmark examples, which includes a simulation of multiple shear bands in a 3D prismatic bar under uniaxial compression.

The shear dynamo problem for small magnetic Reynolds numbers

We study large-scale kinematic dynamo action due to turbulence in the presence of a linear shear flow in the low-conductivity limit. Our treatment is non-perturbative in the shear strength and makes systematic use of both the shearing coordinate transformation and the Galilean invariance of the linear shear flow. The velocity fluctuations are assumed to have low magnetic Reynolds number (Re-m), but could have arbitrary fluid Reynolds number. The equation for the magnetic fluctuations is expanded perturbatively in the small quantity, Re-m. Our principal results are as follows: (i) the magnetic fluctuations are determined to the lowest order in Rem by explicit calculation of the resistive Green's function for the linear shear flow; (ii) the mean electromotive force is then calculated and an integro-differential equation is derived for the time evolution of the mean magnetic field. In this equation, velocity fluctuations contribute to two different kinds of terms, the 'C' and 'D' terms, respectively, in which first and second spatial derivatives of the mean magnetic field, respectively, appear inside the space-time integrals; (iii) the contribution of the D term is such that its contribution to the time evolution of the cross-shear components of the mean field does not depend on any other components except itself. Therefore, to the lowest order in Re-m, but to all orders in the shear strength, the D term cannot give rise to a shear-current-assisted dynamo effect; (iv) casting the integro-differential equation in Fourier space, we show that the normal modes of the theory are a set of shearing waves, labelled by their sheared wavevectors; (v) the integral kernels are expressed in terms of the velocity-spectrum tensor, which is the fundamental dynamical quantity that needs to be specified to complete the integro-differential equation description of the time evolution of the mean magnetic field; (vi) the C term couples different components of the mean magnetic field, so they can, in principle, give rise to a shear-current-type effect. We discuss the application to a slowly varying magnetic field, where it can be shown that forced non-helical velocity dynamics at low fluid Reynolds number does not result in a shear-current-assisted dynamo effect.

High-temperature low-cycle fatigue behavior of a NIMONIC PE-16 superalloy—Correlation with deformation and fracture

Low-cycle fatigue (LCF) responses of NIMONIC PE-16 for various prior microstructures and strain amplitudes have been evaluated and the fatigue behavior has been explained in terms of the operative deformation mechanisms. Total strain-controlled LCF tests were performed at 923 K on samples possessing three different prior microstructures: alloy A in solution-annealed condition (free of γ′ and carbides), alloy B with double aging treatment (spherical γ′ of 18-nm diameter and M23C6), and alloy C with another double aging treatment (γ′ of size 35 nm, MC and M23C6). All three microstructures exhibited an intial cyclic hardening followed by a period of gradual softening at 923 K. Coffin-Manson plots describing the plastic strain amplitudevs number of reversals to failure showed that alloy A had maximum fatigue life while C showed the least. Alloy B exhibited a two-slope behavior in the Coffin-Manson plot over the strain amplitudes investigated. This has been ascribed to the change in the degree of homogeneity of deformation at high and low strain amplitudes. Transmission electron microscopic studies were carried out to characterize the various deformation mechanisms and precipitation reactions occurring during fatigue testign. Fresh precipitation of fine γ′ was confirmed by the development of “mottled contrast” in alloy C. Evidence for the shearing of the ordered γ′ precipitates was revealed by the presence of superdislocations in alloy C. Repeated shearing during cyclic loading led to the reduction in the size of the γ′ and consequent softening. Coarser γ′ precipitates were associated with Orowan loops. The observed fatigue behavior has been rationalized based on the micromechanisms stated above and on the degree of homogenization of slip assessed by slipband spacing measurements on tested samples.

Effect of strain rate on the high-temperature low-cycle fatigue properties of a nimonic PE-16 superalloy

Strain-rate effects on the low-cycle fatigue (LCF) behavior of a NIMONIC PE-16 superalloy have been evaluated in the temperature range of 523 to 923 K. Total-strain-controlled fatigue tests were per-formed at a strain amplitude of +/-0.6 pct on samples possessing two different prior microstructures: microstructure A, in the solution-annealed condition (free of gamma' and carbides); and microstructure B, in a double-aged condition with gamma' of 18-nm diameter and M23C6 carbides. The cyclic stress response behavior of the alloy was found to depend on the prior microstructure, testing temperature, and strain rate. A softening regime was found to be associated with shearing of ordered gamma' that were either formed during testing or present in the prior microstructure. Various manifestations of dynamic strain aging (DSA) included negative strain rate-stress response, serrations on the stress-strain hysteresis loops, and increased work-hardening rate. The calculated activation energy matched well with that for self-diffusion of Al and Ti in the matrix. Fatigue life increased with an increase in strain rate from 3 x 10(-5) to 3 x 10(-3) s-1, but decreased with further increases in strain rate. At 723 and 823 K and low strain rates, DSA influenced the deformation and fracture behavior of the alloy. Dynamic strain aging increased the strain localization in planar slip bands, and impingement of these bands caused internal grain-boundary cracks and reduced fatigue life. However, at 923 K and low strain rates, fatigue crack initiation and propagation were accelerated by high-temperature oxidation, and the reduced fatigue life was attributed to oxidation-fatigue interaction. Fatigue life was maximum at the intermediate strain rates, where strain localization was lower. Strain localization as a function of strain rate and temperature was quantified by optical and scanning electron microscopy and correlated with fatigue life.

Dynamic recrystallisation during hot working of Zr-2 center dot 5Nb: Characterisation using processing maps

The characteristics of the hot deformation of Zr-2.5Nb (wt-%) in the temperature range 650-950 degrees C and in the strain rate range 0.001-100 s(-1) have been studied using hot compression testing. Two different preform microstructures: equiaxed (alpha + beta) and beta transformed have been investigated. For this study, the approach of processing maps has been adopted and their interpretation carried out using the dynamic materials model. The efficiency of power dissipation given by [2m/(m + 1)], where m is the strain rate sensitivity, is plotted as a function of temperature and strain rate to obtain a processing map. A domain of dynamic recrystallisation has been identified in the maps of equiaxed (alpha + beta) and beta transformed preforms. In the case of equiaxed (alpha + beta), the stress-strain curves are steady state and the dynamic recrystallisation domain in the map occurs with a peak efficiency of 45% at 850 degrees C and 0.001 s(-1). On the other hand the beta transformed preform exhibits stress-strain curves with continuous flow softening. The corresponding processing map shows a domain of dynamic recrystallisation occurring by the shearing of alpha platelets followed by globularisation with a peak efficiency of 54% at 750 degrees C and 0.001 s(-1). The characteristics of dynamic recrystallisation are analysed on the basis of a simple model which considers the rates of nucleation and growth of recrystallised gains. Calculations show that these two rates are nearly equal and that the nucleation of dynamic recrystallisation is essentially controlled by mechanical recovery involving the cross-slip of screw dislocations. Analysis of flow instabilities using a continuum criterion revealed that Zi-2.5Nb exhibits flow localisation at temperatures lower than 700 degrees C and strain rates higher than 1 s(-1).

Characterization of hot deformation behaviour of Zr---2.5Nb---0.5Cu using processing maps

The characteristics of hot deformation of beta-quenched Zr-2.5Nb-0.5Cu in the temperature range 650-1050 degrees C and in the strain rate range 0.001-100 s(-1) have been studied using hot compression testing. For this study, the approach of processing maps has been adopted and their interpretation done using the Dynamic Materials Model. The efficiency of power dissipation given by [2m/(m + 1)], where m is strain rate sensitivity, is plotted as a function of temperature and strain rate to obtain a processing map. The processing map for Zr-2.5Nb-0.5Cu within (alpha + beta) phase field showed a domain of dynamic recrystallization, occurring by shearing of alpha-platelets followed by spheroidization, with a peak efficiency of 48% at 750 degrees C and 0.001 s(-1). The stress-strain curves in this domain had features of continuous flow softening and all these are similar to that in Zr-2.5Nb alloy. In the beta-phase field, a second domain with a peak efficiency of 47% occurred at 1050 degrees C and 0.001 s(-1) and this domain is correlated with the superplasticity of beta-phase. The beta-deformation characteristics of this alloy are similar to that observed in pure beta-zirconium with large grain size. Analysis of flow instabilities using a continuum criterion revealed that the Zr-2.5Nb-0.5Cu exhibits flow localization at temperatures higher than 800 degrees C and strain rates higher than about 30 s(-1) and that the addition of copper to Zr-2.5Nb reduces its susceptibility to flow instability, particularly in the (alpha + beta) phase field.

Swing amplification of nonaxisymmetric perturbations in stars and gas in a sheared galactic disk

We present a study of the growth of local, nonaxisymmetric perturbations in gravitationally coupled stars and gas in a differentially rotating galactic disk. The stars and gas are treated as two isothermal fluids of different velocity dispersions, with the stellar velocity dispersion being greater than that for the gas. We examine the physical effects of inclusion of a low-velocity dispersion component (gas) on the growth of non-axisymmetric perturbations in both stars and gas, as done for the axisymmetric case by Jog & Solomon. The amplified perturbations in stars and gas constitute trailing, material, spiral features which may be identified with the local spiral features seen in all spiral galaxies. The formulation of the two-fluid equations closely follows the one-fluid treatment by Goldreich & Lynden-Bell. The local, linearized perturbation equations in the sheared frame are solved to obtain the results for a temporary growth via swing amplification. The problem is formulated in terms of five dimensionless parameters-namely, the Q-factors for stars and gas, respectively; the gas mass fraction; the shearing rate in the galactic disk; and the length scale of perturbation. By using the observed values of these parameters, we obtain the amplifications and the pitch angles for features in stars and gas for dynamically distinct cases, as applicable for different regions of spiral galaxies. A real galaxy consisting of stars and gas may display growth of nonaxisymmetric perturbations even when it is stable against axisymmetric perturbations and/or when either fluid by itself is stable against non-axisymmetric perturbations. Due to its lower velocity dispersion, the gas exhibits a higher amplification than do the stars, and the amplified gas features are slightly more tightly wound than the stellar features. When the gas contribution is high, the stellar amplification and the range of pitch angles over which it can occur are both increased, due to the gravitational coupling between the two fluids. Thus, the two-fluid scheme can explain the origin of the broad spiral arms in the underlying old stellar populations of galaxies, as observed by Schweizer and Elmegreen & Elmegreen. The arms are predicted to be broader in gas-rich galaxies, as is indeed seen for example in M33. In the linear regime studied here, the arm contrast is shown to increase with radius in the inner Galaxy, in agreement with observations of external galaxies by Schweizer. These results follow directly due to the inclusion of gas in the problem.

Effect of substrate on particle arrangement in arrays formed by self-assembly of polymer grafted nanoparticles

We show that the substrate affects the interparticle spacing in monolayer arrays with hexagonal order formed by self-assembly of polymer grafted nanoparticles. Remarkably, arrays with square packing were formed due to convective shearing at a liquid surface induced by miscibility of colloidal solution with the substrate.

The effect of boundaries on the plane Couette flow of granular materials: a bifurcation analysis

The tendency of granular materials in rapid shear flow to form non-uniform structures is well documented in the literature. Through a linear stability analysis of the solution of continuum equations for rapid shear flow of a uniform granular material, performed by Savage (1992) and others subsequently, it has been shown that an infinite plane shearing motion may be unstable in the Lyapunov sense, provided the mean volume fraction of particles is above a critical value. This instability leads to the formation of alternating layers of high and low particle concentrations oriented parallel to the plane of shear. Computer simulations, on the other hand, reveal that non-uniform structures are possible even when the mean volume fraction of particles is small. In the present study, we have examined the structure of fully developed layered solutions, by making use of numerical continuation techniques and bifurcation theory. It is shown that the continuum equations do predict the existence of layered solutions of high amplitude even when the uniform state is linearly stable. An analysis of the effect of bounding walls on the bifurcation structure reveals that the nature of the wall boundary conditions plays a pivotal role in selecting that branch of non-uniform solutions which emerges as the primary branch. This demonstrates unequivocally that the results on the stability of bounded shear how of granular materials presented previously by Wang et al. (1996) are, in general, based on erroneous base states.

On the criteria for reverse transition in a two-dimensional boundary layer flow

Experiments on reverse transition were conducted in two-dimensional accelerated incompressible turbulent boundary layers. Mean velocity profiles, longitudinal velocity fluctuations $\tilde{u}^{\prime}(=(\overline{u^{\prime 2}})^{\frac{1}{2}})$ and the wall-shearing stress (TW) were measured. The mean velocity profiles show that the wall region adjusts itself to laminar conditions earlier than the outer region. During the reverse transition process, increases in the shape parameter (H) are accompanied by a decrease in the skin friction coefficient (Cf). Profiles of turbulent intensity (u’2) exhibit near similarity in the turbulence decay region. The breakdown of the law of the wall is characterized by the parameter \[ \Delta_p (=\nu[dP/dx]/\rho U^{*3}) = - 0.02, \] where U* is the friction velocity. Downstream of this region the decay of $\tilde{u}^{\prime}$ fluctuations occurred when the momentum thickness Reynolds number (R) decreased roughly below 400.

The effect of boundaries on the plane Couette flow of granular materials: a bifurcation analysis

The tendency of granular materials in rapid shear ow to form non-uniform structures is well documented in the literature. Through a linear stability analysis of the solution of continuum equations for rapid shear flow of a uniform granular material, performed by Savage (1992) and others subsequently, it has been shown that an infinite plane shearing motion may be unstable in the Lyapunov sense, provided the mean volume fraction of particles is above a critical value. This instability leads to the formation of alternating layers of high and low particle concentrations oriented parallel to the plane of shear. Computer simulations, on the other hand, reveal that non-uniform structures are possible even when the mean volume fraction of particles is small. In the present study, we have examined the structure of fully developed layered solutions, by making use of numerical continuation techniques and bifurcation theory. It is shown that the continuum equations do predict the existence of layered solutions of high amplitude even when the uniform state is linearly stable. An analysis of the effect of bounding walls on the bifurcation structure reveals that the nature of the wall boundary conditions plays a pivotal role in selecting that branch of non-uniform solutions which emerges as the primary branch. This demonstrates unequivocally that the results on the stability of bounded shear flow of granular materials presented previously by Wang et al. (1996) are, in general, based on erroneous base states.

Model for Predicting the Mean Drop Size in Effervescent Sprays

This paper is focused on the development of a model for predicting the mean drop size in effervescent sprays. A combinatorial approach is followed in this modeling scheme, which is based on energy and entropy principles. The model is implemented in cascade in order to take primary breakup (due to exploding gas bubbles) and secondary breakup (due to shearing action of surrounding medium) into account. The approach in this methodology is to obtain the most probable drop size distribution by maximizing the entropy while satisfying the constraints of mass and energy balance. The comparison of the model predictions with the past experimental data is presented for validation. A careful experimental study is conducted over a wide range of gas-to-liquid ratios, which shows a good agreement with the model predictions: It is observed that the model gives accurate results in bubbly and annular flow regimes. However, discrepancies are observed in the transitional slug flow regime of the atomizer.

Performance analysis for geometrical attack on digital image watermarking

We present a technique for irreversible watermarking approach robust to affine transform attacks in camera, biomedical and satellite images stored in the form of monochrome bitmap images. The watermarking approach is based on image normalisation in which both watermark embedding and extraction are carried out with respect to an image normalised to meet a set of predefined moment criteria. The normalisation procedure is invariant to affine transform attacks. The result of watermarking scheme is suitable for public watermarking applications, where the original image is not available for watermark extraction. Here, direct-sequence code division multiple access approach is used to embed multibit text information in DCT and DWT transform domains. The proposed watermarking schemes are robust against various types of attacks such as Gaussian noise, shearing, scaling, rotation, flipping, affine transform, signal processing and JPEG compression. Performance analysis results are measured using image processing metrics.

Estimation of the shear transformation zone size fin a bulk metallic glass through statistical a analysis of the first pop-in stresses during spherical nanoindentation

The size of the shear transformation zone (STZ) that initiates the elastic to plastic transition in a Zr-based bulk metallic glass was estimated by conducting a statistical analysis of the first pop-in event during spherical nanoindentation. A series of experiments led us to a successful description of the distribution of shear strength for the transition and its dependence on the loading rate. From the activation volume determined by statistical analysis the STZ size was estimated based on a cooperative shearing model. (C) 2012 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

The effect of base roughness on the development of a dense granular flow down an inclined plane

The development of the flow of a granular material down an inclined plane starting from rest is studied as a function of the base roughness. In the simulations, the particles are rough frictional spheres interacting via the Hertz contact law. The rough base is made of a random configuration of fixed spheres with diameter different from the flowing particles, and the base roughness is decreased by decreasing the diameter of the base particles. The transition from an ordered to a disordered flowing state at a critical value of the base particle diameter, first reported by Kumaran and Maheshwari Phys. Fluids 24, 053302 (2012)] for particles with the linear contact model, is observed for the Hertzian contact model as well. The flow development for the ordered and disordered flows is very different. During the development of the disordered flow for the rougher base, there is shearing throughout the height. During the development of the ordered flow for the smoother base, there is a shear layer at the bottom and a plug region with no internal shearing above. In the shear layer, the particles are layered and hexagonally ordered in the plane parallel to the base, and the velocity profile is well approximated by Bagnold law. The flow develops in two phases. In the first phase, the thickness of the shear layer and the maximum velocity increase linearly in time till the shear front reaches the top. In the second phase, after the shear layer encompasses the entire flow, there is a much slower increase in the maximum velocity until the steady state is reached. (C) 2013 AIP Publishing LLC.