Experimental and analytical study of the effect of contact angle on liquid convective heat transfer in microchannels

In this paper we examine the effect of contact angle (or surface wettability) on the convective heat transfer coefficient in microchannels. Slip flow, where the fluid velocity at the wall is non-zero, is most likely to occur in microchannels due to its dependence on shear rate or wall shear stress. We show analytically that for a constant pressure drop, the presence of slip increases the Nusselt number. In a microchannel heat exchanger we modified the surface wettability from a contact angle of 20 degrees-120 degrees using thin film coating technology. Apparent slip flow is implied from pressure and flow rate measurements with a departure from classical laminar friction coefficients above a critical shear rate of approximately 10,000 s(-1). The magnitude of this departure is dependant on the contact angle with higher contact angles surfaces exhibiting larger pressure drop decreases. Similarly, the non-dimensional heat flux is found to decrease relative to laminar non-slip theory, and this decrease is also a function of the contact angle. Depending on the contact angle and the wall shear rate, variations in the heat transfer rate exceeding 10% can be expected. Thus the contact angle is an important consideration in the design of micro, and even more so, nano heat exchangers. (c) 2006 Elsevier Ltd. All rights reserved.

Friction factors and rheological behavior of Australian honey in a straight pipe

The pumping characteristics of four Australian honey samples were investigated in a straight pipe. Six flow rates (100-500 kg h(-1)) were studied at three temperatures (35-50degreesC). The pressure loss increased with an increase in the length of the pipe, as the low rate was increased and as the temperature was reduced. In the 25.4 mm-pipe, the Reynolds number ranged from 0.2-32.0 and are substantially less than the critica value (2040-2180) for laminar condition in the system. The relationship between the wall shear stress and shear rate approximated power-law behaviour, and the power-law index was not significantly (p>0.05) different from 1.0. The honey samples exhibited Newtonian behaviour at all the temperatures and this was confirmed by rheometric studies using Couette geometry. A friction chart was generated independent of temperature and the type of honey. An equation was developed to predict the pressure loss of the honey in a typical pipeline at any temperature once the viscosity-temperature relationship had been established.

Diffusion modeling of percutaneous absorption kinetics: 3. Variable diffusion and partition coefficients, consequences for stratum corneum depth profiles and desorption kinetics

Stratum corneum (SC) desorption experiments have yielded higher calculated steady-state fluxes than those obtained by epidermal penetration studies. A possible explanation of this result is a variable diffusion or partition coefficient across the SC. We therefore developed the diffusion model for percutaneous penetration and desorption to study the effects of either a variable diffusion coefficient or variable partition coefficient in the SC over the diffusion path length. Steady-state flux, lag time, and mean desorption time were obtained from Laplace domain solutions. Numerical inversion of the Laplace domain solutions was used for simulations of solute concentration-distance and amount penetrated (desorbed)-time profiles. Diffusion and partition coefficients heterogeneity were examined using six different models. The effect of heterogeneity on predicted flux from desorption studies was compared with that obtained in permeation studies. Partition coefficient heterogeneity had a more profound effect on predicted fluxes than diffusion coefficient heterogeneity. Concentration-distance profiles show even larger dependence on heterogeneity, which is consistent with experimental tape-stripping data reported for clobetasol propionate and other solutes. The clobetasol propionate tape-stripping data were most consistent with the partition coefficient decreasing exponentially for half the SC and then becoming a constant for the remaining SC. (C) 2004 Wiley-Liss, Inc.

Comparisons of diffusive and viscous contributions to transport coefficients of light gases in single-walled carbon nanotubes

We examine here the relative importance of different contributions to transport of light gases in single walled carbon nanotubes, using methane and hydrogen as examples. Transport coefficients at 298 K are determined using molecular dynamics simulation with atomistic models of the nanotube wall, from which the diffusive and viscous contributions are resolved using a recent approach that provides an explicit expression for the latter. We also exploit an exact theory for the transport of Lennard-Jones fluids at low density considering diffuse reflection at the tube wall, thereby permitting the estimation of Maxwell coefficients for the wall reflection. It is found that reflection from the carbon nanotube wall is nearly specular, as a result of which slip flow dominates, and the viscous contribution is small in comparison, even for a tube as large as 8.1 nm in diameter. The reflection coefficient for hydrogen is 3-6 times as large as that for methane in tubes of 1.36 nm diameter, indicating less specular reflection for hydrogen and greater sensitivity to atomic detail of the surface. This reconciles results showing that transport coefficients for hydrogen and methane, obtained in simulation, are comparable in tubes of this size. With increase in adsorbate density, the reflection coefficient increases, suggesting that adsorbate interactions near the wall serve to roughen the local potential energy landscape perceived by fluid molecules.

