Engineering Papers⌕ Search

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.

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

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Multidimensional Tests of a Finite-Volume Solver for MHD With a Real-Gas Equation of State

This work considers two algorithms of a finite-volume solver for the MHD equations with a real-gas equation of state (EOS). Both algorithms use a multistate form of Harten-Lax-Van Leer approximate Riemann solver as formulated for MHD discontinuities. This solver is modified to use the generalized sound speed from the real-gas EOS. Two methods are tested: EOS evaluation at cell centers and at flux interfaces where the former is more computationally efficient. A battery of 1D and 2D tests are employed: convergence of 1D and 2D linearized waves, shock tube Riemann problems, a 2D nonlinear circularly polarized Alfvén wave, and a 2D magneto-Rayleigh-Taylor instability test. The cell-centered EOS evaluation algorithm produces unresolvable thermodynamic inconsistencies in the intermediate states leading to spurious solutions while the flux-interface EOS evaluation algorithm robustly produces the correct solution. The linearized wave tests show this inconsistency is associated with the magnetosonic waves and the magneto-Rayleigh-Taylor instability test demonstrates simulation findings where the spurious solution leads to an unphysical simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Hybrid large-eddy simulation with adaptive downstream anisotropic eddy viscosity model

Turbulence pervades most flows of engineering interest, and its prediction remains a challenge on both accuracy and cost. One promising predictive approach that reduces cost combines large-eddy simulation (LES) with simulation based on Reynolds averaged Navier–Stokes equations (RANS). This study presents a method to overcome stability and accuracy issues associated with these hybrid LES-RANS methods. The method developed involves extracting the Reynolds stresses from the upstream LES solution and then using that information to convert the downstream RANS model from a scalar eddy-viscosity model to an anisotropic nonlinear eddy-viscosity model. The method developed differs from the downstream tensorial eddy-viscosity model by being independent of the coordinate system. The method developed was evaluated by computing film cooling of a flat plate with the coolant injected through one row of circular holes. Results obtained show instabilities at the LES-to-RANS interface to be eliminated. Also, the method developed yielded solutions that compare reasonably well with those from LES, even though a significant portion of the flow is computed by the adapted anisotropic RANS model instead of LES, which significantly reduced the number of grid points and computational time needed. Since modification of the downstream RANS model is based on information extracted from the upstream LES solution, the method developed is adaptive to the problem being studied.

42 ENGINEERING↗
DOE OSTI2023

$\mathrm{RADICAL}$-Pilot and $\mathrm{PMIx}$/$\mathrm{PRRTE}$: Executing Heterogeneous Workloads at Large Scale on Partitioned $\mathrm{HPC}$ Resources

Execution of heterogeneous workflows on high-performance computing (HPC) platforms present unprecedented resource management and execution coordination challenges for runtime systems. Task heterogeneity increases the complexity of resource and execution management, limiting the scalability and efficiency of workflow execution. Re-source partitioning and distribution of tasks execution over portioned re-sources promises to address those problems but we lack an experimental evaluation of its performance at scale. Here this paper provides a performance evaluation of the Process Management Interface for Exascale (PMIx) and its reference implementation PRRTE on the leadership-class HPC plat-form Summit, when integrated into a pilot-based runtime system called RADICAL-Pilot. We partition resources across multiple PRRTE Distributed Virtual Machine (DVM) environments, responsible for launching tasks via the PMIx interface. We experimentally measure the work-load execution performance in terms of task scheduling/launching rate and distribution of DVM task placement times, DVM startup and termination overheads on the Summit leadership-class HPC platform. Integrated solution with PMIx/PRRTE enables using an abstracted, standardized set of interfaces for orchestrating the launch process, dynamic process management and monitoring capabilities. It extends scaling capabilities allowing to overcome a limitation of other launching mechanisms (e.g., JSM/LSF). Explored different DVM setup configurations provide insights on DVM performance and a layout to leverage it. Our experimental results show that heterogeneous workload of 65,500 tasks on 2048 nodes, and partitioned across 32 DVMs, runs steady with resource utilization not lower than 52%. While having less concurrently executed tasks resource utilization is able to reach up to 85%, based on results of heterogeneous workload of 8200 tasks on 256 nodes and 2 DVMs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Linear stability analysis of two fluid columns of different densities and viscosities in a gravity field

