Modifying peripheral IV catheters with side holes and side slits results in favorable changes in fluid dynamic properties during the injection of iodinated contrast material

OBJECTIVE: The purpose of this study was to compare a standard peripheral end-hole angiocatheter with those modified with side holes or side slits using experimental optical techniques to qualitatively compare the contrast material exit jets and using numeric techniques to provide flow visualization and quantitative comparisons. MATERIALS AND METHODS: A Schlieren imaging system was used to visualize the angiocatheter exit jet fluid dynamics at two different flow rates. Catheters were modified by drilling through-and-through side holes or by cutting slits into the catheters. A commercial computational fluid dynamics package was used to calculate numeric results for various vessel diameters and catheter orientations. RESULTS: Experimental images showed that modifying standard peripheral IV angiocatheters with side holes or side slits qualitatively changed the overall flow field and caused the exiting jet to become less well defined. Numeric calculations showed that the addition of side holes or slits resulted in a 9-30% reduction of the velocity of contrast material exiting the end hole of the angiocatheter. With the catheter tip directed obliquely to the wall, the maximum wall shear stress was always highest for the unmodified catheter and was always lowest for the four-side-slit catheter. CONCLUSION: Modified angiocatheters may have the potential to reduce extravasation events in patients by reducing vessel wall shear stress.

Flow disturbances in stent-related coronary evaginations: a computational fluid-dynamic simulation study

Aims: Angiographic ectasias and aneurysms in stented segments have been associated with late stent thrombosis. Using optical coherence tomography (OCT), some stented segments show coronary evaginations reminiscent of ectasias. The purpose of this study was to explore, using computational fluid-dynamic (CFD) simulations, whether OCT-detected coronary evaginations can induce local changes in blood flow. Methods and results: OCT-detected evaginations are defined as outward bulges in the luminal vessel contour between struts, with the depth of the bulge exceeding the actual strut thickness. Evaginations can be characterised cross ectionally by depth and along the stented segment by total length. Assuming an ellipsoid shape, we modelled 3-D evaginations with different sizes by varying the depth from 0.2-1.0 mm, and the length from 1-9 mm. For the flow simulation we used average flow velocity data from non-diseased coronary arteries. The change in flow with varying evagination sizes was assessed using a particle tracing test where the particle transit time within the segment with evagination was compared with that of a control vessel. The presence of the evagination caused a delayed particle transit time which increased with the evagination size. The change in flow consisted locally of recirculation within the evagination, as well as flow deceleration due to a larger lumen - seen as a deflection of flow towards the evagination. Conclusions: CFD simulation of 3-D evaginations and blood flow suggests that evaginations affect flow locally, with a flow disturbance that increases with increasing evagination size.

Primary Oil Migration Through Buoyancy-Driven Multiple Fracture Propagation: Oil Velocity and Flux

We present a fracture-mechanics-based formulation to investigate primary oil migration through the propagation of an array of periodic, parallel fractures in a sedimentary rock with elevated pore fluid pressure. The rock is assumed to be a linearly elastic medium. The fracture propagation and hence oil migration velocity are determined using a fracture mechanics criterion together with the lubrication theory of fluid mechanics. We find that fracture interactions have profound effects on the primary oil migration behavior. For a given fracture length, the mass flux of oil migration decreases dramatically with an increase in fracture density. The reduced oil flux is due to the decreased fracture propagation velocity as well as the narrowed fracture opening that result from the fracture interactions.

Comparison of morphological and rheological conditions between conventional and eversion carotid endarterectomy using computational fluid dynamics - a pilot study

PURPOSE To compare postoperative morphological and rheological conditions after eversion carotid endarterectomy versus conventional carotid endarterectomy using computational fluid dynamics. BASIC METHODS Hemodynamic metrics (velocity, wall shear stress, time-averaged wall shear stress and temporal gradient wall shear stress) in the carotid arteries were simulated in one patient after conventional carotid endarterectomy and one patient after eversion carotid endarterectomy by computational fluid dynamics analysis based on patient specific data. PRINCIPAL FINDINGS Systolic peak of the eversion carotid endarterectomy model showed a gradually decreased pressure along the stream path, the conventional carotid endarterectomy model revealed high pressure (about 180 Pa) at the carotid bulb. Regions of low wall shear stress in the conventional carotid endarterectomy model were much larger than that in the eversion carotid endarterectomy model and with lower time-averaged wall shear stress values (conventional carotid endarterectomy: 0.03-5.46 Pa vs. eversion carotid endarterectomy: 0.12-5.22 Pa). CONCLUSIONS Computational fluid dynamics after conventional carotid endarterectomy and eversion carotid endarterectomy disclosed differences in hemodynamic patterns. Larger studies are necessary to assess whether these differences are consistent and might explain different rates of restenosis in both techniques.

