Biomedical subjects

Dimos Poulikakos

Publications and source records attributed to Dimos Poulikakos.

10 recordsLinked to original sources

PubMed2006

In-vivo flow simulation in coronary arteries based on computed tomography datasets: feasibility and initial results.

The purpose of this paper was to non-invasively assess hemodynamic parameters such as mass flow, wall shear stress (WSS), and wall pressure with computational fluid dynamics (CFD) in coronary arteries using patient-specific data from computed tomography (CT) angiography. Five patients (two without atherosclerosis, three with atherosclerosis) underwent retrospectively electrocardiogram (ECG) gated 16-detector row CT using ECG-pulsing and geometric models of coronary arteries were reconstructed for CFD analysis. Blood flow was considered laminar, incompressible, Newtonian, and pulsatile. The mass flow, WSS, and wall pressure were quantified and flow patterns were visualized. The wall pressure continuously decreased towards distal segments and showed pressure drops in stenotic segments. In coronary segments without atherosclerotic wall changes, WSS remained low, even during phases of high flow velocity, whereas in atherosclerotic vessels, the WSS was elevated already at low flow velocities. Stenoses and post-stenotic dilatations led to flow acceleration and rapid deceleration, respectively, including a distortion of flow. Areas of high WSS and high flow velocities were found adjacent to plaques, with values correlating with the degree of stenosis. CFD provided detailed mass flow measurements. CFD analysis is feasible in normal and atherosclerotic coronary arteries and provides the rationale for further investigation of the links between hemodynamic parameters and the significance of coronary stenoses.

Aged↗

PubMed2006

Computational investigation of subject-specific cerebrospinal fluid flow in the third ventricle and aqueduct of Sylvius.

The cerebrospinal fluid flow in the third ventricle of the brain and the aqueduct of Sylvius was studied using computational fluid dynamics (CFD) based on subject-specific boundary conditions derived from magnetic resonance imaging (MRI) scans. The flow domain geometry was reconstructed from anatomical MRI scans by manual image segmentation. The movement of the domain boundary was derived from MRI brain motion scans. Velocimetric MRI scans were used to reconstruct the velocity field at the inferior end of the aqueduct of Sylvius based on the theory of pulsatile flow in pipes. A constant pressure boundary condition was assigned at the foramina of Monro. Three main flow features were observed: a fluid jet emerging from the aqueduct of Sylvius, a moderately mobile recirculation zone above the jet and a mobile recirculation below the jet. The flow in the entire domain was laminar with a maximum Reynolds number of 340 in the aqueduct. The findings demonstrate that by combining MRI scans and CFD simulations, subject-specific detailed quantitative information of the flow field in the third ventricle and the aqueduct of Sylvius can be obtained.

Adult↗

PubMed2006

Three-dimensional modeling of mechanical forces in the extracellular matrix during epithelial lumen formation.

Mechanical interactions between cells and extracellular matrix (ECM) mediate epithelial cyst formation. This work relies on the combination of numerical modeling with live cell imaging, to piece together a novel nonintrusive method for determining three-dimensional (3D) mechanical forces caused by shape changes of a multicellular aggregate at the early stages of epithelial cyst formation. We analyzed the evolution of Madin-Darby canine kidney cells in 3D cultures using time-lapse microscopy, with type I collagen gel forming the ECM. The evolving 3D interface between the ECM and the cell aggregate was obtained from microscopy images, and the stress on the surface of a proliferating aggregate and in the surrounding ECM was calculated using the finite element method. The viscoelastic properties of the ECM (a needed input for the finite element method solver) were obtained through oscillatory shear flow experiments on a rheometer. For validation purpose, the forces exerted by an aggregate on a force-sensor array were measured and compared against the computational results.

Animals↗

PubMed2006

Measurement of the thermal conductivity of individual carbon nanotubes by the four-point three-omega method.

The thermal conductivity of individual multiwalled carbon nanotubes was measured by utilizing the four-point-probe third-harmonic method, based on the fact that the third harmonic amplitude and phase as a response to applied alternate current at fundamental frequency, omega, can be expressed in terms of thermal conductivity and diffusivity. To this end, a microfabricated device composed of four metal electrodes was modified to manufacture nanometer-sized wires by using a focused ion beam source. A carbon nanotube could then be suspended over a deep trench milled by the focused ion beam, preventing heat loss to the substrate. Compared with the two-point-probe technique, a significant improvement in accuracy is assured by using four probes, because the contact contribution to the determination of the thermal conductivity is eliminated, making it possible to measure the correct signals of first and third harmonics. The multiwalled carbon nanotube was modeled as a one-dimensional diffusive energy transporter and its thermal conductivity was measured at room temperature under vacuum to be 300 +/- 20 W/mK.

Algorithms↗

PubMed2006

Nanoparticle traffic on helical tracks: thermophoretic mass transport through carbon nanotubes.

Using molecular dynamics simulations, we demonstrate and quantify thermophoretic motion of solid gold nanoparticles inside carbon nanotubes subject to wall temperature gradients ranging from 0.4 to 25 K/nm. For temperature gradients below 1 K/nm, we find that the particles move "on tracks" in a predictable fashion as they follow unique helical orbits depending on the geometry of the carbon nanotubes. These findings markedly advance our knowledge of mass transport mechanisms relevant to nanoscale applications.

Computer Simulation↗

PubMed2005

Reconstruction of cerebrospinal fluid flow in the third ventricle based on MRI data.

