Lagrangian Coherent Structures in Tandem Flapping Wing Hovering

Lagrangian Coherent Structures (LCS) of tandem wings hovering in an inclined stroke plane is studied using ImmersedBoundary Method (IBM) by solving two dimensional (2D) incompressible Navier-Stokes equations.Coherent structures responsible for the force variation are visualized by calculating Finite Time Lyapunov Exponents (FTLE),and vorticity contours.LCS is effective in determining the vortex boundaries,flow separation,and the wing-vortex interactions accurately.The effects of inter-wing distance and phase difference on the force generation are studied.Results show that in-phase stroking generates maximum vertical force and counter-stroking generates the least vertical force.In-phase stroking generates a wake with swirl,and counter stroking generates a wake with predominant vertical velocity.Counter stroking aids the stability of the body in hovering.As the hindwing operates in the wake of the forewing,due to the reduction in the effective Angle of Attack (AoA),the hindwing generates lesser force than that of a single flapping wing.     

Influence of cardiac flow rate on turbulent shear stress from a prosthetic heart valve

Elevated turbulent shear stresses associated with sufficient exposure times are potentially damaging to blood constituents. Since these conditions can be induced by mechanical heart valves, the objectives of this study were to locate the maximum turbulent shear stress in both space and time and to determine how the maximum turbulent shear stress depends on the cardiac flow rate in a pulsatile flow downstream of a tilting disk valve. Two-component, simultaneous, correlated laser velocimeter measurements were recorded at four different axial locations and three different flow rates in a straight tube model of the aorta. All velocity data were ensemble averaged within a 15 ms time window located at approximately peak systolic flow over more than 300 cycles. Shear stresses as high as 992 dynes/cm2 were found 0.92 tube diameters downstream of the monostrut, disk valve. The maximum turbulent shear stress was found to scale with flow rate to the 0.72 power. A repeatable starting vortex was shed from the disk at the beginning of each cycle.     

Pulsatile flow effects on the hemodynamics of intracranial aneurysms

High-resolution numerical simulations are carried out to systematically investigate the effect of the incoming flow waveform on the hemodynamics and wall shear stress patterns of an anatomic sidewall intracranial aneurysm model. Various wave forms are constructed by appropriately scaling a typical human waveform such that the waveform maximum and time-averaged Reynolds numbers, the Womersley number (α), and the pulsatility index (PI) are systematically varied within the human physiologic range. We show that the waveform PI is the key parameter that governs the vortex dynamics across the aneurysm neck and the flow patterns within the dome. At low PI, the flow in the dome is similar to a driven cavity flow and is characterized by a quasi-stationary shear layer that delineates the parent artery flow from the recirculating flow within the dome. At high PI, on the other hand, the flow is dominated by vortex ring formation, transport across the neck, and impingement and breakdown at the distal wall of the aneurysm dome. We further show that the spatial and temporal characteristics of the wall shear stress field on the aneurysm dome are strongly correlated with the vortex dynamics across the neck. We finally argue that the ratio between the characteristic time scale of transport by the mean flow across the neck and the time scale of vortex ring formation can be used to predict for a given sidewall aneurysm model the critical value of the waveform PI for which the hemodynamics will transition from the cavity mode to the vortex ring mode.     

The effect of vortex formation on left ventricular filling and mitral valve efficiency

A new mechanism for quantifying the filling energetics in the left ventricle (LV) and past mechanical heart valves (MHV) is identified and presented. This mechanism is attributed to vortex formation dynamics past MHV leaflets. Recent studies support the conjecture that the natural healthy left ventricle (LV) performs in an optimum, energy-preserving manner by redirecting the flow with high efficiency. Yet to date, no quantitative proof has been presented. The present work provides quantitative results and validation of a theory based on the dynamics of vortex ring formation, which is governed by a critical formation number (FN) that corresponds to the dimensionless time at which the vortex ring has reached its maximum circulation content, in support of this hypothesis. Herein, several parameters (vortex ring circulation, vortex ring energy, critical FN, hydrodynamic efficiencies, vortex ring propagation speed) have been quantified and presented as a means of bridging the physics of vortex formation in the LV. In fact, the diastolic hydrodynamic efficiencies were found to be 60, 41, and 29%, respectively, for the porcine, anti-anatomical, and anatomical valve configurations. This assessment provides quantitative proof of vortex formation, which is dependent of valve design and orientation, being an important flow characteristic and associated to LV energetics. Time resolved digital particle image velocimetry with kilohertz sampling rate was used to study the ejection of fluid into the LV and resolve the spatiotemporal evolution of the flow. The clinical significance of this study is quantifying vortex formation and the critical FN that can potentially serve as a parameter to quantify the LV filling process and the performance of heart valves.     

