PARADISCYL:Cylinder-test: Difference between revisions
No edit summary |
No edit summary |
||
| Line 13: | Line 13: | ||
==TEST1 : Straight edge dislocation along the cylinder axis== |
==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(<math> |
In the first test case, image stress of a straight edge dislocation is considered, as shown in the ''Figure 1'', where dislocation is offset(<math>x0</math>) from the center of the cylinder. We have two cases of <math>x0</math>=0.5<math>R</math> and 0.9<math>R</math>.Here, <math>R</math> is the radius of the cylinder. |
||
{|border="0" align="center" |
{|border="0" align="center" |
||
|[[Image:T1_geom_Isoview.jpg|frameless|400px|caption]] |
|[[Image:T1_geom_Isoview.jpg|frameless|400px|caption]] |
||
Revision as of 10:05, 17 November 2011
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(Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle x0} ) from the center of the cylinder. We have two cases of =0.5Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle R} and 0.9Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle R} .Here, Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle R} is the radius of the cylinder.
|
|
| (a) | (b) |
| Fig.1 (a) Isoview (b) Topview.The Bergurs vector is Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \mathbf{b} = a[1 0 0]}
.
Since the image stress field of this configuration is 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. |
|