Search NASA⌕ Search

SEARCH · Search NASA

Results for “interface structure prediction”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Likelihood-based interactive local docking into cryo-EM maps in ChimeraX

The interpretation of cryo-EM maps often includes the docking of known or predicted structures of the components, which is particularly useful when the map resolution is worse than 4 Å. Although it can be effective to search the entire map to find the best placement of a component, the process can be slow when the maps are large. However, frequently there is a well-founded hypothesis about where particular components are located. In such cases, a local search using a map subvolume will be much faster because the search volume is smaller, and more sensitive because optimizing the search volume for the rotation-search step enhances the signal to noise. A Fourier-space likelihood-based local search approach, based on the previously published em_placement software, has been implemented in the new emplace_local program. Tests confirm that the local search approach enhances the speed and sensitivity of the computations. An interactive graphical interface in the ChimeraX molecular-graphics program provides a convenient way to set up and evaluate docking calculations, particularly in defining the part of the map into which the components should be placed.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Area-based composition predictions of materials fabricated using simultaneous wire-powder-directed energy deposition

Functionally graded materials are an emergent method for designing components with programmable site-specific material properties. These materials are typically fabricated using metal additive manufacturing tools by simultaneously feeding multiple wire and/or powder feedstocks at various rates to achieve spatial composition change. The wire-powder-directed energy deposition (WP-DED) technique is of particular interest for many functionally graded material applications by balancing the low raw materials cost of wire with the high resolution of powder. However, feeding wire and powder are inherently different processes since all extruded wire enters the melt pool, while much of the blown powder is scattered, which makes determining the composition of the build challenging. In this study, we devise a simple area-based measurement method for estimating the composition of WP-DED structures. WP-DED single beads are printed using 309L stainless steel wire and commercially pure Fe powder at five wire feed rates (0.5, 0.75, 1.00, 1.25, 1.50 mm/mm) and five powder feed rates (2, 4, 6, 8, 10 rpm). Characteristic defects including interface gaps and macrosegregation (lack of mixing) tendencies are examined. High powder feed rates (8, 10 rpm) result in interface gaps at all wire feed rates, but smooth deposition and complete mixing is achieved at low powder feed rates, particularly with lower wire feed rates as well. The area-based composition measurement method is within ±20% of energy dispersive x-ray spectroscopy measurements for all samples, showing its effectiveness as a rapid composition estimate for WP-DED materials development.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Substrate-Directed Dimensional and Phase Control of Peptide Assemblies on Two-Dimensional van der Waals Materials

Understanding and controlling biomolecular self-assembly on van der Waals (vdW) materials has the potential to advance hybrid bioelectronic devices by enabling precise tuning of the interface and modulation of the resulting electronic properties of the biomolecule-vdW heterostructure. However, how surface properties of vdW materials direct biomolecule assembly remains poorly understood. To fill this knowledge gap, we investigated the assembly of a peptide known to assemble into two-dimensional (2D) crystalline films on MoS 2 on three representative vdW surfaces: WS 2 , MoS 2 , and highly oriented pyrolytic graphite (HOPG). Using in situ atomic force microscopy (AFM), we find that assembly is substrate-dependent, resulting in multilayers on WS 2 , monolayers on MoS 2 , and multiple coexisting phases on HOPG. WS 2 exhibits a higher negative charge, strong long-range electrostatic interactions, and extensive hydration layering that may promote multilayer stacking. In contrast, MoS 2 has stronger short-range interactions with the peptides but much weaker long-range interactions and hydration structure, which may favor monolayer formation. Molecular dynamics simulations predict a corresponding switch from monolayer to multilayer aggregates of the adsorbed monomers, reflected in their relative mobilities. On hydrophobic HOPG, the peptides bind most strongly and remain as monomers with high surface mobility. The peptide dimers comprising the basic unit of the crystals are more compact on HOPG, which has a smaller lattice constant than WS 2 or MoS 2 , suggesting strain contributes to stabilizing multiple phases. Our results provide mechanistic insights into how surface charge and hydration structure, and the lattice structure of the substrates governs peptide assembly on vdW materials, offering a framework to rationally control the 2D peptide-vdW heterostructures.

