Ralph W. Metcalfe

Secondary instability of a temporally growing mixing layer

The three-dimensional stability of two-dimensional vortical states of planar mixing layers is studied by direct numerical integration of the Navier-Stokes equations. Small-scale instabilities are shown to exist for spanwise scales at which classical linear modes are stable. These modes grow on convective timescales, extract their energy from the mean flow and exist at moderately low Reynolds numbers. Their growth rates are comparable with the most rapidly growing inviscid instability and with the growth rates of two-dimensional subharmonic (pairing) modes. At high amplitudes, they can evolve into pairs of counter-rotating, streamwise vortices, connecting the primary spanwise vortices, which are very similar to the structures observed in laboratory experiments. The three-dimensional modes do not appear to saturate in quasi-steady states as do the purely two-dimensional fundamental and subharmonic modes in the absence of pairing. The subsequent evolution of the flow depends on the relative amplitudes of the pairing modes. Persistent pairings can inhibit threedimensional instability and, hence, keep the flow predominantly two-dimensional. Conversely, suppression of the pairing process can drive the three-dimensional modes to more chaotic, turbulent-like states. An analysis of high-resolution simulations of fully turbulent mixing layers confirms the existence of rib-like structures and that their coherence depends strongly on the presence of the two-dimensional pairing modes.

Direct numerical simulations of homogeneous turbulence in density‐stratified fluids

We report on the results of direct numerical simulations of the decay of homogeneous turbulence in density‐stratified fluids. The effects of stratification are examined by performing simulations for a range of different initial Froude numbers. Also two theoretical approaches are presented, which are applicable in the later stages of decay. For the conditions computed, although the stritification introduced definite wave‐like characteristics into the flow‐fields, the non‐linearity, as expressed for example by the spectral energy transfer, remained strong. However, the selective suppression of certain types of transfer was observed to occur.

Effects of heat release on the large-scale structure in turbulent mixing layers

The effects of chemical heat release on the large-scale structure in a chemically reacting, turbulent mixing layer are investigated using direct numerical simulations. Three-dimensional, time-dependent simulations are performed for a binary, single-step chemical reaction occurring across a temporally developing turbulent mixing layer. It is found that moderate heat release slows the development of the large-scale structures and shifts their wavelengths to larger scales. The resulting entrainment of reactants is reduced, decreasing the overall chemical product formation rate. The simulation results are interpreted in terms of turbulence energetics, vorticity dynamics, and stability theory. The baroclinic torque and thermal expansion in the mixing layer produce changes in the flame vortex structure that result in more diffuse vortices than in the constant-density case, resulting in lower rotation rates of the large-scale structures. Previously unexplained anomalies observed in the mean velocity profiles of reacting jets and mixing layers are shown to result from vorticity generation by baroclinic torques.

Direct numerical simulations of chemically reacting turbulent mixing layers

The results of direct numerical simulations of chemically reacting, turbulent mixing layers are presented. The reaction considered is a binary, irreversible reaction with no heat release, so that only the effect of the turbulence on the chemical reaction is investigated. The simulation results are shown to be consistent with similarity theory, and are found to be in approximate agreement with laboratory data, even though there are no adjustable parameters in the method.

Reproducibility of haemodynamical simulations in a subject-specific stented aneurysm model - A report on the Virtual Intracranial Stenting Challenge 2007

