published_or_final_version / Electrical and Electronic Engineering / Master / Master of Philosophy

No description available.

Thesis (M.S.)--Massachusetts Institute of Technology, Dept. of Civil Engineering, 1980. / MICROFICHE COPY AVAILABLE IN ARCHIVES AND ENGINEERING. / Bibliography: leaves 150-155. / by Paul Christopher Damian Milly. / M.S.

Cloud-wind interactions play an important role in long-lived multiphase flows in various galaxy-related contexts (e.g., galactic fountains and winds, cosmological cold-mode accretion, or multiphase tails of satellites). These interactions occur when a volume-filling hot phase, the wind, moves relative to a cool pressure-confined body of gas, the cloud. The conditions necessary for clouds to survive the destructive effects of mixing and become entrained within the wind (i.e. for the relative velocity to be removed), has been a long-standing problem. This problem has received particular attention in the context of galactic winds: cloud entrainment is expected to play a critical role in explaining observed multiphase structure in these outflows. This thesis investigates a mechanism for facilitating cloud survival in the context of rapid cooling, which we hereafter term TRML (turbulent radiative mixing layer) entrainment. Our investigation leverages numerical (magneto)hydrodynamic ENZO-E simulations of a cool (≲10⁴𝐊) clouds that encounter a hot (≳10⁶𝐊), supersonic winds. We begin by introducing a simple entropy-based formalism to characterize the role of mixing in cloud-wind interactions, and demonstrate example applications using simulations. Under this formalism, the high-dimensional description of the interaction's state at a given time is simplified to the joint distribution of mass over pressure (𝑃 ) and entropy (𝐾=𝑃𝞀^-𝜸). As a result, this approach provides a way for (empirically and analytically) quantifying the impact of different initial conditions and sets of physics on the interaction's evolution. We find that mixing predominantly alters the distribution along the 𝐾 direction and illustrate how the formalism can be used to model mixing and cooling for fluid elements originating in the cloud. We further confirm and generalize a previously suggested survival criterion for clouds undergoing TRML entrainment, and demonstrate that the shape of the cooling curve, particularly at the low temperature end, can play an important role in controlling condensation. Moreover, we discuss the capacity of our approach to generalize such a criterion to apply to additional sets of physics, and to build intuition for the impact of subtle higher order effects not directly addressed by the criterion. Despite the fact that the competition the between turbulent mixing and radiative cooling dictate the outcome of the cloud-wind interaction (as well as many observable properties), turbulence in these interactions remains poorly understood. Thus, we next investigate the turbulence that arises for clouds undergoing TRML entrainment. To obtain robust results, we employ multiple metrics to characterize the turbulent velocity, 𝝂_turb. We find four primary results. First, 𝝂_turb manifests clear temperature dependence. Initially, 𝝂_turb roughly matches the scaling of sound speed on temperature. In gas hotter than the temperature where cooling peaks, this dependence weakens with time until 𝝂_turb is constant. Second, the relative velocity between the cloud and wind initially drives rapid growth of 𝝂_turb. As it drops (from entrainment), 𝝂_turb starts to decay before it stabilizes at roughly half its maximum. At late times cooling flows appear to support turbulence. Third, the magnitude of 𝝂_turb scales with the ratio between the hot phase sound crossing time and the minimum cooling time. Finally, we find tentative evidence for a length-scale associated with resolving turbulence. Under-resolving this scale may cause violent shattering and affect the cloud's large-scale morphological properties. Finally, we propose a new criterion for clouds to survive interactions with the wind in the via TRML entrainment, and validate it with simulations. Properties of TRML entrainment are generally understood to be controlled by ratio between the relevant dynamical and cooling timescales 𝝉_dyn / 𝝉_cool. Previously proposed survival criteria disagree about the size of the smallest surviving cloud by factors of up to ∼100. These criteria primarily differ in their choice of 𝝉_{\rm cool}$; perplexingly, the choices most consistent with the well-modeled micro-scale physics observed in shear-layer studies are associated with less-accurate criteria. We present a new criterion which agrees with previous fitting formulae but is based on a set of simple physical principles. Whereas prior criteria link 𝝉_dyn with the cloud destruction timescale, our new criterion links it to the characteristic cloud-crossing timescale of a hot-phase fluid element. This choice leads to scaling relations that are more physically consistent with shear-layer studies. Additionally, we illustrate that discrepancies among previous criteria primarily emerged due to the choices of simulation conditions, rather than commonly-cited differences in the definition of cloud destruction.

