Correlation equations for average deposition rate coefficients of nanoparticles in a cylindrical pore

Nanoparticle deposition behavior observed at the Darcy scale represents an average of the processes occurring at the pore scale. Hence, the effect of various pore-scale parameters on nanoparticle deposition can be understood by studying nanoparticle transport at pore scale and upscaling the results to the Darcy scale. In this work, correlation equations for the deposition rate coefficients of nanoparticles in a cylindrical pore are developed as a function of nine pore-scale parameters: the pore radius, nanoparticle radius, mean flow velocity, solution ionic strength, viscosity, temperature, solution dielectric constant, and nanoparticle and collector surface potentials. Based on dominant processes, the pore space is divided into three different regions, namely, bulk, diffusion, and potential regions. Advection-diffusion equations for nanoparticle transport are prescribed for the bulk and diffusion regions, while the interaction between the diffusion and potential regions is included as a boundary condition. This interaction is modeled as a first-order reversible kinetic adsorption. The expressions for the mass transfer rate coefficients between the diffusion and the potential regions are derived in terms of the interaction energy profile. Among other effects, we account for nanoparticle-collector interaction forces on nanoparticle deposition. The resulting equations are solved numerically for a range of values of pore-scale parameters. The nanoparticle concentration profile obtained for the cylindrical pore is averaged over a moving averaging volume within the pore in order to get the 1-D concentration field. The latter is fitted to the 1-D advection-dispersion equation with an equilibrium or kinetic adsorption model to determine the values of the average deposition rate coefficients. In this study, pore-scale simulations are performed for three values of Peclet number, Pe = 0.05, 5, and 50. We find that under unfavorable conditions, the nanoparticle deposition at pore scale is best described by an equilibrium model at low Peclet numbers (Pe = 0.05) and by a kinetic model at high Peclet numbers (Pe = 50). But, at an intermediate Pe (e.g., near Pe = 5), both equilibrium and kinetic models fit the 1-D concentration field. Correlation equations for the pore-averaged nanoparticle deposition rate coefficients under unfavorable conditions are derived by performing a multiple-linear regression analysis between the estimated deposition rate coefficients for a single pore and various pore-scale parameters. The correlation equations, which follow a power law relation with nine pore-scale parameters, are found to be consistent with the column-scale and pore-scale experimental results, and qualitatively agree with the colloid filtration theory. These equations can be incorporated into pore network models to study the effect of pore-scale parameters on nanoparticle deposition at larger length scales such as Darcy scale.

Microstructural evolution and strength variability in microwires

Tensile experiments on cold-drawn Ni microwires with diameters from similar to 115 to 50 gm revealed high strengths, with significant strength variability for finer wires with diameters less than similar to 50 gm. The wires showed pronounced necking at fracture. The coarser wires with diameters > 50 mu m exhibited conventional ductile cup-cone fracture, with dimples in the central zone and peripheral shear lips, whereas finer wires failed by shear with knife or chisel-edge fractures. Shear bands were observed in all samples. Further, through- section microscopy of selected fractured samples revealed that the shear bands did not go across the enitre specimen for the coarser wires. The shear bands led to grain fragmention, with a reduction in grain aspect ratio as well as rotations away from the initial < 111 > orientations. The strength data were analysed based on a Weibull approach. The data could be rationalized in terms of failure from volume defects in coarser wires, with a high Weibull modulus, and from surface defects in finer wires, with a low Weibull modulus and greater variability. (C) 2015 Elsevier B.V. All rights reserved.

Structural evolution and second harmonic properties of lithium niobate-tantalate nanocrystals embedded in a borate glass

