Simulations of droplet dispersion in the wakes of cylinders and icing tunnel spray bars

Simulations of droplet dispersion behind cylinder wakes and downstream of icing tunnel spray bars were conducted. In both cases, a range of droplet sizes were investigated numerically with a Lagrangian particle trajectory approach while the turbulent air flow was investigated with a hybrid Reynolds-Averaged Navier-Stokes/Large-Eddy Simulations approach scheme. In the first study, droplets were injected downstream of a cylinder at sub-critical conditions (i.e. with laminar boundary layer separation). A stochastic continuous random walk (CRW) turbulence model was used to capture the effects of sub-grid turbulence. Small inertia droplets (characterized by small Stokes numbers) were affected by both the large-scale and small-scale vortex structures and closely followed the air flow, while exhibiting a dispersion consistent with that of a scalar flow field. Droplets with intermediate Stokes numbers were centrifuged by the vortices to the outer edges of the wake, yielding an increased dispersion. Large Stokes number droplets were found to be less responsive to the vortex structures and exhibited the least dispersion. Particle concentration was also correlated with vorticity distribution which yielded preferential bias effects as a function of different particle sizes. This trend was qualitatively similar to results seen in homogenous isotropic turbulence, though the influence of particle inertia was less pronounced for the cylinder wake case. A similar study was completed for droplet dispersion within the Icing Research Tunnel (IRT) at the NASA Glenn Research Center, where it is important to obtain a nearly uniform liquid water content (LWC) distribution in the test section (to recreate atmospheric icing conditions).. For this goal, droplets are diffused by the mean and turbulent flow generated from the nozzle air jets, from the upstream spray bars, and from the vertical strut wakes. To understand the influence of these three components, a set of simulations was conducted with a sequential inclusion of these components. Firstly, a jet in an otherwise quiescent airflow was simulated to capture the impact of the air jet on flow turbulence and droplet distribution, and the predictions compared well with experimental results. The effects of the spray bar wake and vertical strut wake were then included with two more simulation conditions, for which it was found that the air jets were the primary driving force for droplet dispersion, i.e. that the spray bar and vertical strut wake effects were secondary.

The experimental determination of impact ionization coefficients in gallium arsenide and indium phosphide

U of I Only

Heat dissipation from carbon nano-electronics

The incorporation of graphitic compounds such as carbon nanotubes (CNTs) and graphene into nano-electronic device packaging holds much promise for waste heat management given their high thermal conductivities. However, as these graphitic materials must be used in together with other semiconductor/insulator materials, it is not known how thermal transport is affected by the interaction. Using different simulation techniques, in this thesis, we evaluate the thermal transport properties - thermal boundary conductance (TBC) and thermal conductivity - of CNTs and single-layer graphene in contact with an amorphous SiO2 (a-SiO2) substrate. First, the theoretical methodologies and concepts used in our simulations are presented. In particular, two concepts are described in detail as they are necessary for the understanding of the subsequent chapters. The first is the linear response Green-Kubo (GK) theory of thermal boundary conductance (TBC), which we develop in this thesis, and the second is the spectral energy density method, which we use to directly compute the phonon lifetimes and thermal transport coefficients. After we set the conceptual foundations, the TBC of the CNT-SiO2 interface is computed using non- equilibrium molecular dynamics (MD) simulations and the new Green-Kubo method that we have developed. Its dependence on temperature, the strength of the interaction with the substrate, and tube diameter are evaluated. To gain further insight into the phonon dynamics in supported CNTs, the scattering rates are computed using the spectral energy density (SED) method. With this method, we are able to distinguish the different scattering mechanisms (boundary and CNT-substrate phonon-phonon) and rates. The phonon lifetimes in supported CNTs are found to be reduced by contact with the substrate and we use that lifetime reduction to determine the change in CNT thermal conductivity. Next, we examine thermal transport in graphene supported on SiO2. The phonon contribution to the TBC of the graphene-SiO2 interface is computed from MD simulations and found to agree well with experimentally measured values. We derive the theory of remote phonon scattering of graphene electrons and compute the heat transfer coefficient dependence on doping level and temperature. The thermal boundary conductance from remote phonon scattering is found to be an order of magnitude smaller than that of the phonon contribution. The in-plane thermal conductivity of supported graphene is calculated from MD simulations. The experimentally measured order of magnitude reduction in thermal conductivity is reproduced in our simulations. We show that this reduction is due to the damping of the flexural (ZA) modes. By varying the interaction between graphene and the substrate, the ZA modes hybridize with the substrate Rayleigh modes and the dispersion of the hybridized modes is found to linearize in the strong coupling limit, leading to an increased thermal conductance in the composite structure.

All-regime Lagrangian-Remap numerical schemes for the gas dynamics equations. Applications to the large friction and low Mach coefficients