Microfiltration of recombinant yeast cells using a rotating disk dynamic filtration system

To develop a highly efficient cell harvest step under time constraint, a novel rotating disk dynamic filtration system was studied on the laboratory scale (0.147-ft.(2) nylon membrane) for concentrating recombinant yeast cells containing an intracellular product. The existing cross-flow microfiltration method yielded pseudo-steady state flux values below 25 LMH (L/m(2). h) even at low membrane loadings (10 L/ft.(2)). By creating high shear rates (up to 120,000(-1)) on the membrane surface using a rotating solid disk, this dynamic filter has demonstrated dramatically improved performance, presumably due to minimal cake buildup and reduced membrane fouling. Among the many factors investigated, disk rotating speed, which determines shear rates and flow patterns, was found to be the most important adjustable parameter. Our experimental results have shown that the flux increases with disk rotating speed, increases with transmembrane pressure at higher cell concentrations, and can be sustained at high levels under constant flux mode. At a certain membrane loading level, there was a critical speed below which it behaved similarly to a flat sheet system with equivalent shear. Average flux greater than 200 LMH has been demonstrated at 37-L/ft.(2) loading at maximum speed to complete sixfold concentration and 15-volume diafiltration for less than 100 min. An order of magnitude improvement over the crossflow microfiltration control was projected for large scale production. This superior performance, however, would be achieved at the expense of additional power input and heat dissipation, especially when cell concentration reaches above 80 g dry cell weight (DCW)/L. Although a positive linear relationship between power input and dynamic flux at a certain concentration factor has been established, high cell density associated with high viscosity impacted adversely on effective average shear rates and, eventually, severe membrane fouling, rather than cake formation, would limit the performance of this novel system. (c) 1995 John Wiley & Sons, Inc.     

A detailed fluid mechanics study of tilting disk mechanical heart valve closure and the implications to blood damage

Hemolysis and thrombosis are among the most detrimental effects associated with mechanical heart valves. The strength and structure of the flows generated by the closure of mechanical heart valves can be correlated with the extent of blood damage. In this in vitro study, a tilting disk mechanical heart valve has been modified to measure the flow created within the valve housing during the closing phase. This is the first study to focus on the region just upstream of the mitral valve occluder during this part of the cardiac cycle, where cavitation is known to occur and blood damage is most severe. Closure of the tilting disk valve was studied in a "single shot" chamber driven by a pneumatic pump. Laser Doppler velocimetry was used to measure all three velocity components over a 30 ms period encompassing the initial valve impact and rebound. An acrylic window placed in the housing enabled us to make flow measurements as close as 200 microm away from the closed occluder. Velocity profiles reveal the development of an atrial vortex on the major orifice side of the valve shed off the tip of the leaflet. The vortex strength makes this region susceptible to cavitation. Mean and maximum axial velocities as high as 7 ms and 20 ms were recorded, respectively. At closure, peak wall shear rates of 80,000 s(-1) were calculated close to the valve tip. The region of the flow examined here has been identified as a likely location of hemolysis and thrombosis in tilting disk valves. The results of this first comprehensive study measuring the flow within the housing of a tilting disk valve may be helpful in minimizing the extent of blood damage through the combined efforts of experimental and computational fluid dynamics to improve mechanical heart valve designs.     

Determinants of the occurrence of vortex rings in the left ventricle during diastole

