Life Sciences Applications: Modeling and Simulation for Biomedical - - PowerPoint PPT Presentation
Life Sciences Applications: Modeling and Simulation for Biomedical Device Design Kristian.Debus@cd-adapco.com SGC 2013 Mode deling and S and Simu mulatio ion f for Biomed medical De Device Desig e Design Biomedical device design
Kristian.Debus@cd-adapco.com
– Stents – Pumps – Heartvalves – Artificial Organs – Catheters – Pacemakers – Respiratory Aids
– Lab on a Chip – Implanted sensors – Implanted drug delivery
– MRI/CT Scanners – Ultrasound
– Lung/Heart Machine – Dialysis
– Blood Pressure – ECG, EEG, dissolved gases
– Lasers – Infusion Pumps
for
ing a V&V c V&V commi mmittee t that hat i is appl plication-specif cific t ic to t
medical dical device ice ind indus ustry Some E Example Cas ases: :
– FDA CPI I – Nozzle – Hemolysis Modeling – Drug delivery to the eye, by intravitreal injection – Oscillatory Pipe Flow – Flow in a Flexible Pipe – CFD Challenge – Aneurysms – Porous media modeling (Fiber bundles) – Oxygenator – Particle tracking etc. etc…..
Mouth Cavity
Simpleware was used to obtain the complex geometry from MRI of the human body
8
Ventilation flow and convective cooling as required for MRI/CT scanners, ICU devices
Surface Wrapping, STL Cleanup & Polyhedral Meshing Rapid Turnaround of Complex Geometry
Dissected Aorta Polyhedral Mesh, Geometry Provided by the Methodist DeBakey Heart and Vascular Center, Houston (Dr. Christof Karmonik, Dr. Mark Davies, Dr. Alan Lumsden, Dr. Jean Bismuth)
AAA (Abdominal Aortic Aneurysm) Geometry Provided by Computational Clinical Modeling, New Jersey (Chris Ebeling)
Car ardio diovas ascul cular ar flow w wave form from
applie ied at at the in inlet Materi rial pr properties of bl f blood ( (Newtonian A Appr pproxima mation)
» Density = 1056 kg m-3 » Dynamic Viscosity = 0.0035 Pas
Win indkes essel p par aram ameters t to
define ine t the ou
condit itio ion: n:
» Z = 1.1x107 [kg m-4 s-1] » R = 1.45x108 [kg m-4 s-1] » C = 1.45x10-8 [m4 s2 kg-1]
Lami minar r fl flow mo model
Sim imul ulat atio ion w was r run un for a a nu number of
cycl cles t to e
per erio iod r respons
was ach achie ieved.
Z
2
R
2
C
2
Z
3
R
3 C 3
Z
1
R
1
C
1
Z
4
R
4
C
4
Anal nalyt ytical ical s sol
utio ion [ n [2]
Num umerical ical s solut utio ion n
Analytical Solution
[2] Brown A. G., Patient-Specific Local and Systemic Haemodynamics in the Presence of a Left Ventricular Assist Device, 2012. PhD Thesis, The University of Sheffield.
speed of near 50 m/s when bulk modulus
travel at sound speed of the fluid (1500 m/s).
into radial strain energy in the solid when it travels there is nothing left to push the pulse down the pipe.
expect wave speed, the wave travels one cell down the axis. So the smaller the modulus, the smaller the wave speed, the larger the time step yet still accurate and stable!!
Damon Afkari, Universidad Politécnica de Madrid
Universidad Politécnica de Madrid, Damon Afkari: PhD Student Developed Proprietary Explicit Coupling Methodology
Abaqus
Surgeon Decision Making Aortic Dis ic Dissect ctio ion
CA CAD
–
0.43 mm base size
–
Separate motion for each leaflet
–
Interpolated and mapped
–
“Minimum points in a gap” of 4 nodes used to control proximity
–
Mesh varied from ~700k to 1.8M Valve L Leafl aflet Cr Cros
CFD mo model el showi
ng po polyhedral me mesh.
st Ru
Run Expe periment nt
st Ru
Run Lea eaflet Mo Motion
nd Ru
Run
nd Ru
Run Lea eaflet Mo Motion
Fa Far-Fi Field Pressure
Sid ide v e vie iew Top vi view Bottom v vie iew WSS
Biomedical FSI Applications
Carotid Arteries etc.
FSI & Overset Meshes