In this talk, we propose an all regime Lagrange-Projection like numerical scheme for the gas dynamics equations. By all regime, we mean that the numerical scheme is able to compute accurate approximate solutions with an under-resolved discretization with respect to the Mach number M, i.e. such that the ratio between the Mach number M and the mesh size or the time step is small with respect to 1. The key idea is to decouple acoustic and transport phenomenon and then alter the numerical flux in the acoustic approximation to obtain a uniform truncation error in term of M. This modified scheme is conservative and endowed with good stability properties with respect to the positivity of the density and the internal energy. A discrete entropy inequality under a condition on the modification is obtained thanks to a reinterpretation of the modified scheme in the Harten Lax and van Leer formalism. A natural extension to multi-dimensional problems discretized over unstructured mesh is proposed. Then a simple and efficient semi implicit scheme is also proposed. The resulting scheme is stable under a CFL condition driven by the (slow) material waves and not by the (fast) acoustic waves and so verifies the all regime property. Numerical evidences are proposed and show the ability of the scheme to deal with tests where the flow regime may vary from low to high Mach values.

Analysis of the four wave mixing effect (FWM) in a dispersion decreasing fiber (DDF) for a WDM system

We did a numerical investigation of the propagation of short light pulses in the region of 1.55 mu m and the conversion efficiency (CE) for the four wave mixing generation (FWM) of ordinary and dispersion decreasing fibers for use in wavelength division multiplexing (WDM) systems, Our simulations studies three different profiles, linear, hyperbolic. and constant, One conclude that for all the profiles there is decrease of the conversion efficiency with the increase in the channel separation. The hyperbolic profile present a higher efficiency of around 1000 above in magnitude compared with the others profiles at 0.2 nm of channel separation. We calculate the conversion efficiency versus the fiber length for the three profiles. The conversion efficiency for the hyperbolic profile is higher when compared to the constant and linear profiles. The other interesting point of the hyperbolic profile is that the increase of the CE in the beginning of the fiber does not show my oscillation in the CE value (log eta), which was observed for the constant and linear profiles. For all the profiles there is an increase of the conversion efficiency with the increase of the pump power. The compression factor C-i for the generated FWM signal at omega(3) was measured along the DDF's and the constant profile fibers. One can conclude that with the use of decreasing dispersion profile (DDF) fibers one can have a control of the (CE) conversion efficiency and the compression factor of the four wave mixing (FWM) generation in WDM systems. (c) 2005 Elsevier B.V. All rights reserved.

Motion compensation with minimal residue dispersion matching criteria

Dissertação (mestrado)—Universidade de Brasília, Faculdade de Tecnoloigia, 2016.

The Becker-Döring equations with exponentially size-dependent rate coefficients

This paper is concerned with an analysis of the Becker-Döring equations which lie at the heart of a number of descriptions of non-equilibrium phase transitions and related complex dynamical processes. The Becker-Döring theory describes growth and fragmentation in terms of stepwise addition or removal of single particles to or from clusters of similar particles and has been applied to a wide range of problems of physicochemical and biological interest within recent years. Here we consider the case where the aggregation and fragmentation rates depend exponentially on cluster size. These choices of rate coefficients at least qualitatively correspond to physically realistic molecular clustering scenarios such as occur in, for example, simulations of simple fluids. New similarity solutions for the constant monomer Becker-Döring system are identified, and shown to be generic in the case of aggregation and fragmentation rates that depend exponentially on cluster size.

Fluorinated surfactants in solution: Diffusion coefficients of fluorinated alcohols in water

Intra-diffusion coefficients of three fluorinated alcohols, 2,2,3,3,3-pentafluoropropan-1-ol (PFP), 2,2,3,3,4,4,4-heptafluorobutan-1-ol (HFB) and 2,2,3,3,4,4,5,5,5-nonafluoropentan-1-ol (NFP) in water have been measured by the PFG–NMR spin-echo technique as a function of temperature and composition, focusing on the alcohol dilute region. For comparison, intra-diffusion coefficients of 2,2,2- trifluoroethanol (TFE) and HFB have also been measured in heavy water using the same method and conditions. As far as we know, these are the first experimental measurements of this property for these binary systems. Intra-diffusion coefficients for NFP in water and for TFE and HFB in heavy water have also been obtained by molecular dynamics simulation, complementing those for TFE, PFP and HFB reported in a previous work. The agreement between experimental and simulated results for PFP, HFB and NFP in water is reasonable, although presenting higher deviations than for the TFE/water system. From the dependence of the intra-diffusion coefficients on temperature, diffusion activation energies were estimated for all the solutes in water and heavy water.

Rheological Properties of a Sodium Bentonite and Lactose Aqueous Dispersion

Aim of this paper is show the viscosity measure of a sodium bentonite-water-lactose mixture and your rheological behaviour. This analysis showed the formation of tridimensional structure type and formation of stratified silicate/lactose, this occurred due to different concentrations of organic products into mixture and due to a difference of rotation during viscosity measument. Formation of networks is a consequence of the attraction between the silicate layers in water-lactose mixture. In the present work aqueous solutions of lactose with concentration of 7%, 5%, 3%, 1% and 0% (wt %) were used.

