Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coupled system”

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 OSTI2022

Cumulant methods for electron-phonon problems. I. Perturbative expansions

Here in this work, we investigate the ability of the cumulant expansion (CE) to capture one-particle spectral information in electron-phonon coupled systems at both zero and finite temperatures. In particular, we present a comprehensive study of the second- and fourth-order CEs for the one-dimensional Holstein model as compared with numerically exact methods. We investigate both finite sized systems as well as the approach to the thermodynamic limit, drawing distinctions, and connections between the behavior of systems in and away from the thermodynamic limit that enable a greater understanding of the ability of the CE to capture real-frequency information across the full range of wave vectors. We find that for zero electronic momentum, the spectral function is well described by the second-order CE at low and high temperatures. However, for nonzero electronic momenta, the CE is only accurate at high temperature. We analyze the fourth-order cumulant and find that while it improves the description of the short-time dynamics encoded in the one-particle Green's function, it can introduce divergences in the time domain as well as unphysical negative spectral weight in the spectral function. When well-behaved, the fourth-order CE does provide notable accurate corrections to the second-order CE. Finally, we use our results to comment on the use of the CE as a tool for calculating transport behavior in the realistic ab initio modeling of materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Hybrid Kohn-Sham+Thomas-Fermi scheme for high-temperature density functional theory

We report density functional theory has revolutionised our understanding of dense, strongly coupled systems in the condensed-matter limit. At elevated temperatures, such as those found in warm dense matter and dense plasma, the computational cost of calculating enough states to accurately describe a thermal ensemble of electrons is a major obstacle to practical Kohn-Sham calculations. This problem has recently been tackled from a number of approaches, including a plane-wave approximation for electrons at high energy. In this work, we demonstrate how the Thomas-Fermi density of states, equivalent to a local plane-wave approximation, can yield improved results for thermodynamic quantities, electron density profiles, and interatomic forces while requiring fewer states than existing methods.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Computation of high-frequency magnetoelastic waves in layered materials

Here, the direct calculation of magnetoelastic wave dispersion in layered media is presented using an efficient, accurate computational technique. The governing, coupled equations for elasticity and magnetism, the Navier and Landau-Lifshitz equations, respectively, are linearized to form a quadratic eigenvalue problem that determines a complex web of wave-number–frequency dispersion branches and their corresponding mode profiles. Numerical discretization of the eigenvalue problem via a spectral collocation method (SCM) is employed to determine the complete dispersion maps for both a single, finite-thickness magnetic layer and a finite magnetic-nonmagnetic double-layer. The SCM, previously used to study elastic waves in nonmagnetic media, is fast, accurate, and adaptable to a variety of sample configurations and geometries. Emphasis is placed on the extremely high-frequency regimes being accessed in ultrafast magnetism experiments. The dispersion maps and modes provide insight into how energy propagates through the coupled system, including how energy can be transferred between elastic- and magnetic-dominated waves as well as between different layers. The numerical computations for a single layer are further understood by a simplified analytical calculation in the high-frequency, exchange-dominated regime where the resonance condition required for energy exchange (an anticrossing) between quasi-elastic and quasi-magnetic dispersion branches is determined. Nonresonant interactions are shown to be well approximated by the dispersion of uncoupled elastic and magnetic waves. The methods and results provide fundamental theoretical tools to model and understand current and future magnetic devices powering spintronic innovation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Complex paths around the sign problem

The Monte Carlo evaluation of path integrals is one of a few general purpose methods to approach strongly coupled systems. It is used in all branches of physics, from QCD and nuclear physics to the correlated electron systems. However, many systems of great importance (dense matter inside neutron stars, the repulsive Hubbard model away from half filling, and dynamical and nonequilibrium observables) are not amenable to the Monte Carlo method as it currently stands due to the so-called sign problem. Here, a new set of ideas recently developed to tackle the sign problem based on the complexification of field space and the Picard-Lefshetz theory accompanying it is reviewed. The mathematical ideas underpinning this approach, as well as the algorithms developed thus far, are described together with nontrivial examples where the method has already been proved successful. Directions of future work, including the burgeoning use of machine learning techniques, are delineated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Transient recovery problem in acoustics: A multichannel model-based deconvolution approach

