Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Kernel”

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 109 records · Page 6

A direct numerical simulation of Jet A flame kernel quenching

The safe operation of aeronautical engines requires an understanding of flame ignition, propagation and extinction. In this study, direct numerical simulations are performed using a 29 species reduced chemical mechanism for jet fuel surrogate Jet A to understand the flame quenching process. Here, initially laminar spherical flames of varying sizes and equivalence ratios are subject to an identical periodic domain of decaying and isotropic high intensity turbulence with a turbulent Reynolds number of 2400. All cases become quenched, except for the larger kernel with lower Karlovitz number. An analysis of the flame structure shows broadened preheat zone, flame shortening on the product side, differential species diffusion and partial fuel pyrolysis in the fresh mixture. Two extinction mechanisms are identified arising from flame shortening and high flame stretch. Flame shortening occurs due to turbulence-chemistry interactions that resemble the flame–flame interaction in a laminar counterflow reactant-to-reactant configuration, which contorts and breaks up the ignition kernel. Flame stretch is a local effect that attenuates the heat release rate and causes the flame to retreat towards the product mixtures, similar to what has been observed for reactant-to-product laminar counterflow flames. Chemical explosive mode analysis was also performed to quantify the flame structure and local combustion mode. The diffusion–reaction balance in pinched-off flame islands favors extinction of these smaller structures, while auto-ignition modes are observed within the flame kernel after fresh mixture is engulfed and preheated in the product kernel. Statistics of the density-weighted displacement speed conditional on local combustion mode indicates strong correlation between the local extinction mode and negative displacement speed. The local balance between diffusion and reaction ultimately determines the propensity for local extinction in both laminar and turbulent flames, the extent of which has an impact on global flame propagation.

42 ENGINEERING↗

Simple, low-cost and accurate data-driven geophysical forecasting with learned kernels

Modelling geophysical processes as low-dimensional dynamical systems and regressing their vector field from data is a promising approach for learning emulators of such systems. We show that when the kernel of these emulators is also learned from data (using kernel flows, a variant of cross-validation), then the resulting data-driven models are not only faster than equation-based models but are easier to train than neural networks such as the long short-term memory neural network. In addition, they are also more accurate and predictive than the latter. When trained on geophysical observational data, for example the weekly averaged global sea-surface temperature, considerable gains are also observed by the proposed technique in comparison with classical partial differential equation-based models in terms of forecast computational cost and accuracy. When trained on publicly available re-analysis data for the daily temperature of the North American continent, we observe significant improvements over classical baselines such as climatology and persistence-based forecast techniques. Although our experiments concern specific examples, the proposed approach is general, and our results support the viability of kernel methods (with learned kernels) for interpretable and computationally efficient geophysical forecasting for a large diversity of processes.

58 GEOSCIENCES↗

Kernels of the Direct Collision Integral of the Nonlinear Boltzmann Equation

This paper reviews the properties of kernels of the direct collision integral describing the escape of particles from an element of the phase volume due to collisions, as well as properties of kernels of the integral operators arising in expansion of the direct collision integral in spherical harmonics. It is shown that, for isotropic scattering cross sections depending on the relative velocity by a power law, these kernels can be obtained in an analytical form. Since the analytical expressions for the kernels with large indices become too cumbersome, their asymptotics are found.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Smoothing Lexis diagrams using kernel functions: A contemporary approach

Lexis diagrams are rectangular arrays of event rates indexed by age and period. Analysis of Lexis diagrams is a cornerstone of cancer surveillance research. Typically, population-based descriptive studies analyze multiple Lexis diagrams defined by sex, tumor characteristics, race/ethnicity, geographic region, etc. Inevitably the amount of information per Lexis diminishes with increasing stratification. Several methods have been proposed to smooth observed Lexis diagrams up front to clarify salient patterns and improve summary estimates of averages, gradients, and trends. In this article, we develop a novel bivariate kernel-based smoother that incorporates two key innovations. First, for any given kernel, we calculate its singular values decomposition, and select an optimal truncation point—the number of leading singular vectors to retain—based on the bias-corrected Akaike information criterion. Second, we model-average over a panel of candidate kernels with diverse shapes and bandwidths. The truncated model averaging approach is fast, automatic, has excellent performance, and provides a variance-covariance matrix that takes model selection into account. We present an in-depth case study (invasive estrogen receptor-negative breast cancer incidence among non-Hispanic white women in the United States) and simulate operating characteristics for 20 representative cancers. The truncated model averaging approach consistently outperforms any fixed kernel. Our results support the routine use of the truncated model averaging approach in descriptive studies of cancer.

