Primary Flow Path Research Articles

On Leakage Flows In A Liquid Hydrogen Multi-Stage Pump for Aircraft Engine Applications

Abstract A comprehensive operational characterisation of a representative, Liquid Hydrogen (LH2) aircraft engine pump is presented in this work. The implications of leakage flows are investigated in a 2-stage, high-pressure pump for a wide range of flow rates and rotational speeds, through 3D (U)RANS simulations. Two configurations are compared: a baseline model comprising the primary flow path components - inducers, impellers and volutes and a realisable pump hardware that includes hub, shroud and power unit cavities. Performance metrics, including head changes and efficiencies, are extracted both at a component and system level. Leakage flow rates of 27.6% and up to 92.9% of the overall pump flow rate are recorded at design and lowest flow points respectively. The head loss in the mid to low flow rates does not exceed 4.5%, but the efficiency diminishes by up to 13.5% at off-design operation. The component analysis indicates significant penalties in impeller efficiency. At high flow rates, the presence of leakage flows improves the overall pump performance by 43% and 27% in head rise and efficiency, due to reduced losses in volutes and connecting ducts. The characterisation of pump behaviour described in this work is crucial for developing and designing safe and predictable LH2 aircraft pumps. Contrary to LH2 pumps utilized in rocketry and for cooling in nuclear industry, aircraft pumps are required to operate with wider turn-down ratios and often at low specific speeds. Therefore, this study addresses design considerations for this enabling technology.

Research on the mechanism and transmission path of the oscillation in DFIG-connected flexible DC transmission system

Concerning the oscillations in DFIG (doubly-fed induction generator)-connected flexible DC transmission system caused by large disturbances, an analysis method based on subsystem partitioning and dynamic energy modeling is put forward. Firstly, the complex high-order system is partitioned into several low-order subsystems, and the dynamic energy model of each subsystem is built using the first integration method. Then, by analyzing the directions and variation patterns of the interaction energy flows between different subsystems in the oscillation process, the redistribution of the primary flow paths of the system’s oscillation energy after large disturbances occur is determined. And then, an oscillation tracing method based on the response characteristics of the oscillation energy links of the subsystems is proposed. In the oscillation process, the interaction energy links between different subsystems are categorized into 'oscillation gaining links', 'oscillation maintaining links' and 'oscillation damping links', according to the contribution of each interaction energy link to the system’s oscillation energy. By analyzing the ‘oscillation gaining links’ corresponding to the dominant interaction energy links of different subsystems, the transmission paths by which the dominant oscillation energy links keeps gaining can be obtained, and the energy source that induces oscillation can be located. Finally, hardware-in-loop tests are conducted on the RTLAB semi-physical simulation platform, the results of which verify the effectiveness of the proposed method. This research paper can provide theoretical support for suppressing the oscillation by adding compensation branches between the subsystems (i.e. damping control strategies based on supplementary compensation branches).

Potential groundwater contamination from oil drilling in the Okavango

Canadian oil and gas company, Reconnaissance Energy Africa (ReconAfrica), plans to exploit potential oil reserves in the Cubango Okavango River Basin (CORB) in Namibia and Botswana, where little is known about the local groundwater systems and how contamination and its effects could impact the region. Using borehole data from target oil and gas areas, we calculated hydraulic gradients, flow direction, and flow velocity to map potential groundwater contamination effects from oil and gas drilling in the region. We also plotted the major geological structures, calculated flowpaths, and estimated travel times for contaminated groundwater that may travel along these preferential flowpaths from the drill sites to the Okavango River and Delta. Results indicate that contaminated groundwater from the oil lease areas could take 3–23.5 years to reach the Okavango River system via the shallow sandy aquifer, but in a worst-case scenario, contamination could reach the Delta within four days via structures associated with dykes and faults that serve as primary flowpaths. Such contamination could adversely affect human health and the region's ecosystems and biodiversity. We recommend prohibiting oil exploration and production activities within the CORB until future studies can determine the impacts of hydrocarbon extraction with greater certainty.

