Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical mechanics models”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Theorem on extensive spectral degeneracy for systems with rigid higher symmetries in general dimensions

We establish, in the spirit of the Lieb-Schultz-Mattis theorem, lower bounds on the spectral degeneracy of quantum systems with higher (gaugelike) symmetries with rather generic physical boundary conditions in an arbitrary number of spatial dimensions. Contrary to applying twists or equivalent adiabatic operations, we exploit the effects of modified boundary conditions. We find that when a general choice of boundary geometry is immaterial in determining spectral degeneracies while approaching the thermodynamic limit, systems that exhibit rigid noncommuting gaugelike symmetries, such as the orbital compass model, must have an exponential (in the size of the boundary) degeneracy of each of their spectral levels.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1986

Theoretical analysis of isotope effects on ozone formation in oxygen photochemistry

In situ measurements of stratospheric ozone and laboratory studies of ozone production in electric discharge through oxygen have shown previously that ozone containing heavy isotopes of oxygen (O-17, O-18) may be formed preferentially. In order to assess the relevance of thee latter experiment to the stratospheric measurements, detailed understanding of the effect of isotopic substitution on the O3 formation reaction O + O2 + M yields O3 + M and on the O atom exchange reaction O + O2 + O yields O2 + O is necessary. In this work, an estimate of the effect of isotopic substitution on the recombination rate is made by us of approximate dynamical theories and statistical mechanics. The results indicate the possibility of isotope effects on the O + O2 recombination rate of the order of several percent at stratospheric temperatures. In general, recombination reactions involving heavy (mass 49, 50) O3 formation are found to be slower than the reaction leading to normal (mass 48) O3 formation. The calculated isotope effects are sufficiently small that the uncertainties in the model input and the approximations in the dynamical theories will probably make the quantitative nature of these results subject to considerable uncertainty. This isotope effect should not be observable in the atmosphere given the precision of the current measurements but may be crucial in the understanding of the laboratory experiments, where observed enhancements are only of the order of several percent. Possible reasons for this discrepancy between the observed enhancement and predicted depletion are presented.

Kaye, J. A.↗
NASA NTRS1987

Advanced probabilistic method of development

Advanced structural reliability methods are utilized on the Probabilistic Structural Analysis Methods (PSAM) project to provide a tool for analysis and design of space propulsion system hardware. The role of the effort at the University of Arizona is to provide reliability technology support to this project. PSAM computer programs will provide a design tool for analyzing uncertainty associated with thermal and mechanical loading, material behavior, geometry, and the analysis methods used. Specifically, reliability methods are employed to perform sensitivity analyses, to establish the distribution of a critical response variable (e.g., stress, deflection), to perform reliability assessment, and ultimately to produce a design which will minimize cost and/or weight. Uncertainties in the design factors of space propulsion hardware are described by probability models constructed using statistical analysis of data. Statistical methods are employed to produce a probability model, i.e., a statistical synthesis or summary of each design variable in a format suitable for reliability analysis and ultimately, design decisions.

Wirsching, P. H.↗
DOE OSTI2022

Reduced Dynamics of Full Counting Statistics

We present a theory of modified reduced dynamics in the presence of counting fields. Reduced dynamics techniques are useful for describing open quantum systems at long emergent timescales when the memory timescales are short. However, they can be difficult to formulate for observables spanning the system and its environment, such as those characterizing transport properties. A large variety of mixed system--environment observables, as well as their statistical properties, can be evaluated by considering counting fields. Given a numerical method able to simulate the field-modified dynamics over the memory timescale, we show that the long-lived full counting statistics can be efficiently obtained from the reduced dynamics. We demonstrate the utility of the technique by computing the long-time current in the nonequilibrium Anderson impurity model from short-time Monte Carlo simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Toward a 3d Ising model with a weakly-coupled string theory dual

It has long been expected that the 3d Ising model can be thought of as a string theory, where one interprets the domain walls that separate up spins from down spins as two-dimensional string worldsheets. The usual Ising Hamiltonian measures the area of these domain walls. This theory has string coupling of unit magnitude. We add new local terms to the Ising Hamiltonian that further weight each spin configuration by a factor depending on the genus of the corresponding domain wall, resulting in a new 3d Ising model that has a tunable bare string coupling g_s g s . We use a combination of analytical and numerical methods to analyze the phase structure of this model as g_s g s is varied. We study statistical properties of the topology of worldsheets and discuss the prospects of using this new deformation at weak string coupling to find a worldsheet description of the 3d Ising transition.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRSDate not supplied

