Emphysema is the permanent enlargement of air spaces in the respiratory regions of the lung due to destruction of the inter-alveolar septa. The progressive coalescence of alveoli and alveolar ducts into larger airspaces leads to the disruption of normal airway wall motion and airflow rates within the pulmonary acinus. To contribute to the understanding of the individual effects of emphysema during its earliest stages, computational fluid dynamics (CFD) simulations of airflow in mathematically derived models of the pulmonary acinus were performed. The here generated computational domain consists of two generations of alveolar ducts within the pulmonary acinus, with alveolar geometries approximated as closely packed, 14-sided polygons. Physiologically realistic airflow rates and wall motions were used to study airflow patterns within subsequent generations of alveolar ducts during the inspiratory and expiratory phases of the breathing cycle. The effects of progressive emphysema on the airway wall motion and flow rates were simulated by sequentially removing all alveolar septa within each alveolar duct. Parametric studies were presented to independently assess the relative influence of progressive septal destruction of airway motion and flow rates. The results illustrate that septal destruction lowers the flow resistance through the alveolar ducts but has little influence on the mass transport of oxygen into the alveoli. Septal destruction has a net effect on the flow field by favoring the development of recirculatory flow patterns in individual alveoli.

Availability of essential species like oxygen is critical in shaping the dynamics of tumor growth. When the intracellular oxygen level falls below normal, it initiates major cascades in cellular dynamics leading to tumor cell survival. In a cellular block with cells growing away from the blood vessel, the scenario can be aggravated for the cells further inside the block. In this study, the dynamics of intracellular species inside a colony of tumor cells are investigated by varying the cell-block thickness and cell types in a microfluidic cell culture device. The oxygen transport across the cell block is modeled through diffusion, while ascorbate (AS) transport from the extracellular medium is addressed by a concentration-dependent uptake model. The extracellular and intracellular descriptions were coupled through the consumption and traffic of species from the microchannel to the cell block. Our model shows that the onset of hypoxia is possible in HeLa cell within minutes depending on the cell location, although the nutrient supply inside the channel is maintained in normoxic levels. This eventually leads to total oxygen deprivation inside the cell block in the extreme case, representing the development of a necrotic core that maintains a dynamic balance with growing cells and scarce supply. The numerical model reveals that species concentration and hypoxic response are different for HeLa and HelaS3 cells. Results also indicate that the long-term hypoxic response from a microfluidic cellular block stays within 5% of the values of a tissue with the basal layer. The hybrid model can be very useful in designing microfluidic experiments to satisfactorily predict the tissue-level response in cancer research.

The translation of many tissue engineering/regenerative medicine (TE/RM) therapies that demonstrate promise in vitro are delayed or abandoned due to reduced and inconsistent efficacy when implemented in more complex and clinically relevant preclinical in vivo models. Determining mechanistic reasons for impaired treatment efficacy is challenging after a regenerative therapy is implanted due to technical limitations in longitudinally measuring the progression of key environmental cues in vivo. The ability to acquire real-time measurements of environmental parameters of interest including strain, pressure, pH, temperature, oxygen tension, and specific biomarkers within the regenerative niche in situ would significantly enhance the information available to tissue engineers to monitor and evaluate mechanisms of functional healing or lack thereof. Continued advancements in material and fabrication technologies utilized by microelectromechanical systems (MEMSs) and the unique physical characteristics of passive magnetoelastic sensor platforms have created an opportunity to implant small, flexible, low-power sensors into preclinical in vivo models, and quantitatively measure environmental cues throughout healing. In this perspective article, we discuss the need for longitudinal measurements in TE/RM research, technical progress in MEMS and magnetoelastic approaches to implantable sensors, the potential application of implantable sensors to benefit preclinical TE/RM research, and the future directions of collaborative efforts at the intersection of these two important fields.

