Hi Lammps users
I have simulated shock wave in Lammps with fix wall/atoms and fix append/atoms.
I took the length of thermostated region of fix append/atoms equal to the length of appended region and my fix append/atoms, fix wall/piston and fix nve commands were:
fix walls all wall/piston zlo vel 0.0 0.0 ${VLammps} temp 1.0 10 589999 $l units box
fix 1 all nve
fix 2 all append/atoms zhi size y freq {FREQ} temp $T 100 579946 $y units lattice
But when I saw the temperature profile it seems that the shock on the border of appended region was splited into two waves: one in direction of original wave and the other on the opposite direction. And when I visualized my results in AtomEye when one region is appended, the previously appended region cross section shrinks and it’s seems that the newly appended region have zero temperature.
Can you help me find out what the problem is?
Any help is appreciated in advance
Mahdi Tavakol