Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inverse prediction”

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 1,189 records · Page 66

Prediction of low-Z collinear and noncollinear antiferromagnetic compounds having momentum-dependent spin splitting even without spin-orbit coupling

The over 60 years old Rashba-Dresselhaus effect predicts spin-orbit coupling (SOC) induced momentum-dependent spin splitting and spin polarization in materials with noncentrosymmetric structures. Strong SOC induced effects usually require high-atomic number (Z) elements such as rare-earth elements. It has recently been pointed out that antiferromagnets could hold SOC-independent spin splitting and spin polarization. In the present work we develop the spatial and magnetic symmetry conditions enabling such antiferromagnet (AFM)-induced spin splitting, dividing the 1651 magnetic space groups into seven different spin splitting prototypes (SST-1 to SST-7). This analysis places the physics of AFM spin splitting (SST-4) within the broader context of symmetry conditions that enable the more familiar forms of spin splitting, such as ferromagnetic Zeeman effect (SST-5), nonmagnetic no spin splitting (SST-6), and the nonmagnetic Rashba and Dresselhaus effects (SST-7). The AFM-induced spin splitting and spin polarization do not necessarily require breaking of inversion symmetry or the presence of SOC, hence can exist even in centrosymmetric, low-Z light element compounds, considerably broadening the material base for spin manipulations. We use the “inverse design” approach of first formulating the target property (here, spin splitting in low-Z compounds not restricted to low symmetry structures), then derive the enabling physical design principles—the magnetic symmetry conditions—to search realizable compounds that satisfy these a priori design principles. This process uncovers 422 magnetic space groups (160 centrosymmetric and 262 noncentrosymmetric) that could hold AFM-induced, SOC-independent spin splitting and spin polarization. We then search for stable compounds following such enabling symmetries. We investigate the electronic and spin structures of some selected prototype compounds by density functional theory (DFT) and find spin textures that are different than the traditional Rashba-Dresselhaus patterns and exist even in the absence of SOC effect. We provide the DFT results for all antiferromagnetic spin splitting prototypes (SST-1, SST-2, SST-3, SST-4), and concentrate on two limits of SST-4 that are particularly unusual: When spin splitting is momentum dependent (just like the Rashba effect) but is enabled in antiferromagnets even in the absence of SOC in the Hamiltonian. This includes examples of (a) centrosymmetric SST-4A compounds (e.g., orthorhombic LaMnO 3 illustrating collinear AFM, as well as cubic NiS 2 illustrating noncollinear AFM) and (b) noncentrosymmetric SST-4B compounds (e.g., rhombohedral MnTiO 3 illustrating collinear AFM and hexagonal ScMnO 3 illustrating noncollinear AFM). The symmetry design principles outlined here, along with their transformation into an inverse design material search approach and DFT verification, could open the way to their experimental examination.

36 MATERIALS SCIENCE↗

First high-resolution γ-ray spectroscopy of 41 Si

The first high-resolution in-beam γ -ray spectroscopy is reported for the neutron-rich nucleus 41 Si, a tenant of the N = 28 island of inversion. Excited states were populated in the direct one-proton removal reaction from 42 P projectiles and pn removal from 43 P. Seven γ-ray transitions were observed, only one of which had been reported previously in the literature. Furthermore, this makes 41 Si the most neutron-rich odd-even N = 27 isotone with high-resolution excited-state information. For the one-proton removal, the measured partial cross-section distribution to all observed bound final states is contrasted with results from direct one-proton removal calculations that combine eikonal reaction dynamics with SDPF-MU shell-model spectroscopic factors and assume various possible initial states for the poorly known 42 P projectile. Rather distinct calculated cross-section distributions emerge that, in comparison to the new data, imply that the initial state in 42 P is most likely 3 – or 2 – rather than 1 – or 0 – , the predicted shell-model ground state of 42 P. It is further shown that the level scheme from the novel VS-IMSRG calculation closely agrees with the one of SDPF-MU, the most successful phenomenological shell-model effective interaction in describing the much discussed neighboring isotope 42 Si, perhaps cross-validating these complementary approaches on the quest to model rapid shell evolution away from the valley of β stability.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Geometric Scale-up Experiments on Fluidization of Geldart B Glass Beads

