Coupled Nonlinear Ginzburg-Landau and Mechanics Model for Martensitic Transformations in Polycrystals

Las transformaciones martensíticas (MT) se definen como un cambio en la estructura del cristal para formar una fase coherente o estructuras de dominio multivariante, a partir de la fase inicial con la misma composición, debido a pequeños intercambios o movimientos atómicos cooperativos. En el siglo pasado se han descubierto MT en diferentes materiales partiendo desde los aceros hasta las aleaciones con memoria de forma, materiales cerámicos y materiales inteligentes. Todos muestran propiedades destacables como alta resistencia mecánica, memoria de forma, efectos de superelasticidad o funcionalidades ferroicas como la piezoelectricidad, electro y magneto-estricción etc. Varios modelos/teorías se han desarrollado en sinergia con el desarrollo de la física del estado sólido para entender por qué las MT generan microstructuras muy variadas y ricas que muestran propiedades muy interesantes. Entre las teorías mejor aceptadas se encuentra la Teoría Fenomenológica de la Cristalografía Martensítica (PTMC, por sus siglas en inglés) que predice el plano de hábito y las relaciones de orientación entre la austenita y la martensita. La reinterpretación de la teoría PTMC en un entorno de mecánica del continuo (CM-PTMC) explica la formación de los dominios de estructuras multivariantes, mientras que la teoría de Landau con dinámica de inercia desentraña los mecanismos físicos de los precursores y otros comportamientos dinámicos. La dinámica de red cristalina desvela la reducción de la dureza acústica de las ondas de tensión de red que da lugar a transformaciones débiles de primer orden en el desplazamiento. A pesar de las diferencias entre las teorías estáticas y dinámicas dado su origen en diversas ramas de la física (por ejemplo mecánica continua o dinámica de la red cristalina), estas teorías deben estar inherentemente conectadas entre sí y mostrar ciertos elementos en común en una perspectiva unificada de la física. No obstante las conexiones físicas y diferencias entre las teorías/modelos no se han tratado hasta la fecha, aun siendo de importancia crítica para la mejora de modelos de MT y para el desarrollo integrado de modelos de transformaciones acopladas de desplazamiento-difusión. Por lo tanto, esta tesis comenzó con dos objetivos claros. El primero fue encontrar las conexiones físicas y las diferencias entre los modelos de MT mediante un análisis teórico detallado y simulaciones numéricas. El segundo objetivo fue expandir el modelo de Landau para ser capaz de estudiar MT en policristales, en el caso de transformaciones acopladas de desplazamiento-difusión, y en presencia de dislocaciones. Comenzando con un resumen de los antecedente, en este trabajo se presentan las bases físicas de los modelos actuales de MT. Su capacidad para predecir MT se clarifica mediante el ansis teórico y las simulaciones de la evolución microstructural de MT de cúbicoatetragonal y cúbicoatrigonal en 3D. Este análisis revela que el modelo de Landau con representación irreducible de la deformación transformada es equivalente a la teoría CM-PTMC y al modelo de microelasticidad para predecir los rasgos estáticos durante la MT, pero proporciona una mejor interpretación de los comportamientos dinámicos. Sin embargo, las aplicaciones del modelo de Landau en materiales estructurales están limitadas por su complejidad. Por tanto, el primer resultado de esta tesis es el desarrollo del modelo de Landau nolineal con representación irreducible de deformaciones y de la dinámica de inercia para policristales. La simulación demuestra que el modelo propuesto es consistente fcamente con el CM-PTMC en la descripción estática, y también permite una predicción del diagrama de fases con la clásica forma ’en C’ de los modos de nucleación martensítica activados por la combinación de temperaturas de enfriamiento y las condiciones de tensión aplicada correlacionadas con la transformación de energía de Landau. Posteriomente, el modelo de Landau de MT es integrado con un modelo de transformación de difusión cuantitativa para elucidar la relajación atómica y la difusión de corto alcance de los elementos durante la MT en acero. El modelo de transformaciones de desplazamiento y difusión incluye los efectos de la relajación en borde de grano para la nucleación heterogenea y la evolución espacio-temporal de potenciales de difusión y movilidades químicas mediante el acoplamiento de herramientas de cálculo y bases de datos termo-cinéticos de tipo CALPHAD. El modelo se aplica para estudiar la evolución microstructural de aceros al carbono policristalinos procesados por enfriamiento y partición (Q&P) en 2D. La microstructura y la composición obtenida mediante la simulación se comparan con los datos experimentales disponibles. Los resultados muestran el importante papel jugado por las diferencias en movilidad de difusión entre la fase austenita y martensita en la distibución de carbono en las aceros. Finalmente, un modelo multi-campo es propuesto mediante la incorporación del modelo de dislocación en grano-grueso al modelo desarrollado de Landau para incluir las diferencias morfológicas entre aceros y aleaciones con memoria de forma con la misma ruptura de simetría. La nucleación de dislocaciones, la formación de la martensita ’butterfly’, y la redistribución del carbono después del revenido son bien representadas en las simulaciones 2D del estudio de la evolución de la microstructura en aceros representativos. Con dicha simulación demostramos que incluyendo las dislocaciones obtenemos para dichos aceros, una buena comparación frente a los datos experimentales de la morfología de los bordes de macla, la existencia de austenita retenida dentro de la martensita, etc. Por tanto, basado en un modelo integral y en el desarrollo de códigos durante esta tesis, se ha creado una herramienta de modelización multiescala y multi-campo. Dicha herramienta acopla la termodinámica y la mecánica del continuo en la macroescala con la cinética de difusión y los modelos de campo de fase/Landau en la mesoescala, y también incluye los principios de la cristalografía y de la dinámica de red cristalina en la microescala. ABSTRACT Martensitic transformation (MT), in a narrow sense, is defined as the change of the crystal structure to form a coherent phase, or multi-variant domain structures out from a parent phase with the same composition, by small shuffles or co-operative movements of atoms. Over the past century, MTs have been discovered in different materials from steels to shape memory alloys, ceramics, and smart materials. They lead to remarkable properties such as high strength, shape memory/superelasticity effects or ferroic functionalities including piezoelectricity, electro- and magneto-striction, etc. Various theories/models have been developed, in synergy with development of solid state physics, to understand why MT can generate these rich microstructures and give rise to intriguing properties. Among the well-established theories, the Phenomenological Theory of Martensitic Crystallography (PTMC) is able to predict the habit plane and the orientation relationship between austenite and martensite. The re-interpretation of the PTMC theory within a continuum mechanics framework (CM-PTMC) explains the formation of the multivariant domain structures, while the Landau theory with inertial dynamics unravels the physical origins of precursors and other dynamic behaviors. The crystal lattice dynamics unveils the acoustic softening of the lattice strain waves leading to the weak first-order displacive transformation, etc. Though differing in statics or dynamics due to their origins in different branches of physics (e.g. continuum mechanics or crystal lattice dynamics), these theories should be inherently connected with each other and show certain elements in common within a unified perspective of physics. However, the physical connections and distinctions among the theories/models have not been addressed yet, although they are critical to further improving the models of MTs and to develop integrated models for more complex displacivediffusive coupled transformations. Therefore, this thesis started with two objectives. The first one was to reveal the physical connections and distinctions among the models of MT by means of detailed theoretical analyses and numerical simulations. The second objective was to expand the Landau model to be able to study MTs in polycrystals, in the case of displacive-diffusive coupled transformations, and in the presence of the dislocations. Starting with a comprehensive review, the physical kernels of the current models of MTs are presented. Their ability to predict MTs is clarified by means of theoretical analyses and simulations of the microstructure evolution of cubic-to-tetragonal and cubic-to-trigonal MTs in 3D. This analysis reveals that the Landau model with irreducible representation of the transformed strain is equivalent to the CM-PTMC theory and microelasticity model to predict the static features during MTs but provides better interpretation of the dynamic behaviors. However, the applications of the Landau model in structural materials are limited due its the complexity. Thus, the first result of this thesis is the development of a nonlinear Landau model with irreducible representation of strains and the inertial dynamics for polycrystals. The simulation demonstrates that the updated model is physically consistent with the CM-PTMC in statics, and also permits a prediction of a classical ’C shaped’ phase diagram of martensitic nucleation modes activated by the combination of quenching temperature and applied stress conditions interplaying with Landau transformation energy. Next, the Landau model of MT is further integrated with a quantitative diffusional transformation model to elucidate atomic relaxation and short range diffusion of elements during the MT in steel. The model for displacive-diffusive transformations includes the effects of grain boundary relaxation for heterogeneous nucleation and the spatio-temporal evolution of diffusion potentials and chemical mobility by means of coupling with a CALPHAD-type thermo-kinetic calculation engine and database. The model is applied to study for the microstructure evolution of polycrystalline carbon steels processed by the Quenching and Partitioning (Q&P) process in 2D. The simulated mixed microstructure and composition distribution are compared with available experimental data. The results show that the important role played by the differences in diffusion mobility between austenite and martensite to the partitioning in carbon steels. Finally, a multi-field model is proposed by incorporating the coarse-grained dislocation model to the developed Landau model to account for the morphological difference between steels and shape memory alloys with same symmetry breaking. The dislocation nucleation, the formation of the ’butterfly’ martensite, and the redistribution of carbon after tempering are well represented in the 2D simulations for the microstructure evolution of the representative steels. With the simulation, we demonstrate that the dislocations account for the experimental observation of rough twin boundaries, retained austenite within martensite, etc. in steels. Thus, based on the integrated model and the in-house codes developed in thesis, a preliminary multi-field, multiscale modeling tool is built up. The new tool couples thermodynamics and continuum mechanics at the macroscale with diffusion kinetics and phase field/Landau model at the mesoscale, and also includes the essentials of crystallography and crystal lattice dynamics at microscale.

