Application of standard and refined heat balance integral methods to one-dimensional Stefan problems

The work in this paper concerns the study of conventional and refined heat balance integral methods for a number of phase change problems. These include standard test problems, both with one and two phase changes, which have exact solutions to enable us to test the accuracy of the approximate solutions. We also consider situations where no analytical solution is available and compare these to numerical solutions. It is popular to use a quadratic profile as an approximation of the temperature, but we show that a cubic profile, seldom considered in the literature, is far more accurate in most circumstances. In addition, the refined integral method can give greater improvement still and we develop a variation on this method which turns out to be optimal in some cases. We assess which integral method is better for various problems, showing that it is largely dependent on the specified boundary conditions.

On the origin of the ultradeep East Barents Sea basin

Very large subsidence, with up to 20 km thick sediment layers, is observed in the East Barents Sea basin. Subsidence started in early Paleozoic, accelerated in Permo-Triassic times, finished during the middle Cretaceous, and was followed by moderate uplift in Cenozoic times. The observed gravity signal suggests that the East Barents Sea is at present in isostatic balance and indicates that a mass excess is required in the lithosphere to produce the observed large subsidence. Several origins have been proposed for the mass excess. We use 1-D thermokinematic modeling and 2-D isostatic density models of continental lithosphere to evaluate these competing hypotheses. The crustal density in 2-D thermokinematic models resulting from pressure-, temperature-, and composition-dependent phase change models is computed along transects crossing the East Barents Sea. The results indicate the following. (1) Extension can only explain the observed subsidence provided that a 10 km thick serpentinized mantle lens beneath the basin center is present. We conclude that this is unlikely given that this highly serpentinized layer should be formed below a sedimentary basin with more than 10 km of sediments and crust at least 10 km thick. (2) Phase changes in a compositionally homogeneous crust do not provide enough mass excess to explain the present-day basin geometry. (3) Phase change induced densification of a preexisting lower crustal gabbroic body, interpreted as a mafic magmatic underplate, can explain the basin geometry and observed gravity anomalies. The following model is proposed for the formation of the East Barents Sea basin: (1) Devonian rifting and extension related magmatism resulted in moderate thinning of the crust and a mafic underplate below the central basin area explaining initial late Paleozoic subsidence. (2) East-west shortening during the Permian and Triassic resulted in densification of the previously emplaced mafic underplated body and enhanced subsidence dramatically, explaining the present-day deep basin geometry.

Mecanismos de deformación en aceros inoxidables austeníticos metaestables

La deformación plástica puede inducir a la transformación de la austenita a martensita en los aceros inoxidables austeníticos metaestables. Para analizar este hecho, el inoxidable austenítico metaestable grado AISI 301 LN fue estudiado en dos condiciones diferentes: recocido y laminado en frío. En el primer caso, el acero era completamente austenítico, mientras que después de la laminación presentaba un importante porcentaje de α’-martensita. Se evaluó el cambio de fase cuando el acero es sometido a ensayos monotónicos y cíclicos, así como cuando ha sido modificada la superficie mediante el granallado o se han realizado tratamientos térmicos de reversión. Se utilizaron diferentes técnicas de caracterización microestructural para detectar y cuantificar la martensita, como microscopía óptica, difracción de rayos-X (DRX) y difracción de electrones retrodispersados (EBSD); como también de caracterización mecánica para evaluar el comportamiento de los aceros, trabajo esencial de fractura (TEF), conformabilidad, fatiga de alto número de ciclos (HCF) y nanoindentación. Los resultados mostraron un incremento en la resistencia mecánica del acero laminado en comparación al acero recocido; este hecho está relacionado con la presencia de martensita originada por la laminación. Sin embargo, en términos de deformación y endurecimiento el acero recocido presenta un mejor desempeño como consecuencia del elevado porcentaje de fase austenítica. Así mismo, revertir la martensita de laminación a austenita y refinar la austenita presente permite obtener un acero con una propiedades mecánicas similares a cuando esta en la condición laminado.

One-dimensional solidification of supercooled melts

