DIRECT NUMERICAL SIMULATION OF INCOMPRESSIBLE MULTIPHASE FLOW WITH PHASE CHANGE

Phase change problems arise in many practical applications such as air-conditioning and refrigeration, thermal energy storage systems and thermal management of electronic devices. The physical phenomenon in such applications are complex and are often difficult to be studied in detail with the help of only experimental techniques. The efforts to improve computational techniques for analyzing two-phase flow problems with phase change are therefore gaining momentum. The development of numerical methods for multiphase flow has been motivated generally by the need to account more accurately for (a) large topological changes such as phase breakup and merging, (b) sharp representation of the interface and its discontinuous properties and (c) accurate and mass conserving motion of the interface. In addition to these considerations, numerical simulation of multiphase flow with phase change introduces additional challenges related to discontinuities in the velocity and the temperature fields. Moreover, the velocity field is no longer divergence free. For phase change problems, the focus of developmental efforts has thus been on numerically attaining a proper conservation of energy across the interface in addition to the accurate treatment of fluxes of mass and momentum conservation as well as the associated interface advection. Among the initial efforts related to the simulation of bubble growth in film boiling applications the work in \cite{Welch1995} was based on the interface tracking method using a moving unstructured mesh. That study considered moderate interfacial deformations. A similar problem was subsequently studied using moving, boundary fitted grids \cite{Son1997}, again for regimes of relatively small topological changes. A hybrid interface tracking method with a moving interface grid overlapping a static Eulerian grid was developed \cite{Juric1998} for the computation of a range of phase change problems including, three-dimensional film boiling \cite{esmaeeli2004computations}, multimode two-dimensional pool boiling \cite{Esmaeeli2004} and film boiling on horizontal cylinders \cite{Esmaeeli2004a}. The handling of interface merging and pinch off however remains a challenge with methods that explicitly track the interface. As large topological changes are crucial for phase change problems, attention has turned in recent years to front capturing methods utilizing implicit interfaces that are more effective in treating complex interface deformations. The VOF (Volume of Fluid) method was adopted in \cite{Welch2000} to simulate the one-dimensional Stefan problem and the two-dimensional film boiling problem. The approach employed a specific model for mass transfer across the interface involving a mass source term within cells containing the interface. This VOF based approach was further coupled with the level set method in \cite{Son1998}, employing a smeared-out Heaviside function to avoid the numerical instability related to the source term. The coupled level set, volume of fluid method and the diffused interface approach was used for film boiling with water and R134a at the near critical pressure condition \cite{Tomar2005}. The effect of superheat and saturation pressure on the frequency of bubble formation were analyzed with this approach. The work in \cite{Gibou2007} used the ghost fluid and the level set methods for phase change simulations. A similar approach was adopted in \cite{Son2008} to study various boiling problems including three-dimensional film boiling on a horizontal cylinder, nucleate boiling in microcavity \cite{lee2010numerical} and flow boiling in a finned microchannel \cite{lee2012direct}. The work in \cite{tanguy2007level} also used the ghost fluid method and proposed an improved algorithm based on enforcing continuity and divergence-free condition for the extended velocity field. The work in \cite{sato2013sharp} employed a multiphase model based on volume fraction with interface sharpening scheme and derived a phase change model based on local interface area and mass flux. Among the front capturing methods, sharp interface methods have been found to be particularly effective both for implementing sharp jumps and for resolving the interfacial velocity field. However, sharp velocity jumps render the solution susceptible to erroneous oscillations in pressure and also lead to spurious interface velocities. To implement phase change, the work in \cite{Hardt2008} employed point mass source terms derived from a physical basis for the evaporating mass flux. To avoid numerical instability, the authors smeared the mass source by solving a pseudo time-step diffusion equation. This measure however led to mass conservation issues due to non-symmetric integration over the distributed mass source region. The problem of spurious pressure oscillations related to point mass sources was also investigated by \cite{Schlottke2008}. Although their method is based on the VOF, the large pressure peaks associated with sharp mass source was observed to be similar to that for the interface tracking method. Such spurious fluctuation in pressure are essentially undesirable because the effect is globally transmitted in incompressible flow. Hence, the pressure field formation due to phase change need to be implemented with greater accuracy than is reported in current literature. The accuracy of interface advection in the presence of interfacial mass flux (mass flux conservation) has been discussed in \cite{tanguy2007level,tanguy2014benchmarks}. The authors found that the method of extending one phase velocity to entire domain suggested by Nguyen et al. in \cite{nguyen2001boundary} suffers from a lack of mass flux conservation when the density difference is high. To improve the solution, the authors impose a divergence-free condition for the extended velocity field by solving a constant coefficient Poisson equation. The approach has shown good results with enclosed bubble or droplet but is not general for more complex flow and requires additional solution of the linear system of equations. In current thesis, an improved approach that addresses both the numerical oscillation of pressure and the spurious interface velocity field is presented by featuring (i) continuous velocity and density fields within a thin interfacial region and (ii) temporal velocity correction steps to avoid unphysical pressure source term. Also I propose a general (iii) mass flux projection correction for improved mass flux conservation. The pressure and the temperature gradient jump condition are treated sharply. A series of one-dimensional and two-dimensional problems are solved to verify the performance of the new algorithm. Two-dimensional and cylindrical film boiling problems are also demonstrated and show good qualitative agreement with the experimental observations and heat transfer correlations. Finally, a study on Taylor bubble flow with heat transfer and phase change in a small vertical tube in axisymmetric coordinates is carried out using the new multiphase, phase change method.

