Solution grafting of maleic anhydride onto polyethylene and compatibilization of polyetylene-nylon blends

Maleic anhydride (MAH) has been grafted onto high density polyethylene (HDPE) with benzoyl peroxide (BOP) initiator in toluene solution. Maximum degree of grafting (12%) without crosslinking has been obtained using MAH/HDPE and BOP/HDPE weight ratios of 1.0 and 0.15 respectively, at 110 degrees C. The HDPE-g-MAH compatibilizer is found to drastically reduce the dispersed phase size and also to produce homogeneous blends for relatively low concentrations of dispersed phase in HDPE/nylon blends. Addition of this compatibilizer results in increase of tensile strength and modulus with increasing nylon content of HDPE/nylon blends, while the opposite is found for the blends without any added compatibilizer.

Stability of non-parabolic flow in a flexible tube

Flows with velocity profiles very different from the parabolic velocity profile can occur in the entrance region of a tube as well as in tubes with converging/diverging cross-sections. In this paper, asymptotic and numerical studies are undertaken to analyse the temporal stability of such 'non-parabolic' flows in a flexible tube in the limit of high Reynolds numbers. Two specific cases are considered: (i) developing flow in a flexible tube; (ii) flow in a slightly converging flexible tube. Though the mean velocity profile contains both axial and radial components, the flow is assumed to be locally parallel in the stability analysis. The fluid is Newtonian and incompressible, while the flexible wall is modelled as a viscoelastic solid. A high Reynolds number asymptotic analysis shows that the non-parabolic velocity profiles can become unstable in the inviscid limit. This inviscid instability is qualitatively different from that observed in previous studies on the stability of parabolic flow in a flexible tube, and from the instability of developing flow in a rigid tube. The results of the asymptotic analysis are extended numerically to the moderate Reynolds number regime. The numerical results reveal that the developing flow could be unstable at much lower Reynolds numbers than the parabolic flow, and hence this instability can be important in destabilizing the fluid flow through flexible tubes at moderate and high Reynolds number. For flow in a slightly converging tube, even small deviations from the parabolic profile are found to be sufficient for the present instability mechanism to be operative. The dominant non-parallel effects are incorporated using an asymptotic analysis, and this indicates that non-parallel effects do not significantly affect the neutral stability curves. The viscosity of the wall medium is found to have a stabilizing effect on this instability.

Statistical thermodynamics of lattice models in zeolites: Implications of local versus global mean field interactions

The statistical thermodynamics of adsorption in caged zeolites is developed by treating the zeolite as an ensemble of M identical cages or subsystems. Within each cage adsorption is assumed to occur onto a lattice of n identical sites. Expressions for the average occupancy per cage are obtained by minimizing the Helmholtz free energy in the canonical ensemble subject to the constraints of constant M and constant number of adsorbates N. Adsorbate-adsorbate interactions in the Brag-Williams or mean field approximation are treated in two ways. The local mean field approximation (LMFA) is based on the local cage occupancy and the global mean field approximation (GMFA) is based on the average coverage of the ensemble. The GMFA is shown to be equivalent in formulation to treating the zeolite as a collection of interacting single site subsystems. In contrast, the treatment in the LMFA retains the description of the zeolite as an ensemble of identical cages, whose thermodynamic properties are conveniently derived in the grand canonical ensemble. For a z coordinated lattice within the zeolite cage, with epsilon(aa) as the adsorbate-adsorbate interaction parameter, the comparisons for different values of epsilon(aa)(*)=epsilon(aa)z/2kT, and number of sites per cage, n, illustrate that for -1 0. We compare the isotherms predicted with the LMFA with previous GMFA predictions [K. G. Ayappa, C. R. Kamala, and T. A. Abinandanan, J. Chem. Phys. 110, 8714 (1999)] (which incorporates both the site volume reduction and a coverage-dependent epsilon(aa)) for xenon and methane in zeolite NaA. In all cases the predicted isotherms are very similar, with the exception of a small steplike feature present in the LMFA for xenon at higher coverages. (C) 1999 American Institute of Physics. [S0021-9606(99)70333-8].

Temperature scaling in a dense vibrofluidized granular material

