Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “STATISTICAL MECHANICS”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
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 OSTI2023

Reactor Network Analysis with Various Reaction Mechanisms to Investigate Hydrogen vs. Methane Fuel at Varying Flame Temperatures with Experimental Data

Emissions data were evaluated for a set of test hardware that adapts Collins Aerospace’s aeroengine liquid fueled injector technology to a ground-based turbine. The hardware was specifically developed to operate on 100% hydrogen; however it has the ability to operate on both pure methane and hydrogen and mixtures in between. 16 injector configurations and seven factors (air split, fuel and air swirl, pressure drop, preheat temperature, fuel composition, and flame temperature) were investigated based on a Box Behnken statistical model. Of the 16, configuration 3 was selected for further discussion being in the middle of the statistical design space. A chemical reactor network was developed based on the experimental data obtained and was assessed with three different mechanisms: GRI Mech 3.0, UCSD, and Galway. All factors were held constant except the fuel composition to directly compare chemical pathways for methane and hydrogen. This study was conducted across different flame temperatures of 1500K, 1675K, and 1850K for both hydrogen and methane. Additional lower flame temperatures of 1100K and 1300K were evaluated for hydrogen. Moreover, results for fully premixed and non-premixed results were compared to obtain information on the implications of mixedness on emissions. These results were compared with the Leonard and Stegmaier plot for methane and a similar plot was constructed for hydrogen.

08 HYDROGEN↗
DOE OSTI2022

Instrumentation development and testing of a wind turbine blade for sub-scale wake studies

Wind turbine blade aerodynamics and the resulting wake flow are complex, and wind tunnel testing of these flows can provide critical insight. The data from such tests are also valuable for validation of numerical models. For experiments using sub-scale turbines to be useful, the blade aerodynamics and wakes must exhibit the important features observed in utility-scale turbine flows. In this work, ~1 m blades designed to produce physics relevant to larger scale turbines were manufactured and affixed to an existing turbine. One of the blades was instrumented for blade surface pressure measurements. Time-dependent surface pressure measurements coupled with instantaneous measurements of the inflow with Cobra probes and the wake with particle image velocimetry allowed for characterization of the inflow, blade flow, and near wake. The results demonstrate that the instrumentation was effective in characterizing the blade loading and the flow field. For the one test case discussed in this paper, the measurements of inflow, blade loading, and wake properties facilitate understanding of the wake's behavior.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Advanced Laboratory and Field Arrays: Debris Modeling, Detection,& Mitigation (Task 1)

The statement of project objectives for this task was: develop tools, methods and models to assess, and mitigate the risk of damage to MHK infrastructure from woody debris. Develop the capability to detect woody debris using sonar and/or physical methods for purposes of characterizing debris statistics in river (at UAF’s Tanana River Test Site) and near-shore wave (at Yakutat, AK) environments and to activate debris mitigation measures. Develop debris impact risk maps and tables using statistics on debris size, geometry, type, prevalence, mobility and location. Improve and apply the COUPi discrete element method (DEM) to develop models of debris movement and impact on MHK infrastructure to evaluate risk of damage, and interference, to operations from debris. The proposed final deliverable for the task was a set of tools or techniques for providing estimates of the probability of debris impact, and resulting impact forces, on MHK infrastructure as a function of debris size, type, wave regime, and current velocity. Such estimates are required to assess damage risks to operational MHK infrastructure.

13 HYDRO ENERGY↗
NASA NTRS1973

Sandstorms and eolian erosion on Mars

It is now apparent that windblown sand plays a major role in the Martian environment. The existence of great sand and dust storms implies the lateral transport of large quantities of material on Mars that provide a convenient mechanism for the filling of craters- an event directly implied both by crater frequency-diameter statistics, and by the appearance of Martian craters. Such fast moving sand particles, saltating close to the Martian surface, can produce by collision and subsequent abrasion significant erosion of the Martian terrain.

Sagan, C.↗
NASA NTRS1980

