Search NASA⌕ Search

SEARCH · Search NASA

Results for “Maximum Entropy Model”

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 235 records · Page 13

Hidden first-order phase transitions and large magnetocaloric effects in GdNi 1-x Co x

We report a large magnetocaloric effect tunable between 70 and 115 K in a series of rare earth intermetallic compounds GdNi 1-x Co x with x varying from 0 to 0.15 that manifest rather unconventional second-order phase transitions with concurrent changes in both magnetic and crystallographic sublattices. While the Curie temperature, TC, linearly increases with x(Co), the maximum isothermal entropy change induced by magnetic field varying between 0 and 50 kOe slightly decreases from 14 J/kg K at 71 K when x(Co) = 0 to10 J/kg K at 115 K when x(Co) = 0.15. Additionally, the temperature-averaged entropy change figures of merit calculated at temperature spans of 10 K for all examined compositions are comparable to those of the best known magnetocaloric materials that exhibit second-order phase transitions, including elemental Gd and La(Fe,Co,Si) 13 alloys. Through a detailed analysis of the magnetic field-induced entropy changes along with assessment of critical exponents we explore the role magneto-elastic coupling plays in controlling magnetocaloric properties of GdNi 1-x Co x compounds. We also formulate a model based on the mean field theory approximation to describe both magnetic and magnetocaloric properties of the title materials in qualitative agreement with experimental results.

36 MATERIALS SCIENCE↗

Granular Contact Forces: Proof of "Self-Ergodicity" by Generalizing Boltzmann's Stosszahlansatz and H Theorem

Ergodicity is proved for granular contact forces. To obtain this proof from first principles, this paper generalizes Boltzmann's stosszahlansatz (molecular chaos) so that it maintains the necessary correlations and symmetries of granular packing ensembles. Then it formally counts granular contact force states and thereby defines the proper analog of Boltzmann's H functional. This functional is used to prove that (essentially) all static granular packings must exist at maximum entropy with respect to their contact forces. Therefore, the propagation of granular contact forces through a packing is a truly ergodic process in the Boltzmannian sense, or better, it is self-ergodic. Self-ergodicity refers to the non-dynamic, internal relationships that exist between the layer-by-layer and column-by-column subspaces contained within the phase space locus of any particular granular packing microstate. The generalized H Theorem also produces a recursion equation that may be solved numerically to obtain the density of single particle states and hence the distribution of granular contact forces corresponding to the condition of self-ergodicity. The predictions of the theory are overwhelmingly validated by comparison to empirical data from discrete element modeling.

Metzger, Philip T.↗

Foundational Performance Analyses of Pressure Gain Combustion Thermodynamic Benefits for Gas Turbines

A methodology is described whereby the work extracted by a turbine exposed to the fundamentally nonuniform flowfield from a representative pressure gain combustor (PGC) may be assessed. The method uses an idealized constant volume cycle, often referred to as an Atkinson or Humphrey cycle, to model the PGC. Output from this model is used as input to a scalable turbine efficiency function (i.e., a map), which in turn allows for the calculation of useful work throughout the cycle. Integration over the entire cycle yields mass-averaged work extraction. The unsteady turbine work extraction is compared to steady work extraction calculations based on various averaging techniques for characterizing the combustor exit pressure and temperature. It is found that averages associated with momentum flux (as opposed to entropy or kinetic energy) provide the best match. This result suggests that momentum-based averaging is the most appropriate figure-of-merit to use as a PGC performance metric. Using the mass-averaged work extraction methodology, it is also found that the design turbine pressure ratio for maximum work extraction is significantly higher than that for a turbine fed by a constant pressure combustor with similar inlet conditions and equivalence ratio. Limited results are presented whereby the constant volume cycle is replaced by output from a detonation-based PGC simulation. The results in terms of averaging techniques and design pressure ratio are similar.

Paxson, Daniel E.↗

Statistical Significance Testing for Mixed Priors: A Combined Bayesian and Frequentist Analysis