The objective of this work is to provide a valuable database from controlled experiments for validating computational models. Recently, coarse-grained techniques such as particle-in-cell (PIC) or coarse-grained discrete element modeling (DEM) have gained popularity due to their computational efficiency while modeling large-scale systems; however, the influence of model parameters and their sensitivities at different geometric scales and flow conditions remain to be analyzed. These datasets are critical for the multiphase flow research community to assess predictive capability of modeling techniques as well as elucidate the hydrodynamic behavior in these systems. This study performed fluidization experiments using three different test sections with internal diameters of 2.5, 4, and 6 in. The operating conditions, bed material, and range of flow velocities at the inlet were constant in all the units, which were not hydrodynamically scaled. Glass beads having a Sauter Mean Diameter of 332 μm were used. Superficial velocity was varied from 2.97 to 5.35 times the minimum fluidization velocity and the initial static bed height was 0.1524 m. The order in which the experiments were performed was randomized and replicates were included to provide uncertainty in measurements. Statistics of differential pressure and bed height from these tests were reported. Future plans include validating PIC methodology in the open-source software, MFiX (Multiphase Flow with Interphase Exchanges) using results from this study. This could further be extended to determine optimal model parameters using inverse techniques such as deterministic calibration or Bayesian inference.

20 FOSSIL-FUELED POWER PLANTS↗

Geometric Scale-up Experiments on Fluidization of Geldart B Glass Beads

The objective of this work is to provide a valuable database from controlled experiments for validating computational models. Recently, coarse-grained techniques such as particle-in-cell (PIC) or coarse-grained discrete element modeling (DEM) have gained popularity due to their computational efficiency while modeling large-scale systems; however, the influence of model parameters and their sensitivities at different geometric scales and flow conditions remain to be analyzed. These datasets are critical for the multiphase flow research community to assess predictive capability of modeling techniques as well as elucidate the hydrodynamic behavior in these systems. This study performed fluidization experiments using three different test sections with internal diameters of 2.5, 4, and 6 in. The operating conditions, bed material, and range of flow velocities at the inlet were constant in all the units, which were not hydrodynamically scaled. Glass beads having a Sauter Mean Diameter of 332 μm were used. Superficial velocity was varied from 2.97 to 5.35 times the minimum fluidization velocity and the initial static bed height was 0.1524 m. The order in which the experiments were performed was randomized and replicates were included to provide uncertainty in measurements. Statistics of differential pressure and bed height from these tests were reported. Future plans include validating PIC methodology in the open-source software, MFiX (Multiphase Flow with Interphase Exchanges) using results from this study. This could further be extended to determine optimal model parameters using inverse techniques such as deterministic calibration or Bayesian inference.

42 ENGINEERING↗

Physics Informed Neural Networks as Computational Physics Emulators

This report is a brief overview and evaluation of Physics Informed Neural Networks (PINNs). Karniadakis and co-workers, e.g., Karniadakis et al. (2021) assert that the PINNs approach integrates seamlessly both data and mathematical physics models, even in partially understood, uncertain and high-dimensional contexts. They further claim that PINNs are effective and efficient for ill-posed and inverse problems, and when combined with domain decomposition, are scalable to large problems and a tool to discover hidden physics. While they demonstrate the capabilities in specific academic instances, their overarching claims about PINNs seem to be an overstatement, at least at the current time. We have briefly considered a few of the limitations of PINNs in this investigation. It is not clear to us if the PINNs approach can ever be competitive with approaches that use specialized algorithms to achieve high-accuracy solutions of governing equations and other techniques that can combine observational data with such solutions. As an example of the latter, consistent with the principles of Bayesian inference, data assimilation, or more generally data-model fusion, is a process that fuses observational data typically with a computational model that respects certain constraints such as conservation laws. For example, improvements in observational network combined with data assimilation have been key in improving weather predictions over the past four decades Kalnay (2003).