Grain refinement efficiency and mechanism of aluminium carbide in Mg-Al alloys

Detailed microscopic examination using optical and electron microscopes suggests that Al4C3, often observed in the central regions of magnesium grains on polished sections, is a potent substrate for primary Mg. Calculations of the crystallographic relationships between magnesium and Al4C3 further support the experimental observations. (c) 2005 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Study of the calcification of PHEMA hydrogels using a two compartment permeation cell

To simulate the process of calcification in hydrogel implants, particularly calcification inside hydrogels, in vitro experiments using two compartment permeation cells have been performed. PHEMA hydrogel membranes were synthesized by free radical polymerization in bulk. The permeability and diffusion coefficient for Ca2+ ions at 37 ° C were determined using Fick's laws of diffusion. It was evident that Ca2+ ions either from CaCl2 or SBF solutions may diffuse through PHEMA hydrogel membranes. The fort-nation of calcium phosphate deposits inside the hydrogel was observed and attributed to a heterogeneous nucleation from diffusing calcium and phosphate ions. The morphology of the deposits both on the surface and inside the hydrogels was found to be similar, i.e. spherical aggregates with a diameter of less than one micron. © 2005 Elsevier B.V. All rights reserved.

Precipitation studies on alloys from the system copper chromium zirconium

The precipitation reactions occurring in a series of copper-based alloys selected from the system copper-chromium-zirconium have been studied by resistometric and metallographic techniques. A survey of the factors influencing the development of copper-based alloys for high strength, high conductivity applications is followed by a more general review of contemporary materials, and illustrates that the most promising alloys are those containing chromium and zirconium. The few systematic attempts to study alloys from this system have been collated, discussed, and used as a basis for the selection of four alloy compositions viz:- Cu - 0.4% Cr Cu - 0.24. Zr Cu - 0. 3% Cr - 0.1% Zr Cu - 0.2% Cr - 0.2% Zr A description of the experimental techniques used to study the precipitation behaviour of these materials is preceeded by a discussion of the currently accepted theories relating to precipitate nucleation and growth. The experimental results are presented and discussed for each of the alloys independently, and are then treated jointly to obtain an overall assessment of the way in which the precipitation kinetics, metallography and mechanical properties vary with alloy composition and heat treatment. The metastable solid solution of copper-chromium is found to decompose by the rejection of chromium particles which maintain a coherent interface and a Kurdjumov-Sachs type crystallographic orientation relationship with the copper matrix. The addition of 0.1% zirconium to the alloy retards the rate of transformation by a factor of ten and modifies the dispersion characteristics of the precipitate without markedly altering the morphology. Further additions of zirconium lead to the growth of stacking faults during ageing, which provide favourable nucleation sites for the chromium precipitate. The partial dislocations bounding such stacking faults are also found to provide mobile heterogeneous nucleation sources for the precipitation reactions occurring in copper-zirconium.