In many hypothesis testing applications, we have mixed priors, with well-motivated informative priors for some parameters but not for others. The Bayesian methodology uses the Bayes factor and is helpful for the informative priors, as it incorporates Occam’s razor via the multiplicity or trials factor in the look-elsewhere effect. However, if the prior is not known completely, the frequentist hypothesis test via the false-positive rate is a better approach, as it is less sensitive to the prior choice. We argue that when only partial prior information is available, it is best to combine the two methodologies by using the Bayes factor as a test statistic in the frequentist analysis. We show that the standard frequentist maximum likelihood-ratio test statistic corresponds to the Bayes factor with a non-informative Jeffrey’s prior. We also show that mixed priors increase the statistical power in frequentist analyses over the maximum likelihood test statistic. We develop an analytic formalism that does not require expensive simulations and generalize Wilks’ theorem beyond its usual regime of validity. In specific limits, the formalism reproduces existing expressions, such as the p-value of linear models and periodograms. We apply the formalism to an example of exoplanet transits, where multiplicity can be more than 10 7 . We show that our analytic expressions reproduce the $p$-values derived from numerical simulations. We offer an interpretation of our formalism based on the statistical mechanics. We introduce the counting of states in a continuous parameter space using the uncertainty volume as the quantum of the state. We show that both the $p$-value and Bayes factor can be expressed as an energy versus entropy competition.

97 MATHEMATICS AND COMPUTING↗

Spray Characteristics of a Hybrid Twin-Fluid Pressure-Swirl Atomizer

The spray performance of a fuel injection system applicable for use in main combustion chamber of an oxidizer-rich staged combustion (ORSC) cycles is presented. The experimental data reported here include mean drop size and drop size distribution, spray cone half-angle, and momentum rate (directly related to spray penetration). The maximum entropy formalism, MEF, method to predict drop size distribution is applied and compared to the experimental data. Geometric variables considered include the radius of the injector inlet orifice plate through which oxidizer flows (&) and the exposed length from the fuel inlet to the injector exit plane (L2). Operating conditions that were varied include the liquid mass flow rate and air mass flow rate. For orifices B and C there is a significant dependence of D3Z on both the air and liquid mass flow rates, as well as on L2. For the A orifice, the momentum rate of the air flow appears to exceed a threshold value above which a constant D32 is obtained. Using the MEF method, a semi-analytical process was developed to model the spray distribution using two input parameters (q = 0.4 and Dso). The momentum rate of the spray is directly related to the air and liquid mass flow rates. The cone half angle of the spray ranges from 25 to 17 degrees. The data resulting from this project will eventually be used to develop advanced rocket systems.

Durham, M. J.↗

Fabrication of Engineered Silica-Functionalized-Polyanilines Nanocomposites for Water Decontamination of Cadmium and Lead

In this study, two novel functionalized silica nanocomposites were synthesized via covalent bonding with nanopolyaniline (NPANI) and crosslinked nanopolyaniline (CrossNPANI) to produce [NSi-NPANI] and [NSi-CrossNPANI], respectively. The two nanocomposites were portrayed by SEM, FT-IR, HR-TEM, BET-surface and TGA. The HR-TEM images of [NSi-NPANI] and [NSi-CrossNPANI] confirmed the particle size in the range 14.28–21.43 and 26.19–35.71 nm, respectively and these two nanocomposites were successfully applied to remove divalent cadmium and lead from solutions. As compared to nanosilica [NSi] sorbent, the maximum capacity values of the two nanocomposites for divalent cadmium and lead were identified to increase from 100–250 µmole g{sup −1} (pH 1.0) to 750–800 µmole g{sup −1} (pH 7.0) and from 400–1050 µmole g{sup −1} (pH 1.0) to 650–1350 µmole g{sup −1} (pH 6.0), at the optimum conditions. The adsorption data were compared using two and three parameter equations based on Langmuir, Freundlich, Temkin, Dubinin-Radushkevich, Redlish-Peterson and sips models. The kinetic study of [NSi-NPANI] and [NSi-CrossNPANI] revealed that the pseudo-second order kinetic was the best model to explain the kinetic data for cadmium(II) and lead(II). The free energy of sorption (ΔG{sup o}), enthalpy (ΔH{sup o}), and entropy (ΔS{sup o}) changes were calculated to portend the nature of adsorption. Adsorptive extraction of toxic lead and cadmium from tap water and wastewater was successfully accomplished with the range of percentage recovery values 98.5–100.0 and 91.9–93.0%, respectively.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Discretization and Preconditioning Algorithms for the Euler and Navier-Stokes Equations on Unstructured Meshes

