SEARCH · Engineering Papers
Results for “sequential analysis”
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.
A Strong Gravitational Lens Is Worth a Thousand Dark Matter Halos: Inference on Small-scale Structure Using Sequential Methods
Strong gravitational lenses are a singular probe of the Universe’s small-scale structure—they are sensitive to the gravitational effects of low-mass (<10 10 M ⊙ ) halos even without a luminous counterpart. Recent strong-lensing analyses of dark matter structure rely on simulation-based inference (SBI). Modern SBI methods, which leverage neural networks as density estimators, have shown promise in extracting the halo-population signal. However, it is unclear whether the constraints from these models are limited by the methodology or the data. In this study, we introduce an accelerator-optimized simulation pipeline that can generate lens images with realistic subhalo populations in milliseconds. Leveraging this simulator, we identify the main limitation of our fiducial SBI analysis: training set size. We then adopt a sequential neural posterior estimation (SNPE) approach, allowing us to refine the training distribution to align with the observed data. Using only one-fifth as many mock Hubble Space Telescope images, SNPE matches the constraints on the low-mass halo population produced by our best nonsequential model. Our experiments suggest that an over 3 order-of-magnitude increase in training set size and GPU hours would be required to achieve an equivalent result without sequential methods. While the full potential of the existing lens sample remains to be explored, the notable improvement in constraining power enabled by our sequential approach highlights that current constraints are limited primarily by methodology and not the data itself. Moreover, our results emphasize the need to treat training set generation and model optimization as interconnected stages of any cosmological analysis using SBI.
The advancing wave front on a sloping channel covered by a rod canopy following an instantaneous dam break
The drag coefficient Cd for a rigid and uniformly distributed rod canopy covering a sloping channel following the instantaneous collapse of a dam was examined using flume experiments. The measurements included space x and time t high resolution images of the water surface h(x, t) for multiple channel bed slopes So and water depths behind the dam Ho along with drag estimates provided by sequential load cells. Using these data, an analysis of the Saint-Venant equation (SVE) for the front speed was conducted using the diffusive wave approximation. An inferred Cd=0.4 from the h(x, t) data near the advancing front region, also confirmed by load cell measurements, is much reduced relative to its independently measured steady-uniform flow case. This finding suggests that drag reduction mechanisms associated with transients and flow disturbances are more likely to play a dominant role when compared to conventional sheltering or blocking effects on Cd examined in uniform flow. The increased air volume entrained into the advancing wave front region as determined from an inflow–outflow volume balance partly explains the Cd reduction from unity.
Per-phase Phasor Modeling of GFL and GFM Inverters for Distribution System Dynamic Studies
A significant increase in inverter-based distributed energy resources (DERs) is expected in the near future. Existing electromagnetic models and sequential electromechanical models of DERs supporting dynamic analysis cannot be extended to largescale three-phase distribution systems usually with unbalanced construction. As such, there is a need for appropriate industry models that can support distribution system dynamic studies. This paper describes the latest phasor models of grid-following (GFL) and grid-forming (GFM) inverters that can capture distribution system dynamics under grid disturbances. The overall modeling performance is examined and compared with results from other platforms such as controller hardware-in-the-loop (CHIL), MATLAB/Simulink and PSCAD with electromagnetic transient (EMT) and/or phasor simulation capabilities under balanced and unbalanced faults.
An Integrated Framework for Memory-Centric Analysis: From Trace Collection to Co-Design
The memory wall phenomenon—where advances in processor performance significantly outpace those in memory subsystems—poses a fundamental challenge for contemporary computing systems. In memory-bound applications, memory subsystem behavior dominates performance, yet existing analysis approaches present significant limitations: detailed microarchitectural simulators require days to weeks to simulate modest workloads; hardware performance counters provide only aggregate statistics that obscure temporal and spatial access patterns; and scaled simulation approaches face challenges in capturing certain behaviors that emerge at larger scales. These limitations reflect a processor-centric design philosophy increasingly misaligned with memory-bound workloads where detailed understanding of memory access patterns, cache hierarchy interactions, and contention is critical for effective optimization. This paper presents an integrated framework for memory-centric analysis that enables effective hardware-software co-design. We describe practical trace collection techniques, including hardware-assisted processor tracing with minimal overhead and portable software-based instrumentation with statistical sampling. We present multi-perspective analysis methods that examine memory behavior from temporal, sequential, spatial, and relational viewpoints, revealing distinct optimization opportunities invisible in aggregate metrics. We detail an architectural modeling framework that uses sampled traces with temporal interpolation and confidence-based filtering to evaluate cache and memory configurations. Evaluation on representative benchmarks demonstrates that this framework achieves practical accuracy (L2 cache errors of 2.64\%, confidence-filtered L3 errors of 9.92\%, bandwidth errors of 7.33\%) while providing substantial speedup (26.8×) over cycle-accurate simulation, enabling rapid design space exploration. We demonstrate how this integrated framework enables systematic identification of both hardware optimizations (memory controller tuning, bank partitioning, NUMA configuration) and software optimizations (data layout restructuring, prefetching strategies, memory-aware scheduling). Through this comprehensive treatment of the memory-centric analysis pipeline—from trace collection through architectural modeling to co-design application—we provide researchers and practitioners with practical techniques for addressing memory bottlenecks in contemporary computing systems.
Benchmarking Autonomous Scattering Experiments Illustrated on TAS
With the advancement of artificial intelligence and machine learning methods, autonomous approaches are recognized to have great potential for performing more efficient scattering experiments. In our view, it is crucial for such approaches to provide thorough evidence about respective performance improvements in order to increase acceptance within a scientific community. Therefore, we propose a benchmarking procedure designed as a cost-benefit analysis that is applicable to any scattering method sequentially collecting data during an experiment. For a given approach, the performance assessment is based on how much benefit, given a certain cost budget, it is able to acquire in predefined test cases. Different approaches thus get a chance for comparison and can make their advantages explicit and visible. Key components of the procedure, i.e., cost measures, benefit measures, and test cases, are made precise for the setting of three-axes spectrometry (TAS) as an illustration. Finally, we discuss neglected aspects and possible extensions for the TAS setting and comment on the procedure’s applicability to other scattering methods. A Python implementation of the procedure to simplify its utilization by interested researchers from the field is also provided.
Intricate Genetic Programs Controlling Dormancy in Mycobacterium tuberculosis
Mycobacterium tuberculosis (MTB) displays the remarkable ability to transition in and out of dormancy, a hallmark of the pathogen’s capacity to evade the immune system and exploit susceptible individuals. Uncovering the gene regulatory programs that underlie the phenotypic shifts in MTB during disease latency and reactivation has posed a challenge. We develop an experimental system to precisely control dissolved oxygen levels in MTB cultures in order to capture the transcriptional events that unfold as MTB transitions into and out of hypoxia-induced dormancy. Using a comprehensive genome-wide transcription factor binding map and insights from network topology analysis, we identify regulatory circuits that deterministically drive sequential transitions across six transcriptionally and functionally distinct states encompassing more than three-fifths of the MTB genome. The architecture of the genetic programs explains the transcriptional dynamics underlying synchronous entry of cells into a dormant state that is primed to infect the host upon encountering favorable conditions.
High manganese redox variability and manganate predominance in temperate soil profiles as determined by X-ray absorption spectroscopy
Manganese speciation is a key to understanding the fate of contaminants, nutrients, and organic matter in soils. To date, quantification of Mn species in bulk soils has been performed mainly by sequential extraction methods and rarely supported by spectroscopic analysis. In order to obtain quantitative information on the Mn species inventory of soils, we investigated 46 soil horizons (<2-mm fraction, 45.1–2,280 mg/kg Mn) of nine typical Central European soils (Cambisols, Chernozems, Luvisols, Podzol, Stagnosol) by chemical Mn analyses and Mn K-edge X-ray absorption spectroscopy, and related speciation results to major soil properties. Amounts of Mn 2+ , Mn 3+ , and Mn 4+ and the average oxidation state of Mn were evaluated by linear combination fitting (LCF) of X-ray absorption near edge structure (XANES) spectra. Additionally, we used extended X-ray absorption fine structure (EXAFS) spectroscopy to identify and quantify major Mn species. For this, EXAFS spectra of 20 organic and mineral soil samples from five soils (Cambisols, Chernozem, Luvisol) were analyzed by LCF and shell fitting. XANES analyses revealed a high Mn redox variability in organic surface layers, with Mn 2+ being most abundant (≤100%, $\overline{x}$ = 54%), followed by Mn 3+ (≤80%, $\overline{x}$ = 32%) and Mn 4+ (≤55%, $\overline{x}$ = 14%). Mineral soil horizons contained significantly less Mn 2+ (≤56%, $\overline{x}$ = 23%), about equal quantities of Mn 3+ (≤68%, $\overline{x}$ = 31%), and were enriched in Mn 4+ (≤89%, $\overline{x}$ = 46%). EXAFS analyses implied the presence of six major Mn species groups: manganates, organically complexed Mn, Mn(III) oxyhydroxides, silicate-bound Mn, Mn oxides without tunnel- or layer structure, and physisorbed Mn. In litter horizons, Mn was mainly present in organic complexes (58–91%, $\overline{x}$ = 78%) and as physisorbed Mn (≤15%), but individual horizons also comprised manganates, Mn(III) oxyhydroxides, and silicate-bound Mn. Manganates, likely mixtures of phyllomanganates with hexagonal layer symmetry and tectomanganates, dominated in all mineral soil horizons (37–94%, $\overline{x}$ = 67%). Correlation analysis showed that manganates dissolve completely during dithionite-citrate and acid ammonium oxalate extractions, and suggested that Mn 4+ -rich manganates preferentially form under less acidic soil conditions, partly by oxidation of organically complexed Mn(II), and that they are enriched in the soil clay fraction. Mineral soil horizons also contained minor quantities of organically complexed Mn (≤39%, $\overline{x}$ = 11%), silicate-bound Mn (≤30%, $\overline{x}$ = 8%), Mn(III) oxyhydroxides (≤37%, $\overline{x}$ = 7%), Mn oxides without tunnel- or layer structure (≤18%, $\overline{x}$ = 5%), and physisorbed Mn (≤14%, $\overline{x}$ <1%). The detection of Mn(III) oxyhydroxides such as feitknechtite (β-MnOOH) or groutite (α-MnOOH) as well as the spinel hausmannite (Mn 3 O 4 ) in acidic soils is remarkable, since their formation is normally linked to neutral or alkaline pH conditions. Minor contributions of silicate-bound Mn indicate the release of Mn from primary minerals already at early stages of soil formation, and low concentrations of physisorbed Mn suggest that exchangeable Mn is rapidly converted into manganates in oxic soils. In conclusion, the predominance of manganates in mineral soils has far-reaching implications for the functioning of soils and biogeochemical element cycles, as these minerals play an important role in metal binding, plant nutrition, and redox-related processes.
Competition between Solvation and Intramolecular Hydrogen-Bonding in Microsolvated Protonated Glycine and β-Alanine
Infrared predissociation (IRPD) spectroscopy is used to reveal and compare the microsolvation motifs of GlyH + (H 2 O) n and β-AlaH + (H 2 O) n . The chemical structure of these amino acids differ only in the length of the carbon chain connecting the amine and carboxyl terminals, which nonetheless leads to a significant difference in the strength of the intramolecular C=O···H-N hydrogen bond in the unsolvated ions. This difference makes them useful in our studies of the competition between solvation and internal hydrogen bonding interactions. Analysis of the IRPD results reveals that the sequential addition of water molecules leads to similar effects on the intramolecular interaction in both GlyH + (H 2 O) n and β-AlaH + (H 2 O) n . Solvation of the -NH 3 + group leads to a weakening of the C=O···H-N hydrogen bond, while solvation of the carboxyl -OH leads to a strengthening of this bond. Additionally, we have found that for β-AlaH + , the addition of a H 2 O to the second solvation shell can still influence the strength of the C=O···H-N hydrogen bonding interaction. Lastly because the C=O···H-N interaction in β-AlaH + is stronger than that in GlyH + , more solvent molecules are needed to sufficiently weaken the intramolecular hydrogen bond such that isomers without this bond begin to be energetically competitive; this occurs at n=5 for β-AlaH + and n=1 for GlyH + .
Sequential Learning of Active Subspaces
This study shows that in recent years, active subspace methods (ASMs) have become a popular means of performing subspace sensitivity analysis on black-box functions. Naively applied, however, ASMs require gradient evaluations of the target function. In the event of noisy, expensive, or stochastic simulators, evaluating gradients via finite differencing may be infeasible. In such cases, often a surrogate model is employed, on which finite differencing is performed. When the surrogate model is a Gaussian process (GP), we show that the ASM estimator is available in closed form, rendering the finite-difference approximation unnecessary. We use our closed-form solution to develop acquisition functions focused on sequential learning tailored to sensitivity analysis on top of ASMs. We also show that the traditional ASM estimator may be viewed as a method of moments estimator for a certain class of GPs. We demonstrate how uncertainty on GP hyperparameters may be propagated to uncertainty on the sensitivity analysis, allowing model-based confidence intervals on the active subspace. Our methodological developments are illustrated on several examples.
Arabidopsis GALACTURONOSYLTRANSFERASE (GAUT) 1 synthesizes a homogalacturonan tightly bound to the cell wall and required for cell expansion
Arabidopsis GALACTURONOSYLTRANSFERASE1 (GAUT1) synthesizes homogalacturonan (HG), the most abundant pectin in growing plant cells. GAUT1 has the greatest in vitro enzyme activity of the six confirmed Arabidopsis HG biosynthetic GAUTs, but its biological activity remains elusive. Here we show that Arabidopsis GAUT1 homozygous mutants have a severe dwarfed seedling phenotype, survive several weeks as 2 to 3 mm seedlings, and have severely reduced shoot and root growth and hypocotyl epidermal, cortex and endodermal cell size. gaut1-1 pollen tubes are shorter than WT with increased bursting. Complementation of homozygous gaut1-1 with GAUT1 coding sequence driven by the GAUT1 promoter restored WT-like growth. The extreme dwarf phenotype of homozygous gaut1-1 seedlings precluded their use for detailed cell wall analysis, thus suspensions cultures were produced from callus generated from mutant and WT seedlings. Homozygous gaut1-1 suspension cells were smaller than WT with ∼30% reduced wall GalA content compared to WT. Sequential extraction of the walls with increasingly harsh solvents and sugar composition analysis revealed reduced GalA content in only the 4M KOH post-chlorite fraction, indicating that GAUT1-synthesized HG was held tightly in the wall by direct or indirect hydrogen bonding and/or oxidation-sensitive linkages. Treatment of wall fractions with endopolygalacturonase to hydrolyze HG and gel electrophoretic separation of hydrolysates exposed an HG-associated doublet band markedly downregulated in the homozygous gaut1-1 4M KOH post-chlorite fraction and to a lesser extent in 4M KOH and sodium chlorite fractions. NMR analysis identified the band as rhamnogalacturonan (RG)-II. Super resolution microscopy using anti-HG antibodies showed that, compared to WT, the homozygous gaut1-1 hypocotyl epidermal and callus cells had reduced content and length of HG nanofilaments, HG fibers associated with cell expansion in Arabidopsis. The results demonstrate that GAUT1-synthesized HG resides in a tightly-cell-wall-bound, RG-II-containing polymer required for HG nanofilament formation and seedling cell expansion.
Microalgae Conversion to Biofuels and Biochemical via Sequential Hydrothermal Liquefaction (SEQHTL) and Bioprocessing: 2020 State of Technology
A preliminary techno-economic analysis (TEA) was developed for the fiscal year 2020 state of technology (SOT) assessment to evaluate the benefits and risks for a large-scale microalga hydrothermal liquefaction (HTL) system based on most recent testing results. The focus of the study is directed toward the conversion system, which consists of five processes: two-stage sequential HTL (SEQHTL), biocrude upgrading to final fuels, bioprocessing for co-product generation, hydrogen generation, and steam cycle. In this system, algae biomass with corn stover supplement during the lower algae productivity seasons (winter, fall, and spring) to match the maximum algae seasonal production rate in summer is employed to maintain a constant plant capacity in all the seasons. Algae only (summer season) or algae/corn stover blended feedstock (other seasons) are sent to a two-stage SEQHTL process. In stage I, the carbohydrates in the feedstock are extracted and separated from the residual solid. The residual solid from stage I is further converted to biocrude in the SEQHTL stage II step. The biocrude is upgraded to final fuel products in an upgrading process. The extract stream from HTL stage I is sent to the bioprocessing section for co-product generation via fermentation of carbohydrate. Lactic acid (LA) is assumed to be the co-product based on current bioprocessing testing results.
Two- and three-body fragmentation of multiply charged tribromomethane by ultrafast laser pulses
We investigate the two- and three-body fragmentation of tribromomethane (bromoform, CHBr 3 ) resulting from multiple ionization by 28-femtosecond near-infrared laser pulses with a peak intensity of 6 × 10 14 W cm -2 . The analysis focuses on channels consisting exclusively of ionic fragments, which are measured by coincidence momentum imaging. The dominant two-body fragmentation channel is found to be Br + + CHBr 2 +. Weaker HBr + + CBr 2 +, CHBr + + Br 2 +, CHBr 2+ + Br 2 +, and Br + + CHBr 2 2+ channels, some of which require bond rearrangement prior to or during the fragmentation, are also observed. The dominant three-body fragmentation channel is found to be Br + + Br + + CHBr + . This channel includes both concerted and sequential fragmentation pathways, which we identify using the native frames analysis method. We compare the measured kinetic energy release and momentum correlations with the results of classical Coulomb explosion simulations and discuss the possible isomerization of CHBr 3 to BrCHBr–Br (iso-CHBr 3 ) prior to the fragmentation.
Microscale mechanical modeling of deformable geomaterials with dynamic contacts based on the numerical manifold method
Abstract Micromechanical modeling of geomaterials is challenging because of the complex geometry of discontinuities and potentially large number of deformable material bodies that contact each other dynamically. In this study, we have developed a numerical approach for micromechanical analysis of deformable geomaterials with dynamic contacts. In our approach, we detect contacts among multiple blocks with arbitrary shapes, enforce different contact constraints for three different contact states of separated, bonded, and sliding, and iterate within each time step to ensure convergence of contact states. With these features, we are able to simulate the dynamic contact evolution at the microscale for realistic geomaterials having arbitrary shapes of grains and interfaces. We demonstrate the capability with several examples, including a rough fracture with different geometric surface asperity characteristics, settling of clay aggregates, compaction of a loosely packed sand, and failure of an intact marble sample. With our model, we are able to accurately analyze (1) large displacements and/or deformation, (2) the process of high stress accumulated at contact areas, (3) the failure of a mineral cemented rock samples under high stress, and (4) post-failure fragmentation. The analysis highlights the importance of accurately capturing (1) the sequential evolution of geomaterials responding to stress as motion, deformation, and high stress; (2) large geometric features outside the norms (such as large asperities and sharp corners) as such features can dominate the micromechanical behavior; and (3) different mechanical behavior between loosely packed and tightly packed granular systems.
Metallic Technetium Sequestration in Nickel Core/Shell Microstructure during Fe(OH)2 Transformation with Ni doping
This study investigates the impacts of Ni doping on technetium-99 (Tc) sequestration in aqueous solutions through transformation of Fe(OH)2(s) to iron spinel (magnetite) under alkaline conditions. Extensive solid characterization was performed for the mineral phases produced, as well as the Tc/Ni speciation and distribution within these phases. X-ray diffraction results show that iron spinel was the dominant mineral product without detectable Ni incorporation. The doped Ni ions mainly precipitated as fine Fe/Ni oxide/hydroxide particles, including strongly reduced nanometer?sized spheroidal Ni-rich and metallic Ni phases. High-resolution analytical scanning transmission electron microscopy using energy dispersive X-ray spectroscopy and electron energy loss spectroscopy on the produced solid samples (focused ion beam-prepared specimens) revealed three Tc distribution domains dominated by nanocrystals and, especially, a Tc-rich metallic phase. Instances of metallic Tc were specifically found in spheroidal, Ni-rich and metallic nanoparticles exhibiting a core/shell microstructure that suggests strong reduction and sequential precipitation of Ni-Tc-Ni. Mass balance analysis showed nearly 100% Tc removal from the 4.8 × 10-4 M Tc solutions. The finding of the metallic Tc encapsulation indicates that Tc sequestration through Ni-doped Fe(OH)2(s)?to?iron spinel transformation process likely provides an alternative treatment pathway for Tc removal and could be combined into further waste treatment approaches.
In situ synchrotron diffraction and modeling of non-equilibrium solidification of a MnFeCoNiCu alloy
The solidification mechanism and segregation behavior of laser-melted Mn 35 Fe 5 Co 20 Ni 20 Cu 20 was firstly investigated via in situ synchrotron x-ray diffraction at millisecond temporal resolution. The transient composition evolution of the random solid solution during sequential solidification of dendritic and interdendritic regions complicates the analysis of synchrotron diffraction data via any single conventional tool, such as Rietveld refinement. Therefore, a novel approach combining a hard-sphere approximation model, thermodynamic simulation, thermal expansion measurement and microstructural characterization was developed to assist in a fundamental understanding of the evolution of local composition, lattice parameter, and dendrite volume fraction corresponding to the diffraction data. This methodology yields self-consistent results across different methods. Via this approach, four distinct stages were identified, including: (I) FCC dendrite solidification, (II) solidification of FCC interdendritic region, (III) solid-state interdiffusion and (IV) final cooling with marginal diffusion. It was found out that in Stage I, Cu and Mn were rejected into liquid as Mn 35 Fe 5 Co 20 Ni 20 Cu 20 solidified dendritically. During Stage II, the lattice parameter disparity between dendrite and interdendritic region escalated as Cu and Mn continued segregating into the interdendritic region. After complete solidification, during Stage III, the lattice parameter disparity gradually decreases, demonstrating a degree of composition homogenization. The volume fraction of dendrites slightly grew from 58.3 to 65.5%, based on the evolving composition profile across a dendrite/interdendritic interface in diffusion calculations. Postmortem metallography further confirmed that dendrites have a volume fraction of 64.7 ± 5.3% in the final microstructure.
Unsupervised Segmentation and Clustering Workflow for Efficient Processing of 4D-STEM and 5D-STEM Data
Four-dimensional scanning transmission electron microscopy (4D-STEM) enables mapping of diffraction information with nanometer-scale spatial resolution, offering detailed insight into local structure, orientation, and strain. However, as data dimensionality and sampling density increase, particularly for in situ scanning diffraction experiments (5D-STEM), robust segmentation of structurally consistent behavior across sequential measurements becomes essential for efficient and physically meaningful analysis. Here, we introduce a clustering framework that identifies crystallographically distinct domains from 4D-STEM datasets. By using local diffraction-pattern similarity as a metric, the method extracts closed contours delineating spatially contiguous regions. This approach produces cluster-averaged diffraction patterns that improve signal quality while reducing data volume by orders of magnitude, enabling rapid and accurate orientation, phase, and strain mapping. We demonstrate the applicability of this approach to in situ liquid-cell 4D-STEM data of gold nanoparticle growth. Our method provides a scalable and generalizable route for spatially coherent segmentation, data compression, and quantitative structure–strain mapping across diverse 4D-STEM modalities. The full analysis code and example workflows are publicly available to support reproducibility and reuse.
Custom Accessors: Enabling Scalable Data Ingestion, (Re-)Organization, and Analysis on Distributed Systems
The emerging class of high velocity and high volume data analytic workflows comprise interwoven data ingestion, organization, and processing stages, with ingestion and organization steps often contributing comparable or even higher computational costs than actual processing steps. Since complex workflows consist of a variety of phases that view and use data differently, being able to construct efficient, scalable, distributed data structures (arrays, vectors, sets, maps, and multi-maps) is essential and requires custom methods to extend and shrink containers, analyze and position data, and, maintain globallyconsistent meta-data. In this paper, we propose a novel datastructure access paradigm based on the concept of Accessors. At a high level, accessors are customizable callable objects that can modify the behavior of insert, read, update, and delete operations for distributed containers while preserving atomicity guarantees. Accessors provide a very clean and natural way to implement a variety of programming patterns, e.g., conditional insertion/deletion and cascading computations, which would be otherwise hard (or even impossible) to express in parallel and distributed settings without using locks. We demonstrate the practicality and usefulness of our approach with two representative use cases and study the performance of these applications on a distributed High-Performance Computing system. Our analysis highlights that our proposed abstraction allows for an effective overlapping and concurrent execution of different workflow steps (e.g., data ingestion and analysis), which in a conventional analytics pipeline would execute sequentially, contributing cumulatively to the overall latency.