The multiphase starburst-driven galactic wind in NGC 5394

We present a detailed study of the neutral and ionized gas phases in the galactic wind for the nearby starburst galaxy NGC 5394 based on new integral field spectroscopy obtained with the INTEGRAL fibre system at the William Herschel Telescope. The neutral gas phase in the wind is detected via the interstellar Na I D doublet absorption. After a careful removal of the stellar contribution to these lines, a significant amount of neutral gas (∼10^7 M_⊙) is detected in a central region of ∼1.75 kpc size. This neutral gas is blueshifted by ∼165 km s^−1 with respect to the underlying galaxy. The mass outflow of neutral gas is comparable to the star formation rate of the host galaxy. Simultaneously, several emission lines (Hα, [N II], [S II]) are also analysed looking for the ionized warm phase counterpart of the wind. A careful kinematic decomposition of the line profiles reveals the presence of a secondary, broader, kinematic component. This component is found roughly in the same region where the Na I D absorption is detected. It presents higher [N II]/Hα and [S II]/Hα line ratios than the narrow component at the same locations, indicative of contamination by shock ionization. This secondary component also presents blueshifted velocities, although smaller than those measured for the neutral gas, averaging to ∼−30 km s^−1. The mass and mass outflow rate of the wind is dominated by the neutral gas, of which a small fraction might be able to escape the gravitational potential of the host galaxy. The observations in this system can be readily understood within a bipolar gas flow scenario.

Initiation and early crack growth in VHCF of stainless steels : Experimental and theoretical analysis

Mechanical fatigue is a failure phenomenon that occurs due to repeated application of mechanical loads. Very High Cycle Fatigue (VHCF) is considered as the domain of fatigue life greater than 10 million load cycles. Increasing numbers of structural components have service life in the VHCF regime, for instance in automotive and high speed train transportation, gas turbine disks, and components of paper production machinery. Safe and reliable operation of these components depends on the knowledge of their VHCF properties. In this thesis both experimental tools and theoretical modelling were utilized to develop better understanding of the VHCF phenomena. In the experimental part, ultrasonic fatigue testing at 20 kHz of cold rolled and hot rolled stainless steel grades was conducted and fatigue strengths in the VHCF regime were obtained. The mechanisms for fatigue crack initiation and short crack growth were investigated using electron microscopes. For the cold rolled stainless steels crack initiation and early growth occurred through the formation of the Fine Granular Area (FGA) observed on the fracture surface and in TEM observations of cross-sections. The crack growth in the FGA seems to control more than 90% of the total fatigue life. For the hot rolled duplex stainless steels fatigue crack initiation occurred due to accumulation of plastic fatigue damage at the external surface, and early crack growth proceeded through a crystallographic growth mechanism. Theoretical modelling of complex cracks involving kinks and branches in an elastic half-plane under static loading was carried out by using the Distributed Dislocation Dipole Technique (DDDT). The technique was implemented for 2D crack problems. Both fully open and partially closed crack cases were analyzed. The main aim of the development of the DDDT was to compute the stress intensity factors. Accuracy of 2% in the computations was attainable compared to the solutions obtained by the Finite Element Method.

