SEARCH · Engineering Papers

Results for “computational 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

DOE OSTI2024

On the Statistical Mechanics of Mass Accommodation at Liquid–Vapor Interfaces

Here we propose a framework for describing the dynamics associated with the adsorption of small molecules to liquid-vapor interfaces using an intermediate resolution between traditional continuum theories that are bereft of molecular detail and molecular dynamics simulations that are replete with them. In particular, we develop an effective single particle equation of motion capable of describing the physical processes that determine thermal and mass accommodation probabilities. The effective equation is parametrized with quantities that vary through space away from the liquid-vapor interface. Of particular importance in describing the early time dynamics is the spatially dependent friction, for which we propose a numerical scheme to evaluate from molecular simulation. Taken together with potentials of mean force computable with importance sampling methods, we illustrate how to compute the mass accommodation coefficient and residence time distribution. Throughout, we highlight the case of ozone adsorption in aqueous solutions and its dependence on electrolyte composition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1988

Monolithic ceramic analysis using the SCARE program

The Structural Ceramics Analysis and Reliability Evaluation (SCARE) computer program calculates the fast fracture reliability of monolithic ceramic components. The code is a post-processor to the MSC/NASTRAN general purpose finite element program. The SCARE program automatically accepts the MSC/NASTRAN output necessary to compute reliability. This includes element stresses, temperatures, volumes, and areas. The SCARE program computes two-parameter Weibull strength distributions from input fracture data for both volume and surface flaws. The distributions can then be used to calculate the reliability of geometrically complex components subjected to multiaxial stress states. Several fracture criteria and flaw types are available for selection by the user, including out-of-plane crack extension theories. The theoretical basis for the reliability calculations was proposed by Batdorf. These models combine linear elastic fracture mechanics (LEFM) with Weibull statistics to provide a mechanistic failure criterion. Other fracture theories included in SCARE are the normal stress averaging technique and the principle of independent action. The objective of this presentation is to summarize these theories, including their limitations and advantages, and to provide a general description of the SCARE program, along with example problems.

Manderscheid, Jane M.↗

DOE OSTI2020

Pass-efficient methods for compression of high-dimensional turbulent flow data

The future of high-performance computing, specifically on future Exascale computers, will presumably see memory capacity and bandwidth fail to keep pace with data generated, for instance, from massively parallel partial differential equation (PDE) systems. Current strategies proposed to address this bottleneck entail the omission of large fractions of data, as well as the incorporation of in situ compression algorithms to avoid overuse of memory. To ensure that post-processing operations are successful, this must be done in a way that a sufficiently accurate representation of the solution is stored. Moreover, in situations where the input/output system becomes a bottleneck in analysis, visualization, etc., or the execution of the PDE solver is expensive, the number of passes made over the data must be minimized. In the interest of addressing this problem, this work focuses on the utility of pass-efficient, parallelizable, low-rank, matrix decomposition methods in compressing high-dimensional simulation data from turbulent flows. Additionally, a particular emphasis is placed on using coarse representation of the data – compatible with the PDE discretization grid – to accelerate the construction of the low-rank factorization. This includes the presentation of a novel single-pass matrix decomposition algorithm for computing the so-called interpolative decomposition. The methods are described extensively and numerical experiments on two turbulent channel flow data are performed. In the first (unladen) channel flow case, compression factors exceeding 400 are achieved while maintaining accuracy with respect to first- and second-order flow statistics. In the particle-laden case, compression factors of 100 are achieved and the compressed data is used to recover particle velocities. These results show that these compression methods can enable efficient computation of various quantities of interest in both the carrier and disperse phases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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.↗

DOE OSTI2024

Functional Data Analysis in Wearable Body Sensor Networks

