Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distribution adaptation”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Numerical Modeling & Optimization of the iProTech Pitching Inertial Pump (PIP) Wave Energy Converter (WEC) (CRADA Final Report)

This project represents a continuation of the collaboration between iProTech and NLR to simulate, optimize and design the iProTech Pitching Inertial Pump (PIP) device. The objectives of this TEAMER project are twofold: 1. Refining the physical characteristics of the existing iProTech PIP WEC-Sim model to enhance the model’s fidelity and include controllable components. Key model enhancements target the inclusion of Coulomb friction, the introduction of a controllable bypass valve, and the replacement of traditional check valves with advanced motorized ones. 2. Exploring traditional and advanced control algorithms. From traditional methods like latching control to cutting-edge reinforcement learning (RL) algorithms, the goal is to ensure the PIP device's adaptability and optimal performance across a range of ocean conditions. NLR is tasked with augmenting the WEC-Sim model and implementing the control algorithms, culminating in performance comparison analyses. iProTech will update their existing 3D models, advise on model improvements, and determine crucial system metrics. WEC-Sim, developed in MATLAB/SIMULINK with Simscape Multibody, is the main piece of software that will be used in this project. Coupled with the MATLAB RL Toolbox, it offers a robust platform for in-depth simulation and optimization of the iProTech PIP device. Building on previous work to explore the PIP design space and optimize its geometry, mass distribution, center of gravity and other key parameters, this project aims to refine iProTech’s existing numerical models and develop effective control algorithms that can seamlessly integrate into their future hardware testing campaigns.

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

An Intercomparison of Large‐Eddy Simulations of a Convection Cloud Chamber Using Haze‐Capable Bin and Lagrangian Cloud Microphysics Schemes

Abstract Recent in situ observations show that haze particles exist in a convection cloud chamber. The microphysics schemes previously used for large‐eddy simulations of the cloud chamber could not fully resolve haze particles and the associated processes, including their activation and deactivation. Specifically, cloud droplet activation was modeled based on Twomey‐type parameterizations, wherein cloud droplets were formed when a critical supersaturation for the available cloud condensation nuclei (CCN) was exceeded and haze particles were not explicitly resolved. Here, we develop and adapt haze‐capable bin and Lagrangian microphysics schemes to properly resolve the activation and deactivation processes. Results are compared with the Twomey‐type CCN‐based bin microphysics scheme in which haze particles are not fully resolved. We find that results from the haze‐capable bin microphysics scheme agree well with those from the Lagrangian microphysics scheme. However, both schemes significantly differ from those from a CCN‐based bin microphysics scheme unless CCN recycling is considered. Haze particles from the recycling of deactivated cloud droplets can strongly enhance cloud droplet number concentration due to a positive feedback in haze‐cloud interactions in the cloud chamber. Haze particle size distributions are more realistic when considering solute and curvature effects that enable representing the complete physics of the activation process. Our study suggests that haze particles and their interactions with cloud droplets may have a strong impact on cloud properties when supersaturation fluctuations are comparable to mean supersaturation, as is the case in the cloud chamber and likely is the case in the atmosphere, especially in polluted conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Fluid-Kinetic Coupling: Advanced Discretizations for Simulations on Emerging Heterogeneous Architectures (LDRD FY20-0643)