Radio physics of the sun; Proceedings of the Symposium, University of Maryland, College Park, Md., August 7-10, 1979

Papers are presented in the areas of the radio characteristics of the quiet sun and active regions, the centimeter, meter and decameter wavelength characteristics of solar bursts, space observations of low-frequency bursts, theoretical interpretations of solar active regions and bursts, joint radio, visual and X-ray observations of active regions and bursts, and the similarities of stellar radio characteristics to solar radio phenomena. Specific topics include the centimeter and millimeter wave characteristics of the quiet sun, radio fluctuations arising upon the transit of shock waves through the transition region, microwave, EUV and X-ray observations of active region loops and filaments, interferometric observations of 35-GHz radio bursts, emission mechanisms for radio bursts, the spatial structure of microwave bursts, observations of type III bursts, the statistics of type I bursts, and the numerical simulation of type III bursts. Attention is also given to the theory of type IV decimeter bursts, Voyager observations of type II and III bursts at kilometric wavelengths, radio and whitelight observations of coronal transients, and the possibility of obtaining radio observations of current sheets on the sun.

Kundu, M. R.↗
NASA NTRS1981

A review of transhorizon propagation phenomena

Interference problems underlie the current interest in transhorizon propagation. In particular, statistics of the rare, high-level fields are of interest. This paper reviews the propagation mechanisms which produce the high-level fields and summarizes recent work in the modeling of the transhorizon propagation.

Crane, R. K.↗
NASA NTRS1990

Interplanetary magnetic field enhancements - Evidence for solar wind dust trail interactions

Magnetic disturbances related to the passage of the asteroid 2201 Oljato are described. It is pointed out that the process leading to these disturbances must be associated with material in the orbit of Oljato and not the Oljato asteroid itself. Using data from latest space observations, the statistical association of these disturbances with the orbit of 2201 Oljato is updated, and the possible mechanisms for causing the disturbances are discussed. The phenomenon is considered to be associated with comet-like processes in the solar wind.

Russell, C. T.↗
NASA NTRS1990

Scientific applications of frequency-stabilized laser technology in space

A synoptic investigation of the uses of frequency-stabilized lasers for scientific applications in space is presented. It begins by summarizing properties of lasers, characterizing their frequency stability, and describing limitations and techniques to achieve certain levels of frequency stability. Limits to precision set by laser frequency stability for various kinds of measurements are investigated and compared with other sources of error. These other sources include photon-counting statistics, scattered laser light, fluctuations in laser power, and intensity distribution across the beam, propagation effects, mechanical and thermal noise, and radiation pressure. Methods are explored to improve the sensitivity of laser-based interferometric and range-rate measurements. Several specific types of science experiments that rely on highly precise measurements made with lasers are analyzed, and anticipated errors and overall performance are discussed. Qualitative descriptions are given of a number of other possible science applications involving frequency-stabilized lasers and related laser technology in space. These applications will warrant more careful analysis as technology develops.

Schumaker, Bonny L.↗
NASA NTRS1993

The compositional structure and the effects of exothermicity in a nonpremixed planar jet flame

Results are presented of direct numerical simulation (DNS) of a randomly perturbed compressible, spatially developing two-dimensional planar jet under the influence of a finite rate chemical reaction of the type F + O yields product. The objectives of the simulations are to assess the compositional structure of the flame and to determine the influence of reaction exothermicity by means of statistical sampling of the DNS generated data. It is shown that even with this idealized kinetics model the simulated results exhibit features in accord with experimental data. These results indicate that the Damkohler number is an important parameter in determining the statistical composition of the reacting field and that the results are not very sensitive to the mechanism by which this parameter is varied. It is demonstrated that as the intensity of mixing is increased and the effect of finite rate chemistry is more pronounced, the magnitudes of the ensemble mean and variance of the product mass fraction decrease and those of the reactants' mass fraction increase. Also, at higher mixing rates the joint probability density functions of the reactants' mass fractions shift towards higher values within the composition domain, indicating a lower reactedness.