The linear stability of a vertical interface separating two miscible fluid columns of different densities and viscosities under the influence of gravity is investigated. This flow possesses a time-dependent reference state (each column accelerates at different rates owing to their different densities) and the interface thickness grows as the square root of time (by diffusion). Numerical integration of the linear initial-value problem is carried out and discussed in detail as a function of vertical and spanwise wavenumbers and the flow parameters. Adjoint-based optimization is performed in order to determine initial conditions that lead to maximum growth of disturbances in finite time. Results indicate that the rate of growth of the perturbation energy at small wavenumbers (less affected by viscosity initially) is dominated by two-dimensional modes (no spanwise variation). Substantial transient growth is observed at higher wave modes initially, followed by asymptotic decay of the perturbations at large time. Sensitivity of perturbation growth with respect to initial time, density and viscosity ratios is investigated. This work is complementary to previous inviscid analysis of this configuration, which showed that the interface was unconditionally unstable at all wave modes, even in the presence of surface tension, and that instability grew as the exponential of time squared.

Mechanics↗
DOE OSTI2023

MyCrunchGPT: A LLM Assisted Framework for Scientific Machine Learning

Scientific machine learning (SciML) has advanced recently across many different areas in computational science and engineering. Here, the objective is to integrate data and physics seamlessly without the need of employing elaborate and computationally taxing data assimilation schemes. However, preprocessing, problem formulation, code generation, postprocessing, and analysis are still time- consuming and may prevent SciML from wide applicability in industrial applications and in digital twin frameworks. Here, we integrate the various stages of SciML under the umbrella of ChatGPT, to formulate MyCrunchGPT, which plays the role of a conductor orchestrating the entire workflow of SciML based on simple prompts by the user. Specifically, we present two examples that demonstrate the potential use of MyCrunchGPT in optimizing airfoils in aerodynamics, and in obtaining flow fields in various geometries in interactive mode, with emphasis on the validation stage. To demonstrate the flow of the MyCrunchGPT, and create an infrastructure that can facilitate a broader vision, we built a web app based guided user interface, that includes options for a comprehensive summary report. The overall objective is to extend MyCrunchGPT to handle diverse problems in computational mechanics, design, optimization and controls, and general scientific computing tasks involved in SciML, hence using it as a research assistant tool but also as an educational tool. While here the examples focus on fluid mechanics, future versions will target solid mechanics and materials science, geophysics, systems biology, and bioinformatics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A Domain-Decomposed A-ϕ Formulation Based on Lagrange Multipliers for Low-Frequency Problems

A domain-decomposed A-ϕ formulation based on Lagrange multipliers is proposed to simulate low-frequency elec- tromagnetic problems. This method partitions the computational domain into smaller subdomains, allowing each subdomain to be independently formulated using Lagrange multipliers as Dirichlet boundary conditions, while ensuring continuity of the fields across the interfaces. A mixed finite element method, utilizing both vector and scalar basis functions, is employed to discretize the formulation, resulting in a global system to be solved. The proposed method is validated using TEAM Problem 7 at 50 Hz, demonstrating its effectiveness in handling complex geometries and addressing the low-frequency breakdown issues commonly encountered in traditional finite element methods.

Hossain, Amzad↗
DOE OSTI2025

Corrosion Resistance of an AlCeMg/Stainless-Steel Reactive Bond

A major issue for metal components in many industries is corrosion as it can substantially reduce their lifetime. This issue is especially problematic for materials used in heat exchanger applications. Al–Ce–Mg alloys, which exhibit corrosion resistance and can reactively bond with other metals, may be a viable solution to this problem. This investigation studied the corrosion behavior of Al–2Ce–6Mg (atomic percent)/stainless-steel (SS) reactive bond interfaces after full immersion in nitric, sulfuric, formic, and mixed acid for 267 h. This particular Al–Ce–Mg alloy was chosen due to its good castability. The results of scanning electron microscope characterization showed that reactive bond formations repeatedly occurred throughout the length of the casting in the as-cut samples and that these formations maintained a secure bonding between the alloy and the stainless-steel tubes. Transmission electron microscopy results showed that there was a clear compositional and microstructural transition across the reactive bond. The results of the immersion tests indicated that the nitric, sulfuric, and the mixed acid did not have an observably negative effect on the reactive bond structure. As for the sample exposed to formic acid only, noticable changes were seen in both the microstructural appearance and the elemental profile across the bond, suggesting that oxide formation occurred.

