Hello,

today a new problem has appeared.

Usually, the shear stress tensor Pi is calculated in the form

```
Pi_ij = eta * (dv_i / dx_j + dv_j / dx_i)
```

Obviously, this expression is symmetric. It can be calculated in the LB code in a local way,

```
Pi_ij = (1 - 1 / (2 * tau)) * sum_k f^neq_k * (c_ik * c_jk - delta_ij / 3)
```

in three dimensions. However, this quantity is not physical. Why is that so? If you have an area element in the xy-plane and you want to measure (in an experiment) the force due to shear on that area element PARALLEL to the xy-plane, then you measure

```
Pi'_xz = eta * dv_z / dx
```

and

```
Pi'_yz = eta * dv_z / dy
```

Those two quantities are not the same as Pi_xz and Pi_yz using the definition at the top. The shear tensor Pi_ij has five independent components, since it is symmetric and traceless. The shear tensor Pi’_ij has EIGHT independent components (if the flow is incompressible). The missing three components are the asymmetric parts, which have been removed.

Now, since I want to simulate the physical shear stress tensor, which can be measured in an experiment, I want to calculate the quantity Pi’ and not Pi. It is not possible to extract Pi’ from Pi, since the asymmetric information is gone.

But is it possible to find a local expression usable in a LB code, which directly produces Pi’?

My understanding is that this is not possible, since one cannot find a second moment of the velocities which is asymmetric.

My conclusion is that, if one is interested in Pi’, one has to use a finite difference scheme, which of course is not local.

Does anybody have an idea regarding this issue?

Thanks a lot,

Timm