Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial discretization”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Alkene adsorption and cracking on acidic zeolites – A gradual process of understanding

Cracking of light alkenes on solid acids and, in particular on zeolites, is ubiquitous in refineries. Different zeolite framework topologies in the catalysts used provide spatial constraints for reactants and products, leading to discrete reaction routes and varying selectivities. The conversion of alkenes in cracking proceeds either via monomolecular or bimolecular reaction pathways. In both routes, the reacting alkene adsorbs first on Brønsted acid sites forming the initial state (e.g., π-complex, carbenium ion or alkoxide), followed by β-cleavage in the carbenium transition state, generating a smaller alkene and a carbenium ion (monomolecular route) or forming a dimer by alkene addition that cracks subsequently (bimolecular route). The current state of understanding alkene sorption and cracking is critically discussed in this paper and four questions are formulated that may be the central questions from our perspective in further developing new generations of catalysts and processes for alkene conversion.

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

Realizability-preserving discontinuous Galerkin method for spectral two-moment radiation transport in special relativity

Here we present a realizability-preserving numerical method for solving a spectral two-moment model to simulate the transport of massless, neutral particles interacting with a steady background material moving with relativistic velocities. The model is obtained as the special relativistic limit of a four-momentum-conservative general relativistic two-moment model. Using a maximum-entropy closure, we solve for the Eulerian-frame energy and momentum. The proposed numerical method is designed to preserve moment realizability, which corresponds to moments defined by a nonnegative phase-space density. The realizability-preserving method is achieved with the following key components: (i) a discontinuous Galerkin phase-space discretization with specially constructed numerical fluxes in the spatial and energy dimensions; (ii) a strong stability-preserving implicit-explicit time-integration method; (iii) a realizability-preserving conserved to primitive moment solver; (iv) a realizability-preserving implicit collision solver; and (v) a realizability-enforcing limiter. Component (iii) is necessitated by the closure procedure, which closes higher order moments nonlinearly in terms of primitive moments. The nonlinear conserved to primitive and the implicit collision solves are formulated as fixed-point problems, which are solved with custom iterative solvers designed to preserve the realizability of each iterate. With a series of numerical tests, we demonstrate the accuracy and robustness of this discontinuous-Galerkin-implicit-explicit method.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Imaging Nanoscale Energy Transport and Conversion with Ultrafast Electron Microscopy (Final Technical Report)

Light-matter interactions are ubiquitous in nature and reside at the heart of innumerable technologies. The cascade of processes that occur when a material absorbs a photon of light are exceedingly complex and are interwoven in both space and time, rendering precise determination of the atomic-scale and ultrafast mechanisms immensely challenging. The advent of methods for generating short pulses of light several decades ago led to major advances in understanding the initial moments of light absorption and the resultant effects, though directly interrogating the response of the atoms within the material continued to prove challenging. More recently, methods for generating ultrashort pulses of X-rays and fast electrons have opened the way to probing photoinduced structural dynamics of a wide range of matter in multiple phases. An especially promising laboratory-scale method is ultrafast electron microscopy (UEM), wherein the modalities of conventional transmission electron microscopes (imaging, diffraction, spectroscopy) are extended into the femtosecond temporal regime. Here, the ultrafast imaging and diffraction modalities of UEM were used to study the transient structural aspects of photoexcitation of semiconducting materials (e.g., spatially-resolved electron-phonon coupling, excitation and emission of acoustic phonons, and discrete nanoscale scattering processes). The project had three main objectives: (1) determination of the excitation mechanisms of dense, hypersonic charge-carrier waves and the spatially mediated means by which they couple to the lattice via coherent phonon emission, especially with UEM imaging, wherein effects of nanoscale structural and morphological features on the coupling and relaxation dynamics are expected to dictate nucleation sites and preferred wave vectors, (2) elucidation of photoinduced acoustic-phonon seeding, emergence, propagation, and decay over nanoscale crystal regions and especially with respect to local strain fields and atomic-scale disruptions in lattice order, and (3) interwoven with the first two were efforts aimed at realizing combined angstrom-femtosecond spatiotemporal imaging with UEM. The outcomes and impacts of this project were the generation of new knowledge with respect to fundamental light-matter interactions and, especially, the spatially-mediated excitation and evolution of the structural response of materials following coherent photoexcitation. Importantly, the spatial and temporal resolutions of the UEM imaging modalities used are well-suited for such studies and enable spatially-resolved mechanisms to be determined as a function of atomic order, structural features, and morphology. In addition, ultrafast crystallographic measurements were used to correlate real- and reciprocal-space dynamics in order to determine atomic-scale preferential wave vectors and ultrafast scattering mechanisms, especially as dictated by specimen boundary conditions. The obtained results, as detailed in peer-reviewed publications and presentations, illustrate the importance of ultrafast, angstrom-scale real-space imaging for developing a comprehensive understanding of energy transport and conversion in materials.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Gauging anomalous unitary operators

