Experimentally validated Raman Monte Carlo simulation for a cuboid object to obtain Raman spectroscopic signatures for hidden material

In conventional Raman spectroscopic measurements of liquids or surfaces the preferred geometry for detection of the Raman signal is the backscattering (or reflection) mode. For non-transparent layered materials, sub-surface Raman signals have been retrieved using spatially offset Raman spectroscopy (SORS), usually with light collection in the same plane as the point of excitation. However, as a result of multiple scattering in a turbid medium, Raman photons will be emitted in all directions. In this study, Monte Carlo simulations for a three-dimensional layered sample with finite geometry have been performed to confirm the detectability of Raman signals at all angles and at all sides of the object. We considered a non-transparent cuboid container (high density polyethylene) with explosive material (ammonium nitrate) inside. The simulation results were validated with experimental Raman intensities. Monte Carlo simulation results reveal that the ratio of sub-surface to surface signals improves at geometries other than backscattering. In addition, we demonstrate through simulations the effects of the absorption and scattering coefficients of the layers, and that of the diameter of the excitation beam. The advantage of collecting light from all possible 4 angles, over other collection modes, is that this technique is not geometry specific and molecular identification of layers underneath non-transparent surfaces can be obtained with minimal interference from the surface layer. To what extent all sides of the object will contribute to the total signal will depend on the absorption and scattering coefficients and the physical dimensions. Copyright (c) 2015 John Wiley & Sons, Ltd.

Influence of fuel hydrogen fraction on syngas fueled SI engine: Fuel thermo-physical property analysis and in-cylinder experimental investigations

Hydrogen, either in pure form or as a gaseous fuel mixture specie enhances the fuel conversion efficiency and reduce emissions in an internal combustion engine. This is due to the reduction in combustion duration attributed to higher laminar flame speeds. Hydrogen is also expected to increase the engine convective heat flux, attributed (directly or indirectly) to parameters like higher adiabatic flame temperature, laminar flame speed, thermal conductivity and diffusivity and lower flame quenching distance. These factors (adversely) affect the thermo-kinematic response and offset some of the benefits. The current work addresses the influence of mixture hydrogen fraction in syngas on the engine energy balance and the thermo-kinematic response for close to stoichiometric operating conditions. Four different bio-derived syngas compositions with fuel calorific value varying from 3.14 MJ/kg to 7.55 MJ/kg and air fuel mixture hydrogen fraction varying from 7.1% to 14.2% by volume are used. The analysis comprises of (a) use of chemical kinetics simulation package CHEMKIN for quantifying the thermo-physical properties (b) 0-D model for engine in-cylinder analysis and (c) in-cylinder investigations on a two-cylinder engine in open loop cooling mode for quantifying the thermo-kinematic response and engine energy balance. With lower adiabatic flame temperature for Syngas, the in-cylinder heat transfer analysis suggests that temperature has little effect in terms of increasing the heat flux. For typical engine like conditions (700 K and 25 bar at CR of 10), the laminar flame speed for syngas exceeds that of methane (55.5 cm/s) beyond mixture hydrogen fraction of 11% and is attributed to the increase in H based radicals. This leads to a reduction in the effective Lewis number and laminar flame thickness, potentially inducing flame instability and cellularity. Use of a thermodynamic model to assess the isolated influence of thermal conductivity and diffusivity on heat flux suggests an increase in the peak heat flux between 2% and 15% for the lowest (0.420 MW/m(2)) and highest (0.480 MW/m(2)) hydrogen containing syngas over methane (0.415 MW/m(2)) fueled operation. Experimental investigations indicate the engine cooling load for syngas fueled engine is higher by about 7% and 12% as compared to methane fueled operation; the losses are seen to increase with increasing mixture hydrogen fraction. Increase in the gas to electricity efficiency is observed from 18% to 24% as the mixture hydrogen fraction increases from 7.1% to 9.5%. Further increase in mixture hydrogen fraction to 14.2% results in the reduction of efficiency to 23%; argued due to the changes in the initial and terminal stages of combustion. On doubling of mixture hydrogen fraction, the flame kernel development and fast burn phase duration decrease by about 7% and 10% respectively and the terminal combustion duration, corresponding to 90%-98% mass burn, increases by about 23%. This increase in combustion duration arises from the cooling of the near wall mixture in the boundary layer attributed to the presence of hydrogen. The enhancement in engine cooling load and subsequent reduction in the brake thermal efficiency with increasing hydrogen fraction is evident from the engine energy balance along with the cumulative heat release profiles. Copyright (C) 2015, Hydrogen Energy Publications, LLC. Published by Elsevier Ltd. All rights reserved.

