Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “peak flows”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

Spatially Resolving Ocean Color and Sediment Dispersion in River Plumes, Coastal Systems, and Continental Shelf Waters

Satellite remote sensing of ocean color in dynamic coastal, inland, and nearshorewaters is impeded by high variability in optical constituents, demands specialized atmospheric correction, and is limited by instrument sensitivity. To accurately detect dispersion of bio-optical properties, remote sensors require ample signal-to-noise ratio (SNR) to sense small variations in ocean color without saturating over bright pixels, an atmospheric correction that can accommodate significantwater-leaving radiance in the near infrared (NIR), and spatial and temporal resolution that coincides with the scales of variability in the environment. Several current and historic space-borne sensors have met these requirements with success in the open ocean, but are not optimized for highly red-reflective and heterogeneous waters such as those found near river outflows or in the presence of sediment resuspension. Here we apply analytical approaches for determining optimal spatial resolution, dominant spatial scales of variability ("patches"), and proportions of patch variability that can be resolved from four river plumes around the world between 2008 and 2011. An offshore region in the Sargasso Sea is analyzed for comparison. A method is presented for processing Moderate Resolution Imaging Spectroradiometer (MODIS) Aqua and Terra imagery including cloud detection, stray lightmasking, faulty detector avoidance, and dynamic aerosol correction using short-wave- and near-infrared wavebands in extremely turbid regions which pose distinct optical and technical challenges. Results showthat a pixel size of approx. 520 mor smaller is generally required to resolve spatial heterogeneity in ocean color and total suspended materials in river plumes. Optimal pixel size increases with distance from shore to approx. 630 m in nearshore regions, approx 750 m on the continental shelf, and approx. 1350 m in the open ocean. Greater than 90% of the optical variability within plume regions is resolvable with 500 m resolution, and small, but significant, differences were found between peak and nadir river flow periods in terms of optimal resolution and resolvable proportion of variability.

MODIS↗

Method for Comparing Porous Carbon Electrode Performance in Redox Flow Batteries

In this paper, we present a method for comparing porous electrodes for use in redox flow batteries. By comparing physical and electrochemical properties of six carbon cloth electrodes (CCEs) including hydraulic permeability, Ohmic resistance, and mass transport overpotential, we show how optimal engineering conditions such as compression and electrolyte flow rates can be derived. These results are combined and applied to a quinone-based redox flow battery, leading to a peak power density around 0.4 W cm –2 for several different CCEs. This method can be applied broadly to other electrode architectures and flow battery chemistries, and provides a comprehensive framework for their comparison.

25 ENERGY STORAGE↗

C/NOFS Satellite Electric Field and Plasma Density Observations of Plasma Instabilities Below the Equatorial F-Peak -- Evidence for Approximately 500 km-Scale Spread-F "Precursor" Waves Driven by Zonal Shear Flow and km-Scale, Narrow-Banded Irregularities

As solar activity has increased, the ionosphere F-peak has been elevated on numerous occasions above the C/NOFS satellite perigee of 400km. In particular, during the month of April, 2011, the satellite consistently journeyed below the F-peak whenever the orbit was in the region of the South Atlantic anomaly after sunset. During these passes, data from the electric field and plasma density probes on the satellite have revealed two types of instabilities which had not previously been observed in the C/NOFS data set (to our knowledge): The first is evidence for 400-500km-scale bottomside "undulations" that appear in the density and electric field data. In one case, these large scale waves are associated with a strong shear in the zonal E x B flow, as evidenced by variations in the meridional (outward) electric fields observed above and below the F-peak. These undulations are devoid of smaller scale structures in the early evening, yet appear at later local times along the same orbit associated with fully-developed spread-F with smaller scale structures. This suggests that they may be precursor waves for spread-F, driven by a collisional shear instability, following ideas advanced previously by researchers using data from the Jicamarca radar. A second new result (for C/NOFS) is the appearance of km-scale irregularities that are a common feature in the electric field and plasma density data that also appear when the satellite is below the F -peak at night. The vector electric field instrument on C/NOFS clearly shows that the electric field component of these waves is strongest in the zonal direction. These waves are strongly correlated with simultaneous observations of plasma density oscillations and appear both with, and without, evidence of larger-scale spread-F depletions. These km-scale, quasi-coherent waves strongly resemble the bottomside, sinusoidal irregularities reported in the Atmosphere Explorer satellite data set by Valladares et al. [JGR, 88, 8025, 1983]. We interpret these new observations in terms of fundamental plasma instabilities associated with the unstable, nighttime equatorial ionosphere.

