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.

Homography-based ground plane detection using a single on-board camera

This study presents a robust method for ground plane detection in vision-based systems with a non-stationary camera. The proposed method is based on the reliable estimation of the homography between ground planes in successive images. This homography is computed using a feature matching approach, which in contrast to classical approaches to on-board motion estimation does not require explicit ego-motion calculation. As opposed to it, a novel homography calculation method based on a linear estimation framework is presented. This framework provides predictions of the ground plane transformation matrix that are dynamically updated with new measurements. The method is specially suited for challenging environments, in particular traffic scenarios, in which the information is scarce and the homography computed from the images is usually inaccurate or erroneous. The proposed estimation framework is able to remove erroneous measurements and to correct those that are inaccurate, hence producing a reliable homography estimate at each instant. It is based on the evaluation of the difference between the predicted and the observed transformations, measured according to the spectral norm of the associated matrix of differences. Moreover, an example is provided on how to use the information extracted from ground plane estimation to achieve object detection and tracking. The method has been successfully demonstrated for the detection of moving vehicles in traffic environments.

Selective area growth of a- and c-plane GaN nanocolumns by molecular beam epitaxy using colloidal nanolithography

Selective area growth of a-plane GaN nanocolumns by molecular beam epitaxy was performed for the first time on a-plane GaN templates. Ti masks with 150 nm diameter nanoholes were fabricated by colloidal lithography, an easy, fast and cheap process capable to handle large areas. Even though colloidal lithography does not provide a perfect geometrical arrangement like e-beam lithography, it produces a very homogeneous mask in terms of nanohole diameter and density, and is used here for the first time for the selective area growth of GaN. Selective area growth of a-plane GaN nanocolumns is compared, in terms of anisotropic lateral and vertical growth rates, with GaN nanocolumns grown selectively on the c-plane

Single phase a-plane MgZnO epilayers for UV optoelectronics: substitutional behaviour of Mg at large contents

High quality 1 μm thick a-plane MgxZn1−xO layers were produced by molecular beam epitaxy with Mg contents higher than 50%. Resonant Rutherford backscattering spectrometry combined with ion channeling revealed a uniform growth in both composition and atomic order. The lattice-site location of Mg, Zn and O elements was determined independently, proving the substitutional behaviour of Mg in Zn-sites of the wurtzite lattice. X-Ray diffraction pole figure analysis also confirms the absence of phase separation. Optical properties at such high Mg contents were studied in Schottky photodiodes.

Failure locus of polypropylene nonwoven fabrics under in-plane biaxial deformation

The failure locus, the characteristics of the stress–strain curve and the damage localization patterns were analyzed in a polypropylene nonwoven fabric under in-plane biaxial deformation. The analysis was carried out by means of a homogenization model developed within the context of the finite element method. It provides the constitutive response for a mesodomain of the fabric corresponding to the area associated to a finite element and takes into account the main deformation and damage mechanisms experimentally observed. It was found that the failure locus in the stress space was accurately predicted by the Von Mises criterion and failure took place by the localization of damage into a crack perpendicular to the main loading axis.

A constitutive model for the in-plane mechanical behavior of nonwoven fabrics

A constitutive model is presented for the in-plane mechanical behavior of nonwoven fabrics. The model is developed within the context of the finite element method and provides the constitutive response for a mesodomain of the fabric corresponding to the area associated to a finite element. The model is built upon the ensemble of three blocks, namely fabric, fibers and damage. The continuum tensorial formulation of the fabric response rigorously takes into account the effect of fiber rotation for large strains and includes the nonlinear fiber behavior. In addition, the various damage mechanisms experimentally observed (bond and fiber fracture, interfiber friction and fiber pull-out) are included in a phenomenological way and the random nature of these materials is also taken into account by means of a Monte Carlo lottery to determine the damage thresholds. The model results are validated with recent experimental results on the tensile response of smooth and notched specimens of a polypropylene nonwoven fabric.

Deep levels in a-plane, high Mg-content MgxZn1-xO epitaxial layers grown by molecular beam epitaxy

Deep level defects in n-type unintentionally doped a-plane MgxZn1−xO, grown by molecular beam epitaxy on r-plane sapphire were fully characterized using deep level optical spectroscopy (DLOS) and related methods. Four compositions of MgxZn1−xO were examined with x = 0.31, 0.44, 0.52, and 0.56 together with a control ZnO sample. DLOS measurements revealed the presence of five deep levels in each Mg-containing sample, having energy levels of Ec − 1.4 eV, 2.1 eV, 2.6 V, and Ev + 0.3 eV and 0.6 eV. For all Mg compositions, the activation energies of the first three states were constant with respect to the conduction band edge, whereas the latter two revealed constant activation energies with respect to the valence band edge. In contrast to the ternary materials, only three levels, at Ec − 2.1 eV, Ev + 0.3 eV, and 0.6 eV, were observed for the ZnO control sample in this systematically grown series of samples. Substantially higher concentrations of the deep levels at Ev + 0.3 eV and Ec − 2.1 eV were observed in ZnO compared to the Mg alloyed samples. Moreover, there is a general invariance of trap concentration of the Ev + 0.3 eV and 0.6 eV levels on Mg content, while at least and order of magnitude dependency of the Ec − 1.4 eV and Ec − 2.6 eV levels in Mg alloyed samples.

A canonical UTD solution for electromagnetic scattering by an electrically large impedance circular cylinder illuminated by an obliquely incident plane wave