Boundary theories of static bulk topological phases of matter are obstructed in the sense that they cannot be realized on their own as isolated systems. The obstruction can be quantified/characterized by quantum anomalies, in particular when there is a global symmetry. Similarly, topological Floquet evolutions can realize obstructed unitary operators at their boundaries. In this paper, we discuss the characterization of such obstructions by using quantum anomalies. As a particular example, here we discuss time-reversal symmetric boundary unitary operators in one and two spatial dimensions, where the anomaly emerges as we gauge the so-called Kubo-Martin-Schwinger (KMS) symmetry. We also discuss mixed anomalies between particle number conserving U(1) symmetry and discrete symmetries, such as C and CP, for unitary operators in odd spatial dimensions that can be realized at the boundaries of topological Floquet systems in even spatial dimensions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Improved treatment of multi-material cells in thermal radiation transport codes

High-energy-density physics simulations with non-conformal meshes of materials require a multi-material (MM) closure that affects the thermal radiation transport (TRT). We propose a set of novel closures that work for an arbitrary number of materials, both grey and multigroup energy discretizations, and any angular discretization (such as Sn, IMC, or diffusion). For each spatial cell, our closures let each species (ion and electron) of each material have its own temperature, density, and internal energy, but use a single radiation distribution that interacts with all materials within the cell. Our closures maintain energy conservation, do not incur increased computational cost in the TRT solve itself, are compatible with single-material TRT solvers, and do not make any temperature-equilibrium assumptions. Here we test our closures on a wide range of increasingly realistic problems and find them to be robust.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Capability Demonstration of a 3D CdZnTe Detector on a High-Altitude Balloon Flight

In collaboration between the University of Michigan and Los Alamos National Laboratory, a 3D position-sensing CdZnTe (CZT) detector prototype was built and integrated into a high-altitude balloon platform to evaluate its performance in a space-like mixed-radiation environment. The detector prototype, Orion Eagle, was designed to operate in near-vacuum environments without any temperature regulation. Orion Eagle was hand-launched from NASA’s Columbia Scientific Balloon Facility (CSBF) at Fort Sumner, NM on September 26, 2021, and successfully operated throughout a 9-hour flight, which reached 38.5 km in altitude. The flight met its objectives, successfully detecting atmospheric gamma rays and galactic cosmic rays, and raising the Technical Readiness Level from 4 to 6 for large-volume 3D CZT detector technology for space applications. Ionization tracks produced by charged particles create spatial signatures in the detector that are distinguishable from discrete gamma-ray interactions. Therefore, the 3D position-sensing capabilities using pixelated electrodes on a CZT detector can help enable discrimination of background charged particles from gamma-ray events without an anticoincidence shield. The potential for background rejection capability, ambient-temperature operation, gamma-ray coded-aperture and Compton imaging, and near High Purity Germanium (HPGe) energy resolution motivate the use of large-volume 3D CZT imaging spectrometers in future space missions.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Lectures on statistical mechanics

Presented here is a transcription of the lecture notes from Professor Allan N. Kaufman’s graduate statistical mechanics course Physics 212A and 212B at the University of California Berkeley from the 1972–1973 academic year. 212A addressed equilibrium statistical mechanics with topics: fundamentals (micro-canonical and sub-canonical ensembles, adiabatic law and action conservation, fluctuations, pressure, and virial theorem), classical fluids and other systems (equation of state, deviations from ideality, virial coefficients and van der Waals potential, canonical ensemble and partition function, quasistatic evolution, grand-canonical ensemble and partition function, chemical potential, simple model of a phase transition, quantum virial expansion, numerical simulation of equations of state, and phase transition), chemical equilibrium (systems with multiple species and chemical reactions, law of mass action, Saha equation, chemical equilibrium including ionization and excited states), and long-range interactions (including Coulomb, dipole, and gravitational interactions, Debye–Hückel theory, and shielding). 212B addressed nonequilibrium statistical mechanics with topics: fundamentals (definitions: realizations, moments, characteristic function, and discrete variables), Brownian motion (Langevin equation, fluctuation–dissipation theorem, spatial diffusion, Boltzmann’s H-theorem), Liouville and Klimontovich equations, Landau equation (derivation, elaboration, and H-theorem, and irreversibility), Markov processes and Fokker–Planck equation (derivations of the Fokker–Planck equation and a master equation), linear response and transport theory (linear Boltzmann equation, linear response theory of Kubo and Mori, relation of entropy production to electrical conductivity, transport relations and coefficients, normal mode solutions of the transport equations, sketch of a generalized Langevin equation method for transport theory), and an introduction to nonequilibrium quantum statistical mechanics.

