CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > CFX

FSI Modeling

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

Like Tree1Likes
  • 1 Post By stumpy

Reply
 
LinkBack Thread Tools Display Modes
Old   June 22, 2011, 17:54
Exclamation FSI Modeling
  #1
New Member
 
morteza
Join Date: Jul 2009
Location: United States
Posts: 12
Rep Power: 8
kimiaghalam is on a distinguished road
Hi Everyone

I am trying to simulate the opening motion of a heart valve during systolic phase in an artificial heart valve design. I am using immersed solid model for valves. In this case the opening velocity of the valve should be a function of the torque exerted on the valve by the fluid. I was trying to calculate the torque in each time step using experession torque_z@Locator in CFX Pre and then calcute the angular velocity by another simple expression. But CFX-Solver stops at the beginning of the run without any specific Errors! I am pretty sure that the problem is related to using this torque expression. Do you guys have any idea where the problem comes from ?

Appreciate your help
Morteza
kimiaghalam is offline   Reply With Quote

Old   June 23, 2011, 09:05
Default
  #2
Senior Member
 
Join Date: Apr 2009
Posts: 516
Rep Power: 12
stumpy is on a distinguished road
Rather than calculating the motion yourself use the 6-DOF rigid body solver to calculate the motion, then apply this motion to the immersed solid.
stumpy is offline   Reply With Quote

Old   June 23, 2011, 17:48
Default Hi Stumpy
  #3
New Member
 
morteza
Join Date: Jul 2009
Location: United States
Posts: 12
Rep Power: 8
kimiaghalam is on a distinguished road
Quote:
Originally Posted by stumpy View Post
Rather than calculating the motion yourself use the 6-DOF rigid body solver to calculate the motion, then apply this motion to the immersed solid.
I am using Ansys CFX 12 and I could not find anything about 6-DOF rigid body solver. Could you please help me find some tutorials or some help on this solver?

Thanks for your consideration
kimiaghalam is offline   Reply With Quote

Old   June 23, 2011, 19:00
Default
  #4
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,929
Rep Power: 85
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
The 6DOF solver had has considerable updates and improvements going into V13. In V12 it was just a beta feature. Sounds like it is time to update to the current version.
ghorrocks is offline   Reply With Quote

Old   July 1, 2011, 19:37
Default
  #5
New Member
 
morteza
Join Date: Jul 2009
Location: United States
Posts: 12
Rep Power: 8
kimiaghalam is on a distinguished road
I have updated my Ansys to Ansys 13. Rigid Body solver seems to be a really helpful module. But now there is this Error occuring in ansys solver:

********* WARNING ********* |
| There are 118921 interpolation points ( 36.56%) that are not properly mapped onto a source domain. This will likely introduce incorrect terpolated values on the immersed boundary. If this is due to a slight difference in mesh, please increase the expert parameter "octree bndbox tolerance". Boundary Patch Name : Left_Leaflet Default |

Can you help me with this error. After this warning I run out of memory and I can not increase the memory factor any further than 1.3 and I can not find the "octree bndbox tolerance"

Appreciate your help
kimiaghalam is offline   Reply With Quote

Old   July 2, 2011, 00:19
Default
  #6
New Member
 
morteza
Join Date: Jul 2009
Location: United States
Posts: 12
Rep Power: 8
kimiaghalam is on a distinguished road
I can run the model in 2D mode but I want to know how to make the valves stop moving when they touch the wall ? Right now they just pass the wall !!!
kimiaghalam is offline   Reply With Quote

Old   July 4, 2011, 09:21
Default
  #7
Senior Member
 
Join Date: Apr 2009
Posts: 516
Rep Power: 12
stumpy is on a distinguished road
That's expected, there's no contact modelling in CFX. You can try applying an external force/torque when the contact criteria is met, but that's probably going to produce chatter. you could use an interrupt control to stop the run when contact is made, then restart the run after turning off the rigid body solver (or just set 'solve meshdisp = f' to turn of moving mesh).
kimiaghalam likes this.
stumpy is offline   Reply With Quote

Old   July 6, 2011, 16:50
Default
  #8
New Member
 
morteza
Join Date: Jul 2009
Location: United States
Posts: 12
Rep Power: 8
kimiaghalam is on a distinguished road
Is there any way I can turn of the rigid body solver when it reaches a certain angle !?
kimiaghalam is offline   Reply With Quote

Old   October 12, 2011, 03:16
Default
  #9
Member
 
