Macroscale versus microscale methods for physiological analysis of biofilms formed in 96-well microtiter plates

Microtiter plates with 96 wells have become one of the preferred platforms for biofilm studies mainly because they enable high-throughput assays. In this work, macroscale and microscale methods were used to study the impact of hydrodynamic conditions on the physiology and location of Escherichia coli JM109(DE3) biofilms formed in microtiter plates. Biofilms were formed in shaking and static conditions, and two macroscale parameters were assayed: the total amount of biofilm was measured by the crystal violet assay and the metabolic activity was determined by the resazurin assay. From the macroscale point of view, there were no statistically significant differences between the biofilms formed in static and shaking conditions. However, at a microscale level, the differences between both conditions were revealed using scanning electron microscopy (SEM). It was observed that biofilm morphology and spatial distribution along the wall were different in these conditions. Simulation of the hydrodynamic conditions inside the wells at a microscale was performed by computational fluid dynamics (CFD). These simulations showed that the shear strain rate was unevenly distributed on the walls during shaking conditions and that regions of higher shear strain rate were obtained closer to the air/liquid interface. Additionally, it was shown that wall regions subjected to higher shear strain rates were associated with the formation of biofilms containing cells of smaller size. Conversely, regions with lower shear strain rate were prone to have a more uniform spatial distribution of adhered cells of larger size. The results presented on this work highlight the wealth of information that may be gathered by complementing macroscale approaches with a microscale analysis of the experiments.     

Variations in mass transfer to single endothelial cells

Mass transfer between flowing blood and arterial mural cells (including vascular endothelial cells) may play an important role in atherogenesis. Endothelial cells are known to have an apical surface topography that is not flat, and hence mass transfer patterns to individual endothelial cells are likely affected by the local cellular topography. The purpose of this paper is to investigate the relationship between vascular endothelial cell surface topography and cellular level mass transfer. Confluent porcine endothelial monolayers were cultured under both shear and static conditions and atomic force microscopy was used to measure endothelial cell topography. Using finite element methods and the measured cell topography, flow and concentration fields were calculated for a typical, small, blood-borne solute. A relative Sherwood number was defined as the difference between the computed Sherwood number and that predicted by the Leveque solution for mass transfer over a flat surface: this eliminates the effects of axial location on mass transfer efficiency. The average intracellular relative Sherwood number range was found to be dependent on cell height and not dependent on cell elongation due to shear stress in culture. The mass flux to individual cells reached a maximum at the highest point on the endothelial cell surface, typically corresponding to the nucleus of the cell. Therefore, for small receptor-mediated solutes, increased solute uptake efficiency can be achieved by concentrating receptors near the nucleus. The main conclusion of the work is that although the rate of mass transfer varies greatly over an individual cell, the average mass transfer rate to a cell is close to that predicted for a flat cell. In comparison to other hemodynamic factors, the topography of endothelial cells therefore seems to have little effect on mass transfer rates and is likely physiologically insignificant.     

Characterizing the fluid dynamics in the flow fields of cylindrical orbitally shaken bioreactors with different geometry sizes

Orbitally shaken bioreactors (OSRs) are commonly used for the cultivation of mammalian cells in suspension. To aid the geometry designing and optimizing of OSRs, we conducted a three‐dimensional computational fluid dynamics (CFD) simulation to characterize the flow fields in a 10 L cylindrical OSR with different vessel diameters. The liquid wave shape captured by a camera experimentally validated the CFD models established for the cylindrical OSR. The geometry size effect on volumetric mass transfer coefficient (k_La) and hydromechanical stress was analyzed by varying the ratio of vessel diameter (d) to liquid height at static (h_L), d/h_L. The highest value of k_La about 30 h^?1 was observed in the cylindrical vessel with the d/h_L of 6.35. Moreover, the magnitudes of shear stress and energy dissipation rate in all the vessels tested were below their minimum values causing cells damage separately, which indicated that the hydromechanical‐stress environment in OSRs is suitable for cells cultivation in suspension. Finally, the CFD results suggested that the d/h_L higher than 8.80 should not be adopted for the 10 L cylindrical OSR at the shaking speed of 180 rpm because the “out of phase” state probably will happen there.     

