Search NASA⌕ Search

SEARCH · Search NASA

Results for “CLASSES”

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 775 records · Page 43

Fast meta-solvers for 3D complex-shape scatterers using neural operators trained on a non-scattering problem

Three-dimensional target identification using scattering techniques requires high accuracy solutions and very fast computations for real-time predictions in some critical applications. We first train a deep neural operator (DeepONet) to solve wave propagation problems described by the Helmholtz equation in a domain without scatterers but at different wavenumbers and with a complex absorbing boundary condition. We then design two classes of fast meta-solvers by combining DeepONet with either relaxation methods, such as Jacobi and Gauss-Seidel, or with Krylov methods, such as GMRES and BiCGStab, using the trunk basis of DeepONet as a coarse-scale preconditioner. We leverage the spectral bias of neural networks to account for the lower part of the spectrum in the error distribution while the upper part is handled inexpensively using relaxation methods or fine-scale preconditioners. The meta-solvers are then applied to solve scattering problems with different shape of scatterers, at no extra training cost. We first demonstrate that the resulting meta-solvers are shape-agnostic, fast, and robust, whereas the standard standalone solvers may even fail to converge without the DeepONet. We then apply both classes of meta-solvers to scattering from a submarine, a complex three-dimensional problem. We achieve very fast solutions, especially with the DeepONet-Krylov methods, which require orders of magnitude fewer iterations than any of the standalone solvers.

97 MATHEMATICS AND COMPUTING↗

Finite spatial-grid effects in energy-conserving particle-in-cell algorithms

Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These instabilities originate in the aliasing of interpolation errors between mesh quantities and particles (which live in the space–time continuum). Recently, strictly energy-conserving PIC (EC-PIC) algorithms have been developed that promise enhanced robustness against aliasing instabilities. In this study, we confirm by analysis that EC-PIC is stable against aliasing instabilities for stationary plasmas. For drifting plasmas, we demonstrate by analysis and numerical experiments that, while EC-PIC algorithms are not free from these instabilities in principle, they feature a benign stability threshold for finite-temperature plasmas that make them usable in practice for a large class of problems (featuring ambipolarity and realistic ion-electron mass ratios) without the need to consider the size of the Debye length. We also demonstrate that this threshold is absent for the popular momentum-conserving PIC algorithms, which are therefore unstable for both drifting and stationary plasmas beyond a threshold in cell size compared to Debye length. Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These instabilities originate in the aliasing of interpolation errors between mesh quantities and particles (which live in the space–time continuum). Recently, strictly energy-conserving PIC (EC-PIC) algorithms have been developed that promise enhanced robustness against aliasing instabilities. In this study, we confirm by analysis that EC-PIC is stable against aliasing instabilities for stationary plasmas. For drifting plasmas, we demonstrate by analysis and numerical experiments that, while EC-PIC algorithms are not free from these instabilities in principle, they feature a benign stability threshold for finite-temperature plasmas that make them usable in practice for a large class of problems (featuring ambipolarity and realistic ion-electron mass ratios) without the need to consider the size of the Debye length. Finally, we also demonstrate that this threshold is absent for the popular momentum-conserving PIC algorithms, which are therefore unstable for both drifting and stationary plasmas beyond a threshold in cell size compared to Debye length.

97 MATHEMATICS AND COMPUTING↗

GenASiS Mathematics: Object-oriented manifolds, operations, and solvers for large-scale physics simulations (version 2)

We report GenASiS Mathematics provides modern Fortran classes furnishing extensible object-oriented functionality for the solution of fields governed by selected partial differential equations. The initial release included extensible object-oriented implementations of simple meshes and the evolution of generic conserved currents thereon. This revision - Version 2 of Mathematics - includes significant reorganization and streamlining of these classes, higher-order reconstruction by a different method, a Poisson solver, coarsening to avoid Courant time step limitations near coordinate singularities, and the offloading of computational kernels to GPUs.

97 MATHEMATICS AND COMPUTING↗

Maize specialized metabolome networks reveal organ-preferential mixed glycosides