Nano-crystals of LiNbxTa1 (-) O-x(3) were evolved by subjecting melt-quenched 1.5Li(2)O-2B(2)O(3)-xNb(2)O(5)-(1 - x)Ta2O5 glasses (where x = 0, 0.25, 0.5, 0.75 and 1.00) to a controlled 3-h isothermal heat treatment between 530 and 560 degrees C. Detailed X-ray diffraction and Raman spectral studies confirmed the formation of nano-crystalline LiNbxTa1 (-) O-x(3) along with a minor phase of ferroelectric and non-linear optic Li2B4O7. The sizes of the nanocrystals evolved in the glass were in the range of 19-37 nm for x = 0-0.75 and 23-45 nm for x = 1.00. Electron microscopic studies confirmed a transformation of the morphology of the nano-crystallites from dendritic star-shaped spherulites for x = 0 to rod-shaped structures for x = 1.00 brought about by a coalescence of crystallites. Broad Maker-fringe patterns (recorded at 532 nm) were obtained by subjecting the heat-treated glass plates to 1064 nm fundamental radiation. However, an effective second order non-linear optic coefficient, d(eff), of 0.45 pm/V, which is nearly 1.2 times the d(36) of KDP single crystal, was obtained for a 560 degrees C/3 h heat-treated glass of the representative composition x = 0.50 comprising 37 nm sized crystallites. (C) 2015 Elsevier B.V. All rights reserved.

The time finite element as a robust general scheme for solving nonlinear dynamic equations including chaotic systems

Schemes that can be proven to be unconditionally stable in the linear context can yield unstable solutions when used to solve nonlinear dynamical problems. Hence, the formulation of numerical strategies for nonlinear dynamical problems can be particularly challenging. In this work, we show that time finite element methods because of their inherent energy momentum conserving property (in the case of linear and nonlinear elastodynamics), provide a robust time-stepping method for nonlinear dynamic equations (including chaotic systems). We also show that most of the existing schemes that are known to be robust for parabolic or hyperbolic problems can be derived within the time finite element framework; thus, the time finite element provides a unification of time-stepping schemes used in diverse disciplines. We demonstrate the robust performance of the time finite element method on several challenging examples from the literature where the solution behavior is known to be chaotic. (C) 2015 Elsevier Inc. All rights reserved.

The time finite element as a robust general scheme for solving nonlinear dynamic equations including chaotic systems

Schemes that can be proven to be unconditionally stable in the linear context can yield unstable solutions when used to solve nonlinear dynamical problems. Hence, the formulation of numerical strategies for nonlinear dynamical problems can be particularly challenging. In this work, we show that time finite element methods because of their inherent energy momentum conserving property (in the case of linear and nonlinear elastodynamics), provide a robust time-stepping method for nonlinear dynamic equations (including chaotic systems). We also show that most of the existing schemes that are known to be robust for parabolic or hyperbolic problems can be derived within the time finite element framework; thus, the time finite element provides a unification of time-stepping schemes used in diverse disciplines. We demonstrate the robust performance of the time finite element method on several challenging examples from the literature where the solution behavior is known to be chaotic. (C) 2015 Elsevier Inc. All rights reserved.

Evolution of Jahn-Teller distortion, transport and dielectric properties with doping in perovskite NdFe1-xMnxO3 (0 <= x <= 1) compounds

We have carried out dielectric and transport measurements in NdFe1-xMnxO3 (0 <= x <= 1) series of compounds and studied the variation of activation energy due to a change in Mn concentration. Despite similar ionic radii in Mn3+ and Fe3+, large variation is observed in the lattice parameters and a crossover from dynamic to static Jahn-Teller distortion is discernible. The Fe/Mn-O-Fe/Mn bond angle on the ab plane shows an anomalous change with doping. With an increase in the Mn content, the bond angle decreases until x = 0.6; beyond this, it starts rising until x = 0.8 and again falls after that. A similar trend is observed in activation energies estimated from both transport and dielectric relaxation by assuming a small polaron hopping (SPH) model. Impedance spectroscopy measurements delineate grain and grain boundary contributions separately both of which follow the SPH model. Frequency variation of the dielectric constant is in agreement with the modified Debye law from which relaxation dispersion is estimated.

Structural evolution and second harmonic properties of lithium niobate-tantalate nanocrystals embedded in a borate glass