97 MATHEMATICS AND COMPUTING↗

Spatial and Temporal Variations in Titan's Surface Temperatures from Cassini CIRS Observations

We report a wide-ranging study of Titan's surface temperatures by analysis of the Moon's outgoing radiance through a spectral window in the thermal infrared at 19 mm (530/cm) characterized by lower atmospheric opacity. We begin by modeling Cassini Composite Infrared Spectrometer (CIRS) far infrared spectra collected in the period 2004-2010, using a radiative transfer forward model combined with a non-linear optimal estimation inversion method. At low-latitudes, we agree with the HASI near-surface temperature of about 94 K at 101S (Fulchignoni et al., 2005). We find a systematic decrease from the equator toward the poles, hemispherically asymmetric, of approx. 1 K at 60 deg. south and approx. 3 K at 60 deg. north, in general agreement with a previous analysis of CIRS data and with Voyager results from the previous northern winter. Subdividing the available database, corresponding to about one Titan season, into 3 consecutive periods, small seasonal changes of up to 2 K at 60 deg N became noticeable in the results. In addition, clear evidence of diurnal variations of the surface temperatures near the equator are observed for the first time: we find a trend of slowly increasing temperature from the morning to the early afternoon and a faster decrease during the night. The diurnal change is approx. 1.5 K, in agreement with model predictions for a surface with a thermal inertia between 300 and 600 J/ sq. m s (exp -1/2) / K. These results provide important constraints on coupled surface-atmosphere models of Titan's meteorology and atmospheric dynamic.

titan↗

Evolution and folding of repeat proteins

Repeat proteins are made with tandem copies of similar amino acid stretches that fold into elongated architectures. These proteins constitute excellent model systems to investigate how evolution relates to structure, folding, and function. Here, we propose a scheme to map evolutionary information at the sequence level to a coarse-grained model for repeat-protein folding and use it to investigate the folding of thousands of repeat proteins. We model the energetics by a combination of an inverse Potts-model scheme with an explicit mechanistic model of duplications and deletions of repeats to calculate the evolutionary parameters of the system at the single-residue level. These parameters are used to inform an Ising-like model that allows for the generation of folding curves, apparent domain emergence, and occupation of intermediate states that are highly compatible with experimental data in specific case studies. We analyzed the folding of thousands of natural Ankyrin repeat proteins and found that a multiplicity of folding mechanisms are possible. Fully cooperative all-or-none transitions are obtained for arrays with enough sequence-similar elements and strong interactions between them, while noncooperative element-by-element intermittent folding arose if the elements are dissimilar and the interactions between them are energetically weak. Additionally, we characterized nucleation-propagation and multidomain folding mechanisms. We show that the global stability and cooperativity of the repeating arrays can be predicted from simple sequence scores.

Ezequiel A. Galpern↗

Ensemble Kalman filter for data assimilation coupled with low-resolution computations techniques applied in fluid dynamics

