Modeling of Bubble-Size Distribution in Water and Freon Co-Blown Free Rise Polyurethane Foams

A model has been developed to simulate the foam characteristics obtained, when chemical (water) and physical (Freon) blowing agents are used together for the formation of polyurethane foams. The model considers the rate of reaction, the consequent rise in temperature of the reaction mixture, nucleation of bubbles, and mass transfer of CO2 and Freon to them till the time of gelation. The model is able to explain the experimental results available in literature. It further predicts that the nucleation period gets reduced with increase in water (at constant Freon content), whereas with increase in Freon (at constant water) concentration nucleation period decreases marginally leading to narrower bubble-size distribution. By the use of uniform sized nuclei added initially, the model predicts that the bubble-size distribution can be made independent of the rate of homogeneous nucleation and can, thus, offer an extra parameter for its control. (C) 2014 Wiley Periodicals, Inc.

Confinement enhances dispersion in nanoparticle-polymer blend films

Polymer nanocomposites constitute an important class of materials whose properties depend on the state of dispersion of the nanoparticles in the polymer matrix. Here we report the first observations of confinement-induced enhancement of dispersion in nanoparticle-polymer blend films. Systematic variation in the dispersion of nanoparticles with confinement for various compositions and matrix polymer chain dimensions has been observed. For fixed composition, strong reduction in glass transition temperature, T-g, is observed with decreasing blend-film thickness. The enhanced dispersion occurs without altering the polymer-particle interactions and seems to be driven by enhanced matrix-chain orientation propensity and a tendency to minimize the density gradients within the matrix. This implies the existence of two different mechanisms in polymer nanocomposites, which determines their state of dispersion and glass transition.

Emergent properties of the interferon-signalling network may underlie the success of hepatitis C treatment

Current interferon alpha-based treatment of hepatitis C virus (HCV) infection fails to cure a sizeable fraction of patients treated. The cause of this treatment failure remains unknown. Here using mathematical modelling, we predict treatment failure to be a consequence of the emergent properties of the interferon-signalling network. HCV induces bistability in the network, creating a new steady state where it can persist. Cells that admit the new steady state alone are refractory to interferon. Using a model of viral kinetics, we show that when the fraction of cells refractory to interferon in a patient exceeds a critical value, treatment fails. Direct-acting antivirals that suppress HCV replication can eliminate the new steady state, restoring interferon sensitivity and improving treatment response. Our study thus presents a new conceptual basis of HCV persistence and treatment response, elucidates the origin of the synergy between interferon and direct-acting antivirals, and facilitates rational treatment optimization.

Multimodal shape oscillations of droplets excited by an air stream

The shape dynamics of droplets exposed to an air jet at intermediate droplet Reynolds numbers is investigated. High speed imaging and hot-wire anemometry are employed to examine the mechanism of droplet oscillation. The theory that the vortex shedding behind the droplet induces oscillation is examined. In these experiments, no particular dominant frequency is found in the wake region of the droplet. Hence the inherent free-stream disturbances prove to be driving the droplet oscillations. The modes of droplet oscillation show a band of dominant frequencies near the corresponding natural frequency, further proving that there is no particular forcing frequency involved. In the frequency spectrum of the lowest mode of oscillation for glycerol at the highest Reynolds number, no response is observed below the threshold frequency corresponding to the viscous dissipation time scale. This selective suppression of lower frequencies in the case of glycerol is corroborated by scaling arguments. The influence of surface tension on the droplet oscillation is studied using ethanol as a test fluid. Since a lower surface tension reduces the natural frequency, ethanol shows lower excited frequencies. The oscillation levels of different fluids are quantified using the droplet aspect ratio and correlated in terms of Weber number and Ohnesorge number. (C) 2014 Elsevier Ltd. All rights reserved.

Metal-metal charge transfer and interfacial charge transfer mechanism for the visible light photocatalytic activity of cerium and nitrogen co-doped TiO2