Comparison of the ageing behaviour of PM 2124 Al alloy and Al-SiCp metal-matrix composite

The ageing response of 2124 Al-SiC particulate metal-matrix composite (MMC) and unreinforced alloy has been examined using hardness measurements and Arrhenius analysis. The formation of phases during precipitation has been studied using differential scanning calorimetry (DSC). The MMC exhibits accelerated ageing compared to unreinforced alloy, due to enhanced S′ formation. The activation energy for diffusion is lower in the MMC than in the unreinforced alloy. DSC scans show Guinier-Preston B (GPB) zone nucleation to occur at a lower temperature in the MMC, whilst the total volume of GPB zones formed is smaller than in the unreinforced alloy. A model has been proposed to explain the GPB zone formation behaviour, in which ease of GPB zone nucleation varies within the MMC, as a function of ageing time and of position within the matrix. S′ formation is enhanced in the MMC because of improved diffusion and a large increase in density of heterogeneous nucleation sites compared to the unreinforced alloy. © 1994 Chapman & Hall.

Hydrophobic hydration driven self-assembly of curcumin in water: Similarities to nucleation and growth under large metastability, and an analysis of water dynamics at heterogeneous surfaces

As the beneficial effects of curcumin have often been reported to be limited to its small concentrations, we have undertaken a study to find the aggregation properties of curcumin in water by varying the number of monomers. Our molecular dynamics simulation results show that the equilibrated structure is always an aggregated state with remarkable structural rearrangements as we vary the number of curcumin monomers from 4 to 16 monomers. We find that the curcumin monomers form clusters in a very definite pattern where they tend to aggregate both in parallel and anti-parallel orientation of the phenyl rings, often seen in the formation of beta-sheet in proteins. A considerable enhancement in the population of parallel alignments is observed with increasing the system size from 12 to 16 curcumin monomers. Due to the prevalence of such parallel alignment for large system size, a more closely packed cluster is formed with maximum number of hydrophobic contacts. We also follow the pathway of cluster growth, in particular the transition from the initial segregated to the final aggregated state. We find the existence of a metastable structural intermediate involving a number of intermediate-sized clusters dispersed in the solution. We have constructed a free energy landscape of aggregation where the metatsable state has been identified. The course of aggregation bears similarity to nucleation and growth in highly metastable state. The final aggregated form remains stable with the total exclusion of water from its sequestered hydrophobic core. We also investigate water structure near the cluster surface along with their orientation. We find that water molecules form a distorted tetrahedral geometry in the 1st solvation layer of the cluster, interacting rather strongly with the hydrophilic groups at the surface of the curcumin. The dynamics of such quasi-bound water molecules near the surface of curcumin cluster is considerably slower than the bulk signifying a restricted motion as often found in protein hydration layer. (C) 2014 AIP Publishing LLC.

