Discrete element modeling of a particle heater for energy storage systems
Not Available
Schirck, Jason↗
SEARCH · Search NASA
Results for “low order model”
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.
In situ recording of Mars soundscape
Abstract Before the Perseverance rover landing, the acoustic environment of Mars was unknown. Models predicted that: (1) atmospheric turbulence changes at centimetre scales or smaller at the point where molecular viscosity converts kinetic energy into heat 1 , (2) the speed of sound varies at the surface with frequency 2,3 and (3) high-frequency waves are strongly attenuated with distance in CO 2 (refs. 2–4 ). However, theoretical models were uncertain because of a lack of experimental data at low pressure and the difficulty to characterize turbulence or attenuation in a closed environment. Here, using Perseverance microphone recordings, we present the first characterization of the acoustic environment on Mars and pressure fluctuations in the audible range and beyond, from 20 Hz to 50 kHz. We find that atmospheric sounds extend measurements of pressure variations down to 1,000 times smaller scales than ever observed before, showing a dissipative regime extending over five orders of magnitude in energy. Using point sources of sound (Ingenuity rotorcraft, laser-induced sparks), we highlight two distinct values for the speed of sound that are about 10 m s −1 apart below and above 240 Hz, a unique characteristic of low-pressure CO 2 -dominated atmosphere. We also provide the acoustic attenuation with distance above 2 kHz, allowing us to explain the large contribution of the CO 2 vibrational relaxation in the audible range. These results establish a ground truth for the modelling of acoustic processes, which is critical for studies in atmospheres such as those of Mars and Venus.
79 ASTRONOMY AND ASTROPHYSICS↗
Planar, longitudinal, compressive waves in solids: Thermodynamics and uniaxial strain restrictions
A new tri-axial pressure-based constitutive expression has been found using Cauchy's stress tensor. This stress state emphasizes pressure and shear stress. The description is a pressure plus an effective shear stress allowing for a constitutive law based on atomic solid-state phase changes in crystalline cells due to pressure plus shear-based dislocation motion commonly associated with plasticity. Pressure has a new role in the material's constitutive response as it is separated from plasticity. The thermo-mechanical system describes third-order Gibbs’ expressions without specific volume restrictions placed upon the material. Isothermally, the ratio of heat to shear work in elastic copper is shown to approach zero at a very low temperature and become larger than one as temperature approaches melting. Wave compression models investigated are elastic and plastic: in fully elastic materials, the planar wave is restricted by Poisson's effect although plastic shear changes this constraint. Plastic deformation, dominated by dissipative shear stresses in uniaxial strain, heats the material while excluding phase changes from hydrostatic pressures. The material properties per se across Hugoniot shocks are described with entropy concepts. Shock waves are exceedingly complex since the constitutive laws are linked at extreme temperatures, pressures, and shear stresses. Isothermal, isentropic, isochoric, and iso-shear conditions are used throughout with Jacobian algebra.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Does Bayesian model averaging improve polynomial extrapolations? Two toy problems as tests
We assess the accuracy of Bayesian polynomial extrapolations from small parameter values, x, to large values of x. We consider a set of polynomials of fixed order, intended as a proxy for a fixed-order effective field theory (EFT) description of data. We employ Bayesian model averaging (BMA) to combine results from different order polynomials (EFT orders). Our study considers two 'toy problems' where the underlying function used to generate data sets is known. We use Bayesian parameter estimation to extract the polynomial coefficients that describe these data at low x. A 'naturalness' prior is imposed on the coefficients, so that they are $\mathcal{O}(1)$. We BMA different polynomial degrees by weighting each according to its Bayesian evidence and compare the predictive performance of this BMA with that of the individual polynomials. In conclusion, the credibility intervals on the BMA forecast have the stated coverage properties more consistently than does the highest evidence polynomial, though BMA does not necessarily outperform every polynomial.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Convergent close-coupling calculations of electron scattering on LiH
Cross sections for electron scattering on the ground state of LiH are calculated using the molecular convergent close-coupling method. The fixed-nuclei approximation is utilized, and calculations are performed at the mean internuclear separation of 3.06 𝑎 0 for projectile energies ranging from 0.1 to 500 eV. Here, a model-potential approach is utilized to treat the LiH molecule as a two-electron system, and good agreement is found with the literature for various properties of the structure model. Cross sections are presented for elastic scattering, excitation to several electronic states, and ionization. Comparison is made with existing results, and qualitative agreement is found with small close-coupling calculations at low projectile energies. Good agreement is found with first-order results at high projectile energies.
74 ATOMIC AND MOLECULAR PHYSICS↗
First-principles theory of the Dirac semimetal Cd 3 As 2 under Zeeman magnetic field
Time-reversal broken Weyl semimetals have attracted much attention recently, but certain aspects of their behavior, including the evolution of their Fermi surface topology and anomalous Hall conductivity with Fermi-level position, have remained underexplored. A promising route to obtain such materials may be to start with a nonmagnetic Dirac semimetal and break time-reversal symmetry via magnetic doping or magnetic proximity. Here we explore this scenario in the case of the Dirac semimetal Cd 3 As 2 , based on first-principles density-functional calculations and subsequent low-energy modeling of Cd 3 As 2 in the presence of a Zeeman field applied along the symmetry axis. We clarify how each four-fold degenerate Dirac node splits into four Weyl nodes, two with chirality ±1 and two higher-order nodes with chirality ±2. Using a minimal k · p model Hamiltonian whose parameters are fit to the first-principles calculations, we detail the evolution of the Fermi surfaces and their Chern numbers as the Fermi energy is scanned across the region of the Weyl nodes at fixed Zeeman field. We also compute the intrinsic anomalous Hall conductivity as a function of Fermi-level position, finding a characteristic inverted-dome structure. Cd 3 As 2 is especially well suited to such a study because of its high mobility, but the qualitative behavior revealed here should be applicable to other Dirac semimetals as well.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Correlated metals and unconventional superconductivity in rhombohedral trilayer graphene: A renormalization group analysis
Here, motivated by recent experimental observations of correlated metallic phases and superconductivity in rhombohedral trilayer graphene (RTG), we perform an unbiased study of electronic ordering instabilities in hole-doped RTG. Specifically, we focus on electronic states energetically proximate to Van Hove singularities (VHSs), where a large density of states promotes different interaction-induced symmetry-breaking electronic orders. To resolve the Fermi surface near VHSs, we construct a fermionic hot-spot model and demonstrate that a perpendicular electric field can tune different nesting structures of the Fermi surface. Subsequently, we apply a renormalization group analysis to describe the low-energy phase diagrams of our model under both short-range repulsive interactions as well as realistic (long-range) Coulomb interactions. Our analysis shows instabilities towards either intervalley coherent metallic phases or superconducting phases. The dominant pairing channel depends crucially on the nature of Fermi surface nesting—repulsive Coulomb interaction favors spin-singlet d-wave pairing for relatively small displacement field and spin-singlet i-wave pairing for larger displacement field. We argue that the phase diagram of RTG can be well-understood by modeling the realistic Coulomb interaction as the sum of repulsive density-density interaction and ferromagnetic spin-triplet intervalley coherence (IVC) Hund's coupling, while phonon-mediated electronic interactions have a negligible effect on this system, in sharp contrast to twisted graphene multilayers.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
X 17 boson and the H 3 ( p , e + e – ) He 4 and He 3 ( n , e + e – ) He 4 processes: A theoretical analysis
The present work deals with e + – e – pair production in the four-nucleon system. We first analyze the process as a purely electromagnetic one in the context of a state-of-the-art approach to nuclear strong-interaction dynamics and nuclear electromagnetic currents, derived from chiral effective field theory ( χ EFT ) . Next, we examine how the exchange of a hypothetical low-mass boson would impact the cross section for such a process. We consider several possibilities, that this boson is either a scalar, pseudoscalar, vector, or axial particle. The ab initio calculations use exact hyperspherical-harmonics methods to describe the bound state and low-energy spectrum of the A = 4 continuum, and they fully account for initial state interaction effects in the 3 + 1 clusters. While electromagnetic interactions are treated to high orders in the chiral expansion, the interactions of the hypothetical boson with nucleons are modeled in leading-order χ EFT (albeit, in some instances, selected subleading contributions are also accounted for). We also provide an overview of possible future experiments probing pair production in the A = 4 system at a number of candidate facilities.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Simulation and Development of the Radial Time Projection Chamber For the Bonus12 Experiment in CLAS12
Knowledge of the structure of nucleons (i.e. protons and neutrons) is a central topic of interest to nuclear/particle physicists. Much more is known about the structure of the proton than the neutron due to the lack of high-density free neutron targets. The Barely Off-shell Nucleon Structure experiment (BONuS12) at Jefferson Lab (JLab) is a second generation experiment upgraded/optimized to advance our knowledge of the neutron's structure using the deep-inelastic scattering of electrons off deuterium. Typically, since deuterium is a nuclear target, corrections for off-shell and nuclear binding effects must be taken into account in order to extract results on the neutron. These corrections are model-dependent and therefore have limited our success in extracting neutron information using deuterium targets. In the BONuS12 experiment, 10.6 GeV electrons are scattered off of a deuterium target. By detecting the low momentum spectator proton at backward angles, the uncertainty due to final state interactions is minimized. The goal of the experiment is to measure the ratio of the neutron to proton structure functions ($F^n_2/F^p_2$) as the Bjroken scaling variable x approaches 1. The newly designed Radial Time Projection Chamber (RTPC) for BONuS12 detects the spectator proton in coincidence with the scattered electron, which is detected in the CEBAF Large Acceptance Spectrometer (CLAS12). This work presents the simulation and development of the new BONuS12 RTPC. The design, construction, and testing of the Drift-gas Monitoring Sysytem (DMS) for the BONuS12 experiment is also described. The results of the DMS operation as well as the first preliminary data from the BONuS12 experimental run are given. Because the BONuS12 data analysis depends on CLAS12 working effectively, an effort to verify the CLAS12 operation with the extraction of the inclusive deep inelastic cross section from the first experiment in CLAS12 (Run Group A) will be presented.
Dzbenski, Nathan↗
Gravitational waves from first-order phase transitions in Majoron models of neutrino mass
We show how the generation of right-handed neutrino masses in Majoron models may be associated with a first-order phase transition and accompanied by the production of a stochastic background of gravitational waves (GWs). We explore different energy scales with only renormalizable operators in the effective potential. If the phase transition occurs above the electroweak scale, the signal can be tested by future interferometers. We consider two possible energy scales for phase transitions below the electroweak scale. If the phase transition occurs at a GeV, the signal can be tested at LISA and provide a complementary cosmological probe to right-handed neutrino searches at the FASER detector. If the phase transition occurs below 100 keV, we find that the peak of the GW spectrum is two or more orders of magnitude below the putative NANOGrav GW signal at low frequencies, but well within reach of the SKA and THEIA experiments. We show how searches of very low frequency GWs are motivated by solutions to the Hubble tension in which ordinary neutrinos interact with the dark sector. We also present general calculations of the phase transition temperature and Euclidean action that apply beyond Majoron models.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Kondo Lattice Model of Magic-Angle Twisted-Bilayer Graphene: Hund’s Rule, Local-Moment Fluctuations, and Low-Energy Effective Theory
For this work, we apply a generalized Schrieffer-Wolff transformation to the extended Anderson-like topological heavy fermion (THF) model for the magic-angle (θ=1.05°) twisted bilayer graphene (MATBLG) [Phys. Rev. Lett. 129, 047601 (2022)], to obtain its Kondo lattice limit. In this limit localized $\mathcal{f}$ electrons on a triangular lattice interact with topological conduction $\mathcal{c}$ electrons. By solving the exact limit of the THF model, we show that the integer fillings ν=0, ±1, ±2 are controlled by the heavy $\mathcal{f}$ electrons, while ν=±3 is at the border of a phase transition between two $\mathcal{f}$-electron fillings. For ν=0, ±1, ±2, we then calculate the Ruderman-Kittel-Kasuya-Yosida (RKKY) interactions between the $\mathcal{f}$ moments in the full model and analytically prove the SU(4) Hund’s rule for the ground state which maintains that two $\mathcal{f}$ electrons fill the same valley-spin flavor. Our (ferromagnetic interactions in the) spin model dramatically differ from the usual Heisenberg antiferromagnetic interactions expected at strong coupling. We show the ground state in some limits can be found exactly by employing a positive semidefinite “bond-operators” method. We then compute the excitation spectrum of the $\mathcal{f}$ moments in the ordered ground state, prove the stability of the ground state favored by RKKY interactions, and discuss the properties of the Goldstone modes, the (reason for the accidental) degeneracy of (some of) the excitation modes, and the physics of their phase stiffness. We develop a low-energy effective theory for the $\mathcal{f}$ moments and obtain analytic expressions for the dispersion of the collective modes. We discuss the relevance of our results to the spin-entropy experiments in TBG.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Synthesis and characterization of Pt(Cu 0.67 Sn 0.33 )
Pt(Cu 0.67 Sn 0.33 ) has recently been found in a natural sample. In order to be able to characterize this new ternary compound, we synthesized it from the elements. Samples were characterized by X-ray powder diffraction, differential scanning calorimetry, thermal relaxation calorimetry, and scanning electron microscopy studies. Density functional theory-based model calculations complemented the experimental studies. Pt(Cu 0.67 Sn 0.33 )was already formed at a relatively low temperature of 773 K. Rietveld refinement of Pt(Cu 0.67 Sn 0.33 ) has been carried out in CuAu-type or L1 0 -type structure, space group P4/mmm, with Pt on 0,0,0 and disordered Cu and Sn on 1/2, 1/2, 1/2 and Z = 1. The lattice parameters are a = 2.823(1) Å, c = 3.64(1) Å, and V = 29.00(4) Å which are in good agreement with values obtained earlier on the natural sample and with the results of DFT calculations. The vibrational entropy for Pt(Cu 0.67 Sn 0.33 ) is $S_{298.15}^{vib}$ = 79.9(7) J mol -1 K -1 . The pressure dependence up to 36(2) GPa of the unit-cell volume and the lattice parameters and unit-cell volume have been obtained by synchrotron based powder diffraction using a diamond anvil cell. A fit of a 3rd-order Birch–Murnaghan equation of state to the Pt(Cu 0.67 Sn 0.33 )) (p,V)-data results in a bulk modulus of B 0 = 215(27) GPa and B' = 5(2).
36 MATERIALS SCIENCE↗
Aligned yet large dipoles: a SMEFT study
We study a non-universal flavor scenario at the level of the Standard Model Effective Field Theory, according to which the matrix of Wilson coefficients cuW of an up-type electroweak quark dipole operator is aligned with the up-type Yukawa coupling. Such an alignment usually follows from the assumption of Minimal Flavor Violation (MFV), away from which we step by allowing the entries of cuW to be sizable along the first quark generations. A particular example, which we refer to as “inverse hierarchy MFV”, features Wilson coefficients inversely proportional to quark masses, and arises from BSM models respecting MFV and containing heavy fields that replicate the mass hierarchy of SM quarks. We then analyze the phenomenology driven by cuW at colliders and at lower-energy flavor experiments. We show that precision measurements of the process pp → Wh → γγℓν at FCC-hh could set an upper bound on |cuW| ≲ O (10 -2 )(Λ/TeV) 2 , with Λ the cutoff of the effective field theory. This bound is an order of magnitude stronger than the existing LHC bounds. Moreover, we estimate that Wh → b$\overline{b}$ℓν at HL-LHC could also give competitive bounds. In the low-energy regime, we consider bounds arising from rare kaon decays, which turn out to be loose, |c$_{uW}^{11}$ | < O (1)(Λ/TeV) 2 . We finally demonstrate that our flavor and operator assumptions can be derived from a weakly-coupled UV model, which we choose to simultaneously illustrate the UV origin of inverse hierarchy MFV.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Influence of functional groups on low-temperature combustion chemistry of biofuels
Ongoing progress in synthetic biology, metabolic engineering, and catalysis continues to produce a diverse array of advanced biofuels with complex molecular structure and functional groups. In order to integrate biofuels into existing combustion systems, and to optimize the design of next-generation combustion systems, understanding connections between molecular structure and ignition at low-temperature conditions (< 1000 K) remains a priority that is addressed in part using chemical kinetics modeling. The development of predictive models relies on detailed information, derived from experimental and theoretical studies, on molecular structure and chemical reactivity, both of which influence the balance of chain reactions that occur during combustion – propagation, termination, and branching. In broad context, three main categories of reactions affect ignition behavior: (i) initiation reactions that generate a distribution of organic radicals, ; (ii) competing unimolecular decomposition of and bimolecular reaction of with O2; (iii) decomposition mechanisms of peroxy radical adducts (RO), including isomerization via RO ⇌ OOH. Furthermore, all three categories are influenced by functional groups in different ways, which causes a shift in the balance of chain reactions that unfold over complex temperature- and pressure-dependent mechanisms.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
ECNet is an evolutionary context-integrated deep learning framework for protein engineering
Abstract Machine learning has been increasingly used for protein engineering. However, because the general sequence contexts they capture are not specific to the protein being engineered, the accuracy of existing machine learning algorithms is rather limited. Here, we report ECNet (evolutionary context-integrated neural network), a deep-learning algorithm that exploits evolutionary contexts to predict functional fitness for protein engineering. This algorithm integrates local evolutionary context from homologous sequences that explicitly model residue-residue epistasis for the protein of interest with the global evolutionary context that encodes rich semantic and structural features from the enormous protein sequence universe. As such, it enables accurate mapping from sequence to function and provides generalization from low-order mutants to higher-order mutants. We show that ECNet predicts the sequence-function relationship more accurately as compared to existing machine learning algorithms by using ~50 deep mutational scanning and random mutagenesis datasets. Moreover, we used ECNet to guide the engineering of TEM-1 β-lactamase and identified variants with improved ampicillin resistance with high success rates.
59 BASIC BIOLOGICAL SCIENCES↗
Physics-constrained deep learning of nonlinear normal modes of spatiotemporal fluid flow dynamics
In this study, we present a physics-constrained deep learning method to discover and visualize from data the invariant nonlinear normal modes (NNMs) which contain the spatiotemporal dynamics of the fluid flow potentially containing strong nonlinearity. Specifically, we develop a NNM-physics-constrained convolutional autoencoder (NNM-CNN-AE) integrated with a multi-temporal-step dynamics prediction block to learn the nonlinear modal transformation, the NNMs containing the spatiotemporal dynamics of the flow, and reduced-order reconstruction and long-time future-state prediction of the flow fields, simultaneously. In test cases, we apply the developed method to analyze different flow regimes past a cylinder, including laminar flows with low Reynolds number in transient and steady states (RD = 100) and high Reynolds number flow (RD = 1000), respectively. The results indicate that the identified NNMs are able to reveal the nonlinear spatiotemporal dynamics of these flows, and the NNMs-based reduced-order modeling consistently achieves better accuracy with orders of magnitudes smaller errors in construction and prediction of the nonlinear velocity and vorticity fields, compared to the linear proper orthogonal decomposition (POD) method and the Koopman-constrained-CNN-AE using the same number or dimension of modes. We perform an analysis of the modal energy distribution of NNMs and find that compared to POD modes, the few fundamental NNMs capture a very high level of total energy of the flow, which is advantageous for reduced-order modeling and representation of the complex flows. Finally, we discuss the potentials and limitations of the presented method.
Mechanics↗
Emulating ab initio computations of infinite nucleonic matter
We construct efficient emulators for the computation of the infinite nuclear matter equation of state. These emulators are based on the subspace-projected coupled-cluster method for which we here develop a new algorithm called small-batch voting to eliminate spurious states that might appear when emulating quantum many-body methods based on a non-Hermitian Hamiltonian. The efficiency and accuracy of these emulators facilitate a rigorous statistical analysis within which we explore nuclear matter predictions for > 10 6 different parametrizations of a chiral interaction model with explicit Δ -isobars at next-to-next-to leading order. Constrained by nucleon-nucleon scattering phase shifts and bound-state observables of light nuclei up to He 4 , we use history matching to identify nonimplausible domains for the low-energy coupling constants of the chiral interaction. Within these domains we perform a Bayesian analysis using sampling and importance resampling with different likelihood calibrations and study correlations between interaction parameters, calibration observables in light nuclei, and nuclear matter saturation properties. Published by the American Physical Society 2024
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Low Mach number fluctuating hydrodynamics model for ionic liquids
We present a new mesoscale model for ionic liquids based on a low Mach number fluctuating hydrodynamics formulation for multicomponent charged species. The low Mach number approach eliminates sound waves from the fully compressible equations leading to a computationally efficient incompressible formulation. The model uses a Gibbs free energy functional that includes enthalpy of mixing, interfacial energy, and electrostatic contributions. These lead to a new fourth-order term in the mass equations and a reversible stress in the momentum equations. We calibrate our model using parameters for [DMPI+][F6P-], an extensively-studied room temperature ionic liquid (RTIL), and numerically demonstrate the formation of mesoscopic structuring at equilibrium in two and three dimensions. In simulations with electrode boundaries the measured double layer capacitance decreases with voltage, in agreement with theoretical predictions and experimental measurements for RTILs. Finally, we present a shear electroosmosis example to demonstrate that the methodology can be used to model electrokinetic flows.
42 ENGINEERING↗