Molecular dynamics simulations↗
DOE OSTI2026

Elevating SolTrace's Capabilities for the Next Generation of Concentrating Solar Analysis

SolTrace is an open-source Monte Carlo ray tracing software developed at NREL. SolTrace can characterize concentrating solar thermal (CST) collector optical performance and is CST technology agnostic. Shown in Fig. 1, SolTrace is a foundational tool in NREL's CST system and component modeling suite. SolTrace's generic surface elements can flexibly model novel collector and receiver designs to predict spatial and temporal flux distributions - critical to understand for CST component design, performance prediction, and system integration. Since its initial development, SolTrace has over 1,650 references on Google Scholar, over 9,800 downloads since 2017, and has served the CST research and development community as a benchmark of 3rd party verification. SolTrace provides users with many options for defining surface shape and boundaries. However, SolTrace provides limited documentation which can result in a steep learning curve for new users. Additionally, SolTrace lacks the computational performance required to evaluate optical performance of a CST system over the course of a year and/or iteratively over design parameters in a timely manner. To address this, we are working towards a new release of SolTrace that enables increased computational throughput by implementing ray tracing acceleration structures and enabling GPU parallelization. Additionally, we are working to improve SolTrace's usability, accessibility, and maintainability by (1) automating solar position time-dependent simulation processes, (2) creating general CST collector templates of grouped elements, (3) updating the user interface to better visualize model inputs and outputs, and (4) creating a user support network through forums, "how to" videos, and documentation.

14 SOLAR ENERGY↗
DOE OSTI2026

From bulk to surface: Structure and dynamics of amorphous alumina from deep potential molecular dynamics

Understanding the atomic-scale structure and dynamics of amorphous oxide surfaces is essential for interpreting their chemical reactivity, mechanical stability, and interfacial behavior, yet direct experimental characterization remains challenging. We employ Deep Potential (DP) molecular dynamics to generate large-scale, ab initio -quality models of amorphous Al 2 O 3 bulk glasses and melt-quenched free surfaces, enabling a quantitative analysis of both structure and relaxation dynamics with statistical confidence inaccessible to direct ab initio simulation. The trained DP model reproduces experimental liquid and glass structure, captures the cooling-rate dependence of the bulk glass transition, and corrects systematic biases in the polyhedral populations predicted by widely used classical force fields. At the free surface, mass density recovers to bulk values over ~10 Å, while local coordination requires a slightly wider subsurface region to fully converge. The outermost layer is oxygen-enriched, exhibits altered polyhedral connectivity with contracted Al–O bonds, and hosts a broad population of under-coordinated motifs (notably AlO 3 and OAl 2 ) whose abundances are governed by glass stability. These under-coordinated surface motifs exhibit distinct vibrational signatures and occur as locally paired Lewis acid and Brønsted base sites consistent with bond-valence compensation, yet remain spatially dispersed rather than aggregating into extended clusters. Despite this pronounced structural heterogeneity, surface relaxation and the glass-transition temperature remain comparable to their bulk counterparts, suggesting that the disordered surface is kinetically stable once formed. Together, these results establish a molecular-level picture of amorphous alumina surfaces and demonstrate the capability of machine-learned potentials to resolve structure–property relationships in disordered oxide interfaces.

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

Generalizable, fast, and accurate DeepQSPR with fastprop

