Topological blocking in quantum quench dynamics

We study the nonequilibrium dynamics of quenching through a quantum critical point in topological systems, focusing on one of their defining features: ground-state degeneracies and associated topological sectors. We present the notion of ``topological blocking,'' experienced by the dynamics due to a mismatch in degeneracies between two phases, and we argue that the dynamic evolution of the quench depends strongly on the topological sector being probed. We demonstrate this interplay between quench and topology in models stemming from two extensively studied systems, the transverse Ising chain and the Kitaev honeycomb model. Through nonlocal maps of each of these systems, we effectively study spinless fermionic p-wave paired topological superconductors. Confining the systems to ring and toroidal geometries, respectively, enables us to cleanly address degeneracies, subtle issues of fermion occupation and parity, and mismatches between topological sectors. We show that various features of the quench, which are related to Kibble-Zurek physics, are sensitive to the topological sector being probed, in particular, the overlap between the time-evolved initial ground state and an appropriate low-energy state of the final Hamiltonian. While most of our study is confined to translationally invariant systems, where momentum is a convenient quantum number, we briefly consider the effect of disorder and illustrate how this can influence the quench in a qualitatively different way depending on the topological sector considered.

Edge states of a three-dimensional topological insulator

We use the bulk Hamiltonian for a three-dimensional topological insulator such as Bi-2 Se-3 to study the states which appear on its various surfaces and along the edge between two surfaces. We use both analytical methods based on the surface Hamiltonians (which are derived from the bulk Hamiltonian) and numerical methods based on a lattice discretization of the bulk Hamiltonian. We find that the application of a potential barrier along an edge can give rise to states localized at that edge. These states have an unusual energy-momentum dispersion which can be controlled by applying a potential along the edge; in particular, the velocity of these states can be tuned to zero. The scattering and conductance across the edge is studied as a function of the edge potential. We show that a magnetic field in a particular direction can also give rise to zero energy states on certain edges. We point out possible experimental ways of looking for the various edge states.

Dense shallow granular flows

Simplified equations are derived for a granular flow in the `dense' limit where the volume fraction is close to that for dynamical arrest, and the `shallow' limit where the stream-wise length for flow development (L) is large compared with the cross-stream height (h). The mass and diameter of the particles are set equal to 1 in the analysis without loss of generality. In the dense limit, the equations are simplified by taking advantage of the power-law divergence of the pair distribution function chi proportional to (phi(ad) - phi)(-alpha), and a faster divergence of the derivativ rho(d chi/d rho) similar to (d chi/d phi), where rho and phi are the density and volume fraction, and phi(ad) is the volume fraction for arrested dynamics. When the height h is much larger than the conduction length, the energy equation reduces to an algebraic balance between the rates of production and dissipation of energy, and the stress is proportional to the square of the strain rate (Bagnold law). In the shallow limit, the stress reduces to a simplified Bagnold stress, where all components of the stress are proportional to (partial derivative u(x)/partial derivative y)(2), which is the cross-stream (y) derivative of the stream-wise (x) velocity. In the simplified equations for dense shallow flows, the inertial terms are neglected in the y momentum equation in the shallow limit because the are O(h/L) smaller than the divergence of the stress. The resulting model contains two equations, a mass conservation equations which reduces to a solenoidal condition on the velocity in the incompressible limit, and a stream-wise momentum equation which contains just one parameter B which is a combination of the Bagnold coefficients and their derivatives with respect to volume fraction. The leading-order dense shallow flow equations, as well as the first correction due to density variations, are analysed for two representative flows. The first is the development from a plug flow to a fully developed Bagnold profile for the flow down an inclined plane. The analysis shows that the flow development length is ((rho) over barh(3)/B) , where (rho) over bar is the mean density, and this length is numerically estimated from previous simulation results. The second example is the development of the boundary layer at the base of the flow when a plug flow (with a slip condition at the base) encounters a rough base, in the limit where the momentum boundary layer thickness is small compared with the flow height. Analytical solutions can be found only when the stream-wise velocity far from the surface varies as x(F), where x is the stream-wise distance from the start of the rough base and F is an exponent. The boundary layer thickness increases as (l(2)x)(1/3) for all values of F, where the length scale l = root 2B/(rho) over bar. The analysis reveals important differences between granular flows and the flows of Newtonian fluids. The Reynolds number (ratio of inertial and viscous terms) turns out to depend only on the layer height and Bagnold coefficients, and is independent of the flow velocity, because both the inertial terms in the conservation equations and the divergence of the stress depend on the square of the velocity/velocity gradients. The compressibility number (ratio of the variation in volume fraction and mean volume fraction) is independent of the flow velocity and layer height, and depends only on the volume fraction and Bagnold coefficients.