Characterisation of operation conditions and online monitoring of physiological culture parameters in shaken 24-well microtiter plates

A new online monitoring technique to measure the physiological parameters, dissolved oxygen (DO) and pH of microbial cultures in continuously shaken 24-well microtiter plates (MTP) is introduced. The new technology is based on immobilised fluorophores at the bottom of standard 24-well MTPs. The sensor MTP is installed in a sensor dish reader, which can be fixed on an orbital shaker. This approach allows real online measurements of physiological parameters during continuous shaking of cultures without interrupting mixing and mass transfer like currently available technologies do. The oxygen transfer conditions at one constant shaking frequency (250 1/min) and diameter (25 mm) was examined with the chemical sulphite oxidation method. Varied filling volumes (600–1,200 μL) of Escherichia coli cultures demonstrated the importance of sufficient oxygen transfer to the culture. Cultures with higher filling volumes were subjected to an oxygen limitation, which influenced the cell metabolism and prolongated the cultivation time. The effects could be clearly monitored by online DO and pH measurements. A further study of different media in an E. coli fermentation elucidated the different growth behaviour in response to the medium composition. The MTP fermentations correlated very well with parallel fermentations in shake flasks. The new technique gives valuable new insights into biological processes at a very small scale, thus enabling parallel experimentation and shorter development times in bioprocessing.     

Oxygen transfer phenomena in 48-well microtiter plates: determination by optical monitoring of sulfite oxidation and verification by real-time measurement during microbial growth

Oxygen limitation is one of the most frequent problems associated with the application of shaking bioreactors. The gas-liquid oxygen transfer properties of shaken 48-well microtiter plates (MTPs) were analyzed at different filling volumes, shaking diameters, and shaking frequencies. On the one hand, an optical method based on sulfite oxidation was used as a chemical model system to determine the maximum oxygen transfer capacity (OTR(max)). On the other hand, the Respiration Activity Monitoring System (RAMOS) was applied for online measurement of the oxygen transfer rate (OTR) during growth of the methylotropic yeast Hansenula polymorpha. A proportionality constant between the OTR(max) of the biological system and the OTR(max) of the chemical system were indicated from these data, offering the possibility to transform the whole set of chemical data to biologically relevant conditions. The results exposed "out of phase" shaking conditions at a shaking diameter of 1 mm, which were confirmed by theoretical consideration with the phase number (Ph). At larger shaking diameters (2-50 mm) the oxygen transfer rate in MTPs shaken at high frequencies reached values of up to 0.28 mol/L/h, corresponding to a volumetric mass transfer coefficient (k(L)a) of 1,600 1/h. The specific mass transfer area (a) increases exponentially with the shaking frequency up to values of 2,400 1/m. On the contrary, the mass transfer coefficient (k(L)) is constant at a level of about 0.15 m/h over a wide range of shaking frequencies and shaking diameters. However, at high shaking frequencies, when the complete liquid volume forms a thin film on the cylindric wall of the well, the mass transfer coefficient (k(L)) increases linearly to values of up to 0.76 m/h. Essentially, the present investigation demonstrates that the 48-well plate outperforms the 96-well MTP and shake flasks at widely used operating conditions with respect to oxygen supply. The 48-well plates emerge, therefore, as an excellent alternative for microbial cultivation and expression studies combining the advantages of both the high-throughput 96-well MTP and the classical shaken Erlenmeyer flask.     

Determination of the average shear rate in a stirred and aerated tank bioreactor

