The structural and piezoresponse properties of c-axis-oriented Aurivillius phase Bi5Ti3FeO15 thin films deposited by atomic vapor deposition

The deposition by atomic vapor deposition of highly c-axis-oriented Aurivillius phase Bi 5Ti 3FeO 15 (BTFO) thin films on (100) Si substrates is reported. Partially crystallized BTFO films with c-axis perpendicular to the substrate surface were first deposited at 610°C (8 excess Bi), and subsequently annealed at 820°C to get stoichiometric composition. After annealing, the films were highly c-axis-oriented, showing only (00l) peaks in x-ray diffraction (XRD), up to (0024). Transmission electron microscopy (TEM) confirms the BTFO film has a clear layered structure, and the bismuth oxide layer interleaves the four-block pseudoperovskite layer, indicating the n 4 Aurivillius phase structure. Piezoresponse force microscopy measurements indicate strong in-plane piezoelectric response, consistent with the c-axis layered structure, shown by XRD and TEM.

Room temperature ferroelectric and magnetic investigations and detailed phase analysis of Aurivillius phase Bi 5Ti 3Fe 0.7Co 0.3O 15 thin films

Aurivillius phase Bi 5Ti 3Fe 0.7Co 0.3O 15 (BTF7C3O) thin films on α-quartz substrates were fabricated by a chemical solution deposition method and the room temperature ferroelectric and magnetic properties of this candidate multiferroic were compared with those of thin films of Mn 3 substituted, Bi 5Ti 3Fe 0.7Mn 0.3O 15 (BTF7M3O). Vertical and lateral piezoresponse force microscopy (PFM) measurements of the films conclusively demonstrate that BTF7C3O and BTF7M3O thin films are piezoelectric and ferroelectric at room temperature, with the major polarization vector in the lateral plane of the films. No net magnetization was observed for the in-plane superconducting quantum interference device (SQUID) magnetometry measurements of BTF7M3O thin films. In contrast, SQUID measurements of the BTF7C3O films clearly demonstrated ferromagnetic behavior, with a remanent magnetization, B r, of 6.37 emu/cm 3 (or 804 memu/g), remanent moment 4.99 × 10 -5 emu. The BTF7C3O films were scrutinized by x-ray diffraction, high resolution transmission electron microscopy, scanning transmission electron microscopy, and energy dispersive x-ray analysis mapping to assess the prospect of the observed multiferroic properties being intrinsic to the main phase. The results of extensive micro-structural phase analysis demonstrated that the BTF7C3O films comprised of a 3.95 Fe/Co-rich spinel phase, likely CoFe 2 - xTi xO 4, which would account for the observed magnetic moment in the films. Additionally, x-ray magnetic circular dichroism photoemission electron microscopy (XMCD-PEEM) imaging confirmed that the majority of magnetic response arises from the Fe sites of Fe/Co-rich spinel phase inclusions. While the magnetic contribution from the main phase could not be determined by the XMCD-PEEM images, these data however imply that the Bi 5Ti 3Fe 0.7Co 0.3O 15 thin films are likely not single phase multiferroics at room temperature. The PFM results presented demonstrate that the naturally 2D nanostructured Bi 5Ti 3Fe 0.7Co 0.3O 15 phase is a novel ferroelectric and has potential commercial applications in high temperature piezoelectric and ferroelectric memory technologies. The implications for the conclusive demonstration of ferroelectric and ferromagnetic properties in single-phase materials of this type are discussed.

Magnetic field-induced ferroelectric switching in multiferroic aurivillius phase thin films at room temperature

Single-phase multiferroic materials are of considerable interest for future memory and sensing applications. Thin films of Aurivillius phase Bi 7Ti3Fe3O21 and Bi6Ti 2.8Fe1.52Mn0.68O18 (possessing six and five perovskite units per half-cell, respectively) have been prepared by chemical solution deposition on c-plane sapphire. Superconducting quantum interference device magnetometry reveal Bi7Ti3Fe 3O21 to be antiferromagnetic (TN = 190 K) and weakly ferromagnetic below 35 K, however, Bi6Ti2.8Fe 1.52Mn0.68O18 gives a distinct room-temperature in-plane ferromagnetic signature (Ms = 0.74 emu/g, μ0Hc =7 mT). Microstructural analysis, coupled with the use of a statistical analysis of the data, allows us to conclude that ferromagnetism does not originate from second phase inclusions, with a confidence level of 99.5%. Piezoresponse force microscopy (PFM) demonstrates room-temperature ferroelectricity in both films, whereas PFM observations on Bi6Ti2.8Fe1.52Mn0.68O18 show Aurivillius grains undergo ferroelectric domain polarization switching induced by an applied magnetic field. Here, we show for the first time that Bi6Ti2.8Fe1.52Mn0.68O18 thin films are both ferroelectric and ferromagnetic and, demonstrate magnetic field-induced switching of ferroelectric polarization in individual Aurivillius phase grains at room temperature.