Improving response time of indirect room-size calorimeters is still an outstanding problem in metabolic research. Accurate estimates of instantaneous rates of gaseous exchange require numerical differentiation of measured gaseousgas concentrations. We propose a new method to estimate the instantaneous gaseousgas exchange rates in indirect calorimetry. In contrast to the previously developed techniques, the method addresses the problem of differentiation of gaseous concentrations as an ill-posed problem. By applying the method of regularization, the problem of differentiation is converted into a well-posed problem resulting in smooth and consistent gaseous exchange rates. The validity of the method is tested on a large dataset of calorimeter experiments which included 313 human experiments along with 231 alcohol combustion experiments. It is demonstrated that the method is able to reliably differentiate between the “unphysiological” process of alcohol combustion and physiological variations produced by human metabolism. The method also allowed unraveling the previously unreported relative kinetics of O2 consumption and Respiratory Quotient (RQ) in humans. It was found that the kinetics of oxidative fuel selection lags behind the energy expenditure in humans exhibiting some sort of oxidative inertia. The time lag varies from 2-3 min up to 30 min, depending on particular individual. No such lag was found in alcohol combustion experiments. In addition to the relative kinetics of substrate oxidation, two statistical indexes reflecting variability of minute-by-minute RQ were estimated. The indexes were the RQ’s standard deviation and RQ’s first-order derivative. Both indexes showed statistically significant difference between human experiments and alcohol combustion experiments. We conclude that the proposed method can consistently extract physiologically-relevant information from noisy calorimetry data and the aforesaid information can provide additional insights into the mechanism of metabolic fuel selection in humans.

60 - APPLIED LIFE SCIENCES↗

DOE OSTI2026

Stress and Strain Heterogeneity and Persistence in Uniaxially‐ and Triaxially‐Loaded Sandstone

Two critical questions in brittle rock mechanics are how rocks developed localized strains and to what extent internal stress heterogeneity controls this localization and subsequent macroscopic failure. Definitive answers have not yet been found, but would provide insight into rock fracture mechanics as relevant to hydrocarbon extraction and sequestration. Here, we use synchrotron X‐ray tomography (XRT) and 3D X‐ray diffraction (3DXRD) during uniaxial and triaxial tests on Nugget and Bentheimer sandstones to examine strain and stress localization prior to mechanical failure. 3DXRD was used to measure intra‐granular lattice strains which were used to compute elastic stress tensors of each grain. Digital volume correlation (DVC) was applied to XRT images to determine the strain field in the sample. Both samples featured marked spatial heterogeneity, localization, and temporal persistence of elevated stresses and strains during their mechanical deformation toward failure. Both samples featured a majority of grains with at least one principal stress component that was tensile, a signature of the influence of heterogeneity on stress transmission. Measurements further revealed that compressive stress orientations and statistics evolved in a similar manner to those of inter‐particle forces in loose granular materials, with triaxially‐compressed rock exhibiting enhanced grain stress heterogeneity compared to uniaxially‐compressed rock. Our results complement recent work by others who employed XRT and scanning 3DXRD to study triaxially‐compressed sandstone, but extend those results to uniaxial compression, sandstones of varied porosity, and grain stress measurements throughout the 3D full extent of the samples rather than in a single layer examined with scanning 3DXRD.

58 GEOSCIENCES↗

DOE OSTI2023

The advanced characterization, post-irradiation examination, and materials informatics for the development of ultra high-burnup annular U-10Zr metallic fuel

U-Zr metallic fuel is a promising fuel candidate for Gen Ⅳ fast spectrum reactors. Previous experimental irradiation campaigns showed that the sodium thermal bonded U-10Zr fuel design can achieve a burnup of 10% fissions per initial heavy metal atom (FIMA). Advanced metallic fuel designs are pushing the burnup limit to 20% or even 30% FIMA. To achieve the higher burnup and eliminate the pyrophoric sodium, a prototypical annular fuel has been designed, fabricated, clad with HT-9 in the Materials and Fuels Complex, and irradiated in the Advanced Test Reactors of Idaho National Laboratory (INL) to a peak burnup of 3.3% FIMA. During irradiation, the mechanical contact between fuel and cladding acts as a thermal bond. The irradiation lasted for 132 days in the reactor. Recently, the archived fresh and irradiated fuel samples were characterized using advanced characterization capabilities in the Irradiated Materials Characterization Laboratory (IMCL) of INL. This article summarizes the results of advanced characterization and computer vision-based materials informatics to reveal the irradiation effects on U-Zr metallic fuel. Future work will focus on further implementation of advanced characterization and statistical data mining to improve the fidelity of fuel performance modeling and support U-Zr metallic fuel qualification for fast spectrum reactors.

Yao, Tiankai↗

NASA NTRS2007

Robust Fuzzy Controllers Using FPGAs

Electro-mechanical device controllers typically come in one of three forms, proportional (P), Proportional Derivative (PD), and Proportional Integral Derivative (PID). Two methods of control are discussed in this paper; they are (1) the classical technique that requires an in-depth mathematical use of poles and zeros, and (2) the fuzzy logic (FL) technique that is similar to the way humans think and make decisions. FL controllers are used in multiple industries; examples include control engineering, computer vision, pattern recognition, statistics, and data analysis. Presented is a study on the development of a PD motor controller written in very high speed hardware description language (VHDL), and implemented in FL. Four distinct abstractions compose the FL controller, they are the fuzzifier, the rule-base, the fuzzy inference system (FIS), and the defuzzifier. FL is similar to, but different from, Boolean logic; where the output value may be equal to 0 or 1, but it could also be equal to any decimal value between them. This controller is unique because of its VHDL implementation, which uses integer mathematics. To compensate for VHDL's inability to synthesis floating point numbers, a scale factor equal to 10(sup (N/4) is utilized; where N is equal to data word size. The scaling factor shifts the decimal digits to the left of the decimal point for increased precision. PD controllers are ideal for use with servo motors, where position control is effective. This paper discusses control methods for motion-base platforms where a constant velocity equivalent to a spectral resolution of 0.25 cm(exp -1) is required; however, the control capability of this controller extends to various other platforms.

Monroe, Author Gene S., Jr.↗

DOE OSTI2025

Improved loss functions for machine-learned atomic potentials

Machine learning (ML) has become an invaluable tool across a wide array of domains in science as researchers find new ways to leverage its predictive power. This is especially true in chemistry, where ML is used to fit chemical properties or desirable attributes to the local structure of molecules and materials. In the pursuit of greater accuracy, it is relatively simple to increase the size or complexity of such models, although this often requires simultaneously seeking larger datasets in order to both fit and interpret the larger number of parameters. However, it is equally important to assess the quality and relative importance of the data and how these factors impact the training process. We, therefore, investigate the impact of using different loss functions for training neural network potentials (NNPs), as the loss function defines the error and parameter gradients used to train the NNP. In particular, we test the mean-squared error and Huber loss functions and, using insight from these functions, derive a new loss function based on the Asinh function, which yields significant improvement in the accuracy and generality of NNPs. We show that by discounting/minimizing errors and anomalies in the optimization process, both the Huber and Asinh loss functions improve the training of NNPs, leading to a final potential with a greater effective dimensionality.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Scattering amplitudes for all masses and spins

We introduce a formalism for describing four-dimensional scattering amplitudes for particles of any mass and spin. This naturally extends the familiar spinor-helicity formalism for massless particles to one where these variables carry an extra SU(2) little group index for massive particles, with the amplitudes for spin Sparticles transforming as symmetric rank 2S tensors. We systematically characterise all possible three particle amplitudes compatible with Poincare symmetry. Unitarity, in the form of consistent factorization, imposes algebraic conditions that can be used to construct all possible four-particle tree amplitudes. This also gives us a convenient basis in which to expand all possible four-particle amplitudes in terms of what can be called “spinning polynomials”. Many general results of quantum field theory follow the analysis of four-particle scattering, ranging from the set of all possible consistent theories for massless particles, to spin-statistics, and the Weinberg-Witten theorem. We also find a transparent understanding for why massive particles of sufficiently high spin cannot be “elementary”. The Higgs and Super-Higgs mechanisms are naturally discovered as an infrared unification of many disparate helicity amplitudes into a smaller number of massive amplitudes, with a simple understanding for why this can’t be extended to Higgsing for gravitons. We illustrate a number of applications of the formalism at one-loop, giving few-line computations of the electron (g - 2) as well as the beta function and rational terms in QCD. “Off-shell” observables like correlation functions and form-factors can be thought of as scattering amplitudes with external “probe” particles of general mass and spin, so all these objects — amplitudes, form factors and correlators, can be studied from a common on-shell perspective.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS2019