A mobile robot with a two-degree-of-freedom suspension for traversing uneven terrain with minimal slip: Modeling, simulation and experiments

It is known in literature that a wheeled mobile robot (WMR) with fixed length axle will slip on an uneven terrain. One way to avoid wheel slip is to use a torus-shaped wheel with lateral tilt capability which allows the distance between the wheel-ground contact points to change even with a fixed length axle. Such an arrangement needs a two degree-of-freedom (DOF) suspension for the vertical and lateral tilting motion of the wheel. In this paper modeling, simulation, design and experimentation with a three-wheeled mobile robot, with torus-shaped wheels and a novel two DOF suspension allowing independent lateral tilt and vertical motion, is presented. The suspension is based on a four-bar mechanism and is called the double four-bar (D4Bar) suspension. Numerical simulations show that the three-wheeled mobile robot can traverse uneven terrain with low wheel slip. Experiments with a prototype three-wheeled mobile robot moving on a constructed uneven terrain along a straight line, a circular arc and a path representing a lane change, also illustrate the low slip capability of the three-wheeled mobile robot with the D4Bar suspension. (C) 2015 Elsevier Ltd. All rights reserved.

Kinematic flow patterns in slow deformation of a dense granular material

The kinematic flow pattern in slow deformation of a model dense granular medium is studied at high resolution using in situ imaging, coupled with particle tracking. The deformation configuration is indentation by a flat punch under macroscopic plane-strain conditions. Using a general analysis method, velocity gradients and deformation fields are obtained from the disordered grain arrangement, enabling flow characteristics to be quantified. The key observations are the formation of a stagnation zone, as in dilute granular flow past obstacles; occurrence of vortices in the flow immediately underneath the punch; and formation of distinct shear bands adjoining the stagnation zone. The transient and steady state stagnation zone geometry, as well as the strength of the vortices and strain rates in the shear bands, are obtained from the experimental data. All of these results are well-reproduced in exact-scale non-smooth contact dynamics simulations. Full 3D numerical particle positions from the simulations allow extraction of flow features that are extremely difficult to obtain from experiments. Three examples of these, namely material free surface evolution, deformation of a grain column below the punch and resolution of velocities inside the primary shear band, are highlighted. The variety of flow features observed in this model problem also illustrates the difficulty involved in formulating a complete micromechanical analytical description of the deformation.

Lower Bound Axisymmetric Limit Analysis Using Drucker-Prager Yield Cone and Simulation with Mohr-Coulomb Pyramid

This paper presents a lower bound limit analysis approach for solving an axisymmetric stability problem by using the Drucker-Prager (D-P) yield cone in conjunction with finite elements and nonlinear optimization. In principal stress space, the tip of the yield cone has been smoothened by applying the hyperbolic approximation. The nonlinear optimization has been performed by employing an interior point method based on the logarithmic barrier function. A new proposal has also been given to simulate the D-P yield cone with the Mohr-Coulomb hexagonal yield pyramid. For the sake of illustration, bearing capacity factors N-c, N-q and N-gamma have been computed, as a function of phi, both for smooth and rough circular foundations. The results obtained from the analysis compare quite well with the solutions reported from literature.

Looking for BSM physics using top-quark polarization and decay-lepton kinematic asymmetries