In this paper a one-phase supercooled Stefan problem, with a nonlinear relation between the phase change temperature and front velocity, is analysed. The model with the standard linear approximation, valid for small supercooling, is first examined asymptotically. The nonlinear case is more difficult to analyse and only two simple asymptotic results are found. Then, we apply an accurate heat balance integral method to make further progress. Finally, we compare the results found against numerical solutions. The results show that for large supercooling the linear model may be highly inaccurate and even qualitatively incorrect. Similarly as the Stefan number β → 1&sup&+&/sup& the classic Neumann solution which exists down to β =1 is far from the linear and nonlinear supercooled solutions and can significantly overpredict the solidification rate.

Normative Topographic ERP Analyses of Speed of Speech Processing and Grammar Before and After Grammatical Treatment.

Typically developing (TD) preschoolers and age-matched preschoolers with specific language impairment (SLI) received event-related potentials (ERPs) to four monosyllabic speech sounds prior to treatment and, in the SLI group, after 6 months of grammatical treatment. Before treatment, the TD group processed speech sounds faster than the SLI group. The SLI group increased the speed of their speech processing after treatment. Posttreatment speed of speech processing predicted later impairment in comprehending phrase elaboration in the SLI group. During the treatment phase, change in speed of speech processing predicted growth rate of grammar in the SLI group.

Novas fases estacionárias à base de sílica para cromatografia líquida de alta eficiência

The present work reviews recent advances in the preparation of new reversed phase packing materials such as sterically protected, bidentate, hybrid organic-inorganic and monolithic phases and phases containing embedded polar groups. The bonding chemistry involved in the preparation of these phases as well as their advantages over conventional C8 and C18 reversed phases are discussed. Understanding the reasons behind the development of these newer column packings helps analysts select the best stationary phase for a given application.

Statistical Analysis and Numerics of Heat Exchanger Models

The identifiability of the parameters of a heat exchanger model without phase change was studied in this Master’s thesis using synthetically made data. A fast, two-step Markov chain Monte Carlo method (MCMC) was tested with a couple of case studies and a heat exchanger model. The two-step MCMC-method worked well and decreased the computation time compared to the traditional MCMC-method. The effect of measurement accuracy of certain control variables to the identifiability of parameters was also studied. The accuracy used did not seem to have a remarkable effect to the identifiability of parameters. The use of the posterior distribution of parameters in different heat exchanger geometries was studied. It would be computationally most efficient to use the same posterior distribution among different geometries in the optimisation of heat exchanger networks. According to the results, this was possible in the case when the frontal surface areas were the same among different geometries. In the other cases the same posterior distribution can be used for optimisation too, but that will give a wider predictive distribution as a result. For condensing surface heat exchangers the numerical stability of the simulation model was studied. As a result, a stable algorithm was developed.

Theoretical investigations of the bulk modulus in the tetra-cubic transition of PbTiO3 material

Resulting from ion displacement in a solid under pressure, piezoelectricity is an electrical polarization that can be observed in perovskite-type electronic ceramics, such as PbTiO3, which present cubic and tetragonal symmetries at different pressures. The transition between these crystalline phases is determined theoretically through the bulk modulus from the relationship between material energy and volume. However, the change in the material molecular structure is responsible for the piezoelectric effect. In this study, density functional theory calculations using the Becke 3-Parameter-Lee-Yang-Parr hybrid functional were employed to investigate the structure and properties associated with the transition state of the tetragonal-cubic phase change in PbTiO3 material.

DESENVOLVIMENTO DE MATERIAIS HÍBRIDOS MICRO-MESOPOROSOS DO TIPO ZSM-12/MCM-41

Micro-mesoporous hybrid materials of ZSM-12/MCM-41 type with different micro- and mesoporosity contributions were prepared by a procedure that uses the desilication of the zeolite in an alkaline medium, followed by recrystallization onto the mesostructure, where the zeolite is used as the silica source in the formation of mesoporous phase. The materials were characterized by X-ray diffraction, nitrogen adsorption-desorption at 77 K, scanning electron microscopy and thermal analysis. The results showed that the methodology utilized is efficient for obtaining hybrid materials of ZSM-12/MCM-41 type with optimized micro-and mesoporosity.