This study employs classical inviscid fluid dynamics theory to investigate whether LV diastolic inflow volume and the size of the LV play a role in vortex ring formation. Fluid injection across an orifice into a large container results in the generation of a vortex ring having a constant size and speed. Relations between the vortex size and speed and the injection were obtained by applying conservation laws regarding kinetic energy, impulse and vorticity; the initial state was computed using a bolus injection model, and the final state by using the Kelvin vortex model. An important parameter in the equations is the relative injection length, i.e., the ratio of the length of the injected bolus and the radius of the orifice (L/R). Its estimated highest value in man, L/R = 15, produces a rather thick vortex ring (relative thickness 0.77). Comparable results following from the Hill vortex model convinced us that the Kelvin vortex model can be applied in the whole range of injection lengths in the human left ventricle. In an in vitro model it is shown experimentally that vortex rings can be generated for L/R in the range from 2 to 16. The measured traveling speed of the vortex ring is in fair agreement with the theory, as well as the ring radius for large injections. A vortex ring located in a narrow channel cannot reach its proper traveling speed. The method of images is used to estimate the speed reduction of vortex rings within a cylinder. It turns out that propagation of vortex rings is possible when the ratio of orifice to cylinder radius is less than about 0.5.(ABSTRACT TRUNCATED AT 250 WORDS)     

Computational analysis of confined jet flow and mass transport in a blind tube.

A computational analysis of confined nonimpinging jet flow in a blind tube is performed as an initial investigation of the underlying fluid and mass transport mechanics of tracheal gas insufflation. A two-dimensional axisymmetric model of a laminar steady jet flow into a concentric blind-end tube is put forth and the governing continuity, momentum, and convection-diffusion equations are solved with a finite element code. The effects of the jet diameter based Reynolds number (Re(j)), the ratio of the jet-to-outer tube diameters (epsilon), and the Schmidt number (Sc) are evaluated with the determined velocity and contaminant concentration fields. The normalized penetration depth of the jet is found to increase linearly with increasing Re(j) for epsilon = O(0.1). For a given epsilon, a ring vortex that develops is observed to be displaced downstream and radially outward from the jet tip for increasing Re(j). The axial shear stress profile along the inside wall of the outer tube possesses regions of fixed shear stress in addition to a local minimum and maximum in the vicinity of the jet tip. Corresponding regions of axial shear stress gradients exist between the fixed shear stress regions and the local extrema. Contaminant concentration gradients develop across the ring vortex indicating the inward diffusion of contaminant into the jet flow. For fixed epsilon and Sc and Re(j) approximately 900, normalized contaminant flow rate is observed to be approximately twice that of simple diffusion. This model predicts modest net axial contaminant transport enhancement due to convection-diffusion interaction in the region of the ring vortex.     

Flow of a blood analogue fluid in a compliant abdominal aortic aneurysm model: Experimental modelling

The aim of this work is to develop a unique in vitro set-up in order to analyse the influence of the shear thinning fluid-properties on the flow dynamics within the bulge of an abdominal aortic aneurysm (AAA). From an experimental point of view, the goals are to elaborate an analogue shear thinning fluid mimicking the macroscopic blood behaviour, to characterise its rheology at low shear rates and to propose an experimental device able to manage such an analogue fluid without altering its feature while reproducing physiological flow rate and pressure, through compliant AAA. Once these experimental prerequisites achieved, the results obtained in the present work show that the flow dynamics is highly dependent on the fluid rheology. The main results point out that the propagation of the vortex ring, generated in the AAA bulge, is slower for shear thinning fluids inducing a smaller travelled distance by the vortex ring so that it never impacts the anterior wall in the distal region, in opposition to Newtonian fluids. Moreover, scalar shear rate values are globally lower for shear thinning fluids inducing higher maximum stress values than those for the Newtonian fluids. Consequently, this work highlights that a Newtonian fluid model is finally inadequate to obtain a reliable prediction of the flow dynamics within AAA.     

Characteristics of secondary flow in steady and pulsatile flows through a symmetrical bifurcation

