Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “projection method”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Identification of new marker genes from plant single‐cell RNA‐seq data using interpretable machine learning methods

Summary An essential step in the analysis of single‐cell RNA sequencing data is to classify cells into specific cell types using marker genes. In this study, we have developed a machine learning pipeline called single‐cell predictive marker (SPmarker) to identify novel cell‐type marker genes in the Arabidopsis root. Unlike traditional approaches, our method uses interpretable machine learning models to select marker genes. We have demonstrated that our method can: assign cell types based on cells that were labelled using published methods; project cell types identified by trajectory analysis from one data set to other data sets; and assign cell types based on internal GFP markers. Using SPmarker, we have identified hundreds of new marker genes that were not identified before. As compared to known marker genes, the new marker genes have more orthologous genes identifiable in the corresponding rice single‐cell clusters. The new root hair marker genes also include 172 genes with orthologs expressed in root hair cells in five non‐ Arabidopsis species, which expands the number of marker genes for this cell type by 35–154%. Our results represent a new approach to identifying cell‐type marker genes from scRNA‐seq data and pave the way for cross‐species mapping of scRNA‐seq data in plants.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Topsoil bulk geochemical compositions - An updated harmonized global dataset

Mineral weathering is a key biogeochemical process because of the capacity of minerals to stabilize organic matter. However, predicting soil weathering status across large spatial areas still isn’t possible due to a lack of global data and theoretical frameworks. To address this knowledge gap, multiple global datasets of bulk topsoil geochemical compositions have been harmonized using R. These datasets document topsoil bulk geochemical compositions across five continents (n = ~16,000 observations). Source data for these observations include the EuroGEOSurveys Geochemical Baseline Database (FOREGS), the US Geological Survey National Geochemical Database (NASGLP), the Geochemical Atlas of Australia (GAA), the US Geological Survey Alaska Geochemical Database (AGD84), the National Cooperative Soil Survey (NCSS), the European Geochemical Mapping of Agricultural Soil (GEMAS), Ecorespira-Amazon (ERA), the New Zealand Geochemical Baseline Survey (NZ_GBS), and the African Soil Information Service (AFSIS). Major elements observed include Aluminum (Al), Calcium (Ca), Iron (Fe), Potassium (K), Magnesium (Mg), Sodium (Na), Titanium (Ti), Manganese (Mn), Phosphorus (P), Carbon (C), and Sulfur (S). This data package includes the harmonized dataset itself, and the R scripts necessary to harmonize these datasets, in addition to metadata that describes all columns, files, and databases used in this project. Methods & Sampling Step 1 – Databases of geochemical data identified This study aimed to leverage existing measurements of topsoil geochemical data. Databases were first identified and deemed appropriate for inclusion if they were measuring soils and performed these measurements on the <2mm soil fraction. Databases such as NCSS and AGD84 needed more post processing to include in the database and this was done using the NCSS_datamerge_031626 R file and Alaska_USGSmerge_031626 R file, respectively. Step 2 – Database harmonization Once appropriate databases were identified, they were harmonized for ease of analysis using the R script Database_Harmonization_031826. This included removing columns from original datasets that would not be used in analysis (removed columns are noted in the code). Then, data cleaning procedures specific to each dataset were undertaken. This includes standardizing columns to include units and adding metadata columns regarding procedures for analyzing specific elements. Functions for standardizing measurements and units are outline in R files: calculate element_mg_kg_031626, calculate_oxide_wt_perc_031626, change_oxide_caps_031626, and conv_2_numeric_031626. This also included adding a unique identifier for each sample to identify it with its respective database (see CD_ID in data dictionary). Geographic information: Data reflect a compilation of datasets collected globally. Geographic areas covered by each of the datasets include: - EuroGEOSurveys Geochemical Baseline Database (FOREGS) - European continent - North American Soil Geochemical Landscapes (NASGLP) - continental United States and limited parts of Canada (see database key for more details) - National Geochemical Survey of Australia (GAA) - Australia - Alaska geochemical database (AGDB4) - Alaska - National Cooperative Soil Survey (NCSS) - Global measurements, but concentrated in the continental United States - Geochemical data for arable land and land under permanent grass cover in continental Europe (GEMAS) - continental Europe - Ecorespira-Amazon (ERA) - Geochemical data from the Amazon basin - Geochemical baseline data for New Zealand (NZGBS) - New Zealand - Geochemical data collected across continental Africa (AfSIS) - Measurements across Africa