Measurements of velocity-dependent and slip-dependent frictional strength in laboratory friction experiments on natural samples from IODP Expedition 316 and ODP Hole 190-1174B

Slowslip forms part of the spectrum of fault behaviour between stable creep and destructive earthquakes. Slow slip occurs near the boundaries of large earthquake rupture zones and may sometimes trigger fast earthquakes. It is thought to occur in faults comprised of rocks that strengthen under fast slip rates, preventing rupture as a normal earthquake, or on faults that have elevated pore-fluid pressures. However, the processes that control slow rupture and the relationship between slow and normal earthquakes are enigmatic. Here we use laboratory experiments to simulate faulting in natural rock samples taken from shallow parts of the Nankai subduction zone, Japan, where very low-frequency earthquakes - a form of slow slip - have been observed.We find that the fault rocks exhibit decreasing strength over millimetre-scale slip distances rather than weakening due to increasing velocity. However, the sizes of the slip nucleation patches in our laboratory simulations are similar to those expected for the very lowfrequency earthquakes observed in Nankai. We therefore suggest that this type of fault-weakening behaviour may generate slow earthquakes. Owing to the similarity between the expected behaviour of slow earthquakes based on our data, and that of normal earthquakes during nucleation, we suggest that some types of slow slip may represent prematurely arrested earthquakes.

(Table 1) Bulk and grain density, porosity, and p- and s-wave velocity characteristics of ODP Hole 163-990A sediments

Seismic velocities in rocks are influenced by the properties of the solid, the pore fluid, and the pore space. Cracks dramatically affect seismic velocities in rocks; their influence on the effective elastic moduli of rocks depends on their shape and concentration. Thin cracks (or fractures) substantially lower the moduli of a rock relative to the effect of spherical voids (or vesicles), and lower moduli are reflected by lower P- and S-wave velocities. The objective of this research is to determine the types and concentrations of cracks and their influence on the seismic properties of subaerially erupted basalts drilled from Hole 990A on the Southeast Greenland margin during Ocean Drilling Program Leg 163. Ellipsoidal cracks are used to model the voids in the rocks. The elastic moduli of the solid (grains) are also free parameters in the inverse modeling procedure. The apparent grain moduli reflect a weighted average of the moduli of the constituent minerals (e.g., plagioclase, augite, and clay minerals). The results indicate that (1) there is a strong relationship between P-wave velocity and porosity, suggesting a similarity of pore shape distributions, (2) the distribution of crack types within the massive, central region of aa flows from Hole 990A is independent of total porosity, (3) thin cracks are the first to be effectively sealed by alteration products, and (4) grain densities (an alteration index) and apparent grain moduli of the basalt samples are directly related.

Electrical resistivity, formation factors, and sound velocity at DSDP Holes 75-530A and 75-530B