For newborns diagnosed with pulmonary atresia or severe pulmonary stenosis leading to insufficient pulmonary blood flow, cyanosis can be mitigated with placement of a modified Blalock–Taussig shunt (MBTS) between the innominate and pulmonary arteries. In some clinical scenarios, patients receive two systemic-to-pulmonary connections, either by leaving the patent ductus arteriosus (PDA) open or by adding an additional central shunt (CS) in conjunction with the MBTS. This practice has been motivated by the thinking that an additional source of pulmonary blood flow could beneficially increase pulmonary flow and provide the security of an alternate pathway in case of thrombosis. However, there have been clinical reports of premature shunt occlusion when more than one shunt is employed, leading to speculation that multiple shunts may in fact lead to unfavorable hemodynamics and increased mortality. In this study, we hypothesize that multiple shunts may lead to undesirable flow competition, resulting in increased residence time (RT) and elevated risk of thrombosis, as well as pulmonary overcirculation. Computational fluid dynamics-based multiscale simulations were performed to compare a range of shunt configurations and systematically quantify flow competition, pulmonary circulation, and other clinically relevant parameters. In total, 23 cases were evaluated by systematically changing the PDA/CS diameter, pulmonary vascular resistance (PVR), and MBTS position and compared by quantifying oxygen delivery (OD) to the systemic and coronary beds, wall shear stress (WSS), oscillatory shear index (OSI), WSS gradient (WSSG), and RT in the pulmonary artery (PA), and MBTS. Results showed that smaller PDA/CS diameters can lead to flow conditions consistent with increased thrombus formation due to flow competition in the PA, and larger PDA/CS diameters can lead to insufficient OD due to pulmonary hyperfusion. In the worst case scenario, it was found that multiple shunts can lead to a 160% increase in RT and a 10% decrease in OD. Based on the simulation results presented in this study, clinical outcomes for patients receiving multiple shunts should be critically investigated, as this practice appears to provide no benefit in terms of OD and may actually increase thrombotic risk.

In scaffold-based tissue engineering, sufficient oxygen and nutrient supply into cells within a scaffold is essential to increase cell viability and the proliferation rate. Generally, oxygen and nutrients reach the cells through the media by diffusion in vitro or in vivo, assuming there is no convection flow through a scaffold with small-sized pores. The scaffold diffusion rate depends mainly on the scaffold pore architecture. Thus, understanding the effect of scaffold pore architecture on the diffusion mechanism is necessary to design an efficient scaffold model. This study proposes a computational method to estimate diffusivity using the finite element analysis (FEA). This method can be applied to evaluate and analyze the effective diffusivity of a freeform fabricated 3D scaffold. The diffusion application module of commercial FEA software was used to calculate the spatial oxygen concentration gradient in a scaffold model medium. The effective diffusivities of each scaffold could be calculated from the oxygen concentration data, which revealed that the scaffold pore architecture influences its effective diffusivity. The proposed method has been verified experimentally and can be applied to design pore architectures with efficient diffusion by increasing our understanding of how the diffusion rate within a scaffold is affected by its pore architecture.

Pulmonary arterial hypertension (PAH) is a rapidly fatal disease in which mortality is typically due to right ventricular (RV) failure. An excellent predictor of mortality in PAH is proximal pulmonary artery stiffening, which is mediated by collagen accumulation in hypoxia-induced pulmonary hypertension (HPH) in mice. We sought to investigate the impact of limiting vascular and ventricular collagen accumulation on RV function and the hemodynamic coupling efficiency between the RV and pulmonary vasculature. Inbred mice were exposed to chronic hypoxia for 10 days with either no treatment (HPH) or with treatment with a proline analog that impairs collagen synthesis (CHOP-PEG; HPH + CP). Both groups were compared to control mice (CTL) exposed only to normoxia (no treatment). An admittance catheter was used to measure pressure-volume loops at baseline and during vena cava occlusion, with mice ventilated with either room air or 8% oxygen, from which pulmonary hemodynamics, RV function, and ventricular-vascular coupling efficiency (ηvvc) were calculated. Proline analog treatment limited increases in RV afterload (neither effective arterial elastance Ea nor total pulmonary vascular resistance significantly increased compared to CTL with CHOP-PEG), limited the development of pulmonary hypertension (CHOP-PEG reduced right ventricular systolic pressure by 10% compared to HPH, p < 0.05), and limited RV hypertrophy (CHOP-PEG reduced RV mass by 18% compared to HPH, p < 0.005). In an acutely hypoxic state, treatment improved RV function (CHOP-PEG increased end-systolic elastance Ees by 43%, p < 0.05) and maintained ηvvc at control, room air levels. CHOP-PEG also decreased lung collagen content by 12% measured biochemically compared to HPH (p < 0.01), with differences evident in large and small pulmonary arteries by histology. Our results demonstrate that preventing new collagen synthesis limits pulmonary hypertension development by reducing collagen accumulation in the pulmonary arteries that affect RV afterload. In particular, the proline analog limited structural and functional changes in distal pulmonary arteries in this model of early and somewhat mild pulmonary hypertension. We conclude that collagen plays an important role in small pulmonary artery remodeling and, thereby, affects RV structure and function changes induced by chronic hypoxia.

One of the primary limitations of cell therapy for myocardial infarction is the low survival of transplanted cells, with a loss of up to 80% of cells within 3 days of delivery. The aims of this study were to investigate the distribution of nutrients and oxygen in infarcted myocardium and to quantify how macromolecular transport properties might affect cell survival. Transmural myocardial infarction was created by controlled cryoablation in pigs. At 30 days post-infarction, oxygen and metabolite levels were measured in the peripheral skeletal muscle, normal myocardium, the infarct border zone, and the infarct interior. The diffusion coefficients of fluorescein or FITC-labeled dextran (0.3–70 kD) were measured in these tissues using fluorescence recovery after photobleaching. The vascular density was measured via endogenous alkaline phosphatase staining. To examine the influence of these infarct conditions on cells therapeutically used in vivo, skeletal myoblast survival and differentiation were studied in vitro under the oxygen and glucose concentrations measured in the infarct tissue. Glucose and oxygen concentrations, along with vascular density were significantly reduced in infarct when compared to the uninjured myocardium and infarct border zone, although the degree of decrease differed. The diffusivity of molecules smaller than 40 kD was significantly higher in infarct center and border zone as compared to uninjured heart. Skeletal myoblast differentiation and survival were decreased stepwise from control to hypoxia, starvation, and ischemia conditions. Although oxygen, glucose, and vascular density were significantly reduced in infarcted myocardium, the rate of macromolecular diffusion was significantly increased, suggesting that diffusive transport may not be inhibited in infarct tissue, and thus the supply of nutrients to transplanted cells may be possible. in vitro studies mimicking infarct conditions suggest that increasing nutrients available to transplanted cells may significantly increase their ability to survive in infarct.

Multiphysics modeling is an emerging approach in cellular bioengineering research, used for simulating complex biophysical interactions and their effects on cell viability and function. Our goal in the present study was to integrate cell-specific finite element modeling—which we have developed in previous research to simulate deformation of individual cells subjected to external loading—with oxygen transport in the deformed cells at normoxic and hypoxic environments. We specifically studied individual and combined effects of substrate stretch levels, O 2 concentration in the culture media, and temperature of the culture media on intracellular O 2 levels in cultured myoblasts, in models of two individual cells. We found that (i) O 2 transport became faster with the increasing levels of substrate stretch (ranging from 0 to 24%), and (ii) the effect of a 3 °C temperature drop on slowing down the O 2 transport was milder with respect to the effect that strains had. The changes in cell geometry due to externally applied deformations could, hence, theoretically affect cell respiration, which should be a consideration in cellular mechanics experiments.

An idealized systemic-to-pulmonary shunt anatomy is parameterized and coupled to a closed loop, lumped parameter network (LPN) in a multidomain model of the Norwood surgical anatomy. The LPN approach is essential for obtaining information on global changes in cardiac output and oxygen delivery resulting from changes in local geometry and physiology. The LPN is fully coupled to a custom 3D finite element solver using a semi-implicit approach to model the heart and downstream circulation. This closed loop multidomain model is then integrated with a fully automated derivative-free optimization algorithm to obtain optimal shunt geometries with variable parameters of shunt diameter, anastomosis location, and angles. Three objective functions: (1) systemic; (2) coronary; and (3) combined systemic and coronary oxygen deliveries are maximized. Results show that a smaller shunt diameter with a distal shunt-brachiocephalic anastomosis is optimal for systemic oxygen delivery, whereas a more proximal anastomosis is optimal for coronary oxygen delivery and a shunt between these two anatomies is optimal for both systemic and coronary oxygen deliveries. Results are used to quantify the origin of blood flow going through the shunt and its relationship with shunt geometry. Results show that coronary artery flow is directly related to shunt position.

The aim of this study was to maximize oxygen diffusion within a three-dimensional scaffold in order to improve cell viability and proliferation. To evaluate the effect of pore architecture on oxygen diffusion, we designed a regular channel shape with uniform diameter, referred to as cylinder shaped, and a new channel shape with a channel diameter gradient, referred to as cone shaped. A numerical analysis predicted higher oxygen concentration in the cone-shaped channels than in the cylinder-shaped channels, throughout the scaffold. To confirm these numerical results, we examined cell proliferation and viability in 2D constructs and 3D scaffolds. Cell culture experiments revealed that cell proliferation and viability were superior in the constructs and scaffolds with cone-shaped channels.

