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

Blood Damage Modelling via Particle Tracking in a Centrifugal Heart Pump

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 4, 2011, 10:56
Default Blood Damage Modelling via Particle Tracking in a Centrifugal Heart Pump
  #1
New Member
 
Join Date: Mar 2011
Posts: 3
Rep Power: 15
scatman is on a distinguished road
Hi All!

First of all, I would like to say what an incredible help this forum has been to me already. The knowledge available and shared on this forum is outstanding, it's truly been an invaluable resource.

Now, down to the problem...

I'm investigating centrifugal heart blood pumps (VAD) as my final year engineering project using ANSYS CFX. My supervisor has already performed steady-state simulations (without particle tracking) of the pump I am considering, and the results of my simulations match his closely (also without particle tracking). I desire to model blood damage by introducing particle tracking to my simulations.

Reading the literature available, and other academic papers, one-way coupling (as opposed to fully coupled) is an acceptable implementation of particle tracking for this scenario, which allows me just to re-run my simulations with the results as the initial values (so that only the particle tracks are calculated).

The formula for quantifying red blood cell damage is related to shear stresses and the exposure times of these stresses on the cells, as presented below:

(dHb/Hb) = (1.8e-6)([tau]^1.991)([dT]^0.765)

If we consider one particle track from the inlet of the pump to the outlet, summing the above equation at each position of the particle from the inlet to the outlet will result in the total cell damage for that particle. Thus, the above equation needs to be calculated at each point of each particle track. Furthermore, each particle track needs to be considered separately in this fashion, since multiple particle tracks may cross the some point in space but have different histories and futures.

It would be very appreciated if someone could point me in the correct direction to solve my problem.

ps. I do not need help determining the shear stress, academic papers have suggested the approxiamation of "tau=2*Density*Turbulent Kinetic Energy" is acceptable.

Thanks!!!!

Peter

Last edited by scatman; May 4, 2011 at 11:32.
scatman is offline   Reply With Quote

Old   May 4, 2011, 15:58
Default
  #2
Senior Member
 
Join Date: Apr 2009
Posts: 531
Rep Power: 21
stumpy is on a distinguished road
I wouldn't use particle tracking at all for this. Instead, investigate using an Additional Variable using the transport equation option (i.e. it just follows the flow). Create a source term for the Additional Variable that is a function of the shear stress and timestep. The Additional Variable value represent the total cell damage. Have a search on this forum for using Additional Variables to track fluid residence time, it's the same idea.
stumpy is offline   Reply With Quote

Old   May 4, 2011, 19:17
Default
  #3
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,830
Rep Power: 144
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
I agree with Stumpy. Why do this using a Lagrangian particle tracking approach? The Eularian approach suggested by Stumpy will be much more computationally efficient. It will take a bit of setting up but it sounds like a good method.
ghorrocks is offline   Reply With Quote

Old   May 4, 2011, 23:54
Default
  #4
New Member
 
Join Date: Mar 2011
Posts: 3
Rep Power: 15
scatman is on a distinguished road
Thanks for the quick reply! I'll give your suggestion a shot and see how it goes.

Cheers!
scatman is offline   Reply With Quote

Old   May 5, 2011, 02:36
Default
  #5
New Member
 
Join Date: Mar 2011
Posts: 3
Rep Power: 15
scatman is on a distinguished road
Glenn and Stumpy, I had a meeting with my supervisor today and he insists that a Lagrangian approach be taken. I have a feeling he wants the simulations to have the same method as other simulations in the field of ventricular assist devices, all of which i've seen use Lagrangian tracking.

So, although I realize it's not the most elegant solution, any help with using a Lagrangian approach to this problem would be greatly appreciated.
scatman is offline   Reply With Quote

Old   May 5, 2011, 07:23
Default
  #6
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,830
Rep Power: 144
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
If I was you, until you find a good reason to do it using Lagrangian particle tracking I would do it using both approaches and you will see for yourself the advantages of the approach Stumpy has suggested. I don't find "because everybody else does it that way" a very convincing argument.

Particle tracking is demonstrated in several CFX tutorial examples. That is the best place to start.
ghorrocks is offline   Reply With Quote

Old   January 7, 2018, 21:53
Default
  #7
New Member
 
Someeee
Join Date: Dec 2017
Posts: 4
Rep Power: 8
Someeee is on a distinguished road
Hello,
I am also trying to do blood damage modelling. The blood damage index is function of shear stress and exposure time, as mentioned in this post previously. Can you please give details on how to do this using Eulerian approach without injecting particles? I am uncertain about creating a source term as mentioned in this post.

Any help will be greatly appreciated.
Someeee is offline   Reply With Quote

Old   January 8, 2018, 00:59
Default
  #8
Member
 
usama
Join Date: Jul 2016
Location: Pakistan
Posts: 33
Rep Power: 10
usamaperwez is on a distinguished road
First of all make additional variable:
Variable type: Volumetric
Units: as per your problem
Tensor Type: scalar

Then Add expression:
For example like this "(0.26[s^-1]*(1+tanh(-(T*1[K^-1]-335.3)/5)))/2)"
It need to stable in terms of unit

Then add subdomain:
Tick Sources and add the expression name in that.

Here we go.

I basically did it for DNA Species and it worked excellently (it took me 3 months to model that).
usamaperwez is offline   Reply With Quote

Reply

Tags
centrifugal pump, particle tracking, vad

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
massless particle tracking problem Renold FLUENT 0 January 26, 2011 14:23
Particle Tracking for ion Jun CFX 2 August 31, 2010 08:19
Number density tracking rather than particle tracking Rebecca Main CFD Forum 2 April 23, 2009 12:52
DPM UDF particle position using the macro P_POS(p)[i] dm2747 FLUENT 0 April 17, 2009 01:29
modelling of gap leakage flow in centrifugal pump melvin CFX 3 January 14, 2008 16:40


All times are GMT -4. The time now is 01:57.