Modelización de los Procesos de Cristalización en Poliolefinas

“Por lo tanto, la cristalización de polímeros se supone, y en las teorías se describe a menudo, como un proceso de múltiples pasos con muchos aspectos físico-químicos y estructurales influyendo en él. Debido a la propia estructura de la cadena, es fácil entender que un proceso que es termodinámicamente forzado a aumentar su ordenamiento local, se vea obstaculizado geométricamente y, por tanto, no puede conducirse a un estado de equilibrio final. Como resultado, se forman habitualmente estructuras de no equilibrio con diferentes características dependiendo de la temperatura, presión, cizallamiento y otros parámetros físico-químicos del sistema”. Estas palabras, pronunciadas recientemente por el profesor Bernhard Wunderlich, uno de los mas relevantes fisico-quimicos que han abordado en las ultimas décadas el estudio del estado físico de las macromoléculas, adelantan lo que de alguna manera se explicita en esta memoria y constituyen el “leitmotiv” de este trabajo de tesis. El mecanismo de la cristalización de polímeros esta aun bajo debate en la comunidad de la física de polímeros y la mayoría de los abordajes experimentales se explican a través de la teoría LH. Esta teoría clásica debida a Lauritzen y Hoffman (LH), y que es una generalización de la teoría de cristalización de una molécula pequeña desde la fase de vapor, describe satisfactoriamente muchas observaciones experimentales aunque esta lejos de explicar el complejo fenómeno de la cristalización de polímeros. De hecho, la formulación original de esta teoría en el National Bureau of Standards, a comienzos de la década de los 70, sufrió varias reformulaciones importantes a lo largo de la década de los 80, buscando su adaptación a los hallazgos experimentales. Así nació el régimen III de cristalización que posibilita la creacion de nichos moleculares en la superficie y que dio pie al paradigma ofrecido por Sadler y col., para justificar los experimentos que se obtenian por “scattering” de neutrones y otras técnicas como la técnica de “droplets” o enfriamiento rapido. Por encima de todo, el gran éxito de la teoría radica en que explica la dependencia inversa entre el tamaño del plegado molecular y el subenfriamiento, definido este ultimo como el intervalo de temperatura que media entre la temperatura de equilibrio y la temperatura de cristalización. El problema concreto que aborda esta tesis es el estudio de los procesos de ordenamiento de poliolefinas con distinto grado de ramificacion mediante simulaciones numéricas. Los copolimeros estudiados en esta tesis se consideran materiales modelo de gran homogeneidad molecular desde el punto de vista de la distribución de tamaños y de ramificaciones en la cadena polimérica. Se eligieron estas poliolefinas debido al gran interes experimental en conocer el cambio en las propiedades fisicas de los materiales dependiendo del tipo y cantidad de comonomero utilizado. Además, son modelos sobre los que existen una ingente cantidad de información experimental, que es algo que preocupa siempre al crear una realidad virtual como es la simulación. La experiencia en el grupo Biophym es que los resultados de simulación deben de tener siempre un correlato mas o menos próximo experimental y ese argumento se maneja a lo largo de esta memoria. Empíricamente, se conoce muy bien que las propiedades físicas de las poliolefinas, en suma dependen del tipo y de la cantidad de ramificaciones que presenta el material polimérico. Sin embargo, tal como se ha explicado no existen modelos teóricos adecuados que expliquen los mecanismos subyacentes de los efectos de las ramas. La memoria de este trabajo es amplia por la complejidad del tema. Se inicia con una extensa introducción sobre los conceptos básicos de una macromolecula que son relevantes para entender el contenido del resto de la memoria. Se definen los conceptos de macromolecula flexible, distribuciones y momentos, y su comportamiento en disolución y fundido con los correspondientes parametros caracteristicos. Se pone especial énfasis en el concepto de “entanglement” o enmaranamiento por considerarse clave a la hora de tratar macromoléculas con una longitud superior a la longitud critica de enmaranamiento. Finaliza esta introducción con una reseña sobre el estado del arte en la simulación de los procesos de cristalización. En un segundo capitulo del trabajo se expone detalladamente la metodología usada en cada grupo de casos. En el primer capitulo de resultados, se discuten los estudios de simulación en disolución diluida para sistemas lineales y ramificados de cadena única. Este caso mas simple depende claramente del potencial de torsión elegido tal como se discute a lo largo del texto. La formación de los núcleos “babys” propuestos por Muthukumar parece que son consecuencia del potencial de torsión, ya que este facilita los estados de torsión mas estables. Así que se propone el análisis de otros potenciales que son igualmente utilizados y los resultados obtenidos sobre la cristalización, discutidos en consecuencia. Seguidamente, en un segundo capitulo de resultados se estudian moleculas de alcanos de cadena larga lineales y ramificados en un fundido por simulaciones atomisticas como un modelo de polietileno. Los resultados atomisticos pese a ser de gran detalle no logran captar en su totalidad los efectos experimentales que se observan en los fundidos subenfriados en su etapa previa al estado ordenado. Por esta razon se discuten en los capítulos 3 y 4 de resultados sistemas de cadenas cortas y largas utilizando dos modelos de grano grueso (CG-PVA y CG-PE). El modelo CG-PE se desarrollo durante la tesis. El uso de modelos de grano grueso garantiza una mayor eficiencia computacional con respecto a los modelos atomísticos y son suficientes para mostrar los fenómenos a la escala relevante para la cristalización. En todos estos estudios mencionados se sigue la evolución de los procesos de ordenamiento y de fusión en simulaciones de relajación isoterma y no isoterma. Como resultado de los modelos de simulación, se han evaluado distintas propiedades fisicas como la longitud de segmento ordenado, la cristalinidad, temperaturas de fusion/cristalizacion, etc., lo que permite una comparación con los resultados experimentales. Se demuestra claramente que los sistemas ramificados retrasan y dificultan el orden de la cadena polimérica y por tanto, las regiones cristalinas ordenadas decrecen al crecer las ramas. Como una conclusión general parece mostrarse una tendencia a la formación de estructuras localmente ordenadas que crecen como bloques para completar el espacio de cristalización que puede alcanzarse a una temperatura y a una escala de tiempo determinada. Finalmente hay que señalar que los efectos observados, estan en concordancia con otros resultados tanto teoricos/simulacion como experimentales discutidos a lo largo de esta memoria. Su resumen se muestra en un capitulo de conclusiones y líneas futuras de investigación que se abren como consecuencia de esta memoria. Hay que mencionar que el ritmo de investigación se ha acentuado notablemente en el ultimo ano de trabajo, en parte debido a las ventajas notables obtenidas por el uso de la metodología de grano grueso que pese a ser muy importante para esta memoria no repercute fácilmente en trabajos publicables. Todo ello justifica que gran parte de los resultados esten en fase de publicación. Abstract “Polymer crystallization is therefore assumed, and in theories often described, to be a multi step process with many influencing aspects. Because of the chain structure, it is easy to understand that a process which is thermodynamically forced to increase local ordering but is geometrically hindered cannot proceed into a final equilibrium state. As a result, nonequilibrium structures with different characteristics are usually formed, which depend on temperature, pressure, shearing and other parameters”. These words, recently written by Professor Bernhard Wunderlich, one of the most prominent researchers in polymer physics, put somehow in value the "leitmotiv "of this thesis. The crystallization mechanism of polymers is still under debate in the physics community and most of the experimental findings are still explained by invoking the LH theory. This classical theory, which was initially formulated by Lauritzen and Hoffman (LH), is indeed a generalization of the crystallization theory for small molecules from the vapor phase. Even though it describes satisfactorily many experimental observations, it is far from explaining the complex phenomenon of polymer crystallization. This theory was firstly devised in the early 70s at the National Bureau of Standards. It was successively reformulated along the 80s to fit the experimental findings. Thus, the crystallization regime III was introduced into the theory in order to explain the results found in neutron scattering, droplet or quenching experiments. This concept defines the roughness of the crystallization surface leading to the paradigm proposed by Sadler et al. The great success of this theory is the ability to explain the inverse dependence of the molecular folding size on the supercooling, the latter defined as the temperature interval between the equilibrium temperature and the crystallization temperature. The main scope of this thesis is the study of ordering processes in polyolefins with different degree of branching by using computer simulations. The copolymers studied along this work are considered materials of high molecular homogeneity, from the point of view of both size and branching distributions of the polymer chain. These polyolefins were selected due to the great interest to understand their structure– property relationships. It is important to note that there is a vast amount of experimental data concerning these materials, which is essential to create a virtual reality as is the simulation. The Biophym research group has a wide experience in the correlation between simulation data and experimental results, being this idea highly alive along this work. Empirically, it is well-known that the physical properties of the polyolefins depend on the type and amount of branches presented in the polymeric material. However, there are not suitable models to explain the underlying mechanisms associated to branching. This report is extensive due to the complexity of the topic under study. It begins with a general introduction to the basics concepts of macromolecular physics. This chapter is relevant to understand the content of the present document. Some concepts are defined along this section, among others the flexibility of macromolecules, size distributions and moments, and the behavior in solution and melt along with their corresponding characteristic parameters. Special emphasis is placed on the concept of "entanglement" which is a key item when dealing with macromolecules having a molecular size greater than the critical entanglement length. The introduction finishes with a review of the state of art on the simulation of crystallization processes. The second chapter of the thesis describes, in detail, the computational methodology used in each study. In the first results section, we discuss the simulation studies in dilute solution for linear and short chain branched single chain models. The simplest case is clearly dependent on the selected torsion potential as it is discussed throughout the text. For example, the formation of baby nuclei proposed by Mutukhumar seems to result from the effects of the torsion potential. Thus, we propose the analysis of other torsion potentials that are also used by other research groups. The results obtained on crystallization processes are accordingly discussed. Then, in a second results section, we study linear and branched long-chain alkane molecules in a melt by atomistic simulations as a polyethylene-like model. In spite of the great detail given by atomistic simulations, they are not able to fully capture the experimental facts observed in supercooled melts, in particular the pre-ordered states. For this reason, we discuss short and long chains systems using two coarse-grained models (CG-PVA and CG-PE) in section 3 and 4 of chapter 2. The CG-PE model was developed during the thesis. The use of coarse-grained models ensures greater computational efficiency with respect to atomistic models and is enough to show the relevant scale phenomena for crystallization. In all the analysis we follow the evolution of the ordering and melting processes by both isothermal and non isothermal simulations. During this thesis we have obtained different physical properties such as stem length, crystallinity, melting/crystallization temperatures, and so on. We show that branches in the chains cause a delay in the crystallization and hinder the ordering of the polymer chain. Therefore, crystalline regions decrease in size as branching increases. As a general conclusion, it seems that there is a tendency in the macromolecular systems to form ordered structures, which can grown locally as blocks, occupying the crystallization space at a given temperature and time scale. Finally it should be noted that the observed effects are consistent with both, other theoretical/simulation and experimental results. The summary is provided in the conclusions chapter along with future research lines that open as result of this report. It should be mentioned that the research work has speeded up markedly in the last year, in part because of the remarkable benefits obtained by the use of coarse-grained methodology that despite being very important for this thesis work, is not easily publishable by itself. All this justify that most of the results are still in the publication phase.

