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

Fixing error in averaging

Browse Source
This commit is contained in:
jtclemm
2023-12-01 09:58:40 -07:00
parent 1e26c6d0c5
commit f7aeecd3be
1 changed files with 2 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
src/RHEO/pair_rheo.cpp
View File

@ -248,7 +248,7 @@ void PairRHEO::compute(int eflag, int vflag)
if (harmonic_means_flag) {
kappa_ave = 2.0 * kappai * kappaj / (kappai + kappaj);
} else {
kappa_ave = 0.5 * (kappai * kappaj);
kappa_ave = 0.5 * (kappai + kappaj);
}
dT_prefactor = 2.0 * kappa_ave * (Ti - Tj) * rinv * rinv * voli * volj * 2.0 / (rhoi + rhoj);
@ -270,7 +270,7 @@ void PairRHEO::compute(int eflag, int vflag)
if (harmonic_means_flag) {
eta_ave = 2.0 * etai * etaj / (etai + etaj);
} else {
eta_ave = 0.5 * (etai * etaj);
eta_ave = 0.5 * (etai + etaj);
}
//Add artificial viscous pressure if required