Study of Self Heating Effect in the wake of complete and partial bottom dielectric insertion under 5 nm Stacked Nanosheet Transistor

Beyond 7 nm technology nodes, a Stacked Nanosheet Field Effect Transistor (SNT) is a potential candidate to continue device scaling due to its higher ON/OFF current ratio and less Short Channel Effect than FinFET. However, its gate all-around structure with low thermal conductive HiK oxide material exacerbates the Self-Heating Effect (SHE) issue in the SNT. Further, a dielectric inserted beneath SNT to reduce leakage current disrupts the primary heat flow path through the substrate, aggravating the Self Heating Effect (SHE). This study focuses on exploring the thermal aspect of the SNT with partial and complete dielectric insertion under the device from a performance and reliability perspective. A 3D electrothermal analysis has been performed using Synopsys TCAD tool. The effect of partial and complete dielectric insertion under the device on lattice temperature was investigated in the case of single and double-stack nanosheet transistors. The effects of dielectric insertion on channel and substrate temperature were analyzed. Its variation with width, extension length, interface thermal resistance, and the number of sheets and stack pitch of SNT are discussed. Sheet-dependent heat flow through the substrate and source/drain contact explained for different types of SNT structures and its effect on device performance and reliability discussed. To discuss the effect of multiple parameters variation on the lattice temperature of the device, Analysis of Variance (ANOVA) method was used. It was found that complete dielectric insertion leads to a severe heating issue in the SNT compared to partial and without dielectric insertion.

Fracture Intensity Impacts on Reaction Front Propagation and Mineral Weathering in Three‐Dimensional Fractured Media

AbstractStudying reaction front propagation in heterogeneous natural settings is challenging, but numerical simulations can provide insight into the varying spatial and temporal scales of reaction front propagation. Here, the impact of increasing fracture intensity on mineral dissolution rates, and the extent of reaction front propagation is investigated using reactive transport simulations in upscaled discrete fracture network domains with varied fracture intensity. Domain‐averaged dissolution rates vary less than 0.5 log units regardless of the fracture intensity, but the spatial distribution of reactions is controlled by the location and number of dead‐end fractures and the number of connected flowpaths through the domain. Higher fracture intensities lead to more weathering in the domain because of more available mineral for water‐rock interactions. We find that reaction fronts propagate through the primary flowpaths in the first 10,000 years of the simulation for a 10‐m length domain, then propagate into secondary flowpaths and dead‐end fractures between 10,000 and 100,000 years, and finally into the matrix over timescales of hundreds of thousands of years. The domain‐averaged reaction rates decrease through time corresponding to a transition from dissolution in advection‐dominated, fast‐flowing pathways, to dissolution in transport‐limited zones of disconnected fractures and matrix. Matrix dissolution, or dissolution under transport‐limited conditions, is the dominant process at late times in these simulations. The results of these simulations recreate the observed paradox found in nature where highly fractured hillslopes tend to be more weathered but have slower weathering rates, while hillslopes with fewer fractures, are less weathered but have higher dissolution rates.

Hydro-mechanical coupling characteristics and weakening mechanisms of filling joint resulting from water injection

The injection of fluids into fault gouges in rock formations disturbs the in situ stress conditions, leading to fault slip and increasing the risk of inducing earthquakes. The weakening effect and the permeation of the injected fluid can be influenced significantly by the presence of fault gouges. To investigate this issue, the hydro-mechanical characteristics of fault gouges were evaluated using physical tests to study the combined effects of coupling injecting water and shear deformation. We propose a new experimental apparatus that allows us to measure the spatial distribution of the thickness of a gouge layer sample under combined conditions of shearing and water injection, using 3D scanning technology to evaluate the primary flow path. The test results showed that injecting water had a significant effect in reducing the maximum shear strength, but the degree to which the strength was affected varied with the gouge fill material. The effect of shear deformation is that it will increase the inhomogeneity of the thickness distribution and therefore the distribution of injected water along the fault is not uniformly radial. Additionally, the properties of gouge fill material have an important influence on flow characteristics during fault slipping.