Pulsating emission of droplets from an electrified meniscus.

A numerical description is given for the pulsating emission of droplets from an electrified meniscus of an inviscid liquid of infinite electrical conductivity which is injected at a constant flow rate into a region of uniform, continuous or time periodic, electric field. Under a continuous field, the meniscus attains a periodic regime in which bursts of tiny droplets are emitted from its tip. At low electric fields this regime consists of sequences of emission bursts interspersed with sequences of meniscus oscillations without droplet emission, while at higher fields the bursts occur periodically. These results are in qualitative agreement with experimental results in the literature. Under a time periodic electric field with square waveform, the electric stress that acts on the surface of the liquid while the field is on may generate a tip that emits tiny droplets or may accelerate part of the meniscus and lead to a second emission mode in which a few large droplets are emitted after the electric field is turned off. Conditions under which each emission mode or a combination of the two are realized are discussed for low frequency oscillatory fields. A simplified model is proposed for high electric field frequencies, of the order of the capillary frequency of the meniscus. This model allows computing the average emission rate as a function of the amplitude, duration and bias of the electric field square wave, and shows that droplet emission fails to follow the applied field above a certain frequency

Multi-fluid Eulerian model of an electrospray in a host gas.

An Eulerian multifluid model is used to describe the evolution of an electrospray plume and the flow induced in the surrounding gas by the drag of the electrically charged spray droplets in the space between an injection electrode containing the electrospray source and a collector electrode. The spray is driven by the voltage applied between the two electrodes. Numerical computations and order-of-magnitude estimates for a quiescent gas show that the droplets begin to fly back toward the injection electrode at a certain critical value of the flux of droplets in the spray, which depends very much on the electrical conditions at the injection electrode. As the flux is increased toward its critical value, the electric field induced by the charge of the droplets partially balances the field due to the applied voltage in the vicinity of the injection electrode, leading to a spray that rapidly broadens at a distance from its origin of the order of the stopping distance at which the droplets lose their initial momentum and the effect of their inertia becomes negligible. The axial component of the electric field first changes sign in this region, causing the fly back. The flow induced in the gas significantly changes this picture in the conditions of typical experiments. A gas plume is induced by the drag of the droplets whose entrainment makes the radius of the spray away from the injection electrode smaller than in a quiescent gas, and convects the droplets across the region of negative axial electric field that appears around the origin of the spray when the flux of droplets is increased. This suppresses fly back and allows much higher fluxes to be reached than are possible in a quiescent gas. The limit of large droplet-to-gas mass ratio is discussed. Migration of satellite droplets to the shroud of the spray is reproduced by the Eulerian model, but this process is also affected by the motion of the gas. The gas flow preferentially pushes satellite droplets from the shroud to the core of the spray when the effect of the inertia of the droplets becomes negligible, and thus opposes the well-established electrostatic/inertial mechanism of segregation and may end up concentrating satellite droplets in an intermediate radial region of the spray.

Estudio microestructural y de los procesos de hidratación de cementos con adiciones

