Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “smoothing”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Machine Learning-Assisted Recovery of Delicate Kinetic Information from Transient Reactor Experiments

Identifying active sites and their roles in chemical reaction steps remains a vital challenge in heterogeneous catalysis. Transient experiments offer a unique way to probe active sites and distinguish subtle kinetic features. Although physics-based analysis methods may be well-developed, they can be highly susceptible to experimental noise, and smoothing methods may erase or even distort important features; a smooth curve is not always the best curve. We demonstrate a new workflow for the direct interpretation of intrinsic kinetic information from exit flux curves measured in transient reactor experiments. This workflow contains three artificial neural networks (ANNs), including a noise reducer, a concentration predictor, and a rate predictor to analyze experimental data, followed by the virtual TAP (VTAP) physics-based reactor model and density functional theory (DFT) calculations of adsorption energies on specific sites. We use this workflow to analyze the data from experiments titrating Pt/Al 2 O 3 and Pt/SiO 2 catalysts with carbon monoxide (CO) in the temporal analysis of products (TAP) reactor. Our workflow separates the time-evolving chemical reaction and mass transfer information contained in the TAP pulse response. The existence of strong- and weak-binding sites on the Pt/Al 2 O 3 catalyst is observed in the catalyst titration experiment in the transient reactor. The structures of the strong- and weak-binding sites are then identified by using DFT calculations. We find that the Pt/SiO 2 catalyst has only strong-binding sites, which aligns with the inactive support effect of SiO 2 . We demonstrate how machine learning methods provide unique insights with high-resolution data analysis that cannot be achieved by using state-of-the-art physics-based methods.

Adsorption↗
DOE OSTI2022

Crust and Upper Mantle Structure Beneath the Eastern United States

Abstract The Eastern United States (EUS) has a complex geological history and hosts several seismic active regions. We investigate the subsurface structure beneath the broader EUS. To produce reliable images of the subsurface, we simultaneously invert smoothed P‐wave receiver functions, Rayleigh‐wave phase and group velocity measurements, and Bouguer gravity observations for the 3D shear‐wave speed. Using surface‐wave observations (3–250 s) and spatially smoothed receiver functions, our velocity models are robust, reliable, and rich in detail. The shear‐wave velocity models fit all three types of observations well. The resulting velocity model for the eastern U.S. shows thinner crust beneath New England, the east coast, and the Mississippi Embayment (ME). A relatively thicker crust was found beneath the stable North America craton. A relatively slower upper mantle was imaged beneath New England, the east coast, and western ME. A comparison of crust thickness derived from our model against four recent published models shows first‐order consistency. A relatively small upper mantle low‐speed region correlates with a published P‐wave analysis that has associated the anomaly with a 75 Ma kimberlite volcanic site in Kentucky. We also explored the relationship between the subsurface structure and seismicity in the eastern U.S. We found that earthquakes often locate near regions with seismic velocity variations, but not universally. Not all regions of significant subsurface wave speed changes are loci of seismicity. A weak correlation between upper mantle shear velocity and earthquake focal mechanism has been observed.

58 GEOSCIENCES↗
DOE OSTI2022

Mild-to-wild plastic transition is governed by athermal screw dislocation slip in bcc Nb

Plastic deformation in crystals is mediated by the motion of line defects known as dislocations. For decades, dislocation activity has been treated as a homogeneous, smooth continuous process. However, it is now recognized that plasticity can be determined by long-range correlated and intermittent collective dislocation processes, known as avalanches. Here we demonstrate in body-centered cubic Nb how the long-range and scale-free dynamics at room temperature are progressively quenched out with decreasing temperature, eventually revealing intermittency with a characteristic length scale that approaches the Burgers vector itself. Plasticity is shown to be bimodal across the studied temperature regime, with conventional thermally-activated smooth plastic flow (‘mild’) coexisting with sporadic bursts (‘wild’) controlled by athermal screw dislocation activity, thereby violating the classical notion of temperature-dependent screw dislocation motion at low temperatures. An abrupt increase of the athermal avalanche component is identified at the critical temperature of the material. Our results indicate that plasticity at any scale can be understood in terms of the coexistence of these mild and wild modes of deformation, which could help design better alloys by suppressing one of the two modes in desired temperature windows.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Local spring effect in titanium-based layered oxides