Brechtl, Jamieson [ORNL] (ORCID:0000000217394283)↗
DOE OSTI2021

Improved MPACT energy deposition and explicit heat generation coupling with CTF

The CASL neutronics code MPACT assumes local energy deposition with equilibrium delayed energy for steady-state and transient calculations. These approximations limit the accurate representation of the heat generation in space and its variations over time, both of which are essential for power distribution and thermal–hydraulic coupling analyses. This paper presents an explicit energy deposition model considering the neutron and gamma heating of all regions and the time-dependent delayed energy. The MPACT/CTF interfaces are updated to transfer the explicit heat sources. The new model agrees well with MCNP for problems without thermal–hydraulic feedback. The MPACT/CTF coupled calculation for a hot-full-power quarter-core case using the new model exhibits a reduction of peak pin power by 2.3% and fuel centerline temperature by 17 K. The new model also shows a keff difference up to -100 pcm in depletion and a 30 K peak fuel temperature reduction in a reactivity-initiated-accident problem.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

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

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

Mechanics↗
DOE OSTI2020

Developing high safety Li-metal anodes for future high-energy Li-metal batteries: strategies and perspectives

Developing high-safety Li-metal anodes (LMAs) is extremely important for the application of high-energy Li-metal batteries (LMBs), especially Li–S and Li–O 2 battery systems. However, the notorious Li-dendrite growth problem results in serious safety concerns for any energy storage application. Through a recent combination of interface-based science, nanotechnology-based solutions and characterization methods, the LMA is now primed for a technological boom. In this review, the recent emerging strategies and perspectives on LMAs are summarized, following which the current huge evolution in interfacial chemistry regulation, optimizing electrolyte components, designing a rational ‘host’ for lithium metal, optimizing “solid-state electrolytes” and other emerging strategies for developing high-safety LMAs is highlighted. Furthermore, several state-of-the-art in situ/operando synchrotron-based X-ray techniques for high safety LMB research are introduced. With the further development of LMAs in the future, subsequent application in high energy LMBs is to be expected.

25 ENERGY STORAGE↗
DOE OSTI2025

Quantum Simulators and Applications on Quantum Framework

Simulating quantum circuits is essential for validating quantum algorithms. However, no single simulator consistently performs best - efficiency depends on circuit structure, entanglement, and depth. In this work, we integrate Qiskit-Aer (state-vector and matrix product state) and QTensor, a tree-tensor-network based simulator, into the Quantum Framework (QFw), a modular platform that supports multiple quantum backends via a unified interface. We also enable distributed quantum approximate optimization algorithm (DQAOA) application compatibility with QFw, allowing sub-problems to be solved in parallel at scale. We then benchmark DQAOA and TFIM (transverse field Ising model) circuits across supported simulators, showing how performance varies significantly with problem type. All simulations are deployed on the Frontier supercomputer using QFw's MPI-based orchestration for distributed, multinode execution. These results underscore the need for simulatoragnostic infrastructure to enable systematic evaluation and highperformance scaling of quantum workloads. QFw provides a practical and extensible path toward reproducible quantum algorithm development across diverse application domains.

Chundury, Srikar [ORNL] (ORCID:0009000183359259)↗
DOE OSTI2021

Methods and Usability Enhancements in Shift for Non-LWR Applications