Microscale Constitutive Model Sensitivity on Multiscale Modeling of Fiber Reinforced Composites

Fiber reinforced composites are desirable in applications where low weight and high strength are needed, but are susceptible to microscale variability during manufacturing, making failure predictions difficult. The impact microscale variability has on macroscale mechanical response is difficult to predict due to the computational efficiency needed to simulate many, large, high fidelity, microscale models. In this study, a multiscale approach was taken to model 3-point bend, 4-point bend, and tensile experiments of a unidirectional composite from only having microstructure scans of these samples and constituent properties from literature. These scans were sampled with different sized windows, and statistically equivalent microstructures were generated, then simulated for stiffness, strength, and fracture toughness using an efficient micromechanical. Mesoscale models were created where element sizes equaled microstructure size, and properties were assigned through sampling of microscale simulation results. First, this study showed the effect of using Weibull scaling on constituent matrix strength on macroscale response. Then, a comparison was made between different element sizes and experiments. Finally, model dimensions were fixed, and the effect of randomly distributed local properties alone was examined. Results showed that the scatter of strength and stiffness in the experiments could be predicted well using images of the microscale fiber morphologies and that using stochastic properties produced a 3% coefficient of variation of strength for all experiments.

statistical microstructure↗
DOE OSTI2023

COMPUTATIONAL MODELING OF IGNITION AND PREMIXED FLAME PROPAGATION INITIATED BY A PRE-CHAMBER TURBULENT JET

Addressing the pressing need for reduced carbon emissions, Turbulent Jet Ignition (TJI) emerges as a promising technology for ultra-lean combustion, offering enhanced thermal efficiencies and minimized cyclic variability in spark-ignited engines. To facilitate rapid testing and integration of this technology, a robust computational modeling framework is crucial. This study delves into the predictive capabilities of computational models for main-chamber ignition and premixed flame propagation using a single-cycle TJI rig measured by Biswas et al. (Applied Thermal Engineering, vol 106, 2016). Employing an open-source compressible flow simulation solver with Large Eddy Simulation (LES) for turbulence modeling, the investigation integrates the conventional Laminar Finite Rate Chemistry (LFRC) model alongside the transported Probability Density Method (PDF) for turbulence-chemistry interaction. A fully-consistent Eulerian Monte-Carlo Fields (EMCF) method is utilized to approximate the transported PDF, while Interaction by Exchange with Mean is employed to close micro-mixing terms in stochastic differential equations. A reduced chemical reaction mechanism with 21 species and 84 reactions (DRM-19) is used for solving chemical kinetics, and a double Gaussian energy deposition model is used to approximate the spark ignition in the pre-chamber. An unstructured O-grid mesh with 0.3 million cells in the prechamber and 1 million cells in the main chamber is employed. Results are divided into two phases: pre-chamber initialization and full TJI simulations. Validation of the predicted pre-chamber flame propagation and the lean ignition in the main-chamber is carried out by using available experimental data. Under quiescent conditions, both the LFRC and transported PDF methods largely underestimate the flame speed and subsequent pressure growth in the pre-chamber. A linear momentum forcing technique is applied to investigate the impact of initial turbulence in the pre-chamber, demonstrating a notable influence on flame propagation. Fine-tuning of the forcing coefficient reproduces the sudden pressure growth observed in the experiment. The experimentally validated pre-chamber simulation serves as the initial condition for the full TJI simulations. It is found that the LFRC model fails to predict lean-ignition in the main-chamber, resulting in a misfiring event. Incorporation of turbulence-chemistry interaction using the transported PDF method substantially improves the prediction of the ignition event in the main-chamber, achieving fair qualitative agreement and quantitative validation of combustion parameters within ±10% of the reported experimental data. The rich simulation results consisting of a full set of statistical description of the thermo-chemical states enable us to gain deep insights into the ignition mechanisms in the main chamber, which is limited when done experimentally. A novel dual ignition phenomenon is revealed in the TJI rig for the first time. Initially, a primary ignition kernel is formed at a downstream location which eventually detaches from the main jet. As the jet momentum decreases, a secondary ignition event follows, this time at a more upstream location which eventually combines with the primary ignition kernel to form a single connected flame front. Investigation of these ignition sequences in chemical composition space reveal distinct differences between the two. The primary ignition event in the main-chamber is followed by a large concentration of active radicals from the pre-chamber jet, accelerating the chain-branching steps, characterizing what has been referred to as flame ignition. In contrast, the secondary ignition occurs in the absence of active radicals in the pre-chamber jet, hence characterized as jet ignition. Further analysis of the effect of pre-chamber jet characteristics on lean ignition in the main-chamber is conducted by setting up cases with different initial pressure ratios (por) between the two chambers, a non-dimensional parameter, ranging from 1.2 to 3.2. As the initial pressure ratio increases, jet momentum increases, with dual ignition observed in cases above por= 2.2. Case with por= 3.2 lead to misfiring. The effect of ignition sequence on global combustion characteristics of TJI is analyzed. Dual ignition events lead to non-monotonicity in combustion characteristics such as global reaction progress variable, flame penetration, and global heat release rate. In dual ignition events, although the rate of fuel consumption and global heat release rate is initially lower, the secondary ignition leads to a sudden increase in flame surface area, resulting in a sudden jump and promoting the overall performance of the TJI system.

