Numerical simulation of Ginzburg-Landau-Langevin equations

This work is concerned with non-equilibrium phenomena, with focus on the numerical simulation of the relaxation of non-conserved order parameters described by stochastic kinetic equations known as Ginzburg-Landau-Langevin (GLL) equations. We propose methods for solving numerically these type of equations, with additive and multiplicative noises. Illustrative applications of the methods are presented for different GLL equations, with emphasis on equations incorporating memory effects.

Chemical Nonequilibrium Effects on Flow Field for Reusable Launch Vehicles

High temperature chemical non-equilibrium phenomena have a great effect on the flow field around a reentry vehicle. A set of three dimensional Navier-Stokes equations have been solved by implicit finite volume NND scheme. Both ideal gas viscous flow and chemical non-equilibrium flow are calculated for a spherical-cone at a small angle of attack. The results of the two flows have been compared and the effect of chemical non-equilibrium has been analyzed. The effect of wall material's properties, such as catalysis and radiation were studied. The results are in good agreement with the referenced paper.

Light-matter Interactions in Semiconductors and Metals: From Nitride Optoelectronics to Quantum Plasmonics

This thesis puts forth a theory-directed approach coupled with spectroscopy aimed at the discovery and understanding of light-matter interactions in semiconductors and metals.

The first part of the thesis presents the discovery and development of Zn-IV nitride materials.The commercial prominence in the optoelectronics industry of tunable semiconductor alloy materials based on nitride semiconductor devices, specifically InGaN, motivates the search for earth-abundant alternatives for use in efficient, high-quality optoelectronic devices. II-IV-N2 compounds, which are closely related to the wurtzite-structured III-N semiconductors, have similar electronic and optical properties to InGaN namely direct band gaps, high quantum efficiencies and large optical absorption coefficients. The choice of different group II and group IV elements provides chemical diversity that can be exploited to tune the structural and electronic properties through the series of alloys. The first theoretical and experimental investigation of the ZnSnxGe1−xN2 series as a replacement for III-nitrides is discussed here.

The second half of the thesis shows ab−initio calculations for surface plasmons and plasmonic hot carrier dynamics. Surface plasmons, electromagnetic modes confined to the surface of a conductor-dielectric interface, have sparked renewed interest because of their quantum nature and their broad range of applications. The decay of surface plasmons is usually a detriment in the field of plasmonics, but the possibility to capture the energy normally lost to heat would open new opportunities in photon sensors, energy conversion devices and switching. A theoretical understanding of plasmon-driven hot carrier generation and relaxation dynamics in the ultrafast regime is presented here. Additionally calculations for plasmon-mediated upconversion as well as an energy-dependent transport model for these non-equilibrium carriers are shown.

Finally, this thesis gives an outlook on the potential of non-equilibrium phenomena in metals and semiconductors for future light-based technologies.

Efeitos dinâmicos de programas de estabilização

Com o emprego do modelo de dois setores de acumulação ótima de capital em economia aberta, determina-se o impacto sobre a trajetória do câmbio, dos salários, do investimento, da poupança e, portanto, da dívida externa e do estoque de capital, de uma elevação permanente e não antecipada da produtividade da economia. Em geral, após um choque positivo permanente de produtividade há redução da poupança, piora do balanço de pagamentos em transações correntes e valorização do câmbio. Todos fenômenos que do ponto de vista do modelo são de equilíbrio intertemporal, conseqüência da elevação da renda permanente e do excesso de demanda por bens domésticos que sucede o ganho de produtividade. Supondo que os programas de estabilização elevaram a produtividade da economia é possível com a estrutura analítica construída racionalizar qualitativamente os fenômenos observados após estes planos.

Ajustamento de uma economia após uma elevação da produtividade

Employing the two sector model of capital accumulation in an open economy, the impact on the path of the following variables: exchange rate, wages, investment, saving, and consequently externaI debt and capital stock afier a permanent and non expected elevation of the economy productivity is determinated. Afier this positive shock, saving rate decreases, current transaction deteriorates and the exchange rate appreciates. Those are equilibrium phenomena from 3D intertemporaI point of view due to the permanent income raise and to the domestic good excess demand that follows the productivity increase. Assuming that the stabilization programa augment the economy productivity, the model could rationalize qualitatively the stylized facts witnessed after those programa.

Modelling, diagnostics and experimental analysis of plasma assisted processes for material treatment