3D Time-Lapse Electrical Resistivity Imaging of Rock Damage Patterns and Gas Flow Paths Resulting from Two Underground Chemical Explosions

Rock damage from underground nuclear explosions (UNEs) has a strong influence on sub-surface gas movement and on seismic waveform characteristics, both of which are used to detect UNEs. Although advanced numerical simulation capabilities exist to predict rock damage patterns and corresponding detection signals, those predictions are dependent on (generally) unknown properties of the host rock. For example, the effects of in-situ mechanical heterogeneities on the explosively generated damage/fractures that provide gas flow pathways to the surface are not well understood, due largely to the difficulty in accessing and characterizing the near-source region. In this paper we demonstrate the emerging use of electrical resistivity tomography (ERT) for imaging rock damage and gas flow patterns resulting from two relatively small-scale underground chemical explosions. Pre-explosion ERT and crosshole seismic imaging revealed a natural fracture zone within the test bed. Post-explosion imaging revealed that the damage zone was non-symmetric and was focused primarily within the pre-existing fracture zone, located 10 m above the first explosion and 5 m above the second explosion. Time-lapse ERT imaging of heated air injected into the detonation borehole revealed the primary gas flow paths to be within the upper margin of the same primary damage zone. These results point to the utility of ERT imaging for understanding rock damage and gas flow patterns under experimental conditions, and to the importance of understanding the effects of geologic heterogeneity on UNE detection signals, particularly gas surface breakthrough times.

Inverse design of microreactor flow fields through anisotropic porous media optimization and dehomogenization

Microreactors are widely used in energy storage and conversion systems, where fluid flow fields play a significant role in reaction-fluid performance. While conventional and bio-inspired flow fields have been investigated in literature, they are often limited to a forward design framework. Forward-designed flow fields depend on an initial layout selection with limited design freedom. In this paper, an inverse design and dehomogenization framework is proposed to discover innovative microreactor flow field designs. A gradient-based method is applied to optimize the spatially varying orientation of anisotropic porous media. A dehomogenization method is adopted to recover the optimized porous media performance by way of intricate microchannel structures in a post-processing step. Through a numerical example, we find trade-offs between the reaction performance and fluid flow performance for multiple optimized microreactor flow fields. Parallel flow fields are confirmed as low pressure drop designs. Branching flow fields with primary and secondary flow paths are discovered for enhanced reaction performance. Optimized designs with balanced reaction-fluid performance can also be obtained with input design requirements. The three-dimensional (3D) performance of each flow field design is verified by additional numerical simulation. The optimized flow fields are shown to outperform conventional parallel and serpentine benchmark designs. The findings of the paper will be useful in novel reactor flow field designs for enhanced performance in biomedical, pharmaceutical, and energy applications.

Study on the proportion and flow path of Yangtze River water diversion into Taihu Lake

Abstract The purpose of this study is to quantify the proportion and flow path of the water diversion from Yangtze River (YRD) into Taihu Lake. Based on the analysis of rainfall and data of Taihu basin in the recent 30 years, a 1-D hydrodynamic model of the main inflow river network area of Taihu basin was constructed, coupled the convection-diffusion model with conservative material, the characteristics of YRD and the water inflow into Taihu Lake (WITL) in three typical years were calculated. The results show that the YRD has shown a significant upward trend in the past 30 years, accounting for 26.4, 35.6 and 42% of the total WITL in three typical years of wet, normal and dry conditions respectively. From the perspective of space, Taige River is the largest river in the western part of the lake that is affected by the river diversion (35%–72%), and Wuxi River is the smallest (1–3%). In addition, the primary flow path of YRD to Taihu Lake was through the Wuyi River and Lake Gehu from the water diversion station west of the Zao River.

Characterizing fluid flow paths in the Hellisheidi geothermal field using detailed fault mapping and stress-dependent permeability

