Search NASA⌕ Search

SEARCH · Search NASA

Results for “Positivity preserving”

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 91 records · Page 5

A Computationally Efficient Linear Semi-Lagrangian Scheme for the Advection of Microphysical Variables in Cloud-Resolving Models

An efficient semi-Lagrangian (SL) numerical scheme for the advection of positive-definite microphysical variables in cloud-resolving models is proposed. The scheme is a linear combination of SL schemes of the first and second order. The proposed scheme is compared with two high-order, monotonic, nonlinear Eulerian schemes, in two 2D tests. In the first test, advection of a positive scalar field with large gradients is performed within complicated idealized steering flows. In the second test, a hail storm is simulated using the 2D Hebrew University Cloud Model, describing the interactions between 15 size distributions of different microphysical quantities (each distribution is described by 43 mass bins). The proposed scheme is computationally efficient, has a low numerical diffusion and a high level of mass conservation accuracy, and preserves the sum of multiple advected variables as well as the linear relationships between the advected scalar variables. Although the proposed SL scheme is not exactly positive definite, the negative values and their number, which may appear only in the case of very strong gradients, are negligible. The scheme produces results similar to those of nonlinear Eulerian advection schemes while reducing computation time by approximately a factor of 10.

54 ENVIRONMENTAL SCIENCES↗

Stabilized bases for high-order, interpolation semi-Lagrangian, element-based tracer transport

In a computational fluid model of the atmosphere, the advective transport of trace species, or tracers, can be computationally expensive. For efficiency, models often use semi-Lagrangian advection methods. High-order interpolation semi-Lagrangian (ISL) methods, in particular, can be extremely efficient, if the problem of property preservation specific to them can be addressed. Atmosphere models often use geometrically and logically nonuniform grids for efficiency and, as a result, element-based discretizations. Such grids and discretizations make stability a particular problem for ISL methods. Generally, high-order, element-based ISL methods that use the natural polynomial interpolant associated with a nodal finite-element discretization are unstable. Here, we derive new bases having order of accuracy up to nine, with positive nodal weights, that stabilize the element-based ISL method. We use these bases to construct the linear advection operator in the property-preserving Interpolation Semi-Lagrangian Element-based Transport (Islet) method. Then we discuss key software implementation details. Finally, we show performance results for the Energy Exascale Earth System Model's atmosphere dynamical core, comparing the original and new transport methods. These simulations used up to 27,600 Graphical Processing Units (GPU) on the Oak Ridge Leadership Computing Facility's Summit supercomputer.

97 MATHEMATICS AND COMPUTING↗

CAHS: Context-Aware Homology Search

Protein homology search is foundational to bioinformatics: it supports annotation transfer, structure/function inference, and evolutionary analysis over rapidly expanding sequence repositories (e.g., UniProtKB). Profile hidden Markov models (pHMMs), as implemented in HMMER, remain the most widely trusted approach because they provide statistically calibrated E-values; however, their gap behavior is fixed once a profile is trained, despite biological evidence that insertion/deletion tolerance varies across flexible loops and intrinsically disordered regions. We present CAHS (Context-Aware Homology Search), a lightweight query-time adapter for pHMM search that incorporates learned and biologically motivated signals without changing HMMER's downstream search pipeline or its calibrated E-value reporting. Given a query sequence, CAHS computes per-residue representations from a protein language model and a disorder predictor, maps these to profile coordinates, and modulates only match-state transition rows (gap-open and gap-extension probabilities) while preserving Plan7 constraints. We comprehensively evaluate CAHS across six structurally diverse protein families and multi-domain architectures against a 570k-sequence target corpus. CAHS expands detection capability, retrieving thousands of additional remote homologs at relaxed thresholds by maintaining alignment quality through flexible regions. For multi-domain proteins, context-aware modulation resolves 94% of fragmented alignments. Crucially, CAHS preserves hit-set invariance at stringent operating points (E<10-10), demonstrating increased statistical confidence without inflating false positives. Furthermore, sharper statistical distinction between homologs and background noise during early filter stages yields up to a 3.87× acceleration in end-to-end wall-clock time on high-performance computing clusters. Overall, CAHS illustrates a practical AI-for-science design pattern: augmenting a trusted probabilistic model with query-specific learned signals to improve interpretable, reproducible inference in data-rich biology.

