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 577 records · Page 32

Microwave remote sensing: Active and passive. Volume 2 - Radar remote sensing and surface scattering and emission theory

The fundamental principles of radar backscattering measurements are presented, including measurement statistics, Doppler and pulse discrimination techniques, and associated ambiguity functions. The operation of real and synthetic aperture sidelooking airborne radar systems is described, along with the internal and external calibration techniques employed in scattering measurements. Attention is given to the physical mechanisms responsible for the scattering emission behavior of homogeneous and inhomogeneous media, through a discussion of surface roughness, dielectric properties and inhomogeneity, and penetration depth. Simple semiempirical models are presented. Theoretical models involving greater mathematical sophistication are also given for extended ocean and bare soil surfaces, and the more general case of a vegetation canopy over a rough surface.

Ulaby, F. T.↗

Thermochemical Ablation Analysis of the Orion Heatshield

The Orion Multi-Purpose Crew Vehicle will one day carry astronauts to the Moon and beyond, and Orion's heatshield is a critical component in ensuring their safe return to Earth. The Orion heatshield is the structural component responsible for absorbing the intense heating environment caused by re-entry to Earth's atmosphere. The heatshield is primarily composed of Avcoat, an ablative material that is consumed during the re-entry process. Ablation is primarily characterized by two processes: pyrolysis and recession. The decomposition of in-depth virgin material is known as pyrolysis. Recession occurs when the exposed surface of the heatshield reacts with the surrounding flow. The Orion heatshield design was changed from an individually filled Avcoat honeycomb to a molded block Avcoat design. The molded block Avcoat heatshield relies on an adhesive bond to keep it attached to the capsule. In some locations on the heatshield, the integrity of the adhesive bond cannot be verified. For these locations, a mechanical retention device was proposed. Avcoat ablation was modelled in CHAR and the in-depth virgin material temperatures were used in a Thermal Desktop model of the mechanical retention device. The retention device was analyzed and shown to cause a large increase in the maximum bondline temperature. In order to study the impact of individual ablation modelling parameters on the heatshield sizing process, a Monte Carlo simulation of the sizing process was proposed. The simulation will give the sensitivity of the ablation model to each of its input parameters. As part of the Monte Carlo simulation, statistical uncertainties on material properties were required for Avcoat. Several properties were difficult to acquire uncertainties for: the pyrolysis gas enthalpy, non-dimensional mass loss rate (B´c), and Arrhenius equation parameters. Variability in the elemental composition of Avcoat was used as the basis for determining the statistical uncertainty in pyrolysis gas enthalpy and B´c. A MATLAB program was developed to allow for faster, more accurate and automated computation of Arrhenius reaction parameters. These parameters are required for a material model to be used in the CHAR ablation analysis program. This MATLAB program, along with thermogravimetric analysis (TGA) data, was used to generate uncertainties on the Arrhenius parameters for Avcoat. In addition, the TGA fitting program was developed to provide Arrhenius parameters for the ablation model of the gap filler material, RTV silicone.

Sixel, William↗

A probabilistic-phase field model for the fracture of brittle materials

Abstract We develop a computational method to determine the failure probability of brittle materials under general mechanical loading conditions. The method is a combination of two parts: (1) numerical simulations of materials with multiple cracks using phase field theory, where the complete fracture process is viewed as ‘damage percolation’ along critical paths or clusters of cracks, rather than the traditional weak-link failure mechanism of Weibull, and (2) an extension of the Batdorf statistical theory of fracture to finite domains, where it is implemented within the finite element framework. The results of phase-field simulations at the ‘percolation threshold’ are used as failure data in the Batdorf theory to determine the overall probability of failure. The input to this approach is the size distribution of cracks in a pristine material. An example is shown, where alumina samples that were previously tested by Abe and coworkers (Abeet al2003J. Am. Ceram. Soc.861019–21) in four-point loading are compared to the results of our numerical simulations. The approach developed here has the advantage of being extendable to more complex thermomechanical loading.

Materials Science↗

Mechanisms of Low-Level Jet Formation in the U.S. Mid-Atlantic Offshore