Development, mechanical characterization and qualification in liquid lead of stainless steels for structural applications in the future fission and fusion nuclear reactors

The research activity carried out in the Brasimone Research Center of ENEA concerns the development and mechanical characterization of steels conceived as structural materials for future fission reactors (Heavy Liquid Metal IV Generation reactors: MYRRHA and ALFRED) and for the future fusion reactor DEMO. Within this framework, two parallel lines of research have been carried out: (i) characterization in liquid lead of steels and weldings for the components of the IV Generation fission reactors (GIV) by means of creep and SSRT (Slow Strain Rate Tensile) tests; (ii) development and screening on mechanical properties of RAFM (Reduced Activation Ferritic Martensitic) steels to be employed as structural materials of the future DEMO fusion reactor. The doctoral work represents therefore a comprehensive report of the research carried out on nuclear materials both from the point of view of the qualification of existing (commercial) materials for their application in the typical environmental conditions of 4th generation fission reactors operating with lead as coolant, and from the point of view of the metallurgical study (with annexed microstructural and mechanical characterization of the selected compositions / Thermo Mechanical Treatment (TMT) options) of new compositional variants to be proposed for the “Breeding Blanket” of the future DEMO Fusion Reactor.

Coupled wellbore-reservoir simulation of thermal effects during multiphase CO2 injection

T2Well-ECO2M is a coupled wellbore reservoir simulator still under development at Lawrence Berkeley National Laboratory (USA) with the ability to deal with a mixture of H2O-CO2-NaCl and includes the simulation of CO2 phase transition and multiphase flow. The code was originally developed for the simulation of CO2 injection into deep saline aquifers and the modelling of enhanced geothermal systems; however, the focus of this research was to modify and test T2Well-ECO2M to simulate CO2 injection into depleted gas reservoirs. To this end, the original code was properly changed in a few parts and a dedicated injection case was developed to study CO2 phase transition inside of a wellbore and the corresponding thermal effects. In the first scenario, the injection case was run applying the fully numerical approach of wellbore to formation heat exchange calculation. Results were analysed in terms of wellbore pressure and temperature vertical profiles, wellhead and bottomhole conditions, and characteristic reservoir displacement fronts. Special attention was given to the thorough analysis of bottomhole temperature as the critical parameter for hydrate formation. Besides the expected direct effect of wellbore temperature changes on reservoir conditions, the simulation results indicated also the effect of CO2 phase change in the near wellbore zone on BH pressure distribution. To test the implemented software changes, in a second scenario, the same injection case was reproduced using the improved semi-analytical time-convolution approach for wellbore to formation heat exchange calculation. The comparison of the two scenarios showed that the simulation of wellbore and reservoir parameters after one year of continuous CO2 injection are in good agreement with the computation time to solve the time-convolution semi-analytical reduced. The new updated T2Well-ECO2M version has shown to be a robust and performing wellbore-reservoir simulator that can be also used to simulate the CO2 injection into depleted gas reservoirs.

High-strength steels manufactured via Laser Powder Bed Fusion: effect of post-process heat treatments and surface finishing on microstructure and mechanical properties

Laser Powder Bed Fusion (LPBF) permits the manufacturing of parts with optimized geometry, enabling lightweight design of mechanical components in aerospace and automotive and the production of tools with conformal cooling channels. In order to produce parts with high strength-to-weight ratio, high-strength steels are required. To date, the most diffused high-strength steels for LPBF are hot-work tool steels, maraging and precipitation-hardening stainless steels, featuring different composition, feasibility and properties. Moreover, LPBF parts usually require a proper heat treatment and surface finishing, to develop the desired properties and reduce the high roughness resulting from LPBF. The present PhD thesis investigates the effect of different heat treatments and surface finishing on the microstructure and mechanical properties of a hot-work tool steel and a precipitation-hardening stainless steel manufactured via LPBF. The bibliographic section focuses on the main aspects of LPBF, hot-work tool steels and precipitation-hardening stainless steels. The experimental section is divided in two parts. Part A addresses the effect of different heat treatments and surface finishing on the microstructure, hardness, tensile and fatigue behaviour of a LPBF manufactured hot-work tool steel, to evaluate its feasibility for automotive and racing components. Results indicated the possibility to achieve high hardness and strength, comparable to the conventionally produced steel, but a great sensitivity of fatigue strength on defects and surface roughness resulting from LPBF. Part B investigates the effect of different heat treatments on the microstructure, hardness, tensile and notch-impact behaviour of a LPBF produced precipitation-hardening stainless steel, to assess its feasibility for tooling applications. Results indicated the possibility to achieve high hardness and strength also through a simple Direct Aging, enabling heat treatment simplification by exploiting the microstructural features resulting from LPBF.

