Osmotic flow in porous membranes:effects of electric charge

The electrostatic model for osmotic flow across a porous membrane in our previous study (Akinaga et al. 2008)" was extended to include the streaming potential, for solutes and pores of like charge and fixed surface charge densities. The magnitude of the streaming potential was determined to satisfy zero current condition along the pore axis. It was found that the streaming potential affects the velocity profiles of the pressure driven flow as well as the osmotic flow through the pore, and decreases their flow rates, particularly in the case of large Debye length relative to the pore radius, whereas it has little effect on the reflection coefficients of spherical solutes through cylindrical pores.

Modelling and simulation of biomass gasification in a circulating fluidized bed reactor

Computational Fluid Dynamics (CFD) has found great acceptance among the engineering community as a tool for research and design of processes that are practically difficult or expensive to study experimentally. One of these processes is the biomass gasification in a Circulating Fluidized Bed (CFB). Biomass gasification is the thermo-chemical conversion of biomass at a high temperature and a controlled oxygen amount into fuel gas, also sometime referred to as syngas. Circulating fluidized bed is a type of reactor in which it is possible to maintain a stable and continuous circulation of solids in a gas-solid system. The main objectives of this thesis are four folds: (i) Develop a three-dimensional predictive model of biomass gasification in a CFB riser using advanced Computational Fluid Dynamic (CFD) (ii) Experimentally validate the developed hydrodynamic model using conventional and advanced measuring techniques (iii) Study the complex hydrodynamics, heat transfer and reaction kinetics through modelling and simulation (iv) Study the CFB gasifier performance through parametric analysis and identify the optimum operating condition to maximize the product gas quality. Two different and complimentary experimental techniques were used to validate the hydrodynamic model, namely pressure measurement and particle tracking. The pressure measurement is a very common and widely used technique in fluidized bed studies, while, particle tracking using PEPT, which was originally developed for medical imaging, is a relatively new technique in the engineering field. It is relatively expensive and only available at few research centres around the world. This study started with a simple poly-dispersed single solid phase then moved to binary solid phases. The single solid phase was used for primary validations and eliminating unnecessary options and steps in building the hydrodynamic model. Then the outcomes from the primary validations were applied to the secondary validations of the binary mixture to avoid time consuming computations. Studies on binary solid mixture hydrodynamics is rarely reported in the literature. In this study the binary solid mixture was modelled and validated using experimental data from the both techniques mentioned above. Good agreement was achieved with the both techniques. According to the general gasification steps the developed model has been separated into three main gasification stages; drying, devolatilization and tar cracking, and partial combustion and gasification. The drying was modelled as a mass transfer from the solid phase to the gas phase. The devolatilization and tar cracking model consist of two steps; the devolatilization of the biomass which is used as a single reaction to generate the biomass gases from the volatile materials and tar cracking. The latter is also modelled as one reaction to generate gases with fixed mass fractions. The first reaction was classified as a heterogeneous reaction while the second reaction was classified as homogenous reaction. The partial combustion and gasification model consisted of carbon combustion reactions and carbon and gas phase reactions. The partial combustion considered was for C, CO, H2 and CH4. The carbon gasification reactions used in this study is the Boudouard reaction with CO2, the reaction with H2O and Methanation (Methane forming reaction) reaction to generate methane. The other gas phase reactions considered in this study are the water gas shift reaction, which is modelled as a reversible reaction and the methane steam reforming reaction. The developed gasification model was validated using different experimental data from the literature and for a wide range of operating conditions. Good agreement was observed, thus confirming the capability of the model in predicting biomass gasification in a CFB to a great accuracy. The developed model has been successfully used to carry out sensitivity and parametric analysis. The sensitivity analysis included: study of the effect of inclusion of various combustion reaction; and the effect of radiation in the gasification reaction. The developed model was also used to carry out parametric analysis by changing the following gasifier operating conditions: fuel/air ratio; biomass flow rates; sand (heat carrier) temperatures; sand flow rates; sand and biomass particle sizes; gasifying agent (pure air or pure steam); pyrolysis models used; steam/biomass ratio. Finally, based on these parametric and sensitivity analysis a final model was recommended for the simulation of biomass gasification in a CFB riser.