Abstract Quantitative Structure–Property Relationship studies (QSPR), often referred to interchangeably as QSAR, seek to establish a mapping between molecular structure and an arbitrary target property. Historically this was done on a target-by-target basis with new descriptors being devised to specifically map to a given target. Today software packages exist that calculate thousands of these descriptors, enabling general modeling typically with classical and machine learning methods. Also present today are learned representation methods in which deep learning models generate a target-specific representation during training. The former requires less training data and offers improved speed and interpretability while the latter offers excellent generality, while the intersection of the two remains under-explored. This paper introduces , a software package and general Deep-QSPR framework that combines a cogent set of molecular descriptors with deep learning to achieve state-of-the-art performance on datasets ranging from tens to tens of thousands of molecules. provides both a user-friendly Command Line Interface and highly interoperable set of Python modules for the training and deployment of feedforward neural networks for property prediction. This approach yields improvements in speed and interpretability over existing methods while statistically equaling or exceeding their performance across most of the tested benchmarks. is designed with Research Software Engineering best practices and is free and open source, hosted at github.com/jacksonburns/fastprop.

Burns, Jackson W. (ORCID:0000000206579426)↗
DOE OSTI2024

Elucidating the Transition of 3D Morphological Evolution of Binary Alloys in Molten Salts with Metal Ion Additives

Molten salts serve as effective high-temperature heat transfer fluids and thermal storage media used in a wide range of energy generation and storage facilities, including concentrated solar power plants, molten salt reactors and high-temperature batteries. However, at the salt–metal interfaces, a complex interplay of charge-transfer reactions involving various metal ions, generated either as fission products or through corrosion of structural materials, takes place. Simultaneously, there is a mass transport of ions or atoms within the molten salt and the parent alloys. The precise physical and chemical mechanisms leading to the diverse morphological changes in these materials remain unclear. Here, to address this knowledge gap, this work employed a combination of synchrotron X-ray nanotomography and electron microscopy to study the morphological and chemical evolution of Ni-20Cr in molten KCl-MgCl 2 , while considering the influence of metal ions (Ni 2+ , Ce 3+ , and Eu 3+ ) and variations in salt composition. Our research suggests that the interplay between interfacial diffusivity and reactivity determines the morphological evolution. The summary of the associated mass transport and reaction processes presented in this work is a step forward toward achieving a fundamental comprehension of the interactions between molten salts and alloys. Overall, the findings offer valuable insights for predicting the diverse chemical and structural alterations experienced by alloys in molten salt environments, thus aiding in the development of protective strategies for future applications involving molten salts.

36 - MATERIALS SCIENCE↗
DOE OSTI2026

SAM Code Enhancements for Fission Product Tracking of Noble Gases and Metals in MSRs

This report documents fiscal year 2026 enhancements to the System Analysis Module (SAM) for modeling fission product transport in liquid-fueled molten salt reactors (MSRs). The work advances three principal areas: noble gas transport, noble metal deposition, and user interface improvements. The noble gas transport capability integrates drift-flux gas transport, Henry’s law two-film interphase mass transfer with pressure-based nucleation suppression, Knudsen-regime pore diffusion into porous graphite with a conjugate salt-graphite interface constraint, built-in material properties, five Sherwood-number mass transfer correlations including three derived from high-fidelity NekRS simulations, and xenon-135 reactivity feedback through SAM’s point-kinetics model. This work also presents a comprehensive verification test suite, including new analytically verified cases for pressure-dependent onset of interphase gas transfer in a stagnant vertical pipe, a postulated FLiBe-graphite Xe extraction permeator, a gravity riser with a fission-product source, and a descending pipe with gas redissolution driven by hydrostatic pressure. A machine learning framework for bubble rise velocity prediction in molten salt systems is developed and benchmarked on molten-salt and diverse aqueous bubble datasets. The best-performing fine-tuned transfer-learning networks achieve an 82% reduction in RMSE relative to the Clift correlation, and is implemented directly in SAM. The noble metal transport capability is developed, including a liquid-wall deposition model and a gas-surface flotation mechanism that transfers insoluble particles entrained by sparging gas to wetted structures. Verification tests and demonstration cases cover the surface deposition, flotation efflux, and flotation shedding. Finally, a new [SpeciesTransport] input structure replaces positional global vectors with selfcontained, order-independent, named species blocks, simplifies the specification of multiphase species and decay chains, and remains fully compatible with existing SAM input files. Together, these developments improve the physical fidelity, verification basis, and usability of SAM for system-level analyses of fissionproduct behavior in MSRs.