The Husmuli zone of the SW-Iceland Hellisheidi geothermal field is currently being used for re-injection of geothermal fluids and geothermal CO2 for its permanent storage in the form of carbonate minerals. A fully coupled hydro-thermo-mechanical numerical model was employed to investigate the coupled impacts of these complex processes on the calibration of fluid flow paths, which can have significant implications for the long-term performance of this subsurface reservoir.Employing a combination of high-resolution fault mapping with laboratory measurements of stress dependent permeability coupled into a dual porosity field-scale model, the flow paths were calibrated using results of tracer tests performed at the site using stress-dependent permeability tensors. Although vertically extended faults are the primary fluid flow paths, fractures connecting the faults can play an important role in fluid transport. As the upward flow streamlines manifest, deep geological layers can also deviate the fluid flow towards the shallower layers provoking the vertical flow of geothermal fluids. This highlights the sweet spot for sustainable flow and heat extraction in vicinity of faults intercepting the geological layers at depth of 1100 m.It was also shown that the inclusion of the geomechanical calculations in the history matching of the tracer test could lead to changes in arrival time and peak of the tracer profiles. Results of an independent tracer test were used to validate the model and to demonstrate the predictive capability of the calibrated model. This verifies the consistency of our methodology to incorporate the stress-dependent permeability. The results of this comprehensive modelling study provide insight into the likely fluid flow paths, which can have profound impact on the evaluation of various processes such as CO2 mineralisation taking place in Hellisheidi geothermal reservoir.

Efficient Evaluation of Heat Exchangers Behaviour in Nuclear Power Plants

This paper studies the performance of the primary and secondary flow paths of the Heat Exchangers (HEXs) in the Nuclear Power Plants (NPP). Explicit solutions for the wall temperature, primary and secondary wall temperatures of the HEX are presented. Therefore, closed form expressions for wall temperature, primary and secondary wall temperature are proposed. Moreover, performance analysis of HEX is investigated. Curves for the HEX wall temperature are introduced with special emphasis on the initial wall temperature. The various effects of pipe wall heat capacities, perimeter, fluid Heat Transfer (HT) coefficient and wall HT coefficient on the HEX are evaluated. The obtained results confirm that the surface contact between the two fluid passes reduces the temperature to some extent at 5 seconds. Furthermore, it is limited to the pipe length and contact time. As a final point of observation, the distributed temperature through the HEX structure increases with the decrease of perimeter. Since, the wall area decreases with the decrease of perimeter. The proposed work allows more control of the HT process within the HEX. Consequently, the performance of NPPs can be improved.

Geochemical insights into groundwater movement in alpine karst, Bear River Range, Utah, USA

Alpine karst aquifers control the availability and longevity of some water resources, but are not well understood. A conceptual model of the alpine karst aquifer system in the Bear River Range of northern Utah (USA) has been developed by geochemical analysis (major ions, δ18O, δ2H and δ13C values) of seasonal water samples from seven perennial springs, and residence-time assessment (3H and CFCs) of two low- and two high-discharge springs. All spring data can be explained by reaction paths dominated by the dissolution of calcian dolomite. The δ13C values align well with reaction paths for open-system dissolution. Saturation indices and low Ca:Mg molar ratios indicate that incongruent dissolution exerts a strong control on water–rock interactions, complicating interpretation of natural solute tracers. Values of δ18O and δ2H in springs follow the Utah meteoric water line. Snow δ18O values correlate with elevation, but not with increasing rainout distance, providing qualitative estimates of recharge elevation that generally align with previous dye-traces to five of the seven springs. Concentrations of 3H and CFCs likely are best described by binary mixing of subannual recharge with 60–65-year-old groundwater, suggesting that the alpine karst aquifer system in the Bear River Range is best represented by a double-porosity model. Subannual recharge documented by dye traces implies that caverns are the primary flowpaths to the springs, but the presence of decadal-age water may indicate that lower permeability flowpaths dominate during baseflow. No evidence was found for a longer-residing flow component, suggesting high sensitivity to future climate variability.