We explore beyond-standard-model (BSM) physics signatures in the l + jets channel of the t (t) over bar pair production process at the Tevatron and the LHC. We study the effects of BSM physics scenarios on the top-quark polarization and on the kinematics of the decay leptons. To this end, we construct asymmetries using the lepton energy and angular distributions. Further, we find their correlations with the top polarization, net charge asymmetry and top forward-backward asymmetry. We show that when used together, these observables can help discriminate effectively between SM and different BSM scenarios, which can lead to varying degrees of top polarization at the Tevatron as well as the LHC. We use two types of colored mediator models to demonstrate the effectiveness of proposed observables, an s-channel axigluon and a u-channel diquark.

Simulation of Effective Thermal Conductivity of Metal Hydride Packed Beds

Several mathematical models are available for estimation of effective thermal conductivity of nonreactive packed beds. Keeping in view the salient differences between metal hydride beds in which chemisorption of hydrogen takes place and conventional nonreactive packed beds, modified models are proposed here to predict the effective thermal conductivity. Variation in properties such as solid thermal conductivity and porosity during hydrogen absorption and desorption processes are incorporated. These extended models have been applied to simulate the effective thermal conductivity of the MmNi(4.5)Al(0.5) hydride bed and are compared with the experimental results. Applicability of the extended models for estimation of the effective thermal conductivity at different operating conditions such as pressure, temperature, and hydrogen concentration is discussed.

Computationally efficient models for simulation of non-ideal DC-DC converters operating in continuous and discontinuous conduction modes

This paper discusses dynamic modeling of non-isolated DC-DC converters (buck, boost and buck-boost) under continuous and discontinuous modes of operation. Three types of models are presented for each converter, namely, switching model, average model and harmonic model. These models include significant non-idealities of the converters. The switching model gives the instantaneous currents and voltages of the converter. The average model provides the ripple-free currents and voltages, averaged over a switching cycle. The harmonic model gives the peak to peak values of ripple in currents and voltages. The validity of all these models is established by comparing the simulation results with the experimental results from laboratory prototypes, at different steady state and transient conditions. Simulation based on a combination of average and harmonic models is shown to provide all relevant information as obtained from the switching model, while consuming less computation time than the latter.

Differential dynamics of the serotonin(1A) receptor in membrane bilayers of varying cholesterol content revealed by all atom molecular dynamics simulation

The serotonin(1A) receptor belongs to the superfamily of G protein-coupled receptors (GPCRs) and is a potential drug target in neuropsychiatric disorders. The receptor has been shown to require membrane cholesterol for its organization, dynamics and function. Although recent work suggests a close interaction of cholesterol with the receptor, the structural integrity of the serotonin(1A) receptor in the presence of cholesterol has not been explored. In this work, we have carried out all atom molecular dynamics simulations, totaling to 3s, to analyze the effect of cholesterol on the structure and dynamics of the serotonin(1A) receptor. Our results show that the presence of physiologically relevant concentration of membrane cholesterol alters conformational dynamics of the serotonin(1A) receptor and, on an average lowers conformational fluctuations. Our results show that, in general, transmembrane helix VII is most affected by the absence of membrane cholesterol. These results are in overall agreement with experimental data showing enhancement of GPCR stability in the presence of membrane cholesterol. Our results constitute a molecular level understanding of GPCR-cholesterol interaction, and represent an important step in our overall understanding of GPCR function in health and disease.

New insight into the architecture of oxy-anion pocket in unliganded conformation of GAT domains: A MD-simulation study