GeoGebraTUTOR : développement d’un système tutoriel autonome pour l’accompagnement d’élèves en situation de résolution de problèmes de démonstration en géométrie plane et genèse d’un espace de travail géométrique idoine

Travaux d'études doctorales réalisées conjointement avec les travaux de recherches doctorales de Nicolas Leduc, étudiant au doctorat en génie informatique à l'École Polytechnique de Montréal.

Ultra Compact Inline E-Plane Waveguide Bandpass Filters Using Cross Coupling

This paper presents novel ultra-compact waveguide bandpass filters that exhibit pseudo elliptic responses with ability to place transmission zeros on both sides of the passband to form sharp roll offs. The filters contain E plane extracted pole sections cascaded with cross-coupled filtering blocks. Compactness is achieved by the use of evanescent mode sections and closer arranged resonators modified to shrink in size. The filters containing non-resonating nodes are designed by means of the generalized coupling coefficients (GCC) extraction procedure for the cross-coupled filtering blocks and extracted pole sections. We illustrate the performance of the proposed structures through the design examples of a third and a fourth order filters with center frequencies of 9.2 GHz and 10 GHz respectively. The sizes of the proposed structures suitable for fabricating using the low cost E plane waveguide technology are 38% smaller than ones of the E plane extracted pole filter of the same order.

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 influence of nano-scale second-phase particles on deformation of fine grained calcite mylonites

Grey and white carbonate mylonites were collected along thrust planes of the Helvetic Alps. They are characterised by very small grain sizes and non-random grain shape (SPO) and crystallographic preferred orientation (CPO). Presumably they deformed in the field of grain size sensitive flow by recrystallisation accommodated intracrystalline deformation in combination with granular flow. Both mylonites show a similar mean grain size, but in the grey mylonites the grain size range is larger, the grain shapes are more elongate and the dynamically recrystallised calcite grains are more often twinned. Grey mylonites have an oblique CPO, while the CPO in white mylonites is symmetric with respect to the shear plane. Combustion analysis and TEM investigations revealed that grey mylonites contain a higher amount of highly structured kerogens with particle sizes of a few tens of nanometers, which are finely dispersed at the grain boundaries. During deformation of the rock, nano-scale particles reduced the migration velocity of grain boundaries by Zener drag resulting in slower recrystallisation rates of the calcite aggregate. In the grey mylonites, more strain increments were accommodated by individual grains before they became refreshed by dynamic recrystallisation than in white mylonites, where grain boundary migration was less hindered and recrystallisation cycles were faster. Consequently, grey mylonites represent ‘deformation’ microfabrics while white mylonites are characterised by ‘recrystallisation’ microfabrics. Field geologists must utilise this different deformation behavior when applying the obliquity in CPO and SPO of the respective mylonites as reliable shear sense indicators.

GeoGebraTUTOR : développement d’un système tutoriel autonome pour l’accompagnement d’élèves en situation de résolution de problèmes de démonstration en géométrie plane et genèse d’un espace de travail géométrique idoine

Travaux d'études doctorales réalisées conjointement avec les travaux de recherches doctorales de Nicolas Leduc, étudiant au doctorat en génie informatique à l'École Polytechnique de Montréal.

The distribution of grain boundary and interface plane orientations in transformed microstructures

The properties of interfaces depend not only on the lattice misorientation, but also on the interface plane orientation. Extensive studies of grain boundaries led to the conclusion that in systems evolving by grain growth, the relative areas of different grain boundary planes are inversely correlated to their relative energies. In other words, the low energy grain boundary planes make up a larger part of the population than the higher energy grain boundary planes. The hypothesis of this work is that the interface plane orientation distribution in transformed microstructures depends more on the mechanism of formation than on the relative energy. After a discussion of methods for measuring interface plane orientations, results will be presented for lath martensite in a low carbon steel and for martensite in a Ti-6Al-4V alloy processed in two different ways to promote a displacive transformation in one case and a diffusional transformation in the other.

The role of thermomechanical routes on the distribution of grain boundary and interface plane orientations in transformed microstructures

