Liquid–Liquid Equilibria of Water + Ethanol + 1-Butyl-3-methylimidazolium Bis(trifluoromethanesulfonyl)imide Ternary System: Measurements and Correlation at Different Temperatures

In this work authors present the experimental liquid–liquid equilibria (LLE) data of water + ethanol + 1-butyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide ([bmim][Tf2N]) system at different temperatures. The LLE of the system was obtained in the temperature range from 283.2 to 323.2 K. The nonrandom two liquid (NRTL) and universal quasichemical (UNIQUAC) models were used to correlate ternary systems. The equilibrium compositions were successfully correlated by the interaction parameters from both models, however UNIQUAC gave a more accurate correlation. Finally, a study about the solvent capability of ionic liquid was made in order to evaluate the possibility of separating the mixture formed by ethanol and water using that ionic liquid.

Linear Stability Analysis of Convection in Two-Layer System with an Evaporating Vapor-Liquid Interface

Classical theories have successfully provided an explanation for convection in a liquid layer heated from below without evaporation. However, these theories are inadequate to account for the convective instabilities in an evaporating liquid layer, especially in the case when it is cooled from below. In the present paper, we study the onset of Marangoni convection in a liquid layer being overlain by a vapor layer.A new two-sided model is put forward instead of the one-sided model in previous studies. Marangoni-Bénard instabilities in evaporating liquid thin layers are investigated with a linear instability analysis. We define a new evaporation Biot number, which is different from that in previous studies and discuss the influences of reference evaporating velocity and evaporation Biot number on the vapor-liquid system. At the end, we explain why the instability occurs even when an evaporating liquid layer is cooled from below.

Two fluid numerical model for coastal applications

Wave breaking is an important coastal process, influencing hydro-morphodynamic processes such as turbulence generation and wave energy dissipation, run-up on the beach and overtopping of coastal defence structures. During breaking, waves are complex mixtures of air and water (“white water”) whose properties affect velocity and pressure fields in the vicinity of the free surface and, depending on the breaker characteristics, different mechanisms for air entrainment are usually observed. Several laboratory experiments have been performed to investigate the role of air bubbles in the wave breaking process (Chanson & Cummings, 1994, among others) and in wave loading on vertical wall (Oumeraci et al., 2001; Peregrine et al., 2006, among others), showing that the air phase is not negligible since the turbulent energy dissipation involves air-water mixture. The recent advancement of numerical models has given valuable insights in the knowledge of wave transformation and interaction with coastal structures. Among these models, some solve the RANS equations coupled with a free-surface tracking algorithm and describe velocity, pressure, turbulence and vorticity fields (Lara et al. 2006 a-b, Clementi et al., 2007). The single-phase numerical model, in which the constitutive equations are solved only for the liquid phase, neglects effects induced by air movement and trapped air bubbles in water. Numerical approximations at the free surface may induce errors in predicting breaking point and wave height and moreover, entrapped air bubbles and water splash in air are not properly represented. The aim of the present thesis is to develop a new two-phase model called COBRAS2 (stands for Cornell Breaking waves And Structures 2 phases), that is the enhancement of the single-phase code COBRAS0, originally developed at Cornell University (Lin & Liu, 1998). In the first part of the work, both fluids are considered as incompressible, while the second part will treat air compressibility modelling. The mathematical formulation and the numerical resolution of the governing equations of COBRAS2 are derived and some model-experiment comparisons are shown. In particular, validation tests are performed in order to prove model stability and accuracy. The simulation of the rising of a large air bubble in an otherwise quiescent water pool reveals the model capability to reproduce the process physics in a realistic way. Analytical solutions for stationary and internal waves are compared with corresponding numerical results, in order to test processes involving wide range of density difference. Waves induced by dam-break in different scenarios (on dry and wet beds, as well as on a ramp) are studied, focusing on the role of air as the medium in which the water wave propagates and on the numerical representation of bubble dynamics. Simulations of solitary and regular waves, characterized by both spilling and plunging breakers, are analyzed with comparisons with experimental data and other numerical model in order to investigate air influence on wave breaking mechanisms and underline model capability and accuracy. Finally, modelling of air compressibility is included in the new developed model and is validated, revealing an accurate reproduction of processes. Some preliminary tests on wave impact on vertical walls are performed: since air flow modelling allows to have a more realistic reproduction of breaking wave propagation, the dependence of wave breaker shapes and aeration characteristics on impact pressure values is studied and, on the basis of a qualitative comparison with experimental observations, the numerical simulations achieve good results.

Mass transfer characteristics of two-aqueous-phase liquid-liquid mixtures

