ICME 2017 HW2

From EVOCD
(Difference between revisions)
Jump to: navigation, search
(Step 4: no Ensight)
(Environment Setup)
 
(16 intermediate revisions by 3 users not shown)
Line 8: Line 8:
 
* Dislocation Dynamics (DD) using Multiscale Dislocation Dynamics Plasticity (MDDP)
 
* Dislocation Dynamics (DD) using Multiscale Dislocation Dynamics Plasticity (MDDP)
  
 
+
All the necessary files can be found on the ICME website. Copy them to a folder on your computer.
All necessary input files and scripts are provided here or in /scratch/ICME_2017/Homework2/ . Move any of these files to your own directory (and make a backup copy) before trying to perform any simulations.
+
 
+
 
Use /scratch/"Your Directory" for best results, especially if your job reads/writes much data.
 
Use /scratch/"Your Directory" for best results, especially if your job reads/writes much data.
  
 
Write a full report that follows a journal article manuscript format (include figures and tables in the text).
 
Write a full report that follows a journal article manuscript format (include figures and tables in the text).
  
Upon completion, submit a .pdf and .doc(x) file of your report. Be sure to also include the requested files and plots from each section of the homework.  
+
Upon completion, submit via email a .pdf and .doc(x) file of your report. Be sure to also include the requested files and plots from each section of the homework.
  
 
=Part 1 - Single Dislocation Mobility Calculations using MEAM=
 
=Part 1 - Single Dislocation Mobility Calculations using MEAM=
Line 24: Line 22:
 