60 APPLIED LIFE SCIENCES↗

TEMPORAL EVOLUTION AND SPATIAL DISTRIBUTION OF WHITE-LIGHT FLARE KERNELS IN A SOLAR FLARE

On 2011 September 6, we observed an X2.1-class flare in continuum and H α with a frame rate of about 30 Hz. After processing images of the event by using a speckle-masking image reconstruction, we identified white-light (WL) flare ribbons on opposite sides of the magnetic neutral line. We derive the light curve decay times of the WL flare kernels at each resolution element by assuming that the kernels consist of one or two components that decay exponentially, starting from the peak time. As a result, 42% of the pixels have two decay-time components with average decay times of 15.6 and 587 s, whereas the average decay time is 254 s for WL kernels with only one decay-time component. The peak intensities of the shorter decay-time component exhibit good spatial correlation with the WL intensity, whereas the peak intensities of the long decay-time components tend to be larger in the early phase of the flare at the inner part of the flare ribbons, close to the magnetic neutral line. The average intensity of the longer decay-time components is 1.78 times higher than that of the shorter decay-time components. If the shorter decay time is determined by either the chromospheric cooling time or the nonthermal ionization timescale and the longer decay time is attributed to the coronal cooling time, this result suggests that WL sources from both regions appear in 42% of the WL kernels and that WL emission of the coronal origin is sometimes stronger than that of chromospheric origin.

79 ASTRONOMY AND ASTROPHYSICS↗

StOKeDMD: Streaming Occupation kernel dynamic mode decomposition

Dynamic mode decomposition (DMD) has become a common technique for constructing surrogate models for dynamical systems from observed system states. The Occupation Kernel DMD (OKDMD) method proposed in (Rosenfeld et al., 2022) and (Rosenfeld et al., 2024) is a Liouville operator based method that builds surrogate models from system state trajectories. Here, this paper proposes an extension of OKDMD to the case when the system states are observed in a streaming fashion, i.e., only a small fraction of the state trajectory is available at a given time. The developed method, Streaming Occupation Kernel DMD (StOKeDMD), accommodates the streaming data input by leveraging properties of specific choices of kernel functions and occupation kernels. We apply the StoKeDMD method as a compression method for streaming data, analyze the memory complexity, and demonstrate the performance of StoKeDMD in the compression of streaming data generated from a Lorenz system and a fluid flow simulation.

97 MATHEMATICS AND COMPUTING↗

Leveraging a Neural Network-Enhanced Reproducing Kernel Particle Method for Multiphysics Degradation Modeling of Energy Storage Materials

Energy storage materials exhibit strong electro-chemo-mechanical coupling and highly anisotropic material properties, contributing to the formation and propagation of micro-cracking during charge/discharge cycling and resulting in reduced performance and service life. A coupled electro-chemo-mechanical reproducing kernel particle method (RKPM) formulation has been developed to analyze this system. With microstructural images supplied by the National Renewable Energy Laboratory (NREL), pixel-based model construction by RKPM is used to represent the complex material microstructures that dictate the coupled physics of these systems. Traditional electro-chemo-mechanical models rely on mesh-based finite element methods, which can lead to difficulties in meshing such complex geometries and capturing crack propagation due to mesh dependency. Here, a neural network-enhanced reproducing kernel particle method (NN-RKPM) [1, 2] is introduced to effectively model damage and crack propagation in the material microstructures; the location, orientation, and solution transition near a localization are automatically captured by superimposed block-level NN optimizations. This NN enrichment approach allows for effective modeling of localizations via a fixed background discretization, relieving tedious efforts for adaptive refinement in traditional mesh-based methods. Applications to the heterogeneous microstructures of Li-ion battery cathodes will be presented to demonstrate the effectiveness of the proposed methods. NN-RKPM is additionally used to inform how crack opening and closure in turn affect the coupled chemical equations and material microstructure. Reference: [1] Baek, J., Chen, J. S., Susuki, K., "Neural Network enhanced Reproducing Kernel Particle Method for Modeling Localizations," International Journal for Numerical Methods in Engineering, Vol. 123, pp 4422-4454, https://doi.org/10.1002/nme.7040, 2022. [2] Baek, J., Chen, J. S., "A Neural Network-Based Enrichment of Reproducing Kernel Approximation for Modeling Brittle Fracture", Computer Methods in Applied Mechanics and Engineering Vol. 410, 116590, 2024.