Mui, Travis (ORCID:0000000303736470)↗
DOE OSTI2024

Effect of Solvent Composition on Non-DLVO Forces and Oriented Attachment of Zinc Oxide Nanoparticles

Oriented attachment (OA) occurs when nanoparticles in solution align their crystallographic axes prior to colliding and subsequently fuse into single crystals. Traditional colloidal theories such as DLVO provide a framework for evaluating OA but fail to capture key particle interactions due to the atomistic details of both the crystal structure and the interfacial solution structure. Using zinc oxide as a model system, we investigated the effect of the solvent on short-ranged and long-ranged particle interactions and the resulting OA mechanism. In situ TEM imaging showed that ZnO nanocrystals in toluene undergo long-range attraction comparable to 1kT at separations of 10 nm and 3kT near particle contact. These observations were rationalized by considering non-DLVO interactions, namely dipole-dipole forces and torques between the polar ZnO nanocrystals. Langevin dynamics simulations showed stronger interactions in toluene compared to methanol solvents, consistent with the experimental results. Concurrently, we performed atomic force microscopy measurements using ZnO-coated probes for the short-ranged interaction. Our data provided valuable insights into another type of non-DLVO interaction, namely the repulsive solvation force. Specifically, the solvation force was stronger in water compared to ethanol and methanol, due to the stronger hydrogen bonding and denser packing of water molecules at the interface. In conclusion, our results highlight the importance of non-DLVO forces in a general framework for understanding and predicting particle aggregation and attachment.

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

Rigor and Reproducibility in Electrocatalysis: Best Practices for Operando Studies

Operando measurements have rapidly expanded the scope of electrocatalysis by enabling direct observation of catalytic interfaces under working conditions and by linking structural, compositional, and spectroscopic observables to activity and selectivity. However, the growth of operando methods has outpaced the adoption of broadly shared experimental standards, creating persistent challenges in reproducibility, interpretation, and comparison across laboratories and platforms. This perspective synthesizes discussions from the 2025 National Science Foundation Workshop on Rigor and Reproducibility in Electrocatalysis and outlines a practical framework for the rigorous use of operando measurements in electrocatalysis. We highlight three recurring needs: careful implementation of complex methods to avoid overinterpretation; recognition that (subtle) differences in reactor architecture, hydrodynamics, and electrical boundary conditions can alter apparent kinetics and selectivity; and transparent reporting standards that enable meaningful cross-comparison without constraining measurement-specific cell innovation. Focusing on widely used techniques (including X-ray and vibrational spectroscopies, mass spectrometry, and electron microscopy), we discuss technique-specific pitfalls, cross-validation strategies, and recurring platform-agnostic considerations such as mass transport, current distribution, temporal-resolution mismatches, and catalyst evolution. This Perspective aims to strengthen the mechanistic inference and improve the reproducibility, comparability, and predictive value of operando electrocatalysis research.

X-ray absorption spectroscopy↗
DOE OSTI2024

Equilipy: a python package for calculating phase equilibria