Nano-crystals of LiNbxTa1 (-) O-x(3) were evolved by subjecting melt-quenched 1.5Li(2)O-2B(2)O(3)-xNb(2)O(5)-(1 - x)Ta2O5 glasses (where x = 0, 0.25, 0.5, 0.75 and 1.00) to a controlled 3-h isothermal heat treatment between 530 and 560 degrees C. Detailed X-ray diffraction and Raman spectral studies confirmed the formation of nano-crystalline LiNbxTa1 (-) O-x(3) along with a minor phase of ferroelectric and non-linear optic Li2B4O7. The sizes of the nanocrystals evolved in the glass were in the range of 19-37 nm for x = 0-0.75 and 23-45 nm for x = 1.00. Electron microscopic studies confirmed a transformation of the morphology of the nano-crystallites from dendritic star-shaped spherulites for x = 0 to rod-shaped structures for x = 1.00 brought about by a coalescence of crystallites. Broad Maker-fringe patterns (recorded at 532 nm) were obtained by subjecting the heat-treated glass plates to 1064 nm fundamental radiation. However, an effective second order non-linear optic coefficient, d(eff), of 0.45 pm/V, which is nearly 1.2 times the d(36) of KDP single crystal, was obtained for a 560 degrees C/3 h heat-treated glass of the representative composition x = 0.50 comprising 37 nm sized crystallites. (C) 2015 Elsevier B.V. All rights reserved.

Orientation Selection and Microstructural Evolution in Directionally Solidified Tb0.3Dy0.7Fe1.95

Tb0.3Dy0.7Fe1.95 alloy was directionally solidified by using a modified Bridgman technique at a wide range of growth rates of 5 to 100 cm/h. The directionally grown samples exhibited plane front solidification morphology up to a growth rate of 90 cm/h. Typical island banding feature was observed closer to the chilled end, which eventually gave rise to irregular peritectic coupled growth (PCG). The PCG gained prominence with an increase in the growth rate. The texture study revealed formation of strong aOE (c) 311 > texture in a lower growth rate regime, aOE (c) 110 > and ``rotated aOE (c) 110 > aEuroe in an intermediate growth regime, and aOE (c) 112 > in a higher growth rate regime. In-depth analysis of the atomic configuration of a solid-liquid interface revealed that the growth texture is influenced by the kinetics of atomic attachment to the solid-liquid interface, which is intimately related to a planar packing fraction and an atomic stacking sequence of the interfacial plane. The mechanism proposed in this article is novel and will be useful in addressing the orientation selection mechanism of topologically closed packed intermetallic systems. The samples grown at a higher growth rate exhibit larger magnetostriction (lambda) and d lambda/dH owing to the absence of pro-peritectic (Tb,Dy)Fe-3 and formation of aOE (c) 112 > texture, which lies closer to the easy magnetization direction (EMD).

Archean tectonics and crustal evolution of the Biligiri Rangan Block, southern India

The Southern Granulite Terrain in India is a collage of crustal blocks ranging in age from Archean to Neoproterozoic. This study investigate the tectonic evolution of one of the northernmost block- the Biligiri Block (BRB) through a multidisciplinary approach involving field investigation, petrographic studies, LA-ICPMS zircon U-Pb geochronology, Hf isotopic analyses, metamorphic P-T phase diagram computations, and crustal thickness modeling. The garnet bearing quartzofeldspathic gneiss from the central BRB preserve Mesoarchean magmatic zircons with ages between 3207 and 2806 Ma and positive epsilon Hf value (+2.7) which possibly indicates vestiges of a Mesoarchean primitive continental crust. The occurrence of quartzite-iron formation intercalation as well as ultramafic lenses along the western boundary of the BRB is interpreted to indicate that the Kollegal structural lineament is a possible paleo-suture. Phase diagram computation of a metagabbro from the southwestern periphery of the Kollegal suture zone reveals high-pressure (similar to 18.5 kbar) and medium-temperature (similar to 840 degrees C) metamorphism, likely during eastward subduction of the Western Dharwar oceanic crust beneath the Mesoarchean BRB. In the model presented here, slab subduction, melting and underplating processes generated arc magmatism and subsequent charnockitization within the BRB between ca. 2650 Ma and ca. 2498 Ma. These results thus reveal Meso- to Neoarchean tectonic evolution of the BRB. The spatial variation of crustal thickness, derived from flexure inversion technique, provides additional constraints on the tectonic linkage of the BRB with its surrounding terrains. In conjunction with published data, the Moyar and the Kollegal suture zones are considered to mark the trace of ocean closure along which the Nilgiri and Biligiri Rangan Blocks accreted on to the Western Dharwar Craton. (C) 2016 Elsevier B.V. All rights reserved.

