Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “root”

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 937 records · Page 52

Assessing horizontal gene transfer in the rhizosphere of Brachypodium distachyon using fabricated ecosystems (EcoFABs)

ABSTRACT Horizontal gene transfer (HGT) is a major process by which genes are transferred between microbes in the rhizosphere. However, examining HGT remains challenging due to the complexity of mimicking conditions within the rhizosphere. Fabricated ecosystems (EcoFABs) have been used to investigate several complex processes in plant-associated environments. Here we show that EcoFABs are efficient tools to examine and measure HGT frequency in the rhizosphere. We provide the first demonstration of gene transfer via a triparental conjugation system in the Brachypodium distachyon rhizosphere in an EcoFAB using Pseudomonas putida KT2440 as both donor and recipient bacterial strain with the donor containing a mobilizable and non-self-transmissible plasmid. We observed that the frequency of plasmid transfer in the rhizosphere is potentially dependent on the plant developmental stage and the composition and amount of root exudates. The frequency of plasmid transfer also increased with higher numbers of donor cells. We demonstrate the transfer of plasmid from P. putida to another B. distachyon root colonizer, Burkholderia sp. OAS925, showing HGT within a rhizosphere microbial community. Environmental stresses also influenced the rate and efficiency of HGT in the rhizosphere between different species and genera. This study provides a robust workflow to evaluate transfer of engineered plasmids in the rhizosphere when such plasmids are potentially introduced in a field or other plant-associated environments. IMPORTANCE We report the use of EcoFABs to investigate the HGT process in a rhizosphere environment. It highlights the potential of EcoFABs in recapitulating the dynamic rhizosphere conditions as well as their versatility in studying plant-microbe interactions. This study also emphasizes the importance of studying the parameters impacting the HGT frequency. Several factors such as plant developmental stages, nutrient conditions, number of donor cells, and environmental stresses influence gene transfer within the rhizosphere microbial community. This study paves the way for future investigations into understanding the fate and movement of engineered plasmids in a field environment.

Priya, Shweta↗

Global Diversity and Biogeography of the Zostera marina Mycobiome

Seagrasses are marine flowering plants that provide critical ecosystem services in coastal environments worldwide. Marine fungi are often overlooked in microbiome and seagrass studies, despite terrestrial fungi having critical functional roles as decomposers, pathogens, or endophytes in global ecosystems. Here, we characterize the distribution of fungi associated with the seagrass Zostera marina, using leaves, roots, and rhizosphere sediment from 16 locations across its full biogeographic range. Using high-throughput sequencing of the ribosomal internal transcribed spacer (ITS) region and 18S rRNA gene, we first measured fungal community composition and diversity. We then tested hypotheses of neutral community assembly theory and the degree to which deviations suggested that amplicon sequence variants (ASVs) were plant selected or dispersal limited. Finally, we identified a core mycobiome and investigated the global distribution of differentially abundant ASVs. We found that the fungal community is significantly different between sites and that the leaf mycobiome follows a weak but significant pattern of distance decay in the Pacific Ocean. Generally, there was evidence for both deterministic and stochastic factors contributing to community assembly of the mycobiome, with most taxa assembling through stochastic processes. The Z. marina core leaf and root mycobiomes were dominated by unclassified Sordariomycetes spp., unclassified Chytridiomycota lineages (including Lobulomycetaceae spp.), unclassified Capnodiales spp., and Saccharomyces sp. It is clear from the many unclassified fungal ASVs and fungal functional guilds that knowledge of marine fungi is still rudimentary. Further studies characterizing seagrass-associated fungi are needed to understand the roles of these microorganisms generally and when associated with seagrasses. Fungi have important functional roles when associated with land plants, yet very little is known about the roles of fungi associated with marine plants, like seagrasses. In this study, we report the results of a global effort to characterize the fungi associated with the seagrass Zostera marina across its full biogeographic range. Although we defined a putative global core fungal community, it is apparent from the many fungal sequences and predicted functional guilds that had no matches to existing databases that general knowledge of seagrass-associated fungi and marine fungi is lacking. This work serves as an important foundational step toward future work investigating the functional ramifications of fungi in the marine ecosystem.