Mass transfer rates were studied using the falling drop method. Cibacron Blue 3 GA dye was the transferring solute from the salt phase to the PEG phase. Measurements were undertaken for several concentrations of the dye and the phase-forming solutes and with a range of different drop sizes, e.g. 2.8, 3.0 and 3.7 mm. The dye was observed to be present in the salt phase as finely dispersed solids but a model confirmed that the mass transfer process could still be described by an equation based upon the Whitman two-film model. The overall mass transfer coefficient increased with increasing concentration of the dye. The apparent mass transfer coefficient ranged from 1 x 10-5 to 2 x 10 -4 m/s. Further experiments suggested that mass transfer was enhanced at high concentration by several mechanisms. The dye was found to change the equilibrium composition of the two phases, leading to transfer of salt between the drop and continuous phases. It also lowered the interfacial tension (i.e. from 1.43 x 10-4 N/m for 0.01% w/w dye concentration to 1.07 x 10-4 N/m for 0.2% w/w dye concentration) between the two phases, which could have caused interfacial instabilities (Marangoni effects). The largest drops were deformable, which resulted in a significant increase in the mass transfer rate. Drop size distribution and Sauter mean drop diameter were studied on-line in a 1 litre agitated vessel using a laser diffraction technique. The effects of phase concentration, dispersed phase hold-up and impeller speed were investigated for the salt-PEG system. An increase in agitation speed in the range 300 rpm to 1000 rpm caused a decrease in mean drop diameter, e.g. from 50 m to 15 m. A characteristic bimodal drop size distribution was established within a very short time. An increase in agitation rate caused a shift of the larger drop size peak to a smaller size.

Dynamics of transformation and fragmentation of composite liquid nano-particles

Recent research on particle size distributions and particle concentrations near a busy road cannot be explained by the conventional mechanisms for particle evolution of combustion aerosols. Specifically they appear to be inadequate to explain the experimental observations of particle transformation and the evolution of the total number concentration. This resulted in the development of a new mechanism based on their thermal fragmentation, for the evolution of combustion aerosol nano-particles. A complex and comprehensive pattern of evolution of combustion aerosols, involving particle fragmentation, was then proposed and justified. In that model it was suggested that thermal fragmentation occurs in aggregates of primary particles each of which contains a solid graphite/carbon core surrounded by volatile molecules bonded to the core by strong covalent bonds. Due to the presence of strong covalent bonds between the core and the volatile (frill) molecules, such primary composite particles can be regarded as solid, despite the presence of significant (possibly, dominant) volatile component. Fragmentation occurs when weak van der Waals forces between such primary particles are overcome by their thermal (Brownian) motion. In this work, the accepted concept of thermal fragmentation is advanced to determine whether fragmentation is likely in liquid composite nano-particles. It has been demonstrated that at least at some stages of evolution, combustion aerosols contain a large number of composite liquid particles containing presumably several components such as water, oil, volatile compounds, and minerals. It is possible that such composite liquid particles may also experience thermal fragmentation and thus contribute to, for example, the evolution of the total number concentration as a function of distance from the source. Therefore, the aim of this project is to examine theoretically the possibility of thermal fragmentation of composite liquid nano-particles consisting of immiscible liquid v components. The specific focus is on ternary systems which include two immiscible liquid droplets surrounded by another medium (e.g., air). The analysis shows that three different structures are possible, the complete encapsulation of one liquid by the other, partial encapsulation of the two liquids in a composite particle, and the two droplets separated from each other. The probability of thermal fragmentation of two coagulated liquid droplets is discussed and examined for different volumes of the immiscible fluids in a composite liquid particle and their surface and interfacial tensions through the determination of the Gibbs free energy difference between the coagulated and fragmented states, and comparison of this energy difference with the typical thermal energy kT. The analysis reveals that fragmentation was found to be much more likely for a partially encapsulated particle than a completely encapsulated particle. In particular, it was found that thermal fragmentation was much more likely when the volume ratio of the two liquid droplets that constitute the composite particle are very different. Conversely, when the two liquid droplets are of similar volumes, the probability of thermal fragmentation is small. It is also demonstrated that the Gibbs free energy difference between the coagulated and fragmented states is not the only important factor determining the probability of thermal fragmentation of composite liquid particles. The second essential factor is the actual structure of the composite particle. It is shown that the probability of thermal fragmentation is also strongly dependent on the distance that each of the liquid droplets should travel to reach the fragmented state. In particular, if this distance is larger than the mean free path for the considered droplets in the air, the probability of thermal fragmentation should be negligible. In particular, it follows form here that fragmentation of the composite particle in the state with complete encapsulation is highly unlikely because of the larger distance that the two droplets must travel in order to separate. The analysis of composite liquid particles with the interfacial parameters that are expected in combustion aerosols demonstrates that thermal fragmentation of these vi particles may occur, and this mechanism may play a role in the evolution of combustion aerosols. Conditions for thermal fragmentation to play a significant role (for aerosol particles other than those from motor vehicle exhaust) are determined and examined theoretically. Conditions for spontaneous transformation between the states of composite particles with complete and partial encapsulation are also examined, demonstrating the possibility of such transformation in combustion aerosols. Indeed it was shown that for some typical components found in aerosols that transformation could take place on time scales less than 20 s. The analysis showed that factors that influenced surface and interfacial tension played an important role in this transformation process. It is suggested that such transformation may, for example, result in a delayed evaporation of composite particles with significant water component, leading to observable effects in evolution of combustion aerosols (including possible local humidity maximums near a source, such as a busy road). The obtained results will be important for further development and understanding of aerosol physics and technologies, including combustion aerosols and their evolution near a source.

