Case 1: One threading dislocation in the thinfilm code

Benchmarks to check PBC in thin film code

ILL RYU

May 31, 2010

From these benchmarks, we like to confirm if the periodic boundary condition works properly. To do that, two models are prepared as following figures. In the first model, one threading dislocation is located at the center of the simulation domain. In the other model, threading dislocation is split and is put into both right and the left surface.

Model geometries

• 

  [One threading dislocaiton at the center]

• 

  [Splited dislocaitons at the right and left surfaces]

To check if PBC works, forces at the nodes are printed out in each cases. If forces output in both cases are the same, we can conclude PBC works well in the thinfilm code.

Test cases

In the thinfilm code, there are several flags which can turn on'off the spectral method (DEFS += -D_TFIMGSTRESS), with or without virtual segments (DEFS += -D_NOVIRTUALSEG) and with or without the Yoffe image stress (DEFS += -D_NOYOFFESTRESS). In order to test these different cases, force outputs are compared as the following table.

Test cases	Spectral method	Virtural segment	Yoffe stress
1.threading_force_NoImg_NoVS_NoYoffe	OFF	OFF	OFF
2.threading_force_NoImg_NoVS_YesYoffe	OFF	OFF	ON
3.threading_force_NoImg_YesVS_NoYoffe	OFF	ON	OFF
4.threading_force_NoImg_YesVS_YesYoffe	OFF	ON	ON
5.threading_force_YesImg_NoVS_NoYoffe	ON	OFF	OFF
6.threading_force_YesImg_YesVS_YesYoffe	ON	ON	ON

RESULTS

case 1:threading_force_NoImg_NoVS_NoYoffe

• 

case 2: threading_force_NoImg_NoVS_YesYoffe

• 

case 3:threading_force_NoImg_YesVS_NoYoffe

• 

case 4:threading_force_NoImg_YesVS_YesYoffe

• 

case 5:threading_force_YesImg_NoVS_NoYoffe

• 

case 6:threading_force_YesImg_YesVS_YesYoffe

• 

Conclusion

From the plot, we can conclude that PBC works well in thinfilm code, because the force output fits well to each other.

Input & plot files

Media: TBD
Media: TBD