<p>Continuum atomistic and mesoscopic models are developed and utilized in the context of studying microstructural evolution and phase selection in materials systems. Numerous phenomena are examined, ranging from defect-solute interaction in solid state systems to microstructural evolution under external driving conditions. Emphasis is placed on the derivation and development of models capable of self consistently describing the intricate mechanisms at work in the systems undergoing these phenomena.</p> <p>Namely, grain growth dynamics are studied in nanocrystalline systems under external driving conditions using a newly developed phase-field-crystal model, which couples an additional free energy source term to the standard phase-field-crystal model. Such external driving can be attributed to incident energetic particles. The nanocrystalline system is found to be susceptible to enhanced grain growth as a function of the intensity/flux associated with the external driving and the energy of driving. Static kinetic phase diagram calculations also seem to confirm that systems under external driving conditions can be forced into long metastable states.</p> <p>Early stage solute clustering and precipitation in Al alloys is also examined with a variant of the phase-field-crystal method, so-called structural phase-field-crystal models for multi-component alloys developed as part of this thesis. We find that clustering is aided by quenched-in defects (dislocations), whereby the nucleation barrier is reduced and at times eliminated, a mechanism proposed by Cahn for a single dislocation for spinodal systems. In a three-component system, we predict a multi-step mechanism for clustering, where the nature and amount of the third species plays an important role in relieving stresses caused by the quenched-in dislocations before clustering, i.e., segregation as predicted by the equilibrium phase diagram, can occur.</p> <p>Finally, we present a new coarse-graining procedure for generating complex amplitude models, i.e., complex order-parameter phase-field models, derived from phase-field-crystal models. They retain many salient atomistic features and behaviours of the original phase-field-crystal model, however is now capable of describing mesoscopic length scales like the phase-field model. We demonstrate the scheme by generating an amplitude model of the two-dimensional structural phase-fieldcrystal model, which allows multiple crystal structures to be stable in equilibrium, a crucial aspect of proper multi-scale modelling of materials systems. The dynamics are demonstrated by examining solidification and coarsening, peritectic growth, along with grain growth and the emergence of secondary phases.</p> / Doctor of Science (PhD)

Martensitic phase transformation in NiTi shape memory alloys (SMA) occurs over a hierarchy of spatial scales, as evidenced from observed multiscale patterns of the martensitic phase fraction, which depend on the material microstructure and on the size of the SMA specimen. This paper presents a methodology for the multiscale tracking of the thermally induced martensitic phase transformation process in NiTi SMA. Fine scale stochastic phase field simulations are coupled to macroscale experimental measurements through the compound wavelet matrix method (CWM). A novel process for obtaining CWM fine scale wavelet coefficients is used that enhances the effectiveness of the method in transferring uncertainties from fine to coarse scales, and also ensures the preservation of spatial correlations in the phase fraction pattern. Size effects, well-documented in the literature, play an important role in designing the multiscale tracking methodology. Molecular dynamics (MD) simulations are employed to verify the phase field simulations in terms of different statistical measures and to demonstrate size effects at the nanometer scale. The effects of thermally induced martensite phase fraction uncertainties on the constitutive response of NiTi SMA is demonstrated.