Steady and pulsatile flow in a glass model simulating an arterial bifurcation was investigated by flow visualization techniques. Secondary flow generated at the bifurcation has a similar pattern to a vortex, called the horseshoe vortex, produced around a wall-based protuberance in a circular tube. The same flow disturbance was clearly observed during the decelerating phase of pulsatile flow. The vortex produces a stagnation point on the top and bottom wall just upstream from the bifurcation apex. When aluminium dust was suspended in the test fluid perfusing the blood vessel model, particles deposited over an area spreading from the stagnation point to the lateral corners of the bifurcation. Comparison between the present results and topographical patterns of atherosclerosis reported in the literature suggests that it is in such low shear regions that lipid deposition tends to occur most.     

Chirality generated by flows in pseudocyanine dye J-aggregates: revisiting 40 years old reports

Spontaneous symmetry breaking in J-aggregates of cyanine dyes has a long history in chemical literature. In 1976, Honda and Hada claimed that they had achieved chiral induction (CD) by stirring J-aggregates of pseudocyanine. However, this report is controversial, as the combinations of linear dichroism and birefringence can lead to artifactual circular dichroic signals that are unrelated to molecular chirality. A Mueller matrix spectroscopy study, with an approach for the application of a gradient of the shear rate (solution layer between a rotating and a fixed disk) that differs from the simple vortex stirring used in the original report, shows that true CD can be induced in the sample. The phenomenon is discussed, taking into account the flow dynamics that allows the alignment of the aggregate particles and the gradient of shear rates that determines their folding/torsion, which leads to a chiral excitonic transition.     

A formulation for calculating the translational velocity of a vortex ring or pair

Cephalopods, among other marine animals, use jet propulsion for swimming. A simple actuator is designed to loosely mimic pulsatile jet formation in squid and jellyfish. The actuator is comprised of a cavity with an oscillating diaphragm and an exit orifice. Periodic oscillation of the diaphragm results in the formation of an array of vortex rings and eventually could generate a periodic pulsatile jet. A general formulation for calculating the velocity of a steadily translating vortical structure in two-dimensional and axi-symmetric shear flows is presented. This technique is based on taking the variational derivative of an energetic function at its critical point. This technique is general, applicable to vortices in liquid and gas media, with no limitation on the relative size of the vortex core. The technique is then implemented to estimate the translational velocity of a vortex ring in a Helmholtz vortex ring generator.     

Flow on the symmetry plane of a total cavo-pulmonary connection

The flow inside a total cavo-pulmonary connection, a bypass operation of the right heart adopted in the presence of congenital malformation, is here studied for a specific geometry which has been recently introduced in clinics. The analysis has been performed by preliminary experimental observation and a novel Navier-Stokes formulation on the symmetry plane. This method, once some basic hypotheses are verified, allows to reproduce the flow on the symmetry plane of a three-dimensional field by using an extension of the two-dimensional approach. The analysis has confirmed the existence of a central vortex showing that it is not a real vortex (i.e. a place with accumulation of vorticity) but, rather, a weakly dissipative recirculating zone. It is surrounded by a shear layer that becomes spontaneously unsteady at moderately high Reynolds number. The topological changes and energy dissipation have been analysed in both cases of unbalanced and of balanced pulmonary artery and caval flows.     

Vortex Formation with a Snapping Shrimp Claw

Snapping shrimp use one oversized claw to generate a cavitating high speed water jet for hunting, defence and communication. This work is an experimental investigation about the jet generation. Snapping shrimp (Alpheus-bellulus) were investigated by using an enlarged transparent model reproducing the closure of the snapper claw. Flow inside the model was studied using both High-Speed Particle Image Velocimetry (HS-PIV) and flow visualization. During claw closure a channel-like cavity was formed between the plunger and the socket featuring a nozzle-type contour at the orifice. Closing the mechanism led to the formation of a leading vortex ring with a dimensionless formation number of approximate ΔT*≈4. This indicates that the claw might work at maximum efficiency, i.e. maximum vortex strength was achieved by a minimum of fluid volume ejected. The subsequent vortex cavitation with the formation of an axial reentrant jet is a reasonable explanation for the large penetration depth of the water jet. That snapping shrimp can reach with their claw-induced flow. Within such a cavitation process, an axial reentrant jet is generated in the hollow cylindrical core of the cavitated vortex that pushes the front further downstream and whose length can exceed the initial jet penetration depth by several times.     

The vortex wake of a 'hovering' model hawkmoth