Stress Effects on Flow and Transport in Three‐Dimensional Fracture Networks

AbstractWe investigate the effects of various external stress regimes on fracture apertures, fluid flow, and solute transport in three‐dimensional fracture networks. We use well‐established geomechanics equations coupled with discrete fracture network modeling to characterize changes in primary flow paths within a complex network as a function of stress magnitude and orientation. These changes manifest in the alterations of the fluid flow field and are measured in terms of Eulerian and Lagrangian flow observables including solute transport, which is a key problem in many hydrologic applications. Changes in primary flow paths affect the solute transport in the network by promoting anomalously early arrival or long tailing behavior. However, early time arrival is not ubiquitous in anisotropically stressed networks and in most cases there is a delayed arrival of solute, which is attributed to (i) the presence of low‐velocity zones normal to the flow direction, (ii) the angle between flow direction and major compressive principal stress directions, and (iii) changes in the primary flow paths (i.e., increases in tortuosity and active network structure). Overall, flows become more channelized in anisotropically stressed fracture networks than in unstressed and isotropically stressed networks.

Multiphase DMA Design Methodology Based on Graph Theory and Many-Objective Optimization

Partitioning a water distribution system (WDS) into district metered areas (DMAs) is a difficult task due to the complex WDS structure and simultaneous consideration of multiple constraints. Further, alternative metrics can define DMA goals, and expressing and quantifying those objectives in an efficient algorithm is also challenging. To address the multifaceted set of objectives, a multiphase DMA design method is developed. DMA feed pipes are identified as the primary flow paths in a branched network. In the methodology presented here, the feed pipe network is first laid out by determining node clusters and boundary pipes that maximize the dissimilarity of pressures and modularities between clusters and minimize the number of cuts between clusters while defining the number of DMAs closest to the desired number. In the methodology’s second phase, secondary DMA feed pipes are identified by minimizing the number of secondary feed pipes and maximizing the nodal excess energy while maintaining desired pressure. Finally, a postoptimization analysis compares the performance of the Pareto solutions based on their availability, water quality, and daily leakage. System availability is calculated based on the minimum cut-set method combined with a new pressure-driven analysis method. To accelerate the optimization algorithm, two strategies are applied: step-by-step optimization and reducing the decision variable searching space by considering desirable DMA characteristics. The effectiveness of the methods is examined by applying it to the C-Town and real B-Town water distribution networks. Results demonstrate that the search space reduction method effectively decomposes the full network into DMAs in the face of multiple hydraulic and water quality metrics.

RELAP5-3D Three-Dimensional Analysis Based on PHÉNIX Dissymmetric Transient Test

Abstract Before the final shutdown of the PHÉNIX fast reactor, the CEA carried out a final set of experimental tests to gather data and additional knowledge on relevant sodium fast reactors (SFR) operation and safety aspects. One of these experiments conducted was the dissymmetrical configuration test, which was selected as benchmark transient within the H2020 SESAME project. ENEA and Sapienza University of Rome are participating in the benchmark using the RELAP5-3D© code. The thermal hydraulic analysis focuses on adequate core cooling prediction in accidental scenario. With the goal of investigating asymmetric thermal hydraulic behavior inside of the reactor pool, two different nodalization approaches have been applied for the RELAP5-3D model, which adopt the same geometrical scheme for the primary flow path, with the exception of the hot and cold pools and the core bypass. The first scheme has been developed using vertical parallel pipes with cross junctions for the hot and cold pools and an equivalent pipe to reproduce the core bypass. The second model includes a multidimensional (MULTID) component, which simulates the pools and provides a detailed nodalization of the core bypass. This study aims at assessing whether the two modeling approaches are equally capable to predict the asymmetrical temperature evolution over the test, caused by the azimuthal asymmetry of the boundary conditions. Blind calculation results are presented and discussed. The paper will be a first step toward the RELAP5-3D code assessment against the experimental results collected as part of the PHÉNIX dissymmetric test.