A sinterização do cobre ocorre através do transporte difusivo de matéria de regiões de alto para as de baixo potencial químico. Esse processo tem como força motriz a minimização da energia associada às interfaces. Na tentativa de quantificar o processo de sinterizaçao, diversos modelos analíticos foram desenvolvidos desde 1945. O presente trabalho teve como intuito implementar e validar um modelo matemático baseado no modelo de campo de fases (\"phase field method\") para simular o primeiro estágio do processo de sinterização. Para isso, um estudo termodinâmico detalhado foi realizado de modo a definir as equações a serem empregadas no modelo. Foi feita uma análise quantitativa (análise estatística) e qualitativamente (análise gráfica e pelo expoente do tempo) dos modelos analíticos teóricos comparando-os com os valores experimentais publicados em cinco artigos de grande relevância na área. A partir dos resultados estatísticos observou-se que o melhor mecanismo para descrever o processo de sinterização do cobre é o modelo combinado entre quatro dos seis principais modelos individuais. Os mecanismos de transporte de fase fictícia via evaporação-condensação e via difusão gasosa contribuem de maneira irrisória na sinterização do cobre, sendo por muitos autores desconsiderados. Foi verificado que a configuração inicial do metal, se na forma de esferas ou cilindros, modifica o processo de difusão dominante. Sendo que o efeito da difusão superficial é mais predominante nas esferas do que nos cilindros, consequentemente, o mecanismo combinado para a esferas inclui a difusão superficial, enquanto que o dos cilindros não. Na simulação em condições unidimensionais, o modelo foi capaz de impor as condições de equilíbrio termodinâmico local e de movimentar a interface no sentido contrário ao fluxo de lacunas. Sob condições bidimensionais, o mesmo impôs automaticamente a fração de lacunas de equilíbrio sob o efeito do raio de curvatura, responsável pela expansão ou retração do sólido/poro cilíndrico. Na simulação da formação do pescoço entre dois cilindros de cobre puro, observou-se um comportamento qualitativo consistente com o comportamento físico. A principal dificuldade encontrada na modelagem foi o tamanho da malha e o tempo de processamento computacional necessário. Para resolução destes dois aspectos, usou-se uma malha adaptativa e foi feita a paralelização em placa de vídeo do código computacional. / Copper sintering occurs through a diffusive transport of matter from regions with high chemical potential to regions of low chemical potential. The driving force of this process is the minimization of the energy associated with the interfaces of the system. In an attempt to quantify the sintering process, several analytical models have been developed since 1945. The aim of the present work was to implement and validate a mathematical model based on the phase field model to simulate the first stage of sintering. A very detailed thermodynamic study was done in order to define which equations should me used in the computational model. As well as, the use of a quantitative (statistical analysis) and a qualitative analysis (graphical analysis and by the exponent of time) to compare the theoretical models with the experimental values published in five articles of great relevance in the area. From the statistical results it was observed that the best mechanism to describe the copper sintering is the combined model between the main individual models (lattice diffusion from surface, lattice diffusion from grain boundary, surface diffusion and grain boundary diffusion). The mechanisms of gas-phase transport via evaporation-condensation and gas diffusion contribute in a negligible way in copper sintering, considered irrelevant by many authors . It has been found that the initial configuration of the metal, whether in the form of spheres or cylinders, modifies the dominant diffusion mechanism. Since the effect of surface diffusion is more predominant in the spheres than in the cylinders, therefore the combined mechanism for the spheres includes surface diffusion, while the cylinders do not. With the computational modeling, some important mechanisms that occur during stage I of sintering were simulated. Simulations performed under unidimensional conditions indicated that the model is able to impose local thermodynamic equilibrium conditions and to move the interface in the opposite direction of the vacancies flow. When used to simulate the transport of vacancies under two-dimensional conditions, the model automatically imposed the fraction of equilibrium vacancies under the effect of the radius of curvature. This fraction results in a flow that causes the expansion or retraction of the solid / cylindrical pore, which was reproduced by the implemented model. The main difficulties found in the computational modeling were the size of the mesh and the computational processing time required. To solve these two aspects, an adaptive mesh was used and the parallelization of the computational code was done, which resulted in a significant reduction in the simulation time.

La ferrite allotriomorphe est une des morphologies de la ferrite dont la répartition spatiale influe fortement sur les propriétés mécaniques dans les aciers dual-phase. En fonction des traitements qu'ils subissent, la ferrite peut s'y répartir suivant les bandes de ségrégation en manganèse, issues de l'étape de solidification. Pour établir le rôle que joue le processus de croissance de la ferrite allotriomorphe sur la mise en place de la structure en bandes, nous avons développé un modèle de champ de phases possédant deux spécificités originales, imposées par le problème. D'une part, ce modèle est capable de reproduire les différents régimes cinétiques observés dans les alliages ternaires Fe-C-X, pilotés par la présence concomittante du carbone diffusant rapidement,et d'un élément substitutionnel X diffusant lentement. Nous avons ainsi mis en évidence la transition d'un régime initial rapide de paraéquilibre vers une croissance lente en orthoéquilibre, en bon accord avec des résultats expérimentaux de la littérature. D'autre part, notre modèle incorpore de manière économe la présence des joints de grains austénitiques, dont le rôle dans l'élimination des structures en bande est souligné par nos calculs. Nous observons ainsi qu'il existe un seuil d'intensité deségrégation en manganèse en dessous duquel le mouillage de la ferrite le long des joints de grain de plus grande énergie peut contrecarrer la croissance dans les bandes ségrégées négativement / The growth of allotriomorphic ferrite plays a major role in the formation of martensite bands in Dual-Phase steels. We have thus developed a phase field model to study the ferritic growth in different ternary Fe-C-X alloys, incorporating two necessary features. First, we have paid a particular attention to recover the different growth regimes due to the huge difference between the diffusion rates of Cand X substitutional species. Our calculations have exhibited a transition from fast paraequilibrium to slow orthoequilibrium in good agreement with experimental measurements in the literature. Second, austenite grain boundaries have been included in the model because they can conterbalance the manganese segregation bands, as shown in our calculations. Indeed, our results show that the bands can be broken bythe wetting of ferrite along the austenite grain boundaries, provided that the segregation is below a threshold value, and provided that the grain boundary energies are sufficiently high

The primary objective of this work is to propose a state-of-the-art physics based multiscale modeling framework for simulating material phase change problems. Both ice melting and copper crystallization problems are selected to demonstrate this multiscale modeling and simulation. The computational methods employed in this thesis include: classical molecular dynamics, finite element method, phase-field method, and multiscale (nano/micro coupling) methods. Classical molecular dynamics (MD) is a well-known method to study material behaviors at atomic level. Due to the limit of MD, it is not realistic to provide a complete molecular model for simulations at large length and time scales. Continuum methods, including finite element methods, should be employed in this case. In this thesis, MD is employed to study phase change problems at the nanoscale. In order to study material phase change problems at the microscale, a thermal wave method one-way coupling with the MD and a phase-field method one-way coupling with MD are proposed. The thermal wave method is more accurate than classical thermal diffusion for the study of heat transfer problems especially in crystal based structures. The second model is based on the well-known phase-field method. It is modified to respond to the thermal propagation in the crystal matrix by the thermal wave method, as well as modified to respond to temperature gradients and heat fluxes by employing the Dual-Phase-Lag method. Both methods are coupled with MD to obtain realistic results. It should be noted that MD simulations can be conducted to obtain material/thermal properties for microscopic and/or macroscopic simulations for the purpose of hierarchical/sequential multiscale modeling. These material parameters include thermal conductivity, specific heat, latent heat, and relaxation time. Other type of interfacial parameters that occur during the phase change process, such as nucleus shape, interfacial energy, interfacial thickness, etc., are also obtained by MD simulation since these have so far been too difficult to measure experimentally. I consider two common phase change phenomena, ice melting and copper crystallization, in this thesis. For the case of ice melting, MD is first employed to study its phase change process and obtain thermal properties of ice and water. Several potential models are used. I conduct simulations of both bulk ice and ice/water contacting cases. It is found that various potential models result in similar melting phenomena, especially melting speed. Size effects are also studied and it is found that the melting time is longer for larger bulk ice segments but that the average melting speed is size dependent. There is no size effect for the melting speed at ice/water interface at the nanoscale if the same temperature gradient is applied. The melting speed of ice should depend on the temperature gradient. To study ice melting at the microscale, the thermal wave model is employed with parameters obtained from MD simulations. It is found that ice melting speed is scale, for both length scale and time scale, dependent. For the case of copper crystallization, an EAM potential is first employed to conduct MD simulations for studying the copper crystallization process at the nanoscale. I obtain thermal properties and interfacial parameters, including thermal diffusion coefficient, latent heat, relaxation time, interfacial thickness, interfacial energy and the anisotropy coefficients, and nucleus shape etc. A central symmetry parameter is used to identify an atom in solid state or liquid state. And then an initial nucleus shape is obtained and used as the input for microscale simulation, in which the phase-field method is used to study copper crystallization at the microscale.

Martensite is considered to be the backbone of the high strength of many commercial steels. Martensite is formed by a rapid diffusionless phase transformation, which has been the subject of extensive research studies for more than a century. Despite such extensive studies, martensitic transformation is still considered to be intriguing due to its complex nature. Phase-field method, a computational technique used to simulate phase transformations, could be an aid in understanding the transformation. Moreover, due to the growing interest in the field of “Integrated computational materials engineering (ICME)”, the possibilities to couple the phase-field method with other computational techniques need to be explored. In the present work a three dimensional elastoplastic phase-field model, based on the works of Khachaturyan et al. and Yamanaka et al., is developed to study the athermal and the stress-assisted martensitic transformations occurring in single crystal and polycrystalline steels. The material parameters corresponding to the carbon steels and stainless steels are considered as input data for the simulations. The input data for the simulations is acquired from computational as well as from experimental works. Thus an attempt is made to create a multi-length scale model by coupling the ab-initio method, phase-field method, CALPHAD method, as well as experimental works. The model is used to simulate the microstructure evolution as well as to study various physical concepts associated with the martensitic transformation. The simulation results depict several experimentally observed aspects associated with the martensitic transformation, such as twinned microstructure and autocatalysis. The results indicate that plastic deformation and autocatalysis play a significant role in the martensitic microstructure evolution. The results indicate that the phase-field simulations can be used as tools to study some of the physical concepts associated with martensitic transformation, e.g. embryo potency, driving forces, plastic deformation as well as some aspects of crystallography. The results obtained are in agreement with the experimental results. The effect of stress-states on the stress-assisted martensitic microstructure evolution is studied by performing different simulations under different loading conditions. The results indicate that the microstructure is significantly affected by the loading conditions. The simulations are also used to study several important aspects, such as TRIP effect and Magee effect. The model is also used to predict some of the practically important parameters such as Ms temperature as well as the volume fraction of martensite formed. The results also indicate that it is feasible to build physically based multi-length scale model to study the martensitic transformation. Finally, it is concluded that the phase-field method can be used as a qualitative aid in understanding the complex, yet intriguing, martensitic transformations. / QC 20120525 / Hero-m

In this thesis, we have used a phase field model for studying the effect of atomic mobility inside the precipitate phase on coarsening behaviour in two dimensional (2D) systems. In all the available coarsening theories, the diffusivity inside the precipitate phase is not explicitly taken into account; this would imply that there is no chemical potential gradient inside the precipitate. This assumption is valid if (a) the atomic mobility inside the precipitate is much higher than that in the matrix, or (b) the precipitate volume fraction is small (i.e. the interparticle spacing is far higher than the average particle size). We undertook this study to evaluate the potential effect of diffusivity in the precipitate on coarsening in situations where conditions (a) and (b), above, do not hold, by studying systems with moderate volume fractions (20% and 30%) and with low atomic mobilities in the precipitate. In our study, we have fixed the atomic mobility in the matrix at a constant value. We have used the well known Cahn-Hilliard model in which the microstructure is described in terms of a composition field variable. The evolution of microstructure is studied by numerically solving a non-classical diffusion equation known as the Cahn-Hilliard equation. We have used a semi-implicit Fourier spectral technique for solving the CH equation using periodic boundary conditions. The coarsening behaviour is tracked and analyzed using number density of particles, their average size and their size distribution. The main conclusion from this study is that, contrary to expectations, the atomic mobility in the precipitate phase has only a small effect on coarsening behavior. Specifically, with decreasing atomic mobility in the precipitate phase, we report a small increase in the number density, a slightly wider size distribution and a slightly smaller coarsening rate. We also add that these effects are too small to allow experimental verification. These results indicate that the need for chemical potential equilibration within each precipitate is not an important factor during coarsening.

In this thesis two kinds of phase change i.e., solid state phase transformation in steels and solid-to-liquid phase transformation in paraffin, have been modeled and numerically simulated. The solid state phase transformation is modeled using the phase field theory while the solid-to-liquid phase transformation is modeled using the Stokes equation and exploiting the viscous nature of the paraffin, by treating it as a liquid in both states.The theoretical base of the solid state, diffusionless phase transformation or the martensitic transformation comes from the Khachaturyan's phase field microelasticity theory. The time evolution of the variable describing the phase transformation is computed using the time dependent Ginzburg-Landau equation. Plasticity is also incorporated into the model by solving another time dependent equation. Simulations are performed both in 2D and 3D, for a single crystal and a polycrystal. Although the model is valid for most iron-carbon alloys, in this research an Fe-0.3\%C alloy is chosen.In order to simulate martensitic transformation in a polycrystal, it is necessary to include the effect of the grain boundary to correctly capture the morphology of the microstructure. One of the important achievements of this research is the incorporation of the grain boundary effect in the Khachaturyan's phase field model. The developed model is also employed to analyze the effect of external stresses on the martensitic transformation, both in 2D and 3D. Results obtained from the numerical simulations show good qualitative agreement with the empirical observations found in the literature.The microactuators are generally used as a micropump or microvalve in various miniaturized industrial and engineering applications. The phase transformation in a paraffin based thermohydraulic membrane microactuator is modeled by treating paraffin as a highly viscous liquid, instead of a solid, below its melting point.  The fluid-solid interaction between paraffin and the enclosing membrane is governed by the ALE technique. The thing which sets apart the presented model from the previous models, is the use of geometry independent and realistic thermal and mechanical properties. Numerical results obtained by treating paraffin as a liquid in both states show better conformity with the experiments, performed on a similar microactuator. The developed model is further employed to analyze the time response of the system, for different input powers and geometries of the microactuator. / <p>QC 20130219</p>

No description available.

No description available.