The generalized Onsager model for a binary gas mixture

The Onsager model for the secondary flow field in a high-speed rotating cylinder is extended to incorporate the difference in mass of the two species in a binary gas mixture. The base flow is an isothermal solid-body rotation in which there is a balance between the radial pressure gradient and the centrifugal force density for each species. Explicit expressions for the radial variation of the pressure, mass/mole fractions, and from these the radial variation of the viscosity, thermal conductivity and diffusion coefficient, are derived, and these are used in the computation of the secondary flow. For the secondary flow, the mass, momentum and energy equations in axisymmetric coordinates are expanded in an asymptotic series in a parameter epsilon = (Delta m/m(av)), where Delta m is the difference in the molecular masses of the two species, and the average molecular mass m(av) is defined as m(av) = (rho(w1)m(1) + rho(w2)m(2))/rho(w), where rho(w1) and rho(w2) are the mass densities of the two species at the wall, and rho(w) = rho(w1) + rho(w2). The equation for the master potential and the boundary conditions are derived correct to O(epsilon(2)). The leading-order equation for the master potential contains a self-adjoint sixth-order operator in the radial direction, which is different from the generalized Onsager model (Pradhan & Kumaran, J. Fluid Mech., vol. 686, 2011, pp. 109-159), since the species mass difference is included in the computation of the density, viscosity and thermal conductivity in the base state. This is solved, subject to boundary conditions, to obtain the leading approximation for the secondary flow, followed by a solution of the diffusion equation for the leading correction to the species mole fractions. The O(epsilon) and O(epsilon(2)) equations contain inhomogeneous terms that depend on the lower-order solutions, and these are solved in a hierarchical manner to obtain the O(epsilon) and O(epsilon(2)) corrections to the master potential. A similar hierarchical procedure is used for the Carrier-Maslen model for the end-cap secondary flow. The results of the Onsager hierarchy, up to O(epsilon(2)), are compared with the results of direct simulation Monte Carlo simulations for a binary hard-sphere gas mixture for secondary flow due to a wall temperature gradient, inflow/outflow of gas along the axis, as well as mass and momentum sources in the flow. There is excellent agreement between the solutions for the secondary flow correct to O(epsilon(2)) and the simulations, to within 15 %, even at a Reynolds number as low as 100, and length/diameter ratio as low as 2, for a low stratification parameter A of 0.707, and when the secondary flow velocity is as high as 0.2 times the maximum base flow velocity, and the ratio 2 Delta m/(m(1) + m(2)) is as high as 0.5. Here, the Reynolds number Re = rho(w)Omega R-2/mu, the stratification parameter A = root m Omega R-2(2)/(2k(B)T), R and Omega are the cylinder radius and angular velocity, m is the molecular mass, rho(w) is the wall density, mu is the viscosity and T is the temperature. The leading-order solutions do capture the qualitative trends, but are not in quantitative agreement.

Breakup processes of pressure swirl spray in gaseous cross-flow

