Skip to content

compute stress

New to LBM and OpenLB? Attend our upcoming one-week hands-on Spring School in Marseille/France from 19.05.2025 – 23.05.2025.

OpenLB – Open Source Lattice Boltzmann Code › Forums › on OpenLB › General Topics › compute stress

  • This topic has 2 replies, 2 voices, and was last updated 1 month ago by liu.
Viewing 3 posts - 1 through 3 (of 3 total)
  • Author
    Posts
  • February 13, 2025 at 9:50 am #9876
    liu
    Participant

    dear community:
    In the three-dimensional case, what does pi[2] represent, and in the two-dimensional case, what does pi[2] represent? thanks

    template <typename CELL, typename V=typename CELL::value_t>
    static V computePiNeqNormSqr(CELL& cell) any_platform
    {
    V rho, u[DESCRIPTOR::d], pi[util::TensorVal<DESCRIPTOR>::n];
    computeAllMomenta(cell, rho, u, pi);
    V PiNeqNormSqr = pi[0]*pi[0] + 2.*pi[1]*pi[1] + pi[2]*pi[2];
    if constexpr (util::TensorVal<DESCRIPTOR >::n == 6) {
    PiNeqNormSqr += pi[2]*pi[2] + pi[3]*pi[3] + 2.*pi[4]*pi[4] +pi[5]*pi[5];
    }
    return PiNeqNormSqr;
    }

    February 13, 2025 at 9:58 am #9877
    Adrian
    Keymaster

    pi[2] is the xz component of the strain rate tensor in 3D. In 2D pi[2] is the yy component.

    February 13, 2025 at 1:44 pm #9879
    liu
    Participant

    thank you

  • Author
    Posts
Viewing 3 posts - 1 through 3 (of 3 total)
  • You must be logged in to reply to this topic.