SEARCH · Engineering Papers
Results for “Directional derivatives”
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.
Shortwave Array Spectroradiometer-Hemispheric (SAS-He): design and evaluation
A novel ground-based radiometer, referred to as the Shortwave Array Spectroradiometer-Hemispheric (SAS-He), is introduced. This radiometer uses the shadow-band technique to report total irradiance and its direct and diffuse components frequently (every 30 s) with continuous spectral coverage (350–1700 nm) and moderate spectral (~2.5 nm ultraviolet–visible and ~6 nm shortwave-infrared) resolution. The SAS-He's performance is evaluated using integrated datasets collected over coastal regions during three field campaigns supported by the US Department of Energy's Atmospheric Radiation Measurement (ARM) program, namely the (1) Two-Column Aerosol Project (TCAP; Cape Cod, Massachusetts), (2) Tracking Aerosol Convection Interactions Experiment (TRACER; in and around Houston, Texas), and (3) Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE; La Jolla, California). We compare (i) aerosol optical depth (AOD) and total optical depth (TOD) derived from the direct irradiance, as well as (ii) the diffuse irradiance and direct-to-diffuse ratio (DDR) calculated from two components of the total irradiance. As part of the evaluation, both AOD and TOD derived from the SAS-He direct irradiance are compared to those provided by a collocated Cimel sunphotometer (CSPHOT) at five (380, 440, 500, 675, 870 nm) and two (1020, 1640 nm) wavelengths, respectively. Additionally, the SAS-He diffuse irradiance and DDR are contrasted with their counterparts offered by a collocated multifilter rotating shadowband radiometer (MFRSR) at six (415, 500, 615, 675, 870, 1625 nm) wavelengths. Overall, reasonable agreement is demonstrated between the compared products despite the challenging observational conditions associated with varying aerosol loadings and diverse types of aerosols and clouds. For example, the AOD- and TOD-related values of root mean square error remain within 0.021 at 380, 440, 500, 675, 870, 1020, and 1640 nm wavelengths during the three field campaigns.
JetGP: A derivative enhanced Gaussian process library
Derivative enhanced Gaussian Processes (DEGPs) can significantly improve surrogate model accuracy over standard Gaussian Process (GP) formulations by incorporating derivative information. However, standard implementations scale poorly with dimension, limiting their use in high dimensional engineering problems. JetGP is a Python framework that unifies existing derivative enhanced GP methodologies into a single library and extends them to support arbitrary order derivative information. The library implements four complementary formulations: standard derivative enhanced Gaussian Processes (DEGP), directional DEGP (DDEGP), generalized directional DEGP (GDDEGP), and weighted DEGP (WDEGP). By unifying these approaches in a consistent interface with robust numerical implementations, JetGP enables practitioners to balance predictive accuracy and computational efficiency for high dimensional optimization, uncertainty quantification, and sensitivity analysis in engineering design.
Recommendations for developing, documenting, and distributing data products derived from NEON data
The National Ecological Observatory Network (NEON) provides over 180 distinct data products from 81 sites (47 terrestrial and 34 freshwater aquatic sites) within the United States and Puerto Rico. These data products include both field and remote sensing data collected using standardized protocols and sampling schema, with centralized quality assurance and quality control (QA/QC) provided by NEON staff. Such breadth of data creates opportunities for the research community to extend basic and applied research while also extending the impact and reach of NEON data through the creation of derived data products—higher level data products derived by the user community from NEON data. Derived data products are curated, documented, reproducibly-generated datasets created by applying various processing steps to one or more lower level data products—including interpolation, extrapolation, integration, statistical analysis, modeling, or transformations. Derived data products directly benefit the research community and increase the impact of NEON data by broadening the size and diversity of the user base, decreasing the time and effort needed for working with NEON data, providing primary research foci through the development via the derivation process, and helping users address multidisciplinary questions. Creating derived data products also promotes personal career advancement to those involved through publications, citations, and future grant proposals. However, the creation of derived data products is a nontrivial task. Here we provide an overview of the process of creating derived data products while outlining the advantages, challenges, and major considerations.
Excited-State Densities from Time-Dependent Density Functional Response Theory
While the variational principle for excited-state energies leads to a route to obtaining excited-state densities from time-dependent density functional theory, relatively little attention has been paid to the quality of the resulting densities in real space obtained with different exchange-correlation functional approximations or how nonadiabatic approximations developed for energies of states of double-excitation character perform for their densities. Here we derive an expression directly in real space for the excited-state density, which includes the case of nonadiabatic kernels and consequently is able, for the first time, to yield densities of states of double-excitation character. Under some well-defined simplifications, we compare the performance of the local-density approximation and exact-exchange approximation, which are in a sense at the opposite extremes of the fundamental functional approximations, on local and charge-transfer excitations in one-dimensional model systems and show that the dressed Time-Dependent Density Functional Theory (TDDFT) approach gives good densities of double excitations.
Hofmann Stability Charts Revisited for PIP-II: From Classical Theory to Assumption-Free and ML-Driven Maps
The Hofmann stability chart remains a standard for visualizing parametric resonances in space-charge–dominated linacs, but its use typically relies on non-oscillatory Vlasov dispersion relations with simplifying assumptions (continuous focusing, KV phase space, linear optics, limited transverse–longitudinal coupling). We revisit the chart for the PIP-II linac along three tracks. (1) We reproduce the conventional maps in the (νz/νx, νx/ν0x) plane for relevant εz/εx, providing a validated reference. (2) We remove key assumptions by deriving stability surfaces directly from multi-particle tracking with realistic lattice discreteness, RF defocusing, solenoid/quad optics, and bunched-beam dynamics; local tunes and early-time growth rates are estimated from envelope oscillations and projected to the same coordinates. These assumption-reduced maps recover the canonical stopbands while revealing shifts and broadenings driven by tune modulation, non-KV distributions, and transverse–longitudinal coupling at PIP-II intensities. (3) We train a compact machine-learning surrogate that emulates the growth surface from zero-current optics, tune depression, emittance ratio, bunching factor, and selected lattice descriptors, enabling rapid scans and online working-point selection. We compare the three representations on representative PIP-II sections and discuss implications for commissioning guard bands, resonance avoidance, and routine operations.
Nyström type exponential integrators for strongly magnetized charged particle dynamics
Solving for charged particle motion in electromagnetic fields (i.e. the particle pushing problem) is a computationally intensive component of particle-in-cell (PIC) methods for plasma physics simulations. This task is especially challenging when the plasma is strongly magnetized due numerical stiffness arising from the wide range of time scales between highly oscillatory gyromotion and long term macroscopic behavior. A promising approach to solve these problems is by a class of methods known as exponential integrators that can solve linear problems exactly and are A-stable. This work extends the standard exponential integration framework to derive Nyström-type exponential integrators that integrates the Newtonian equations of motion as a second-order differential equation directly. In particular, we derive second-order and third-order Nyström-type exponential integrators for strongly magnetized particle pushing problems. Numerical experiments show that the Nyström-type exponential integrators exhibit significant improvement in computation speed over the standard exponential integrators.
BLOC Site - Radar Wind Profiler / Derived Data Reformatted
A radar wind profiler measures the Doppler shift of electromagnetic energy scattered back from atmospheric turbulence and hydrometeors along 3-5 vertical and off-vertical point beam directions. Back-scattered signal strength and radial-component velocities are remotely sensed along all beam directions and combined to derive the horizontal wind field over the radar. These data are typically sampled and averaged hourly or sub-hourly (15-min) and usually have 60-m and/or 100-m vertical resolutions up to 4 km for the 915 MHz system. Both a high-resolution, lower height coverage mode and a low-resolution, higher height coverage mode are used to collect the data.
NANT Site - Radar Wind Profiler / Derived Data Reformatted
A radar wind profiler measures the Doppler shift of electromagnetic energy scattered back from atmospheric turbulence and hydrometeors along 3-5 vertical and off-vertical point beam directions. Back-scattered signal strength and radial-component velocities are remotely sensed along all beam directions and combined to derive the horizontal wind field over the radar. These data are typically sampled and averaged hourly or sub-hourly (15-min) and usually have 60-m and/or 100-m vertical resolutions up to 4 km for the 915 MHz system. Both a high-resolution, lower height coverage mode and a low-resolution, higher height coverage mode are used to collect the data.
Net-Zero Ethylene: On the Sustainability, Economics, and Scalability of Synthetic and Fossil Production Pathways
The ethylene industry has contributed over 260 million tons of CO 2 annually, warranting a more sustainable approach. The conversion of CO 2 and H 2 O into ethylene is an appealing technology capable of decoupling chemical production from fossil fuels. However, the large energy demand from this process can potentially lead to adverse environmental impacts. Here, in this article, we critically analyze the economic viability, environmental impact, and scalability of the conversion of CO 2 to ethylene via electrochemical reduction (CO 2 R) and compare this with those of CO 2 -neutral fossil routes utilizing carbon capture and direct air capture. Ethylene derived from CO 2 may be economically competitive under optimistic conditions; however, its large energy requirements pose environmental and scalability challenges. Meeting forecast 2050 ethylene demand using CO 2 R would require half of all electricity produced globally today, and, if powered by solar PV, may have greater CO 2 emissions than current petrochemical ethylene production, negating the purpose of this technology. Using Carbon Capture and Storage and Direct Air Capture to decarbonize petrochemical pathways would require roughly an order of magnitude less energy but would have disproportionate health and climate impacts. Lastly, the analysis highlights the importance of low-carbon energy sources to ensure sustainable CO 2 R ethylene production.
All-organic self-separating three-dimensionally nanoarchitected electrochemical energy storage devices
This work realizes a three-dimensionally (3D) nanoarchitected, all organic, "self-separating" lithium-ion electrochemical energy storage (EES) device that is cycled as a solid-state full cell. The device is enabled by a monolithic carbon anode with a co-continuous pore network, derived from the structure direction of resols by an ultra-large molar mass block copolymer (BCP), poly(styrene-block-2-dimethylaminoethyl methacrylate) (SA). Electropolymerization of a single-phase conductive and redox-active material, poly((2,3-dihydrothieno[3,4-b][1,4]dioxin-2-yl)methyl 9,10-dioxo-9,10-dihydroanthracene-2-carboxylate) (PAQEDOT), into the pore space provides the cathode of the cell. The device is electronically contacted to the relevant electrode network enabled by the co-continuous nature of each electrode. Electrochemical processing via cycling against external lithium in an electrolyte generates a solid electrolyte interphase (SEI) as a separator and lithiates the cell electrodes, after which the EES device is cycled in the solid state. While the full cell does not demonstrate high cyclability, the best full cell demonstrates a discharge capacity of 267 milliamp hours per gram (mAh/g). This work marks, to the best of knowledge of the authors, the first example of an all-organic materials derived 3D nanoarchitected EES device, as well as the first design of "self-separating" cell fabrication. Furthermore, generalization of the design to another co-continuous carbon form factor is demonstrated.
Life Cycle Analysis of Growing Canola for Biofuel Production in the United States
This study quantifies and compares the life cycle greenhouse gas (GHG) emissions of renewable diesel (RD), sustainable aviation fuel (SAF), and biodiesel (BD) produced from two U.S. canola production systems: 1) emerging intermediate winter canola, typically grown in double- or relay-cropping systems between the growing seasons of main crops, and 2) main canola, mostly spring canola but also including winter canola, which are grown as primary crops occupying the field for a full growing season. Using the Research and Development version of the Greenhouse gases, Regulated Emissions, and Energy use in Technologies (R&D GREET) model and the most up-to-date life cycle inventory data─field trial data for intermediate winter canola (>37,000 acres) and recent national survey data for spring canola─this life cycle analysis (LCA) estimates the direct emissions from canola cultivation and harvest, the conversion of canola into fuels, fuel transportation, and combustion. In addition, we account for market-mediated emissions associated with a scenario of 0.5 billion gallons per year of spring canola-based biofuels, including induced land use change (ILUC), induced other crop (nonfeedstock) production changes, and induced livestock production changes. For intermediate winter canola, these market-mediated effects were not modeled, as ILUC is expected to be negligible due to its integration into existing rotations, and data are currently insufficient to reliably quantify other market-mediated changes. The estimated life cycle direct emissions of RD/SAF derived from intermediate winter canola and main spring canola are about 32 and 33 g of CO2-equivalent per megajoule of fuel (g CO 2 e/MJ), respectively. Corresponding emissions for BD from intermediate winter canola and main spring canola are about 30 and 31 g of CO 2 e/MJ, respectively. Farming is the dominant emissions source for both canola systems, with intermediate winter canola and main spring canola emitting about 19 and 20 g of CO 2 e/MJ, respectively. ILUC and other induced changes increase emissions of main spring canola-derived RD/SAF and BD by about 18 and 17 g of CO 2 e/MJ, respectively. These results indicate that the GHG emissions of biofuels produced from the two canola systems may differ substantially due to the different land use dynamics of the systems.
Photoluminescence line shapes of nanocrystals: Contributions from first- and second-order vibronic couplings
Here, we present a microscopic, parameter-free approach for computing the photoluminescence spectra of a single semiconductor nanocrystal. The method derives exciton–phonon coupling directly from the semi-empirical pseudopotential framework and systematically incorporates both diagonal and off-diagonal exciton-phonon interactions, expanded to second-order in the phonon coordinates. The dipole–dipole correlation function was calculated using a Dyson expansion within the Kubo–Toyozawa formalism, enabling a consistent description of the role of pure dephasing and population transfer on the photoluminescence spectral features. Applied to CdSe/CdS core–shell nanocrystals, the approach quantitatively reproduces experimental photoluminescence spectra over a wide temperature range, revealing that quadratic phonon couplings account for nearly half of the homogeneous linewidth above ≈ 100−150 K, while off-diagonal couplings leading to exciton thermalization play only a minor role and only as T → 300 K.
Fast and accurate calculation of EXAFS Debye-Waller factors in UO2 using the dynamical matrix method
Theoretical modeling of bonding dynamics in metal oxides is required for predicting their thermal conductivity, catalytic activity, and mechanical properties. A primary challenge is the scarcity of experimental methods for validating theoretical predictions of these atomic-scale dynamics. This work presents a workflow that uses experimental extended x-ray absorption fine structure (EXAFS) data collected at high temperatures to validate an interatomic force field for uranium dioxide (UO2), an important model material. The validated force field is then used to drive computationally intensive molecular dynamics (MD) simulations and as input for the much faster dynamical matrix Debye-Waller (DMDW) method. The predicted values of the Debye-Waller factors from the DMDW calculations are in good agreement with those obtained from the MD simulations, with residual pair-specific differences attributable to quantum zero-point motion at low temperatures and lattice anharmonicity at high temperatures. We further show that theoretical EXAFS spectra constructed directly from DMDW-derived Debye-Waller factors reproduce the experimental data (at relatively low temperatures) with accuracy comparable to full MD-EXAFS, providing an additional validation of the choice of the potential. This study establishes a validated, rapid computational pathway for modeling bond dynamics, naturally incorporating quantum nuclear\\\\r\\\\nstatistics absent in classical simulations, which are essential for the mechanistic understanding of complex oxide materials.
SENTRA: A Modular Computational Graph Framework for Critical Mineral and Materials Supply Chains: Part I: Network Construction Latent-Quantity Estimation, and Temporal Graph Forecasting
Global supply chains for critical minerals and materials are complex, evolving networks of countries, products, production stages, and trade relationships. Existing analytical approaches are limited by fragmented data and static network representations that do not capture the dynamic production dependencies linking raw materials, intermediate products, and final goods across multiple countries. Trade and production statistics provide only a partial view of domestic production, inventories, and material flows, making it difficult to identify indirect sourcing pathways, hidden dependencies, and embedded foreign exposures. This paper introduces the Supply Chain Exposure Network Tracking and Risk Assessment (SENTRA) framework, a modular graph-based computational framework for constructing, analyzing, and forecasting dynamic supply chain networks. As the first paper in a three-part methodological series, it establishes the computational foundation of SENTRA by constructing a temporal attributed multi-relational graph whose nodes represent product–country pairs and whose edges encode observed trade and within-country value-chain relationships. Statistical estimation and constrained optimization recover latent production, final demand, and product input dependency coefficients while enforcing economic accounting constraints. Graph-derived exposure measures quantify direct, transshipment, value-chain, and multi-hop supply chain dependencies independently of the forecasting model. A temporal graph forecasting architecture based on a relational graph neural network then forecasts the evolution of the graph under mass-balance constraints with distribution-free conformal uncertainty quantification. Validation on the global aluminum supply chain shows that the learned graph representations recover economically meaningful supply chain structure, accurately forecast out-of-sample trade relationships, and produce well-calibrated prediction intervals. Subsequent papers apply this computational foundation to exposure assessment, disruption analysis, and scenario-based policy analysis, and extend the framework to multimaterial supply chain modeling and decision support.
Determining Magnetic Field Directions and Omni-directional Fluxes of Particles from STPSat-6 ZPS Plasma Measurements
This report summarizes our recent data derivation efforts on STPSat-6 ZPS plasma measurements, presenting final outcomes. We begin by outlining the methodology developed to determine local magnetic field directions from measured directional intensities of ~20 keV ions, leveraging symmetry in particles’ pitch-angle distributions, and the determined field directions over a storm period are shown with errors analyzed through comparison with NOAA GOES measurements. Next, we further assess the validity of the method and quantify its overall high performance while acknowledging certain caveats. Finally, given local magnetic field directions and ZPS data, we explain how omni-directional fluxes of particles can be derived from ZPS measurements with limited directional coverage, and estimate the related errors under several theoretical scenarios. The methods developed in this study can be applied for processing and augmenting ZPS data for the next, and the insights gained here can also inform instrument design in the future.
Mechanism of Molecular Polariton Decoherence in the Collective Light–Matter Couplings Regime
Molecular polaritons, the hybridization of electronic states in molecules with photonic excitation inside a cavity, play an important role in fundamental quantum science and technology. Understanding the decoherence mechanism of molecular polaritons is among the most significant fundamental questions. We theoretically demonstrate that hybridizing many molecular excitons in a cavity protects the overall quantum coherence from phononinduced decoherence. The polariton coherence time can be prolonged up to 100 fs with a realistic collective Rabi splitting and quality factor at room temperature, compared to the typical electronic coherence time which is around 15 fs. Our numerically exact simulations and analytic theory suggest that the dominant decoherence mechanism is the population transfer from the upper polariton state to the dark state manifold. Increasing the collective coupling strength will increase the energy gap between these two sets of states and thus prolong the coherence lifetime. We further derived valuable scaling relations that directly indicate how polariton coherence depends on the number of molecules, Rabi splittings, and light–matter detunings.
Global River Topology (GRIT): A Bifurcating River Hydrography
Existing global river networks underpin a wide range of hydrological applications but do not represent channels with divergent river flows (bifurcations, multi‐threaded channels, canals), as these features defy the convergent flow assumption that elevation‐derived networks (e.g., HydroSHEDS, MERIT Hydro) are based on. Yet, bifurcations are important features of the global river drainage system, especially on large floodplains and river deltas, and are also often found in densely populated regions. Here we developed the first raster and vector‐based Global RIver Topology that not only represents the tributaries of the global drainage network but also the distributaries, including multi‐threaded rivers, canals and deltas. We achieve this by merging a 30 m Landsat‐based river mask with elevation‐generated streams to ensure a homogeneous drainage density outside of the river mask for rivers narrower than approximately 30 m. Crucially, we employ the new 30 m digital terrain model, FABDEM, based on TanDEM‐X, which shows greater accuracy over the traditionally used SRTM derivatives. After vectorization and pruning, directionality is assigned by a series of elevation, flow angle and continuity approaches. The new global network and its attributes are validated using gauging stations, comparison with existing networks, and randomized manual checks. The new network represents 19.6 million km of streams and rivers with drainage areas greater than 50 km 2 and includes 67,495 bifurcations. With the advent of hyper‐resolution modeling and artificial intelligence, GRIT is expected to greatly improve the accuracy of many river‐based applications such as flood forecasting, water availability and quality simulations, or riverine habitat mapping.