Post-test simulation of a PLOFA transient test in the CIRCE-HERO facility

CIRCE is a lead–bismuth eutectic alloy (LBE) pool facility aimed to simulate the primary system of a heavy liquid metal (HLM) cooled pool-type fast reactor. The experimental facility was implemented with a new test section, called HERO (Heavy liquid mEtal pRessurized water cOoled tubes), which consists of a steam generator composed of seven double-wall bayonet tubes (DWBT) with an active length of six meters. The experimental campaign aims to investigate HERO behavior, which is representative of the tubes that will compose ALFRED SG. In the framework of the Horizon 2020 SESAME project, a transient test was selected for the realization of a validation benchmark. The test consists of a protected loss of flow accident (PLOFA) simulating the shutdown of primary pumps, the reactor scram and the activation of the DHR system.A RELAP5-3D© nodalization scheme was developed in the pre-test phase at DIAEE of “Sapienza” University of Rome, providing useful information to the experimentalists. The model consisted to a mono-dimensional scheme of the primary flow path and the SG secondary side, and a multi-dimensional component simulating the large LBE pool. The analysis of experimental data, provided by ENEA, has suggested to improve the thermal–hydraulic model with a more detailed nodalization scheme of the secondary loop, looking to reproduce the asymmetries observed on the DWBTs operation. The paper summarizes the post-test activity performed in the frame of the H2020 SESAME project as a contribution of the benchmark activity, highlighting a global agreement between simulations and experiment for all the primary circuit physical quantities monitored. Then, the attention is focused on the secondary system operation, where uncertainties related to the boundary conditions affect the computational results.

Predictive Erosion Model for Mixed Flow Centrifugal Pump

Electrical submersible pumps (ESPs) are widely used in upstream oil production. The presence of a low concentration solid phase, particle-laden flow, in the production fluid may cause severe damage in the internal sections of the pump which reduces its operating lifetime. To better understand the ESP pump's endurance, two different designs of commonly used mixed flow ESPs were studied numerically to determine the pump's flow behavior at its best efficiency point. Computational fluid dynamics (CFD) analysis was conducted on two stages of one design type of pump's primary flow path employing Eulerian–Granular scheme in ANSYS FLUENT. The key parameters affecting the erosion phenomena within the pump such as turbulence kinetic energy, local sand concentration, and near wall relative sand velocity were identified. The predictive erosion model applicable to pumps was developed by correlating the erosion key parameters with available experimental results. It is concluded that the use of an erosion model on the second design of ESP proves the model's versatility to predict the erosion on different designs of ESPs.

Determination of Primary and Secondary Lahar Flow Paths of the Fuego Volcano (Guatemala) Using Morphometric Parameters

On 3 June 2018, a strong eruption of the Fuego volcano in Guatemala produced a dense cloud of 10-km-high volcanic ash and destructive pyroclastic flows that caused nearly 200 deaths and huge economic losses in the region. Subsequently, due to heavy rains, destructive secondary lahars were produced, which were not plotted on the hazard maps using the LAHAR Z software. In this work we propose to complement the mapping of this type of lahars using remote-sensing (Differential Interferometry, DINSAR) in Sentinel images 1A and 2A, to locate areas of deformation of the relief on the flanks of the volcano, areas that are possibly origin of these lahars. To determine the trajectory of the lahars, parameters and morphological indices were analyzed with the software System for Automated Geoscientific Analysis (SAGA). The parameters and morphological indices used were the accumulation of flow (FCC), the topographic wetness index (TWI), the length-magnitude factor of the slope (LS). Finally, a slope stability analysis was performed using the Shallow Landslide Susceptibility software (SHALSTAB) based on the Mohr–Coulomb theory and its parameters: internal soil saturation degree and effective precipitation, parameters required to destabilize a hillside. In this case, the application of this complementary methodology provided a more accurate response of the areas destroyed by primary and secondary lahars in the vicinity of the volcano.

Multiscale simulations of shale gas transport in micro/nano-porous shale matrix considering pore structure influence