EARTH SCIENCE > LAND SURFACE > SOILS↗
DOE OSTI2022

Development of the uncertainty quantification toolkit's python interface and surrogate construction tutorial

The uncertainty quantification toolkit (UQTk) is a collection of c++ libraries that assess the confidence of numerical models. Surrogate approximations, often polynomial chaos expansions (PCEs), lessen the computational cost of these assessments. I developed a Python interface in UQTk for regression and Bayesian compressive sensing to add to the existing Galerkin projection method. These methods receive an object containing the polynomial basis information and NumPy arrays of sample points, call c++ methods, and return the PCE coefficients in a NumPy array. To demonstrate these methods, I wrote a tutorial in which I use them to construct surrogates for Genz functions and calculate the resulting error.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Electron Spin Polarization in Large Electric Fields

This project investigated the effects of large electric fields on electron spin polarization in semiconductors, which is of interest for spin-based electronics (“spintronics”) and quantum information processing. It has been theoretically proposed that large electric fields could potentially preserve or amplify electron spin polarization, which would make it easier to generate, transport, and detect spin polarization in spin-based devices; however, previous measurements had not measured the response of electron spin polarization at large electric fields. During this project, methods to perform measurements of electron spin polarization at large electric fields were developed and performed using ultrafast time- and spatially-resolved magneto-optical spectroscopy, high-voltage pulses, and lock-in detection. The device geometry was engineered to reduce the threshold voltage required to reach the Gunn threshold, and the effect of optical illumination on the sample conductivity and Gunn oscillations was investigated. Electron spin polarization measurements were performed at electric fields up to and beyond the Gunn threshold for the device. The measurements showed that the electron spin polarization magnitude decreased with increasing electric field in the devices measured, and a spin amplification effect was not observed. However, these effects could potentially exist in a different material or in a material with a different doping concentration or in a different device geometry.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Rigid-Mode Limit of the Yokoya Matrix Formalism and the Burov-Lebedev Dispersion Equation

Transverse single-bunch instabilities of space-charge-dominated coasting beams with round and flat transverse geometries are studied using a unified dispersion-relation framework. The analysis combines the Burov-Lebedev formalism, which captures space-charge tune spread, Landau damping, and instability threshold behavior, with Yokoya’s projection method for representing coherent transverse mode structure and its dependence on beam aspect ratio. In the rigid-beam limit, the formulation reduces to a scalar dispersion relation of Burov-Lebedev paper. For non-rigid transverse oscillations, truncation of Yokoya’s Hermite-based expansion yields a finite-dimensional matrix eigenvalue problem in which space-charge and coupling impedance effects enter through Burov-Lebedev–type denominators. This approach provides a consistent basis for comparing rigid and non-rigid instability behavior in round and flat beams and for assessing the role of beam ellipticity in modifying coherent mode structure and stability thresholds.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Implementation of Manifold-Based Combustion Models in a Highly Scalable Low Mach Number Reacting Flow Solver: Preprint

Manifold-based representations of the thermochemistry are often employed in conjunction with large eddy eimulation (LES) to lower the cost of combustion simulations. This work describes steps taken to implement this modeling approach in PeleLM, a scalable and performance-portable low Mach number flow solver. Most significantly, this includes adapting the projection method used by PeleLM to satisfy the mass conservation constraint for use with manifold-based models. The implementation is designed to be general across manifold-based models, including both those that employ traditional tabulation and those that employ neural networks. An initial demonstration for simple test cases is presented and will be used for performance assessment.

high-performance computing↗
DOE OSTI2025

Symmetry-projected spin-AGP methods applied to spin systems

Symmetry-projected wave function methods capture static correlation by breaking and restoring the symmetries of a system. In this article, we present the symmetry-projected spin antisymmetrized geminal power (spin-AGP) state projected onto space group symmetry as well as complex conjugation, spin-flip, and time-reversal symmetries. The method is benchmarked on the 1D XXZ model and the 2D J 1 − J 2 model with square and triangular lattices. Our results indicate that symmetry projection methods provide a powerful tool for frustrated spin systems.

Antisymmetrized geminal power↗
DOE OSTI2024

An Evaluation of Dynamical Downscaling Methods Used to Project Regional Climate Change

