SEARCH · Engineering Papers

Results for “Interface Problems”

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 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Ionic Conductive and Highly-Stable Interface for Alkali Metal Anodes

Alkali metals are regarded as the most promising candidates for advanced anode for the next-generation batteries due to their high specific capacity, low electrochemical potential, and lightweight. However, critical problems of the alkali metal anodes, especially dendrite formation and interface stabilization, remain challenging to overcome. The solid electrolyte interphase (SEI) is a key factor affecting Li and Na deposition behavior and electrochemical performances. Herein, a facile and universal approach is successfully developed to fabricate ionic conductive interfaces for Li and Na metal anodes by modified atomic layer deposition (ALD). In this process, the Li metal (or Na metal) plays the role of Li (or Na) source without any additional Li (or Na) precursor during ALD. Moreover, the key questions about the influence of ALD deposition temperature on the compositions and structure of the coatings are addressed. The optimized ionic conductive coatings have significantly improved the electrochemical performances. In addition, the electrochemical phase-field model is performed to prove that the ionic conductive coating is very effective in promoting uniform electrodeposition. Further, this approach is universal and can be potentially applied to other different metal anodes. At the same time, it can be extended to other types of coatings or other deposition techniques.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Partitioning of Complex Fluids at Mineral Surfaces

This report summarizes the results obtained during the LDRD project entitled "Partitioning of Complex Fluids at Mineral Interfaces." This research addressed fundamental aspects of such interfaces, which are relevant to energy-water applications in the subsurface, including fossil energy extraction and carbon sequestration. This project directly addresses the problem of selectivity of complex fluid components at mineral-fluid interfaces, where complex fluids are defined as a mixture of hydrophobic and hydrophilic components: e.g., water, aqueous ions, polar/nonpolar organic compounds. Specifically, this project investigates how adsorption selectivity varies with surface properties and fluid composition. Both experimental and molecular modeling techniques were used to better understand trends in surface wettability on mineral surfaces. The experimental techniques spanned the macroscale (contact angle measurements) to the nanoscale (cryogenic electronic microscopy and vibrational spectroscopy). We focused on an anionic surfactant and a well-characterized mineral phase representative of clay phases present in oil- and gas-producing shale deposits. Collectively, the results consistently demonstrate that the presence of surfactant in the aqueous fluid significantly affects the mineral-fluid interfacial structure. Experimental and molecular modeling results reveal details of the surfactant structure at the interface, and how this structure varies with surfactant coverage and fluid composition.

04 OIL SHALES AND TAR SANDS↗

DOE OSTI2023

Interface flux recovery framework for constructing partitioned heterogeneous time-integration methods

A common approach for the development of partitioned schemes employing different time integrators on different subdomains is to lag the coupling terms in time. This can lead to accuracy issues, especially in multistage methods. Here, in this article, we present a novel framework for partitioned heterogeneous time-integration methods, which allows the coupling of arbitrary multistage and multistep methods without reducing their order of accuracy. At the core of our approach are accurate estimates of the interface flux obtained from the Schur complement of an auxiliary monolithic system. We use these estimates to construct a polynomial-in-time approximation of the interface flux over the current time coupling window. This approximation provides the interface boundary conditions necessary to decouple the subdomain problems at any point within the coupling window. In so doing our framework enables a flexible choice of time-integrators for the individual subproblems without compromising the time-accuracy at the coupled problem level. This feature is the main distinction between our framework and other approaches. To demonstrate the framework, we construct a family of partitioned heterogeneous time-integration methods, combining multistage and multistep methods, for a simplified tracer transport component of the coupled air-sea system in Earth system models. We report numerical tests evaluating accuracy and flux conservation for different pairs of time-integrators from the explicit Runge-Kutta and Adams-Moulton families.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

A mineral precipitation model based on the volume of fluid method