Several stabilized demoralization procedures for conservation law equations on triangulated domains will be considered. Specifically, numerical schemes based on upwind finite volume, fluctuation splitting, Galerkin least-squares, and space discontinuous Galerkin demoralization will be considered in detail. A standard energy analysis for several of these methods will be given via entropy symmetrization. Next, we will present some relatively new theoretical results concerning congruence relationships for left or right symmetrized equations. These results suggest new variants of existing FV, DG, GLS, and FS methods which are computationally more efficient while retaining the pleasant theoretical properties achieved by entropy symmetrization. In addition, the task of Jacobean linearization of these schemes for use in Newton's method is greatly simplified owing to exploitation of exact symmetries which exist in the system. The FV, FS and DG schemes also permit discrete maximum principle analysis and enforcement which greatly adds to the robustness of the methods. Discrete maximum principle theory will be presented for general finite volume approximations on unstructured meshes. Next, we consider embedding these nonlinear space discretizations into exact and inexact Newton solvers which are preconditioned using a nonoverlapping (Schur complement) domain decomposition technique. Elements of nonoverlapping domain decomposition for elliptic problems will be reviewed followed by the present extension to hyperbolic and elliptic-hyperbolic problems. Other issues of practical relevance such the meshing of geometries, code implementation, turbulence modeling, global convergence, etc, will. be addressed as needed.

Barth, Timothy J.↗

Water anomalous thermodynamics, attraction, repulsion, and hydrophobic hydration

A model composed of van der Waals-like and hydrogen bonding contributions that simulates the low-temperature anomalous thermodynamics of pure water while exhibiting a second, liquid-liquid critical point [P. H. Poole et al., Phys. Rev. Lett. 73, 1632 (1994)] is extended to dilute solutions of nonionic species. Critical lines emanating from such second critical point are calculated. While one infers that the smallness of the water molecule may be a relevant factor for those critical lines to move towards experimentally accessible regions, attention is mainly focused on the picture our model draws for the hydration thermodynamics of purely hydrophobic and amphiphilic non-electrolyte solutes. We first focus on differentiating solvation at constant volume from the corresponding isobaric process. Both processes provide the same viewpoint for the low solubility of hydrophobic solutes: it originates from the combination of weak solute-solvent attractive interactions and the specific excluded-volume effects associated with the small molecular size of water. However, a sharp distinction is found when exploring the temperature dependence of hydration phenomena since, in contrast to the situation for the constant-V process, the properties of pure water play a crucial role at isobaric conditions. Specifically, the solubility minimum as well as enthalpy and entropy convergence phenomena, exclusively ascribed to isobaric solvation, are closely related to water’s density maximum. Furthermore, the behavior of the partial molecular volume and the partial molecular isobaric heat capacity highlights the interplay between water anomalies, attraction, and repulsion. The overall picture presented here is supported by experimental observations, simulations, and previous theoretical results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Lattice NRQCD study on in-medium bottomonium spectra using a novel Bayesian reconstruction approach

