Dear LAMMPS/ATC users:
I guess my questions should be addressed to Drs. Zimmerman, Jones, Templeton or other AtC folks.
I’m playing with AtC package trying to find the pros and cons of using a coupled FE/MD simulation for predicting the mechanical and thermal properties of solids. Most of the tests and input scripts provided with the package work fine and I am truly impressed with the ATC code, but sometimes I face with lack of some useful functionalities. Here are two examples:
-
After a few days of playing around with the simple input files found in the examples/USER/atc/thermal sub-directory of the LAMMPS public release (10 Aug 2015), I was confident that I could produce any script with minimal effort. Unfortunately I cannot visualize my results when I use the thermal fix and the binary output (Ensight *.case file): the temperature field is present in the text .DATA file, however, it is not displayed by Paraview. How can I sort out this problem? In particular, how can I generate with Paraview/EnSight software the temperature fields like those shown in Fig.7 in the paper of Wagner et al. (Comput. Methods Appl. Mech. Engrg. 2008, 197, 3351)?
-
My second question is related to the in.eshelby_static input file from the examples/USER/atc/hardy sub-directory. I understand that the FE region can be larger than the MD region, but how can I visualize the coarse-grained stress fields outside the MD domain, not only inside it? Rephrasing my question, is it possible to analyze the whole FE+MD domain when FE >> MD? In this regard I’d like to refer to a nice picture from Zimmerman’s presentation (2013) where J-integral calculations are discussed (see e.g. slide #9). I would be grateful for examples of the corresponding input scripts.
Thanks and best regards,
Helmut Grauer
The University of Ulm
89069 Ulm,
Germany