Traditional CFD Boundary Conditions Applied to Blood Analog Flow through a Patient-Specific Aortic Coarctation

Author(s):  
Xiao Wang ◽  
D. Keith Walters ◽  
Greg W. Burgreen ◽  
David S. Thompson
2008 ◽  
Vol 130 (2) ◽  
Author(s):  
Matthew D. Ford ◽  
Hristo N. Nikolov ◽  
Jaques S. Milner ◽  
Stephen P. Lownie ◽  
Edwin M. DeMont ◽  
...  

Computational fluid dynamics (CFD) modeling of nominally patient-specific cerebral aneurysms is increasingly being used as a research tool to further understand the development, prognosis, and treatment of brain aneurysms. We have previously developed virtual angiography to indirectly validate CFD-predicted gross flow dynamics against the routinely acquired digital subtraction angiograms. Toward a more direct validation, here we compare detailed, CFD-predicted velocity fields against those measured using particle imaging velocimetry (PIV). Two anatomically realistic flow-through phantoms, one a giant internal carotid artery (ICA) aneurysm and the other a basilar artery (BA) tip aneurysm, were constructed of a clear silicone elastomer. The phantoms were placed within a computer-controlled flow loop, programed with representative flow rate waveforms. PIV images were collected on several anterior-posterior (AP) and lateral (LAT) planes. CFD simulations were then carried out using a well-validated, in-house solver, based on micro-CT reconstructions of the geometries of the flow-through phantoms and inlet/outlet boundary conditions derived from flow rates measured during the PIV experiments. PIV and CFD results from the central AP plane of the ICA aneurysm showed a large stable vortex throughout the cardiac cycle. Complex vortex dynamics, captured by PIV and CFD, persisted throughout the cardiac cycle on the central LAT plane. Velocity vector fields showed good overall agreement. For the BA, aneurysm agreement was more compelling, with both PIV and CFD similarly resolving the dynamics of counter-rotating vortices on both AP and LAT planes. Despite the imposition of periodic flow boundary conditions for the CFD simulations, cycle-to-cycle fluctuations were evident in the BA aneurysm simulations, which agreed well, in terms of both amplitudes and spatial distributions, with cycle-to-cycle fluctuations measured by PIV in the same geometry. The overall good agreement between PIV and CFD suggests that CFD can reliably predict the details of the intra-aneurysmal flow dynamics observed in anatomically realistic in vitro models. Nevertheless, given the various modeling assumptions, this does not prove that they are mimicking the actual in vivo hemodynamics, and so validations against in vivo data are encouraged whenever possible.


Author(s):  
Maria Nicole Antonuccio ◽  
Alessandro Mariotti ◽  
Benigno Marco Fanni ◽  
Katia Capellini ◽  
Claudio Capelli ◽  
...  

AbstractComputational Fluid Dynamics (CFD) simulations of blood flow are widely used to compute a variety of hemodynamic indicators such as velocity, time-varying wall shear stress, pressure drop, and energy losses. One of the major advances of this approach is that it is non-invasive. The accuracy of the cardiovascular simulations depends directly on the level of certainty on input parameters due to the modelling assumptions or computational settings. Physiologically suitable boundary conditions at the inlet and outlet of the computational domain are needed to perform a patient-specific CFD analysis. These conditions are often affected by uncertainties, whose impact can be quantified through a stochastic approach. A methodology based on a full propagation of the uncertainty from clinical data to model results is proposed here. It was possible to estimate the confidence associated with model predictions, differently than by deterministic simulations. We evaluated the effect of using three-element Windkessel models as the outflow boundary conditions of a patient-specific aortic coarctation model. A parameter was introduced to calibrate the resistances of the Windkessel model at the outlets. The generalized Polynomial Chaos method was adopted to perform the stochastic analysis, starting from a few deterministic simulations. Our results show that the uncertainty of the input parameter gave a remarkable variability on the volume flow rate waveform at the systolic peak simulating the conditions before the treatment. The same uncertain parameter had a slighter effect on other quantities of interest, such as the pressure gradient. Furthermore, the results highlight that the fine-tuning of Windkessel resistances is not necessary to simulate the post-stenting scenario.


2018 ◽  
Vol 13 (3) ◽  
pp. 1-10 ◽  
Author(s):  
I.Sh. Nasibullayev ◽  
E.Sh Nasibullaeva ◽  
O.V. Darintsev