Steinberger, C. J.↗
NASA NTRS1992

Mechanical properties of silicate glasses exposed to a low-Earth orbit

The effects of a 5.8 year exposure to low earth orbit environment upon the mechanical properties of commercial optical fused silica, low iron soda-lime-silica, Pyrex 7740, Vycor 7913, BK-7, and the glass ceramic Zerodur were examined. Mechanical testing employed the ASTM-F-394 piston on 3-ball method in a liquid nitrogen environment. Samples were exposed on the Long Duration Exposure Facility (LDEF) in two locations. Impacts were observed on all specimens except Vycor. Weibull analysis as well as a standard statistical evaluation were conducted. The Weibull analysis revealed no differences between control samples and the two exposed samples. We thus concluded that radiation components of the Earth orbital environment did not degrade the mechanical strength of the samples examined within the limits of experimental error. The upper bound of strength degradation for meteorite impacted samples based upon statistical analysis and observation was 50 percent.

Wiedlocher, David E.↗
NASA NTRS2020

The Roughness Properties of Small Ice-Bearing Craters at the South Pole of the Moon: Implications for Accessing Fresh Water Ice in Future Surface Operations

The lunar poles provide a fascinating thermal environment capable of cold-trapping water ice on geologic timescales [1]. While there have been many observations indicating the presence of water ice at the lunar surface [e.g., 2–4], it is still not clear when this ice was delivered to the Moon. The timing of volatile dep-osition provides important constraints on the origin of lunar ice because different delivery mechanisms have been active at different times throughout lunar history. We previously found that some small (<10 km) cra-ters at the south pole of the Moon have morphologies suggestive of relatively young ages, on the basis of crisp crater rims [5]. These craters are too small to date with robust cratering statistics [5], but the possibility of ice in young craters is intriguing because it suggests that there is some recent and perhaps ongoing mechanism that is delivering or redistributing water to polar cold traps. Therefore, understanding if these small, ice-bear-ing craters are indeed young is essential in understand-ing the age and source of volatiles on the Moon. Here we take a new approach to understand the ages of these small polar cold traps: analyzing the roughness properties of small ice-bearing craters. It is well under-stood that impact crater properties (e.g., morphology, rock abundance, and roughness) evolve with time due to a variety of geologic and space-weathering processes [6–11]. Topographic roughness is a measurement of the local deviation from the mean topography, providing a measurement of surface texture, and is a powerful tool for evaluating surface evolution over geologic time [e.g., 11–14]. In this study we analyze the roughness of southern lunar craters (40°S–90°S) from all geologic eras, and determine how the roughness of small (<10 km) ice-bearing craters compare. We discuss the implications of the ages of ice-bearing south polar craters, and potential strategies for accessing fresh ice on the Moon.

Deutsch, Ariel N.↗
NASA NTRS1999

Statistical physics and physiology: monofractal and multifractal approaches

Even under healthy, basal conditions, physiologic systems show erratic fluctuations resembling those found in dynamical systems driven away from a single equilibrium state. Do such "nonequilibrium" fluctuations simply reflect the fact that physiologic systems are being constantly perturbed by external and intrinsic noise? Or, do these fluctuations actually, contain useful, "hidden" information about the underlying nonequilibrium control mechanisms? We report some recent attempts to understand the dynamics of complex physiologic fluctuations by adapting and extending concepts and methods developed very recently in statistical physics. Specifically, we focus on interbeat interval variability as an important quantity to help elucidate possibly non-homeostatic physiologic variability because (i) the heart rate is under direct neuroautonomic control, (ii) interbeat interval variability is readily measured by noninvasive means, and (iii) analysis of these heart rate dynamics may provide important practical diagnostic and prognostic information not obtainable with current approaches. The analytic tools we discuss may be used on a wider range of physiologic signals. We first review recent progress using two analysis methods--detrended fluctuation analysis and wavelets--sufficient for quantifying monofractual structures. We then describe recent work that quantifies multifractal features of interbeat interval series, and the discovery that the multifractal structure of healthy subjects is different than that of diseased subjects.