A method for evaluating the average shear rate () in a stirred and aerated tank bioreactor has been proposed for non-Newtonian fluids. The volumetric oxygen transfer coefficient (k _L a) was chosen as the appropriate characteristic parameter to evaluate the average shear rate (). The correlations for the average shear rate as a function of N and rheological properties of the fluid (K and n) were obtained for two airflow rate conditions (ϕ_air). The shear rate values estimated by the proposed methodology lay within the range of the values calculated by classical correlations. The proposed correlations were utilized to predict the during the Streptomyces clavuligerus cultivations carried out at 0.5 vvm and four different rotational impeller speeds. The results show that the values of the average shear rate () varied from 437 to 2,693 s⁻¹ by increasing with N and flow index (n) and decreasing with the fluid consistency index (K).     

Role of diastolic vortices in flow and energy dynamics during systolic ejection

MRI-based computational fluid dynamics simulations were performed in the left ventricles of two adult porcine subjects with varying physiological states (before and after an induced infarction). The hypothesis that diastolic vortices store kinetic energy and assist systolic ejection was tested, by performing systolic simulations in the presence and absence of diastolic vortices. The latter was achieved by reinitializing the entire velocity field to be zero at the beginning of systole. A rudimentary prescribed motion model of a mitral valve was included in the simulations to direct the incoming mitral jet towards the apex. Results showed that the presence or absence of diastolic vortex rings had insignificant impact on the energy expended by walls of the left ventricles for systolic ejection for both the porcine subjects, under all physiological conditions. Although substantial kinetic energy was stored in diastolic vortices by end diastole, it provided no appreciable savings during systolic ejection, and most likely continued to complete dissipation during systole. The role of diastolic vortices in apical washout was investigated by studying the cumulative mass fraction of passive dye that was ejected during systole in the presence and absence of vortices. Results indicated that the diastolic vortices play a crucial role in ensuring efficient washout of apical blood during systolic ejection.     

Development and evaluation of a 24 well microtitre plate method for isolation of Listeria spp. or Listeria monocytogenes from foods

A rapid and simple method (24M) using 24 well microtitre plates was developed to determine the presence of Listeria monocytogenes or Listeria spp. in food samples. The 24M was composed of two 24 well microtitre plates connected with a yellow tip. The 24M was evaluated with pathogen cocktails and ground beef samples and compared with the conventional method for presumptive identification of Listeria spp. Only food-borne pathogen cocktails and ground beef samples containing L. monocytogenes or Listeria spp. showed a positive reaction in 24M after 24 h incubation at 35 degrees C. Test results were the same with the conventional method and the 24M method and showed high efficiency for recovery of Listeria spp. from foods. This new, convenient and economical method can isolate Listeria spp. simultaneously from 24 different food samples.     

The physical characterisation of a microscale parallel bioreactor platform with an industrial CHO cell line expressing an IgG4

There is a growing body of evidence that the ambr™ workstation from TAP Biosystems performs well in terms of helping to select appropriate clones for scale-up studies. Here we have investigated the physical characteristics of this microscale bioreactor system and found that these are quite different from those that exist in larger scale stirred bioreactors. For example, the flow regime in the ambr™ vessel is transitional rather than turbulent and the sparged air/oxygen superficial gas velocity is relatively very low whilst the specific power input is much higher (~400 W/m³) when compared to that used at larger scales (typically ~20 W/m³). This specific power input is necessary in order to achieve k_La values sufficiently high to satisfy the oxygen demand of the cells and control of dO₂. In line with other studies, we find that the culture of CHO cells in a 15 mL ambr™ bioreactor gave similar cell growth and productivity to that achieved in a 5 L stirred bioreactor whilst the results from shake flasks were significantly different. Given the differences in physical characteristics between the ambr™ and larger stirred bioreactors, we suggest that this similarity in biological performance is due to their similar control capabilities and the ‘equivalence of the stress parameters’ across the scales when compared with shake flasks.     

Computational fluid dynamics modeling of mass-transfer behavior in a bioreactor for hairy root culture. II. Analysis of ultrasound-intensified process

