ICODE-ADC/lammps
4
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Update smd_material_models.cpp

Browse Source 
The pressure pFinal using ShockEOS() was of the wrong sign when Gamma = 0.
This commit is contained in:
abbatux
2017-10-24 11:04:11 +11:00
committed by GitHub
parent 39df9f5d94
commit 41c25877e8
1 changed files with 1 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
src/USER-SMD/smd_material_models.cpp
View File

@ -75,7 +75,7 @@ void ShockEOS(double rho, double rho0, double e, double e0, double c0, double S,
double mu = rho / rho0 - 1.0;
double pH = rho0 * square(c0) * mu * (1.0 + mu) / square(1.0 - (S - 1.0) * mu);
pFinal = (pH + rho * Gamma * (e - e0));
pFinal = (-pH + rho * Gamma * (e - e0));
//printf("shock EOS: rho = %g, rho0 = %g, Gamma=%f, c0=%f, S=%f, e=%f, e0=%f\n", rho, rho0, Gamma, c0, S, e, e0);
//printf("pFinal = %f\n", pFinal);