En los últimos años, las sociedades industrializadas han tomado una mayor conciencia sobre el problema que suponen las emisiones indiscriminadas de gases de efecto invernadero a la atmósfera. El hormigón, cuyo principal componente es el cemento, es probablemente el material más utilizado en construcción. En la actualidad, las emisiones globales de CO2 debidas a la combustión del CaCO3 del cemento Pórtland representan entre el 5% y el 10% respecto del total. Estos valores son de gran interés si se considera que el compromiso aceptado al firmar el Protocolo de Kioto es de una reducción del 5% antes del año 2020, sobre el total de gases producidos. El principal objetivo del presente trabajo es el estudio microestructural y de los procesos de hidratación de los cementos con adiciones. Para ello se propone contribuir a la investigación sobre nuevos productos cementicios basados en micropartículas esféricas vítreas que pueden adicionarse al cemento antes del proceso de amasado. Los resultados obtenidos se han contrastado con las adiciones convencionales de más uso en la actualidad. El nuevo material basa su composición en la química del aluminio y el silicio. Al disminuir la cantidad de CaCO3, se contribuye al desarrollo sostenible y a la reducción de emisiones de CO2. La patente creada por el Grupo Cementos Pórtland Valderrivas (GCPV), describe el proceso de producción de las cemesferas (WO 2009/007470, 2010). Los productos que forman la materia prima para la elaboración de las cemesferas son arcillas, calizas, margas o productos o subproductos industriales, que tras su molienda, son fundidos mediante un fluido gaseoso a elevada temperatura (entre 1250ºC y 1600ºC). Este proceso permite obtener un producto final en forma esférica maciza o microesfera, que tras estabilizarse mediante un enfriamiento rápido, consigue una alta vitrificación idónea para su reactividad química, con una mínima superficie específica en relación a su masa. El producto final obtenido presenta prácticamente la finura requerida y no precisa ser molido, lo que reduce las emisiones de CO2 por el ahorro de combustible durante el proceso de molienda. El proceso descrito permite obtener un amplio abanico de materiales cementantes que, no solo pueden dar respuesta a los problemas generados por las emisiones de CO2, sino también a la disponibilidad de materiales en países donde hasta el momento no se puede fabricar cemento debido a la falta de calizas. Complementariamente se ha optimizado el método de cálculo del grado de hidratación a partir de los resultados del ensayo de ATD-TG en base a los modelos de cálculo de Bhatty y Pane. El método propuesto permite interpretar el comportamiento futuro del material a partir de la interpolación numérica de la cantidad de agua químicamente enlazada. La evolución del grado de hidratación tiene una relación directa con el desarrollo de la resistencia mecánica del material. Con el fin de caracterizar los materiales de base cemento, se ha llevado a cabo una amplia campaña experimental en pasta de cemento, mortero y hormigón. La investigación abarca tres niveles: caracterización microestructural, macroestructural y caracterización del comportamiento a largo plazo, fundamentalmente durabilidad. En total se han evaluado ocho adiciones diferentes: cuatro adiciones convencionales y cuatro tipos de cemesferas con diferente composición química. Los ensayos a escala microscópica comprenden la caracterización química, granulométrica y de la superficie específica BET de los materiales anhidros, análisis térmico diferencial y termogravimétrico en pasta de cemento y mortero, resonancia magnética de silicio en pasta de cemento, difracción de rayos X de los materiales anhidros y de las probetas de pasta, microscopía electrónica de barrido con analizador de energía dispersiva por rayos X en pasta y mortero, y porosimetría por intrusión de mercurio en mortero. La caracterización macroscópica del material comprende ensayos de determinación del agua de consistencia normal y de los tiempos de inicio y fin de fraguado en pasta de cemento, ensayos de resistencia mecánica a flexión y compresión en probetas prismáticas de mortero, y ensayos de resistencia a compresión en probetas de hormigón. Para caracterizar la durabilidad se han desarrollado ensayos de determinación del coeficiente de migración de cloruros y ensayos de resistividad eléctrica en probetas de mortero. Todos los ensayos enumerados permiten clarificar el comportamiento de las cemesferas y compararlo con las distintas adiciones de uso convencional. Los resultados reflejan un buen comportamiento resistente y durable de los materiales con adición de cemesferas. La caracterización microscópica refleja su relación con las propiedades mesoscópicas y permite comprender mejor la evolución en los procesos de hidratación de las cemesferas. In recent years industrialised societies have become increasingly aware of the problem posed by indiscriminate emission of greenhouse gases into the atmosphere. Concrete, with a main component being cement, is arguably the most widely used construction material. At present, global emissions of CO2 due to the combustion of CaCO3 from Portland cement represent between 5% and 10% of the total. If the requirement of the Kyoto Protocol of a reduction of 5% of the total gas produced before 2020 is considered, then such values are of significant interest. The main objective of this work is the assessment of the microstructure and the hydration processes of cements with additions. Such an examination proposes research into new cementitious products based on vitreous spherical microparticles that may be added to the cement before the mixing process. The results are compared with the most commonly used conventional additions. The new material bases its composition on the chemistry of aluminium and silicates. By decreasing the amount of CaCO3, it is possible both to contribute to sustainable development and reduce CO2 emissions. The patent created by Grupo Cementos Portland Valderrivas (GCPV) describes the production process of microspheres (WO 2009/007470, 2010). The products that form the raw material for manufacture are clays, lime-stone, marl and industrial products or by-products that melt after being ground and fed into a gaseous fluid at high temperatures (1250°C and 1600°C). This process allows the obtaining of a product with a solid-spherical or micro-spherical shape and which, after being stabilised in a solid state by rapid cooling, obtains a high vitrification suitable for chemical reactivity, having a minimal surface in relation to its mass. Given that the final product has the fineness required, it prevents grinding that reduces CO2 emissions by saving fuel during this process. The process, which allows a wide range of cementitious materials to be obtained, not only addresses the problems caused by CO2 emissions but also enhances the availability of materials in countries that until the time of writing have not produced cement due to a lack of limestone. In addition, the calculation of the degree of hydration from the test results of DTA-TG is optimised and based on Bhatty and Pane calculation models. The proposed method allows prediction of the performance of the material from numerical interpolation of the amount of chemically bound water. The degree of hydration has a direct relationship with the development of material mechanical strength. In order to characterise the cement-based materials, an extensive experimental campaign in cement paste, concrete and mortar is conducted. The research comprises three levels: micro-structural characterisation, macro-structural and long-term behaviour (mainly durability). In total, eight additions are assessed: four conventional additions and four types of microspheres with different chemical compositions. The micro-scale tests include characterisation of chemical composition, particle size distribution and the BET specific surface area of anhydrous material, differential thermal and thermogravimetric analysis in cement paste and mortar, silicon-29 nuclear magnetic resonance in cement paste, X-ray diffraction of the anhydrous materials and paste specimens, scanning of electron microscopy with energy dispersive X-ray analyser in cement paste and mortar, and mercury intrusion porosimetry in mortar. The macroscopic material characterisation entails determination of water demand for normal consistency, and initial and final setting times of cement paste, flexural and compressive mechanical strength tests in prismatic mortar specimens, and compressive strength tests in concrete specimens. Tests for determining the chloride migration coefficient are performed to characterise durability, together with electrical resistivity tests in mortar specimens. All the tests listed allow clarification of the behaviour of the microspheres and comparison with the various additions of conventional use. The results show good resistance and durable behaviour of materials with a microsphere addition. Microscopic characterisation reflects their relationship with mesoscopic properties and provides insights into the hydration processes of the microspheres.

Sheath vaporization of a monodisperse fuel-spray jet

The group vaporization of a monodisperse fuel-spray jet discharging into a hot coflowing gaseous stream is investigated for steady flow by numerical and asymptotic methods with a two-continua formulation used for the description of the gas and liquid phases. The jet is assumed to be slender and laminar, as occurs when the Reynolds number is moderately large, so that the boundary-layer form of the conservation equations can be employed in the analysis. Two dimensionless parameters are found to control the flow structure, namely the spray dilution parameter 1, defined as the mass of liquid fuel per unit mass of gas in the spray stream, and the group vaporization parameter e, defined as the ratio of the characteristic time of spray evolution due to droplet vaporization to the characteristic diffusion time across the jet. It is observed that, for the small values of e often encountered in applications, vaporization occurs only in a thin layer separating the spray from the outer droplet-free stream. This regime of sheath vaporization, which is controlled by heat conduction, is amenable to a simplified asymptotic description, independent of ε,in which the location of the vaporization layer is determined numerically as a free boundary in a parabolic problem involving matching of the separate solutions in the external streams, with appropriate jump conditions obtained from analysis of the quasi-steady vaporization front. Separate consideration of dilute and dense sprays, corresponding, respectively, to the asymptotic limits λ<<1 and λ>>1, enables simplified descriptions to be obtained for the different flow variables, including explicit analytic expressions for the spray penetration distance.

Dynamics of thermal ignition of spray flames in mixing layers