From 0 to 277 m at Site 530 are found Holocene to Miocene diatom ooze, nannofossil ooze, marl, clay, and debrisflow deposits; from 277 to 467 m are Miocene to Oligocene mud; from 467 to 1103 m are Eocene to late Albian Cenomanian interbedded mudstone, marlstone, chalk, clastic limestone, sandstone, and black shale in the lower portion; from 1103 to 1121 m are basalts. In the interval from 0 to 467 m, in Holocene to Oligocene pelagic oozes, marl, clay, debris flows, and mud, velocities are 1.5 to 1.8 km/s; below 200 m velocities increase irregularly with increasing depth. From 0 to 100 m, in Holocene to Pleistocene diatom and nannofossil oozes (excluding debris flows), velocities are approximately equivalent to that of the interstitial seawater, and thus acoustic reflections in the upper 100 m are primarily caused by variations in density and porosity. Below 100 or 200 m, acoustic reflections are caused by variations in both velocity and density. From 100 to 467 m, in Miocene-Oligocene nannofossil ooze, clay, marl, debris flows, and mud, acoustic anisotropy irregularly increases to 10%, with 2 to 5% being typical. From 467 to 1103 m in Paleocene to late Albian Cenomanian interbedded mudstone, marlstone, chalk, clastic limestone, and black shale in the lower portion of the hole, velocities range from 1.6 to 5.48 km/s, and acoustic anisotropies are as great as 47% (1.0 km/s) faster horizontally. Mudstone and uncemented sandstone have anisotropies which irregularly increase with increasing depth from 5 to 10% (0.2 km/s). Calcareous mudstones have the greatest anisotropies, typically 35% (0.6 km/s). Below 1103 m, basalt velocities ranged from 4.68 to 4.98 km/s. A typical value is about 4.8 km/s. In situ velocities are calculated from velocity data obtained in the laboratory. These are corrected for in situ temperature, hydrostatic pressure, and porosity rebound (expansion when the overburden pressure is released). These corrections do not include rigidity variations caused by overburden pressures. These corrections affect semiconsolidated sedimentary rocks the most (up to 0.25 km/s faster). These laboratory velocities appear to be greater than the velocities from the sonic log. Reflection coefficients derived from the laboratory data, in general, agree with the major features on the seismic profiles. These indicate more potential reflectors than indicated from the reflection coefficients derived using the Gearhart-Owen Sonic Log from 625 to 940 m, because the Sonic Log data average thin beds. Porosity-density data versus depth for mud, mudstone, and pelagic oozes agree with data for similar sediments as summarized in Hamilton (1976). At depths of about 400 m and about 850 m are zones of relatively higher porosity mudstones, which may suggest anomalously high pore pressure; however, they are more probably caused by variations in grain-size distribution and lithology. Electrical resistivity (horizontal) from 625 to 950 m ranged from about 1.0 to 4.0 ohm-m, in Maestrichtian to Santonian- Coniacian mudstone, marlstone, chalk, clastic limestone, and sandstone. An interstitial-water resistivity curve did not indicate any unexpected lithology or unusual fluid or gas in the pores of the rock. These logs were above the black shale beds. From 0 to 100 m at Sites 530 and 532, the vane shear strength on undisturbed samples of Holocene-Pleistocene diatom and nannofossil ooze uniformly increases from about 80 g/cm**2 to about 800 g/cm**2. From 100 to 300 m, vane shear strength of Pleistocene-Miocene nannofossil ooze, clay, and marl are irregular versus depth with a range of 500 to 2300 g/cm**2; and at Site 532 the vane shear strength appears to decrease irregularly and slightly with increasing depth (gassy zone). Vane shear strength values of gassy samples may not be valid, for the samples may be disturbed as gas evolves, and the sediments may not be gassy at in situ depths.

Effective stress, porosity, p-wave velocity and mineral composition of ODP Hole 174A-1073A sediments

Porosity, permeability, and compressional (P-wave) velocity were measured as a function of stress on sediments from Ocean Drilling Program Site 1073, U.S. Mid-Atlantic continental slope. Thin sections, scanning electron microscopy, and X-ray diffraction analyses provided mineralogical characteristics of the samples. Uniaxial strain boundary conditions were imposed on the samples during consolidation tests with the maximum effective axial stress reaching 13 MPa. The maximum effective radial stress necessary to maintain uniaxial strain was 7.6 MPa. Over an effective axial stress interval of 0 to 5.2 MPa, Sample 174A-1073A-26X-2, 82-89 cm (226.65 meters below seafloor [mbsf]), exhibited the largest decrease in porosity (51% to 41%), whereas Sample 71X-1, 2-8 cm (644.70 mbsf), exhibited the smallest decrease in porosity (48% to 45%). All samples showed negligible porosity increases during unloading. The permeability (on the order of 1 x 10-17 m**2) of Sample 174A-1073A-71X-1, 2-8 cm, was twice that measured on Sample 8H-1, 23-26 cm (63.75 mbsf), even though the former was considerably deeper and older. The differences in porosity-stress behavior and permeability between shallow and deep samples is related to lithologic, mineralogic, and diagenetic differences between the sediments above and below the Pliocene-Pleistocene to Miocene unconformity. P-wave velocity for Samples 174A-1073A-41X-5, 97-103 cm (372.35 mbsf), and 71X-1, 2-8 cm, increased with decreasing porosity, but did not change significantly during unloading.

Exact solution for the conjugate fluid-fluid problem in the thermal entrance region of laminar counterflow heat exchangers