plasma dynamics↗
DOE OSTI2023

11-th order of accuracy for numerical solution of 3-D Poisson equation with irregular interfaces on unfitted Cartesian meshes

For the first time the optimal local truncation error method (OLTEM) with 125-point stencils and unfitted Cartesian meshes has been developed in the general 3-D case for the Poisson equation for heterogeneous materials with smooth irregular interfaces. The 125-point stencils equations that are similar to those for quadratic finite elements are used for OLTEM. The interface conditions for OLTEM are imposed as constraints at a small number of interface points and do not require the introduction of additional unknowns, i.e., the sparse structure of global discrete equations of OLTEM is the same for homogeneous and heterogeneous materials. The stencils coefficients of OLTEM are calculated by the minimization of the local truncation error of the stencil equations. These derivations include the use of the Poisson equation for the relationship between the different spatial derivatives. Such a procedure provides the maximum possible accuracy of the discrete equations of OLTEM. In contrast to known numerical techniques with quadratic elements and third order of accuracy on conforming and unfitted meshes, OLTEM with the 125-point stencils provides 11-th order of accuracy, i.e., an extremely large increase in accuracy by 8 orders for similar stencils. The numerical results show that OLTEM yields much more accurate results than high-order finite elements with much wider stencils. The increased numerical accuracy of OLTEM leads to an extremely large increase in computational efficiency. Additionally, a new post-processing procedure with the 125-point stencil has been developed for the calculation of the spatial derivatives of the primary function. The post-processing procedure includes the minimization of the local truncation error and the use of the Poisson equation. It is demonstrated that the use of the partial differential equation (PDE) for the 125-point stencils improves the accuracy of the spatial derivatives by 6 orders compared to post-processing without the use of PDE as in existing numerical techniques. At an accuracy of 0.1% for the spatial derivatives, OLTEM reduces the number of degrees of freedom by 900 - 4∙10 6 times compared to quadratic finite elements. The developed post-processing procedure can be easily extended to unstructured meshes and can be independently used with existing post-processing techniques (e.g., with finite elements).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

An atomistic survey of shear coupling in asymmetric tilt grain boundaries and interpretation using the disconnections framework

Grain Boundaries (GB) play an important role in determining the behavior of polycrystalline materials. While the mechanisms of motion and associated shear response for symmetric tilt grain boundaries (STGBs) are well studied, the same is not true for asymmetric tilt grain boundaries (ATGBs) despite their greater prevalence in polycrystals. Here, this study aims to investigate the shear response of a large collection of asymmetric tilt grain boundaries (ATGBs) using molecular dynamics (MD) simulations and interpret the data using a discrete disconnections model that works remarkably well for STGBs. MD simulations of shear-driven ATGBs show that the plastic shear (shear coupling factor) in the region swept by a GB exhibits a complex dependence on the inclination angle, and this dependence changes with the misorientation of the GB. In addition, the shear response was observed to be highly sensitive to the applied shear rate and temperature. Recognizing the spatial and temporal scale limitations of MD simulations, we extended the discrete disconnections mesoscale model of Khateret al. (2012) to calculate the nucleation barriers of disconnection modes and predict the effective shear coupling of an ATGB. We observed that the mesoscale model’s predictions of the shear coupling factor of ATGBs do not agree with those observed in MD simulations. Finally, we examine the hypotheses of our mesoscale model that contribute to disagreements between MD simulations and the mesoscale model and propose improvements to the mesoscale model for future work.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Advances in hexagon mesh-based flow direction modeling

Watershed delineation and flow direction representation are the foundations of streamflow routing in spatially distributed hydrologic modeling. A recent study showed that hexagon-based watershed discretization has several advantages compared to the traditional Cartesian (latitude–longitude) discretization, such as uniform connectivity and compatibility with other Earth system model components based on unstructured mesh systems (e.g., oceanic models). Despite these advantages, hexagon-based discretization has not been widely adopted by the current generation of hydrologic models. One major reason is that there is no existing model that can delineate hexagon-based watersheds while maintaining accurate representations of flow direction across various spatial resolutions. In this study, we explored approaches such as spatial resampling and hybrid breaching-filling stream burning techniques to improve watershed delineation and flow direction representation using a newly developed hexagonal mesh watershed delineation model (HexWatershed). We applied these improvements to the Columbia River basin and performed 16 simulations with different configurations. The results show that (1) spatial resampling modulates flow direction around headwaters and provides an opportunity to extract subgrid information; and (2) stream burning corrects the flow directions in mountainous areas with complex terrain features.