The CALPHAD (CALculation of PHAse Diagram) approach (Nigel Saunders & Miodownik, 1998) provides predictions for thermodynamically stable phases in multicomponent-multiphase materials across a wide range of temperatures. Consequently, the CALPHAD calculations became an essential tool in materials and process design (Luo, 2015). Such design tasks frequently require navigating a high-dimensional space due to multiple components involved in the system. This increasing complexity demands high-throughput CALPHAD calculations, especially in the rapidly evolving field of alloy design. In response to the need, we developed Equilipy an open-source Python package designed for calculating phase equilibria of multicomponent-multiphase systems. Equilipy is specifically tailored for high-throughput CALPHAD calculations, offering parallel computations across multiple processors and nodes with the given NPT input conditions namely elemental compositions (N), pressure (P), and temperature (T). Equilipy utilizes the program structure and Gibbs energy functions from the Fortran-based program, Thermochimica (Piro et al., 2013), with incorporating a new Gibbs energy minimization algorithm. This algorithm, originally developed by Capitani and Brown in 1987 (Capitani & Brown, 1987), has been revised and implemented to enhance the stability and performance of calculations. The Fortran codes are precompiled and interfaced with Python via F2PY, ensuring high computation speed. Benchmark tests shown in Figure 1 demonstrate that Equilipy’s computation speed is comparable to those of established commercial software, TC-Python and PanPython. This result highlights its efficiency and potential applications in various scientific and industrial fields.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Convergent Manufacturing of Large-Scale Components for Nuclear Applications, via Additive Manufacturing and Powder Metallurgy Hot Isostatic Pressing

Powder metallurgy (PM)–hot isostatic pressing (PM-HIP) has long been recognized as a powerful route for producing fully dense, near net shape metallic components. By consolidating powders under high temperature and pressure, HIP provides isotropic properties, uniform microstructures, and scalability to complex geometries that are vital for sectors such as aerospace, energy, and nuclear power. Yet despite these advantages, the technology has remained constrained by costly trial and error canister fabrication, limitations of conventional forging, and incomplete knowledge about how the canister design influences final part properties. Additive manufacturing (AM), by contrast, thrives on design freedom and geometric flexibility but struggles with speed, scalability, and cost when applied to very large structures. The research presented in this report investigated how a convergent manufacturing approach, combining AM with PM-HIP, can merge the strengths of both technologies, leveraging AM’s flexibility for canister design and HIP’s consolidation capability to deliver reliable, large, and complex parts. The work progressed through three case studies that built on one another in scale and complexity. Small cylindrical canisters fabricated by conventional methods, laser powder bed fusion, and directed energy deposition were filled with stainless steel powders and subjected to HIP. The resulting parts demonstrated near-full density and mechanical properties on par with wrought stainless steel, showing for the first time that AM canisters can be a direct substitute for conventional ones without sacrificing quality. The next step involved a medium-scale, noncentrosymmetric T-valve, which is an enclosed, multibranch geometry that tested the limits of AM + PM-HIP integration. The T-valve achieved predictable shrinkage and uniform densification, confirming feasibility for enclosed designs. However, this study also revealed oxide inclusions and interfacial challenges at the AM + HIP boundary, underscoring the critical importance of controlling interface chemistry and employing robust, in situ strategies, such as melt pool monitoring and thermal monitoring, coupled with nondestructive evaluation techniques such as x-ray computed tomography. Finally, the effort culminated in fabricating a large-scale impeller weighing nearly 2000 lb and spanning 5 ft in diameter. Produced via multirobot wire arc AM and hot isostatic pressed to near-full density, the impeller validated industrial-scale feasibility. Predictive models closely matched experimental shrinkage, tensile properties were spatially uniform across the component, and the AM + PM-HIP interface proved mechanically sound despite the presence of oxide-decorated prior particle boundaries. This large-scale demonstration is a major milestone, showing that hybrid AM + PM‑HIP can reliably deliver components at reactor-relevant scales. Collectively, these studies charted a logical pathway: small-scale work built scientific confidence, medium-scale work highlighted opportunities and challenges, and large-scale work proved industrial impact. The overarching conclusion of this report is that AM + PM-HIP should not be seen as a replacement for forging but as a complementary pathway that provides the US with flexibility, resilience, and new options for manufacturing nuclear-grade components. Looking ahead, several directions emerge as critical to sustaining progress. Predictive modeling must become faster, more accessible, and more accurate, with digital twins and machine learning reducing reliance on trial and error. Powders and alloys must be optimized for HIP, with improved cleanliness, reduced oxides, and tailored chemistries that enhance creep, fatigue, and irradiation resistance. Interfaces between AM and HIP regions must be better engineered through coatings, machining strategies, and surface treatments to mitigate oxide formation and ensure reliable bonding to explore opportunities for HIP of targeted compositional parts, as well as multimaterial HIP cladding applications. Monitoring and nondestructive evaluation need to expand, incorporating multimodal sensors, x-ray computed tomography, and real-time data integration through platforms such as Pelican. At the same time, the pathway to industrial adoption requires techno-economic analysis, machinability studies, and qualification frameworks aligned with industry and regulatory standards. Finally, workforce and academic engagement must be strengthened. Programs that train technicians and engineers for US Navy and US Department of Energy manufacturing challenges should be paired with academic partnerships to support fundamental research, with open sharing of non-export-controlled data to accelerate innovation and build the next generation of experts. In conclusion, this report demonstrates that hybrid AM + PM-HIP is scientifically viable and strategically important. By combining the design agility of AM with the consolidation strength of HIP and embedding modeling, monitoring, and workforce development, this approach provided a transformative new capability for US manufacturing. The path forward is clear: hybrid AM + PM-HIP is not just a promising research direction but is also potentially an industrially relevant pathway that can reshape how nuclear-grade components are designed, qualified, and deployed.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A multiphysics coupling framework for exascale simulation of fracture evolution in subsurface energy applications