NASA Discipline Cardiopulmonary↗
NASA NTRS1993

On-board closed-loop congestion control for satellite based packet switching networks

NASA LeRC is currently investigating a satellite architecture that incorporates on-board packet switching capability. Because of the statistical nature of packet switching, arrival traffic may fluctuate and thus it is necessary to integrate congestion control mechanism as part of the on-board processing unit. This study focuses on the closed-loop reactive control. We investigate the impact of the long propagation delay on the performance and propose a scheme to overcome the problem. The scheme uses a global feedback signal to regulate the packet arrival rate of ground stations. In this scheme, the satellite continuously broadcasts the status of its output buffer and the ground stations respond by selectively discarding packets or by tagging the excessive packets as low-priority. The two schemes are evaluated by theoretical queuing analysis and simulation. The former is used to analyze the simplified model and to determine the basic trends and bounds, and the later is used to assess the performance of a more realistic system and to evaluate the effectiveness of more sophisticated control schemes. The results show that the long propagation delay makes the closed-loop congestion control less responsive. The broadcasted information can only be used to extract statistical information. The discarding scheme needs carefully-chosen status information and reduction function, and normally requires a significant amount of ground discarding to reduce the on-board packet loss probability. The tagging scheme is more effective since it tolerates more uncertainties and allows a larger margin of error in status information. It can protect the high-priority packets from excessive loss and fully utilize the downlink bandwidth at the same time.

Chu, Pong P.↗
DOE OSTI2022

Physics-Informed Gaussian Process Regression for States Estimation and Forecasting in Power Grids

Real-time state estimation and forecasting are critical for the efficient operation of power grids. In this paper, a physics-informed Gaussian process regression (PhI-GPR) method is presented and used for forecasting and estimating the phase angle, angular speed, and wind mechanical power of a three-generator power grid system using sparse measurements. In standard data-driven Gaussian process regression (GPR), parameterized models for the prior statistics are fit by maximizing the marginal likelihood of observed data. In the PhI-GPR method, we propose to compute the prior statistics offline by solving stochastic differential equations (SDEs) governing the power grid dynamics. The short-term forecast of a power grid system dominated by wind generation is complicated by the stochastic nature of the wind and the resulting uncertainty in wind mechanical power. Here, we assume that the power grid dynamics are governed by swing equations, with the wind mechanical power fluctuating randomly in time. We solve these equations for the mean and covariances of the power grid states using the Monte Carlo simulation method. We demonstrate that the proposed PhI-GPR method can accurately forecast and estimate observed and unobserved states. For the considered problem, PhI-GPR has computational advantages over the ensemble Kalman filter (EnKF) method: In PhI-GPR, ensembles are computed offline and independently of the data acquisition process, whereas for EnFK, ensembles are computed online with data acquisition, rendering real-time forecast more challenging. We also demonstrate that the PhI-GPR forecast is more accurate than the EnKF forecast when the random mechanical wind power is non-Markovian. In contrast, the two methods produce similar forecasts for the Markovian mechanical wind power. For observed states, we show that PhI-GPR provides a forecast comparable to the standard data-driven GPR; both forecasts are significantly more accurate than the autoregressive integrated moving average (ARIMA) forecast. We also show that the ARIMA forecast is more sensitive to observation frequency and measurement errors than the PhI-GPR forecast.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Thermo-mechanical characterization of shale using nanoindentation