Despite the scientific and economic importance of maize, little is known about its specialized metabolism. Here, five maize organs were profiled using different reversed-phase liquid chromatography-mass spectrometry methods. The resulting spectral metadata, combined with candidate substrate-product pair (CSPP) networks, allowed the structural characterization of 427 of the 5,420 profiled compounds, including phenylpropanoids, flavonoids, benzoxazinoids, and auxin-related compounds, among others. Only 75 of the 427 compounds were already described in maize. Analysis of the CSPP networks showed that phenylpropanoids are present in all organs, whereas other metabolic classes are rather organ-enriched. Frequently occurring CSPP mass differences often corresponded with glycosyl- and acyltransferase reactions. The interplay of glycosylations and acylations yields a wide variety of mixed glycosides, bearing substructures corresponding to the different biochemical classes. For example, in the tassel, many phenylpropanoid and flavonoid-bearing glycosides also contain auxin-derived moieties. The characterized compounds and mass differences are an important step forward in metabolic pathway discovery and systems biology research.

59 BASIC BIOLOGICAL SCIENCES↗

Evaluation of riparian enhancement actions in the Columbia River Basin

Riparian enhancement is a common restoration technique in the Columbia River Basin (CRB) and the Pacific Northwest. However, relatively few studies have evaluated its success and even fewer studies include long-term monitoring. Forty-one riparian planting projects located in the CRB, each with a paired treatment and control reach, were evaluated in the summer of 2018 and 2019 using an extensive post-treatment (EPT) design. At each reach, we quantified woody plant abundance, richness, diversity, and vegetation and canopy cover as riparian response variables, and measured terrace height and documented individuals with bud browse, deceased, and with predator protection as potential explanatory variables. Species richness, woody plant abundance (all height classes combined and within shrub height class), and the proportion of woody plants with bud browse or deceased were all higher in treatment than control reaches. However, no significant improvements were observed for any other riparian response or explanatory variables. Implementation techniques, including invasive species removal efforts, follow-up planting, watering, and the type of restoration project (floodplain restoration and planting, in-stream restoration and planting, or solely planting) were also considered as potential explanatory variables and were found to have significant impacts on restoration response for abundance and cover variables. Further, our results suggest that planting implementation methods, site level physical factors, and time since planting all influence the success of riparian planting projects and if not addressed in restoration design and implementation can contribute to a lack of detectable response.

54 ENVIRONMENTAL SCIENCES↗

Agent-based modeling to evaluate the effects of harvesting biomass and hunting on ring-necked pheasant ( Phasianus colchicus ) populations

As the US strives to divest from fossil fuels, dedicated energy crops have increased in popularity because of their potential to serve as a renewable source of transportation fuel and biomass-generated electricity. Whereas the impacts of crop agriculture on biodiversity have been well-documented, less is known about the potential impacts that growing bioenergy crops could have on wildlife. Here we look specifically at the ring-necked pheasant (Phasianus colchicus) and how pheasant populations may be influenced by growing and harvesting bioenergy crops. We explored effects of temporal harvest strategies at the field scale on biomass yield and pheasant population size. To investigate, we developed an agent-based model (ABM) that simulates ring-necked pheasants, tractors, hunters, and vegetation classes. One such vegetation class is land that is enrolled in a Conservation Reserve Program (CRP) plan specifically targeting pheasant conservation. Using this ABM, we assessed four different landscapes—corn-dominated, CRP-dominated, grassland-dominated, and mixed landscape—under strategies that varied the time of harvest. We also used ecological valuation to compare scenario outcomes from an economic perspective. We determined that biomass yield and pheasant population size were sensitive to harvesting times. Our scenarios totaled between ~$931,000 (minimum) and $3.8 million (maximum) over all landscapes, though on average harvesting in late spring generated the most value ($2.45 million) and harvesting in fall was a close second in terms of profitability ($2.39 million).

59 BASIC BIOLOGICAL SCIENCES↗

Coulombic friction in metamaterials to dissipate mechanical energy

