Womersley pulsatile boundary condition for blood flow.
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? 
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. 
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! 
Quote:

1 Attachment(s)
I forgot to say is that the Womersley number is known in my experiment.

1 Attachment(s)
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 
Quote:

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! 
All times are GMT 4. The time now is 23:32. 