Eric Paterson - Profile on Academia.edu (original) (raw)

Papers by Eric Paterson

Decomposition of Periodic Eddies with Varying Energy in a Turbulent Flow Using a Directional Averaging Technique

Ocean Engineering, Dec 1, 2015

Fluids, Jul 11, 2019

A hybrid Reynolds-averaged Navier Stokes/large-eddy simulation (RANS/LES) turbulence model integr... more A hybrid Reynolds-averaged Navier Stokes/large-eddy simulation (RANS/LES) turbulence model integrated with a transition formulation is developed and tested on a surrogate model problem through a joint experimental and computational fluid dynamic approach. The model problem consists of a circular cylinder for generating coherent unsteadiness and a downstream airfoil in the cylinder wake. The cylinder flow is subcritical, with a Reynolds number of 64,000 based upon the cylinder diameter. The quantitative dynamics of vortex shedding and Reynolds stresses in the cylinder near wake are well captured, owing to the turbulence-resolving large eddy simulation mode that was activated in the wake. The hybrid model switched between RANS and LES modes outside the boundary layers, as expected. According to the experimental and simulation results, the airfoil encountered local flow angle variations up to ±50 • . Further analysis through a phase-averaging technique found phase lags in the airfoil boundary layer along the chordwise locations, and both the phase-averaged and mean velocity profiles collapsed into the Law-of-the-wall in the range of 0 < y + < 50. The features of high blade-loading fluctuations due to unsteadiness and transitional boundary layers are of interest in the aerodynamic studies of full-scale wind turbine blades, making the current model problem a comprehensive benchmark case for future model development and validation.

Assessment of DES Models for Separated Flow From a Hump in a Turbulent Boundary Layer

Turbulent flow past the Glauert-Goldschmied body, a flow-control hump in a turbulent boundary lay... more Turbulent flow past the Glauert-Goldschmied body, a flow-control hump in a turbulent boundary layer, is studied numerically using detached-eddy simulation (DES), zonal detached-eddy simulation (ZDES), delayed detached-eddy simulation (DDES), and Reynolds-Averaged Navier-Stokes (RANS) modeling. The geometry is smooth so the downstream separation point is not set by facets of the geometry but is a function of the pressure gradient, a challenging condition for turbulence models. Comparisons to experimental data show that RANS with the Spalart-Allmaras turbulence model predicts the mean-field statistics well. The ZDES and DDES methods perform better than the DES formulation and are comparable to RANS in most statistics. An analysis of model behavior indicates that modeled stress depletion in the detached shear layer shortly after separation leads to loss of accuracy in the DES variants.

Development and Validation of a Computational Fluid Dynamics Methodology for Simulation of Pulsatile Left Ventricular Assist Devices

Asaio Journal, Mar 1, 2007

ABSTRACT An unsteady computational fluid dynamic methodology was developed so that design analyse... more ABSTRACT An unsteady computational fluid dynamic methodology was developed so that design analyses could be undertaken for devices such as the 50cc Penn State positive-displacement left ventricular assist device (LVAD). The piston motion observed in vitro was modeled, yielding the physiologic flow waveform observed during pulsatile experiments. Valve closure was modeled numerically by locally increasing fluid viscosity during the closed phase. Computational geometry contained Bjork-Shiley Monostrut mechanical heart valves in mitral and aortic positions. Cases for computational analysis included LVAD operation under steady-flow and pulsatile-flow conditions. Computations were validated by comparing simulation results with previously obtained in vitro particle image velocimetry (PIV) measurements. The steady portion of the analysis studied effects of mitral valve orientation, comparing the computational results with in vitro data obtained from mock circulatory loop experiments. The velocity field showed good qualitative agreement with the in vitro PIV data. The pulsatile flow simulations modeled the unsteady flow phenomena associated with a positive-displacement LVAD operating through several beat cycles. Flow velocity gradients allowed computation of the scalar wall strain rate, an important factor for determining hemodynamics of the device. Velocity magnitude contours compared well with PIV data throughout the cycle. Computational wall shear rates over the pulsatile cycle were found to be in the same range as wall shear rates observed in vitro.