In the current study, a series of thermomechanical routes were used to produce different microstructures (i.e., ferrite and martensite) in low-carbon low alloy steels. The five-parameter grain boundary character distribution was measured for all microstructures. The thermomechanical processing route altered the texture of the fully ferritic microstructure and significantly influenced the anisotropy of the grain boundary character distribution. Generally, the population of (111) planes increased with an increase in the γ-fiber texture for the ferritic microstructure, but it did not change the shape of the grain boundary plane distribution at specific misorientations. The most commonly observed boundaries in the fully ferritic structures produced through different routes were {112} symmetric tilt boundaries with the Σ3 = 60 deg/[111] misorientation; this boundary also had a low energy. However, the grain boundary plane distribution was significantly changed by the phase transformation path (i.e., ferrite vs martensite) for a given misorientation. In the martensitic steel, the most populous Σ3 boundary was the {110} symmetric tilt boundary. This results from the crystallographic constraints associated with the shear transformation (i.e., martensite) rather than the low-energy interface that dominates in the diffusional phase transformation (i.e., ferrite).

Liquid phase adsorption kinetics and equilibrium of toluene by novel modified-diatomite

The adsorption equilibria of toluene from aqueous solutions on natural and modified diatomite were examined at different operation parameters such as pH, contact time, initial toluene concentration was evaluated and optimum experimental conditions were identified. The surface area and morphology of the nanoparticles were characterized by SEM, BET, XRD, FTIR and EDX analysis. It was found that in order to obtain the highest possible removal of toluene, the experiments can be carried out at pH 6, temperature 25°C, an agitation speed of 200 rpm, an initial toluene concentration of 150 mg/L, a centrifugal rate of 4000 rpm, adsorbent dosage = 0.1 g and a process time of 90 min. The results of this work show that the maximum percentage removal of toluene from aqueous solution in the optimum conditions for NONMD was 96.91% (145.36 mg/g). Furthermore, under same conditions, the maximum adsorption of natural diatomite was 71.45% (107.18 mg/g). Both adsorption kinetic and isotherm experiments were carried out. The experimental data showed that the adsorption follows the Langmuir model and Freundlich model on natural and modified diatomite respectively. The kinetics results were found to conform well to pseudo-second order kinetics model with good correlation. Thus, this study demonstrated that the modified diatomite could be used as potential adsorbent for removal of toluene from aqueous solution.

Direct visualization of magnetic-field-induced magnetoelectric switching in multiferroic Aurivillius phase thin films

Multiferroic materials displaying coupled ferroelectric and ferromagnetic order parameters could provide a means for data storage whereby bits could be written electrically and read magnetically, or vice versa. Thin films of Aurivillius phase Bi6Ti2.8Fe1.52Mn0.68O18, previously prepared by a chemical solution deposition (CSD) technique, are multiferroics demonstrating magnetoelectric coupling at room temperature. Here, we demonstrate the growth of a similar composition, Bi6Ti2.99Fe1.46Mn0.55O18, via the liquid injection chemical vapor deposition technique. High-resolution magnetic measurements reveal a considerably higher in-plane ferromagnetic signature than CSD grown films (MS = 24.25 emu/g (215 emu/cm3), MR = 9.916 emu/g (81.5 emu/cm3), HC = 170 Oe). A statistical analysis of the results from a thorough microstructural examination of the samples, allows us to conclude that the ferromagnetic signature can be attributed to the Aurivillius phase, with a confidence level of 99.95%. In addition, we report the direct piezoresponse force microscopy visualization of ferroelectric switching while going through a full in-plane magnetic field cycle, where increased volumes (8.6 to 14% compared with 4 to 7% for the CSD-grown films) of the film engage in magnetoelectric coupling and demonstrate both irreversible and reversible magnetoelectric domain switching.

Condensed phase combustion travelling waves with sequential exothermic or endothermic reactions

The one-dimensional propagation of a combustion wave through a premixed solid fuel for two-stage kinetics is studied. We re-examine the analysis of a single reaction travelling-wave and extend it to the case of two-stage reactions. We derive an expression for the travelling wave speed in the limit of large activation energy for both reactions. The analysis shows that when both reactions are exothermic, the wave structure is similar to the single reaction case. However, when the second reaction is endothermic, the wave structure can be significantly different from single reaction case. In particular, as might be expected, a travelling wave does not necessarily exist in this case. We establish conditions in the limiting large activation energy limit for the non-existence, and for monotonicity of the temperature profile in the travelling wave.