Plasma physics simulations are vital for a host of Sandia mission concerns, for fundamental science, and for clean energy in the form of fusion power. Sandia's most mature plasma physics simulation capabilities come in the form of particle-in-cell (PIC) models and magnetohydrodynamics (MHD) models. MHD models for a plasma work well in denser plasma regimes when there is enough material that the plasma approximates a fluid. PIC models, on the other hand, work well in lower-density regimes, in which there is not too much to simulate; error in PIC scales as the square root of the number of particles, making high-accuracy simulations expensive. Real-world applications, however, almost always involve a transition region between the high-density regimes where MHD is appropriate, and the low-density regimes for PIC. In such a transition region, a direct discretization of Vlasov is appropriate. Such discretizations come with their own computational costs, however; the phase-space mesh for Vlasov can involve up to six dimensions (seven if time is included), and to apply appropriate homogeneous boundary conditions in velocity space requires meshing a substantial padding region to ensure that the distribution remains sufficiently close to zero at the velocity boundaries. Moreover, for collisional plasmas, the right-hand side of the Vlasov equation is a collision operator, which is non-local in velocity space, and which may dominate the cost of the Vlasov solver. The present LDRD project endeavors to develop modern, foundational tools for the development of continuum-kinetic Vlasov solvers, using the discontinuous Petrov-Galerkin (DPG) methodology, for discretization of Vlasov, and machine-learning (ML) models to enable efficient evaluation of collision operators. DPG affords several key advantages. First, it has a built-in, robust error indicator, allowing us to adapt the mesh in a very natural way, enabling a coarse velocity-space mesh near the homogeneous boundaries, and a fine mesh where the solution has fine features. Second, it is an inherently high-order, high-intensity method, requiring extra local computations to determine so-called optimal test functions, which makes it particularly suited to modern hardware in which floating-point throughput is increasing at a faster rate than memory bandwidth. Finally, DPG is a residual-minimizing method, which enables high-accuracy computation: in typical cases, the method delivers something very close to the $L^2$ projection of the exact solution. Meanwhile, the ML-based collision model we adopt affords a cost structure that scales as the square root of a standard direct evaluation. Moreover, we design our model to conserve mass, momentum, and energy by construction, and our approach to training is highly flexible, in that it can incorporate not only synthetic data from direct-simulation Monte Carlo (DSMC) codes, but also experimental data. We have developed two DPG formulations for Vlasov-Poisson: a time-marching, backward-Euler discretization and a space-time discretization. We have conducted a number of numerical experiments to verify the approach in a 1D1V setting. In this report, we detail these formulations and experiments. We also summarize some new theoretical results developed as part of this project (published as papers previously): some new analysis of DPG for the convection-reaction problem (of which the Vlasov equation is an instance), a new exponential integrator for DPG, and some numerical exploration of various DPG-based time-marching approaches to the heat equation. As part of this work, we have contributed extensively to the Camellia open-source library; we also describe the new capabilities and their usage. We have also developed a well-documented methodology for single-species collision operators, which we applied to argon and demonstrated with numerical experiments. We summarize those results here, as well as describing at a high level a design extending the methodology to multi-species operators. We have released a new open-source library, MLC, under a BSD license; we include a summary of its capabilities as well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Multiscale Electricity Modeling for Evaluating Carbon Capture and Sequestration Technologies (MEME-CCS)

This effort employs and adapts a rigorous multiscale electricity modeling platform at the National Renewable Energy Laboratory (NREL) to evaluate carbon capture and sequestration (CCS) and negative emissions technologies (NET) from the ARPA-E FLECCS program. NREL's modeling platform includes the Regional Energy Deployment System (ReEDS) electric sector capacity expansion model, which projects future electricity generation mixes at sub-state-level resolution that are downscaled to the unit-level to enable hourly, zonal or nodal electricity production cost modeling in the PLEXOS model. The ReEDS-PLEXOS modeling suite is well-established for examining electric sector futures with high renewable energy penetrations, energy storage, electrification, and distributed generation. The key advancement proposed herein utilizes collaboration with CCS experts at the University of Wyoming (U.WY) and the FLECCS technology development teams to create innovative methods for representing CCS and NET in the ReEDS and PLEXOS models. Expanded technology options and new operational parameterizations are integrated into these models to allow an unprecedented combination of scope and resolution for exploring the future of CCS and NET. The resulting capabilities permit wide-ranging scenario analysis to assess CCS and NET deployment under alternative scenarios of CO2 prices and competitiveness of flexible CCS and NET technologies. We demonstrate sample deployment and operational outcomes to show how these models are being used to assess the future potential for FLECCS technologies and their impacts on the U.S. electricity system. These products will help an emerging CCS/NET industry in the United States by providing economics-driven guidance to technology developers while informing policy and investment decisions in the public and private sectors.

capacity expansion↗
DOE OSTI2024

In-Process Monitoring and Structural Health Monitoring of Large-Scale Additive Manufacturing Using Acoustic Emission Technique

