SEARCH · Engineering Papers

Results for “reversible computing”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 919 records · Page 51

All sources NASA NTRS DOE OSTI

NASA NTRS2012

Microwave Properties of Ice-Phase Hydrometeors for Radar and Radiometers: Sensitivity to Model Assumptions

A simplied framework is presented for assessing the qualitative sensitivities of computed microwave properties, satellite brightness temperatures, and radar reflectivities to assumptions concerning the physical properties of ice-phase hydrometeors. Properties considered included the shape parameter of a gamma size distribution andthe melted-equivalent mass median diameter D0, the particle density, dielectric mixing formula, and the choice of complex index of refraction for ice. We examine these properties at selected radiometer frequencies of 18.7, 36.5, 89.0, and 150.0 GHz; and radar frequencies at 2.8, 13.4, 35.6, and 94.0 GHz consistent with existing and planned remote sensing instruments. Passive and active microwave observables of ice particles arefound to be extremely sensitive to the melted-equivalent mass median diameter D0 ofthe size distribution. Similar large sensitivities are found for variations in the ice vol-ume fraction whenever the geometric mass median diameter exceeds approximately 1/8th of the wavelength. At 94 GHz the two-way path integrated attenuation is potentially large for dense compact particles. The distribution parameter mu has a relatively weak effect on any observable: less than 1-2 K in brightness temperature and up to 2.7 dB difference in the effective radar reflectivity. Reversal of the roles of ice and air in the MaxwellGarnett dielectric mixing formula leads to a signicant change in both microwave brightness temperature (10 K) and radar reflectivity (2 dB). The choice of Warren (1984) or Warren and Brandt (2008) for the complex index of refraction of ice can produce a 3%-4% change in the brightness temperature depression.

Microwave properties↗

NASA NTRS2021

Analysis of Exercise Loads to Inform Vibration Isolation System Design

BACKGROUND: This study was conducted with the primary interest of providing data that would inform Vibration Isolation and Stabilization (VIS) system design and performance for the European Enhanced Exploration Exercise Device (E4D). In preparation for the International Space Station (ISS) in-flight demonstration, a list of critical Human Health Countermeasures (HHC) exercises was compiled [1]. The goal of this study was to assess the ground reaction forces and moments imposed by an exercising subject in each of the six VIS Degrees of Freedom (DOFs) during a comprehensive set of these critical exercises performed on the E4D. METHODS AND RESULTS: The ISS in-flight demonstration list of critical exercises included seated aerobic rowing, bent-over rowing, cycling, front squats, back squats, conventional deadlifts, Romanian deadlifts, heel raises, overhead presses, reverse chops, and power clean presses. At the NASA Johnson Space Center (JSC) Prototype Immersive Technology (PIT) laboratory, motion capture data were collected on critical E4D exercises for six subjects. At the NASA JSC Active Response Gravity Offload System (ARGOS) facility, additional motion capture and load cell data were collected on offloaded trials for four subjects. Select data were extrapolated to represent a 5th percentile female subject and a 95th percentile male subject. A previous investigation comparing the forces obtained from the load cell and from motion capture based data found a satisfactory level of agreement between the two measurements [2]. The motion capture based data were analyzed for this study since it is driven by the subject’s trajectory alone, automatically excluding any forces exerted on the subject by the ARGOS offloading harness. The OpenSim [3, 4] biomechanical simulation inverse kinematics tool was used to calculate the joint angles based on the locations of motion capture markers placed at key positions on the subject’s body. An OpenSim plugin was then used to obtain the forces and moments generated by the subject during each trial, with the moments computed relative to the equilibrium location of the subject’s feet [5]. The force of gravity was also removed to simulate the loads generated by the exercise when performed in microgravity. The load plots for each trial were generated and visually analyzed to obtain the magnitudes of the peak loads for each exercise in each DOF. The typical period of exercise for each trial was also estimated and used to calculate the frequency for each trial. The exercise loads data was then organized in multiple ways to capture different aspects of the data. As a result of this study, we present a summary of the load magnitudes observed during these critical exercises utilizing the E4D.

C A Bell↗

DOE OSTI2025

Risks and Benefits of Pressurized Water Reactor Coupling to Energy Storage and Water Desalination

