|
[Sponsors] |
Blood Damage Modelling via Particle Tracking in a Centrifugal Heart Pump |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
May 4, 2011, 10:56 |
Blood Damage Modelling via Particle Tracking in a Centrifugal Heart Pump
|
#1 |
New Member
Join Date: Mar 2011
Posts: 3
Rep Power: 15 |
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. |
|
May 4, 2011, 15:58 |
|
#2 |
Senior Member
Join Date: Apr 2009
Posts: 531
Rep Power: 21 |
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.
|
|
May 4, 2011, 19:17 |
|
#3 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,830
Rep Power: 144 |
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.
|
|
May 4, 2011, 23:54 |
|
#4 |
New Member
Join Date: Mar 2011
Posts: 3
Rep Power: 15 |
Thanks for the quick reply! I'll give your suggestion a shot and see how it goes.
Cheers! |
|
May 5, 2011, 02:36 |
|
#5 |
New Member
Join Date: Mar 2011
Posts: 3
Rep Power: 15 |
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. |
|
May 5, 2011, 07:23 |
|
#6 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,830
Rep Power: 144 |
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. |
|
January 7, 2018, 21:53 |
|
#7 |
New Member
Someeee
Join Date: Dec 2017
Posts: 4
Rep Power: 8 |
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. |
|
January 8, 2018, 00:59 |
|
#8 |
Member
usama
Join Date: Jul 2016
Location: Pakistan
Posts: 33
Rep Power: 10 |
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). |
|
Tags |
centrifugal pump, particle tracking, vad |
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 |