Fluids, Dec 11, 2020

Electrochemical deposition (ECD) is a common method used in the field of microelectronics to grow... more Electrochemical deposition (ECD) is a common method used in the field of microelectronics to grow metallic coatings on an electrode. The deposition process occurs in an electrolyte bath where dissolved ions of the depositing material are suspended in an acid while an electric current is applied to the electrodes. The proposed computational model uses the finite volume method and the finite area method to predict copper growth on the plating surface without the use of a level set method or deforming mesh because the amount of copper layer growth is not expected to impact the fluid motion. The finite area method enables the solver to track the growth of the copper layer and uses the current density as a forcing function for an electric potential field on the plating surface. The current density at the electrolyte-plating surface interface is converged within each PISO (Pressure Implicit with Splitting Operator) loop iteration and incorporates the variance of the electrical resistance that occurs via the growth of the copper layer. This paper demonstrates the application of the finite area method for an ECD problem and additionally incorporates coupling between fluid mechanics, ionic diffusion, and electrochemistry.

Numerical and experimental study of the unsteady transitional boundary layer on a wind turbine airfoil

35th Wind Energy Symposium, 2017

System-Level Simulation of Floating Platform and Wind Turbine Using High-Fidelity and Engineering Models

Numerical Study of Transitional Unsteady Boundary Layer on Wind Turbine Airfoil Using Hybrid RANS/LES Turbulence Model

Fluids, 2020

The problem of simulating wakes in a stratified oceanic environment with active background turbul... more The problem of simulating wakes in a stratified oceanic environment with active background turbulence is considered. Anisotropic RANS turbulence models are tested against laboratory and eddy-resolving models of the problem. An important aspect of our work is to acknowledge that the environment is not quiescent; therefore, additional sources are included in the models to provide a non-zero background turbulence. The RANS models are found to reproduce some key features from the eddy-resolving and laboratory descriptions of the problem. Tests using the freestream sources show the intuitive result that background turbulence causes more rapid wake growth and decay.

Fluids, 2020

A modified set of governing differential equations for geophysical fluid flows is derived. All of... more A modified set of governing differential equations for geophysical fluid flows is derived. All of the simulation fields are decomposed into a nominal large-scale background state and a small-scale perturbation from this background, and the new system is closed by the assumption that the perturbation is one-way coupled to the background. The decomposition method, termed the multi-scale localized perturbation method (MSLPM), is then applied to the governing equations of stratified fluid flows, implemented in OpenFOAM, and exercised in order to simulate the interaction of a vertically-varying background shear flow with an axisymmetric perturbation in a turbulent ocean environment. The results demonstrate that the MSLPM can be useful in visualizing the evolution of a perturbation within a complex background while retaining the complex physics that are associated with the original governing equations. The simulation setup may also be simplified under the MSLPM framework. Further applicat...

Experiments in Fluids, 2018

Wind turbines with thick blade profiles experience turbulent, periodic approach flow, leading to ... more Wind turbines with thick blade profiles experience turbulent, periodic approach flow, leading to unsteady blade loading and large torque fluctuations on the turbine drive shaft. Presented here is an experimental study of a surrogate problem representing some key aspects of the wind turbine unsteady fluid mechanics. This experiment has been designed through joint consideration by experiment and computation, with the ultimate goal of numerical model development for aerodynamics in unsteady and turbulent flows. A cylinder at diameter Reynolds number of 65,000 and Strouhal number of 0.184 is placed 10.67 diameters upstream of a NACA 63215b airfoil with chord Reynolds number of 170,000 and chord-reduced frequency of k = 2 f c 2 ∕V = 1.5 . Extensive flow field measurements using particle image velocimetry provide a number of insights about this flow, as well as data for model validation and development. Velocity contours on the airfoil suction side in the presence of the upstream cylinder indicate a redistribution of turbulent normal stresses from transverse to streamwise, consistent with rapid distortion theory predictions. A study of the boundary layer over the suction side of the airfoil reveals very low Reynolds number turbulent mean streamwise velocity profiles. The dominance of the high amplitude large eddy passages results in a phase lag in streamwise velocity as a function of distance from the wall. The results and accompanying description provide a new test case incorporating moderate-reduced frequency inflow for computational model validation and development.