Abstract Shale can be a potential buffer for high-level radioactive nuclear wastes. To be an effective buffer while subject to waste heat, shale's mechanical response at elevated temperature must be known. Many researchers have experimentally characterized the mechanical behavior of various shales at different length scales in adiabatic conditions. However, its mechanical performance at elevated temperatures at the nano-scale remains unknown. To investigate the temperature dependency of nanomechanical properties of shale, we conducted both experimental and numerical studies. In this study, we measured mechanical and fracture properties of shale, such as hardness, elastic modulus, anisotropy, and fracture toughness from 25 °C up to 300 °C at different bedding planes. Statistical analysis of the results suggests that hardness and fracture toughness significantly increased at temperatures from 100 to 300 °C; while, temperature does not have a significant impact on elastic modulus. Data also shows that the bedding plane orientations have a substantial impact on both mechanical and fracture properties of shale at the nano-scale leading to distinct anisotropic behavior at elevated temperature below 100 °C. Additionally, we numerically investigated the mechanical performance of the shale samples at room temperature to gain an insight into its mechanical response through the thickness. Numerical results were validated against the experimental results, confirming the simulation can be used to predict shale deformation at the nano-scale or potentially be used in multi-scale simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Comments on foliated gauge theories and dualities in 3+1d

We investigate the properties of foliated gauge fields and construct several foliated field theories in 3+1d that describe foliated fracton orders both with and without matter, including the recent hybrid fracton models. These field theories describe Abelian or non-Abelian gauge theories coupled to foliated gauge fields, and they fall into two classes of models that we call the electric models and the magnetic models. We show that these two classes of foliated field theories enjoy a duality. We also construct a model (using foliated gauge fields and an exactly solvable lattice Hamiltonian model) for a subsystem-symmetry protected topological (SSPT) phase, which is analogous to a one-form symmetry protected topological phase, with the subsystem symmetry acting on codimension-two subregions. We construct the corresponding gauged SSPT phase as a foliated two-form gauge theory. Some instances of the gauged SSPT phase are a variant of the X-cube model with the same ground state degeneracy and the same fusion, but different particle statistics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Benchmarking Microscale Ductility Measurements (Final Report of the Project DE-NE0008799)

Conventional macroscale experimentation is generally considered to be straightforward with few limitations. Conversely, micro/nanoscale experimentation presents numerous challenges in loading device design, sample preparation and handling, as well as accurate understanding of grain size and local texture effects on recorded measurements. Despite these challenges, nanopillar compression, MEMs based micro-tension, and nanoindentation approaches have been able to provide fundamental contributions to the understanding of material behavior at small lengthscales. However, the overarching shortcoming of these micro/nanoscale experimentation approaches, is the inability to directly translate measurements evaluated at the nm and µm length scales (e.g., hardness) to macroscale tensile material behavior (i.e., elastic modulus, yield strength, and ductility). The objectives of the proposed study are, 1) to establish best practices for obtaining tensile microscale ductility measurements, and 2) to validate methodologies to for comparing microscale ductility measurements to macroscale ductility measurements. In order to achieve these objectives, a multi-lengthscale, multi-temperature testing protocol and simulation framework are executed first on copper as a model material to validate the following approach, and second on reactor grade Zircaloy-2. Experiments are conducted on specimens extracted from the same test piece to ensure nominally identical grain size and texture from specimens to specimen. Motivated by the need to isolate the contribution of size-effects on obtained mechanical property measurements, specimens are manufactured with thicknesses at the micro- (1-10 µm), meso- (10-100's µm), and macroscales (sub-sized ASTM E8). In-situ full-field deformation techniques (scanning electron microscopy (SEM) grid methods and optical DIC) are incorporated into testing at each specimen length-scale to capture plasticity localization and evolution. Experimental testing for all specimens is conducted at both room temperature and elevated temperatures to probe the role of thermal activation on plastic deformation accommodation processes. Simulation efforts focus on examining the mechanical behavior of microscale specimens using a finite element approach with explicitly resolved grain morphologies, and an embedded crystal plasticity model. The cost-efficient implementation method allows for the modeling of a statistically significant number of both real (i.e., digital twin) and generated microstructures to obtain an understanding of the interrelationships between specimen microstructure and geometric variables (grain size, texture, specimen geometry, etc.) on microscale mechanical behavior.

36 MATERIALS SCIENCE↗