Analytical and numerical study of the electro-osmotic annular flow of viscoelastic fluids

In this work we present semi-analytical solutions for the electro-osmotic annular flow of viscoelastic fluids modeled by the Linear and Exponential PTT models. The viscoelastic fluid flows in the axial direction between two concentric cylinders under the combined influences of electrokinetic and pressure forcings. The analysis invokes the Debye-Hückel approximation and includes the limit case of pure electro-osmotic flow. The solution is valid for both no slip and slip velocity at the walls and the chosen slip boundary condition is the linear Navier slip velocity model. The combined effects of fluid rheology, electro-osmotic and pressure gradient forcings on the fluid velocity distribution are also discussed.

Driven fragmentation of granular gases

The dynamics of homogeneously heated granular gases which fragment due to particle collisions is analyzed. We introduce a kinetic model which accounts for correlations induced at the grain collisions and analyze both the kinetics and relevant distribution functions these systems develop. The work combines analytical and numerical studies based on direct simulation Monte Carlo calculations. A broad family of fragmentation probabilities is considered, and its implications for the system kinetics are discussed. We show that generically these driven materials evolve asymptotically into a dynamical scaling regime. If the fragmentation probability tends to a constant, the grain number diverges at a finite time, leading to a shattering singularity. If the fragmentation probability vanishes, then the number of grains grows monotonously as a power law. We consider different homogeneous thermostats and show that the kinetics of these systems depends weakly on both the grain inelasticity and driving. We observe that fragmentation plays a relevant role in the shape of the velocity distribution of the particles. When the fragmentation is driven by local stochastic events, the longvelocity tail is essentially exponential independently of the heating frequency and the breaking rule. However, for a Lowe-Andersen thermostat, numerical evidence strongly supports the conjecture that the scaled velocity distribution follows a generalized exponential behavior f (c)~exp (−cⁿ), with n ≈1.2, regarding less the fragmentation mechanisms

Aspereza equivalente de ríos con contorno granular grueso

El concepto de aspereza equivalente (ks) permite relacionar la resistencia al flujo de cauces de contorno granular con el tamaño de las partículas que lo conforman, mediante la ley logarítmica de distribución vertical de velocidad de la corriente. En este artículo se revisa el estado del conocimiento acerca de los criterios para la predicción de dicha relación con el objetivo de brindar una guía para la selección de los mismos. Pese a que el análisis de las experiencias compiladas desvela la dispersión de resultados, a efectos prácticos se recomienda para ríos de grava y de montaña: ks 3•d90; ks 3,4•d84 y ks 7•d50. Dicha dispersión y el hecho de que ks sea varias veces superior al diámetro medio del sedimento puede atribuirse en diferente grado: a la heterogeneidad del sedimento en lechos naturales, a las formas de fondo, al transporte sólido de fondo,así como a las limitaciones del modelo logarítmico de distribución de velocidad bajo ciertas condiciones de flujo. La predicción de ks en función de variables estadísticas del campo de elevaciones del lecho se ha revelado como una alternativa con una elevada capacidad explicativa, por lo que en la medida que avancen las técnicas microtopográficas se consolidará esta vía como el método del futuro.

Selluradan leijutus ja hallinta ilmakuivaimessa

Diplomityön tavoitteena oli tutkia selluradan leijutusta sekä hallintaa kapealla ja täysleveällä selluradalla ilmakuivaimessa. Työn sisältö on jaettu neljään osaan teoriaan, numeeriseen virtauslaskentaan, kokeelliseen osuuteen ja johtopäätöksiin. Työn alussa kiinnitettiin huomiota ilmakuivaimen ja selluradan rakenteeseen, selluradan lujuusominaisuuksiin ja niiden testaamiseen laboratorio-oloissa. Teoria osuudessa tarkasteltiin selluradan leijutusteoriaa ja siihen liittyvää virtaustekniikkaa. Leijutusta kuvattiin mekaniikasta tutun vapaakappalekuvan muodossa. Vapaakappalekuva auttaa ymmärtämään sellurataan vaikuttavia voimia ja niiden aiheuttajia. Numeerisessa virtauslaskennassa rakennettiin Fluent -virtauslaskentamalli. Malli kuvaa leijutussuuttimen alusta 0,2 metrin pituista pätkää. Laskennassa simuloitiin selluradan lei-jutusta 2, 5 ja 7 mm etäisyyksillä. Tuloksia tarkasteltiin paine- ja nopeusjakaumien muodossa. Kokeellisessa osuudessa tehtiin mittauksia yksittäisellä ylä- ja alasuuttimella sekä suutinpaketilla. Yksittäisistä suuttimista mitattiin ilmasuihkujen ulostulokulmaa, suuttimensisäistä painejakaumaa ja selluradan pintaan vaikuttavaa staattista painetta suuttimen pituuden funktiona. Suutinpaketilla tutkittiin ilmasuihkujen ja tasomaisen levyn välistä vuorovaikutusta. Ylä- ja alasuuttimen välissä leijutettiin levyä, josta mitattiin ilmasuihkujen aiheuttama z-suuntainen voima ja leijutuskorkeus. Työn lopussa arvioitiin teorian, numeerisen virtauslaskennan ja kokeellisen osuuden tuloksia. Näiden pohjalta tehtiin teknisiä parannusehdotuksia ilmakuivaimen toiminnan tehostamiseksi ja selluradan hallinnan parantamiseksi kapealla ja täysleveällä selluradalla.