Predicting the evolution of fractured media is challenging due to coupled thermal, hydrological, chemical and mechanical processes that occur over a broad range of spatial scales, from the microscopic pore scale to field scale. We present a software framework and scientific workflow that couples the pore scale flow and reactive transport simulator Chombo-Crunch with the field scale geomechanics solver in GEOS to simulate fracture evolution in subsurface fluid-rock systems. This new multiphysics coupling capability comprises several novel features. An HDF5 data schema for coupling fracture positions between the two codes is employed and leverages the coarse resolution of the GEOS mechanics solver which limits the size of data coupled, and is, thus, not taxed by data resulting from the high resolution pore scale Chombo-Crunch solver. The coupling framework requires tracking of both before and after coarse nodal positions in GEOS as well as the resolved embedded boundary in Chombo-Crunch. We accomplished this by developing an approach to geometry generation that tracks the fracture interface between the two different methodologies. The GEOS quadrilateral mesh is converted to triangles which are organized into bins and an accessible tree structure; the nodes are then mapped to the Chombo representation using a continuous signed distance function that determines locations inside, on and outside of the fracture boundary. The GEOS positions are retained in memory on the Chombo-Crunch side of the coupling. The time stepping cadence for coupled multiphysics processes of flow, transport, reactions and mechanics is stable and demonstrates temporal reach to experimental time scales. The approach is validated by demonstration of 9 days of simulated time of a core flood experiment with fracture aperture evolution due to invasion of carbonated brine in wellbore-cement and sandstone. We also demonstrate usage of exascale computing resources by simulating a high resolution version of the validation problem on OLCF Frontier.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Exploring Domain-Wall Pinning in Ferroelectrics via Automated High-Throughput Atomic Force Microscopy

Domain-wall dynamics in ferroelectric materials are strongly position-dependent, since each polar interface is locked into a unique local microstructure. This necessitates spatially resolved studies of wall pinning using scanning-probe microscopy techniques. The pinning centers and pre-existing domain walls are usually sparse within the image plane, precluding the use of dense hyperspectral imaging modes and requiring time-consuming human experimentation. Here, a large-area epitaxial PbTiO 3 film on cubic KTaO 3 was investigated to quantify the electric-field-driven dynamics of the polar–strain domain structures using ML-controlled automated piezoresponse force microscopy. Analysis of 1500 switching events reveals that domain-wall displacement depends not only on field parameters but also on the local ferroelectric–ferroelastic configuration. For example, twin boundaries in polydomains regions, like a 1 – /c+ ∥ a 2 – /c – , stay pinned up to a certain level of bias magnitude and change only marginally as the bias increases from 20 to 30 V, whereas single-variant boundaries, like the a 2 + /c + ∥ a 2 – /c – stack, are already activated at 20 V. These statistics on the possible ferroelectric and ferroelastic wall orientations, together with the automated high-throughput AFM workflow, can be distilled into a predictive map that links domain configurations to pulse parameters. Here, this microstructure-specific rule set forms the foundation for the design of ferroelectric memories.

