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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Effects of Normal Aging on Visuo-Motor Plasticity

Normal aging is associated with declines in neurologic function. Uncompensated visual and vestibular problems may have dire consequences including dangerous falls. Visuomotor plasticity is a form of behavioral neural plasticity which is important in the process of adapting to visual or vestibular alteration, including those changes due to pathology, pharmacotherapy, surgery or even entry into a microgravity or underwater environment. In order to determine the effects of aging on visuomotor plasticity, we chose the simple and easily measured paradigm of visual-motor re-arrangement created by using visual displacement prisms while throwing small balls at a target. Subjects threw balls before, during and after wearing a set of prisms which displace the visual scene by twenty degrees to the right. Data obtained during adaptation were modeled using multilevel analyses for 73 subjects aged 20 to 80 years. We found no statistically significant difference in measures of visuomotor plasticity with advancing age. Further studies are underway examining variable practice training as a potential mechanism for enhancing this form of behavioral neural plasticity.

Roller, Carrie A.↗
DOE OSTI2022

Time-dependent THMC properties and microstructural evolution of damaged rocks in excavation damage zone

Modeling coupled thermo-hydro-mechanical-chemical (THMC) processes in host rocks near high-level nuclear waste (HLW) repositories at various time scales is an extremely challenging task. The current study integrates experimental, theoretical, and numerical methods in assessing the evolution of excavation damage zone (EDZ) over time and its implication on the long-term migration of hazardous species. Argillite and rock salt and are the focus of this study. The first part of the report presents a novel time-dependent directional microcrack damage theory for generic brittle rocks. It features detailed statistical description of the microcracks within a damaged solid, permitting a direct upscaling of microscale processes such as crack growth kinetics, crack closure/opening, sliding friction to explain the macroscopic creep, nonlinear elasticity, shear dilation, and loading-unloading hysteresis. This provides a basic platform for describing the anisotropic mechanical and transport properties of damaged rocks during excavation and subsequent THMC loadings. The model is validated and numerically implemented to Finite Element (FE) package ABAQUS through the user-defined material (UMAT) interface and have demonstrated great potential in resolving the time-dependent and anisotropic evolution of damage in the EDZ. The second part of the report presents a multi-scale experimental effort in characterizing the thermal, hydraulic, and mechanical properties of Mancos shale and Avery Island salt. For the Mancos shale, triaxial compression tests are performed at different confining pressures and temperatures to probe its thermomechanical properties relevant to HLW repositories. The obtained stress-strain data are interpreted using the proposed directional damage theory. Post-test specimens are subjected to gas permeability tests to reveal the correlation between permeability and the degree of microcracking. At microscale, temperature-controlled nanoindentation tests are performed and found a linear correlation between fracture toughness and elastic modulus from 25°C to 300°C. For Avery Island salt, we have designed and manufactured a novel relative-humidity controlled uniaxial creep device. Long-term creep tests at low stresses (< 5 MPa) are performed at different levels of relative humidity (RH). Besides confirming the much higher creep rates as one would expect through extrapolating the high-stress creep data, the results reveal that the steady-state creep rate of rock salt is strongly dependent on the ambient RH, an aspect that is often neglected in the literature. Both behaviors can be attributed to the pressure-solution creep mechanism which dominates at low stress and high RH levels. The third part of the report explores a set of numerical strategies in modeling the THMC behavior of porous geomaterials. A fully implicit, monolithic FE solution that can flexibly interface with different material models and coupling mechanisms for THMC problems is developed and verified through the ABAQUS user-defined element (UEL) interface. The scheme is used to study the THM response of a hypothetical HLW storage site with reference to an existing in-situ heater test. Strategies for integrating the UEL and the microcrack UMAT are suggested. Another numerical endeavor of this study is to implement a higher-order asymptotic homogenization method to account for the heterogeneous porous structures. The same method is then extended to perform microstructure-informed thermo-mechanical modeling of generalized continua. The above outcomes of this project provide a strong thrust towards enhancing the fundamental understanding and modeling capability of the long-term evolution of host rocks in EDZ, thus helping achieve the design goal of 1-million-year isolation of high-level nuclear wastes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

A coupled vacancy diffusion-dislocation dynamics model for the climb-glide motion of jogged screw dislocations