Concentration polarization model of spiral-wound membrane modules with application to batch-mode RO desalination of brackish water

Batch-mode reverse osmosis (batch-RO) operation is considered a promising desalination method due to its low energy requirement compared to other RO system arrangements. To improve and predict batch-RO performance, studies on concentration polarization (CP) are carried out. The Kimura-Sourirajan mass-transfer model is applied and validated by experimentation with two different spiral-wound RO elements. Explicit analytical Sherwood correlations are derived based on experimental results. For batch-RO operation, a new genetic algorithm method is developed to estimate the Sherwood correlation parameters, taking into account the effects of variation in operating parameters. Analytical procedures are presented, then the mass transfer coefficient models are developed for different operation processes, i.e., batch-RO and continuous RO. The CP related energy loss in batch-RO operation is quantified based on the resulting relationship between feed flow rates and mass transfer coefficients. It is found that CP increases energy consumption in batch-RO by about 25% compared to the ideal case in which CP is absent. For continuous RO process, the derived Sherwood correlation predicted CP accurately. In addition, we determined the optimum feed flow rate of our batch-RO system.

A comparative study of life history and gene flow in darters (Pisces: Percidae)

Gene flow, or the exchange of genes between populations, is important because it determines the evolutionary trajectory of a species, including the relative influences of genetic drift and natural selection in the process of population differentiation. Gene flow differs among species because of variation in dispersal capability and abundances across taxa, and historical forces related to geological or lineage history. Both history and ecology influence gene flow in potentially complicated ways, and accounting for their effects remains an important problem in evolutionary biology. This research is a comparative study of gene flow and life-history in a monophyletic group of stream fishes, the darters. As a first step in disentangling historical and ecological effects, I reconstructed the phylogenetic relationships of the study species from nucleotide sequences in the mtDNA control region. I then used this phylogeny and regional glaciation history to infer historical effects on life-history evolution and gene flow in 15 species of darters. Gene flow was estimated indirectly, using information from 20 resolvable and polymorphic allozyme loci. When I accounted for historical effects, comparisons across taxa revealed that gene flow rates were closely associated with differences in clutch sizes and reproductive investment patterns. I hypothesized that differences in larval dispersal among taxa explained this relationship. Results from a field study of larval drift were consistent with this hypothesis. Finally, I asked whether there was an interaction between species' ecology and genetic differentiation across biogeographically distinct regions. Information from allozymes and mtDNA sequences revealed that life history plays an important role in the magnitude of species divergence across biogeographic boundaries. These results suggested an important association between life histories and rates of speciation following an allopatric isolation event. This research, along with other studies from the literature, further illustrates the enormous potential of North American freshwater fishes as a system for studying speciation processes. ^

Lattice Boltzmann method for flow and heat transfer in microgeometries

Recent technological developments have made it possible to design various microdevices where fluid flow and heat transfer are involved. For the proper design of such systems, the governing physics needs to be investigated. Due to the difficulty to study complex geometries in micro scales using experimental techniques, computational tools are developed to analyze and simulate flow and heat transfer in microgeometries. However, conventional numerical methods using the Navier-Stokes equations fail to predict some aspects of microflows such as nonlinear pressure distribution, increase mass flow rate, slip flow and temperature jump at the solid boundaries. This necessitates the development of new computational methods which depend on the kinetic theory that are both accurate and computationally efficient. In this study, lattice Boltzmann method (LBM) was used to investigate the flow and heat transfer in micro sized geometries. The LBM depends on the Boltzmann equation which is valid in the whole rarefaction regime that can be observed in micro flows. Results were obtained for isothermal channel flows at Knudsen numbers higher than 0.01 at different pressure ratios. LBM solutions for micro-Couette and micro-Poiseuille flow were found to be in good agreement with the analytical solutions valid in the slip flow regime (0.01 < Kn < 0.1) and direct simulation Monte Carlo solutions that are valid in the transition regime (0.1 < Kn < 10) for pressure distribution and velocity field. The isothermal LBM was further extended to simulate flows including heat transfer. The method was first validated for continuum channel flows with and without constrictions by comparing the thermal LBM results against accurate solutions obtained from analytical equations and finite element method. Finally, the capability of thermal LBM was improved by adding the effect of rarefaction and the method was used to analyze the behavior of gas flow in microchannels. The major finding of this research is that, the newly developed particle-based method described here can be used as an alternative numerical tool in order to study non-continuum effects observed in micro-electro-mechanical-systems (MEMS).

