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.

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.

Temperature-Induced Reversible First-Order Single Crystal to Single Crystal Phase Transition in Boc-gamma(4)(R)Val-Val-OH: Interplay of Enthalpy and Entropy

Crystals of Boc-gamma y(4)(R)Val-Val-OH undergo a reversible first-order single crystal to single crystal phase transition at T-c approximate to 205 K from the orthorhombic space group P22(1)2(1) (Z' = 1) to the monoclinic space group P2(1) (Z' = 2) with a hysteresis of similar to 2.1 K. The low-temperature monoclinic form is best described as a nonmerohedral twin with similar to 50% contributions from its two components. The thermal behavior of the dipeptide crystals was characterized by differential scanning calorimetry experiments. Visual changes in birefringence of the sample during heating and cooling cycles on a hot-stage microscope with polarized light supported the phase transition. Variable-temperature unit cell check measurements from 300 to 100 K showed discontinuity in the volume and cell parameters near the transition temperature, supporting the first-order behavior. A detailed comparison of the room-temperature orthorhombic form with the low-temperature (100 K) monoclinic form revealed that the strong hydrogen-bonding motif is retained in both crystal systems, whereas the non-covalent interactions involving side chains of the dipeptide differ significantly, leading to a small change in molecular conformation in the monoclinic form as well as a small reorientation of the molecules along the ac plane. A rigid-body thermal motion analysis (translation, libration, screw; correlation of translation and libration) was performed to study the crystal entropy. The reversible nature of the phase transition is probably the result of an interplay between enthalpy and entropy: the low-temperature monoclinic form is enthalpically favored, whereas the room-temperature orthorhombic form is entropically favored.

Probing a spin-glass state in SrRuO3 thin films through higher-order statistics of resistance fluctuations

The complex perovskite oxide SrRuO3 shows intriguing transport properties at low temperatures due to the interplay of spin, charge, and orbital degrees of freedom. One of the open questions in this system is regarding the origin and nature of the low-temperature glassy state. In this paper we report on measurements of higher-order statistics of resistance fluctuations performed in epitaxial thin films of SrRuO3 to probe this issue. We observe large low-frequency non-Gaussian resistance fluctuations over a certain temperature range. Our observations are compatible with that of a spin-glass system with properties described by hierarchical dynamics rather than with that of a simple ferromagnet with a large coercivity.

Hydrogen bonding, halogen bonding and lithium bonding: an atoms in molecules and natural bond orbital perspective towards conservation of total bond order, inter- and intra-molecular bonding

One hundred complexes have been investigated exhibiting D-X center dot center dot center dot A interactions, where X = H, Cl or Li and DX is the `X bond' donor and A is the acceptor. The optimized structures of all these complexes have been used to propose a generalized `Legon-Millen rule' for the angular geometry in all these interactions. A detailed Atoms in Molecules (AIM) theoretical analysis confirms an important conclusion, known in the literature: there is a strong correlation between the electron density at the X center dot center dot center dot A bond critical point (BCP) and the interaction energy for all these interactions. In addition, we show that extrapolation of the fitted line leads to the ionic bond for Li-bonding (electrostatic) while for hydrogen and chlorine bonding, it leads to the covalent bond. Further, we observe a strong correlation between the change in electron density at the D-X BCP and that at the X center dot center dot center dot A BCP, suggesting conservation of the bond order. The correlation found between penetration and electron density at BCP can be very useful for crystal structure analysis, which relies on arbitrary van der Waals radii for estimating penetration. Various criteria proposed for shared-and closed-shell interactions based on electron density topology have been tested for H/Cl/Li bonded complexes. Finally, using the natural bond orbital (NBO) analysis it is shown that the D-X bond weakens upon X bond formation, whether it is ionic (DLi) or covalent (DH/DCl) and the respective indices such as ionicity or covalent bond order decrease. Clearly, one can think of conservation of bond order that includes ionic and covalent contributions to both D-X and X center dot center dot center dot A bonds, for not only X = H/Cl/Li investigated here but also any atom involved in intermolecular bonding.

Assignment of the C-13 NMR spectrum by correlation to dipolar coupled proton-pairs and estimation of order parameters of a thiophene based liquid crystal