A novel volume of fluid method is presented for mineral precipitation coupled with fluid flow and reactive transport. The approach describes the fluid-solid interface as a smooth transitional region, which is designed to provide the same precipitation rate and viscous drag force as a sharp interface. Specifically, the governing equation of mineral precipitation is discretized by an upwind scheme, and a rigorous effective viscosity model is derived around the interface. The model is validated against analytical solutions for mineral precipitation in channel and ring-shaped structures. It also compares well with interface tracking simulations of advection-diffusion-reaction problems. Here, the methodology is finally employed to model mineral precipitation in fracture networks, which is challenging due to the low porosity and complex geometry. Compared to other approaches, the proposed model has a concise algorithm and contains no free parameters. In the modeling, only the pore space requires meshing, which improves the computational effciency especially for low-porosity media.

58 GEOSCIENCES↗

DOE OSTI2025

MTUQ: a framework for estimating moment tensors, point forces, and their uncertainties

SUMMARY We introduce MTUQ, an open-source Python package for seismic source estimation and uncertainty quantification, emphasizing flexibility and operational scalability. MTUQ provides MPI-parallelized grid search and global optimization capabilities, compatibility with 1-D and 3-D Green’s function database formats, customizable data processing, C-accelerated waveform and first-motion polarity misfit functions, and utilities for plotting seismic waveforms and visualizing misfit and likelihood surfaces. Applicability to a range of full- and constrained-moment tensor, point force, and centroid inversion problems is possible via a documented application programming interface, accompanied by example scripts and integration tests. We demonstrate the software using three different types of seismic events: (1) a 2009 intraslab earthquake near Anchorage, Alaska; (2) an episode of the 2021 Barry Arm landslide in Alaska; and (3) the 2017 Democratic People’s Republic of Korea underground nuclear test. With these events, we illustrate the well-known complementary character of body waves, surface waves, and polarities for constraining source parameters. We also convey the distinct misfit patterns that arise from each individual data type, the importance of uncertainty quantification for detecting multimodal or otherwise poorly constrained solutions, and the software’s flexible, modular design.

58 GEOSCIENCES↗

DOE OSTI2025

Gate-tunable enhancement of supercurrent in hybrid planar Josephson junctions

Planar Josephson junctions (JJs) have emerged as a promising platform for the realization of topological superconductivity and Majorana zero modes. To obtain robust quasi one-dimensional (1D) topological superconducting states using planar JJs, limiting the number of 1D Andreev bound states’ subbands that can be present, and increasing the size of the topological superconducting gap are two fundamental challenges. It has been suggested that both problems can be addressed by properly designing the interfaces between the JJ’s normal region and the superconducting leads. We fabricated Josephson junctions with periodic hole structures on the superconducting contact leads on InAs heterostructures with epitaxial superconducting Al. By depleting the chemical potential inside the holes region with a top gate, we observed an enhancement of the supercurrent across the junction. The theoretical analysis shows that the enhancement of the JJ’s critical current is achieved when the depletion of the holes is such to optimize the matching of quasiparticles’ wave function at the normal/superconductor interface. Furthermore, these results show how the combination of carefully designed patterns for the Al coverage, and external gates, can be successfully used to tune the density and wave functions’ profiles in the normal region of the JJ, and therefore open an avenue to tune some of the critical properties, such as number of subbands and size of the topological gap, that must be optimized to obtain robust quasi-1D superconducting states supporting Majorana bound states.

Bogoliubov-de Gennes equations↗

DOE OSTI2020

Stress Waves in Solids [Slides]

The outline for SSDDS Lecture 3 includes: stress and strain, stress waves, waves at interfaces and boundaries, formation of shock waves, and problems.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Development of Modified Perturbation Solutions to the One-Phase Stefan Problems With a Convective Boundary

The classical Stefan problem is used to track the moving solid-liquid interface during the freezing process. Perturbation theory has often been applied to find an approximate analytical solution due to the nonlinearity of the moving interface. However, the Stefan number (i.e., the sensible over latent heat) must be small and usually less than 0.01 to assume the perturbation expansion, which in turn limits the thermal engineering applications. In this study, a modified perturbation solution is developed by adding a correction term after the leading-order solution to be valid for a much wider range of Stefan numbers (i.e., 0.01 less than or equal to Ste less than or equal to 1). Specifically, a one-phase Stefan problem is first formulated subjected to a convective boundary in the Cartesian, cylindrical, and spherical coordinate systems for inward solidification. The leading-order solution is calculated based on the regular perturbation theory, while the correction term is obtained using the Monte-Carlo method and a multi-variant regression. Results show that the correction term has a linear relationship with the Stefan number and is not significantly influenced by the Biot number. The proposed modified perturbation solution can accurately and rapidly predict the nonlinear moving interface motion for the freezing process.