Oxygen is essential to maintaining normal brain function. A large body of evidence suggests that the partial pressure of oxygen ( p O 2 ) in brain tissue is physiologically maintained within a narrow range in accordance with region-specific brain activity. Since the transportation of oxygen in the brain tissue is mainly driven by a diffusion process caused by a concentration gradient of oxygen from blood to cells, the spatial organization of the vascular system, in which the oxygen content is higher than in tissue, is a key factor for maintaining effective transportation. In addition, a local mechanism that controls energy demand and blood flow supply plays a critical role in moment-to-moment adjustment of tissue p O 2 in response to dynamically varying brain activity. In this review, we discuss the spatiotemporal structures of brain tissue oxygen transport in relation to local brain activity based on recent reports of tissue p O 2 measurements with polarographic oxygen microsensors in combination with simultaneous recordings of neural activity and local cerebral blood flow in anesthetized animal models. Although a physiological mechanism of oxygen level sensing and control of oxygen transport remains largely unknown, theoretical models of oxygen transport are a powerful tool for better understanding the short-term and long-term effects of local changes in oxygen demand and supply. Finally, emerging new techniques for three-dimensional imaging of the spatiotemporal dynamics of p O 2 map may enable us to provide a whole picture of how the physiological system controls the balance between demand and supply of oxygen during both normal and pathological brain activity.

The intervertebral disk (IVD) is the largest avascular structure in the human body. Transport of small molecules in IVD is mainly through diffusion from the endplates and the peripheral blood vessels surrounding IVD. Studies have investigated the structure, chemical components, and water content in IVD, but to our knowledge no study has investigated the effect of mechanical loading on oxygen transport in IVD. The objective of this study was to determine the strain-dependent behavior of oxygen diffusivity in IVD tissue. A one-dimensional steady-state diffusion experiment was designed and performed to determine the oxygen diffusivity in bovine annulus fibrosus (AF). The oxygen diffusivity was calculated using equation derived from Fick’s law. A total of 20 AF specimens ( d = 6 mm , h ∼ 0.5 mm ) from bovine coccygeal IVD were used to determine oxygen diffusivity at three levels of compressive strain. The average oxygen diffusivity ( mean ± SD ) of bovine AF in the axial direction was 1.43 ± 0.242 × 10 − 5 cm 2 / s ( n = 20 ) at 4.68 ± 1.67 % compressive strain level, 1.05 ± 0.282 × 10 − 5 cm 2 / s ( n = 20 ) at 14.2 ± 1.50 % strain level, and 7.71 ± 1.63 × 10 − 6 cm 2 / s ( n = 20 ) at 23.7 ± 1.34 % strain level. There was a significant decrease in oxygen diffusivity with increasing level of compressive strain (ANOVA, p < 0.05 ). Oxygen diffusivity of bovine AF in the axial direction has been determined. The mechanical loading has a significant effect on oxygen transport in IVD tissues. This study is important in understanding nutritional transport in IVD tissues and related disk degeneration.

Various cell populations have been shown to respond to hydrostatic pressure; however, many of the culture systems suffer from shortcomings in design or methodology. Of particular interest to us is the potential role of pressure and other environmental factors in modulating stem cell behavior in intervertebral disk repair. A system was developed for the growth of cells and tissues under intermittent hydrostatic pressure. The system was validated with NIH 3T3 fibroblasts for sterilizability and cytotoxicity. Further experiments were conducted with canine mesenchymal stem cells under various levels of pressure, oxygen, glucose, and conditioned medium. The culture system showed no cytotoxicity and was able to demonstrate that the proliferation and metabolism of mesenchymal stem cells are sensitive to medium glucose and oxygen concentration and hydrostatic pressure. The cells exposed to hydrostatic pressure differed in their morphology from nonexposed cells. The system is capable of supporting long-term cell culture and examining the role of mechanical and environmental stimulation in vivo.

The influence of time-dependent flows on oxygen transport from hollow fibers was computationally and experimentally investigated. The fluid average pressure drop, a measure of resistance, and the work required by the heart to drive fluid past the hollow fibers were also computationally explored. This study has particular relevance to the development of an artificial lung, which is perfused by blood leaving the right ventricle and in some cases passing through a compliance chamber before entering the device. Computational studies modeled the fiber bundle using cylindrical fiber arrays arranged in in-line and staggered rectangular configurations. The flow leaving the compliance chamber was modeled as dampened pulsatile and consisted of a sinusoidal perturbation superimposed on a steady flow. The right ventricular flow was modeled to depict the period of rapid flow acceleration and then deceleration during systole followed by zero flow during diastole. Experimental studies examined oxygen transfer across a fiber bundle with either steady, dampened pulsatile, or right ventricular flow. It was observed that the dampened pulsatile flow yielded similar oxygen transport efficiency to the steady flow, while the right ventricular flow resulted in smaller oxygen transport efficiency, with the decrease increasing with Re. Both computations and experiments yielded qualitatively similar results. In the computational modeling, the average pressure drop was similar for steady and dampened pulsatile flows and larger for right ventricular flow while the pump work required of the heart was greatest for right ventricular flow followed by dampened pulsatile flow and then steady flow. In conclusion, dampening the artificial lung inlet flow would be expected to maximize oxygen transport, minimize work, and thus improve performance.

Flow analysis at microvascular bifurcation after partial replacement of red blood cell (RBC) with liposome-encapsulated hemoglobin (LEH) was performed using the lattice Boltzmann method. A two-dimensional symmetric Y bifurcation model with a parent vessel diameter of 20 μ m and daughter branch diameters of 20 μ m was considered, and the distributions of the RBC, LEH, and oxygen fluxes were calculated. When only RBCs flow into the daughter branches with unevenly distributed flows, plasma separation occurred and the RBC flow to the lower-flow branch was disproportionately decreased. On the other hand, when half of RBC are replaced by LEH, the biasing of RBC flow was enhanced whereas LEH flowed favorably into the lower-flow branch, because many LEH within the parent vessel are suspended in the plasma layer, where no RBCs exist. Consequently, the branched oxygen fluxes became nearly proportional to flows. These results indicate that LEH facilitates oxygen supply to branches that are inaccessible to RBCs.

The effect of blood viscosity on oxygen transport in a stenosed coronary artery during the postangioplasty scenario is studied. In addition to incorporating varying blood viscosity using different hematocrit (Hct) concentrations, oxygen consumption by the avascular wall and its supply from vasa vasorum, nonlinear oxygen binding capacity of the hemoglobin, and basal to hyperemic flow rate changes are included in the calculation of oxygen transport in both the lumen and the avascular wall. The results of this study show that oxygen transport in the postangioplasty residual stenosed artery is affected by non-Newtonian shear-thinning property of the blood viscosity having variable Hct concentration. As Hct increases from 25% to 65%, the diminished recirculation zone for the increased Hct causes the commencement of p O 2 decrease to shift radially outward by ∼ 20 % from the center of the artery for the basal flow, but by ∼ 10 % for the hyperemic flow at the end of the diverging section. Oxygen concentration increases from a minimum value at the core of the recirculation zone to over 90 mm Hg before the lumen-wall interface at the diverging section for the hyperemic flow, which is attributed to increased shear rate and thinner lumen boundary layer for the hyperemic flow, and below 90 mm Hg for the basal flow. As Hct increases from 25% to 65%, the average of p O 2 , min beyond the diverging section drops by ∼ 25 % for the basal flow, whereas it increases by ∼ 15 % for the hyperemic flow. Thus, current results with the moderate stenosed artery indicate that reducing Hct might be favorable in terms of increasing O 2 flux and p O 2 , min , in the medial region of the wall for the basal flow, while higher Hct is advantageous for the hyperemic flow beyond the diverging section. The results of this study not only provide significant details of oxygen transport under varying pathophysiologic blood conditions such as unusually high blood viscosity and flow rate, but might also be extended to offer implications for drug therapy related to blood-thinning medication and for blood transfusion and hemorrhage.