Human Guanine Monophosphate Synthetase (hGMPS) converts XMP to GMP, and acts as a bifunctional enzyme with N-terminal ``glutaminase'' (GAT) and C-terminal ``synthetase'' domain. The enzyme is identified as a potential target for anticancer and immunosuppressive therapies. GAT domain of enzyme plays central role in metabolism, and contains conserved catalytic residues Cys104, His190, and Glu192. MD simulation studies on GAT domain suggest that position of oxyanion in unliganded conformation is occupied by one conserved water molecule (W1), which also stabilizes that pocket. This position is occupied by a negatively charged atom of the substrate or ligand in ligand bound crystal structures. In fact, MD simulation study of Ser75 to Val indicates that W1 conserved water molecule is stabilized by Ser75, while Thr152, and His190 also act as anchor residues to maintain appropriate architecture of oxyanion pocket through water mediated H-bond interactions. Possibly, four conserved water molecules stabilize oxyanion hole in unliganded state, but they vacate these positions when the enzyme (hGMPS)-substrate complex is formed. Thus this study not only reveals functionally important role of conserved water molecules in GAT domain, but also highlights essential role of other non-catalytic residues such as Ser75 and Thr152 in this enzymatic domain. The results from this computational study could be of interest to experimental community and provide a testable hypothesis for experimental validation. Conserved sites of water molecules near and at oxyanion hole highlight structural importance of water molecules and suggest a rethink of the conventional definition of chemical geometry of inhibitor binding site.

Atomistic Simulation of Stacked Nucleosome Core Particles: Tail Bridging, the H4 Tail, and Effect of Hydrophobic Forces

We report the first atomistic simulation of two stacked nucleosome core particles (NCPs), with an aim to understand, in molecular detail, how they interact, the effect of salt concentration, and how different histone tails contribute to their interaction, with a special emphasis on the H4 tail, known to have the largest stabilizing effect on the NCP-NCP interaction. We do not observe specific K16-mediated interaction between the H4 tail and the H2A-H2B acidic patch, in contrast with the findings from crystallographic studies, but find that the stacking was stable even in the absence of this interaction. We perform simulations with the H4 tail (partially/completely) removed and find that the region between LYS-16 and LYS-20 of the H4 tail holds special importance in mediating the inter-NCP interaction. Performing similar tail-clipped simulations with the H3 tail removed, we compare the roles of the H3 and H4 tails in maintaining the stacking. We discuss the relevance of our simulation results to the bilayer and other liquid-crystalline phases exhibited by NCPs in vitro and, through an analysis of the histone-histone interface, identify the interactions that could possibly stabilize the inter-NCP interaction in these columnar mesophases. Through the mechanical disruption of the stacked nucleosome system using steered molecular dynamics, we quantify the strength of inter-NCP stacking in the presence and absence of salt. We disrupt the stacking at some specific sites of internucleosomal tail-DNA contact and perform a comparative quantification of the binding strengths of various tails in stabilizing the stacking. We also examine how hydrophobic interactions may contribute to the overall stability of the stacking and find a marked difference in the role of hydrophobic forces as compared with electrostatic forces in determining the stability of the stacked nucleosome system.

Efficient simulation and rendering of realistic motion of one-dimensional flexible objects

In gross motion of flexible one-dimensional (1D) objects such as cables, ropes, chains, ribbons and hair, the assumption of constant length is realistic and reasonable. The motion of the object also appears more natural if the motion or disturbance given at one end attenuates along the length of the object. In an earlier work, variational calculus was used to derive natural and length-preserving transformation of planar and spatial curves and implemented for flexible 1D objects discretized with a large number of straight segments. This paper proposes a novel idea to reduce computational effort and enable real-time and realistic simulation of the motion of flexible 1D objects. The key idea is to represent the flexible 1D object as a spline and move the underlying control polygon with much smaller number of segments. To preserve the length of the curve to within a prescribed tolerance as the control polygon is moved, the control polygon is adaptively modified by subdivision and merging. New theoretical results relating the length of the curve and the angle between the adjacent segments of the control polygon are derived for quadratic and cubic splines. Depending on the prescribed tolerance on length error, the theoretical results are used to obtain threshold angles for subdivision and merging. Simulation results for arbitrarily chosen planar and spatial curves whose one end is subjected to generic input motions are provided to illustrate the approach. (C) 2016 Elsevier Ltd. All rights reserved.