We present recent results on the in-medium modification of S- and P-wave bottomonium states around the deconfinement transition. Our study uses lattice QCD with N{sub f} = 2 + 1 light quark flavors to describe the non-perturbative thermal QCD medium between 140MeV < T < 249MeV and deploys lattice regularized non-relativistic QCD (NRQCD) effective field theory to capture the physics of heavy quark bound states immersed therein. The spectral functions of the {sup 3}S{sub 1} (ϒ) and {sup 3}P{sub 1} (χ{sub b1}) bottomonium states are extracted from Euclidean time Monte Carlo simulations using a novel Bayesian prescription, which provides higher accuracy than the Maximum Entropy Method. Based on a systematic comparison of interacting and free spectral functions we conclude that the ground states of both the S-wave (ϒ) and P-wave (χ{sub b1}) channel survive up to T = 249MeV. Stringent upper limits on the size of the in-medium modification of bottomonium masses and widths are provided.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Heat capacity and magnetic properties of fluoride CsFe{sup 2+}Fe{sup 3+}F{sub 6} with defect pyrochlore structure

Heat capacity, Mössbauer and Raman spectra as well as magnetic properties of fluoride CsFe{sub 2}F{sub 6} with defect pyrochlore structure were studied. In addition to recently found above room temperature three successive structural transformations Pnma-Imma-I4{sub 1}amd-Fd-3m, phase transition of antiferromagnetic nature with the 13.7 K Neel temperature and a broad heat capacity anomaly with a maximum at about 30 K were observed. The room temperature symmetry Pnma is unchanged at least down to 7 K. Simple model of indirect bond used to estimate the exchange interactions and to propose a magnetic structure model. - Graphical abstract: The ordered arrangement of Fe{sup 2+} and Fe{sup 3+} ions in high-spin states as well as antiferromagnetic phase transition followed by significant magnetic frustrations were found in pyrocholore-related CsFe{sup 2+}Fe{sup 3+}F{sub 6}. A magnetic structure was proposed using a simple model of indirect bonds. - Highlights: • The Pnma structure in pyrocholore CsFe{sub 2}F{sub 6} is stable down to helium temperature. • Mössbauer spectra confirmed the ordering of Fe{sup 2+} and Fe{sup 3+} ions. • Antiferromagnetic transformation and significant magnetic frustrations are found. • Experimental magnetic entropy agrees with entropy for Fe ions in high-spin state. • Superexchange interactions were calculated and a magnetic structure was proposed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Discretization and Preconditioning Algorithms for the Euler and Navier-Stokes Equations on Unstructured Meshes

Several stabilized discretization procedures for conservation law equations on triangulated domains will be considered. Specifically, numerical schemes based on upwind finite volume, fluctuation splitting, Galerkin least-squares, and space discontinuous Galerkin discretization will be considered in detail. A standard energy analysis for several of these methods will be given via entropy symmetrization. Next, we will present some relatively new theoretical results concerning congruence relationships for left or right symmetrized equations. These results suggest new variants of existing FV, DG, GLS and FS methods which are computationally more efficient while retaining the pleasant theoretical properties achieved by entropy symmetrization. In addition, the task of Jacobian linearization of these schemes for use in Newton's method is greatly simplified owing to exploitation of exact symmetries which exist in the system. These variants have been implemented in the "ELF" library for which example calculations will be shown. The FV, FS and DG schemes also permit discrete maximum principle analysis and enforcement which greatly adds to the robustness of the methods. Some prevalent limiting strategies will be reviewed. Next, we consider embedding these nonlinear space discretizations into exact and inexact Newton solvers which are preconditioned using a nonoverlapping (Schur complement) domain decomposition technique. Elements of nonoverlapping domain decomposition for elliptic problems will be reviewed followed by the present extension to hyperbolic and elliptic-hyperbolic problems. Other issues of practical relevance such the meshing of geometries, code implementation, turbulence modeling, global convergence, etc. will be addressed as needed.

Barth, Timothy↗

A crisis for the verification and validation of turbulence simulations