The leading order "temperature" of a dense two-dimensional granular material fluidized by external vibrations is determined. The grain interactions are characterized by inelastic collisions, but the coefficient of restitution is considered to be close to 1, so that the dissipation of energy during a collision is small compared to the average energy of a particle. An asymptotic solution is obtained where the particles are considered to be elastic in the leading approximation. The velocity distribution is a Maxwell-Boltzmann distribution in the leading approximation,. The density profile is determined by solving the momentum balance equation in the vertical direction, where the relation between the pressure and density is provided by the virial equation of state. The temperature is determined by relating the source of energy due to the vibrating surface and the energy dissipation due to inelastic collisions. The predictions of the present analysis show good agreement with simulation results at higher densities where theories for a dilute vibrated granular material, with the pressure-density relation provided by the ideal gas law, sire in error. [:S1063-651X(99)04408-6].

The effect of boundaries on the plane Couette flow of granular materials: a bifurcation analysis

The tendency of granular materials in rapid shear flow to form non-uniform structures is well documented in the literature. Through a linear stability analysis of the solution of continuum equations for rapid shear flow of a uniform granular material, performed by Savage (1992) and others subsequently, it has been shown that an infinite plane shearing motion may be unstable in the Lyapunov sense, provided the mean volume fraction of particles is above a critical value. This instability leads to the formation of alternating layers of high and low particle concentrations oriented parallel to the plane of shear. Computer simulations, on the other hand, reveal that non-uniform structures are possible even when the mean volume fraction of particles is small. In the present study, we have examined the structure of fully developed layered solutions, by making use of numerical continuation techniques and bifurcation theory. It is shown that the continuum equations do predict the existence of layered solutions of high amplitude even when the uniform state is linearly stable. An analysis of the effect of bounding walls on the bifurcation structure reveals that the nature of the wall boundary conditions plays a pivotal role in selecting that branch of non-uniform solutions which emerges as the primary branch. This demonstrates unequivocally that the results on the stability of bounded shear how of granular materials presented previously by Wang et al. (1996) are, in general, based on erroneous base states.

Synthesis and Degradation of Sorbitol-Based Polymers

A new class of biodegradable copolyesters was synthesized by the catalyst-free melt condensation of sorbitol with citric acid, tartaric acid, and sebacic acid. The resulting polymers were designated as poly(sorbitol citric sebacate) p(SCS)] and poly(sorbitol tartaric sebacate) p(STS)]. The synthesized polymers were characterized by Fourier transform infrared spectroscopy, H-1-NMR spectroscopy, and differential scanning calorimetry analysis. Porous spongelike scaffolds were prepared with a salt-leaching technique and characterized with scanning electron microscopy. Tensile testing of the p(SCS) and p(STS) polymers showed that they exhibited a wide range of mechanical properties. The Young's modulus and tensile strengths of the polymers ranged from 1.06 +/- 0.12 to 462.65 +/- 34.21 MPa and from 0.45 +/- 0.04 to 20.32 +/- 2.54 MPa, respectively. In vitro degradation studies were performed on disc-shaped polymer samples. The half-life of the polymers ranged from 0.54 to 38.52 days. The percentage hydration of the polymers was in the range 9.36 +/- 1.26 to 78.25 +/- 1.91, with sol contents of 2-14%. At any given polymer composition, the Young's modulus and tensile strength of p(SCS) was higher than that of p(STS), whereas the degradation rates of p(SCS) was lower than that of p(STS). This was attributed to the structural difference between the citric and tartaric monomers and to the degree of crosslinking. (C) 2011 Wiley Periodicals, Inc. J Appl Polym Sci 121: 2861-2869, 2011

Effect of Zr4+-ion substitution in CeO2 on H2O2-assisted degradation of orange G

A series of Pd ion-substituted CeO2-ZrO2 solid solutions were synthesized using the solution combustion technique. H2O2-assisted degradation of orange G was carried out in the presence of the catalysts. The activity of the catalysts was found to increase with the introduction of the second component in the solid solution, as signified by an increase in the rate constants and lowering of activation energy. The study showed the involvement of lattice oxygen and the importance of reducibility of the compound for the reaction. (C) 2011 Elsevier B.V. All rights reserved.

Velocity distribution for a dilute vibrated granular material