Pfaff, R.↗

Radiative Extinction of Gaseous Spherical Diffusion Flames in Microgravity

Radiative extinction of spherical diffusion flames was investigated experimentally and numerically. The experiments involved microgravity spherical diffusion flames burning ethylene and propane at 0.98 bar. Both normal (fuel flowing into oxidizer) and inverse (oxidizer flowing into fuel) flames were studied, with nitrogen supplied to either the fuel or the oxygen. Flame conditions were chosen to ensure that the flames extinguished within the 2.2 s of available test time; thus extinction occurred during unsteady flame conditions. Diagnostics included color video and thin-filament pyrometry. The computations, which simulated flow from a porous sphere into a quiescent environment, included detailed chemistry, transport and radiation, and yielded transient results. Radiative extinction was observed experimentally and simulated numerically. Extinction time, peak temperature, and radiative loss fraction were found to be independent of flow rate except at very low flow rates. Radiative heat loss was dominated by the combustion products downstream of the flame and was found to scale with flame surface area, not volume. For large transient flames the heat release rate also scaled with surface area and thus the radiative loss fraction was largely independent of flow rate. Peak temperatures at extinction onset were about 1100 K, which is significantly lower than for kinetic extinction. One observation of this work is that while radiative heat losses can drive transient extinction, this is not because radiative losses are increasing with time (flame size) but rather because the heat release rate is falling off as the temperature drops.

Santa, K. J.↗

Discharge Oscillations in a Permanent Magnet Cylindrical Hall-Effect Thruster

Measurements of the discharge current in a cylindrical Hall thruster are presented to quantify plasma oscillations and instabilities without introducing an intrusive probe into the plasma. The time-varying component of the discharge current is measured using a current monitor that possesses a wide frequency bandwidth and the signal is Fourier transformed to yield the frequency spectra present, allowing for the identification of plasma oscillations. The data show that the discharge current oscillations become generally greater in amplitude and complexity as the voltage is increased, and are reduced in severity with increasing flow rate. The breathing mode ionization instability is identified, with frequency as a function of discharge voltage not increasing with discharge voltage as has been observed in some traditional Hall thruster geometries, but instead following a scaling similar to a large-amplitude, nonlinear oscillation mode recently predicted in for annular Hall thrusters. A transition from lower amplitude oscillations to large relative fluctuations in the oscillating discharge current is observed at low flow rates and is suppressed as the mass flow rate is increased. A second set of peaks in the frequency spectra are observed at the highest propellant flow rate tested. Possible mechanisms that might give rise to these peaks include ionization instabilities and interactions between various oscillatory modes.

Polzin, K. A.↗

Rheoencephalographic (REG) Assessment of Head and Neck Cooling for use with Multiple Sclerosis Patients

We have prepared a computer program (RHEOSYS:RHEOencephalographic impedance trace scanning SyStem) that can be used to automate the analysis of segmental impedance blood flow waveforms. This program was developed to assist in the post test analysis of recorded impedance traces from multiple segments of the body. It incorporates many of the blood flow, segmental volume, and vascular state indices reported in the world literature. As it is currently programmed, seven points are selected from each blood flow pulse and associated ECG waveforrn: 1. peak of the first ECG QRS complex, 2. start of systolic slope on the blood flow trace, 3. maximum amplitude of the impedance pulse, 4. position of the dicrotic notch, 5. maximum amplitude of the postdicrotic segment, 6. peak of the second ECG QRS complex, and 7. start of the next blood flow pulse. These points we used to calculate various geometric, area, and time-related values associated with the impedance pulse morphology. RHEOSYS then calculates a series of 34 impedance and cardiac cycle parameters which include pulse amplitudes; areas; pulse propagation times; cardiac cycle times; and various measures of arterial and various tone, contractility, and pulse volume. We used this program to calculate the scalp and intracranial blood flow responses to head and neck cooling as it may be applied to lower the body temperatures of multiple sclerosis patients. Twelve women and twelve men were tested using a commercially available head and neck cooling system operated at its maximum cooling capacity for a period of 30 minutes. Head and neck cooling produced a transient change in scalp blood flow and a significant, (P<0.05) decrease of approx. 30% in intracranial blood flow. Results of this experiment will illustrate how REG and RHEOSYS can be used in biomedical applications.

Montogomery, Leslie D.↗