Critical acoustical systems operating in complex environments contaminated with disturbances and noise offer an extreme challenge when excited by out-of-the-ordinary, impulsive, transient events that can be undetected and seriously affect their overall performance. Transient impulse excitations must be detected, extracted, and evaluated to determine any potential system damage that could have been imposed; therefore, the problem of recovering the excitation in an uncertain measurement environment becomes one of multichannel deconvolution. Recovering a transient and its initial energy has not been solved satisfactorily, especially when the measurement has been truncated and only a small segment of response data is available. The development of multichannel deconvolution techniques for both complete and incomplete excitation data is discussed, employing a model-based approach based on the state-space representation of an identified acoustical system coupled to a forward modeling solution and a Kalman-type processor for enhancement and extraction. In conclusion, synthesized data are utilized to assess the feasibility of the various approaches, demonstrating that reasonable performance can be achieved even in noisy environments.

42 ENGINEERING↗
DOE OSTI2021

On a nonlocal Cahn–Hilliard model permitting sharp interfaces

A nonlocal Cahn–Hilliard model with a non-smooth potential of double-well obstacle type that promotes sharp interfaces in the solution is presented. To capture long-range interactions between particles, a nonlocal Ginzburg–Landau energy functional is defined which recovers the classical (local) model as the extent of nonlocal interactions vanish. In contrast to the local Cahn–Hilliard problem that always leads to diffuse interfaces, the proposed nonlocal model can lead to a strict separation into pure phases of the substance. In this work, the lack of smoothness of the potential is essential to guarantee the aforementioned sharp-interface property. Mathematically, this introduces additional inequality constraints that, in a weak formulation, lead to a coupled system of variational inequalities which at each time instance can be restated as a constrained optimization problem. We prove the well-posedness and regularity of the semi-discrete and continuous in time weak solutions, and derive the conditions under which pure phases are admitted. Moreover, we develop discretizations of the problem based on finite element methods and implicit–explicit time-stepping methods that can be realized efficiently. Finally, we illustrate our theoretical findings through several numerical experiments in one and two spatial dimensions that highlight the differences in features of local and nonlocal solutions and also the sharp interface properties of the nonlocal model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Sparsified Time-dependent PDEs FNO (STFNO) v1.0.0

STFNO (Sparsified Time-dependent PDEs FNO code) is an extension of the popular Fourier Neural Operator (FNO) architecture to the solution of coupled systems of time-dependent partial differential equations. STFNO leverages the sparsified dependencies on the field quantities based on the semi-discretiezed form of the PDEs, enabling significant reduction in the number of model parameters. STFNO has been extensively tested on two fusion simulation codes, NIMROD and GTC, and can be easily tailored to other systems of PDEs.

