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

Blood flow with FSI in Ansys 12

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 15, 2010, 10:28
Default Blood flow with FSI in Ansys 12
  #1
New Member
 
Marten
Join Date: Jun 2010
Posts: 6
Rep Power: 7
Marten is on a distinguished road
Hi everyone,

for my diplom thesis I'm trying to set up a pulsatile blood flow in a vessel including fluid structur interaction. For this I use CFX and ansys under workbench enviroment. I created the vessel in the design modeler and also used the meshing tools included in workbench. The set up is almost the same like in the cfx tutorial with the oscillating plate. But when I come to run the solver I get an error after some timesteps, that the mesh gets folded. I think I have to do a remeshing everytime the mesh quality gets too low. I'm not able to do this with the documentation provided with cfx. can anybody help me with this? I also tried to set the stiffness of the mesh to 1, which seems to be the lowest, but I would rather like to solve the problem with remeshing.
Do I have to remesh the structure and the fluid domain seperatly, 'cause this is how I generated the meshes in the first place?
I'm confused and need help.


kind regards
Marten

Last edited by Marten; June 15, 2010 at 10:49.
Marten is offline   Reply With Quote

Old   July 25, 2010, 08:57
Default mesh type
  #2
New Member
 
Mehran Saeedi
Join Date: Jan 2010
Posts: 8
Rep Power: 7
Mehran is on a distinguished road
Hi,

Can I ask you what is the type of mesh you are using?
Mehran is offline   Reply With Quote

Old   July 27, 2010, 05:25
Default
  #3
New Member
 
Join Date: Apr 2009
Posts: 12
Rep Power: 8
spatialtime is on a distinguished road
You can't do an automatic remesh when using FSI. CFX won't allow you to add a "configuration" sub module.
spatialtime is offline   Reply With Quote

Old   March 23, 2012, 13:58
Default help to set the inflow BC in CFX
  #4
Member
 
Hesam Moghaddam
Join Date: Mar 2012
Posts: 49
Rep Power: 5
hesamking is on a distinguished road
i am doing a research on an aneurysm, I want to redo an
already done work to learn details of this work.
I don't know how I should set up the inlet velocity. I have the inlet velocity waveform graph from the paper, but I do not have any expression or equation which i can input into CFX. will you please guide me what I should do?
I am doing my study for an aneurysm at a bifurcation of a basilar artery. I have one inlet and two outlets? I also wanted to ask what B.C i should use for the outlet. In paper, I think it mentions Newmann B.C.
I have attached the picture of inflow waveform. I would really appreciate your kind help.
Regards,
Heasam
Attached Images
File Type: jpg untitled.JPG (40.2 KB, 27 views)
hesamking is offline   Reply With Quote

Old   March 23, 2012, 15:13
Default
  #5
New Member
 
Marten
Join Date: Jun 2010
Posts: 6
Rep Power: 7
Marten is on a distinguished road
you can use a "user function" (for example called 'veloprofile') and make a series of time points connecting to the specific velocity. this user function can be used in an expression .in the example type 'veloprofile(time)' . this makes your user function dependent of the simulation time. now use your expression as inlet velocity by selecting inlet velocity normal to boundary or something like that. I would suggest that you use an opening boundary condition at your to outlets and set the difference pressure to 0 Pa. Neumann B.C. simply means that you specify a value for the BC and not a gradient.

many greetings
Marten
Marten is offline   Reply With Quote

Old   March 23, 2012, 15:31
Default
  #6
Member
 
Hesam Moghaddam
Join Date: Mar 2012
Posts: 49
Rep Power: 5
hesamking is on a distinguished road
Thanks Marten,
The problem is finding the function to input as UDF. how should i do that? i have already extracted the data points out of the curve I showed in the picture. I have it as CVS file. Any Idea for getting the function?
Best,
Hesam
hesamking is offline   Reply With Quote

Old   March 23, 2012, 18:19
Default
  #7
Senior Member
 
Rikio
Join Date: Mar 2009
Location: SH, China
Posts: 182
Blog Entries: 1
Rep Power: 8
rikio is on a distinguished road
Send a message via Skype™ to rikio
Hesam, UDF works in Fluent, and CEL in CFX.

You can get this profile function from Excel, for a simply expression, or obtain it from a profesional mathmatical software, like Matlab.
rikio is offline   Reply With Quote

Old   March 24, 2012, 07:52
Default
  #8
New Member
 
Marten
Join Date: Jun 2010
Posts: 6
Rep Power: 7
Marten is on a distinguished road
In CFX you can simply use the data points and it will interpolate between these. You can import your data points for example as a .txt-file (each data pair in a line devided by komma, I think) by clicking in the user function area under the unit selection and select 'import data'. in the unit field type for example [s] and below [m s^-1]. This is a very practical tool in CFX. If you want to solve it as a mathematical expression you might need to split the function into several parts (and produce polynomes in Excel for example) and use the step function in CEL. If you use Fluent you need to write a UDF file, which is a bit more complex.

greetings
Marten
Marten is offline   Reply With Quote

Old   March 24, 2012, 08:02
Default
  #9
New Member
 
Marten
Join Date: Jun 2010
Posts: 6
Rep Power: 7
Marten is on a distinguished road
This is an example of FLUENT UDF for an unsteady flow with sinus-form

/************************************************** ********************
unsteady.c
UDF for specifying a transient velocity profile boundary condition
************************************************** *********************/



#include "udf.h"

DEFINE_PROFILE(unsteady_velocity, thread, position)
{
face_t f;
real t = CURRENT_TIME;
real coeff;

coeff = 2.*3.14159265357/(60./60.);
begin_f_loop(f, thread)
{
F_PROFILE(f, thread, position) = 0.6 + 0.4*sin(t*coeff);
}
end_f_loop(f, thread)
}
Marten is offline   Reply With Quote

Old   April 14, 2012, 16:33
Default
  #10
Member
 
Hesam Moghaddam
Join Date: Mar 2012
Posts: 49
Rep Power: 5
hesamking is on a distinguished road
Dear Marten thank you so much.
Actually as I figured out the best way is to define a periodic function for
your inlet velocity cause the user profile data may not do the interpolation good and end up with some wrong results as the input. Anyhow, I created a simple expression like 1+sin(t) and then for my inlet velocity I took cylinderical coordinated with just axial velocity components and the other 2 one equal to zero.
So for my axial component I took this:
u= (1+sin(t)*(1-r^2)
and my steady state run converges good but my transient one is fluctuating
and after the run is complete I can see any results, for example for pressure on the artery wall.
I hope I could explain it good.
Thanks
hesamking is offline   Reply With Quote

Reply

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
Pulsatile blood flow Duncan FLUENT 11 April 10, 2013 12:17
fsi model of water flow in a pipe alinik CFX 11 November 5, 2009 12:29
ccl for fsi modeling of blood vessel smn CFX 0 July 19, 2009 05:36
Modeling blood flow - FloWorks mcneelyd FloEFD, FloWorks & FloTHERM 2 June 15, 2009 12:53
modelling blood flow in arteries with starcd sara CD-adapco 5 April 10, 2007 09:17


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