automated scanning probe microscopy↗
DOE OSTI2024

Myna

The additive manufacturing (AM) community has been developing digital factory tools over the past decade to better leverage the multi-modal process data coming out of the advanced manufacturing process. As a result, numerous databases of additive manufacturing process data exist in the literature and in the archival storage of disparate research groups. While some efforts have been made to create a standard ontology for storing and sharing AM data, in practice a variety of data structures are used to store AM build data, even within a single institution. This causes many problems for maintainability and extensibility when attempting to integrate computational modeling tools with experimental data to either validate models or to provide further insight into results and trends. Myna is a Python-based framework that aims to decrease the effort needed to connect individual computational models to the variety of AM process data that exist in different research groups and institutions. This type of software is sometimes referred to as "middleware" or “glueware,” in that it connects disparate databases and applications into a single computational ecosystem. Instead of maintaining unique interfaces between each application and each database, developers can create a single interface from each application to Myna and thereby gain access to the implemented database connections. Similarly, developing a database connection in Myna provides access to the developed simulation applications. This framework greatly simplifies the maintainability of model applications that rely on experimental data. Using external simulation tools, users will also be able to run pre-configured workflows using the built-in workflow manager. Several examples of input files are provided with Myna for different workflows, including melt pool geometry predictions and detailed melt pool and solidification microstructure predictions.