42 ENGINEERING↗
DOE OSTI2024

Robust finite-temperature many-body scarring on a quantum computer

Mechanisms for suppressing thermalization in disorder-free many-body systems, such as Hilbert space fragmentation and quantum many-body scars, have recently attracted much interest in foundations of quantum statistical physics and potential quantum information processing applications. However, their sensitivity to realistic effects such as finite temperature remains largely unexplored. Here, we have utilized IBM's Kolkata quantum processor to demonstrate an unexpected robustness of quantum many-body scars at finite temperatures when the system is prepared in a thermal Gibbs ensemble. We identify such robustness in the PXP model, which describes quantum many-body scars in experimental systems of Rydberg atom arrays and ultracold atoms in tilted Bose-Hubbard optical lattices. By contrast, other theoretical models which host exact quantum many-body scars are found to lack such robustness and their scarring properties quickly decay with temperature. Our study sheds light on the important differences between scarred models in terms of their algebraic structures, which impacts their resilience to finite temperature. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Zentropy Theory for Transformative Functionalities of Magnetic and Superconducting Materials

The proposed research developed the zentropy theory through applications to complex magnetic materials and superconductors under the hypothesis that the emergent properties of complex magnetic materials and superconductors can be predicted by statistical mechanics of ergodic microstates with their partition functions computed from DFT-predicted free energies. The key objective is to develop approaches to systematically determine the types and number of microstates and the supercell size in DFT-based calculations through convergency of macroscopic functionalities, with the incorporation of our mixed-space approach accounting for the interactions between periodic supercells. In addition to use scientific intuitions to guide the design of important microstates, the key innovation of the proposed research is to integrate the domain knowledge and the material-property-descriptor database (MPDD) with 4 million microstates, which is supported by our deep neural network machine learning models (SIPFENN: structure-informed prediction of formation energy using neural networks) and integrated with our high throughput DFT Tool Kit (DFTTK). For complex magnetic materials, one of the objectives is to develop approaches to calculate short-range ordering from the statistical distribution of each microstate. For superconductors, the divergency of quasiparticle effective mass at a quantum critical point will be investigated, and the superconducting and non-superconducting microstates will be delineated through analysis of electronic band structure, density of states, charge density, and Fermi surface.

36 MATERIALS SCIENCE↗
NASA NTRS1987

Cyclic hardening mechanisms in Nimonic 80A

A nickel base superalloy was fatigued under constant plastic strain range control. The hardening response was investigated as a function of plastic strain range and particle size of the gamma prime phase. Hardening was found to be a function of the slip band spacing. Numerous measurements of the slip band spacing and other statistical data on the slip band structures were obtained. Interactions between intersecting slip systems were shown to influence hardening. A Petch-Hall model was found to describe best this relationship between the response stress and the slip band spacing.

Lerch, B. A.↗
DOE OSTI2022

Quantum wake dynamics in Heisenberg antiferromagnetic chains

Traditional spectroscopy, by its very nature, characterizes physical system properties in the momentum and frequency domains. However, the most interesting and potentially practically useful quantum many-body effects emerge from local, short-time correlations. Here, using inelastic neutron scattering and methods of integrability, we experimentally observe and theoretically describe a local, coherent, long-lived, quasiperiodically oscillating magnetic state emerging out of the distillation of propagating excitations following a local quantum quench in a Heisenberg antiferromagnetic chain. This “quantum wake” displays similarities to Floquet states, discrete time crystals and nonlinear Luttinger liquids. We also show how this technique reveals the non-commutativity of spin operators, and is thus a model-agnostic measure of a magnetic system’s “quantumness.”

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1993

Development of a model to assess orthostatic responses

A major change for crewmembers during weightlessness in microgravity is the redistribution of body fluids from the legs into the abdomen, thorax, and head. The fluids continue to be sequestered in these areas throughout the flight. Upon reentry into gravity on landing, these same body fluids are displaced again to their normal locations, however, not without hazardous incidence to the crewmembers. The problem remains that upon landing, crewmembers are subject to orthostasis, that is, the blood flowing into the legs reduces the blood supply to the brain and may result in the crewmember fainting. The purpose of this study was to develop a model of testing orthostatic responses of blood pressure regulating mechanisms of the cardiovascular system, when challenged, to maintain blood pressure to the brain. To accomplish this, subjects' responses were assessed as they proceeded from the supine position of progressive head-up tilt positions of 30 deg, 60 deg, and 90 deg angles. A convenience sample consisted of 21 subjects, females (N=11) and males (N=10), selected from a list of potential subjects available through the NASA subject screening office. The methodology included all non-invasive measurements of blood pressure, heart rate, echocardiograms, cardiac output, cardiac stroke volume, fluid shifts in the thorax, ventricular ejection and velocity times, and skin blood perfusion. The Fischer statistical analysis was done of all data with the significance level at .05. Significant differences were demonstrated in many instances of changes of posture for all variables. Based on the significance of the findings of this study, this model for assessing orthostatic responses does provide an adequate challenge to the blood pressure regulatory systems. While individuals may use different adaptations to incremental changes in gravity, the subjects, in aggregate, demonstrated significant adaptive cardiovascular changes to orthostatic challenges which were presented to them.

Rubin, Marilyn↗
DOE OSTI2021

A model for the fragmentation kinetics of crumpled thin sheets

As a confined thin sheet crumples, it spontaneously segments into flat facets delimited by a network of ridges. Despite the apparent disorder of this process, statistical properties of crumpled sheets exhibit striking reproducibility. Experiments have shown that the total crease length accrues logarithmically when repeatedly compacting and unfolding a sheet of paper. Here, we offer insight to this unexpected result by exploring the correspondence between crumpling and fragmentation processes. We identify a physical model for the evolution of facet area and ridge length distributions of crumpled sheets, and propose a mechanism for re-fragmentation driven by geometric frustration. This mechanism establishes a feedback loop in which the facet size distribution informs the subsequent rate of fragmentation under repeated confinement, thereby producing a new size distribution. We then demonstrate the capacity of this model to reproduce the characteristic logarithmic scaling of total crease length, thereby supplying a missing physical basis for the observed phenomenon.

42 ENGINEERING↗
DOE OSTI2021

The metabolic origins of non-photorespiratory CO 2 release during photosynthesis: a metabolic flux analysis

Respiration in the light (R L ) releases CO 2 in photosynthesizing leaves and is a phenomenon that occurs independently from photorespiration. Since R L lowers net carbon fixation, understanding R L could help improve plant carbon-use efficiency and models of crop photosynthesis. Although R L was identified more than 75 years ago, its biochemical mechanisms remain unclear. To identify reactions contributing to R L , we mapped metabolic fluxes in photosynthesizing source leaves of the oilseed crop and model plant camelina (Camelina sativa). We performed a flux analysis using isotopic labeling patterns of central metabolites during 13 CO 2 labeling time course, gas exchange, and carbohydrate production rate experiments. To quantify the contributions of multiple potential CO 2 sources with statistical and biological confidence, we increased the number of metabolites measured and reduced biological and technical heterogeneity by using single mature source leaves and quickly quenching metabolism by directly injecting liquid N 2 ; we then compared the goodness-of-fit between these data and data from models with alternative metabolic network structures and constraints. Our analysis predicted that R L releases 5.2 μmol CO2 g –1 FW h –1 of CO 2 , which is relatively consistent with a value of 9.3 μmol CO 2 g –1 FW h –1 measured by CO 2 gas exchange. The results indicated that ≤10% of R L results from TCA cycle reactions, which are widely considered to dominate R L . Further analysis of the results indicated that oxidation of glucose-6-phosphate to pentose phosphate via 6-phosphogluconate (the G6P/OPP shunt) can account for >% of CO 2 released by R L .

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Brittle failure analysis and modeling of high-burnup PWR fuel cladding alloys