Abstract Low-level jets (LLJs), in which the wind speed attains a local maximum at low altitudes, have been found to occur in the U.S. mid-Atlantic offshore, a region of active wind energy deployment as of 2023. In contrast to widely studied regions such as the U.S. southern Great Plains and the California coastline, the mechanisms that underlie LLJs in the U.S. mid-Atlantic are poorly understood. This work analyzes floating lidar data from buoys deployed in the New York Bight to understand the characteristics and causes of coastal LLJs in the region. Application of the Hilbert–Huang transform, a frequency analysis technique, to LLJ case studies reveals that mid-Atlantic jets frequently occur during times of adjustment in synoptic-scale motions, such as large-scale temperature and pressure gradients or frontal passages, and that they do not coincide with motions at the native inertial oscillation frequency. Subsequent analysis with theoretical models of inertial oscillation and thermal winds further reveals that these jets can form in the stationary geostrophic wind profile from horizontal temperature gradients alone—in contrast to canonical LLJs, which arise from low-level inertial motions. Here, inertial oscillation can further modulate the intensity and altitude of the wind speed maximum. Statistical evidence indicates that these oscillations arise from stable stratification and the associated frictional decoupling due to warmer air flowing over a cold sea surface during the springtime land–sea breeze. These results improve our conceptual understanding of mid-Atlantic jets and may be used to better predict low-level wind speed maxima. Significance Statement The purpose of this work is to identify and characterize the atmospheric mechanisms that result in an occasional low-level maximum in the wind speed off the U.S. mid-Atlantic coastline. Our findings show that these low-level jets form due to horizontal temperature gradients arising from fronts and synoptic systems, as well as from the land–sea breeze that forces warmer air over the cold ocean surface. This work aids predictability of such jets, improves our understanding of this coastal environment, and has implications for future deployment of offshore wind energy in this region.

17 WIND ENERGY↗

A data driven approach for cross-slip modelling in continuum dislocation dynamics

Cross-slip is a thermally activated process by which screw dislocation changes its glide plane to another slip plane sharing the same Burgers vector. The rate at which this process happens is determined by a Boltzmann type expression that is a function of the screw segment length and the stress acting on the dislocation. In continuum dislocation dynamics (CDD), the information regarding the length of the screw dislocation segment and local stress state on dislocations are lost due to the coarse-grained representation of the density. Here, in this work, a data driven approach to characterize the lost information by analyzing the discrete dislocation configurations is proposed to enable cross-slip modeling in the CDD framework in terms of the coarse-grained dislocation density and stress fields. The analysis showed that the screw segment length follows an exponential distribution, and the stress fluctuations, defined as the difference between the stress on the dislocations and the mean field stress in CDD, follows a Lorentzian distribution. A novel approach for cross slip implementation in CDD employing the screw segment length and stress fluctuation statistics was proposed and rigorously tested by comparing the CDD cross-slip rates with discrete dislocation dynamics (DDD) rates. This approach has been applied in conjunction with three cross-slip models used in DDD simulations differing mainly in the functional form of cross slip activation energy. It was found that different cross-slip activation energy formulations yielded different cross-slip rates, yet the effect on mechanical stress-strain response and dislocation density evolution was minimal for the [001] type loading.

42 ENGINEERING↗

Possible Martian brines - Radar observations and models

The 1971 and 1973 Goldstone 12.6-cm radar observations of Mars are separate data sets which include reflectivity as a function of latitude, longitude, and season. It has been argued that secular reflectivity variations of Mars' surface are indicated by the data and that shallow subsurface melting is the causal mechanism most compatible with the observations; however, the melting hypothesis conflicts with accepted notions of the state and distribution of water on Mars. The data are examined to identify temporal and spatial domains within which statistically significant changes in measured reflectivity are clustered. Brines which might satisfy the best supported reflectivity variations are out of equilibrium with the chemical megaenvironment. It is unclear whether such a brine, if emplaced in the Martian regolith at a depth shallow enough to affect the radar reflectivity, could survive even a single freeze-thaw cycle. Some combination of unique scattering properties or some as yet unidentified process other than melting is responsible for any genuine reflectivity variations.

Zent, Aaron P.↗

Impact of Land Surface Initialization Approach on Subseasonal Forecast Skill: a Regional Analysis in the Southern Hemisphere

