| 87 |
|
HydroProps cr; |
| 88 |
|
HydroProps cd; |
| 89 |
|
calcHydroPropsAtCR(beads_, viscosity, temperature, cr); |
| 90 |
< |
calcHydroPropsAtCD(beads_, viscosity, temperature, cd); |
| 90 |
> |
//calcHydroPropsAtCD(beads_, viscosity, temperature, cd); |
| 91 |
|
setCR(cr); |
| 92 |
|
setCD(cd); |
| 93 |
|
|
| 156 |
|
|
| 157 |
|
Xiott += Cij; |
| 158 |
|
Xiotr += U[i] * Cij; |
| 159 |
< |
Xiorr += -U[i] * Cij * U[j] + (6 * viscosity * volume) * I; |
| 159 |
> |
//Xiorr += -U[i] * Cij * U[j] + (6 * viscosity * volume) * I; |
| 160 |
> |
Xiorr += -U[i] * Cij * U[j]; |
| 161 |
|
} |
| 162 |
|
} |
| 163 |
|
|
| 320 |
|
|
| 321 |
|
Xitt += Cij; |
| 322 |
|
Xitr += U[i] * Cij; |
| 323 |
< |
Xirr += -U[i] * Cij * U[j] + (6 * viscosity * volume) * I; |
| 323 |
> |
//Xirr += -U[i] * Cij * U[j] + (6 * viscosity * volume) * I; |
| 324 |
> |
Xirr += -U[i] * Cij * U[j]; |
| 325 |
|
} |
| 326 |
|
} |
| 327 |
|
|