Dynamics and Heterogeneity of Pb(II) Binding by SiO2 Nanoparticles in an Aqueous Dispersion

Pb(II) binding by SiO2 nanoparticles in an aqueous dispersion was investigated under conditions where the concentrations of Pb2+ ions and nanoparticles are of similar magnitude. Conditional stability constants (log K) obtained at different values of pH and ionic strength varied from 4.4 at pH 5.5 and I = 0.1 M to 6.4 at pH 6.5 and I = 0.0015 M. In the range of metal to nanoparticle ratios from 1.6 to 0.3, log K strongly increases, which is shown to be due to heterogeneity in Pb(II) binding. For an ionic strength of 0.1 M the Pb2+/SiO2 nanoparticle system is labile, whereas for lower ionic strengths there is loss of lability with increasing pH and decreasing ionic strength. Theoretical calculations on the basis of Eigen-type complex formation kinetics seem to support the loss of lability. This is related to the nanoparticulate nature of the system, where complexation rate constants become increasingly diffusion controlled. The ion binding heterogeneity and chemodynamics of oxidic nanoparticles clearly need further detailed research.

Assessing the activity coefficients of water in cholinium-based ionic liquids: experimental measurements and COSMO-RS modeling

The vapor liquid-equilibrium of water + ionic liquids is relevant for a wide range of applications of these compounds. It is usually measured by ebulliometric techniques, but these are time consuming and expensive. In this work it is shown that the activity coefficients of water in a series of cholinium-based ionic liquids can be reliably and quickly estimated at 298.15K using a humidity meter instrument. The cholinium based ionic liquids were chosen to test this experimental methodology since data for water activities of quaternary ammonium salts are available in the literature allowing the validation of the proposed technique. The COSMO-RS method provides a reliable description of the data and was also used to understand the molecular interactions occurring on these binary systems. The estimated excess enthalpies indicate that hydrogen bonding between water and ionic liquid anion is the dominant interaction that governs the behavior of water and cholinium-based ionic liquids systems, while the electrostatic-misfit and van der Walls forces have a minor contribution to the total excess enthalpies.

Modélisation hydro-sédimentaire de la rade de Toulon : application à l’étude de la dispersion de radionucléides

L’objectif de ce travail est mettre en place un modèle hydro-sédimentaire de dispersion des radionucléides, comme outil d’aide à la décision suite à une hypothétique contamination marine accidentelle. Appliqué à la Rade de Toulon, il utilise le modèle MARS-3D pour l’étude hydrodynamique et le module MIXSED pour la dynamique des différentes classes de sédiments. Les cas d’application présentés concernent le césium et le plutonium, qui ont des affinités très différentes avec les sédiments fins. Les simulations de rejet, effectuées en fond de baie, montrent une variabilité saisonnière marquée, due aux conditions météorologiques typiques, mais aussi aux conditions de débit et de charge sédimentaire associée du fleuve Las, qui se jette dans la Rade. Ainsi, les simulations par forts vents, qui favorisent la circulation et les échanges de masses d’eau, présentent une diminution rapide des activités dissoutes. Par ailleurs, le Vent d’Est hivernal, qui engendre d’importantes crues du Las et un apport de sédiments à la Rade, favorise le piégeage des radionucléides dans les sédiments de fond, où les deuxtiers de la contamination initiale sont piégés après une simulation de deux mois.

Measurements of activity coefficients at infinite dilution of organic solutes and water on polar imidazolium-based ionic liquids

The activity coefficients at infinite dilution, gamma(infinity)(13), of 55 organic solutes and water in three ionic liquids with the common cation 1-butyl-3-methylimidazolium and the polar anions Cl--,Cl- [CH3SO3](-) and [(CH3)(2)PO4](-), were determined by (gas + liquid) chromatography at four temperatures in the range (358.15 to 388.15) K for alcohols and water, and T = (398.15 to 428.15) K for the other organic solutes including alkanes, cycloalkanes, alkenes, cycloalkenes, alkynes, ketones, ethers, cyclic ethers, aromatic hydrocarbons, esters, butyraldehyde, acetonitrile, pyridine, 1-nitropropane and thiophene. From the experimental gamma(infinity)(13) values, the partial molar excess Gibbs free energy, (G) over bar (E infinity)(m), enthalpy (H) over bar (E infinity)(m), and entropy (S) over bar (E infinity)(m), at infinite dilution, were estimated in order to provide more information about the interactions between the solutes and the ILs. Moreover, densities were measured and (gas + liquid) partition coefficients (KL) calculated. Selectivities at infinite dilution for some separation problems such as octane/benzene, cyclohexane/benzene and cyclohexane/thiophene were calculated using the measured gamma(infinity)(13), and compared with literature values for N-methyl-2-pyrrolidinone (NMP), sulfolane, and other ionic liquids with a common cation or anion of the ILs here studied. From the obtained infinite dilution selectivities and capacities, it can be concluded that the ILs studied may replace conventional entrainers applied for the separation processes of aliphatic/aromatic hydrocarbons.