The internal wave field in Sau reservoir : Observation and modeling of a third vertical mode

Water withdrawal from Mediterranean reservoirs in summer is usually very high. Because of this, stratification is often continuous and far from the typical two-layered structure, favoring the excitation of higher vertical modes. The analysis of wind, temperature, and current data from Sau reservoir (Spain) shows that the third vertical mode of the internal seiche (baroclinic mode) dominated the internal wave field at the beginning of September 2003. We used a continuous stratification two-dimensional model to calculate the period and velocity distribution of the various modes of the internal seiche, and we calculated that the period of the third vertical mode is ;24 h, which coincides with the period of the dominating winds. As a result of the resonance between the third mode and the wind, the other oscillation modes were not excited during this period

Synthesis of percarboxylic acids in microreactor

Percarboxylic acids are commonly used as disinfection and bleaching agents in textile, paper, and fine chemical industries. All of these applications are based on the oxidative potential of these compounds. In spite of high interest in these chemicals, they are unstable and explosive chemicals, which increase the risk of synthesis processes and transportation. Therefore, the safety criteria in the production process should be considered. Microreactors represent a technology that efficiently utilizes safety advantages resulting from small scale. Therefore, microreactor technology was used in the synthesis of peracetic acid and performic acid. These percarboxylic acids were produced at different temperatures, residence times and catalyst i.e. sulfuric acid concentrations. Both synthesis reactions seemed to be rather fast because with performic acid equilibrium was reached in 4 min at 313 K and with peracetic acid in 10 min at 343 K. In addition, the experimental results were used to study the kinetics of the formation of performic acid and peracetic acid. The advantages of the microreactors in this study were the efficient temperature control even in very exothermic reaction and good mixing due to the short diffusion distances. Therefore, reaction rates were determined with high accuracy. Three different models were considered in order to estimate the kinetic parameters such as reaction rate constants and activation energies. From these three models, the laminar flow model with radial velocity distribution gave most precise parameters. However, sulfuric acid creates many drawbacks in this synthesis process. Therefore, a ´´greener´´ way to use heterogeneous catalyst in the synthesis of performic acid in microreactor was studied. The cation exchange resin, Dowex 50 Wx8, presented very high activity and a long life time in this reaction. In the presence of this catalyst, the equilibrium was reached in 120 second at 313 K which indicates a rather fast reaction. In addition, the safety advantages of microreactors were investigated in this study. Four different conventional methods were used. Production of peracetic acid was used as a test case, and the safety of one conventional batch process was compared with an on-site continuous microprocess. It was found that the conventional methods for the analysis of process safety might not be reliable and adequate for radically novel technology, such as microreactors. This is understandable because the conventional methods are partly based on experience, which is very limited in connection with totally novel technology. Therefore, one checklist-based method was developed to study the safety of intensified and novel processes at the early stage of process development. The checklist was formulated using the concept of layers of protection for a chemical process. The traditional and three intensified processes of hydrogen peroxide synthesis were selected as test cases. With these real cases, it was shown that several positive and negative effects on safety can be detected in process intensification. The general claim that safety is always improved by process intensification was questioned.

Emission characteristics and dynamics of C2 from laser produced graphite plasma

The emission features of laser ablated graphite plume generated in a helium ambient atmosphere have been investigated with time and space resolved plasma diagnostic technique. Time resolved optical emission spectroscopy is employed to reveal the velocity distribution of different species ejected during ablation. At lower values of laser fluences only a slowly propagating component of C2 is seen. At high fluences emission from C2 shows a twin peak distribution in time. The formation of an emission peak with diminished time delay giving an energetic peak at higher laser fluences is attributed to many body recombination. It is also observed that these double peaks get modified into triple peak time of flight distribution at distances greater than 16 mm from the target. The occurrence of multiple peaks in the C2 emission is mainly due to the delays caused from the different formation mechanism of C2 species. The velocity distribution of the faster peak exhibits an oscillating character with distance from the target surface.

Numerical Simulation of Electroosmotic Flow with Step Change in Zeta Potential