Ambiguous Stability of Glaciers at Bed Peaks

Increasing ice flux from glaciers retreating over deepening (retrograde) bed topography has been implicated in the recent acceleration of mass loss from the Greenland and Antarctic ice sheets. We show in observations that some glaciers have remained at peaks in bed topography without retreating despite enduring significant changes in climate. Observations also indicate that some glaciers which persist at bed peaks undergo sudden retreat years or decades after the onset of local ocean or atmospheric warming. Using model simulations, we show that persistence of a glacier at a bed peak is caused by ice slowing as it flows up a reverse-sloping bed to the peak. Persistence at bed peaks may lead to two very different future behaviors for a glacier: one where it persists at a bed peak indefinitely, and another where it retreats from the bed peak after potentially long delays following climate forcing. However, it is nearly impossible to distinguish which of these two future behaviors will occur from current observations. We conclude that inferring glacier stability from observations of persistence obscures our true commitment to future sea-level rise under climate change. We recommend that further research is needed on seemingly stable glaciers to determine their likely future.

ice dynamics↗

Interferometric investigations of compressible dynamic stall over a transiently pitching airfoil

The dynamic stall flow field over NACA 0012 airfoil pitching transiently from 0 - 60 at a constant rate under compressible flow conditions has been studied using the real-time technique of point diffraction interferometry. This investigation using nonintrusive diagnostics provides a quantitative description of the overall flow field, including the finer details of dynamic stall vortex formation, growth and the concomitant changes in the pressure distribution. Analysis of several hundred interferograms obtained for a range of flow conditions shows that the peak leading edge suction pressure coefficient that stall is nearly constant for a given free stream Mach number at all nondimensional pitch rates. Also, this value is below that seen in steady flow at static stall for the same Mach number, indicating that dynamic effects significantly effect the separation behavior. Further, for a given Mach number, the dynamic stall vortex seems to form rapidly at nearly the same angle of attack for all pitch rates studied. As the vortex is shed, it induces an anti-clockwise trailing edge vortex, which grows in a manner similar to that of a starting vortex. The measured peak suction pressure coefficient drops as the free stream Mach number increases. For free stream Mach numbers above 0.4, small multiple shocks appear near the leading edge.

Chandrasekhara, M. S.↗

Experimental Investigation of Unsteady Shock Wave Turbulent Boundary Layer Interactions About a Blunt Fin

A series of experiments were performed to investigate the effects of Mach number variation on the characteristics of the unsteady shock wave/turbulent boundary layer interaction generated by a blunt fin. A single blunt fin hemicylindrical leading edge diameter size was used in all of the experiments which covered the Mach number range from 2.0 to 5.0. The measurements in this investigation included surface flow visualization, static and dynamic pressure measurements, both on centerline and off-centerline of the blunt fin axis. Surface flow visualization and static pressure measurements showed that the spatial extent of the shock wave/turbulent boundary layer interaction increased with increasing Mach number. The maximum static pressure, normalized by the incoming static pressure, measured at the peak location in the separated flow region ahead of the blunt fin was found to increase with increasing Mach number. The mean and standard deviations of the fluctuating pressure signals from the dynamic pressure transducers were found to collapse to self-similar distributions as a function of the distance perpendicular to the separation line. The standard deviation of the pressure signals showed initial peaked distribution, with the maximum standard deviation point corresponding to the location of the separation line at Mach number 3.0 to 5.0. At Mach 2.0 the maximum standard deviation point was found to occur significantly upstream of the separation line. The intermittency distributions of the separation shock wave motion were found to be self-similar profiles for all Mach numbers. The intermittent region length was found to increase with Mach number and decrease with interaction sweepback angle. For Mach numbers 3.0 to 5.0 the separation line was found to correspond to high intermittencies or equivalently to the downstream locus of the separation shock wave motion. The Mach 2.0 tests, however, showed that the intermittent region occurs significantly upstream of the separation line. Power spectral densities measured in the intermittent regions were found to have self-similar frequency distributions when compared as functions of a Strouhal number for all Mach numbers and interaction sweepback angles. The maximum zero-crossing frequencies were found to correspond with the peak frequencies in the power spectra measured in the intermittent region.

Barnhart, Paul J.↗

DASSH-F: Subchannel Based Thermal Analysis

