From e0dbf8daf27a4779411dbe9f0213226f46d890cf Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Fabien=20P=C3=A9an?= Date: Fri, 22 Nov 2024 11:00:28 +0100 Subject: [PATCH 1/2] Fix speed of sound calculation for elastic solids --- src/shared/materials/elastic_solid.cpp | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/shared/materials/elastic_solid.cpp b/src/shared/materials/elastic_solid.cpp index 6776a6cb23..734a5f85b6 100644 --- a/src/shared/materials/elastic_solid.cpp +++ b/src/shared/materials/elastic_solid.cpp @@ -10,8 +10,8 @@ namespace SPH //=================================================================================================// void ElasticSolid::setSoundSpeeds() { - c0_ = sqrt(K0_ / rho0_); - ct0_ = sqrt(E0_ / rho0_); + c0_ = sqrt((K0_ + 4.*G/3.) / rho0_); // https://en.wikipedia.org/wiki/Speed_of_sound#Three-dimensional_solids + ct0_ = sqrt(E0_ / rho0_); // ⚠️ Calculation valid only for 1D solid https://en.wikipedia.org/wiki/Speed_of_sound#One-dimensional_solids cs0_ = sqrt(G0_ / rho0_); }; //=================================================================================================// From 546c251f9d73e14d0239f4be82b95a1d5caaf901 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Fabien=20P=C3=A9an?= Date: Fri, 22 Nov 2024 12:09:16 +0100 Subject: [PATCH 2/2] Fix typo --- src/shared/materials/elastic_solid.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/shared/materials/elastic_solid.cpp b/src/shared/materials/elastic_solid.cpp index 734a5f85b6..c2eb5aa5c9 100644 --- a/src/shared/materials/elastic_solid.cpp +++ b/src/shared/materials/elastic_solid.cpp @@ -10,7 +10,7 @@ namespace SPH //=================================================================================================// void ElasticSolid::setSoundSpeeds() { - c0_ = sqrt((K0_ + 4.*G/3.) / rho0_); // https://en.wikipedia.org/wiki/Speed_of_sound#Three-dimensional_solids + c0_ = sqrt((K0_ + 4.*G0_/3.) / rho0_); // https://en.wikipedia.org/wiki/Speed_of_sound#Three-dimensional_solids ct0_ = sqrt(E0_ / rho0_); // ⚠️ Calculation valid only for 1D solid https://en.wikipedia.org/wiki/Speed_of_sound#One-dimensional_solids cs0_ = sqrt(G0_ / rho0_); };