Several development and analysis tasks were undertaken in FY21 under the Nuclear Energy Advanced Modeling and Simulation program to enhance modeling of non light–water reactors (LWRs) with Shift. Specifically, these efforts targeted enhancements for tristructural isotropic (TRISO) fuel modeling. A new Shift user interface was developed that allows for much better usability and ease of modeling for non-LWR problems and TRISO fuel. Performance studies were conducted using an HTR-10 fuel pebble model by comparing different geometry packages in Shift, KENO-VI, and Serpent. These studies showed that the new geometry package in Shift performs well compared to Serpent for TRISO fuel modeling with consistent tracking options between both packages. The studies also identified the most critical areas of improvement for more efficiently performing Monte Carlo transport on TRISO fuel models with Shift. Tally calculations in Shift were optimized for non-LWR cross section generation and depletion calculations, and areas for further optimization and accuracy improvements were identified. Finally, initial collaboration efforts were formed between Idaho National Laboratory, Argonne National Laboratory, and the Nuclear Regulatory Commission to use Shift for Comprehensive Reactor Analysis Bundle support.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Progress in Genetic Algorithm Fitting of X-Ray Fluorescence Data for Absorption Spectroscopy of Liquid Surfaces

X-ray fluorescence (XRF) spectroscopy of liquid surfaces is a promising approach to elemental analysis of forensic samples. Elements in the mixture can be identified from the XRF data because chemical elements are associated with unique spectra. This report describes progress in the development of automated identification of unknown mixtures of elements from XRF spectra using a genetic algorithm (GA). The XRF spectra are measured at the Advanced Photon Source (APS) at the liquid-liquid interface during solvent extraction process. The usage of the GA is demonstrated and future implementation for related problems are then discussed.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

A Scalable Parallel Hypergraph Generator (HyGen)

Graphs are extensively used to model real-world complex systems. An edge in a graph can model pairwise relationships. However, multiway relationships (connections between three or more vertices) are common in many complex systems such as cellular process, image segmentation, and circuit design. A graph edge cannot model multiway relationships. A hypergraph, which can connect more than two vertices, is thus a better option to model multiway relationships. A large-scale hypergraph analysis has the potential to find useful insights from a complex system and assist in knowledge discovery. Currently a limited number of hypergraphs exists that are representative of real-world datasets. Moreover, real-world hypergraph datasets are small in size and inadequate to incorporate future needs. A graph generator that can produce large-scale synthetic hypergraphs can solve the above mentioned problems. In this paper, we present a scalable parallel hypergraph generator (HyGen) based on the Message Passing Interface (MPI) standard. To generate hypergraphs, HyGen takes the following parameter values as inputs: i) number of vertices, ii) number of hyperedges, iii) number of clusters, iv) vertex distribution, v) hyperedge distribution, vi) local cluster cardinality, and vii) global cluster cardinality. We have demonstrated that HyGen can generate hypergraphs of various sizes in a scalable fashion. HyGen takes approximately four minutes to generate a hypergraph with 4.8 million vertices, 1.6 million hyperedges, and 800 clusters using 1,024 processes on a leadership class computing platform. Our strong and weak scaling experiments on supercomputers demonstrate that HyGen can quickly create large-scale hypergraphs in a parallel manner, thus providing a useful capability for hypergraph analysis.

Hasan, S M Shamimul↗
DOE OSTI2022

An OpenMP GPU-offload implementation of a non-equilibrium solidification cellular automata model for additive manufacturing

Here, in this paper, performance strategies on GPU-based HPC platforms of a cellular automata (CA) simulation code for non-equilibrium solidification, including nucleation, grain growth, solute partitioning and transport for the metal additive manufacturing (AM) process are investigated using OpenMP 4.5. To accurately report the speed-up for multicore CPUs and GPUs, a rigorous performance analysis employed optimizations appropriate for both CPU-only code (baseline) and GPU offload codes for an isothermal test problem. The performance results on Summit at the Oak Ridge Leadership Computing Facility indicate that using a precomputed list of interface cells significantly decreased the wall-clock time on GPUs. The speedup due to GPU acceleration was evaluated for a full Summit node and measured to be 1.8X when comparing a 6 MPI tasks run with 6 GPUs versus 36 MPI tasks on the CPU only. That speed-up was found to be 7.9X when comparing 6 MPI tasks with 6 GPUs versus the 6 MPI tasks running on the CPU only. Performance measurements showed that system total time is almost constant for runs with more than 96 MPI tasks (or GPUs), indicating that the GPU-accelerated code showed an excellent weak scaling performance. Finally, a rapid directional solidification problem was considered to demonstrate the CA code capability on Summit. It was found that a mesh size of at least 0.05 μm is recommended for the AM-like simulations in order to obtain accurate elongated grain microstructure and elongated subgrain features, which are in qualitative good agreement with experimental data. The results presented in this study indicate that the performance strategies on GPU-based HPC platforms for the CA code are appropriate for novel HPC exascale platforms.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A moving discontinuous Galerkin finite element method with interface condition enforcement for compressible flows

