# FSI Modeling

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 June 22, 2011, 17:54 FSI Modeling #1 New Member   morteza Join Date: Jul 2009 Location: United States Posts: 12 Rep Power: 10 Sponsored Links 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

 June 23, 2011, 09:05 #2 Senior Member   Join Date: Apr 2009 Posts: 532 Rep Power: 14 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.

June 23, 2011, 17:48
Hi Stumpy
#3
New Member

morteza
Join Date: Jul 2009
Location: United States
Posts: 12
Rep Power: 10
Quote:
 Originally Posted by stumpy 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

 June 23, 2011, 19:00 #4 Super Moderator   Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 13,725 Rep Power: 106 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.

 July 1, 2011, 19:37 #5 New Member   morteza Join Date: Jul 2009 Location: United States Posts: 12 Rep Power: 10 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

 July 2, 2011, 00:19 #6 New Member   morteza Join Date: Jul 2009 Location: United States Posts: 12 Rep Power: 10 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 !!!

 July 4, 2011, 09:21 #7 Senior Member   Join Date: Apr 2009 Posts: 532 Rep Power: 14 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.

 July 6, 2011, 16:50 #8 New Member   morteza Join Date: Jul 2009 Location: United States Posts: 12 Rep Power: 10 Is there any way I can turn of the rigid body solver when it reaches a certain angle !?

October 12, 2011, 03:16
#9
Member

Join Date: May 2010
Posts: 42
Rep Power: 9
Quote:
 Originally Posted by kimiaghalam 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?

 October 12, 2011, 16:52 #10 Senior Member   Join Date: Apr 2009 Posts: 532 Rep Power: 14 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.

 October 14, 2011, 05:47 #11 Member   Join Date: May 2010 Posts: 42 Rep Power: 9 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.

 October 14, 2011, 14:23 #12 Senior Member   Join Date: Apr 2009 Posts: 532 Rep Power: 14 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.

 October 17, 2011, 11:03 #13 Member   Join Date: May 2010 Posts: 42 Rep Power: 9 Can you please explain this in more detail? Where exactly should I add the expression?

 October 27, 2011, 04:25 #14 Member   Join Date: May 2010 Posts: 42 Rep Power: 9 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?

 May 22, 2012, 08:56 #15 New Member   belgacem Join Date: Jan 2012 Posts: 22 Rep Power: 7 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!

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

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post smn CFX 18 April 27, 2015 23:50 smn CFX 10 January 26, 2010 17:26 mortazavi CFX 0 July 21, 2009 02:06 smn CFX 0 July 19, 2009 05:36 Eugene CFX 7 March 27, 2005 11:19