Recently, cichoric acid production from hairy roots of Echinacea purpurea was significantly improved by ultrasound stimulation in an airlift bioreactor. In this article, the possible mechanism on ultrasound-intensified hairy root culture of E. purpurea in the bioreactor was elucidated with the help of computational fluid dynamics (CFD) simulation, membrane permeability detection, dissolved oxygen concentration detection, confocal laser-scanning microscopy (LSM) observation, and phenylalanine ammonium lyase (PAL) activity analysis. The CFD model developed in Part I was used to simulate the hydrodynamics and oxygen mass transfer in hairy root bioreactor culture stimulated by ultrasound. A dynamic mesh model combined with a changing Schmidt number method was used for the simulation of the ultrasound field. Simulation results and experimental data illustrated that ultrasound intensified oxygen mass transfer in the hairy root clump, which subsequently stimulated root growth and cichoric acid biosynthesis. Ultrasound increased the hairy root membrane permeability, and a high root membrane permeability of 0.359 h(-1) was observed at the bottom region in the bioreactor. LSM observation showed that the change in the membrane permeability recovered to normal in the further culture after ultrasound stimulation. PAL activity in the hairy roots was stimulated by ultrasound increase and was correlated well to cichoric acid accumulation in the hairy roots of E. purpurea.     

Using CFD simulations and statistical analysis to correlate oxygen mass transfer coefficient to both geometrical parameters and operating conditions in a stirred-tank bioreactor

Optimization of a bioreactor design can be an especially challenging process. For instance, testing different bioreactor vessel geometries and different impeller and sparger types, locations, and dimensions can lead to an exceedingly large number of configurations and necessary experiments. Computational fluid dynamics (CFD), therefore, has been widely used to model multiphase flow in stirred-tank bioreactors to minimize the number of optimization experiments. In this study, a multiphase CFD model with population balance equations are used to model gas–liquid mixing, as well as gas bubble distribution, in a 50 L single-use bioreactor vessel. The vessel is the larger chamber in an early prototype of a multichamber bioreactor for mammalian cell culture. The model results are validated with oxygen mass transfer coefficient (k_La) measurements within the prototype. The validated model is projected to predict the effect of using ring or pipe spargers of different sizes and the effect of varying the impeller diameter on k_La. The simulations show that ring spargers result in a superior k_La compared to pipe spargers, with an optimum sparger-to-impeller diameter ratio of 0.8. In addition, larger impellers are shown to improve k_La. A correlation of k_La is presented as a function of both the reactor geometry (i.e., sparger-to-impeller diameter ratio and impeller-to-vessel diameter ratio) and operating conditions (i.e., Reynolds number and gas flow rate). The resulting correlation can be used to predict k_La in a bioreactor and to optimize its design, geometry, and operating conditions.     

Torsional motion of the left ventricle does not affect ventricular fluid dynamics of both foetal and adult hearts

Left ventricular torsion is caused by shortening and relaxation of the helical fibres in the myocardium, and is thought to be an optimal configuration for minimizing myocardial tissue strains. Characteristics of torsional motion has also been proposed to be markers for cardiac dysfunction. However, its effects on fluid and energy dynamics in the left ventricle have not been comprehensively investigated. To investigate this, we performed image-based flow simulations on five healthy adult porcine and two healthy human foetal left ventricles (representing two different length scales) at different degrees of torsional motions. In the adult porcine ventricles, cardiac features such as papillary muscles and mitral valves, and cardiac conditions such as myocardial infarctions, were also included to investigate the effect of twist. The results showed that, for all conditions investigated, ventricular torsional motion caused minimal changes to flow patterns, and consistently accounted for less than 2% of the energy losses, wall shear stresses, and ejection momentum energy. In contrast, physiological characteristics such as chamber size, stroke volume and heart rate had a much greater influence on flow patterns and energy dynamics. The results thus suggested that it might not be necessary to model the torsional motion to study the flow and energy dynamics in left ventricles.     