This paper presents the results of the Virtual Intracranial Stenting Challenge (VISC) 2007, an international initiative whose aim was to establish the reproducibility of state-of-the-art haemodynamical simulation techniques in subject-specific stented models of intracranial aneurysms (IAs). IAs are pathological dilatations of the cerebral artery walls, which are associated with high mortality and morbidity rates due to subarachnoid haemorrhage following rupture. The deployment of a stent as flow diverter has recently been indicated as a promising treatment option, which has the potential to protect the aneurysm by reducing the action of haemodynamical forces and facilitating aneurysm thrombosis. The direct assessment of changes in aneurysm haemodynamics after stent deployment is hampered by limitations in existing imaging techniques and currently requires resorting to numerical simulations. Numerical simulations also have the potential to assist in the personalized selection of an optimal stent design prior to intervention. However, from the current literature it is difficult to assess the level of technological advancement and the reproducibility of haemodynamical predictions in stented patient-specific models. The VISC 2007 initiative engaged in the development of a multicentre-controlled benchmark to analyse differences induced by diverse grid generation and computational fluid dynamics (CFD) technologies. The challenge also represented an opportunity to provide a survey of available technologies currently adopted by international teams from both academic and industrial institutions for constructing computational models of stented aneurysms. The results demonstrate the ability of current strategies in consistently quantifying the performance of three commercial intracranial stents, and contribute to reinforce the confidence in haemodynamical simulation, thus taking a step forward towards the introduction of simulation tools to support diagnostics and interventional planning.

A new mechanism of small-scale transition in a plane mixing layer : core dynamics of spanwise vortices

We present a new mechanism of small-scale transition via core dynamics instability (CDI) in an incompressible plane mixing layer, a transition which is not reliant on the presence of longitudinal vortices (‘ribs’) and which can originate much earlier than ribinduced transition. Both linear stability analysis and direct numerical simulation are used to describe CDI growth and subsequent transition in terms of vortex dynamics and vortex line topology. CDI is characterized by amplifying oscillations of core size non-uniformity and meridional flow within spanwise vortices (‘rolls’), produced by a coupling of roll swirl and meridional flow that is manifested by helical twisting and untwisting of roll vortex lines. We find that energetic CDI is excited by subharmonic oblique modes of shear layer instability after roll pairing, when adjacent rolls with out-of-phase undulations merge. Starting from moderate initial disturbance amplitudes, twisting of roll vortex lines generates within the paired roll opposing spanwise flows which even exceed the free-stream velocity. These flows collide to form a nearly irrotational bubble surrounded by a thin vorticity sheath of a large diameter, accompanied by folding and reconnection of roll vortex lines and local transition. We find that accelerated energy transfer to high wavenumbers precedes the development of roll internal intermittency; this transfer, inferred from increased energy at high wavenumbers and an intensification of roll vorticity, occurs prior to the development of strong opposite-signed (to the mean) spanwise vorticity and granularity of the roll vorticity distribution. We demonstrate that these core dynamics are not reliant upon special symmetries and also occur in the presence of moderate-strength ribs, despite entrapment of ribs within pairing rolls. In fact, the roll vorticity dynamics are dominated by CDI if ribs are not sufficiently strong to first initiate transition; thus CDI may govern small-scale transition for moderate initial 3D disturbances, typical of practical situations. Results suggest that CDI constitutes a new generic mechanism for transition to turbulence in shear flows.

Stabilization of laminar boundary layers by compliant membranes

Direct numerical simulations have been performed to study boundary layer flow over a compliant membrane. While coupling directly to the fluid dynamics code through the wall boundary condition, this membrane model displays many of the important features of other compliant coatings. While membrane parameters are identified that increase the critical Reynolds number by a factor of 2 compared with the rigid wall, this is mitigated by a significant increase in the growth rates in the unstable region. A detailed analysis of the temporal evolution and spatial structure of the terms in the kinetic energy balance equation shows significant differences in the behavior of two different classes of modes of the system: class A waves, which are destabilized by increased dissipation in the membrane, and class B waves, which are stabilized by membrane dissipation. For class A waves, the dominant dissipative mechanism is viscous damping in the fluid augmented by negative energy production, the principal stabilization mech...

Digital thermal monitoring of vascular function: a novel tool to improve cardiovascular risk assessment

