An Optimizing Code Generator for a Class of Lattice-Boltzmann Computations

The Lattice-Boltzmann method (LBM), a promising new particle-based simulation technique for complex and multiscale fluid flows, has seen tremendous adoption in recent years in computational fluid dynamics. Even with a state-of-the-art LBM solver such as Palabos, a user has to still manually write the program using library-supplied primitives. We propose an automated code generator for a class of LBM computations with the objective to achieve high performance on modern architectures. Few studies have looked at time tiling for LBM codes. We exploit a key similarity between stencils and LBM to enable polyhedral optimizations and in turn time tiling for LBM. We also characterize the performance of LBM with the Roofline performance model. Experimental results for standard LBM simulations like Lid Driven Cavity, Flow Past Cylinder, and Poiseuille Flow show that our scheme consistently outperforms Palabos-on average by up to 3x while running on 16 cores of an Intel Xeon (Sandybridge). We also obtain an improvement of 2.47x on the SPEC LBM benchmark.

Field induced domain switching as the origin of anomalous lattice strain along non-polar direction in rhombohedral BiScO3-PbTiO3 close to the morphotropic phase boundary

The lattice strain and domain switching behavior of xBiScO(3)-(1-x) PbTiO3 (x = 0.40) was investigated as a function of cyclic field and grain orientation by in situ X-ray diffraction during application of electric fields. The electric field induced 200 lattice strain was measured to be five times larger than the 111 lattice strain in pseudorhombohedral xBiScO(3)-(1-x) PbTiO3 (x = 0.40). It is shown that the anomalous 200 lattice strain is not an intrinsic phenomenon, but arises primarily due to stress associated with the reorientation of the 111 domains in dense polycrystalline ceramic. (C) 2015 AIP Publishing LLC.

