Numerical study of the impact of inoculant and grain transport on macrosegregation and microstructure formation during solidification of an Al-22%Cu alloy

We investigate the impact of the nucleation law for nucleation on Al-Ti-B inoculant particles, of the motion of inoculant particles and of the motion of grains on the predicted macrosegregation and microstructure in a grain-refined Al-22 wt.% Cu alloy casting. We conduct the study by numerical simulations of a casting experiment in a side-cooled 76×76×254 mm sand mould. Macrosegregation and microstructure formation are studied with a volume-averaged two-phase model accounting for macroscopic heat and solute transport, melt convection, and transport of inoculant particles and equiaxed grains. On the microscopic scale it accounts for nucleation on inoculant particles with a given size distribution (and corresponding activation undercooling distribution)and for the growth of globular solid grains. The growth kinetics is described by accounting for limited solute diffusion in both liquid and solid phases and for convective effects. We show that the consideration of a size distribution of the inoculants has a strong impact on the microstructure(final grain size) prediction. The transport of inoculants significantly increases the microstructure heterogeneities and the grain motion refines the microstructure and reduces the microstructure heterogeneities.

Simulation of inhomogeneous distributions of ultracold atoms in an optical lattice via a massively parallel implementation of nonequilibrium strong-coupling perturbation theory

We present a nonequilibrium strong-coupling approach to inhomogeneous systems of ultracold atoms in optical lattices. We demonstrate its application to the Mott-insulating phase of a two-dimensional Fermi-Hubbard model in the presence of a trap potential. Since the theory is formulated self-consistently, the numerical implementation relies on a massively parallel evaluation of the self-energy and the Green's function at each lattice site, employing thousands of CPUs. While the computation of the self-energy is straightforward to parallelize, the evaluation of the Green's function requires the inversion of a large sparse 10(d) x 10(d) matrix, with d > 6. As a crucial ingredient, our solution heavily relies on the smallness of the hopping as compared to the interaction strength and yields a widely scalable realization of a rapidly converging iterative algorithm which evaluates all elements of the Green's function. Results are validated by comparing with the homogeneous case via the local-density approximation. These calculations also show that the local-density approximation is valid in nonequilibrium setups without mass transport.

Large deviation function and fluctuation theorem for classical particle transport

We analytically evaluate the large deviation function in a simple model of classical particle transfer between two reservoirs. We illustrate how the asymptotic long-time regime is reached starting from a special propagating initial condition. We show that the steady-state fluctuation theorem holds provided that the distribution of the particle number decays faster than an exponential, implying analyticity of the generating function and a discrete spectrum for its evolution operator.

Electronic band structure and Fermi surfaces of the quasi-two-dimensional monophosphate tungsten bronze, P4W12O44

The electronic structure of quasi-two-dimensional monophosphate tungsten bronze, P4W12O44, has been investigated by high-resolution angle-resolved photoemission spectroscopy and density functional theoretical calculations. Experimental electron-like bands around Gamma point and Fermi surfaces have similar shapes as predicted by calculations. Fermi surface mapping at different temperatures shows a depletion of density of states at low temperature in certain flat portions of the Fermi surfaces. These flat portions of the Fermi surfaces satisfy the partial nesting condition with incommensurate nesting vectors q(1) and q(2), which leads to the formation of charge density waves in this phosphate tungsten bronzes. The setting up of charge density wave in these bronzes can well explain the anomaly observed in its transport properties. Copyright (C) EPLA, 2014

IS A DEEP ONE-CELL MERIDIONAL CIRCULATION ESSENTIAL FOR THE FLUX TRANSPORT SOLAR DYNAMO?

The solar activity cycle is successfully modeled by the flux transport dynamo, in which the meridional circulation of the Sun plays an important role. Most of the kinematic dynamo simulations assume a one-cell structure of the meridional circulation within the convection zone, with the equatorward return flow at its bottom. In view of the recent claims that the return flow occurs at a much shallower depth, we explore whether a meridional circulation with such a shallow return flow can still retain the attractive features of the flux transport dynamo (such as a proper butterfly diagram, the proper phase relation between the toroidal and poloidal fields). We consider additional cells of the meridional circulation below the shallow return flow-both the case of multiple cells radially stacked above one another and the case of more complicated cell patterns. As long as there is an equatorward flow in low latitudes at the bottom of the convection zone, we find that the solar behavior is approximately reproduced. However, if there is either no flow or a poleward flow at the bottom of the convection zone, then we cannot reproduce solar behavior. On making the turbulent diffusivity low, we still find periodic behavior, although the period of the cycle becomes unrealistically large. In addition, with a low diffusivity, we do not get the observed correlation between the polar field at the sunspot minimum and the strength of the next cycle, which is reproduced when diffusivity is high. On introducing radially downward pumping, we get a more reasonable period and more solar-like behavior even with low diffusivity.

