Skip to content

Bug in HLLC Jacobian calculation #2198

@ashwaniassam

Description

@ashwaniassam

https://github.com/su2code/SU2/blob/60270690cd749bd190e16b4c34cf724c42edbeb1/SU2_CFD/src/numerics/flow/convection/hllc.cpp#L543C3-L543C3

As per the https://www.sciencedirect.com/science/article/pii/S0021999197957934 (original paper), multiplying the Jacobian with \kappa = 0.5 is not right. It is leading to non-convergence.

Further, dpStar_dU[iVar] calculations need to be relooked. As per the original paper, density should be at "j" when it is being used at "i". Though I have not seen a major difference, but currently testing the same more for strong shock cases.

Metadata

Metadata

Assignees

No one assigned

    Labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions