Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical Physics”

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 973 records · Page 54

Fast strength model characterization using Bayesian statistics

A variety of flow stress models exist with new models constantly being developed. These models aim to approximate the strength of materials in a variety of regimes from quasistatic loading through shock scenarios. All models contain an array of parameters which need to be tuned to the material under study. Some models perform well under limited conditions, requiring adjustment of the parameters when venturing outside of those predefined ranges. Other models perform well over a wide range of conditions with a set of parameters, but may be outperformed by other models optimized on a tighter range of conditions. Recent research by Los Alamos demonstrated the ability to optimize the Johnson Cook (JC) model using a set of 3 plate-impact experiments on Aluminum. They utilized Bayesian statistics and emulation to determine optimal parameters for the model with a quantification of parameter uncertainty. We present an extension of this capability to incorporate velocimetry from plate-impact tests, stress-strain data from split Hopkinson pressure bar and quasistatic compression tests, plus profiles from Taylor cylinders in a unified fashion. Statistically robust comparisons of the performance and uncertainty of different realizations of the JC flow stress model were carried out based on calibration to several possible combinations of these three different experiment types.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Evaluation of statistical variation of microstructural properties and temperature effects on creep fracture of Grade 91

This report describes three improvements made to a physically-based model for creep and creep rupture in Grade 91 steel developed as part of the Advanced Reactor Technologies program: 1. Progress on transitioning the model framework to the MOOSE finite element package to improve parallel scalability, increase the physical size of the simulations, and incorporate multiphysics effects into the model. 2. The development of a physically-based method of scaling the model parameters to accurately capture creep in Grade 91 for temperatures in the range of 450° to 500° C. 3. Extending the model to account for microstructural statistical variations, specifically capturing the effect of grain boundary energy on the physical parameters underling the model for grain boundary void nucleation and cavitation. Put together, these improvements result in an accurate, predictive model for the physical response of Grade 91 steel over the expected use temperature range for the material in future liquid metal cooled fast reactors. The model can be used to more accurately predict engineering properties of the material for very long service lives, which could lead to safer, more economical future advanced reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Impact of horizontal resolution on simulation of precipitation extremes in an aqua-planet version of Community Atmospheric Model (CAM3)

One key question regarding current climate models is whether the projection of climate extremes converges to a realistic representation as the spatial and temporal resolutions of the model are increased. Ideally the model extreme statistics should approach a fixed distribution once the resolutions are commensurate with the characteristic length and time scales of the processes governing the formation of the extreme phenomena of interest. In this study, a series of AGCM runs with idealized 'aquaplanet-steady-state' boundary conditions have been performed with the Community Atmosphere Model CAM3 to investigate the effect of horizontal resolution on climate extreme simulations. The use of the aquaplanet framework highlights the roles of model physics and dynamics and removes any apparent convergence in extreme statistics due to better resolution of surface boundary conditions and other external inputs. Assessed at a same large spatial scale, the results show that the horizontal resolution and time step have strong effects on the simulations of precipitation extremes. The horizontal resolution has a much stronger impact on precipitation extremes than on mean precipitation. Updrafts are strongly correlated with extreme precipitation at tropics at all the resolutions, while positive low-tropospheric temperature anomalies are associated with extreme precipitation at mid-latitudes.

54 ENVIRONMENTAL SCIENCES↗

Peer-to-Peer Energy Trading and Grid Impact Studies in Smart Communities

The rise of peer-to-peer (P2P) marketplace paradigms has transformed existing marketplace models, but the extent to which this approach can be applied to the energy marketplace has yet to be considered. In this paper, we examine existing approaches taken in the application of a P2P paradigm to the energy marketplace, further presenting an approach towards facilitating an online P2P energy marketplace, implementing a prototype P2P web application named SolTrade. Furthermore, we submit initial statistics based on simulated transactions facilitated through the platform, which illustrate the physical impact of marketplace transactions on the energy grid. In particular, these results show that, as the number of users rises, the chance of overloading the grid rises, but the chance of the grid being unable to sustain itself without an external source of energy falls.

peer-to-peer economy, web development, energy shar↗

Path integral Monte Carlo approach to the structural properties and collective excitations of liquid $$^3{\text {He}}$$ without fixed nodes

Abstract Due to its nature as a strongly correlated quantum liquid, ultracold helium is characterized by the nontrivial interplay of different physical effects. Bosonic $$^4{\text {He}}$$ 4 He exhibits superfluidity and Bose-Einstein condensation. Its physical properties have been accurately determined on the basis of ab initio path integral Monte Carlo (PIMC) simulations. In contrast, the corresponding theoretical description of fermionic $$^3{\text {He}}$$ 3 He is severely hampered by the notorious fermion sign problem, and previous PIMC results have been derived by introducing the uncontrolled fixed-node approximation. In this work, we present extensive new PIMC simulations of normal liquid $$^3{\text {He}}$$ 3 He without any nodal constraints. This allows us to to unambiguously quantify the impact of Fermi statistics and to study the effects of temperature on different physical properties like the static structure factor $$S({\mathbf {q}})$$ S ( q ) , the momentum distribution $$n({\mathbf {q}})$$ n ( q ) , and the static density response function $$\chi ({\mathbf {q}})$$ χ ( q ) . In addition, the dynamic structure factor $$S({\mathbf {q}},\omega )$$ S ( q , ω ) is rigorously reconstructed from imaginary-time PIMC data. From simulations of $$^3{\text {He}}$$ 3 He , we derived the familiar phonon–maxon–roton dispersion function that is well-known for $$^4{\text {He}}$$ 4 He and has been reported previously for two-dimensional $$^3{\text {He}}$$ 3 He films (Nature 483:576–579 (2012)). The comparison of our new results for both $$S({\mathbf {q}})$$ S ( q ) and $$S({\mathbf {q}},\omega )$$ S ( q , ω ) with neutron scattering measurements reveals an excellent agreement between theory and experiment.

97 MATHEMATICS AND COMPUTING↗

PROPERTIES OF SUPERSONIC EVERSHED DOWNFLOWS

We study supersonic Evershed downflows in a sunspot penumbra by means of high spatial resolution spectropolarimetric data acquired in the Fe i 617.3 nm line with the CRISP instrument at the Swedish 1 m Solar Telescope. Physical observables, such as Dopplergrams calculated from line bisectors and Stokes V zero-crossing wavelengths, and Stokes V maps in the far red-wing, are used to find regions where supersonic Evershed downflows may exist. We retrieve the line-of-sight velocity and the magnetic field vector in these regions using two-component inversions of the observed Stokes profiles with the help of the SIR code. We follow these regions during their lifetime to study their temporal behavior. Finally, we carry out a statistical analysis of the detected supersonic downflows to characterize their physical properties. Supersonic downflows are contained in compact patches moving outward, which are located in the mid- and outer penumbra. They are observed as bright, roundish structures at the outer end of penumbral filaments that resemble penumbral grains. The patches may undergo fragmentations and mergings during their lifetime; some of them are recurrent. Supersonic downflows are associated with strong and rather vertical magnetic fields with a reversed polarity compared to that of the sunspot. Our results suggest that downflows returning back to the solar surface with supersonic velocities are abruptly stopped in dense deep layers and produce a shock. Consequently, this shock enhances the temperature and is detected as a bright grain in the continuum filtergrams, which could explain the existence of outward-moving grains in the mid- and outer penumbra.

79 ASTRONOMY AND ASTROPHYSICS↗

Secondary carbides in high chromium cast irons: An alternative approach to their morphological and spatial distribution characterization

Highlights: • Systematic characterization of second phase particles by combining a variety of techniques • Comparison of morphological parameters obtained from SEM, BSE, etching+CLSM and EBSD • Quantitative analysis and degree of distribution of secondary carbides • Particle size, distribution and morphology from image analysis of CLSM, SEM and EBS images - Abstract: Secondary carbides precipitated in high chromium cast irons during thermal treatments were characterized by means of different characterization techniques, including scanning electron microscopy, energy dispersive X-ray spectroscopy, electron backscattered diffraction and a combination of chemical etching with confocal scanning laser microscopy. This set of techniques provides a full morphological, chemical and crystallographic description of the analysed phases. This work evaluated different methods for optimizing the image acquisition for a further image analysis (IA) based on the threshold binarization. Finally, the carbide size, distribution and morphology were determined after IA of the images acquired by aforementioned characterization techniques. Although the different techniques report some dispersion in the value for the average particle size, the particle inter-spacing and aspect ratio meet within the error value. The proposed characterization methodology provides statistically reliable data for a further evaluation of related physical properties in composites.

36 MATERIALS SCIENCE↗

Phase transitions in few-monolayer spin ice films

Vertex models are an important class of statistical mechanical system that admit exact solutions and exotic physics. Applications include water ice, ferro- and antiferro-electrics, spin ice and artificial spin ice. Here we show that it is possible to engineer spin ice films with atomic-layer precision down to the monolayer limit. Specific heat measurements show that these films, which have a fundamentally different symmetry to bulk spin ice, realise systems close to the two-dimensional F-model, with exotic phase transitions on topologicallyconstrained configurational manifolds. Our results show how spin ice thin films can release the celebrated Pauling entropy of spin ice without an anomaly in the specific heat. They also significantly expand the class of vertex models available to experiment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Hydrodynamic analysis of laser-driven cylindrical implosions

Three-dimensional hydrodynamic simulations are performed to study laser-driven cylindrical implosions in the context of experiments (F. Perez et al., Plasma Phys. Controlled Fusion 51, 124035 (2009)) carried out at the Rutherford Appleton Laboratory in the framework of the HiPER project. The analysis is carried out by using the 3D version of the hydrocode MULTI (R. Ramis et al., Comput. Phys. Commun. 49, 475-505 (1988)). The influence of the main laser parameters on implosion performance and symmetry is consistently studied and compared with the results of 2D analysis. Furthermore, the effects of uncertainties in laser irradiation (pointing, focusing, power balance, and time jitter) on implosion performance (average peak density and temperature) are studied by means of statistical analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Understanding glass-forming ability through sluggish crystallization of atomically thin metallic glassy films

The glass-forming ability (GFA) of an alloy, closely related to its ability to resist crystallization, is a crucial issue in condensed matter physics. So far, the studies on GFA are mostly statistical and empirical guides. Benefiting from the ultrahigh thermal stability of ultrathin metallic glassy film and high resolution spherical aberration-corrected transmission electron microscope, the crystallization of atomically thin ZrCu and its microalloyed ZrCuAl glasses with markedly different GFA was investigated at the atomic scale. We find the Zr diffusivity estimated from the density of nuclei is dramatically decreased by adding of Al, which is the major reason for the much better GFA of the ZrCuAl metallic glass.

36 MATERIALS SCIENCE↗

A tale of two modes: neutrino free-streaming in the early universe

We present updated constraints on the free-streaming nature of cosmological neutrinos from cosmic microwave background (CMB) temperature and polarization power spectra, baryonic acoustic oscillation data, and distance ladder measurements of the Hubble constant. Specifically, we consider a Fermi-like four-fermion interaction between massless neutrinos, characterized by an effective coupling constant G {sub eff}, and resulting in a neutrino opacity τ-dot {sub ν∝} G {sub eff}{sup 2} T {sub ν}{sup 5}. Using a conservative flat prior on the parameter log{sub 10}( G {sub eff} MeV{sup 2}), we find a bimodal posterior distribution with two clearly separated regions of high probability. The first of these modes is consistent with the standard ΛCDM cosmology and corresponds to neutrinos decoupling at redshift z {sub ν,dec} > 1.3×10{sup 5}, that is before the Fourier modes probed by the CMB damping tail enter the causal horizon. The other mode of the posterior, dubbed the 'interacting neutrino mode', corresponds to neutrino decoupling occurring within a narrow redshift window centered around z {sub ν,dec}∼8300. This mode is characterized by a high value of the effective neutrino coupling constant, log{sub 10}( G {sub eff} MeV{sup 2}) = −1.72 ± 0.10 (68% C.L.), together with a lower value of the scalar spectral index and amplitude of fluctuations, and a higher value of the Hubble parameter. Using both a maximum likelihood analysis and the ratio of the two mode's Bayesian evidence, we find the interacting neutrino mode to be statistically disfavored compared to the standard ΛCDM cosmology, and determine this result to be largely driven by the low- l CMB temperature data. Interestingly, the addition of CMB polarization and direct Hubble constant measurements significantly raises the statistical significance of this secondary mode, indicating that new physics in the neutrino sector could help explain the difference between local measurements of H {sub 0}, and those inferred from CMB data. A robust consequence of our results is that neutrinos must be free streaming long before the epoch of matter-radiation equality in order to fit current cosmological data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Good and proper: self-similarity of N -body simulations with proper force softening

ABSTRACT Analysis of self-similarity in scale-free N-body simulations reveals the spatial and temporal scales for which statistics measured in cosmological simulations are converged to the physical continuum limit. We examine how the range of scales in which the two-point correlation function is converged depends on the force softening length and whether it is held constant in comoving or proper coordinates. We find that a proper softening that reaches roughly 1/30th of the inter-particle spacing by the end of the simulation resolves the same spatial and temporal scales as a comoving softening of the same length while using a third fewer time-steps, for a range of scale factors typical to Lambda cold dark matter (ΛCDM) simulations. We additionally infer an inherent resolution limit, set by the particle mass and scaling as a−1/2, beyond which reducing the softening does not improve the resolution. We postulate a mapping of these results with spectral index n = −2 to ΛCDM simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Searching for axionlike particles from core-collapse supernovae with Fermi LAT’s low-energy technique

Light axionlike particles (ALPs) are expected to be abundantly produced in core-collapse supernovae (CCSNe), resulting in a ~ 10 -second long burst of ALPs. These particles subsequently undergo conversion into gamma rays in external magnetic fields to produce a long gamma-ray burst (GRB) with a characteristic spectrum peaking in the 30–100-MeV energy range. At the same time, CCSNe are invoked as progenitors of ordinary long GRBs, rendering it relevant to conduct a comprehensive search for ALP spectral signatures using the observations of long GRBs with the Fermi Large Area Telescope (LAT). We perform a data-driven sensitivity analysis to determine CCSN distances for which a detection of an ALP signal is possible with the LAT’s low-energy technique which, in contrast to the standard LAT analysis, allows for a a larger effective area for energies down to 30 MeV. Assuming an ALP mass m a ≲ 10 - 10 eV and ALP-photon coupling g a γ = 5.3 × 10 - 12 GeV - 1 , values considered and deduced in ALP searches from SN1987A, we find that the distance limit ranges from ~ 0.5 to ~ 10 Mpc , depending on the sky location and the CCSN progenitor mass. Furthermore, we select a candidate sample of 24 GRBs and carry out a model comparison analysis in which we consider different GRB spectral models with and without an ALP signal component. We find that the inclusion of an ALP contribution does not result in any statistically significant improvement of the fits to the data. Finally, we discuss the statistical method used in our analysis and the underlying physical assumptions, the feasibility of setting upper limits on the ALP-photon coupling, and give an outlook on future telescopes in the context of ALP searches.

79 ASTRONOMY AND ASTROPHYSICS↗

Improving neutrino energy estimation of charged-current interaction events with recurrent neural networks in MicroBooNE

We present a deep learning-based method for estimating the neutrino energy of charged-current neutrino-argon interactions. We employ a recurrent neural network (RNN) architecture for neutrino energy estimation in the MicroBooNE experiment, utilizing liquid argon time projection chamber (LArTPC) detector technology. Traditional energy estimation approaches in LArTPCs, which largely rely on reconstructing and summing visible energies, often experience sizable biases and resolution smearing because of the complex nature of neutrino interactions and the detector response. The estimation of neutrino energy can be improved after considering the kinematics information of reconstructed final-state particles. Utilizing kinematic information of reconstructed particles, the deep learning-based approach shows improved resolution and reduced bias for the muon neutrino Monte Carlo simulation sample compared to the traditional approach. In order to address the common concern about the effectiveness of this method on experimental data, the RNN-based energy estimator is further examined and validated with dedicated data-simulation consistency tests using MicroBooNE data. We also assess its potential impact on a neutrino oscillation study after accounting for all statistical and systematic uncertainties and show that it enhances physics sensitivity. This method has good potential to improve the performance of other physics analyses. Published by the American Physical Society 2024

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Phenomenology of large scale structure in scalar-tensor theories: Joint prior covariance of w DE , Σ, and μ in Horndeski theories