Challenges in Qualification of Thermal Protection Systems in Extreme Entry Environments

Planetary entry vehicles employ ablative TPS materials to shield the aeroshell from entry aeroheating environments. To ensure mission success, it must be demonstrated that the heat shield system, including local features such as seams, does not fail at conditions that are suitably margined beyond those expected in flight. Furthermore, its thermal response must be predictable, with acceptable fidelity, by computational tools used in heat shield design. Mission assurance is accomplished through a combination of ground testing and material response modelling. A material's robustness to failure is verified through arcjet testing while its thermal response is predicted by analytical tools that are verified against experimental data. Due to limitations in flight-like ground testing capability and lack of validated high-fidelity computational models, qualification of heat shield materials is often achieved by piecing together evidence from multiple ground tests and analytical simulations, none of which fully bound the flight conditions and vehicle configuration. Extreme heating environments (>2000 W/sq. cm heat flux and >2 atm pressure), experienced during entries at Venus, Saturn and Ice Giants, further stretch the current testing and modelling capabilities for applicable TPS materials. Fully-dense Carbon Phenolic was the material of choice for these applications; however, since heritage raw materials are no longer available, future uses of re-created Carbon Phenolic will require re-qualification. To address this sustainability challenge, NASA is developing a new dual-layer material based on 3D weaving technology called Heat shield for Extreme Entry Environments (HEEET). Regardless of TPS material, extreme environments pose additional certification challenges beyond what has been typical in recent NASA missions. Scope of this presentation: This presentation will give an overview of challenges faced in verifying TPS performance at extreme heating conditions. Examples include: (1) Bounding aeroheating parameters (heat flux, pressure, shear and enthalpy) in ground facilities. How to certify TPS if environments can't be bounded or aeroheating parameters can't be simultaneously achieved. (2) Higher uncertainties in ground test environments (facility calibration and analytical predictions) at extreme conditions. (3) Testing in flows similar to planetary atmosphere composition (H2/He for Gas and Ice Giants). (4) Test sample size limitations for qualifying seam designs. (5) Lack of computational tools capable of simulating all significant aspects of TPS performance (including initiation and propagation of failures). This presentation will provide recommendations on how the EDL community can address these challenges and mitigate some of the risks involved in flying TPS materials at extreme conditions. Examples include: (1) Dedicated activity to understanding TPS failure modes. Develop computational tools capable of modelling fluid interaction with material's thermostructural response. Validate these tools through failure testing. A better understanding of failure mechanisms may eliminate the need to fully bound all aeroheating parameters in ground testing. (2) Enhancements to current testing facilities to simulate flight-like ablation mechanism (ex. testing in Nitrogen at Ames Interaction Heating Facility to limit oxidation in favor of more sublimation). (3) Improved characterization of test conditions with new diagnostic methods and determination of environment uncertainty through rigorous statistical analysis of available data. (4) Design margin policies that are directly tied to uncertainties in ground test environments and modelling fidelity

Mahzari, Milad↗

DOE OSTI2020

Computing Potential of the Mean Force Profiles for Ion Permeation Through Channelrhodopsin Chimera, C1C2

Umbrella sampling, coupled with a weighted histogram analysis method (US-WHAM), can be used to construct potentials of mean force (PMFs) for studying the complex ion permeation pathways of membrane transport proteins. Despite the widespread use of US-WHAM, obtaining a physically meaningful PMF can be challenging. Here, we provide a protocol to resolve that issue. Then, we apply that protocol to compute a meaningful PMF for sodium ion permeation through channelrhodopsin chimera, C1C2, for illustration.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2026

