Application of the N-alkane molecular alloys to thermally protected containers for catering

A thermally controlled transport device was designed and tested. As hot food needs to be transported at temperatures between 60 and 70ºC in order to avoid contamination by microorganisms, the use of Molecular Alloy Phase Change Materials (MAPCM) can lead to improvements in this field of application. A heat transfer numerical simulation of the box used for transporting the food was conducted. Despite obvious simplifications, a good agreement between numerical simulation and experimental results was obtained. Furthermore, we compared our experimental results with those from other experiments related to the transport of hot food. Here, pizza is taken as the example, and it is shown that delivering time can be increased three-fold.

Puumassapesurin suotautumisolojen kartoitus ja simulointi

Työn tavoitteena oli tutkia sellumassan suotautumiseen liittyviä tekijöitä Ahlstrom Machineryn Drum Displacer™ -puumassapesurissa (DD-pesuri). Teoriaosassa tarkasteltiin aluksi suotautumisen teoriaa kuitu-vesi-suspensiossa, minkä jälkeen esiteltiin suotautumisnopeuteen vaikuttavia fysikaalisia ja kemiallisia vaikutusmekanismeja. Seuraavaksi kuvattiin massan pesun yleisiä perusteita sekä teoriaa puumassapesureissa. Lopuksi tarkasteltiin pesurien kytkeytymistä muuhun kuitulinjaan sekä prosessista johtuvia pesun toiminnan ulkoisia häiriötekijöitä. Kokeellisen osan aluksi tarkasteltiin paine- ja lämpötilamittauksien avulla massapesurissa vallitsevia prosessioloja. Mittaustulosten perusteella pumppausolot pesurin suodoslinjoissa ovat vaikeahkot ja häiriötilanteita voi esiintyä, mutta käytäntö on osoittanut tästä olevan vain harvoin haittaa prosessin toiminnalle. Pesureissa toteutuneet syrjäytysnopeudet laskettiin ja niitä verrattiin syrjäytystestien antamiin tuloksiin. Kuitulinjasta riippuen testin vastaavuus tehdasprosessiin vaihteli suuresti. Syrjäytystesteillä kokeiltiin myös tehdasprosesseissa usein esiintyvien muuttujien vaikutusta sellukakun syrjäytettävyyteen. Kakun paksuus ja syrjäytyslämpötila vaikuttivat syrjäytysnopeuteen Darcyn lain mukaisesti. Alipaineen massakakun alapuolella havaittiin huonontavan syrjäytysnopeutta verrattuna tilanteeseen, jossa kakun alla vallitsi ilmanpaine. Tämä havainto on selvästi ristiriidassa suotautumisen teorian kanssa. Massakakun muodostumis-pH osoittautui ratkaisevaksi lopulliselle syrjäytysnopeudelle, sillä alkalisissa oloissa muodostetun kuitukakun syrjäytysnopeus ei enää parantunut happamalla syrjäytysnesteellä. Happamissa oloissa muodostetun kakun syrjäytysnopeus oli alkalisista parempi, mutta se alkoi hitaasti alentua, kun syrjäytysneste vaihtui alkaliseen. Massan laimentaminen ennen syrjäytystä alkalisella tehdassuodoksella puhtaan veden sijasta alensi ligniinipitoisella massalla lopullista syrjäytysnopeutta. Shirato-Tillerin mallilla ja Jönssonin staattisella mallilla simuloitiin numeerisesti syrjäytystestiä kahdessa eri pH:ssa, ja simulointituloksia verrattiin vastaavissa oloissa tehtyihin syrjäytystesteihin. Shirato-Tillerin mallin antamien syrjäytysnopeuksien havaittiin olevan lähellä syrjäytystestien nopeuksia, kun Jönssonin mallin antamat tulokset jäivät huomattavasti testituloksia alemmiksi. Herkkyystarkastelussa havaittiin mallien olevan varsin herkkiä parametrien virheille. Hajonta vaadittavien kuituparametrien määrityksissä ja menetelmien työläys rajoittavat numeerisen simuloinnin käytettävyyttä, sillä kuituparametrien määrityksen vaatima työmäärä on ainakin toistaiseksi syrjäytystestiä suurempi. Lopuksi todettiin, että oikeiden syrjäytysolosuhteiden käyttö on ensiarvoisen tärkeää oikeiden tulosten saamiseksi sekä kokeellisessa että numeerisessa simuloinnissa. Nykyinen syrjäytystestilaitteisto on pienin muutoksin käyttökelpoinen, kun massan testaus prosessioloissa tulee rutiininomaiseksi.

