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.

Influence of HZSM-5 catalyst on the thermal degradation of poly(vinyl chloride) in solution

The kinetics of thermal degradation of poly(vinyl chloride) (PVC) in solution was investigated at various temperatures (210-250degreesC). The degradation rate coefficients were determined from the time evolution of the molecular weight distribution (MWD). The energy of activation, determined from the temperature dependence of the rate coefficient, was 26.6 kcal/mol. The degradation of PVC was also studied in the presence of a catalyst (HZSM-5 zeolite). The results indicated that increase of the degradation rate of PVC is first order with the HZSM-5 concentration up to 50 g/L and zero order at higher concentrations. The thermal degradation kinetics of PVC in the presence of 50 g/L of the catalyst was studied at various temperatures. The temperature dependency of the rate coefficient was used to calculate the activation energy (21.5 kcal/mol). This is consistent with the observation that the presence of a catalyst generally decreases the activation energy and promotes degradation. (C) 2002 John Wiley Sons, Inc.

Kinetics of thermal degradation of vinyl polyperoxides in solution

The thermal degradation of vinyl polyperoxides, poly(styrene peroxide, (PSP), poly(alpha-methylstyrene peroxide) (PAMSP) and poly(alpha-phenylstyrene pet-oxide) (PAPSP), was carried out at different temperatures in toluene. The time evolution of molecular weight distributions (MWDs) was determined by gel permeation chromatography (GPC). A continuous distribution model was used to evaluate the random chain degradation rate coefficients. The activation energies, determined from the temperature dependence of the rate coefficients, suggest that thermal degradation of polyperoxides is controlled by the dissociation of the O-O bonds in the backbone of the polymer chain. Among the three polyperoxides investigated, the thermal stability is the highest for PAPSP and the lowest for PAMSP. (C) 2002 Elsevier Science Ltd. All rights reserved.

Mechanical properties of Nylon-6/LDPE blends containing an epoxy resin compatibilizer

Binary and ternary blends of nylon-6/low density polyethylene (nylon-6/LDPE) and Nylon-6/LDPE/poly(ethylene-co-glycidyl methacrylate) were prepared by melt mixing. The blends exhibit two phase morphology with LDPE dispersed in the form of spherical domains in the nylon-6 matrix. The mechanical properties of the blends were measured by standard methods. It is shown that the use of the epoxy copolymer as a compatibilizer improves the impact strength of the blend as compared to nylon-6, which is attributed to better stress transfer across the interface due to the compatibilizer. The data for each mechanical property were also fitted into a best fit model equation and the method of steepest ascent was applied to arrive at the optimum composition of the blend for that property.

Microwave heating in multiphase systems: evaluation of series solutions

The 1D electric field and heat-conduction equations are solved for a slab where the dielectric properties vary spatially in the sample. Series solutions to the electric field are obtained for systems where the spatial variation in the dielectric properties can be expressed as polynomials. The series solution is used to obtain electric-field distributions for a binary oil-water system where the dielectric properties are assumed to vary linearly within the sample. Using the finite-element method temperature distributions are computed in a three-phase oil, water and rock system where the dielectric properties vary due to the changing oil saturation in the rock. Temperature distributions predicted using a linear variation in the dielectric properties are compared with those obtained using the exact nonlinear variation.

Nonlinear dynamics of eucaryotic pyruvate dehydrogenase multienzyme complex: Decarboxylation rate, oscillations, and multiplicity

Pyruvate conversion to acetyl-CoA by the pyruvate dehydrogenase (PDH) multienzyme complex is known as a key node in affecting the metabolic fluxes of animal cell culture. However, its possible role in causing possible nonlinear dynamic behavior such as oscillations and multiplicity of animal cells has received little attention. In this work, the kinetic and dynamic behavior of PDH of eucaryotic cells has been analyzed by using both in vitro and simplified in vivo models. With the in vitro model the overall reaction rate (v(1)) of PDH is shown to be a nonlinear function of pyruvate concentration, leading to oscillations under certain conditions. All enzyme components affect v, and the nonlinearity of PDH significantly, the protein X and the core enzyme dihydrolipoamide acyltransferase (E2) being mostly predominant. By considering the synthesis rates of pyruvate and PDH components the in vitro model is expanded to emulate in vivo conditions. Analysis using the in vivo model reveals another interesting kinetic feature of the PDH system, namely, multiple steady states. Depending on the pyruvate and enzyme levels or the operation mode, either a steady state with high pyruvate decarboxylation rate or a steady state with significantly lower decarboxylation rate can be achieved under otherwise identical conditions. In general, the more efficient steady state is associated with a lower pyruvate concentration. A possible time delay in the substrate supply and enzyme synthesis can also affect the steady state to be achieved and lead's to oscillations under certain conditions. Overall, the predictions of multiplicity for the PDH system agree qualitatively well with recent experimental observations in animal cell cultures. The model analysis gives some hints for improving pyruavte metabolism in animal cell culture.

Adsorption of 1,1,1,2-Tetrafluoroethane on Activated Charcoal

This paper presents adsorption isotherms for HFC-134a on activated charcoal, in the temperature range of 273-353 K and for pressures up to 0.65 MPa, measured using the volumetric method. Three samples of charcoals with widely varying surface areas were chosen. The shapes of the isotherms,obtained from the experimental data were similar in all cases and comparable to those reported in the literature. Adsorption parameters were evaluated from the isotherms using the Dubinin-Astakhov (DA) equation. The concentration dependence of the isosteric enthalpies of adsorption is extracted from the data.

Role of length scales on microwave thawing dynamics in 2D cylinders

Microwave (MW) thawing of 2D frozen cylinders exposed to uniform plane waves from one face, is modeled using the effective heat capacity formulation with the MW power obtained from the electric field equations. Computations are illustrated for tylose (23% methyl cellulose gel) which melts over a range of temperatures giving rise to a mushy zone. Within the mushy region the dielectric properties are functions of the liquid volume fraction. The resulting coupled, time dependent non-linear equations are solved using the Galerkin finite element method with a fixed mesh. Our method efficiently captures the multiple connected thawed domains that arise due to the penetration of MWs in the sample. For a cylinder of diameter D, the two length scales that control the thawing dynamics are D/D-p and D/lambda(m), where D-p and lambda(m) are the penetration depth and wavelength of radiation in the sample respectively. For D/D-p, D/lambda(m) much less than 1 power absorption is uniform and thawing occurs almost simultaneously across the sample (Regime I). For D/D-p much greater than 1 thawing is seen to occur from the incident face, since the power decays exponentially into the sample (Regime III). At intermediate values, 0.2 < D/D-p, D/lambda(m) < 2.0 (Regime II) thawing occurs from the unexposed face at smaller diameters, from both faces at intermediate diameters and from the exposed and central regions at larger diameters. Average power absorption during thawing indicates a monotonic rise in Regime I and a monotonic decrease in Regime III. Local maxima in the average power observed for samples in Regime II are due to internal resonances within the sample. Thawing time increases monotonically with sample diameter and temperature gradients in the sample generally increase from Regime I to Regime III. (C) 2002 Elsevier Science Ltd. All rights reserved.

Stability of fluid flow past a membrane

The stability of fluid flow past a membrane of infinitesimal thickness is analysed in the limit of zero Reynolds number using linear and weakly nonlinear analyses. The system consists of two Newtonian fluids of thickness R* and H R*, separated by an infinitesimally thick membrane, which is flat in the unperturbed state. The dynamics of the membrane is described by its normal displacement from the flat state, as well as a surface displacement field which provides the displacement of material points from their steady-state positions due to the tangential stress exerted by the fluid flow. The surface stress in the membrane (force per unit length) contains an elastic component proportional to the strain along the surface of the membrane, and a viscous component proportional to the strain rate. The linear analysis reveals that the fluctuations become unstable in the long-wave (alpha --> 0) limit when the non-dimensional strain rate in the fluid exceeds a critical value Lambda(t), and this critical value increases proportional to alpha(2) in this limit. Here, alpha is the dimensionless wavenumber of the perturbations scaled by the inverse of the fluid thickness R*(-1), and the dimensionless strain rate is given by Lambda(t) = ((gamma) over dot* R*eta*/Gamma*), where eta* is the fluid viscosity, Gamma* is the tension of the membrane and (gamma) over dot* is the strain rate in the fluid. The weakly nonlinear stability analysis shows that perturbations are supercritically stable in the alpha --> 0 limit.