This paper presents an innovative Reduced-order model (ROM) for merging experimental and simulation data using data assimilation (DA) to estimate the "True" state of a fluid dynamics system, leading to more accurate predictions. Our methodology introduces a novel approach by implementing the ensemble Kalman filter (EnKF) within a reduced-dimensional framework, grounded in a robust theoretical foundation and applied to fluid dynamics. To address the substantial computational demands of DA, the proposed ROM employs low-resolution (LR) techniques to drastically reduce computational costs. This innovative approach involves downsampling datasets for DA computations, followed by an advanced reconstruction technique based on low-cost singular value decomposition (lcSVD). The lcSVD method, a key innovation in this paper, has never been applied to DA before and offers a highly efficient way to enhance resolution with minimal computational resources. Our results demonstrate significant reductions in both computation time and RAM usage through these LR techniques without compromising the accuracy of the estimations. For instance, in a turbulent test case, for a data compression rate of 15.9, the LR approach can achieve a speed-up of 13.7 and a RAM compression of 90.9% while maintaining a low relative root mean square error (RRMSE) of 2.6%, compared to 0.8% in the high-resolution (HR) reference. Furthermore, we highlight the effectiveness of the EnKF in estimating and predicting the state of fluid flow systems based on limited observations and given low-fidelity numerical data. This paper highlights the potential of the proposed DA method in fluid dynamics applications, particularly for improving computational efficiency in CFD and related fields. Its ability to balance accuracy with low computational and memory costs makes it especially suitable for large-scale and real-time applications, such as environmental monitoring or engineering design. This method will be incorporated into ModelFLOWs-app.

Data Assimilation↗

Investigating the Variation in Solvation Interactions of Choline Chloride-Based Deep Eutectic Solvents Formed Using Different Hydrogen Bond Donors

Choline chloride-based deep eutectic solvents (DESs) have gained widespread recognition as green solvents in catalysis and separation science as their designer properties permit solvation properties to be modulated by the choice of the hydrogen bond donor (HBD) as well as its molar ratio. Despite being one of the most popular classes of DESs, very little is currently understood regarding how their solvation characteristics differ among themselves as well as from other classes of DESs. Previous studies have shown that the catalytic activity, extraction efficiency, and solubilities of natural compounds can be influenced by the type of choline salt and HBD, but no study has systematically related these parameters with its solvation properties. In this study, inverse gas chromatography was employed for the first time to study the individual solvation characteristics of a diverse range of choline chloride and acetate-based DESs composed of amides, diols, organic acids, and carbohydrates. Furthermore, hydrogen bond acidity was found to be a dominant interaction for all DESs, especially those composed of organic acid HBDs. DESs composed of isomers of butane diol and hexane diol as HBDs exhibited higher hydrogen bond basicity and dispersive-type interactions compared to those composed of urea, acetamide, and organic acids with the position of the hydroxyl functional group influencing their solvation interactions. Choline acetate-based DESs exhibited lower dipolarity and stronger hydrogen bond basicity and dispersive-type interactions compared to choline chloride DESs. The solvation models developed in this study were used to interpret and explain DES behavior in previously reported studies involving catalysis, organic synthesis, and extractions demonstrating that the measured solute–solvent interactions can serve to predict the performance of choline chloride-based DESs when used in various applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Observation of paramagnetic spin-degeneracy lifting in EuZn 2 Sb 2

Taken together, time-reversal and spatial inversion symmetries impose a twofold spin degeneracy of the electronic states in crystals. In centrosymmetric materials, this degeneracy can be lifted by introducing magnetism, either via an externally applied field or through internal magnetization. However, a correlated alignment of spins, even in the paramagnetic phase, can lift the spin degeneracy of electronic states. Here, we report an in-depth study of the electronic band structure of the Eu-ternary pnictide EuZn 2 Sb 2 through a combination of high-resolution angle-resolved photoemission spectroscopy measurements and first-principles calculations. An analysis of the photoemission line shapes over a range of incident photon energies and sample temperatures is shown to reveal the presence of band spin-degeneracy lifting in the paramagnetic phase. Our angle-resolved photoemission spectroscopy results are in good agreement with theoretical ferromagnetic-phase calculations, which indicates the importance of ferromagnetic fluctuations in the system. Through our calculations, we predict that spin-polarized bands in EuZn 2 Sb 2 generate a single pair of Weyl nodes. Our observation of band splittingin EuZn 2 Sb 2 provides a key step toward realizing time-reversal symmetry breaking physics in the absence of long-range magnetic order.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Evaluating Northern Hemisphere Growing Season Net Carbon Flux in Climate Models Using Aircraft Observations