Numerical and experimental study of horizontal round turbulent forced plume in a static homogeneous environment

Heat and Mass Transfer, 2015

Detached-Eddy Simulation of High Reynolds Number Beveled-Trailing-Edge Flows and Wakes

Communications in Computational Physics, 2019

Feature identification is an important task in many fluid dynamics applications and diverse metho... more Feature identification is an important task in many fluid dynamics applications and diverse methods have been developed for this purpose. These methods are based on a physical understanding of the underlying behavior of the flow in the vicinity of the feature. Particularly, they rely on definition of suitable criteria (i.e. point-based or neighborhood-based derived properties) and proper selection of thresholds. For instance, among other techniques, vortex identification can be done through computing the Q-criterion or by considering the center of looping streamlines. However, these methods rely on creative visualization of physical idiosyncrasies of specific features and flow regimes, making them non-universal and requiring significant effort to develop. Here we present a physics-based, data-driven method capable of identifying any flow feature it is trained to. We use convolutional neural networks, a machine learning approach developed for image recognition, and adapt it to the problem of identifying flow features. The method was tested using mean flow fields from numerical simulations, where the recirculation region and boundary layer were identified in a two-dimensional flow through a convergent-divergent channel, and the horseshoe vortex was identified in three-dimensional flow over a wing-body junction. The novelty of the method is its ability to identify any type of feature, even distinguish between similar ones, without the need to explicitly define the physics (i.e. through development of suitable criterion and tunning of threshold). This provides a general method and removes the large burden placed on identifying new features. We expect this method can supplement existing techniques and allow for more automatic and discerning feature detection. The method can be easily extended to time-dependent flows, where it could be particularly impactful. For instance, it could be used in the identification of coherent structures in turbulent flows, a hindrance in the ongoing effort to establish a link between coherent structures and turbulence statistics.

Physical Review Fluids, 2017

Turbulence modeling is a critical component in numerical simulations of industrial flows based on... more Turbulence modeling is a critical component in numerical simulations of industrial flows based on Reynolds-averaged Navier-Stokes (RANS) equations. However, after decades of efforts in the turbulence modeling community, universally applicable RANS models with predictive capabilities are still lacking. Large discrepancies in the RANS-modeled Reynolds stresses are the main source that limits the predictive accuracy of RANS models. Identifying these discrepancies is of significance to possibly improve the RANS modeling. In this work, we propose a data-driven, physics-informed machine learning approach for reconstructing discrepancies in RANS modeled Reynolds stresses. The discrepancies are formulated as functions of the mean flow features. By using a modern machine learning technique based on random forests, the discrepancy functions are trained by existing DNS databases and then used to predict Reynolds stress discrepancies in different flows where data are not available. The proposed method is evaluated by two classes of flows: (1) fully developed turbulent flows in a square duct at various Reynolds numbers and (2) flows with massive separations. In separated flows, two training flow scenarios of increasing difficulties are considered: (1) the flow in the same periodic hills geometry yet at a lower Reynolds number, and (2) the flow in a different hill geometry with a similar recirculation zone. Excellent predictive performances were observed in both scenarios, demonstrating the merits of the proposed method.

Verification and Validation of CFD Simulation of Pulsating Laminar Flow in a Straight Pipe

17th AIAA Computational Fluid Dynamics Conference, 2005

