Laminar heat transfer to non-Newtonian fluids in pipes with abrupt changes in diameter

In this thesis the results of experimental work performed to determine local heat transfer coefficients for non-Newtonian fluids in laminar flow through pipes with abrupt discontinuities are reported. The fluids investigated were water-based polymeric solutiorrs of time-indpendent, pseudoplastic materials, with flow indices "n" ranging from 0.39 to 0.9.The tube configurations were a 3.3 :1 sudden convergence, and a 1: 3.3 sudden divergence.The condition of a prescribed uniform wall heat flux was considered, with both upstream and downstream tube sections heated. Radial temperature traverses were also under­ taken primarily to justify the procedures used in estimating the tube wall and bulk fluid temperatures and secondly to give further insight into the mechanism of heat transfer beyond a sudden tube expansion. A theoretical assessment of the influence of viscous dissipation on a non-Newtonian pseudoplastic fluid of' arbitrary index "n" was carried out. The effects of other secondary factors such as free convection and temperature-dependent consistency were evaluated empirically. In the present investigations, the test conditions were chosen to minimise the effects of natural convection and the estimates of viscous heat generation showed the effect to be insignificant with the polymeric concentrations tested here. The final results have been presented as the relationships between local heat transfer coef'ficient and axial distance downstream of the discontinuities and relationships between dimensionless wall temperature and reduced radius. The influence of Reynolds number, Prandtl number, non-Newtonian index and heat flux have been indicated.

Investigating viscous fluid flow in an internal mixer using computational fluid dynamics

This thesis presents an effective methodology for the generation of a simulation which can be used to increase the understanding of viscous fluid processing equipment and aid in their development, design and optimisation. The Hampden RAPRA Torque Rheometer internal batch twin rotor mixer has been simulated with a view to establishing model accuracies, limitations, practicalities and uses. As this research progressed, via the analyses several 'snap-shot' analysis of several rotor configurations using the commercial code Polyflow, it was evident that the model was of some worth and its predictions are in good agreement with the validation experiments, however, several major restrictions were identified. These included poor element form, high man-hour requirements for the construction of each geometry and the absence of the transient term in these models. All, or at least some, of these limitations apply to the numerous attempts to model internal mixes by other researchers and it was clear that there was no generally accepted methodology to provide a practical three-dimensional model which has been adequately validated. This research, unlike others, presents a full complex three-dimensional, transient, non-isothermal, generalised non-Newtonian simulation with wall slip which overcomes these limitations using unmatched ridding and sliding mesh technology adapted from CFX codes. This method yields good element form and, since only one geometry has to be constructed to represent the entire rotor cycle, is extremely beneficial for detailed flow field analysis when used in conjunction with user defined programmes and automatic geometry parameterisation (AGP), and improves accuracy for investigating equipment design and operation conditions. Model validation has been identified as an area which has been neglected by other researchers in this field, especially for time dependent geometries, and has been rigorously pursued in terms of qualitative and quantitative velocity vector analysis of the isothermal, full fill mixing of generalised non-Newtonian fluids, as well as torque comparison, with a relatively high degree of success. This indicates that CFD models of this type can be accurate and perhaps have not been validated to this extent previously because of the inherent difficulties arising from most real processes.

The atomisation of detergent slurries

A study was made to determine the conditions under which the optimum droplet size distribution (ie., narrowest size range with a minimum of fines and over-sized agglomerates), is generated in sprays from centrifugal pressure nozzles. A range of non-Newtonian detergent slurries were tested but the results are of wider application and parallel work was undertaken with water, ionic solutions and chalk slurries. Six centrifugal pressure nozzles were used and the drop-size distributions correlated as a function of fluid properties, pressure, fiowrate, feed temperature, and nozzle characteristics. Measurements were made using a Malvern Particle Size Anayser slung across a specially-designed transparent tower section of approximately 1.2m diameter in order to reduce obscuration caused by the spray and improve existing droplet sizing techniques. The results obtained were based upon the Rosin-Rammler distribution model and the Size Analyser provided a direct print-out of size distribution and the parameters characterising it. A Spraying System nozzle, AAASSTC8-8, produced the optimum spray distribution with the detergent slurry at a temperature of 60°C whilst operating at 1200 psi. With other fluids the Delevan 2.2SJ nozzle produced the optimum spray distribution operating at 1200 psi but with the Spraying Systems nozzles there was no clear-cut optimum set of conditions, ie. the nozzle and pressure varied depending upon the fluid being sprayed. The mechanisms of liquid sheet break-up and droplet dispersion were investigated in specially-constructed, scaled-up, transparent nozzles. A mathematical model of centrifugal pressure nozzle atomisation was developed based upon fundamental operating parameters rather than resorting to empirical correlations. This enabled theoretical predictions to be made over a wide range of operating conditions and nozzle types. The model predictions for volumetric fiowrate, liquid sheet length and air core diameter showed good agreement with the experimentally determined results. However, the model predicted smaller droplet sizes than were produced experimentally due to inaccuracies identified in the initial assumptions.