degradation↗

The effective number of parameters in kernel density estimation

We devise a new formula for measuring the effective degrees of freedom (EDoF) in kernel density estimation (KDE). Starting from the orthogonal polynomial sequence (OPS) expansion for the ratio of the empirical to the oracle density, we show how convolution with the kernel leads to a new OPS with respect to which one may express the resulting KDE. The expansion coefficients of the two OPS systems can then be related via a kernel sensitivity matrix, which leads to a natural oracle definition of EDoF through the trace operator. Asymptotic properties of the (empirical) plug-in EDoF are worked out through influence functions, and connections with other empirical EDoFs are established. Minimization of Kullback-Leibler divergence is investigated as an alternative to integrated squared error based bandwidth selection rules, yielding a new normal scale rule. The methodology, which arises from a proper oracle formulation and is not restricted to convolution kernels, suggests the possibility of a new bandwidth selection rule based on an information criterion such as AIC.

bandwidth selection↗

Microstructure and fission products in the UCO kernel of an AGR-1 TRISO fuel particle after post irradiation safety testing

Tristructural isotropic (TRISO) coated particle fuel as a key fuel concept for the High Temperature Gas-cooled Reactors (HTGRs) and a candidate accident tolerant fuel (ATF) has been currently investigated under the US-DOE Advanced Gas Reactor (AGR) fuel development and qualification program. Over the past decades, extensive studies have been conducted to evaluate the fission product release, diffusion on Ag, Pd and Cd in the SiC layer and the TRISO coating system, and post-safety-test performance. However, to date, there are limited reported results on the fuel kernels’ response to irradiation with or without post irradiation high temperature annealing. To incrementally fulfill this knowledge gap extensive studies using transmission electron microscopes (TEM) and atom probe tomography (APT) were conducted on a TRISO fuel particle kernel with a 19.74% 235U enrichment, irradiated to 18.63% FIMA and subsequently subjected to post-safety testing at 1600ºC for 300 hours.. Microstructure characterizations, elemental analysis, and phase identification were conducted using conventional TEM and scanning TEM (STEM) imaging, energy-dispersive X-ray spectroscopy (EDS), selected-area electron diffraction (SAED) and atom probe tomography (APT). The following findings were made: (1) significant reconstructions and phase evolutions occurred in the irradiated and post-safety annealed fuel kernel, and its microstructure consists of two primary phases, namely a “high-Z (atomic mass)” UC phase and a “low-Z” UO phase, (2) no fission gas bubbles are identified within the fuel kernel, that can be attributed to the high temperature post safety annealing, (3) fission products Zr, Nb, Mo, Ru, Tc and Rh were found to preferentially segregate into UC phase or to form metallic precipitates, while the lanthanide fission products tend to stay in the solution of UO phase, (4) Pd was detected in the rod-shaped precipitates, (5) APT data reveals segregation of Te, Pd, and Ag at the UO and UC phase boundary, and (6) Pd is detected in both the UO and UC phases but with an enrichment in the former one.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Physicochemical evolution of uranium nitride kernel microstructure with varying carbon distribution for advanced TRISO fuel forms

Uranium nitride (UN) has emerged as a fuel candidate for advanced nuclear reactor concepts due to its superior uranium density, thermal conductivity, and high melting temperature. However, the fabrication route for converting UO 2 to UN is complex and difficult to standardize. Although the chemistry of this conversion process is well-studied, more insight into the physicochemical dynamics of this conversion using advanced characterization techniques can help further our understanding of this material system. This work leveraged thermogravimetric analysis (TGA), X-ray diffraction (XRD), and nondestructive 3D X-ray computed tomography (XCT) to characterize dynamic microstructural changes in the UO 2 → UCO → UN fabrication pathway for two kernels with a varying carbon distribution in the starting composition. TGA and XRD were used to quantify changes in the mass, density, and chemical composition of the two kernels, while three-dimensional image processing and segmentation of XCT data were used to quantify the volume, surface area, and spatial distribution of features within each kernel for multiple steps along the fabrication pathway. The analysis indicates distinct differences between the two kernels that are correlated to downstream conversion efficiency. In conclusion, this work is among the first to perform 3D quantification of physicochemical evolution during UN conversion, providing quantitative correlation between processing, properties, and expected fuel performance.

