Multiphase modeling of melting / solidification with high density variations using XFEM

La modélisation de la cryolite, utilisée dans la fabrication de l’aluminium, implique plusieurs défis, notament la présence de discontinuités dans la solution et l’inclusion de la difference de densité entre les phases solide et liquide. Pour surmonter ces défis, plusieurs éléments novateurs ont été développés dans cette thèse. En premier lieu, le problème du changement de phase, communément appelé problème de Stefan, a été résolu en deux dimensions en utilisant la méthode des éléments finis étendue. Une formulation utilisant un multiplicateur de Lagrange stable spécialement développée et une interpolation enrichie a été utilisée pour imposer la température de fusion à l’interface. La vitesse de l’interface est déterminée par le saut dans le flux de chaleur à travers l’interface et a été calculée en utilisant la solution du multiplicateur de Lagrange. En second lieu, les effets convectifs ont été inclus par la résolution des équations de Stokes dans la phase liquide en utilisant la méthode des éléments finis étendue aussi. Troisièmement, le changement de densité entre les phases solide et liquide, généralement négligé dans la littérature, a été pris en compte par l’ajout d’une condition aux limites de vitesse non nulle à l’interface solide-liquide pour respecter la conservation de la masse dans le système. Des problèmes analytiques et numériques ont été résolus pour valider les divers composants du modèle et le système d’équations couplés. Les solutions aux problèmes numériques ont été comparées aux solutions obtenues avec l’algorithme de déplacement de maillage de Comsol. Ces comparaisons démontrent que le modèle par éléments finis étendue reproduit correctement le problème de changement phase avec densités variables.

Estabilidade de taludes em escavação: um caso de estudo

Esta dissertação aborda a análise de estabilidade de taludes em escavação, para infraestruturas de engenharia civil, nomeadamente obras rodoviárias. Este trabalho encontra-se dividido em quatro partes. Na primeira parte é abordado o papel das descontinuidades no maciço. A segunda parte, na representação dos dados adquiridos. A terceira parte centra-se nos fundamentos teóricos sobre os taludes. Por último, a quarta parte, onde se desenvolve o caso de estudo, no âmbito de um estágio curricular na Geoma. O estudo refere-se a análise de estabilidade de taludes para uma obra rodoviária denominada Sublanço – Nó de Ligação do IP4/Túnel do Marão.

Anomaly Detection in Noisy Images

Finding rare events in multidimensional data is an important detection problem that has applications in many fields, such as risk estimation in insurance industry, finance, flood prediction, medical diagnosis, quality assurance, security, or safety in transportation. The occurrence of such anomalies is so infrequent that there is usually not enough training data to learn an accurate statistical model of the anomaly class. In some cases, such events may have never been observed, so the only information that is available is a set of normal samples and an assumed pairwise similarity function. Such metric may only be known up to a certain number of unspecified parameters, which would either need to be learned from training data, or fixed by a domain expert. Sometimes, the anomalous condition may be formulated algebraically, such as a measure exceeding a predefined threshold, but nuisance variables may complicate the estimation of such a measure. Change detection methods used in time series analysis are not easily extendable to the multidimensional case, where discontinuities are not localized to a single point. On the other hand, in higher dimensions, data exhibits more complex interdependencies, and there is redundancy that could be exploited to adaptively model the normal data. In the first part of this dissertation, we review the theoretical framework for anomaly detection in images and previous anomaly detection work done in the context of crack detection and detection of anomalous components in railway tracks. In the second part, we propose new anomaly detection algorithms. The fact that curvilinear discontinuities in images are sparse with respect to the frame of shearlets, allows us to pose this anomaly detection problem as basis pursuit optimization. Therefore, we pose the problem of detecting curvilinear anomalies in noisy textured images as a blind source separation problem under sparsity constraints, and propose an iterative shrinkage algorithm to solve it. Taking advantage of the parallel nature of this algorithm, we describe how this method can be accelerated using graphical processing units (GPU). Then, we propose a new method for finding defective components on railway tracks using cameras mounted on a train. We describe how to extract features and use a combination of classifiers to solve this problem. Then, we scale anomaly detection to bigger datasets with complex interdependencies. We show that the anomaly detection problem naturally fits in the multitask learning framework. The first task consists of learning a compact representation of the good samples, while the second task consists of learning the anomaly detector. Using deep convolutional neural networks, we show that it is possible to train a deep model with a limited number of anomalous examples. In sequential detection problems, the presence of time-variant nuisance parameters affect the detection performance. In the last part of this dissertation, we present a method for adaptively estimating the threshold of sequential detectors using Extreme Value Theory on a Bayesian framework. Finally, conclusions on the results obtained are provided, followed by a discussion of possible future work.

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.