In the past decade, dynamical downscaling using “pseudo-global-warming” (PGW) techniques has been applied frequently to project regional climate change. Such techniques generate signals by adding mean global climate model (GCM)-simulated climate change signals in temperature, moisture, and circulation to lateral and surface boundary conditions derived from reanalysis. An alternative to PGW is to downscale GCM data directly. This technique should be advantageous, especially for simulation of extremes, since it incorporates the GCM's full spectrum of changing synoptic-scale dynamics in the regional solution. Here, we test this assumption, by comparing simulations in Europe and Western North America. We find that for warming and changes in temperature extremes, PGW often produces similar results to direct downscaling in both regions. For mean and extreme precipitation changes, PGW generally also performs surprisingly well in many cases. Moisture budget analysis in the Western North America domain reveals why. Large fractions of the downscaled hydroclimate changes arise from mean changes in large-scale thermodynamics and circulation, that is, increases in temperature, moisture, and winds, included in PGW by design. The one component PGW may have difficulty with is the contribution from changes in synoptic-scale variability. When this component is large, PGW performance could be degraded. Global analysis of GCM data shows there are regions where it is large or dominant. Hence, our results provide a road map to identify, through GCM analyses, the circumstances when PGW would not be expected to accurately regionalize GCM climate signals.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Adaptive Mesh Refinement for Parallel in Time Methods

The project applied the multigrid-reduction-in-time (MGRIT) algorithm to an existing sub-cycled adaptive mesh refinement (AMR) code to investigate the performance of flows dominated by inertial physics. Previous work demonstrated good performance from MGRIT+AMR applied to flows dominated by diffusive physics. Consistent with previous experience, inertial physics negatively affected convergence rates and performance. Efforts to circumvent this issue by appealing to the physics of turbulence were investigated. It has been demonstrated that scales can be effectively transferred between multigrid levels for a turbulent flow resulting in a) partial convergence observed and b) nearly identical results to sequential time-stepping. Performance improvements have not yet been demonstrated - attempts at coarsening the grid on coarser MG levels compromises the solution quality and leads to divergence. This report summarize the accomplishments for the time-frame from 10/5/2020 to 12/31/2020 with an informal no-cost-extension to 05/20/2021

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Tracking Photochemical and Photophysical Processes for Solar Energy Conversion Via Multidimensional Visible and Vibrational Spectroscopic Methods

This project was motivated by an overarching goal to elucidate the mechanism of energy and electron transfer that governs the efficient charge separation in photosystem I (PSI) complexes. PSI is a natural light harvesting complex that drives oxygenic photosynthesis in plants, algae, and cyanobacteria. It uses ~300 tightly packed chlorophylls (Chls) to absorb photons, transfer the excitation energy to the reaction center (RC), and generate a charge separated state with near unity quantum efficiency (QE). A better understanding of the mechanism of energy transfer and charge separation in PSI is required for understanding the high QE of natural light harvesting complexes, and it could lead to the further development of artificial photosynthetic systems for solar energy conversion. We applied two-dimensional electronic spectroscopy (2DES) to different cyanobacterial photosystem I complexes to map energy transfer pathways to gain insight into the efficient light harvesting of PSI. We interpreted the 2DES spectra through a global analysis procedure, finding two ultrafast equilibration processes involving red-shifted Chl states. The analysis procedure is general and can be readily applied to other natural and artificial light harvesting systems. In addition to investigating PSI complexes, we also performed measurements on model systems to establish general data analysis procedures for interpreting 2D spectra.

14 SOLAR ENERGY↗
DOE OSTI2022

Accelerated Discovery of Solar Thermochemical Hydrogen Production Materials via High-Throughput Computational and Experimental Methods