In this paper, verification and validation analysis for pulsatile flow in a straight pipe is pres... more In this paper, verification and validation analysis for pulsatile flow in a straight pipe is presented. Numerical results were obtained using the finite volume method, and optimum grid density was determined using the Grid Convergence Index (GCI) calculation for the steady flow case. Based on the GCI calculation procedure we have developed a similar calculation for Time Convergence Index for the unsteady simulations. Transient computations were performed to obtain time dependent solutions of pulsating flow using three dierent time-steps and TCI was calculated to check time-step convergence. The numerical solution for the unsteady case was verified by comparing with accurate analytical solutions of Navier-Stokes equations for the pulsating flow in a pipe. Finally the numerical results were validated with experimental data. The results indicate that there is good agreement between the present results and analytical and experimental solutions. This paper stands as a good example for verification studies of unsteady CFD simulations where we have shown that TCI calculation can be applied in a similar procedure to grid convergence index calculation.

Comparison of LVAD Design Variants Using CFD Computed Thrombus Susceptibility Parameter

ASME 2008 Summer Bioengineering Conference, Parts A and B, 2008

The 70cc LionHeart was developed for use in patients weighing more than 70 kg, which prevents its... more The 70cc LionHeart was developed for use in patients weighing more than 70 kg, which prevents its use in many women, smaller men, and children. Scaled-down devices have been prone to thrombus formation during animal testing [1, 2]. It is suspected that the increased level of thrombosis is related to changes in the flow field when scaling from the 70cc to the 50cc device. Low fluid shear stresses and high residence times are traits that are believed to influence thrombus formation and deposition within artificial blood pumps [1].

Notices of the American Mathematical Society, 2014

ASAIO Journal, 2006

We developed a real-time remote managing system for an artificial heart using CDMA-based PDA phon... more We developed a real-time remote managing system for an artificial heart using CDMA-based PDA phone. It can both telemonitor and telecontrol the operating status of an artificial heart. The system consists of an artificial heart controller that contains cellular phone module, a PDA phone, and a gateway server computer that connects an heart controller and a PDA phone via TCP/IP network. Using this managing system, a medical staff can connect to the artificial heart controller remotely anytime, anywhere, see the operating status of an artificial heart, and adjust control parameters of the heart controller if needed. For security and reliability, several techniques are involved that checks communication error and user identification. Developed system showed satisfactory performance at in vitro and animal experiments. Using this mobile-based remote managing technique, more efficient outpatient treatment are possible. And also, the technique used for artificial heart remote managing can be applied to various medical fields.

Decomposition of Periodic Eddies with Varying Energy in a Turbulent Flow Using a Directional Averaging Technique

Ocean Engineering, Dec 1, 2015

Fluids, Jul 11, 2019

A hybrid Reynolds-averaged Navier Stokes/large-eddy simulation (RANS/LES) turbulence model integr... more A hybrid Reynolds-averaged Navier Stokes/large-eddy simulation (RANS/LES) turbulence model integrated with a transition formulation is developed and tested on a surrogate model problem through a joint experimental and computational fluid dynamic approach. The model problem consists of a circular cylinder for generating coherent unsteadiness and a downstream airfoil in the cylinder wake. The cylinder flow is subcritical, with a Reynolds number of 64,000 based upon the cylinder diameter. The quantitative dynamics of vortex shedding and Reynolds stresses in the cylinder near wake are well captured, owing to the turbulence-resolving large eddy simulation mode that was activated in the wake. The hybrid model switched between RANS and LES modes outside the boundary layers, as expected. According to the experimental and simulation results, the airfoil encountered local flow angle variations up to ±50 • . Further analysis through a phase-averaging technique found phase lags in the airfoil boundary layer along the chordwise locations, and both the phase-averaged and mean velocity profiles collapsed into the Law-of-the-wall in the range of 0 < y + < 50. The features of high blade-loading fluctuations due to unsteadiness and transitional boundary layers are of interest in the aerodynamic studies of full-scale wind turbine blades, making the current model problem a comprehensive benchmark case for future model development and validation.

Assessment of DES Models for Separated Flow From a Hump in a Turbulent Boundary Layer

