Metodología de simulación numérica de la transferencia de calor por ebullición en película en torno a superficies esféricas refrigeradas por líquido saturado en convección libre: aplicación a reactores nucleares de lecho de bolas (PBR)

En la presente Tesis se realizó un análisis numérico, usando el código comercial Ansys-Fluent, de la refrigeración de una bola de combustible de un reactor de lecho de bolas (PBR, por sus siglas en inglés), ante un escenario de emergencia en el cual el núcleo sea desensamblado y las bolas se dejen caer en una piscina de agua, donde el mecanismo de transferencia de calor inicialmente sería la ebullición en película, implicando la convección y la radiación al fluido. Previamente se realizaron pruebas de validación, comparando los resultados numéricos con datos experimentales disponibles en la literatura para tres geometrías diferentes, lo cual permitió seleccionar los esquemas y modelos numéricos con mejor precisión y menor costo computacional. Una vez identificada la metodología numérica, todas las pruebas de validación fueron ampliamente satisfactorias, encontrándose muy buena concordancia en el flujo de calor promedio con los datos experimentales. Durante estas pruebas de validación se lograron caracterizar numéricamente algunos parámetros importantes en la ebullición en película con los cuales existen ciertos niveles de incertidumbre, como son el factor de acoplamiento entre convección y radiación, y el factor de corrección del calor latente de vaporización. El análisis térmico de la refrigeración de la bola del reactor por ebullición en película mostró que la misma se enfría, a pesar del calor de decaimiento, con una temperatura superficial de la bola que desciende de forma oscilatoria, debido al comportamiento inestable de la película de vapor. Sin embargo, la temperatura de esta superficie tiene una buena uniformidad, notándose que las áreas mejor y peor refrigeradas están localizadas en la parte superior de la bola. Se observó la formación de múltiples domos de vapor en diferentes posiciones circunferenciales, lo cual causa que el área más caliente de la superficie se localice donde se forman los domos más grandes. La separación entre los domos de vapor fue consistente con la teoría hidrodinámica, con la adición de que la separación entre domos se reduce a medida que evolucionan y crecen, debido a la curvatura de la superficie. ABSTRACT A numerical cooling analysis of a PBR fuel pebble, after an emergency scenario in which the nucleus disassembly is made and the pebbles are dropped into a water pool, transmitting heat by film boiling, involving convection and radiation to the fluid, is carried out in this Thesis. First, were performed validation tests comparing the numerical results with experimental works available for three different geometries, which allowed the selection of numerical models and schemes with better precision and lower computational cost. Once identified the numerical methodology, all validation tests were widely satisfactory, finding very good agreement with experimental works in average heat flux. During these validation tests were achieved numerically characterize some important parameters in film boiling with which there are certain levels of uncertainty, such as the coupling factor between convection and radiation, and the correction factor of the latent heat of vaporization. The thermal analysis of pebble cooling by film boiling shows that despite its decay heat, cooling occurs, with pebble surface temperature descending from an oscillatory manner, due to the instability of the vapor film. However, the temperature of this surface has a good uniformity, noting that the best and worst refrigerated area is located at the top of the pebble. The formation of multiple vapor domes at different circumferential positions is observed, which cause that the hottest area of the surface was located where biggest vapor domes were formed. The separation between vapor domes was consistent with the hydrodynamic theory, with the addition that the separation is reduced as the vapor dome evolves and grows, due to the surface curvature.

Modelo matemático homogéneo e isótropo para la operación y optimización de Sistemas CAES