Lattice Boltzmann Method for Flow and Heat Transfer in Microgeometries

Recent technological developments have made it possible to design various microdevices where fluid flow and heat transfer are involved. For the proper design of such systems, the governing physics needs to be investigated. Due to the difficulty to study complex geometries in micro scales using experimental techniques, computational tools are developed to analyze and simulate flow and heat transfer in microgeometries. However, conventional numerical methods using the Navier-Stokes equations fail to predict some aspects of microflows such as nonlinear pressure distribution, increase mass flow rate, slip flow and temperature jump at the solid boundaries. This necessitates the development of new computational methods which depend on the kinetic theory that are both accurate and computationally efficient. In this study, lattice Boltzmann method (LBM) was used to investigate the flow and heat transfer in micro sized geometries. The LBM depends on the Boltzmann equation which is valid in the whole rarefaction regime that can be observed in micro flows. Results were obtained for isothermal channel flows at Knudsen numbers higher than 0.01 at different pressure ratios. LBM solutions for micro-Couette and micro-Poiseuille flow were found to be in good agreement with the analytical solutions valid in the slip flow regime (0.01 < Kn < 0.1) and direct simulation Monte Carlo solutions that are valid in the transition regime (0.1 < Kn < 10) for pressure distribution and velocity field. The isothermal LBM was further extended to simulate flows including heat transfer. The method was first validated for continuum channel flows with and without constrictions by comparing the thermal LBM results against accurate solutions obtained from analytical equations and finite element method. Finally, the capability of thermal LBM was improved by adding the effect of rarefaction and the method was used to analyze the behavior of gas flow in microchannels. The major finding of this research is that, the newly developed particle-based method described here can be used as an alternative numerical tool in order to study non-continuum effects observed in micro-electro-mechanical-systems (MEMS).

Radionuclide analyses, and sedimentation and accumulation rates of sediment core PS2319-1

There is increasing evidence indicating that syndepositional redistribution of sediment on the seafloor by bottom currents is common and can significantly affect sediment mass accumulation rates. Notwithstanding its common incidence, this process (generally referred to as sediment focusing) is often difficult to recognize. If redistribution is near synchronous to deposition, the stratigraphy of the sediment is not disturbed and sediment focusing can easily be overlooked. Ignoring it, however, can lead to serious misinterpretations of sedimentary fluxes, particularly when past changes in export flux from the overlying water are inferred. In many instances, this problem can be resolved, at least for sediments deposited during the late Quaternary, by normalizing to the flux of 230Th scavenged from seawater, which is nearly constant and equivalent to the known rate of production of 230Th from the decay of dissolved 234U. We review the principle, advantages and limitations of this method. Notwithstanding its limitations, it is clear that 230Th normalization does provide a means of achieving more accurate interpretations of sedimentary fluxes and eliminates the risk of serious misinterpretations of sediment mass accumulation rates.

Volume averaging of multiphase flows with hydrate formation in subsea pipelines

In oil and gas pipeline operations, the gas, oil, and water phases simultaneously move through pipe systems. The mixture cools as it flows through subsea pipelines, and forms a hydrate formation region, where the hydrate crystals start to grow and may eventually block the pipeline. The potential of pipe blockage due to hydrate formation is one of the most significant flow-assurance problems in deep-water subsea operations. Due to the catastrophic safety and economic implications of hydrate blockage, it is important to accurately predict the simultaneous flow of gas, water, and hydrate particles in flowlines. Currently, there are few or no studies that account for the simultaneous effects of hydrate growth and heat transfer on flow characteristics within pipelines. This thesis presents new and more accurate predictive models of multiphase flows in undersea pipelines to describe the simultaneous flow of gas, water, and hydrate particles through a pipeline. A growth rate model for the hydrate phase is presented and then used in the development of a new three-phase model. The conservation equations of mass, momentum, and energy are formulated to describe the physical phenomena of momentum and heat transfer between the fluid and the wall. The governing equations are solved based on an analytical-numerical approach using a Newton-Raphson method for the nonlinear equations. An algorithm was developed in Matlab software to solve the equations from the inlet to the outlet of the pipeline. The developed models are validated against a single-phase model with mixture properties, and the results of comparative studies show close agreement. The new model predicts the volume fraction and velocity of each phase, as well as the mixture pressure and temperature profiles along the length of the pipeline. The results from the hydrate growth model reveal the growth rate and location where the initial hydrates start to form. Finally, to assess the impact of certain parameters on the flow characteristics, parametric studies have been conducted. The results show the effect of a variation in the pipe diameter, mass flow rate, inlet pressure, and inlet temperature on the flow characteristics and hydrate growth rates.

