A simple lifestyle score predicts survival in healthy elderly men

Background. Although improvements in life expectancy have been attributed in part to the adoption or a more prudent lifestyle, few studies have examined the association of lifestyle with survival, using several lifestyle factors simultaneously, in a healthy elderly population. Methods. We investigated the association of health related behaviors with mortality in 7989 men aged 65 to 83 years participating in a population-based trial in Perth, Western Australia, by calculating a lifestyle score as a simple tally of flow many or eight prudent behaviors each individual followed. Results. Invitations to screening produced a corrected response of 70.5%. Out of a possible score of 8.46% of men had a score of less than 5. Within 5 years, a total of 703 men (9%) had died from any cause. The hazard ratio in men with a low lifestyle score was 1.3 [95% confidence interval (CI): 1.1-1.5] compared with men with a score of 5 or more. Conclusions. Lifestyle remains an important predictor of mortality even in old age. Survival in older men without a history of cardiovascular disease can potentially be enhanced by promoting a healthy lifestyle. © 2004 Elsevier Inc, All rights reserved.

Atmospheric boundary layer development over a narrow coastal plain during onshore flow

A study of the structure of the daytime atmospheric boundary layer during onshore flow over a narrow coastal plain is presented. The main emphasis of the study is on the nature and causes of heating and cooling observed in the boundary layer temperature profiles. Measurements included vertical temperature profiles above at least two sites derived from radiosondes and aircraft, as well as surface estimates of radiative and sensible heat fluxes. Surface meteorological and pilot balloon data were also available, providing further evidence of short-term changes in atmospheric boundary layer structure. The Manawatu case was representative of autumnal anticyclonic conditions with weak pressure gradients, and illustrated typical diurnal development of a convective boundary layer over a coastal plain bordered by mountain ranges, with a transition from a stable nocturnal situation to a well-mixed profile in the afternoon. The profiles show surface input of heat propagating upwards through the boundary layer during the day, as well as entrainment of heat at the top associated with shear induced turbulence and/or penetrative convection. Applying a one-dimensional model, estimates of boundary layer heat budget components were obtained for four time periods during the day. Later periods were affected by cumulus cloud development at the top of the boundary layer, resulting in significant changes in individual components. Input of sensible heat from the surface decreased, while the addition of heat to the boundary layer from both cloud condensation and advection increased.

Pendant drop formation of shear-thinning and yield stress fluids

A volume-of-fluid numerical method is used to predict the dynamics of shear-thinning liquid drop formation in air from a circular orifice. The validity of the numerical calculation is confirmed for a Newtonian liquid by comparison with experimental measurements. For particular values of Weber number and Froude number, predictions show a more rapid pinch-off, and a reduced number of secondary droplets, with increasing shear-thinning. Also a minimum in the limiting drop length occurs for the smallest Weber number as the zero-shear viscosity is varied. At the highest viscosity, the drop length is reduced due to shear-thinning, whereas at lower viscosities there is little effect of shear-thinning. The evolution of predicted drop shape, drop thickness and length, and the configuration at pinch-off are discussed for shear-thinning drops. The evolution of a drop of Bingham yield stress liquid is also considered as a limiting case. In contrast to the shear-thinning cases, it exhibits a plug flow prior to necking, an almost step-change approach to pinch-off of a torpedo shaped drop following the onset of necking, and a much smaller neck length; no secondary drops are formed. The results demonstrate the potential of the numerical model as a design tool in tailoring the fluid rheology for controlling drop formation behaviour. (c) 2006 Elsevier Inc. All rights reserved.

Application of a coupled ground-surface water flow model to simulate periodic groundwater flow influenced by a sloping boundary, capillarity and vertical flows

Theoretical developments as well as field and laboratory data have shown the influence of the capillary fringe on water table fluctuations to increase with the fluctuation frequency. The numerical solution of a full, partially saturated flow equation can be computationally expensive. In this paper, the influence of the capillary fringe on water table fluctuations is simplified through its parameterisation into the storage coefficient of a fully-saturated groundwater flow model using the complex effective porosity concept [Nielsen, P., Perrochet, P., 2000. Water table dynamics under capillary fringes: experiments and modelling. Advances in Water Resources 23 (1), 503-515; Nielsen, P., Perrochet, P., 2000. ERRATA: water table dynamics under capillary fringes: experiments and modelling (Advances in Water Resources 23 (2000) 503-515). Advances in Water Resources 23, 907-908]. The model is applied to sand flume observations of periodic water table fluctuations induced by simple harmonic forcing across a sloping boundary, analogous to many beach groundwater systems. While not providing information on the moisture distribution within the aquifer, this approach can reasonably predict the water table fluctuations in response to periodic forcing across a sloping boundary. Furthermore, he coupled ground-surface water model accurately predicts the extent of the seepage face formed at the sloping boundary. (C) 2005 Elsevier Ltd. All rights reserved.