La operación cíclica de carga y descarga del almacenamiento en los sistemas CAES induce variaciones de presión y temperatura del aire almacenado. Se requiere pues de modelos que predigan la magnitud de dichas actuaciones para lograr una óptima configuración de la planta de potencia. El modelo se deduce directamente de las ecuaciones de conservación de la masa y de la energía para un sistema homogéneo e isótropo. Es resuelto analítica, para ciertos casos límite de interés, y numéricamente. Se somete al modelo a un análisis de sensibilidad para la determinación de los parámetros característicos más influyentes que afectan a las actuaciones de presión y al volumen de almacenamiento requerido. Se observó que el flujo de calor a través del medio circundante al almacenamiento presenta una influencia sobre éste poco despreciable que debe tenerse en cuenta. Se comparan los casos para sistemas adiabáticos e isotermos para calibrar dicha desviación. Además, se esboza un modelo aproximado que describa el comportamiento del almacenamiento con un almacén térmico (AT) acoplado que regule las variaciones de la temperatura del aire inyectado. Abstract Cyclic charge and discharge operation in CAES systems induce pressure and temperatura fluctuations in the storaged air. It is then required a model that properly describes the magnitud of such fluctutions in order to achieve an optimum power plant configuration. This model is directly deduced from the general conservation principles: conservation of mass and energy, applied into homogeneous and isotropic systems. It was solved analytically for certain specific limiting cases of interest as well as numerically. The model was subjected to a sensibility analysis in order to determine the intrinsic parameters that affect the pressure fluctuations as well as the required storage volume the most. Furthermore it was observed that the heat flux through the sourrounding media of the storage has a non-negligible effect that must be taken into account. Adiabatic and isothermal cases are then compared to each other to determine such a deviation. What is more, it is given a general outline about the model which describes the behaviour of the storage coupled with a thermal storage system(TES) that regulates the inyected air temperature.

A model intercomparison of changes in the Atlantic thermohaline circulation in response to increasing atmospheric CO2 concentration

As part of the Coupled Model Intercomparison Project, integrations with a common design have been undertaken with eleven different climate models to compare the response of the Atlantic thermohaline circulation ( THC) to time-dependent climate change caused by increasing atmospheric CO2 concentration. Over 140 years, during which the CO2 concentration quadruples, the circulation strength declines gradually in all models, by between 10 and 50%. No model shows a rapid or complete collapse, despite the fairly rapid increase and high final concentration of CO2. The models having the strongest overturning in the control climate tend to show the largest THC reductions. In all models, the THC weakening is caused more by changes in surface heat flux than by changes in surface water flux. No model shows a cooling anywhere, because the greenhouse warming is dominant.

Polymer Fluid Dynamics: Continuum and Molecular Appoaches

To solve problems in polymer fluid dynamics, one needs the equation of continuity, motion, and energy. The last two equations contain the stress tensor and the heat-flux vector for the material. There are two ways to formulate the stress tensor: (1) one can write a continuum expression for the stress tensor in terms of kinematic tensors, or (2) one can select a molecular model that represents the polymer molecule, and then develop an expression for the stress tensor from kinetic theory. The advantage of the kinetic theory approach is that one gets information about the relation between the molecular structure of the polymers and the rheological properties. In this review, we restrict the discussion primarily to the simplest stress tensor expressions or “constitutive equations” containing from two to four adjustable parameters, although we do indicate how these formulations may be extended to give more complicated expressions. We also explore how these simplest expressions are recovered as special cases of a more general framework, the Oldroyd 8-constant model. The virtue of studying the simplest models is that we can discover some general notions as to which types of empiricisms or which types of molecular models seem to be worth investigating further. We also explore equivalences between continuum and molecular approaches. We restrict the discussion to several types of simple flows, such as shearing flows and extensional flows. These are the flows that are of greatest importance in industrial operations. Furthermore, if these simple flows cannot be well described by continuum or molecular models, then it is not necessary to lavish time and energy to apply them to more complex flow problems.

Geochemistry of sediments of stage 5 of the Nordic seas