Structure-rheology relationship in a sheared lamellar fluid

The structure-rheology relationship in the shear alignment of a lamellar fluid is studied using a mesoscale model which provides access to the lamellar configurations and the rheology. Based on the equations and free energy functional, the complete set of dimensionless groups that characterize the system are the Reynolds number (rho gamma L-2/mu), the Schmidt number (mu/rho D), the Ericksen number (mu(gamma)/B), the interface sharpness parameter r, the ratio of the viscosities of the hydrophilic and hydrophobic parts mu(r), and the ratio of the system size and layer spacing (L/lambda). Here, rho and mu are the fluid density and average viscosity, (gamma) over dot is the applied strain rate, D is the coefficient of diffusion, B is the compression modulus, mu(r) is the maximum difference in the viscosity of the hydrophilic and hydrophobic parts divided by the average viscosity, and L is the system size in the cross-stream direction. The lattice Boltzmann method is used to solve the concentration and momentum equations for a two dimensional system of moderate size (L/lambda = 32) and for a low Reynolds number, and the other parameters are systematically varied to examine the qualitative features of the structure and viscosity evolution in different regimes. At low Schmidt numbers where mass diffusion is faster than momentum diffusion, there is fast local formation of randomly aligned domains with ``grain boundaries,'' which are rotated by the shear flow to align along the extensional axis as time increases. This configuration offers a high resistance to flow, and the layers do not align in the flow direction even after 1000 strain units, resulting in a viscosity higher than that for an aligned lamellar phase. At high Schmidt numbers where momentum diffusion is fast, the shear flow disrupts layers before they are fully formed by diffusion, and alignment takes place by the breakage and reformation of layers by shear, resulting in defects (edge dislocations) embedded in a background of nearly aligned layers. At high Ericksen number where the viscous forces are large compared to the restoring forces due to layer compression and bending, shear tends to homogenize the concentration field, and the viscosity decreases significantly. At very high Ericksen number, shear even disrupts the layering of the lamellar phase. At low Ericksen number, shear results in the formation of well aligned layers with edge dislocations. However, these edge dislocations take a long time to anneal; the relatively small misalignment due to the defects results in a large increase in viscosity due to high layer stiffness and due to shear localization, because the layers between defects get pinned and move as a plug with no shear. An increase in the viscosity contrast between the hydrophilic and hydrophobic parts does not alter the structural characteristics during alignment. However, there is a significant increase in the viscosity, due to pinning of the layers between defects, which results in a plug flow between defects and a localization of the shear to a part of the domain.

Electrochemical deposition of manganese oxide-phosphate-reduced graphene oxide composite and electrocatalysis of the oxygen evolution reaction

A composite of manganese oxide and reduced graphene oxide (rGO) is prepared in a single step electrochemical reduction process in a phosphate buffer solution for studying as an electrocatalyst for the oxygen evolution reaction (OER). The novel composite catalyst, namely, MnOx-Pi-rGO, is electrodeposited from a suspension of graphene oxide (GO) in a neutral phosphate buffer solution containing KMnO4. The manganese oxide incorporates phosphate ions and deposits on the rGO sheet, which in turn is formed on the substrate electrode by electrochemical reduction of GO in the suspension. The OER is studied with the MnOx-Pi-rGO catalyst in a neutral phosphate electrolyte by linear sweep voltammetry. The results indicate a positive influence of rGO in the catalyst. By varying the ratio of KMnO4 and GO in the deposition medium and performing linear sweep voltammetry for the OER, the optimum composition of the deposition medium is obtained as 20 mM KMnO4 + 6.5% GO in 0.1 M phosphate buffer solution of pH 7. Under identical conditions, the MnOx-Pi-rGO catalyst exhibits 6.2 mA cm(-2) OER current against 2.9 mA cm(-2) by MnOx-Pi catalyst at 2.05 V in neutral phosphate solution. The Tafel slopes measured for OER at MnOx-Pi and MnOx-Pi-rGO are similar in magnitude at about 0.180 V decade(-1). The high Tafel slopes are attributed to partial dissolution of the catalyst during oxygen evolution. The O-2 evolved at the catalyst is measured by the water displacement method and the positive role of rGO on catalytic activity of MnOx-Pi is demonstrated.