Injection of liquid fuel in cross flowing air has been a strategy for future aircraft engines in order to control the emissions. In this context, breakup of a pressure swirl spray in gaseous cross-flow is investigated experimentally. The atomizer discharges a conical swirling sheet of liquid that interacts with cross-flowing air. This complex interaction and the resulting spray structures at various flow conditions are studied through flow visualization using still as well as high speed photography. Experiments are performed over a wide range of aerodynamic Weber number (2-300) and liquid-to-air momentum flux ratio (5-150). Various breakup regimes exhibiting different breakup processes are mapped on a parameter space based on flow conditions. This map shows significant variations from breakup regime map for a plain liquid jet in cross-flow. It is observed that the breakup of leeward side of the sheet is dominated by bag breakup and the windward side of the sheet undergoes breakup through surface waves. Similarities and differences between bag breakup present in plain liquid jet in cross-flow and swirl spray in cross-flow are explained. Multimodal drop size distribution from bag breakup, frequency of bag breakup, wavelength of surface waves and trajectory of spray in cross-flow are measured by analyzing the spray images and parametric study of their variations is also presented. (C) 2014 Elsevier Ltd. All rights reserved.

Tuning photoinduced terahertz conductivity in monolayer graphene: Optical-pump terahertz-probe spectroscopy

Optical-pump terahertz-probe differential transmission measurements of as-prepared single layer graphene (AG) (unintentionally hole dopedwith Fermi energy E-F at similar to -180 meV), nitrogen doping compensated graphene (NDG) with E-F similar to -10 meV, and thermally annealed doped graphene (TAG) are examined quantitatively to understand the opposite signs of photoinduced dynamic terahertz conductivity Delta sigma. It is negative for AG and TAG but positive for NDG. We show that the recently proposed mechanism of multiple generations of secondary hot carriers due to Coulomb interaction of photoexcited carriers with the existing carriers together with the intraband scattering can explain the change of photoinduced conductivity sign and its magnitude. We give a quantitative estimate of Delta sigma in terms of controlling parameters-the Fermi energy E-F and momentum relaxation time tau. Furthermore, the cooling of photoexcited carriers is analyzed using a supercollision model which involves a defect mediated collision of the hot carriers with the acoustic phonons, thus giving an estimate of the deformation potential.

Large deformation dynamic finite element analysis of delaminated composite plates using contact-impact conditions

A new C-0 composite plate finite element based on Reddy's third order theory is used for large deformation dynamic analysis of delaminated composite plates. The inter-laminar contact is modeled with an augmented Lagrangian approach. Numerical results show that the widely used ``unconditionally stable'' beta-Newmark method presents instability problems in the transient simulation of delaminated composite plate structures with large deformation. To overcome this instability issue, an energy and momentum conserving composite implicit time integration scheme presented by Bathe and Baig is used. It is found that a proper selection of the penalty parameter is very crucial in the contact simulation. (C) 2014 Elsevier Ltd. All rights reserved.

Time resolved terahertz spectroscopy of low frequency electronic resonances and optical pump-induced terahertz photoconductivity in reduced graphene oxide membrane

Towards ultrafast optoelectronic applications of single and a few layer reduced graphene oxide (RGO), we study time domain terahertz spectroscopy and optical pump induced changes in terahertz conductivity of self-supported RGO membrane in the spectral window of 0.5-3.5 THz. The real and imaginary parts of conductivity spectra clearly reveal low frequency resonances, attributed to the energy gaps due to the van Hove singularities in the density of states flanking the Dirac points arising due to the relative rotation of the graphene layers. Further, optical pump induced terahertz conductivity is positive, pointing to the dominance of intraband scattering processes. The relaxation dynamics of the photo-excited carriers consists of three cooling pathways: the faster (similar to 450 fs) one due to optical phonon emission followed by disorder mediated large momentum and large energy acoustic phonon emission with a time constant of a few ps (called the super-collision mechanism) and a very large time (similar to 100 ps) arising from the deep trap states. The frequency dependence of the dynamic conductivity at different delay times is analyzed in term of Drude-Smith model. (C) 2014 Published by Elsevier Ltd.

Penetration Characteristics of Air, Carbon Dioxide and Helium Transverse Sonic Jets in Mach 5 Cross Flow