Diffusion in an elastic medium: A model for macromolecule transport across the nuclear pore complex

Nuclear pore complexes (NPCs) are very selective filters that sit on the membrane of the nucleus and monitor the transport between the cytoplasm and the nucleoplasm. For the central plug of NPC two models have been suggested in the literature. The first suggests that the plug is a reversible hydrogel while the other suggests that it is a polymer brush. Here we propose a model for the transport of a protein through the plug, which is general enough to cover both the models. The protein stretches the plug and creates a local deformation, which together with the protein, we refer to as the bubble. We start with the free energy for creation of the bubble and consider its motion within the plug. The relevant coordinate is the center of the bubble which executes random walk. We find that for faster relaxation of the gel, the diffusion of the bubble is greater. (C) 2014 Elsevier-B.V. All rights reserved.

Study of liquid fuel transport in a small carburetted engine in the context of cold-start HC emission control

In the present study, a detailed visualization of the transport of fuel film has been performed in a small carburetted engine with a transparent manifold at the exit of the carburettor. The presence of fuel film is observed significantly on the lower half of the manifold at idling, while at load conditions, the film is found to be distributed all throughout the manifold walls. Quantitative measurement of the fuel film in a specially-designed manifold of square cross section has also been performed using the planar laser-induced fluorescence (PLIF) technique. The measured fuel film thickness is observed to be of the order of 1 nun at idling, and in the range of 0.1 to 0.4 mm over the range of load and speed studied. These engine studies are complemented by experiments conducted in a carburettor rig to study the state of the fuel exiting the carburettor. Laser-based Particle/Droplet Image Analysis (PDIA) technique is used to identify fuel droplets and ligaments and estimate droplet diameters. At a throttle position corresponding to idling, the fuel exiting the carburettor is found to consist of very fine droplets of size less than 15 mu m and large fuel ligaments associated with length scales of the order of 500 mu m and higher. For a constant pressure difference across the carburettor, the fuel consists of droplets with an SMD of the order of 30 mu m. Also, the effect of liquid fuel film on the cold start HC emissions is studied. Based on the understanding obtained from these studies, strategies such as manifold heating and varying carburettor main jet nozzle diameter are implemented. These are observed to reduce emissions under both idling and varying load conditions.

Monte Carlo simulation of light transport in turbid medium with embedded object-spherical, cylindrical, ellipsoidal, or cuboidal objects embedded within multilayered tissues

Monte Carlo modeling of light transport in multilayered tissue (MCML) is modified to incorporate objects of various shapes (sphere, ellipsoid, cylinder, or cuboid) with a refractive-index mismatched boundary. These geometries would be useful for modeling lymph nodes, tumors, blood vessels, capillaries, bones, the head, and other body parts. Mesh-based Monte Carlo (MMC) has also been used to compare the results from the MCML with embedded objects (MCML-EO). Our simulation assumes a realistic tissue model and can also handle the transmission/reflection at the object-tissue boundary due to the mismatch of the refractive index. Simulation of MCML-EO takes a few seconds, whereas MMC takes nearly an hour for the same geometry and optical properties. Contour plots of fluence distribution from MCML-EO and MMC correlate well. This study assists one to decide on the tool to use for modeling light propagation in biological tissue with objects of regular shapes embedded in it. For irregular inhomogeneity in the model (tissue), MMC has to be used. If the embedded objects (inhomogeneity) are of regular geometry (shapes), then MCML-EO is a better option, as simulations like Raman scattering, fluorescent imaging, and optical coherence tomography are currently possible only with MCML. (C) 2014 Society of Photo-Optical Instrumentation Engineers (SPIE)

Selenium in Diketopyrrolopyrrole-based Polymers: Influence on Electronic Properties and Charge Carrier Mobilities