Stable isotope, foraminifera and ice rafted detritus (IRD) records covering the last interglacial (the Eemian) from 7 sediment cores in a transect from the Norwegian to the Greenland Sea are presented. The percentages of Neogloboquadrina pachyderma (s.) and Globigerina quinqueloba, foraminiferal content, and to some extent planktonic stable isotope records, demonstrate marked, regional changes in surface water conditions. Importantly, the variability in the abundances of subpolar foraminifera and foraminiferal content are not coherent, implying that these two types of proxies fluctuated independently of each other and most likely reflect changes in sea surface temperature and surface water carbonate productivity, respectively. Paleoceanographic reconstructions demonstrate significant movements of the oceanographic fronts. At the warmest periods, the Arctic front was located far west of the present-day location, at least within the Iceland Sea region. At 126-125 ka, this was most probably due to a stronger or more westerly located Norwegian current. Within the later warm intervals, higher heat flux to the western part of the basin reflects a combination of a stronger Irminger current and/or a weaker east Greenland current. During the main cold spell at ~124 ka, a diffuse Arctic front had a more southeasterly location than today, and intrusion of Atlantic surface waters was probably limited to a narrow corridor in the Eastern Norwegian Sea. A general correspondence between minima in sea surface temperatures and light benthic delta18O may indicate enhanced influx of freshwater to the basin within the cold events. At least in the Norwegian Sea, we find some evidence that the changes in surface water conditions are associated with changes in deep water ventilation. The majority of the fluctuations may be related to occasional breakdown or reduction of the thermohaline circulation within the Nordic seas. In the earliest Eemian, this could result from meltwater forcing. During the remaining part of the last interglacial the fine balance between temperature and salinity, which the deep water formation is depending on, may have been disturbed by periodic increases in fresh water supply or variable influx of warm Atlantic surface waters.

Stable oxygen isotope and Mg/Ca ratios of planktonic foraminifera from the Caribbean and North Atlantic

Variations in the strength of the North Atlantic Ocean thermohaline circulation have been linked to rapid climate changes during the last glacial cycle through oscillations in North Atlantic Deep Water formation and northward oceanic heat flux. The strength of the thermohaline circulation depends on the supply of warm, salty water to the North Atlantic, which, after losing heat to the atmosphere, produces the dense water masses that sink to great depths and circulate back south. Here we analyse two Caribbean Sea sediment cores, combining Mg/Ca palaeothermometry with measurements of oxygen isotopes in foraminiferal calcite in order to reconstruct tropical Atlantic surface salinity during the last glacial cycle. We find that Caribbean salinity oscillated between saltier conditions during the cold oxygen isotope stages 2, 4 and 6, and lower salinities during the warm stages 3 and 5, covarying with the strength of North Atlantic Deep Water formation. At the initiation of the Bølling/Allerød warm interval, Caribbean surface salinity decreased abruptly, suggesting that the advection of salty tropical waters into the North Atlantic amplified thermohaline circulation and contributed to high-latitude warming.

Colour reflectance and core descriptions from sediment cores from the Weddell Sea

The timing of the last maximum extent of the Antarctic ice sheets relative to those in the Northern Hemisphere remains poorly understood. We develop a chronology for the Weddell Sea sector of the East Antarctic Ice Sheet that, combined with ages from other Antarctic ice-sheet sectors, indicates that the advance to and retreat from their maximum extent was within dating uncertainties synchronous with most sectors of Northern Hemisphere ice sheets. Surface climate forcing of Antarctic mass balance would probably cause an opposite response, whereby a warming climate would increase accumulation but not surface melting. Our new data support teleconnections involving sea-level forcing from Northern Hemisphere ice sheets and changes in North Atlantic deep-water formation and attendant heat flux to Antarctic grounding lines to synchronize the hemispheric ice sheets.

Modeling collisional processes in plasmas using discontinuous numerical methods

Thesis (Ph.D.)--University of Washington, 2016-06

Mixing Regimes in the Amundsen Basin of the Arctic Ocean

Thesis (Ph.D.)--University of Washington, 2016-06

Methods of detecting fouling caused by heating of milk

Fouling is the deposition of milk solids on heat transfer sur aces, particularly heat exchangers. It is a major industrial problem, which causes a decrease in heat transfer efficiency and shortens run times. The resultant effect is a decrease in process efficiency and economy. For studying and monitoring deposit formation, suitable fouling detectors or methods of measuring the deposit are required. This can be achieved through direct means, whereby the deposit is analyzed after a certain time, or indirectly through instrumentation for monitoring parameters such as temperature, pressure, flow rate, overall heat transfer coefficient, heat flux, and other physical properties. This article reviews the various reported fouling detection methods.

Optimization of pipeline transport for CO2 sequestration