Estudo por espectroscopia micro-Raman dos mecanismos de separação de fase em vidros fosfatos de metais de transição

Glass-ceramics are prepared by controlled separation of crystal phases in glasses, leading to uniform and dense grain structures. On the other hand, chemical leaching of soluble crystal phases yields porous glass-ceramics with important applications. Here, glass/ceramic interfaces of niobo-, vanado- and titano-phosphate glasses were studied by micro-Raman spectroscopy, whose spatial resolution revealed the multiphase structures. Phase-separation mechanisms were also determined by this technique, revealing that interface composition remained unchanged as the crystallization front advanced for niobo- and vanadophosphate glasses (interface-controlled crystallization). For titanophosphate glasses, phase composition changed continuously with time up to the equilibrium composition, indicating a spinodal-type phase separation.

Bone response to a Ca- and P-enriched titanium surface obtained by anodization

This study evaluated bone response to a Ca- and P- enriched titanium (Ti) surface treated by a multiphase anodic spark deposition coating (BSP-AK). Two mongrel dogs received bilateral implantation of 3 Ti cylinders (4.1 x 12 mm) in the humerus, being either BSP-AK treated or untreated (machined - control). At 8 weeks postimplantation, bone fragments containing the implants were harvested and processed for histologic and histomorphometric analyses. Bone formation was observed in cortical area and towards the medullary canal associated to approximately 1/3 of implant extension. In most cases, in the medullary area, collagen fiber bundles were detected adjacent and oriented parallel to Ti surfaces. Such connective tissue formation exhibited focal areas of mineralized matrix lined by active osteoblasts. The mean percentages of bone-to-implant contact were 2.3 (0.0-7.2 range) for BSP-AK and 0.4 (0.0-1.3 range) for control. Although the Mann-Whitney test did not detect statistically significant differences between groups, these results indicate a trend of BSP-AK treated surfaces to support contact osteogenesis in an experimental model that produces low bone-to-implant contact values.

Estudo do mecanismo de corrosão por pites em água do mar de aços inoxidáveis supermartensíticos microligados com Nb e Ti

Os aços inoxidáveis supermartensíticos (SMSS) são usados em ambientes agressivos devido à sua boa soldabilidade, boas propriedades mecânicas em temperaturas elevadas e superior resistência à corrosão sob tensão. Aplicações na exploração de petróleo demandam superior combinação de propriedades e os aços inoxidáveis duplex e superduplex têm sido aplicados nessa área, a despeito de seus custos elevados. Os SMSS consistem numa alternativa técnica e econômica ao uso daqueles aços. Nesse trabalho, adições de Nb e Ti foram realizadas com o intuito de minimizar o efeito de sensitização, promover o refino de grãos e foram estudados os aspectos microestruturais e a resistência à corrosão por pites em água do mar. A formação e a evolução dos pites foram acompanhadas por ensaios de corrosão, microscopia óptica e eletrônica, focalizando suas morfologias. O aço com Ti apresentou o melhor desempenho quanto à corrosão, com o maior potencial de corrosão e menor potencial de pite entre os aços em estudo. O aço com Nb, apesar de apresentar potencial de corrosão superior ao do aço sem adição, teve um potencial de pite inferior ao do mesmo.

Distinção entre cachaças destiladas em alambiques e em colunas usando quimiometria

One hundred fifteen cachaça samples derived from distillation in copper stills (73) or in stainless steels (42) were analyzed for thirty five itens by chromatography and inductively coupled plasma optical emission spectrometry. The analytical data were treated through Factor Analysis (FA), Partial Least Square Discriminant Analysis (PLS-DA) and Quadratic Discriminant Analysis (QDA). The FA explained 66.0% of the database variance. PLS-DA showed that it is possible to distinguish between the two groups of cachaças with 52.8% of the database variance. QDA was used to build up a classification model using acetaldehyde, ethyl carbamate, isobutyl alcohol, benzaldehyde, acetic acid and formaldehyde as chemical descriptors. The model presented 91.7% of accuracy on predicting the apparatus in which unknown samples were distilled.