Conditions are identified under which analyses of laminar mixing layers can shed light on aspects of turbulent spray combustion. With this in mind, laminar spray-combustion models are formulated for both non-premixed and partially premixed systems. The laminar mixing layer separating a hot-air stream from a monodisperse spray carried by either an inert gas or air is investigated numerically and analytically in an effort to increase understanding of the ignition process leading to stabilization of high-speed spray combustion. The problem is formulated in an Eulerian framework, with the conservation equations written in the boundary-layer approximation and with a one-step Arrhenius model adopted for the chemistry description. The numerical integrations unveil two different types of ignition behaviour depending on the fuel availability in the reaction kernel, which in turn depends on the rates of droplet vaporization and fuel-vapour diffusion. When sufficient fuel is available near the hot boundary, as occurs when the thermochemical properties of heptane are employed for the fuel in the integrations, combustion is established through a precipitous temperature increase at a well-defined thermal-runaway location, a phenomenon that is amenable to a theoretical analysis based on activation-energy asymptotics, presented here, following earlier ideas developed in describing unsteady gaseous ignition in mixing layers. By way of contrast, when the amount of fuel vapour reaching the hot boundary is small, as is observed in the computations employing the thermochemical properties of methanol, the incipient chemical reaction gives rise to a slowly developing lean deflagration that consumes the available fuel as it propagates across the mixing layer towards the spray. The flame structure that develops downstream from the ignition point depends on the fuel considered and also on the spray carrier gas, with fuel sprays carried by air displaying either a lean deflagration bounding a region of distributed reaction or a distinct double-flame structure with a rich premixed flame on the spray side and a diffusion flame on the air side. Results are calculated for the distributions of mixture fraction and scalar dissipation rate across the mixing layer that reveal complexities that serve to identify differences between spray-flamelet and gaseous-flamelet problems.

Source terms for calculations of vaporizing and burning fuel sprays with non-unity Lewis numbers in gases with temperature-dependent thermal conductivities

Liquid-fueled burners are used in a number of propulsion devices ranging from internal combustion engines to gas turbines. The structure of spray flames is quite complex and involves a wide range of time and spatial scales in both premixed and non-premixed modes (Williams 1965; Luo et al. 2011). A number of spray-combustion regimes can be observed experimentally in canonical scenarios of practical relevance such as counterflow diffusion flames (Li 1997), as sketched in figure 1, and for which different microscalemodelling strategies are needed. In this study, source terms for the conservation equations are calculated for heating, vaporizing and burning sprays in the single-droplet combustion regime. The present analysis provides extended formulation for source terms, which include non-unity Lewis numbers and variable thermal conductivities.

Flamelet structures in spray ignition

In typical liquid-fueled burners the fuel is injected as a high-velocity liquid jet that breaks up to form the spray. The initial heating and vaporization of the liquid fuel rely on the relatively large temperatures of the sourrounding gas, which may include hot combustion products and preheated air. The heat exchange between the liquid and the gas phases is enhanced by droplet dispersion arising from the turbulent motion. Chemical reaction takes place once molecular mixing between the fuel vapor and the oxidizer has occurred in mixing layers separating the spray flow from the hot air stream. Since in most applications the injection velocities are much larger than the premixed-flame propagation velocity, combustion stabilization relies on autoignition of the fuel-oxygen mixture, with the combustion stand-off distance being controlled by the interaction of turbulent transport, droplet heating and vaporization, and gas-phase chemical reactions. In this study, conditions are identified under which analyses of laminar flamelets canshed light on aspects of turbulent spray ignition. This study extends earlier fundamental work by Liñan & Crespo (1976) on ignition in gaseous mixing layers to ignition of sprays. Studies of laminar mixing layers have been found to be instrumental in developing un-derstanding of turbulent combustion (Peters 2000), including the ignition of turbulent gaseous diffusion flames (Mastorakos 2009). For the spray problem at hand, the configuration selected, shown in Figure 1, involves a coflow mixing layer formed between a stream of hot air moving at velocity UA and a monodisperse spray moving at velocity USUA. The boundary-layer approximation will be used below to describe the resulting sl ender flow, which exhibits different igniting behaviors depending on the characteristics of t he fuel. In this approximation, consideration of the case U A = U S enables laminar ignition distances to be related to ignition times of unstrained spray flamelets, thereby pro viding quantitative information of direct applicability in regions of low scala r dissipation-rate in turbulent reactive flows (see the discussion in pp. 181–186 of Peters (2000)) . This report is organized as follows. Effects of droplet dispersion dynamics on ignition of sprays in turbulent mixing layers are discussed in Section 2. The formulation f or ignition in laminar mixing layers is outlined in Sections 3 and 4. The results are presented in Section 5. In Section 6, the mixture-fraction field and associated scalar dissipat ion rates for spray ignition are discussed. Finally, some brief conclusions are drawn in Section 7.

Análisis y simulación numérica de problemas electrohidrodinámicos relacionados con el electrospray