Product designs from a wide range of industries such as aerospace, automotive, biomedical, and others can benefit from new metamaterials for mechanical energy dissipation. In this study, we explore a novel new class of metamaterials with unit cells that absorb energy via sliding Coulombic friction. Remarkably, even materials such as metals and ceramics, which typically have no intrinsic reversible energy dissipation, can be architected to provide dissipation akin to elastomers. The concept is demonstrated at different scales (centimeter to micrometer), with different materials (metal and polymer), and in different operating environments (high and low temperatures), all showing substantial dissipative improvements over conventional non-contacting lattice unit cells. Further, as with other ‘programmable’ metamaterials, the degree of Coulombic absorption can be tailored for a given application. Additionally, an analytic expression is derived to allow rapid first-order optimization. This new class of Coulombic friction energy absorbers can apply broadly to many industrial sectors such as transportation (e.g. monolithic shock absorbers), biomedical (e.g. prosthetics), athletic equipment (e.g. skis, bicycles, etc.), defense (e.g. vibration tolerant structures), and energy (e.g. survivable electrical grid components).

36 MATERIALS SCIENCE↗

Profiling overstory survival trends following varying thinning and burning disturbance regimes in a mixed pine-hardwood forest in the US South

Prescribed burning as a silvicultural tool may be effective for achieving various management objectives, yet some practitioners have concerns about inadvertently increasing overstory tree mortality. Using extensive data from a long-term, ongoing study in northcentral Alabama, USA, that systematically applies a variety of thinning and prescribed burning disturbances on mixed pine-hardwood stands, we assess survival trends for different groups of trees. The primary research question is whether more frequent prescribed burns adversely affects overstory survival versus infrequent or no burns; secondary questions explore whether overstory survival trends differ based on thinning level, species group, or size class. Interest is in broad groups of trees, not individual tree survival or mortality; consequently, survival analysis methods are used, including nonparametric Kaplan-Meier (KM) techniques for examining single grouping factors, and parametric accelerated failure time models for analyzing simultaneous effects of multiple covariates. Leveraging relatively recent methodological advances, the statistical techniques are adjusted for interval censored data. KM survival curves showed that more frequent prescribed fires did not result in differing survival trends. Also, trees in unthinned stands had the lowest 14-year survival compared to thinned stands, oaks had the highest survival compared to pines and the others species group, and the smallest size class of overstory trees had the lowest survival. In conclusion, these results support managing transitional mixed pine-hardwoods using thinning and multiple prescribed fires to restore specific species composition and structure.

54 ENVIRONMENTAL SCIENCES↗

Signatures of rare earth element distributions in fly ash derived from the combustion of Central Appalachian, Illinois, and Powder River basin coals

We report the distribution of Rare Earth elements (REE) in coal-derived fly ashes can have distinctive patterns when fly ashes are produced from different coals within or between basins, such as the Pennsylvanian Class F fly ashes from the Illinois and Central Appalachian basins. Both the Fire Clay coal and a blend of a number of eastern Kentucky coals show strong Gd peaks and an H-type distribution in the Upper Continental Crust-corrected plots. The Fire Clay coal-derived ash has a higher heavy REE concentration than the blended coal-derived ash. The Illinois Basin-derived fly as has an overall lower REE concentration than the latter ashes. Class C fly ash derived from Powder River Basin coals has, with the exception of an Eu peak, a flatter distribution of REE and an overall L-type or indistinct H- versus L-type distribution. The signatures of the REE in fly ashes may be useful in predicting their behavior in the extraction of the REE; simple extrapolations from the basic concentrations and the predicted extraction percentages for ashes from different basins are not necessarily indicative of the actual distribution of the extracted REE.

01 COAL, LIGNITE, AND PEAT↗

Parallel sorting algorithm classification: is manual instrumentation necessary?