Materials with widely varying molecular topologies and exhibiting liquid crystalline properties have attracted considerable attention in recent years. C-13 NMR spectroscopy is a convenient method for studying such novel systems. In this approach the assignment of the spectrum is the first step which is a non-trivial problem. Towards this end, we propose here a method that enables the carbon skeleton of the different sub-units of the molecule to be traced unambiguously. The proposed method uses a heteronuclear correlation experiment to detect pairs of nearby carbons with attached protons in the liquid crystalline core through correlation of the carbon chemical shifts to the double-quantum coherences of protons generated through the dipolar coupling between them. Supplemented by experiments that identify non-protonated carbons, the method leads to a complete assignment of the spectrum. We initially apply this method for assigning the C-13 spectrum of the liquid crystal 4-n-pentyl-4'-cyanobiphenyl oriented in the magnetic field. We then utilize the method to assign the aromatic carbon signals of a thiophene based liquid crystal thereby enabling the local order-parameters of the molecule to be estimated and the mutual orientation of the different sub-units to be obtained.

Formation Flying with Nonlinear Partial Integrated Guidance and Control

Using the dynamic inversion philosophy, a nonlinear partial integrated guidance and control approach is presented in this paper for formation flying. It is based on the evolving philosophy of integrated guidance and control. However, it also retains the advantages of the conventional guidance then control philosophy by retaining the timescale separation between translational and rotational dynamics explicitly. Simulation studies demonstrate that the proposed technique is effective in bringing the vehicles into formation quickly and maintaining the formation.

A 5th and 7th order harmonic suppression scheme using capacitive filtering for 2-level VSI fed IM drive

In this paper, a 5th and 7th harmonic suppression technique for a 2-level VSI fed IM drive, by using capacitive filtering is proposed. A capacitor fed 2-level inverter is used on an open-end winding induction motor to suppress all 5th and 7th order harmonics. A PWM scheme that maintains the capacitor voltage, while suppressing the harmonics is also proposed. The proposed scheme is valid for the entire modulation range, including overmodulation and six-step mode of operation of the main inverter.

Adaptive mesh generation for singularly perturbed fourth-order ordinary differential equations

In this paper, we consider a singularly perturbed boundary-value problem for fourth-order ordinary differential equation (ODE) whose highest-order derivative is multiplied by a small perturbation parameter. To solve this ODE, we transform the differential equation into a coupled system of two singularly perturbed ODEs. The classical central difference scheme is used to discretize the system of ODEs on a nonuniform mesh which is generated by equidistribution of a positive monitor function. We have shown that the proposed technique provides first-order accuracy independent of the perturbation parameter. Numerical experiments are provided to validate the theoretical results.

Wave propagation analysis in adhesively bonded composite joints using the wavelet spectral finite element method

A wavelet spectral finite element (WSFE) model is developed for studying transient dynamics and wave propagation in adhesively bonded composite joints. The adherands are formulated as shear deformable beams using the first order shear deformation theory (FSDT) to obtain accurate results for high frequency wave propagation. Equations of motion governing wave motion in the bonded beams are derived using Hamilton's principle. The adhesive layer is modeled as a line of continuously distributed tension/compression and shear springs. Daubechies compactly supported wavelet scaling functions are used to transform the governing partial differential equations from time domain to frequency domain. The dynamic stiffness matrix is derived under the spectral finite element framework relating the nodal forces and displacements in the transformed frequency domain. Time domain results for wave propagation in a lap joint are validated with conventional finite element simulations using Abaqus. Frequency domain spectrum and dispersion relation results are presented and discussed. The developed WSFE model yields efficient and accurate analysis of wave propagation in adhesively-bonded composite joints. (C) 2014 Elsevier Ltd. All rights reserved.

Transitory second-order reciprocal connection for two surfaces in point contact

This paper investigates the instantaneous spatial higher pair to lower pair substitute-connection which is kinematically equivalent up to acceleration analysis for two smooth surfaces in point contact. The existing first-order equivalent substitute-connection consisting of a Hooke's joint (U-joint) and a spherical joint (S-joint) connected by an additional link is extended up to second-order. A two step procedure is chalked out for achieving this equivalence. First, the existing method is employed for velocity equivalence. In the second step, the two centers of substitution are obtained as a conjugate relationship involving the principal normal curvatures of the surfaces at the contact point and the screw coordinates of the instantaneous screw axis (ISA) of the first-order relative motion. Unlike the classical planar replacement, this particular substitution cannot be done by merely examining the profiles of the contacting surfaces. An illustrative example of a three-link direct-contact mechanism is presented. (C) 2014 Elsevier Ltd. All rights reserved.

REDUCED ORDER MODELLING OF COMBUSTION INSTABILITY IN A BACKWARD FACING STEP COMBUSTOR