[Truncated abstract] As gravitational wave detectors approach sensitivities that will allow observations to become routine, astrophysics lies on the cusp of an exciting new era. Potential sources will include transients such as merging neutron stars and black holes, supernova explosions or the engines that power gamma-ray bursts. This thesis will be devoted to the astrophysical gravitational wave background signal produced by cosmological populations of such transient signals. Particular attention will be devoted to the observation-time dependence imposed on the individual sources that accumulate to produce a gravitational wave background signal. The ultimate aim is to determine what information is encoded in the temporal evolution of such a signal. To lay the foundations for further investigation, the stochastic gravitational wave background signal from neutron star birth throughout the Universe has been calculated. In view of the uncertainties in both the single-source emissions and source rate histories, several models of each are employed. The results show that that the resulting signals are only weakly dependent on the source-rate evolution model and that prominent features in the single-source spectra can be related to the background spectra. In comparison with previous studies, the use of relativistic single-source gravitational wave waveforms rather than Newtonian models and a more slowly evolving source-rate density results in a 1 { 2 order of magnitude reduction in signal. ... A comparison with the more commonly used brightness distribution of events shows that when applying both methods to a data stream containing a background of Gaussian distributed false alarms, the brightness distribution yielded lower standard errors, but was biased by the false alarms. In comparison, a fitting procedure based on the time evolution of events was less prone to errors resulting from false alarms, but as fewer events contributed to the data, had a lower resolution. In further support of the time dependent signature of transient events, an alternative technique is fiapplied to the same source population. In this case, the local rate density is probed by measuring the statistical compatibility of the filtered data against synthetic time dependent data. Although this method is not as compact as the fitting procedure, the rate estimates are compatible. To further investigate how the observation time dependence of transient populations can be used to constrain global parameters, the method is applied to Swift long gamma-ray burst data. By considering a distribution in peak °ux rather than a gravitational wave amplitude, gamma-ray bursts can be considered as a surrogate for resolved gravitational wave transients. For this application a peak °ux{observation time relation is described that takes the form of a power law that is invariant to the luminosity distribution of the sources. Additionally, the method is enhanced by invoking time reversal invariance and the temporal cosmological principle. Results are presented to show that the peak °ux{observation time relation is in good agreement with recent estimates of source parameters. Additionally, to show that the intrinsic time dependence allows the method to be used as a predictive tool, projections are made to determine the upper limits in peak °ux of future gamma-ray burst detections for Swift.

The atmospheres of exoplanets are being characterised in increasing detail by observational facilities and will be examined with even greater clarity with upcoming space based missions such as the James Webb Space Telescope (JWST) and the Wide Field InfraRed Survey Telescope (WFIRST). A major component of exoplanet atmospheres is the presence of cloud particles which produce characteristic observational signatures in transit spectra and influence the geometric albedo of exoplanets. Despite a decade of observational evidence, the formation, dynamics and radiative-transport of exoplanet atmospheric cloud particles remains an open question in the exoplanet community. In this thesis, we investigate the kinetic chemistry of cloud formation in hot Jupiter exoplanets, their effect on the atmospheric dynamics and observable properties. We use a static 1D cloud formation code to investigate the cloud formation properties of the hot Jupiter HD 189733b. We couple a time-dependent kinetic cloud formation to a 3D radiative-hydrodynamic simulation of the atmosphere of HD 189733b and investigate the dynamical properties of cloud particles in the atmosphere. We develop a 3D multiple-scattering Monte Carlo radiative-transfer code to post-process the results of the cloudy HD 189733b RHD simulation and compare the results to observational results. We find that the cloud structures of the hot Jupiter HD 189733b are likely to be highly inhomogeneous, with differences in cloud particle sizes, number density and composition with longitude, latitude and depth. Cloud structures are most divergent between the dayside and nightside faces of the planet due to the instability of silicate materials on the hotter dayside. We find that the HD 189733b simulation in post-processing is consistent with geometric albedo observations of the planet. Due to the scattering properties of the cloud particles we predict that HD 189733b will be brighter in the upcoming space missions CHaracterising ExOPlanet Satellite (CHEOPS) bandpass compared to the Transiting Exoplanet Space Survey (TESS) bandpass.

This thesis is primarily focused on problems in extremal combinatorics, although we will also consider some questions of analytic and algorithmic nature. The d-dimensional hypercube is the graph with vertex set {0,1}^d where two vertices are adjacent if they differ in exactly one coordinate. In Chapter 2 we obtain an upper bound on the 'saturation number' of Q_m in Q_d. Specifically, we show that for m ≥ 2 fixed and d large there exists a subgraph G of Q_d of bounded average degree such that G does not contain a copy of Q_m but, for every G' such that G ⊊ G' ⊆ Q_d, the graph G' contains a copy of Q_m. This result answers a question of Johnson and Pinto and is best possible up to a factor of O(m). In Chapter 3, we show that there exists ε > 0 such that for all k and for n sufficiently large there is a collection of at most 2^{(1-ε)k} subsets of [n] which does not contain a chain of length k+1 under inclusion and is maximal subject to this property. This disproves a conjecture of Gerbner, Keszegh, Lemons, Palmer, Pálvölgyi and Patkós. We also prove that there exists a constant c ∈ (0,1) such that the smallest such collection is of cardinality 2^{(1+o(1))^ck} for all k. In Chapter 4, we obtain an exact expression for the 'weak saturation number' of Q_m in Q_d. That is, we determine the minimum number of edges in a spanning subgraph G of Q_d such that the edges of E(Q_d)\E(G) can be added to G, one edge at a time, such that each new edge completes a copy of Q_m. This answers another question of Johnson and Pinto. We also obtain a more general result for the weak saturation of 'axis aligned' copies of a multidimensional grid in a larger grid. In the r-neighbour bootstrap process, one begins with a set A₀ of 'infected' vertices in a graph G and, at each step, a 'healthy' vertex becomes infected if it has at least r infected neighbours. If every vertex of G is eventually infected, then we say that A₀ percolates. In Chapter 5, we apply ideas from weak saturation to prove that, for fixed r ≥ 2, every percolating set in Q_d has cardinality at least (1+o(1))(d choose r-1)/r. This confirms a conjecture of Balogh and Bollobás and is asymptotically best possible. In addition, we determine the minimum cardinality exactly in the case r=3 (the minimum cardinality in the case r=2 was already known). In Chapter 6, we provide a framework for proving lower bounds on the number of comparable pairs in a subset S of a partially ordered set (poset) of prescribed size. We apply this framework to obtain an explicit bound of this type for the poset 𝒱(q,n) consisting of all subspaces of 𝔽_qⁿordered by inclusion which is best possible when S is not too large. In Chapter 7, we apply the result from Chapter 6 along with the recently developed 'container method,' to obtain an upper bound on the number of antichains in 𝒱(q,n) and a bound on the size of the largest antichain in a p-random subset of 𝒱(q,n) which holds with high probability for p in a certain range. In Chapter 8, we construct a 'finitely forcible graphon' W for which there exists a sequence (ε_i)^∞_i=1 tending to zero such that, for all i ≥ 1, every weak ε_i-regular partition of W has at least exp(ε_i^-2/2^{5log∗ε_i^-2}) parts. This result shows that the structure of a finitely forcible graphon can be much more complex than was anticipated in a paper of Lovász and Szegedy. For positive integers p,q with p/q ❘≥ 2, a circular (p,q)-colouring of a graph G is a mapping V(G) → ℤ_p such that any two adjacent vertices are mapped to elements of ℤ_p at distance at least q from one another. The reconfiguration problem for circular colourings asks, given two (p,q)-colourings f and g of G, is it possible to transform f into g by recolouring one vertex at a time so that every intermediate mapping is a p,q-colouring? In Chapter 9, we show that this question can be answered in polynomial time for 2 ≤ p/q < 4 and is PSPACE-complete for p/q ≥ 4.

Le bicarbonate de soude raffiné, produit industriellement par la société Solvay, est fabriqué dans des colonnes à bulles de grande taille, appelées les colonnes BIR.<p>Dans ces colonnes, une phase gazeuse contenant un mélange d’air et dioxyde de carbone (CO2) est dispersée sous forme de bulles dans une solution aqueuse de carbonate et de bicarbonate de sodium (respectivement Na2CO3 et NaHCO3). Cette dispersion donne lieu à un transfert de CO2 des bulles vers la phase liquide. Au sein des colonnes, la phase gazeuse se répartit dans deux populations de bulles :des petites bulles (diamètre de quelques mm) et des grandes bulles (diamètre de quelques cm). Le transfert bulle-liquide de CO2 est couplé à des réactions chimiques prenant place en phase liquide, qui conduisent à la conversion du Na2CO3 en NaHCO3. Une fois la concentration de saturation dépassée le NaHCO3 précipite sous forme de cristaux et un mélange liquide-solide est recueilli à la sortie de ces colonnes.<p>Ce travail, réalisé en collaboration avec la société Solvay, porte sur l’étude et la modélisation mathématique des phénomènes de transfert de matière entre phases, couplés à des réactions chimiques, prenant place au sein d’une colonne BIR. L’association d’études sur des colonnes à bulles à l’échelle industrielle ou réduite (pilote) et d’études plus fondamentales sur des dispositifs de laboratoire permet de développer une meilleure compréhension du fonctionnement des colonnes BIR et d’en construire un modèle mathématique détaillé.<p>L’objectif appliqué de ce travail est la mise au point d’un modèle mathématique complet et opérationnel d’une colonne BIR. Cet objectif est supporté par trois blocs de travail, dans lesquels différents outils sont développés et exploités.<p><p>Le premier bloc est consacré à la modélisation mathématique du transfert bulle-liquide de CO2 dans une solution aqueuse de NaHCO3 et de Na2CO3. Ce transfert est couplé à des réactions chimiques en phase liquide qui influencent sa vitesse. Dans un premier temps, des modèles sont développés selon des approches unidimensionnelles classiquement rencontrées dans la littérature. Ces approches passent par une idéalisation de l’écoulement du liquide autour des bulles. Une expression simplifiée de la vitesse du transfert bulle-liquide de CO2, est également développée et validée pour le modèle de colonne BIR.<p>Dans un second temps, une modélisation complète des phénomènes de transport (convection et diffusion), couplés à des réactions chimiques, est réalisée en suivant une approche bidimensionnelle axisymétrique. L’influence de la vitesse de réactions sur la vitesse de transfert est étudiée et les résultats des deux approches sont également comparés.<p><p>Le deuxième bloc est consacré à l’étude expérimentale du transfert gaz-liquide de CO2 dans des solutions aqueuses de NaHCO3 et de Na2CO3. A cette fin, un dispositif expérimental est développé et présenté. Du CO2 est mis en contact avec des solutions aqueuses de NaHCO3 et de Na2CO3 dans une cellule transparente. Les phénomènes provoqués en phase liquide par le transfert de CO2 sont observés à l’aide d’un interféromètre de Mach-Zehnder.<p>Les résultats expérimentaux sont comparés à des résultats de simulation obtenus avec un des modèles unidimensionnels développés dans le premier bloc. De cette comparaison, il apparaît qu’une mauvaise estimation de la valeur de certains paramètres physico-chimiques apparaissant dans les équations de ce modèle conduit à des écarts significatifs entre les grandeurs observées expérimentalement et les grandeurs estimées par simulation des équations du modèle.<p>C’est pourquoi une méthode d’estimation paramétrique est également développée afin d’identifier les valeurs numériques de ces paramètres physico-chimiques sur base des résultats expérimentaux. Ces dernières sont également discutées.<p><p>Dans le troisième bloc, nous apportons une contribution à l’étude des cinétiques de précipitation du NaHCO3 dans un cristallisoir à cuve agitée. Cette partie du travail est réalisée en collaboration avec Vanessa Gutierrez (du service Matières et Matériaux de l’ULB).<p>Nous contribuons à cette étude par le développement de trois outils :une table de calcul Excel permettant de synthétiser les résultats expérimentaux, un ensemble de simulations de l’écoulement au sein du cristallisoir par mécanique des fluides numérique et une nouvelle méthode d’extraction des cinétiques de précipitation du NaHCO3 à partir des résultats expérimentaux. Ces trois outils sont également utilisés de façon combinée pour estimer les influences de la fraction massique de solide et de l’agitation sur la cinétique de germination secondaire du NaHCO3.<p><p>Enfin, la synthèse de l’ensemble des résultats de ces études est réalisée. Le résultat final est le développement d’un modèle mathématique complet et opérationnel des colonnes BIR. Ce modèle est développé en suivant l’approche de modélisation en compartiments, développée au cours du travail de Benoît Haut. Ce modèle synthétise les trois blocs d’études réalisées dans ce travail, ainsi que les travaux d’Aurélie Larcy (du service Transferts, Interfaces et Procédés de l’ULB) et de Vanessa Gutierrez. Les équations modélisant les différents phénomènes sont présentées, ainsi que la méthode utilisée pour résoudre ces équations. Des simulations des équations du modèle sont réalisées et discutées. Les résultats de simulation sont également comparés à des mesures effectuées sur une colonne BIR. Un accord raisonnable est observé.<p>A l’issue de ce travail, nous disposons donc d’un modèle opérationnel de colonne BIR. Bien que ce modèle doive encore être optimisé et validé, il peut déjà être utilisé pour étudier l’effet des caractéristiques géométriques des colonnes BIR et des conditions appliquées à ces colonnes sur le comportement des simulations des équations du modèle et pour identifier des tendances.<p>//<p>The refined sodium bicarbonate is produced by the Solvay company using large size bubble columns, called the BIR columns.<p>In these columns, a gaseous phase containing an air-carbon dioxyde mixture (CO2) is dispersed under the form of bubbles in an aqueous solution of sodium carbonate and sodium bicarbonate (Na2CO3 and NaHCO3, respectively). This dispersion leads to a CO2 transfer from the bubbles to the liquid phase. Inside these columns, the gaseous phase is distributed in two bubbles populations :small bubbles (a few mm of diameter) and large bubbles (a few cm of diameter).<p>The bubble-liquid CO2 transfer is coupled with chemical reactions taking places in the liquid phase that leads to the conversion of Na2CO3 to NaHCO3. When the solution is supersaturated in NaHCO3, the NaHCO3 precipitates under the form of crystals and a liquid-solid mixture is extracted at the outlet of the BIR columns.<p>This work, realized in collaboration with Solvay, aims to study and to model mathematically the mass transport phenomena between the phases, coupled with chemical reactions, taking places inside a BIR column. Study of bubble columns at the industrial and the pilot scale is combined to a more fundamental study at laboratory scale to improve the understanding of the BIR columns functioning and to develop a detailed mathematical modeling.<p>The applied objective of this work is to develop a complete and operational mathematical modeling of a BIR column. This objective is supported by three blocks of work. In each block, several tools are developed and used.<p><p>The first block is devoted to the mathematical modeling of the bubble-liquid CO2 transfer in an NaHCO3 and Na2CO3 aqueous solution. This transfer is coupled with chemical reactions in liquid phase, which affect the transfer rate.<p>In a first time, mathematical models are developed following the classical one-dimensional approaches of the literature. These approaches idealize the liquid flow around the bubbles. A simplified expression of the bubble-liquid CO2 transfer rate is equally developed and validated for the BIR column model.<p>In a second time, a complete modeling of the transport phenomena (convection and diffusion) coupled with chemical reactions is developed, following an axisymmetrical twodimensional approach. The chemical reaction rate influence on the bubble-liquid transfer rate is studied and the results of the two approaches are then compared.<p><p>The second block is devoted to the experimental study of the gas-liquid CO2 transfer to NaHCO3 and Na2CO3 aqueous solutions. An experimental set-up is developed and presented. CO2 is put in contact with NaHCO3 and Na2CO3 aqueous solutions in a transparent cell. The phenomena induced in liquid phase by the CO2 transfer are observed using a Mach-Zehnder interferometer.<p>The experimental results are compared to simulation results that are obtained using one of the one-dimensional model developed in the first block. From this comparison, it appears that a wrong estimation of some physico-chemical parameter values leads to significative differences between the experimentally observed quantities and those estimated by simulation of the model equations. Therefore, a parametric estimation method is developed in order to estimate those parameters numerical values from the experimental results. The found values are then discussed.<p><p>In the third block is presented a contribution to the NaHCO3 precipitation kinetic study in a stirred-tank crystallizer. This part of the work is realized in collaboration with Vanessa Gutierrez (Chemicals and Materials Department of ULB).<p>Three tools are developed :tables in Excel sheet to synthetize the experimental results, a set of simulations of the flow inside the crystallizer by Computational Fluid Dynamic (CFD) and a new method to extract the NaHCO3 precipitation kinetics from the experimental measurements. These three tools are combined to estimate the influences of the solid mass fraction and the flow on the NaHCO3 secondary nucleation rate.<p><p>Finally, the synthesis of all these results is realized. The final result is the development of a complete and operational mathematical model of BIR columns. This model is developed following the compartmental modeling approach, developed in the PhD thesis of Benoît Haut. This model synthetizes the three block of study realized in this work and the studies of Aurélie Larcy (Transfers, Interfaces and Processes Department of ULB) and those of Vanessa Gutierrez. The equations modeling the phenomena taking place in a BIR column are presented as the used method to solve these equations. The equations of the model are simulated and the results are discussed. The results are equally compared to experimental measurement realized on a BIR column. A reasonable agreement is observed.<p>At the end of this work, an operational model of a BIR column is thus developed. Although this model have to be optimized and validated, it can already be used to study the influences of the geometrical characteristics of the BIR columns and of the conditions applied to these columns on the behaviour of the model equation simulations and to identity tendencies. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished