Gastric flow and mixing studied using computer simulation

The fed human stomach displays regular peristaltic contraction waves that originate in the proximal antrum and propagate to the pylorus. High-resolution concurrent manometry and magnetic resonance imaging (MRI) studies of the stomach suggest a primary function of antral contraction wave (ACW) activity unrelated to gastric emptying. Detailed evaluation is difficult, however, in vivo. Here we analyse the role of ACW activity on intragastric fluid motions, pressure, and mixing with computer simulation. A two-dimensional computer model of the stomach was developed with the 'lattice-Boltzmann' numerical method from the laws of physics, and stomach geometry modelled from MRI. Time changes in gastric volume were specified to match global physiological rates of nutrient liquid emptying. The simulations predicted two basic fluid motions: retrograde 'jets' through ACWs, and circulatory flow between ACWs, both of which contribute to mixing. A well-defined 'zone of mixing', confined to the antrum, was created by the ACWs, with mixing motions enhanced by multiple and narrower ACWs. The simulations also predicted contraction-induced peristaltic pressure waves in the distal antrum consistent with manometric measurements, but with a much lower pressure amplitude than manometric data, indicating that manometric pressure amplitudes reflect direct contact of the catheter with the gastric wall. We conclude that the ACWs are central to gastric mixing, and may also play an indirect role in gastric emptying through local alterations in common cavity pressure.     

Solid Loss of Carrots During Simulated Gastric Digestion

The knowledge of solid loss kinetics of foods during digestion is crucial for understanding the factors that constrain the release of nutrients from the food matrix and their fate of digestion. The objective of this study was to investigate the solid loss of carrots during simulated gastric digestion as affected by pH, temperature, viscosity of gastric fluids, mechanical force present in stomach, and cooking. Cylindrical carrot samples were tested by static soaking method and using a model stomach system. The weight retention, moisture, and loss of dry mass were determined. The results indicated that acid hydrolysis is critical for an efficient mass transfer and carrot digestion. Internal resistance rather than external resistance is dominant in the transfer of soluble solids from carrot to gastric fluid. Increase in viscosity of gastric fluid by adding 0.5% gum (w/w) significantly increased the external resistance and decreased mass transfer rate of carrots in static soaking. When mechanical force was not present, 61% of the solids in the raw carrot samples were released into gastric fluid after 4 h of static soaking in simulated gastric juice. Mechanical force significantly increased solid loss by causing surface erosion. Boiling increased the disintegration of carrot during digestion that may favor the loss of solids meanwhile reducing the amount of solids available for loss in gastric juice. Weibull function was successfully used to describe the solid loss of carrot during simulated digestion. The effective diffusion coefficients of solids were calculated using the Fick’s second law of diffusion for an infinite cylinder, which are between 0.75 × 10⁻¹¹ and 8.72 × 10⁻¹¹ m²/s, depending on the pH of the gastric fluid.     

The effect of compliance on wall shear in casts of a human aortic bifurcation

Rigid and compliant casts of a human aortic bifurcation were subjected to physiologically realistic pulsatile fluid flows. At a number of sites near the wall in the approximate median plane of the bifurcation of these models, fluid velocity was measured with a laser Doppler velocimeter, and wall motion (in the case of the compliant cast) was determined with a Reticon linescan camera. The velocity and wall motion data were combined to estimate the instantaneous shear rates at the cast wall. Analysis showed that at the outer walls the cast compliance reduced shear rates, while at the walls of the flow divider the shear rate was increased.     

Tortuosity triggers platelet activation and thrombus formation in microvessels

Tortuous blood vessels are often seen in humans in association with thrombosis, atherosclerosis, hypertension, and aging. Vessel tortuosity can cause high fluid shear stress, likely promoting thrombosis. However, the underlying physical mechanisms and microscale processes are poorly understood. Accordingly, the objectives of this study were to develop and use a new computational approach to determine the effects of venule tortuosity and fluid velocity on thrombus initiation. The transport, collision, shear-induced activation, and receptor-ligand adhesion of individual platelets in thrombus formation were simulated using discrete element method. The shear-induced activation model assumed that a platelet became activated if it experienced a shear stress above a relative critical shear stress or if it contacted an activated platelet. Venules of various levels of tortuosity were simulated for a mean flow velocity of 0.10?cm s(-1), and a tortuous arteriole was simulated for a mean velocity of 0.47?cm s(-1). Our results showed that thrombus was initiated at inner walls in curved regions due to platelet activation in agreement with experimental studies. Increased venule tortuosity modified fluid flow to hasten thrombus initiation. Compared to the same sized venule, flow in the arteriole generated a higher amount of mural thrombi and platelet activation rate. The results suggest that the extent of tortuosity is an important factor in thrombus initiation in microvessels.     

Numerical simulation of blood flow in the left ventricle and aortic sinus using magnetic resonance imaging and computational fluid dynamics

Understanding cardiac blood flow patterns has many applications in analysing haemodynamics and for the clinical assessment of heart function. In this study, numerical simulations of blood flow in a patient-specific anatomical model of the left ventricle (LV) and the aortic sinus are presented. The realistic 3D geometry of both LV and aortic sinus is extracted from the processing of magnetic resonance imaging (MRI). Furthermore, motion of inner walls of LV and aortic sinus is obtained from cine-MR image analysis and is used as a constraint to a numerical computational fluid dynamics (CFD) model based on the moving boundary approach. Arbitrary Lagrangian–Eulerian finite element method formulation is used for the numerical solution of the transient dynamic equations of the fluid domain. Simulation results include detailed flow characteristics such as velocity, pressure and wall shear stress for the whole domain. The aortic outflow is compared with data obtained by phase-contrast MRI. Good agreement was found between simulation results and these measurements.     

Shear stress gradient over endothelial cells in a curved microchannel system.

Our purpose was to test a scale model of the microcirculation by measuring the shear forces to which endothelial cells were exposed, and comparing this to computer simulations. In vitro experiments were performed to measure the 2-dimensional projected velocity profile along endothelial cell lined microchannels (D-shaped, 10-30 microns radius, n = 15), or in microchannels without endothelial cells (n = 18). Microchannels were perfused with fluorescently labeled microspheres (0.5 micron dia., < 1%) suspended in cell culture media. The velocity of individual microspheres was obtained off-line (videorecording), using an interactive software program; velocity was determined as the distance traveled in one video field (1/60 s). Mass balance was verified in the microchannels by comparing the microsphere velocities to the perfusion pump rate. In confluent endothelial cell lined microchannels, a velocity profile was obtained as microspheres passed an endothelial cell nucleus (identified by fluorescent dye), and again, for a paired region 100 microns away without nuclei (cytoplasm region). The velocity profile was significantly shifted and sharpened by the endothelial cell nucleus, as anticipated. Over the nucleus, data are consistent with a normal sized nucleus extending into the lumen, further confirming that this scale model can be used to determine the wall shear stress to which endothelial cells are exposed. Using the experimental bulk phase fluid parameters as boundary conditions, we used computational fluid dynamics (CFD) to predict the expected wall shear stress gradient along an endothelial cell lined D-shaped tube. The wall shear stress gradient over the nucleus was 2-fold greater in the radial versus axial directions, and was sensitive to lateral versus midline positioned nuclei.     

On the kinematics and efficiency of advective mixing during gastric digestion – A numerical analysis

The mixing performance of gastric contents during digestion is expected to have a major role on the rate and final bioavailability of nutrients within the body. The aim of this study was to characterize the ability of the human stomach to advect gastric contents with different rheological properties. The flow behavior of two Newtonian fluids (10⁻³ Pa s, 1 Pa s) and a pseudoplastic solution (K=0.223 Pa s^0.59) during gastric digestion were numerically characterized within a simplified 3D model of the stomach geometry and motility during the process (ANSYS-FLUENT). The advective performances of each of these gastric flows were determined by analyzing the spatial distribution and temporal history of their stretching abilities (Lagrangian analysis). Results illustrate the limited influence that large retropulsive and vortex structures have on the overall dynamics of gastric flows. Even within the distal region, more than 50% of the flow experienced velocity and shear values lower than 10% of their respective maximums. While chaotic, gastric advection was always relatively poor (with Lyapunov exponents an order of magnitude lower than those of a laminar stirred tank). Contrary to expectations, gastric rheology had only a minor role on the advective properties of the flow (particularly within the distal region). As viscosity increased above 1 St, the role of fluid viscosity became largely negligible. By characterizing the fluid dynamic and mixing conditions that develop during digestion, this work will inform the design of novel in vitro systems of enhanced biomechanical performance and facilitate a more accurate diagnosis of gastric digestion processes.     

An Analytic Study on the Effect of Alginate on the Velocity Profiles of Blood in Rectangular Microchannels Using Microparticle Image Velocimetry

It is desired to understand the effect of alginic acid sodium salt from brown algae (alginate) as a viscosity modifier on the behavior of blood in vitro using a micro-particle image velocimetry (µPIV) system. The effect of alginate on the shape of the velocity profile, the flow rate and the maximum velocity achieved in rectangular microchannels channels are measured. The channels were constructed of polydimethylsiloxane (PDMS), a biocompatible silicone. Porcine blood cells suspended in saline was used as the working fluid at twenty percent hematocrit (H = 20). While alginate was only found to have minimal effect on the maximum velocity and the flow rate achieved, it was found to significantly affect the shear rate at the wall by between eight to a hundred percent.     

Pulsatile flow in a coronary artery using multiphase kinetic theory

Pulsatile flow in a model of a right coronary artery (RCA) was previously modeled as a single-phase fluid and as a two-phase fluid using experimental rheological data for blood as a function of hematocrit and shear rate. Here we present a multiphase kinetic theory model which has been shown to compute correctly the viscosity of red blood cells (RBCs) and their migration away from vessel walls: the Fahraeus–Lindqvist effect. The computed RBC viscosity decreases with shear rate and vessel size, consistent with measurements. The pulsatile computations were performed using a typical cardiac waveform until a limit cycle was well established. The RBC volume fractions, shear stresses, shear stress gradients, granular temperatures, viscosities, and phase velocities varied with time and position during each cardiac cycle. Steady-state computations were also performed and were found to compare well with time-averaged transient results. The wall shear stress and wall shear stress gradients (both spatial and temporal) were found to be highest on the inside area of maximum curvature. Potential atherosclerosis sites are identified using these computational results.     

A lubrication analysis of pharyngeal peristalsis: Application to flavour release

After eating a liquid or a semi-liquid food product, a thin film responsible for the dynamic profile of aroma release coats the pharyngeal mucosa. The aim of this article was to analyse the fluid mechanics of pharyngeal peristalsis and to develop a simple biomechanical model in order to understand the role of saliva and food bolus viscosity on the coating of pharyngeal mucosa. We began by analysing the physiology and the biomechanics of swallowing in order to determine relevant model assumptions. This analysis of the literature clarified the types of mechanical solicitations applied on the food bolus. Moreover, we showed that the pharyngeal peristalsis in the most occluded region is equivalent to a forward roll coating process, the originality of which is lubrication by a film of saliva. A model based on the lubrication theory for Newtonian liquids was developed in dimensionless form. The parametric study showed the strong influence of relative saliva thickness on the food bolus coating. A specific experimental device was designed that confirms the model predictions. Two sets of conditions that depend on the relative thickness of saliva were distinguished. The first is characterised by a relatively thin film of saliva: food bolus viscosity has a strong impact on mucosa coating. These phenomena are well represented by the model developed here. The second is obtained when the saliva film is relatively thick: hydrodynamic mixing with saliva, interdiffusion or instabilities may govern mucosa coating. Finally, these results were extrapolated to determine the influence of food bolus viscosity on the dynamic profile of flavour release according to physiological parameters.     

Motion and deformation of a red blood cell in a shear flow: a two-dimensional simulation of the wall effect

The motion and deformation of a single red blood cell in a simple shear flow between two parallel walls is studied theoretically. A two-dimensional deformable microcapsule is adopted as a model for the cell, which has a thin moving membrane, like a tank-tread, around the interior and is deformed into an elliptical shape with a constant area. Applying the finite element method to the Stokes equations, the tank-tread motion and deformation is determined in a stationary motion, under fluid dynamic interaction between the cell and the walls. It is shown that the motion and deformation of the microcapsule crucially depends on the channel width between the two walls. As the width decreases, the microcapsule is more elongated and the frequency of tank-tread motion decreases at a constant shear rate. In addition, the angle of inclination decreases at the low range of the viscosity ratio of internal to external fluids and increases at the high range. The results obtained are compared with experimental observations and applied to the behavior of cells under mutual interaction.     

Computational study of the effect of geometric and flow parameters on the steady flow field at the rabbit aorto-celiac bifurcation

Arterial hemodynamic forces may play a role in the localization of early atherosclerotic lesions. We have been developing numerical techniques based on overset or "Chimera" type formulations to solve the Navier-Stokes equations in complex geometries simulating arterial bifurcations. This paper presents three-dimensional steady flow computations in a model of the rabbit aorto-celiac bifurcation. The computational methods were validated by comparing the numerical results to previously-obtained flow visualization data. Once validated, the numerical algorithms were used to investigate the sensitivity of the computed flow field and resulting wall shear stress distribution to various geometric and hemodynamic parameters. The results demonstrated that a decrease in the extent of aortic taper downstream of the celiac artery induced looping fluid motion along the lateral walls of the aorta and shifted the peak wall shear stress from downstream of the celiac artery to upstream. Increasing the flow Reynolds number led to a sharp increase in spatial gradients of wall shear stress. The flow field was highly sensitive to the flow division ratio, i.e., the fraction of total flow rate that enters the celiac artery, with larger values of this ratio leading to the occurrence of flow separation along the dorsal wall of the aorta. Finally, skewness of the inlet velocity profile had a profound impact on the wall shear stress distribution near the celiac artery. While not physiological due to the assumption of steady flow, these results provide valuable insight into the fluid physics at geometries simulating arterial bifurcations.     

Flow-pressure drop measurement and calculation in a tapered femoral artery of a dog

This study describes the in vivo measurement of pressure drop and flow during the cardiac cycle in the femoral artery of a dog, and the computer simulation of the experiment based on the use of the measured flow, vessel dimensions and blood viscosity. In view of the experimental uncertainty in obtaining the accurate velocity profile at the wall region, the velocity pulse at the center was measured and numerical calculations were performed for the center Une instantaneous velocity and within the two limits of spatial distribution of inlet flow conditions: uniform and parabolic. Temporal and spatial variations of flow parameters, i.e., velocity profile, shear rate, non-Newtonian viscosity, wall shear stress, and pressure drop were calculated. There existed both positive and negative shear rates during a pulse cycle, i.e., the arterial wall experiences zero shear three times during a cardiac cycle. For the parabolic inlet condition, the taper of the artery not only increased the magnitude of the positive and negative shear rates, but caused a steep gradient in shear rate, a phenomenon which in turn affects wall shear stress and pressure. In contrast, for the uniform inlet condition, the flow through the tapered artery was predominantly the developing type, which resulted in reduction in magnitude of wall shear rate along the axial direction.     

Gastric secretion in dogs following three day fasting and resumed feeding

The investigation into the influence of a three-day starvation on the gastric secretion in dogs with Pavlov pouches stimulated by meat, histamine and pentagastrin, was carried out. A 72-hour starvation did not change the summary volume of the gastric juice, debit of the gastric acid, and quantity of pepsin. At the same time the starvation decreased the average rate of gastric juice secretion, gastric acid and pepsin secretion in response to histamine and decreased the pepsin secretion in response to pentagastrin. In this way re-feeding enhanced the average rate of gastric juice secretion and gastric acid secretion on 3-day and pepsin on 5-day in response to meat. The average rate of gastric juice secretion increased on the 5-day after refeeding in response to histamine and the average rate of gastric juice, gastric acid and pepsin secretion in response to pentagastrin.     