La presente tesis es un estudio analítico y numérico del electrospray. En la configuración más sencilla, un caudal constante del líquido a atomizar, que debe tener una cierta conductividad eléctrica, se inyecta en un medio dieléctrico (un gas u otro líquido inmiscible con el primero) a través de un tubo capilar metálico. Entre este tubo y un electrodo lejano se aplica un voltaje continuo que origina un campo eléctrico en el líquido conductor y en el espacio que lo rodea. El campo eléctrico induce una corriente eléctrica en el líquido, que acumula carga en su superficie, y da lugar a un esfuerzo eléctrico sobre la superficie, que tiende a alargarla en la dirección del campo eléctrico. El líquido forma un menisco en el extremo del tubo capilar cuando el campo eléctrico es suficientemente intenso y el caudal suficientemente pequeño. Las variaciones de presión y los esfuerzos viscosos asociados al movimiento del líquido son despreciables en la mayor parte de este menisco, siendo dominantes los esfuerzos eléctrico y de tensión superficial que actúan sobre la superficie del líquido. En el modo de funcionamiento llamado de conochorro, el balance de estos esfuerzos hace que el menisco adopte una forma cónica (el cono de Taylor) en una región intermedia entre el extremo del tubo y la punta del menisco. La velocidad del líquido aumenta al acercarse al vértice del cono, lo cual propicia que las variaciones de la presión en el líquido generadas por la inercia o por la viscosidad entren en juego, desequilibrando el balance de esfuerzos mencionado antes. Como consecuencia, del vértice del cono sale un delgado chorro de líquido, que transporta la carga eléctrica que se acumula en la superficie. La acción del campo eléctrico tangente a la superficie sobre esta carga origina una tracción eléctrica que tiende a alargar el chorro. Esta tracción no es relevante en el menisco, donde el campo eléctrico tangente a la superficie es muy pequeño, pero se hace importante en el chorro, donde es la causa del movimiento del líquido. Lejos del cono, el chorro puede o bien desarrollar una inestabilidad asimétrica que lo transforma en una espiral (whipping) o bien romperse en un spray de gotas prácticamente monodispersas cargadas eléctricamente. La corriente eléctrica transportada por el líquido es la suma de la corriente de conducción en el interior del líquido y la corriente debida a la convección de la carga acumulada en su superficie. La primera domina en el menisco y la segunda en el chorro lejano, mientras que las dos son comparables en una región intermedia de transferencia de corriente situada al comienzo del chorro aunque aguas abajo de la región de transición cono-chorro, en la que el menisco deja de ser un cono de Taylor. Para un campo exterior dado, la acumulación de carga eléctrica en la superficie del líquido reduce el campo eléctrico en el interior del mismo, que llega a anularse cuando la carga alcanza un estado final de equilibrio. El tiempo característico de este proceso es el tiempo de relajación dieléctrica, que es una propiedad del líquido. Cuando el tiempo de residencia del líquido en la región de transición cono-chorro (o en otra región del campo fluido) es grande frente al tiempo de relajación dieléctrica, la carga superficial sigue una sucesión de estados de equilibrio y apantalla al líquido del campo exterior. Cuando esta condición deja de cumplirse, aparecen efectos de relajación de carga, que se traducen en que el campo exterior penetra en el líquido, a no ser que su constante dieléctrica sea muy alta, en cuyo caso el campo inducido por la carga de polarización evita la entrada del campo exterior en el menisco y en una cierta región del chorro. La carga eléctrica en equilibrio en la superficie de un menisco cónico intensifica el campo eléctrico y determina su variación espacial hasta distancias aguas abajo del menisco del orden de su tamaño. Este campo, calculado por Taylor, es independiente del voltaje aplicado, por lo que las condiciones locales del flujo y el valor de la corriente eléctrica son también independientes del voltaje en tanto los tamaños de las regiones que determinan estas propiedades sean pequeños frente al tamaño del menisco. Los resultados experimentales publicados en la literatura muestran que existe un caudal mínimo para el que el modo cono-chorro que acabamos de describir deja de existir. El valor medio y la desviación típica de la distribución de tamaños de las gotas generadas por un electrospray son mínimos cuando se opera cerca del caudal mínimo. A pesar de que los mecanismos responsables del caudal mínimo han sido muy estudiados, no hay aún una teoría completa del mismo, si bien su existencia parece estar ligada a la aparición de efectos de relajación de carga en la región de transición cono-chorro. En esta tesis, se presentan estimaciones de orden de magnitud, algunas existentes y otras nuevas, que muestran los balances dominantes responsables de las distintas regiones de la estructura asintótica de la solución en varios casos de interés. Cuando la inercia del líquido juega un papel en la transición cono-chorro, los resultados muestran que la región de transferencia de corriente, donde la mayor parte de la corriente pasa a la superficie, está en el chorro aguas abajo de la región de transición cono-chorro. Los efectos de relajación de carga aparecen de forma simultánea en el chorro y la región de transición cuando el caudal se disminuye hasta valores de un cierto orden. Para caudales aún menores, los efectos de relajación de carga se notan en el menisco, en una región grande comparada con la de transición cono-chorro. Cuando el efecto de las fuerzas de viscosidad es dominante en la región de transición, la región de transferencia de corriente está en el chorro pero muy próxima a la región de transición cono-chorro. Al ir disminuyendo el caudal, los efectos de relajación de carga aparecen progresivamente en el chorro, en la región de transición y por último en el menisco. Cuando el caudal es mucho mayor que el mínimo del modo cono-chorro, el menisco deja de ser cónico. El campo eléctrico debido al voltaje aplicado domina en la región de transferencia de corriente, y tanto la corriente eléctrica como el tamaño de las diferentes regiones del problema pasan a depender del voltaje aplicado. Como resultado de esta dependencia, el plano caudal-voltaje se divide en diferentes regiones que se analizan separadamente. Para caudales suficientemente grandes, la inercia del líquido termina dominando frente a las fuerzas de la viscosidad. Estos resultados teóricos se han validado con simulaciones numéricas. Para ello se ha formulado un modelo simplificado del flujo, el campo eléctrico y el transporte de carga en el menisco y el chorro del electrospray. El movimiento del líquido se supone casi unidireccional y se describe usando la aproximación de Cosserat para un chorro esbelto. Esta aproximación, ampliamente usada en la literatura, permite simular con relativa facilidad múltiples casos y cubrir amplios rangos de valores de los parámetros reteniendo los efectos de la viscosidad y la inercia del líquido. Los campos eléctricos dentro y fuera del liquido están acoplados y se calculan sin simplificación alguna usando un método de elementos de contorno. La solución estacionaria del problema se calcula mediante un método iterativo. Para explorar el espacio de los parámetros, se comienza calculando una solución para valores fijos de las propiedades del líquido, el voltaje aplicado y el caudal. A continuación, se usa un método de continuación que permite delinear la frontera del dominio de existencia del modo cono-chorro, donde el método iterativo deja de converger. Cuando el efecto de la inercia del líquido domina en la región de transición cono-chorro, el caudal mínimo para el cual el método iterativo deja de converger es del orden del valor estimado del caudal para el que comienza a haber efectos de relajación de carga en el chorro y el cono. Aunque las simulaciones no convergen por debajo de dicho caudal, el valor de la corriente eléctrica para valores del caudal ligeramente mayores parece ajustarse a las estimaciones para caudales menores, reflejando un posible cambio en los balances aplicables. Por el contrario, cuando las fuerzas viscosas dominan en la región de transición, se pueden obtener soluciones estacionarias para caudales bastante menores que aquel para el que aparecen efectos de relajación de carga en la región de transición cono-chorro. Los resultados numéricos obtenidos para estos pequeños caudales se ajustan perfectamente a las estimaciones de orden de magnitud que se describen en la memoria. Por último, se incluyen como anexos dos estudios teóricos que han surgido de forma natural durante el desarrollo de la tesis. El primero hace referencia a la singularidad en el campo eléctrico que aparece en la línea de contacto entre el líquido y el tubo capilar en la mayoría de las simulaciones. Primero se estudia en qué situaciones el campo eléctrico tiende a infinito en la línea de contacto. Después, se comprueba que dicha singularidad no supone un fallo en la descripción del problema y que además no afecta a la solución lejos de la línea de contacto. También se analiza si los esfuerzos eléctricos infinitamente grandes a los que da lugar dicha singularidad pueden ser compensados por el resto de esfuerzos que actúan en la superficie del líquido. El segundo estudio busca determinar el tamaño de la región de apantallamiento en un chorro de líquido dieléctrico sin carga superficial. En esta región, el campo exterior es compensado parcialmente por el campo que induce la carga de polarización en la superficie del líquido, de forma que en el interior del líquido el campo eléctrico es mucho menor que en el exterior. Una región como ésta aparece en las estimaciones cuando los efectos de relajación de carga son importantes en la región de transferencia de corriente en el chorro. ABSTRACT This aim of this dissertation is a theoretical and numerical analysis of an electrospray. In its most simple configuration, a constant flow rate of the liquid to be atomized, which has to be an electrical conductor, is injected into a dielectric medium (a gas or another inmiscible fluid) through a metallic capillary tube. A constant voltage is applied between this tube and a distant electrode that produces an electric field in the liquid and the surrounding medium. This electric field induces an electric current in the liquid that accumulates charge at its surface and leads to electric stresses that stretch the surface in the direction of the electric field. A meniscus appears on the end of the capillary tube when the electric field is sufficiently high and the flow rate is small. Pressure variations and viscous stresses due to the motion of the liquid are negligible in most of the meniscus, where normal electric and surface tension stresses acting on the surface are dominant. In the so-called cone-jet mode, the balance of these stresses forces the surface to adopt a conical shape -Taylor cone- in a intermediate region between the end of the tube and the tip of the meniscus. When approaching the cone apex, the velocity of the liquid increases and leads to pressure variations that eventually disturb the balance of surfaces tension and electric stresses. A thin jet emerges then from the tip of the meniscus that transports the charge accumulated at its surface. The electric field tangent to the surface of the jet acts on this charge and continuously stretches the jet. This electric force is negligible in the meniscus, where the component of the electric field tangent to the surface is small, but becomes very important in the jet. Far from the cone, the jet can either develop an asymmetrical instability named “whipping”, whereby the jet winds into a spiral, or break into a spray of small, nearly monodisperse, charged droplets. The electric current transported by the liquid has two components, the conduction current in the bulk of the liquid and the convection current due to the transport of the surface charge by the flow. The first component dominates in the meniscus, the second one in the far jet, and both are comparable in a current transfer region located in the jet downstream of the cone-jet transition region where the meniscus ceases to be a Taylor cone. Given an external electric field, the charge that accumulates at the surface of the liquid reduces the electric field inside the liquid, until an equilibrium is reached in which the electric field induced by the surface charge counters the external electric field and shields the liquid from this field. The characteristic time of this process is the electric relaxation time, which is a property of the liquid. When the residence time of the liquid in the cone-jet transition region (or in other region of the flow) is greater than the electric relaxation time, the surface charge follows a succession of equilibrium states and continuously shield the liquid from the external field. When this condition is not satisfied, charge relaxation effects appear and the external field penetrates into the liquid unless the liquid permittivity is large. For very polar liquids, the field due to the polarization charge at the surface prevents the external field from entering the liquid in the cone and in certain region of the jet. The charge at the surface of a conical meniscus intensifies the electric field around the cone, determining its spatial variation up to distances downstream of the apex of the order of the size of the meniscus. This electric field, first computed by Taylor, is independent of the applied voltage. Therefore local flow characteristics and the electric current carried by the jet are also independent of the applied voltage provided the size of the regions that determine these magnitudes are small compared with the size of the meniscus. Many experiments in the literature show the existence of a minimum flow rate below which the cone-jet mode cannot be established. The mean value and the standard deviation of the electrospray droplet size distribution are minimum when the device is operated near the minimum flow rate. There is no complete explanation of the minimum flow rate, even though possible mechanisms have been extensively studied. The existence of a minimum flow rate seems to be connected with the appearance of charge relaxation effects in the transition region. In this dissertation, order of magnitude estimations are worked out that show the dominant balances in the different regions of the asymptotic structure of the solution for different conditions of interest. When the inertia of the liquid plays a role in the cone-jet transition region, the region where most of the electric current is transfered to the surface lies in the jet downstream the cone-jet transition region. When the flow rate decreases to a certain value, charge relaxation effects appear simultaneously in the jet and in the transition region. For smaller values of the flow rate, charge relaxation effects are important in a region of the meniscus larger than the transition region. When viscous forces dominate in the flow in the cone-jet transition region, the current transfer region is located in the jet immediately after the transition region. When flow rate is decreased, charge relaxation effects appears gradually, first in the jet, then in the transition region, and finally in the meniscus. When flow rate is much larger than the cone-jet mode minimum, the meniscus ceases to be a cone. The electric current and the structure of the solution begin to depend on the applied voltage. The flow rate-voltage plane splits into different regions that are analyzed separately. For sufficiently large flow rates, the effect of the inertia of the liquid always becomes greater than the effect of the viscous forces. A set of numerical simulations have been carried out in order to validate the theoretical results. A simplified model of the problem has been devised to compute the flow, the electric field and the surface charge in the meniscus and the jet of an electrospray. The motion of the liquid is assumed to be quasi-unidirectional and described by Cosserat’s approximation for a slender jet. This widely used approximation allows to easily compute multiple configurations and to explore wide ranges of values of the governing parameters, retaining the effects of the viscosity and the inertia of the liquid. Electric fields inside and outside the liquid are coupled and are computed without any simplification using a boundary elements method. The stationary solution of the problem is obtained by means of an iterative method. To explore the parameter space, a solution is first computed for a set of values of the liquid properties, the flow rate and the applied voltage, an then a continuation method is used to find the boundaries of the cone-jet mode domain of existence, where the iterative method ceases to converge. When the inertia of the liquid dominates in the cone-jet transition region, the iterative method ceases to converge for values of the flow rate for which order-of-magnitude estimates first predict charge relaxation effects to be important in the cone and the jet. The electric current computed for values of the flow rate slightly above the minimum for which convergence is obtained seems to agree with estimates worked out for lower flow rates. When viscous forces dominate in the transition region, stationary solutions can be obtained for flow rates significantly smaller than the one for which charge relaxation effects first appear in the transition region. Numerical results obtained for those small values of the flow rate agree with our order of magnitude estimates. Theoretical analyses of two issues that have arisen naturally during the thesis are summarized in two appendices. The first appendix contains a study of the singularity of the electric field that most of the simulations show at the contact line between the liquid and the capillary tube. The electric field near the contact line is analyzed to determine the ranges of geometrical configurations and liquid permittivity where a singularity appears. Further estimates show that this singularity does not entail a failure in the description of the problem and does not affect the solution far from the contact line. The infinite electric stresses that appear at the contact line can be effectively balanced by surface tension. The second appendix contains an analysis of the size and slenderness of the shielded region of a dielectric liquid in the absence of free surface charge. In this region, the external electric field is partially offset by the polarization charge so that the inner electric field is much lower than the outer one. A similar region appears in the estimates when charge relaxation effects are important in the current transfer region.