Tooth wear in Irish teenagers: a laboratory and epidemiological study

Aim: To investigate the characteristics, development and determinants of toothwear among Irish schoolchildren. Methods: A cross-sectional (examination at 16-years-old) and longitudinal (examinations at 5-,12-,14-years) study were conducted. Two indices were used to measure toothwear, children/parents completed a demographic profile and questionnaire on oral hygiene and dietary practices, health, and lifestyle in both studies. Saliva was collected from consenting 16-year-olds. The explanatory variables for the cross-sectional and longitudinal study were derived from children/parents responses. Differences in salivary profiles were determined for subsets; the protein concentration was determined with Bradford protein assay and protein carbonyl concentration (a protein oxidation marker) was determined spectrophotometrically. Gel-electrophoresis and mass spectrometry determined proteins and ion chromatography inorganic ions. Statistical significance was accepted at p<0.05. Results: At 16-years-old the prevalence of toothwear with dentine visible was 44%. No difference in salivary flow rates existed. In unstimulated saliva a higher mean, protein carbonyl (p<0.0001) and total calcium concentration (p<0.002) existed for the group with moderate toothwear. In stimulated saliva the moderate toothwear group had a lower mean protein concentration(p<0.0001). The 2-DE protein spots prepared for a sub-group differed between those with toothwear and without. Mass spectrometry, identified one of the different proteins as IgA. For 16-year-olds, the self-reported factors indicated that brushing after breakfast was associated with lower toothwear scores(p<0.03). Nail-biting, being asthmatic or reporting a dry mouth were associated with higher toothwear scores(all p<0.05). Eating an apple daily or less was associated with less toothwear(p<0.002). In the longitudinal study toothwear into dentine at age five or 12-years was associated with more toothwear at age 14(all p<0.05). Discussion: The results illustrate the multifactorial aetiology of toothwear. The biochemical and physical correlates of saliva with toothwear requires further research. Conclusion: The impact of previous toothwear, salivary, dietary and personal factors on toothwear in the early permanent dentition is demonstrated.

Age determinaton and age model of sediments from the Bay of Biscay

Continuous high-resolution mass accumulation rates (MAR) and X-ray fluorescence (XRF) measurements from marine sediment records in the Bay of Biscay (NE Atlantic) have allowed the determination of the timing and the amplitude of the 'Fleuve Manche' (Channel River) discharges during glacial stages MIS 10, MIS 8, MIS 6 and MIS 4-2. These results have yielded detailed insight into the Middle and Late Pleistocene glaciations in Europe and the drainage network of the western and central European rivers over the last 350 kyr. This study provides clear evidence that the 'Fleuve Manche' connected the southern North Sea basin with the Bay of Biscay during each glacial period and reveals that 'Fleuve Manche' activity during the glaciations MIS 10 and MIS 8 was significantly less than during MIS 6 and MIS 2. We correlate the significant 'Fleuve Manche' activity, detected during MIS 6 and MIS 2, with the extensive Saalian (Drenthe Substage) and the Weichselian glaciations, respectively, confirming that the major Elsterian glaciation precedes the glacial MIS 10. In detail, massive 'Fleuve Manche' discharges occurred at ca 155 ka (mid-MIS 6) and during Termination I, while no significant discharges are found during Termination II. It is assumed that a substantial retreat of the European ice sheet at ca 155 kyr, followed by the formation of ice-free conditions between the British Isles and Scandinavia until Termination II, allowed meltwater to flow northwards through the North Sea basin during the second part of the MIS 6. We assume that this glacial pattern corresponds to the Warthe Substage glacial maximum, therefore indicating that the data presented here equates to the Drenthe and the Warthe glacial advances at ca 175-160 ka and ca 150-140 ka, respectively. Finally, the correlation of our records with ODP site 980 reveals that massive 'Fleuve Manche' discharges, related to partial or complete melting of the European ice masses, were synchronous with strong decreases in both the rate of deep-water formation and the strength of the Atlantic thermohaline circulation. 'Fleuve Manche' discharges over the last 350 kyr probably participated, with other meltwater sources, in the collapse of the thermohaline circulation by freshening the northern Atlantic surface water.