A finite-volume model of the cerebrospinal fluid (CSF) system encompassing the third ventricle and the aqueduct of Sylvius was used to reconstruct CSF velocity and pressure fields based on MRI data. The flow domain geometry was obtained through segmentation of MRI brain anatomy scans. The movement of the domain walls was interpolated from brain motion MRI scans. A constant pressure boundary condition (BC) was specified at the foramina of Monro. A transient velocity BC reconstructed from velocimetric MRI scans was employed at the inferior end of the aqueduct of Sylvius. It could be shown that a combination of MRI scans and computational fluid dynamics (CFD) simulation can be used to reconstruct the flow field in the third ventricle. Pre-interventional knowledge of patient-specific CSF flow has the potential to improve neurosurgical interventions such as shunt placement in case of hydrocephalus.

Adult↗

PubMed2005

Computational simulation of a non-newtonian model of the blood separation process.

The aim of this work is to construct a computational fluid dynamics model capable of simulating the transient non-Newtonian process of apheresis. A Lagrangian-Eulerian model has been developed which tracks the blood particles within a two-dimensional flow configuration. Within the Eulerian method, the fluid mass and momentum conservation equations within the separator are solved using the density and the viscosity is calculated from the blood particle concentrations. Subsequently, the displacement of the blood particles is calculated with a Lagrangian method. Hawksley's model for the density of supensions is used in the variable density calculation. The viscosity is calculated with two models based on Vand's rigid particle suspension viscosity concepts, followed by the flow field calculation in the separator. Simulations were performed for various inlet hematocrit values and separator lengths. The simulations are in satisfactory agreement with experimental results reported in literature, indicating a complete separation of plasma and red blood cells (RBCs), as well as nearly complete separation of red blood cells and platelets. No hemolysis was observed in the simulations because the shear rate remained under the critical value of 150 N/m2.

Biophysical Phenomena↗

PubMed2005

Computational simulation of the blood separation process.

The aim of this work is to construct a computational fluid dynamics model capable of simulating the quasitransient process of apheresis. To this end a Lagrangian-Eulerian model has been developed which tracks the blood particles within a delineated two-dimensional flow domain. Within the Eulerian method, the fluid flow conservation equations within the separator are solved. Taking the calculated values of the flow field and using a Lagrangian method, the displacement of the blood particles is calculated. Thus, the local blood density within the separator at a given time step is known. Subsequently, the flow field in the separator is recalculated. This process continues until a quasisteady behavior is reached. The simulations show good agreement with experimental results. They shows a complete separation of plasma and red blood cells, as well as nearly complete separation of red blood cells and platelets. The white blood cells build clusters in the low concentrate cell bed.

Animals↗

PubMed2005

Computational modeling of the mechanical behavior of the cerebrospinal fluid system.

A computational fluid dynamics (CFD) model of the cerebrospinal fluid system was constructed based on a simplified geometry of the brain ventricles and their connecting pathways. The flow is driven by a prescribed sinusoidal motion of the third ventricle lateral walls, with all other boundaries being rigid. The pressure propagation between the third and lateral ventricles was examined and compared to data obtained from a similar geometry with a stenosed aqueduct. It could be shown that the pressure amplitude in the lateral ventricles increases in the presence of aqueduct stenosis. No difference in phase shift between the motion of the third ventricle walls and the pressure in the lateral ventricles because of the aqueduct stenosis could be observed. It is deduced that CFD can be used to analyze the pressure propagation and its phase shift relative to the ventricle wall motion. It is further deduced that only models that take into account the coupling between ventricles, which feature a representation of the original geometry that is as accurate as possible and which represent the ventricle boundary motion realistically, should be used to make quantitative statements on flow and pressure in the ventricular space.

Animals↗

PubMed2004

Computational simulation of intracoronary flow based on real coronary geometry.

OBJECTIVE: To assess the feasibility of computationally simulating intracoronary blood flow based on real coronary artery geometry and to graphically depict various mechanical characteristics of this flow. METHODS: Explanted fresh pig hearts were fixed using a continuous perfusion of 4% formaldehyde at physiological pressures. Omnipaque dye added to lead rubber solution was titrated to an optimum proportion of 1:25, to cast the coronary arterial tree. The heart was stabilized in a phantom model so as to suspend the base and the apex without causing external deformation. High resolution computerized tomography scans of this model were utilized to reconstruct the three-dimensional coronary artery geometry, which in turn was used to generate several volumetric tetrahedral meshes of sufficient density needed for numerical accuracy. The transient equations of momentum and mass conservation were numerically solved by employing methods of computational fluid dynamics under realistic pulsatile inflow boundary conditions. RESULTS: The simulations have yielded graphic distributions of intracoronary flow stream lines, static pressure drop, wall shear stress, bifurcation mass flow ratios and velocity profiles. The variability of these quantities within the cardiac cycle has been investigated at a temporal resolution of 1/100th of a second and a spatial resolution of about 10 microm. The areas of amplified variations in wall shear stress, mostly evident in the neighborhoods of arterial branching, seem to correlate well with clinically observed increased atherogenesis. The intracoronary flow lines showed stasis and extreme vorticity during the phase of minimum coronary flow in contrast to streamlined undisturbed flow during the phase of maximum flow. CONCLUSIONS: Computational tools of this kind along with a state-of-the-art multislice computerized tomography or magnetic resonance-based non-invasive coronary imaging, could enable realistic, repetitive, non-invasive and multidimensional quantifications of the effects of stenosis on distal hemodynamics, and thus help in precise surgical/interventional planning. It could also add insights into coronary and bypass graft atherogenesis.

Animals↗