Time-frequency and time-scale analysis of Barkhausen noise signals

Carrying out information about the microstructure and stress behaviour of ferromagnetic steels, magnetic Barkhausen noise (MBN) has been used as a basis for effective non-destructive testing methods, opening new areas in industrial applications. One of the factors that determines the quality and reliability of the MBN analysis is the way information is extracted from the signal. Commonly, simple scalar parameters are used to characterize the information content, such as amplitude maxima and signal root mean square. This paper presents a new approach based on the time-frequency analysis. The experimental test case relates the use of MBN signals to characterize hardness gradients in a AISI4140 steel. To that purpose different time-frequency (TFR) and time-scale (TSR) representations such as the spectrogram, the Wigner-Ville distribution, the Capongram, the ARgram obtained from an AutoRegressive model, the scalogram, and the Mellingram obtained from a Mellin transform are assessed. It is shown that, due to nonstationary characteristics of the MBN, TFRs can provide a rich and new panorama of these signals. Extraction techniques of some time-frequency parameters are used to allow a diagnostic process. Comparison with results obtained by the classical method highlights the improvement on the diagnosis provided by the method proposed.

Modeling of the Barkhausen jump in low carbon steel

This work presents a model for the magnetic Barkhausen jump in low carbon content steels. The outcomes of the model evidence that the Barkhausen jump height depends on the coercive field of the pinning site and on the mean free path of the domain wall between pinning sites. These results are used to deduce the influence of the microstructural features and of the magnetizing parameters on the amplitude and duration of the Barkhausen jumps. In particular, a theoretical expression, establishing the dependence of the Barkbausen jump height on the carbon content and grain size, is obtained. The model also reveals the dependence of the Barkhausen jump on the applied frequency and amplitude. Theoretical and experimental results are presented and compared, being in good agreement. (C) 2008 American Institute of Physics.

Morphological analysis of pitting corrosion in AISI 310S austenitic stainless steel after salt spray exposition

This paper aims to study evolution of increase, distribution and classification of pits in 310S austenitic stainless steels obtained in the state as-received and heat-treated under different exposure times in saline. This work applicability has been based on a technique development for morphologic characterization of localized corrosion associated with description aspects of shapes, size and population-specific parameters. Methodology has been consisted in the following steps: specimens preparation, corrosion tests via salt spray in different conditions, microstructural analysis, pits profiles analysis and images analysis, digital processing and image analysis in order to characterize the pits distribution, morphology and size. Results obtained in digital processing and profiles image analysis have been subjected to statistical analysis using median as parameter in the alloy as received and treated. The alloy as received displays the following morphology: hemispheric pits> transition region A> transition region B> irregular> conic. The pits amount in the treated alloy at each exposure time is: transition region B> hemispherical> transition region A> conic> irregular.

ECAP of Fe. Experiments and simulations of the in-elbow textures

The study of deformation properties of low carbon steels is of particular interest because of their many technological applications. Obtaining fine grained Fe based materials can be approached by one of the several available Severe Plastic Deformation (SPD) techniques. The current paper shows experimental data and simulations of the deformation process of iron samples by Equal Channel Angular Extrusion (ECAE). The samples were extruded in a 120 degrees channel die either by one or a few passes. The heterogeneity and local development of the deformation on the elbow of the channel has been studied by X-ray measuring and simulation of the texture evolution. The Self Consistent models used for simulation allowed the calculation of the spin of the main texture components which agreed pretty well with the experiments.

Annealing effects on microstructure and coercive field of ferritic-martensitic ODS Eurofer steel

Oxide dispersion strengthened reduced-activation ferritic-martensitic steels are promising candidates for applications in future fusion power plants. Samples of a reduced activation ferritic-martensitic 9 wt.%Cr-oxide dispersion strengthened Eurofer steel were cold rolled to 80% reduction in thickness and annealed in vacuum for 1 h from 200 to 1350 degrees C to evaluate its thermal stability. Vickers microhardness testing and electron backscatter diffraction (EBSD) were used to characterize the microstructure. The microstructural changes were also followed by magnetic measurements, in particular the corresponding variation of the coercive field (H(c)), as a function of the annealing treatment. Results show that magnetic measurements were sensitive to detect the changes, in particular the martensitic transformation, in samples annealed above 850 degrees C (austenitic regime). (C) 2010 Elsevier B.V. All rights reserved.