Dear Ray Shan and lammps users,
I am attaching the video as well as plot (velocity_time_ region wise) got from my shock simulation through porous medium , I was trying to use your logic to build the script initially. Thanks for helping me by sharing the script , from that I really have learnt how to use wall/reflect to generate a shock in classical method.
After a long steps my system behaves in peculiar way, the last region behind the shock front approaches opposite to the shock direction while the system dimension along shock direction increases gradually, finally my simulation stops showing error related to memory allocation.
my particle velocity was 2km/s. I had used your logic almost directly only change was system topology and force field information,here my script portion
variable SHOCKVEL equal 2 # km/s
variable xleftwall equal -0.0749979
variable xrightwall equal 812.925
variable Nshock equal 1.10*(v_xrightwall/0.00005)/(v_SHOCKVEL/0.10)
.
.
.
fix 4 all wall/reflect xlo {xleftwall} xhi {xrightwall} units box
.
variable vel equal ${SHOCKVEL}*-10.000 # A/ps
.
velocity all set {vel} 0.0 0.0 sum yes units box fix 1 all nve fix 4 all wall/reflect xlo {xleftwall} xhi ${xrightwall} units box
run ${Nshock}
Why the shock front cant reach the right wall before that the last portion behaves in such a manner. Please need your valuable comment. Please give some insights into this regarding the video and as well as plot.
Thanks in advance.
Thanks
&
Regards
Oxem Selesta
shock_result_animation.avi (3.33 MB)