# shock simulation

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

&

Regards

Oxem Selesta

shock_result_animation.avi (3.33 MB)

Selesta,

I don't understand your question, if there even was one. This has
become a non LAMMPS-related question, and the mailing list is not
suitable for interpreting your shock results.

Anyhow, your velocity unit is wrong - default velocity is in A/fs for
real units, not in A/ps. You were actually using a 2000 km/s piston
velocity. Might explain your run ended with a seg fault?

Ray

Dear Ray,
Thanks for your reply,I had used metal unit, so my velocity unit is A/ps not default A/fs. Here particle velocity was 2km/s.

The shock front was propagating from left to right(along +ve X dirn) but before reaching the right-wall, it seems that it have stopped, in this time the last 20% of the material near the right wall start motion towards -ve X dirn.