Bhattaram, Swethasree [Georgia Institute of Techno↗

Preservation of a Dust Crystal as it Falls in an Afterglow Plasma

In an experiment, the power that sustains a plasma was extinguished, so that microspheres, which had been levitated, fell downward toward a lower electrode. At the beginning of their fall, the microspheres were self-organized with a crystalline structure. This structure was found to be preserved as the microspheres accelerated all the way to the lower electrode. Although microspheres had, in this afterglow plasma, large positive charges of 12,500 e, their interparticle repulsion was unable to significantly alter the crystalline arrangement of the microspheres, as they fell. After their impact on the lower electrode, the microspheres bounced upward, and only then was the crystalline structure lost.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Can increasing the size and flexibility of a molecule reduce decoherence and prolong charge migration?

Coherent superposition of electronic states, created by ionizing a molecule, can initiate ultrafast dynamics of the electron density. Correlation between nuclear and electron motions, however, typically dissipates the electronic coherence in only a few femtoseconds, especially in larger and more flexible molecules. We, therefore, use ab initio semiclassical dynamics to study decoherence in a sequence of analogous organic molecules of increasing size and find, surprisingly, that extending the carbon skeleton in propynal analogs slows down decoherence and prolongs charge migration. To elucidate this observation, we decompose the overall decoherence into contributions from individual vibrational modes and show that: 1) The initial decay of electronic coherence is caused by high- and intermediate-frequency vibrations via momentum separation of nuclear wavepackets evolving on different electronic surfaces. 2) At later times, the coherence disappears completely due to the increasing position separation in the low-frequency modes. 3) In agreement with another study, we observe that only normal modes that preserve the symmetry of the molecule induce decoherence. All together, we justify the enhanced charge migration by a combination of increased hole-mixing and the disappearance of decoherence contributions from specific vibrational modes—CO stretching in butynal and various H rockings in pentynal.

Science & Technology - Other Topics↗

Sensitivity of the pseudo-global warming method under flood conditions: a case study from the northeastern US

Abstract. Intensified extreme precipitation and corresponding floods are the most relevant consequences of climate change over the northeastern US (NEUS). To evaluate the impacts of climate change or certain climate perturbations on future extreme weather events which are dynamically similar to historic analogs, the pseudo-global warming (PGW) method has been frequently employed; however, this method lacks precise definition and guidelines, thus limiting its application. More specifically, three key questions related to the application of the PGW method remain unanswered: at what spatial scale should climate perturbations be applied? Among the different meteorological variables available, which ones should be perturbed? And will PGW projections vary significantly when different perturbations are applied? To address these questions, we examine the sensitivity and robustness of conclusions drawn from the PGW method over the NEUS by conducting multiple PGW experiments with varied perturbation spatial scales and choice of perturbed meteorological variables. The results show that the projections of precipitation and other essential variables at the regional mean scale are consistent across the PGW simulations, with a relative difference of much less than 10 %; however, different perturbation modifications can cause significant displacements of the storm events being simulated. Several previously assumed advantages of modifying only the temperature at regional mean scale, such as the preservation of geostrophic balance, do not appear to hold. Also, for these experiments, we find the regional mean perturbation produces a positive precipitation bias because it ignores the land–ocean warming contrast, which is a robust regional response to global warming. Overall, PGW experiments with perturbations from temperature or the combination of temperature and wind at the grid point scale are both recommended, depending on the research questions. The first approach can isolate the spatially dependent thermodynamic impact, and the latter incorporates both the thermodynamic and dynamic impacts.