The fundamental study of blood flow past a circular cylinder filled with an oxygen source is investigated as a building block for an artificial lung. The Casson constitutive equation is used to describe the shear-thinning and yield stress properties of blood. The presence of hemoglobin is also considered. Far from the cylinder, a pulsatile blood flow in the x direction is prescribed, represented by a time periodic (sinusoidal) component superimposed on a steady velocity. The dimensionless parameters of interest for the characterization of the flow and transport are the steady Reynolds number (Re), Womersley parameter ( α ) , pulsation amplitude ( A ) , and the Schmidt number (Sc). The Hill equation is used to describe the saturation curve of hemoglobin with oxygen. Two different feed-gas mixtures were considered: pure O 2 and air. The flow and concentration fields were computed for Re=5, 10, and 40, 0 ≤ A ≤ 0.75 , α = 0.25 , 0.4, and Schmidt number, Sc=1000. The Casson fluid properties result in reduced recirculations (when present) downstream of the cylinder as compared to a Newtonian fluid. These vortices oscillate in size and strength as A and α are varied. Hemoglobin enhances mass transport and is especially important for an air feed which is dominated by oxyhemoglobin dispersion near the cylinder. For a pure O 2 feed, oxygen transport in the plasma dominates near the cylinder. Maximum oxygen transport is achieved by operating near steady flow (small A ) for both feed-gas mixtures. The time averaged Sherwood number, Sh ̿ , is found to be largely influenced by the steady Reynolds number, increasing as Re increases and decreasing with A . Little change is observed with varying α for the ranges investigated. The effect of pulsatility on Sh ̿ is greater at larger Re. Increasing Re aids transport, but yields a higher cylinder drag force and shear stresses on the cylinder surface which are potentially undesirable.

The pulsatile blood flow and gas transport of oxygen and carbon dioxide through a cylindrical array of microfibers are numerically simulated. Blood is modeled as a homogeneous Casson fluid, and hemoglobin molecules in blood are assumed to be in local equilibrium with oxygen and carbon dioxide. It is shown that flow pulsatility enhances gas transport and the amount of gas exchange is sensitive to the blood flow field across the fibers. The steady Sherwood number dependence on Reynolds number was shown to have a linear relation consistent with experimental findings. For most cases, an enhancement in gas transport is accompanied with an increase in flow resistance. Maximum local shear stress is provided as a possible indicator of thrombosis, and the computed shear stress is shown to be below the threshold value for thrombosis formation for all cases evaluated.

The Circle of Willis (CoW) is a ringlike structure of blood vessels found at the base of the brain. Its main function is to distribute oxygen-rich arterial blood to the cerebral mass. In a previous study, a one-dimensional (1D) model of the CoW was created to simulate a series of possible clinical scenarios such as occlusions in afferent arteries, absent or stringlike circulus vessels, or arterial infarctions ( Moorhead et al., 2004, Comput. Methods Biomech. Biomed. Eng., 7(3), pp. 121–130 ). The model captured cerebral haemodynamic autoregulation by using a proportional-integral-derivative (PID) controller to modify efferent artery resistances. Although some good results and correlations were achieved, the model was too simple to capture all the transient dynamics of autoregulation. Hence a more physiologically accurate model has been created that additionally includes the oxygen dynamics that drive the autoregulatory response. Results very closely match accepted physiological response and limited clinical data. In addition, a set of boundary conditions and geometry is presented for which the autoregulated system cannot provide sufficient perfusion, representing a condition with increased risk of stroke and highlighting the importance of modeling the haemodynamics of the CoW. The system model created is computationally simple so it can be used to identify at-risk cerebral arterial geometries and conditions prior to surgery or other clinical procedures.

The pulsatile flow and gas transport of a Newtonian passive fluid across an array of cylindrical microfibers are numerically investigated. It is related to an implantable, artificial lung where the blood flow is driven by the right heart. The fibers are modeled as either squared or staggered arrays. The pulsatile flow inputs considered in this study are a steady flow with a sinusoidal perturbation and a cardiac flow. The aims of this study are twofold: identifying favorable array geometry/spacing and system conditions that enhance gas transport; and providing pressure drop data that indicate the degree of flow resistance or the demand on the right heart in driving the flow through the fiber bundle. The results show that pulsatile flow improves the gas transfer to the fluid compared to steady flow. The degree of enhancement is found to be significant when the oscillation frequency is large, when the void fraction of the fiber bundle is decreased, and when the Reynolds number is increased; the use of a cardiac flow input can also improve gas transfer. In terms of array geometry, the staggered array gives both a better gas transfer per fiber (for relatively large void fraction) and a smaller pressure drop (for all cases). For most cases shown, an increase in gas transfer is accompanied by a higher pressure drop required to power the flow through the device.