Synthesis of highly focused fields with circular polarization at any transverse plane

We develop a method for generating focused vector beams with circular polarization at any transverse plane. Based on the Richards-Wolf vector model, we derive analytical expressions to describe the propagation of these set of beams near the focal area. Since the polarization and the amplitude of the input beam are not uniform, an interferometric system capable of generating spatially-variant polarized beams has to be used. In particular, this wavefront is manipulated by means of spatial light modulators displaying computer generated holograms and subsequently focused using a high numerical aperture objective lens. Experimental results using a NA=0.85 system are provided: irradiance and Stokes images of the focused field at different planes near the focal plane are presented and compared with those obtained by numerical simulation.

Propagação da poluição atmosférica por flúor nas águas subterrâneas e solos de regiões próximas às indústrias de fertilizantes (Rio Grande, RS)

Fluoride concentration was determined in rainwater, ground water and soil in the zone of fertilizer industry in the city of Rio Grande. In contaminated rainwater fluorine concentration was registered up to the value of 4,4 mg.L-1. Fluorine concentration in the shallow ground water in general reflects its distribution in the atmosphere, but cannot be used as marker of atmosphere contamination in the urban area due to dissolving influence of residential effluents. The 0,01% HCl extracts from the set of surface soil samples demonstrates fluoride distribution in the zone of influence of industrial emissions, which coincides to the numerical simulation of fluorides dispersion in the air.

Chaotic dynamics of electric-field domains in periodically driven superlattices

Self-sustained time-dependent current oscillations under dc voltage bias have been observed in recent experiments on n-doped semiconductor superlattices with sequential resonant tunneling. The current oscillations are caused by the motion and recycling of the domain wall separating low- and high-electric-field regions of the superlattice, as the analysis of a discrete drift model shows and experimental evidence supports. Numerical simulation shows that different nonlinear dynamical regimes of the domain wall appear when an external microwave signal is superimposed on the dc bias and its driving frequency and driving amplitude vary. On the frequency-amplitude parameter plane, there are regions of entrainment and quasiperiodicity forming Arnold tongues. Chaos is demonstrated to appear at the boundaries of the tongues and in the regions where they overlap. Coexistence of up to four electric-field domains randomly nucleated in space is detected under ac+dc driving.

Application of a non isotropic turbulence model to stable atmospheric flows and dispersion over 3D topography

A non isotropic turbulence model is extended and applied to three dimensional stably stratified flows and dispersion calculations. The model is derived from the algebraic stress model (including wall proximity effects), but it retains the simplicity of the "eddy viscosity" concept of first order models. The "modified k-epsilon" is implemented in a three dimensional numerical code. Once the flow is resolved, the predicted velocity and turbulence fields are interpolated into a second grid and used to solve the concentration equation. To evaluate the model, various steady state numerical solutions are compared with small scale dispersion experiments which were conducted at the wind tunnel of Mitsubishi Heavy Industries, in Japan. Stably stratified flows and plume dispersion over three distinct idealized complex topographies (flat and hilly terrain) are studied. Vertical profiles of velocity and pollutant concentration are shown and discussed. Also, comparisons are made against the results obtained with the standard k-epsilon model.

Thermalhydraulic behavior of electrically heated rod during a critical heat flux transient

In nuclear reactors, the occurrence of critical heat flux leads to fuel rod overheating with clad fusion and radioactive products leakage. To predict the effects of such phenomenon, experiments are performed using electrically heated rods to simulate operational and accidental conditions of nuclear fuel rods. In the present work, it is performed a theoretical analysis of the drying and rewetting front propagation during a critical heat flux experiment, starting with the application of an electrical power step from steady state condition. After the occurrence of critical heat flux, the drying front propagation is predicted. After a few seconds, a power cut is considered and the rewetting front behavior is analytically observed. Studies performed with various values of coolant mass flow rate show that this variable has more influence on the drying front velocity than on the rewetting one.

Surface Tension Implementation for Gensmac2d

In the present work we describe a method which allows the incorporation of surface tension into the GENSMAC2D code. This is achieved on two scales. First on the scale of a cell, the surface tension effects are incorporated into the free surface boundary conditions through the computation of the capillary pressure. The required curvature is estimated by fitting a least square circle to the free surface using the tracking particles in the cell and in its close neighbors. On a sub-cell scale, short wavelength perturbations are filtered out using a local 4-point stencil which is mass conservative. An efficient implementation is obtained through a dual representation of the cell data, using both a matrix representation, for ease at identifying neighbouring cells, and also a tree data structure, which permits the representation of specific groups of cells with additional information pertaining to that group. The resulting code is shown to be robust, and to produce accurate results when compared with exact solutions of selected fluid dynamic problems involving surface tension.