Knapp, GerryL. [Oak Ridge National Laboratory (ORN↗
DOE OSTI2024

The Radiation-Induced Fate of Fission Product Iodine in Molten Salts

Understand and Predict Radiation-Induced Iodine Speciation, Chemistry, and Transport in High-Temperature Molten Salts Award Number: DE-AC07-05ID1451 Gregory P. Holmbeck (Gregory.Holmbeck@inl.gov), Center for Radiation Chemistry Research, Idaho National Laboratory, 1955 N. Fremont Avenue Idaho Falls, ID, 83415, USA. Project Scope The goal of this Chemical and Materials Sciences to Advance Clean Energy Technologies and Low-Carbon Manufacturing project is to understand and predict the radiation-induced speciation, chemistry, and transport of fission product iodine in the triumvirate extremes of high-temperature, ionizing radiation, and corrosive molten salts. This missing fundamental information is critical for the accelerated development and deployment of safe, clean nuclear energy based on molten salt reactor (MSR) and pyrochemical reprocessing technologies. The central hypothesis driving this research is, the radiation-induced conversion of iodide will yield an extensive suite of transient and steady-state iodine radiolysis products that will alter the bulk chemical and physical properties of the irradiated molten salt system—the speciation, distribution, and chemical transport of which will be dictated by the composition and the availability of multivalent metal cations and metal alloy interfaces. To test this hypothesis, this project initiated three synergistic Research Objectives: (1) determine how the inclusion of iodine/iodide influences the chemical and physical properties of complex molten salt mixtures; (2) elucidate the speciation and fundamental chemical behavior of transient and steady-state iodine/iodide species formed by the irradiation of molten and solid salt mixtures; and (3) understand the influence of interfacial processes on determining the final disposition of iodine in high temperature molten salts, notably the structure and chemical speciation of iodine at metal-salt interfaces.

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

Discoveries in Blast-Driven Turbulence of Astrophysical Relevance

The fluid mixing caused by variable-density instabilities is important in a wide variety of scenarios from ocean mixing and astrophysical phenomena to nuclear fusion techniques and atomic weapons. This thesis explores the mixing resulting from a specific instability known as the Blast Driven Instability (BDI). This work investigates the variable density mixing in an explosively driven environment due to the fluid instabilities at the material interfaces. Specifically, diverging Richtmyer-Meshkov (impulsive-acceleration environment) and Rayleigh-Taylor (variable-acceleration environment) instabilities (present in supernova and inertial confinement fusion) are studied using advanced high-speed diagnostics in carefully designed laboratory experiments. The BDI morphology is presented through a time development of Mie scattering images, and steps through the parameter space (varying density ratio and driver speed), highlighting the development of the structures that form during mixing. A scaling criterion is used to relate the two systems of vastly different spatiotemporal scales. Velocity fields in the BDI have been captured for the first time using the high temporal resolution PIV technique. Subsequent analysis of the dynamics of the instability from the velocity fields illustrates the distribution of kinetic energy, the transition to turbulence, and the characteristic growth of the instability are discussed. This study furthers understanding of how blast-driven instability pertains to supernova and inertial confinement fusion science. The morphology of the BDI has been characterized for the first time. This work steps through the parameter space covered in Mie scattering experiments, and how the different parameters contribute to development of structures and mixing. It also examines a scaling of the Atwood number for expanding predictive capabilities to other experimental conditions and simulations. The first collection of velocity fields acquired for the BDI are recorded, and subsequent analysis evaluating the distribution of kinetic energy throughout space and time for two density ratios from the overall parameter space, as well as the transition to turbulence, estimated from a Reynolds number calculated based on momentum mixing are all presented. This information is useful in advancing the development of models to predict physics of high energy density applications where experiments are not always readily available. This research has successfully demonstrated understanding for the time criteria defining regimes where the shock driven (Richtmyer-Meshkov instability) and the buoyancy driven (Rayleigh-Taylor instability) dominates through a parametric study of density variation (Atwood number) and driver speed (Mach number). All this furthers understanding of how the BDI pertains to SN and ICF.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

A Finite Element Method for Compressible and Turbulent Multiphase Flow Instabilities with Heat Transfer

We present a new finite element framework for modeling compressible, turbulent multiphase flows with heat transfer. For two-fluid systems with a free surface, the Volume of Fluid (VOF) method is implemented without the need for interface reconstruction, while turbulence is resolved using a dynamic Vreman large eddy simulation (LES) model. Unlike most two-phase VOF studies, which neglect heat transfer, the present approach incorporates energy transport equations within the VOF formulation to account for heat exchange, an effect particularly important in turbulent flows. Conjugate heat transfer is often challenging in finite volume methods, which require explicit specification of heat fluxes at the solid–fluid interface, limiting accuracy and predictive capability. By contrast, the finite element formulation does not require heat flux inputs, allowing more accurate and robust simulation of heat transfer between solids and fluids. The method is demonstrated through three representative cases. First, a two-fluid instability with a single-mode perturbation is simulated and validated against analytical growth rates. Second, conjugate heat transfer is examined in a high-temperature flow over a cold metal cylinder, with validation performed both quantitatively—via pressure coefficient comparisons with experimental data—and qualitatively using vector field topology. Finally, compressible spray injection and breakup are modeled, demonstrating the ability of the framework to capture interfacial dynamics and atomization under turbulent, high-speed conditions. In the compressible spray injection and breakup case, the results indicate that the finite element formulation achieved higher predictive accuracy and robustness than the finite-volume method. With the same mesh resolution, the FEM reduced the root mean square error (RMSE) and mean absolute percentage error (MAPE) from 6.96 mm and 26.0% (for the FVM) to 4.85 mm and 12.7%, respectively, demonstrating improved accuracy and robustness in capturing interfacial dynamics and heat transfer. The study also introduced vector field topology to visualize and interpret coherent flow structures and instabilities, offering insights beyond conventional scalar-field analyses.

97 MATHEMATICS AND COMPUTING↗