109 |
|
straightforward. The major modifications to the method are the addition of a rotational shearing term and the use of more versatile hot / cold regions instead |
110 |
|
of rectangular slabs. A temperature profile along the $r$ coordinate is created by recording the average temperature of concentric spherical shells. |
111 |
|
|
112 |
< |
\begin{figure} |
113 |
< |
\center{\includegraphics[width=7in]{figures/VSS}} |
114 |
< |
\caption{Schematics of periodic (left) and nonperiodic (right) Velocity Shearing and Scaling RNEMD. A kinetic energy or momentum flux is applied from region B to region A. Thermal gradients are depicted by a color gradient. Linear or angular velocity gradients are shown as arrows.} |
115 |
< |
\label{fig:VSS} |
116 |
< |
\end{figure} |
112 |
> |
% \begin{figure} |
113 |
> |
% \center{\includegraphics[width=7in]{figures/VSS}} |
114 |
> |
% \caption{Schematics of periodic (left) and nonperiodic (right) Velocity Shearing and Scaling RNEMD. A kinetic energy or momentum flux is applied from region B to region A. Thermal gradients are depicted by a color gradient. Linear or angular velocity gradients are shown as arrows.} |
115 |
> |
% \label{fig:VSS} |
116 |
> |
% \end{figure} |
117 |
|
|
118 |
|
At each time interval, the particle velocities ($\mathbf{v}_i$ and |
119 |
|
$\mathbf{v}_j$) in the cold and hot shells ($C$ and $H$) are modified by a |