Electroosmotic flow is a convenient mechanism for transporting polar fluid in a microfluidic device. The flow is generated through the application of an external electric field that acts on the free charges that exists in a thin Debye layer at the channel walls. The charge on the wall is due to the chemistry of the solid-fluid interface, and it can vary along the channel, e.g. due to modification of the wall. This investigation focuses on the simulation of the electroosmotic flow (EOF) profile in a cylindrical microchannel with step change in zeta potential. The modified Navier-Stoke equation governing the velocity field and a non-linear two-dimensional Poisson-Boltzmann equation governing the electrical double-layer (EDL) field distribution are solved numerically using finite control-volume method. Continuities of flow rate and electric current are enforced resulting in a non-uniform electrical field and pressure gradient distribution along the channel. The resulting parabolic velocity distribution at the junction of the step change in zeta potential, which is more typical of a pressure-driven velocity flow profile, is obtained.

Adaptation-induced collective dynamics of a single-cell protozoan

We investigate the behavior of a single-cell protozoan in a narrow tubular ring. This environment forces them to swim under a one-dimensional periodic boundary condition. Above a critical density, single-cell protozoa aggregate spontaneously without external stimulation. The high-density zone of swimming cells exhibits a characteristic collective dynamics including translation and boundary fluctuation. We analyzed the velocity distribution and turn rate of swimming cells and found that the regulation of the turing rate leads to a stable aggregation and that acceleration of velocity triggers instability of aggregation. These two opposing effects may help to explain the spontaneous dynamics of collective behavior. We also propose a stochastic model for the mechanism underlying the collective behavior of swimming cells.

Location and characteristics of the reconnection X line deduced from low-altitude satellite and ground-based observations: 2. Defense Meteorological Satellite Program and European Incoherent Scatter data

We present an analysis of a cusp ion step, observed by the Defense Meteorological Satellite Program (DMSP) F10 spacecraft, between two poleward moving events of enhanced ionospheric electron temperature, observed by the European Incoherent Scatter (EISCAT) radar. From the ions detected by the satellite, the variation of the reconnection rate is computed for assumed distances along the open-closed field line separatrix from the satellite to the X line, do. Comparison with the onset times of the associated ionospheric events allows this distance to be estimated, but with an uncertainty due to the determination of the low-energy cutoff of the ion velocity distribution function, ƒ(ν). Nevertheless, the reconnection site is shown to be on the dayside magnetopause, consistent with the reconnection model of the cusp during southward interplanetary magnetic field (IMF). Analysis of the time series of distribution function at constant energies, ƒ(ts), shows that the best estimate of the distance do is 14.5±2 RE. This is consistent with various magnetopause observations of the signatures of reconnection for southward IMF. The ion precipitation is used to reconstruct the field-parallel part of the Cowley D ion distribution function injected into the open low-latitude boundary layer in the vicinity of the X line. From this reconstruction, the field-aligned component of the magnetosheath flow is found to be only −55±65 km s−1 near the X line, which means either that the reconnection X line is near the stagnation region at the nose of the magnetosphere, or that it is closely aligned with the magnetosheath flow streamline which is orthogonal to the magnetosheath field, or both. In addition, the sheath Alfvén speed at the X line is found to be 220±45 km s−1, and the speed with which newly opened field lines are ejected from the X line is 165±30 km s−1. We show that the inferred magnetic field, plasma density, and temperature of the sheath near the X line are consistent with a near-subsolar reconnection site and confirm that the magnetosheath field makes a large angle (>58°) with the X line.

EISCAT observations of ion composition and temperature anisotropy in the high-latitude F-region

The papers by Winser et al. [(1990) J. atmos. terr. Phys.52, 501] and Häggström and Collis [(1990) J. atmos. terr. Phys.52, 519] used plasma flows and ion temperatures, as measured by the EISCAT tristatic incoherent scatter radar, to investigate changes in the ion composition of the ionospheric F-layer at high latitudes, in response to increases in the speed of plasma convection. These studies reported that the ion composition rapidly changed from mainly O+ to almost completely (>90%) molecular ions, following rapid increases in ion drift speed by >1 km s−1. These changes appeared inconsisent with theoretical considerations of the ion chemistry, which could not account for the large fractions of molecular ions inferred from the obsevations. In this paper, we discuss two causes of this discrepancy. First, we reevaluate the theoretical calculations for chemical equilibrium and show that, if we correct the derived temperatures for the effect of the molecular ions, and if we employ more realistic dependences of the reaction rates on the ion temperature, the composition changes derived for the faster convection speeds can be explained. For the Winser et al. observations with the radar beam at an aspect angle of ϕ = 54.7° to the geomagnetic field, we now compute a change to 89% molecular ions in < 2 min, in response to the 3 km s−1 drift. This is broadly consistent with the observations. But for the two cases considered by Häggström and Collis, looking along the field line (ϕ = 0°), we compute the proportion of molecular ions to be only 4 and 16% for the observed plasma drifts of 1.2 and 1.6 km s−1, respectively. These computed proportions are much smaller than those derived experimentally (70 and 90%). We attribute the differences to the effects of non-Maxwellian, anisotropic ion velocity distribution functions. We also discuss the effect of ion composition changes on the various radar observations that report anisotropies of ion temperature.