ITS2↗

Comparative genomics and stable isotope analysis reveal the saprotrophic-pathogenic lifestyle of a neotropical fungus

In terrestrial forested ecosystems, fungi may interact with trees in at least three distinct ways: (i) associated with roots as symbionts; (ii) as pathogens in roots, trunks, leaves, flowers, and fruits; or (iii) decomposing dead tree tissues on soil or even on dead tissues in living trees. Distinguishing the latter two nutrition modes is rather difficult in Hymenochaetaceae (Basidiomycota) species. Herein, we have used an integrative approach of comparative genomics, stable isotopes, host tree association, and bioclimatic data to investigate the lifestyle ecology of the scarcely known neotropical genus Phellinotus, focusing on the unique species Phellinotus piptadeniae. This species is strongly associated with living Piptadenia gonoacantha (Fabaceae) trees in the Atlantic Forest domain on a relatively high precipitation gradient. Phylogenomics resolved P. piptadeniae in a clade that also includes both plant pathogens and typical wood saprotrophs. Furthermore, both genome-predicted Carbohydrate-Active Enzymes (CAZy) and stable isotopes (δ 13 C and δ 15 N) revealed a rather flexible lifestyle for the species. Altogether, our findings suggest that P. piptadeniae has been undergoing a pathotrophic specialization in a particular tree species while maintaining all the metabolic repertoire of a wood saprothroph.

59 BASIC BIOLOGICAL SCIENCES↗

Fungal-Bacterial Cooccurrence Patterns Differ between Arbuscular Mycorrhizal Fungi and Nonmycorrhizal Fungi across Soil Niches

Soil bacteria and fungi are known to form niche-specific communities that differ between actively growing and decaying roots. Yet almost nothing is known about the cross-kingdom interactions that frame these communities and the environmental filtering that defines these potentially friendly or competing neighbors. We explored the temporal and spatial patterns of soil fungal (mycorrhizal and nonmycorrhizal) and bacterial cooccurrence near roots of wild oat grass, Avena fatua, growing in its naturalized soil in a greenhouse experiment. Amplicon sequences of the fungal internal transcribed spacer (ITS) and bacterial 16S rRNA genes from rhizosphere and bulk soils collected at multiple plant growth stages were used to construct covariation-based networks as a step toward identifying fungal-bacterial associations. Corresponding stable-isotope-enabled metagenome-assembled genomes (MAGs) of bacteria identified in cooccurrence networks were used to inform potential mechanisms underlying the observed links. Bacterial-fungal networks were significantly different in rhizosphere versus bulk soils and between arbuscular mycorrhizal fungi (AMF) and nonmycorrhizal fungi. Over 12 weeks of plant growth, nonmycorrhizal fungi formed increasingly complex networks with bacteria in rhizosphere soils, while AMF more frequently formed networks with bacteria in bulk soils. Analysis of network-associated bacterial MAGs suggests that some of the fungal-bacterial links that we identified are potential indicators of bacterial breakdown and consumption of fungal biomass, while others intimate shared ecological niches.

54 ENVIRONMENTAL SCIENCES↗

Ectomycorrhizal fungi alter soil food webs and the functional potential of bacterial communities

Most of Earth’s trees rely on critical soil nutrients that ectomycorrhizal fungi (EcMF) liberate and provide, and all of Earth’s land plants associate with bacteria that help them survive in nature. Yet, our understanding of how the presence of EcMF modifies soil bacterial communities, soil food webs, and root chemistry requires direct experimental evidence to comprehend the effects that EcMF may generate in the belowground plant microbiome. To this end, we grew Pinus muricata plants in soils that were either inoculated with EcMF and native forest bacterial communities or only native bacterial communities. We then profiled the soil bacterial communities, applied metabolomics and lipidomics, and linked omics data sets to understand how the presence of EcMF modifies belowground biogeochemistry, bacterial community structure, and their functional potential. We found that the presence of EcMF (i) enriches soil bacteria linked to enhanced plant growth in nature, (ii) alters the quantity and composition of lipid and non-lipid soil metabolites, and (iii) modifies plant root chemistry toward pathogen suppression, enzymatic conservation, and reactive oxygen species scavenging. Using this multi-omic approach, we therefore show that this widespread fungal symbiosis may be a common factor for structuring soil food webs.