This work presents results from experimental investigations of several different atmospheric pressure plasmas applications, such as Metal Inert Gas (MIG) welding and Plasma Arc Cutting (PAC) and Welding (PAW) sources, as well as Inductively Coupled Plasma (ICP) torches. The main diagnostic tool that has been used is High Speed Imaging (HSI), often assisted by Schlieren imaging to analyse non-visible phenomena. Furthermore, starting from thermo-fluid-dynamic models developed by the University of Bologna group, such plasma processes have been studied also with new advanced models, focusing for instance on the interaction between a melting metal wire and a plasma, or considering non-equilibrium phenomena for diagnostics of plasma arcs. Additionally, the experimental diagnostic tools that have been developed for industrial thermal plasmas have been used also for the characterization of innovative low temperature atmospheric pressure non equilibrium plasmas, such as dielectric barrier discharges (DBD) and Plasma Jets. These sources are controlled by few kV voltage pulses with pulse rise time of few nanoseconds to avoid the formation of a plasma arc, with interesting applications in surface functionalization of thermosensitive materials. In order to investigate also bio-medical applications of thermal plasma, a self-developed quenching device has been connected to an ICP torch. Such device has allowed inactivation of several kinds of bacteria spread on petri dishes, by keeping the substrate temperature lower than 40 degrees, which is a strict requirement in order to allow the treatment of living tissues.

Optische Mikroskopie an elektrostatisch stabilisierten kolloidalen Suspensionen unter dem Einfluss geladener Wände - ausgewählte Gleichgewichts- und Nichtgleichgewichtsphänomene

Das Verhalten kolloidaler Suspensionen unter räumlich beschränkter Geometrie ist von großer Bedeutung für die statistische Physik wie auch für die Technologie. Von speziellem Interesse sind Modellsysteme geladener kolloidaler Sphären aufgrund ihrer langreichweitigen und veränderbaren Wechselwirkungen. In dieser Arbeit wurde ein experimenteller Aufbau für die optische mikroskopische Untersuchung solcher, zwischen ebenen Wänden beschränkter Systeme realisiert. Anhand von Piezo-Aktuatoren kann die Zellgeometrie flexibel und präzise eingestellt werden. Unter Verwendung eines Pumpkreislaufs mit einer Ionentauschersäule können kolloidale Suspensionen unter stark entsalzten Bedingungen effizient präpariert werden. Anhand dieses Aufbaus wurde zunächst das Gleichgewichtsphasendiagramm monodisperser geladener kolloidaler Sphären zwischen parallelen Wänden untersucht. Es wurden quantitative Resultate für den Grenzfall starker Entsalzung erzielt, welche mit theoretischen Grundzustandsvorhersagen übereinstimmen. In Doppellagensystemen konnte die Existenz transienter kolloidaler Moiré-Rotationsmuster demonstriert werden, welche besondere zweidimensionale Kristallstrukturen mit komplexer Basis darstellen. Es wurden ferner Nichtgleichgewichtsphänomene untersucht, welche durch Gradienten von lokal freigesetzten Elektrolyten verursacht werden. Durch hauptsächlich diffusioosmotischen Partikeltransport entlang einer einzelnen geladenen Substratoberfläche konnten die Bildung kristalliner Ordnung sowie komplexe, selbstorganisierte Bewegungszustände in einem verdünnten kolloidalen Monolagenfluid bei kleinen Reynolds-Zahlen induziert werden. Interessante Perspektiven für die zukünftige Verwendung des experimentellen Aufbaus ergeben sich aus Beobachtungen verschiedener weiterer Phänomene.

Critical point phenomena in the electrical resistivity of binary liquid systems: Cs2 + CH3CN and CS2 + CH3NO2

Electrical resistance measurements are reported on the binary liquid mixtures CS2 + CH3CN and CS2 + CH3NO2 with special reference to the critical region. Impurity conduction seems to be the dominant mechanism for charge transport. For the liquid mixture filled at the critical composition, the resistance of the system aboveT c follows the relationR=R c−A(T−T c) b withb=0·6±0·1. BelowT c the conductivities of the two phases obey a relation σ2−σ1=B(T c−T)β with β=0·34±0·02, the exponent of the transport coefficient being the same as the exponent of the order parameter, an equilibrium property.

Reversible shear-induced crystallization above equilibrium freezing temperature in a lyotropic surfactant system