The authors use a sophisticated coupled land-atmosphere modeling system for a Southern Hemisphere subdomain centered over southeastern Australia to evaluate differences in simulation skill from two different land surface initialization approaches. The first approach uses equilibrated land surface states obtained from offline simulations of the land surface model, and the second uses land surface states obtained from reanalyses. The authors find that land surface initialization using prior offline simulations contribute to relative gains in subseasonal forecast skill. In particular, relative gains in forecast skill for temperature of 10%-20% within the first 30 days of the forecast can be attributed to the land surface initialization method using offline states. For precipitation there is no distinct preference for the land surface initialization method, with limited gains in forecast skill irrespective of the lead time. The authors evaluated the asymmetry between maximum and minimum temperatures and found that maximum temperatures had the largest gains in relative forecast skill, exceeding 20% in some regions. These results were statistically significant at the 98% confidence level at up to 60 days into the forecast period. For minimum temperature, using reanalyses to initialize the land surface contributed to relative gains in forecast skill, reaching 40% in parts of the domain that were statistically significant at the 98% confidence level. The contrasting impact of the land surface initialization method between maximum and minimum temperature was associated with different soil moisture coupling mechanisms. Therefore, land surface initialization from prior offline simulations does improve predictability for temperature, particularly maximum temperature, but with less obvious improvements for precipitation and minimum temperature over southeastern Australia.

Precipitation↗

Seasonal Water Transport in the Atmosphere of Mars: Applications of a Mars General Circulation Model Using Mars Global Surveyor Data

This is a Final Report for a Joint Research Interchange (JRI) between NASA Ames Research Center and San Jose State University, Department of Meteorology. We present below a summary of progress made during the duration of this JRI. The focus of this JRI has been to investigate seasonal water vapor transport in the atmosphere of Mars and its effects on the planet's present climate. To this end, the primary task has been to adapt a new dynamical processor for the adiabatic tendencies of the atmospheric circulation into the NASA Ames Mars general circulation model (MGCM). Using identical boundary and initial conditions, several comparative tests between the new and old MGCMs have been performed and the nature of the simulated circulations have been diagnosed. With confidence that the updated version of the Ames MGCM produces quite similar mean and eddy circulation statistics, the new climate model is well poised as a tool to pursue fundamental questions related to the spatial and seasonal variations of atmospheric water vapor on Mars, and to explore exchanges of water with non-atmospheric reservoirs and transport within its atmosphere. In particular, the role of surface sources and sinks can be explored, the range of water-vapor saturation altitudes can be investigated, and plausible precipitation mechanisms can be studied, for a range of atmospheric dust loadings, such future investigations can contribute to a comprehensive study of surface inventories, exchange mechanisms, and the relative importance of atmospheric transport Mars' water cycle. A listing of presentations made and manuscripts submitted during the course of this project is provided.

Hollingsworth, Jeffery L.↗

Seasonal Water Transport in the Atmosphere of Mars: Applications of a Mars General Circulation Model Using Mars Global Surveyor Data

This is a Final Report for a Joint Research Interchange (JRI) between NASA Ames Research Center and San Jose State University, Department of Meteorology. We present below a summary of progress made during the duration of this JRI. The focus of this JRI has been to investigate seasonal water vapor transport in the atmosphere of Mars and its effects on the planet's present climate. To this end, the primary task has been to adapt a new dynamical processor for the adiabatic tendencies of the atmospheric circulation into the NASA Ames Mars general circulation model (MGCM). Using identical boundary and initial conditions, several comparative tests between the new and old MGCMs have been performed and the nature of the simulated circulations have been diagnosed. With confidence that the updated version of the Ames MGCM produces quite similar mean and eddy circulation statistics, the new climate model is well poised as a tool to pursue fundamental questions related to the spatial and seasonal variations of atmospheric water vapor on Mars, and to explore exchanges of water with non-atmospheric reservoirs and transport within its atmosphere. In particular, the role of surface sources and sinks can be explored, the range of water-vapor saturation altitudes can be investigated, and plausible precipitation mechanisms can be studied, for a range of atmospheric dust loadings. Such future investigations can contribute to a comprehensive study of surface inventories, exchange mechanisms, and the relative importance of atmospheric transport Mars' water cycle. A listing of presentations made and manuscripts submitted during the course of this project is provided.

Hollingsworth, Jeffery L.↗

Equibiaxial flexural strength determination of UO 2 using a ball-on-ring test