54 ENVIRONMENTAL SCIENCES↗

ETG turbulence in a tokamak pedestal

This paper explores the fundamental characteristics of electron-temperature-gradient (ETG)-driven turbulence in the tokamak pedestal. The extreme gradients in the pedestal produce linear instabilities and nonlinear turbulence that are distinct from the corresponding ETG phenomenology in the core plasma. The linear system exhibits multiple (greater than ten) unstable eigenmodes at each perpendicular wave vector, representing different toroidal and slab branches of the ETG instability. Proper orthogonal decomposition of the nonlinear fluctuations reveals no clear one-to-one correspondence between the linear and nonlinear modes for most wave vectors. Moreover, nonlinear frequencies deviate strongly from those of the linear instabilities, with spectra peaking at positive frequencies, which is opposite the sign of the ETG instability. The picture that emerges is one in which the linear properties are preserved only in a narrow range of k-space. Outside this range, nonlinear processes produce strong deviations from both the linear frequencies and eigenmode structures. This is interpreted in the context of critical balance, which enforces alignment between the parallel scales and fluctuation frequencies. We also investigate the nonlinear saturation processes. We observe a direct energy cascade from the injection scale to smaller scales in both perpendicular directions. However, in the bi-normal direction, there is also nonlocal inverse energy transfer to larger scales. Neither streamers nor zonal flows dominate the saturation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fast automated energy changes at synchrotron radiation beamlines equipped with transfocator or focusing mirrors

Algorithms and procedures to fully automate retuning of synchrotron radiation beamlines over wide energy ranges are discussed. The discussion is based on the implementation at the National Institute of General Medical Sciences and the National Cancer Institute Structural Biology Facility at the Advanced Photon Source. When a user selects a new beamline energy, software synchronously controls the beamline monochromator and undulator to maintain the X-ray beam flux after the monochromator, preserves beam attenuation by determining a new set of attenuator foils, changes, as needed, mirror reflecting stripes and the undulator harmonic, preserves beam focal distance of compound refractive lens focusing by changing the In/Out combination of lenses in the transfocator, and, finally, restores beam position at the sample by on-the-fly scanning of either the Kirkpatrick–Baez mirror angles or the transfocator up/down and inboard/outboard positions. The sample is protected from radiation damage by automatically moving it out of the beam during the energy change and optimization.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Towards Three-Dimensional Neutron Imaging with Light-field Technology

A broad range of applications in nuclear safeguards and security can benefit from compact, fast neutron imagers with large angular acceptance. However, accurate 3D multi-vertex reconstruction in a monolithic detector remains a significant barrier to realization. Current approaches, such as the centroiding approach and the use of shadow masks, have not yet resulted in a successful demonstration in a light-starved environment. This project explored the use of cutting-edge technology—a light field camera—which inherently preserves both the spatial and directional information of incident light. Such technology can, in principle, resolve multi-vertex events and offer accurate position reconstruction in 3D with relatively simple readout electronics. Since the project's inception in May 2023, we have built an experimental setup to calibrate and characterize a commercial light-field camera (Lytro Illum). We assessed its 3D event reconstruction capability in a relatively low-light intensity environment by analyzing the cross-correlation of a series of 2D images of a characterized tunable light source taken at various distances. A sub-cm resolution, lower than typical neutron interaction separations in a compact scintillator volume, was observed with the LED light source, suggesting the promising capability of its nominal optics design to provide the adequate resolution required for a compact monolithic neutron directional detector. However, using a light-field-based readout for a neutron camera requires incorporating an ultra-low-noise light sensor, which is beyond the scope of this project. We also initiated the development of a 3D light-field-based reconstruction algorithm tailored to sparse scenes, as expected from particle interactions in a scintillator medium. Additionally, we demonstrated the capability of the algorithm to replicate the ground truth. Finally, we began the development of a neutron directional detector simulation to determine the performance criteria for light-field-based reconstruction that allows for good neutron directional reconstruction. This Feasibility Study has led to a successful follow-up project under the DNN R&D innovation portfolio starting May 2024.