Here, we develop a novel model to study the climb/glide motion of jogged screw dislocations within the discrete dislocation dynamics (DDD) framework. We present results for the dependence of the climb velocity on the applied stress and on the jog size and distribution statistics. We show that the model predictions are consistent with experimental data in both γ-TiAl and Zircaloy-4. The ranges of the applied stress and jog spacing that determine the dominance of one of three dislocation mechanisms are identified. These are the jog dragging, dipole dragging, and dipole bypass mechanisms, respectively. The overall dislocation motion in the jog dragging regime is composed of glide of screw segments and climb of jogs, controlling the plastic strain and the creep rate, respectively. Based on current simulations and on a detailed examination of the predicted jog heights compared to experiments, we advance the hypothesis that a combination of jog dragging and dipole bypass mechanisms is necessary to reproduce the high creep rate observed in some experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Towards Quantum Computing Phase Diagrams of Gauge Theories with Thermal Pure Quantum States

The phase diagram of strong interactions in nature at finite temperature and chemical potential remains largely theoretically unexplored due to inadequacy of Monte-Carlo–based computational techniques in overcoming a sign problem. Quantum computing offers a sign-problem-free approach, but evaluating thermal expectation values is generally resource intensive on quantum computers. To facilitate thermodynamic studies of gauge theories, we propose a generalization of the thermal-pure-quantum-state formulation of statistical mechanics applied to constrained gauge-theory dynamics, and numerically demonstrate that the phase diagram of a simple low-dimensional gauge theory is robustly determined using this approach, including mapping a chiral phase transition in the model at finite temperature and chemical potential. Quantum algorithms, resource requirements, and algorithmic and hardware error analysis are further discussed to motivate future implementations. Thermal pure quantum states, therefore, may present a suitable candidate for efficient thermal simulations of gauge theories in the era of quantum computing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1980

Ion chemistry in interstellar space

The temperature dependence of the bimolecular reactions of Cl(+), HCl(+), CH(+), CH2(+), N(+), NH(+), and NH2(+) reacting with H2 have been investigated. For Cl(+) and HCl(+) rate constants have been determined over the temperature range of 150 to 400 K. Preliminary data have been obtained for the other systems. The Cl(+)/H2 system shows a weak, but significant positive temperature dependence that could be important in interstellar modeling studies. The HCl(+)/H2 system shows a substantial negative temperature dependence with the rate constant approaching the collision rate at low temperatures. The association reaction CH3(+) + HCN yields CH3 HCN(+) has been theoretically modeled using statistical phase space theory. Both radiative and collisional stabilization have been included. The results are compared with experiment with good agreement obtained over wide variations in T and p. This system is potentially important in the mechanism of large molecule synthesis in interstellar space.

Bowers, M. T.↗
DOE OSTI2020

Mapping the kinetic evolution of metastable grain boundaries under non-equilibrium processing

The kinetic evolution of a multiplicity of metastable grain boundaries (GBs) under fast driving conditions are studied by atomistic modeling. Assisted with an enhanced statistical analysis, the energetic evolution of GBs over a broad metastability-temperature space is mapped out, wherein two distinct regimes—an ageing regime and a rejuvenating regime—are retrieved with high fidelity. By comparing the results under various conditions (e.g. random perturbations, isothermal annealing, and fast heating/cooling), it is shown that such ageing/rejuvenating mechanism map is universal, irrespective of the actual stimuli used to elicit the metastable GBs. Here, the ageing/rejuvenating phenomena are demonstrated to stem from the energy imbalance of uphill climbing and downhill dropping during sequential transitions in the system's potential energy landscape. Without the necessity of introducing free parameters, such model can reconcile experimentally measured hardness variation of nanocrystalline metals subjected to femto-second laser irradiation, and it therefore provides a novel perspective on achieving a plurality of interfacial states and facilitating previously inaccessible property regimes.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Direct Monte Carlo Simulations of Chemical Equilibrium Composition of Molecular Fluid Mixtures Under Shock Conditions

A new Monte Carlo simulation method is implemented which holds the number and kind of atoms constant rather than the number and kind of molecules for a molecular fluid mixture. Chemical reactions are allowed as a correlated move of atoms such as an interchange of atoms between molecules. Equilibrium chemical composition is determined as an average over the simulation along with equation of state quantities. Results are presented for N 2 + O 2 → 2NO at 30GPa and 3000K. The equilibrium composition is shown to be sensitive to the potentials between unlike species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2002