Influence of sparger on energy dissipation, shear rate, and mass transfer to sea water in a concentric-tube airlift bioreactor

Data on volumetric mass-transfer coefficient, K_La_L, in a 12 × 10⁻³ m³ airlift bioreactor are reported. Measurements were made in sea water. The superficial gas velocity ranged up to 0.21 m/s. Four cylindrical spargers (60–1000 μm pore size) were tested. In bubbly flow, the sparger pore size strongly influenced the K_La_L; the highest K_La_L values were obtained with the smallest pore size. In contrast, in the transition and heterogeneous flow regimes, the pore size had little influence on K_La_L. The best correlation of the mass transfer data was obtained when both gas holdup and liquid superficial velocity were taken as independent variables. Shear rates were estimated in the different zones of the reactor. The highest values were found in the bottom zone of the reactor and in the gas-liquid separator. The penetration and isotropic turbulence models were used to develop a semi-theoretical equation relating the volumetric mass-transfer coefficient to shear rate; hence providing a better understanding of how the operational variables may be manipulated to attain a moderate shear rate and an appropriate level of mass transfer, two extremely important parameters for the growth of sensible microorganisms as those used in marine biotechnology.     

Using 96-well tissue culture polystyrene plates and a fluorescence plate reader as tools to study the survival and inactivation of viruses on surfaces

A method for studying the behavior of viruses on surfaces has been developed and is illustrated by determining the temperatures that inactivate adsorbed viral hemorrhagic septicemia virus (VHSV) and the concentration of 1-propanol that disinfected surfaces with adsorbed VHSV and chum salmon virus (CSV). VHSV is a rhabdovirus; CSV, a reovirus, and they were detected with two fish cell lines, EPC and CHSE-214, respectively. When polystyrene tissue culture surfaces were incubated with virus, rinsed, and left to dry, they still supported the attachment and spreading of cell lines and after 7 days these cells showed the characteristic CPE of the viruses. Thus cells appeared to be infected directly from surfaces on which viruses had been adsorbed. Applying this property to 96-well plates allowed duplicate surfaces to be examined for their infectiousness or support of CPE. For each treatment 80 replicate surfaces in a 96-well plate were tested at one time and the results expressed as the number of wells showing CPE. VHSV adsorbed to polystyrene was inactivated by drying in the dark at temperatures above 14 °C, but remained infectious for at least 15 days of drying at 4 °C. For chemical sterilization of polystyrene surfaces with adsorbed virus, disinfection was achieved with 1-propanol at 40% for VHSV and at 60% for CSV. As CPE can be conveniently monitored in 96-well plates with a fluorescence plate reader, this method can be used to rapidly evaluate a variety of treatments for their ability to inactivate surface-bound viruses.     

Haemodynamic effects of stent diameter and compaction ratio on flow-diversion treatment of intracranial aneurysms: A numerical study of a successful and an unsuccessful case

Background: Compacting a flow-diverting (FD) stent is an emerging technique to create a denser configuration of wires across the aneurysm ostium. However, quantitative analyses of post-stenting haemodynamics affected by the compaction level of different stent sizes remain inconclusive.Objective: To compare the aneurysmal haemodynamic alterations after virtual FD treatments with different device diameters at different compaction ratios.Methods: We virtually implanted three sizes of FD stent, with each size deployed at four compaction ratios, into two patient aneurysms previously treated with the Silk + FD—one successful case and the other unsuccessful. Wire configurations of the FD in the 24 treatment scenarios were examined, and aneurysmal haemodynamic alterations were resolved by computational fluid dynamics (CFD) simulations. We investigated the aneurysmal flow patterns, aneurysmal average velocity (AAV), mass flowrate (MF), and energy loss (EL) in each scenario.Results: Compactions of the stent in the successful case resulted in a greater metal coverage rate than that achieved in the unsuccessful one. A 25% increment in compaction ratio further decreased the AAV (12%), MF (11%), and EL (9%) in both cases (average values). The averaged maximum differences attributable to device size were 10% (AAV), 8% (MF), and 9% (EL).Conclusions: Both stent size and compaction level could markedly affect the FD treatment outcomes. It is therefore important to individualise the treatment plan by selecting the optimal stent size and deployment procedure. CFD simulation can be used to investigate the treatment outcomes, thereby assisting doctors in choosing a favourable treatment plan.     