Diketopyrrolopyrrole (DPP)-based pi-conjugated copolymers with thiophene have exceptionally high electron mobilities. This paper investigates electronic properties and charge carrier mobilities of selenophene containing analogues. Two new copolymers, with alternating thiophene DPP (TDPP) and selenophene DPP (SeDPP) units, were synthesized. Two side-chains, hexyl (Hex) and triethylene glycol (TEG) were employed, yielding polymers designated as PTDPPSeDPP-Hex and PTDPPSeDPP-TEG. Selenophene systems have smaller band gaps, with concomitant enhancement of the stability of the reduced state. For both polymers, ambipolar mobilities were observed in organic field-effect transistors (OFET). Grazing incidence X-ray diffraction (GIXD) data indicates preferential edge-on orientation of PTDPPSeDPP-TEG, which leads to superior charge transport properties of the TEG substituted polymer, as compared to its Hex analogue. Time-dependent-density functional theory (TDDFT) calculations corroborate the decrease in the optical band gap with the inclusion of selenophene. Ambipolar charge transport is rationalized by exceptionally wide conduction bands. Delta SCF calculations confirm the larger electron affinity, and therefore the greater stability, of the reduced form of the selenophene-containing DPP polymer in presence of chloroform.

A DYNAMO MODEL OF MAGNETIC ACTIVITY IN SOLAR-LIKE STARS WITH DIFFERENT ROTATIONAL VELOCITIES

We attempt to provide a quantitative theoretical explanation for the observations that Ca II H/K emission and X-ray emission from solar-like stars increase with decreasing Rossby number (i.e., with faster rotation). Assuming that these emissions are caused by magnetic cycles similar to the sunspot cycle, we construct flux transport dynamo models of 1M(circle dot) stars rotating with different rotation periods. We first compute the differential rotation and the meridional circulation inside these stars from a mean-field hydrodynamics model. Then these are substituted in our dynamo code to produce periodic solutions. We find that the dimensionless amplitude f(m) of the toroidal flux through the star increases with decreasing rotation period. The observational data can be matched if we assume the emissions to go as the power 3-4 of f(m). Assuming that the Babcock-Leighton mechanism saturates with increasing rotation, we can provide an explanation for the observed saturation of emission at low Rossby numbers. The main failure of our model is that it predicts an increase of the magnetic cycle period with increasing rotation rate, which is the opposite of what is found observationally. Much of our calculations are based on the assumption that the magnetic buoyancy makes the magnetic flux tubes rise radially from the bottom of the convection zone. Taking into account the fact that the Coriolis force diverts the magnetic flux tubes to rise parallel to the rotation axis in rapidly rotating stars, the results do not change qualitatively.

Virus-sized colloid transport in a single pore: Model development and sensitivity analysis

A mathematical model is developed to simulate the transport and deposition of virus-sized colloids in a cylindrical pore throat considering various processes such as advection, diffusion, colloid-collector surface interactions and hydrodynamic wall effects. The pore space is divided into three different regions, namely, bulk, diffusion and potential regions, based on the dominant processes acting in each of these regions. In the bulk region, colloid transport is governed by advection and diffusion whereas in the diffusion region, colloid mobility due to diffusion is retarded by hydrodynamic wall effects. Colloid-collector interaction forces dominate the transport in the potential region where colloid deposition occurs. The governing equations are non-dimensionalized and solved numerically. A sensitivity analysis indicates that the virus-sized colloid transport and deposition is significantly affected by various pore-scale parameters such as the surface potentials on colloid and collector, ionic strength of the solution, flow velocity, pore size and colloid size. The adsorbed concentration and hence, the favorability of the surface for adsorption increases with: (i) decreasing magnitude and ratio of surface potentials on colloid and collector, (ii) increasing ionic strength and (iii) increasing pore radius. The adsorbed concentration increases with increasing Pe, reaching a maximum value at Pe = 0.1 and then decreases thereafter. Also, the colloid size significantly affects particle deposition with the adsorbed concentration increasing with increasing particle radius, reaching a maximum value at a particle radius of 100 nm and then decreasing with increasing radius. System hydrodynamics is found to have a greater effect on larger particles than on smaller ones. The secondary minimum contribution to particle deposition has been found to increase as the favorability of the surface for adsorption decreases. The sensitivity of the model to a given parameter will be high if the conditions are favorable for adsorption. The results agree qualitatively with the column-scale experimental observations available in the literature. The current model forms the building block in upscaling colloid transport from pore scale to Darcy scale using Pore-Network Modeling. (C) 2014 Elsevier By. All rights reserved.

Modelling of optical transport behavior of organic photovoltaic devices with nano-pillar transparent conducting electrodes