To increase nuclear fuel performance and reliability, their mechanical properties require an accurate and statistically relevant assessment. Biaxial flexural strength tests provide an alternative to bend bar techniques for assessing mechanical behavior; namely, the transverse rupture strength (TRS) of ceramic samples. Biaxial test samples require simple geometries and minimal surface preparation, reducing fabrication costs and handling hazards. This study investigated the TRS of polycrystalline UO 2 fuel forms at room temperature using a ball-on-ring test fixture. Here, pellets were fabricated from UO 2 powder using conventional powder processing and sintering techniques. The TRS and Weibull parameters were obtained through Weibull statistics on over 60 UO 2 samples tested under equibiaxial flexure. The larger sample size in this study enabled a more robust Weibull statistical analysis than alternative test methods, which may not capture the stochastic failure of typical ceramics. Furthermore, two different loading ball diameters were employed to assess the impact of contact damage on fracture strength. While Hertzian contact damage was observed with the smaller loading ball, the fracture strength remained unaffected. A fracture analysis of the tested UO 2 samples indicated a mixture of intergranular and transgranular fracture that transitions to transgranular fracture with increasing distance from the fracture origin. The characteristic strength of the combined data sets was determined to be 148 MPa, and the Weibull modulus was determined to be 9.1. The TRS values and Weibull parameters were close to values found in the literature for alternative testing techniques using samples with similar microstructure and density. The findings in this study validate the ball-on-ring method used to obtain the TRS of UO 2 with a sample geometry more representative of nuclear fuels. Additionally, experimental TRS results from this study can be implemented in modeling codes to predict fuel performance, which is critical to fuel burnup extension and advanced nuclear fuel technologies.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Exotic Higgs Decays

Exotic decays of the Standard Model (SM)-like Higgs boson into beyond-the-SM particles are predicted in a wide range of well-motivated theories. The enormous samples of Higgs bosons that have been and will be produced at the Large Hadron Collider thus constitute one of the key discovery opportunities at that facility, particularly in the upcoming high-statistics, high-luminosity run. Here we review recent theoretical work on models that predict or accommodate exotic Higgs decays, survey the status of current experimental searches, and look forward to future capabilities at dedicated Higgs factories and beyond.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Single-Particle Studies Reveal a Nanoscale Mechanism for Elastic, Bright, and Repeatable ZnS:Mn Mechanoluminescence in a Low-Pressure Regime

Mechanoluminescent materials, which emit light in response to elastic deformation, are demanded for use as in situ stress sensors. ZnS doped with Mn is known to exhibit one of the lowest reported thresholds for appearance of mechanoluminescence, with repeatable light emission under contact pressure <10 MPa. Here, the physical basis for such behavior remains as yet unclear. Here, reliable microscopic detection of mechanoluminescence of single ZnS:Mn microparticles, in combination with nanoscale structural characterization, provides evidence that the mechanoluminescent properties of these particles result from interplay between a non-centrosymmetric crystal lattice and its defects, viz., dislocations and stacking faults. Statistical analysis of the distributions of mechanoluminescence energy release trajectories reveals two distinct mechanisms of excitation: one attributable to a piezo-phototronic effect and the other due to dislocation motion. At pressures below 8.1 MPa, both mechanisms contribute to mechanoluminescent output, with a dominant contribution from the piezo-phototronic mechanism. In contrast, above 8.1 MPa, dislocation motion is the primary excitation source. For the piezo-phototronic mechanism, we propose a specific model that accounts for elastic ZnS:Mn mechanoluminescence under very low pressure. The charged interfaces in stacking faults lead to the presence of filled traps, which otherwise would be empty in the absence of the built-in electric field. Upon application of external stress, local enhancement of the piezoelectric field at the stacking faults’ interfaces facilitates release of the trapped carriers and subsequent luminescence. This field enhancement explains how <10 MPa pressure produces thousands of photons.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Calculation of self-diffusion coefficients in supercritical carbon dioxide using mean force kinetic theory