ORNL collaborated with MISTRAS Group, Inc. to investigate acoustic emission (AE) as a structural health monitoring (SHM) method for large-scale additive manufacturing (AM). Large-scale AM is being adapted as method of producing large structures in a short lead time and cost-effective way. With the growing advancement in AM techniques and application, machine monitoring and part qualification is highly needed. There has been leading research focused on the manufacturing, feedstock material but minimum research on the SHM, defect detection, and nondestructive evaluation (NDE) for AM. Scanning large structure using conventional nondestructive testing (NDT) techniques, such as ultrasound or X-ray, and searching for potential defects can be very time consuming, challenging and cost prohibitive. AE is a passive technique that can be used to monitor and locate defect progression in large structure by distributing group of sensors around the part. This project utilized AE technique and system manufactured/designed by MISTRAS Group to monitor large-scale AM equipment (i.e. Big Area Additive Manufacturing (BAAM) system located at the Oak Ridge National Laboratory – Manufacturing Demonstration Facility (ORNL-MDF) and the printed parts it produces. The AE system provided valuable insight on defect development/progression during and post-printing process.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Enhanced Feedstock Characterization and Modeling to Facilitate Optimal Preprocessing and Deconstruction of Corn Stover (Final Report)

This project addresses the challenge of processing corn stover by fractionating this biomass feedstock to both streamline processing and generate new potential co-products. Additionally, the project developed new field-deployable analytical tools that can be coupled with empirical models that were used to predict feedstock properties and processing performance. The overall scope of this project was: (1) identify conditions for optimal corn stover fractionation using a two- stage physical fractionation, (2) assess how physical fractionation impacts properties, partitioning of biomass, and response to processing, (3) further adapt, develop, and validate several advanced characterization tools for assessing biomass properties that can be linked to processing behavior, and (4) develop and validate predictive models based on measurements that can be performed “in the field” or “at the biorefinery gate” to predict feedstock processing behavior (preprocessing and deconstruction). The first objective employed pre-separation processing (size reduction) which was next subjected to enhanced separations to yield fractions enriched or depleted in select compositional components or properties. For the second objective, fractions were screened for their response to post-separation processing (pretreatment and enzymatic hydrolysis). Detailed characterization profiles were developed and dynamic image analysis to assess distribution of particle size and morphology. For the final objective, we utilized these tools to develop empirical models to assess the relative abundance of tissue type in order to assess fractionation efficacy and to predict fraction performance during pretreatment and enzymatic hydrolysis.

09 BIOMASS FUELS↗
DOE OSTI2026

HydroDCM: Hydrological Domain-Conditioned Modulation for Cross-Reservoir Inflow Prediction

Deep learning models have shown promise in reservoir inflow prediction, yet their performance often deteriorates when applied to different reservoirs due to distributional differences, referred to as the domain shift problem. Domain generalization (DG) solutions aim to address this issue by extracting domain-invariant representations that mitigate errors in unseen domains. However, in hydrological settings, each reservoir exhibits unique inflow patterns, while some metadata beyond observations like spatial information exerts indirect but significant influence. This mismatch limits the applicability of conventional DG techniques to many-domain hydrological systems. To overcome these challenges, we propose HydroDCM, a scalable DG framework for cross-reservoir inflow forecasting. Spatial metadata of reservoirs is used to construct pseudo-domain labels that guide adversarial learning of invariant temporal features. During inference, HydroDCM adapts these features through light-weight conditioning layers informed by the target reservoir’s metadata, reconciling DG’s invariance with location-specific adaptation. Experiment results on 30 real-world reservoirs in the Upper Colorado River Basin demonstrate that our method substantially outperforms state-of-the-art DG baselines under many-domain conditions and remains computationally efficient.

Hu, Pengfei [ORNL] (ORCID:0009000367130950)↗
DOE OSTI2026

Exploring life’s hidden majority: microbial dark matter symposium highlights

The Microbial Dark Matter Symposium held on August 28–29, 2025, in Laguna Beach, Orange County, CA, convened a multidisciplinary group of scientists to address the vast unknowns in microbial life—from uncultured taxa and uncharacterized proteins to elusive viruses and spacefaring microbes. Set against a scenic coastal backdrop, the symposium highlighted advances in single-cell genomics, proximity ligation sequencing, and artificial intelligence-ready bioinformatics, while also probing the limits of microbial persistence, metabolism, and ecological distribution. Sessions explored microbial dark matter from multiple dimensions: cultivability, where new strategies are enabling recovery of elusive microbes; functional ambiguity, where metagenomic dark zones are illuminated by computational annotation; and genomic representation, where single-cell methods bridge gaps left by shotgun community sequencing. Researchers shared breakthroughs in identifying atmospheric microbiomes, “dark oxygen” production in groundwater ecosystems, and microbial survival on the International Space Station. The symposium emphasized integration of methods, disciplines, and ecosystems, advancing a collective push to illuminate the microbial dark matter on Earth and beyond. By highlighting emerging tools, pressing questions, and cross-domain insights, the symposium underscored the need for collaborative, open, and adaptive approaches to study the microbial unknown. The meeting marks a pivotal moment in microbiology, where cultivating knowledge of the uncultivated promises transformative understanding of life, everywhere.

Podar, Mircea [ORNL] (ORCID:0000000327760205)↗
DOE OSTI2023

Envisioning urban environments resilient to vector-borne diseases: a protocol to study dengue in Vietnam

Transmitted primarily by Aedes aegypti (Ae. aegypti) and Aedes albopictus (Ae. albopictus), arboviral diseases pose a major global public health threat. Dengue, chikungunya, and zika are increasingly prevalent in Southeast Asia. Among other arboviruses, dengue and zika are becoming more common in Central and South America. Given human encroachment into previously uninhabited, often deforested areas, to provide new housing in regions of population expansion, conceptualizing built urban environments in a novel way is urgently needed to safeguard against the growing climate change-driven threat of vector-borne diseases. By understanding the spread from a One Health perspective, enhanced control and prevention can be achieved. This is particularly important considering that climate change is likely to significantly impact the persistence of ponded water where mosquitoes breed due to increasing temperature and shifting rainfall patterns with regard to magnitude, duration, frequency, and season. Models can incorporate aquatic mosquito stages and adult spatial dynamics when habitats are heterogeneously available, thereby including dispersal and susceptible-exposed-infected-recovered (SEIR) epidemiology. Coupled with human population distribution (density, locations), atmospheric conditions (air temperature, precipitation), and hydrological conditions (soil moisture distribution, ponding persistence in topographic depressions), modeling has improved predictive ability for infection rates. However, it has not informed interventional approaches from an urban environment perspective which considers the role of ponds/lakes that support green spaces, the density of population that enables rapid spread of disease, and varying micro-habitats for various mosquito stages under climate change. Here, for an example of dengue in Vietnam, a preventive and predictive approach to design resilient urban environments is proposed, which uses data from rapidly expanding metropolitan communities to learn continually. This protocol deploys computational approaches including simulation and machine learning/artificial intelligence, underpinned by surveillance and medical data for validation and adaptive learning. Its application may best inform urban planning in low-middle income countries in tropical zones where arboviral pathogens are prevalent.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Self-oscillating synchronematic colloids

Self-oscillators that sustain periodic dynamics under constant input are ubiquitous in natural and engineered systems, where their interactions enable spatiotemporal coordination among many individual units. New forms of organization can emerge when these self-oscillating units are free to move and rotate, coupling their spatial arrangement and alignment with their oscillation frequencies and phases. Here, we report experiments and simulations on populations of Quincke colloids that behave as self-oscillating units with position, orientation, frequency, and phase. Depending on the initial distribution, these active oscillators spontaneously organize into distinct collective states characterized by temporal synchronization and directional alignment, which we term synchronematic order. In fluid-like clusters, this order is short-ranged and decays over a length scale set by the competition between hydrodynamic interactions and athermal noise. In crystalline clusters, these interactions drive flobal synchronization and circular alignment-synchronematic crystals-whose collective frequency increases with cluster size due to non-reciprocal interactions. Our results establish self-oscillating colloids as a model system for active oscillatory matter and reveal fundamental principles by which synchronization, alignment, and structure co-emerge, offering new pathways for designing adaptive, frequency-tunable materials.

colloids↗
DOE OSTI2021

Optimal Bayesian supervised domain adaptation for RNA sequencing data

Abstract Motivation When learning to subtype complex disease based on next-generation sequencing data, the amount of available data is often limited. Recent works have tried to leverage data from other domains to design better predictors in the target domain of interest with varying degrees of success. But they are either limited to the cases requiring the outcome label correspondence across domains or cannot leverage the label information at all. Moreover, the existing methods cannot usually benefit from other information available a priori such as gene interaction networks. Results In this article, we develop a generative optimal Bayesian supervised domain adaptation (OBSDA) model that can integrate RNA sequencing (RNA-Seq) data from different domains along with their labels for improving prediction accuracy in the target domain. Our model can be applied in cases where different domains share the same labels or have different ones. OBSDA is based on a hierarchical Bayesian negative binomial model with parameter factorization, for which the optimal predictor can be derived by marginalization of likelihood over the posterior of the parameters. We first provide an efficient Gibbs sampler for parameter inference in OBSDA. Then, we leverage the gene-gene network prior information and construct an informed and flexible variational family to infer the posterior distributions of model parameters. Comprehensive experiments on real-world RNA-Seq data demonstrate the superior performance of OBSDA, in terms of accuracy in identifying cancer subtypes by utilizing data from different domains. Moreover, we show that by taking advantage of the prior network information we can further improve the performance. Availability and implementation The source code for implementations of OBSDA and SI-OBSDA are available at the following link. https://github.com/SHBLK/BSDA. Supplementary information Supplementary data are available at Bioinformatics online.

Biochemistry & Molecular Biology↗
DOE OSTI2021

Butterfly Factorization Via Randomized Matrix-Vector Multiplications

This paper presents an adaptive randomized algorithm for computing the butterfly factorization of an m × n matrix with m ≈ n provided that both the matrix and its transpose can be rapidly applied to arbitrary vectors. The resulting factorization is composed of O(log n) sparse factors, each containing O(n) nonzero entries. The factorization can be attained using O(n 3/2 log n) computation and O(n log n) memory resources. Furthermore, the proposed algorithm can be implemented in parallel and can apply to matrices with strong or weak admissibility conditions arising from surface integral equation solvers as well as multi-frontal-based finite-difference, finite-element, or finite-volume solvers. A distributed-memory parallel implementation of the algorithm demonstrates excellent scaling behavior.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Ultrasmall amorphous zirconia nanoparticles catalyse polyolefin hydrogenolysis

Carbon–carbon bond cleavage reactions, adapted to deconstruct aliphatic hydrocarbon polymers and recover the intrinsic energy and carbon value in plastic waste, have typically been catalysed by metal nanoparticles or air-sensitive organometallics. Metal oxides that serve as supports for these catalysts are typically considered to be inert. Here we show that Earth-abundant, non-reducible zirconia catalyses the hydrogenolysis of polyolefins with activity rivalling that of precious metal nanoparticles. To harness this unusual reactivity, our catalytic architecture localizes ultrasmall amorphous zirconia nanoparticles between two fused platelets of mesoporous silica. Macromolecules translocate from bulk through radial mesopores to the highly active zirconia particles, where the chains undergo selective hydrogenolytic cleavage into a narrow, C18-centred distribution. Calculations indicated that C–H bond heterolysis across a Zr–O bond of a Zr(O) 2 adatom model for unsaturated surface sites gives a zirconium hydrocarbyl, which cleaves a C–C bond via β-alkyl elimination.

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

Knowledge Oriented Graph Unified Transformer (KOGUT) v0.1

KOGUT — Knowledge Oriented Graph Unified Transformer KOGUT implements the Relational Graph Transformer (RelGT) architecture for knowledge graph link prediction in biological domains, with a primary focus on microbial growth media prediction. While the original RelGT (arXiv:2505.10960) targets relational tables, time series, and multi-table databases, KOGUT adapts this architecture for heterogeneous biological knowledge graphs, providing first-in-class AI predictive models for microbial cultivation. Key Adaptations Beyond Original RelGT: - Knowledge Graph Focus: Applied to biological KGs with semantic node types (taxa, chemicals, media, phenotypes, environments) versus generic relational database tables, trained on the KG-Microbe knowledge graph (1.3M entities, 2.9M edges, 24 relation types). - Multimodal Node Encoding: Integrates node labels, categories, descriptions, and synonyms from KG metadata through learned embedding layers—adapting relational column features to graph node attributes with textual semantics. - Extended K-Hop Subgraph Strategy: Optimized neighborhood sampling (3-hop default, configurable up to 200 nodes) tuned for sparse biological networks, building on the original local-global attention framework with biological relation preservation. - Biolink Predicate Preservation: Type-specific transformations for 24 biological edge semantics (occurs_in, consumes, produces, has_phenotype, subclass_of) beyond standard relational foreign keys, enabling multi-relation link prediction. - Inductive Learning Support: Enables zero-shot predictions for novel taxa through feature-based embeddings (temperature, oxygen requirements, gram stain, cell shape), extending the original transductive relational benchmark scope to uncultured microorganisms. CheapSOTA Performance Optimizations (This Distribution): - VQ-EMA Centroid Attention: Vector quantization with exponential moving average for improved global context modeling (+5-10% MRR improvement). - HDF5 Precomputed Data Loading: One-time preprocessing of k-hop subgraphs to eliminate redundant graph traversals (2-5× training speedup). - Distributed Data Parallel Training: Multi-GPU support for scaling to larger knowledge graphs (tested on 4× NVIDIA A100 GPUs at NERSC Perlmutter). - Mixed Precision Training: Automatic mixed precision (AMP) for memory efficiency and faster training. Advantages Over Standard Knowledge Graph Embedding Models: Combines RelGT's proven multi-element tokenization (features, type, hop, structure) with graph-native biological representations, enabling interpretable link prediction across heterogeneous entities that standard embedding models (TransE, RotatE, ComplEx) and table-based transformers cannot directly model. Achieves near-perfect performance on microbial growth media prediction (MRR: 0.9966, Precision@1: 0.9932, Hit@10: 1.0000) while maintaining explainability through attention-based reasoning over biological pathways. Training Data: - KG-Microbe merged knowledge graph: 1,379,337 nodes, 2,960,472 edges - 24 biological relation types including taxonomic hierarchies, metabolic interactions, phenotype associations, and environmental relationships - Primary prediction task: Growth media suitability for microbial taxa (biolink:occurs_in, 50K edges) - Multi-relation capability: Predicts links for any of the 24 relation types, including chemical consumption/production, phenotype associations, and taxonomic classification Citation: Original RelGT Architecture: Dwivedi et al., "Relational Graph Transformer", arXiv:2505.10960, 2025 KOGUT Implementation: Knowledge Oriented Graph Unified Transformer for Microbial Growth Media Prediction Developed at Lawrence Berkeley National Laboratory (LBNL) Trained on NERSC Perlmutter supercomputer

Joachimiak, Marcin [Lawrence Berkeley National Lab↗
DOE OSTI2025

Enhancing Security and Resiliency in Operational Technology Environments Through Network Slicing and Federated Learning

The growing convergence of Information Technology (IT) and Operational Technology (OT) within Industry 4.0 environments has introduced new demands on industrial network infrastructure. As cyber-physical systems become increasingly interconnected, ensuring the secure, timely, and efficient exchange of critical data is essential. This thesis explores how network slicing, a method of creating isolated virtual network segments, can be applied within OT environments to address challenges such as latency, security, and resource allocation. The first research question addressed in this thesis is: How can OT networks take advantage of NFV and SDN technology to become cyber resilient? This study examines the operational, security, and architectural implications of introducing network slicing into traditionally static OT infrastructures such as Industrial Control Systems (ICS) and SCADA. Through simulated deployments and case studies, the research demonstrates how slicing enables better isolation between critical and non-critical services, thereby improving response time, throughput, and security in sensitive environments. The second question considers: How to dynamically implement network slicing and take advantage of network resources towards integrating decentralized machine learning? In response, this thesis proposes a framework that combines Software-Defined Networking (SDN), Network Function Virtualization (NFV), and Federated Learning (FL) to enable real-time analytics while maintaining data locality. The proposed approach reduces the burden on centralized infrastructure and minimizes privacy risks by supporting on-site training of models across distributed OT nodes, coordinated through dynamically allocated network slices. The third focus explores: How slicing helps to increase the resiliency of OT networks through the orchestration of a dynamic DMZ? To answer this, the thesis presents a method for creating and managing Dynamic Demilitarized Zones (DMZs) using network slicing. This enables flexible and automated isolation of sensitive subsystems during threat scenarios or high-risk operations. Coupled with intelligent orchestration and containerized security services, the dynamic DMZ significantly enhances the system's ability to respond to cyber incidents without halting production. Ultimately, this thesis contributes a comprehensive architecture that blends network slicing with machine learning, secure segmentation, and automation, paving the way for resilient, adaptive, and intelligent OT environments. Performance evaluations across multiple scenarios show improvements in system reliability, threat response time, model accuracy, and resource utilization, providing a strong foundation for future industrial automation systems.

Rodiles Delgado, Brian G↗
DOE OSTI2020

P SI4 1.4: Open-source software for high-throughput quantum chemistry

PSI4 is a free and open-source ab initio electronic structure program providing implementations of Hartree–Fock, density functional theory, many-body perturbation theory, configuration interaction, density cumulant theory, symmetry-adapted perturbation theory, and coupled-cluster theory. Most of the methods are quite efficient, thanks to density fitting and multi-core parallelism. The program is a hybrid of C++ and Python, and calculations may be run with very simple text files or using the Python API, facilitating post-processing and complex workflows; method developers also have access to most of PSI4 's core functionalities via Python. Job specification may be passed using The Molecular Sciences Software Institute (MolSSI) QCSCHEMA data format, facilitating interoperability. We announce a rewrite of our top-level computation driver, and concomitant adoption of the MolSSI QCARCHIVE INFRASTRUCTURE project, makes the latest version of PSI4 well suited to distributed computation of large numbers of independent tasks. The project has fostered the development of independent software components that may be reused in other quantum chemistry programs.

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

Simulating single-particle dynamics in magnetized plasmas: The RMF code

The RMF (Rotating Magnetic Field) code is designed to calculate the motion of a charged particle in a given electromagnetic field. It integrates Hamilton’s equations in cylindrical coordinates using an adaptive predictor-corrector double-precision variable-coefficient ordinary differential equation solver for speed and accuracy. RMF has multiple capabilities for the field. Particle motion is initialized by specifying the position and velocity vectors. Here, the six-dimensional state vector and derived quantities are saved as functions of time. A post-processing graphics code, XDRAW, is used on the stored output to plot up to 12 windows of any two quantities using different colors to denote successive time intervals. Multiple cases of RMF may be run in parallel and perform data mining on the results. Recent features are a synthetic diagnostic for simulating the observations of charge-exchange-neutral energy distributions and RF grids to explore a Fermi acceleration parallel to static magnetic fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Towards fully predictive gyrokinetic full- f simulations: validation and triangularity studies in TCV

Designing economical magnetic confinement fusion power plants motivates computational tools that can estimate plasma behavior from engineering parameters without direct reliance on experimental measurement of the plasma profiles. In this work, we present full-f global long-wavelength gyrokinetic simulations of edge and scrape-off layer turbulence in tokamaks that use only magnetic geometry, heating power, and particle inventory as inputs. Unlike many modeling approaches that employ free parameters fitted to experimental data, raising uncertainties when extrapolating to reactor scales. This approach directly simulates turbulence and resulting profiles through gyrokinetics without such empirical adjustments. This is achieved via an adaptive sourcing algorithm in Gkeyll that strictly controls energy injection and emulates particle sourcing due to neutral recycling. We show that the simulated kinetic profiles compare reasonably well with Thomson scattering and Langmuir probe data for Tokamak á Configuration Variable (TCV) discharge #65125, and that the simulations reproduce characteristic features such as blob transport and self-organized electric fields. Applying the same framework to study triangularity effects suggests mechanisms contributing to the improved confinement reported for negative triangularity (NT). Simulations of TCV discharges #65125 and #65130 indicate that NT increases the E x B flow shear (by about 20% in these cases), which correlates with reduced turbulent losses and a modest change in the distribution of power exhaust to the vessel wall. While the physical models contain approximations that can be refined in future work, the predictive capability demonstrated here, evolving multiple profile relaxation times with kinetic electron and ion models in hundreds of GPU hours, indicates the feasibility of using Gkeyll to support design studies of fusion devices.

Hoffmann, Antoine Cyril David [Princeton Plasma Ph↗