On the solution of bivariate population balance equations for aggregation: Pivotwise expansion of solution domain

The solution of a bivariate population balance equation (PBE) for aggregation of particles necessitates a large 2-d domain to be covered. A correspondingly large number of discretized equations for particle populations on pivots (representative sizes for bins) are solved, although at the end only a relatively small number of pivots are found to participate in the evolution process. In the present work, we initiate solution of the governing PBE on a small set of pivots that can represent the initial size distribution. New pivots are added to expand the computational domain in directions in which the evolving size distribution advances. A self-sufficient set of rules is developed to automate the addition of pivots, taken from an underlying X-grid formed by intersection of the lines of constant composition and constant particle mass. In order to test the robustness of the rule-set, simulations carried out with pivotwise expansion of X-grid are compared with those obtained using sufficiently large fixed X-grids for a number of composition independent and composition dependent aggregation kernels and initial conditions. The two techniques lead to identical predictions, with the former requiring only a fraction of the computational effort. The rule-set automatically reduces aggregation of particles of same composition to a 1-d problem. A midway change in the direction of expansion of domain, effected by the addition of particles of different mean composition, is captured correctly by the rule-set. The evolving shape of a computational domain carries with it the signature of the aggregation process, which can be insightful in complex and time dependent aggregation conditions. (c) 2012 Elsevier Ltd. All rights reserved.

Microstructure and texture evolution during accumulative roll bonding of aluminium alloys AA2219/AA5086 composite laminates

Accumulative roll bonding of two aluminium alloys, AA2219 and AA5086 was carried out up to 8 passes. During the course of ARB, the deformation inhomogeneity between the two alloy layers results in interfacial instability after the 4th pass, necking of the AA5086 layers after the 6th pass and fracture along the necked regions after the 7th and 8th pass. The EBSD analysis shows deformation bands along the interfaces after 8 passes of ARB. The ARB-processed materials predominantly show characteristic deformation texture components. The weak texture after the 2nd pass results from the combination of a weakly-textured starting AA2219 layer and a strongly-textured starting AA5086 layer. A strong deformation texture forms due to the high imposed strain after a higher number of ARB passes. Subgrain formation and related shear banding induces copper/S components in the case of the small elongated grains, while planar slip leads to the formation of brass component in the large elongated grains.

An Accurate Fractional Period Delay Generation System

An all-digital technique is proposed for generating an accurate delay irrespective of the inaccuracies of a controllable delay line. A subsampling technique-based delay measurement unit (DMU) capable of measuring delays accurately for the full period range is used as the feedback element to build accurate fractional period delays based on input digital control bits. The proposed delay generation system periodically measures and corrects the error and maintains it at the minimum value without requiring any special calibration phase. Up to 40x improvement in accuracy is demonstrated for a commercial programmable delay generator chip. The time-precision trade-off feature of the DMU is utilized to reduce the locking time. Loop dynamics are adjusted to stabilize the delay after the minimum error is achieved, thus avoiding additional jitter. Measurement results from a high-end oscilloscope also validate the effectiveness of the proposed system in improving accuracy.

Practical fully three-dimensional reconstruction algorithms for diffuse optical tomography

We have developed an efficient fully three-dimensional (3D) reconstruction algorithm for diffuse optical tomography (DOT). The 3D DOT, a severely ill-posed problem, is tackled through a pseudodynamic (PD) approach wherein an ordinary differential equation representing the evolution of the solution on pseudotime is integrated that bypasses an explicit inversion of the associated, ill-conditioned system matrix. One of the most computationally expensive parts of the iterative DOT algorithm, the reevaluation of the Jacobian in each of the iterations, is avoided by using the adjoint-Broyden update formula to provide low rank updates to the Jacobian. In addition, wherever feasible, we have also made the algorithm efficient by integrating along the quadratic path provided by the perturbation equation containing the Hessian. These algorithms are then proven by reconstruction, using simulated and experimental data and verifying the PD results with those from the popular Gauss-Newton scheme. The major findings of this work are as follows: (i) the PD reconstructions are comparatively artifact free, providing superior absorption coefficient maps in terms of quantitative accuracy and contrast recovery; (ii) the scaling of computation time with the dimension of the measurement set is much less steep with the Jacobian update formula in place than without it; and (iii) an increase in the data dimension, even though it renders the reconstruction problem less ill conditioned and thus provides relatively artifact-free reconstructions, does not necessarily provide better contrast property recovery. For the latter, one should also take care to uniformly distribute the measurement points, avoiding regions close to the source so that the relative strength of the derivatives for measurements away from the source does not become insignificant. (c) 2012 Optical Society of America

Asymmetric and symmetric rolling of magnesium: Evolution of microstructure, texture and mechanical properties