Nanosized cerium and nitrogen co-doped TiO2 (Ce-TiO2-xNx) was synthesized by sol gel method and characterized by powder X-ray diffraction (PXRD), X-ray photoelectron spectroscopy (XPS), FESEM, Fourier transform infrared, N-2 adsorption and desorption methods, photoluminescence and ultraviolet-visible (UV-vis) DRS techniques. PXRD analysis shows the dopant decreases the crystallite sizes and slows the crystallization of the titania matrix. XPS confirm the existence of cerium ion in +3 or +4 state, and nitrogen in -3 state in Ce-TiO2-xNx. The modified surface of TiO2 provides highly active sites for the dyes at the periphery of the Ce-O-Ti interface and also inhibits Ce particles from sintering. UV-visible DRS studies show that the metal-metal charge transfer (MMCT) of Ti/Ce assembly (Ti4+/Ce3+ -> Ti3+/Ce4+) is responsible for the visible light photocatalytic activity. Photoluminescence was used to determine the effect of cerium ion on the electron-hole pair separation between the two interfaces Ce-TiO2-xNx and Ce2O3. This separation increases with the increase of cerium and nitrogen ion concentrations of doped samples. The degradation kinetics of methylene blue and methyl violet dyes in the presence of sol gel TiO2, Ce-TiO2-xNx and commercial Degussa P25 was determined. The higher visible light activity of Ce-TiO2-xNx was due to the participation of MMCT and interfacial charge transfer mechanism.

Effect of adsorbate loading on selectivity during adsorption of C/C and C/C n-alkane binary mixtures in silicalite

Adsorption experiments of mixtures of long chain alkanes into silicalite under liquid phase conditions show selectivity inversion and azeotrope formation. These effects are due to the subtle interplay between the size of the adsorbed molecules and pore topology of the adsorbent. In this study, the selective uptake of lighter component during liquid phase adsorption of C/C and C/C n-alkane binary mixtures in the zeolite silicalite is understood through configurational bias grand-canonical Monte Carlo molecular simulation technique and a coarse-grained siting analysis. The simulations are conducted under conditions of low and intermediate levels of loading. The siting pattern of the adsorbates inside the zeolite pores explain the selectivity as seen in experiments.

Evaluating methane storage targets: from powder samples to onboard storage systems

The development of a viable adsorbed natural gas onboard fuel system involves synthesizing materials that meet specific storage target requirements. We assess the impact on natural gas storage due to intermediate processes involved in taking a laboratory powder sample to an onboard packed or adsorbent bed module. We illustrate that reporting the V/V (volume of gas/volume of container) capacities based on powder adsorption data without accounting for losses due to pelletization and bed porosity, grossly overestimates the working storage capacity for a given material. Using data typically found for adsorbent materials that are carbon and MOF based materials, we show that in order to meet the Department of Energy targets of 180 V/V (equivalent STP) loading at 3.5 MPa and 298 K at the onboard packed bed level, the volumetric capacity of the pelletized sample should be at least 245 V/V and the corresponding gravimetric loading varies from 0.175 to 0.38 kg/kg for pellet densities ranging from 461.5 to 1,000 . With recent revision of the DOE target to 263 V/V at the onboard packed bed level, the volumetric loadings for the pelletized sample should be about 373 V/V.

In situ formation of silver nanowire networks on paper

Simple, universally adaptable techniques for fabricating conductive patterns are required to translate laboratory-scale innovations into low-cost solutions for the developing world. Silver nanostructures have emerged as attractive candidates for forming such conductive patterns. We report here the in situ formation of conductive silver-nanowire networks on paper, thereby eliminating the need for either cost-intensive ink formulation or substrate preparation or complex post-deposition sintering steps. Reminiscent of the photographic process of `salt printing', a desktop office printer was used to deposit desired patterns of silver bromide on paper, which were subsequently exposed to light and then immersed in a photographic developer. Percolating silver nanowire networks that conformally coated the paper fibres were formed after 10 min of exposure to light from a commercial halogen lamp. Thus, conductive and patterned films with sheet resistances of the order of 4 Omega/rectangle can be easily formed by combining two widely used processes - inkjet printing and photographic development.

Calculation of three-phase methane-ethane binary clathrate hydrate phase equilibrium from Monte Carlo molecular simulations

Methane and ethane are the simplest hydrocarbon molecules that can form clathrate hydrates. Previous studies have reported methods for calculating the three-phase equilibrium using Monte Carlo simulation methods in systems with a single component in the gas phase. Here we extend those methods to a binary gas mixture of methane and ethane. Methane-ethane system is an interesting one in that the pure components form sII clathrate hydrate whereas a binary mixture of the two can form the sII clathrate. The phase equilibria computed from Monte Carlo simulations show a good agreement with experimental data and are also able to predict the sI-sII structural transition in the clathrate hydrate. This is attributed to the quality of the TIP4P/Ice and TRaPPE models used in the simulations. (C) 2014 Elsevier B.V. All rights reserved.

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.

Shear-imposed falling film

The study of a film falling down an inclined plane is revisited in the presence of imposed shear stress. Earlier studies regarding this topic (Smith, J. Fluid Mech., vol. 217, 1990, pp. 469-485; Wei, Phys. Fluids, vol. 17, 2005a, 012103), developed on the basis of a low Reynolds number, are extended up to moderate values of the Reynolds number. The mechanism of the primary instability is provided under the framework of a two-wave structure, which is normally a combination of kinematic and dynamic waves. In general, the primary instability appears when the kinematic wave speed exceeds the speed of dynamic waves. An equality criterion between their speeds yields the neutral stability condition. Similarly, it is revealed that the nonlinear travelling wave solutions also depend on the kinematic and dynamic wave speeds, and an equality criterion between the speeds leads to an analytical expression for the speed of a family of travelling waves as a function of the Froude number. This new analytical result is compared with numerical prediction, and an excellent agreement is achieved. Direct numerical simulations of the low-dimensional model have been performed in order to analyse the spatiotemporal behaviour of nonlinear waves by applying a constant shear stress in the upstream and downstream directions. It is noticed that the presence of imposed shear stress in the upstream (downstream) direction makes the evolution of spatially growing waves weaker (stronger).

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.

Modelling of Effect of Non-Uniform Current Density on the Performance of Soluble Lead Redox Flow Batteries

Soluble lead acid redox flow battery (SLRFB) offers a number of advantages. These advantages can be harnessed after problems associated with buildup of active material on. electrodes (residue) are resolved. A mathematical model is developed to understand residue formation in SLRFB. The model incorporates fluid flow, ion transport, electrode reactions, and non-uniform current distribution on electrode surfaces. A number of limiting cases are studied to conclude that ion transport and electrode reaction on anode simultaneously control battery performance. The model fits the reported cell voltage vs. time profiles very well. During the discharge cycle, the model predicts complete dissolution of deposited material from trailing edge side of the electrodes. With time, the active surface area of electrodes decreases rapidly. The corresponding increase in current density leads to precipitous decrease in cell potential before all the deposited material is dissolved. The successive charge-discharge cycles add to the residue. The model correctly captures the marginal effect of flow rate on cell voltage profiles, and identifies flow rate and flow direction as new variables for controlling residue buildup. Simulations carried out with alternating flow direction and a SLRFB with cylindrical electrodes show improved performance with respect to energy efficiency and residue buildup. (C) 2014 The Electrochemical Society. All rights reserved.

Oscillation dynamics of sessile droplets subjected to substrate vibration

Sessile droplets on a vibrating substrate are investigated focusing on axisymmetric oscillations with pinned contact line. Proper orthogonal decomposition is employed to identify the different modes of droplet shape oscillation and quantitatively assess the droplet oscillation and spectral response. We offer the first experimental evidence for the analogy of an oscillating sessile droplet with a non-linear spring mass damper system. The qualitative and quantitative agreement of amplitude response and phase response curves and limit cycles of the model dynamical system with that observed experimentally suggest that the bulk oscillations in the fundamental mode of a sessile droplet can be very well modeled by a Duffing oscillator with a hard spring, especially near the resonance. The red shift of the resonance peak with an increase in the glycerol concentration is clearly evidenced by both the experimental and predicted amplitude response curves. The influence of various operational parameters such as excitation frequency and amplitude and fluid properties on the droplet oscillation characteristics is adequately captured by the model. (C) 2014 Elsevier Ltd. All rights reserved.

A new microscopic insight into membrane penetration and reorganization by PETIM dendrimers

Dendrimers are highly branched polymeric nanoparticles whose structure and topology, largely, have determined their efficacy in a wide range of studies performed so far. An area of immense interest is their potential as drug and gene delivery vectors. Realizing this potential, depending on the nature of cell surface-dendrimer interactions, here we report controlled model membrane penetration and reorganization, using a model supported lipid bilayer and poly(ether imine) (PETIM) dendrimers of two generations. By systematically varying the areal density of the lipid bilayers, we provide a microscopic insight, through a combination of high resolution scattering, atomic force microscopy and atomistic molecular dynamics simulations, into the mechanism of PETIM dendrimer membrane penetration, pore formation and membrane re-organization induced by such interactions. Our work represents the first systematic observation of a regular barrel-like membrane spanning pore formation by dendrimers, tunable through lipid bilayer packing, without membrane disruption.