Peristaltic transport in a channel with a porous peripheral layer: model of a flow in gastrointestinal tract

Peristaltic transport in a two dimensional channel, filled with a porous medium in the peripheral region and a Newtonian fluid in the core region, is studied under the assumptions of long wavelength and low Reynolds number. The fluid flow is investigated in the waveframe of reference moving with the velocity of the peristaltic wave. Brinkman extended Darcy equation is utilized to model the flow in the porous layer. The interface is determined as a part of the solution using the conservation of mass in both the porous and fluid regions independently. A shear-stress jump boundary condition is used at the interface. The physical quantities of importance in peristaltic transport like pumping, trapping, reflux and axial velocity are discussed for various parameters of interest governing the flow like Darcy number, porosity, permeability, effective viscosity etc. It is observed that the peristalsis works as a pump against greater pressure in two-layered model with a porous medium compared with a viscous fluid in the peripheral layer. Increasing Darcy number Da decreases the pumping and increasing shear stress jump constant beta results in increasing the pumping. The limits on the time averaged flux Q for trapping in the core layer are obtained. The discussion on pumping, trapping and reflux may be helpful in understanding some of the fluid dynamic aspects of the transport of chyme in gastrointestinal tract.     

Computational and experimental investigation of flow and fluid mixing in the roller bottle bioreactor

The fully three-dimensional velocity field in a roller bottle bioreactor is simulated for two systems (creeping flow and inertial flow conditions) using a control volume-finite element method, and validated experimentally using particle imaging velocimetry. The velocity fields and flow patterns are described in detail using velocity contour plots and tracer particle pathline computations. Bulk fluid mixing in the roller bottle is then examined using a computational fluid tracer program and flow visualization experiments. It is shown that the velocity fields and flow patterns are substantially different for each of these flow cases. For creeping flow conditions the flow streamlines consist of symmetric, closed three-dimensional loops; and for inertial flow conditions, streamlines consist of asymmetric toroidal surfaces. Fluid tracers remain trapped on these streamlines and are unable to contact other regions of the flow domain. As a result, fluid mixing is greatly hindered, especially in the axial direction. The lack of efficient axial mixing is verified computationally and experimentally. Such mixing limitations, however, are readily overcome by introducing a small-amplitude vertical rocking motion that disrupts both symmetry and recirculation, leading to much faster and complete axial mixing. The frequency of such motion is shown to have a significant effect on mixing rate, which is a critical parameter in the overall performance of roller bottles.     

Pressure flow patterns in a cylinder with reabsorbing walls

The problem of fluid motion in renal tubules, in contrast to ordinary flow through cylinders with impermeable walls, is complicated by the existence of radial velocities generated by reabsorption processes. As a first approach to this problem, the Navier Stokes equations for axially symmetric, slow flow in an infinite cylinder whose walls reabsorb fluid are integrated. If the rate of reabsorption is constant, the solutions resemble the conventional Poiseuille flow, i.e., the longitudinal velocity profile is parabolic. In addition the drop in mean pressure is proportional to the mean axial flow, the length of tube between reference points, and inversely proportional to the fourth power of the radius. If the rate of reabsorption is a linear function of the distance from the origin, the presence of an additive term alters these relations. If, for example, the gradient in reabsorption is positive, the axial velocity profile tends to flatten and when the gradient is sufficiently large, the maximum velocity moves from the center of the stream toward the periphery, leaving a relative minimum at the center. In passing from the center of the tube to the walls, the radial velocity passes through a miximum, regardless of the reabsorption properties of the wall.     

Computational fluid dynamics modeling of the paddle dissolution apparatus: Agitation rate,mixing patterns,and fluid velocities