Intermittent Dynamics, Stochastic Resonance and Dynamical Heterogeneity in Supercooled Liquid Water

Here we find through computer simulations and theoretical analysis that the low temperature thermodynamic anomalies of liquid water arises from the intermittent fluctuation between its high density and low density forms, consisting largely of 5-coordinated and 4-coordinated water molecules, respectively. The fluctuations exhibit strong dynamic heterogeneity (defined by the four point time correlation function), accompanied by a divergence like growth of the dynamic correlation length, of the type encountered in fragile supercooled liquids. The intermittency has been explained by invoking a two state model often employed to understand stochastic resonance, with the relevant periodic perturbation provided here by the fluctuation of the total volume of the system.

Constructing a new individual-based model of phosphine resistance in lesser grain borer (Rhyzopertha dominica): do we need to include two loci rather than one?

In this article, we describe and compare two individual-based models constructed to investigate how genetic factors influence the development of phosphine resistance in lesser grain borer (R. dominica). One model is based on the simplifying assumption that resistance is conferred by alleles at a single locus, while the other is based on the more realistic assumption that resistance is conferred by alleles at two separate loci. We simulated the population dynamic of R. dominica in the absence of phosphine fumigation, and under high and low dose phosphine treatments, and found important differences between the predictions of the two models in all three cases. In the absence of fumigation, starting from the same initial frequencies of genotypes, the two models tended to different stable frequencies, although both reached Hardy-Weinberg equilibrium. The one-locus model exaggerated the equilibrium proportion of strongly resistant beetles by 3.6 times, compared to the aggregated predictions of the two-locus model. Under a low dose treatment the one-locus model overestimated the proportion of strongly resistant individuals within the population and underestimated the total population numbers compared to the two-locus model. These results show the importance of basing resistance evolution models on realistic genetics and that using oversimplified one-locus models to develop pest control strategies runs the risk of not correctly identifying tactics to minimise the incidence of pest infestation.

Real-time kinetic analysis of hCG-monoclonal antibody interaction using radiolabeled hCG probe: presence of two forms of Ag-mAb complex as revealed by solid phase dissociation studies

Real-time kinetics of ligand-ligate interaction has predominantly been studied by either fluorescence or surface plasmon resonance based methods. Almost all such studies are based on association between the ligand and the ligate. This paper reports our analysis of dissociation data of monoclonal antibody-antigen (hCG) system using radio-iodinated hCG as a probe and nitrocellulose as a solid support to immobilize mAb. The data was analyzed quantitatively for a one-step and a two-step model. The data fits well into the two-step model. We also found that a fraction of what is bound is non-dissociable (tight-binding portion (TBP)). The TBP was neither an artifact of immobilization nor does it interfere with analysis. It was present when the reaction was carried out in homogeneous solution in liquid phase. The rate constants obtained from the two methods were comparable. The work reported here shows that real-time kinetics of other ligand-ligate interaction can be studied using nitrocellulose as a solid support. (C) 2002 Elsevier Science B.V. All rights reserved.

Vaporization and collision modeling of liquid fuel sprays in a co-axial fuel and air pre-mixer