Acoustic Analogy and Alternative Theories for Jet Noise Prediction

Several methods for the prediction of jet noise are described. All but one of the noise prediction schemes are based on Lighthill's or Lilley's acoustic analogy, whereas the other is the jet noise generation model recently proposed by Tam and Auriault. In all of the approaches, some assumptions must be made concerning the statistical properties of the turbulent sources. In each case the characteristic scales of the turbulence are obtained from a solution of the Reynolds-averaged Navier-Stokes equation using a kappa-sigma turbulence model. It is shown that, for the same level of empiricism, Tam and Auriault's model yields better agreement with experimental noise measurements than the acoustic analogy. It is then shown that this result is not because of some fundamental flaw in the acoustic analogy approach, but instead is associated with the assumptions made in the approximation of the turbulent source statistics. If consistent assumptions are made, both the acoustic analogy and Tam and Auriault's model yield identical noise predictions. In conclusion, a proposal is presented for an acoustic analogy that provides a clearer identification of the equivalent source mechanisms, as is a discussion of noise prediction issues that remain to be resolved.

Morris, Philip J.↗
NASA NTRS2002

The Acoustic Analogy and Alternative Theories for Jet Noise Prediction

This paper describes several methods for the prediction of jet noise. All but one of the noise prediction schemes are based on Lighthill's or Lilley's acoustic analogy while the other is the jet noise generation model recently proposed by Tam and Auriault. In all the approaches some assumptions must be made concerning the statistical properties of the turbulent sources. In each case the characteristic scales of the turbulence are obtained from a solution of the Reynolds-averaged Navier Stokes equation using a k - epsilon turbulence model. It is shown that, for the same level of empiricism, Tam and Auriault's model yields better agreement with experimental noise measurements than the acoustic analogy. It is then shown that this result is not because of some fundamental flaw in the acoustic analogy approach: but, is associated with the assumptions made in the approximation of the turbulent source statistics. If consistent assumptions are made, both the acoustic analogy and Tam and Auriault's model yield identical noise predictions. The paper concludes with a proposal for an acoustic analogy that provides a clearer identification of the equivalent source mechanisms and a discussion of noise prediction issues that remain to be resolved.

Morris, Philip J.↗
NASA NTRS2002

The Acoustic Analogy and Alternative Theories for Jet Noise Prediction

This paper describes several methods for the prediction of jet noise. All but one of the noise prediction schemes are based on Lighthill's or Lilley's acoustic analogy while the other is the jet noise generation model recently proposed by Tam and Auriault. In all the approaches some assumptions must be made concerning the statistical properties of the turbulent sources. In each case the characteristic scales of the turbulence are obtained from a solution of the Reynolds-averaged Navier Stokes equation using a k-epsilon turbulence model. It is shown that, for the same level of empiricism, Tam and Auriault's model yields better agreement with experimental noise measurements than the acoustic analogy. It is then shown that this result is not because of some fundamental flaw in the acoustic analogy approach: but, is associated with the assumptions made in the approximation of the turbulent source statistics. If consistent assumptions are made, both the acoustic analogy and Tam and Auriault's model yield identical noise predictions. The paper concludes with a proposal for an acoustic analogy that provides a clearer identification of the equivalent source mechanisms and a discussion of noise prediction issues that remain to be resolved.

Morris, Philip J.↗
DOE OSTI2025

Fuel performance analysis of fully-resolved TRISO compact

The TRi-structural ISOtropic (TRISO) fuel multilayered coating structure offers multiple barriers to fission product release, enhancing safety and performance. The heterogeneous nature of TRISO fuel compacts, comprising thousands of randomly distributed coated fuel particles embedded in a graphite matrix, creates intricate stress fields and thermal gradients that cannot be accurately modeled using simplified one-dimensional or homogenized approaches. Consequently, three-dimensional modeling enables the prediction of fuel compact dimensional changes, internal pressure buildup, and fission product transport pathways under diverse irradiation and thermal conditions. This capability facilitates detailed analysis of particle-to-particle interactions, matrix cracking mechanisms, and the statistical distribution of coating failures, which directly impact fuel performance and safety margins. This capability is particularly critical for advanced reactors, such as high-temperature gas-cooled reactors and other Generation IV reactor designs where TRISO fuel operates at elevated temperatures and burn-up levels. This work introduces a novel method to generate an optimized packing of TRISO compacts and a complete 3D mesh with random distribution of TRISO particles, which are discretized into each coating component layer.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1983