Turbulent flow past the Glauert-Goldschmied body, a flow-control hump in a turbulent boundary lay... more Turbulent flow past the Glauert-Goldschmied body, a flow-control hump in a turbulent boundary layer, is studied numerically using detached-eddy simulation (DES), zonal detached-eddy simulation (ZDES), delayed detached-eddy simulation (DDES), and Reynolds-Averaged Navier-Stokes (RANS) modeling. The geometry is smooth so the downstream separation point is not set by facets of the geometry but is a function of the pressure gradient, a challenging condition for turbulence models. Comparisons to experimental data show that RANS with the Spalart-Allmaras turbulence model predicts the mean-field statistics well. The ZDES and DDES methods perform better than the DES formulation and are comparable to RANS in most statistics. An analysis of model behavior indicates that modeled stress depletion in the detached shear layer shortly after separation leads to loss of accuracy in the DES variants.

Development and Validation of a Computational Fluid Dynamics Methodology for Simulation of Pulsatile Left Ventricular Assist Devices

Asaio Journal, Mar 1, 2007

ABSTRACT An unsteady computational fluid dynamic methodology was developed so that design analyse... more ABSTRACT An unsteady computational fluid dynamic methodology was developed so that design analyses could be undertaken for devices such as the 50cc Penn State positive-displacement left ventricular assist device (LVAD). The piston motion observed in vitro was modeled, yielding the physiologic flow waveform observed during pulsatile experiments. Valve closure was modeled numerically by locally increasing fluid viscosity during the closed phase. Computational geometry contained Bjork-Shiley Monostrut mechanical heart valves in mitral and aortic positions. Cases for computational analysis included LVAD operation under steady-flow and pulsatile-flow conditions. Computations were validated by comparing simulation results with previously obtained in vitro particle image velocimetry (PIV) measurements. The steady portion of the analysis studied effects of mitral valve orientation, comparing the computational results with in vitro data obtained from mock circulatory loop experiments. The velocity field showed good qualitative agreement with the in vitro PIV data. The pulsatile flow simulations modeled the unsteady flow phenomena associated with a positive-displacement LVAD operating through several beat cycles. Flow velocity gradients allowed computation of the scalar wall strain rate, an important factor for determining hemodynamics of the device. Velocity magnitude contours compared well with PIV data throughout the cycle. Computational wall shear rates over the pulsatile cycle were found to be in the same range as wall shear rates observed in vitro.

Fluids, Dec 11, 2020

Electrochemical deposition (ECD) is a common method used in the field of microelectronics to grow... more Electrochemical deposition (ECD) is a common method used in the field of microelectronics to grow metallic coatings on an electrode. The deposition process occurs in an electrolyte bath where dissolved ions of the depositing material are suspended in an acid while an electric current is applied to the electrodes. The proposed computational model uses the finite volume method and the finite area method to predict copper growth on the plating surface without the use of a level set method or deforming mesh because the amount of copper layer growth is not expected to impact the fluid motion. The finite area method enables the solver to track the growth of the copper layer and uses the current density as a forcing function for an electric potential field on the plating surface. The current density at the electrolyte-plating surface interface is converged within each PISO (Pressure Implicit with Splitting Operator) loop iteration and incorporates the variance of the electrical resistance that occurs via the growth of the copper layer. This paper demonstrates the application of the finite area method for an ECD problem and additionally incorporates coupling between fluid mechanics, ionic diffusion, and electrochemistry.

Numerical and experimental study of the unsteady transitional boundary layer on a wind turbine airfoil

35th Wind Energy Symposium, 2017

System-Level Simulation of Floating Platform and Wind Turbine Using High-Fidelity and Engineering Models

Numerical Study of Transitional Unsteady Boundary Layer on Wind Turbine Airfoil Using Hybrid RANS/LES Turbulence Model

Fluids, 2020

The problem of simulating wakes in a stratified oceanic environment with active background turbul... more The problem of simulating wakes in a stratified oceanic environment with active background turbulence is considered. Anisotropic RANS turbulence models are tested against laboratory and eddy-resolving models of the problem. An important aspect of our work is to acknowledge that the environment is not quiescent; therefore, additional sources are included in the models to provide a non-zero background turbulence. The RANS models are found to reproduce some key features from the eddy-resolving and laboratory descriptions of the problem. Tests using the freestream sources show the intuitive result that background turbulence causes more rapid wake growth and decay.

