Particle dynamics in the channel flow of a turbulent particle-gas suspension at high Stokes number. Part 1. DNS and fluctuating force model

The fluctuating force model is developed and applied to the turbulent flow of a gas-particle suspension in a channel in the limit of high Stokes number, where the particle relaxation time is large compared to the fluid correlation time, and low particle Reynolds number where the Stokes drag law can be used to describe the interaction between the particles and fluid. In contrast to the Couette flow, the fluid velocity variances in the different directions in the channel are highly non-homogeneous, and they exhibit significant variation across the channel. First, we analyse the fluctuating particle velocity and acceleration distributions at different locations across the channel. The distributions are found to be non-Gaussian near the centre of the channel, and they exhibit significant skewness and flatness. However, acceleration distributions are closer to Gaussian at locations away from the channel centre, especially in regions where the variances of the fluid velocity fluctuations are at a maximum. The time correlations for the fluid velocity fluctuations and particle acceleration fluctuations are evaluated, and it is found that the time correlation of the particle acceleration fluctuations is close to the time correlations of the fluid velocity in a `moving Eulerian' reference, moving with the mean fluid velocity. The variances of the fluctuating force distributions in the Langevin simulations are determined from the time correlations of the fluid velocity fluctuations and the results are compared with direct numerical simulations. Quantitative agreement between the two simulations are obtained provided the particle viscous relaxation time is at least five times larger than the fluid integral time.

Particle dynamics in the channel flow of a turbulent particle-gas suspension at high Stokes number. Part 2. Comparison of fluctuating force simulations and experiments

The particle and fluid velocity fluctuations in a turbulent gas-particle suspension are studied experimentally using two-dimensional particle image velocimetry with the objective of comparing the experiments with the predictions of fluctuating force simulations. Since the fluctuating force simulations employ force distributions which do not incorporate the modification of fluid turbulence due to the particles, it is of importance to quantify the turbulence modification in the experiments. For experiments carried out at a low volume fraction of 9.15 x 10(-5) (mass loading is 0.19), where the viscous relaxation time is small compared with the time between collisions, it is found that the gas-phase turbulence is not significantly modified by the presence of particles. Owing to this, quantitative agreement is obtained between the results of experiments and fluctuating force simulations for the mean velocity and the root mean square of the fluctuating velocity, provided that the polydispersity in the particle size is incorporated in the simulations. This is because the polydispersity results in a variation in the terminal velocity of the particles which could induce collisions and generate fluctuations; this mechanism is absent if all of the particles are of equal size. It is found that there is some variation in the particle mean velocity very close to the wall depending on the wall-collision model used in the simulations, and agreement with experiments is obtained only when the tangential wall-particle coefficient of restitution is 0.7. The mean particle velocity is in quantitative agreement for locations more than 10 wall units from the wall of the channel. However, there are systematic differences between the simulations and theory for the particle concentrations, possibly due to inadequate control over the particle feeding at the entrance. The particle velocity distributions are compared both at the centre of the channel and near the wall, and the shape of the distribution function near the wall obtained in experiments is accurately predicted by the simulations. At the centre, there is some discrepancy between simulations and experiment for the distribution of the fluctuating velocity in the flow direction, where the simulations predict a bi-modal distribution whereas only a single maximum is observed in the experiments, although both distributions are skewed towards negative fluctuating velocities. At a much higher particle mass loading of 1.7, where the time between collisions is smaller than the viscous relaxation time, there is a significant increase in the turbulent velocity fluctuations by similar to 1-2 orders of magnitude. Therefore, it becomes necessary to incorporate the modified fluid-phase intensity in the fluctuating force simulation; with this modification, the mean and mean-square fluctuating velocities are within 20-30% of the experimental values.

Temperature variation of the grüneisen parameter in magnesium oxide

The thermal expansion of magnesium oxide has been measured below room temperature from 140°K to 284.5°K, using an interferometric method. The accuracy of measurement is better than 3% in the temperature range studied. The agreement of these results with Durand's is quite good, but consistently higher over most of the range by 2 or 3%, for the most part within the estimated experimental error. The Grüneisen parameter remains constant at about 1.51 over the present experimental range; but an isolated measurement of Durand at 85°K suggests that at lower temperatures it rises quite sharply above this value. This possibility is therefore investigated theoretically. With a non-central force model to represent MgO, γ(−3) and γ(2) are calculated and it is found that γ(−3) > γ(2), again suggesting that the Grüneisen parameter increases with falling temperature. Of the two reported experimental values for the infra-red absorption frequency, correlation with the heat capacity strongly indicates a wavelength of 25.26μm rather than 17.3μm. Thermal expansion measurements at still lower temperatures must be carried out to confirm definitely the rise in the Grüneisen parameter.

Force measurements on hypersonic waveriders in the IISc hypersonic shock tunnel HST2

The drag and lift coefficients for a viscous optimized Mach 6 conical waverider has been measured using an accelerometer force balance system in the IISc hypersonic shock tunnel. A rubber bush placed in between the waverider model and the steel sting ensures unrestrained motion to the model during shock tunnel testing (500 mu s). Two accelerometers mounted on the model are used to measure the model accelerations in the axial and normal directions. The measured value of lift to drag ratio at zero angle of incidence for the IISc conical waverider with viscous optimized leading edge is 2.149, which compares well with the value reported in the open literature (Anderson et al 1991) for similar class of waveriders designed for a flight Mach number of 6. The details of the experimental study along with illustrative numerical results are discussed in this paper.

A new 6-component accelerometer force balance for short duration ground testing facilities

A new six-component accelerometer force balance is developed and used in the HST2 shock tunnel of Indian Institute of Science. Aerodynamic forces and moments for a hypersonic slender body measured using this balance system at a free stream Mach number of 5.75 and Reynolds number of 1.5 million and stagnation enthalpy of 1.5 and 2 MJ/kg are presented. These measured values compare well with the theoretical values estimated using modified Newtonian theory.

A novel free floating accelerometer force balance system for shock tunnel applications

In order to overcome the interference of the model mounting system with the external aerodynamics of the body during shock tunnel testing, a new free floating internally mountable balance system that ensures unrestrained model motion during testing has been designed, fabricated and tested. Minimal friction ball bearings are used for ensuring the free floating condition of the model during tunnel testing. The drag force acting on a blunt leading edge flat plate at hypersonic Mach number has been measured using the new balance system. Finite element model (FEM) and CFD are exhaustively used in the design as well as for calibrating the new balance system. The experimentally measured drag force on the blunt leading edge flat plate at stagnation enthalpy of 0.7 and 1.2 MJ/kg and nominal Mach number of 5.75 matches well with FEM results. The concept can also be extended for measuring all the three fundamental aerodynamic forces in short duration test facilities like free piston driven shock tunnels.

Tensile force distribution along the reinforcement for reinforced soil foundations

In a reinforced soil bed system reinforcement layer is usually placed with or without end anchorage. Since soil is weak in tension reinforcement develop tension under the applied load or the displacement of the footing. This tensile force is distributed along the length of the reinforcement subjected to the end condition. The reinforccement tension helps in distributing the load over a wider area, and becomes more effective at large induced settlements. As a result, vertical componenent of tensile force generated becomes effective in reducing applied load. However, very few studies to quantify the tensile force along the reinforcement have been reported in the literature. In this paper an attempt has been made to obtain the true nature of tensile force distribution along the reinforcement. For a reinforced soil bed below a strip footing this paper brings out induced tensile force distribution along the reinforcement at different load levels and for different types of reinforcements.

Studies on Field Dependent Domain Structures in Multi-Grained 0.85PbMg(1/3)Nb(2/3)O3-0.15PbTiO(3) Thin Films by Scanning Force Microscopy

0.85PbMg(1/3)Nb(2/3)O(3)-0.15PbTiO(3) (0.85PMN-0.15PT) ferroelectric relaxor thin films have been deposited on La0.5Sr0.5CoO3/(111) Pt/TiO2/SiO2/Si by pulsed laser ablation by varying the oxygen partial pressures from 50 mTorr to 400 mTorr. The X-ray diffraction pattern reveals a pyrochlore free polycrystalline film. The grain morphology of the deposited films was studied using scanning electron microscopy and was found to be affected by oxygen pressure. By employing dynamic contact-electrostatic force microscopy we found that the distribution of polar nanoregions is majorly affected by oxygen pressure. Finally, the electric field induced switching in these films is discussed in terms of domain wall pinning.

Online Computation of Hysteresis Boundary for Constant Switching Frequency Current-Error Space-Vector-Based Hysteresis Controller for VSI-Fed IM Drives

This paper proposes a current-error space-vector-based hysteresis controller with online computation of boundary for two-level inverter-fed induction motor (IM) drives. The proposed hysteresis controller has got all advantages of conventional current-error space-vector-based hysteresis controllers like quick transient response, simplicity, adjacent voltage vector switching, etc. Major advantage of the proposed controller-based voltage-source-inverters-fed drive is that phase voltage frequency spectrum produced is exactly similar to that of a constant switching frequency space-vector pulsewidth modulated (SVPWM) inverter. In this proposed hysteresis controller, stator voltages along alpha- and beta-axes are estimated during zero and active voltage vector periods using current errors along alpha- and beta-axes and steady-state model of IM. Online computation of hysteresis boundary is carried out using estimated stator voltages in the proposed hysteresis controller. The proposed scheme is simple and capable of taking inverter upto six-step-mode operation, if demanded by drive system. The proposed hysteresis-controller-based inverter-fed drive scheme is experimentally verified. The steady state and transient performance of the proposed scheme is extensively tested. The experimental results are giving constant frequency spectrum for phase voltage similar to that of constant frequency SVPWM inverter-fed drive.

A MATLAB Code for Three Dimensional Linear Elastostatics using Constant Boundary Elements

Present work presents a code written in the very simple programming language MATLAB, for three dimensional linear elastostatics, using constant boundary elements. The code, in full or in part, is not a translation or a copy of any of the existing codes. Present paper explains how the code is written, and lists all the formulae used. Code is verified by using the code to solve a simple problem which has the well known approximate analytical solution. Of course, present work does not make any contribution to research on boundary elements, in terms of theory. But the work is justified by the fact that, to the best of author’s knowledge, as of now, one cannot find an open access MATLAB code for three dimensional linear elastostatics using constant boundary elements. Author hopes this paper to be of help to beginners who wish to understand how a simple but complete boundary element code works, so that they can build upon and modify the present open access code to solve complex engineering problems quickly and easily. The code is available online for open access (as supplementary file for the present paper), and may be downloaded from the website for the present journal.

Morphology dependent electrochemical performance of sputter deposited Sn thin films

This study deals with tailoring of the surface morphology, microstructure, and electrochemical properties of Sn thin films deposited by magnetron sputtering with different deposition rates. Scanning electron microscopy and atomic force microscopy are used to characterize the film surface morphology. Electrochemical properties of Sn thin film are measured and compared by cyclic voltammetry and charge-discharge cycle data at a constant current density. Sn thin film fabricated with a higher deposition rate exhibited an initial discharge capacity of 798 mAh g(-1) but reduced to 94 mAh g(-1) at 30th cycle. Film deposited with lower deposition rate delivered 770 mAh g(-1) during 1st cycle with improved capacity retention of 521 mAh g(-1) on 30th cycle. Comparison of electrochemical performances of these films has revealed important distinctions, which are associated with the surface morphology and hence on rate of deposition. (C) 2012 Elsevier Ltd. All rights reserved.

Study of n-ZnO/p-Si (100) thin film heterojunctions by pulsed laser deposition without buffer layer

ZnO/Si heterojunctions were fabricated by growing ZnO thin films on p-type Si (100) substrate by pulsed laser deposition without buffer layers. The crystallinity of the heterojunction was analyzed by high resolution X-ray diffraction and atomic force microscopy. The optical quality of the film was analyzed by room temperature (RT) photoluminescence measurements. The high intense band to band emission confirmed the high quality of the ZnO thin films on Si. The electrical properties of the junction were studied by temperature dependent current-voltage measurements and RT capacitance-voltage (C-V) analysis. The charge carrier concentration and the barrier height (BH) were calculated, to be 5.6x10(19) cm(-3) and 0.6 eV respectively from the C-V plot. The BH and ideality factor, calculated using the thermionic emission (TE) model, were found to be highly temperature dependent. We observed a much lower value in Richardson constant, 5.19x10(-7)A/cm(2) K-2 than the theoretical value (32 A/cm(2) K-2) for ZnO. This analysis revealed the existence of a Gaussian distribution (GD) with a standard deviation of sigma(2)=0.035 V. By implementing the GD to the TE, the values of BH and Richardson constant were obtained as 1.3 eV and 39.97 A/cm(2) K-2 respectively from the modified Richardson plot. The obtained Richardson constant value is close to the theoretical value for n-ZnO. These high quality heterojunctions can be used for solar cell applications. (C) 2012 Elsevier B.V. All rights reserved.

Conformational dynamics of HIV-1 protease: a comparative molecular dynamics simulation study with multiple amber force fields

Flap dynamics of HIV-1 protease (HIV-pr) controls the entry of inhibitors and substrates to the active site. Dynamical models from previous simulations are not all consistent with each other and not all are supported by the NMR results. In the present work, the er effect of force field on the dynamics of HIV-pr is investigated by MD simulations using three AMBER force fields ff99, ff99SB, and ff03. The generalized order parameters for amide backbone are calculated from the three force fields and compared with the NMR S2 values. We found that the ff99SB and ff03 force field calculated order parameters agree reasonably well with the NMR S2 values, whereas ff99 calculated values deviate most from the NMR order parameters. Stereochemical geometry of protein models from each force field also agrees well with the remarks from NMR S2 values. However, between ff99SB and ff03, there are several differences, most notably in the loop regions. It is found that these loops are, in general, more flexible in the ff03 force field. This results in a larger active site cavity in the simulation with the ff03 force field. The effect of this difference in computer-aided drug design against flexible receptors is discussed.

A Compliant End-Effector to Passively Limit the Force in Tele-Operated Tissue-Cutting

This paper presents a compliant end-effector that cuts soft tissues and senses the cutting forces. The end-effector is designed to have an upper threshold on cutting forces to facilitate safe handling of tissue during automated cutting. This is demonstrated with nonlinear finite element analysis and experimental results obtained by cutting inhomogeneous phantom tissue. The cutting forces are estimated using a vision-based technique that uses amplified elastic deformation of the compliant end-effector. We also demonstrate an immersive tele-operated tissue-cutting system together with a haptic device that gives real-time force feedback to the user. DOI: 10.1115/1.4007638]

Minimum power consumption for drag reduction on a circular cylinder by tangential surface motion

We investigate the effect of a prescribed tangential velocity on the drag force on a circular cylinder in a spanwise uniform cross flow. Using a combination of theoretical and numerical techniques we make an attempt at determining the optimal tangential velocity profiles which will reduce the drag force acting on the cylindrical body while minimizing the net power consumption characterized through a non-dimensional power loss coefficient (C-PL). A striking conclusion of our analysis is that the tangential velocity associated with the potential flow, which completely suppresses the drag force, is not optimal for both small and large, but finite Reynolds number. When inertial effects are negligible (R e << 1), theoretical analysis based on two-dimensional Oseen equations gives us the optimal tangential velocity profile which leads to energetically efficient drag reduction. Furthermore, in the limit of zero Reynolds number (Re -> 0), minimum power loss is achieved for a tangential velocity profile corresponding to a shear-free perfect slip boundary. At finite Re, results from numerical simulations indicate that perfect slip is not optimum and a further reduction in drag can be achieved for reduced power consumption. A gradual increase in the strength of a tangential velocity which involves only the first reflectionally symmetric mode leads to a monotonic reduction in drag and eventual thrust production. Simulations reveal the existence of an optimal strength for which the power consumption attains a minima. At a Reynolds number of 100, minimum value of the power loss coefficient (C-PL = 0.37) is obtained when the maximum in tangential surface velocity is about one and a half times the free stream uniform velocity corresponding to a percentage drag reduction of approximately 77 %; C-PL = 0.42 and 0.50 for perfect slip and potential flow cases, respectively. Our results suggest that potential flow tangential velocity enables energetically efficient propulsion at all Reynolds numbers but optimal drag reduction only for Re -> infinity. The two-dimensional strategy of reducing drag while minimizing net power consumption is shown to be effective in three dimensions via numerical simulation of flow past an infinite circular cylinder at a Reynolds number of 300. Finally a strategy of reducing drag, suitable for practical implementation and amenable to experimental testing, through piecewise constant tangential velocities distributed along the cylinder periphery is proposed and analysed.