Three algorithms have been proposed for solution of the Rayleigh–Taylor turbulent mixing problem. They are based upon three different physical principles governing the Euler equations for fluid flow. The principles serve to select the physically relevant solution from among many nonunique solutions. The admissibility principle is in dispute. The three different algorithms, expressing the three physical admissibility principles can be formulated in terms of the three energy dissipation rates or the entropy production rates, as selected by the size of the sub grid scale coefficients. These have maximal values or less than maximal values. The resulting solutions are markedly different. We find strong validation evidence that supports the maximum rate principle, based on a review of prior results and on new results presented here. We review experimental data used for validation and sufficient to discriminate among the three. We present a new analysis of this data. We show that the hypothesized long wave length perturbations in the initial conditions are not significant, so that validation can be based on this data in a straight forward manner. One of the algorithms is labeled direct numerical simulation, but is not, and as a consequence, the two algorithms with less than maximal SGS coefficients are variants of one another. Finally recommendations for the numerical modeling of the deflagration to detonation transition in type Ia supernova are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

An entropy-based debiasing approach to quantifying experimental coverage for novel applications of interest in the nuclear community

This manuscript proposes a novel information-theoretic approach to the quantification of experimental relevance, i.e., coverage, to achieve optimal data assimilation results for nuclear engineering applications. Specifically, this work posits the need for a new metric, called coverage (q C ) of an application’s quantity of interest, i.e., eigenvalue or power peaking for an advanced reactor concept, defined herein as the theoretically maximum achievable reduction in the quantity’s uncertainty given measurements from a pool of experiments in a manner that is independent of the data assimilation procedure employed. Currently, reduction in a quantity’s uncertainty is strongly biased by the underlying assumptions of the assimilation procedure to account for the under-determined nature of such problems and the similarity criterion employed to identify relevant experiments. To address this challenge, this work has developed a coverage metric, q C , based on mutual information, which establishes a new conceptual framework for assessing coverage, one that is independent of the model parameters and responses degree of variations in both the experimental and application domains, i.e., linear vs non-linear, and their prior uncertainty distributions, i.e., Gaussian vs. non-Gaussian. The q C is an entropic measure capable of addressing coverage for general nonlinear problems with non-Gaussian uncertainties and inclusive of the measurement uncertainties from multiple experiments. Numerical experiments from manufactured analytical problems as well as a set of benchmarks from the ICSBEP handbook are employed to demonstrate its theoretical and practical performance as compared to the c k -based experiment selection methodology, commonly employed in the neutronic community. The manuscript then employs other well-known adaptations to existing data assimilation methodologies for nonlinear and non-Gaussian problems capable of achieving the coverage posited by q C .

Bayesian data assimilation↗

Absolute entropy and the observer’s no-boundary state

We investigate the no-boundary proposal for closed universes with an observer. We argue that the observer’s no-boundary state is the identity operator on the physical Hilbert space, i.e., the maximum entropy state and show this explicitly in Jackiw-Teitelboim gravity. Geometrically, the no-boundary state is a bra-ket wormhole. Expectation values in the no-boundary state provide a trace for the observer’s algebra, which allows one to define von Neumann entropy for observers in different universes as the relative entropy with respect to the no-boundary state. This result is consistent with all previously discussed cases of traces for invariantly defined regions.

Cosmological models↗

Enhancement and maximum in the isobaric specific-heat capacity measurements of deeply supercooled water using ultrafast calorimetry

Knowledge of the temperature dependence of the isobaric specific heat (C p ) upon deep supercooling can give insights regarding the anomalous properties of water. If a maximum in C p exists at a specific temperature, as in the isothermal compressibility, it would further validate the liquid–liquid critical point model that can explain the anomalous increase in thermodynamic response functions. The challenge is that the relevant temperature range falls in the region where ice crystallization becomes rapid, which has previously excluded experiments. Here, we have utilized a methodology of ultrafast calorimetry by determining the temperature jump from femtosecond X-ray pulses after heating with an infrared laser pulse and with a sufficiently long time delay between the pulses to allow measurements at constant pressure. Evaporative cooling of ~15-µm diameter droplets in vacuum enabled us to reach a temperature down to ~228 K with a small fraction of the droplets remaining unfrozen. We observed a sharp increase in C p , from 88 J/mol/K at 244 K to about 218 J/mol/K at 229 K where a maximum is seen. The C p maximum is at a similar temperature as the maxima of the isothermal compressibility and correlation length. From the C p measurement, we estimated the excess entropy and self-diffusion coefficient of water and these properties decrease rapidly below 235 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