CFD Modelling of Direct Contact Condensation in Suppression Pools by Applying Condensation Models of Separated Flow

The condensation rate has to be high in the safety pressure suppression pool systems of Boiling Water Reactors (BWR) in order to fulfill their safety function. The phenomena due to such a high direct contact condensation (DCC) rate turn out to be very challenging to be analysed either with experiments or numerical simulations. In this thesis, the suppression pool experiments carried out in the POOLEX facility of Lappeenranta University of Technology were simulated. Two different condensation modes were modelled by using the 2-phase CFD codes NEPTUNE CFD and TransAT. The DCC models applied were the typical ones to be used for separated flows in channels, and their applicability to the rapidly condensing flow in the condensation pool context had not been tested earlier. A low Reynolds number case was the first to be simulated. The POOLEX experiment STB-31 was operated near the conditions between the ’quasi-steady oscillatory interface condensation’ mode and the ’condensation within the blowdown pipe’ mode. The condensation models of Lakehal et al. and Coste & Lavi´eville predicted the condensation rate quite accurately, while the other tested ones overestimated it. It was possible to get the direct phase change solution to settle near to the measured values, but a very high resolution of calculation grid was needed. Secondly, a high Reynolds number case corresponding to the ’chugging’ mode was simulated. The POOLEX experiment STB-28 was chosen, because various standard and highspeed video samples of bubbles were recorded during it. In order to extract numerical information from the video material, a pattern recognition procedure was programmed. The bubble size distributions and the frequencies of chugging were calculated with this procedure. With the statistical data of the bubble sizes and temporal data of the bubble/jet appearance, it was possible to compare the condensation rates between the experiment and the CFD simulations. In the chugging simulations, a spherically curvilinear calculation grid at the blowdown pipe exit improved the convergence and decreased the required cell count. The compressible flow solver with complete steam-tables was beneficial for the numerical success of the simulations. The Hughes-Duffey model and, to some extent, the Coste & Lavi´eville model produced realistic chugging behavior. The initial level of the steam/water interface was an important factor to determine the initiation of the chugging. If the interface was initialized with a water level high enough inside the blowdown pipe, the vigorous penetration of a water plug into the pool created a turbulent wake which invoked the chugging that was self-sustaining. A 3D simulation with a suitable DCC model produced qualitatively very realistic shapes of the chugging bubbles and jets. The comparative FFT analysis of the bubble size data and the pool bottom pressure data gave useful information to distinguish the eigenmodes of chugging, bubbling, and pool structure oscillations.

Impact of considering the effects of latent heat on caldera-hosted hydrothermal systems

To investigate the thennal effects of latent heat in hydrothennal settings, an extension was made to the existing finite-element numerical modelling software, Aquarius. The latent heat algorithm was validated using a series of column models, which analysed the effects of penneability (flow rate), thennal gradient, and position along the two-phase curve (pressure). Increasing the flow rate and pressure increases displacement of the liquid-steam boundary from an initial position detennined without accounting for latent heat while increasing the thennal gradient decreases that displacement. Application to a regional scale model of a caldera-hosted hydrothennal system based on a representative suite of calderas (e.g., Yellowstone, Creede, Valles Grande) led to oscillations in the model solution. Oscillations can be reduced or eliminated by mesh refinement, which requires greater computation effort. Results indicate that latent heat should be accounted for to accurately model phase change conditions in hydrothennal settings.

Variation des biomarqueurs dans le spectre visible non résolu de la Terre