Ongoing and upcoming cosmological surveys will significantly improve our ability to probe the equation of state of dark energy, w DE , and the phenomenology of large scale structure. They will allow us to constrain deviations from the Λ cold dark matter predictions for the relations between the matter density contrast and the weak lensing and the Newtonian potential, described by the functions Σ and μ, respectively. In this work, we derive the theoretical prior for the joint covariance of w DE , Σ and μ, expected in general scalar-tensor theories with second order equations of motion (Horndeski gravity), focusing on their time-dependence at certain representative scales. We employ Monte Carlo methods to generate large ensembles of statistically independent Horndeski models, focusing on those that are physically viable and in broad agreement with local tests of gravity, the observed cosmic expansion history and the measurement of the speed of gravitational waves from a binary neutron star merger. Here, we identify several interesting features and trends in the distribution functions of w DE , Σ and μ, as well as in their covariances; we confirm the high degree of correlation between Σ and μ in scalar-tensor theories. The derived prior covariance matrices will allow us to reconstruct jointly w DE , Σ and μ in a nonparametric way.

79 ASTRONOMY AND ASTROPHYSICS↗

Supersymmetric Dynamics and Zeta-Functions

Boson, fermion, and super oscillators and (statistical) mechanism of cosmological constant; finite approximation of the zeta-function and fermion factorization of the bosonic statistical sum concidered.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Rapid Electrochemical Diagnosis of Battery Health and Safety from Cells to Modules

Rapid electrochemical diagnosis of battery health and failure is critical for ensuring reliable battery performance and battery safety. Traditional battery health diagnostics such as capacity measurements and DC pulse tests are reliable and well-understood, however, these measurements of battery capacity and resistance do not capture all aspects of battery degradation. Other aspects of degradation, such as electrolyte decomposition, lithium-plating, and particle cracking are difficult to detect electrochemically but are crucial to measure to get a full picture of battery safety and flag out potential failures. In this work, lab- and field-aged commercial lithium-ion batteries and modules of various chemistries and formats are tested using a variety of traditional electrochemical characterization methods as well as using 2-minute pseudo-random DC pulse sequences at rest and during charge/discharge. The electrochemical measurements are compared to physical cell measurements, cell efficiency, drive cycle performance, physical and thermal heterogeneity, and qualitative safety metrics using statistical and machine-learning methods to discover if a comprehensive "battery health map" can be accurately identified using only rapid DC measurements.

ADVANCED PROPULSION SYSTEMS,ENERGY STORAGE↗

Robust parallel iterative solvers for linear and least-squares problems (Final Technical Report)

The primary goal of this project is to study and develop robust iterative methods for solving linear systems of equations and least squares systems. The focus of the Minnesota team is on algorithms development, robustness issues, and on tests and validation of the methods on realistic problems. 1. The project begun with an investigation on how to practically update a preconditioner obtained from an ILU-type factorization, when the coefficient matrix changes. 2. We investigated strategies to improve robustness in parallel preconditioners in a specific case of a PDE with discontinuous coefficients. 3. We explored ways to adapt standard preconditioners for solving linear systems arising from the Helmholtz equation. These are often difficult linear systems to solve by iterative methods. 4. We have also worked on purely theoretical issues related to the analysis of Krylov subspace methods for linear systems. 5. We developed an effective strategy for performing ILU factorizations for the case when the matrix is highly indefinite. The strategy uses shifting in some optimal way. The method was extended to the solution of Helmholtz equations by using complex shifts, yielding very good results in many cases. 6. We addressed the difficult problem of preconditioning sparse systems of equations on GPUs. 7. A by-product of the above work is a software package consisting of an iterative solver library for GPUs based on CUDA. This was made publicly available. It was the first such library that offers complete iterative solvers for GPUs. 8. We considered another form of ILU which blends coarsening techniques from Multigrid with algebraic multilevel methods. 9. We have released a new version on our parallel solver - called pARMS [new version is version 3]. As part of this we have tested the code in complex settings - including the solution of Maxwell and Helmholtz equations and for a problem of crystal growth.10. As an application of polynomial preconditioning we considered the problem of evaluating f(A)v which arises in statistical sampling. 11. As an application to the methods we developed, we tackled the problem of computing the diagonal of the inverse of a matrix. This arises in statistical applications as well as in many applications in physics. We explored probing methods as well as domain-decomposition type methods. 12. A collaboration with researchers from Toulouse, France, considered the important problem of computing the Schur complement in a domain-decomposition approach. 13. We explored new ways of preconditioning linear systems, based on low-rank approximations.

97 MATHEMATICS AND COMPUTING↗