Understanding parallel algorithms is crucial for accelerating scientific simulations on complex, distributed memory, high-performance computers. Modern algorithm classification approaches learn semantics directly from source code to differentiate between algorithms, however, accessing source code is not always possible. We can learn about parallel algorithms from observing their performance, as programs running the same algorithms and using the same hardware should exhibit similar performance characteristics. We present an approach to learn algorithm classes from parallel performance data directly in order to classify algorithms without access to the source code. We extend previous work to enable classifying parallel sorting algorithms using automatic instrumentation instead of requiring manual region annotations in the source code. In this work, we design and demonstrate a study for classification of parallel sorting algorithms using parallel performance data collected from automatic instrumentation, and evaluate the performance of our new methodology on classification. We leverage Caliper to collect the performance data, Thicket for our exploratory data analysis (EDA), and PyTorch and Scikit-learn to evaluate the effectiveness of random forests, support vector machines (SVMs), decision trees, neural networks, and logistic regressions on parallel performance data. Additionally, we study noise in parallel performance data, whether the removal of noise and pre-processing of the data is necessary to accurately classify parallel sorting algorithms, and determine the effectiveness of features created from performance data. In conclusion, we demonstrate classification accuracy for these five different models of up to 97.7% across four different parallel algorithm classes.

Algorithm Classification↗

On the convexity of phase-field fracture formulations: Analytical study and comparison of various degradation functions

Efficient and accurate fracture modeling is of great importance in applications where catastrophic outcomes under extreme scenarios are possible. The phase-field (PF) approach to fracture received significant attention over the past decade, due to its capability to capture complicated fracture patterns (e.g., crack merging and branching). Specifically, crack initiation and propagation are modeled via minimization of the total energy functional, which is regularized with the aid of a phase field. Despite the promising results and modeling capabilities of the PF method in many applications, the solution of fracture problems remains computationally challenging mainly due to the non-convexity of the total energy functional with respect to the combined unknown (phase field and displacement) fields. Understanding the effects of their coupling on convexity is crucial in order to address frequently encountered hurdles in fracture modeling (e.g., inefficient solvers and non-physical crack nucleation). In this paper, we develop convexity criteria for a wide class of PF fracture formulations. For this class of formulations, the second variation of the total energy functional is expressed in terms of Hessian matrices (evaluated at individual material points). Depending on the choice of geometric crack functions and degradation functions, we classify the formulations into three categories and analytically study each one separately. To study the sign of the second variation, we derive inequalities which are satisfied at material points when the Hessian matrix is locally positive semi-definite. These inequalities provide objective criteria for comparing degradation functions. Finally, the applicability of the proposed convexity criteria is demonstrated in the context of a one-dimensional problem, solved using a conventional monolithic solver.

97 MATHEMATICS AND COMPUTING↗

Neutralizing monoclonal antibodies elicited by mosaic RBD nanoparticles bind conserved sarbecovirus epitopes

Increased immune evasion by SARS-CoV-2 variants of concern highlights the need for new therapeutic neutralizing antibodies. Immunization with nanoparticles co-displaying spike receptor-binding domains (RBDs) from eight sarbecoviruses (mosaic-8 RBD-nanoparticles) efficiently elicits cross-reactive polyclonal antibodies against conserved sarbecovirus RBD epitopes. Here, we identified monoclonal antibodies (mAbs) capable of cross-reactive binding and neutralization of animal sarbecoviruses and SARS-CoV-2 variants by screening single mouse B cells secreting IgGs that bind two or more sarbecovirus RBDs. Single-particle cryo-EM structures of antibody-spike complexes, including a Fab-Omicron complex, mapped neutralizing mAbs to conserved class 1/4 RBD epitopes. Structural analyses revealed neutralization mechanisms, potentials for intra-spike trimer cross-linking by IgGs, and induced changes in trimer upon Fab binding. In addition, we identified a mAb-resembling Bebtelovimab, an EUA-approved human class 3 anti-RBD mAb. These results support using mosaic RBD-nanoparticle vaccination to generate and identify therapeutic pan-sarbecovirus and pan-variant mAbs.

59 BASIC BIOLOGICAL SCIENCES↗

Stabilization of glucose-6-phosphate dehydrogenase oligomers enhances catalytic activity and stability of clinical variants

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is a genetic trait that can cause hemolytic anemia. To date, over 150 nonsynonymous mutations have been identified in G6PD, with pathogenic mutations clustering near the dimer and/or tetramer interface and the allosteric NADP + -binding site. Recently, our lab identified a small molecule that activates G6PD variants by stabilizing the allosteric NADP + and dimer complex, suggesting therapeutics that target these regions may improve structural defects. Here, we elucidated the connection between allosteric NADP + binding, oligomerization, and pathogenicity to determine whether oligomer stabilization can be used as a therapeutic strategy for G6PD deficiency (G6PD def ). We first solved the crystal structure for G6PD K403Q , a mutant that mimics the physiological acetylation of wild-type G6PD in erythrocytes and demonstrated that loss of allosteric NADP + binding induces conformational changes in the dimer. These structural changes prevent tetramerization, are unique to Class I variants (the most severe form of G6PD def ), and cause the deactivation and destabilization of G6PD. We also introduced nonnative cysteines at the oligomer interfaces and found that the tetramer complex is more catalytically active and stable than the dimer. Furthermore, stabilizing the dimer and tetramer improved protein stability in clinical variants, regardless of clinical classification, with tetramerization also improving the activity of G6PD K403Q and Class I variants. These findings were validated using enzyme activity and thermostability assays, analytical size-exclusion chromatography (SEC), and SEC coupled with small-angle X-ray scattering (SEC-SAXS). Taken together, our findings suggest a potential therapeutic strategy for G6PD def and provide a foundation for future drug discovery efforts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The crystal structure of Grindelia robusta 7,13-copalyl diphosphate synthase reveals active site features controlling catalytic specificity

Diterpenoid natural products serve critical functions in plant development and ecological adaptation and many diterpenoids have economic value as bioproducts. The family of class II diterpene synthases catalyzes the committed reactions in diterpenoid biosynthesis, converting a common geranylgeranyl diphosphate precursor into different bicyclic prenyl diphosphate scaffolds. Enzymatic rearrangement and modification of these precursors generate the diversity of bioactive diterpenoids. We report the crystal structure of Grindelia robusta 7,13-copalyl diphosphate synthase, GrTPS2, at 2.1 Å of resolution. GrTPS2 catalyzes the committed reaction in the biosynthesis of grindelic acid, which represents the signature metabolite in species of gumweed (Grindelia spp., Asteraceae). Grindelic acid has been explored as a potential source for drug leads and biofuel production. The GrTPS2 crystal structure adopts the conserved three-domain fold of class II diterpene synthases featuring a functional active site in the γβ-domain and a vestigial α-domain. Substrate docking into the active site of the GrTPS2 apo protein structure predicted catalytic amino acids. Biochemical characterization of protein variants identified residues with impact on enzyme activity and catalytic specificity. Specifically, mutagenesis of Y457 provided mechanistic insight into the position-specific deprotonation of the intermediary carbocation to form the characteristic 7,13 double bond of 7,13-copalyl diphosphate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Improved lipid analysis using a 2D-LC-MS system with a novel injection procedure

The aim of this study was to improve analysis of nonpolar lipidomics sample extracts using reversed phase (RP) chromatography. A 4/3/3 (v/v/v) mixture of methanol/methyl tert-butyl ether/chloroform (MeOH/MTBE/CHCl 3 , MMC) was chosen for sample extraction solvent based on its proven extraction capability for several lipid classes. To avoid carry over, loss of analytes and peak distortion the loops and all capillaries of the presented LC system were flushed and filled up with methanol until the analytical column. Here, the choice of methanol was due to its weak elution strength and being infinitely miscible with MMC and several other nonpolar solvents. This allowed injection of a 100 μl sample that was 20 μl nonpolar extraction solvent diluted fivefold with methanol. All lipids of 25 lipid classes were transferred quantitatively to the column head where the online dilution of methanol was carried out with aqueous eluent for focusing the lipid analytes. The weak elution strength of methanol prevented peak distortions. The consecutive reversed phase elution resulted in remarkably narrow peaks (full width at half maximum was 0.07–0.08 min typically) and enhanced sensitivity (limit of detection usually in sub nM region) because of increased sample injection volume and narrow peaks. Calibration and quality control samples made by diluting commercial lipid standards 200–50000 times confirmed the applicability of this approach both for targeted lipid quantification and for untargeted quantitative comparison of lipids from different sources.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Goal-oriented a-posteriori estimation of model error as an aid to parameter estimation