Unravelling the transport characteristic of shale gas in actual extraction process is significantly important to improve the gas recovery efficiency and production. But the geological complexity of shale formations in particular makes it challenging to comprehensively cognize the transport mechanisms at multiscale. Traditional shale gas transport model based on circle nanopores system would be failed to capture the realistic transport process of shale gas in the matrix consisting of abundant micro/nano-pores with complex pore structures. Towards this end, herein multiscale simulations were performed to investigate the transport characteristic and mechanism of shale gas in micro/nano-porous shale matrix combining the molecular dynamics (MD) simulations, analytical model and pore network model. MD simulations of shale gas (methane) transport in nanopores demonstrated that the transport behavior is determined by the competition between gas-wall interaction (gas diffusion) and gas-gas intermolecular interaction (viscous flow). Considering different pore structures (e.g., circle, square, triangular, and slit), we proposed a multiscale analytical model with the coupling of continuum flow theory and diffusion effect (Knudsen diffusion and surface diffusion), which is well verified by results from MD simulations and exhibits the practicability of predicting shale gas transport from nanoscale to macroscale. Furthermore, with the aid of pore network model, micro/nano-porous structure was constructed to simulate gas transport in shale matrix. It was found that traditional simulations of shale gas transport in matrix based on circular nanopores would underestimate the transport capacity, for example, that of slit nanopore system is enhanced by as much as 126%. In particular, shale gas transport in shale matrix can be remarkably influenced by the reservoir pressure. At early stage of exploitation, the shale gas transport in matrix is mainly contributed by macropores (H > 100 nm) with the dominant transport characteristic of viscous flow. With the decrease of pressure in continuous exploitation, nanopores (H < 10 nm) will become the primary flow paths due to the enormously enhanced diffusion effect. The intrinsic tendency of our simulated multiscale transport may be useful for understanding the transport behavior and extraction of shale gas in porous shale formations.

An acoustic method for systematic ventricular assist device thrombus evaluation with a novel artificial thrombus model.

Pump thrombosis (PT) is still one of the major adverse events in patients supported with left ventricular assist devices. Nowadays, thrombus detection relies on clinical parameters like reoccurring heart failure symptoms, on changes in pump power consumption, and on laboratory parameters such as increased LDH and hemolysis. Once detected PT is most often persistent and refractory to medical therapy. We therefore designed a novel, non-invasive acoustic method for early pump thrombus detection in an in vitro artificial thrombus model. The study was performed in vitro using a mock circulation loop, artificial blood (water-glycerin) and artificial thrombus material (silicon) allowing for repeatable and defined testing. Tested ventricular assist device (VAD) type was HVAD (Medtronic). Three different thrombus locations were evaluated: on the tilted pad of the rotor, in the primary flow path, and in the secondary flow path beneath the rotor. After evaluating baseline parameters (no thrombus, n=20 for each pump), the influence of thrombi of seven different masses (no thrombus, 0.5-5.0 mg) on pump power consumption and acoustic emission of four HVAD devices was investigated via a microphone system (Sennheiser) and subsequent frequency spectrum analysis (n=12). The acoustic analysis algorithm included the number of frequency peaks recorded. Measurements with thrombi on the tilted pad showed an increased number of frequency peaks with all thrombus sizes compared to baseline measurements without any thrombus (baseline: 32.7±7.4; 0.5 mg: 45.3±10.4 up to 5 mg: 80.4±5.5). Power consumption was relevantly elevated in 5mg thrombus measurement only (6.3±1.29 W compared to 4.9±0.14 W at baseline). Measurements with thrombi in the primary and secondary showed no relevant alteration in power consumption and frequency peak count. We present an acoustic method that detects pump thrombi located on the tilted pad of the HVAD rotor requiring ten times less mass compared to thrombi detected by power consumption alterations used in current detection algorithms. Assuming that pump thrombi are growing over several days, the presented method may detect PT much earlier thereby increasing efficacy of medical therapy and helping to avoid pump exchange.