(Table 1) Radiolarian accumulation rates and geochemistry of ODP Site 184-1143

Late Miocene-Recent micropaleontological and geochemical records from Ocean Drilling Program (ODP) Site 1143 in the southern South China Sea (SCS) indicate that increase and decrease in abundance of siliceous plankton may be controlled mainly by the input of nutrients derived from land and provided by upwelling. A high export production event - a "biogenic bloom" event - occurred in the southern SCS between 12 and 6 Ma. During this period, high ratios of smectite/(illite + chlorite), smectite/quartz and Al/K indicate a high weathering intensity of the Asian continent, possibly due to the intensification of the East Asian Summer Monsoon (EASM), which may have increased the net flux of nutrients to the ocean, both directly through terrestrial input and indirectly through upwelling activity. A drop in Ba/Ti, Al/Ti and Ca/Ti values around 6 Ma may indicate a lowering of productivity, possibly due to the large consumption of sea surface nutrients by the "biogenic bloom". Alternatively, it may indicate a shift in terrigenous input source area. At about 5.4 Ma, a decrease in weathering intensity, as indicated by a sudden decrease in the values of smectite/(illite + chlorite), smectite/quartz and Al/K, might have led to a sudden decrease of terrestrial nutrient input to the SCS. We suggest that the biogenic bloom ended when nutrients in surface waters were exhausted, because of a decrease in supply as well as a decrease in upwelling intensity due to weakening of the EASM. As a result, radiolarians were absent in the studied area between ~6 and 3.2 Ma. At ~3.2 Ma, radiolarians began to recover, possibly because the start of Northern Hemispheric glaciation and the rapid uplift of the Tibet Plateau led to intensification of the East Asian monsoon. After the Mid-Pleistocene Climate Transition at 0.9 Ma, the abundance and mass accumulation rates of radiolarians increased, probably as a result of increased upwelling activity driven by the increasing intensity of the summer monsoon.

Computer simulations of soft particles in flow

