CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Journal

Coe et al. - 2023 - IKA-FLOW : A Flexible Body Overset Mesh Impl. for Fish Swimming

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 6, 2023, 02:20
Default Coe et al. - 2023 - IKA-FLOW : A Flexible Body Overset Mesh Impl. for Fish Swimming
  #1
New Member
 
Jozsef Nagy
Join Date: Dec 2018
Posts: 14
Rep Power: 7
jnagy is on a distinguished road
Comments and questions are welcomed on the article Coe et al., IKA-FLOW : A Flexible Body Overset Mesh Implementation for Fish Swimming, OpenFOAM Journal, https://doi.org/10.51560/ofj.v3.89.
jnagy is offline   Reply With Quote

Old   July 7, 2023, 05:08
Default
  #2
mco
New Member
 
Michael Coe
Join Date: Jan 2021
Posts: 7
Rep Power: 5
mco is on a distinguished road
Please let me know if you have any questions or comments, or need any clarification about the code.
mco is offline   Reply With Quote

Old   August 10, 2023, 14:14
Default
  #3
New Member
 
Rithwik Kandukuri
Join Date: Mar 2023
Posts: 4
Rep Power: 3
rkandukuri is on a distinguished road
in figures.31-34 , does the background mesh only part refer to a body fitted mesh or the overset assembly with the component meshes blanked out ? If it's the latter, whats the rationale behind including that ? I'm sorry If I missed that bit in the article ?
rkandukuri is offline   Reply With Quote

Old   August 10, 2023, 17:12
Default
  #4
mco
New Member
 
Michael Coe
Join Date: Jan 2021
Posts: 7
Rep Power: 5
mco is on a distinguished road
Hi rkandukuri,

The background mesh only refers to the background mesh with the overset component mesh not rendered (using a threshold). This is to show the hole and the corresponding flow field for just the background mesh.

The rationale is that we found it interesting comparing the flow fields of the two and thought a reader might also be interested.

A comparison with a body fitted mesh would be good as well, but it was quite difficult to get the diffusion functions to work with the inflation layer requirements we had (35 layers). Without going into too much detail, we were looking at how fish propulsion scales from Reynolds 1e3 to 1e7, which is why the high number of inflation layers was needed. This is what prompted the development of the overset meshing route.

Please let me know if you have any other questions.

Michael
mco is offline   Reply With Quote

Old   August 26, 2023, 11:57
Default
  #5
New Member
 
Manchester
Join Date: Aug 2022
Posts: 27
Rep Power: 3
Tianyz is on a distinguished road
Quote:
Originally Posted by mco View Post
Please let me know if you have any questions or comments, or need any clarification about the code.
Hi Michael,

Your work is amazing, I'm using overset mesh based on OpenFOAM. But it's only for a rigid body, your contribution fills this empty. I still reading your code and want to talk with you further. Since your solver is focusing on prescribed motion (one way) currently, I'm wondering if you plan to extend it to two-way coupling in the future?

Best,
Tian
Tianyz is offline   Reply With Quote

Old   August 27, 2023, 18:18
Default
  #6
mco
New Member
 
Michael Coe
Join Date: Jan 2021
Posts: 7
Rep Power: 5
mco is on a distinguished road
Hi Tianyz,

Thank you for the comment. I did have a plan to implement the forces class inside the solver to get the average forces and then change the amplitude variable to balance the thrust and drag. This is outlined in the great work done by Audrey Maertens in Maertens 2017.

Currently, I'm working on another project developing heat transfer solvers. If you would like to implement something as described above, send me a private message.

Hope that helps,

Michael
mco is offline   Reply With Quote

Old   October 11, 2023, 02:10
Default
  #7
New Member
 
SUN WEIYI
Join Date: Jul 2022
Posts: 10
Rep Power: 3
sunweiyi2020 is on a distinguished road
Hi jnagy,

About the overset mesh, I have one question. Can it be used to deal with the deformation of elastic structure connected with the bottom boundary? I try to use overset mesh on this situation. But the simulation diverged quickly.
sunweiyi2020 is offline   Reply With Quote

Old   October 11, 2023, 03:25
Default
  #8
mco
New Member
 
Michael Coe
Join Date: Jan 2021
Posts: 7
Rep Power: 5
mco is on a distinguished road
Hi sunweiyi2020,


With an overset mesh, you need to make sure there are enough cells between the hole, overset mesh, and the boundary of your background mesh. I believe a minimum of 3 is the required. So you would want at least 3 overset mesh cells between your hole and any boundary in you background mesh.


I'm not sure what your case is exactly, but it sounds like the hole mesh would be right at the boundary. This is not going to really work. If you can give a screenshot of your mesh, it would be easier to give a more exact answer.
mco is offline   Reply With Quote

Old   October 12, 2023, 00:32
Default
  #9
New Member
 
SUN WEIYI
Join Date: Jul 2022
Posts: 10
Rep Power: 3
sunweiyi2020 is on a distinguished road
Hi mco,