Droplet collision occurs frequently in regions where the droplet number density is high. Even for Lean Premixed and Pre-vaporized (LPP) liquid sprays, the collision effects can be very high on the droplet size distributions, which will in turn affect the droplet vaporization process. Hence, in conjunction with vaporization modeling, collision modeling for such spray systems is also essential. The standard O'Rourke's collision model, usually implemented in CFD codes, tends to generate unphysical numerical artifact when simulations are performed on Cartesian grid and the results are not grid independent. Thus, a new collision modeling approach based on no-time-counter method (NTC) proposed by Schmidt and Rutland is implemented to replace O'Rourke's collision algorithm to solve a spray injection problem in a cylindrical coflow premixer. The so called ``four-leaf clover'' numerical artifacts are eliminated by the new collision algorithm and results from a diesel spray show very good grid independence. Next, the dispersion and vaporization processes for liquid fuel sprays are simulated in a coflow premixer. Two liquid fuels under investigation are jet-A and Rapeseed Methyl Esters (RME). Results show very good grid independence in terms of SMD distribution, droplet number distribution and fuel vapor mass flow rate. A baseline test is first established with a spray cone angle of 90 degrees and injection velocity of 3 m/s and jet-A achieves much better vaporization performance than RME due to its higher vapor pressure. To improve the vaporization performance for both fuels, a series of simulations have been done at several different combinations of spray cone angle and injection velocity. At relatively low spray cone angle and injection velocity, the collision effect on the average droplet size and the vaporization performance are very high due to relatively high coalescence rate induced by droplet collisions. Thus, at higher spray cone angle and injection velocity, the results expectedly show improvement in fuel vaporization performance since smaller droplet has a higher vaporization rate. The vaporization performance and the level of homogeneity of fuel-air mixture can be significantly improved when the dispersion level is high, which can be achieved by increasing the spray cone angle and injection velocity. (C) 2012 Elsevier Ltd. All rights reserved.

Homogeneous isotropic superfluid turbulence in two dimensions: Inverse and forward cascades in the Hall-Vinen-Bekharevich-Khalatnikov model

We present the first direct-numerical-simulation study of the statistical properties of two-dimensional superfluid turbulence in the simplified, Hall-Vinen-Bekharevich-Khalatnikov two-fluid model. We show that both normalfluid and superfluid energy spectra can exhibit two power-law regimes, the first associated with an inverse cascade of energy and the second with the forward cascade of enstrophy. We quantify the mutual-friction-induced alignment of normal and superfluid velocities by obtaining probability distribution functions of the angle between them and the ratio of their moduli.

Marangoni-Benard instability with the exchange of evaporation at liquid-vapour interface

A new two-sided model rather than the one-sided model in previous works is put forward. The linear instability analysis is performed on the Marangoni-Benard convection in the two-layer system with an evaporation interface. We define a new evaporation Biot number which is different from that in the one-sided model, and obtain the curves of critical Marangoni number versus wavenumber. The influence of evaporation velocity and Biot number on the system is discussed and a new phenomenon uninterpreted before is now explained from our numerical results.

Hydrothermal Wave in a Shallow Liquid Layer

The oscillatory thermocapillary convection and hydrothermal wave in a shallow liquid layer, where a temperature difference is applied between two parallel sidewalls, have been numerically investigated in a two-dimensional model. The oscillatory thermocapillary convection and hydrothermal wave appear if the Marangoni number is larger than a critical value. The critical phase speed and critical wave number of the hydrothermal wave agree with the ones given analytically by Smith and Davis in the microgravity environment, and it travels in the direction opposed to the surface flow. Another wave traveled downstream in addition to the hydrothermal wave traveled upstream was observed in the case of earth gravity condition.

Modeling and Simulation of Liquid Pulsed Particulate Fluidized Beds

Pulsed fluidization is of considerable interest in process engineering for improving fluidization quality. Quantitative understanding of the pulsed two-phase flow behaviors is very important for proper design and optimum operation of such contactors. The

Numerical Simulation of Pulsed Liquid Fluidized Bed and Its Experimental Validation

The pulsed liquid fluidized bed was studied using numerical simulation and experimental methods, The area-averaged two-fluid model (TFM) was used to simulate the pulsed fluidization. The bed expansion and collapse processes were simulated first and the phenomena obtained from the calculation were consistent with our previous experiments and observations. In the pulsed fluidization, the variation of bed height, the variations of particle velocity and concentration distribution were obtained and analyzed. Experiments were carried out to validate the simulation results. The pressure variation with time at different locations was measured using pressure transducers and compared with the simulated results. The variations of bed height and particle concentration distribution were recorded using a digital video camera recorder. The results were consistent with the simulation results as a whole.

Oil-Water Two-Phase Flow Inside T-Junction

The oil/water two-phase flow inside T-junctions was numerically simulated with a 3-D two-fluid model, and the turbulence was described using the mixture k - epsilon model. Some experiments of oil/water flow inside a single T-junction were conducted in the laboratory. The results show that the separating performance of T-junction largely depends oil the inlet volumetric fraction and flow patterns. A reasonable agreement is reached between the numerical simulation and the experiments for both the oil fraction distribution and the separation efficiency.