In this work, a Bayesian model calibration framework is presented that utilizes goal-oriented a-posterior error estimates in quantities of interest (QoIs) for classes of high-fidelity models characterized by PDEs. It is shown that for a large class of computational models, it is possible to develop a computationally inexpensive procedure for calibrating parameters of high-fidelity models of physical events when the parameters of low-fidelity (surrogate) models are known with acceptable accuracy. The main ingredients in the proposed model calibration scheme are goal-oriented a-posteriori estimates of error in QoIs computed using a so-called lower fidelity model compared to those of an uncalibrated higher fidelity model. The estimates of error in QoIs are used to define likelihood functions in Bayesian inversion analysis. A standard Bayesian approach is employed to compute the posterior distribution of model parameters of high-fidelity models. As applications, parameters in a quasi-linear second-order elliptic boundary-value problem (BVP) are calibrated using a second-order linear elliptic BVP. In a second application, parameters of a tumor growth model involving nonlinear time-dependent PDEs are calibrated using a lower fidelity linear tumor growth model with known parameter values.

A-posterior estimates↗

Investigating uncertainties in human adaptation and their impacts on water scarcity in the Colorado river Basin, United States

The Colorado River Basin (CRB) supports the water supply for seven states and forty million people in the Western United States (US) and has been suffering an extensive drought for more than two decades. As climate change continues to reshape water resources distribution in the CRB, its impact can differ in intensity and location, resulting in variations in human adaptation behaviors. The feedback from human systems in response to the environmental changes and the associated uncertainty is critical to water resources management, especially for water-stressed basins. This paper investigates how human adaptation affects water scarcity uncertainty in the CRB and highlights the uncertainties in human behavior modeling. Our focus is on agricultural water consumption, as approximately 80% of the water consumption in the CRB is used in agriculture. We adopted a coupled agent-based and water resources modeling approach for exploring human-water system dynamics, in which an agent is a human behavior model that simulates a farmer’s water consumption decisions. We examined uncertainties at the system, agent, and parameter levels through uncertainty, clustering, and sensitivity analyses. The uncertainty analysis results suggest that the CRB water system may experience 13 to 30 years of water shortage during the 2019–2060 simulation period, depending on the paths of farmers’ adaptation. The clustering analysis identified three decision-making classes: bold, prudent, and forward-looking, and quantified the probabilities of an agent belonging to each class. The sensitivity analysis results indicated agents whose decision-making models require further investigation and the parameters with the higher uncertainty reduction potentials. Here, by conducting numerical experiments with the coupled model, this paper presents quantitative and qualitative information about farmers’ adaptation, water scarcity uncertainties, and future research directions for improving human behavior modeling.

Agent-based modeling↗

Structural Insights into the Mechanisms and Pharmacology of K 2P Potassium Channels

Leak currents, defined as voltage and time independent flows of ions across cell membranes, are central to cellular electrical excitability control. The K 2P (KCNK) potassium channel class comprises an ion channel family that produces potassium leak currents that oppose excitation and stabilize the resting membrane potential in cells in the brain, cardiovascular system, immune system, and sensory organs. Due to their widespread tissue distribution, K 2P s contribute to many physiological and pathophysiological processes including anesthesia, pain, arrythmias, ischemia, hypertension, migraine, intraocular pressure regulation, and lung injury responses. Structural studies of six homomeric K 2P s have established the basic architecture of this channel family, revealed key moving parts involved in K 2P function, uncovered the importance of asymmetric pinching and dilation motions in the K 2P selectivity filter (SF) C-type gate, and defined two K 2P structural classes based on the absence or presence of an intracellular gate. Further, a series of structures characterizing K 2P :modulator interactions have revealed a striking polysite pharmacology housed within a relatively modestly sized (~70 kDa) channel. Binding sites for small molecules or lipids that control channel function are found at every layer of the channel structure, starting from its extracellular side through the portion that interacts with the membrane bilayer inner leaflet. This framework provides the basis for understanding how gating cues sensed by different channel parts control function and how small molecules and lipids modulate K 2P activity. Such knowledge should catalyze development of new K 2P modulators to probe function and treat a wide range of disorders.

59 BASIC BIOLOGICAL SCIENCES↗