Accessing chain length dependent termination rate coefficients of methyl methacrylate (MMA) via the reversible addition fragmentation chain transfer (RAFT) process

The RAFT-CLD-T methodology is demonstrated to be not only applicable to 1-substituted monomers such as styrene and acrylates, but also to 1,1-disubstituted monomers such as MMA. The chain length of the terminating macromolecules is controlled by CPDB in MMA bulk free radical polymerization at 80 degrees C. The evolution of the chain length dependent termination rate coefficient, k(t)(i,i), was constructed in a step-wise fashion, since the MMA/CPDB system displays hybrid behavior (between conventional and living free radical polymerization) resulting in initial high molecular weight polymers formed at low RAFT agent concentrations. The obtained CLD of k(t) in MMA polymerizations is compatible with the composite model for chain length dependent termination. For the initial chain-length regime, up to a degree of polymerization of 100, k(t) decreases with alpha (in the expression k(t)(i,i) = k(t)(0) . i(-alpha)) being close to 0.65 at 80 degrees C. At chain lengths exceeding 100, the decrease is less pronounced (affording an alpha of 0.15 at 80 degrees C). However, the data are best represented by a continuously decreasing nonlinear functionality implying a chain length dependent alpha.

Control of hypersonic turbulent skin friction by boundary-layer combustion of hydrogen

Shvab-Zeldovich coupling of flow variables has been used to extend Van Driest's theory of turbulent boundary-layer skin friction to include injection and combustion of hydrogen in the boundary layer. The resulting theory is used to make predictions of skin friction and heat transfer that are found to be consistent with experimental and numerical results. Using the theory to extrapolate to larger downstream distances at the same experimental conditions, it is found that the reduction in skin-friction drag with hydrogen mixing and combustion is three times that with mixing alone. In application to flow on a flat plate at mainstream velocities of 2, 4, and 6 knits, and Reynolds numbers from 3 X 10(6) to 1 x 10(8), injection and combustion of hydrogen yielded values of skin-friction drag that were less than one-half of the no-injection skin-friction drag, together with a net reduction in heat transfer when the combustion heat release in air was less than the stagnation enthalpy. The mass efficiency of hydrogen injection, as measured by effective specific impulse values, was approximately 2000 s.

Customised birthweight: Coefficients for an Australian population and validation of the model

Background: Published birthweight references in Australia do not fully take into account constitutional factors that influence birthweight and therefore may not provide an accurate reference to identify the infant with abnormal growth. Furthermore, studies in other regions that have derived adjusted (customised) birthweight references have applied untested assumptions in the statistical modelling. Aims: To validate the customised birthweight model and to produce a reference set of coefficients for estimating a customised birthweight that may be useful for maternity care in Australia and for future research. Methods: De-identified data were extracted from the clinical database for all births at the Mater Mother's Hospital, Brisbane, Australia, between January 1997 and June 2005. Births with missing data for the variables under study were excluded. In addition the following were excluded: multiple pregnancies, births less than 37 completed week's gestation, stillbirths, and major congenital abnormalities. Multivariate analysis was undertaken. A double cross-validation procedure was used to validate the model. Results: The study of 42 206 births demonstrated that, for statistical purposes, birthweight is normally distributed. Coefficients for the derivation of customised birthweight in an Australian population were developed and the statistical model is demonstrably robust. Conclusions: This study provides empirical data as to the robustness of the model to determine customised birthweight. Further research is required to define where normal physiology ends and pathology begins, and which segments of the population should be included in the construction of a customised birthweight standard.

Negative second virial coefficients as predictors of protein crystal growth: Evidence from sedimentation equilibrium studies that refutes the designation of those light scattering parameters as osmotic virial coefficients

The effects of ammonium sulphate concentration on the osmotic second virial coefficient (B-AA/M-A) for equine serum albumin (pH 5.6, 20 degrees C) have been examined by sedimentation equilibrium. After an initial steep decrease with increasing ammonium sulphate concentration, B-AA/M-A assumes an essentially concentration-independent magnitude of 8-9 ml/g. Such behaviour conforms with the statistical-mechanical prediction that a sufficient increase in ionic strength should effectively eliminate the contributions of charge interactions to B-AA/M-A but have no effect on the covolume contribution (8.4 ml/g for serum albumin). A similar situation is shown to apply to published sedimentation equilibrium data for lysozyme (pH 4.5). Although termed osmotic second virial coefficients and designated as such (B-22), the negative values obtained in published light scattering studies of both systems have been described incorrectly because of the concomitant inclusion of the protein-salt contribution to thermodynamic nonideality of the protein. Those negative values are still valid predictors of conditions conducive to crystal growth inasmuch as they do reflect situations in which there is net attraction between protein molecules. However, the source of attraction responsible for the negative virial coefficient stems from the protein-salt rather than the protein-protein contribution, which is necessarily positive. (c) 2005 Elsevier B.V. All rights reserved.