Screening and qualification methodology for SiC end plug processing methods

Deployment of SiC-ceramic-based fuel cladding for light water reactors requires a hermetic end plug–to–cladding joint that can withstand neutron irradiation during normal operation and maintain integrity during design-basis accidents. Reactor experiments have shown that some SiC composite tubes with SiC end plugs can retain hermeticity after irradiation. However, achieving consistent joint performance under irradiation remains a key challenge. Resolving this issue is essential to enable integral irradiation testing and to demonstrate fuel integrity under commercial-reactor irradiation conditions. This report aims to: (1) provide guidance for designing radiation-tolerant end plug joints for SiC cladding; (2) demonstrate experimental methods to detect processing defects that are unstable under neutron irradiation at light-water-reactor-relevant temperatures and doses; and (3) outline a step-by-step approach for designing and conducting reactor experiments to screen joining methods. The resulting data will be used to improve joint processing and to define critical defect types and sizes that must be detected and eliminated through non-destructive evaluation for quality assurance. Based on prior irradiation experiments at the High Flux Isotope Reactor, differential swelling among the cladding, bonding layer, and end plug was identified as an underlying mechanism for irradiation-induced joint degradation. Accordingly, this effect must be considered in the design of radiation-tolerant joining techniques. In this work, miniature SiC end plug joint specimens irradiated during the previous project were analyzed using X-ray computed tomography to characterize the joint microstructure. Digital volume correlation of the tomography data quantified radiation-induced microstructural changes and enabled evaluation of defect-related risks. Finally, ongoing neutron irradiation efforts using larger specimen volumes are presented. These efforts aim to statistically assess joint performance and to build a microstructure–performance (e.g., leak-tightness) dataset to inform processing improvements and quality control.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS2002

Risk-Based Probabilistic Approach to Aeropropulsion System Assessment

In an era of shrinking development budgets and resources, where there is also an emphasis on reducing the product development cycle, the role of system assessment, performed in the early stages of an engine development program, becomes very critical to the successful development of new aeropropulsion systems. A reliable system assessment not only helps to identify the best propulsion system concept among several candidates, it can also identify which technologies are worth pursuing. This is particularly important for advanced aeropropulsion technology development programs, which require an enormous amount of resources. In the current practice of deterministic, or point-design, approaches, the uncertainties of design variables are either unaccounted for or accounted for by safety factors. This could often result in an assessment with unknown and unquantifiable reliability. Consequently, it would fail to provide additional insight into the risks associated with the new technologies, which are often needed by decision makers to determine the feasibility and return-on-investment of a new aircraft engine. In this work, an alternative approach based on the probabilistic method was described for a comprehensive assessment of an aeropropulsion system. The statistical approach quantifies the design uncertainties inherent in a new aeropropulsion system and their influences on engine performance. Because of this, it enhances the reliability of a system assessment. A technical assessment of a wave-rotor-enhanced gas turbine engine was performed to demonstrate the methodology. The assessment used probability distributions to account for the uncertainties that occur in component efficiencies and flows and in mechanical design variables. The approach taken in this effort was to integrate the thermodynamic cycle analysis embedded in the computer code NEPP (NASA Engine Performance Program) and the engine weight analysis embedded in the computer code WATE (Weight Analysis of Turbine Engines) with the fast probability integration technique (FPI). FPI was developed by Southwest Research Institute under contract with the NASA Glenn Research Center. The results were plotted in the form of cumulative distribution functions and sensitivity analyses and were compared with results from the traditional deterministic approach. The comparison showed that the probabilistic approach provides a more realistic and systematic way to assess an aeropropulsion system. The current work addressed the application of the probabilistic approach to assess specific fuel consumption, engine thrust, and weight. Similarly, the approach can be used to assess other aspects of aeropropulsion system performance, such as cost, acoustic noise, and emissions. Additional information is included in the original extended abstract.

Tong, Michael T.↗