Understanding the dynamics of blood cells is a crucial element to discover biological mechanisms, to develop new efficient drugs, design sophisticated microfluidic devices, for diagnostics. In this work, we focus on the dynamics of red blood cells in microvascular flow. Microvascular blood flow resistance has a strong impact on cardiovascular function and tissue perfusion. The flow resistance in microcirculation is governed by flow behavior of blood through a complex network of vessels, where the distribution of red blood cells across vessel cross-sections may be significantly distorted at vessel bifurcations and junctions. We investigate the development of blood flow and its resistance starting from a dispersed configuration of red blood cells in simulations for different hematocrits, flow rates, vessel diameters, and aggregation interactions between red blood cells. Initially dispersed red blood cells migrate toward the vessel center leading to the formation of a cell-free layer near the wall and to a decrease of the flow resistance. The development of cell-free layer appears to be nearly universal when scaled with a characteristic shear rate of the flow, which allows an estimation of the length of a vessel required for full flow development, $l_c \approx 25D$, with vessel diameter $D$. Thus, the potential effect of red blood cell dispersion at vessel bifurcations and junctions on the flow resistance may be significant in vessels which are shorter or comparable to the length $l_c$. The presence of aggregation interactions between red blood cells lead in general to a reduction of blood flow resistance. The development of the cell-free layer thickness looks similar for both cases with and without aggregation interactions. Although, attractive interactions result in a larger cell-free layer plateau values. However, because the aggregation forces are short-ranged at high enough shear rates ($\bar{\dot{\gamma}} \gtrsim 50~\text{s}^{-1}$) aggregation of red blood cells does not bring a significant change to the blood flow properties. Also, we develop a simple theoretical model which is able to describe the converged cell-free-layer thickness with respect to flow rate assuming steady-state flow. The model is based on the balance between a lift force on red blood cells due to cell-wall hydrodynamic interactions and shear-induced effective pressure due to cell-cell interactions in flow. We expect that these results can also be used to better understand the flow behavior of other suspensions of deformable particles such as vesicles, capsules, and cells. Finally, we investigate segregation phenomena in blood as a two-component suspension under Poiseuille flow, consisting of red blood cells and target cells. The spatial distribution of particles in blood flow is very important. For example, in case of nanoparticle drug delivery, the particles need to come closer to microvessel walls, in order to adhere and bring the drug to a target position within the microvasculature. Here we consider that segregation can be described as a competition between shear-induced diffusion and the lift force that pushes every soft particle in a flow away from the wall. In order to investigate the segregation, on one hand, we have 2D DPD simulations of red blood cells and target cell of different sizes, on the other hand the Fokker-Planck equation for steady state. For the equation we measure force profile, particle distribution and diffusion constant across the channel. We compare simulation results with those from the Fokker-Planck equation and find a very good correspondence between the two approaches. Moreover, we investigate the diffusion behavior of target particles for different hematocrit values and shear rates. Our simulation results indicate that diffusion constant increases with increasing hematocrit and depends linearly on shear rate. The third part of the study describes development of a simulation model of complex vascular geometries. The development of the model is important to reproduce vascular systems of small pieces of tissues which might be gotten from MRI or microscope images. The simulation model of the complex vascular systems might be divided into three parts: modeling the geometry, developing in- and outflow boundary conditions, and simulation domain decomposition for an efficient computation. We have found that for the in- and outflow boundary conditions it is better to use the SDPD fluid than DPD one because of the density fluctuations along the channel of the latter. During the flow in a straight channel, it is difficult to control the density of the DPD fluid. However, the SDPD fluid has not that shortcoming even in more complex channels with many branches and in- and outflows because the force acting on particles is calculated also depending on the local density of the fluid.

Salivary and serum cortisol levels, salivary alpha-amylase and unstimulated whole saliva flow rate in pregnant and non-pregnant

OBJETIVO: Comparar os níveis de cortisol sérico e salivar, alfa-amilase salivar (sAA) e fluxo de saliva não estimulada (UWS) em gestantes e não gestantes. MÉTODOS: Trata-se de um estudo longitudinal realizado no centro de promoção da saúde de um hospital universitário. Nove gestantes e 12 não gestantes participaram do estudo. Foram coletados e analisados soro e UWS nos três trimestres gestacionais e duas vezes por mês durante o ciclo menstrual. A análise do cortisol salivar e sérico foi realizada com o uso de quimiluminescência e a atividade da sAA foi determinada por meio de analisador automático para bioquímica. RESULTADOS: Foi verificado que a mediana (intervalo interquartil) dos níveis de cortisol sérico no grupo de gestantes foi maior que 23,8 µL/dL (19,4-29,4) quando comparado ao grupo de não gestantes, que teve média de 12,3 (9,6-16,8; p<0,001). Os níveis de sAA seguiram o mesmo padrão, com médias de 56,7 U/L (30,9-82,2) e 31,8 (18,1-53,2; p<0,001), respectivamente. Foram observadas diferenças dos níveis de cortisol sérico e salivar (µL/dL) e de sAA entre a fase folicular versus a fase lútea (p<0,001). As medianas dos fluxos salivares (UWS) foram semelhantes em gestantes (0,26 [0,15-0,30] mL/min) e não gestantes (0,23 [0,20-0,32] mL/min). Foram encontradas correlações significativas entre o cortisol salivar e o sérico (p=0,02) e entre o cortisol salivar e a sAA (p=0,01). CONCLUSÕES: Os níveis de cortisol sérico de sAA durante a gestação elevam-se. Na fase lútea do ciclo ovariano, os níveis de cortisol salivar aumentam ao passo que os níveis de cortisol sérico e sAA diminuem. _______________________________________________________________________________________ ABSTRACT

Periodic emission of droplets from an oscillating electrified meniscus of a low-viscosity, highly conductive liquid