You can obtain LAMMPS and find the user manual [http://lammps.sandia.gov/doc/Manual.html here].
 
You can obtain LAMMPS and find the user manual [http://lammps.sandia.gov/doc/Manual.html here].
 
Visit the [[LAMMPS tutorials]] page for a wide range of examples.
 
Visit the [[LAMMPS tutorials]] page for a wide range of examples.
 +
LAMMPS video tutorials - [https://youtu.be/GXA2PyqKYdY Intro] [https://www.youtube.com/watch?v=7Ila18g8zSY Input file Part 1] [https://www.youtube.com/watch?v=BOJPl9A7-K8 Input file Part 2]
 +
[https://www.youtube.com/watch?v=TTDXJXJJi18/ Molecular Dynamics Theory and Examples]
  
 
==Environment Setup==
 
==Environment Setup==
Line 39: Line 39:
  
 
===Walkthrough===
 
===Walkthrough===
 +
 +
Video tutorials: [https://www.youtube.com/watch?v=S7w2wFHzX7o Part 1] [https://youtu.be/oZwmPV5mK0s Part 2]
  
 
====Step 1====
 
====Step 1====
Generate the atomic structure using [[Atomistic Dislocation Generation]]. Once you have downloaded and compiled the fortran routine, simply run it at the command line as explained on the script page.
+
Generate the atomic structure using [[Atomistic Dislocation Generation]]. Once you have downloaded and compiled the fortran routine, as directed, simply run it at the command line as explained on the script page.
For this calculation, you want make an edge dislocation in a PAD geometry. It will create the required data file as atoms.*.edge.pad, where * is either fcc or bcc, depending on your material. You will need to copy the data file to the directory where you will run the MD simulation.
+
For this calculation, you want to make an edge dislocation in a PAD geometry. It will create the required data file as atoms.*.edge.pad, where * is either fcc or bcc, depending on your material. You will need to copy the data file to the directory where you will run the MD simulation.
  
 
====Step 2====
 
====Step 2====
Line 63: Line 65:
  
 
====Step 3====
 
====Step 3====
 +
Run the simulation:
 +
<pre>
 +
lmp_version -in DislocationVelocity.in
 +
</pre>
  
 
LAMMPS will create several output files from the simulation. You can load each of them in [[Ovito]] to see what each includes. You will run the post-processing script, [[Single Defect Velocity in Ovitos]], on the <code>dump.shear.unwrap</code> file, following the directions on the script page.
 
LAMMPS will create several output files from the simulation. You can load each of them in [[Ovito]] to see what each includes. You will run the post-processing script, [[Single Defect Velocity in Ovitos]], on the <code>dump.shear.unwrap</code> file, following the directions on the script page.
Line 94: Line 100:
 
==Environment Setup==
 
==Environment Setup==
  
Source code and example inputs can be downloaded from this [[Media:MDDP_HW2.zip|link]]. <br>
+
Source code and example inputs can be downloaded from this [[Media:MDDP_BCC_HW2.zip|link]]. <br>
Post processing instructions for MDDP are found [[MDDP Post Processing|here]].  
+
Post processing instructions for MDDP are found [[MDDP Post Processing|here]]. <br>
 +
A video tutorial on running MDDP can be found [https://youtu.be/3g_MX2bKDes here]. <br>
 +
A second video on running MDDP can be found [https://www.youtube.com/watch?v=VsqUBnpqJu0&feature=youtu.be here].
  
 
==Simulating the Motion of a Dislocation==
 
==Simulating the Motion of a Dislocation==
Line 113: Line 121:
 
*Shear modulus (shr)
 
*Shear modulus (shr)
 
*Poisson's ratio (both elastic moduli can be upscaled from DFT)
 
*Poisson's ratio (both elastic moduli can be upscaled from DFT)
*Mobility (from the first homework)
+
*Mobility (from the first part of this homework)
 
*Burger's vector (Ba)
 
*Burger's vector (Ba)
 
*Intrinsic stacking fault energy (stkFE)
 
*Intrinsic stacking fault energy (stkFE)
Line 132: Line 140:
  
 
====Step 4====
 
====Step 4====
Post-process using TecPlot or Paraview
+
Post-process using [[MDDP Post Processing|TecPlot]] or Paraview
  
 
===Homework Assignment===
 
===Homework Assignment===
1. Run [[MDDP|MDDP]] using the single Frank-Read source (SFRS) input. Be sure to change the data file to reflect the properties of your material as determined from LAMMPS.  
+
1. Run [[MDDP|MDDP]] using a single Frank-Read source (SFRS). Be sure to change the data file to reflect the properties of your material as determined from LAMMPS.  
:a. Generate stress-strain curves using a minimum of three (3) different mobilities.  
+
:a. Generate stress-strain curves using a minimum of three (3) different mobilities.  
:b. Illustrate the SFRS at several intervals as the dislocation loop propagates.
+
:b. Illustrate the SFRS at several intervals as the dislocation loop propagates.
 
+
2. Run [[MDDP|MDDP]] using the multiple Frank-Read sources (MFRS) input. Be sure to change the data file to reflect the properties of your material as determined from LAMMPS.
+
:a. Generate stress-strain curves using a minimum of three (3) different mobilities. These will be used for upscaling to crystal plasticity.
+
:b. Illustrate the MFRS at several intervals as the dislocation loops propagate.
+
  
 
{| class="wikitable" style="width: 40%; height: 100px; margin: 1em auto 1em auto;"
 
{| class="wikitable" style="width: 40%; height: 100px; margin: 1em auto 1em auto;"

Latest revision as of 21:40, 19 April 2019

< ICME 2017 Overview

Contents

[edit] Overview

In this homework, we will bridge information from the nanoscale to the microscale by calculating the dislocation mobility drag coefficient using Molecular Dynamics, and using it to run Dislocation Dynamics simulations. There are two parts to this homework,

  • Molecular Dynamics (MD) using Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS)
  • Dislocation Dynamics (DD) using Multiscale Dislocation Dynamics Plasticity (MDDP)

All the necessary files can be found on the ICME website. Copy them to a folder on your computer. Use /scratch/"Your Directory" for best results, especially if your job reads/writes much data.

Write a full report that follows a journal article manuscript format (include figures and tables in the text).

Upon completion, submit via email a .pdf and .doc(x) file of your report. Be sure to also include the requested files and plots from each section of the homework.

[edit] Part 1 - Single Dislocation Mobility Calculations using MEAM

[edit] Objective

For this section, you will use the Modified Embedded Atom Method (MEAM) in Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) to acquire the dislocation mobility drag coefficient for your material.

You can obtain LAMMPS and find the user manual here. Visit the LAMMPS tutorials page for a wide range of examples. LAMMPS video tutorials - Intro Input file Part 1 Input file Part 2 Molecular Dynamics Theory and Examples

[edit] Environment Setup

LAMMPS is an open source code and can be downloaded here. The source code can be easily compiled, or binary distributions are available for easy installation.

In addition to the software, you will need an input file, an atom position file, a post-processing script, and two files for your atomistic potential.

[edit] Calculating Dislocation Velocity

[edit] Walkthrough

Video tutorials: Part 1 Part 2

[edit] Step 1

Generate the atomic structure using Atomistic Dislocation Generation. Once you have downloaded and compiled the fortran routine, as directed, simply run it at the command line as explained on the script page. For this calculation, you want to make an edge dislocation in a PAD geometry. It will create the required data file as atoms.*.edge.pad, where * is either fcc or bcc, depending on your material. You will need to copy the data file to the directory where you will run the MD simulation.

[edit] Step 2

Edit the LAMMPS input file for your material, applied stress, and desired temperature. These can all be set by variables at the top of the input file.

# Variable definitions
variable initTemp equal 300.          	# desired temperature
variable sigma equal 15000.    		# applied stress in bar
variable material string Ta		# material symbol
variable atom_file string atoms.bcc.edge.pad # the configuration was generated by SG with the preprocessor dislocation.f90

initTemp will determine the temperature of the system.

sigma will determine the applied shear stress in bar.

material will determine the name of your material used in your MEAM potential files. It expects to find a library file, *.library.meam, and a parameter file, *.meam where the * represents the material symbol.

atom_file is the name of the data file that contains the initial atom positions.

[edit] Step 3

Run the simulation:

lmp_version -in DislocationVelocity.in

LAMMPS will create several output files from the simulation. You can load each of them in Ovito to see what each includes. You will run the post-processing script, Single Defect Velocity in Ovitos, on the dump.shear.unwrap file, following the directions on the script page.

It will output a time versus position file which you can use to calculate the velocity of the dislocation.

[edit] Homework Assignment

  1. Generate the atom positions file to be used for studying the mobility of an edge dislocation for your FCC or BCC material. A unit cell size of 100 x 60 x 2 will produce a simulation box containing ~70,000 atoms for an FCC structure.
  2. Run LAMMPS using the atom positions file generated in the previous step along with the LAMMPS input file for each of the following:
    1. Show the atom positions before the calculation illustrating the dislocation by looking at the dump.equilibration file.
    2. Use a minimum of three (3) different MEAM parameter sets based on the sensitivity analysis from HW1. Compare the position vs. time curves for each set.
    3. Study the effects of the applied shear stress on the dislocation velocity in your material compared to aluminum as in Figure 9.7 (a) in the ICME for Metals textbook.
    4. Determine the drag coefficient using Equation 9.2 in the ICME for Metals textbook from the study in Part (c).
Dislocation Mobility

[edit] Part 2 - Dislocation Dynamics

[edit] Objective

In this section, you will use Multiscale Dislocation Dynamics Plasticity (MDDP) to simulate the motion of a dislocation, thereby upscaling the dislocation mobility from the nanoscale.

[edit] Environment Setup

Source code and example inputs can be downloaded from this link.
Post processing instructions for MDDP are found here.
A video tutorial on running MDDP can be found here.
A second video on running MDDP can be found here.

[edit] Simulating the Motion of a Dislocation

[edit] Walkthrough

[edit] Step 1

Edit datain and data.

The first file, datain, will be used by the utility that creates the dislocation geometry. You will need to make sure that the two files are consistent - specifically that they both have the same system size (sidex,sidey,sidez).

The second file, data, will be used by MDDP to define material and system parameters.

For your material, edit:

  • Crystal structure
  • Density
  • Shear modulus (shr)
  • Poisson's ratio (both elastic moduli can be upscaled from DFT)
  • Mobility (from the first part of this homework)
  • Burger's vector (Ba)
  • Intrinsic stacking fault energy (stkFE)

Make sure that you maintain the format of the file

[edit] Step 2

Create the dislocation geometry.

Depending on your crystal structure, run FCCdata or BCCdata. The utility will prompt you for several quantities, including the type of structures to create. Specify Frank-Read Source pairs as your geometry.

[edit] Step 3

Run MDDP.

Several output files will be created, but only two kinds that we care about:

  • DDtimeResults.out - This is a text file with the history data of dislocation density, stress, and strain.
  • tech* - These will be a series of files that can be used to get the visual output including the dislocation lines.

[edit] Step 4

Post-process using TecPlot or Paraview

[edit] Homework Assignment

1. Run MDDP using a single Frank-Read source (SFRS). Be sure to change the data file to reflect the properties of your material as determined from LAMMPS.

a. Generate stress-strain curves using a minimum of three (3) different mobilities.
b. Illustrate the SFRS at several intervals as the dislocation loop propagates.
Frank Read Source Operation

[edit] Part 3 - Room for Improvement

Improve the instructions and/or tutorials for running LAMMPS/MDDP using your experience gained from Parts 1 and 2.

[edit] Additional Guides

ICME 2012 HW2

ICME 2013 HW2

ICME 2015 HW2

LAMMPS tutorials

MDDP Post Processing

[edit] License

By using the codes provided here you accept the the Mississippi State University's license agreement. Please read the agreement carefully before usage.

Personal tools
Namespaces

Variants
Actions
home
Materials
Material Models
Design
Resources
Projects
Education
Toolbox