Integrating life design counseling and psychotherapy: Possibilities and practices.

The debate on the relation between career counseling and psychotherapy revealed the close relation between psychosocial and career issues. The connection between these two approaches paves the way for the integration of career counseling into psychotherapy. One of the key challenges of integrating career counseling into psychotherapy is to avoid discontinuities and inconsistencies in the psychotherapeutic intervention. The aim of this article is, therefore, to present the Systematic Treatment Selection perspective as a framework for an effective integration of psychotherapy and career counseling. Then, Life Design Counseling is described, highlighting its possibilities to fit the psychotherapeutic approach. Finally, an illustrative case is presented and discussed.

Implementação de programas de educação ambiental nas escolas municipais de Curitiba, Paraná, Brasil

O presente trabalho pretende analisar a implementação dos programas de Educação Ambiental nas escolas de ensino fundamental, da Rede Municipal de Ensino Público, do Município de Curitiba, no estado do Paraná, Brasil. Esta análise tem por objeto de estudo as ações desenvolvidas pelas escolas, no que se refere às condições de saúde e bem-estar, a partir do que estabelecem os Parâmetros Curriculares Nacionais. Buscou-se compreender como os conteúdos definidos por ‘temas transversais' propiciam a constituição de um espaço de transformação em nível municipal. Este trabalho se constrói com base em discussões sobre o conhecimento produzido na área de Educação Ambiental e pelo levantamento de informações por meio de entrevistas e questionários junto aos profissionais ligados aos órgãos municipais responsáveis e às escolas. A análise realizada demonstrou que os projetos implementados vêm surtindo efeito positivo sobre as transformações dos pensamentos relativos às questões ambientais, sofrendo descontinuidades devido às mudanças nas gestões políticas. ABSTRACT; This study aims to analyze the implementation of the Environment Education Programs in the schools of the Municipal Public Education Network, in the City of Curitiba, in the state of Paraná, Brazil. The object study on this analysis is the actions developed by the schools, referring the conditions of health and well-being, from what is established by the National Curricular Parameters. We try to understand how the contents defined for 'transversal subjects' can propitiate the constitution of a space of transformation in municipal level. This work is built on discussion about the knowledge produced in the area of Environment Education besides the survey of informations through interviews and questionnaires with professionals connected to the responsible offices and the schools. This analysis demonstrated that the implemented projects is causing positive effects on the transformations of the relative thoughts to the environment questions, come suffering discontinuities due to the changes in the managements politics. RESUMEN; Este proyecto tiene por objetivo examinar la implementación de los programas de educación ambiental en las escuelas de educación básica de la Red Municipal de Educación Pública, en la ciudad de Curitiba, Paraná, Brasil. El objeto de estudio de la análisis son las medidas adoptadas por las escuelas, con respecto a la salud y el bienestar, de acuerdo con lo establecido por los Parámetros Curriculares Nacionales. Tratamos de comprender como el contenido, tal como se define por "temas transversales" puede ofertar un espacio de transformación en el nivel municipal. Este trabajo se basa en los debates de los conocimientos producidos en el ámbito de la educación ambiental y por las informaciones planteadas a través de entrevistas y cuestionarios con los profesionales de los departamentos responsables y en las escuelas. El análisis mostró que los proyectos ejecutados han tenido efecto positiva en la transformación de ideas sobre cuestiones ambientales, sufriendo discontinuidades debidas a cambios en las gestiones políticas.

Controle de sistema não linear de separação gás/líquido sob ação de efeitos transientes fluídicos por sliding control

Dissertação (mestrado)—Universidade de Brasília, Faculdade de Tecnologia, Departamento de Engenharia Mecânica, 2015.

Intraplate versus ridge volcanism on the Pacific-Antarctic Ridge near 37°S -111°W