Wildland inventory and resource modeling for Douglas and Carson City Counties, Nevada, using LANDSAT and digital terrain data

The potential of using LANDSAT satellite imagery to map and inventory pinyon-juniper desert forest types in Douglas and Carson City Counties, Nevada was demonstrated. Specific map and statistical products produced include land cover, mechanical operations capability, big game winter range habitat, fire hazard, and forest harvestability. The Nevada Division of Forestry determined that LANDSAT can produce a reliable and low-cost resource data. Added benefits become apparent when the data are linked to a geographical information system (GIS) containing existing ownership, planning, elevation, slope, and aspect information.

Brass, J. A.↗
DOE OSTI2021

Chemical Understanding of Actinide Separations

This work provides a broad overview of the separation processes used to isolate actinides and the experimentally and computationally determined chemical characteristics that define those separations. The redox chemistry of the actinides plays a pivotal role in both aqueous and pyrochemical processing separations. The near-overlapping energies of the 6d and 5f orbitals in the light actinides allow for facile adjustment of actinide oxidation states, which is used in many established separation methods. In contrast, the stable, generally 3+ oxidation states of the mid- and heavy actinides can make it difficult to separate them from the similarly lanthanides(III). In aqueous separations, the tendency of the actinides to form anionic and neutral aqueous complexes with a variety of complexants (especially soft donors) is used to achieve high separation factors between chemically similar elements in both solid–liquid separations and liquid–liquid extraction. This selectivity can be further tuned through the use of specialized organic or solid-phase ligands. Pyroprocessing separations utilize the unique redox behavior of the actinides to adjust their distribution between a molten salt electrolyte and either a solid electrode or molten metal phase. Atomic-level insights into the mechanisms underlying actinide separation processes, with the ultimate goal of predicting separation behavior, can be provided by electronic structure and statistical mechanical-based calculation methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2000

The Turbulent/Non-Turbulent Interface Bounding a Far-Wake

The velocity fields of a turbulent wake behind a flat plate obtained from the direct numerical simulations of Moser et al. are used to study the structure of the flow in the intermittent zone where there are, alternately, regions of fully turbulent flow and non-turbulent velocity fluctuations either side of a thin randomly moving interface. Comparisons are made with a wake that is 'forced' by amplifying initial velocity fluctuations. There is also a random temperature field T in the flow; T varies between constant values of 0.0 and 1.0 on the sides of the wake. The value of the Reynolds number based on the centreplane mean velocity defect and halfwidth b of the wake is Re approx. = 2000. It is found that the thickness of the continuous interface is about equal to 0.07b, whereas the amplitude of fluctuations of the instantaneous interface displacement y(sub I)(t) is an order of magnitude larger, being about 0.5b. This explains why the mean statistics of vorticity in the intermittent zone can be calculated in terms of the probability distribution of y(sub I) and the instantaneous discontinuity in vorticity across the interface. When plotted as functions of y - y(sub I), the conditional mean velocity (U) and temperature (T) profiles show sharp jumps Delta(U) and Delta(T) at the interface adjacent to a thick zone where (U) and (T) vary much more slowly. Statistics for the vorticity and velocity variances, available in such detail only from DNS data, show how streamwise and spanwise components of vorticity are generated by vortex stretching in the bulges of the interface. Flow fields around the interface, analyzed in terms of the local streamline pattern, confirm previous results that the advancement of the vortical interface into the irrotational flow is driven by large-scale eddy motion. It is argued that because this is an inviscid mechanism the entrainment process is not sensitive to the value of Re, and that small-scale nibbling only plays a subsidiary role. While mean Reynolds stresses decrease gradually in the intermittent zone, conditional stresses are found to decrease sharply towards zero at the interface. Using one-point turbulence models applied to either unconditional or conditional statistics for the turbulent region and then averaged, the entrainment rate E(sub b) would, if calculated exactly, be zero. But if computed with standard computational methods, E(sub b) would be non-zero because of numerical diffusion. It is concluded that the current practice in statistical models of approximating entrainment by a diffusion process is computationally arbitrary and physically incorrect. An analysis shows how E(sub b) is related to Delta(U) and the jump in shear stress at the interface, and correspondingly to Delta(T) and the heat flux.