Join Date: May 2010
Posts: 42
Rep Power: 7
Zymon is on a distinguished road
Quote:
Originally Posted by kimiaghalam View Post
Is there any way I can turn of the rigid body solver when it reaches a certain angle !?
I have the same question right now. Has anybody solved this problem yet? I am also dealing with artificial heart valves. Due to the large deformations the immersed solid approach is quite promising. But the motion of the valves should stop after about 90 degrees. Furthermore, I still have the problem that the flow passes through the rigid body (streamlines just change direction at the valve slightly). I already set the Momentum Source Scaling Factor to 400 (together with the command 'smooth inside ims = t'). Is there any other option to handle this?
Zymon is offline   Reply With Quote

Old   October 12, 2011, 16:52
Default
  #10
Senior Member
 
Join Date: Apr 2009
Posts: 516
Rep Power: 12
stumpy is on a distinguished road
You can stop the simulation when a certain angle is reached using the solver interrupt controls together with an expression like rbstate(Euler Angle X)@RigidBody. You can't turn off the rigid body solver without stopping, but I guess your boundary condition that uses the displacements calculated by the rigid body solver could contain some logic to ignore the displacement if rbstate(Euler Angle X)@RigidBody (for example) exceeded some value.
stumpy is offline   Reply With Quote

Old   October 14, 2011, 05:47
Default
  #11
Member
 
Join Date: May 2010
Posts: 42
Rep Power: 7
Zymon is on a distinguished road
Hi stumpy, thanks a lot! I think this is working.

The flow still passes through the solid. At least it seems so when I look at the streamlines. I guess one problem is that the valves are quite thin. I already use a very fine mesh of the flow field, but maybe this is not enough. Furthermore, the flow is partly relative slow. So the momentum scaling factor has little effect. Still it looks quite nice already. Most likely with the mesh deformation approach this is not possible yet due to the large deformations.
Zymon is offline   Reply With Quote

Old   October 14, 2011, 14:23
Default
  #12
Senior Member
 
Join Date: Apr 2009
Posts: 516
Rep Power: 12
stumpy is on a distinguished road
Yes, I've seen that problem with slow/stationary flows - if you scale zero by a large number it is still zero! You could try adding a fixed momentum sink via a subdomain and make use of an expression like:
if(inside()@ImmersedSolidDomain,1,0)*MomentumSrc
so that you are getting a larger momentum sink in the immersed solid.
stumpy is offline   Reply With Quote

Old   October 17, 2011, 11:03
Default
  #13
Member
 
Join Date: May 2010
Posts: 42
Rep Power: 7
Zymon is on a distinguished road
Can you please explain this in more detail? Where exactly should I add the expression?
Zymon is offline   Reply With Quote

Old   October 27, 2011, 04:25
Default
  #14
Member
 
Join Date: May 2010
Posts: 42
Rep Power: 7
Zymon is on a distinguished road
I have a question concerning the restart of the simulation. When I restart the simulation with switching off the rigid body solver (immersed solid to stationary) the flow is in its initial condition but the rigid bodies are in its very first position of the first simulation. Do I have to rotate them manually for the restart? Or can I somehow initialize them as well?
Zymon is offline   Reply With Quote

Old   May 22, 2012, 08:56
Default
  #15
New Member
 
belgacem
Join Date: Jan 2012
Posts: 22
Rep Power: 5
belgacem is on a distinguished road
Hi friend
First, I have the same problem right now. Please, can you explain more how to added a fixed momentum sink via subdimain and to make this expression if(inside()@ImmersedSolidDomain,1,0)*MomentumSrc
Also, I am studying immersed boundary method and and try to simulate a block falling in the water. I am using "immersed solid" then rigid body 6DOF and I let it fall freely but it can't stoped in the bottom where the velocity must be zero. I give a density to the block and i let it fall freely under gravity. Noted that the rigid body is defined as an immersed solid. i have specified a stationary coordinate frame that has its origin at the center of mass of the physical rigid body. Another fixed coordinate frame was specified related to the water at rest.
What can I do to stopped the rigid body in the bottom where the potentiel energy must be zero?

thank you!
belgacem is offline   Reply With Quote

Reply

Tags
force, fsi, hear valve, immersed solid, torque

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
FSI modeling of blood vessel??? smn CFX 18 April 27, 2015 23:50
FSI modeling of a surface effect ship!!! smn CFX 10 January 26, 2010 17:26
initial file for fsi modeling? mortazavi CFX 0 July 21, 2009 02:06
ccl for fsi modeling of blood vessel smn CFX 0 July 19, 2009 05:36
Urgent..Needs for FSI modeling using ANSYS Eugene CFX 7 March 27, 2005 11:19


All times are GMT -4. The time now is 12:19.