Transport of simple fluids in nanopores: Theory and simulation

A theory is discussed of single-component transport in nanopores, recently developed by Bhatia and coworkers. The theory considers the oscillatory motion of molecules between diffuse wall collisions, arising from the fluid-wall interaction, along with superimposed viscous flow due to fluid-fluid interaction. The theory is tested against molecular dynamics simulations for hydrogen, methane, and carbon tetrafluoride flow in cylindrical nanopores in silica. Although exact at low densities, the theory performs well even at high densities, with the density dependency of the transport coefficient arising from viscous effects. Such viscous effects are reduced at high densities because of the large increase in viscosity, which explains the maximum in the transport coefficient with increase in density. Further, it is seen that in narrow pore sizes of less than two molecular diameters, where a complete monolayer cannot form on the surface, the mutual interference of molecules on opposite sides of the cross section can reduce the transport coefficient, and lead to a maximum in the transport coefficient with increasing density. The theory is also tested for the case of partially diffuse reflection and shows the viscous contribution to be negligible when the reflection is nearly specular. (c) 2005 American Institute of Chemical Engineers AIChE J, 52: 29-38, 2006.

Measuring and predicting the dynamics of linear monodisperse entangled polymers in rapid flow through an abrupt contraction. A small angle neutron scattering study

Small-angle neutron scattering measurements on a series of monodisperse linear entangled polystyrene melts in nonlinear flow through an abrupt 4:1 contraction have been made. Clear signatures of melt deformation and subsequent relaxation can be observed in the scattering patterns, which were taken along the centerline. These data are compared with the predictions of a recently derived molecular theory. Two levels of molecular theory are used: a detailed equation describing the evolution of molecular structure over all length scales relevant to the scattering data and a simplified version of the model, which is suitable for finite element computations. The velocity field for the complex melt flow is computed using the simplified model and scattering predictions are made by feeding these flow histories into the detailed model. The modeling quantitatively captures the full scattering intensity patterns over a broad range of data with independent variation of position within the contraction geometry, bulk flow rate and melt molecular weight. The study provides a strong, quantitative validation of current theoretical ideas concerning the microscopic dynamics of entangled polymers which builds upon existing comparisons with nonlinear mechanical stress data. Furthermore, we are able to confirm the appreciable length scale dependence of relaxation in polymer melts and highlight some wider implications of this phenomenon.

Differential diffusion: Often a finite-mising length effect

Many instances of differential diffusion, i e, different species having different turbulent diffusion coefficients in the same flow, can be explained as a finite mixing length effect. That is, in a simple mixing length scenario, the turbulent diffusion coefficient has the form 1 ( m )2 m m c l K w l OL   =  +    where, wm is the mixing velocity, lm the mixing length and Lc the overall distribution scale for a particular species. The first term represents the familiar gradient diffusion while the second term becomes important when lm/Lc is finite. This second term shows that different species will have different diffusion coefficients if they have different overall distribution scales. Such different Lcs may come about due to different boundary conditions and different intrinsic properties (molecular diffusivity, settling velocity etc) for different species. For momentum transfer in turbulent oscillatory boundary layers the second term is imaginary and explains observed phase leads of shear stresses ahead of velocity gradients.

Direct numerical simulation of turbulent flow and mass transfer around a rotating circular cylinder

Turbulent flow around a rotating circular cylinder has numerous applications including wall shear stress and mass-transfer measurement related to the corrosion studies. It is also of interest in the context of flow over convex surfaces where standard turbulence models perform poorly. The main purpose of this paper is to elucidate the basic turbulence mechanism around a rotating cylinder at low Reynolds numbers to provide a better understanding of flow fundamentals. Direct numerical simulation (DNS) has been performed in a reference frame rotating at constant angular velocity with the cylinder. The governing equations are discretized by using a finite-volume method. As for fully developed channel, pipe, and boundary layer flows, a laminar sublayer, buffer layer, and logarithmic outer region were observed. The level of mean velocity is lower in the buffer and outer regions but the logarithmic region still has a slope equal to the inverse of the von Karman constant. Instantaneous flow visualization revealed that the turbulence length scale typically decreases as the Reynolds number increases. Wavelet analysis provided some insight into the dependence of structural characteristics on wave number. The budget of the turbulent kinetic energy was computed and found to be similar to that in plane channel flow as well as in pipe and zero pressure gradient boundary layer flows. Coriolis effects show as an equivalent production for the azimuthal and radial velocity fluctuations leading to their ratio being lowered relative to similar nonrotating boundary layer flows.

