Blood flow simulations in unicorn

Asked by Martin Vymazal

Dear FEniCS developers,

   the domain of my interest are blood flow simulations in bypass geometries. I would like to ask if the solver Unicorn can currently be used for fluid-structure interaction simulations. Could you also tell me if it is possible to prescribe flux boundary conditions for the fsi simulation and if it is possible to use meshes generated by gmsh?
 There used to be a bit of information about Unicorn before the main FEniCS website was attacked. I would appreciate any links to general description of Unicorn, papers, tutorials etc. Thank you.

Best regards,

   Martin Vymazal

Question information

Language:
English Edit question
Status:
Answered
For:
DOLFIN Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:
Revision history for this message
Jeannette Spühler (j-spuehler) said :
#1

Hi Martin,

Yes, you can use Unicorn to simulate FSI problems.
Please look at the examples you can find under unicorn/ucsolver/fsi to see how to set boundary conditions.

You have to convert the gmsh in to DOLFIN XML format. Please check the following link:
http://manpages.ubuntu.com/manpages/lucid/man1/dolfin-convert.1.html

There is no tutorial about unicorn which can be found online.
However they are working with a chapter about Unicorn for the FEniCS book.

What kind of problem do you want to solve? In my project we are simulating the blood flow in the heart. So I might help you in more details if you need.

Best regards
Jeannette Spühler

Revision history for this message
Martin Vymazal (martin-vymazal-vki) said :
#2

Hi Jeannette,

 thank you for your reply. I'm interested in simulating flow in (lower limb) arteries and in artery bypasses. Could you please also tell me if it's possible to set 'flux' boundary conditions? Imagine I can measure the total flow rate at the inlet of my geometry, but I don't know the velocity profile. Is there a way of prescribing only the flow rate?

 Thank you for the link to the dolfin converter. In fact, the mesh format is not a big issue - I could eventually write a script to convert Gmsh->Dolfin. I was more curious about what's the status of Unicorn at this moment (which you answered) and available boundary condition types.

 I remember reading somewhere that the FSI solver is monolithic. Can you confirm this?

 Best regards,

   Martin Vymazal

Revision history for this message
Jeannette Spühler (j-spuehler) said :
#3

Hi Martin,
sorry again for the delay.
Yes, the FSI solver is monolithic.
About the flow rate: There is no class taking care of this at the moment. Can't you express it with pressure condition at the inlet?

Best regards
Jeannette

Can you help with this problem?

Provide an answer of your own, or ask Martin Vymazal for more information if necessary.

To post a message you must log in.