Dynamic Analysis Model of Spherical Roller Bearings with Defects

Rolling element bearings are essential components of rotating machinery. The spherical roller bearing (SRB) is one variant seeing increasing use, because it is self-aligning and can support high loads. It is becoming increasingly important to understand how the SRB responds dynamically under a variety of conditions. This doctoral dissertation introduces a computationally efficient, three-degree-of-freedom, SRB model that was developed to predict the transient dynamic behaviors of a rotor-SRB system. In the model, bearing forces and deflections were calculated as a function of contact deformation and bearing geometry parameters according to nonlinear Hertzian contact theory. The results reveal how some of the more important parameters; such as diametral clearance, the number of rollers, and osculation number; influence ultimate bearing performance. Distributed defects, such as the waviness of the inner and outer ring, and localized defects, such as inner and outer ring defects, are taken into consideration in the proposed model. Simulation results were verified with results obtained by applying the formula for the spherical roller bearing radial deflection and the commercial bearing analysis software. Following model verification, a numerical simulation was carried out successfully for a full rotor-bearing system to demonstrate the application of this newly developed SRB model in a typical real world analysis. Accuracy of the model was verified by comparing measured to predicted behaviors for equivalent systems.

SEPARATION OF SATURED AND UNSATURATED FATTY ACIDS FROM PALM FATTY ACIDS DISTILLATES IN CONTINUOUS MULTISTAGE COUNTERCURRENT COLUMNS WITH SUPERCRITICAL CARBON DIOXIDE AS SOLVENT: A PROCESS DESIGN METHODOLOGY

In this work the separation of multicomponent mixtures in counter-current columns with supercritical carbon dioxide has been investigated using a process design methodology. First the separation task must be defined, then phase equilibria experiments are carried out, and the data obtained are correlated with thermodynamic models or empirical functions. Mutual solubilities, Ki-values, and separation factors aij are determined. Based on this data possible operating conditions for further extraction experiments can be determined. Separation analysis using graphical methods are performed to optimize the process parameters. Hydrodynamic experiments are carried out to determine the flow capacity diagram. Extraction experiments in laboratory scale are planned and carried out in order to determine HETP values, to validate the simulation results, and to provide new materials for additional phase equilibria experiments, needed to determine the dependence of separation factors on concetration. Numerical simulation of the separation process and auxiliary systems is carried out to optimize the number of stages, solvent-to-feed ratio, product purity, yield, and energy consumption. Scale-up and cost analysis close the process design. The separation of palmitic acid and (oleic+linoleic) acids from PFAD-Palm Fatty Acids Distillates was used as a case study.

Interaction between vortices and impurities in a d-wave superconductor /

High temperature superconductors were discovered in 1986, but despite considerable research efforts, both experimental and theoretical, these materials remain poorly understood. Because their electronic structure is both inhomogeneous and highly correlated, a full understanding will require knowledge of quasiparticle properties both in real space and momentum space. In this thesis, we will present a theoretical analysis of the scanning tunneling microscopy (STM) data in BSCCO. We introduce the Bogoliubov-De Gennes Hamiltonian and solve it numerically on a two-dimensional 20 x 20 lattice under a magnetic field perpendicular to the surface. We consider a vortex at the center of our model. We introduce a Zn impurity in our lattice as a microscopic probe of the physical properties of BSCCO. By direct numerical diagonalization of the lattice BogoliubovDe Gennes Hamiltonian for different positions of the impurity, we can calculate the interaction between the vortex and the impurity in a d-wave superconductor.

Modélisation MHD tridimensionnelle de tubes de flux coronaux utilisant l'assimilation des donnés 4D-VAR

Mémoire numérisé par la Division de la gestion de documents et des archives de l'Université de Montréal

Simulations magnétohydrodynamiques en régime idéal