We demonstrate a unique shear-induced crystallization phenomenon above the equilibrium freezing temperature (T-K(o)) in weakly swollen isotropic (L-i) and lamellar (L-alpha) mesophases with bilayers formed in a cationic-anionic mixed surfactant system. Synchrotron rheological X-ray diffraction study reveals the crystallization transition to be reversible under shear (i.e., on stopping the shear, the nonequilibrium crystalline phase L-c melts back to the equilibrium mesophase). This is different from the shear-driven crystallization below T-K(o), which is irreversible. Rheological optical observations show that the growth of the crystalline phase occurs through a preordering of the L-i phase to an L-alpha phase induced by shear flow, before the nucleation of the Lc phase. Shear diagram of the L-i phase constructed in the parameter space of shear rate ((gamma)) over dot vs. temperature exhibits L-i -> L-c and L-i -> L-alpha transitions above the equilibrium crystallization temperature (T-K(o)), in addition to the irreversible shear-driven nucleation of L-c in the L-i phase below T-K(o). In addition to revealing a unique class of nonequilibrium phase transition, the present study urges a unique approach toward understanding shear-induced phenomena in concentrated mesophases of mixed amphiphilic systems.

I. The effect of droplet solidification upon two-phase flow in a rocket nozzle. II. A kinetic theory investigation of some condensation-evaporation phenomena by a moment method

Part I:

The perturbation technique developed by Rannie and Marble is used to study the effect of droplet solidification upon two-phase flow in a rocket nozzle. It is shown that under certain conditions an equilibrium flow exists, where the gas and particle phases have the same velocity and temperature at each section of the nozzle. The flow is divided into three regions: the first region, where the particles are all in the form of liquid droplets; a second region, over which the droplets solidify at constant freezing temperature; and a third region, where the particles are all solid. By a perturbation about the equilibrium flow, a solution is obtained for small particle slip velocities using the Stokes drag law and the corresponding approximation for heat transfer between the particle and gas phases. Singular perturbation procedure is required to handle the problem at points where solidification first starts and where it is complete. The effects of solidification are noticeable.

Part II:

When a liquid surface, in contact with only its pure vapor, is not in the thermodynamic equilibrium with it, a net condensation or evaporation of fluid occurs. This phenomenon is studied from a kinetic theory viewpoint by means of moment method developed by Lees. The evaporation-condensation rate is calculated for a spherical droplet and for a liquid sheet, when the temperatures and pressures are not too far removed from their equilibrium values. The solutions are valid for the whole range of Knudsen numbers from the free molecule to the continuum limit. In the continuum limit, the mass flux rate is proportional to the pressure difference alone.

Modelling the fluid dynamics and coupled phenomena in arc weld pools

A 3D model of melt pool created by a moving arc type heat sources has been developed. The model solves the equations of turbulent fluid flow, heat transfer and electromagnetic field to demonstrate the flow behaviour phase-change in the pool. The coupled effects of buoyancy, capillary (Marangoni) and electromagnetic (Lorentz) forces are included within an unstructured finite volume mesh environment. The movement of the welding arc along the workpiece is accomplished via a moving co-ordinator system. Additionally a method enabling movement of the weld pool surface by fluid convection is presented whereby the mesh in the liquid region is allowed to move through a free surface. The surface grid lines move to restore equilibrium at the end of each computational time step and interior grid points then adjust following the solution of a Laplace equation.

Neutral and Charged 1-Butyl-3-methylimidazolium Triflate Clusters: Equilibrium Concentration in the Vapor Phase and Thermal Properties of Nanometric Droplets

Ground state energy, structure, and harmonic vibrational modes of 1-butyl-3-methylimidazolium triflate ([bmim][Tf]) clusters have been computed using an all-atom empirical potential model. Neutral and charged species have been considered up to a size (30 [bmim][Tf] pairs) well into the nanometric range. Free energy computations and thermodynamic modeling have been used to predict the equilibrium composition of the vapor phase as a function of temperature and density. The results point to a nonnegligible concentration of very small charged species at pressures (P ~ 0.01 Pa) and temperatures (T 600 K) at the boundary of the stability range of [bmim][Tf]. Thermal properties of nanometric neutral droplets have been investigated in the 0 T 700 K range. A near-continuous transition between a liquidlike phase at high T and a solidlike phase at low T takes place at T ~ 190 K in close correspondence with the bulk glass point Tg ~ 200 K. Solidification is accompanied by a transition in the dielectric properties of the droplet, giving rise to a small permanent dipole embedded into the solid cluster. The simulation results highlight the molecular precursors of several macroscopic properties and phenomena and point to the close competition of Coulomb and dispersion forces as their common origin.

