
[Sponsors] 
August 5, 2011, 02:03 
Womersley pulsatile boundary condition for blood flow.

#1 
Member
Guifan Li
Join Date: Apr 2011
Location: New York City, U.S.
Posts: 96
Rep Power: 8 
Dear all,
I know this topic is discussed before, but I didn't find solutions. So please bear with me. I try to use groovyBC to implement this function but: The womersley formula include a Bessel function which need to be passed into a imaginary complex number. To my knowledge, Bessel function under OpenFOAM J0(s) can only solve the real number. And I don't know how to transfer womersley function to real domain. I write a matlab code to simulate the womersley pulsatile flow that I want. (besselj can resolve imaginary complex number). But that's it, no idea how to corporate my matlab data with OpenFOAM since the data is dependent on both t (time) and y (position). I have beening struggling with this for a while and quite frustrated... Can anyone give some help or hint please? 

August 5, 2011, 03:47 

#2 
Senior Member

Hello,
have you tried to play around with pythonflu? http://pythonflu.wikidot.com/ maybe it gives you some opportunities that might be of value to you. 

August 5, 2011, 06:02 

#3  
Member
Guifan Li
Join Date: Apr 2011
Location: New York City, U.S.
Posts: 96
Rep Power: 8 
Quote:
Thanks for your reply, I went to pythonflu, is seems like can integrate with Matlab results. But I didn't find any manual of pythonflu. Can you please suggest me in some details how to solve womersley inlet in pythonflu. Do I need to import data from matlab or I can just write some function in pythonflu? Cheers! 

August 5, 2011, 07:06 

#4  
Assistant Moderator
Bernhard Gschaider
Join Date: Mar 2009
Posts: 4,017
Rep Power: 43 
Quote:


August 5, 2011, 07:32 

#5 
Member
Guifan Li
Join Date: Apr 2011
Location: New York City, U.S.
Posts: 96
Rep Power: 8 
I forgot to say is that the Womersley number is known in my experiment.
Last edited by liguifan; August 5, 2011 at 09:56. 

August 5, 2011, 09:55 

#6  
Member
Guifan Li
Join Date: Apr 2011
Location: New York City, U.S.
Posts: 96
Rep Power: 8 
Quote:
Do you think I can use the following formula which which is in real domain to implement the Womersley solution so that I can pass into groovyBC. Many thanks 

August 8, 2011, 05:40 

#7 
Assistant Moderator
Bernhard Gschaider
Join Date: Mar 2009
Posts: 4,017
Rep Power: 43 
If the computer algebra program of your choice (or you calculating by hand) gives back a reasonable expression for the closed part (like "Re[u[y,t]]" in Mathematica) there should be no problem. I only glanced at your formula and I don't quite get what sould make it complex anyway (if t and y are real)


August 10, 2011, 06:59 

#8  
Member
Guifan Li
Join Date: Apr 2011
Location: New York City, U.S.
Posts: 96
Rep Power: 8 
Quote:
"FullSimplify[Re[u[y,t]]]" , the return value is still in the form of "Re[...]". I will put more time on this see whether can I get a nice derivation. Thanks again for your help! 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Wind turbine simulation  Saturn  CFX  45  February 8, 2016 05:42 
Need Help on boundary condition in MultiPhase flow  creddy_trddc  CFX  0  December 9, 2010 08:53 
Swirling flow as a boundary condition  hami9293  FLUENT  0  November 2, 2010 08:09 
mass flow in is not equal to mass flow out  saii  CFX  2  September 18, 2009 08:07 
Convective Heat Transfer  Heat Exchanger  Mark  CFX  6  November 15, 2004 16:55 