Rahman, Mustafa [Lawrence Berkeley National Labora↗
DOE OSTI2025

Photon-magnon coupling using gain-assisted spoof-localized surface plasmons

Improving the photon-magnon coupling strength can be done by tuning the structure of microwave resonators to better interact with the magnon counterpart. Planar resonators accommodating unconventional photon modes beyond the half- and quarter-wavelength designs have been explored due to their optimized mode profiles and potentials for on-chip integration. Here, we designed and fabricated an actively controlled ring resonator supporting the spoof localized surface plasmons (LSPs), and implemented it in the investigation of photon-magnon coupling for hybrid magnonic applications. We demonstrated gain-assisted photon-magnon coupling with the YIG magnon mode under several different sample geometries. The achieved coupling amplification largely benefits from the high quality factor (Q-factor) due to the additional gain provided by a semiconductor amplifier, which effectively increases the Q-factor from a nearly null state (passive resonance) to more than 1000 for a quadrupole LSP mode. Our results suggest an additional control knob for manipulating photon-magnon coupled systems exploiting external controls of gain and loss.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Snow thermal conductivity and conductive flux in the Central Arctic: Estimates from observations and implications for models

During the Arctic winter, the conductive heat flux through the sea ice and snow balances the radiative and turbulent heat fluxes at the surface. Snow on sea ice is a thermal insulator that reduces the magnitude of the conductive flux. The thermal conductivity of snow, that is, how readily energy is conducted, is known to vary significantly in time and space from observations, but most forecast and climate models use a constant value. This work begins with a demonstration of the importance of snow thermal conductivity in a regional coupled forecast model. Varying snow thermal conductivity impacts the magnitudes of all surface fluxes, not just conduction, and their responses to atmospheric forcing. Given the importance of snow thermal conductivity in models, we use observations from sea ice mass balance buoys installed during the Multidisciplinary drifting Observatory for the Study of Arctic Climate expedition to derive the profiles of thermal conductivity, density, and conductive flux. From 13 sites, median snow thermal conductivity ranges from 0.33 W m −1 K −1 to 0.47 W m −1 K −1 with a median from all data of 0.39 W m −1 K −1 from October to February. In terms of surface energy budget closure, estimated conductive fluxes are generally smaller than the net atmospheric flux by as much as 20 W m −2 , but the average residual during winter is −6 W m −2 , which is within the uncertainties. The spatial variability of conductive heat flux is highest during clear and cold time periods. Higher surface temperature, which often occurs during cloudy conditions, and thicker snowpacks reduce temporal and spatial variability. These relationships are compared between observations and the coupled forecast model, emphasizing both the importance and challenge of describing thermodynamic parameters of snow cover for modeling the Arctic as a coupled system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Entanglement islands in higher dimensions

It has been suggested in recent work that the Page curve of Hawking radiation can be recovered using computations in semi-classical gravity provided one allows for ``islands" in the gravity region of quantum systems coupled to gravity. The explicit computations so far have been restricted to black holes in two-dimensional Jackiw-Teitelboim gravity. In this note, we numerically construct a five-dimensional asymptotically AdS geometry whose boundary realizes a four-dimensional Hartle-Hawking state on an eternal AdS black hole in equilibrium with a bath. We also numerically find two types of extremal surfaces: ones that correspond to having or not having an island. The version of the information paradox involving the eternal black hole exists in this setup, and it is avoided by the presence of islands. Thus, recent computations exhibiting islands in two-dimensional gravity generalize to higher dimensions as well.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Design Optimization of A generic Fissile Solution for Mo99 production using Electron Beam-based Neutron Generator using MCNP+CFD

Molybdenum-99 ( 99 Mo) is a medical isotope used in 80% of all medical imaging procedures today. However, 99 Mo is not directly administered to patients for imaging; its decay product, technetium-99m ( 99 mTc), is a metastable isotope with a half-life of ~6 hours. Such a short half-life does not allow for production, separation, and application to a patient before decaying further. Therefore, the solution is to produce the precursor to 99 mTc, 99 Mo. 99 Mo has a half-life of 66 hours, long enough to produce the material and send it to the desired location before it fully decays. Under National Nuclear Security Administration’s (NNSA) highly enriched uranium (HEU) minimization mission, the Material Management and Minimization (M3) office leads the Molybdenum-99 (Mo-99) program. A major goal of the Mo-99 program is to develop methods of producing 99 Mo using low enriched uranium (LEU) opposed to the traditionally used HEU domestically in the U.S. Internationally, LEU is the new standard for use in research reactors and isotope production facility. In this study, we introduce a generic fissile solution system for 99 Mo production facility meeting the design requirement of 1) LEU usage as fissile material and 2) subcritical fissile system coupled with a horizontal electron beam (E-beam) accelerator. The focus of this research is to optimize a fissile solution configuration for maximum 99 Mo production yield using MCNP model, and ultimately implement the optimized design with computational fluid dynamics (CFD) model to understand thermal hydraulic behavior and solution convection characteristic. The report consists of two parts: 1) isotope yield calculation and criticality analysis with a wide range of parametric test matrix, and 2) thermal hydraulic analysis for solution mixing and cooling assessment.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

CO 2 to Bioplastics: Beneficial Re-use of Carbon Emissions from Coal-fired Power Plants using Microalgae