In this project, combinatorial synthesis and testing methods were combined with high-throughput materials theory calculations to greatly accelerate the discovery of thermodynamically suitable candidates for green hydrogen production via a two-stage solar thermochemical water splitting (STCH) process. Over the course of the project, more than 8000 quinary and higher oxide compositions were computationally screened for STCH viability, and detailed stability calculations were performed for more than 30 of the most promising identified compositional archetypes. As a result, three new STCH capable compositional families were discovered and experimentally verified. The first, Ce x Sr 2-x MnO 4 (CSM), represents the first known Ruddlesden-Popper compound to show STCH activity, and thus demonstrates that perovskite-related structures may hold promise for this application. The second family, Sr 1-x Ce x MnO 3 (SCM), is the simple perovskite sister-analog to CSM. Sr 0.7 Ce 0.3 MnO 3 (SCM30), a member of this compositional family, was found to produce the highest hydrogen yields of any compound tested in this project, exceeding the end of project milestone target of > 150 μmol H 2 /gram oxide at a reduction temperature of 1350 °C, although only at steam-to-hydrogen ratios greater than 1000:1. Finally, we proved that a third novel Sr-and Mn-containing family, Sr 1-x Ca x Ti 1-y Mn y O 3 (SCTM), which was identified by Materials Project tools, also splits water. The behavior of the SCTM system was found to be similar to the previously discovered Sr 1-x La x Al 1-y Mn y O 3 (SLMA) family, albeit with lower H 2 yields. Across the three thrusts of the project (computational, combinatorial, and bulk testing), five journal articles were published. As part of Program End Analysis and Data Dissemination, relevant data used for the publications was uploaded to the HydroGEN Data Hub for public access, and in certain cases, results were added to public materials databases.

08 HYDROGEN↗
DOE OSTI2022

Training and projecting: A reduced basis method emulator for many-body physics

Here, we present the reduced basis method as a tool for developing emulators for equations with tun able parameters within the context of the nuclear many-body problem. The method uses a basis expansion informed by a set of solutions for a few values of the model parameters and then projects the equations over a well-chosen low-dimensional subspace. We connect some of the results in the eigenvector continuation literature to the formalism of reduced basis methods and show how these methods can be applied to a broad set of problems. As we illustrate, the possible success of the formalism on such problems can be diagnosed beforehand by a principal component analysis. We apply the reduced basis method to the one-dimensional Gross-Pitaevskii equation with a harmonic trap ping potential and to nuclear density functional theory for 48 Ca, achieving speed-ups of more than x150 and x250, respectively, when compared to traditional solvers. The outstanding performance of the approach, together with its straightforward implementation, show promise for its application to the emulation of computationally demanding calculations, including uncertainty quantification.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Using projection operators with maximum overlap methods to simplify challenging self-consistent field optimization

Maximum overlap methods are effective tools for optimizing challenging ground- and excited-state wave functions using self-consistent field models such as Hartree-Fock and Kohn-Sham density functional theory. Nevertheless, such models have shown significant sensitivity to the user-defined initial guess of the target wave function. Here in this work, a projection operator framework is defined and used to provide a metric for non-aufbau orbital selection in maximum-overlap-methods. The resulting algorithms, termed the Projection-based Maximum Overlap Method (PMOM) and Projection-based Initial Maximum Overlap Method (PIMOM), are shown to perform exceptionally well when using simple user-defined target solutions based on occupied/virtual molecular orbital permutations. This work also presents a new metric that provides a simple and conceptually convenient measure of agreement between the desired target and the current or final SCF results during a calculation employing a maximum-overlap method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Atmospheric River Detection Under Changing Seasonality and Mean-State Climate: ARTMIP Tier 2 Paleoclimate Experiments

Atmospheric rivers (ARs) are filamentary structures within the atmosphere that account for a substantial portion of poleward moisture transport and play an important role in Earth's hydroclimate. However, there is no one quantitative definition for what constitutes an atmospheric river, leading to uncertainty in quantifying how these systems respond to global change. This study seeks to better understand how different AR detection tools (ARDTs) respond to changes in climate states utilizing single-forcing climate model experiments under the aegis of the Atmospheric River Tracking Method Intercomparison Project (ARTMIP). We compare a simulation with an early Holocene orbital configuration and another with CO2 levels of the Last Glacial Maximum to a preindustrial control simulation to test how the ARDTs respond to changes in seasonality and mean climate state, respectively. We find good agreement among the algorithms in the AR response to the changing orbital configuration, with a poleward shift in AR frequency that tracks seasonal poleward shifts in atmospheric water vapor and zonal winds. In the low CO2 simulation, the algorithms generally agree on the sign of AR changes, but there is substantial spread in their magnitude, indicating that mean-state changes lead to larger uncertainty. This disagreement likely arises primarily from differences between algorithms in their thresholds for water vapor and its transport used for identifying ARs. These findings warrant caution in ARDT selection for paleoclimate and climate change studies in which there is a change to the mean climate state, as ARDT selection contributes substantial uncertainty in such cases.

Atmospheric river, paleoclimate↗