Scramjet lift, thrust and pitching-moment characteristics measured in a shock tunnel

Lift, pitching moment, and thrust/drag on a supersonic combustion ramjet were measured in the T4 free-piston shock tunnel using a three-component stress-wave force balance. The scramjet model was 0.567 m long and weighed approximately 6 kg. Combustion occurred at a nozzle-supply enthalpy of 3.3 MJ/kg and nozzle-supply pressure of 32 MPa at Mach 6.6 for equivalence ratios up to 1.4. The force coefficients varied approximately linearly with equivalence ratio. The location of the center of pressure changed by 10% of the chord of the model over the range of equivalence ratios tested. Lift and pitching-moment coefficients remained constant when the nozzle-supply enthalpy was increased to 4.9 MJ/kg at an equivalence ratio of 0.8, but the thrust coefficient decreased rapidly. When the nozzle-supply pressure was reduced at a nozzle-supply enthalpy of 3.3 MJ/kg and an equivalence ratio of 0.8, the combustion-generated increment of lift and thrust was maintained at 26 MPa, but disappeared at 16 MPa. Measured lift and thrust forces agreed well with calculations made using a simplified force prediction model, but the measured pitching moment substantially exceeded predictions. Choking occurred at nozzle-supply enthalpies of less than 3.0 MJ/kg with an equivalence ratio of 0.8. The tests failed to yield a positive thrust because of the skin-friction drag that accounted for up to 50% of the fuel-off drag.

Testing predictions of macroscopic binary diffusion coefficients using lattice models with site heterogeneity

Quantitatively predicting mass transport rates for chemical mixtures in porous materials is important in applications of materials such as adsorbents, membranes, and catalysts. Because directly assessing mixture transport experimentally is challenging, theoretical models that can predict mixture diffusion coefficients using Only single-component information would have many uses. One such model was proposed by Skoulidas, Sholl, and Krishna (Langmuir, 2003, 19, 7977), and applications of this model to a variety of chemical mixtures in nanoporous materials have yielded promising results. In this paper, the accuracy of this model for predicting mixture diffusion coefficients in materials that exhibit a heterogeneous distribution of local binding energies is examined. To examine this issue, single-component and binary mixture diffusion coefficients are computed using kinetic Monte Carlo for a two-dimensional lattice model over a wide range of lattice occupancies and compositions. The approach suggested by Skoulidas, Sholl, and Krishna is found to be accurate in situations where the spatial distribution of binding site energies is relatively homogeneous, but is considerably less accurate for strongly heterogeneous energy distributions.

A unified friction description and its application to the simulation of frictional instability using the finite element method

This article first summarizes some available experimental results on the frictional behaviour of contact interfaces, and briefly recalls typical frictional experiments and relationships, which are applicable for rock mechanics, and then a unified description is obtained to describe the entire frictional behaviour. It is formulated based on the experimental results and applied with a stick and slip decomposition algorithm to describe the stick-slip instability phenomena, which can describe the effects observed in rock experiments without using the so-called state variable, thus avoiding related numerical difficulties. This has been implemented to our finite element code, which uses the node-to-point contact element strategy proposed by the authors to handle the frictional contact between multiple finite-deformation bodies with stick and finite frictional slip, and applied here to simulate the frictional behaviour of rocks to show its usefulness and efficiency.

The lattice solid model to simulate the physics of rocks and earthquakes: Incorporation of friction

The particle-based lattice solid model developed to study the physics of rocks and the nonlinear dynamics of earthquakes is refined by incorporating intrinsic friction between particles. The model provides a means for studying the causes of seismic wave attenuation, as well as frictional heat generation, fault zone evolution, and localisation phenomena. A modified velocity-Verlat scheme that allows friction to be precisely modelled is developed. This is a difficult computational problem given that a discontinuity must be accurately simulated by the numerical approach (i.e., the transition from static to dynamical frictional behaviour). This is achieved using a half time step integration scheme. At each half time step, a nonlinear system is solved to compute the static frictional forces and states of touching particle-pairs. Improved efficiency is achieved by adaptively adjusting the time step increment, depending on the particle velocities in the system. The total energy is calculated and verified to remain constant to a high precision during simulations. Numerical experiments show that the model can be applied to the study of earthquake dynamics, the stick-slip instability, heat generation, and fault zone evolution. Such experiments may lead to a conclusive resolution of the heat flow paradox and improved understanding of earthquake precursory phenomena and dynamics. (C) 1999 Academic Press.