In the present study, asymmetric rolling was carried out for incorporating a shear component during the rolling at different temperatures, and was compared with conventional (symmetric) rolling. The microstructures were investigated using electron back-scatter diffraction (EBSD). The strain incorporated was compared with the help of grain orientation spread (GOS). GOS was eventually used as a criterion to partition the microstructure for separating the deformed and the dynamically recrystallized (DRX) grains. The texture of the partitioned DRX grains was shifted by similar to 30 degrees along the c-axis from the deformed grains. The mechanism of dynamic recrystallization (DRX) has been identified as continuous dynamic recovery and recrystallization (CDRR). The partitioned deformed grains for the higher temperature rolled specimens exhibited a texture similar to the room temperature rolled specimen. The asymmetric rolling introduces a shear component which shifts the texture fibre by similar to 5-10 degrees from the conventional rolling texture. This led to an increase in ductility with little compromise on strength. (c) 2012 Elsevier B.V. All rights reserved.

Effect of addition of aluminum on the evolution of microstructure in HITPERM class Fe44Co44Zr7B4Cu1 alloy

The paper reports the effect of the addition of small amount of Al on the microstructure and properties of HITPERM class rapidly solidified Fe44Co44Zr7B4Cu1 glassy alloy. Using three dimensional atom probe measurements we present evidence for the formation of Cu clusters on annealing in the metallic glass matrix of the Al containing alloy Fe43Co43Al2Zr7B4Cu1. Such clusters are otherwise absent in the parent alloy under similar conditions. The Cu clusters provides heterogeneous nucleation sites for the formation of bcc alpha'-FeCo phase leading to an increase in number density of this nanocrystalline phase and thereby enhancing the magnetic properties. (C) 2012 Elsevier B.V. All rights reserved.

An experimental study on cracking evolution in concrete and cement mortar by the b-value analysis of acoustic emission technique

