SEARCH · Engineering Papers
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.
A multiphysics model is presented to study potential losses in Cu-MEAs and how various physical phenomena impact the product distribution of CO 2 reduction.
ABSTRACT Methane-producing archaea are key organisms in the anaerobic carbon cycle. These organisms, also called methanogens, grow by converting substrate to methane gas in a process called methanogenesis. Previous research showed that the reduction of the terminal electron acceptor is the rate-limiting step in methanogenesis by Methanosarcina acetivorans . In order to gain insight into how the cells sense and respond to the availability of the terminal electron acceptor, we designed an experiment to deplete cells of the essential terminal oxidase enzyme, HdrED. We found that the depletion of HdrED in vivo results in a higher abundance of transcripts for methyltransferases ( mtaC2, mtaB3, mtaC3 ), coenzyme B biosynthesis, C1 metabolism, and pyrimidine compounds. In most cases, these changes were distinct from transcript abundance changes observed during the transition from exponential growth to stationary phase cultures. These data implicate the methylotrophic methanogenesis regulator MsrC (MA4383) in CoM-S-S-CoB heterodisulfide sensing and indicate cells have a specific mechanism to sense intracellular ratio of CoM-S-S-CoB, coenzyme M, and coenzyme B thiols and further suggest transcripts encoding translation and methanogenesis functions are controlled by feed-forward regulation depending on substrate availability. IMPORTANCE Methanosarcina is an emerging model archaeon and synthetic biology platform for the production of renewable energy and sustainable chemicals to reduce dependence on petroleum. Research into metabolic networks and gene regulation in this organism and other methanogens will inform genome-scale metabolic modeling and microbial function prediction in uncultured or non-model anaerobes and archaea. This study suggests methanogens use unknown mechanisms to efficiently couple methanogenesis to gene regulation via CoM-S-S-CoB and ATP availability.
In this work we introduce a manifold learning-based method for uncertainty quantification (UQ) in systems describing complex spatiotemporal processes. Our first objective is to identify the embedding of a set of high-dimensional data representing quantities of interest of the computational or analytical model. For this purpose, we employ Grassmannian diffusion maps, a two-step nonlinear dimension reduction technique which allows us to reduce the dimensionality of the data and identify meaningful geometric descriptions in a parsimonious and inexpensive manner. Polynomial chaos expansion is then used to construct a mapping between the stochastic input parameters and the diffusion coordinates of the reduced space. An adaptive clustering technique is proposed to identify an optimal number of clusters of points in the latent space. The similarity of points allows us to construct a number of geometric harmonic emulators which are finally utilized as a set of inexpensive pretrained models to perform an inverse map of realizations of latent features to the ambient space and thus perform accurate out-of-sample predictions. Thus, the proposed method acts as an encoder-decoder system which is able to automatically handle very high-dimensional data while simultaneously operating successfully in the small-data regime. The method is demonstrated on two benchmark problems and on a system of advection-diffusion-reaction equations which model a first-order chemical reaction between two species. In all test cases, the proposed method is able to achieve highly accurate approximations which ultimately lead to the significant acceleration of UQ tasks.
Turbulent flows over a large surface area (S) covered by n obstacles experience an overall drag due to the presence of the ground and the protruding obstacles into the flow. The drag partition between the roughness obstacles and the ground is analyzed using an analytical model proposed by Raupach (Boundary-Layer Meteorol 60:375-395, 1992) and is hereafter referred to as R92. The R92 is based on the premise that the wake behind an isolated roughness element can be described by a shelter area A and a shelter volume V. The individual sizes of A and V without any interference from other obstacles can be determined from scaling analysis for the spread of wakes. To upscale from an individual roughness element to n/S elements where wakes may interact, R92 adopted a background stress re-normalizing instead of reducing A or V with each element addition. This work demonstrates that R92’s approach results in a linear background stress reduction in A and V only when the ratio of n/S is small, due to a low probability of wake interactions. This probabilistic nature suggests that up-scaling from individual to multiple roughness elements can be re-formulated using stochastic averaging methods proposed here. The two approaches are shown to recover R92 under plausible conditions. An alternative scaling for the shelter volume is also proposed here using thermodynamic arguments of work and dissipation though the final outcome remains similar to R92. Comparisons between R92 and available data spanning more than two decades after R92 on blocks and vegetation-like roughness elements confirm the practical utility of R92. The agreement between R92 and this updated databases of experiments and simulations confirm the potential use of R92 in large-scale models provided that the relevant parameters accommodate certain features of the roughness element type (cube versus vegetation-like) and, to a lesser extent, their configuration throughout S. Last, a comparison between R92 and models based on first-order closure principles with constant mixing length suggests that R92 can outperform such models when evaluated across a wide range of roughness densities.
Oceans cover two-thirds of the Earth and understanding the interactions of aerosols with clouds in these large marine regions requires quantifying the man-made contributions to the budget of cloud-drop forming particles (known as cloud condensation nuclei, or CCN) relative to the non-manmade “baseline” conditions and understanding the meteorology of boundary layer clouds. Modeling studies have shown substantial uncertainties and sensitivities to natural marine CCN sources, meaning that to reduce uncertainties in indirect effects we must be able to better quantify the CCN budget in ocean regions. While models provide important constraints on these uncertainties, actually reducing uncertainties requires substantial observations in open-ocean and coastal regions in order to establish the baseline on which manmade emissions are added. The tropical South Atlantic Ocean is one of the least-sampled regions of the planet, making the comprehensive measurements of the Department of Energy Atmospheric Radiation Measurement Layered Atlantic Smoke Interactions with Clouds (LASIC) campaign provides the longest record of aerosol size distribution measurements from a differential mobility analyzer in a cloud-influenced marine location in the ARM database, with 17 months of ground-based aerosol size distribution measurements. This project addressed the research topic of Aerosol-Cloud Interactions (ACI) by supporting three publications from the LASIC measurements: 1. The ARM measurements were combined with a new value-added technique that was pioneered by the Russell group for quantifying sea salt. This quantification of sea salt uses supermicron scattering measurements for retrieving reasonable sea spray mass concentrations, providing the best-available, observationally-constrained estimate of the sea spray mode properties when supermicron size distribution measurements are not available. 2. The fitted modes of the distribution were used to investigate the signatures of cloud processing for very clean to very smoky aerosol conditions, revealing not only differences in the particles that activate in clouds but also in the mechanisms that control that droplet formation process. In clean air, the size required to form a cloud droplet is influenced by the number of particles, as well as how quickly particles take up water during growth in cloud. 3. Building on this aerosol characterization, cloud and meteorological conditions were used to evaluate aerosol-related changes in cloud albedo and optical depth. We introduced a new method of decomposing the impact of aerosols on clouds known as the Twomey effect by incorporating retrieved supersaturation from two independent sets of observations to constrain the feedback of aerosol particles on cloud properties. The method quantifies the reduction of the Twomey effect at high aerosol concentrations, which has never been explained quantitatively by observations. In addition, the results provide the first direct validation for the conditions observed in the tropical South Atlantic of a parcel-based approach that is embedded in many climate models. Together these findings illustrate how ARM extended field campaigns in stratocumulus-covered regions can be used to constrain ACI processes with direct observations, providing specific radiative effects without models. By making such process-specific constraints available to improve ACI in global climate models, ARM observations play a key role in supporting model development.
We report that recently, there has been a significant interest in detonation-based combustion systems, such as rotating detonation engines (RDE), due to potential performance advantages in propulsion and energy applications. This includes improvements in thermodynamic efficiency, inherently high gas velocities, and the possibility of achieving an increase in total pressure (i.e. pressure gain) through the combustion process. Due to approximating a constant volume combustion process, these devices also tend to exhibit extremely high local gas temperatures relative to a comparable constant pressure combustor. These advantages overlap well with desirable performance characteristics of direct power extraction (DPE) technologies, such as a magnetohydrodynamic (MHD) generator. Typically, in DPE systems hot combustion products are seeded with an easily ionizable material such as potassium carbonate (K 2 CO 3 ) in order to boost the electrical conductivity. However, due to the short gas residence times within an RDE, it was unclear whether forming an electrically conductive combustion plasma would be feasible for integration with a downstream MHD generator. A model is presented which describes the heating, decomposition, and ionization of solid K 2 CO 3 particles and aqueous solutions of K 2 CO 3 in water, for a given initial particle diameter. This model was combined with available computational fluid dynamics (CFD) data for an oxygen-methane RDE in a one-way coupled Eulerian-Lagrangian framework to predict particle trajectories and the corresponding heating, decomposition, and ionization histories. Electrical conductivities were computed using a previously developed model, and a method was proposed to determine an equivalent average electrical conductivity. Results show that particle sizes below ~30 µm are able to fully decompose before reaching the exit of the RDE. While the one-way coupled nature of the simulations precluded rigorous evaluation of the effects of seed material on detonability, a substantial temperature reduction is expected at the detonation wave front due to heating and decomposition. A preliminary comparison is presented between the RDE and an equivalent constant pressure adiabatic combustor, showing a potential performance advantage for the RDE.
Experiments on ASDEX Upgrade (AUG) in 2021 and 2022 have addressed a number of critical issues for ITER and EU DEMO. A major objective of the AUG programme is to shed light on the underlying physics of confinement, stability, and plasma exhaust in order to allow reliable extrapolation of results obtained on present day machines to these reactor-grade devices. Concerning pedestal physics, the mitigation of edge localised modes (ELMs) using resonant magnetic perturbations (RMPs) was found to be consistent with a reduction of the linear peeling-ballooning stability threshold due to the helical deformation of the plasma. Conversely, ELM suppression by RMPs is ascribed to an increased pedestal transport that keeps the plasma away from this boundary. Candidates for this increased transport are locally enhanced turbulence and a locked magnetic island in the pedestal. The enhanced D-alpha (EDA) and quasi-continuous exhaust (QCE) regimes have been established as promising ELM-free scenarios. Here, the pressure gradient at the foot of the H-mode pedestal is reduced by a quasi-coherent mode, consistent with violation of the high-n ballooning mode stability limit there. This is suggestive that the EDA and QCE regimes have a common underlying physics origin. In the area of transport physics, full radius models for both L- and H-modes have been developed. These models predict energy confinement in AUG better than the commonly used global scaling laws, representing a large step towards the goal of predictive capability. A new momentum transport analysis framework has been developed that provides access to the intrinsic torque in the plasma core. In the field of exhaust, the X-Point Radiator (XPR), a cold and dense plasma region on closed flux surfaces close to the X-point, was described by an analytical model that provides an understanding of its formation as well as its stability, i.e., the conditions under which it transitions into a deleterious MARFE with the potential to result in a disruptive termination. With the XPR close to the divertor target, a new detached divertor concept, the compact radiative divertor, was developed. Here, the exhaust power is radiated before reaching the target, allowing close proximity of the X-point to the target. No limitations by the shallow field line angle due to the large flux expansion were observed, and sufficient compression of neutral density was demonstrated. With respect to the pumping of non-recycling impurities, the divertor enrichment was found to mainly depend on the ionisation energy of the impurity under consideration. In the area of MHD physics, analysis of the hot plasma core motion in sawtooth crashes showed good agreement with nonlinear 2-fluid simulations. This indicates that the fast reconnection observed in these events is adequately described including the pressure gradient and the electron inertia in the parallel Ohm’s law. Concerning disruption physics, a shattered pellet injection system was installed in collaboration with the ITER International Organisation. Thanks to the ability to vary the shard size distribution independently of the injection velocity, as well as its impurity admixture, it was possible to tailor the current quench rate, which is an important requirement for future large devices such as ITER. Progress was also made modelling the force reduction of VDEs induced by massive gas injection on AUG. The H-mode density limit was characterised in terms of safe operational space with a newly developed active feedback control method that allowed the stability boundary to be probed several times within a single discharge without inducing a disruptive termination. Regarding integrated operation scenarios, the role of density peaking in the confinement of the ITER baseline scenario (high plasma current) was clarified. The usual energy confinement scaling ITER98(p,y) does not capture this effect, but the more recent H20 scaling does, highlighting again the importance of developing adequate physics based models. Advanced tokamak scenarios, aiming at large non-inductive current fraction due to non-standard profiles of the safety factor in combination with high normalised plasma pressure were studied with a focus on their access conditions. A method to guide the approach of the targeted safety factor profiles was developed, and the conditions for achieving good confinement were clarified. Based on this, two types of advanced scenarios (‘hybrid’ and ‘elevated’
The effective weight reduction in the automotive industry by the wide adoption of lightweight magnesium (Mg) alloys demands high-quality joint between magnesium alloys and massively-used steels in order to wring the excess weight with strength and safety assurance. However, Mg-steel joint is difficult to achieve because there is no mutual solubility between magnesium and steel and huge disparity in physical properties. An impact-based welding method recently showed successful Mg-steel joining. In this work, the characteristics of Mg-steel interface joined by the impact welding method were investigated. Synchrotron high-energy X-ray computed tomography and diffraction were applied to characterize the microstructure across Mg-steel interface. Results revealed a deposit layer formed at the joint interface where Fe-rich particles spread deep into the Mg matrix. High-resolution 3D morphology of Mg-steel interface demonstrated the trapped pores and cracks inside the deposit layer. Finally, the formation of the deposit layer and the void/cracking evolution were analyzed by using finite element models. These findings provide insights into the immiscible Mg-steel joining process.
This work presents the validation of the local structural load modeling capability in OpenFAST for floating substructures based on data from the FOCAL experimental campaign. Previously, OpenFAST could only represent the floating substructure as a rigid body, and though this approach can model the global response of the floater in most cases, it is not able to capture the structural loads within the floater's individual members. Consideration of local substructure loads is important for some floating designs, because the pursuit of cost reduction often results in lighter and more flexible structures. To address this limitation, the HydroDyn (hydrodynamics) and SubDyn (substructure dynamics) modules of OpenFAST have been recently extended to account for the flexibility of floating substructures. To validate this new capability, we compare the results obtained by OpenFAST with data measured during the FOCAL experimental campaign, which analyzed a 1:70 scale performance-matched model of the IEA 15-MW reference turbine atop a modified University of Maine VolturnUS-S semisubmersible in a wave basin under the action of both wind and waves. For the purposes of the present work, the most important feature of the experiment is the presence of load cells at the root of each pontoon, and our objective is to assess how well those loads are reproduced by OpenFAST. To model the distributed hydrodynamic and hydrostatic loads along the floating substructure, we adopt a strip-theory approach based on the Morison equation, and we discuss the impact of different hydrodynamic modeling options (wave stretching, MacCamy-Fuchs correction, and second-order wave kinematics) on both motions and loads. For simplicity, we focus on wave-only conditions, both regular and irregular. The results demonstrate good overall agreement for the loads at the root of the pontoons for the waves analyzed in this work, especially given the assumptions and simplifications inherent to a simple strip-theory model.
In this paper, we study real-time classical matrix mechanics of a simplified 2 × 2 matrix model inspired by the black hole evaporation problem. This is a step towards making a quantitative model of real-time evaporation of a black hole, which is realized as a bound state of D0-branes in string theory. The model we study is the reduction of Yang–Mills in 2 + 1 dimension to 0 + 1 dimensions, which has been corrected with an additional potential that can be interpreted as a zero-point energy for fermions. Our goal is to understand the lifetime of such a classical bound state object in the classical regime. To do so, we pay particular attention to when D-particles separate to check that the “off-diagonal modes” of the matrices become adiabatic and use that information to improve on existing models of evaporation. It turns out that the naive expectation value of the lifetime with the fermionic correction is infinite. This is a logarithmic divergence that arises from very large excursions in the separation between the branes near the threshold for classical evaporation. The adiabatic behavior lets us get some analytic control of the dynamics in this regime to get this estimate. This divergence is cutoff in the quantum theory due to quantization of the adiabatic parameter, resulting in a long lifetime of the bound state, with a parametric dependence of order log(1/ℏ).
This work presents the validation of the local structural load modeling capability in OpenFAST for floating substructures based on data from the FOCAL experimental campaign. Previously, OpenFAST could only represent the floating substructure as a rigid body, and though this approach can model the global response of the floater in most cases, it is not able to capture the structural loads within its individual members. Consideration of local substructure loads is important for some floating designs, as the pursuit of cost reduction often results in lighter and more flexible structures. To address this limitation, the HydroDyn (hydrodynamics) and SubDyn (substructure dynamics) modules of OpenFAST have been recently extended to account for the flexibility of floating substructures. To validate this new capability, we compare the results obtained by OpenFAST with data measured during the FOCAL experimental campaign, which analyzed a 1:70 scale performance-matched model of the IEA 15-MW reference turbine atop a modified University of Maine VolturnUS-S semisubmersible in a wave basin under the action of both wind and waves. For the purposes of the present work, the most important feature of the experiment is the presence of load cells at the root of each pontoon, and our objective is to assess how well those loads are reproduced by OpenFAST. To model the distributed hydrodynamic and hydrostatic loads along the floating substructure, we adopt a strip-theory approach based on the Morison equation, and we discuss the impact of different hydrodynamic modeling options (wave stretching, MacCamy-Fuchs correction, and second-order wave kinematics) on both motions and loads. For simplicity, we focus on wave-only conditions, both regular and irregular. The results demonstrate good overall agreement for the loads at the root of the pontoons for the waves analyzed in this work, especially given the assumptions and simplifications inherent to a simple strip-theory model.
Spectroscopic measurements in model 1D correlated systems offer insights for understanding their two-dimensional counterparts, which include the cuprate and iron pnictide/chalcogenide superconductors. A major challenge is the identification of such correlated systems with dominantly 1D physics. Here, in this Letter, inelastic neutron scattering measurements on NaFe 0.53 Cu 0.47 As single crystal directly reveal quasi-1D spin excitations, resulting from atomic order that leads to magnetic Fe and nonmagnetic Cu chains. The dominant exchange interaction is antiferromagnetic along the chain [𝑆𝐽 ∥ ≈ 90.1(3) meV], whereas the inter-chain couplings are much weaker [𝑆𝐽 ⊥ ≈ −2.4(1) meV and 𝑆𝐽 c ≈ 0.15(5) meV]. The quasi-1D spin excitations in NaFe 0.53 Cu 0.47 As stem from both the Néel and stripe vectors, with Néel excitations sensitive to Fe impurities on the Cu site. The spin excitations in quasi-1D NaFe 0.53 Cu 0.47 As and quasi-2D FeSe exhibit a striking resemblance, suggesting a common origin for their coexistent stripe and Néel excitations. Our findings demonstrate magnetic dilution in NaFeAs leads to dimension reduction of its magnetic degree of freedom, presenting a strategy for discovering low-dimensional quantum materials.
Patterned functional materials offer improved properties (electrical, thermal, etc.) over their bulk counterparts in many applications, including energy storage, flexible electronics, and sensors. However, manufacturing approaches for patterning materials over large areas with features on the order of hundreds of microns or less are limited. Acoustophoresis, which uses acoustic forces to control particle arrangement in a fluid medium, is a pathway to address this challenge. This process is dependent on particle and fluid properties and enables patterning of a broad range of materials. Herein, a model with experimental validation is presented to demonstrate that acoustophoresis can be combined with direct-ink writing (DIW) to fabricate line patterns over large cm-scale areas. An in-nozzle particle interaction model was developed to investigate the impact of processing conditions on multi-nodal acoustophoretic DIW. The model predicts patterned line widths within a factor of two relative to experimental results for a high viscosity case study. Here, the model was used to investigate the impact of frequency, particle loading, particle radius, and acoustic pressure on line width and patterning time, providing critical feedback regarding the processing conditions suitable for a target application. Model results illustrate that frequency has the greatest impact on line patterns: increasing from 1 to 3 MHz resulted in a greater than 65% reduction in line width and a greater than 85% reduction in patterning time. Additionally, experiments were conducted with an alumina-epoxy ink and a ~21 cm 2 area pattern was rastered in ~5.5 minutes, demonstrating a path towards large-area line-patterned composite fabrication.
Filamentous fungi produce a wide variety of enzymes in order to efficiently degrade plant cell wall polysaccharides. The production of these enzymes is controlled by transcriptional regulators, which also control the catabolic pathways that convert the released monosaccharides. Two transcriptional regulators, GalX and GalR, control d-galactose utilization in the model filamentous fungus Aspergillus nidulans, while the arabinanolytic regulator AraR regulates l-arabinose catabolism. d-Galactose and l-arabinose are commonly found together in polysaccharides, such as arabinogalactan, xylan and rhamnogalacturonan I. Therefore, the catabolic pathways that convert d-galactose and l-arabinose are often also likely to be active simultaneously. In this study, we investigated the interaction between GalX, GalR and AraR in d-galactose and l-arabinose catabolism. For this, we generated single, double and triple mutants of the three regulators, and analysed their growth and enzyme and gene expression profiles. Our results clearly demonstrated that GalX, GalR and AraR co-regulate d-galactose catabolism in A. nidulans. GalX has a prominent role on the regulation of genes of d-galactose oxido-reductive pathway, while AraR can compensate for the absence of GalR and/or GalX.
Taking a nuclear fuel concept through the research, development, and qualification stages has historically taken on the order of 20 to 25 years because of extensive irradiation tests required for a variety of conditions. The concept of accelerated fuel qualification (AFQ) has been proposed to increase the innovation pace for nuclear fuels. The goal of AFQ is not to replace the traditional qualification approach but rather to reduce the total number of experiments required to ensure approval from the regulatory authority. Of the many AFQ approaches being explored, advanced modeling—and, in particular, mechanistic modeling—is in a uniquely cross-cutting position to reduce the number of required integral tests through the inclusion of separate-effects testing, while helping to extrapolate reactor performance during rare events. We make the case that propagation of uncertainty through various computational length scales helps contextualize mechanistic modeling. We will utilize UO 2 fission gas diffusion predictions from the atomistically informed cluster dynamics code Centipede to inform fuel performance rodlet simulations using the BISON finite element code as the metric for showing how multiscale mechanistic uncertainty quantification can help reduce uncertainty in fuel performance. In conclusion, by quantifying uncertainty and its reduction through multiscale modeling, the qualification process may be accelerated through the reduction of costly irradiation experiments.
Non-negative Matrix Factorization (NMF) has proven to be a powerful unsupervised learning method for uncovering hidden features in complex and noisy datasets with applications in data mining, text recognition, dimension reduction, face recognition, anomaly detection, blind source separation, and many other fields. An important input for NMF is the latent dimensionality of the data, that is, the number of hidden features, K, present in the explored dataset. Unfortunately, and this quantity is rarely known a priori. The existing methods for determining latent dimensionality, such as Automatic Relevance Determination (ARD), are mostly heuristic and utilize different characteristics to estimate the number of hidden features. However, all of them require human presence to make a final determination of K. Here we utilize a supervised machine learning approach in combination with a recent method for model determination, called NMFk, to determine the number of hidden features automatically. NMFk performs a set of NMF simulations on an ensemble of matrices, obtained by bootstrapping the initial dataset, and estimates which K produces stable groups of latent features that reconstruct the initial dataset well. We then train a Multi-Layter Perceptron (MLP) classifier network to determine the correct number of latent features utilizing the statistics and characteristics of the NMF solution, obtained from NMFk. In order to train the MLP classifier, a training set of 58,660 matrices with predetermined latent features were factorized with NMFk. The MLP classifier in conjunction with NMFk maintains a greater than 95% success rate when applied to a held out test set. Additionally, when applied to two well-known benchmark datasets, the swimmer and MIT face data, NMFk/MLP correctly recovers the established number of hidden features. Finally, we compare the accuracy of our method to the ARD, AIC and Stability-based methods.
Multiple anomalous features in electronic spectra of metals with a kagome lattice structure—van Hove singularities, Dirac points, and flat bands—imply that materials containing this structural motif may lie at a nexus of topological and correlated electron physics. Due to the prospects of such exceptional electronic behavior, the recent discovery of superconductivity coexisting with charge-density wave (CDW) order in the layered kagome metals 𝐴V 3 Sb 5 (𝐴 = K,Rb,Cs) has attracted considerable attention. Notably, these archetypal kagome metals express unconventional magnetotransport behavior, including an unexpected linear-in-𝐻 diagonal resistivity at low fields, and an even more peculiar, nonmonotonic sign-changing behavior of the Hall resistivity, which has been speculated to arise from a chiral CDW. We argue here that this unusual magnetotransport derives not from such unconventional phenomena, but rather from the unique fermiology of the 𝐴V 3 Sb 5 materials. Specifically, it is caused by a large, concave hexagonal Fermi surface sheet formed in the close proximity to the van Hove singularities, which is backfolded into a small hexagonal sheet and two large triangular sheets in the CDW state. We introduce and analyze a model of the electronic structure of these Fermi surface sheets that allows for a full analytical treatment within Boltzmann kinetic theory and that enables semi-quantitative fits of our transport data. Specifically, we find that the anomalous magnetotransport behavior is caused by the confluence of strong reduction of the Fermi velocity near the van Hove singularities located near the vertices of the hexagonal sheet and sharp corners in Fermi surface generated by the CDW reconstruction. In conclusion, our analytical approach not only explains the anomalous magnetotransport in the kagome superconductors but also can be extended to a variety of metallic systems hosting singular features in their Fermi surfaces.
Dilute combustion with exhaust gas recirculation (EGR) in spark-ignition engines presents a cost-effective method for achieving higher levels of engine efficiency. At high levels of EGR, however, cycle-to-cycle variability (CCV) of the combustion process is exacerbated by sporadic occurrences of misfires and partial burns. Previous studies have shown that temporal deterministic patterns emerge at such conditions and certain combustion cycles have a significant influence over future events. Due to the complexity of the combustion process and the nature of CCV, harnessing all the deterministic information for control purposes has remained challenging even with physics based 0-D, 1-D, and high-fidelity computational fluid dynamics (CFD) models. In this study, we present a data-driven approach to optimize the combustion process by controlling CCV adjusting the cycle-to-cycle fuel injection quantity. Readily available data from in-cylinder pressure was used to train a spiking neural network (SNN) which learns the optimal way to manage fuel injection in order to reduce CCV while maintaining acceptable levels of fuel consumption. SNNs are particularly well suited for powertrain control applications due to their ability to be deployed on FPGA-based neuromorphic hardware which are small, inexpensive, and have a low power demand. The high-performance computing (HPC) resources of Oak Ridge National Laboratory were used to run an evolutionary-based training approach for choosing the best SNN configuration that minimizes the size of the network while achieving the desired goal. The neuromorphic hardware with the optimized SNN deployed was connected to the rapid prototyping engine control system for real-time control implementation and tested on a single cylinder version of a GM LNF 4-cylinder engine. The results show a significant reduction of CCV with a small percentage of additional fuel used to stabilize the charge.