Exploration of the Foundation Volcanic Chain (33 degrees S-131 degrees W; 37 degrees S-111 degrees W) revealed the existence of different magmatic provinces with relation to their geological settings. (1) The Pacific-Antarctic Ridge (PAR) is made up of several en echelon segments where both glassy midocean ridge basalts (MORBs) with low incompatible elements (K2O<200 ppm, Zr<120 ppm and Ce <20 ppm) as well as andesites and dacites have erupted, (2) Oblique Ridges located up to 300 lan from the PAR axis are topped with seamounts made up essentially of transitional (T) and enriched (E) MORBs with intermediate incompatible elements (K2O=0.11-0.40 %, Zr=70-140 ppm and Ce=15-30 ppm), (3) the Foundation Seamounts (FS) consisting essentially of isolated volcanoes which have erupted alkalic lavas (alkali basalt, trachybasalt and trachyandesite) with high incompatible elements (K2O (0.50-1.1 %, Zr (>150 ppm) and Ce (>48 ppm)) at about 306-1300 km from the PAR axis, (4) The Old Pacific Seamounts built on a crust older than 23 m. y. located west of longitude 124 degrees W (> 1300 km from the PAR axis) consist of T and EMORB. On the PAR axis, extensive crystal fractionation (>65%) produced the silicic lavas. On the basis of Pacific plate reconstruction using a half spreading rate of about 50 mm/yr and integrating the observed compositional changes with respect to the structural settings, it is inferred that the last volcanic events giving rise to the FS took place at about 110 km from the PAR axis about 5 m. y. ago. The Oblique Ridges built between 5 m. y. and <1 m. y. are believed to represent ancient leaky transforms and/or large discontinuities between accreting ridge segments filled by volcanic cones during the interaction (mixing) of the enriched plume components of the FS with PAR depleted (MORB type) magmatism. The Old Pacific Seamounts built on ancient crust (>23 m. y.) with MORB volcanics comparable to those of the the Oblique Ridge-PAR provinces, could also have been formed by an interaction between the Foundation Seamount (dredge site 28) hotspot magmatism and that of an ancient accreting ridge magmatism precursor of the PAR.

LHC main dipole magnet circuits: sustaining near-nominal beam energies

Crossing the Franco-Swiss border, the Large Hadron Collider (LHC), designed to collide 7 TeV proton beams, is the world's largest and most powerful particle accelerator the operation of which was originally intended to commence in 2008. Unfortunately, due to an interconnect discontinuity in one of the main dipole circuit's 13 kA superconducting busbars, a catastrophic quench event occurred during initial magnet training, causing significant physical system damage. Furthermore, investigation into the cause found that such discontinuities were not only present in the circuit in question, but throughout the entire LHC. This prevented further magnet training and ultimately resulted in the maximum sustainable beam energy being limited to approximately half that of the design nominal, 3.5-4 TeV, for the first three years of operation (Run 1, 2009-2012) and a major consolidation campaign being scheduled for the first long shutdown (LS 1, 2012-2014). Throughout Run 1, a series of studies attempted to predict the amount of post-installation training quenches still required to qualify each circuit to nominal-energy current levels. With predictions in excess of 80 quenches (each having a recovery time of 8-12+ hours) just to achieve 6.5 TeV and close to 1000 quenches for 7 TeV, it was decided that for Run 2, all systems be at least qualified for 6.5 TeV operation. However, even with all interconnect discontinuities scheduled to be repaired during LS 1, numerous other concerns regarding circuit stability arose. In particular, observations of an erratic behaviour of magnet bypass diodes and the degradation of other potentially weak busbar sections, as well as observations of seemingly random millisecond spikes in beam losses, known as unidentified falling object (UFO) events, which, if persist at 6.5 TeV, may eventually deposit sufficient energy to quench adjacent magnets. In light of the above, the thesis hypothesis states that, even with the observed issues, the LHC main dipole circuits can safely support and sustain near-nominal proton beam energies of at least 6.5 TeV. Research into minimising the risk of magnet training led to the development and implementation of a new qualification method, capable of providing conclusive evidence that all aspects of all circuits, other than the magnets and their internal joints, can safely withstand a quench event at near-nominal current levels, allowing for magnet training to be carried out both systematically and without risk. This method has become known as the Copper Stabiliser Continuity Measurement (CSCM). Results were a success, with all circuits eventually being subject to a full current decay from 6.5 TeV equivalent current levels, with no measurable damage occurring. Research into UFO events led to the development of a numerical model capable of simulating typical UFO events, reproducing entire Run 1 measured event data sets and extrapolating to 6.5 TeV, predicting the likelihood of UFO-induced magnet quenches. Results provided interesting insights into the involved phenomena as well as confirming the possibility of UFO-induced magnet quenches. The model was also capable of predicting that such events, if left unaccounted for, are likely to be commonplace or not, resulting in significant long-term issues for 6.5+ TeV operation. Addressing the thesis hypothesis, the following written works detail the development and results of all CSCM qualification tests and subsequent magnet training as well as the development and simulation results of both 4 TeV and 6.5 TeV UFO event modelling. The thesis concludes, post-LS 1, with the LHC successfully sustaining 6.5 TeV proton beams, but with UFO events, as predicted, resulting in otherwise uninitiated magnet quenches and being at the forefront of system availability issues.

Análise hidromecânica do problema de fraturamento hidráulico

Dissertação (mestrado)—Universidade de Brasília, Faculdade de Tecnologia, Departamento de Engenharia Civil e Ambiental, 2016.

Spectral-based mesh segmentation

In design and manufacturing, mesh segmentation is required for FACE construction in boundary representation (BRep), which in turn is central for featurebased design, machining, parametric CAD and reverse engineering, among others -- Although mesh segmentation is dictated by geometry and topology, this article focuses on the topological aspect (graph spectrum), as we consider that this tool has not been fully exploited -- We preprocess the mesh to obtain a edgelength homogeneous triangle set and its Graph Laplacian is calculated -- We then produce a monotonically increasing permutation of the Fiedler vector (2nd eigenvector of Graph Laplacian) for encoding the connectivity among part feature submeshes -- Within the mutated vector, discontinuities larger than a threshold (interactively set by a human) determine the partition of the original mesh -- We present tests of our method on large complex meshes, which show results which mostly adjust to BRep FACE partition -- The achieved segmentations properly locate most manufacturing features, although it requires human interaction to avoid over segmentation -- Future work includes an iterative application of this algorithm to progressively sever features of the mesh left from previous submesh removals

Comportamento magnético de compósitos de matriz polimérica com adição de ferritas

Polymer matrix composites offer advantages for many applications due their combination of properties, which includes low density, high specific strength and modulus of elasticity and corrosion resistance. However, the application of non-destructive techniques using magnetic sensors for the evaluation these materials is not possible since the materials are non-magnetizable. Ferrites are materials with excellent magnetic properties, chemical stability and corrosion resistance. Due to these properties, these materials are promising for the development of polymer composites with magnetic properties. In this work, glass fiber / epoxy circular plates were produced with 10 wt% of cobalt or barium ferrite particles. The cobalt ferrite was synthesized by the Pechini method. The commercial barium ferrite was subjected to a milling process to study the effect of particle size on the magnetic properties of the material. The characterization of the ferrites was carried out by x-ray diffraction (XRD), field emission gun scanning electron microscopy (FEG-SEM) and vibrating sample magnetometry (VSM). Circular notches of 1, 5 and 10 mm diameter were introduced in the composite plates using a drill bit for the non-destructive evaluation by the technique of magnetic flux leakage (MFL). The results indicated that the magnetic signals measured in plates with barium ferrite without milling and cobalt ferrite showed good correlation with the presence of notches. The milling process for 12 h and 20 h did not contribute to improve the identification of smaller size notches (1 mm). However, the smaller particle size produced smoother magnetic curves, with fewer discontinuities and improved signal-to-noise ratio. In summary, the results suggest that the proposed approach has great potential for the detection of damage in polymer composites structures

Espaços descontínuos

Dissertação para obtenção do grau de Mestre em Arquitectura com especialização em Urbanismo, apresentada na Universidade de Lisboa - Faculdade de Arquitetura.

Ontologia do espaço : critica da critica da entificação social do ser enquanto pressuposto a uma teoria espacial interpenetrada a “ontologia do ser social”, de Gyorgy Lukacs

Tese (doutorado)—Universidade de Brasília, Instituto de Ciências Humanas, Departamento de Geografia, Programa de Pós Graduação em Geografia, 2015.

2ª Circular

Dissertação para obtenção do grau de Mestre em Arquitectura com Especialização em Gestão Urbanística, apresentada na Universidade de Lisboa - Faculdade de Arquitectura.