This paper presents an application of mean force kinetic theory (MFT) to the calculation of the self-diffusivity of CO 2 in the supercritical fluid regime. Two modifications to the typical application of MFT are employed to allow its application to a system of molecular species. Therefore, the first is the assumption that the inter-particle potential of mean force can be obtained from the molecule center-of-mass pair correlation function, which in the case of CO 2 is the C–C pair correlation function. The second is a new definition of the Enskog factor that describes the effect of correlations at the surface of the collision volume. The new definition retains the physical picture that this quantity represents a local density increase, resulting from particle correlations, relative to that in the zero density homogeneous fluid limit. These calculations are facilitated by the calculation of pair correlation functions from molecular dynamics (MD) simulations using the FEPM2 molecular CO2 model. The self-diffusivity calculated from theory is in good agreement with that from MD simulations up to and slightly beyond the density at the location of the Frenkel line. The calculation is compared with and is found to perform similarly well to other commonly used models but has a greater potential for application to systems of mixed species and to systems of particles with long range interatomic potentials due to electrostatic interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Heat current fluctuations and anomalous transport in low-dimensional carbon lattices

Molecular dynamics simulations and methods of importance sampling are used to study the heat transport of low-dimensional carbon lattices. For both carbon nanotubes and graphene sheets, heat transport is found to be anomalous, violating Fourier’s law of conduction with a system size dependent thermal conductivity and concomitant nonlinear temperature profiles. For carbon nanotubes, the thermal conductivity is found to increase as the square root of the length of the nanotube, while for graphene sheets the thermal conductivity is found to increase as the logarithm of the length of the sheet over the system sizes considered. Furthermore, the particular length dependence and nonlinear temperature profiles place carbon lattices into a universality class with nonlinear lattice models, and suggest that heat transport through carbon nanostructures is better described by a Levy walk rather than simple diffusion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Microstructurally based model of fatigue initiation and growth

A model to calculate fatigue life is developed based on the assumption that fatigue life is entirely composed of crack growth from an initial microstructural inhomogeneity. Specifically, growth is considered to start from either an ellipsoidal void, a cracked particle, or a debonded particle. The capability of predicting fatigue life from material microstructure is based on linear elastic fracture mechanics principles, the sizes of the crack-initiating microstructural inhomogeneities, and an initiation parameter that is proportional to the cyclic plastic zone size. A key aspect of this modeling approach is that it is linked with a general purpose probability program to analyze the effect of the distribution of controlling microstructural features within the material. This enables prediction of fatigue stress versus life curves for various specimen geometries using distributional statistics obtained from characterizations of the microstructure. Results are compared to experimental fatigue data from an aluminum alloy.

Brockenbrough, J. R.↗

Mechanical Properties of Irradiated U-10wt%Mo Alloy Degraded by Porosity Development

A plate-type nuclear fuel consisting of a solid monolithic foil of U-10wt%Mo is under development for use in the United States’ high performance research reactors. In support of developing this fuel, the fuel has been fabricated for the first time by a commercial fuel vendor and subsequently irradiated in a test reactor. This provides an opportunity to evaluate post-irradiation mechanical properties of commercially fabricated fuel. Four-point bend testing was conducted on the irradiated U-10Mo fuel and the data produced includes bending strength and Young’s modulus. Although the material behaves in a brittle manner due to the developed porosity, a general trend of strength and modulus reduction are found as fission density increases. The data produced is evaluated using both Weibull statistics and a modulus degradation model with recommendations provided.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Probabilistic Fracture Mechanics Analysis of the Orbiter's LH2 Feedline Flowliner

Work performed by Southwest Research Institute (SwRI) as part of an Independent Technical Assessment (ITA) for the NASA Engineering and Safety Center (NESC) is summarized. The ITA goal was to establish a flight rationale in light of a history of fatigue cracking due to flow induced vibrations in the feedline flowliners that supply liquid hydrogen to the space shuttle main engines. Prior deterministic analyses using worst-case assumptions predicted failure in a single flight. The current work formulated statistical models for dynamic loading and cryogenic fatigue crack growth properties, instead of using worst-case assumptions. Weight function solutions for bivariant stressing were developed to determine accurate crack "driving-forces". Monte Carlo simulations showed that low flowliner probabilities of failure (POF = 0.001 to 0.0001) are achievable, provided pre-flight inspections for cracks are performed with adequate probability of detection (POD)-specifically, 20/75 mils with 50%/99% POD. Measurements to confirm assumed POD curves are recommended. Since the computed POFs are very sensitive to the cyclic loads/stresses and the analysis of strain gage data revealed inconsistencies with the previous assumption of a single dominant vibrant mode, further work to reconcile this difference is recommended. It is possible that the unaccounted vibrational modes in the flight spectra could increase the computed POFs.

Bonacuse, Peter J.↗