An experimental investigation of sonic air, CO2 and Helium transverse jets in Mach 5 cross flow was carried out over a flat plate. The jet to freestream momentum flux ratio, J, was kept the same for all gases. The unsteady flow topology was examined using high speed schlieren visualisation and PIV. Schlieren visualisation provided information regarding oscillating jet shear layer structures and bow shock, Mach disc and barrel shocks. Two-component PIV measurements at the centreline, provided information regarding jet penetration trajectories. Barrel shocks and Mach disc forming the jet boundary were visualised/quantified also jet penetration boundaries were determined. Even though J is kept the same for all gases, the penetration patterns were found to be remarkably different both at the nearfield and the farfield. Air and CO2 jet resulted similar nearfield and farfield penetration pattern whereas Helium jet spread minimal in the nearfield.

Optical diagnostics of fuel-air mixing and vortex formation in a cavity combustor

The current work reports optical diagnostic measurements of fuel-air mixing and vortex structure in a single cavity trapped vortex combustor (TVC). Specifically, the mixture fraction using acetone PLIF technique in the non-reacting flow, and PIV measurements in the reacting flow are reported for the first time in trapped vortex combustors. The fuel-air momentum flux ratio, where the air momentum corresponds to that entering the cavity through a specially-incorporated flow guide vane, is used to characterize the mixing. The acetone PLIF experiments show that at high momentum flux ratios, the fuel-air mixing in the cavity is very minimal and is enhanced as the momentum flux ratio reduces, due to a favourable vortex formation in the cavity. Stoichiometric mixture fraction surfaces show that the mixing causes the reaction surfaces to shift from non-premixed to partially-premixed stratified mixtures. PIV measurements conducted in the non-reacting flow in the cavity further reinforce this observation. The scalar dissipation rates of mixture fraction were compared with the contours of RMS of fluctuating velocity and showed very good agreement. The regions of maximum mixing are observed to be along the fuel air interface. Reacting flow Ply measurements which differ substantially from the non-reacting cases primarily because of the heat release from combustion and the resulting gas expansion show that the vortex is displaced from the centre of the cavity towards the guide vane. Overall, the measurements show interesting features of the flow including the presence of the dual cavity structure and lead to a clear understanding of the underlying physics of the cavity flow highlighting the importance of the fuel-air momentum ratio parameter. (C) 2014 Elsevier Inc. All rights reserved.

Transverse single spin asymmetry in e plus p up arrow -> e plus J/psi plus X and Q(2) evolution of Sivers function-II

We present estimates of single spin asymmetry in the electroproduction of J/psi taking into account the transverse momentum-dependent (TMD) evolution of the gluon Sivers function. We estimate single spin asymmetry for JLab, HERMES, COMPASS and eRHIC energies using the color evaporation model of J/psi. We have calculated the asymmetry using recent parameters extracted by Echevarria et al. using the Collins-Soper-Sterman approach to TMD evolution. These recent TMD evolution fits are based on the evolution kernel in which the perturbative part is resummed up to next-to-leading logarithmic accuracy. We have also estimated the asymmetry by using parameters which had been obtained by a fit by Anselmino et al., using both an exact numerical and an approximate analytical solution of the TMD evolution equations. We find that the variation among the different estimates obtained using TMD evolution is much smaller than between these on one hand and the estimates obtained using DGLAP evolution on the other. Even though the use of TMD evolution causes an overall reduction in asymmetries compared to the ones obtained without it, they remain sizable. Overall, upon use of TMD evolution, predictions for asymmetries stabilize.

Transverse Single Spin Asymmetries and Charmonium Production

We estimate transverse spin single spin asymmetry(TSSA) in the process e + p(up arrow) -> J/psi + X using color evaporation model of charmonium production. We take into account transverse momentum dependent(TMD) evolution of Sivers function and parton distribution function and show that the there is a reduction in the asymmetry as compared to our earlier estimates wherein the Q(2) - evolution was implemented only through DGLAP evolution of unpolarized gluon densities.

Associated jet and subjet rates in light-quark and gluon jet discrimination