Coal fired power generation will continue to provide energy to the world for the foreseeable future. However, this energy use is a significant contributor to increased atmospheric CO2 concentration and, hence, global warming. Capture and disposal Of CO2 has received increased R&D attention in the last decade as the technology promises to be the most cost effective for large scale reductions in CO2 emissions. This paper addresses CO2 transport via pipeline from capture site to disposal site, in terms of system optimization, energy efficiency and overall economics. Technically, CO2 can be transported through pipelines in the form of a gas, a supercritical. fluid or in the subcooled liquid state. Operationally, most CO2 pipelines used for enhanced oil recovery transport CO2 as a supercritical fluid. In this paper, supercritical fluid and subcooled liquid transport are examined and compared, including their impacts on energy efficiency and cost. Using a commercially available process simulator, ASPEN PLUS 10.1, the results show that subcooled liquid transport maximizes the energy efficiency and minimizes the Cost Of CO2 transport over long distances under both isothermal and adiabatic conditions. Pipeline transport of subcooled liquid CO2 can be ideally used in areas of cold climate or by burying and insulating the pipeline. In very warm climates, periodic refrigeration to cool the CO2 below its critical point of 31.1 degrees C, may prove economical. Simulations have been used to determine the maximum safe pipeline distances to subsequent booster stations as a function of inlet pressure, environmental temperature and ground level heat flux conditions. (c) 2005 Published by Elsevier Ltd.

On the numerical solution of a Cauchy problem for the Laplace equation via a direct integral equation approach

We investigate the problem of determining the stationary temperature field on an inclusion from given Cauchy data on an accessible exterior boundary. On this accessible part the temperature (or the heat flux) is known, and, additionally, on a portion of this exterior boundary the heat flux (or temperature) is also given. We propose a direct boundary integral approach in combination with Tikhonov regularization for the stable determination of the temperature and flux on the inclusion. To determine these quantities on the inclusion, boundary integral equations are derived using Green’s functions, and properties of these equations are shown in an L2-setting. An effective way of discretizing these boundary integral equations based on the Nystr¨om method and trigonometric approximations, is outlined. Numerical examples are included, both with exact and noisy data, showing that accurate approximations can be obtained with small computational effort, and the accuracy is increasing with the length of the portion of the boundary where the additionally data is given.

Determining the temperature from Cauchy data in corner domains

We consider the problem of reconstruction of the temperature from knowledge of the temperature and heat flux on a part of the boundary of a bounded planar domain containing corner points. An iterative method is proposed involving the solution of mixed boundary value problems for the heat equation (with time-dependent conductivity). These mixed problems are shown to be well-posed in a weighted Sobolev space.

The method of fundamental solutions for free surface Stefan problems

In this paper, free surface problems of Stefan-type for the parabolic heat equation are investigated using the method of fundamental solutions. The additional measurement necessary to determine the free surface could be a boundary temperature, a heat flux or an energy measurement. Both one- and two-phase flows are investigated. Numerical results are presented and discussed.

An iterative method based on boundary integrals for elliptic Cauchy problems in semi-infinite domains

In this study, we investigate the problem of reconstruction of a stationary temperature field from given temperature and heat flux on a part of the boundary of a semi-infinite region containing an inclusion. This situation can be modelled as a Cauchy problem for the Laplace operator and it is an ill-posed problem in the sense of Hadamard. We propose and investigate a Landweber-Fridman type iterative method, which preserve the (stationary) heat operator, for the stable reconstruction of the temperature field on the boundary of the inclusion. In each iteration step, mixed boundary value problems for the Laplace operator are solved in the semi-infinite region. Well-posedness of these problems is investigated and convergence of the procedures is discussed. For the numerical implementation of these mixed problems an efficient boundary integral method is proposed which is based on the indirect variant of the boundary integral approach. Using this approach the mixed problems are reduced to integral equations over the (bounded) boundary of the inclusion. Numerical examples are included showing that stable and accurate reconstructions of the temperature field on the boundary of the inclusion can be obtained also in the case of noisy data. These results are compared with those obtained with the alternating iterative method.