NASA NTRS1977

Flight inspection data and crack initiation times

Lockheed C-130 service-flight inspection data is reduced by means of fracture mechanics and statistical analysis into a form from which the probability of time to crack initiation and the distribution of initial flaw sizes for various locations on the aircraft can be determined. Crack sizes are normalized by extrapolating a growing crack backwards from its first recorded size to .03 inches (crack initiation size) yielding the time to crack initiation. In a similar fashion, the crack is grown backwards to a time equal to zero to yield an initial flaw size. In order to perform the discussed computations, crack growth constants were computed for each individual location and were analyzed statistically. Statistical distributions were fitted to times to crack initiation and to initial flaw sizes in order to describe their expected behavior. Weibull and Johnson distributions have been found to fit the data reasonably well but for a better fit convolution integrals of Poisson load distributions and Weibull strength distributions are needed.

Johnson, W. S.↗

NASA NTRS1990

Models of space averaged energetics of plates

The analysis of high frequency vibrations in plates is of particular interest in the study of structure borne noise in aircrafts. The current methods of analysis are either too expensive (finite element method) or may have a confidence band wider than desirable (Statistical Energy Analysis). An alternative technique to model the space and time averaged response of structural acoustics problems with enough detail to include all significant mechanisms of energy generation, transmission, and absorption is highly desirable. The focus of this paper is the development of a set of equations which govern the space and time averaged energy density in plates. To solve this equation, a new type of boundary value problem must be treated in terms of energy density variables using energy and intensity boundary conditions. A computer simulation verification study of the energy governing equation is performed. A finite element formulation of the new equations is also implemented and several test cases are analyzed and compared to analytical solutions.

Bouthier, O. M.↗

DOE OSTI2022

Computing Thermodynamic Properties of Fluids Augmented by Nanoconfinement: Application to Pressurized Methane

Nanoconfined fluids exhibit remarkably different thermodynamic behavior compared to the bulk phase. These confinement effects render predictions of thermodynamic quantities of nanoconfined fluids challenging. In particular, confinement creates a spatially varying density profile near the wall that is primarily responsible for adsorption and capillary condensation behavior. Significant fluctuations in thermodynamic quantities, inherent in such nanoscale systems, coupled to strong fluid–wall interactions give rise to this near-wall density profile. Empirical models have been proposed to explain and model these effects, yet no first-principles based formulation has been developed. We present a statistical mechanics framework that embeds such a coupling to describe the effect of the fluid–wall interaction in amplifying the near-wall density behavior for compressible gases at elevated pressures such as pressurized methane in confinement. We show that the proposed theory predicts accurately the adsorbed layer thickness as obtained with small-angle neutron scattering measurements. Furthermore, the predictions of density under confinement from the proposed theory are shown to be in excellent agreement with available experimental and atomistic simulations data for a range of temperatures for nanoconfined methane. While the framework is presented for evaluating the near-wall density, owing to its rigorous foundation in statistical mechanics, the proposed theory can also be generalized for predicting phase-transition and nonequilibrium transport of nanoconfined fluids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1985

Statistical mechanics of light elements at high pressure. VII - A perturbative free energy for arbitrary mixtures of H and He

A model free energy is presented which accurately represents results from 45 high-precision Monte Carlo calculations of the thermodynamics of hydrogen-helium mixtures at pressures of astrophysical and planetophysical interest. The free energy is calculated using free-electron perturbation theory (dielectric function theory), and is an extension of the expression given in an earlier paper in this series. However, it fits the Monte Carlo results more accurately, and is valid for the full range of compositions from pure hydrogen to pure helium. Using the new free energy, the phase diagram of mixtures of liquid metallic hydrogen and helium is calculated and compared with earlier results. Sample results for mixing volumes are also presented, and the new free energy expression is used to compute a theoretical Jovian adiabat and compare the adiabat with results from three-dimensional Thomas-Fermi-Dirac theory. The present theory gives slightly higher densities at pressures of about 10 megabars.

Hubbard, W. B.↗