59 BASIC BIOLOGICAL SCIENCES↗

Nutrient limitation shapes functional traits of mycorrhizal fungi and phosphorus-cycling bacteria across an elevation gradient

In nutrient-limited high-elevation ecosystems, plants rely on arbuscular mycorrhizal (AM) fungi to provide mineral phosphorus (P) in the form of phosphate (PO43-). AM fungi gather these nutrients from phosphorus-cycling bacteria (PCBs) that can mineralize PO43- from organic matter and solubilize mineral-bound P. How climate, soil factors, and nutrient limitation influence AM fungi and PCB assembly remains unclear. We collected soil from montane meadows across a 1,000-m elevation gradient on three replicate mountainsides and analyzed AM fungal marker genes, P-cycling genes from shotgun metagenomes, and edaphic measurements. High-elevation soils had nearly 50-fold less soil PO₄³⁻ and 60% more AM fungal hyphae than low-elevation soils. AM fungal turnover was linked to changes in pH, organic carbon, and PO₄³-. The composition of 198 P-cycling genes was influenced by the AM fungal community structure. Drivers of individual PCB functional genes, including pH and organic carbon, varied with gene phylogeny. We found a trade-off in P-cycling strategies across elevation: P-rich, low-elevation soils supported root-colonizing AM fungi and organic P-mineralizing bacteria. P-poor, high-elevation soils were dominated by stress-tolerant AM fungi and mineral P-solubilizing bacteria. Our results suggest that AM fungi and PCB community turnover across elevation are both shaped by pH, organic carbon, and P availability. With continued climate warming, the structure and function of mountaintop ecosystems might shift to resemble lower elevations, disrupting long-established and specialized microbial assemblages, with consequences for P-cycling dynamics and the total P available to plant communities.IMPORTANCEPhosphorus (P) limits plant productivity in high-elevation ecosystems, yet the microbial networks that mobilize P, including arbuscular mycorrhizal (AM) fungi and phosphorus-cycling bacteria (PCBs), remain under-characterized in these nutrient-poor soils. We show that across a 10,00-m elevation gradient, AM fungi and P-cycling gene assemblages shift predictably with pH, organic carbon, and phosphate availability. Higher elevations, with less available P, select for stress-tolerant AM fungal taxa and PCB strategies geared toward mineral solubilization, while low-elevation sites favor root colonization by AM fungi and organic P mineralization. These results suggest that nutrient limitation can constrain microbial community assembly in consistent ways across landscapes. High mountain soils are low in P and rely on a network of underground AM fungi and PCB to deliver nutrients to plants. This study shows how those underground relationships reorganize with elevation and how climate change could collapse long-standing microbial strategies by pushing high-elevation ecosystems toward lowland conditions. As soils warm and dry, the microbial scaffolding that supports alpine plant life may become increasingly unstable.

arbuscular mycorrhizal fungi↗

Agricultural Management Affects the Active Rhizosphere Bacterial Community Composition and Nitrification