Structural Evolution and Property Changes in Nd_{60}Al_{10}Fe_{20}Co_{10} Bulk Metallic Glass During Crystallization

The structural evolution and property changes in Nd60Al10Fe20Co10 bulk metallic glass (BMG) upon crystallization are investigated by the ultrasonic method, x-ray diffraction, density measurement, and differential scanning calorimetry. The elastic constants and Debye temperature of the BMG are obtained as a function of annealing temperature. Anomalous changes in ultrasonic velocities, elastic constants, and density are observed between 600–750 K, corresponding to the formation of metastable phases as an intermediate product in the crystallization process. The changes in acoustic velocities, elastic constants, density, and Debye temperature of the BMG relative to its fully crystallized state are much smaller, compared with those of other known BMGs, the differences being attributed to the microstructural feature of the BMG.

Spallation analysis with a closed trans-scale formulation damage evolution

A closed, trans-scale formulation of damage evolution based on the statistical microdamage mechanics is summarized in this paper. The dynamic function of damage bridges the mesoscopic and macroscopic evolution of damage. The spallation in an aluminium plate is studied with this formulation. It is found that the damage evolution is governed by several dimensionless parameters, i.e., imposed Deborah numbers De* and De, Mach number M and damage number S. In particular, the most critical mode of the macroscopic damage evolution, i.e., the damage localization, is deter-mined by Deborah number De+. Deborah number De* reflects the coupling and competition between the macroscopic loading and the microdamage growth. Therefore, our results reveal the multi-scale nature of spallation. In fact, the damage localization results from the nonlinearity of the microdamage growth. In addition, the dependence of the damage rate on imposed Deborah numbers De* and De, Mach number M and damage number S is discussed.

Numerical Investigation of Richtmyer-Meshkov Instability Driven by Cylindrical Shocks

In this paper, a numerical method with high order accuracy and high resolution was developed to simulate the Richtmyer-Meshkov(RM) instability driven by cylindrical shock waves. Compressible Euler equations in cylindrical coordinate were adopted for the cylindrical geometry and a third order accurate group control scheme was adopted to discretize the equations. Moreover, an adaptive grid technique was developed to refine the grid near the moving interface to improve the resolution of numerical solutions. The results of simulation exhibited the evolution process of RM instability, and the effect of Atwood number was studied. The larger the absolute value of Atwood number, the larger the perturbation amplitude. The nonlinear effect manifests more evidently in cylindrical geometry. The shock reflected from the pole center accelerates the interface for the second time, considerably complicating the interface evolution process, and such phenomena of reshock and secondary shock were studied.

Chaotic dynamic analysis of viscoelastic plates

The dynamic buckling of viscoelastic plates with large deflection is investigated in this paper by using chaotic and fractal theory. The material behavior is given in terms of the Boltzmann superposition principle. in order to obtain accurate computation results, the nonlinear integro-differential dynamic equation is changed into an autonomic four-dimensional dynamical system. The numerical time integrations of equations are performed by using the fourth-order Runge-Kutta method. And the Lyapunov exponent spectrum, the fractal dimension of strange attractors and the time evolution of deflection are obtained. The influence of geometry nonlinearity and viscoelastic parameter on the dynamic buckling of viscoelastic plates is discussed.