The DASSH thermal analysis code is designed to rapidly allow a reactor design engineer to obtain flow rates requirements that satisfy peak temperature constraints in the domain. The advantage of using DASSH over a hand calculation is that it has a more rigorous treatment of the pin power distribution and coolant heat transfer within an assembly and between assemblies. The advantage of using DASSH over a conventional 3D subchannel code or a computational fluid dynamics code (CFD) is that it can obtain the desired solution in a matter of minutes in serial with minor computer memory needs. The DASSH methodology is virtually identical to SUPERENERGY-2 with additional functionalities taken from follow on work to SUPERENERGY-2 done at ANL in the 1980s. DASSH today is an integral component of the Argonne Fast Reactor analysis suite for reactor design work. DASSH obtains the power distribution from a coupled neutron-gamma heating calculation in GAMSOR (including DIF3D) at each time point of a companion fuel cycle analysis calculation with REBUS. The domain in DASSH assumes a hexagonal grid typical for fast reactors with much of the geometry information taken from the DIF3D model. DASSH assumes the assemblies that are loaded into each grid position are ducted to control the coolant flow. Considerable detail is given on the subchannel formulation of DASSH in this document. Much of the formulation and design of the code builds upon research done by previous authors with little new investigation. Thus the decisions made in developing the subchannel model used in DASSH have their origins over 50 years ago. Much of the heat transfer methodology in DASSH is built upon correlations for both the coolant mixing and heat transfer coefficients for pins and ducts. DASSH is thus not a rigorous treatment of a given problem, but a rapid assessment of the temperature field that has known limitations with respect to an experimental measurement or CFD calculation. The DASSH input and output are detailed along with usage of the software. The DASSH output provides tables of evaluated material properties and key coolant and pin temperature results. DASSH can create Python scripts that generate domain summary pictures. DASSH can also generate assembly temperature maps and VTK output files which allow the DASSH solution to be visualized. As the primary purpose of the DASSH software is to compute the coolant and fuel pin temperature distribution for a given model of a reactor, much of the output focus is giving the user quick summary tables needed to assess the performance of a given orifice flow specification. The present version of DASSH has a crude orifice search capability and an efficient orifice flow search capability. The flow search tries to meet user specified constraints for 1) peak 2-sigma clad midwall temperature, 2) peak coolant temperature, and 3) desired bulk outlet temperature. There are many development shortcomings in DASSH detailed in this document, but this version is functional for modern analysis needs. This document serves as the manual for the Fortran based DASSH software that was developed to replace the Python version of DASSH developed as part of the VTR program.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DASSH-F: Subchannel Based Thermal Analysis

The DASSH thermal analysis code is designed to rapidly allow a reactor design engineer to obtain flow rates requirements that satisfy peak temperature constraints in the domain. The advantage of using DASSH over a hand calculation is that it has a more rigorous treatment of the pin power distribution and coolant heat transfer within an assembly and between assemblies. The advantage of using DASSH over a conventional 3D subchannel code or a computational fluid dynamics code (CFD) is that it can obtain the desired solution in a matter of minutes in serial with minor computer memory needs. The DASSH methodology for pin lattice models is virtually identical to SUPERENERGY-2 with additional functionalities taken from follow on work to SUPERENERGY-2 done at ANL in the 1980s. DASSH today is an integral component of the Argonne Fast Reactor analysis suite for reactor design work. DASSH obtains the power distribution from a coupled neutron-gamma heating calculation in GAMSOR (including DIF3D) at each time point of a companion fuel cycle analysis calculation with REBUS. The domain in DASSH assumes a hexagonal grid typical for fast reactors with much of the geometry information taken from the DIF3D model. DASSH assumes the assemblies that are loaded into each grid position are ducted to control the coolant flow. The user can alternatively provide their own geometry and power profile instead of inheriting it from DIF3D. Considerable detail is given on the subchannel formulation of DASSH in this document. Much of the formulation and design of the code builds upon research done by previous authors with little new investigation. Thus the decisions made in developing the subchannel model used in DASSH have their origins over 50 years ago. Much of the heat transfer methodology in DASSH is built upon correlations for both the coolant mixing and heat transfer coefficients for pins and ducts. DASSH is thus not a rigorous treatment of a given problem, but a rapid assessment of the temperature field that has known limitations with respect to an experimental measurement or CFD calculation. The DASSH input and output are detailed along with usage of the software. The DASSH output provides tables of evaluated material properties and key coolant and pin temperature results. DASSH can create Python scripts that generate domain summary pictures. DASSH can also generate assembly temperature maps and VTK output files which allow the DASSH solution to be visualized. As the primary purpose of the DASSH software is to compute the coolant and fuel pin temperature distribution for a given model of a reactor, much of the output focus is giving the user quick summary tables needed to assess the performance of a given orifice flow specification. The present version of DASSH has a crude orifice search capability and a sufficient orifice flow search capability. The flow search tries to meet user specified constraints for 1) peak 2-sigma clad temperature, 2) peak coolant temperature, and 3) desired bulk outlet temperature. This document serves as the manual for the Fortran based DASSH software that was developed to replace the Python version of DASSH developed as part of the VTR program.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Performance of transonic fan stage with weight flow per unit annulus area of 178 kilograms per second per square meter (6.5(lb/sec)/(sq ft))