A generalized Lévêque solution is presented for the conjugate fluid–fluid problem that arises in the thermal entrance region of laminar counterflow heat exchangers. The analysis, carried out for constant property fluids, assumes that the Prandtl and Peclet numbers are both large compared to unity, and neglects axial conduction both in the fluids and in the plate, assumed to be thermally thin. Under these conditions, the thermal entrance region admits an asymptotic self-similar description where the temperature varies as a power ϳ of the axial distance, with the particularity that the self-similarity exponent must be determined as an eigenvalue by solving a transcendental equation arising from the requirement of continuity of heat fluxes at the heat conducting wall. Specifically, the analysis reveals that j depends only on the lumped parameter ƙ = (A2/A1)1/3 (α1/α2)1/3(k2/k1), defined in terms of the ratios of the wall velocity gradients, A, thermal diffusivities, α i, and thermal conductivities,k i, of the fluids entering, 1, and exiting, 2, the heat exchanger. Moreover, it is shown that for large (small) values of K solution reduces to the classical first (second) Lévêque solution. Closed-form analytical expressions for the asymptotic temperature distributions and local heat-transfer rate in the thermal entrance region are given and compared with numerical results in the counterflow parallel-plate configuration, showing very good agreement in all cases.

Numerical issues in Lagrangian tracking and topological evolution of fluid particles in wallbounded turbulent flows

The determination of the local Lagrangian evolution of the flow topology in wall-bounded turbulence, and of the Lagrangian evolution associated with entrainment across the turbulent / non-turbulent interface into a turbulent boundary layer, require accurate tracking of a fluid particle and its local velocity gradients. This paper addresses the implementation of fluid-particle tracking in both a turbulent boundary layer direct numerical simulation and in a fully developed channel flow simulation. Determination of the sub-grid particle velocity is performed using both cubic B-spline, four-point Hermite spline and higher-order Hermite spline interpolation. Both wall-bounded flows show similar oscillations in the Lagrangian tracers of both velocity and velocity gradients, corresponding to the movement of particles across the boundaries of computational cells. While these oscillation in the particle velocity are relatively small and have negligible effect on the particle trajectories for time-steps of the order of CFL = 0.1, they appear to be the cause of significant oscillations in the evolution of the invariants of the velocity gradient tensor.

Measurement of velocity in rotational flows using ultrasonic anemometry: the flowmeter

In this paper a previously developed theoretical model of the measurement process performed by a transit-time ultrasonic anemometer is applied to a fluid flowing through a circular section pipe. This model considers the influence of the shift of the acoustic pulse trajectory from straight propagation due to the flow on the measured speed. The aim of this work is to estimate the errors induced in the measured velocity by the shift of the acoustic pulse trajectory. Using different duct’s flow models, laminar and turbulent regimes have been analyzed. The results show that neglecting the effect of shift of the acoustic pulse trajectory leads to flow rate measurement underestimation.

Discrete linear stability and sensitivity analysis of fluid flows