The flow of a liquid through a tube deformed by a piezoelectric cell under a harmonic law is studied in this paper. Linear deformations are compared for the Dirichlet and Neumann boundary conditions on the contact surface of the tube and piezoelectric element. The flow of fluid through a deformed channel for two flow regimes is investigated: in a tube with one closed end due to deformation of the tube; for a tube with two open ends due to deformation of the tube and the differential pressure applied to the channel. The flow rate of the liquid is calculated as a function of the frequency of the deformations, the pressure drop and the physical parameters of the liquid.


Dynamics ◽  
2021 ◽  
Vol 1 (1) ◽  
pp. 9-17
Author(s):  
Andrea Natale Impiombato ◽  
Giorgio La Civita ◽  
Francesco Orlandi ◽  
Flavia Schwarz Franceschini Zinani ◽  
Luiz Alberto Oliveira Rocha ◽  
...  

As it is known, the Womersley function models velocity as a function of radius and time. It has been widely used to simulate the pulsatile blood flow through circular ducts. In this context, the present study is focused on the introduction of a simple function as an approximation of the Womersley function in order to evaluate its accuracy. This approximation consists of a simple quadratic function, suitable to be implemented in most commercial and non-commercial computational fluid dynamics codes, without the aid of external mathematical libraries. The Womersley function and the new function have been implemented here as boundary conditions in OpenFOAM ESI software (v.1906). The discrepancy between the obtained results proved to be within 0.7%, which fully validates the calculation approach implemented here. This approach is valid when a simplified analysis of the system is pointed out, in which flow reversals are not contemplated.


Author(s):  
D. Keith Walters ◽  
Greg W. Burgreen ◽  
Robert L. Hester ◽  
David S. Thompson ◽  
David M. Lavallee ◽  
...  

Computational fluid dynamics (CFD) simulations were performed for unsteady periodic breathing conditions, using large-scale models of the human lung airway. The computational domain included fully coupled representations of the orotracheal region and large conducting zone up to generation four (G4) obtained from patient-specific CT data, and the small conducting zone (to G16) obtained from a stochastically generated airway tree with statistically realistic geometrical characteristics. A reduced-order geometry was used, in which several airway branches in each generation were truncated, and only select flow paths were retained to G16. The inlet and outlet flow boundaries corresponded to the oronasal opening (superior), the inlet/outlet planes in terminal bronchioles (distal), and the unresolved airway boundaries arising from the truncation procedure (intermediate). The cyclic flow was specified according to the predicted ventilation patterns for a healthy adult male at three different activity levels, supplied by the whole-body modeling software HumMod. The CFD simulations were performed using Ansys FLUENT. The mass flow distribution at the distal boundaries was prescribed using a previously documented methodology, in which the percentage of the total flow for each boundary was first determined from a steady-state simulation with an applied flow rate equal to the average during the inhalation phase of the breathing cycle. The distal pressure boundary conditions for the steady-state simulation were set using a stochastic coupling procedure to ensure physiologically realistic flow conditions. The results show that: 1) physiologically realistic flow is obtained in the model, in terms of cyclic mass conservation and approximately uniform pressure distribution in the distal airways; 2) the predicted alveolar pressure is in good agreement with previously documented values; and 3) the use of reduced-order geometry modeling allows accurate and efficient simulation of large-scale breathing lung flow, provided care is taken to use a physiologically realistic geometry and to properly address the unsteady boundary conditions.


2013 ◽  
Vol 8-9 ◽  
pp. 225-234
Author(s):  
Dalia Sabina Cimpean

The present study is focused on the mixed convection fluid flow through a porous medium, when a different amount of nanoparticles is added in the base fluid. The nanofluid saturates the porous matrix and different situations of the flow between two walls are presented and discussed. Alternatively mathematical models are presented and discussed. A solution of a system which contains the momentum, Darcy and energy equations, together with the boundary conditions involved, is given. The behavior of different nanofluids, such thatAu-water, Ag-waterandFe-wateris graphically illustrated and compared with the previous results.The research target is to observe the substantial increase of the thermophysical fluid properties, when the porous medium issaturated by a nanofluid instead of a classical Newtonian fluid.


Sign in / Sign up

Export Citation Format

Share Document