The velocity distribution for a vibrated granular material is determined in the dilute limit where the frequency of particle collisions with the vibrating surface is large compared to the frequency of binary collisions. The particle motion is driven by the source of energy due to particle collisions with the vibrating surface, and two dissipation mechanisms-inelastic collisions and air drag-are considered. In the latter case, a general form for the drag force is assumed. First, the distribution function for the vertical velocity for a single particle colliding with a vibrating surface is determined in the limit where the dissipation during a collision due to inelasticity or between successive collisions due to drag is small compared to the energy of a particle. In addition, two types of amplitude functions for the velocity of the surface, symmetric and asymmetric about zero velocity, are considered. In all cases, differential equations for the distribution of velocities at the vibrating surface are obtained using a flux balance condition in velocity space, and these are solved to determine the distribution function. It is found that the distribution function is a Gaussian distribution when the dissipation is due to inelastic collisions and the amplitude function is symmetric, and the mean square velocity scales as [[U-2](s)/(1 - e(2))], where [U-2](s) is the mean square velocity of the vibrating surface and e is the coefficient of restitution. The distribution function is very different from a Gaussian when the dissipation is due to air drag and the amplitude function is symmetric, and the mean square velocity scales as ([U-2](s)g/mu(m))(1/(m+2)) when the acceleration due to the fluid drag is -mu(m)u(y)\u(y)\(m-1), where g is the acceleration due to gravity. For an asymmetric amplitude function, the distribution function at the vibrating surface is found to be sharply peaked around [+/-2[U](s)/(1-e)] when the dissipation is due to inelastic collisions, and around +/-[(m +2)[U](s)g/mu(m)](1/(m+1)) when the dissipation is due to fluid drag, where [U](s) is the mean velocity of the surface. The distribution functions are compared with numerical simulations of a particle colliding with a vibrating surface, and excellent agreement is found with no adjustable parameters. The distribution function for a two-dimensional vibrated granular material that includes the first effect of binary collisions is determined for the system with dissipation due to inelastic collisions and the amplitude function for the velocity of the vibrating surface is symmetric in the limit delta(I)=(2nr)/(1 - e)much less than 1. Here, n is the number of particles per unit width and r is the particle radius. In this Limit, an asymptotic analysis is used about the Limit where there are no binary collisions. It is found that the distribution function has a power-law divergence proportional to \u(x)\((c delta l-1)) in the limit u(x)-->0, where u(x) is the horizontal velocity. The constant c and the moments of the distribution function are evaluated from the conservation equation in velocity space. It is found that the mean square velocity in the horizontal direction scales as O(delta(I)T), and the nontrivial third moments of the velocity distribution scale as O(delta(I)epsilon(I)T(3/2)) where epsilon(I) = (1 - e)(1/2). Here, T = [2[U2](s)/(1 - e)] is the mean square velocity of the particles.

A fixed-grid finite element based enthalpy formulation for generalized phase change problems: role of superficial mushy region