L’évolution rapide des technologies de détection et de caractérisation des exoplanètes depuis le début des années 1990 permet de croire que de nouveaux instruments du type Terrestrial Planet Finder (TPF) pourront prendre les premiers spectres d’exoplanètes semblables à la Terre d’ici une ou deux décennies. Dans ce contexte, l’étude du spectre de la seule planète habitée connue, la Terre, est essentielle pour concevoir ces instruments et analyser leurs résultats. Cette recherche présente les spectres de la Terre dans le visible (390-900 nm), acquis lors de 8 nuits d’observation étalées sur plus d’un an. Ces spectres ont été obtenus en observant la lumière cendrée de la Lune avec le télescope de 1.6 m de l’Observatoire du Mont-Mégantic (OMM). La surface de la Lune réfléchissant de manière diffuse la lumière provenant d’une portion de la Terre, ces spectres sont non résolus spatialement. L’évolution de ces spectres en fonction de la lumière réfléchie à différentes phases de Terre est analogue à celle du spectre d’une exoplanète, dont la phase change selon sa position autour de l’étoile. L'eau, l'oxygène et l'ozone de l’atmosphère, détectés dans tous nos spectres, sont des biomarqueurs dont la présence suggère l’habitabilité de la planète et/ou la présence d’une activité biologique. Le Vegetation Red Edge (VRE), une autre biosignature spectrale, dû aux organismes photosynthétiques à la surface, est caractérisé par l’augmentation de la réflectivité autour de 700 nm. Pour les spectres de 5 nuits, cette augmentation a été évaluée entre -5 et 15% ±~5%, après que les contributions de la diffusion de Rayleigh, des aérosols et d’une large bande moléculaire de l’ozone aient été enlevées. Les valeurs mesurées sont cohérentes avec la présence de végétation dans la phase de la Terre contribuant au spectre, mais s’étendent sur une plage de variations plus large que celles trouvées dans la littérature (0-10%). Cela pourrait s’expliquer par des choix faits lors de la réduction des données et du calcul du VRE, ou encore par la présence d’autres éléments de surface ou de l’atmosphère dont la contribution spectrale autour de 700 nm serait variable.

Development of solar air heaters & thermal energy storage system for drying applications in food processing industries

In the present work, the author has designed and developed all types of solar air heaters called porous and nonporous collectors. The developed solar air heaters were subjected to different air mass flow rates in order to standardize the flow per unit area of the collector. Much attention was given to investigate the performance of the solar air heaters fitted with baffles. The output obtained from the experiments on pilot models, helped the installation of solar air heating system for industrial drying applications also. Apart from these, various types of solar dryers, for small and medium scale drying applications, were also built up. The feasibility of ‘latent heat thermal energy storage system’ based on Phase Change Material was also undertaken. The application of solar greenhouse for drying industrial effluent was analyzed in the present study and a solar greenhouse was developed. The effectiveness of Computational Fluid Dynamics (CFD) in the field of solar air heaters was also analyzed. The thesis is divided into eight chapters.

Quasi-biennial oscillation and tracer distributions in a coupled chemistry-climate model

[ 1] We have used a fully coupled chemistry-climate model (CCM), which generates its own wind and temperature quasi-biennial oscillation (QBO), to study the effect of coupling on the QBO and to examine the QBO signals in stratospheric trace gases, particularly ozone. Radiative coupling of the interactive chemistry to the underlying general circulation model tends to prolong the QBO period and to increase the QBO amplitude in the equatorial zonal wind in the lower and middle stratosphere. The model ozone QBO agrees well with Stratospheric Aerosol and Gas Experiment II and Total Ozone Mapping Spectrometer satellite observations in terms of vertical and latitudinal structure. The model captures the ozone QBO phase change near 28 km over the equator and the column phase change near +/- 15 degrees latitude. Diagnosis of the model chemical terms shows that variations in NOx are the main chemical driver of the O-3 QBO around 35 km, i.e., above the O-3 phase change.

Similarity solutions describing the melting of a mushy layer

A model of the melting of a mushy region in the absence of fluid flow is presented. Similarity solutions are obtained which are used to describe melting from a hot plate with and without the generation of a completely molten region. These solutions are extended to describe the melting of a mushy region in contact with a hot liquid. A significant feature of melting mushy regions is that the phase change occurs internally by dissolution. Our solutions for melting of a mushy region are used to investigate this internal phase change and are compared with the classical Neumann solutions for melting of a pure substance.