Este trabajo presenta un método discreto para el cálculo de estabilidad hidrodinámica y análisis de sensibilidad a perturbaciones externas para ecuaciones diferenciales y en particular para las ecuaciones de Navier-Stokes compressible. Se utiliza una aproximación con variable compleja para obtener una precisión analítica en la evaluación de la matriz Jacobiana. Además, mapas de sensibilidad para la sensibilidad a las modificaciones del flujo de base y a una fuerza constante permiten identificar las regiones del campo fluido donde una modificacin (ej. fuerza puntual) tiene un efecto estabilizador del flujo. Se presentan cuatro casos de prueba: (1) un caso analítico para comprobar la derivación discreta, (2) una cavidad cerrada a bajo Reynolds para mostrar la mayor precisión en el cálculo de los valores propios con la aproximación de paso complejo, (3) flujo 2D en un cilindro circular para validar la metodología, y (4) flujo en un cavidad abierta, presentado para validar el método en casos de inestabilidades convectivamente inestables. Los tres últimos casos mencionados (2-4) se resolvieron con las ecuaciones de Navier-Stokes compresibles, utilizando un método Discontinuous Galerkin Spectral Element Method. Se obtuvo una buena concordancia para el caso de validación (3), cuando se comparó el nuevo método con resultados de la literatura. Además, este trabajo muestra que para el cálculo de los modos propios directos y adjuntos, así como para los mapas de sensibilidad, el uso de variables complejas es de suprema importancia para obtener una predicción precisa. El método descrito es aplicado al análisis para la estabilización de la estela generada por un disco actuador, que representa un modelo sencillo para hélices, rotores de helicópteros o turbinas eólicas. Se explora la primera bifurcación del flujo para un disco actuador, y se sugiere que está asociada a una inestabilidad de tipo Kelvin-Helmholtz, cuya estabilidad se controla con en el número de Reynolds y en la resistencia del disco actuador (o fuerza resistente). En primer lugar, se verifica que la disminución de la resistencia del disco tiene un efecto estabilizador parecido a una disminución del Reynolds. En segundo lugar, el análisis hidrodinmico discreto identifica dos regiones para la colocación de una fuerza puntual que controle las inestabilidades, una cerca del disco y otra en una zona aguas abajo. En tercer lugar, se muestra que la inclusión de un forzamiento localizado cerca del actuador produce una estabilización más eficiente que al forzar aguas abajo. El análisis de los campos de flujo controlados confirma que modificando el gradiente de velocidad cerca del actuador es más eficiente para estabilizar la estela. Estos resultados podrían proporcionar nuevas directrices para la estabilización de la estela de turbinas de viento o de marea cuando estén instaladas en un parque eólico y minimizar las interacciones no estacionarias entre turbinas. ABSTRACT A discrete framework for computing the global stability and sensitivity analysis to external perturbations for any set of partial differential equations is presented. In particular, a complex-step approximation is used to achieve near analytical accuracy for the evaluation of the Jacobian matrix. Sensitivity maps for the sensitivity to base flow modifications and to a steady force are computed to identify regions of the flow field where an input could have a stabilising effect. Four test cases are presented: (1) an analytical test case to prove the theory of the discrete framework, (2) a lid-driven cavity at low Reynolds case to show the improved accuracy in the calculation of the eigenvalues when using the complex-step approximation, (3) the 2D flow past a circular cylinder at just below the critical Reynolds number is used to validate the methodology, and finally, (4) the flow past an open cavity is presented to give an example of the discrete method applied to a convectively unstable case. The latter three (2–4) of the aforementioned cases were solved with the 2D compressible Navier–Stokes equations using a Discontinuous Galerkin Spectral Element Method. Good agreement was obtained for the validation test case, (3), with appropriate results in the literature. Furthermore, it is shown that for the calculation of the direct and adjoint eigenmodes and their sensitivity maps to external perturbations, the use of complex variables is paramount for obtaining an accurate prediction. An analysis for stabilising the wake past an actuator disc, which represents a simple model for propellers, helicopter rotors or wind turbines is also presented. We explore the first flow bifurcation for an actuator disc and it suggests that it is associated to a Kelvin- Helmholtz type instability whose stability relies on the Reynolds number and the flow resistance applied through the disc (or actuator forcing). First, we report that decreasing the disc resistance has a similar stabilising effect to an decrease in the Reynolds number. Second, a discrete sensitivity analysis identifies two regions for suitable placement of flow control forcing, one close to the disc and one far downstream where the instability originates. Third, we show that adding a localised forcing close to the actuator provides more stabilisation that forcing far downstream. The analysis of the controlled flow fields, confirms that modifying the velocity gradient close to the actuator is more efficient to stabilise the wake than controlling the sheared flow far downstream. An interesting application of these results is to provide guidelines for stabilising the wake of wind or tidal turbines when placed in an energy farm to minimise unsteady interactions.

Behavior of several two-dimensional fluid equations in singular scenarios

We give conditions that rule out formation of sharp fronts for certain two-dimensional incompressible flows. We show that a necessary condition of having a sharp front is that the flow has to have uncontrolled velocity growth. In the case of the quasi-geostrophic equation and two-dimensional Euler equation, we obtain estimates on the formation of semi-uniform fronts.

Chemical evolution and estimated flow velocity of water in the Trinity aquifer, south-central Texas /

Shipping list no.: 98-0242-P.

Prediction of slurry transport in SAG mills using SPH fluid flow in a dynamic DEM based porous media

DEM modelling of the motion of coarse fractions of the charge inside SAG mills has now been well established for more than a decade. In these models the effect of slurry has broadly been ignored due to its complexity. Smoothed particle hydrodynamics (SPH) provides a particle based method for modelling complex free surface fluid flows and is well suited to modelling fluid flow in mills. Previous modelling has demonstrated the powerful ability of SPH to capture dynamic fluid flow effects such as lifters crashing into slurry pools, fluid draining from lifters, flow through grates and pulp lifter discharge. However, all these examples were limited by the ability to model only the slurry in the mill without the charge. In this paper, we represent the charge as a dynamic porous media through which the SPH fluid is then able to flow. The porous media properties (specifically the spatial distribution of porosity and velocity) are predicted by time averaging the mill charge predicted using a large scale DEM model. This allows prediction of transient and steady state slurry distributions in the mill and allows its variation with operating parameters, slurry viscosity and slurry volume, to be explored. (C) 2006 Published by Elsevier Ltd.