This research explored potential options to improve nuclear pressurized water reactors’ (PWRs) operational flexibility through integration with thermal energy storage (TES) and water desalination technologies. Benefits and challenges of TES technologies were reviewed, which includes sensible heat storage systems, latent heat storage systems, and thermochemical storage systems. Water desalination technologies were reviewed and compared, which includes multi-stage flash, reverse osmosis, and multi-effect distillation. The multi-effect distillation technology was selected to be coupled with a generic PWR. Risk of this cogeneration was quantified using the probabilistic risk assessment (PRA) methodology. Results suggest that nuclear-desalination cogeneration can be done safely because the additional risk to the PWR is trivial, and even in certain cases the PWR risk is decreased. This cogeneration operation is expected to improve efficiency of thermal utilization and the nuclear plant’s revenues. In summary, this study supports the potential of nuclear-desalination cogeneration systems to address water scarcity while optimizing nuclear energy use and managing safety.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Modification and analysis of context-specific genome-scale metabolic models: methane-utilizing microbial chassis as a case study

ABSTRACT Context-specific genome-scale model (CS-GSM) reconstruction is becoming an efficient strategy for integrating and cross-comparing experimental multi-scale data to explore the relationship between cellular genotypes, facilitating fundamental or applied research discoveries. However, the application of CS modeling for non-conventional microbes is still challenging. Here, we present a graphical user interface that integrates COBRApy, EscherPy, and RIPTiDe, Python-based tools within the BioUML platform, and streamlines the reconstruction and interrogation of the CS genome-scale metabolic frameworks via Jupyter Notebook. The approach was tested using -omics data collected for Methylotuvimicrobium alcaliphilum 20Z R , a prominent microbial chassis for methane capturing and valorization. We optimized the previously reconstructed whole genome-scale metabolic network by adjusting the flux distribution using gene expression data. The outputs of the automatically reconstructed CS metabolic network were comparable to manually optimized i IA409 models for Ca-growth conditions. However, the CS model questions the reversibility of the phosphoketolase pathway and suggests higher flux via primary oxidation pathways. The model also highlighted unresolved carbon partitioning between assimilatory and catabolic pathways at the formaldehyde-formate node. Only a very few genes and only one enzyme with a predicted function in C1 metabolism, a homolog of the formaldehyde oxidation enzyme ( fae1-2 ), showed a significant change in expression in La-growth conditions. The CS-GSM predictions agreed with the experimental measurements under the assumption that the Fae1-2 is a part of the tetrahydrofolate-linked pathway. The cellular roles of the tungsten (W)-dependent formate dehydrogenase ( fdhAB ) and fae homologs ( fae1-2 and fae3 ) were investigated via mutagenesis. The phenotype of the f dhAB mutant followed the model prediction. Furthermore, a more significant reduction of the biomass yield was observed during growth in La-supplemented media, confirming a higher flux through formate. M. alcaliphilum 20Z R mutants lacking fae1-2 did not display any significant defects in methane or methanol-dependent growth. However, contrary to fae1, the fae1-2 homolog failed to restore the formaldehyde-activating enzyme function in complementation tests. Overall, the presented data suggest that the developed computational workflow supports the reconstruction and validation of CS-GSM networks of non-model microbes. IMPORTANCE The interrogation of various types of data is a routine strategy to explore the relationship between genotype and phenotype. An efficient approach for integrating and cross-comparing experimental multi-scale data in the context of whole-genome-based metabolic network reconstruction becomes a powerful tool that facilitates fundamental and applied research discoveries. The present study describes the reconstruction of a context-specific (CS) model for the methane-utilizing bacterium, Methylotuvimicrobium alcaliphilum 20Z R . M. alcaliphilum 20Z R is becoming an attractive microbial platform for the production of biofuels, chemicals, pharmaceuticals, and bio-sorbents for capturing atmospheric methane. We demonstrate that this pipeline can help reconstruct metabolic models that are similar to manually curated networks. Furthermore, the model is able to highlight previously overlooked pathways, thus advancing fundamental knowledge of non-model microbial systems or promoting their development toward biotechnological or environmental implementations.

Kulyashov, M. A.↗

DOE OSTI2023

Octahedral to tetrahedral bonding transitions in the local structure of phase change optical media Ge 2 Sb 2 Se 5 x Te 5-5 x with Se doping

