PARADISCYL:Cylinder-test

Test Cases in cylinder code

ill Ryu and Wei Cai

This tutorial describes test cases to check if cylinder code works well, especially for the image stress calculation. We provide Matlab files which generates ParaDiS inputs and plot the results. The theoretical background is published in Computing Image Stress in an Elastic Cylinder(Journal of the Mechanics and Physics of Solids, 55, 2027 (2007)(PDF)) and Dislocation Dynamics Simulations in a Cylinder(Proceedings of the Dislocations 2008 International Conference, IOP Conference Series: Materials Science and Engineering, vol 3,012007 (2009).(PDF)).

---

TEST1 : Straight edge dislocation along the cylinder axis

In the first test case, image stress of a straight edge dislocation is considered, as shown in the Figure 1, where dislocation is offset($x0$) from the center of the cylinder. We have two cases of $x0$=0.5$R$ and 0.9$R$.Here, $R$ is the radius of the cylinder.

(a)	(b)
Fig.1 (a) Isoview (b) Topview.The Bergurs vector is $\mathbf {b} =a[100]$. Since the image stress field of this configuration had been solved analytically[Eshelby, J.D., 1979],we compare the image force between the simulation result and analytic solution.The number of grids in circumferential direction is same as the one along the cylinder axis.To consider periodic images along the cylinder axis under PBC, analytical stress field of edge dislocation is implemented in the function of "AllSegmentStress_no_cell_test1".To use this function, modify makefile in ParaDiS/cylinder directory so that the following line is active. DEFS += -D_CYL_TEST1 Using Method I,the relative error in the image force is computed from the image stress calcualtion. Right plot clearly show exponential decay of the relative error with increasing number of grids. The convergence is slower when the location of the edge dislocation is close to the surface. TEST2 : Circular prismatic dislocation loop In the second test case, we consider a circular prismatic dislocation loop concentric with the cylindrical axis, discretized with 12 equally spaced nodes connected by straight segments, as shown in Fig. 3, where the dislocation loop radius is 0.8 times to the cylinder radius. (a) (b) Fig.3 (a) Isoview (b) Topview.The Bergurs vector is $\mathbf {b} =a[001]$. For simplicity, cell is not used for the force calculatoin. To do this, modify makefile in ParaDiS/cylinder directory so that the following line is active. DEFS += -D_CYL_TEST23 Due to the symmetry of this problem, the image forces on all nodes point to the radial direction and have the same magnitude. An estimate of their relative error is plotted as a function of number of grid in the circumferential direction($n_{\theta }$) in the following figure. Because this problem does not have an analytic solution, the reference value is taken to be the value obtained using Method I (Bessel) with $n_{\theta }=181$. As expected, the relative error in radial image force decreases with increasing number of grids. TEST3 : Circular inclined glide dislocation loop In the third test case, we consider more general shaped dislocation loop, discretized with 10 equally spaced nodes connected by straight segments, as shown in Fig. 5. (a) (b) Fig.5 (a) Isoview (b) Topview.The Bergurs vector is $\mathbf {b} =\mathbf {a} /{\sqrt {2}}[{\bar {1}}10]$ and the slip plane is $\mathbf {n} =1/{\sqrt {6}}[{\bar {1}}{\bar {1}}2]$ Likewise, cell is not used for the force calculatoin. To do this, modify makefile in ParaDiS/cylinder directory so that the following line is active. DEFS += -D_CYL_TEST23 The orientation and magnitude of the image force on every node is different. Fig. 10 plots the component of the image force in the direction away from the center of the loop for all 10 nodes with different values of nq using Method I. With more than 45 number of grids, they converge to almost same value.