Understanding terrestrial ecosystems and their response to anthropogenic climate change requires quantification of land-atmosphere carbon exchange. However, top-down and bottom-up estimates of large-scale land-atmosphere fluxes, including the northern extratropical growing season net flux (GSNF), show significant discrepancies. We developed a data-driven metric for the GSNF using atmospheric carbon dioxide concentration observations collected during the High-Performance Instrumented Airborne Platform for Environmental Research Pole-to-Pole Observations and Atmospheric Tomography Mission flight campaigns. This aircraft-derived metric is bias-corrected using three independent atmospheric inversion systems. We estimate the northern extratropical GSNF to be 5.7 ± 0.3 Pg C and use it to evaluate net biosphere productivity from the Coupled Model Intercomparison Project phase 5 and 6 (CMIP5 and CMIP6) models. While the model-to-model spread in the GSNF has decreased in the CMIP6 models relative to that of the CMIP5 models, there is still disagreement on the magnitude and timing of seasonal carbon uptake with most models underestimating the GSNF and overestimating the length of the growing season relative to the observations. We also use an emergent constraint approach to estimate annual northern extratropical gross primary productivity to be 56 ± 17 Pg C, heterotrophic respiration to be 25 ± 13 Pg C, and net primary productivity to be 28 ± 12 Pg C. The flux inferred from these aircraft observations provides an additional constraint on large-scale gross fluxes in prognostic Earth system models that may ultimately improve our ability to accurately predict carbon-climate feedbacks.

54 ENVIRONMENTAL SCIENCES↗

Assessing Relay Communications for Mars Sample Return Surface Mission Concepts

The Mars Sample Return (MSR) Campaign would be a 3-mission campaign concept supported by NASA and ESA to return samples from the Mars surface. MSR would, for the first time ever, present a need to communicate with multiple surface assets that are co-located on Mars in a coordinated effort to accomplish the unified objective of fetching, transporting, and returning samples from Mars. Currently, Mars surface assets relay data to and from Earth using a number of orbiters in what’s known as the Mars Relay Network (MRN). This network is characterized by a small number of surface assets distributed across the Martian globe and a larger number of orbiters to provide relay services. As of June 2020, there are two surface assets for which five orbiters are providing relay. During the MSR Campaign, there would be two rovers and a lander that all would require relay communication from a small number of Mars orbiters to meet the aggressive MSR timeline. The inversion of the current MRN paradigm, a system of many surface assets requiring relay and few orbiters to provide relay, necessitates the unique challenge of optimally allocating relay passes to maximize the operational capability of all assets. The allocation must consider a large number of trade variables including Mars asset operational requirements and Earth ground system constraints, including staffing schedules, operations planning across time zones, and more. To address these telecommunication challenges, the Mars Asset Relay Mission Link Allocation Design Environment (MARMLADE) tool was developed. It is a MATLAB-based tool to assign orbiter passes or Direct-From-Earth (DFE) links to each of the three surface assets and quantify the operational efficiency of each surface asset.MARMLADE uses a data set of simulated Mars relay orbiter geometry and telecommunication capabilities provided by JPL’s Telecom Orbit Analysis and Simulation Tool (TOAST) software to compute which asset should get each pass based on a series of heuristics and predictions of all assets’ states. Within MARMLADE, the user can provide inputs including the option for time-based pass splitting, fixed FWD data rate capabilities, DFE communication capabilities, and link parameters allowing for the assessment of complex operations and hardware trades using surface mission operational efficiency as a primary figure of merit. As the MSR mission concepts continue to mature, MARMLADE is being used to assess ability of all MSR elements to meet the surface mission timeline requirements and to provide relay link allocations to each of the MSR surface assets.This paper will describe the motivation and design of the MARMLADE tool and how it is being used to perform campaign and mission level trades, generate requirements, and support development of the MSR surface mission scenarios.

Lee, Charles↗

An Assessment of Molten Metal Detachment Hazards During Electron Beam Welding in Space