Cropping system diversity provides yield benefits that may result from shifts in the composition of root-associated bacterial and fungal communities, which either enhance nutrient availability or limit nutrient loss. We investigated whether temporal diversity of annual cropping systems (four versus two crops in rotation) influences the composition and metabolic activities of root-associated microbial communities in maize at a developmental stage when the peak rate of nitrogen uptake occurs. We monitored total (DNA-based) and potentially active (RNA-based) bacterial communities and total (DNA-based) fungal communities in the soil, rhizosphere, and endosphere. Cropping system diversity strongly influenced the composition of the soil microbial communities, which influenced the recruitment of the resident microbial communities and, in particular, the potentially active rhizosphere and endosphere bacterial communities. The diversified cropping system rhizosphere recruited a more diverse bacterial community (species richness), even though there was little difference in soil species richness between the two cropping systems. In contrast, fungal species richness was greater in the conventional rhizosphere, which was enriched in fungal pathogens; the diversified rhizosphere, however, was enriched in Glomeromycetes. While cropping system influenced endosphere community composition, greater correspondence between DNA- and RNA-based profiles suggests a higher representation of active bacterial populations. Cropping system diversity influenced the composition of ammonia oxidizers, which coincided with diminished potential nitrification activity and gross nitrate production rates, particularly in the rhizosphere. The results of our study suggest that diversified cropping systems shift the composition of the rhizosphere’s active bacterial and total fungal communities, resulting in tighter coupling between plants and microbial processes that influence nitrogen acquisition and retention.

60 APPLIED LIFE SCIENCES↗

Visualizing the Hidden Half: Plant-Microbe Interactions in the Rhizosphere

Plant roots and the associated rhizosphere constitute a dynamic environment that fosters numerous intra- and interkingdom interactions, including metabolite exchange between plants and soil mediated by root exudates and the rhizosphere microbiome. These interactions affect plant fitness and performance, soil health, and the belowground carbon budget.

59 BASIC BIOLOGICAL SCIENCES↗

Challenging Practices of Algebraic Battery Life Models through Statistical Validation and Model Identification via Machine-Learning

Various modeling techniques are used to predict the capacity fade of Li-ion batteries. Algebraic reduced-order models, which are inherently interpretable and computationally fast, are ideal for use in battery controllers, technoeconomic models, and multi-objective optimizations. For Li-ion batteries with graphite anodes, solid-electrolyte-interphase (SEI) growth on the graphite surface dominates fade. This fade is often modeled using physically informed equations, such as square-root of time for predicting solvent-diffusion limited SEI growth, and Arrhenius and Tafel-like equations predicting the temperature and state-of-charge rate dependencies. In some cases, completely empirical relationships are proposed. However, statistical validation is rarely conducted to evaluate model optimality, and only a handful of possible models are usually investigated. This article demonstrates a novel procedure for automatically identifying reduced-order degradation models from millions of algorithmically generated equations via bi-level optimization and symbolic regression. Identified models are statistically validated using cross-validation, sensitivity analysis, and uncertainty quantification via bootstrapping. On a LiFePO 4 /Graphite cell calendar aging data set, automatically identified models utilizing square-root, power law, stretched exponential, and sigmoidal functions result in greater accuracy and lower uncertainty than models identified by human experts, and demonstrate that previously known physical relationships can be empirically "rediscovered" using machine learning.

25 ENERGY STORAGE↗

Resilience Measurement Framework For Post-deployment Artificial Intelligence (ai) Integrated Systems

Resilience is largely defined as the ability to adapt or recover from adverse conditions, stresses, attacks, or compromises on systems that use or are enabled by digital resources. In Artificial Intelligence Management and Research for Advanced Networked Testbed Hub (AMARANTH), resilience is measured in the amount of time it took from the beginning of a testing period for the model to reach predictions outside of the original 95% confidence interval or using the Kullback-Leibler (KL) divergence theorem, the Population Stability Index (PSI), and traditional methods such as root mean squared error (RMSE) threshold. Artificial Intelligence (AI) model drift is of significant concern when deploying AI-integrated systems into critical and/or secure environments. Drift can impact resilience of the AI-integrated system post-deployment and requires consistent maintenance and upkeep to ensure the model is accurate and precise. To quantify model drift and predict the point when a model's drift becomes unacceptable, we describe using Kullback-Leibler (KL) divergence, Population Stability Index (PSI) and/or confidence interval width estimations to determine the point of failure and time to failure of a model post-deployment. Through simple code functions, the KL-divergence, PSI, confidence interval, and root mean squared (RMSE) point of failures can be used to derive when a model needs to be maintained as well as the impact of adversarial action through statistical means.