Titanium-based layered oxides (TLOs) are one of the most promising electrode material families for sodium-ion batteries (NIBs) due to their smooth charge/discharge profiles and excellent cycle performance. However, the reaction mechanism of these materials, especially the reason for the disappearance of multiple voltage plateaus, is still not clear. Herein, two representative TLOs (Na 2/3 Ni 1/3 Ti 2/3 O 2 and Na 2/3 Co 1/3 Ti 2/3 O 2 ) with the same P2 crystal structure have been studied to scrutinize those unexplained issues. In situ synchrotron high-energy X-ray diffraction revealed a solid solution reaction mechanism for both, suggesting the absence of rigid phase transitions upon electrochemical cycling. An interesting “spring effect” of the TiO 6 octahedron, i.e., the reversible vibration of the central Ti atom inside the local octahedron upon electrochemical redox, was demonstrated by advanced X-ray absorption spectroscopy and theoretical calculations. Such an effect could suppress the rigid phase transitions, and result in smooth charge/discharge profiles and enhanced cycle stability. Finally, this work not only accounts for the disappearance of multiple voltage plateaus of TLOs for NIBs, but also provides an effective local-structure viewpoint to increase the cycle stability of electrode materials for other advanced battery systems.

25 ENERGY STORAGE↗
DOE OSTI2021

Computational general relativistic force-free electrodynamics: II. Characterization of numerical diffusivity

Scientific codes are an indispensable link between theory and experiment; in (astro-)plasma physics, such numerical tools are one window into the universe’s most extreme flows of energy. The discretization of Maxwell’s equations – needed to make highly magnetized (astro)physical plasma amenable to its numerical modeling – introduces numerical diffusion. It acts as a source of dissipation independent of the system’s physical constituents. Understanding the numerical diffusion of scientific codes is the key to classifying their reliability. It gives specific limits in which the results of numerical experiments are physical. We aim at quantifying and characterizing the numerical diffusion properties of our recently developed numerical tool for the simulation of general relativistic force-free electrodynamics by calibrating and comparing it with other strategies found in the literature. Our code correctly models smooth waves of highly magnetized plasma. We evaluate the limits of general relativistic force-free electrodynamics in the context of current sheets and tearing mode instabilities. We identify that the current parallel to the magnetic field (j ∥ ), in combination with the breakdown of general relativistic force-free electrodynamics across current sheets, impairs the physical modeling of resistive instabilities. We find that at least eight numerical cells per characteristic size of interest (e.g., the wavelength in plasma waves or the transverse width of a current sheet) are needed to find consistency between resistivity of numerical and of physical origins. High-order discretization of the force-free current allows us to provide almost ideal orders of convergence for (smooth) plasma wave dynamics. The physical modeling of resistive layers requires suitable current prescriptions or a sub-grid modeling for the evolution of j ∥ .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Anomalies in the topology of the temperature fluctuations in the cosmic microwave background: An analysis of the NPIPE and FFP10 data releases

