phase field modelling
play

Phase field modelling Phase field modelling Current challenges and - PowerPoint PPT Presentation

Phase field modelling Phase field modelling Current challenges and opportunities for high performance M P Gururajan computing Preamble Microstructure A few movies P G Tennyson , D Mohan , G Phanikumar , M P Gururajan 1 Phase


  1. Phase field modelling Phase field modelling Current challenges and opportunities for high performance M P Gururajan computing Preamble Microstructure A few movies P G Tennyson † , D Mohan ‡ , G Phanikumar ‡ , M P Gururajan 1 Phase field modelling Department of Metallurgical Engineering and Materials Science, Indian Institute of Technology Bombay, Mumbai 400076 INDIA † : Tata Research Development and Design Center (TRDDC), Pune ‡ : Department of Metallurgical and Materials Engineering, IIT Madras, Chennai December, 2018 1 guru.mp@iitb.ac.in,gururajan.mp@gmail.com

  2. Acknowledgements Phase field modelling M P Gururajan Preamble Microstructure Funding: IRCC, IIT Bombay, DST, Government of India, DRDO A few movies and SASE, Ministry of Defence, Government of India, Phase field modelling DST-DAAD, Tata Steel, GE India Computational resources: Spinode, Dendrite, Nebula / Space-Time, Param Yuva (C-DAC, Pune) Organisers, specifically, Dr. Shenoy, C-DAC, Pune Teachers, collaborators, students

  3. Outline Phase field modelling M P Gururajan Preamble Microstructure A few movies Phase field Microstructure and its evolution modelling Phase field modelling Examples: Six-fold anisotropy on morphology / Solidification Way forward!

  4. Computational Materials Engineering Group Phase field modelling M P Gururajan Preamble Microstructure A few movies Phase field modelling Figure: CMEG: part of materials and process modelling lab

  5. The problem Phase field modelling M P Gururajan Preamble Microstructure A few movies Phase field modelling Figure: New material development cycle: 10 to 20 years. Can we bring it down to less than 5 years? Xiong and Olson, npj Computational Materials, 2016

  6. ICME Phase field modelling M P Gururajan Preamble Microstructure A few movies Phase field modelling Figure: ICME: The minerals, metals and materials society (TMS) study, 2013

  7. Tools and techniques Phase field modelling M P Gururajan Preamble Microstructure A few movies Phase field modelling Figure: Computational materials science: tools and techniques

  8. What is microstructure? Phase field modelling M P Gururajan Preamble Microstructure A few movies Phase field modelling Figure: Microstructure (a Ni-base superalloy). Xu et al, Met. Mat. Trans. A, 1998 Structure, shapes, sizes and distribution of interfaces

  9. Microstructural evolution Phase field modelling M P Gururajan Preamble Microstructure A few movies Phase field modelling Figure: Effect of heat treatment. Xu et al, Met. Mat. Trans. A, 1998

  10. Microstructural evolution Phase field modelling M P Gururajan Preamble Microstructure A few movies Phase field modelling Figure: Dendrites during solidification. David et al, JOM, 2003

  11. Spinodal decomposition Phase field modelling M P Gururajan Preamble Microstructure Homogeneous system with noise A few movies Phase field modelling 0.5 Figure: A homogeneous alloy with a slighlty noisy composition profile

  12. Pure material solidification Phase field modelling M P Gururajan Preamble Undercooled melt Microstructure A few movies Phase field modelling Nucleus Figure: An undercooled melt with insulated sides and nucleus on one of the walls. The interfacial energy is 4-fold anisotropic.

  13. Six-fold dendrites Phase field modelling M P Gururajan Preamble Undercooled melt Microstructure A few movies Phase field modelling Nucleus Figure: An undercooled melt with insulated sides and nucleus at the centre. The interfacial energy is 6-fold anisotropic.

  14. Spinodal decomposition Phase field modelling M P Gururajan Preamble Microstructure A few movies Phase field modelling Figure: Regions rich in A (B) become richer in A (B) with time. Microstructures at times 0, 100 and 1000 units.

  15. Issue 1 Phase field modelling M P Gururajan The phase field method, like many other modeling ap- Preamble proaches, is practically limited by the computational expense Microstructure A few movies entailed in running large simulations. The challenge stems Phase field from the need to resolve a diffuse interface that has a dif- modelling fuseness that is on a much smaller length scale than a typical microstructural evolution length scale. –Modeling Across Scales: A Roadmapping Study for Con- necting Materials Models and Simulations Across Length and Time Scales, TMS study report, 2015

  16. Issue 2 Phase field modelling M P Gururajan Preamble Microstructure A few movies Common software: Micress TM , FiPy TM , OpenPhase TM , and MOOSE Phase field modelling (Marmot) TM Compare with VASP, LAMMPS, ParaDIS, ... Phase field: an approach and not a set methodology (like FEM) pfHUB: maintained by NIST

  17. Phase field models Phase field modelling M P Gururajan Preamble Microstructure A few movies ∂ c Phase field ∂ t = ∇ M ∇ µ = ∇ M ∇ [ g ( c ) − κ ∇ 2 c ] (1) modelling ∂φ ∂ t = − L µ = L [ κ ∇ 2 φ − g ( φ )] (2) Ginzburg-Landau, Alan Turing (Chemical morphogenesis), ...

  18. Characteristics of phase field models Phase field modelling M P Gururajan Preamble Microstructure Interfaces are not sharp; diffuse interface model A few movies No tracking of interface: numerical solutions are easier Phase field modelling Gradient energy coefficient: interfacial energy contributions (Gibbs-Thomson, for example) are automatically accounted for Topological singularities (splitting or disappearance of interfaces): naturally taken care of Elastic stress, magnetic and electric field: can be coupled by adding the relevant free energy term!

  19. What is phase field modelling? Phase field modelling M P Gururajan Some representative viewpoints: Preamble Microstructure An approach to obtain solutions of PDEs that are hard to solve A few movies – by introducing artificial regions of continuity where there are Phase field modelling discontinuities (Mathematical) Non-linear partial differential equations that lead to solutions which are interesting patterns (Biology) Continuum equations derived from statistical mechanics that lead (as solution) to interesting patterns (Physics) Partial differential equations that describe diffusion (of atoms and heat) as well as phase transformations (Materials science)

  20. Spectral technique Phase field modelling M P Gururajan ∂ c ∂ t = D ∇ 2 c (3) Preamble Microstructure Spatial Fourier transform of c: ˜ c = � c ( x ) exp [ − i k · r ] dV A few movies Turns the PDE into ODE: Phase field modelling d ˜ c dt = − Dk 2 ˜ c (4) Semi-implicit Fourier spectral technique ∂ c ∂ t = ∇ M ∇ µ = ∇ M ∇ [ g ( c ) − κ ∇ 2 c ] (5) ∂φ ∂ t = − L µ = L [ κ ∇ 2 φ − g ( φ )] (6)

  21. Advantages of FFT Phase field modelling M P Gururajan Preamble Microstructure A few movies Phase field modelling Periodic boundary conditions: representative volume elements Semi-implicit Fourier spectral technique Good, fast, open source FFT codes: FFTW

  22. Extended Cahn-Hilliard free energy: anisotropic interfacial energy Phase field modelling M P Gururajan Preamble Microstructure µ = g ( c ) A few movies 2 κ I Phase field − ij c ij modelling 12 β I ijkl c ij c k c l + 2 β III − ijkl c ijkl 30 α I ijklmn c ij c k c l c m c n − 2 α VII (7) − ijklmn c ijklmn For details: E S Nani and M P Gururajan, Philosophical Magazine Letters (2014)

  23. Six fold anisotropy Phase field modelling M P Gururajan Preamble Microstructure A few movies Phase field modelling Interfacial energy anisotropy / Point effect of diffusion / FG to CG Att. kinetics anisotropy / SG to CG / Noise and Point Effect of Diff From unpublished M Tech thesis of Mr. Abhinav Soni

  24. Profiling on NIVIDIA � - K40C GPUs (Ternary alloy code) Phase field Strong scaling modelling M P Gururajan Preamble Microstructure A few movies Phase field modelling Weak scaling

  25. Profiling Phase field modelling M P Gururajan Preamble Microstructure A few movies Phase field modelling Figure: Ni - 19 Cr - 5 Nb (wt.%) alloy - 3D isothermal dendrite at ∆ T = 8 . 0 K, ∆ t = 58 . 0 ns for ∆ x = 50 . 0 nm. Figure: Profiling of 3D phase-field code (384 × 384 × 1024) Mohan and Phanikumar, Unpublished

  26. Performance Phase field modelling 5000 M P Gururajan ATI HD5870 MPI 16 Preamble MPI 32 4000 Microstructure MPI 64 A few movies Time in seconds Phase field 3000 Time taken for 10000 timesteps modelling 2000 1000 0 1024 1536 2048 2560 3672 4096 Domain length P G Tennyson, G M Karthik, and G Phanikumar, Computer Physics Communications, 2015.

  27. Data Visualization Phase field modelling M P Gururajan 4000 Preamble Data files from each processor 0.16 3500 0.14 Microstructure written at specific time intervals as 3000 0.12 2500 A few movies 0.1 unformatted .bin files 2000 0.08 Phase field 1500 modelling Data files collated and converted to 0.06 1000 0.04 500 .mat files by Mat I/O library by 0.02 500 1000 1500 2000 2500 3000 3500 4000 Christopher Hulbert Visualization of data output was done in Matlab R � MayaVi, created by Prabhu Ramachandran, was used for 3- D data visualization P G Tennyson, G M Karthik, and G Phanikumar, Computer Physics Communications, 2015.

Recommend


More recommend