Laser induced forward transfer of silver pastes for printing of fingers in c-si cells.

The main objective of this work is to adapt the Laser Induced Forward Techniques (LIFT), a well- known laser direct writing technique for material transfer, to define metallic contacts (fingers and busbars) onto c-Si cells. The silver paste (with viscosity around 30-50 kcPs) is applied over a glass substrate using a coater. The thickness of the paste can be control changing the deposit parameters. The glass with the silver paste is set at a controlled gap over the c-Si cell. A solid state pulsed laser (532 nm) is focused at the glass/silver interface producing a droplet of silver that it is transferred to the c-Si cell. A scanner is used to print lines. The process parameters (silver paste thickness, gap and laser parameters -spot size, pulse energy and overlapping of pulses) are modified and the morphology of the lines is studied using confocal microscopy. Long lines are printed and the uniformity (in thickness and height) is studied. Some examples of metallization of larger areas (up to 10 cm x 10 cm) are presented.

Conservación de plantas regeneradas in vitro y análisis de la variación somaclonal de Cinchona officinalis, Linneo

Cinchona officinalis (Rubiaceae), especie endémica del Valle de Loja, ubicado en la región sur del Ecuador, es un recurso forestal de importancia medicinal y ecológica, además la especie ha sido catalogada como planta nacional y es un ícono de la región sur por su aporte a la farmacopea mundial. Esta especie, entre los siglos XVII-XIX sufrió una gran presión en sus poblaciones debido a la extracción masiva de la corteza para la cura del paludismo. Aunque la actividad extractiva generó grandes ingresos a la Corona Española y a la región Sur del Ecuador, ésta fue poco o nada sustentable ecológicamente, provocando la desaparición de la especie en muchos sitios de la provincia, pues, en su momento, no se consideraron alternativas de recuperación de las poblaciones naturales. Actualmente la extracción y consumo de la corteza en la zona de origen es baja o nula, sin embargo esta zona enfrenta nuevas amenazas. La deforestación a causa de proyectos de desarrollo en infraestructuras, la práctica de actividades agrícolas y de ganadería, y los efectos del cambio climático han ocasionado, en estos últimos años, la fragmentación de los ecosistemas. La mayoría de los bosques del sur del Ecuador se han convertido en parches aislados (los bosques en los que se distribuye C. officinalis no son la excepción) siendo esta la principal causa para que la especie se encuentre en estado de amenaza. Los individuos de la especie tienen una alta capacidad de rebrote y producen semillas durante todo el año; sin embargo la capacidad germinativa y la tasa de sobrevivencia son bajas, además de estas dificultades la especie requiere de la asociación con otras especies vegetales para su desarrollo, lo cual ha limitado su distribución en pequeños parches aislados. Con esta problemática, la recuperación natural de las poblaciones es una necesidad evidente. Varios trabajos y esfuerzos previos se han realizado a nivel local: i. Identificación de la distribución actual y potencial; ii. Determinación de la fenología y fructificación iii. Programas de educación ambiental, iv. Análisis moleculares para determinar la diversidad genética. v. Ensayos de propagación vegetativa; y otras acciones de tipo cultural. No obstante, el estado de conservación y manejo de las poblaciones naturales no ha mejorado significativamente, siendo necesaria la aplicación de estrategias integradas de conservación in situ y ex situ, que permitan la recuperación y permanencia de las poblaciones naturales a largo plazo. El presente trabajo tiene como fin dar alternativas para el cultivo de tejidos in vitro de Cinchona officinalis centrados en la propagación masiva a partir de semillas, análisis de la fidelidad genética y alternativas de conservación de tejidos. Los objetivos específicos que se plantean son: i. Analizar el proceso de germinación y proliferación in vitro. ii. Evaluar la estabilidad genética en explantes cultivados in vitro, mediante marcadores ISSR. iii. Establecer protocolos de conservación in vitro mediante limitación del crecimiento y criopreservación de segmentos nodales y yemas. Los resultados más significativos de esta investigación fueron: i. El desarrollo de protocolos eficientes para mejorar los porcentajes de germinación y la proliferación de brotes en explantos cultivados in vitro. Para evaluar el efecto de los fenoles sobre la germinación, se determinó el contenido total de fenoles y el porcentaje de germinación en semillas de C. officinalis comparados con una especie de control, C. pubescens. Para inducir a proliferación, se utilizaron segmentos nodales de plántulas germinadas in vitro en medio Gamborg (1968) suplementado con diferentes combinaciones de reguladores de crecimiento (auxinas y citoquininas). Los resultados obtenidos sugieren que el contenido de compuestos fenólicos es alto en las semillas de C. officinalis en comparación con las semillas de C. pubescens. Estos fenoles pueden eliminarse con peróxido de hidrógeno o con lavados de agua para estimular la germinación. La formación de nuevos brotes y callos en la mayoría de las combinaciones de reguladores de crecimiento se observó en un período de 45 días. El mayor porcentaje de proliferación de brotes, formación de callos y presencia de brotes adventicios se obtuvo en medio Gamborg (B5) suplementado con 5.0 mg/l 6-bencil-aminopurina y 3.0 mg/l de ácido indol-3-butírico. ii. La evaluación de la fidelidad genética de los explantes obtenidos con distintas combinaciones de reguladores de crecimiento vegetal y diversos subcultivos. Se realizó el seguimiento a los explantes obtenidos de la fase anterior, determinando el índice de multiplicación y analizando la fidelidad genética de los tejidos obtenidos por las dos vías regenerativas: brotación directa y regeneración de brotes a partir de callos. Este análisis se realizó por amplificación mediante PCR de las secuencias ubicadas entre microsatélites-ISSR (Inter simple sequence repeat). El medio Gamborg (B5) con 3.0 mg/l de AIB y 5.0 mg/l de BAP usado como medio de inducción en la primera etapa de cultivo generó el mayor índice de proliferación (11.5). Un total de 13 marcadores ISSR fueron analizados, 6 de éstos fueron polimórficos. El mayor porcentaje de variación somaclonal fue inducido en presencia de 1.0 mg/l 2,4-D combinado con 0.2 mg/l Kin con un 1.8% en el segundo sub-cultivo de regeneración, la cual incrementó a 3.6% en el tercer sub-cultivo. Todas las combinaciones con presencia de 2,4-D produjeron la formación de callos y presentaron variación genética. Por su parte la fidelidad genética se mantuvo en los sistemas de propagación directa a través de la formación de brotes a partir de meristemos preformados. iii. El establecimiento de protocolos de conservación in vitro y crioconservación de segmentos nodales y yemas. Para la conservación limitando el crecimiento, se cultivaron segmentos nodales en los medios MS y B5 en tres concentraciones de sus componentes (25, 50 y 100%); y en medio B5 más agentes osmóticos como el manitol, sorbitol y sacarosa en diferentes concentraciones (2, 4 y 8%); los cultivos se mantuvieron por 12 meses sin subcultivos. Para el establecimiento de protocolos para la crioconservación (paralización del metabolismo) se usaron yemas axilares y apicales a las cuales se les aplicaron los métodos de encapsulación-deshidratación y vitrificación. La efectividad de los protocolos usados se determinó en función de la sobrevivencia, reducción del crecimiento y regeneración. Los resultados obtenidos en este apartado reflejan que un crecimiento limitado puede mantener tejidos durante 12 meses de almacenamiento, usando medio B5 más manitol entre 2 y 8%. En los protocolos de crioconservación, se obtuvo el mayor porcentaje de recuperación tras la congelación en NL en el tratamiento control seguido por el método crioprotector de encapsulación-deshidratación. Este trabajo brinda alternativas para la propagación de C. officinalis bajo condiciones in vitro, partiendo de material vegetal con alta diversidad genética. El material propagado puede ser fuente de germoplasma para la recuperación y reforzamiento de las poblaciones naturales así como una alternativa de producción para las comunidades locales debido a la demanda actual de corteza de la zona de origen para la elaboración de agua tónica. ABSTRACT Cinchona officinalis (Rubiaceae) is endemic to the Loja Valley, located in the southern area of Ecuador. The importance of this plant as medical and ecological resource is so great that it has been designated as the national flower and is an icon of the southern region for its contribution to the world pharmacopoeia. Between XVII-XIX centuries its population suffered great reduction due to massive harvesting of the bark to cure malaria. Although extraction activity generated large revenues to the Spanish Crown and the southern region of Ecuador, this was not ecologically sustainable, causing the disappearance of the species in many areas of the province, because during that time alternatives to prevent extinction and recover natural populations were not taken in account. Currently the extraction and consumption of bark in the area of origin is almost absent, but this species faces new threats. Deforestation due to infrastructure development, the practice of farming and ranching, and the effects of climate change had led to the fragmentation of ecosystems during the recent years. Most of the forests of southern Ecuador have become isolated patches, including those where C. officinalis is diffused. The lack of suitable habitat is today the main threat for the species. The species has a high capacity for regeneration and produces seeds throughout the year, but the germination rate is low and the growth is slow. In addition, the species requires the association with other plant species to develop. All these factors had limited its distribution to small isolated patches. The natural recovery of populations is essential to face this problem. Several studies and previous efforts had been made at local level: i. Identification of current and potential distribution; ii. Phenology determination. iii. Environmental education programs, iv. Molecular analisis to determine the genetic diversity. v. Testing of vegetative propagation; and other actions of cultural nature. Despite these efforts, the state of conservation and management of natural populations has not improved significantly. Implementation of integrated in situ and ex situ conservation strategies for the recovery and permanence of long-term natural populations is still needed. This work aims to provide alternatives for in vitro culture of tissue of Cinchona officinalis focused on mass propagation from seeds, genetic fidelity analysis and tissue conservation alternatives. The specific aims are: i. Analyze the process of germination and proliferation in vitro. ii. To evaluate the genetic stability of the explants cultured in vitro by ISSR markers. iii. Establish protocols for in vitro conservation by limiting growth and cryopreservation of nodal segments and buds. The most significant results of this research were: i. The development of efficient protocols to improve germination rates and proliferation of buds in explants cultured in vitro. To study the effect of phenols on germination, the total phenolic content and percentage germination was measured in C. officinalis and in a control species, C. pubescens, for comparison. The content of phenolic compounds in C. officinalis seeds is higher than in C. pubescens. These phenols can be removed with hydrogen peroxide or water washes to stimulate germination. To analyze the regeneration, we used nodal explants from seedlings germinated in vitro on Gamborg medium (1968) supplemented with different combinations of growth regulators (auxins and cytokinins) to induce proliferation. The formation of new shoots and calluses was observed within a period of 45 days in most combinations of growth regulators. The highest percentage of shoot proliferation, callus formation and adventitious buds were obtained in B5 medium supplemented with 5.0 mg/l 6-benzyl-aminopurine and 3.0 mg/l indole-3-butyric acid. ii. Evaluating genetic fidelity explants obtained with various combinations of plant growth regulators and different subcultures. The genetic fidelity was analyzed in tissues obtained by the two regenerative pathways: direct sprouting and shoot regeneration from callus. This analysis was performed by PCR amplification of the sequences located between microsatellite-ISSR (Inter Simple Sequence Repeat). Among a total of 13 ISSR markers analyzed, 6 were polymorphic. The highest percentage of somaclonal variation was induced in the presence of 1.0 mg/l 2,4-D combined with 0.2 mg/l Kin with 1.8% in the second round of regeneration, and increased to 3.6% in the third round. The presence of 2,4-D induced genetic variation in all the combinations of growth regulators. Meanwhile genetic fidelity remained systems propagation through direct shoot formation from meristems preformed. iii. Establishing conservation protocols in vitro and cryoconservation of nodal segments and buds. For medium-term conservation (limited growth) nodal segments were cultured in MS and B5 media at three concentrations (25, 50 and 100%); we tested B5 medium with different concentrations of osmotic agents such as mannitol, sorbitol and sucrose (2, 4 and 8%); cultures were maintained for 12 months with regular subculturing. To establish protocols for cryoconservation (cessation of metabolism) different methods of encapsulation-dehydration and vitrification were applied to axillary and apical buds. The effectiveness of the used protocols is determined based on the survival, growth and regeneration success. The results show that these tissues can be maintained in storage for 12 months, using B5 medium plus mannitol between 2 and 8%. The cryoconservation protocol with highest percentage of recovery was obtained by contral treatment, followed by freezing in NL with encapsulation-dehydration method. This work provides alternatives for the propagation in vitro of C. officinalis, starting from plant material with high genetic diversity. The obtained material represents a source of germplasm to support the recovery and strengthening of natural populations as well as a creation of alternative sources for local communities due to the current demand of bark for the preparation of tonic water.