Abstract Digital thermal monitoring (DTM) of vascular function during cuff-occlusive reactive hyperemia relies on the premise that changes in fingertip temperature during and after an ischemic stimulus reflect changes in blood flow. To determine its utility in individuals with and without known coronary heart disease (CHD), 133 consecutive individuals (age 54 ± 10 years, 50% male, 19 with known CHD) underwent DTM during and after 2 minutes of supra-systolic arm cuff inflation. Fingertip temperatures of the occluded and non-occluded fingertips were measured simultaneously. Post-cuff deflation temperature rebound (TR) was lower in the CHD patients and in those with an increased Framingham risk score (FRS) compared to the normal group. After adjustment for age, sex, and cardiac risk factors, TR was significantly lower in those with CHD compared to those without CHD (p < 0.05). This study demonstrates that vascular dysfunction measured by DTM is associated with CHD and an increased FRS, and could potentially be used to identify high-risk patients.

Spatial Resolution Limits of Multislice Computed Tomography (MS-CT), C-arm-CT, and Flat Panel-CT (FP-CT) Compared to MicroCT for Visualization of a Small Metallic Stent

RATIONALE AND OBJECTIVES Small metallic stents are increasingly used in the treatment of cerebral aneurysms and for revascularization in ischemic strokes. Realistic three-dimensional datasets of a stent were obtained by using three x-ray-based imaging methods in current clinical use. Multislice-CT (MS-CT), C-arm flat detector-CT (C-arm CT, ACT), and flat panel-CT (FP-CT) were compared with high-resolution laboratory MicroCT scans that served as a reference standard. The purpose was to assess and compare the quality and accuracy of current clinical three-dimensional reconstructions of a vascular stents. MATERIAL & METHODS A 3 × 20 mm Cypher stent was deployed in a straight polytetrafluoroethylene tube and filled with nondiluted iodine contrast and BaSO(4). MS-CT images of the static tube phantom and stent were acquired using GE LightSpeed VCT Series, C-arm CT images were obtained using Artis (DynaCT, Siemens), FP-CT were obtained using a preclinical research CT (GE), and MicroCT images were obtained using eXplore Locus SP (GE). DICOM datasets were analyzed using Amira and Matlab. RESULTS Because of blooming effects, the maximum intensity projections (MIPs) and volume renderings generated from MS-CT showed significantly increased strut dimensions with no distinction between the regular struts and connector struts while the lumen diameter is artificially reduced. The shape of the reconstructed stent surface differed remarkably from the real stent. C-arm CT and FP-CT volume renderings more accurately represented the struts. Consistently capturing the structure of the connectors and the strut shape definition was highly threshold dependent. The stent lumen was about 30% underestimated by MS-CT when compared to MicroCT. CONCLUSION The spatial resolution of current clinical CT for imaging of small metallic stents is insufficient to visualize fine geometrical details. Further improvement in the spatial resolution of clinical imaging technologies combined with better software and hardware for image postprocessing will be necessary for detailed structural analysis, evaluation of the stent lumen in vivo, and to permit accurate assessment of stent patency and early detection potential in-stent stenosis.

Preload sensitivity of the Jarvik 2000 and HeartMate II left ventricular assist devices.

The in vitro sensitivity of continuous flow pumps to preload and afterload pressure has been well characterized. We compared flow in the Jarvik 2000 and HeartMate II continuous flow left ventricular assist devices (LVADs) at different inflow and outflow pressures and different pump speeds. This allowed us to measure the impact of a changing inflow pressure on the pump flow rate at different speeds but against a constant afterload. The resulting preload sensitivity curves showed that, overall, both LVADs have a mean preload sensitivity of 0.07 L/min/mm Hg in the physiologic ranges of pressures and flows encountered during normal operation. The HeartMate II pump had an increased preload sensitivity (up to approximately 0.1 L/min/mm Hg) as the preload was increased. The preload sensitivity of the Jarvik 2000 LVAD was more variable, having several peaks and troughs as the preload was increased. In future LVADs, improved preload sensitivity may allow passive regulation of pump output, optimize ventricular unloading, and decrease the risk of ventricular suction by the pump.