Yockey, Patience [Idaho National Laboratory (INL),↗

Reducing Southern Ocean Shortwave Radiation Errors in the ERA5 Reanalysis with Machine Learning and 25 Years of Surface Observations

Earth system models struggle to simulate clouds and their radiative effects over the Southern Ocean, partly due to a lack of measurements and targeted cloud microphysics knowledge. We have evaluated biases of downwelling shortwave radiation in the ERA5 climate reanalysis using 25 years (1995–2019) of summertime surface measurements, collected on the Research and Supply Vessel (RSV) Aurora Australis, the Research Vessel (R/V) Investigator, and at Macquarie Island. During October–March daylight hours, the ERA5 simulation of SW down exhibited large errors (mean bias = 54 W m -2 , mean absolute error = 82 W m -2 , root-mean-square error = 132 W m -2 , and R 2 = 0.71). To determine whether we could improve these statistics, we bypassed ERA5’s radiative transfer model for SW down with machine learning–based models using a number of ERA5’s gridscale meteorological variables as predictors. These models were trained and tested with the surface measurements of SW down using a 10-fold shuffle split. An extreme gradient boosting (XGBoost) and a random forest–based model setup had the best performance relative to ERA5, both with a near complete reduction of the mean bias error, a decrease in the mean absolute error and root-mean-square error by 25% ± 3%, and an increase in the R 2 value of 5% ± 1% over the 10 splits. Large improvements occurred at higher latitudes and cyclone cold sectors, where ERA5 performed most poorly. We further interpret our methods using Shapley additive explanations. Our results indicate that data-driven techniques could have an important role in simulating surface radiation fluxes and in improving reanalysis products.

54 ENVIRONMENTAL SCIENCES↗

The Influence of Shear on Deep Convection Initiation. Part I: Theory

Abstract This article introduces a novel hypothesis for the role of vertical wind shear (“shear”) in deep convection initiation (DCI). In this hypothesis, initial moist updrafts that exceed a width and shear threshold will “root” within a progressively deeper steering current with time, increase their low-level cloud-relative flow and inflow, widen, and subsequently reduce their susceptibility to entrainment-driven dilution, evolving toward a quasi-steady self-sustaining state. In contrast, initial updrafts that do not exceed the aforementioned thresholds experience suppressed growth by shear-induced downward pressure gradient accelerations, will not root in a deep-enough steering current to increase their inflow, will narrow with time, and will succumb to entrainment-driven dilution. In the latter case, an externally driven lifting mechanism is required to sustain deep convection, and deep convection will not persist in the absence of such lifting mechanism. A theoretical model is developed from the equations of motion to further explore this hypothesis. The model indicates that shear generally suppresses DCI, raising the initial subcloud updraft width that is necessary for it to occur. However, there is a pronounced bifurcation in updraft growth in the model after the onset of convection. Sufficiently wide initial updrafts grow and eventually achieve a steady state. In contrast, insufficiently wide initial updrafts shrink with time and eventually decay completely without external support. A sharp initial updraft radius threshold discriminates between these two outcomes. Thus, consistent with our hypothesis and observations, shear inhibits DCI in some situations, but facilitates it in others.

42 ENGINEERING↗

Data for Influence of Particle Size on NIR Spectroscopic Characterization of Sorghum Biomass for the Biofuel Industry

NIR spectroscopy is a rapid and accurate green technology for high-throughput biomass characterization, including sorghum ( Sorghum bicolor ), a promising energy crop for the biofuel industry. This study assessed the influence of particle size on NIR spectroscopic analysis (wavelength range: 867–2535 nm) of sorghum biomass composition. Grown under field conditions, a total of 113 types of genetically diverse sorghum accessions were dried, ground, and sieved (<250, 250–600, 600–850, and > 850 µm particle size) for developing partial least square regression (PLSR) prediction models for moisture, ash, extractive, glucan, xylan, acid-soluble lignin (ASL), acid-insoluble lignin (AIL), and total lignin (ASL + AIL). Overall, smaller particle sizes provided better model performance, while no single particle size provided the best performance for all the selected components. With only 9 selected bands and 4 latent variables (LVs), the best PLSR model was obtained for moisture with particle size of 600–850 µm with the square root of the coefficient of determination (R) of 0.85, the ratio of prediction to deviation (RPD) of 2.2, and the root mean square error (RMSE) of 0.46 % in external validation. Similar model performances were also obtained for ash, extractive, glucan, and xylan. This study showed that size reduction could effectively improve NIR spectroscopic analysis for lipid-producing sorghum biomass for the biofuel industry.

Biomass Analytics↗

Data for The Stem Cell-Type Transcriptome of Bioenergy Sorghum Reveals the Spatial Regulation of Secondary Cell Wall Networks

Bioenergy sorghum is a low-input, drought-resilient, deep-rooting annual crop that has high biomass yield potential enabling the sustainable production of biofuels, biopower, and bioproducts. Bioenergy sorghum’s 4-5 m stems account for ~80% of the harvested biomass. Stems accumulate high levels of sucrose that could be used to synthesize bioethanol and useful biopolymers if information about stem cell-type gene expression and regulation was available to enable engineering. To obtain this information, Laser Capture Microdissection (LCM) was used to isolate and collect transcriptome profiles from five major cell types that are present in stems of the sweet sorghum Wray. Transcriptome analysis identified genes with cell-type specific and cell-preferred expression patterns that reflect the distinct metabolic, transport, and regulatory functions of each cell type. Analysis of cell-type specific gene regulatory networks (GRNs) revealed that unique TF families contribute to distinct regulatory landscapes, where regulation is organized through various modes and identifiable network motifs. Cell-specific transcriptome data was combined with a stem developmental transcriptome dataset to identify the GRN that differentially activates the secondary cell wall (SCW) formation in stem xylem sclerenchyma and epidermal cells. The cell-type transcriptomic dataset provides a valuable source of information about the function of sorghum stem cell types and GRNs that will enable the engineering of bioenergy sorghum stems.

Software↗

Data for Impacts of Legacy and Contemporary Nitrogen Inputs on N2O and CO2 Emissions in Miscanthus and Maize Cultivated Soils

Nutrient inputs influence the sustainability of bioenergy crop production through contemporary (shortly after addition) and legacy effects (persisting over years) on microbial nitrogen (N) and carbon cycling, which contribute to greenhouse gas emissions. However, the relative importance of contemporary and legacy effects and how that could vary by crop functional types is poorly understood. Considering its rhizomatous roots and perennial growth, we hypothesized that Miscanthus × giganteu s ( M × g ) would be more sensitive to legacy N fertilization and the historical context of its environment than an annual crop like maize. To test this hypothesis, we examined the effects of legacy and contemporary N inputs on nitrous oxide (N2O) and carbon dioxide (CO2) emissions, as well as key N cycling genes in soils where M × g and maize were grown. A 150-day soil incubation experiment was conducted using soils from a long-term M × g and maize fertility experiment with three historic N fertilization rates (0, 112, and 336 kg N ha−1 year−1) and a contemporary amendment (60 mg N kg−1) with negative control (0 mg N kg−1). We observed significant increases in cumulative N2O emissions in M × g soils relative to maize soils, particularly at higher legacy fertilization rates, while contemporary N had no significant effect. Bacterial amoA gene abundance, which plays a significant role in nitrification in nutrient-rich soils, also increased with higher legacy fertilization rates in M × g soils but was unaffected by the contemporary N. In maize soils, legacy and contemporary N did not significantly affect N2O emissions, but cumulative CO2 emissions and amoA gene abundance significantly increased. The abundances of norB genes were not significantly influenced by either legacy fertilization or contemporary N amendments in either soil. Our findings demonstrate the greater importance of fertilization history over contemporary N in mediating soil N2O emissions, particularly for perennial bioenergy crops.

Carbon↗

A glycan receptor kinase facilitates intracellular accommodation of arbuscular mycorrhiza and symbiotic rhizobia in the legume Lotus japonicus

Receptors that distinguish the multitude of microbes surrounding plants in the environment enable dynamic responses to the biotic and abiotic conditions encountered. In this study, we identify and characterise a glycan receptor kinase, EPR3a, closely related to the exopolysaccharide receptor EPR3. Epr3a is up-regulated in roots colonised by arbuscular mycorrhizal (AM) fungi and is able to bind glucans with a branching pattern characteristic of surface-exposed fungal glucans. Expression studies with cellular resolution show localised activation of the Epr3a promoter in cortical root cells containing arbuscules. Fungal infection and intracellular arbuscule formation are reduced in epr3a mutants. In vitro , the EPR3a ectodomain binds cell wall glucans in affinity gel electrophoresis assays. In microscale thermophoresis (MST) assays, rhizobial exopolysaccharide binding is detected with affinities comparable to those observed for EPR3, and both EPR3a and EPR3 bind a well-defined β-1,3/β-1,6 decasaccharide derived from exopolysaccharides of endophytic and pathogenic fungi. Both EPR3a and EPR3 function in the intracellular accommodation of microbes. However, contrasting expression patterns and divergent ligand affinities result in distinct functions in AM colonisation and rhizobial infection in Lotus japonicus . The presence of Epr3a and Epr3 genes in both eudicot and monocot plant genomes suggest a conserved function of these receptor kinases in glycan perception.

59 BASIC BIOLOGICAL SCIENCES↗

Design and execution of a verification, validation, and uncertainty quantification plan for a numerical model of left ventricular flow after LVAD implantation

Left ventricular assist devices (LVADs) are implantable pumps that act as a life support therapy for patients with severe heart failure. Despite improving the survival rate, LVAD therapy can carry major complications. Particularly, the flow distortion introduced by the LVAD in the left ventricle (LV) may induce thrombus formation. While previous works have used numerical models to study the impact of multiple variables in the intra-LV stagnation regions, a comprehensive validation analysis has never been executed. The main goal of this work is to present a model of the LV-LVAD system and to design and follow a verification, validation and uncertainty quantification (VVUQ) plan based on the ASME V&V40 and V&V20 standards to ensure credible predictions. The experiment used to validate the simulation is the SDSU cardiac simulator, a bench mock-up of the cardiovascular system that allows mimicking multiple operation conditions for the heart-LVAD system. The numerical model is based on Alya, the BSC’s in-house platform for numerical modelling. Alya solves the Navier-Stokes equation with an Arbitrary Lagrangian-Eulerian (ALE) formulation in a deformable ventricle and includes pressure-driven valves, a 0D Windkessel model for the arterial output and a LVAD boundary condition modeled through a dynamic pressure-flow performance curve. The designed VVUQ plan involves: (a) a risk analysis and the associated credibility goals; (b) a verification stage to ensure correctness in the numerical solution procedure; (c) a sensitivity analysis to quantify the impact of the inputs on the four quantities of interest (QoIs) (average aortic root flow $Q^{avg}_{Ao}$, maximum aortic root flow $Q^{max}_{Ao}$, average LVAD flow $Q^{avg}_{VAD}$, and maximum LVAD flow $Q^{max}_{VAD}$; (d) an uncertainty quantification using six validation experiments that include extreme operating conditions. Numerical code verification tests ensured correctness of the solution procedure and numerical calculation verification showed a grid convergence index (GCI)95% <3.3%. The total Sobol indices obtained during the sensitivity analysis demonstrated that the ejection fraction, the heart rate, and the pump performance curve coefficients are the most impactful inputs for the analysed QoIs. The Minkowski norm is used as validation metric for the uncertainty quantification. It shows that the midpoint cases have more accurate results when compared to the extreme cases. The total computational cost of the simulations was above 100 [core-years] executed in around three weeks time span in Marenostrum IV supercomputer. This work details a novel numerical model for the LV-LVAD system, that is supported by the design and execution of a VVUQ plan created following recognised international standards. We present a methodology demonstrating that stringent VVUQ according to ASME standards is feasible but computationally expensive.

59 BASIC BIOLOGICAL SCIENCES↗