How to simulate flow around multiple bodies?

243 views
Skip to first unread message

Li-Ming Chao

unread,
Nov 7, 2020, 2:25:59 AM11/7/20
to IBAMR Users
Dear all,

I want to simulate flow around multiple bodies using ConstraintIB, and I hope to find some existed examples to have a good beginning. Are there any examples conducted on the similar topic in IBAMR/examples folder?

With best wishes!

L-Ming

Amneet Bhalla

unread,
Nov 18, 2020, 6:15:07 PM11/18/20
to ibamr...@googlegroups.com, darv...@gmail.com, Nishant Nangia
Nishant/Thomas can one of you guys attach the multiple bodies using ConstraintIBMethod example here for Li-Ming. 

--
You received this message because you are subscribed to the Google Groups "IBAMR Users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to ibamr-users...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/ibamr-users/d99be87d-6064-49fd-a6c3-27c9a6365931n%40googlegroups.com.


--
--Amneet 



Nishant Nangia

unread,
Nov 18, 2020, 8:28:24 PM11/18/20
to ibamr...@googlegroups.com, darv...@gmail.com
Hi Li-Ming,

Please see attached source files for an example of how to simulate multiple bodies in ConstraintIBMethod.

Nishant
SinglePhase_DKT.zip

Li-Ming Chao

unread,
Nov 18, 2020, 8:33:43 PM11/18/20
to ibamr...@googlegroups.com, darv...@gmail.com
Dear Bhalla, Nangia,

Thanks so much for your help and examples!

With best regards!

Li-Ming

Nishant Nangia <nishantn...@gmail.com> 于2020年11月19日周四 上午9:28写道:

Li-Ming Chao

unread,
Nov 21, 2020, 2:40:37 AM11/21/20
to IBAMR Users
Dear  Nangia,  

Thanks so much for your multiple bodies example!

Following your example,  I wrote the code to simulate two oscillating foils. However, I found that the parameters of the second structure  have not been used during the simulation. I have upload the document, would you mind helping me to find the error?

With best wishes!

Li-Ming
test.zip
Message has been deleted

Haotian Hang

unread,
Nov 22, 2020, 10:13:59 PM11/22/20
to IBAMR Users
Hello all,
    I have a further question on the similar problem. If the structures are connected with some constrains, such as they are hinged, could you please give me some advise on how to implement this? Thank you so much!
Yours sincerely,
HH

Boyce Griffith

unread,
Dec 14, 2020, 1:29:37 AM12/14/20
to IBAMR Users
Which IB formulation are you using for handling rigidity constraints?

On Nov 22, 2020, at 1:20 AM, Haotian Hang <haot...@usc.edu> wrote:

Hello all,
    I have a further question on the same issue. If I want to simulate two structures which are connected, such as they are hinged, can you give me an advise on doing this? Thank you so much!

On Friday, November 20, 2020 at 11:40:37 PM UTC-8 clm0...@gmail.com wrote:

Haotian Hang

unread,
Dec 14, 2020, 1:35:33 PM12/14/20
to ibamr...@googlegroups.com
I'm trying to use ConstraintIB to solve this problem, I don't know whether it is appropriate. Thank you!

You received this message because you are subscribed to a topic in the Google Groups "IBAMR Users" group.
To unsubscribe from this topic, visit https://groups.google.com/d/topic/ibamr-users/6Zs0y7zS6nc/unsubscribe.
To unsubscribe from this group and all its topics, send an email to ibamr-users...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/ibamr-users/6EE60AD1-FA3F-4EE0-8359-BD08DDC13814%40gmail.com.


--
Haotian Hang

Amneet Bhalla

unread,
Mar 28, 2021, 11:45:00 AM3/28/21
to ibamr...@googlegroups.com, yangd...@gmail.com
Hi Dewu,

Please follow this example to include multiple bodies in the code. In the driver code, you would need to create a vector of kinematics object.

Thanks,
—Amneet

--
--Amneet 



SinglePhase_DKT.zip

Dewu Yang

unread,
Mar 29, 2021, 1:38:20 AM3/29/21
to IBAMR Users
Dear professor Amneet,

     Thanks for your reply. I followed the example which you gave me and changed the code of eel2d, but the error shown in the figure appeared when it was running. And I still have a question, shouldn't the eel2d.vertex file be a file output by executing code? Why do I give eel2d_0.vertex and eel2d_1.vertex in the code, and it shows that the eel2d_1.vertex file cannot be found when I called 'make examples'. I added the eel2d_1.vertex file by copying the eel2d_0.vertex file and renaming it,and manually changed the value. Why the eel2d_1.vertex file not automatically produced? I would appreciate it if you can give me some advice about it.

Best,
—Dewu Yang





sendpix0.jpg

Yang Zhang

unread,
May 4, 2022, 9:38:41 AM5/4/22
to IBAMR Users
Dear Professor,

I have the similar issue. In IBFE,  can we implement two structures which are connected at some nodes?

Best,
Yang

Boyce Griffith

unread,
May 4, 2022, 12:45:37 PM5/4/22
to noreply-spamdigest via IBAMR Users
On May 4, 2022, at 9:38 AM, Yang Zhang <yang.zha...@gmail.com> wrote:

Dear Professor,

I have the similar issue. In IBFE,  can we implement two structures which are connected at some nodes?

With IBFEMethod, if you use nodal coupling (as described here: https://arxiv.org/abs/2111.09958), then if two structures have overlapping nodes, those nodes will always move together. Will that work for your use case?

It also is possible to use both IBFEMethod and IBMethod in the same simulation (using class IBStrategySet). That would allow you to add springs (using IBMethod) to “tie” the IBFEMethod-based structures together. Again, you will want to use nodal coupling for the IBFEMethod-based model.

To setup IBFEMethod to use nodal coupling, you would include the following settings in the IBFEMethod configuration:

IBFEMethod {
   IB_use_nodal_quadrature = TRUE
   use_consistent_mass_matrix = FALSE
   // … other settings …
}

If you don’t use nodal interaction, then overlapping parts of the IBFE structures will tend to move together, but it is possible for them to detach from each other.

Amneet Bhalla

unread,
May 4, 2022, 1:09:32 PM5/4/22
to ibamr...@googlegroups.com
The eel vertex file is generated using the MATLAB script eel2d_straightswimmer.m and not through the cpp code. 

--
You received this message because you are subscribed to the Google Groups "IBAMR Users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to ibamr-users...@googlegroups.com.


--
--Amneet 



Yang Zhang

unread,
May 4, 2022, 10:13:55 PM5/4/22
to IBAMR Users
Prof. Griffith,

Thank you for your reply. I will follow your instructions and have a try.

Best,
Yang
Reply all
Reply to author
Forward
0 new messages