Dynamos driven by weak thermal convection and heterogeneous outer boundary heat flux

We use numerical dynamo models with heterogeneous core-mantle boundary (CMB) heat flux to show that lower mantle lateral thermal variability may help support a dynamo under weak thermal convection. In our reference models with homogeneous CMB heat flux, convection is either marginally supercritical or absent, always below the threshold for dynamo onset. We find that lateral CMB heat flux variations organize the flow in the core into patterns that favour the growth of an early magnetic field. Heat flux patterns symmetric about the equator produce non-reversing magnetic fields, whereas anti-symmetric patterns produce polarity reversals. Our results may explain the existence of the geodynamo prior to inner core nucleation under a tight energy budget. Furthermore, in order to sustain a strong geomagnetic field, the lower mantle thermal distribution was likely dominantly symmetric about the equator. (C) 2015 Elsevier B.V. All rights reserved.

A three-dimensional analytical solution of the microheater temperature before bubble nucleation

A three-dimensional analytical solution of the microheater temperature based on heat diffusion equation is developed and compared with experimental results. Dimensionless parameters are introduced to analyze the temperature rise time and the distribution under steady state. To study the microheater temperatures before bubble nucleation, a set of working fluids and microheaters are considered. It is shown that the dimensionless time xi(-)(0) required for the temperature rise from room to 95% of the steady state temperature is about 75, not dependent on working fluids and microheaters. Heat transfer to the surrounding liquid is mainly caused by conduction, not by convection and radiation mechanisms. The microheater length affects the surface temperature uniformity, while its width influences the steady temperatures significantly, yielding the transition from heterogeneous to homogeneous nucleation mechanism from square microheaters to narrow line microheaters. 