Bisset, David K.↗
NASA NTRS1994

A time-dependent gyro-kinetic model of thermal ion upflows in the high-latitude F region

Ample evidence supports the significance of the high-latitude ionospheric contribution to magnetospheric plasma. Assuming flux conservation along a flux tube, the upward field-aligned ion flows observed in the magnetosphere require high-latitude ionospheric field-aligned ion upflows of the order of 10(exp 8) to 10(exp 9)/sq cm/s. Since radar and satellite observations of high-latitude F region flows at times exceed this flux requirement by an order of magnitude, the thermal ionospheric upflows are not simply the ionospheric response to a magnetospheric flux requirement. Several ionospheric ion upflow mechanisms have been proposed, but simulations based on fluid theory do not reproduce all the observed features of ionospheric ion upflows. Certain asymmetries in the statistical morphology of high-latitude F region ion upflows suggest that the ion upflows may be generated by ion-neutral frictional heating. We developed a single-component (O(+)), time-dependent gyro-kinetic model of the high-latitude F region response to frictional heating in which the neutral exobase is a discontinuous boundary between fully collisional and collisionless plasmas. The concept of a discontinuous neutreal exobase and the assumption of a constant and uniform polarization electric field reduce the ion velocity distribution function, from which we can compute the ion density, parallel velocity, parallel and perpendicular temperature, and parallel flux. Using our model, we simulated the response of a convecting flux tube between 500 km and 2500 km to various frictional heating inputs; the results were both qualitatively and quantitatively different from fluid model results, which may indicate an inadequacy of the fluid theory approach. The gyro-kinetic frictional heating model responses to the various simulations were qualitatively similar: (1) initial perturbations of all the modeled parameters propagated rapidly up the flux tube, (2) transient values of the ion parallel velocity, temperature, and flux exceeded 3 km/s, 2 x 10(exp 4) K, and 10(exp 9)/sq cm/s, respectively, (3) a second transient regime developed wherein the parallel temperature drops to very low values (a few hundred Kelvins), and (4) well after heating ceased, large parallel temperatures and large downward parallel velocities and fluxes developed as the flux tube slowly returned to diffusive equilibrium. The ion velocity distributions during the simulation are often non-Maxwellian and are sometimes composed of two distinct ion populations.

Loranc, M.↗
DOE OSTI2020

Quantifying the effect of tow architecture variability on the performance of biaxially braided composite tubes

Targeted experiments are performed with stereoscopic digital image correlation to quantify the effect of irregular braid structure on the local surface mechanical response of carbon-fiber reinforced epoxy composite tubes. Virtual models comparing both nonuniform and uniform tow structure of the outer braid of each specimen are integrated into numerical simulations similar to experiments. Experimental and numerical results agree well. Statistical tests are used to support correlations between as-manufactured braid variation ranges to local increases in stress up to 18% compared to an ideal, uniform structure in the elastic regime. Results support improved manufacturing and modeling efforts of defect sensitive braided composites.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Learning turbulent flows with generative models for super resolution and sparse flow reconstruction

Neural operators are promising surrogates for dynamical systems but when trained with standard L 2 losses they tend to oversmooth fine-scale turbulent structures. Here, we show that combining operator learning with generative modeling overcomes this limitation. We consider three practical turbulent-flow challenges where conventional neural operators fail: spatio-temporal super-resolution, forecasting, and sparse flow reconstruction. For Schlieren jet super-resolution, an adversarially trained neural operator (adv-NO) reduces the energy-spectrum error by 15 × while preserving sharp gradients at neural operator-like inference cost. For 3D homogeneous isotropic turbulence, adv-NO trained on only 160 timesteps from a single trajectory forecasts accurately for five eddy-turnover times and offers 114 × wall-clock speed-up at inference than the baseline diffusion-based forecasters, enabling near-real-time rollouts. For reconstructing cylinder wake flows from highly sparse Particle Tracking Velocimetry-like inputs, a conditional generative model infers full 3D velocity and pressure fields with correct phase alignment and statistics. These advances enable accurate reconstruction and forecasting at low compute cost, bringing near-real-time analysis and control within reach in experimental and computational fluid mechanics.

Fluid dynamics↗