Visualization experiments with Manduca sexta have revealed the presence of a leading-edge vortex and a highly three-dimensional flow pattern. To further investigate this important discovery, a scaled-up robotic insect was built (the ''flapper'') which could mimic the complex movements of the wings of a hovering hawkmoth. Smoke released from the leading edge of the flapper wing revealed a small but strong leading-edge vortex on the downstroke. This vortex had a high axial flow velocity and was stable, separating from the wing at approximately 75 per cent of the wing length. It connected to a large, tangled tip vortex, extending back to a combining stopping and starting vortex from pronation. At the end of the downstroke, the wake could be approximated as one vortex ring per wing. Based on the size and velocity of the vortex rings, the mean lift force during the downstroke was estimated to be about 1.5 times the body weight of a hawkmoth, confirming that the downstroke is the main provider of lift force.     

The scaling and structure of aquatic animal wakes

Animal generated water movements are visualized and quantifiedusing two-dimensional particle image velocimetry (PIV). Theresulting vector flow fields allow for the study of the distributionof velocity, vorticity and vortices. Structural and temporalaspects of animal-induced flows covering a range of Reynolds(Re) numbers between less than 1 to more than 10⁴ are presented. Maps of flow induced by continuous foraging and intermittentescape responses of tethered nauplius and copepodid stages ofthe marine copepod Temora longicornis offer insight in viscosity-dominatedflow regimes. Fast escape responses of the equally sized largestnauplius stage and the smallest copepodid stage are compared.The nauplius moves by generating a viscous flow pattern withhigh velocities and vorticity; the copepodid moves by usinginertial effects to produce a vortex ring with a rearward jetthrough the center. Larvae and small adult fish (zebra danio) use a burst-and-coast-swimmingmode at Re numbers up to 6,000, shedding a vortex ring withthe associated jet at the tail during the burst phase. Flowpatterns during the coasting phase differ between the smalllarvae and larger adults due to the changes in importance ofviscosity. A 12 cm long mullet swimming in a continuous mode generatesa chain of vortex rings with a backward undulating jet throughthe centers of the rings at Re numbers of 4 x 10⁴ in inertia-dominatedregimes. Our empirical results provide realistic insight in the scaleeffects determining the morphology of the interactions betweenanimals and water.     

Two-component laser velocimeter measurements downstream of heart valve prostheses in pulsatile flow

Elevated turbulent shear stresses resulting from disturbed blood flow through prosthetic heart valves can cause damage to red blood cells and platelets. The purpose of this study was to measure the turbulent shear stresses occurring downstream of aortic prosthetic valves during in-vitro pulsatile flow. By matching the indices of refraction of the blood analog fluid and model aorta, correlated, simultaneous two-component laser velocimeter measurements of the axial and radial velocity components were made immediately downstream of two aortic prosthetic valves. Velocity data were ensemble averaged over 200 or more cycles for a 15-ms window opened at peak systolic flow. The systolic duration for cardiac flows of 8.4 L/min was 200 ms. Ensemble-averaged total shear stress levels of 2820 dynes/cm2 and 2070 dynes/cm2 were found downstream of a trileaflet valve and a tilting disk valve, respectively. These shear stress levels decreased with axial distance downstream much faster for the tilting disk valve than for the trileaflet valve.     

Experimental technique of measuring dynamic fluid shear stress on the aortic surface of the aortic valve leaflet