DIRECT ENERGY CONVERSION,MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Challenges and Prospects of All-Solid-State Electrodes for Solid-State Lithium Batteries

In the development of all-solid-state lithium batteries (ASSLB), progress is made with solid-state electrolytes; however, challenges regarding compatibility and stability still exist with solid electrodes. These issues result in a low battery capacity and short cycle life, which limit the commercial application of ASSLBs. Here, this review summarizes the recent research progress on solid-state electrodes in ASSLBs including the solid–solid interface phenomena such as the interface between electrode materials and electrolytes. The mechanical stability problems in solid electrodes, including fracture, brittleness, and deformation of electrode materials, are also discussed, and corresponding methods to measure the solid electrode stress are provided. In addition, strategies for mitigating stress-related issues are examined. Finally, the fabrication process of solid electrodes is introduced and their future developments, including the exploration of new electrode materials and the design of more intelligent electrode structures, are proposed.

25 ENERGY STORAGE↗

DOE OSTI2025

Data for: Subsurface Interface Structure Controlling Local Electronic Properties of Epitaxial Graphene on SiC(0001)

Recently realized high-mobility semiconducting epitaxial graphene on silicon carbide, provided an important step towards integration of the graphene-based system into active components in post-silicon micro- and nano-electronics. However, the exact atomic-scale structure and the complex bonding configurations of the first epitaxial graphene carbon layer remain an open problem. Our recent report has shed new light on understanding this interface, where the external transverse electric field-dependent dynamic switching behavior of the Cbuffer-SiC bonds was observed. Here, using scanning tunneling microscopy and spectroscopy (STM and STS), we present the direct evidence of silicon (Si) vacancies at the interface and provide their distribution at the topmost reconstructed SiC(0001) layer. Experimental STM and density functional theory modeling data were used in the preparation of figures in a published article in the Journal of Physical Chemistry Letters. Files related to the figures and supplementary materials in the article are present in this dataset in .txt format.

Condensed matter imaging↗

DOE OSTI2025

Data for: Subsurface Interface Structure Controlling Local Electronic Properties of Epitaxial Graphene on SiC(0001)

Condensed matter imaging↗

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

A time-parallel multiple-shooting method for large-scale quantum optimal control

Quantum optimal control plays a crucial role in quantum computing by providing the interface between compiler and hardware. Solving the optimal control problem is particularly challenging for multi-qubit gates, due to the exponential growth in computational complexity with the system's dimensionality and the deterioration of optimization convergence. To ameliorate the computational complexity of time-integration, this paper introduces a multiple-shooting approach in which the time domain is divided into multiple windows and the intermediate states at window boundaries are treated as additional optimization variables. Further, this enables parallel computation of state evolution across time-windows, significantly accelerating objective function and gradient evaluations. Since the initial state matrix in each window is only guaranteed to be unitary upon convergence of the optimization algorithm, the conventional gate trace infidelity is replaced by a generalized infidelity that is convex for non-unitary state matrices. Continuity of the state across window boundaries is enforced by equality constraints. A quadratic penalty optimization method is used to solve the constrained optimal control problem, and an efficient adjoint technique is employed to calculate the gradients in each iteration. We demonstrate the effectiveness of the proposed method through numerical experiments on quantum Fourier transform gates in systems with 2, 3, and 4 qubits, noting a speedup of 80x for evaluating the gradient in the 4-qubit case, highlighting the method's potential for optimizing control pulses in multi-qubit quantum systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

On the Abuse and Detection of Polyglot Files