The generation of identical droplets of controllable size in the micrometer range is a problem of much interest owing to the numerous technological applications of such droplets. This work reports an investigation of the regime of periodic emission of droplets from an electrified oscillating meniscus of a liquid of low viscosity and high electrical conductivity attached to the end of a capillary tube, which may be used to produce droplets more than ten times smaller than the diameter of the tube. To attain this periodic microdripping regime, termed axial spray mode II by Juraschek and Röllgen [R. Juraschek and F. W. Röllgen, Int. J. Mass Spectrom. 177, 1 (1998)], liquid is continuously supplied through the tube at a given constant flow rate, while a dc voltage is applied between the tube and a nearby counter electrode. The resulting electric field induces a stress at the surface of the liquid that stretches the meniscus until, in certain ranges of voltage and flow rate, it develops a ligament that eventually detaches, forming a single droplet, in a process that repeats itself periodically. While it is being stretched, the ligament develops a conical tip that emits ultrafine droplets, but the total mass emitted is practically contained in the main droplet. In the parametrical domain studied, we find that the process depends on two main dimensionless parameters, the flow rate nondimensionalized with the diameter of the tube and the capillary time, q, and the electric Bond number BE, which is a nondimensional measure of the square of the applied voltage. The meniscus oscillation frequency made nondimensional with the capillary time, f, is of order unity for very small flow rates and tends to decrease as the inverse of the square root of q for larger values of this parameter. The product of the meniscus mean volume times the oscillation frequency is nearly constant. The characteristic length and width of the liquid ligament immediately before its detachment approximately scale as powers of the flow rate and depend only weakly on the applied voltage. The diameter of the main droplets nondimensionalized with the diameter of the tube satisfies dd≈(6/π)1/3(q/f)1/3, from mass conservation, while the electric charge of these droplets is about 1/4 of the Rayleigh charge. At the minimum flow rate compatible with the periodic regimen, the dimensionless diameter of the droplets is smaller than one-tenth, which presents a way to use electrohydrodynamic atomization to generate droplets of highly conducting liquids in the micron-size range, in marked contrast with the cone-jet electrospray whose typical droplet size is in the nanometric regime for these liquids. In contrast with other microdripping regimes where the mass is emitted upon the periodic formation of a narrow capillary jet, the present regime gives one single droplet per oscillation, except for the almost massless fine aerosol emitted in the form of an electrospray.

WATER TRANSPORT THROUGH POROUS MEDIA AND IN FLOW CHANNELS OF PROTON EXCHANGE MEMBRANE FUEL CELL

Proton exchange membrane (PEM) fuel cell has been known as a promising power source for different applications such as automotive, residential and stationary. During the operation of a PEM fuel cell, hydrogen is oxidized in anode and oxygen is reduced in the cathode to produce the intended power. Water and heat are inevitable byproducts of these reactions. The water produced in the cathode should be properly removed from inside the cell. Otherwise, it may block the path of reactants passing through the gas channels and/or gas diffusion layer (GDL). This deteriorates the performance of the cell and eventually can cease the operation of the cell. Water transport in PEM fuel cell has been the subject of this PhD study. Water transport on the surface of the GDL, through the gas flow channels, and through GDL has been studied in details. For water transport on the surface of the GDL, droplet detachment has been measured for different GDL conditions and for anode and cathode gas flow channels. Water transport through gas flow channels has been investigated by measuring the two-phase flow pressure drop along the gas flow channels. As accumulated liquid water within gas flow channels resists the gas flow, the pressure drop increases along the flow channels. The two-phase flow pressure drop can reveal useful information about the amount of liquid water accumulated within gas flow channels. Liquid water transport though GDL has also been investigated by measuring the liquid water breakthrough pressure for the region between the capillary fingering and the stable displacement on the drainage phase diagram. The breakthrough pressure has been measured for different variables such as GDL thickness, PTFE/Nafion content within the GDL, GDL compression, the inclusion of a micro-porous layer (MPL), and different water flow rates through the GDL. Prior to all these studies, GDL microstructural properties have been studied. GDL microstructural properties such as mean pore diameter, pore diameter distribution, and pore roundness distribution have been investigated by analyzing SEM images of GDL samples.