We present a topological analysis of the temperature fluctuation maps from the Planck 2020 Data Release 4 NPIPE dataset and the Planck 2018 Data Release 3 FFP10 dataset. We performed a multiscale analysis in terms of the homology characteristics of the maps, invoking relative homology to account for the analysis in the presence of masks. We performed our analysis for a range of smoothing scales spanning sub- and super-horizon scales corresponding to a full width at half maximum (FWHM) of 5',10',20',40',80',160',320', and 640', and employed simulations based on the standard model for comparison, which assumes the initial fluctuation field to be an isotropic and homogeneous Gaussian random field. Examining the behavior of topological components, represented by the 0D homology group, we find the observations to be approximately 2σ or less deviant from the simulations for all resolutions and scales for the NPIPE dataset. For the FFP10 dataset, we detect a 2.96σ deviation between the observations and simulations at N = 128, FWHM = 80'. For the topological loops, represented by the first homology group, the simulations and observations are consistent within 2σ for most resolutions and scales for both the datasets. However, for the NPIPE dataset, we observe a high deviation between the observation and simulations in the number of loops at FWHM = 320', but at a low dimensionless threshold ν = –2.5. Under a Gaussian assumption, this would amount to a deviation of ~4σ. However, the distribution in this bin is manifestly non-Gaussian and does not obey Poisson statistics either. In the absence of a true theoretical understanding, we simply note that the significance is higher than what may be resolved by 600 simulations, yielding an empirical p-value of at most 0.0016. Specifically in this case, our tests indicate that the numbers arise from a statistically stable regime, despite being based on small numbers. For the FFP10 dataset, the differences are not as strong as for the NPIPE dataset, indicating a 2.77σ deviation at this resolution and threshold. The Euler characteristic, which is the alternating sum of the ranks of relative homology groups, reflects the deviations in the components and loops. To assess the significance of combined levels for a given scale, we employed the empirical and theoretical versions of the χ 2 test as well as the nonparametric Tukey depth test. Although all statistics exhibit a stable distribution, we favor the empirical version of the χ 2 test in the final interpretation, as it indicates the most conservative differences. For the NPIPE dataset, we find that the components and loops differ at more than 95%, but agree within the 99% confidence level with respect to the base model at N = 32, FWHM = 320'. The Euler characteristic at this resolution displays a per mil deviation. In contrast, the FFP10 dataset shows that the observations are consistent with the base model within the 95% confidence level, at this and smaller scales. This is consistent with the observations of the Planck analysis pipeline via Minkowski functionals. For the largest smoothing scale, N = 16, FWHM = 640', both datasets exhibit an anomalous behavior of the loops, where FFP10 data exhibit a deviation that is larger by an order of magnitude than that of the NPIPE dataset. In contrast, the values for the topological components and the Euler characteristic agree between observations and model to within a confidence level of 99%. However, for the largest scales, the statistics are based on low numbers and may have to be regarded with caution. Even though both datasets exhibit mild to significant discrepancies, they also exhibit contrasting behaviors at various instances. Therefore, we do not find it feasible to convincingly accept or reject the null hypothesis. Disregarding the large-scale anomalies that persist at similar scales in WMAP and Planck, observations of the cosmic microwave background are largely consistent with the standard cosmological model within 2σ.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Image segmentation for analyzing galaxy-galaxy strong lensing systems

The goal of this Letter is to develop a machine learning model to analyze the main gravitational lens and detect dark substructure (subhalos) within simulated images of strongly lensed galaxies. Using the technique of image segmentation, we turn the task of identifying subhalos into a classification problem, where we label each pixel in an image as coming from the main lens, a subhalo within a binned mass range, or neither. Our network is only trained on images with a single smooth lens and either zero or one subhalo near the Einstein ring. On an independent test set with lenses with large ellipticities, quadrupole and octopole moments, and for source apparent magnitudes between 17–25, the area of the main lens is recovered accurately. On average, only 1.3% of the true area is missed and 1.2% of the true area is added to another part of the lens. In addition, subhalos as light as 10 8.5 M ⊙ can be detected if they lie in bright pixels along the Einstein ring. Furthermore, the model is able to generalize to new contexts it has not been trained on, such as locating multiple subhalos with varying masses or more than one large smooth lens.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

BeyondPlanck: VII. Bayesian estimation of gain and absolute calibration for cosmic microwave background experiments