Random access memories utilize fast, reversible switching between ordered and disordered states of matter in phase change materials (PCMs) such as Ge2Sb2Te5-5x. The short-range structure in the disordered phase has been described either as (i) a network of Ge tetrahedra or (ii) Peierls distorted Ge/Sb octahedra. The PCM transition was investigated in bulk Ge2Sb2Se5xTe5-5x (GSST), in which amorphization sets in with Se doping (x ≈ 0.85) upon quenching. GSST has a hexagonal crystalline ground state with Ge/Sb octahedral coordination, but the phase change transition to the amorphous state that is only observed when the system is quenched brings a short-range structure with sharp, tetrahedrally coordinated Ge/Sb correlations and shortened bonds that are distinctly different from the expected octahedral pairing.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

RAMSeS: Rapid Analysis of Mission Software Systems

Over the past few decades, software has become ubiquitous as it has been integrated into nearly every aspect of society, including household appliances, consumer electronics, industrial control systems, public utilities, government operations, and military systems. Consequently, many critical national security questions can no longer be answered convincingly without understanding software, including its purpose, its capabilities, its flaws, its communication, or how it processes and stores data. As software continues to become larger, more complex, and more widespread, our ability to answer important mission questions and reason about software in a timely way is falling behind. Today, to achieve such understanding of third-party software, we rely predominantly on the ability of reverse engineering experts to manually answer each particular mission question for every software system of interest. This approach often requires heroic human effort that nevertheless fails to meet current mission needs and will never scale to meet future needs. The result is an emerging crisis: a massive and expanding gap between the national security need to answer mission questions about software and our ability to do so. Sandia National Laboratories has established the Rapid Analysis of Mission Software Systems (RAMSeS) effort, a collaborative long-term effort aimed at dramatically improving our nation’s ability to answer mission questions about third-party software by growing an ecosystem of tools that augment the human reverse engineer through automation, interoperability, and reuse. Focusing on static analysis of binary programs, we are attempting to identify reusable software analysis components that advance our ability to reason about software, to automate useful aspects of the software analysis process, and to integrate new methodologies and capabilities into a working ecosystem of tools and experts. We aim to integrate existing tools where possible, adapt tools when modest modifications will enable them to interoperate, and implement missing capability when necessary. Although we do hope to automate a growing set of analysis tasks, we will approach this goal incrementally by assisting the human in an ever-widening range of tasks.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2026

A New Route Toward Atomically Flat and Defect-Free Ge/SiGe Planar Heterostructures

Germanium-based planar heterostructures are emerging as versatile platforms for realizing quantum devices. In particular, planar Ge/SiGe quantum wells (QWs) host hole states with exceptionally high mobility and strong, electrically tunable spin–orbit interactions, enabling full electrical control of quantum information. A key requirement for these systems is the growth of microscopic, defect-free, atomically flat Ge QW heterostructures on relaxed or reverse-graded SiGe buffer layers, as well as on commercial Ge substrates. While several physical deposition techniques have demonstrated high-quality planar Ge/SiGe QWs, a major challenge remains: minimizing defect density typically requires high growth temperatures, which are incompatible with standard CMOS process flows. Here, we present a convenient low-temperature process for realizing high-quality planar SiGe/Ge heterostructures using a combination of thermal and electron-beam evaporation. We systematically map the effects of ex-situ and in-situ substrate preparation protocols, growth temperature, and post-deposition annealing conditions, and correlate these parameters with surface roughness and defect density. We find that in-situ oxide desorption conditions and post-annealing parameters have the most pronounced impact on improving surface quality. Under optimized conditions, we achieve atomically smooth surfaces with root-mean-square roughness σrms ≤ 10 Å and negligible defect density. Interestingly, thermally evaporated Ge layers exhibit oriented triangular crystallites that elongate upon post-annealing in the presence of high Ge vapor pressure. These results demonstrate that this simple, low-temperature deposition approach is a viable and effective route for achieving high-quality Ge QW heterostructures, with strong potential for scalable quantum computing and sensing applications.

