chimad hackathon 2: team mcgill
TRANSCRIPT
![Page 1: CHiMaD Hackathon 2: Team mcgill](https://reader036.vdocument.in/reader036/viewer/2022081605/5872f2361a28ab8c718b4cd3/html5/thumbnails/1.jpg)
TEAM 2N. SMITHT. PINOMAAK.N. SHAMPURMcGill
CHiMaD Phase Field Methods Workshop III & HackathonMAY 3-5, 2016
![Page 2: CHiMaD Hackathon 2: Team mcgill](https://reader036.vdocument.in/reader036/viewer/2022081605/5872f2361a28ab8c718b4cd3/html5/thumbnails/2.jpg)
Question #no. 1.
![Page 3: CHiMaD Hackathon 2: Team mcgill](https://reader036.vdocument.in/reader036/viewer/2022081605/5872f2361a28ab8c718b4cd3/html5/thumbnails/3.jpg)
CHiMaD Phase Field Methods Workshop III & Hackathon – MAY 3-5, 2016
Q1 – GIF of EvolutionPF contour superimposed on U field
![Page 4: CHiMaD Hackathon 2: Team mcgill](https://reader036.vdocument.in/reader036/viewer/2022081605/5872f2361a28ab8c718b4cd3/html5/thumbnails/4.jpg)
CHiMaD Phase Field Methods Workshop III & Hackathon – MAY 3-5, 2016
What final mesh would look like
![Page 5: CHiMaD Hackathon 2: Team mcgill](https://reader036.vdocument.in/reader036/viewer/2022081605/5872f2361a28ab8c718b4cd3/html5/thumbnails/5.jpg)
CHiMaD Phase Field Methods Workshop III & Hackathon – MAY 3-5, 2016
Q1 – Tip Velocity vs Time (part a) Vnum – Interface position per time stepVsharp - Tracking interface, then taking local numerical gradients in liquid and solid phase.
• V_num = 0.089• V_sharp =
0.031
• V_num = 0.014• V_sharp =
0.018
• V_num = 0.027• V_sharp =
0.025L = 20 L = 100 L = 200
![Page 6: CHiMaD Hackathon 2: Team mcgill](https://reader036.vdocument.in/reader036/viewer/2022081605/5872f2361a28ab8c718b4cd3/html5/thumbnails/6.jpg)
CHiMaD Phase Field Methods Workshop III & Hackathon – MAY 3-5, 2016
Q1 – Details
• dx = 20/32 for L = 20, 100/128 for L = 100• dt < dx*dx/(4*D) (linear stability criterion)• # cores: 1• Wall clock time: 10 min for L = 100, <1m for L
= 20;• Intel(R) Xeon(R) CPU E5620 @
2.40GHz• Total 65 GB RAM• https://www.dropbox.com/sh/pyolv5zju1fheqd/
AACHZf_6RcTs83nsF4mxO4mya?dl=0
![Page 7: CHiMaD Hackathon 2: Team mcgill](https://reader036.vdocument.in/reader036/viewer/2022081605/5872f2361a28ab8c718b4cd3/html5/thumbnails/7.jpg)
Question 2
![Page 8: CHiMaD Hackathon 2: Team mcgill](https://reader036.vdocument.in/reader036/viewer/2022081605/5872f2361a28ab8c718b4cd3/html5/thumbnails/8.jpg)
CHiMaD Phase Field Methods Workshop III & Hackathon – MAY 3-5, 2016
Q2 – Energy vs “Time” (Plot)
• Stress free • Cahn-Hilliard dynamics
final free energy = 1.2681020e-14 J
![Page 9: CHiMaD Hackathon 2: Team mcgill](https://reader036.vdocument.in/reader036/viewer/2022081605/5872f2361a28ab8c718b4cd3/html5/thumbnails/9.jpg)
CHiMaD Phase Field Methods Workshop III & Hackathon – MAY 3-5, 2016
Q2 – 3D representation of result
Initialized cube relaxing to a sphere
![Page 10: CHiMaD Hackathon 2: Team mcgill](https://reader036.vdocument.in/reader036/viewer/2022081605/5872f2361a28ab8c718b4cd3/html5/thumbnails/10.jpg)
CHiMaD Phase Field Methods Workshop III & Hackathon – MAY 3-5, 2016
Q2 – Details
• Implemented in Julia• Iterations to converge: 600• Cores: 4 (8 threads)• CPU: i7-4710MQ• Wall clock time: 28 seconds• Source available at
github.com/nsmith/chimadQ2
![Page 11: CHiMaD Hackathon 2: Team mcgill](https://reader036.vdocument.in/reader036/viewer/2022081605/5872f2361a28ab8c718b4cd3/html5/thumbnails/11.jpg)
CHiMaD Phase Field Methods Workshop III & Hackathon – MAY 3-5, 2016
thoughts• U field reaches boundary very fast• Large dx leads to not so round initial
seed• Rotated aniso & n-fold symmetry
were buggy! (in house code)• Not familiar with elasticity models
![Page 12: CHiMaD Hackathon 2: Team mcgill](https://reader036.vdocument.in/reader036/viewer/2022081605/5872f2361a28ab8c718b4cd3/html5/thumbnails/12.jpg)
CHiMaD Phase Field Methods Workshop III & Hackathon – MAY 3-5, 2016
Thanks.