A polyglot is a file that is valid in two or more formats. Polyglot files pose a problem for file-upload and generative AI web interfaces that rely on format identification to determine how to securely handle incoming files. In this work we found that existing file-format and embedded-file detection tools, even those developed specifically for polyglot files, fail to reliably detect polyglot files used in the wild. To address this issue, we studied the use of polyglot files by malicious actors in the wild, finding 30 polyglot samples and 15 attack chains that leveraged polyglot files. Using knowledge from our survey of polyglot usage in the wild---the first of its kind---we created a novel data set based on adversary techniques. We then trained a machine learning detection solution, PolyConv, using this data set. PolyConv achieves a precision-recall area-under-curve score of 0.999 with an F1 score of 99.20% for polyglot detection and 99.47% for file-format identification, significantly outperforming all other tools tested. We developed a content disarmament and reconstruction tool, ImSan, that successfully sanitized 100% of the tested image-based polyglots, which were the most common type found via the survey. Our work provides concrete tools and suggestions to enable defenders to better defend themselves against polyglot files, as well as directions for future work to create more robust file specifications and methods of disarmament.

Oesch, T [ORNL] (ORCID:0000000269091022)↗

DOE OSTI2020

Innovations in optimization and control of accelerators using methods of differential geometry and genetic algorithms (Final Report)

Online tuning of particle accelerators is necessary in order to achieve optimal machine performance. However, it is also a major challenge due to the large parameter space which must be searched and the fact that many of the desir- able objectives compete with one another, and so will not reach their optimal values simultaneously. In order to mitigate these issues, we have explored using dimension-reduction techniques to reduce the size of the parameter space which must be searched and multi-objective genetic algorithms to obtain the sets of tuning parameters which provide pareto-optimal values for the objectives. These methods have enabled us to obtain improved values of the vertical emittance at the Cornell Electron Storage Ring (CESR), and to do so with greater control of orbit errors. Algorithmic tuning is a multidisciplinary endeavour, requiring expertise in beam dynamics, diagnostics, control systems and computer science, and thus a key practical problem is to formulate a common language in which experts with different specialties can communicate. We developed a solution to this problem in the form of a generic accelerator software interface that allows for rapid prototyping of optimization and control algorithms. Our interface is built on the Experimental Physics and Industrial Control Systems (EPICS) and makes possible testing control code in simulation before deployment on real accelerators, as well as deployment of third-party optimization code.

43 PARTICLE ACCELERATORS↗

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 OSTI2024

Options for Achieving Cost Reduction in Advanced Reactors through Open Architecture

A key contributor to high capital costs and schedule overruns for new nuclear power plants is lack of standardization, driven by site-specific customization and construction of multiple designs by competing vendors rather than commitment to a single standardized program. While advanced reactor vendors typically individually target repeat construction of standardized units, the many competing designs could exacerbate the problem. “Open Architecture”, the open specification of requirements and interfaces for structures, systems and components (SSCs), has been proposed as a means of promoting standardization, by facilitating existing non-nuclear suppliers to enter the industry and/or allowing SSCs to be configured for more than one reactor within the same technology type. Contracting mechanisms that facilitate information sharing and alignment of incentives between stakeholders may complement such an approach. A preliminary scheme is presented for selection of SSCs for which such strategies could be adopted, based on a vendor make/buy decision model and stakeholder interviews. SSCs are categorized according to number of suppliers and their contribution to the reactor’s competitive edge. SSCs with many potential suppliers and a high contribution to competitive edge may be attractive for widening the supply chain via open specification of system requirements and interfaces, e.g., SSCs in the power island. SSCs with few suppliers and low contribution to competitive edge may be potential avenues for common system specification between vendors, e.g., some of the auxiliary SSCs. Potential cost reductions from such strategies will depend upon the size of the build program and the reactor type.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

Mass-conserving implicit–explicit methods for coupled compressible Navier–Stokes equations

Earth system models are composed of coupled components that separately model systems such as the global atmosphere, ocean, and land surface. While these components are well developed, coupling them in a single system can be a significant challenge. Computational efficiency, accuracy, and stability are principal concerns. In this study we focus on these issues. In particular, implicit–explicit (IMEX) tight and loose coupling strategies are explored for handling different time scales. For a simplified model for the air–sea interaction problem, we consider coupled compressible Navier–Stokes equations with an interface condition. Under the rigid-lid assumption, horizontal momentum and heat flux are exchanged through the interface. Several numerical experiments are presented to demonstrate the stability of the coupling schemes. Here, we show both numerically and theoretically that our IMEX coupling methods are mass conservative for a coupled compressible Navier–Stokes system with the rigid-lid condition

42 ENGINEERING↗