Analysis of wave phenomena in a morphogenetic mechanochemical model and an application to post-fertilization waves on eggs

Wave solutions to a mechanochemical model for cytoskeletal activity are studied and the results applied to the waves of chemical and mechanical activity that sweep over an egg shortly after fertilization. The model takes into account the calcium-controlled presence of actively contractile units in the cytoplasm, and consists of a viscoelastic force equilibrium equation and a conservation equation for calcium. Using piecewise linear caricatures, we obtain analytic solutions for travelling waves on a strip and demonstrate uiat the full nonlinear system behaves as predicted by the analytic solutions. The equations are solved on a sphere and the numerical results are similar to the analytic solutions. We indicate how the speed of the waves can be used as a diagnostic tool with which the chemical reactivity of the egg surface can be measured.

An equilibrium approach to modelling social interaction

The aim of this work is to put forward a statistical mechanics theory of social interaction, generalizing econometric discrete choice models. After showing the formal equivalence linking econometric multinomial logit models to equilibrium statical mechanics, a multi- population generalization of the Curie-Weiss model for ferromagnets is considered as a starting point in developing a model capable of describing sudden shifts in aggregate human behaviour. Existence of the thermodynamic limit for the model is shown by an asymptotic sub-additivity method and factorization of correlation functions is proved almost everywhere. The exact solution for the model is provided in the thermodynamical limit by nding converging upper and lower bounds for the system's pressure, and the solution is used to prove an analytic result regarding the number of possible equilibrium states of a two-population system. The work stresses the importance of linking regimes predicted by the model to real phenomena, and to this end it proposes two possible procedures to estimate the model's parameters starting from micro-level data. These are applied to three case studies based on census type data: though these studies are found to be ultimately inconclusive on an empirical level, considerations are drawn that encourage further refinements of the chosen modelling approach, to be considered in future work.

Modeling of complex chemical reactions and macromolecular orientation phenomena in confined geometries

The cooperative motion algorithm was applied on the molecular simulation of complex chemical reactions and macromolecular orientation phenomena in confined geometries. First, we investigated the case of equilibrium step-growth polymerization in lamellae, pores and droplets. In such systems, confinement was quantified as the area/volume ratio. Results showed that, as confinement increases, polymerization becomes slower and the average molecular weight (MW) at equilibrium decreases. This is caused by the sterical hindrance imposed by the walls since chain growth reactions in their close vicinity have less realization possibilities. For reactions inside droplets at surfaces, contact angles usually increased after polymerization to compensate conformation restrictions imposed by confinement upon growing chains. In a second investigation, we considered monodisperse and chemically inert chains and focused on the effect of confinement on chain orientation. Simulations of thin polymer films showed that chains are preferably oriented parallel to the surface. Orientation increases as MW increases or as film thickness d decreases, in qualitative agreement with experiments with low MW polystyrene. It is demonstrated that the orientation of simulated chains results from a size effect, being a function of the ratio between chain end-to-end distance and d. This study was complemented by experiments with thin films of pi-conjugated polymers like MEH-PPV. Anisotropic refractive index measurements were used to analyze chain orientation. With increasing MW, orientation is enhanced. However, for MEH-PPV, orientation does not depend on d even at thicknesses much larger than the chain contour length. This contradiction with simulations was discussed by considering additional causes for orientation, for instance the appearance of nematic-like ordering in polymer films. In another investigation, we simulated droplet evaporation at soluble surfaces and reproduced the formation of wells surrounded by ringlike deposits at the surface, as observed experimentally. In our simulations, swollen substrate particles migrate to the border of the droplet to minimize the contact between solvent and vacuum, which costs the most energy. Deposit formation in the beginning of evaporation results in pinning of the droplet. When polymer chains at the substrate surface have strong uniaxial orientation, the resulting pattern is no longer similar to a ring but to a pair of half-moons. In a final stage, as an extension for the model developed for polymerization in nanoreactors, we studied the effect of geometrical confinement on a hypothetical oscillating reaction following the mechanism of the so called periodically forced Brusselator. It was shown that a reaction which is chaotic in the bulk may be driven to periodicity by confinement and vice-versa, opening new perspectives for chaos control.