Variations of local heat transfer coefficient in piped flow of viscous liquids

Measurements were carried out to determine local coefficients of heat transfer in short lengths of horizontal pipe, and in the region of an discontinuity in pipe diameter. Laminar, transitional and turbulent flow regimes were investigated, and mixtures of propylene glycol and water were used in the experiments to give a range of viscous fluids. Theoretical and empirical analyses were implemented to find how the fundamental mechanism of forced convection was modified by the secondary effects of free convection, temperature dependent viscosity, and viscous dissipation. From experiments with the short tube it was possible to determine simple empirical relationships describing the axial distribution of the local 1usselt number and its dependence on the Reynolds and Prandtl numbers. Small corrections were made to account for the secondary effects mentioned above. Two different entrance configurations were investigated to demonstrate how conditions upstream could influence the heat transfer coefficients measured downstream In experiments with a sudden contraction in pipe diameter the distribution of local 1u3se1t number depended on the Prandtl number of the fluid in a complicated way. Graphical data is presented describing this dependence for a range of fluids indicating how the local Nusselt number varied with the diameter-ratio. Ratios up to 3.34:1 were considered. With a sudden divergence in pipe diameter, it was possible to derive the axial distribution of the local Nusse1t number for a range of Reynolds and Prandtl numbers in a similar way to the convergence experiments. Difficulty was encountered in explaining some of the measurements obtained at low Reynolds numbers, and flow visualization techniques wore used to determine the complex flow patterns which could lead to the anomalous results mentioned. Tests were carried out with divergences up to 1:3.34 to find the way in which the local Nusselt number varied with the diameter ratio, and a few experiments were carried out with very large ratios up .to 14.4. A limited amount of theoretical analysis of the 'divergence' system was carried out to substantiate certain explanations of the heat transfer mechanisms postulated.

Transition in homogeneously heated inclined plane parallel shear flows

The transition of internally heated inclined plane parallel shear flows is examined numerically for the case of finite values of the Prandtl number Pr. We show that as the strength of the homogeneously distributed heat source is increased the basic flow loses stability to two-dimensional perturbations of the transverse roll type in a Hopf bifurcation for the vertical orientation of the fluid layer, whereas perturbations of the longitudinal roll type are most dangerous for a wide range of the value of the angle of inclination. In the case of the horizontal inclination transverse roll and longitudinal roll perturbations share the responsibility for the prime instability. Following the linear stability analysis for the general inclination of the fluid layer our attention is focused on a numerical study of the finite amplitude secondary travelling-wave solutions (TW) that develop from the perturbations of the transverse roll type for the vertical inclination of the fluid layer. The stability of the secondary TW against three-dimensional perturbations is also examined and our study shows that for Pr=0.71 the secondary instability sets in as a quasi-periodic mode, while for Pr=7 it is phase-locked to the secondary TW. The present study complements and extends the recent study by Nagata and Generalis (2002) in the case of vertical inclination for Pr=0.

Modelling and measurements on a condensing counter-flow heat exchanger

Experimental investigations and computer modelling studies have been made on the refrigerant-water counterflow condenser section of a small air to water heat pump. The main object of the investigation was a comparative study between the computer modelling predictions and the experimental observations for a range of operating conditions but other characteristics of a counterflow heat exchanger are also discussed. The counterflow condenser consisted of 15 metres of a thermally coupled pair of copper pipes, one containing the R12 working fluid and the other water flowing in the opposite direction. This condenser was mounted horizontally and folded into 0.5 metre straight sections. Thermocouples were inserted in both pipes at one metre intervals and transducers for pressure and flow measurement were also included. Data acquisition, storage and analysis was carried out by a micro-computer suitably interfaced with the transducers and thermocouples. Many sets of readings were taken under a variety of conditions, with air temperature ranging from 18 to 26 degrees Celsius, water inlet from 13.5 to 21.7 degrees, R12 inlet temperature from 61.2 to 81.7 degrees and water mass flow rate from 6.7 to 32.9 grammes per second. A Fortran computer model of the condenser (originally prepared by Carrington[1]) has been modified to match the information available from experimental work. This program uses iterative segmental integration over the desuperheating, mixed phase and subcooled regions for the R12 working fluid, the water always being in the liquid phase. Methods of estimating the inlet and exit fluid conditions from the available experimental data have been developed for application to the model. Temperature profiles and other parameters have been predicted and compared with experimental values for the condenser for a range of evaporator conditions and have shown that the model gives a satisfactory prediction of the physical behaviour of a simple counterflow heat exchanger in both single phase and two phase regions.

Hairpin vortex solution in planar Couette flow: a tapestry of knotted vortices

A numerical continuation method has been carried out seeking solutions for two distinct flow configurations, planar Couette flow (PCF) and laterally heated flow in a vertical slot (LHF). We found that the spanwise vortex solution in LHF identifies a new solution in PCF. The vortical structure of our new solution has the shape of a hairpin observed ubiquitously in high-Reynolds-number turbulent flow, and we believe this discovery may provide the paradigm for a hierarchical organization of coherent structures in turbulent shear layers.

Symmetry of coherent vortices in plane couette flow

A numerical continuation method is carried out in a homotopy space connecting two different flows, the Plane Couette Flow (PCF) and the Laterally Heated Flow in a vertical slot (LHF). This numerical continuation method enables us to obtain an exact steady solution in PCF. The new solution has the shape of hairpin vortices (HVS: hairpin vortex solution), which is observed ubiquitously in turbulent shear flows.

Flow field analysis of some mixing and conveying screw element regions, within a closely intermeshing, co-rotating twin-screw extruder.

Grafting of antioxidants and other modifiers onto polymers by reactive extrusion, has been performed successfully by the Polymer Processing and Performance Group at Aston University. Traditionally the optimum conditions for the grafting process have been established within a Brabender internal mixer. Transfer of this batch process to a continuous processor, such as an extruder, has, typically, been empirical. To have more confidence in the success of direct transfer of the process requires knowledge of, and comparison between, residence times, mixing intensities, shear rates and flow regimes in the internal mixer and in the continuous processor.The continuous processor chosen for the current work in the closely intermeshing, co-rotating twin-screw extruder (CICo-TSE). CICo-TSEs contain screw elements that convey material with a self-wiping action and are widely used for polymer compounding and blending. Of the different mixing modules contained within the CICo-TSE, the trilobal elements, which impose intensive mixing, and the mixing discs, which impose extensive mixing, are of importance when establishing the intensity of mixing. In this thesis, the flow patterns within the various regions of the single-flighted conveying screw elements and within both the trilobal element and mixing disc zones of a Betol BTS40 CICo-TSE, have been modelled using the computational fluid dynamics package Polyflow. A major obstacle encountered when solving the flow problem within all of these sets of elements, arises from both the complex geometry and the time-dependent flow boundaries as the elements rotate about their fixed axes. Simulation of the time dependent boundaries was overcome by selecting a number of sequential 2D and 3D geometries, used to represent partial mixing cycles. The flow fields were simulated using the ideal rheological properties of polypropylene and characterised in terms of velocity vectors, shear stresses generated and a parameter known as the mixing efficiency. The majority of the large 3D simulations were performed on the Cray J90 supercomputer situated at the Rutherford-Appleton laboratories, with pre- and postprocessing operations achieved via a Silicon Graphics Indy workstation. A mechanical model was constructed consisting of various CICo-TSE elements rotating within a transparent outer barrel. A technique has been developed using coloured viscous clays whereby the flow patterns and mixing characteristics within the CICo-TSE may be visualised. In order to test and verify the simulated predictions, the patterns observed within the mechanical model were compared with the flow patterns predicted by the computational model. The flow patterns within the single-flighted conveying screw elements in particular, showed good agreement between the experimental and simulated results.

Purification of plasmid DNA for use in human gene therapy

