Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Bug fixes for ustruct Guccione model and follower pressure load #201

Merged
merged 5 commits into from
Mar 6, 2024
Merged

Bug fixes for ustruct Guccione model and follower pressure load #201

Show file tree
Hide file tree
Changes from 1 commit
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
751b966
Bug fixes for ustruct Gucci and follower pressure
aabrown100-git Feb 28, 2024
b70edc0
Merge branch 'SimVascular:main' into ustruct_bug_fix
aabrown100-git Mar 1, 2024
4e3da2d
Reverting ustruct/LV_Guccione_active max_iterations
aabrown100-git Mar 1, 2024
91d1945
Updating LV_Guccione_active reference solution
aabrown100-git Mar 1, 2024
93ba1ef
Replacing Sb = Sb + with Sb +=
aabrown100-git Mar 6, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
4 changes: 2 additions & 2 deletions Code/Source/svFSI/mat_models.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -690,10 +690,10 @@ void get_pk2cc_dev(const ComMod& com_mod, const CepMod& cep_mod, const dmnType&
2.0*g3*(Es(0,1)*RmRm.slice(3) + Es(0,2)*RmRm.slice(5));

auto CCb = 2.0*ten_dyad_prod(Sb, Sb, nsd);
Sb += Sb * r2;
Sb = Sb * r2;

// Fiber reinforcement/active stress
Sb += Sb + Tfa*mat_dyad_prod(fl.col(0), fl.col(0), nsd);
Sb = Sb + Tfa*mat_dyad_prod(fl.col(0), fl.col(0), nsd);
aabrown100-git marked this conversation as resolved.
Show resolved Hide resolved

double r1 = J2d*mat_ddot(C, Sb, nsd) / nd;
S = J2d*Sb - r1*Ci;
Expand Down
2 changes: 1 addition & 1 deletion Code/Source/svFSI/ustruct.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -142,7 +142,7 @@ void b_ustruct_3d(const ComMod& com_mod, const int eNoN, const double w, const V
for (int a = 0; a < eNoN; a++) {
NxFi(0,a) = Nx(0,a)*Fi(0,0) + Nx(1,a)*Fi(1,0) + Nx(2,a)*Fi(2,0);
NxFi(1,a) = Nx(0,a)*Fi(0,1) + Nx(1,a)*Fi(1,1) + Nx(2,a)*Fi(2,1);
NxFi(2,a) = Nx(1,a)*Fi(1,2) + Nx(1,a)*Fi(1,2) + Nx(2,a)*Fi(2,2);
NxFi(2,a) = Nx(0,a)*Fi(0,2) + Nx(1,a)*Fi(1,2) + Nx(2,a)*Fi(2,2);
}

nFi(0) = nV(0)*Fi(0,0) + nV(1)*Fi(1,0) + nV(2)*Fi(2,0);
Expand Down
2 changes: 1 addition & 1 deletion tests/cases/ustruct/LV_Guccione_active/svFSI.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -39,7 +39,7 @@
<Add_equation type="ustruct" >
<Coupled> true </Coupled>
<Min_iterations> 1 </Min_iterations>
<Max_iterations> 19 </Max_iterations>
<Max_iterations> 10 </Max_iterations>
<Tolerance> 1e-12 </Tolerance>
<Use_taylor_hood_type_basis> false </Use_taylor_hood_type_basis>

Expand Down