42 ENGINEERING↗

Impact of storage and blending of algae and forest product residue on fuel blendstock production

Seasonal impacts on the production of algae biomass require blending with other feedstocks such as wood to maintain consistent annual conversion capacity. Idaho National Laboratory (INL) has developed a long-term storage strategy for algae biomass using ensiling, or anaerobic wet storage, that was tested for blends of algae and wood to stabilize the feedstock supply for processing to fuels. Additionally, blending biomass (algae and wood) leverages existing biomass storage approaches commonly used by the feed and forage industry. Earlier research also has demonstrated positive results that indicate ensiling is more effective and economical in normalizing biomass feedstock supply prior to conversions such as HTL (hydrothermal liquefaction) than preserving biomass by drying. By assessing impacts beginning with upstream logistic operations and proceeding through each conversion step, this work focuses on conversion of ensiled and blended biomass to fungible liquid transportation fuel blendstocks. This is accomplished by conversion of the blended biomass (62% Chlorella sp. blended with 38% loblolly pine forest product residues [FPR]) to bio-crude through HTL and subsequent upgrading through HT (hydrotreating) with a commercial refinery catalyst analogous to a refinery process to produce hydrocarbon fuel. In this work, results indicate that carbon retention and quality are preserved in both the bio-crude and upgraded fuel from the ensiled blend, indicating the potential of this approach for managing seasonal variations in algae biomass productivity.

09 BIOMASS FUELS↗

A performant energy-conserving particle reweighting method for Particle-in-Cell simulations

A new particle-based reweighting method is developed and demonstrated in the Aleph Particle-in-Cell with Direct Simulation Monte Carlo (PIC-DSMC) program. Novel splitting and merging algorithms ensure that modified particles maintain physically consistent positions and velocities. This method allows a single reweighting simulation to efficiently model plasma evolution over orders of magnitude variation in density, while accurately preserving energy distribution functions (EDFs). Demonstrations on electrostatic sheath and collisional rate dynamics show that reweighting simulations achieve accuracy comparable to fixed weight simulations with substantial computational time savings. This highly performant reweighting method is recommended for modeling plasma applications that require accurate resolution of EDFs or exhibit significant density variations in time or space.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The DOE CSSG: 2021 the Year in Review [Slides]

2021 accomplishments include completed taskings (2021-01, CSSG Review of Draft Revision to DOE-STD-1134 amd 2021-02, CSSG Review of FY22 NCSP Proposals) and a panel presentation at the 2021 ANS Winter Meeting. Projected taskings for 2022 include CSSG Guidance to NNSA on Contract Transition Regarding Criticality Safety, self-assessment of the CSSG, and a focus on improving the visibility and relevance of the CSSG. The CSSG is active and always seeking opportunities to assist. The CSSG is striving to be a strong, unified voice, bringing reason to and/or keeping reason in Criticality Safety. The CSSG is well positioned to continue their support role into the future.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

Antarctic permafrost processes and antiphase dynamics of cold-based glaciers in the McMurdo Dry Valleys inferred from 10 Be and 26 Al cosmogenic nuclides