Estudo fenomenológico e numérico do escoamento estratificado óleo-água ondulado e com mistura na interface

Escoamentos bifásicos estão presentes em diversos processos naturais e industriais, como na indústria de petróleo. Podem apresentar-se em diferentes configurações topológicas, ou, padrões de escoamento, entre eles o escoamento estratificado ondulado e o estratificado com mistura na interface. Os escoamentos bifásicos estratificados óleo-água têm sido utilizados como uma forma conveniente de evitar a formação de emulsões de água em óleo em oleodutos e possuem uma ocorrência comum em poços de petróleo direcionais. Quando a onda interfacial ultrapassa determinado limite geométrico e cinemático, surge o fenômeno do entranhamento de gotas, representado por misturas entre as fases junto à interface que promovem um aumento na queda de pressão. Modelos têm sido apresentados pela literatura na tentativa de descrever o fenômeno do entranhamento de gotas. Neste trabalho é apresentada uma nova proposta de modelagem matemática unidimensional para o entranhamento de gotas com o objetivo de melhorar a previsão dos parâmetros envolvidos, em especial, da fração volumétrica de óleo e da queda de pressão bifásica. Também foi utilizada simulação numérica computacional, CFD (Computational Fluid Dynamics), com o uso de software comercial para obtenção dos valores dos parâmetros do escoamento estratificado ondulado óleo-água (fração volumétrica de óleo, queda de pressão, amplitude e comprimento da onda interfacial). Os resultados da modelagem fenomenológica para entranhamento e os de CFD foram comparados com bancos de dados experimentais. Os resultados em CFD mostram concordância com os resultados experimentais, tanto na análise qualitativa das propriedades geométricas das ondas interfaciais, quanto na comparação direta com os dados para fração volumétrica e queda de pressão. Os resultados numéricos da modelagem fenomenológica para fatores de entranhamento apresentam boa concordância com dados da literatura.