The safety issue has been raised with regards to potential molten metal detachments from the weld pool and cold filler wire during electron beam welding in space. This investigation was undertaken to evaluate if molten metal could detach and come in contact with astronauts and burn through the fabric of the astronauts' Extravehicular Mobility Unit (EMU) during electron beam welding in space. Molten metal detachments from either the weld/cut substrate or weld wire could present harm to a astronaut if the detachment was to burn through the fabric of the EMU. Theoretical models were developed to predict the possibility and size of the molten metal detachment hazards during the electron beam welding exercises at Low Earth Orbit (LEO). The primary molten metal detachment concerns were those cases of molten metal separation from the metal surface due to metal cutting, weld pool splashing, entrainment and release of molten metal due to filler wire snap-out from the weld puddle, and molten metal accumulation and release from the end of the weld wire. Some possible ways of obtaining molten metal drop detachments would include an impulse force, or bump, to the weld sample, cut surface, or filler wire. Theoretical models were developed for these detachment concerns from principles of impact and kinetic energies, surface tension, drop geometry, surface energies, and particle dynamics. The surface tension represents the force opposing the liquid metal drop from detaching whereas the weight of the liquid metal droplet represents a force that is tending to detach the molten metal drop. Theoretical calculations have indicated that only a small amount of energy is required to detach a liquid metal drop; however, much of the energy of an impact is absorbed in the sample or weld plate before it reaches the metal drop on the cut edge or surface. The tendency for detachment is directly proportional to the weld pool radius and metal density and inversely proportional to the surface tension of the liquid metal. For a detachment the initial kinetic energy of the weld pool with respect to the plate has to exceed the energy to form the extra surface required for the detachment of the pool. The difficulty is in transferring the energy from the point of impact through the plate and sample to the cut edge. It is likely that not all of the kinetic energy is available for detaching the pool; some may be sequestered in weld pool oscillations. The coefficient of restitution for the collision will be lower than one if irreversible deformation, for example plastic flow deformation, takes place during the collision. Thus determining the amount of energy from an impact that actually reaches the molten metal droplet is critical. Various molten metal detachment scenarios were tested experimentally in an enclosed vacuum chamber using the Ukrainian Universal Hand Tool, an electron beam welder designed for space welding. The experimental testing was performed in a 4 ft. X 4 ft. vacuum chamber at Marshall Space Flight Center, evacuated to vacuum levels of at least 50 microTorr, and also some welding garment material was utilized to observe the effect of the molten metal detachments on the material. A "carillon" apparatus consisting of four pendulum hammer strikers, each weighing approximately 3.65 lbs, raised to predetermined specific heights was used to apply an impact force to the weld sample/plate during electron beam welding and cutting exercises. The strikers were released by switching on an electric motor to rotate a pin holding wires retaining the strikers at desired heights. The specimens were suspended so as to be free to respond to the blows with a sudden velocity increment. The specimens were mounted on a hinged plate for minimizing effective mass with the option to fasten it down so as to raise its effective mass closer to that anticipated for an actual space welding scenario. Measurements were made of the impact energy and the horizontal fling distances of the detached metal drops. It was not particularly easy to generate the detachments for this experiment. This document presents the details of the theoretical modeling effort and a summary of the experimental effort to measure molten metal drop detachments from terrestrial electron beam welding in the enclosed vacuum chamber. The results of the experimental effort have shown that molten metal detachments can occur from the sample/weld plate only if a sufficiently large impact force is applied to the weld plate. A "weld pool detachment parameter" was determined to indicate whether detachment would occur. Detachment can be either full or partial (dripping), Partial detachment means that the weld pool detached from one side of the liquid-solid boundary so as to leave a hole at the puddle site but remained attached over part of the liquid-solid boundary and dripped down the plate with no fully detached material detected. Full detachment, however, does not necessarily mean that the whole pool fully detached; in some cases only a smaller portion of the pool detached, the remainder dripping down the plate. The weld pool detachment parameter according to theory and according to the empirical data allows a determination of whether full detachments might occur. Theoretical calculations indicated titanium alloy would be the most difficult from which to detach molten metal droplets followed by stainless steel and then by aluminum. The experimental results were for the most part consistent with the theoretical analysis and predictions. The above theory is applicable to other situations as desired for assessing the potential for molten metal detachments.

