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.

Estudo da solidificação equiaxial utilizando o modelo do campo de fases tridimensional.

Este trabalho apresenta um estudo da solidificação de metais puros utilizando o modelo de campo de fases. O modelo é utilizado para simular a solidificação com o intuito de obter a morfologia da interface sólido-líquido sob diversas condições de transferência de calor. Foram realizados testes de validação comparando as morfologias da interface sólido-líquido obtida com as morfologias apresentadas em trabalhos anteriores para os casos bi e tridimensionais. O modelo do campo de fases adotado consiste principalmente de duas equações diferenciais: uma para calcular a variável de campo de fases e outra para calcular o campo de temperaturas. As equações foram solucionadas numericamente para um oitavo do domínio devido a simetria do problema. Os cálculos do modelo indicam que um sólido esférico com um raio inicial menor que o raio crítico de nucleação refunde. Entretanto uma esfera de raio maior cresce. Quando o sólido inicial cresce em uma malha numérica relativamente grosseira, a forma do sólido desvia da forma esférica devido perturbações na interface sólido-líquido. Quando a malha é refinada, as perturbações não são detectadas; contudo, quando introduzidas artificialmente as perturbações crescem e distorcem o formato esférico.

Description of manganese nodules and crust collected from the Hakurei Maru Cruise GH76-1, January-March, 1976, in the Central Pacific Basin