This project sought to address the technical and economic barriers to carbon dioxide (CO 2 ) capture and utilization using microalgae. Specifically, a dual photobioreactor (PBR)/open raceway pond (ORP) cultivation system was evaluated with respect to capital and operational costs, productivity, and culture health, and compared to an ORP-only cultivation system. In the dual (or “hybrid”) system, a cyclic flow PBR was used to provide inoculum for two open raceway ponds after each harvest, with the hypothesis that this PBR + pond system should result in improved performance compared to traditional raceway ponds. Two other raceway ponds were operated conventionally as a control, the experiments being performed at Duke Energy’s East Bend Station in northern KY using coal-derived flue gas as the CO 2 source. Although algae productivity achieved was mediocre due to the fact that the experiments had to be performed late in the growing season when climatic conditions were not optimal for high growth, the hybrid cultivation system showed a higher level of algae productivity (statistically significant) compared to the traditional ponds. In order to realize the maximum value of the algal biomass produced, fractionation of the biomass was examined. Bioplastic compounding and material characterization was subsequently performed using three feed stocks: whole biomass, lipid-extracted biomass and the proteinaceous residue from full fractionation. In general, the whole and lipid-extracted biomass gave similar results in terms of the properties of the resulting bioplastics, while the proteinaceous residue showed promise for the production of a polybutylene adipate terephthalate (PBAT) blend. Finally, sustainability assessment of a putative biorefinery system was conducted, including techno-economic and life cycle impact assessment. Nine different production scenarios were considered, comprising combinations of the three different biomass growth architectures (ORP, PBR and dual PBR-ORP systems) coupled with three different algae biomass processing pathways (drying only, lipid extraction and fractionation). Results show that the minimum selling price of the bioplastic feed stock (BPFS) is within the realm of economic competition with prices as low as $970 USD tonne -1 . Additionally, life cycle impact assessment results indicate drastic improvements in performance of the produced BPFS, with reductions in greenhouse gas emissions ranging between 67 and 116% compared to a petroleum based plastic feedstock.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2020

Parallel-in-Time Simulation of Lindblad's Equation

Constructing fast quantum logic gates is critical to building a scalable quantum computer. We consider a qudit, a quantum version of a bit that can take an arbitrary number of states, coupled with a cavity. In this project, we wish to force the qudit to reach the 0-state, for any possible initial state. The coupled system changes in time according to Lindblad’s equation, an ordinary differential equation on the density matrix of the quantum system. Lindblad’s equation contains some parameters that we can control, so-called control functions. We seek control functions which force the qudit to the 0-state within 2 microseconds, which is much faster than what is currently done in practice. The search method is gradient descent, a numerical optimization method that uses gradient information to iteratively improve the control parameters. My contribution to this project is an attempt to speed up the computation of the gradient. It currently takes about 40 seconds to compute the gradient which involves solving a set of ODEs sequentially. Current supercomputers have thousands of cores, but sequential computations can only make use of 1 core at a time. We wish to divide up the work better, so that we can use many more cores at once. To this end, we have implemented the Multigrid Reduction in Time (MGRIT) algorithm. We perform a systematic parameter search on how to best apply this algorithm. Results indicate a 25 percent speed up for solving Lindblad’s equation and determining how close the final state is the 0-state.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Thermal analysis of non-actinide bearing salt

The Netzsch STA 409 CD 403/5/G STA-MS Skimmer Coupling System was used to perform thermal analysis on FLiNaK to establish a method for characterizing vapor pressures over molten salts. The measurements include differential scanning calorimetry, thermogravimetry, and mass spectrometry to generate comprehensive thermochemical property data as stand-alone information to aid in MSR design and licensing as well as inputs into the Molten Salt Thermodynamic Database (MSTDB). The FY20 focus is to establish procedures and quantify uncertainty levels for measurements of non-actinide compositions for moving toward actinide bearing salt measurements in FY21.

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

A New Approach to Fundamental Mechanism Discovery in Polymer Upcycling