The overall and blade-element performances are presented over the stable flow operating range from 50 to 100 percent of design speed. Stage peak efficiency of 0.834 was obtained at a weight flow of 26.4 kg/sec (58.3 lb/sec) and a pressure ratio of 1.581. The stall margin for the stage was 7.5 percent based on weight flow and pressure ratio at stall and peak efficiency conditions. The rotor minimum losses were approximately equal to design except in the blade vibration damper region. Stator minimum losses were less than design except in the tip and damper regions.

Moore, R. D.↗

Numerical study of unsteady viscous hypersonic blunt body flows with an impinging shock

A complex two-dimensional, unsteady, viscous hypersonic shock wave interaction is numerically simulated by a high-resolution, second-order fully implicit shock-capturing scheme. The physical model consists of a nonstationary oblique shock impinging on the bow shock of a blunt body. Studies indicated that the unsteady flow patterns are slightly different from their steady counterparts. However, for the sample cases investigated the peak surface pressures for the unsteady flows seem to occur at very different impingement locations than for the steady flow cases.

Klopfer, G. H.↗

Numerical study of unsteady viscous hypersonic blunt body flows with an impinging shock

A complex two-dimensional, unsteady, viscous hypersonic shock wave interaction is numerically simulated by a high-resolution, second-order fully implicit shock-capturing scheme. The physical model consists of a nonstationary oblique shock impinging on the bow shock of a blunt body. Studies indicate that the unsteady flow patterns are slightly different from their steady counterparts. However, for the sample cases investigated the peak surface pressures for the unsteady flows seem to occur at very different impingement locations than for the steady flow cases.

Klopfer, G. H.↗

Requirements Description of DASSH-F

This report reviews the modeling and simulation capabilities of Argonne National Laboratory’s DASSH code that is used in present reactor analysis activities. These capabilities will be used to establish the set of verification tasks necessary to verify DASSH for use on commercial projects. A similar approach was taken for the PERSENT, REBUS and DIF3D software packages. The DASSH program is a thermal analysis code designed to rapidly allow a reactor design engineer to obtain flow rates requirements that satisfy peak temperature constraints in the domain. DASSH is a follow-on development to the SE2-ANL software and SUPERENERGY-2 software that it is based upon. DASSH was designed to account for both neutron and gamma heating and is inherently connected to the GAMSOR part of the ARC suite of fast reactor analysis software. SE2-ANL is a developed piece of software from the 1980s while DASSH is a modern implementation with notable improvements in geometry handling. The most important upgrade of DASSH relative to SE2-ANL is that it can analyze multiple time points in a single run where SE2-ANL can only treat a single time point. This allows the user to understand the impact of and search the flow distribution for the entire operational period of a reactor design considering pressure drop, peak coolant and fuel temperatures, and thermal striping. DASSH has three input paths that have to be verified. The first input path builds the geometry and power distribution based upon the DIF3D model but ignores the gamma heating aspects of the problem. The second input path also builds the geometry from the DIF3D model but it takes the neutron and gamma heating distributions from GAMSOR. The third input path is to take the geometry and power distribution directly from user input (i.e. not coupled to DIF3D or GAMSOR). DASSH also has many built in correlations for material properties along with a user defined specification of the fuel, structure, and coolant properties. There are correlations for flow split, mixing, pressure drop, and heat transfer coefficients (subchannel rather than a direct methodology). In total, verification of DASSH will require an extensive testing to cover all possible user features of the software.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Experimental investigation of turbulent wake flows in a helically wrapped rod bundle in presence of localized blockages