Observation of Nucleation and Growth of CdS Nanocrystals in a Two-phase System

The nucleation and growth kinetics of CdS nanocrystals in a two-phase synthesis system have been investigated. It was found that the nucleation process is quite lengthy and overlapped with the growth process; nevertheless, as formed nanocrystals show extremely narrow size distribution owing to the unique heterogeneous reacting environment and Ostwald ripening growth. The nucleation and growth kinetics of the nanocrystals were also influenced strongly by the monomer concentration, capping agent concentration, and solvent polarity. It was also found that a high monomer concentration, a low capping agent concentration, and low solvent polarity lead to a higher maximum nucleus concentration and nanocrystal concentration, while high polarity solvents are favorable for the formation of nanocrystals with narrower size distribution and higher photoluminescence quantum yield.

Analytical model for heterogeneous crystallization kinetics of spherical glass particles

An analytical model developed to describe the crystallization kinetics of spherical glass particles has been derived in this work. A continuous phase transition from three-dimensional (3D)-like to 1D-like crystal growth has been considered and a procedure for the quantitative evaluation of the critical time for this 3D-1D transition is proposed. This model also allows straightforward determination of the density of surface nucleation sites on glass powders using differential scanning calorimetry data obtained under different thermal conditions. © 2009 The American Ceramic Society.

On the nucleation and growth of amyloid beta-protein fibrils: detection of nuclei and quantitation of rate constants.

We have studied the fibrillogenesis of synthetic amyloid beta-protein-(1-40) fragment (A beta) in 0.1 M HCl. At low pH, A beta formed fibrils at a rate amenable to detailed monitoring by quasi-elastic light-scattering spectroscopy. Examination of the fibrils with circular dichroism spectroscopy and electron microscopy showed them to be highly similar to those found in amyloid plaques. We determined the hydrodynamic radii of A beta aggregates during the entire process of fibril nucleation and growth. Above an A beta concentration of approximately 0.1 mM, the initial rate of elongation and the final size of fibrils were independent of A beta concentration. Below an A beta concentration of 0.1 mM, the initial elongation rate was proportional to the peptide concentration, and the resulting fibrils were significantly longer than those formed at higher concentration. We also found that the surfactant n-dodecylhexaoxyethylene glycol monoether (C12E6) slowed nucleation and elongation of fibrils in a concentration-dependent manner. Our observations are consistent with a model of A beta fibrillogenesis that includes the following key steps: (i) peptide micelles form above a certain critical A beta concentration, (ii) fibrils nucleate within these micelles or on heterogeneous nuclei (seeds), and (iii) fibrils grow by irreversible binding of monomers to fibril ends. Interpretation of our data enabled us to determine the sizes of fibril nuclei and A beta micelles and the rates of fibril nucleation (from micelles) and fibril elongation. Our approach provides a powerful means for the quantitative assay of A beta fibrillogenesis.