We present a Bayesian calibration algorithm for cosmic microwave background (CMB) observations as implemented within the global end-to-end BEYONDPLANCK framework and applied to the Planck Low Frequency Instrument (LFI) data. Following the most recent Planck analysis, we decomposed the full time-dependent gain into a sum of three nearly orthogonal components: one absolute calibration term, common to all detectors, one time-independent term that can vary between detectors, and one time-dependent component that was allowed to vary between one-hour pointing periods. Each term was then sampled conditionally on all other parameters in the global signal model through Gibbs sampling. The absolute calibration is sampled using only the orbital dipole as a reference source, while the two relative gain components were sampled using the full sky signal, including the orbital and Solar CMB dipoles, CMB fluctuations, and foreground contributions. We discuss various aspects of the data that influence gain estimation, including the dipole-polarization quadrupole degeneracy and processing masks. Comparing our solution to previous pipelines, we find good agreement in general, with relative deviations of -0.67% (-0.84%) for 30 GHz, 0.12% (-0.04%) for 44 GHz and -0.03% (-0.64%) for 70 GHz, compared to Planck PR4 and Planck 2018, respectively. We note that the BEYONDPLANCK calibration was performed globally, which results in better inter-frequency consistency than previous estimates. Additionally, WMAP observations were used actively in the BEYONDPLANCK analysis, which both breaks internal degeneracies in the Planck data set and results in an overall better agreement with WMAP. Finally, we used a Wiener filtering approach to smoothing the gain estimates. We show that this method avoids artifacts in the correlated noise maps as a result of oversmoothing the gain solution, which is difficult to avoid with methods like boxcar smoothing, as Wiener filtering by construction maintains a balance between data fidelity and prior knowledge. Although our presentation and algorithm are currently oriented toward LFI processing, the general procedure is fully generalizable to other experiments, as long as the Solar dipole signal is available to be used for calibration.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Evaluation of Light-Element Reactions in the Resolved Resonance Region

Light-element reactions at low energies in the resolved resonance region are important for a range of applications in basic and applied sciences including nuclear reactors, nonproliferation, cultural heritage, forensics and environmental control, rare event investigations and nuclear astrophysics. In this paper, we report on an effort to evaluate charged-particle cross sections in the resolved resonance region and produce evaluated nuclear data files for further processing and inclusion in evaluated data libraries. We discuss the open issues in R-matrix calculations as we extend to higher energies, such as dealing with the rapidly growing number of open channels and merging with the regime of smooth cross sections described by the statistical model, and present attempts to address these issues in neutron-induced reactions relevant to nuclear reactor applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Driven electrostatic phase space vortices in a 1D weakly dissipative Vlasov–Poisson system

In this paper, the effect of collisions on driven electrostatic phase space vortices is analyzed by means of Eulerian simulation for two different collision models. It was demonstrated recently that in the absence of collisions, at late times, steady state phase space vortices manifest to form a plateau in the resonant region of the particle velocity distribution function, due to trapping of particles supporting multiextrema giant phase space vortices (PSVs). In the presence of collisions, over long time, this multiextrema plateau are found to smooth out, since collisions drive the velocity distribution toward Maxwellian, irrespective of how weak the collisions are as long as they are non-zero. In these conditions, kinetic processes and collisionality are found to be in competition, and the evolution of the plasma is found, therefore, to be a result of nontrivial combination of these two effects. An attempt has been made by means of numerical simulations to study the effect of weak collisionality on the electrostatic driven phase space vortices with two types of collision operators: (1) Bhatnagar–Gross–Krook (Krook) collision operator, where the colliding particles can be treated as isolated pairs and, (2) Fokker–Planck (FP) type collision operator (Zakharov–Karpman) in one dimension, where many weak collisions lead to particle diffusion in velocity space. It is shown that depending on the collision model used, the nature of smoothing in velocity space of giant PSVs results in qualitatively very different phase space structures. However, irrespective of the collision model used, excess density fractions over 10% are retained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Generalized Grad–Shafranov equation for non-axisymmetric MHD equilibria