Cette thèse s’intéresse à la modélisation magnétohydrodynamique des écoulements de fluides conducteurs d’électricité multi-échelles en mettant l’emphase sur deux applications particulières de la physique solaire: la modélisation des mécanismes des variations de l’irradiance via la simulation de la dynamo globale et la reconnexion magnétique. Les variations de l’irradiance sur les périodes des jours, des mois et du cycle solaire de 11 ans sont très bien expliquées par le passage des régions actives à la surface du Soleil. Cependant, l’origine ultime des variations se déroulant sur les périodes décadales et multi-décadales demeure un sujet controversé. En particulier, une certaine école de pensée affirme qu’une partie de ces variations à long-terme doit provenir d’une modulation de la structure thermodynamique globale de l’étoile, et que les seuls effets de surface sont incapables d’expliquer la totalité des fluctuations. Nous présentons une simulation globale de la convection solaire produisant un cycle magnétique similaire en plusieurs aspects à celui du Soleil, dans laquelle le flux thermique convectif varie en phase avec l’ ́energie magnétique. La corrélation positive entre le flux convectif et l’énergie magnétique supporte donc l’idée qu’une modulation de la structure thermodynamique puisse contribuer aux variations à long-terme de l’irradiance. Nous analysons cette simulation dans le but d’identifier le mécanisme physique responsable de la corrélation en question et pour prédire de potentiels effets observationnels résultant de la modulation structurelle. La reconnexion magnétique est au coeur du mécanisme de plusieurs phénomènes de la physique solaire dont les éruptions et les éjections de masse, et pourrait expliquer les températures extrêmes caractérisant la couronne. Une correction aux trajectoires du schéma semi-Lagrangien classique est présentée, qui est basée sur la solution à une équation aux dérivées partielles nonlinéaire du second ordre: l’équation de Monge-Ampère. Celle-ci prévient l’intersection des trajectoires et assure la stabilité numérique des simulations de reconnexion magnétique pour un cas de magnéto-fluide relaxant vers un état d’équilibre.

Patient-specific anisotropic model of human trunk based on MR data

There are many ways to generate geometrical models for numerical simulation, and most of them start with a segmentation step to extract the boundaries of the regions of interest. This paper presents an algorithm to generate a patient-specific three-dimensional geometric model, based on a tetrahedral mesh, without an initial extraction of contours from the volumetric data. Using the information directly available in the data, such as gray levels, we built a metric to drive a mesh adaptation process. The metric is used to specify the size and orientation of the tetrahedral elements everywhere in the mesh. Our method, which produces anisotropic meshes, gives good results with synthetic and real MRI data. The resulting model quality has been evaluated qualitatively and quantitatively by comparing it with an analytical solution and with a segmentation made by an expert. Results show that our method gives, in 90% of the cases, as good or better meshes as a similar isotropic method, based on the accuracy of the volume reconstruction for a given mesh size. Moreover, a comparison of the Hausdorff distances between adapted meshes of both methods and ground-truth volumes shows that our method decreases reconstruction errors faster. Copyright © 2015 John Wiley & Sons, Ltd.

Studies on Some Aspects of Light Beam Propagation Through Certain Nonlinear Media

This thesis deals with the study of light beam propagation through different nonlinear media. Analytical and numerical methods are used to show the formation of solitonS in these media. Basic experiments have also been performed to show the formation of a self-written waveguide in a photopolymer. The variational method is used for the analytical analysis throughout the thesis. Numerical method based on the finite-difference forms of the original partial differential equation is used for the numerical analysis.In Chapter 2, we have studied two kinds of solitons, the (2 + 1) D spatial solitons and the (3 + l)D spatio-temporal solitons in a cubic-quintic medium in the presence of multiphoton ionization.In Chapter 3, we have studied the evolution of light beam through a different kind of nonlinear media, the photorcfractive polymer. We study modulational instability and beam propagation through a photorefractive polymer in the presence of absorption losses. The one dimensional beam propagation through the nonlinear medium is studied using variational and numerical methods. Stable soliton propagation is observed both analytically and numerically.Chapter 4 deals with the study of modulational instability in a photorefractive crystal in the presence of wave mixing effects. Modulational instability in a photorefractive medium is studied in the presence of two wave mixing. We then propose and derive a model for forward four wave mixing in the photorefractive medium and investigate the modulational instability induced by four wave mixing effects. By using the standard linear stability analysis the instability gain is obtained.Chapter 5 deals with the study of self-written waveguides. Besides the usual analytical analysis, basic experiments were done showing the formation of self-written waveguide in a photopolymer system. The formation of a directional coupler in a photopolymer system is studied theoretically in Chapter 6. We propose and study, using the variational approximation as well as numerical simulation, the evolution of a probe beam through a directional coupler formed in a photopolymer system.