This paper develops a fully coupled time domain Reduced Order Modelling (ROM) approach to model unsteady combustion dynamics in a backward facing step combustor The acoustic field equations are projected onto the canonical acoustic eigenmodes of the systems to obtain a coupled system of modal evolution equations. The heat release response of the flame is modelled using the G-equation approach. Vortical velocity fluctuations that arise due to shear layer rollup downstream of the step are modelled using a simplified 1D-advection equation whose phase speed is determined from a linear, local, temporal stability analysis of the shear layer just downstream of the step. The hydrodynamic stability analysis reveals a abrupt change in the value of disturbance phase speed from unity for Re < Re-crit to 0.5 for Re > Re-crit, where Remit for the present geometry was found to be approximate to 10425. The results for self-excited flame response show highly wrinkled flame shapes that are qualitatively similar to those seen in prior experiments of acoustically forced flames. The effect of constructive and destructive interference between the two contributions to flame surface wrinkling results in high amplitude wrinkles for the case when K-c -> 1.

Orientational order as the origin of the long-range hydrophobic effect

The long range attractive force between two hydrophobic surfaces immersed in water is observed to decrease exponentially with their separation-this distance-dependence of effective force is known as the hydrophobic force law (HFL). We explore the microscopic origin of HFL by studying distance-dependent attraction between two parallel rods immersed in 2D Mercedes Benz model of water. This model is found to exhibit a well-defined HFL. Although the phenomenon is conventionally explained by density-dependent theories, we identify orientation, rather than density, as the relevant order parameter. The range of density variation is noticeably shorter than that of orientational heterogeneity. The latter is comparable to the observed distances of hydrophobic force. At large separation, attraction between the rods arises primarily from a destructive interference among the inwardly propagating oppositely oriented heterogeneity generated in water by the two rods. As the rods are brought closer, the interference increases leading to a decrease in heterogeneity and concomitant decrease in free energy of the system, giving rise to the effective attraction. We notice formation of hexagonal ice-like structures at the onset of attractive region which suggests that metastable free energy minimum may play a role in the origin of HFL. (C) 2015 AIP Publishing LLC.

Direct measurements of growing amorphous order and non-monotonic dynamic correlations in a colloidal glass-former

The transformation of flowing liquids into rigid glasses is thought to involve increasingly cooperative relaxation dynamics as the temperature approaches that of the glass transition. However, the precise nature of this motion is unclear, and a complete understanding of vitrification thus remains elusive. Of the numerous theoretical perspectives(1-4) devised to explain the process, random first-order theory (RFOT; refs 2,5) is a well-developed thermodynamic approach, which predicts a change in the shape of relaxing regions as the temperature is lowered. However, the existence of an underlying `ideal' glass transition predicted by RFOT remains debatable, largely because the key microscopic predictions concerning the growth of amorphous order and the nature of dynamic correlations lack experimental verification. Here, using holographic optical tweezers, we freeze a wall of particles in a two-dimensional colloidal glass-forming liquid and provide direct evidence for growing amorphous order in the form of a static point-to-set length. We uncover the non-monotonic dependence of dynamic correlations on area fraction and show that this non-monotonicity follows directly from the change in morphology and internal structure of cooperatively rearranging regions(6,7). Our findings support RFOT and thereby constitute a crucial step in distinguishing between competing theories of glass formation.

Wireless scheduling with partial channel state information: large deviations and optimality

We consider a server serving a time-slotted queued system of multiple packet-based flows, where not more than one flow can be serviced in a single time slot. The flows have exogenous packet arrivals and time-varying service rates. At each time, the server can observe instantaneous service rates for only a subset of flows ( selected from a fixed collection of observable subsets) before scheduling a flow in the subset for service. We are interested in queue length aware scheduling to keep the queues short. The limited availability of instantaneous service rate information requires the scheduler to make a careful choice of which subset of service rates to sample. We develop scheduling algorithms that use only partial service rate information from subsets of channels, and that minimize the likelihood of queue overflow in the system. Specifically, we present a new joint subset-sampling and scheduling algorithm called Max-Exp that uses only the current queue lengths to pick a subset of flows, and subsequently schedules a flow using the Exponential rule. When the collection of observable subsets is disjoint, we show that Max-Exp achieves the best exponential decay rate, among all scheduling algorithms that base their decision on the current ( or any finite past history of) system state, of the tail of the longest queue. To accomplish this, we employ novel analytical techniques for studying the performance of scheduling algorithms using partial state, which may be of independent interest. These include new sample-path large deviations results for processes obtained by non-random, predictable sampling of sequences of independent and identically distributed random variables. A consequence of these results is that scheduling with partial state information yields a rate function significantly different from scheduling with full channel information. In the special case when the observable subsets are singleton flows, i.e., when there is effectively no a priori channel state information, Max-Exp reduces to simply serving the flow with the longest queue; thus, our results show that to always serve the longest queue in the absence of any channel state information is large deviations optimal.