We present a new experimental methodology for detailed experimental investigations of depolymerization reactions over solid catalysts. This project aims to address a critical need in fundamental research on chemical upcycling of polymers – the lack of rapid, sensitive, isomerselective probing techniques for the detection of reaction intermediates and products. Our method combines a heterogeneous catalysis reactor for the study of multiphase (gas/polymer melt/solid) systems, coupled to a vacuum UV photoionization time-of-flight mass spectrometer. This apparatus draws on our expertise in probing complex gas-phase chemistry and enables highthroughput, detailed chemical speciation measurements of the gas phase above the catalyst, providing valuable information on the heterogeneous catalytic reactions. Using this approach, we investigated the depolymerization of high-density polyethylene (HDPE) over Ir-doped zeolite catalysts. We showed that the product distribution was dominated by low-molecular weight alkenes with terminal C=C double bonds and revealed the presence of many methyl-substituted alkenes and alkanes, suggesting extensive methyl radical chemistry. In addition, we investigated the fundamental reactivity of model oligomer molecules n-butane and isobutane over ZSM-5 zeolites. We demonstrated the first direct detection of methyl radical intermediates, confirming the key role of methyl in zeolite-catalyzed activation of alkanes. Our results show the potential of this experimental method to achieve deep insight into the complex depolymerization reactions and pave the way for detailed mechanistic studies, leading to increased fundamental understanding of key processes in chemical upcycling of polymers.

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

HYSPLIT/MACCS Atmospheric Dispersion Model Technical Documentation and Benchmark Analysis

The nuclear accident consequence analysis code MACCS has traditionally modeled dispersion during downwind transport using a Gaussian plume segment model. MACCS is designed to estimate consequence measures such as air concentrations and ground depositions, radiological doses, and health and economic impacts on a statistical basis over the course of a year to produce annualaveraged output measures. The objective of this work is to supplement the Gaussian atmospheric transport and diffusion (ATD) model currently in MACCS with a new option using the HYSPLIT model. HYSPLIT/MACCS coupling has been implemented, with HYSPLIT as an alternative ATD option. The subsequent calculations in MACCS use the HYSPLIT-generated air concentration, and ground deposition values to calculate the same range of output quantities (dose, health effects, risks, etc.) that can be generated when using the MACCS Gaussian ATD model. Based on the results from the verification test cases, the implementation of the HYSPLIT/MACCS coupling is confirmed. This report contains technical details of the HYSPLIT/MACCS coupling and presents a benchmark analysis using the HYSPLIT/MACCS coupling system. The benchmark analysis, which involves running specific scenarios and sensitivity studies designed to examine how the results generated by the traditional MACCS Gaussian plume segment model compare to the new, higher fidelity HYSPLIT/MACCS modeling option, demonstrates the modeling results that can be obtained by using this new option. The comparisons provided herein can also help decision-makers evaluate the potential benefit of using results based on higher fidelity modeling with the additional computational burden needed to perform the calculations. Three sensitivity studies to investigate the potential impact of alternative modeling options, regarding 1) input meteorological data set, 2) method to estimate stability class, and 3) plume dispersion model for larger distances, on consequence results were also performed. The results of these analyses are provided and discussed in this report.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Identification and Quantification of Photosynthetic Pigments in Algae (Laboratory Analytical Procedure (LAP))

The Laboratory Analytical Procedure (LAP) outlined here describes a method to quantitatively extract phytopigments from microalgae biomass, as well as to identify and quantify individual pigments based on separation and detection with a High-Performance Liquid Chromatography (HPLC) system coupled to a Diode Array Detector (DAD). Pigments were extracted with greater than 95 % extraction efficiency. Chromatographic conditions allowed for isomeric resolution between pigments and identification based on UV/Vis spectra and comparison to analytical standards.

09 BIOMASS FUELS↗
DOE OSTI2023

Block smoothers and generalized ideal interpolation in AMG (Final Report)

The Pennsylvania State University (“Subcontractor”) worked on developing new parallel algebraic multilevel methods suitable for solving PDEs. Specifically, work on the design of multigrid solvers for coupled systems of partial differential equations arising in numerical modeling of various applications was completed. A main emphasis was on the design of new ideal algebraic multigrid interpolation for problems such as Maxwell’s equations where block smoothers are needed and the standard form of ideal interpolation is not an effective choice.

97 MATHEMATICS AND COMPUTING↗