A variation of moving discontinuous Galerkin finite element method with interface condition enforcement (MDG-ICE) is developed for solving the compressible Euler equations. The MDG-ICE method, originating from the work of Corrigan et al. [1], [2], [3], [4], is based on the space-time DG formulation, where both flow field and grid geometry are considered as independent variables and the conservation laws are enforced both on discrete elements and element interfaces. The element conservation laws are solved in the standard discontinuous solution space to determine conservative quantities, while the interface conservation is enforced using a variational formulation in a continuous space to determine discrete grid geometry. The resulting over-determined system of nonlinear equations arising from the MDG-ICE formulation can then be solved in a least-squares sense, leading to an unconstrained nonlinear least-squares problem that is regularized and solved by Levenberg-Marquardt method. A number of numerical experiments for both 1D unsteady and 2D steady state compressible flow problems are conducted to assess the accuracy and robustness of the MDG-ICE method. Numerical results obtained indicate that the MDG-ICE method is able to implicitly detect and track all types of discontinuities via interface conservation enforcement and satisfy the conservation law on both elements and interfaces via grid movement and grid management, demonstrating that an exponential rate of convergence for Sod and Lax-Harden shock tube problems can be achieved and highly accurate solutions without overheating to both double-rarefaction wave and Noh problems can be obtained.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Generalized master equation for particle transport in binary random media with renewal statistics

Particle transport in binary stochastic mixtures is classically modeled assuming Markovian or exponential mixing statistics but in many applications material memory invalidates the Markov assumption. For non-Markovian mixing characterized by alternating renewal processes, a transport-theoretic framework is presented that provides an exact description of transport in nonscattering random binary media with general non-exponential statistics. Our approach is to Markovianize the problem by augmenting the {material type, particle flux} state space with the age or distance from the last interface. A Chapman-Kolmogorov equation is formulated for the joint probability density of the material type, particle flux, and age, and subsequently reduced to a generalized Master equation (GME) in differential form. This constitutes the primary result of this work. A state-updating Monte Carlo algorithm consistent with the GME is developed and benchmarked against analytical solutions for multiple chord-length laws. For purely absorbing renewal statistical media, the GME reproduces analytical benchmarks for the equilibrium age distribution, interior mean/variance of material-conditioned fluxes, and boundary transmittance. Simulations further demonstrate that a Markov (exponential) approximation of non-exponential statistics can introduce large errors in transmittance and interior flux profiles. Lastly, the reintroduction of memory due to scattering is briefly addressed through heuristic considerations.

Fluctuations & noise↗
DOE OSTI2026

Nonequilibrium dynamics of two-level systems directly after cryogenic alternating bias

Two-level systems (TLSs) are a dominant decoherence mechanism in superconducting qubits, microwave kinetic inductance detectors, and quantum dots. TLSs are tunneling states commonly found in amorphous materials and interfaces, which interact electromechanically, leading to energy loss. Fluctuations in the frequency of the TLSs are a significant problem for qubit operation and force recalibration every few hours. In this study, we probe the effect of alternating bias at cryogenic temperatures on TLSs in a purpose-built LC oscillator. When an in situ alternating bias is applied, the steady-state spectra of the TLSs disappear. Spectroscopy reveals transient behavior, in which the TLS frequency fluctuates on the order of minutes. Thermal cycling above 10 K reverses these effects, restoring the TLS spectrum to its original state. Importantly, the intrinsic loss tangent of the LC oscillator remains unchanged before and after the application of the alternating bias. We propose that the disappearance of the steady-state spectrum is caused by nonequilibrium energy buildup from strain in the oxide film introduced by the pulsed voltage-bias sequence. Understanding this nonequilibrium energy could shed light on TLS fluctuations and reduce the calibration requirements for qubits.

general physics↗