Two key issues defined the focus of this research in manufacturing plasmid DNA for use In human gene therapy. First, the processing of E.coli bacterial cells to effect the separation of therapeutic plasmid DNA from cellular debris and adventitious material. Second, the affinity purification of the plasmid DNA in a Simple one-stage process. The need arises when considering the concerns that have been recently voiced by the FDA concerning the scalability and reproducibility of the current manufacturing processes in meeting the quality criteria of purity, potency, efficacy, and safety for a recombinant drug substance for use in humans. To develop a preliminary purification procedure, an EFD cross-flow micro-filtration module was assessed for its ability to effect the 20-fold concentration, 6-time diafiltration, and final clarification of the plasmid DNA from the subsequent cell lysate that is derived from a 1 liter E.coli bacterial cell culture. Historically, the employment of cross-flow filtration modules within procedures for harvesting cells from bacterial cultures have failed to reach the required standards dictated by existing continuous centrifuge technologies, frequently resulting in the rapid blinding of the membrane with bacterial cells that substantially reduces the permeate flux. By challenging the EFD module, containing six helical wound tubular membranes promoting centrifugal instabilities known as Dean vortices, with distilled water between the Dean number's of 187Dn and 818Dn,and the transmembrane pressures (TMP) of 0 to 5 psi. The data demonstrated that the fluid dynamics significantly influenced the permeation rate, displaying a maximum at 227Dn (312 Imh) and minimum at 818Dn (130 Imh) for a transmembrane pressure of 1 psi. Numerical studies indicated that the initial increase and subsequent decrease resulted from a competition between the centrifugal and viscous forces that create the Dean vortices. At Dean numbers between 187Dn and 227Dn , the forces combine constructively to increase the apparent strength and influence of the Dean vortices. However, as the Dean number in increases above 227 On the centrifugal force dominates the viscous forces, compressing the Dean vortices into the membrane walls and reducing their influence on the radial transmembrane pressure i.e. the permeate flux reduced. When investigating the action of the Dean vortices in controlling tile fouling rate of E.coli bacterial cells, it was demonstrated that the optimum cross-flow rate at which to effect the concentration of a bacterial cell culture was 579Dn and 3 psi TMP, processing in excess of 400 Imh for 20 minutes (i.e., concentrating a 1L culture to 50 ml in 10 minutes at an average of 450 Imh). The data demonstrated that there was a conflict between the Dean number at which the shear rate could control the cell fouling, and the Dean number at which tile optimum flux enhancement was found. Hence, the internal geometry of the EFD module was shown to sub-optimal for this application. At 579Dn and 3 psi TMP, the 6-fold diafiltration was shown to occupy 3.6 minutes of process time, processing at an average flux of 400 Imh. Again, at 579Dn and 3 psi TMP the clarification of the plasmid from tile resulting freeze-thaw cell lysate was achieved at 120 Iml1, passing 83% (2,5 mg) of the plasmid DNA (6,3 ng μ-1 10.8 mg of genomic DNA (∼23,00 Obp, 36 ng μ-1 ), and 7.2 mg of cellular proteins (5-100 kDa, 21.4 ngμ-1 ) into the post-EFD process stream. Hence the EFD module was shown to be effective, achieving the desired objectives in approximately 25 minutes. On the basis of its ability to intercalate into low molecular weight dsDNA present in dilute cell lysates, and be electrophoresed through agarose, the fluorophore PicoGreen was selected for the development of a suitable dsDNA assay. It was assesseel for its accuracy, and reliability, In determining the concentration and identity of DNA present in samples that were eleclrophoresed through agarose gels. The signal emitted by intercalated PicoGreen was shown to be constant and linear, and that the mobility of the PicaGreen-DNA complex was not affected by the intercalation. Concerning the secondary purification procedure, various anion-exchange membranes were assessed for their ability to capture plasmid DNA from the post-EFD process stream. For a commercially available Sartorius Sartobind Q15 membrane, the reduction in the equilibriumbinding capacity for  ctDNA in buffer of increasing ionic demonstrated that DNA was being.adsorbed by electrostatic  interactions only. However, the problems associated with fluid distribution across the membrane demonstrated that the membrane housing was the predominant cause of the .erratic breakthrough curves. Consequently, this would need to be rectified before such a membrane could be integrated into the current system, or indeed be scaled beyond laboratory scale. However, when challenged with the process material, the data showed that considerable quantities of protein (1150 μg) were adsorbed preferentially to the plasmid DNA (44 μg). This was also shown for derived Pall Gelman UltraBind US450 membranes that had been functionalised by varying molecular weight poly-L~lysine and polyethyleneimine ligands. Hence the anion-exchange membranes were shown to be ineffective in capturing plasmid DNA from the process stream. Finally, work was performed to integrate a sequence-specific DNA·binding protein into a single-stage DNA chromatography, isolating plasmid DNA from E.coli cells whilst minimising the contamination from genomic DNA and cellular protein. Preliminary work demonstrated that the fusion protein was capable of isolating pUC19 DNA into which the recognition sequence for the fusion-protein had been inserted (pTS DNA) when in the presence of the conditioned process material. Althougth the pTS recognition sequence differs from native pUC19 sequences by only 2 bp, the fusion protein was shown to act as a highly selective affinity ligand for pTS DNA alone. Subsequently, the scale of the process was scaled 25-fold and positioned directly following the EFD system. In conclusion, the integration of the EFD micro-filtration system and zinc-finger affinity purification technique resulted in the capture of approximately 1 mg of plasmid DNA was purified from 1L of E.coli  culture in a simple two stage process, resulting in the complete removal of genomic DNA and 96.7% of cellular protein in less than 1 hour of process time.