Abstract. Soil and sediment mixing and associated permafrost processes are not widely studied or understood in the McMurdo Dry Valleys of Antarctica. In this study, we investigate the stability and depositional history of near-surface permafrost sediments to ∼ 3 m depth in the Pearse and lower Wright valleys using measured cosmogenic 10Be and 26Al depth profiles. In Pearse Valley, we estimate a minimum depositional age of ∼ 74 ka for the active layer and paleoactive-layer sediments (< 0.65 m). Combined depth profile modelling of 10Be and 26Al gives a depositional age for near-surface (< 1.65 m) permafrost in Pearse Valley of 180 +20/-40 ka, implying that the deposition of permafrost sediments predates MIS 5 advances of Taylor Glacier. Deeper permafrost sediments (> 2.09 m) in Pearse Valley are thus inferred to have a depositional age of > 180 ka. At a coastal, lower-elevation site in neighbouring lower Wright Valley, 10Be and 26Al depth profiles from a second permafrost core exhibit near-constant concentrations with depth and indicate the sediments are either vertically mixed after deposition or sufficiently young so that post-depositional nuclide production is negligible relative to inheritance. 26Al/10Be concentration ratios for both depth profiles range between 4.0 and 5.2 and are all lower than the nominal surface production rate ratio of 6.75, indicating that prior to deposition, these sediments experienced complex, yet similar, exposure–burial histories. Assuming a single-cycle exposure–burial scenario, the observed 26Al/10Be ratios are equivalent to a total minimum exposure–burial history of ∼ 1.2 Myr. In proximity to the depth profile core site, we measured cosmogenic 10Be and 26Al in three granite cobbles from thin, patchy drift (Taylor 2 Drift) in Pearse Valley to constrain the timing of retreat of Taylor Glacier. Assuming simple continuous exposure, our minimum, zero-erosion exposure ages suggest Taylor Glacier partially retreated from Pearse Valley no later than 65–74 ka. The timing of retreat after 65 ka and until the Last Glacial Maximum (LGM) when Taylor Glacier was at a minimum position remains unresolved. The surface cobble ages and permafrost processes reveal Taylor Glacier advances during MIS 5 were non-erosive or mildly erosive, preserving the underlying permafrost sediments and peppering boulders and cobbles upon an older, relict surface. Our results are consistent with U/Th ages from central Taylor Valley and suggest changes in moisture delivery over Taylor Dome during MIS 5e, 5c, and 5a appear to be associated with the extent of the Ross Ice Shelf and sea ice in the Ross Sea. These data provide further evidence of antiphase behaviour through retreat of a peripheral lobe of Taylor Glacier in Pearse Valley, a region that was glaciated during MIS 5. We suggest a causal relationship of cold-based glacier advance and retreat that is controlled by an increase in moisture availability during retreat of sea ice and perhaps the Ross Ice Shelf, as well as, conversely, a decrease during times of sea ice and Ross Ice Shelf expansion in the Ross Sea.

58 GEOSCIENCES↗

Dynamics of oil–water interface at the beginning of the ultrasonic emulsification process

A lot of effort has been dedicated in recent years towards understanding the basics of cavitation induced emulsification, mainly in the form of single cavitation bubbles. Regarding bulk acoustic emulsification, a lot less research has been done. In our here presented work we utilize advanced high-speed observation techniques in visible light and X-Rays to build upon that knowledge and advance the understanding of bulk emulsion preparation. During research we discovered that emulsion formation has an acute impact on the behavior of the interface and more importantly on its position relative to the horn, hence their interdependence must be carefully studied. We did this by observing bulk emulsification with 2 cameras simultaneously and corroborating these measurements with observation under X-Rays. Since the ultrasonic horns location also influences interface behavior, we shifted its initial position to different locations nearer to and further away from the oil–water interface in both phases. We found that a few millimeters distance between the horn and interface is not enough for fine emulsion formation, but that they must be completely adjacent to each other, with the horn being located inside the oil–water interface. We also observed some previously undiscovered phenomena, such as the splitting of the interface to preserve continuous emulsion formation, climbing of the interface up the horn and circular interface protrusions towards the horn forming vertical emulsion streams. Interestingly, no visible W/O emulsion was ever formed during our experiments, only O/W regardless of initial horn position.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Automatic Determination of the Weak-Beam Condition in Dark Field X-ray Microscopy