Fluids, 2020

A modified set of governing differential equations for geophysical fluid flows is derived. All of... more A modified set of governing differential equations for geophysical fluid flows is derived. All of the simulation fields are decomposed into a nominal large-scale background state and a small-scale perturbation from this background, and the new system is closed by the assumption that the perturbation is one-way coupled to the background. The decomposition method, termed the multi-scale localized perturbation method (MSLPM), is then applied to the governing equations of stratified fluid flows, implemented in OpenFOAM, and exercised in order to simulate the interaction of a vertically-varying background shear flow with an axisymmetric perturbation in a turbulent ocean environment. The results demonstrate that the MSLPM can be useful in visualizing the evolution of a perturbation within a complex background while retaining the complex physics that are associated with the original governing equations. The simulation setup may also be simplified under the MSLPM framework. Further applicat...

Experiments in Fluids, 2018

Wind turbines with thick blade profiles experience turbulent, periodic approach flow, leading to ... more Wind turbines with thick blade profiles experience turbulent, periodic approach flow, leading to unsteady blade loading and large torque fluctuations on the turbine drive shaft. Presented here is an experimental study of a surrogate problem representing some key aspects of the wind turbine unsteady fluid mechanics. This experiment has been designed through joint consideration by experiment and computation, with the ultimate goal of numerical model development for aerodynamics in unsteady and turbulent flows. A cylinder at diameter Reynolds number of 65,000 and Strouhal number of 0.184 is placed 10.67 diameters upstream of a NACA 63215b airfoil with chord Reynolds number of 170,000 and chord-reduced frequency of k = 2 f c 2 ∕V = 1.5 . Extensive flow field measurements using particle image velocimetry provide a number of insights about this flow, as well as data for model validation and development. Velocity contours on the airfoil suction side in the presence of the upstream cylinder indicate a redistribution of turbulent normal stresses from transverse to streamwise, consistent with rapid distortion theory predictions. A study of the boundary layer over the suction side of the airfoil reveals very low Reynolds number turbulent mean streamwise velocity profiles. The dominance of the high amplitude large eddy passages results in a phase lag in streamwise velocity as a function of distance from the wall. The results and accompanying description provide a new test case incorporating moderate-reduced frequency inflow for computational model validation and development.

Numerical and experimental study of horizontal round turbulent forced plume in a static homogeneous environment

Heat and Mass Transfer, 2015

Detached-Eddy Simulation of High Reynolds Number Beveled-Trailing-Edge Flows and Wakes

Communications in Computational Physics, 2019

Feature identification is an important task in many fluid dynamics applications and diverse metho... more Feature identification is an important task in many fluid dynamics applications and diverse methods have been developed for this purpose. These methods are based on a physical understanding of the underlying behavior of the flow in the vicinity of the feature. Particularly, they rely on definition of suitable criteria (i.e. point-based or neighborhood-based derived properties) and proper selection of thresholds. For instance, among other techniques, vortex identification can be done through computing the Q-criterion or by considering the center of looping streamlines. However, these methods rely on creative visualization of physical idiosyncrasies of specific features and flow regimes, making them non-universal and requiring significant effort to develop. Here we present a physics-based, data-driven method capable of identifying any flow feature it is trained to. We use convolutional neural networks, a machine learning approach developed for image recognition, and adapt it to the problem of identifying flow features. The method was tested using mean flow fields from numerical simulations, where the recirculation region and boundary layer were identified in a two-dimensional flow through a convergent-divergent channel, and the horseshoe vortex was identified in three-dimensional flow over a wing-body junction. The novelty of the method is its ability to identify any type of feature, even distinguish between similar ones, without the need to explicitly define the physics (i.e. through development of suitable criterion and tunning of threshold). This provides a general method and removes the large burden placed on identifying new features. We expect this method can supplement existing techniques and allow for more automatic and discerning feature detection. The method can be easily extended to time-dependent flows, where it could be particularly impactful. For instance, it could be used in the identification of coherent structures in turbulent flows, a hindrance in the ongoing effort to establish a link between coherent structures and turbulence statistics.