Bifurcation study for a vertical channel with constant flux and large aspect ratio

Using suitable coupled Navier-Stokes Equations for an incompressible Newtonian fluid we investigate the linear and non-linear steady state solutions for both a homogeneously and a laterally heated fluid with finite Prandtl Number (Pr=7) in the vertical orientation of the channel. Both models are studied within the Large Aspect Ratio narrow-gap and under constant flux conditions with the channel closed. We use direct numerics to identify the linear stability criterion in parametric terms as a function of Grashof Number (Gr) and streamwise infinitesimal perturbation wavenumber (making use of the generalised Squire’s Theorem). We find higher harmonic solutions at lower wavenumbers with a resonance of 1:3exist, for both of the heating models considered. We proceed to identify 2D secondary steady state solutions, which bifurcate from the laminar state. Our studies show that 2D solutions are found not to exist in certain regions of the pure manifold, where we find that 1:3 resonant mode 2D solutions exist, for low wavenumber perturbations. For the homogeneously heated fluid, we notice a jump phenomenon existing between the pure and resonant mode secondary solutions for very specific wavenumbers .We attempt to verify whether mixed mode solutions are present for this model by considering the laterally heated model with the same geometry. We find mixed mode solutions for the laterally heated model showing that a bridge exists between the pure and 1:3 resonant mode 2D solutions, of which some are stationary and some travelling. Further, we show that for the homogeneously heated fluid that the 2D solutions bifurcate in hopf bifurcations and there exists a manifold where the 2D solutions are stable to Eckhaus criterion, within this manifold we proceed to identify 3D tertiary solutions and find that the stability for said 3D bifurcations is not phase locked to the 2D state. For the homogeneously heated model we identify a closed loop within the neutral stability curve for higher perturbation wavenumubers and analyse the nature of the multiple 2D bifurcations around this loop for identical wavenumber and find that a temperature inversion occurs within this loop. We conclude that for a homogeneously heated fluid it is possible to have abrup ttransitions between the pure and resonant 2D solutions, and that for the laterally heated model there exist a transient bifurcation via mixed mode solutions.

Development and evaluation of a novel intranasal spray for the delivery of amantadine

The aim of this study was to develop and characterize an intranasal delivery system for amantadine hydrochloride (AMT). Optimal formulations consisted of a thermosensitive polymer Pluronic® 127 and either carboxymethyl cellulose or chitosan which demonstrated gel transition at nasal cavity temperatures (34 ± 1°C). Rheologically, the loss tangent (Tan δ) confirmed a 3-stage gelation phenomena at 34 ± 1°C and non-Newtonian behavior. Storage of optimized formulation carboxymethyl cellulose and optimal formulation chitosan at 4°C for 8 weeks resulted in repeatable release profiles at 34°C when sampled, with a Fickian mechanism earlier on but moving toward anomalous transport by week 8. Polymers (Pluronic® 127, carboxymethyl cellulose, and chitosan) demonstrated no significant cellular toxicity to human nasal epithelial cells up to 4 mg/mL and up to 1 mM for AMT (IC50: 4.5 ± 0.05 mM). Optimized formulation carboxymethyl cellulose and optimal formulation chitosan demonstrated slower release across an in vitro human nasal airway model (43%-44% vs 79 ± 4.58% for AMT). Using a human nasal cast model, deposition into the olfactory regions (potential nose-to-brain) was demonstrated on nozzle insertion (5 mm), whereas tilting of the head forward (15°) resulted in greater deposition in the bulk of the nasal cavity.

Alginates as therapeutic and drug delivery systems

Alginate is widely used as a viscosity enhancer in many different pharmaceutical formulations. The aim of this thesis is to quantitatively describe the functions of this polyelectrolyte in pharmaceutical systems. To do this the techniques used were Viscometry, Light Scattering, Continuous and Oscillatory Shear Rheometry, Numerical Analysis and Diffusion. Molecular characterization of the Alginate was carried out using Viscometry and Light Scattering to determine the molecular weight, the radius of gyration, the second virial coefficient and the Kuhn statistical segment length. The results showed good agreement with similar parameters obtained in previous studies. By blending Alginate with other polyelectrolytes, Xanthan Gum and 'Carbopol', in various proportions and with various methods of low and high shear preparation, a very wide range of dynamic rheological properties was found. Using oscillatory testing, the parameters often varied over several decades of magnitude. It was shown that the determination of the viscous and elastic components is particularly useful in describing the rheological 'profiles' of suspending agent blends and provides a step towards the non-empirical formulation of pharmaceutical disperse systems. Using numerical analysis of equations describing planar diffusion, it was shown that the analysis of drug release profiles alone does not provide unambiguous information about the mechanism of rate control. These principles were applied to the diffusion of Ibuprofen in Calcium Alginate gels. For diffusion in such non-Newtonian systems, emphasis was placed on the use of the elastic as well as the viscous component of viscoelasticity. It was found that the diffusion coefficients were relatively unaffected by increases in polymer concentration up to 5 per cent, yet the elasticities measured by oscillatory shear rheometry were increased. This was interpreted in the light of several theories of diffusion in gels.

Non-aqueous silicone elastomer gels as a vaginal microbicide delivery system for the HIV-1 entry inhibitor maraviroc

Aqueous semi-solid polymeric gels, such as those based on hydroxyethylcellulose (HEC) and polyacrylic acid (e.g. Carbopol®), have a long history of use in vaginal drug delivery. However, despite their ubiquity, they often provide sub-optimal clinical performance, due to poor mucosal retention and limited solubility for poorly water-soluble actives. These issues are particularly pertinent for vaginal HIV microbicides, since many lead candidates are poorly water-soluble and where a major goal is the development of a coitally independent, once daily gel product. In this study, we report the use of a non-aqueous silicone elastomer gel for vaginal delivery of the HIV-1 entry inhibitor maraviroc. In vitro rheological, syringeability and retention studies demonstrated enhanced performance for silicone gels compared with a conventional aqueous HEC gel, while testing of the gels in the slug model confirmed a lack of mucosal irritancy. Pharmacokinetic studies following single dose vaginal administration of a maraviroc silicone gel in rhesus macaques showed higher and sustained MVC levels in vaginal fluid, vaginal tissue and plasma compared with a HEC gel containing the same maraviroc loading. The results demonstrate that non-aqueous silicone gels have potential as a formulation platform for coitally independent vaginal HIV microbicides.

Transition in inclined internally heated fluid layers

Non-linear solutions and studies of their stability are presented for flows in a homogeneously heated fluid layer under the influence of a constant pressure gradient or when the mass flux across any lateral cross-section of the channel is required to vanish. The critical Grashof number is determined by a linear stability analysis of the basic state which depends only on the z-coordinate perpendicular to the boundary. Bifurcating longitudinal rolls as well as secondary solutions depending on the streamwise x-coordinate are investigated and their amplitudes are determined as functions of the supercritical Grashof number for various Prandtl numbers and angles of inclination of the layer. Solutions that emerge from a Hopf bifurcation assume the form of propagating waves and can thus be considered as steady flows relative to an appropriately moving frame of reference. The stability of these solutions with respect to three-dimensional disturbances is also analyzed in order to identify possible bifurcation points for evolving tertiary flows.

Range of validity of weakly non-linear theory in Rayleigh-Bénard problem

In this paper we examine the equilibrium states of finite amplitude flow in a horizontal fluid layer with differential heating between the two rigid boundaries. The solutions to the Navier-Stokes equations are obtained by means of a perturbation method for evaluating the Landau constants and through a Newton-Raphson iterative method that results from the Fourier expansion of the solutions that bifurcate above the linear stability threshold of infinitesimal disturbances. The results obtained from these two different methods of evaluating the convective flow are compared in the neighborhood of the critical Rayleigh number. We find that for small Prandtl numbers the discrepancy of the two methods is noticeable. © 2009 The Physical Society of Japan.

Design and characterization of zeolite membranes for fluid separation