Optical transport behavior of organic photo-voltaic devices with nano-pillar transparent electrodes is investigated in this paper in order to understand possible enhancement of their charge-collection efficiency. Modeling and simulations of optical transport due to this architecture show an interesting regime of length-scale dependent optical characteristics. An electromagnetic wave propagation model is employed with simulation objectives toward understanding the mechanism of optical scattering and waveguide effects due to the nano-pillars and effective transmission through the active layer. Partial filling of gaps between the nano-pillars due to the nano-fabrication process is taken into consideration. Observations made in this paper will facilitate appropriate design rules for nano-pillar electrodes. (C) 2014 AIP Publishing LLC.

Sensitivity of meteorological input and soil properties in simulating aerosols (dust, PM10, and BC) using CHIMERE chemistry transport model

The objective of this study is to evaluate the ability of a European chemistry transport model, `CHIMERE' driven by the US meteorological model MM5, in simulating aerosol concentrations dust, PM10 and black carbon (BC)] over the Indian region. An evaluation of a meteorological event (dust storm); impact of change in soil-related parameters and meteorological input grid resolution on these aerosol concentrations has been performed. Dust storm simulation over Indo-Gangetic basin indicates ability of the model to capture dust storm events. Measured (AERONET data) and simulated parameters such as aerosol optical depth (AOD) and Angstrom exponent are used to evaluate the performance of the model to capture the dust storm event. A sensitivity study is performed to investigate the impact of change in soil characteristics (thickness of the soil layer in contact with air, volumetric water, and air content of the soil) and meteorological input grid resolution on the aerosol (dust, PM10, BC) distribution. Results show that soil parameters and meteorological input grid resolution have an important impact on spatial distribution of aerosol (dust, PM10, BC) concentrations.

Improved quasiparticle wave functions and mean field for G(0)W(0) calculations: Initialization with the COHSEX operator

The GW approximation to the electron self-energy has become a standard method for ab initio calculation of excited-state properties of condensed-matter systems. In many calculations, the G W self-energy operator, E, is taken to be diagonal in the density functional theory (DFT) Kohn-Sham basis within the G0 W0 scheme. However, there are known situations in which this diagonal Go Wo approximation starting from DFT is inadequate. We present two schemes to resolve such problems. The first, which we called sc-COHSEX-PG W, involves construction of an improved mean field using the static limit of GW, known as COHSEX (Coulomb hole and screened exchange), which is significantly simpler to treat than GW W. In this scheme, frequency-dependent self energy E(N), is constructed and taken to be diagonal in the COHSEX orbitals after the system is solved self-consistently within this formalism. The second method is called off diagonal-COHSEX G W (od-COHSEX-PG W). In this method, one does not self-consistently change the mean-field starting point but diagonalizes the COHSEX Hamiltonian within the Kohn-Sham basis to obtain quasiparticle wave functions and uses the resulting orbitals to construct the G W E in the diagonal form. We apply both methods to a molecular system, silane, and to two bulk systems, Si and Ge under pressure. For silane, both methods give good quasiparticle wave functions and energies. Both methods give good band gaps for bulk silicon and maintain good agreement with experiment. Further, the sc-COHSEX-PGW method solves the qualitatively incorrect DFT mean-field starting point (having a band overlap) in bulk Ge under pressure.

Angular momentum transport and evolution of lopsided galaxies

The surface brightness distribution in the majority of stellar galactic discs falls off exponentially. Often what lies beyond such a stellar disc is the neutral hydrogen gas whose distribution also follows a nearly exponential profile at least for a number of nearby disc galaxies. Both the stars and gas are commonly known to host lopsided asymmetry especially in the outer parts of a galaxy. The role of such asymmetry in the dynamical evolution of a galaxy has not been explored so far. Following Lindblad's original idea of kinematic density waves, we show that the outer part of an exponential disc is ideally suitable for hosting lopsided asymmetry. Further, we compute the transport of angular momentum in the combined stars and gas disc embedded in a dark matter halo. We show that in a pure star and gas disc, there is a transition point where the free precession frequency of a lopsided mode, Omega - kappa, changes from retrograde to prograde and this in turn reverses the direction of angular momentum flow in the disc leading to an unphysical behaviour. We show that this problem is overcome in the presence of a dark matter halo, which sets the angular momentum flow outwards as required for disc evolution, provided the lopsidedness is leading in nature. This, plus the well-known angular momentum transport in the inner parts due to spiral arms, can facilitate an inflow of gas from outside perhaps through the cosmic filaments.