Polimerização em emulsão de estireno em microrreator.

A polimerização em emulsão de estireno em um microrreator Syrris de 250 µL com misturador estático junção \"T\" foi estudada em duas etapas. Primeiro somente a fluidodinâmica deste dispositivo não convencional foi avaliada, depois, foi desenvolvida a reação de polimerização de forma a observar como este fator influencia no sistema. Os experimentos foram realizados procurando se atingir maiores conversões, mas mantendo a estabilidade da emulsão. Foi um trabalho exploratório, portanto se assemelha mais a um processo de evolução (evolutionary process). Foram verificados a partir de qual relação das vazões dos dois fluidos ocorre a formação de gotas, e que com o aumento da vazão da fase contínua, aquosa (Qc), mantendo constante a vazão da fase dispersa (Qd), foi verificado uma diminuição do diâmetro das gotas e um regime de fluxo laminar. Posteriormente, realizou-se a polimerização em emulsão do estireno no microrreator, porém com restrições para altas vazões. Os parâmetros de processo testados foram a proporção Qc e Qd, a temperatura e a concentração do iniciador para então verificar o efeito que a variação destas ocasionam na conversão de monômero, no diâmetro e número de partículas e nas massas moleculares médias. A polimerização foi feita para soma das vazões Qc e Qd da ordem de 100 µ L/min, com 15% de monômero na formulação e com o maior tempo de residência possível de 2,5 minutos. Para maiores concentrações de monômero, acima de 15% foi verificado entupimento do canal do microrreator. A taxa de conversão de monômero aumentou com o aumento da temperatura e com o aumento da concentração do iniciador, mas o maior valor atingido foi de apenas 37% devido ao baixo tempo de residência. Nos casos de maiores taxas de conversão, as massas moleculares obtidas foram as menores conforme o esperado pela teoria. Finalmente, os índices de polidispersão (PDI), obtidos foram da ordem de 2,5 a 3,5.

Estudo comparativo de técnicas numéricas de inversão para obtenção de distribuição de tamanho de gotas em emulsões.

O desenvolvimento de algoritmos computacionais para a obtenção de distribuições de tamanho de partícula em dispersões e que utilizam dados espectroscópicos em tempo real e in-line a partir de sensores permitirá uma variedade de aplicações, como o monitoramento de propriedades em fluidos de corte industriais, acompanhamento de processos de polimerização, tratamento de efluentes e sensoriamento atmosférico. O presente estudo tem como objetivo a implementação e comparação de técnicas para resolução de problemas de inversão, desenvolvendo algoritmos que forneçam distribuição de tamanho de partículas em dispersões a partir de dados de espectroscopia UV-Vis-Nir (Ultravioleta, Visível e Infravermelho próximo). Foram implementadas quatro técnicas, sendo uma delas um método alternativo sem a presença de etapas de inversão. Os métodos que utilizaram alguma técnica de inversão evidenciaram a dificuldade em se obter distribuições de tamanho de gotas (DTG) de boa qualidade, enquanto o método alternativo foi aquele que se mostrou mais eficiente e confiável. Este estudo é parte de um programa cooperativo entre a Universidade de São Paulo e a Universidade de Bremen chamado programa BRAGECRIM (Brazilian German Cooperative Research Initiative in Manufacturing) e é financiado pela FAPESP, CAPES, FINEP e CNPq (Brasil) e DFG (Alemanha).

Fluorescent Biological Aerosol Particles in Coastal Canada and the Link to Atmospheric Ice Nuclei

Bioaerosols are a subgroup of atmospheric aerosols and are often linked to the spread of human, animal and plant diseases. Bioaerosols also may play an indirect effect on environmental processes, including the formation of precipitation and alteration of the global climate through their role as nuclei for cloud droplet formation. Several types of biological organisms (e.g., fungi and bacteria) have been shown to be effective ice nuclei (IN) and cloud condensation nuclei (CCN). During 21 days in August 2013 we participated in a collaborative international campaign at a rural, coastal site near the village of Ucluelet on the west coast of Vancouver Island, British Columbia, Canada. The experiments were conducted as part of the NETCARE project (the NETwork on Climate and Aerosols: Addressing Key Uncertainties in Remote Canadian Environments), in part to examine cloud nuclei properties of marine aerosol. The study was conducted from a mobile trailer located approximately 100 m from the coast. A suite of aerosol instrumentation was operated for approximately one month. Key instruments utilized as a part of this thesis include the wideband integrated bioaerosol sensor (WIBS-4A) and the multiple orifice uniform deposition impactor (MOUDI) coupled with an off-line droplet freezing technique (DFT) for the measurement of ice nucleation activity of particles in immersion mode. The WIBS measures the concentration and properties of individual fluorescent particles suspended in the air, which can serve as a proxy for airborne biological particle content. Particles shown to be fluorescent by the WIBS instrument were divided into seven categories based on the pattern of fluorescence each particle exhibited in the three fluorescent channels. Results of the WIBS analysis show that the fluorescent particle concentration in the region correlated well with IN number. The fluorescent particle concentration correlated well with the number of particles shown to be ice active as a function of both particle size and freezing temperature. Correlations involving marine aerosols and marine biological activity indicate that the majority of IN measured at the coastal site likely are not from have marine sources.