Notched three point bend (TPB) specimens made with plain concrete and cement mortar were tested under crack mouth opening displacement (CMOD) control at a rate of 0.0004 mm/s and simultaneously acoustic emissions (AE) released were recorded during the experiments. Amplitude distribution analysis of AE released during concrete was carried out to study the development of fracture process in concrete and mortar specimens. The slope of the log-linear frequency-amplitude distribution of AE is known as the AE based b-value. The AE based b-value was computed in terms of physical process of time varying applied load using cumulative frequency distribution (Gutenberg-Richter relationship) and discrete frequency distribution (Aki's method) of AE released during concrete fracture. AE characteristics of plain concrete and cement mortar were studied and discussed and it was observed that the AE based b-value analysis serves as a tool to identify the damage in concrete structural members. (C) 2012 Elsevier Ltd. All rights reserved.

Gibbs Energy of Formation of Ca3Ti8Al12O37 and Phase Relations and Chemical Potentials in the System Al2O3-TiO2-CaO

The quaternary oxide in the system Al2O3-CaO-TiO2 is found to have the composition Ca3Ti8Al12O37 rather than CaTi3Al8O19 as reported in the literature. The standard Gibbs energy of formation of Ca3Ti8Al12O37 from component binary oxides is measured in the temperature range from 900 to 1250 K using a solid-state electrochemical cell incorporating single crystal CaF2 as the solid electrolyte. The results can be represented by the equation: delta G(f(0x))(0) (+/- 70)/J mol(-1) = -248474 - 15.706(T/K). Combining this information with thermodynamic data on calcium aluminates and titanates available in the literature, subsolidus phase relations in the pseudo-ternary system Al2O3-CaO-TiO2 are computed and presented as isothermal sections. The evolution of phase relations with temperature is highlighted. Chemical potential diagrams are computed at 1200 K, showing the stability domains of the various phases in the chemical potential-composition space. In each chemical potential diagram, chemical potential of one component is plotted against the cationic fraction of the other two components. The diagrams are valid at relatively high oxygen potentials where Ti is present in its four-valent state in all the oxide phases.

Evolution of quantum discord and its stability in two-qubit NMR systems

We investigate evolution of quantum correlations in ensembles of two-qubit nuclear spin systems via nuclear magnetic resonance techniques. We use discord as a measure of quantum correlations and the Werner state as an explicit example. We, first, introduce different ways of measuring discord and geometric discord in two-qubit systems and then describe the following experimental studies: (a) We quantitatively measure discord for Werner-like states prepared using an entangling pulse sequence. An initial thermal state with zero discord is gradually and periodically transformed into a mixed state with maximum discord. The experimental and simulated behavior of rise and fall of discord agree fairly well. (b) We examine the efficiency of dynamical decoupling sequences in preserving quantum correlations. In our experimental setup, the dynamical decoupling sequences preserved the traceless parts of the density matrices at high fidelity. But they could not maintain the purity of the quantum states and so were unable to keep the discord from decaying. (c) We observe the evolution of discord for a singlet-triplet mixed state during a radio-frequency spin-lock. A simple relaxation model describes the evolution of discord, and the accompanying evolution of fidelity of the long-lived singlet state, reasonably well.

Computational algorithms inspired by biological processes and evolution

In recent times computational algorithms inspired by biological processes and evolution are gaining much popularity for solving science and engineering problems. These algorithms are broadly classified into evolutionary computation and swarm intelligence algorithms, which are derived based on the analogy of natural evolution and biological activities. These include genetic algorithms, genetic programming, differential evolution, particle swarm optimization, ant colony optimization, artificial neural networks, etc. The algorithms being random-search techniques, use some heuristics to guide the search towards optimal solution and speed-up the convergence to obtain the global optimal solutions. The bio-inspired methods have several attractive features and advantages compared to conventional optimization solvers. They also facilitate the advantage of simulation and optimization environment simultaneously to solve hard-to-define (in simple expressions), real-world problems. These biologically inspired methods have provided novel ways of problem-solving for practical problems in traffic routing, networking, games, industry, robotics, economics, mechanical, chemical, electrical, civil, water resources and others fields. This article discusses the key features and development of bio-inspired computational algorithms, and their scope for application in science and engineering fields.

Microstructure and Texture Evolution in Interstitial-free (IF) Steel processed by Multi-Axial Forging

In this study, severe plastic deformation (SPD) of Ti-bearing interstitial-free steel was carried out by multi-axial forging (MAF) technique. The grain refinement achieved was comparable to that by other SPD techniques. A considerable heterogeneity was observed in the microstructure and texture. Texture of multi-axially forged steels has been evaluated and reported for the first time. The material exhibited a six-fold increase in the yield strength after four cycles of MAF.

Cannibalism can drive the evolution of behavioural phase polyphenism in locusts

During outbreaks, locust swarms can contain millions of insects travelling thousands of kilometers while devastating vegetation and crops. Such large-scale spatial organization is preceded locally by a dramatic density-dependent phenotypic transition in multiple traits. Behaviourally, low-density solitarious individuals avoid contact with one another; above a critical local density, they undergo a rapid behavioural transition to the gregarious phase whereby they exhibit mutual attraction. Although proximate causes of this phase polyphenism have been widely studied, the ultimate driving factors remain unclear. Using an individual-based evolutionary model, we reveal that cannibalism, a striking feature of locust ecology, could lead to the evolution of density-dependent behavioural phase-change in juvenile locusts. We show that this behavioural strategy minimizes risk associated with cannibalistic interactions and may account for the empirically observed persistence of locust groups during outbreaks. Our results provide a parsimonious explanation for the evolution of behavioural plasticity in locusts.

Markov chain modeling of evolution of strains in reinforced concrete flexural beams

From the analysis of experimentally observed variations in surface strains with loading in reinforced concrete beams, it is noted that there is a need to consider the evolution of strains (with loading) as a stochastic process. Use of Markov Chains for modeling stochastic evolution of strains with loading in reinforced concrete flexural beams is studied in this paper. A simple, yet practically useful, bi-level homogeneous Gaussian Markov Chain (BLHGMC) model is proposed for determining the state of strain in reinforced concrete beams. The BLHGMC model will be useful for predicting behavior/response of reinforced concrete beams leading to more rational design.

Effect of isothermal aging and electromigration on the microstructural evolution of solder interconnections during thermomechanical loading

The effect of different pre-aging treatments on the microstructural evolution of lead-free solder and growth of interfacial intermetallic compound layers under thermal cycling has been investigated in this work. The results show that there are distinct differences in the microstructural changes between samples with no pretreatment, samples that have experienced thermal annealing at 125A degrees C for 750 h before thermal cycling, and those that have had direct current (DC) stressing for 750 h as pretreatment. The microstructural evolution of the solder matrix is rationalized by utilizing the science of microstructures and analysis of the influence of electron flow on the precipitation phenomena. The finite-element method is utilized to understand the loading conditions imposed on the solder interconnections during cyclic stressing. The growth of intermetallic reaction layers is further analyzed by utilizing quantitative thermodynamic calculations coupled with kinetic analysis. The latter is based on the changes in the intrinsic diffusion fluxes of elements induced by current flow and alloying elements present in the system. With this concurrent approach the differences seen in thermal cycling behavior between the different pre-aging treatments can be explained.

A finite population model of mlecular evolution: theory and computation

This article is concerned with the evolution of haploid organisms that reproduce asexually. In a seminal piece of work, Eigen and coauthors proposed the quasispecies model in an attempt to understand such an evolutionary process. Their work has impacted antiviral treatment and vaccine design strategies. Yet, predictions of the quasispecies model are at best viewed as a guideline, primarily because it assumes an infinite population size, whereas realistic population sizes can be quite small. In this paper we consider a population genetics-based model aimed at understanding the evolution of such organisms with finite population sizes and present a rigorous study of the convergence and computational issues that arise therein. Our first result is structural and shows that, at any time during the evolution, as the population size tends to infinity, the distribution of genomes predicted by our model converges to that predicted by the quasispecies model. This justifies the continued use of the quasispecies model to derive guidelines for intervention. While the stationary state in the quasispecies model is readily obtained, due to the explosion of the state space in our model, exact computations are prohibitive. Our second set of results are computational in nature and address this issue. We derive conditions on the parameters of evolution under which our stochastic model mixes rapidly. Further, for a class of widely used fitness landscapes we give a fast deterministic algorithm which computes the stationary distribution of our model. These computational tools are expected to serve as a framework for the modeling of strategies for the deployment of mutagenic drugs.