This is my mesh. Please check it. The arc is the structure, a thin membrane. It is connected with bottom boundary. I try it with overset mesh but the result is so bad. [IMG]file:///C:/Users/Coastal%20Eng%20Lab/Desktop/mesh.png[/IMG]
If this link is not available, please get the mesh figure in this link https://github.com/NUGUY2022/mesh.git.
sunweiyi2020 is offline   Reply With Quote

Old   October 12, 2023, 21:21
Default
  #10
mco
New Member
 
Michael Coe
Join Date: Jan 2021
Posts: 7
Rep Power: 5
mco is on a distinguished road
In this case, it might be easier to just use a regular morphing mesh. If you look under the src/fvMotionSolver/pointPatchFields/derived folder in OpenFOAM, you will see motion solvers for patches. You can convert one of these to suite your needs.


Another option would be to code it all in the dynamicMeshDict using the new coding feature. An example of this is in the /tutorials/mesh/moveDynamicMesh/twistingColumn folder.


Are you calculating forces on the patch and moving it as an FSI simulation or are you prescribing the motion?
mco is offline   Reply With Quote

Old   October 14, 2023, 02:26
Default
  #11
New Member
 
SUN WEIYI
Join Date: Jul 2022
Posts: 10
Rep Power: 3
sunweiyi2020 is on a distinguished road
Sorry for late, I just check the forum.
Yes, you are right. I am calculating forces on the patch and moving it as an FSI simulation. But large mesh deformation occurred on my case, resulting in divergency. I am not familiar with these motion solvers in src/fvMotionSolver/pointPatchFields/derived. Could you give me any advice for selecting the most suitable solvers? Regarding others' suggestion, they recommend RBFMeshMotionSolver to me. But it seems that wave2foam doesn't support it. Do you have experience to integrate it into interFoam or wave2foam?
sunweiyi2020 is offline   Reply With Quote

Old   October 15, 2023, 16:42
Default
  #12
mco
New Member
 
Michael Coe
Join Date: Jan 2021
Posts: 7
Rep Power: 5
mco is on a distinguished road
Sorry, I'm not familiar with those two solvers. No reason not to try them out. You can play around with the diffusion equations in the dynamic mesh dict and lowering your time step.



What is happening with the mesh when it diverges? Is it a mesh problem or a numerics problem? If the forces are diverging, where on the patch is this happening?


Michael
mco is offline   Reply With Quote

Old   October 15, 2023, 23:05
Default
  #13
New Member
 
SUN WEIYI
Join Date: Jul 2022
Posts: 10
Rep Power: 3
sunweiyi2020 is on a distinguished road
I have uploaded the crashed mesh on the link https://github.com/NUGUY2022/mesh.git.
It is obvious that the mesh diverges rather than a numerical problem, because the case with small deformation can run without any problems. Now, I use the displaceLaplacian solver for mesh motion. I will try different mesh motion solvers for this case. In addition, do you think it can be solved by optimizing mesh?
sunweiyi2020 is offline   Reply With Quote

Old   November 1, 2023, 13:47
Default
  #14
New Member
 
Manchester
Join Date: Aug 2022
Posts: 27
Rep Power: 3
Tianyz is on a distinguished road
Quote:
Originally Posted by mco View Post
Please let me know if you have any questions or comments, or need any clarification about the code.
Hi, Coe

Sorry to reply to you late. I was tightening some other things. And I got one question for you. Form your paper, you said that "Besides fish motion, this library can be used to deform an overset mesh to accommodate any motion desired", what if I want to simulate a FSI case with a deformable structure with other motionsolver, say displacementLaplacian or 6DOF Motion Solver? should I recompile it? Thank you for any hints.

Best,
Tian
Tianyz is offline   Reply With Quote

Old   April 24, 2024, 03:10
Default
  #15
New Member
 
Forrest Abbott
Join Date: Apr 2024
Posts: 1
Rep Power: 0
Forrest Abbott is on a distinguished road
Another approach is to use mesh adaptation techniques, such as mesh refinement or coarsening, to dynamically adjust the resolution of the mesh near the deforming boundaries.
__________________
Forrest Abbott
Forrest Abbott is offline   Reply With Quote

Reply
Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
sliding mesh problem in CFX Saima CFX 46 September 11, 2021 07:38
Issues on the simulation of high-speed compressible flow within turbomachinery dowlee OpenFOAM Running, Solving & CFD 11 August 6, 2021 06:40
[ANSYS Meshing] Not able to mesh flow region only surface of body FredRaw ANSYS Meshing & Geometry 2 April 19, 2013 05:06
[ANSYS Meshing] Simple Symmetry Mesh Question - Flow around Bluff Body Matlab69 ANSYS Meshing & Geometry 20 April 23, 2012 09:55
fluid flow fundas ram Main CFD Forum 5 June 17, 2000 21:31


All times are GMT -4. The time now is 22:34.