Mechanical properties in crystals are strongly correlated to the arrangement of 1D line defects, termed dislocations. Recently, Dark field X-ray Microscopy (DFXM) has emerged as a new tool to image and interpret dislocations within crystals using multidimensional scans. However, the methods required to reconstruct meaningful dislocation information from high-dimensional DFXM scans are still nascent and require significant manual oversight (i.e. supervision). In this work, we present a new relatively unsupervised method that extracts dislocation-specific information (features) from a 3D dataset (x, y, Φ) using Gram-Schmidt orthogonalization to represent the large dataset as an array of 3-component feature vectors for each position, corresponding to the weak-beam conditions and the strong-beam condition. Finally, this method offers key opportunities to significantly reduce dataset size while preserving only the crystallographic information that is important for data reconstruction.

36 MATERIALS SCIENCE↗

Electronic Structure Modification of Rectangular Phosphorene Quantum Dots Via Edge Passivation

Phosphorene is an elemental two-dimensional material that shows great promise as a metal-free photocatalyst owing to its high carrier mobility and direct band gap. Doping and edge functionalization of phosphorene can be exploited to tune its band gap and alter other properties. Herein, the electronic properties of a series of edge-passivated monolayer phosphorene quantum dots are investigated using an array of functional groups (OH, SH, NH 2 , SCH 3 , OCN, CN, and Cl) that are known to have electron-donating or electron-withdrawing capabilities. By employing density functional theory calculations, we demonstrate that edge functionalization can perturbatively shift the HOMO (highest occupied molecular orbital) and LUMO (lowest unoccupied molecular orbital) energy levels while simultaneously preserving the geometric structure and other favorable properties of pristine phosphorene. We use the average overall electrostatic potential of these systems to successfully predict the relative positioning of the HOMO and LUMO energy levels while presenting qualitative explanations on potential jumps introduced by the dipole layers formed due to edge functionalization. Furthermore, our results suggest that selective functionalization permits considerable control over the band edge positions, which could be utilized for applications in energy and optoelectronic devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exploring Atomic-Scale Interactions at the Interface of Reducible Oxide and Ruthenium Nanocatalyst for Ammonia Decomposition

Restructuring atomic-scale interfaces between noble metal nanoparticles and metal oxides provides a promising approach to enhancing catalytic properties. In this study, Ru nanoparticles were supported on CeO 2 , ZrO 2 , and HfO 2 via pyrolysis of MOFs under oxygen-suppressed high-temperature conditions. Despite its nearly ideal particle size (∼2.3 nm), the Ru–CeO 2 /C catalyst exhibited the lowest ammonia decomposition activity. X-ray absorption spectroscopy and DFT calculations revealed that electron transfer from Ru to CeO 2 formed positively charged Ru and partially reduced Ce 4+ , weakening the catalytic performance. In contrast, oxygen-deficient ZrO 2 and HfO 2 support donated electrons to Ru, preserving its metallic state. These findings demonstrate that the reducibility of oxide supports governs the direction and magnitude of interfacial charge transfer, directly tuning catalytic behavior. This study provides insights into the design of oxide-supported Ru catalysts for ammonia decomposition.

Catalysts↗

Tunneling between Multiple Histories as a Solution to the Information Loss Paradox

The information loss paradox associated with black hole Hawking evaporation is an unresolved problem in modern theoretical physics. In a recent brief essay, we revisited the evolution of the black hole entanglement entropy via the Euclidean path integral (EPI) of the quantum state and allow for the branching of semi-classical histories along the Lorentzian evolution. We posited that there exist at least two histories that contribute to EPI, where one is an information-losing history, while the other is an information-preserving one. At early times, the former dominates EPI, while at the late times, the latter becomes dominant. By doing so, we recovered the essence of the Page curve, and thus, the unitarity, albeit with the turning point, i.e., the Page time, much shifted toward the late time. In this full-length paper, we fill in the details of our arguments and calculations to strengthen our notion. One implication of this modified Page curve is that the entropy bound may thus be violated. We comment on the similarity and difference between our approach and that of the replica wormholes and the islands’ conjectures.

79 ASTRONOMY AND ASTROPHYSICS↗