The purpose of this research was to further investigate the hydrodynamics of the United States Pharmacopeia (USP) paddle dissolution apparatus using a previously generated computational fluid dynamics (CFD) model. The influence of paddle rotational speed on the hydrodynamics in the dissolution vessel was simulated. The maximum velocity magnitude for axial and tangential velocities at different locations in the vessel was found to increase linearly with the paddle rotational speed. Path-lines of fluid mixing, which were examined from a central region at the base of the vessel, did not reveal a region of poor mixing between the upper cylin-drical and lower hemispherical volumes, as previously speculated. Considerable differences in the resulting flow patterns were observed for paddle rotational speeds between 25 and 150 rpm. The approximate time required to achieve complete mixing varied between 2 to 5 seconds at 150 rpm and 40 to 60 seconds at 25 rpm, although complete mixing was achievable for each speed examined. An analysis of CFD-generated velocities above the top surface of a cylindrical compact positioned at the base of the vessel, below the center of the rotating paddle, revealed that the fluid in this region was undergoing solid body rotation. An examination of the velocity boundary layers adjacent to the curved surface of the compact revealed large peaks in the shear rates for a region within∼3 mm from the base of the compact, consistent with a ‘grooving’ effect, which had been previously seen on the surface of compacts following dissolution, associated with a higher dissolution rate in this region.     

A dynamical study of the mechanical stimuli and tissue differentiation within a CaP scaffold based on micro-CT finite element models

The control of the mechanical stimuli transmitted to the cells is critical for the design of functional scaffolds for tissue engineering. The objective of this study was to investigate the dynamics of the mechanical stimuli transmitted to the cells during tissue differentiation in an irregular morphology scaffold under compressive load and perfusion flow. A calcium phosphate-based glass porous scaffold was used. The solid phase and the fluid flow within the pores were modeled as linear elastic solid material and Newtonian fluid, respectively. In the fluid model, different levels of viscosity were used to simulate tissue differentiation. Compressive strain of 0.5% and fluid flow with constant inlet velocity of 10 μm/s or constant inlet pressure of 3 Pa were applied. Octahedral shear strain and fluid shear stress were used as mechano-regulatory stimuli. For constant inlet velocity, stimuli equivalent to bone were predicted in 80% of pore volume for the case of low tissue viscosity. For the cases of high viscosity, fluctuations between stimuli equivalent to tissue formation and cell death were predicted due to the increase in the fluid shear stress when tissue started to fill pores. When constant pressure was applied, stimuli equivalent to bone were predicted in 62% of pore volume when low tissue viscosity was used and 42% when high tissue viscosity was used. This study predicted critical variations of fluid shear stress when cells differentiated. If these variations are not controlled in vitro, they can impede the formation of new matured tissue.     

Three-dimensional cerebrospinal fluid flow within the human ventricular system

Cerebrospinal fluid (CSF) is a Newtonian fluid and can, therefore, be modelled using computational fluid dynamics (CFD). Previous modelling of the CSF has been limited to simplified geometric models. This work describes a geometrically accurate three dimensional (3D) computational model of the human ventricular system (HVS) constructed from magnetic resonance images (MRI) of the human brain. It is an accurate and full representation of the HVS and includes appropriately positioned CSF production and drainage locations. It was used to investigate the pulsatile motion of CSF within the human brain. During this investigation CSF flow rate was set at a constant 500 ml/day, to mimic real life secretion of CSF into the system, and a pulsing velocity profile was added to the inlets to incorporate the effect of cardiac pulsations on the choroid plexus and their subsequent influence on CSF motion in the HVS. Boundary conditions for the CSF exits from the ventricles (foramina of Magendie and Lushka) were found using a “nesting” approach, in which a simplified model of the entire central nervous system (CNS) was used to examine the effects of the CSF surrounding the ventricular system (VS). This model provided time varying pressure data for the exits from the VS nested within it. The fastest flow was found in the cerebral aqueduct, where a maximum velocity of 11.38 mm/s was observed over five cycles. The maximum Reynolds number recorded during the simulation was 15 with an average Reynolds number of the order of 0.39, indicating that CSF motion is creeping flow in most of the computational domain and consequently will follow the geometry of the model. CSF pressure also varies with geometry with a maximum pressure drop of 1.14 Pa occurring through the cerebral aqueduct. CSF flow velocity is substantially slower in the areas that are furthest away from the inlets; in some areas flow is nearly stagnant.