A new class of high strength high temperature Cobalt based gamma-gamma ` Co-Mo-Al alloys stabilized with Ta addition

The present paper reports a new class of Co based superalloys that has gamma-gamma' microstructure and exhibits much lower density compared to other commercially available Co superalloys including Co-Al-W based alloys. The basic composition is Co-10Al-5Mo (at%) with addition of 2 at% Ta for stabilization of gamma' phase. The gamma-gamma' microstructure evolves through solutionising and aging treatment. Using first principles calculations, we observe that Ta plays a crucial role in stabilizing gamma' phase. By addition of Ta in the basic stoichiometric composition Co-3(Al, Mo), the enthalpy of formation (Delta H-f) of L1(2) structure (gamma' phase) becomes more negative in comparison to DO19 structure. The All of the L12 structure becomes further more negative by the occupancy of Ni and Ti atoms in the lattice suggesting an increase in the stability of the gamma' precipitates. Among large number of alloys studied experimentally, the paper presents results of detailed investigations on Co-10Al-5Mo-2Ta, Co-30Ni-10Al-5Mo-2Ta and Co-30Ni-10Al-5Mo-2Ta-2Ti. To evaluate the role alloying elements, atom probe tomography investigations were carried out to obtain partition coefficients for the constituent elements. The results show strong partitioning of Ni, Al, Ta and Ti in ordered gamma' precipitates. 2015 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Anomalous Dielectric Features and Structural Implications in Gallic Acid-lsoniazid Cocrystal Mediated by Lattice Water

An organic supramolecular ternary salt (gallic acid:isoniazid:water; GINZH) examined earlier for its proton conducting characteristics is observed to display step-like dielectric behavior across the structural phase transition mediated by loss of water of hydration at 389 K. The presence of hydration in the crystal lattice along with proton mobility between acid base pairs controls the ``ferroelectric like'' behavior until the phase transition temperature.

Attenuating microwave radiation by absorption through controlled nanoparticle localization in PC/PVDF blends

Nanoscale ordering in a polymer blend structure is indispensable to obtain materials with tailored properties. It was established here that controlling the arrangement of nanoparticles, with different characteristics, in co-continuous PC/PVDF (polycarbonate/poly(vinylidene fluoride)) blends can result in outstanding microwave absorption (ca. 90%). An excellent reflection loss (RL) of ca. -71 dB was obtained for a model blend structure wherein the conducting (multiwall carbon nanotubes, MWNTs) and the magnetic inclusions (Fe3O4) are localized in PVDF and the dielectric inclusion (barium titanate, BT) is in PC. The MWNTs were modified using polyaniline, which facilitates better charge transport in the blends. Furthermore, by introducing surface active groups on BT nanoparticles and changing the macroscopic processing conditions, the localization of BT nanoparticles can be tailored, otherwise BT nanoparticles would localize in the preferred phase (PVDF). In this study, we have shown that by ordered arrangement of nanoparticles, the incoming EM radiation can be attenuated. For instance, when PANI-MWNTs were localized in PVDF, the shielding was mainly through reflection. Now by localizing the conducting inclusion and the magnetic lossy materials in PVDF and the dielectric materials in PC, an outstanding shielding effectiveness of ca. -37 dB was achieved where shielding was mainly through absorption (ca. 90%). Thus, this study clearly demonstrates that lightweight microwave absorbers can be designed using polymer blends as a tool.

Ordered Pd2Ge Intermetallic Nanoparticles as Highly Efficient and Robust Catalyst for Ethanol Oxidation

Pd2Ge nanoparticles were synthesized by superhydride reduction of K2PdCl4 and GeCl4. The syntheses were performed using a solvothermal method in the absence of surfactants, and the size of the nanoparticles was controlled by varying the reaction time. The powder X-ray diffraction (PXRD) and transmission electron microscopy data suggest that Pd2Ge nanoparticles were formed as an ordered intermetallic phase. In the crystal structure, Pd and Ge atoms occupy two different crystallographic positions with a vacancy in one of the Ge sites, which was proved by PXRD and energy-dispersive X-ray analysis. The catalyst is highly efficient for the electrochemical oxidation of ethanol and is stable up to the 250th cycle in alkaline medium. The electrochemical active surface area and current density values obtained, 1.41 cm(2) and 4.1 mA cm(-2), respectively, are superior to those of the commercial Pd on carbon. The experimentally observed data were interpreted in terms of the combined effect of adsorption energies of CH3CO and OH radical, d-band center model, and work function of the corresponding catalyst surfaces.

Effect of Super-exchange Interaction on Ground state Magnetic Properties of Spin-dependent Falicov-Kimball Model on a Triangular Lattice

Ground state magnetic properties are studied by incorporating the super-exchange interaction (J(se)) in the spin-dependent Falicov-Kimball model (FKM) between localized (f-) electrons on a triangular lattice for half filled case. Numerical diagonalization and Monte-Carlo simulation are used to study the ground state magnetic properties. We have found that the magnetic moment of (d-) and (f-) electrons strongly depend on the value of Hund's exchange (J), super-exchange interaction (J(se)) and also depends on the number of (d-) electrons (N-d). The ground state changes from antiferromagnetic (AFM) to ferromagnetic (FM) state as we decrease (N-d). Also the density of d electrons at each site depends on the value of J and J(se).

The role of intra- and inter-site exchange correlations in the extended Falicov-Kimball model on a triangular lattice

Ground state magnetic properties of the spin-dependent Falicov-Kimball model (FKM) are studied by incorporating the intrasite exchange correlation J (between itinerant d- and localized f-electrons) and intersite (superexchange) correlation J (between localized f-electrons) on a triangular lattice for two different fillings. Numerical diagonalization and Monte-Carlo techniques are used to determine the ground state magnetic properties. Transitions from antiferromagnetic to ferromagnetic and again to re-entrant antiferromagnetic phase is observed in a wide range of parameter space. The magnetic moments of d- and f-electrons are observed to depend strongly on the value off, J and also on the total number of d-electrons (N-d). (C) 2015 Elsevier Ltd. All rights reserved.

Groups constrain the use of risky habitat by individuals: a new cost to sociality?

Predation risk can strongly constrain how individuals use time and space. Grouping is known to reduce an individual's time investment in costly antipredator behaviours. Whether grouping might similarly provide a spatial release from antipredator behaviour and allow individuals to use risky habitat more and, thus, improve their access to resources is poorly known. We used mosquito larvae, Aedes aegypti, to test the hypothesis that grouping facilitates the use of high-risk habitat. We provided two habitats, one darker, low-risk and one lighter, high-risk, and measured the relative time spent in the latter by solitary larvae versus larvae in small groups. We tested larvae reared under different resource levels, and thus presumed to vary in body condition, because condition is known to influence risk taking. We also varied the degree of contrast in habitat structure. We predicted that individuals in groups should use high-risk habitat more than solitary individuals allowing for influences of body condition and contrast in habitat structure. Grouping strongly influenced the time spent in the high-risk habitat, but, contrary to our expectation, individuals in groups spent less time in the high-risk habitat than solitary individuals. Furthermore, solitary individuals considerably increased the proportion of time spent in the high-risk habitat over time, whereas individuals in groups did not. Both solitary individuals and those in groups showed a small increase over time in their use of riskier locations within each habitat. The differences between solitary individuals and those in groups held across all resource and contrast conditions. Grouping may, thus, carry a poorly understood cost of constraining habitat use. This cost may arise because movement traits important for maintaining group cohesion (a result of strong selection on grouping) can act to exaggerate an individual preference for low-risk habitat. Further research is needed to examine the interplay between grouping, individual movement and habitat use traits in environments heterogeneous in risk and resources. (C) 2015 The Association for the Study of Animal Behaviour. Published by Elsevier Ltd. All rights reserved.

NH3 adsorption on PtM (Fe, Co, Ni) surfaces: Cooperating effects of charge transfer, magnetic ordering and lattice strain

Adsorption of a molecule or group with an atom which is less electronegative than oxygen (0) and directly interacting with the surface is very relevant to development of PtM (M = 3d-transition metal) catalysts with high activity. Here, we present theoretical analysis of the adsorption of NH3 molecule (N being less electronegative than 0) on (111) surfaces of PtM (Fe, Co, Ni) alloys using the first principles density functional approach. We find that, while NH3-Pt interaction is stronger than that of NH3 with the elemental M-surfaces, it is weaker than the strength of interaction of NH3 with M-site on the surface of PtM alloy. (C) 2016 Published by Elsevier B.V.

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.

Social optimum in social groups with give-and-take criterion

We consider a Social Group' of networked nodes, seeking a universe' of segments. Each node has a subset of the universe and access to an expensive resource for downloading data. Nodes can also acquire the universe by exchanging copies of segments among themselves, at low cost, using inter-node links. While exchanges over inter-node links ensure minimum cost, some nodes in the group try to exploit the system. We term such nodes as non-reciprocating nodes' and prohibit such behavior by proposing the give-and-take' criterion, where exchange is allowed if each node has segments unavailable with the other. Under this criterion, we consider the problem of maximizing the number of nodes with the universe at the end of local exchanges. First, we present a randomized algorithm that is shown to be optimal in the asymptotic regime. Then, we present greedy links algorithm, which performs well for most of the scenarios and yields an optimal result when the number of nodes is four. The polygon algorithm is proposed, which yields an optimal result when each of the nodes has a unique segment. After presenting some intuitive algorithms (e.g., greedy incremental algorithm and rarest first algorithm), we compare the performances of all proposed algorithms with the optimal. Copyright (c) 2015 John Wiley & Sons, Ltd.

Traveling waves in trimer granular lattice I: Bifurcation structure of traveling waves in the unit-cell model

Present paper is the first one in the series devoted to the dynamics of traveling waves emerging in the uncompressed, tri-atomic granular crystals. This work is primarily concerned with the dynamics of one-dimensional periodic granular trimer (tri-atomic) chains in the state of acoustic vacuum. Each unit cell consists of three spherical particles of different masses subject to periodic boundary conditions. Hertzian interaction law governs the mutual interaction of these particles. Under the assumption of zero pre-compression, this interaction is modeled as purely nonlinear, which means the absence of linear force component. The dynamics of such chains is governed by the two system parameters that scale the mass ratios between the particles of the unit cell. Such a system supports two different classes of periodic solutions namely the traveling and standing waves. The primary objective of the present study is the numerical analysis of the bifurcation structure of these solutions with emphasis on the dynamics of traveling waves. In fact, understanding of the bifurcation structure of the traveling wave solutions emerging in the unit-cell granular trimer is rather important and can shed light on the more complex nonlinear wave phenomena emerging in semi-infinite trimer chains. (c) 2016 Elsevier B.V. All rights reserved.