The structure of static MHD equilibria that admit continuous families of Euclidean symmetries is well understood. Such field configurations are governed by the classical Grad–Shafranov equation, which is a single elliptic partial differential equation in two space dimensions. By revealing a hidden symmetry, here we show that in fact all smooth solutions of the equilibrium equations with non-vanishing pressure gradients away from the magnetic axis satisfy a generalization of the Grad–Shafranov equation. In contrast to solutions of the classical Grad–Shafranov equation, solutions of the generalized equation are not automatically equilibria, but instead only satisfy force balance averaged over the one-parameter hidden symmetry. We then explain how the generalized Grad–Shafranov equation can be used to reformulate the problem of finding exact three-dimensional smooth solutions of the equilibrium equations as finding an optimal volume-preserving symmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Enabling complete multichannel nonadiabatic dynamics: A global representation of the two-channel coupled, 1,2 1 A and 1 3 A states of NH 3 using neural networks

Global coupled three-state two-channel potential energy and property/interaction (dipole and spin–orbit coupling) surfaces for the dissociation of NH 3 (Ã) into NH + H 2 and NH 2 + H are reported. The permutational invariant polynomial-neural network approach is used to simultaneously fit and diabatize the electronic Hamiltonian by fitting the energies, energy gradients, and derivative couplings of the two coupled lowest-lying singlet states as well as fitting the energy and energy gradients of the lowest-lying triplet state. The key issue in fitting property matrix elements in the diabatic basis is that the diabatic surfaces must be smooth, that is, the diabatization must remove spikes in the original adiabatic property surfaces attributable to the switch of electronic wavefunctions at the conical intersection seam. Here, we employ the fit potential energy matrix to transform properties in the adiabatic representation to a quasi-diabatic representation and remove the discontinuity near the conical intersection seam. The property matrix elements can then be fit with smooth neural network functions. The coupled potential energy surfaces along with the dipole and spin–orbit coupling surfaces will enable more accurate and complete treatment of optical transitions, as well as nonadiabatic internal conversion and intersystem crossing.

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

The Area Localized Coupled Model for Analytical Mean Flow Prediction in Arbitrary Wind Farm Geometries

This work introduces the area localized coupled (ALC) model, which extends the applicability of approaches that couple classical wake superposition models and atmospheric boundary layer models to wind farms with arbitrary layouts. Coupling wake and top–down boundary layer models is particularly challenging since the latter requires averaging over planform areas associated with turbine-specific regions of the flow that need to be specified. The ALC model uses Voronoi tessellation to define this local area around each turbine. A top–down description of a developing internal boundary layer is then applied over Voronoi cells upstream of each turbine to estimate the local mean velocity profile. Coupling between the velocity at hub-height based on this localized top–down model and a wake model is achieved by enforcing a minimum least-square-error in mean velocity in each cell. The wake model in the present implementation takes into account variations in wind farm inflow velocity and represents the wake profile behind each turbine as a super-Gaussian function that smoothly transitions between a top-hat shape in the region immediately following the turbine to a Gaussian profile downstream. Detailed comparisons to large-eddy simulation (LES) data from two different wind farms demonstrate the efficacy of the model in accurately predicting both wind farm power output and local turbine hub-height velocity for different wind farm geometries. These validations using data generated from two different LES codes demonstrate the model's versatility with respect to capturing results from different simulation setups and wind farm configurations.

