Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Moving interface”

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.

71 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The Elastic Analysis Facility's (EAF's) Contribution to the Future of Analysis at Multi-Experiment Institutions and Future Colliders

The Elastic Analysis Facility (EAF) hosted at Fermi National Accelerator Laboratory (Fermilab) is a platform being developed with the goal of providing a fast and efficient facility for physics analysis. As high-energy physics moves towards collecting larger datasets, such as those from the High-Luminosity LHC, the EAF strives to provide a powerful and adaptable framework for future colliders and multi-experiment institutions. Currently, the EAF supports several experiments including CMS, NOvA, and DUNE as well as serving accelerator physicists and beam line operations through integrated software and secure connections to Fermilab's computing resources. In addition, the EAF was designed with a user-friendly interface, intended to be more intuitive for emerging generations of physicists, that is still accessible for established styles of analysis. The EAF can also achieve better analysis efficiency due to the modernization of software and tools that can better utilize Fermilab's computing power. Furthermore, its design incorporates industry standards whenever possible, enhancing its sustainability and making it a possible template for other national or international laboratories and research facilities. Overall, the EAF is a forward-looking solution that will meet the evolving needs of particle physics, ensuring readiness for future colliders and multi-experiment research institutions.

Chavez, Elise [Wisconsin U., Madison]↗
DOE OSTI2025

Commercialization of the Transportation-Security, Tracking, and Reporting System (T-STAR)

The Transportation-Security, Tracking, and Reporting System (T-STAR) was developed by the National Nuclear Security Administration, NA-21, Office of Radiological Security (ORS) to provide a transportation security system for detection and tracking during transport of Category 1 and Category 2 radiological material. Few off-the-shelf systems for conveyance tracking offer detection of a cargo compartment breach or a removal of the cargo. Systems that do offer this capability often require permanent installation through modifying of the conveyance itself. This is not sustainable in many countries where ORS is building use, storage, and transport security capacity. The development of T-STAR has moved from fielding robust prototypes deployed in countries ranging from North America, Latin America and Central Asia to a commercially produced product that can now be deployed to provide enhanced security during transit. Each prototype deployment resulted in important lessons learned, which informed the requirements for the final commercial product. T-STAR uses both cellular and Iridium satellite modems to provide redundant communications to provide the configuration, status, and alerts to a server monitoring the shipment, which is accessible using a multilanguage browser-based user interface. A wireless security system employing using Z-wave sensors for intrusion detection located in the conveyance provide low cost but effective solution for a wide range of conveyance types. Additional capabilities include the ability to monitor a vehicles’ CANBUS (Controller Area Network) system, an ethernet port for high throughput sensor information such as video cameras, and the ability to power and use advanced external sensor payloads. These features make the T-STAR a capable and expandable security gateway that can be deployed on a variety of conveyances from box trucks to open trailers. The ability to provide tracking, monitoring, and detection provide a key component in overall best practices designed to protect shipments of radioactive material.

Schultze, Michael [ORNL] (ORCID:0000000283205671)↗
DOE OSTI2023

Activity-induced migration of viscous droplets on a solid substrate

Active matter exploits motion to induce changes in shape and conformation via external input. Here, in this paper, we establish theoretically that viscous liquid droplets containing magnetic nanoparticles with frozen-in magnetic moments, sitting on a solid substrate and surrounded by an ambient gas phase, can deform and migrate under the influence of a magnetic torque. The effect arises because the collective rotation of the magnetic nanoparticles at the liquid–gas interface tilts the droplet away from a symmetric configuration, breaks the reflection symmetry with respect to the centre axis, and leads to a left–right asymmetry of the contact angles. A sufficiently strong magnetic torque leads the contact angles to overcome hysteresis effects leading the droplet to migrate. We develop a general framework to explain how symmetry-breaking affects droplet migration. Thus previous results of droplet spreading and migration can be recovered as special cases. Such droplets can be employed as agents in active surfaces and can move against gravity, chemical and thermal gradients, providing a mechanism that could be utilized by both industry and medicine.