Fragomeni, James M.↗

Quadrature Based Neural Network Learning of Stochastic Hamiltonian Systems

Hamiltonian Neural Networks (HNNs) provide structure-preserving learning of Hamiltonian systems. In this paper, we extend HNNs to structure-preserving inversion of stochastic Hamiltonian systems (SHSs) from observational data. We propose the quadrature-based models according to the integral form of the SHSs’ solutions, where we denoise the loss-by-moment calculations of the solutions. The integral pattern of the models transforms the source of the essential learning error from the discrepancy between the modified Hamiltonian and the true Hamiltonian in the classical HNN models into that between the integrals and their quadrature approximations. This transforms the challenging task of deriving the relation between the modified and the true Hamiltonians from the (stochastic) Hamilton–Jacobi PDEs, into the one that only requires invoking results from the numerical quadrature theory. Meanwhile, denoising via moments calculations gives a simpler data fitting method than, e.g., via probability density fitting, which may imply better generalization ability in certain circumstances. Numerical experiments validate the proposed learning strategy on several concrete Hamiltonian systems. The experimental results show that both the learned Hamiltonian function and the predicted solution of our quadrature-based model are more accurate than that of the corrected symplectic HNN method on a harmonic oscillator, and the three-point Gaussian quadrature-based model produces higher accuracy in long-time prediction than the Kramers–Moyal method and the numerics-informed likelihood method on the stochastic Kubo oscillator as well as other two stochastic systems with non-polynomial Hamiltonian functions. Moreover, the Hamiltonian learning error εH arising from the Gaussian quadrature-based model is lower than that from Simpson’s quadrature-based model. These demonstrate the superiority of our approach in learning accuracy and long-time prediction ability compared to certain existing methods and exhibit its potential to improve learning accuracy via applying precise quadrature formulae.

Mathematics↗

Relation between Hertz Stress-Life Exponent, Ball-Race Conformity, and Ball Bearing Life

ANSI/ABMA and ISO standards based on Lundberg-Palmgren bearing life theory are normalized for ball bearings having inner- and outer-race conformities of 52 percent (0.52) and made from pre-1940 bearing steel. The Lundberg-Palmgren theory incorporates an inverse 9th power relation between Hertz stress and fatigue life for ball bearings. The effect of race conformity on ball set life independent of race life is not incorporated into the Lundberg-Palmgren theory. In addition, post-1960 vacuum-processed bearing steel exhibits a 12th power relation between Hertz stress and life. The work reported extends the previous work of Zaretsky, Poplawski, and Root to calculate changes in bearing life, that includes the life of the ball set, caused by race conformity, Hertz stress-life exponent, ball bearing type and bearing series. The bearing fatigue life in actual application will usually be equal to or greater than that calculated using the ANSI/ABMA and ISO standards that incorporate the Lundberg-Palmgren theory. The relative fatigue life of an individual race is more sensitive to changes in race conformity for Hertz stress-life exponent n of 12 than where n = 9. However, when the effects are combined to predict actual bearing life for a specified set of conditions and bearing geometry, the predicted life of the bearing will be greater for a value of n = 12 than n = 9.

Zaretsky, Erwin V.↗

Relation Between Hertz Stress-Life Exponent, Ball-Race Conformity, and Ball Bearing Life