Nuclear fuel↗

Transfer learning of memory kernels for transferable coarse-graining of polymer dynamics

The present work concerns the transferability of coarse-grained (CG) modeling in reproducing the dynamic properties of the reference atomistic systems across a range of parameters. In particular, we focus on implicit-solvent CG modeling of polymer solutions. The CG model is based on the generalized Langevin equation, where the memory kernel plays the critical role in determining the dynamics in all time scales. Thus, we propose methods for transfer learning of memory kernels. The key ingredient of our methods is Gaussian process regression. By integration with the model order reduction via proper orthogonal decomposition and the active learning technique, the transfer learning can be practically efficient and requires minimum training data. Through two example polymer solution systems, we demonstrate the accuracy and efficiency of the proposed transfer learning methods in the construction of transferable memory kernels. The transferability allows for out-of-sample predictions, even in the extrapolated domain of parameters. Built on the transferable memory kernels, the CG models can reproduce the dynamic properties of polymers in all time scales at different thermodynamic conditions (such as temperature and solvent viscosity) and for different systems with varying concentrations and lengths of polymers.

Ma, Zhan↗

On the memory kernel and the reduced system propagator

We relate the memory kernel in the Nakajima–Zwanzig–Mori time-convolution approach to the reduced system propagator which is often used to obtain the kernel in the Tokuyama–Mori time-convolutionless approach. The connection provides a robust and simple formalism to compute the memory kernel for a generalized system-bath model circumventing the need to compute high order system-bath observables, thus streamlining the use of numerically exact solvers for calculating the memory kernel. We illustrate this for a model system with electron-electron and electron-phonon couplings, driven away from equilibrium.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Efficient construction of generalized master equation memory kernels for multi-state systems from nonadiabatic quantum-classical dynamics

Techniques derived from the generalized quantum master equation (GQME) framework have provided the basis for elucidating energy and charge transfer in systems ranging from molecular solids to photosynthetic complexes. Recently, the nonperturbative combination of the GQME with quantum-classical methods has resulted in approaches whose accuracy and efficiency exceed those of the original quantumclassical schemes while offering significant accuracy improvements over perturbative expansions of the GQME. In this work, we show that, while the non-Markovian memory kernel required to propagate the GQME scales quartically with the number of subsystem states, the number of trajectories required scales at most quadratically when using quantum-classical methods to construct the kernel. We then present an algorithm that allows further acceleration of the quantum-classical GQME by offering a way to selectively sample the kernel matrix elements that are most important to the process of interest. We demonstrate the utility of these advances by applying the combination of Ehrenfest mean field theory with the GQME (MF-GQME) to models of the Fenna-Matthews-Olson (FMO) complex and the light harvesting complex II (LHCII), with 7 and 14 states, respectively. This allows us to show that the MF-GQME is able to accurately capture all the relevant dynamical time scales in LHCII: the initial nonequilibrium population transfer on the femtosecond time scale, the steady state-type trapping on the picosecond time scale, and the long time population relaxation. Remarkably, all of these physical effects spanning tens of picoseconds can be encoded in a memory kernel that decays only after ~65 fs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Pyrolytic Carbon Coating Effects on Oxide and Carbide Kernels Intended for Nuclear Fuel Applications

The coating of nuclear fuel kernels with pyrolytic carbon (PyC) is a well-understood practice dating back over half a century. In spite of decades of studies related to these coatings, no study has yet investigated the effect of the PyC deposition coating process on the kernels themselves. In this study, the composition and crystallographic phase of kernel materials were observed to change after exposure to the thermal and chemical environment of PyC coating process. Specifically, the coating process increased the fraction of high carbon content phase within carbide microsphere kernels, with W 2 C containing microspheres driven towards WC, and UC containing microspheres driven towards UC 2 . Oxide microspheres consisted of a mixture of two crystalline phases. The monoclinic phase within yttria-stabilized zirconia microspheres was eliminated by the coating process resulting in a purely tetragonal phase. Hafnium oxide microspheres were more stable showing no detectable change in composition or crystal structure after coating.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Hot Hydrogen Testing of W-Coated UN Kernels in a Mo30W Matrix

