Linear Elastic model enabled

Author: divyaadil23

Code/Source/svFSI/fsi.cpp

@@ -38,6 +38,7 @@
 #include "nn.h"
 #include "sv_struct.h"
 #include "utils.h"
+#include "l_elas.h"
 
 #include <array>
 #include <iomanip>
@@ -224,11 +225,12 @@ void construct_fsi(ComMod& com_mod, CepMod& cep_mod, const mshType& lM, const Ar
             auto N0 = fs_1[0].N.col(g);
             struct_ns::struct_3d(com_mod, cep_mod, fs_1[0].eNoN, nFn, w, N0, Nwx, al, yl, dl, bfl, fN, pS0l, pSl, ya_l, lR, lK);
           } break;
-          case Equation_lElas:
-            throw std::runtime_error("[construct_fsi] LELAS3D not implemented");
+          case Equation_lElas:{
+            // throw std::runtime_error("[construct_fsi] LELAS3D not implemented");
             //CALL LELAS3D(fs(1).eNoN, w, fs(1).N(:,g), Nwx, al, dl, bfl, pS0l, pSl, lR, lK)
-          break;
-
+            auto N0 = fs_1[0].N.col(g);
+            l_elas::l_elas_3d(com_mod, fs_1[0].eNoN, w, N0, Nwx, al, dl, bfl, pS0l, pSl, lR, lK);
+          } break;
           case Equation_ustruct:
             throw std::runtime_error("[construct_fsi] USTRUCT3D_M not implemented");
             //CALL USTRUCT3D_M(vmsStab, fs(1).eNoN, fs(2).eNoN, nFn, w, Jac, fs(1).N(:,g), fs(2).N(:,g), Nwx, al, yl, 

tests/cases/fsi/pipe_3d_LinElastic/README.md

@@ -0,0 +1,19 @@
+
+# **Problem Description**
+
+Simulate pressure wave propagation in an arterial model using the Arbitrary Lagrangian-Eulerian method [1]. The problem set-up is as follows.
+
+<p align="center">
+   <img src="./configuration.png" width="600">
+</p>
+
+And the results are
+
+<p align="center">
+   <img src="./results.gif" width="600">
+</p>
+
+
+## References
+
+1.  Liu, Ju, and Alison L. Marsden.  A Unified Continuum and Variational Multiscale Formulation for Fluids, Solids, and Fluid Structure Interaction.  *Computer Methods in Applied Mechanics and Engineering* 337 (August 2018): 549 97. https://doi.org/10.1016/j.cma.2018.03.045.

tests/cases/fsi/pipe_3d_LinElastic/configuration.png

@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:bff21b54094fe528e6f7cd9b57dafc5863b530e112e08f8a62a04431ec14a7a3
+size 67969

tests/cases/fsi/pipe_3d_LinElastic/mesh/fluid/mesh-complete.mesh.vtu

@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:57511367796bab30c07b3b5dbd3b6ec596e433cb1aae2afe763f05ab357655bf
+size 13661

tests/cases/fsi/pipe_3d_LinElastic/mesh/fluid/mesh-surfaces/end.vtp

@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:4c4a25e331d2c3ccbfa2ffcdd52e0ebbd37a8045989bec9abd7b2c0ed775e2b8
+size 28724

tests/cases/fsi/pipe_3d_LinElastic/mesh/fluid/mesh-surfaces/interface.vtp

@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:d90d5fb5ce21d0889b27fd0a592aafa6121d19047bebd801710caf5cba6fb48a
+size 13685

tests/cases/fsi/pipe_3d_LinElastic/mesh/fluid/mesh-surfaces/start.vtp

@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:2527472b08ece18c95175c0661ca838f90ac8b97406d0c1c77459e94df620fe7
+size 36980

tests/cases/fsi/pipe_3d_LinElastic/mesh/solid/mesh-complete.mesh.vtu

@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:abb10cc3c374742b481031050fd6ceb2af30b8a14bd5627d89f930a4df01708a
+size 12351

tests/cases/fsi/pipe_3d_LinElastic/mesh/solid/mesh-surfaces/end.vtp

@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:9302a80f392f30ae07b98d0048cb7ad85f2fd3fcb0a86dbe50def69ba5a6cafc
+size 28900