The enthalpy method is primarily developed for studying phase change in a multicomponent material, characterized by a continuous liquid volume fraction (phi(1)) vs temperature (T) relationship. Using the Galerkin finite element method we obtain solutions to the enthalpy formulation for phase change in 1D slabs of pure material, by assuming a superficial phase change region (linear (phi(1) vs T) around the discontinuity at the melting point. Errors between the computed and analytical solutions are evaluated for the fluxes at, and positions of, the freezing front, for different widths of the superficial phase change region and spatial discretizations with linear and quadratic basis functions. For Stefan number (St) varying between 0.1 and 10 the method is relatively insensitive to spatial discretization and widths of the superficial phase change region. Greater sensitivity is observed at St = 0.01, where the variation in the enthalpy is large. In general the width of the superficial phase change region should span at least 2-3 Gauss quadrature points for the enthalpy to be computed accurately. The method is applied to study conventional melting of slabs of frozen brine and ice. Regardless of the forms for the phi(1) vs T relationships, the thawing times were found to scale as the square of the slab thickness. The ability of the method to efficiently capture multiple thawing fronts which may originate at any spatial location within the sample, is illustrated with the microwave thawing of slabs and 2D cylinders. (C) 2002 Elsevier Science Ltd. All rights reserved.

Bubble size measurement in a gas-liquid ejector for a sodium chloride-air system

A high speed photographic technique has been employed to measure the Sauter mean diameter of bubbles experimentally in a gas liquid ejector using a sodium chloride-air system. The measured values are compared with the theoretically predicted maximum bubble size diameter using Sprow's correlation. Bubble size as a function of the liquid flow rate and also of its distance from the throat of the ejector has been reported in this paper. The results obtained for this non-reactive system are also compared with those obtained earlier for the air-water system.

Continuous distribution kinetics for the thermal degradation of LDPE in solution

Polymer degradation in solution has several advantages over melt pyrolysis, The degradation of low-density polyethylene (LDPE) occurs at much lower temperatures in solution (280-360degreesC) than in conventional melt pyrolysis (400-450degreesC). The thermal degradation kinetics of LDPE in solution was investigated in this work. LDPE was dissolved in liquid paraffin and degraded for 3 h at various temperatures (280-360degreesC). Samples were taken at specific times and analyzed with high-pressure liquid chromatography/gel permeation chromatography for the molecular weight distribution (MWD), The time evolution of the MWD was modeled with continuous distribution kinetics. Data indicated that LDPE followed random-chain-scission degradation. The rapid initial drop in molecular weight, observed up to 45 min, was attributed to the presence of weak links in the polymer. The rate coefficients for the breakage of weak and strong links were determined, and the corresponding average activation energies were calculated to be 88 and 24 kJ/mol, respectively. (C) 2002 John Wiley Sons, Inc.

Boundary conditions at a rigid wall for rough granular gases

We derive boundary conditions at a rigid wall for a granular material comprising rough, inelastic particles. Our analysis is confined to the rapid flow, or granular gas, regime in which grains interact by impulsive collisions. We use the Chapman-Enskog expansion in the kinetic theory of dense gases, extended for inelastic and rough particles, to determine the relevant fluxes to the wall. As in previous studies, we assume that the particles are spheres, and that the wall is corrugated by hemispheres rigidly attached to it. Collisions between the particles and the wall hemispheres are characterized by coefficients of restitution and roughness. We derive boundary conditions for the two limiting cases of nearly smooth and nearly perfectly rough spheres, as a hydrodynamic description of granular gases comprising rough spheres is appropriate only in these limits. The results are illustrated by applying the equations of motion and boundary conditions to the problem of plane Couette flow.

Dynamics of crystal size distributions with size-dependent rates

The growth and dissolution dynamics of nonequilibrium crystal size distributions (CSDs) can be determined by solving the governing population balance equations (PBEs) representing reversible addition or dissociation. New PBEs are considered that intrinsically incorporate growth dispersion and yield complete CSDs. We present two approaches to solving the PBEs, a moment method and a numerical scheme. The results of the numerical scheme agree with the moment technique, which can be solved exactly when powers on mass-dependent growth and dissolution rate coefficients are either zero or one. The numerical scheme is more general and can be applied when the powers of the rate coefficients are non-integers or greater than unity. The influence of the size dependent rates on the time variation of the CSDs indicates that as equilibrium is approached, the CSDs become narrow when the exponent on the growth rate is less than the exponent on the dissolution rate. If the exponent on the growth rate is greater than the exponent on the dissolution rate, then the polydispersity continues to broaden. The computation method applies for crystals large enough that interfacial stability issues, such as ripening, can be neglected. (C) 2002 Elsevier Science B.V. All rights reserved.

Effect of convection on domain growth during demixing transitions in fluids

Current analytical work on the effect of convection on the late stages of spinodal decomposition in liquids is briefly described. The morphology formed during the spinodal decomposition process depends on the relative composition of the two species. Droplet spinodal decomposition occurs when the concentration of one of the species is small. Convective transport has a significant effect on the scaling laws in the late-stage coarsening of droplets in translational or shear flows. In addition, convective transport could result in an attractive interaction between non-Brownian droplets which could lead to coalescence. The effect of convective transport for the growth of random interfaces in a near-symmetric quench was analysed using an area distribution function, which gives the distribution of surface area of the interface in curvature space. It was found that the curvature of the interface decreases proportional to time t in the late stages of spinodal decomposition, and the surface area also decreases proportional to t.

Enzymatic synthesis of flavors in supercritical carbon dioxide

Commercially important flavor esters of isoamyl alcohol, catalyzed by crude hog pancreas lipase (HPL), were synthesized under solvent-free conditions and in supercritical carbon dioxide. The esters synthesized were isoamyl acetate, isoamyl propionate, isoamyl butyrate, and isoamyl octanoate. Very low yields (3-4%) of isoamyl acetate were obtained, but high yields for the other three esters were obtained under both supercritical and solvent-free conditions. The yields of esters of the even-carbon acids, isoamyl acetate, butyrate, and octanoate, increased with increasing chain length, whereas the yield of isoamyl propionate was higher than that of isoamyl butyrate. The optimum temperature of the reaction was higher under supercritical conditions (45 degreesC) than under solvent-free conditions (35-40 degreesC). The effects of other parameters such as alcohol concentration, water concentration, and enzyme loading were investigated. An increase in the water concentration decreased the conversion significantly in supercritical carbon dioxide but not under solvent-free conditions. The optimum ratio of alcohol to acid was dependent on the extent of inhibition by the acid. Although providing a higher apparent yield by being run in a highly concentrated medium, the overall conversion under solvent-free conditions was lower than that under supercritical conditions for similar enzyme concentrations, indicating that the synthesis of esters in supercritical carbon dioxide might be a viable option.