Aggarwal, A. [Northwestern Univ., Evanston, IL (Un↗
DOE OSTI2025

Generalizing the compressible pairwise interaction extended point-particle model

Ejecta physics plays an important role in material interfaces that are impacted by a strong shock wave. When a shock impacts a rough surface of solid material and melts it, the Richtmyer–Meshkov instability grows perturbations on the surface, which can eject particles. After release, the ejecta travel through the post-shock compressible flow. To accurately simulate a large number of ejecta particles, an Euler–Lagrange approach is preferred, which requires modeling the subgrid-scale physics involved with fluid–particle interactions. We generalize the previous work from Hsiao et al. (2023) to consider systems of moving particles subject to any loading shock. The following improvements were made: (1) Particles are allowed to move relative to each other (2) Non-planar shocks are accounted for along with allowing for variable shock speeds. As a result, the generalized algorithm was tested with particle-resolved simulations for canonical test cases. The results of these tests are discussed and analyzed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Building Artificial Layered Solids from the Bottom-up: Materials by Design to Enable New Energy Technologies

This final technical report summarizes the key accomplishments on this DOE Early Career Program award received by PI Guihua Yu at the University of Texas at Austin. The main goal of this ECP award was to develop synthetic methodologies, self-assembly approaches towards structurally controlled nanosheets-like solids from the bottom up, and to understand and characterize their chemical/physical characteristics for the design of future-generation energy devices with novel functions and properties, that will have direct implications for energy science and technologies such as energy storage and conversion, and catalysis. The project accomplished these goals by completing the following objectives: • Rational design, synthesis, and self-assembly of structurally-controlled molecular ‘nanosheets’ materials. Moving beyond conventional van der Waals two-dimensional crystals, this project enables new synthesis and processing strategies to develop ultrathin nanosheets materials, from inorganic transition-metal oxides-based, to hybrid organic-inorganic nanomaterials, with structural factors such as facet, thickness and pore structure being well controlled during synthesis or assembly processes. • Fundamental understanding and electrochemical characterization of these assembled molecular ‘nanosheets’ materials via charge/mass transport studies through electrochemical intercalation of different metal ions for next-generation energy storage, as well as electrocatalytic studies using this new material platform for understanding catalytic reaction kinetics and the roles of surface functionalization and interface interactions owing to atomic thickness, nanoscale porosity, and other structural factors.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Enhancing the Optically Detected Magnetic Resonance Signal of Organic Molecular Qubits

In quantum information science and sensing, electron spins are often purified into a specific polarization through an optical-spin interface, a process known as optically detected magnetic resonance (ODMR). Diamond-NV centers and transition metals are both excellent platforms for these so-called color centers, while metal-free molecular analogues are also gaining popularity for their extended polarization lifetimes, milder environmental impacts, and reduced costs. In our earlier attempt at designing such organic high-spin π-diradicals, we proposed to spin-polarize by shelving triplet M S = ±1 populations as singlets. This was recently verified by experiments albeit with low ODMR contrasts of <1% at temperatures above 5 K. In this work, we propose to improve the ODMR signal by moving singlet populations back into the triplet M S = 0 sublevel, designing a true carbon-based molecular analogue to the NV center. Our proposal is based upon transition-orbital and group-theoretical analyses of beyond-nearest-neighbor spin–orbit couplings, which are further confirmed by ab initio calculations of a realistic trityl-based radical dimer. Microkinetic analyses point toward high ODMR contrasts of around 30% under experimentally feasible conditions, a stark improvement from previous works. Finally, in our quest toward ground-state optically addressable molecular spin qubits, we exemplify how our symmetry-based design avoids Zeeman-induced singlet–triplet mixings, setting the scene for realizing electron spin qubit gates.

Group theory↗
DOE OSTI2025

Electrical and Thermal Load Impacts of Three District Heating and Cooling Designs for an Existing Community in Washington, DC

District energy systems that provide building heating and cooling are a promising option to provide low-cost, energy efficient heating and cooling solutions for communities. As some buildings move to electrified heating designs, there may be significant increases in electric grid demand, particularly in the winter. District energy systems may be able to help reduce these high electrical demands on the electric grid, while also providing decreased overall energy consumption with an increased flexibility in electrical energy usage. To evaluate and compare district options for district-based heating and cooling, an existing neighborhood of Washington, DC consisting of 35 existing buildings was selected as a case study. This study models and compares a fourth generation (4G) district heating and cooling systems that provide hot and chilled water from a central plant directly to each prosumer and a fifth generation (5G) district heating and cooling systems that provide near-ambient water via geothermal boreholes to interface with an energy transfer station with a heat pump at each prosumer. These district systems are compared against a baseline of the non-connected buildings with their self-contained and current HVAC systems. The neighborhood was selected as a good candidate for a district system for its diversity of loads (mix of different commercial and property types) and the availability of building characteristics. URBANopt District Energy Systems (DES) was used to model the buildings and create initial Modelica models for the district systems. Finally, the DES models were tuned and simulated in Dymola, and the outputs were post-processed for verification and to calculate the electrical and thermal grid impact metrics. This is the first documentation of this workflow and its complete analysis. The three designed systems are compared via their grid metrics including: total energy consumption, electrical demand peak loads, daily peak-to-valley ratios, and system ramping. A detailed discussion is provided about how each system impacts and interacts with the electric grid.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Evanescent and inertial-like waves in rigidly rotating odd viscous liquids

Three-dimensional non-rotating odd viscous liquids give rise to Taylor columns and support axisymmetric inertial-like waves ( J. Fluid Mech ., vol. 973, 2023, A30). When an odd viscous liquid is subjected to rigid-body rotation however, there arise in addition a plethora of other phenomena that need to be clarified. In this paper, we show that three-dimensional incompressible or two-dimensional compressible odd viscous liquids, rotating rigidly with angular velocity 𝛺, give rise to both oscillatory and evanescent inertial-like waves or a combination thereof (which we call of mixed type) that can be non-axisymmetric. By evanescent, we mean that along the radial direction, typically when moving away from a solid boundary, the velocity field decreases exponentially. These waves precess in a prograde or retrograde manner with respect to the rotating frame. The oscillatory and evanescent waves resemble respectively the body and wall-modes observed in (non-odd) rotating Rayleigh–Bénard convection ( J. Fluid Mech ., vol. 248, 1993, pp. 583–604). We show that the three types of waves (wall, body or mixed) can be classified with respect to pairs of planar wavenumbers 𝜅 which are complex, real or a combination, respectively. Experimentally, by observing the precession rate of the patterns, it would be possible to determine the largely unknown values of the odd viscosity coefficients. This formulation recovers as special cases recent studies of equatorial or topological waves in two-dimensional odd viscous liquids which provided examples of the bulk–interface correspondence at frequencies 𝜔 < 2⁢𝛺. We finally point out that the two- and three-dimensional problems are formally equivalent. Their difference then lies in the way data propagate along characteristic rays in three dimensions, which we demonstrate by classifying the resulting Poincaré–Cartan equations.

Mechanics↗
DOE OSTI2026

Alchemy: A Model-Based Approach for 2D to 3D Autonomous Nuclear System Design

Engineering design of nuclear power plant (NPP) piping and equipment systems frequently bypasses crucial 2D system planning, instead moving straight to 3D modeling. This often leads to designs that exceed building envelope constraints, forcing expensive and time-consuming redesigns. When 2D modeling is employed, it typically involves labor-intensive manual workflows that convert 2D drawings into 3D models, resulting in inefficiencies and errors across design iterations. These workflows further suffer from poor software interoperability and dependence on proprietary software ecosystems, thereby contributing to schedule delays and cost overruns. This paper presents Alchemy, an autonomous framework that transforms 2D system definitions into Industry Foundation Classes (IFC)-compliant 3D building information models (BIMs) for expediting nuclear facility design at the conceptual preliminary phase. Using a model-based approach, the framework treats the 2D system diagram as the central reference model employed to automatically generate all subsequent outputs, ensuring consistency between the system definition and the resulting physical design. A web-based interface enables engineers to define hierarchical system topologies including associated equipment, geometric properties, and connectivity requirements. A two-phase equipment layout optimization algorithm automatically computes collision-free spatial configurations within predefined building envelopes. An artificial intelligence (AI)-assisted pipe routing module then generates orthogonal, collision-free routing paths, allowing the user to select either an A* search-based method or an Ant Colony Optimization (ACO)-based method. All outputs are authored natively in IFC format, relying on open-source technologies and standardized formats in order to ensure extensibility and eliminate proprietary software dependencies. The proposed framework is validated on two representative pressurized-water reactor (PWR)-based case studies, for which it autonomously generates IFC-compliant 3D models in minutes, drastically reducing workflows that typically require hours of manual effort. The generated model demonstrates topologically correct equipment placement, physically plausible spatial relationships, and collision-free pipe routing consistent with known PWR loop configurations. This work represents a foundational step toward digital engineering for nuclear facility preliminary design, with future ongoing development targeting design code compliance and expanded system complexity.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2025

LDRD Abbreviated report: High-Order General-Discrete-Ordinates Method Enabling Efficient Deterministic Transport in Hydrodynamic Simulations

Deterministic transport simulations for national-security and energy applications often operate in high-dimensional phase-space, where accuracy and cost both become major challenges. A common numerical artifact in such problems is the “ray-effect,” which appears as unphysical streaks. Beyond misinterpretation, these artifacts can contaminate tightly coupled physics, such as fluid dynamics, radiation-hydrodynamics, and laser-plasma interactions, eroding the predictive capability of entire multiphysics workflows. Our objective was to make high-dimension studies practical on modern hardware while mitigating the ray-effect without relying on prohibitively expensive sampling approaches such as Monte Carlo methods. We developed the Generic Discretization Library (GenDiL), a Graphics Processing Unit (GPU)-first framework that uses high-order Discontinuous Galerkin (DG) methods and matrix-free algorithms to reduce memory usage and improve computational efficiency, critical for phase-space simulations. GenDiL supports phase-space adaptivity in both mesh size and polynomial order (hp-adaptivity) to place resolution only where it is needed. A central capability is Local Dimensional Refinement (LDR), which couples lower-dimension continuum models to higher-dimension kinetic models through stable and conservative interfaces, so that high-fidelity physics is applied only in regions where it is essential. Building on the GenDiL framework, we developed the General SN (GSN) family of algorithms as a true generalization of the polar SN approach (discrete ordinates, often denoted SN). Rather than tying discrete ordinates to a specific polar change of coordinates, GSN formulates transport on an arbitrary change of coordinates chosen to reduce ray-effect. We studied two complementary variants: an analytic variant, where the coordinate map is prescribed in advance by a closed-form function; and a data-driven variant, where a quantity of interest, such as the net flux, guides the coordinate system. GenDiL provides the library infrastructure for efficient GPU execution, but the GSN concept is algorithmic and independent of any one library. Across representative high-dimension tests, including non-symmetric solutions, both variants delivered strong ray-effect mitigation at practical cost, moving four- to six-dimensional analysis toward repeatable, routine studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Richtmyer–Meshkov instability when a shock is reflected for fluids with arbitrary equation of state

First predicted by Richtmyer in 1960 and experimentally confirmed by Meshkov in 1969, the Richtmyer–Meshkov instability (RMI) is crucial in fields such as physics, astrophysics, inertial confinement fusion and high-energy-density physics. These disciplines often deal with strong shocks moving through condensed materials or high-pressure plasmas that exhibit non-ideal equations of state (EoS), thus requiring theoretical models with realistic fluid EoS for accurate RMI simulations. Approximate formulae for asymptotic growth rates, like those proposed by Richtmyer, are helpful but rely on heuristic prescriptions for compressible materials. These prescriptions can sometimes approximate the RMI growth rate well, but their accuracy remains uncertain without exact solutions, as the fully compressible RMI growth rate is influenced by both vorticity deposited during shock refraction and multiple sonic wave refractions. This study advances previous work by presenting an analytic, fully compressible theory of RMI for reflected shocks with arbitrary EoS. It compares theoretical predictions with heuristic prescriptions using ideal gas, van der Waals gas and three-term constitutive equations for simple metals, the latter being analysed with detailed and simplified ideal-gas-like EoS. We additionally offer an alternative explicit approximate formula for the asymptotic growth rate. The comprehensive model also incorporates the effects of constant-amplitude acoustic waves at the interface, associated with the D'yakov–Kontorovich instability in shocks.

Napieralski, Mario (ORCID:0009000692344901)↗
DOE OSTI2026

Integration of Online Cross-Section Generation Capability with Depletion and Transient Solvers in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE)-based reactor multiphysics analysis application jointly developed by Argonne and Idaho National Laboratories under the DOENE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In FY25, an online crosssection generation capability based on the Self-Shielding Application Programming Interface (SSAPI) was demonstrated for TRISO-fueled reactor problems under steady-state conditions. This fiscal year, that capability was extended to support depletion and transient multiphysics calculations, enabling high-fidelity analyses that generate self-shielded cross sections on the fly from the actual evolving composition and temperature states rather than from pre-tabulated libraries. For depletion, a two-way coupling was established in which SSAPI computes compact-averaged self-shielded cross sections that the depletion solver then uses to advance the Bateman equations, with the updated compositions returned to SSAPI at each step; the depletion module was refactored to support both library-based and SSAPI-based cross sections, and additional logic was added to track daughter isotopes and to exclude minor isotopes for efficiency. For transient analysis, the SSAPI multigroup library was extended with the kinetics data required for time-dependent calculations, the Improved Quasi-Static (IQS) scheme was coupled with SSAPI, and several supporting capabilities were implemented, including a self-shielding treatment that lets control rods and drums move within a self-shielded model, which had previously been impossible and had ruled out rod- and drum-movement transients with on-the-fly cross sections altogether, a new mixing scheme for delayed-neutron precursor decay constants, a checkpoint-based restart workflow, and performance improvements such as pointwise cross-section interpolation and the bypassing of unnecessary Dancoff factor calculations. The implemented capabilities were verified against Serpent Monte Carlo solutions. For depletion, a prismatic pin-cell problem based on a Next Generation Nuclear Plant (NGNP) Very High Temperature Reactor benchmark showed excellent agreement, with eigenvalue differences within 200 pcm over the entire burnup range (up to 140 MWD/kgU) and fission-product and actinide inventories agreeing to within 0.8% and 2.5%, respectively; a heat-pipe microreactor assembly problem with a much higher fuel loading confirmed the same behavior and quantified the bias introduced when the multigroup equivalence effect is neglected. For transient analysis, a pin-cell problem with a step reactivity insertion and temperature feedback reproduced the analytically expected asymptotic power and showed close agreement between the direct and IQS solutions, and a two-dimensional microreactor core problem with control-drum rotation exercised the new moving-drum self-shielding treatment and demonstrated successful coupling of the online crosssection generation with both the direct and IQS transient methods. The capability was further exercised on a full-core pebble-bed problem, in which Griffin was coupled with the System Analysis Module (SAM) to simulate load-following operation of the gPBR with the Doppler feedback resolved at the TRISO fuel kernel temperature. These developments in Griffin provide a convenient, high-fidelity approach to cross-section generation for advanced thermal reactors with geometrically complex and highly heterogeneous configurations, including TRISO-fueled prismatic and pebble-bed systems, and support steady-state, depletion, and transient multiphysics calculations. They also enable self-shielded cross sections to be evaluated directly at the actual coupled state of the system, thereby establishing a foundation for high-fidelity, fully coupled multiphysics analysis of advanced reactors

Park, H.↗
DOE OSTI2023

Interface Design for High‐Performance All‐Solid‐State Lithium Batteries

All‐solid‐state batteries suffer from high interface resistance and lithium dendrite growth leading to low Li plating/stripping Coulombic efficiency (CE) of <90% and low critical current density at high capacity. Here, in this work, both challenges are simultaneously addressed and the Li plating/stripping CE is significantly increased to 99.6% at 0.2 mA cm −2 /0.2 mAh cm −2 , and critical current density (CCD) of > 3.0 mA cm −2 /3.0 mAh cm −2 by inserting a mixed ionic‐electronic conductive (MIEC) and lithiophobic LiF‐C‐Li 3 N‐Bi nanocomposite interlayer between Li 6 PS 5 Cl electrolyte and Li anode. The highly lithiophobic LiF‐C‐Li 3 N‐Bi interlayer with high ionic conductivity (10 −5 S cm −1 ) and low electronic conductivity (3.4×10 −7 S cm −1 ) enables Li to plate on the current collector (CC) surface rather than on Li 6 PS 5 Cl surface avoiding Li 6 PS 5 Cl electrolyte reduction. During initial Li plating on CC, Li penetrates into porous LiF‐C‐Li 3 N‐Bi interlayer and lithiates Bi nanoparticles into Li 3 Bi. The lithiophilic Li 3 Bi and Li 3 N nanoparticles in LiF‐C‐Li 3 N‐Li 3 Bi sub‐interlayer will move to CC along with plated Li, forming LiF‐C/Li 3 N‐Li 3 Bi lithiophobic/lithiophilic sublayer during the following Li stripping. This interlayer enables Co 0.1 Fe 0.9 S 2 /Li 6 PS 5 Cl/Li cell with an areal capacity of 1.4 mAh cm −2 to achieve a cycle life of >850 cycles at 150 mA g −1 . The lithiophobic/lithiophilic interlayer enables solid‐state metal batteries to simultaneously achieve high energy and long cycle life.

25 ENERGY STORAGE↗
DOE OSTI2026

Toward Trustworthy Autonomous Science: A Two-Year Community Roadmap

One year ago, the AISLE roadmap argued that autonomous laboratories operated as isolated islands and proposed a grassroots network organized around five critical dimensions. The field has since moved faster than that roadmap anticipated: multi-agent systems have produced experimentally validated hypotheses, self-driving laboratories have grown more interoperable and orchestrated, reasoning-trained and domain foundation models have raised the capability ceiling, and the Genesis Mission has placed autonomous experimentation at the center of U.S. federal science strategy, with industry emerging as a primary actor. Progress has met a sobering counter-current, including a corrected flagship discovery result, benchmarks showing that agents which rival experts on closed-ended questions still complete only a fraction of open-ended research, and fabricated citations surfacing at leading venues. We read this as the defining tension of the field: producing a candidate discovery is no longer the hard part, but verifying it is, and this asymmetry now limits autonomous science more than raw model capability. Accordingly, we update the roadmap around seven dimensions, revisiting the original five and elevating two former cross-cutting concerns, trust, verification, and reproducibility, and safety, security, and governance, to first-class status. We assess the original milestones (M1 through M14) as achieved, partially achieved, reframed, or open, add four new milestones (M15 through M18) for the elevated dimensions, and scope the path forward to a two-year horizon, with the first year concentrating on interfaces, protocol adoption, and the scaffolding of verification, and the second targeting federation, zero-trust coordination, and governance. Throughout, we position the grassroots network as the interoperability fabric that lets national programs, international initiatives, and commercial platforms connect rather than re-silo.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2026

Enhancing photoionization rate calculations in low-temperature plasmas using spectral methods

Photoionization plays a central role in the development of streamer discharges and other non-equilibrium plasma phenomena. It creates seed electrons, which are essential for positive streamer propagation, allowing the ionization front to move forward. Because of this, accurate modeling of photoionization is very important for predicting streamer behavior and plasma evolution. The photoionization process in air (N 2 – O 2 mixture) is often described by the Zheleznyak model (1982). This model is usually solved through Helmholtz-type equations that approximate the Zheleznyak photoionization model (Zheleznyak et al. 1982) as Partial Differential Equations (PDEs). Conventional numerical methods, such as the Finite Difference Method (FDM) or Finite Volume Method (FVM), are widely used to solve these equations. Although they are prevalent, the computational cost of these methods due to their need for matrix operations and iterative solver is demanding. To address this challenge, this work develops a spectral solver based on the Fast Fourier Transform (FFT) combined with Discrete Cosine Transform (DCT) and Discrete Sine Transform (DST) to calculate the photoionization rate efficiently in an axisymmetric cylindrical domain. This method naturally satisfies the boundary conditions used in the model and converts the PDE into algebraic ones in spectral space. Thus, avoids the need for iterative matrix solvers. When compared with FDM results, it is demonstrated that the new solver not only maintains accuracy, but also reduces the computational cost, showing a performance increase of approximately 100 compared to FDM over a wide range of problem sizes. The method is parallelized using Message Passing Interface (MPI) and has been integrated into a fluid plasma model for streamer simulation. Here, this FFT-based approach provides a fast and reliable alternative for calculating photoionization in fluid models, helping large-scale plasma simulations run faster and efficiently, and allows higher-resolution simulation without extra computational cost.

Axisymmetric system↗
DOE OSTI2025

MontePy: a Python library for reading, editing, and writing MCNP input files.

The Monte Carlo N-Particle (MCNP) radiation transport code is a highly capable and accurate code with a long legacy. MCNP uses the Monte Carlo simulation process to simulate the path of particles (e.g., neutrons, photons, charged particles, etc.), and their interaction with materials. It is widely used in nuclear engineering, high-energy physics, and other fields. Its origins in the mid-twentieth century predate many modern software conventions. MCNP users provide an input file to MCNP, which it then uses to create an internal representation of the simulation problem. These input files originally had to be stored as punchcard decks, and the user manual still uses the terminology of cards and decks, despite moving beyond punchcards. MCNP predates nearly all modern human readable markup or data serialization languages, such as the extensible Markup Language (XML), the Standard Generalized Markup Language (SGML), YAML (YAML Ain’t Markup Language), and Javascript Object Notation (JSON). Due to this, MCNP uses an entirely custom defined syntax language for its input, making off-the-shelf libraries for XML, YAML, and JSON impossible to use for scripting various operations on MCNP input files (Kulesza et al., 2022).

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Innovating the next generation of commercial smart building software

Nearly 30% of commercial building energy use is wasted due to equipment faults and HVAC controls problems. The result is increased emissions, compromised comfort and productivity, and less reliable coordination of building power needs with a clean grid. The energy impact alone represents $17 billion in potential savings. Today’s smart building software provides a robust solution to address these operational deficiencies. Energy management and information systems (EMIS) are saving up to 9% on average, with two-year paybacks. They are being incorporated into energy management processes, commissioning services, and utility programs. As effective as they are, two barriers prevent even deeper benefits; limited personnel to fix problems once they are identified, and the expense and time to manually implement changes in control systems. In partnership with the research community, the EMIS industry is developing new capabilities to overcome these barriers. Moving beyond siloed products for either fault detection and diagnostics, or optimal control, these new capabilities empower users to not only automatically identify faults, but also to push corrective action, and control improvements to their buildings. In this paper, several areas for enhancements are documented: ‘one-time’ correction of faults such as setpoints, schedules, and economizer lockouts; short-term active testing for automated proportional integral derivative (PID) loop tuning and functional testing; and continuous supervisory control for demand flexibility and year-round efficiency. Results are presented from a pair of partner implementations out of a dozen providers integrating these enhancements into their products, including field tests from across the country, and insights into operator acceptance and integration into operations and maintenance practices.

Casillas, Armando↗