In nuclear sodium fast reactors, bundles of rods are tightly packed into a triangular lattice, enclosed in a hexagonal duct, and each pin is spirally wrapped with a thin wire. Flow blockages can potentially impact the local flow characteristics and heat transfer mechanisms in the bundle due to its small subchannel size. The effects of the blockage on the flow structures and heat transfer mechanisms are important aspects that require an accurate investigation. In this study, the flow-field characteristics in the vicinity of a blockage located in the exterior subchannel of rod bundles with helically wrapped wires were experimentally investigated. The velocity fields in the exterior subchannel were acquired by applying matched-index-of-refraction and time-resolved particle image velocimetry (TR-PIV) techniques for Reynolds numbers of Re 1 = 4000 and Re 2 = 17 000, i.e., equivalent to Rew 1 = 19 600 and Rew 2 = 83 200, respectively, based on the blockage width. The results from the TR-PIV measurements revealed an arch-shaped vortex with a large flow recirculation and a pair of counter-rotating vortices in the wake region downstream of the blockage, which is commonly observed in the wake flow of bluff bodies. The relative lateral distance and angle between the two vortices decreased when the Reynolds numbers increased. Profiles of maximum turbulence intensity along the shear layers illustrated the transition process including the growth, peak, and decay along the flow direction. From the spectral analysis of the turbulent velocities extracted at points along the shear layer, the Strouhal numbers ( St ) representing the vortex shedding frequency were found to be St = 0.25 and St = 0.56 for the left and right shear layers, respectively. Characteristics of shear layers generated by the blockage in the exterior subchannel were investigated via the two-point cross correlation of fluctuating velocities. The patiotemporal cross correlations of turbulent velocities, computed at points in the region where the left shear layer exhibited rolling effects and vortex breakdowns, were considerably wider and longer. The convection velocity U c was estimated to be ~0.82 U m to 0.93 U m . Proper orthogonal decomposition (POD) analysis was applied to the instantaneous velocity fields to extract the statistically dominant flow structures. It was found that POD modes 2–3 and 4–5 formed the pair modes when the corresponding POD temporal coefficients depicted sinusoidal shapes and exhibited nearly circular orbits in the phase space. Spectral analysis of the POD temporal coefficients confirmed the vortex shedding frequencies detected in the analysis of turbulent velocities.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Performance of inlet stage of transonic compressor

The overall and blade-element performances are presented over the stable flow operating range of the stage at the design tip speed of 426 m/sec. Stage peak efficiency of 0.83 was obtained at a weight flow of 28.8 kg/sec and a pressure ratio of 1.52. The stall margin for the stage was 8 percent based on weight flow and pressure ratio at peak efficiency and stall. The rotor appears to be stalling prematurely as evidenced by high rotor tip losses.

Urasek, D. C.↗

Assessment of internal mammary artery and saphenous vein graft patency and flow reserve using transthoracic Doppler echocardiography

OBJECTIVE: To investigate transthoracic Doppler echocardiography in the identification of coronary artery bypass graft (CABG) flow for assessing graft patency. DESIGN: The initial study group comprised 45 consecutive patients with previous CABG undergoing elective cardiac catheterisation for recurrent ischaemia. The Doppler variables best correlated with angiographic graft patency were then tested prospectively in a further 84 patients (test group). SETTING: Three tertiary referral centres. INTERVENTIONS: Flow velocities in grafts were recorded at rest and during hyperaemia induced by dipyridamole (0.56 mg/kg/4 min), under the guidance of transthoracic colour Doppler flow mapping. Findings on transthoracic Doppler were compared with angiography. MAIN OUTCOME MEASURES: Feasibility of identifying open grafts by Doppler and diagnostic accuracy for Doppler detection of significant (>/= 70%) graft stenosis. RESULTS: In the test group the identification rate for mammary artery grafts was 100%, for saphenous vein grafts to left anterior descending coronary artery 91%, for vein grafts to right coronary artery 96%, and for vein grafts to circumflex artery 90%. Coronary flow reserve (the ratio between peak diastolic velocity under hyperaemia and at baseline) of < 1.9 (95% confidence interval 1.83 to 2.08) had 100% sensitivity, 98% specificity, 87.5% positive predictive value, and 100% negative predictive value for mammary artery graft stenosis. Coronary flow reserve of < 1.6 (95% CI 1.51 to 1.73) had 91% sensitivity, 87% specificity, 85.4% positive predictive value, and 92.3% negative predictive value for significant vein graft stenosis. CONCLUSIONS: Transthoracic Doppler can provide non-invasive assessment of CABG patency.

Multicenter Study↗