ZEEMAN DOPPLER MAPS: ALWAYS UNIQUE, NEVER SPURIOUS?

Numerical models of atomic diffusion in magnetic atmospheres of ApBp stars predict abundance structures that differ from the empirical maps derived with (Zeeman) Doppler mapping (ZDM). An in-depth analysis of this apparent disagreement investigates the detectability by means of ZDM of a variety of abundance structures, including (warped) rings predicted by theory, but also complex spot-like structures. Even when spectra of high signal-to-noise ratio are available, it can prove difficult or altogether impossible to correctly recover shapes, positions, and abundances of a mere handful of spots, notwithstanding the use of all four Stokes parameters and an exactly known field geometry; the recovery of (warped) rings can be equally challenging. Inversions of complex abundance maps that are based on just one or two spectral lines usually permit multiple solutions. It turns out that it can by no means be guaranteed that any of the regularization functions in general use for ZDM (maximum entropy or Tikhonov) will lead to a true abundance map instead of some spurious one. Attention is drawn to the need for a study that would elucidate the relation between the stratified, field-dependent abundance structures predicted by diffusion theory on the one hand, and empirical maps obtained by means of “canonical” ZDM, i.e., with mean atmospheres and unstratified abundances, on the other hand. Finally, we point out difficulties arising from the three-dimensional nature of the atomic diffusion process in magnetic ApBp star atmospheres.

79 ASTRONOMY AND ASTROPHYSICS↗

Experimental X-ray Charge-Density Studies–A Suitable Probe for Superconductivity? A Case Study on MgB 2

Case studies of 1T-TiSe 2 and YBa 2 Cu 3 O 7-δ have demonstrated that X-ray diffraction (XRD) studies can be used to trace even subtle structural phase transitions which are inherently connected with the onset of superconductivity in these benchmark systems. However, the utility of XRD in the investigation of superconductors like MgB 2 lacking an additional symmetry-breaking structural phase transition is not immediately evident. Nevertheless, high-resolution powder XRD experiments on MgB 2 in combination with maximum entropy method analyses hinted at differences between the electron density distributions at room temperature and 15 K, that is, below the T c of approx. 39 K. The high-resolution single-crystal XRD experiments in combination with multipolar refinements presented here can reproduce these results but show that the observed temperature-dependent density changes are almost entirely due to a decrease of atomic displacement parameters as a natural consequence of a reduced thermal vibration amplitude with decreasing temperature. Our investigations also shed new light on the presence or absence of magnesium vacancies in MgB 2 samples–a defect type claimed to control the superconducting properties of the compound. Here, we propose that previous reports on the tendency of MgB 2 to form non-stoichiometric Mg 1–x B 2 phases (1 – x ~ 0.95) during high-temperature (HT) synthesis might result from the interpretation of XRD data of insufficient resolution and/or usage of inflexible refinement models. Indeed, advanced refinements based on an Extended Hansen–Coppens multipolar model and high-resolution X-ray data, which consider explicitly the contraction of core and valence shells of the magnesium cations, do not provide any significant evidence for the formation of non-stoichiometric Mg 1–x B 2 phases during HT synthesis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Communication: Improved ab initio molecular dynamics by minimally biasing with experimental data

Accounting for electrons and nuclei simultaneously is a powerful capability of ab initio molecular dynamics (AIMD). However, AIMD is often unable to accurately reproduce properties of systems such as water due to inaccuracies in the underlying electronic density functionals. This shortcoming is often addressed by added empirical corrections and/or increasing the simulation temperature. We present here a maximum-entropy approach to directly incorporate limited experimental data via a minimal bias. Biased AIMD simulations of water and an excess proton in water are shown to give significantly improved properties both for observables which were biased to match experimental data and for unbiased observables. Here, this approach also yields new physical insight into inaccuracies in the underlying density functional theory as utilized in the unbiased AIMD.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