Physical Review Fluids, 2017

Turbulence modeling is a critical component in numerical simulations of industrial flows based on... more Turbulence modeling is a critical component in numerical simulations of industrial flows based on Reynolds-averaged Navier-Stokes (RANS) equations. However, after decades of efforts in the turbulence modeling community, universally applicable RANS models with predictive capabilities are still lacking. Large discrepancies in the RANS-modeled Reynolds stresses are the main source that limits the predictive accuracy of RANS models. Identifying these discrepancies is of significance to possibly improve the RANS modeling. In this work, we propose a data-driven, physics-informed machine learning approach for reconstructing discrepancies in RANS modeled Reynolds stresses. The discrepancies are formulated as functions of the mean flow features. By using a modern machine learning technique based on random forests, the discrepancy functions are trained by existing DNS databases and then used to predict Reynolds stress discrepancies in different flows where data are not available. The proposed method is evaluated by two classes of flows: (1) fully developed turbulent flows in a square duct at various Reynolds numbers and (2) flows with massive separations. In separated flows, two training flow scenarios of increasing difficulties are considered: (1) the flow in the same periodic hills geometry yet at a lower Reynolds number, and (2) the flow in a different hill geometry with a similar recirculation zone. Excellent predictive performances were observed in both scenarios, demonstrating the merits of the proposed method.

Verification and Validation of CFD Simulation of Pulsating Laminar Flow in a Straight Pipe

17th AIAA Computational Fluid Dynamics Conference, 2005

In this paper, verification and validation analysis for pulsatile flow in a straight pipe is pres... more In this paper, verification and validation analysis for pulsatile flow in a straight pipe is presented. Numerical results were obtained using the finite volume method, and optimum grid density was determined using the Grid Convergence Index (GCI) calculation for the steady flow case. Based on the GCI calculation procedure we have developed a similar calculation for Time Convergence Index for the unsteady simulations. Transient computations were performed to obtain time dependent solutions of pulsating flow using three dierent time-steps and TCI was calculated to check time-step convergence. The numerical solution for the unsteady case was verified by comparing with accurate analytical solutions of Navier-Stokes equations for the pulsating flow in a pipe. Finally the numerical results were validated with experimental data. The results indicate that there is good agreement between the present results and analytical and experimental solutions. This paper stands as a good example for verification studies of unsteady CFD simulations where we have shown that TCI calculation can be applied in a similar procedure to grid convergence index calculation.

Comparison of LVAD Design Variants Using CFD Computed Thrombus Susceptibility Parameter

ASME 2008 Summer Bioengineering Conference, Parts A and B, 2008

The 70cc LionHeart was developed for use in patients weighing more than 70 kg, which prevents its... more The 70cc LionHeart was developed for use in patients weighing more than 70 kg, which prevents its use in many women, smaller men, and children. Scaled-down devices have been prone to thrombus formation during animal testing [1, 2]. It is suspected that the increased level of thrombosis is related to changes in the flow field when scaling from the 70cc to the 50cc device. Low fluid shear stresses and high residence times are traits that are believed to influence thrombus formation and deposition within artificial blood pumps [1].

Notices of the American Mathematical Society, 2014

ASAIO Journal, 2006

We developed a real-time remote managing system for an artificial heart using CDMA-based PDA phon... more We developed a real-time remote managing system for an artificial heart using CDMA-based PDA phone. It can both telemonitor and telecontrol the operating status of an artificial heart. The system consists of an artificial heart controller that contains cellular phone module, a PDA phone, and a gateway server computer that connects an heart controller and a PDA phone via TCP/IP network. Using this managing system, a medical staff can connect to the artificial heart controller remotely anytime, anywhere, see the operating status of an artificial heart, and adjust control parameters of the heart controller if needed. For security and reliability, several techniques are involved that checks communication error and user identification. Developed system showed satisfactory performance at in vitro and animal experiments. Using this mobile-based remote managing technique, more efficient outpatient treatment are possible. And also, the technique used for artificial heart remote managing can be applied to various medical fields.