Experimental and theoretical methods have been used to study zeolite structures, properties and applications as membranes for separation purposes. Thin layers of silicalite-1 and Na-LTA zeolites have been synthesised onto carbon-graphite supports using a hydrothermal synthesis procedure. The separation behaviour of the composite membranes was characterized by gas permeation studies of pure, binary and ternary mixtures of methane, ethane and propane. The influence of temperature and feed gas mixture composition on the separation and selectivity performance of the membranes was also investigated. It was found that the silicalite-1 composite membranes synthesised onto the 4 hour oxidized carbon-graphite supports showed the most promising separation behaviour of all the composite membranes investigated. Molecular simulation methods were used to gain an understanding of how hydrocarbon molecules behave both within the pores and on the surfaces of silicalite-1, mordenite and LTA zeolites. Molecular dynamic simulations were used to investigate the influence of temperature and molecular loadings on the diffusional behaviour of hydrocarbons in zeolites. Both hydroxylated (surface termination with hydroxyl groups) and non-hydroxylated silicalite-1 and Na-mordenite surfaces were generated. For both zeolites the most stable surfaces correspond to the {010} surface. For the silicalite-1 {010} surface the adsorption of hydrocarbons and molecular water onto the hydroxylated surface showed a favourable exothermic adsorption process compared to adsorption on the non-hydroxylated surface. With the Na-mordenite {010} surface the adsorption of hydrocarbons onto both the hydroxylated and non-hydroxylated surfaces had a combination of favourable and non-favourable adsorption energies, while the adsorption of molecular water onto both types of surface was found to be a favourable adsorption process.

Phytoestrogens modulate breast cancer resistance protein expression and function at the blood-cerebrospinal fluid barrier

PURPOSE: Breast cancer resistance protein (BCRP/ABCG2) is a drug efflux transporter expressed at the blood cerebrospinal fluid barrier (BCSFB), and influences distribution of drugs into the central nervous systems (CNS). Current inhibitors have failed clinically due to neurotoxicity. Novel approaches are needed to identify new modulators to enhance CNS delivery. This study examines 18 compounds (mainly phytoestrogens) as modulators of the expression/function of BCRP in an in vitro rat choroid plexus BCSFB model. METHODS: Modulators were initially subject to cytotoxicity (MTT) assessment to determine optimal non-toxic concentrations. Reverse-transcriptase PCR and confocal microscopy were used to identify the presence of BCRP in Z310 cells. Thereafter modulation of the intracellular accumulation of the fluorescent BCRP probe substrate Hoechst 33342 (H33342), changes in protein expression of BCRP (western blotting) and the functional activity of BCRP (membrane insert model) were assessed under modulator exposure. RESULTS: A 24 hour cytotoxicity assay (0.001 µM-1000 µM) demonstrated the majority of modulators possessed a cellular viability IC50 > 148 µM. Intracellular accumulation of H33342 was significantly increased in the presence of the known BCRP inhibitor Ko143 and, following a 24 hour pre-incubation, all modulators demonstrated statistically significant increases in H33342 accumulation (P < 0.001), when compared to control and Ko143. After a 24 hour pre-incubation with modulators alone, a 0.16-2.5-fold change in BCRP expression was observed for test compounds. The functional consequences of this were confirmed in a permeable insert model of the BCSFB which demonstrated that 17-β-estradiol, naringin and silymarin (down-regulators) and baicalin (up-regulator) can modulate BCRP-mediated transport function at the BCSFB. CONCLUSION: We have successfully confirmed the gene and protein expression of BCRP in Z310 cells and demonstrated the potential for phytoestrogen modulators to influence the functionality of BCRP at the BCSFB and thereby potentially allowing manipulation of CNS drug disposition.

Comprehensive secondary pyrolysis in fluidized-bed fast pyrolysis of biomass, a fluid dynamics based modelling effort

Homogenous secondary pyrolysis is category of reactions following the primary pyrolysis and presumed important for fast pyrolysis. For the comprehensive chemistry and fluid dynamics, a probability density functional (PDF) approach is used; with a kinetic scheme comprising 134 species and 4169 reactions being implemented. With aid of acceleration techniques, most importantly Dimension Reduction, Chemistry Agglomeration and In-situ Tabulation (ISAT), a solution within reasonable time was obtained. More work is required; however, a solution for levoglucosan (C6H10O5) being fed through the inlet with fluidizing gas at 500 °C, has been obtained. 88.6% of the levoglucosan remained non-decomposed, and 19 different decomposition product species were found above 0.01% by weight. A homogenous secondary pyrolysis scheme proposed can thus be implemented in a CFD environment and acceleration techniques can speed-up the calculation for application in engineering settings.