58 GEOSCIENCES↗
DOE OSTI2026

A fourth order sharp immersed method for the incompressible Navier-Stokes equations with stationary and moving boundaries and interfaces

We propose a fourth order Navier-Stokes solver based on the immersed interface method (IIM), for flow problems with stationary and one-way coupled moving boundaries and interfaces. Our algorithm employs a Runge-Kutta-based projection method that maintains high-order temporal accuracy in both velocity and pressure for steady and unsteady velocity boundary conditions. Fourth order spatial accuracy is achieved through a novel fifth order IIM discretization scheme for the advection term, as well as existing high-order interface-corrected finite difference schemes for the other differential operators. Using a set of manufactured flow problems with stationary and moving boundaries, we demonstrate fourth order convergence of velocity and pressure in the infinity norm, both inside the domain and on the immersed boundaries. The solver’s performance is further validated through a range of practical flow simulations, highlighting its efficiency over a second order scheme. Finally, we showcase the ability of our immersed discretization scheme to handle interface-coupled multiphysics problems by solving a conjugate heat transfer problem with multiple immersed solids. Overall, the proposed approach robustly combines the efficiency of high order discretization schemes with the flexibility of immersed discretizations for flow problems with complex, moving boundaries and interfaces.

42 ENGINEERING↗
DOE OSTI2025

Quasinormal Coupled-Mode Analysis of Dynamic Gain in Exceptional-Point Lasers

One of the key features of lasers operating near exceptional points (EPs) is that the gain medium can support an oscillating population inversion above a pump threshold, leading to self-modulated laser dynamics. This unusual behavior opens up new possibilities for frequency comb generation and temporal modulation. However, the dynamic population inversion couples signals with different frequencies and is difficult to capture using conventional temporal coupled-mode theory (TCMT) based on stationary saturable gain. In this paper, we develop a perturbative coupled-mode analysis framework to capture the spatial-temporal dynamics of near-EP lasers. By decomposing discrete frequency generation into multiple excitations of resonant modes, our analysis establishes a minimal physical model that translates the local distribution of dynamic population-inversion into a resonant modal interpretation of laser gain. Furthermore, this work enables the exploration of unique properties in this self-time-modulated systems, such as time-varying scattering and nonreciprocal transmission.

Absorption↗
DOE OSTI2025

Spatial String Tension and Its Effects on Screening Correlators in a Thermal QCD Plasma

We calculate the spatial Wilson line correlator for 2 +1 flavor QCD using highly improved staggered quark discretization for fermions and in quenched QCD for a wide range of temperatures, from the chiral crossover temperature 𝑇 pc ≃ 156 MeV or the deconfinement temperature ≃300 MeV, respectively, up to 2 GeV. Extracting the spatial string tension for different lattice cutoffs and by performing a continuum extrapolation of this observable, we show that the soft (magnetic) gluons interact nonperturbatively even at temperatures ≳1 GeV. We provide incriminating evidences to demonstrate that dimensionally reduced effective theories can describe these soft quark and gluon quasi-particles for both quenched and 2 +1 flavor QCD, at temperatures 𝑇 ≳ 5⁢𝑇 pc . We also show for the first time the imprints of the nonperturbative pseudopotential in the properties of mesonic screening masses for temperatures ranging from 0.8 to 164 GeV in the quark-gluon plasma.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

On the three-dimensional spatial correlations of curved dislocation systems

Abstract Coarse-grained descriptions of dislocation motion in crystalline metals inherently represent a loss of information regarding dislocation-dislocation interactions. In the present work, we consider a coarse-graining framework capable of re-capturing these interactions by means of the dislocation-dislocation correlation functions. The framework depends on a convolution length to define slip-system-specific dislocation densities. Following a statistical definition of this coarse-graining process, we define a spatial correlation function which will allow the arrangement of the discrete line system at two points—and thus the strength of their interactions at short range—to be recaptured into a mean field description of dislocation dynamics. Through a statistical homogeneity argument, we present a method of evaluating this correlation function from discrete dislocation dynamics simulations. Finally, results of this evaluation are shown in the form of the correlation of dislocation densities on the same slip-system. These correlation functions are seen to depend weakly on plastic strain, and in turn, the dislocation density, but are seen to depend strongly on the convolution length. Implications of these correlation functions in regard to continuum dislocation dynamics as well as future directions of investigation are also discussed.