ANSI/ABMA and ISO standards based on Lundberg-Palmgren bearing life theory are normalized for ball bearings having inner- and outerrace conformities of 52 percent (0.52) and made from pre-1940 bearing steel. The Lundberg-Palmgren theory incorporates an inverse 9th power relation between Hertz stress and fatigue life for ball bearings. The effect of race conformity on ball set life independent of race life is not incorporated into the Lundberg-Palmgren theory. In addition, post-1960 vacuum-processed bearing steel exhibits a 12th power relation between Hertz stress and life. The work reported extends the previous work of Zaretsky, Poplawski, and Root to calculate changes in bearing life--that includes the life of the ball set--caused by race conformity, Hertz stress-life exponent, ball bearing type and bearing series. The bearing fatigue life in actual application will usually be equal to or greater than that calculated using the ANSI/ABMA and ISO standards that incorporate the Lundberg-Palmgren theory. The relative fatigue life of an individual race is more sensitive to changes in race conformity for Hertz stress-life exponent n of 12 than where n = 9. However, when the effects are combined to predict actual bearing life for a specified set of conditions and bearing geometry, the predicted life of the bearing will be greater for a value of n = 12 than n = 9.

Zaretsky, Erwin V.↗

Spectroscopic Probe Molecule Selection Using Quantum Theory, First-Principles Calculations, and Machine Learning

Probe molecule vibrational spectra have a long history of being used to characterize materials including metals, oxides, metal-organic frameworks, and even human proteins. Furthermore, recent advances in machine learning have enabled computationally generated spectra to aid in detailed characterization of complex surfaces with probe molecules. Despite widespread use of probe molecules, the science of probe molecule selection is underdeveloped. Here, we develop physical concepts, including orbital interaction energy and the energy overlap integral, to explain and predict the ability of probe molecules to discriminate structural descriptors. We resolve the crystal orbital overlap population (COOP) to specific molecular orbitals and quantify their bonding character, which directly influences vibrational frequencies. Using only a single adsorbate calculation from density function theory (DFT), we compute the interaction energy of individual adsorbate molecular orbitals with adsorption site atomic orbitals across many different sites. Combining the molecular orbital resolved COOP and changes in orbital interaction energy enables probe molecule selection for improved discrimination of various sites. We demonstrate these concepts by comparing the predicted effectiveness of carbon monoxide (CO), nitric oxide (NO), and ethylene (C 2 H 4 ) to probe Pt adsorption sites. Finally, using a previously developed machine learning framework, we show that models trained on hundreds of thousand C 2 H 4 spectra, computed from DFT, which regress surface binding-type and generalized coordination number (GCN), outperform those trained using CO and NO spectra. Lastly, a python package, pDOS_overlap, for implementing the electron density based analysis on any combination of adsorbates and materials, is also made available.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

The Langdon effect in laser plasmas: Absorption and conduction

A plasma heated by inverse bremsstrahlung absorption of laser light develops a non-Maxwellian electron distribution function, called the Langdon effect [A. B. Langdon, Phys. Rev. Lett. 44, 575 (1980)]. These non-Maxwellian distributions are sufficiently long-lived to impact the absorption processes itself as well as the transport of heat by electrons. The theory of the Langdon effect in a homogeneous plasma is reviewed to clarify some aspects of Langdon's derivation as well as to confirm that the widely used super-Gaussian approximation works fairly well to describe the shape of the distribution function and reduction of the absorption rate. The Langdon effect on thermal conduction in an inhomogeneous plasma is developed by considering perturbations in a homogeneous absorbing plasma, which develops a heat flux due to both temperature and density gradients. A practical theory of the heat flux is developed by fitting the results of Vlasov–Fokker–Planck simulations, which avoids several approximations that compromised the usefulness of past theoretical predictions, most critically, the effect of electron–electron collisions on the fluxes. The present fits parameterize the coefficients of the temperature gradient (thermal conductivity) and the density gradient for a plasma of any ionization state and for any laser intensity where the theory of the Langdon effect remains locally valid. It is expected that this generalized theory of heat flow in an absorbing plasma will improve the predictive capability of radiation-hydrodynamics simulations of laser-produced plasmas, especially those formed in inertial confinement fusion experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