Use of computational fluid dynamics as a tool for establishing process design space for mixing in a bioreactor

The concept of "design space" plays an integral part in implementation of quality by design for pharmaceutical products. ICH Q8 defines design space as "the multidimensional combination and interaction of input variables (e.g., material attributes) and process parameters that have been demonstrated to provide assurance of quality. Working within the design space is not considered as a change. Movement out of the design space is considered to be a change and would normally initiate a regulatory post-approval change process. Design space is proposed by the applicant and is subject to regulatory assessment and approval." Computational fluid dynamics (CFD) is increasingly being used as a tool for modeling of hydrodynamics and mass transfer. In this study, a laboratory-scale aerated bioreactor is modeled using CFD. Eulerian-Eulerian multiphase model is used along with dispersed k-ε turbulent model. Population balance model is incorporated to account for bubble breakage and coalescence. Multiple reference frame model is used for the rotating region. We demonstrate the usefulness of CFD modeling for evaluating the effects of typical process parameters like impeller speed, gas flow rate, and liquid height on the mass transfer coefficient (k(L)a). Design of experiments is utilized to establish a design space for the above mentioned parameters for a given permissible range of k(L)a.     

Mathematical modeling and study of mass transfer parameters in supercritical fluid extraction of fatty acids from Trout powder

Mathematical modeling of fatty acids (FAs) extraction from Trout powder by supercritical carbon dioxide was performed in the present work. Trout powder with its low cost contains high amount of essential FAs and it is commonly available as a proper source of FAs. The effect of process parameters, such as pressure (25, 28, 31, 34 and 37 MPa) and temperature (310, 318 and 326 K) of extraction and void fraction of the bed (0.25, 0.35 and 0.45, v/v) on the yield of FAs extraction was examined in a series of experiments conducted in a laboratory scale apparatus. The results indicated a significant increase of extraction yield with an increase of pressure from 25 to 34 MPa, but working at the higher pressure (37 MPa) caused reduction of the extract. Increasing the temperature higher than 318 K revealed significant reduction of the FAs yield and increasing the bed void fraction from 0.25 to 0.45 showed enhancement of the extraction.

The mathematical model was developed considering diffusion-controlled regime in the particle and film mass transfer resistance around the particle with axial dispersion of the bulk phase at dynamic conditions. Henry law was used to describe the equilibrium state of solid and fluid phases. The proposed mass balance equations were numerically solved using implicit finite difference method and the model parameters were correlated using the experimental results of the outlet FAs concentration in the oil extracted at dynamic conditions. Well-known Nelder–Mead method was applied to estimate the four parameters of the model, namely, mass transfer coefficient (k_f), axial dispersion coefficient (D_ax) in the bulk phase, effective diffusivity (D_eff) into the pores and Henry coefficient (H). In the range of studied conditions, the higher extraction efficiency with higher pressure resulted lower correlated H, although the temperature increasing which showed a retrograde phenomena in the FAs yield, revealed H passing though a minimum.     

Environmentally low-temperature kinetic and thermodynamic study of lactate dehydrogenase from Atlantic cod (G. morhua) using a 96-well microplate technique