Anderson, Joseph Pierre (ORCID:0000000329760903)↗
DOE OSTI2022

SCALE Modeling of the Fast Spectrum Heat Pipe Reactor

As part of the severe accident analysis collaboration with Sandia National Laboratories (SNL) and the Nuclear Regulatory Commission (NRC), SCALE models were developed for a fast-spectrum heat pipe reactor. These models were based on the Idaho National Laboratory (INL) Design A concept, which is an alternative design to the Los Alamos National Laboratory (LANL) Special Purpose Reactor (SPR), also known as the Megapower reactor. The model contains 1,134 heat pipes, surrounded by hexagonal fuel elements, with a potassium working fluid; the fuel is UO 2 with 19.75 wt% 235 U enrichment. The model contains axial beryllium oxide (BeO) reflectors above and below the active fuel region along with a radial alumina reflector containing 12 B 4 C control drums. The center of the core is left unfueled to make room for two shutdown control rods, one annular and one solid. The active region of the core was discretized into twenty axial and five radial zones to analyze spatial variations in power and burnup. Infinite lattice unit cell sensitivity studies were used to perform verification between the SCALE and INL models. The eigenvalue results agreed well with the reported results to within roughly 50 percent mille (pcm). Full-core model verification was performed by analyzing system eigenvalues with differing configurations of control drum and shutdown rod positions. These full core results all had eigenvalue differences less than 310 pcm. Control drum and shutdown rod worths were also compared, with differences of 3.2% or less. Using the verified model, the isotopic inventory and decay heat, as well as temperature feedback coefficients, were calculated and provided to SNL as input to the MELCOR severe accident code to analyze potential releases from this class of reactor. The results of the MELCOR analysis are provided in a different report.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Journey Over Destination: Differentiable Sensor Placement Enhances Generalization [Poster]

The challenge of reconstructing spatial fields that change over time from limited sensor data has been a focal point for many research studies. Various machine learning methods have been used in attempts to address this complex issue, including convolutional neural networks. All of the proposed methods share a common requirement that the user needs to manually determine the sensor positions. This requirement remains a limiting factor in the ongoing quest for efficient learning and accurate field reconstruction. This study aims to present a method that enables a model to optimize sensor positions via backpropagation, thereby facilitating the model’s exploration of the spatial domain and enhancing sensor positioning effectively. Indexing naturally incorporates discrete decisions. This operation is nondifferentiable which is a requirement for the application of gradient-based optimization methods. We showcased its effectiveness by training an attention-based neural network, which achieved top-tier performance on two separate datasets. To our knowledge, this represents the first fully end-to-end differentiable workflow for enhancing sensor placement within a neural network model.

58 GEOSCIENCES↗
DOE OSTI2022

Development of a concurrent coupled Atomistic - Continuum model to predict the defect and grain structure for Additive Manufacturing process [Slides]

Development of a coupled atomistic-continuum model for metal AM process. Two-way coupling between atomistic (heat flux) and continuum (temperature) domain. Identification of grain and dislocation structure at the atomistic domain. Conversion of dislocation data from discrete atomistic to density form at the continuum. Obtain the spatial distribution of grain structure and dislocation density in additively manufactured metal polycrystal material.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Windowed least-squares model reduction for dynamical systems

Here we propose a windowed least-squares (WLS) approach for model reduction of dynamical systems. The proposed approach sequentially minimizes the time-continuous full-order-model residual within a low-dimensional space–time trial subspace over time windows. The approach comprises a generalization of existing model reduction approaches, as particular instances of the methodology recover Galerkin, least-squares Petrov–Galerkin (LSPG), and space–time LSPG projection. In addition, the approach addresses key deficiencies in existing model reduction techniques, e.g., the dependence of LSPG and space–time LSPG projection on the time discretization and the exponential growth in time exhibited by a posteriori error bounds for both Galerkin and LSPG projection. We consider two types of space–time trial subspaces within the proposed approach: one that reduces only the spatial dimension of the full-order model, and one that reduces both the spatial and temporal dimensions of the full-order model. For each type of trial subspace, we consider two different solution techniques: direct (i.e., discretize then optimize) and indirect (i.e., optimize then discretize). Numerical experiments conducted using trial subspaces characterized by spatial dimension reduction demonstrate that the WLS approach can yield more accurate solutions with lower space–time residuals than Galerkin and LSPG projection.

97 MATHEMATICS AND COMPUTING↗