This report of the GH76-1 cruise mainly includes the results of the on-board observations in the survey area of the medial-eastern part of Central Pacific Basin (5 degree -10 degree N, 170 degree -175 degree W) and partly of analytical work at the on-shore laboratory. In addition, the results of some on-board optical and geophysical works along the tracks of Japan-Ogasawara-survey area-Hawaii, are described in appendices. The GH76-1 cruise of the R/V Hakurei-maru was carried out from the 10th January to the 9th March, 1976 as the second phase field work of the Geological Survey of Japan five-year research program of study on the manganese nodule deposits of the Central Pacific Basin and also as a part of the National Research Institute for Pollution and Resources research program of technological study on the exploitation of deep sea mineral resources. The GSJ research program (F.Y. 1974-F.Y. 1978) aims at providing basic information on the manganese nodule distribution and their origin on the deep sea floor of the Central Pacific Basin bounded by the Marshall Ridge to the west, the Christmas Ridge to the east, and the Mid-Pacific Mountains to the north. The first phase of investigation was carried out during the GH74-5 cruise in the eastern part of the area (6 degree -10 degree 30'W, 164 degree 30'-171 degree 30'N)(Mizuno and Chujo, eds., 1975), and the present second phase covered an areas of 5 degree square, just west of the GH74-5 area.

Numerical methods for solar tomography in the STEREO era

In this thesis, we propose several advances in the numerical and computational algorithms that are used to determine tomographic estimates of physical parameters in the solar corona. We focus on methods for both global dynamic estimation of the coronal electron density and estimation of local transient phenomena, such as coronal mass ejections, from empirical observations acquired by instruments onboard the STEREO spacecraft. We present a first look at tomographic reconstructions of the solar corona from multiple points-of-view, which motivates the developments in this thesis. In particular, we propose a method for linear equality constrained state estimation that leads toward more physical global dynamic solar tomography estimates. We also present a formulation of the local static estimation problem, i.e., the tomographic estimation of local events and structures like coronal mass ejections, that couples the tomographic imaging problem to a phase field based level set method. This formulation will render feasible the 3D tomography of coronal mass ejections from limited observations. Finally, we develop a scalable algorithm for ray tracing dense meshes, which allows efficient computation of many of the tomographic projection matrices needed for the applications in this thesis.

Spin dynamics of NiCl(2)-4SC(NH(2))(2) in the field-induced ordered phase

NiCl(2)-4SC(NH(2))(2) (known as DTN) is a spin-1 material with a strong single-ion anisotropy that is regarded as a new candidate for Bose-Einstein condensation (BEC) of spin degrees of freedom. We present a systematic study of the low-energy excitation spectrum of DTN in the field-induced magnetically ordered phase by means of high-field electron spin resonance measurements at temperatures down to 0.45 K. We argue that two gapped modes observed in the experiment can be consistently interpreted within a four-sublattice antiferromagnet model with a finite interaction between two tetragonal subsystems and unbroken axial symmetry. The latter is crucial for the interpretation of the field-induced ordering in DTN in terms of BEC.

Sublattice magnetism in sigma-phase Fe(100-x)Vx (x=34.4, 39.9, and 47.9) studied via zero-field (51)V NMR

The successful measurements of a sublattice magnetism with (51)V NMR techniques in the sigma-phase Fe(100-x)V(x) alloys with x=34.4, 39.9, and 47.9 are reported. Vanadium atoms, which were revealed to be present on all five crystallographic sites, are found to be under the action of the hyperfine magnetic fields produced by the neighboring Fe atoms, which allow the observation of (51)V NMR signals. Their nuclear magnetic properties are characteristic of a given site, which strongly depend on the composition. Site A exhibits the strongest magnetism while site D is the weakest. The estimated average magnetic moment per V atom decreases from 0.36 mu(B) for x=34.4 to 0.20 mu(B) for x=47.9. The magnetism revealed at V atoms is linearly correlated with the magnetic moment of Fe atoms, which implies that the former is induced by the latter.

Effect of disorder on quantum phase transitions in anisotropic XY spin chains in a transverse field

We present some exact results for the effect of disorder on the critical properties of an anisotropic XY spin chain in a transverse held. The continuum limit of the corresponding fermion model is taken and in various cases results in a Dirac equation with a random mass. Exact analytic techniques can then be used to evaluate the density of states and the localization length. In the presence of disorder the ferromagnetic-paramagnetic or Ising transition of the model is in the same universality class as the random transverse field Ising model solved by Fisher using a real-space renormalization-group decimation technique (RSRGDT). If there is only randomness in the anisotropy of the magnetic exchange then the anisotropy transition (from a ferromagnet in the x direction to a ferromagnet in the y direction) is also in this universality class. However, if there is randomness in the isotropic part of the exchange or in the transverse held then in a nonzero transverse field the anisotropy transition is destroyed by the disorder. We show that in the Griffiths' phase near the Ising transition that the ground-state energy has an essential singularity. The results obtained for the dynamical critical exponent, typical correlation length, and for the temperature dependence of the specific heat near the Ising transition agree with the results of the RSRODT and numerical work. [S0163-1829(99)07125-8].

High-field (3T) magnetic resonance defecography with functional assessment of the evacuation phase : a pictorial essay : P10

Purpose: Dynamic high-field magnetic resonance (MR) defecography including the evacuation phase is a promising tool for the assessment of functional pelvic disorders, nowadays seen with increasing frequency in elderly women in particular. Learning objectives: 1. To describe the adequate technique of dynamic high-field MRI (3T) in assessing pelvic floor disorders. 2. To provide an overview of the most common pathologies occurring during the evacuation phase, especially in comparison with results of conventional defecography. Methods and materials: After description of the ideal technical parameters of MR defecography performed in supine position after gel rectal filling with a 3 Tesla unit and including the evacuation phase we stress the importance of using a standardized evaluation system for the exact assessment of pelvic floor pathophysiology. Results: The typical pelvic floor disorders occurring before and/or during the evacuation phase, such as sphincter insufficiency, vaginal vault and/or uterine prolapse, cystourethrocele, peritoneo-/ entero-/ sigmoïdocele or rectal prolapse, are demonstrated. The difference between the terms "pelvic floor descent" and "pelvic floor relaxation" are pictorially outlined. MR results are compared with these of conventional defecography. Conclusion: Exact knowledge about the correct technique including the evacuation phase and the use of a standardized evaluation system in assessing pelvic floor disorders by dynamic high-field MRI is mandatory for accurate and reproducible diagnosis.

Advantages of digital holographic microscopy for real-time full field absolute phase imaging

Different interferometric techniques were developed last decade to obtain full field, quantitative, and absolute phase imaging, such as phase-shifting, Fourier phase microscopy, Hilbert phase microscopy or digital holographic microscopy (DHM). Although, these techniques are very similar, DHM combines several advantages. In contrast, to phase shifting, DHM is indeed capable of single-shot hologram recording allowing a real-time absolute phase imaging. On the other hand, unlike to Fourier phase or Hilbert phase microscopy, DHM does not require to record in focus images of the specimen on the digital detector (CCD or CMOS camera), because a numerical focalization adjustment can be performed by a numerical wavefront propagation. Consequently, the depth of view of high NA microscope objectives is numerically extended. For example, two different biological cells, floating at different depths in a liquid, can be focalized numerically from the same digital hologram. Moreover, the numerical propagation associated to digital optics and automatic fitting procedures, permits vibrations insensitive full- field phase imaging and the complete compensation for a priori any image distortion or/and phase aberrations introduced for example by imperfections of holders or perfusion chamber. Examples of real-time full-field phase images of biological cells have been demonstrated. ©2008 COPYRIGHT SPIE

Noise-induced phase separation: Mean-field results

We present a study of a phase-separation process induced by the presence of spatially correlated multiplicative noise. We develop a mean-field approach suitable for conserved-order-parameter systems and use it to obtain the phase diagram of the model. Mean-field results are compared with numerical simulations of the complete model in two dimensions. Additionally, a comparison between the noise-driven dynamics of conserved and nonconserved systems is made at the level of the mean-field approximation.

Field Evaluation of Compaction Monitoring Technology: Phase II Volumes 1 and 2, TR-495, 2006

This report documents an extensive field program carried out to identify the relationships between soil engineering properties, as measured by various in situ devices, and the results of machine compaction monitoring using prototype compaction monitoring technology developed by Caterpillar Inc. Primary research tasks for this study include the following: (1) experimental testing and statistical analyses to evaluate machine power in terms of the engineering properties of the compacted soil (e.g., density, strength, stiffness) and (2) recommendations for using the compaction monitoring technology in practice. The compaction monitoring technology includes sensors that monitor the power consumption used to move the compaction machine, an on-board computer and display screen, and a GPS system to map the spatial location of the machine. In situ soil density, strength, and stiffness data characterized the soil at various stages of compaction. For each test strip or test area, in situ soil properties were compared directly to machine power values to establish statistical relationships. Statistical models were developed to predict soil density, strength, and stiffness from the machine power values. Field data for multiple test strips were evaluated. The R2 correlation coefficient was generally used to assess the quality of the regressions. Strong correlations were observed between averaged machine power and field measurement data. The relationships are based on the compaction model derived from laboratory data. Correlation coefficients (R2) were consistently higher for thicker lifts than for thin lifts, indicating that the depth influencing machine power response exceeds the representative lift thickness encountered under field conditions. Caterpillar Inc. compaction monitoring technology also identified localized areas of an earthwork project with weak or poorly compacted soil. The soil properties at these locations were verified using in situ test devices. This report also documents the steps required to implement the compaction monitoring technology evaluated.

In vivo assessment of myelination by phase imaging at high magnetic field.

The present study evaluated the potential of using the phase of T2* weighted MR images to characterize myelination during brain development and pathology in rodents at 9.4 T. Phase contrast correlated with myelin content assessed by histology and suggests that most contrast between white and cortical gray matter is modulated by myelin. Ex vivo experiments showed that gray-white matter phase contrast remains unchanged after iron extraction. In dysmyelinated shiverer mice, phase imaging correlated strongly with myelin staining, showing reduced contrast between white and gray matter when compared to healthy controls. We conclude that high-resolution phase images, acquired at high field, allow assessment of myelination and dysmyelination.