This page is intended to developers of the ParaDiS code. We share a few programming tips we discover in ParaDiS.
The Peach-Koehler force at a given node of a dislocation coming from the stress of a dislocation segment σ is given by
The way this force is implemented in ParaDiS is as follows (see ComputeSegSigbRem in NodeForce.c). Consider two neighboring nodes Node 1 and Node 2 where we want to compute the force. Suppose these nodes have coordinates x1 and x2 respectively.
is computed for the pair Node 1 and Node 2. Each node gets the same , where here b1 is the Burgers vector of Node 1. Note that nodes are ordered in ParaDiS and only the pair (Node 1, Node 2) is considered and not (Node 2, Node 1) in this case. 2. This dot product is then use to calculate the Peach-Koehler force. Node 1 gets the force and Node 2 gets the force .