A uniform geometrical theory of diffraction (UTD) solution is developed for the canonical problem of the electromagnetic (EM) scattering by an electrically large circular cylinder with a uniform impedance boundary condition (IBC), when it is illuminated by an obliquely incident high frequency plane wave. A solution to this canonical problem is first constructed in terms of an exact formulation involving a radially propagating eigenfunction expansion. The latter is converted into a circumferentially propagating eigenfunction expansion suited for large cylinders, via the Watson transform, which is expressed as an integral that is subsequently evaluated asymptotically, for high frequencies, in a uniform manner. The resulting solution is then expressed in the desired UTD ray form. This solution is uniform in the sense that it has the important property that it remains continuous across the transition region on either side of the surface shadow boundary. Outside the shadow boundary transition region it recovers the purely ray optical incident and reflected ray fields on the deep lit side of the shadow boundary and to the modal surface diffracted ray fields on the deep shadow side. The scattered field is seen to have a cross-polarized component due to the coupling between the TEz and TMz waves (where z is the cylinder axis) resulting from the IBC. Such cross-polarization vanishes for normal incidence on the cylinder, and also in the deep lit region for oblique incidence where it properly reduces to the geometrical optics (GO) or ray optical solution. This UTD solution is shown to be very accurate by a numerical comparison with an exact reference solution.

HABIT o cómo medir la subjetividad en la movilidad en transporte urbano

Tradicionalmente en los modelos de transporte, las variables explicativas de la elección modal del usuario se reducen a las variables objetivas (coste del transporte, tiempo de viaje, etc.). Sin embargo, la literatura más reciente reconoce que el perfil del usuario así como sus percepciones, influyen también en la elección modal. Por ello, antes de formular un modelo de demanda de transporte, es conveniente integrar factores más subjetivos (percepciones y actitudes) de los usuarios del sistema de transporte. Para evaluar dicha influencia se ha optado por realizar una encuesta en las nuevas estaciones de la línea 2 de Metro, al Sureste de Madrid. Será una manera de medir la consistencia de los hábitos cuando una nueva infraestructura de transporte como el Metro se pone en marcha. La toma de datos se realiza a través de encuestas de preferencias reveladas a los usuarios del nuevo Metro. Gracias a dicha encuesta se puede conocer el anterior modo de viaje utilizado por los usuarios para realizar el mismo trayecto, por lo que se dispone de datos para evaluar el cambio de modo. En concreto, el encuestado responde a preguntas sobre el motivo del viaje, la disponibilidad de vehículo propio y el motivo principal para cambiar de modo. La encuesta de movilidad realizada en la nueva infraestructura permite recoger para desarrollar una metodología capaz de estimar los efectos de shock e inercia del cambio modal. Este artículo utiliza modelos de elección discreta integrando variables cuantitativas (coste del viaje, por ejemplo) y cualitativas (motivo del cambio modal) para explicar la elección modal. Los primeros resultados muestran que los viajes obligados (trabajo y estudios) que antes del Metro se realizaban en coche, presentan poca resistencia al cambio.

A new UTD for the electromagnetic plane wave scattering by a canonical circular cylinder with an impedance boundary condition

A UTD solution is developed for describing the scattering by a circular cylinder with an impedance boundary condition (IBC), when it is illuminated by an obliquely incident electromagnetic (EM) plane wave. The solution to this canonical problem will be crucial for the construction of a more general UTD solution valid for an arbitrary smooth convex surface with an IBC, when it is illuminated by an arbitrary EM ray optical field. The canonical solution is uniformly valid across the surface shadow boundary that is tangent to the surface at grazing incidence. This canonical solution contains cross polarized terms in the scattered fields, which arise from a coupling of the TEz and TMz waves at the impedance boundary on the cylinder. Here, z is the cylinder axis. Numerical results show very good accuracy for the simpler and efficient UTD solution, when compared to exact but very slowly convergent eigenfunction solution.

Plano Horizontal Plano. Del plano horizontal como límite entre lo estereotómico y lo tectónico = Flat Horizontal Plane. On Horizon

Plano Horizontal Plano. Del plano horizontal como límite entre lo estereotómico y lo tectónico

Flat Horizontal Plane - On Horizon

The Flat Horizontal Plane, the platform, is more than just one of the most basic mechanisms of Architecture. In this essay, I would like to move towards understanding this Flat Horizontal Plane not only as the primary mechanism of Architecture, but also, when it is erected, as the spatial limit between the stereotomic and the tectonic.

Light polarization sensitive photodetectors with m- and r-plane homoepitaxial ZnO/ZnMgO quantum wells

Homoepitaxial ZnO/(Zn,Mg)O multiple quantum wells (MQWs) grown with m- and r-plane orientations are used to demonstrate Schottky photodiodes sensitive to the polarization state of light. In both orientations, the spectral photoresponse of the MQW photodiodes shows a sharp excitonic absorption edge at 3.48 eV with a very low Urbach tail, allowing the observation of the absorption from the A, B and C excitonic transitions. The absorption edge energy is shifted by ∼30 and ∼15 meV for the m- and r-plane MQW photodiodes, respectively, in full agreement with the calculated polarization of the A, B, and C excitonic transitions. The best figures of merit are obtained for the m-plane photodiodes, which present a quantum efficiency of ∼11%, and a specific detectivity D* of ∼6.4 × 1010 cm Hz1/2/W. In these photodiodes, the absorption polarization sensitivity contrast between the two orthogonal in-plane axes yields a maximum value of (R⊥/R||)max ∼ 9.9 with a narrow bandwidth of ∼33 meV.