Measuring ion temperatures and studying the ion energy balance in the high-latitude ionosphere

Data are presented for a nighttime ion heating event observed by the EISCAT radar on 16 December 1988. In the experiment, the aspect angle between the radar beam and the geomagnetic field was fixed at 54.7°, which avoids any ambiguity in derived ion temperature caused by anisotropy in the ion velocity distribution function. The data were analyzed with an algorithm which takes account of the non-Maxwellian line-of-sight ion velocity distribution. During the heating event, the derived spectral distortion parameter (D∗) indicated that the distribution function was highly distorted from a Maxwellian form when the ion drift increased to 4 km s−1. The true three-dimensional ion temperature was used in the simplified ion balance equation to compute the ion mass during the heating event. The ion composition was found to change from predominantly O4 to mainly molecular ions. A theoretical analysis of the ion composition, using the MSIS86 model and published values of the chemical rate coefficients, accounts for the order-of-magnitude increase in the atomic/molecular ion ratio during the event, but does not successfully explain the very high proportion of molecular ions that was observed.

Analysis of incoherent scatter radar data from non-thermal F-region plasma

A procedure is presented for fitting incoherent scatter radar data from non-thermal F-region ionospheric plasma, using theoretical spectra previously predicted. It is found that values of the shape distortion factor D∗, associated with deviations of the ion velocity distribution from a Maxwellian distribution, and ion temperatures can be deduced (the results being independent of the path of iteration) if the angle between the line-of-sight and the geomagnetic field is larger than about 15–20°. The procedure can be used with one or both of two sets of assumptions. These concern the validity of the adopted model for the line-of-sight ion velocity distribution in the one case or for the full three-dimensional ion velocity distribution function in the other. The distribution function employed was developed to describe the line-of-sight velocity distribution for large aspect angles, but both experimental data and Monte Carlo simulations indicate that the form of the field-perpendicular distribution can also describe the distribution at more general aspect angles. The assumption of this form for the line-of-sight velocity distribution at a general aspect angle enables rigorous derivation of values of the one-dimensional, line-of-sight ion temperature. With some additional assumptions (principally that the field-parallel distribution is always Maxwellian and there is a simple relationship between the ion temperature anisotropy and the distortion of the field-perpendicular distribution from a Maxwellian), fits to data for large aspect angles enable determination of line-of-sight temperatures at all aspect angles and hence, of the average ion temperature and the ion temperature anisotropy. For small aspect angles, the analysis is restricted to the determination of the line-of-sight ion temperature because the theoretical spectrum is insensitive to non-thermal effects when the plasma is viewed along directions almost parallel to the magnetic field. This limitation is expected to apply to any realistic model of the ion velocity distribution function and its consequences are discussed. Fit strategies which allow for mixed ion composition are also considered. Examples of fits to data from various EISCAT observing programmes are presented.

Radar observations of non-thermal plasmas at different aspect angles

Data are presented from the EISCAT CP-3-E experiment which show the presence of non-thermal plasma over a range of latitudes. The O+ ion-velocity distribution function is almost toroidal when the electric field reaches values of 125 mV m−1. The ion temperature derived from such data assuming a Maxwellian distribution function will overestimate the true ion temperature when the observing angle is large with respect to the magnetic field, and underestimate the temperature when the aspect angle is small. When the expressions for the distribution function are extended to include mixed ion composition, an improvement is sometimes found in fitting the observed data, and estimates of the composition can be made. Such an analysis suggests that N2+ can occasionally form a significant part of the total ion density in a narrow height region centred at 275 km.

Analysis of incoherent scatter spectra from non-Maxwellian plasma

Incoherent scatter data from non-thermal F-region ionospheric plasma are analysed, using theoretical spectra predicted by Raman et al. It is found that values of the semi-empirical drift parameter D∗, associated with deviations of the ion velocity distribution from a Maxwellian, and the plasma temperatures can be rigorously deduced (the results being independent of the path of iteration) if the angle between the line-of-sight and the geomagnetic field is larger than about 15–20 degrees. For small aspect angles, the deduced value of the average (or 3-D) ion temperature remains ambiguous and the analysis is restricted to the determination of the line-of-sight temperature because the theoretical spectrum is insensitive to non-thermal effects when the plasma is viewed along directions almost parallel to the magnetic field. This limitation is expected to apply to any realistic model of the ion velocity distribution, and its consequences are discussed. Fit strategies which allow for mixed ion composition are also considered. Examples of fits to data from various EISCAT observing programmes are presented.