The aim of this research is the development of methods for predicting mechanical behavior and identification of limiting conditions to prevent brittle failure of high-burnup (HBU) pressure water reactor (PWR) fuel cladding alloys. A finite element (FE) model of the ring compression test (RCT) was created to analyze the failure behavior of zirconium-based alloys with radial hydrides during the RCT. An elastic-plastic material model describes the zirconium alloy. The stress-strain curve needed for the elastic-plastic material model was derived by inverse finite element analyses. Cohesive zone modeling is used to reproduce sudden load drops during RCT loading. Based on the failure mechanism in non-irradiated ZIRLO (R) claddings, a micro-mechanical model was developed that distinguishes between brittle failure along hydrides and ductile failure of the zirconium matrix. Two different cohesive laws representing these types of failure are present in the same cohesive interface. The key differences between these constitutive laws are the cohesive strength, the stress at which damage initiates, and the cohesive energy, which is the damage energy dissipated by the cohesive zone. Statistically generated matrix-hydride distributions were mapped onto the cohesive elements and simulations with focus on the first load drop were performed. Computational results are in good agreement with the RCT results conducted on high-burnup M5 (R) samples. It could be shown that crack initiation and propagation strongly depend on the specific configuration of hydrides and matrix material in the fracture area.

Simbruner, Kai↗
NASA NTRS1975

The simulated depth history of dust grains in the lunar regolith

Trajectories giving the individual depth variations of lunar dust grains with time are randomly generated by a Monte Carlo code, where the variables are the mass and speed distribution of meteorites at the lunar surface and the geometrical shape of impact craters. A statistical analysis of a great number of such trajectories is then used to define the 'average' depth history of lunar dust grains for two grain radii: 1 and 50 microns. This yields: (1) estimates for time constants involved during the dynamic evolution of the regolith, (2) a model for the layering of the regolith, and (3) a better understanding of the basic dust-grain mechanisms responsible for the formation of the most mature lunar soil samples. The validity of various soil models proposed for the dynamic evolution of the regolith is discussed in terms of experimental constraints on the models.

Duraud, J. P.↗
NASA NTRS1990

The Venus nitric oxide night airglow - Model calculations based on the Venus Thermospheric General Circulation Model

The mechanism responsible for the Venus nitric oxide (0,1) delta band nightglow observed in the Pioneer Venus Orbiter UV spectrometer (OUVS) images was investigated using the Venus Thermospheric General Circulation Model (Dickinson et al., 1984), modified to include simple odd nitrogen chemistry. Results obtained for the solar maximum conditions indicate that the recently revised dark-disk average NO intensity at 198.0 nm, based on statistically averaged OUVS measurements, can be reproduced with minor modifications in chemical rate coefficients. The results imply a nightside hemispheric downward N flux of (2.5-3) x 10 to the 9th/sq cm sec, corresponding to the dayside net production of N atoms needed for transport.

Bougher, S. W.↗
NASA NTRS2007

Prospects for Pulsar Studies with the GLAST Large Area Telescope

The Large Area Telescope (LAT) on the Gamma-ray Large Area Space Telescope (GLAST), due to launch in November 2007, will have unprecedented sensitivity and energy resolution for gamma-rays in the range of 30 MeV to 200 GeV. GLAST is therefore expected to provide major advances in the understanding of high-energy emission from rotation-powered p ulsars. As the only presently known galactic GeV source class, pulsar s will be one of the most important sources for study with GLAST. The main science goals of the LAT for pulsar studies include an increase in the number of detected radio-loud and radio-quiet gamma-ray pulsar s, including millisecond pulsars, giving much better statistics for e lucidating population characteristics, measurement of the high-energy spectrum and the shape of spectral cutoffs and determining pulse profiles for a variety of pulsars of different age. Further, measurement of phase-resolved spectra and energy dependent pulse profiles of the brighter pulsars should allow detailed tests of magnetospheric partic le acceleration and radiation mechanisms, by comparing data with theo retical models that have been developed. Additionally, the LAT will have the sensitivity to allow blind pulsation searches of nearly all un identified EGRET sources, to possibly uncover more radio-quiet Geming a-like pulsars.

Harding, Alice K.↗