Analyses of temperature-dependent kinetic parameters in enzymes extracted from tissues of ectothermic animals are usually carried out within the range of physiological temperatures (0-40 degrees C). However, multisample spectrophotometers (so-called microplate readers) with efficient wide-range temperature control (including cooling) have previously been unavailable. This limits the statistical quality of the measurements. A temperature-controlled microplate was designed for a 96-well microplate reader to overcome this limitation. This so-called T-microplate is able to control assay temperature between the freezing point of a liquid sample and 60 degrees C with high stability and accuracy in any data acquisition mode. At 4 degrees C the accuracy of the temperature control was +/-0.1 degrees C and temperature homogeneity across the microplate was +/-0.3 degrees C. As examples, analyses of the temperature dependence of Michaelis-Menten (K'(PYR)(m) and substrate inhibition (K'(PYR)(si) constants for pyruvate, of the maximal rate of reaction (V'(max), of the apparent Arrhenius activation energy (E(A), and of the Gibbs free-energy change (deltaG) of lactate dehydrogenases from muscle of Atlantic cod Gadus morhua acclimated to 4 degrees C are described. The large dataset obtained allowed evaluation of a new mechanism of metabolic compensation in response to seasonal temperature change.     

CFD of mixing of multi‐phase flow in a bioreactor using population balance model

Mixing in bioreactors is known to be crucial for achieving efficient mass and heat transfer, both of which thereby impact not only growth of cells but also product quality. In a typical bioreactor, the rate of transport of oxygen from air is the limiting factor. While higher impeller speeds can enhance mixing, they can also cause severe cell damage. Hence, it is crucial to understand the hydrodynamics in a bioreactor to achieve optimal performance. This article presents a novel approach involving use of computational fluid dynamics (CFD) to model the hydrodynamics of an aerated stirred bioreactor for production of a monoclonal antibody therapeutic via mammalian cell culture. This is achieved by estimating the volume averaged mass transfer coefficient (k_La) under varying conditions of the process parameters. The process parameters that have been examined include the impeller rotational speed and the flow rate of the incoming gas through the sparger inlet. To undermine the two‐phase flow and turbulence, an Eulerian‐Eulerian multiphase model and k‐ε turbulence model have been used, respectively. These have further been coupled with population balance model to incorporate the various interphase interactions that lead to coalescence and breakage of bubbles. We have successfully demonstrated the utility of CFD as a tool to predict size distribution of bubbles as a function of process parameters and an efficient approach for obtaining optimized mixing conditions in the reactor. The proposed approach is significantly time and resource efficient when compared to the hit and trial, all experimental approach that is presently used. © 2016 American Institute of Chemical Engineers Biotechnol. Prog., 32:613–628, 2016     

Computational simulations of flow dynamics and blood damage through a bileaflet mechanical heart valve scaled to pediatric size and flow

Despite pressing needs, there are currently no FDA approved prosthetic valves available for use in the pediatric population. This study is performed for predictive assessment of blood damage in bileaflet mechanical heart valves (BMHVs) with pediatric sizing and flow conditions. A model of an adult-sized 23 mm St. Jude Medical (SJM) Regent^™ valve is selected for use in simulations, which is scaled in size for a 5-year old child and 6-month old infant. A previously validated lattice-Boltzmann method (LBM) is used to simulate pulsatile flow with thousands of suspended platelets for cases of adult, child, and infant BMHV flows. Adult BMHV flows demonstrate more disorganized small-scale flow features, but pediatric flows are associated with higher fluid shear stresses. Platelet damage in the pediatric cases is higher than in adult flow, highlighting thrombus complication dangers of pediatric BMHV flows. This does not necessarily suggest clinically important differences in thromboembolic potential. Highly damaged platelets in pediatric flows are primarily found far downstream of the valve, as there is less flow recirculation in pediatric flows. In addition, damage levels are well below expected thresholds for platelet activation. The extent of differences here documented between the pediatric and adult cases is of concern, demanding particular attention when pediatric valves are designed and manufactured. However, the differences between the pediatric and adult cases are not such that development of pediatric sized valves is untenable. This study may push for eventual approval of prosthetic valves resized for the pediatric population. Further studies will be necessary to determine the validity and potential thrombotic and clinical implications of these findings.