Tripathi, Malvika [Fermilab] (ORCID:00000001989251↗

DOE OSTI2026

A New Route Toward Atomically Flat and Defect-Free Ge/SiGe Planar Heterostructures

Tripathi, Malvika [Fermilab] (ORCID:00000001989251↗

DOE OSTI2025

Emergence of Moiré Dirac Fermions at the Interface of Topological and 2D Magnetic Insulators

Dirac Fermions on the surface of the topological insulator are spin-momentum locked and topologically protected, making them interesting for spintronics and quantum computing applications. When in proximity to magnetism and superconductivity, these electronic states could result in quantum anomalous Hall effect and Majorana Fermions, respectively. An even more dramatic enrichment of the topological insulators’ physics is expected for moiré superlattices, where, analogously to the twisted graphene layers, electronic correlations could be strongly enhanced, a task previously notoriously difficult to achieve in topological matter. Until now, the experimental confirmation of such moiré properties has remained elusive. Here, we grow the two-dimensional van der Waals magnetic insulators FeX 2 (where X = Cl or Br) on top of the topological insulator Bi 2 Se 3 and establish a moiré superlattice formation at the interface. By means of scanning tunneling microscopy and angle-resolved photoemission spectroscopy, we investigate the electronic properties of the formed moiré superlattice and demonstrate its tunability via the film choice. We reveal replicated Dirac cones and focus on their intersections, which, in the case of FeBr 2 /Bi 2 Se 3 , occur below the Fermi level. We identify the signatures of small gaps at the intersections around the M̅ i points that we attribute to the moiré interaction. These findings point to the specific type of magnetic moiré potential that breaks the time-reversal symmetry at these points but not at the $\barΓ$ point. Our observations provide an intriguing scenario of correlated topological phases induced by moiré superlattice that may result in topological superconductivity, high Chern number phases, and exotic noncollinear magnetic textures.

2D magnets↗

DOE OSTI2023

Spatially Distributed Ramp Reversal Memory in VO 2

Ramp‐reversal memory has recently been discovered in several insulator‐to‐metal transition materials where a non‐volatile resistance change can be set by repeatedly driving the material partway through the transition. This study uses optical microscopy to track the location and internal structure of accumulated memory as a thin film of VO 2 is temperature cycled through multiple training subloops. These measurements reveal that the gain of insulator phase fraction between consecutive subloops occurs primarily through front propagation at the insulator‐metal boundaries. By analyzing transition temperature maps, it is found, surprisingly, that the memory is also stored deep inside both insulating and metallic clusters throughout the entire sample, making the metal‐insulator coexistence landscape more rugged. This non‐volatile memory is reset after heating the sample to higher temperatures, as expected. Diffusion of point defects is proposed to account for the observed memory writing and subsequent erasing over the entire sample surface. By spatially mapping the location and character of non‐volatile memory encoding in VO 2 , this study results enable the targeting of specific local regions in the film where the full insulator‐to‐metal resistivity change can be harnessed in order to maximize the working range of memory elements for conventional and neuromorphic computing applications.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Bead–Spring Simulation of Ionomer Melts—Studying the Effects of Chain-Length and Associating Group Fraction on Equilibrium Structure and Extensional Flow Behavior

Ionomers are associative polymers with diverse applications ranging from selective membranes and high-performance adhesives to abrasion- and chemical-resistant coatings, insulation layers, vacuum packaging, and foamed sheets. Within equilibrium melt, the ionic or associating groups are known to form thermally reversible, associative clusters whose presence can significantly affect the system’s mechanical, viscoelastic, and transport properties. It is, thus, of great interest to understand how to control such clusters’ size distribution, shape, and stability through the designed choice of polymer architecture and the ionic groups’ fraction, arrangement, and interaction strength. In this work, we represent linear associating polymers using a Kremer–Grest type bead–spring model and perform large-scale MD simulations to explore the effect of polymer chain-length (l) and fraction (f s ) of randomly placed associating groups on the size distribution and stability of formed clusters. We consider different chain-lengths (below and above entanglement), varying fractions of associating groups (represented by ‘sticky’ beads) between 5 and 20%, and a fixed sticky–sticky nonbond interaction strength of four times that between regular non-associating beads. For all melts containing associating groups the equilibrium structure factor S(q) displays a signature ionomer peak at low wave vector q whose intensity increases with increasing f s and l. The average cluster size Nc increases with f s . However, the effect of chain-length on Nc appears to be pronounced only at higher values of f s . Under extensional flows, the computed stress (and viscosity) is higher at higher f s and l regardless of strain rate. Beyond a critical strain rate, we observe fragmentation of the associative clusters, which has interesting effects on the stress/viscous response.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Redox Poise during Rhodospirillum rubrum Phototrophic Growth Drives Large-scale Changes in Macromolecular Synthesis Pathways

During photoheterotrophic growth on organic substrates, purple nonsulfur photosynthetic bacteria like Rhodospirillum rubrum can acquire electrons by multiple means, including oxidation of organic substrates, oxidation of inorganic electron donors (e.g., H2), and by reverse electron flow from the photosynthetic electron transport chain. These electrons are stored as reduced electron-carrying cofactors (e.g., NAD(P)H and ferredoxin). The overall ratio of oxidized to reduced cofactors (e.g., NAD(P)+:NAD(P)H), or ’redox poise’, is difficult to understand or predict, as are the cellular processes for dissipating these reducing equivalents. Using physics-based models that capture mass action kinetics consistent with the thermodynamics of reactions and pathways, a range of redox conditions for heterophototrophic growth are evaluated, from conditions in which the NADP+/NADPH levels approach thermodynamic equilibrium to conditions in which the NADP+/NADPH ratio is far above the typical physiological values. Modeling predictions together with experimental measurements indicate that the redox poise of the cell results in large-scale changes in the activity of biosynthetic pathways and, thus, changes in cell macromolecule levels (DNA, RNA, proteins, and fatty acids). Furthermore, modeling predictions indicate that during phototrophic growth, reverse electron flow from the quinone pool is a minor contributor to the production of reduced cofactors (e.g., NAD(P)H) compared to other oxidative processes (H2 and carbon substrate oxidation). Instead, the quinone pool primarily operates to aid ATP production. The high level of ATP, in turn, drives reduction processes even when NADPH levels are relatively low compared to NADP+ by coupling ATP hydrolysis to the reductive processes. The model, in agreement with experimental measurements of macromolecule ratios of cells growing on different carbon substrates, indicates that the dynamics of nucleotide versus lipid and protein production is likely a significant mechanism of balancing oxidation and reduction in the cell.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS2001

Experimental and Computational Studies of the Control of Convection of Non-Conducting Liquids During solidification by Means of a Magnetic Field Gradient

The elimination of convection is essential in experimental investigations of diffusive transport (of heat and matter) during solidification. One classical approach to damping convection in a conducting liquid is the application of a magnetic field. The damping phenomenon is the induction, by the motion of a conductor in a magnetic field, of currents which interact with the field to produce Lorentz forces that oppose the flow. However, there are many liquids which are not sufficiently conducting to exploit this phenomenon; examples include the transparent liquids (such as succinonitrile-acetone) that are used as "model alloys" in fundamental solidification studies. There have been several investigations of the solidification of these liquids that have been carried out in orbiting laboratories to eliminate natural convection. The paper describes an investigation of an alternative approach whereby a magnetic field gradient is applied to the liquid. A magnetic body force then arises which is dependent on the susceptibility of the liquid and thereby on the temperature and or concentration. With the field gradient aligned vertically and of correct magnitude, the variation of gravitational body force due to temperature/concentration dependent density can be counterbalanced by a variation in magnetic body force. Experiments have been carried out in a super-conducting magnet at Marshall Space Flight Center to measure velocities in an aqueous manganese chloride solution. The solution was contained in a chamber with temperature controlled end walls and glass side walls. Velocities were measured by particle image velocimetry. Starting from zero current in the magnet (zero field gradient) flow driven by the temperature difference between the end walls was measured. At a critical current the flow was halted. At higher currents the normal convection was reversed. The experiments included ones where the solution was solidified and were accompanied by solution of the flow/transport equations using the software package FLUENT.

Seybert, C. D.↗

NASA NTRSDate not supplied

Capturing, Analyzing, Maintaining, and Disseminating Shape Memory Material Data Between Information Management Systems

With an increased demand on reducing the time, cost, and effort to develop new materials, Integrated Computational Materials Engineering (ICME) has received widespread attention in various engineering disciplines as a catalyst for significantly reducing experimental testing during the material design process. An ICME approach to design can enable ‘fit-for-purpose’ materials to be realized in engineering applications by incorporating well-understood process-property-performance relationships between the various length and time scales in a material’s structure, enabling material optimization. However, such an approach requires validated multiscale models at the various length scales for a material, which in turn requires a large amount of data, a robust means of storing the data, and the ability to link data to developed material models. The NASA Vision 2040 [1] has identified nine key elements to enabling ICME approaches in system level design, with one being “Data, Information, and Visualization”, thus outlining the importance of a robust information management system for ICME. As the relationship between microstructure, properties, and material performance become better understood and incorporated into multiscale models that can be leveraged in application design, the emergence of new materials with application-driven properties can be realized. One such new material class that has seen growing attention are shape memory materials (SMM), in which a material can transition between a deformed and undeformed state via a reversible phase transformation when subject to a thermal, mechanical, or magnetic load [2]. SMMs have been used widely in aerospace and biomedical industries, including applications such as actuators, low-shock mechanisms, medical staples, braces, and stents [3, 4]. These materials exhibit unique behavior due to their ability to transition between phases, and thus the mechanisms that enable this transition must be captured in a data information management system and incorporated into SMM material models. At NASA Glenn Research Center, the Shape Memory Materials Database (SMMD) Tool has been developed to capture the necessary information that governs SMM material behavior and provide users the ability to select and visualize various SMMs for a specific application [5]. The database contains point-wise data for published SMM materials, along with the pedigree metadata for traceability necessary for a robust information management system. The database is also capable of storing in-house test data performed at NASA GRC by interacting with the developed Shape Memory Alloy (SMA) Analytics tool to extract the necessary point-wise values and populate the database. Although the SMMD Tool offers its users a single, authoritative source for SMM material data that is critical for model development and material design, the full material pedigree of the in-house test data for SMMs is not currently captured and is out of the scope for the SMMD tool. In this work, the schema for capturing SMM test data within the larger NASA GRC ICME Schema [6, 7, 8, 9] will be developed and implemented for thermomechanical tests conducted at NASA GRC. The developed schema will not only store the relevant data needed for the SMMD tool, but also the material pedigree (i.e., production of the bulk material, bulk material analysis, sample cut-out diagrams, sample fabrication procedure, etc.), test pedigree (i.e., test equipment used, measurement systems used, raw test data), and analysis pedigree (i.e., how the data in the SMMD tool is calculated). Furthermore, a Python-based framework will be developed to seamlessly interact between the SMA Analytics and SMMD tools, which will write the full dataset and associated metadata to the GRC Information Management System before passing the required point-wise data to the SMMD tool. Data informatics is a key element of the NASA Vision 2040, which requires not only that data is stored and maintained throughout the material lifecycle, but that the data is also accessible and reusable such that material development efforts can be minimized. Therefore, for an ICME design approach to be realized, a centralized information management system that drives the ICME process must be able to communicate with other databases. The work that will be presented in this presentation will therefore not only demonstrate the ability of NASA GRC’s information management system to capture SMM data, but also its ability to interact with pre-existing tools specialized for such materials.

Data management↗

DOE OSTI2024

Validation of a Hybrid Domain Overlapping Coupling Between SAM and CFD Against the TALL-3D Transients

The System Thermal Hydraulics (STH) code SAM has been coupled to the Computational Fluid Dynamics (CFD) code Simcenter STAR-CCM+ utilizing a hybrid domain overlapping method with an explicit coupling in time. The coupling aims to extend the STH code’s applicability to scenarios where local momentum and energy transfers are important yet difficult for STH codes to capture, such as three-dimensional (3D) mixing. The coupling method’s numerical stability was verified in the past against two closed-loop configurations, and it was validated against a double T-junction experiment with 3D scalar mixing. In the present work, the coupling method is validated against the TALL-3D STH/CFD coupling benchmark facility. TALL-3D is a three-legged, liquid-metal facility with a large, pool-type enclosure (test section) that exhibits 3D flow effects to be modeled by a CFD code. The rest of the system exhibits approximately 1D behavior well-predicted by an STH code. First, the present STAR-CCM+ CFD model of the 3D test section is validated against experimental data. Then, the SAM-STARCCM+ coupled model is validated against six different TALL-3D steady states, including SAM standalone model results for comparison. Lastly, the SAM-STARCCM+ coupled model is validated against two TALL-3D transients, one exhibiting flow reversal in the test section and one exhibiting nonlinear, Limit Cycle Oscillations (LCO). For the first transient, the SAM-STARCCM+ coupled model properly predicts an increase in the test section’s inlet temperature during flow reversal, and this is not predicted by the SAM standalone model. Following flow reversal, the SAM-STARCCM+ coupled model better-predicts the initial flow recovery and following oscillations as the system approaches a final natural circulation state. For the second transient, no true final steady state is observed due to LCO. Neither the SAM-STARCCM+ coupled model nor the SAM standalone model can perfectly capture the experiment’s changing oscillation frequency during the transient. However, the SAM-STARCCM+ coupled model does reproduce the oscillatory feedback observed in the system. This is a significant achievement as the SAM-STARCCM+ coupled model only uses an explicit coupling in time, as opposed to a semi-implicit coupling. In comparison, previous STH/CFD coupling efforts of the TALL-3D facility required semi-implicit coupling to obtain similar results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Efficient ptychographic phase retrieval via a matrix-free Levenberg-Marquardt algorithm

The phase retrieval problem, where one aims to recover a complex-valued image from far-field intensity measurements, is a classic problem encountered in a range of imaging applications. Modern phase retrieval approaches usually rely on gradient descent methods in a nonlinear minimization framework. Calculating closed-form gradients for use in these methods is tedious work, and formulating second order derivatives is even more laborious. Additionally, second order techniques often require the storage and inversion of large matrices of partial derivatives, with memory requirements that can be prohibitive for data-rich imaging modalities. We use a reverse-mode automatic differentiation (AD) framework to implement an efficient matrix-free version of the Levenberg-Marquardt (LM) algorithm, a longstanding method that finds popular use in nonlinear least-square minimization problems but which has seen little use in phase retrieval. Furthermore, we extend the basic LM algorithm so that it can be applied for more general constrained optimization problems (including phase retrieval problems) beyond just the least-square applications. Since we use AD, we only need to specify the physics-based forward model for a specific imaging application; the first and second-order derivative terms are calculated automatically through matrix-vector products, without explicitly forming the large Jacobian or Gauss-Newton matrices typically required for the LM method. We demonstrate that this algorithm can be used to solve both the unconstrained ptychographic object retrieval problem and the constrained “blind” ptychographic object and probe retrieval problems, under the popular Gaussian noise model as well as the Poisson noise model. We compare this algorithm to state-of-the-art first order ptychographic reconstruction methods to demonstrate empirically that this method outperforms best-in-class first-order methods: it provides excellent convergence guarantees with (in many cases) a superlinear rate of convergence, all with a computational cost comparable to, or lower than, the tested first-order algorithms.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Virtual Time III, Part 1: Unified Virtual Time Synchronization for Parallel Discrete Event Simulation

Algorithms for synchronization of parallel discrete event simulation have historically been divided between conservative methods that require lookahead but not rollback, and optimistic methods that require rollback but not lookahead. In this paper we present a new approach in the form of a framework called Unified Virtual Time (UVT) that unifies the two approaches, combining the advantages of both within a single synchronization theory. Whenever timely lookahead information is available, a logical process (LP) executes conservatively using an irreversible event handler. When lookahead information is not available the LP does not block, as it would in a classical conservative execution, but instead executes optimistically using a reversible event handler. The switch from conservative to optimistic synchronization and back is decided on an event-by-event basis by the simulator, transparently to the model code. UVT treats conservative synchronization algorithms as optional accelerators for an underlying optimistic synchronization algorithm, enabling the speed of conservative execution whenever it is applicable, but otherwise falling back on the generality of optimistic execution. We describe UVT in a novel way, based on fundamental invariants, monotonicity requirements, and synchronization rules. UVT permits zero-delay messages and pays careful attention to tie-handling using superposition. We prove that under fairly general conditions a UVT simulation always makes progress in virtual time. This is Part 1 of a trio of papers describing the UVT framework for PDES, mixing conservative and optimistic synchronization and integrating throttling control.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Superstructure Optimization for Brine Valorization from Brackish Water Desalination

This poster presents preliminary results from a superstructure optimization framework developed to identify cost-optimal brine valorization configurations for brackish water desalination plants across diverse U.S. regional feed chemistries. The study uses brackish groundwater compositions from Arizona, California, Florida, New Mexico, and Texas. Using Pyomo Generalized Disjunctive Programming (GDP) within the WaterTAP modeling environment, the optimization framework simultaneously evaluates thousands of candidate treatment configurations, spanning nanofiltration, reverse osmosis, and chemical precipitation, to minimize the levelized cost of water (LCOW) while meeting water recovery targets and product recovery constraints. Results across eight representative feed clusters demonstrate water recovery rates of 57–87% and net LCOW values ranging from -$0.032/m³ (net revenue-positive) to $0.80/m. Notably, no single process configuration was optimal across all feed types, underscoring the necessity of feed-specific optimization. Products targeted include calcium carbonate (CaCO₃) at $0.01/kg and sodium chloride (NaCl) at $0.10/kg, both at 95% purity, with product revenues offsetting treatment costs in several scenarios. The work advances NAWI's process systems engineering capabilities for multi-configuration screening.

58 GEOSCIENCES↗