Ceramic uranium mononitride (UN) is being considered as a reactor fuel for nuclear thermal propulsion. To avoid or reduce the dissociation of UN at the high temperatures needed, embedding it in a metallic matrix (cermet) has been proposed. To assess the viability of this concept, hot hydrogen testing of tungsten-coated UN kernels embedded in a Mo-30 wt% W (Mo30W) alloy matrix has been performed at temperatures from 1800°C to 2300°C. Both the isolated kernels and kernels consolidated by spark plasma sintering in the Mo30W matrix were tested. In addition to direct observations and mass loss measurements, the samples were analyzed by X-ray diffraction (XRD) and scanning electron microscopy (SEM)/energy dispersive X-ray spectroscopy (EDS) after each run. The decomposition of UN started at 1800°C despite the coating and matrix, and increased at 2000°C. Uranium seeped through the tungsten grain boundaries of the coating at all temperatures. The consolidated sample expanded irregularly at 2000°C through the formation of voids, and SEM/EDS analysis showed uranium-containing veins in the matrix consisting of U 2 Mo according to the XRD data. The observed pore generation at 2000°C was explained by the formation of water vapor from residual oxides and diffused hydrogen. At 2200°C and above, both the kernels and the consolidated samples melted through the formation of uranium or low–melting point uranium-molybdenum alloys.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Collins-Soper kernel in the QCD instanton vacuum

We outline a general framework for evaluating the nonperturbative soft functions in the quantum chromodynamics (QCD) instanton vacuum. In particular, from the soft function we derive the Collins-Soper (CS) kernel, which drives the rapidity evolution of the transverse-momentum-dependent parton distributions. The resulting CS kernel, when supplemented with the perturbative contribution, agrees well with recent lattice results and some phenomenological parametrizations. Moreover, our CS kernel depends logarithmically on the large quark transverse separation, providing a key constraint on its phenomenological parametrization. Finally, a lattice calculation can be directly compared to our generic results in Euclidean signature, thus providing a new approach for evalulating the soft function and extracting the CS kernel by analytical continuation.

QCD phenomenology↗

MiniApp for Density Matrix Renormalization Group Hamiltonian Application Kernel

We present two miniapps that implement the core computational kernel of the DMRG++ application, a generic C++ code that implements the Density Matrix Renormalization Group (DMRG) algorithm. The DMRG++ core Kronecker multiplication kernel is formulated using a batched BLAS approach, with implementation that targets both multi-core CPUs using OpenMP and GPGPU using the MAGMA library. The kernel evaluates the matrix-vector multiplication of the target Hamiltonian matrix used in Lanczos algorithm for computing the lowest eigenvalue and eigenvector. The Hamiltonian matrix is expressed compactly as sums of Kronecker products of small dense matrices. We demonstrate improved performance of the miniapp on synthetic problem, and show the performance of the DMRG++ application using a plugin based on the miniapp. We also present an OpenMP miniapp that explores the use of nested parallel constructs to implement the Kronecker multiplication kernel, exploring the use of nested OpenMP worksharing and tasking abstractions to implement the multi-level parallel multiplication algorithm. The miniapp has been used as a co-design vehicle for evaluating features in the OpenMP-4.5 and upcoming OpenMP-5.0 standards.

Elwasif, Wael↗

Scalable Extended Dynamic Mode Decomposition Using Random Kernel Approximation

The Koopman operator is a linear, infinite-dimensional operator that governs the dynamics of system observables; extended dynamic mode decomposition (EDMD) is a data-driven method for approximating the Koopman operator using functions (features) of the system state snapshots. This paper investigates an approach to EDMD in which the features used provide random approximations to a particular kernel function. The objective of this is computational economy for large data sets: EDMD is generally ill-suited for problems with large state dimension, and its dual kernel formulation (KDMD) is well-suited for such problems only if the number of data snapshots is relatively small. We discuss two specific methods for generating features: random Fourier features and the Nÿstrom method. The first method is a data-independent method for translation-invariant kernels only and involves random sampling in feature space; the second method is a data-dependent empirical method that may be used for any kernel and involves random sampling of data. We first discuss how these ideas may be applied in an EDMD context, as well as a means for adaptively adding random Fourier features. As a result, we demonstrate these methods on two example problems and conclude with an analysis of the relative benefits and drawbacks of each method.

59 BASIC BIOLOGICAL SCIENCES↗