Aortic valve (AV) calcification is a highly prevalent disease with serious impact on mortality and morbidity. The exact cause and mechanism of the progression of AV calcification is unknown, although mechanical forces have been known to play a role. It is thus important to characterize the mechanical environment of the AV. In the current study, we establish a methodology of measuring shear stresses experienced by the aortic surface of the AV leaflets using an in vitro valve model and adapting the laser Doppler velocimetry (LDV) technique. The valve model was constructed from a fresh porcine aortic valve, which was trimmed and sutured onto a plastic stented ring, and inserted into an idealized three-lobed sinus acrylic chamber. Valve leaflet location was measured by obtaining the location of highest back-scattered LDV laser light intensity. The technique of performing LDV measurements near to biological surfaces as well as the leaflet locating technique was first validated in two phantom flow systems: (1) steady flow within a straight tube with AV leaflet adhered to the wall, and (2) steady flow within the actual valve model. Dynamic shear stresses were then obtained by applying the techniques on the valve model in a physiologic pulsatile flow loop. Results show that aortic surface shear stresses are low during early systole (<5 dyn/cm2) but elevated to its peak during mid to late systole at about 18-20 dyn/cm2. Low magnitude shear stress (<5 dyn/cm2) was observed during early diastole and dissipated to zero over the diastolic duration. Systolic shear stress was observed to elevate only with the formation of sinus vortex flow. The presented technique can also be used on other in vitro valve models such as congenitally geometrically malformed valves, or to investigate effects of hemodynamics on valve shear stress. Shear stress data can be used for further experiments investigating effects of fluid shear stress on valve biology, for conditioning tissue engineered AV, and to validate numerical simulations.     

Microfiltration and ultrafiltration of polysaccharides produced by fermentation using a rotating disk dynamic filtration system

The recovery of exopolysaccharides (EPS) produced by Sinorhizobium meliloti bacteria by dynamic microfiltration was investigated using a rotating disk device designed in our laboratory, equipped with a 0.2 microm nylon membrane. This system differs from commercially available systems by the presence of vanes on the disk which produce a very important increase in permeate flux while yielding excellent EPS transmission. For polymers produced under standard fermentation conditions (70 h at 30 degrees C), the mass flux rose to 650 g h(-1) m(-2) using a disk equipped with 2 mm vanes rotating at 2000 rpm against 380 g h(-1) m(-2) with a smooth disk at the same speed. The maximum flux observed was 1560 g h(-1) m(-2) with a 6-mm vanes disk rotating at 3000 rpm and a 36 degrees C broth. An interesting finding was that the permeate flux J(f) for various disks can be correlated by the same function of the mean shear stress at the membrane tau(wm) according to J(f) = 4.6 tau(wm) (0.717) for a 30 degrees C broth, showing that the effect of vanes is merely to increase the shear stress by raising the fluid core velocity between the membrane and the disk. With 6-mm vanes the core angular velocity was found to be 84% of disk velocity vs. 45% for a smooth disk. When the fermentation temperature was increased to 36 degrees C to produce a lower molecular weight polymer, the permeate flux rose by about 250%, much more than what could be expected from the reduction in permeate viscosity and followed the same power law with membrane shear stress as for 30 degrees C. The same device was equipped with a PES 50 kDa membrane to concentrate EPS by ultrafiltration. Permeate fluxes were of the order of 160 L h(-1) m(-2) at 2000 rpm and 30 degrees C with nearly complete EPS rejection. Finally, the net electrical power consumed by the disk was measured by subtracting the power consumed without fluid from the power during filtration at the same speed. This power increases with speed and with the presence of vanes, but since the gain provided by the vanes is very high, the specific energy per m(3) of permeate is minimal with the highest vanes tested (6 mm) and maximal for smooth disks.     

The three-dimensional leading-edge vortex of a 'hovering' model hawkmoth

Recent flow visualisation experiments with the hawkmoth, Manduca sexta, revealed small but clear leading-edge vortex and a pronounced three-dimensional flow. Details of this flow pattern were studied with a scaled-up, robotic insect (''the flapper'') that accurately mimicked the wing movements of a hovering hawkmoth. Smoke released from the leading edge of the flapper wing confirmed the existence of a small, strong and stable leading-edge vortex, increasing in size from wingbase to wingtip. Between 25 and 75 per cent of the wing length, its diameter increased approximately from 10 to 50 per cent of the wing chord. The leading-edge vortex had a strong axial flow veolocity, which stabilized it and reduced its diamater. The vortex separated from the wing at approximately 75 per cent of the wing length and thus fed vorticity into a large, tangled tip vortex. If the circulation of the leading-edge vortex were fully used for lift generation, it could support up to two-thirds of the hawkmoth''s weight during the downstroke. The growth of this circulation with time and spanwise position clearly identify dynamic stall as the unsteady aerodynamic mechanism responsible for high lift production by hovering hawkmoths and possibly also by many other insect species.