We show that in studies of light quark- and gluon-initiated jet discrimination, it is important to include the information on softer reconstructed jets (associated jets) around a primary hard jet. This is particularly relevant while adopting a small radius parameter for reconstructing hadronic jets. The probability of having an associated jet as a function of the primary jet transverse momentum (PT) and radius, the minimum associated jet pi, and the association radius is computed up to next-to-double logarithmic accuracy (NDLA), and the predictions are compared with results from Herwig++, Pythia6 and Pythia8 Monte Carlos (MC). We demonstrate the improvement in quark-gluon discrimination on using the associated jet rate variable with the help of a multivariate analysis. The associated jet rates are found to be only mildly sensitive to the choice of parton shower and hadronization algorithms, as well as to the effects of initial state radiation and underlying event. In addition, the number of k(t) subjets of an anti-k(t) jet is found to be an observable that leads to a rather uniform prediction across different MC's, broadly being in agreement with predictions in NDLA, as compared to the often used number of charged tracks observable.

Airblast spray in crossflow - Structure, trajectory and droplet sizing

This study reports results of an experimental investigation of airblast spray of water and ethanol in crossflow. Laser shadowgraphy and Particle/Droplet Imaging Analysis (PDIA) are used to derive spray trajectory and drop size information while Particle Tracking Velocimetry (PTV) is used to measure droplet velocities. A new phenomenon of spray bifurcation is observed for low Gas to Liquid Ratio (GLR) cases. The reasons for the spatial bifurcation can be attributed to a combination of reasons. These are (a) presence of large ligaments and droplets in the near-nozzle region for low GLRs (b) secondary breakup experienced by ligaments/droplets leading to formation of a large number of small droplets, and (c) the crossflow causing differential dispersion of the small and large droplets. A novel correlation for spray trajectory is proposed incorporating the momentum ratio and liquid surface tension. This correlation is shown to be effective in predicting the non-linear spray trajectory over a large range of conditions for not only water but ethanol and Jet-A also. It is observed that the larger droplets penetrate further into the crossflow, in the direction of injection. Thus, with increase in height of the measurement location from the injection plane, the droplet Sauter Mean Diameter (SMD) is found to increase. Moreover, as the droplets travel downstream in the crossflow direction, the droplet SMD is observed to decrease. The effect of drag is assessed by comparing velocity of different sizes of droplets at various locations. Smaller droplets are entrained into the crossflow at much lower elevations, whereas larger droplets tend to penetrate further into the crossflow. (C) 2015 Elsevier Ltd. All rights reserved.

Numerical study of solidification of A356 aluminum alloy flowing on an oblique plate with experimental validation

Molten A356 aluminum alloy flowing on an oblique plate is water cooled from underneath. The melt partially solidifies on plate wall with continuous formation of columnar dendrites. These dendrites are continuously sheared off into equiaxed/fragmented grains and carried away with the melt by producing semisolid slurry collected at plate exit. Melt pouring temperature provides required solidification whereas plate inclination enables necessary shear for producing slurry of desired solid fraction. A numerical model concerning transport equations of mass, momentum, energy and species is developed for predicting velocity, temperature, macrosegregation and solid fraction. The model uses FVM with phase change algorithm, VOF and variable viscosity. The model introduces solid phase movement with gravity effect as well. Effects of melt pouring temperature and plate inclination on hydrodynamic and thermo-solutal behaviors are studied subsequently. Slurry solid fractions at plate exit are 27%, 22%, 16%, and 10% for pouring temperatures of 620 degrees C, 625 degrees C, 630 degrees C, and 635 degrees C, respectively. And, are 27%, 25%, 22%, and 18% for plate inclinations of 30, 45, 60, and 75, respectively. Melt pouring temperature of 625 degrees C with plate inclination of 60 generates appropriate quality of slurry and is the optimum. Both numerical and experimental results are in good agreement with each other. (C) 2015 Taiwan Institute of Chemical Engineers. Published by Elsevier B.V. All rights reserved.