49 EE - Wind and Water Power Program - Wind (EE-4W↗
DOE OSTI2021

A simple model of a strong shock driven by a spherical or cylindrical piston

A simple model of piston-driven spherical and cylindrical shocks is suggested. The model is based on a consistent use of two factors: (a) an almost uniform pressure across the shocked layer and, (b) continuous geometrical stretching of the surface elements of the expanding piston. It turns out that for a uniform pre-shock medium the gas between the piston and the shock behaves essentially as an incompressible fluid. An algebraic equation for the shock vs piston position is obtained. Detailed evaluation of the accuracy of the proposed solution shows that its accuracy is a few percent for the adiabatic index γ = 5/3 (as in ideal plasma). A closed-form solution describing enhancement of a weak ambient magnetic field by the shock is presented. The proposed model of piston-driven shocks goes beyond the classical self-similar solutions in that it: (1) naturally covers an early, non-asymptotic dynamics and its transition to asymptotic regime; (2) allows for smooth radial density variation of an ambient gas of the form of bumps, dimples or ramps between two constant values; and (3) allows for smooth temporal variation of piston velocity of the form of bumps, dimples, or ramps. Furthermore, this simple and versatile model provides some new insights into a classical hydrodynamical problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A semi-analytical model for laser induced densification in fused silica at elevated temperatures

A semi-analytical model for deterministic surface smoothing of fused silica below the evaporation temperature threshold is presented. Using this model, a method for calculating deterministic track depths for laser polishing fused silica is derived and validated with experiments. In this work, model predictions were in good agreement with experiment across varying laser preheating temperature and power. With the semi-analytical method presented here, deterministic laser smoothing can be achieved in the densification regime to close the process loop of laser polishing.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Construction of meta-GGA functionals through restoration of exact constraint adherence to regularized SCAN functionals

The strongly constrained and appropriately normed (SCAN) meta-GGA exchange–correlation functional is constructed as a chemical environment-determined interpolation between two separate energy densities: one describes single-orbital electron densities accurately and another describes slowly varying densities accurately. To conserve constraints known for the exact exchange–correlation functional, the derivatives of this interpolation vanish in the slowly varying limit. While theoretically convenient, this choice introduces numerical challenges that degrade the functional’s efficiency. We have recently reported a modification to the SCAN meta-GGA, termed restored-regularized-SCAN (r 2 SCAN) , that introduces two regularizations into SCAN, which improve its numerical performance at the expense of not recovering the fourth order term of the slowly varying density gradient expansion for exchange. Here, we show the derivation of a progression of density functional approximations [regularized SCAN (rSCAN), r++SCAN, r 2 SCAN, and r 4 SCAN] with increasing adherence to exact conditions while maintaining a smooth interpolation. The greater smoothness of r 2 SCAN seems to lead to better general accuracy than the additional exact constraint of SCAN or r 4 SCAN does.

, correlation energy↗
DOE OSTI2022

Modeling atomically mixed graded density impactors

Graded density impactors (GDIs) are multi-material composite impactors used in gas gun experiments to tailor the drive conditions imparted to a sample test material. Previous graded density impactors generally rely on thin, but discrete, layers of different materials. The thinner and the greater number of layers will result in smoother compression. Taken to the limit of very thin layers would be pure material 1 at one surface, such as the front surface of an impactor, smoothly transitioning at the atomic scale to pure material 2 on the back surface. Such an impactor can initially shock, then smoothly compresses a material during a dynamic experiment. This type of experiment can serve to explore a larger region of thermodynamic space than a single or even multi-shock experiments. An overview of how graded density impactors are made is reviewed and sample results are given. Additionally, a strategy for modeling these kinds of impactors is presented. The length scales of constituent mixing are given from the experimental build through electrochemical-deposition. Equation of state models for pure constituents and their subsequent mixtures are presented. It is demonstrated that the time scales for pressure and temperature equilibration, for atomically mixed GDIs, are short enough to be a justifiable closure for the resulting multiphase flow. Furthermore, we present simulation results of dynamic shock followed by a ramp compression, utilizing a silver/gold graded density impactor, onto a tantalum sample.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Dynamic control of the spatial frequency content of an intense laser via intra-beam energy transfer

We present theory and simulations demonstrating that bandwidth-inducing optical smoothing techniques such as smoothing by spectral dispersion (SSD) can initiate intra-beam energy transfer (IBET) within a single beam. IBET consists of energy transfer between the different frequency components of the beam and results in a simultaneous and correlated redistribution of the temporal and spatial frequencies of the laser's electric field. We identify a threshold in bandwidth for IBET to occur and validate the theory with simulations. The bandwidth and the color cycling of SSD can result in IBET in conditions relevant for high-energy density experiments at the National Ignition Facility and Omega Laser. We show that IBET could be used to manipulate a laser's spatial frequency content within a plasma. In particular, this technique could be used to increase a laser's effective f-number, with potential applications for self-guiding in plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