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 |
|
|