Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Inference mechanisms”

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 757 records · Page 42

Mineralogical Characterization of Fe-Bearing AGB and Supernova Silicate Grains From the Queen Alexandra Range 99177 Meteorite

Spectroscopic observations of the circumstellar envelopes of evolved O-rich stars indicate the dust is mostly amorphous silicate with olivine-like compositions. Spectral modeling suggests these grains are Fe-rich [Mg/(Mg+Fe) ~0.5], but it is not known whether the Fe is distributed within the silicate matrix or exists as metal inclusions. In contrast, the crystalline silicates are inferred to be extremely Mg-rich [Mg/(Mg+Fe) > 0.95]. The mineralogies and chemical compositions of dust in supernova (SN) remnants are not as well constrained, but abundant silicates of olivine-like and enstatite-like compositions have been fit to the infrared emission features. Silicates in the interstellar medium (ISM) are >99% amor-phous and Fe-bearing. The dearth of crystalline silicates in the ISM requires that some amorphization or destruction mechanisms process these grains.

Nguyen, A. N.↗

Optimal experimental design for prediction based on push-forward probability measures

Incorporating experimental data is essential for increasing the credibility of simulation-aided decision making and design. This paper presents a method which uses a computational model to guide the optimal acquisition of experimental data to produce data-informed predictions of quantities of interest (QoI). Many strategies for optimal experimental design (OED) select data that maximize some utility that measures the reduction in uncertainty of uncertain model parameters, for example the expected information gain between prior and posterior distributions of these parameters. In this paper, we seek to maximize the expected information gained from the pushforward of an initial (prior) density to the push-forward of the updated (posterior) density through the parameter-to-prediction map. The formulation presented is based upon the solution of a specific class of stochastic inverse problems which seeks a probability density that is consistent with the model and the data in the sense that the push-forward of this density through the parameter-to-observable map matches a target density on the observable data. While this stochastic inverse problem forms the mathematical basis for our approach, we develop a one-step algorithm, focused on push-forward probability measures, that leverages inference-for-prediction to bypass constructing the solution to the stochastic inverse problem. A number of numerical results are presented to demonstrate the utility of this optimal experimental design for prediction and facilitate comparison of our approach with traditional OED.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

nPINNs: nonlocal Physics-Informed Neural Networks for a parametrized nonlocal universal Laplacian operator. Algorithms and Applications

Physics-informed neural networks (PINNs) are effective in solving inverse problems based on differential and integro-differential equations with sparse, noisy, unstructured, and multifidelity data. PINNs incorporate all available information, including governing equations (reflecting physical laws), initial-boundary conditions, and observations of quantities of interest, into a loss function to be minimized, thus recasting the original problem into an optimization problem. In this paper, we extend PINNs to parameter and function inference for integral equations such as nonlocal Poisson and nonlocal turbulence models, and we refer to them as nonlocal PINNs (nPINNs). The contribution of the paper is three-fold. First, we propose a unified nonlocal Laplace operator, which converges to the classical Laplacian as one of the operator parameters, the nonlocal interaction radius δ goes to zero, and to the fractional Laplacian as δ goes to infinity. This universal operator forms a super-set of classical Laplacian and fractional Laplacian operators and, thus, has the potential to fit a broad spectrum of data sets. We provide theoretical convergence rates with respect to δ and verify them via numerical experiments. Second, we use nPINNs to estimate the two parameters, δ and α, characterizing the kernel of the unified operator. The strong non-convexity of the loss function yielding multiple (good) local minima reveals the occurrence of the operator mimicking phenomenon, that is, different pairs of estimated parameters could produce multiple solutions of comparable accuracy. Third, we propose another nonlocal operator with spatially variable order α(γ), which is more suitable for modeling turbulent Couette flow. Our results show that nPINNs can jointly infer this function as well as δ. More importantly, these parameters exhibit a universal behavior with respect to the Reynolds number, a finding that contributes to our understanding of nonlocal interactions in wall-bounded turbulence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Interstitial hydrogen enhances the mobility of some grain boundaries in tungsten

Segregation of interstitials at a grain boundary (GB) is known to generally lower its mobility. This phenomenon, called 'solute-drag', has important ramifications on the process of recrystallization and microstructural evolution. Here, we present predictions from molecular dynamics (MD) simulations which demonstrate that interstitial hydrogen in tungsten can in fact increase the mobility of some GBs which exhibit shear coupling. Assuming a disconnection-based mechanism, activation energies and pre-factors for disconnection nucleation are predicted from simulations of shear-coupled motion. In GBs where enhanced mobility is predicted, interstitial H reduces both the activation energy and the pre-factor for disconnection nucleation, thus effectively increasing the mobility. For GBs with diminished mobility, MD predicts that presence of interstitial H reduces the pre-factor and, in some cases, increases the activation energy. The reduction in the activation energy inferred from MD simulations are confirmed by nudged elastic band calculations. Temperature-dependent structural transitions are observed for some GBs, and the effect of interstitial H is found to change with the changes in structure. The effect of interstitial H is predicted to be complex and highly variable, providing some plausible explanations for experimental observations on the recrystallization of tungsten in presence of H-loaded plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Generic Hilbert space fragmentation in Kogut-Susskind lattice gauge theories

At the heart of quantum many-body physics lies the understanding of mechanisms that avoid quantum thermalization in an isolated system quenched far from equilibrium. A prominent example is Hilbert space fragmentation, which has recently emerged as an ergodicity-breaking mechanism in constrained spin models. Here, we show that Kogut-Susskind formulations of lattice gauge theories in d+1D (d spatial and one temporal dimensions) give rise to Hilbert space fragmentation, and discuss possible implications for understanding continuum physics. Lastly, our findings not only prove that lattice gauge theories are a natural platform for Hilbert space fragmentation, they also serve as a guide to the conditions under which these models can be faithfully used to infer the thermalization properties of quantum chromodynamics.

Eigenstate thermalization↗

Environmental controls on the light use efficiency of terrestrial gross primary production

Abstract Gross primary production (GPP) by terrestrial ecosystems is a key quantity in the global carbon cycle. The instantaneous controls of leaf‐level photosynthesis are well established, but there is still no consensus on the mechanisms by which canopy‐level GPP depends on spatial and temporal variation in the environment. The standard model of photosynthesis provides a robust mechanistic representation for C 3 species; however, additional assumptions are required to “scale up” from leaf to canopy. As a consequence, competing models make inconsistent predictions about how GPP will respond to continuing environmental change. This problem is addressed here by means of an empirical analysis of the light use efficiency (LUE) of GPP inferred from eddy covariance carbon dioxide flux measurements, in situ measurements of photosynthetically active radiation (PAR), and remotely sensed estimates of the fraction of PAR (fAPAR) absorbed by the vegetation canopy. Focusing on LUE allows potential drivers of GPP to be separated from its overriding dependence on light. GPP data from over 100 sites, collated over 20 years and located in a range of biomes and climate zones, were extracted from the FLUXNET2015 database and combined with remotely sensed fAPAR data to estimate daily LUE. Daytime air temperature, vapor pressure deficit, diffuse fraction of solar radiation, and soil moisture were shown to be salient predictors of LUE in a generalized linear mixed‐effects model. The same model design was fitted to site‐based LUE estimates generated by 16 terrestrial ecosystem models. The published models showed wide variation in the shape, the strength, and even the sign of the environmental effects on modeled LUE. These findings highlight important model deficiencies and suggest a need to progress beyond simple “goodness of fit” comparisons of inferred and predicted carbon fluxes toward an approach focused on the functional responses of the underlying dependencies.

54 ENVIRONMENTAL SCIENCES↗

Flexural analysis of uplifted rift flanks on Venus

Knowledge of the thermal structure of a planet is vital to a thorough understanding of its general scheme of tectonics. Since no direct measurements of heat flow or thermal gradient are available for Venus, most estimates have been derived from theoretical considerations or by analog with the Earth. The flexural response of the lithosphere to applied loads is sensitive to regional thermal structure. Under the assumption that the yield strength as a function of depth can be specified, the temperature gradient can be inferred from the effective elastic plate thickness. Previous estimates of the effective elastic plate thickness of Venus range from 11-18 km for the foredeep north of Uorsar Rupes to 30-60 km for the annular troughs around several coronae. Thermal gradients inferred for these regions are 14-23 K km(exp -1) and 4-9 K km(exp -1) respectively. In this study, we apply the same techniques to investigate the uplifted flanks of an extensional rift. Hypotheses for the origin of uplifted rift flanks on Earth include lateral transport of heat from the center of the rift, vertical transport of heat by small-scale convection, differential thinning of the lithosphere, dynamical uplift, and isostatic response to mechanical uploading of the lithosphere. The 1st hypothesis is considered the dominant contributor to terrestrial rift flanks lacking evidence for volcanic activity, particularly for rift structures that are no longer active. In this study, we model the uplifted flanks of a venusian rift as the flexural response to a vertical end load.

Evans, Susan A.↗

Chemical transitions for interstellar C2 and CN in cloud envelopes

Observations were made of absorption from CH, C2, and CN toward moderately reddened stars in Sco, OB2, Ceo OB3, and Taurus/Auriga. For these directions, most of the reddening is associated with a single cloud complex, for example, the rho Ophiuchus molecular cloud, and as a result, the observations probe moderately dense material. When combined with avaliable data for nearby directions, the survey provides the basis for a comprehensive analysis of the chemistry for these species. The chemical transitions affecting C2 and CN in cloud envelopes were analyzed. The depth into a cloud at which a transition takes place was characterized by tau(sub uv), the grain optical depth at 1000 A. One transition at tau(sub uv) approx. = 2, which arises from, the conversion of C(+) into CO, affects the chemistries for both molecules because of the key role this ion plays. A second one involving production terms in the CN chemistry occurs at tau(sub uv) of approx. = 3; neutral reactions which C2 and CH is more important at larger values for tau(sub uv). The transition from photodissociation to chemical destruction takes place at tau(sub uv) approx. = 4.5 for C2 and CN. The observational data for stars in Sco OB2, Cep OB3, and Taurus/Auriga were studied with chemical rate equations containing the most important production and destruction mechanisms. Because the sample of stars in Sco OB2 includes sight lines with A(sub v) ranging from 1-4 mag, sight lines dominated by photochemistry could be analyzed separately from those controlled by gas-phase destruction. The analysis yielded values for two poorly known rate constants for reactions involved in the production of CN; the reactions are C2 + N yields CN + C and C(+) + NH yields all products. The other directions were analyzed with the inferred values. The predicted column densities for C2 and CN agree with the observed values to better than 50%, and in most instances 20%. When combining the estimates for density and temperature derived from chemical modeling and molecular excitation for a specific cloud, such as the rho Ophiuchus molecular cloud, the portion of the cloud envelope probed by C2 and CN absorption was found to be in pressure equilibrium.

Federman, S. R.↗

Structure and Early Soot Oxidation Properties of Laminar Diffusion Flames

Soot is an important unsolved problem of combustion science because it is present in most hydrocarbon-fueled flames and current understanding of the reactive and physical properties of soot in flame environments is limited. This lack of understanding affects progress toward developing reliable predictions of flame radiation properties, reliable predictions of flame pollutant emission properties and reliable methods of computational combustion, among others. Motivated by these observations, the present investigation extended past studies of soot formation in this laboratory, to consider soot oxidation in laminar diffusion flames using similar methods. Early work showed that O2 was responsible for soot oxidation in high temperature O2-rich environments. Subsequent work in high temperature flame environments having small O2 concentrations, however, showed that soot oxidation rates substantially exceeded estimates based on the classical O2 oxidation rates of Nagle and Strickland-Constable and suggests that radicals such as O and OH might be strong contributors to soot oxidation for such conditions. Neoh et al. subsequently made observations in premixed flames, supported by later work, that showed that OH was responsible for soot oxidation at these conditions with a very reasonable collision efficiency of 0.13. Subsequent studies in diffusion flames, however, were not in agreement with the premixed flame studies: they agreed that OH played a dominant role in soot oxidation in flames, but found collision efficiencies that varied with flame conditions and were not in good agreement with each other or with Neoh et al. One explanation for these discrepancies is that optical scattering and extinction properties were used to infer soot structure properties for the studies that have not been very successful for representing the optical properties of soot. Whatever the source of the problem, however, these differences among observations of soot oxidation in premixed and diffusion flames clearly must be resolved. Motivated by these findings, the present study undertook measurements of soot and flame properties within the soot oxidation region of some typical laminar diffusion flames and exploited the new measurements to identify soot oxidation mechanisms for these conditions. Present considerations were limited to the early stages of soot oxidation (carbon consumption less than 70%) where reactions at the surface of primary soot particles dominate the process, rather than the later stages when particle porosity and internal particle oxidation become important as discussed by Neoh et al.

El-Leathy, A. M.↗

Radiative Losses in the Chromosphere During A C-Class Flare

Context. Solar flares release an enormous amount of energy (~10(exp 32) erg) into the corona. A substantial fraction of this energy is transported to the lower atmosphere, which results in chromospheric heating. The mechanisms that transport energy to the lower solar atmosphere during a flare are still not fully understood. Aims. We aim to estimate the temporal evolution of the radiative losses in the chromosphere at the footpoints of a C-class flare, in order to set observational constraints on the electron beam parameters of a RADYN flare simulation. Methods. We estimated the radiative losses from hydrogen, and singly ionized Ca and Mg using semiempirical model atmospheres, which were inferred from a multiline inversion of observed Stokes profiles obtained with the CRISP and CHROMIS instruments on the Swedish 1-m Solar Telescope. The radiative losses were computed taking into account the effect of partial redistribution and non-local thermodynamic equilibrium. To estimate the integrated radiative losses in the chromosphere, the net cooling rates were integrated between the temperature minimum and the height where the temperature reaches 10 kK. We also compared our time series of radiative losses with those from the RADYN flare simulations. Results. We obtained a high spatial-resolution map of integrated radiative losses around the flare peak time. The stratification of the net cooling rate suggests that the Ca IR triplet lines are responsible for most of the radiative losses in the flaring atmosphere. During the flare peak time, the contribution from Ca II H and K and Mgii h and k lines are strong and comparable to the Ca IR triplet (~32kW m(exp -2)). Since our flare is a relatively weak event, the chromosphere is not heated above 11 kK, which in turn yields a subdued Lyα contribution (~7kW m(exp -2)) in the selected limits of the chromosphere. The temporal evolution of total integrated radiative losses exhibits sharply rising losses (0.4kW m(exp -2) (s(exp -1)) and a relatively slow decay (0.23kW m(exp -2) s(exp -1)). The maximum value of total radiative losses is reached around the flare peak time and can go up to 175kWm􀀀2 for a single pixel located at footpoint. After a small parameter study, we find the best model-data consistency in terms of the amplitude of radiative losses and the overall atmospheric structure with a RADYN flare simulation in the injected energy flux of 5 × 10(exp 10) erg s(exp -1) cm(exp -2).

Chromosphere↗

The attenuation of gamma-ray emission in strongly-magnetized pulsars

Gamma rays from pulsars can be efficiently attenuated in their magnetospheres via the mechanism of single photon pair production and the exotic quantum electrodynamics (QED) process of photon splitting. The modeling of strongly magnetized gamma ray pulsars focusing on the escape or attenuation of photons emitted near the pole at the neutron star surface in dipole fields in a Schwarzschild metric is considered. It was found that pair production and splitting totally inhibit emission above a value of between 10 and 30 MeV in PSR 1509-58 whose surface field is inferred as being high. The principle predictions of the attenuation analysis are reviewed and the observational diagnostic capabilities of the model are considered. The diagnostics include the energy of the gamma ray turnover and the spectral polarization, which constrain the estimated polar cap size and field strength and can determine the relative strength of splitting and pair creation.

Baring, Matthew G.↗

Constraints on the Opening Rate of Bands on Europa

The opening rates of two bands on Europa, inferred to be sites of spreading of the icy lithosphere, are constrained based on a mid-ocean ridge analog model. Estimates of brittle-ductile transition depth combined with a conductive cooling model limit active band lifetimes to 0.24 - 35 Myr and strain rates of 8.1 x 10(exp -13) - 8.2 x 10(exp -15)/s. These values suggest tensile strengths for ice on Europa of 0.46 - 2.3 MPa, consistent with nonsynchronous rotation as the dominant driving mechanism for band opening.

Stempel, M. M.↗

The heliospheric ambipolar potential inferred from sunward-propagating halo electrons

ABSTRACT We provide evidence that the sunward-propagating half of the solar wind electron halo distribution evolves without scattering in the inner heliosphere. We assume the particles conserve their total energy and magnetic moment, and perform a ‘Liouville mapping’ on electron pitch angle distributions measured by the Parker Solar Probe SPAN-E instrument. Namely, we show that the distributions are consistent with Liouville’s theorem if an appropriate interplanetary potential is chosen. This potential, an outcome of our fitting method, is compared against the radial profiles of proton bulk flow energy. We find that the inferred potential is responsible for nearly 100 per cent of the proton acceleration in the solar wind at heliocentric distances 0.18-0.79 AU. These observations combine to form a coherent physical picture: the same interplanetary potential accounts for the acceleration of the solar wind protons as well as the evolution of the electron halo. In this picture the halo is formed from a sunward-propagating population that originates somewhere in the outer heliosphere by a yet-unknown mechanism.

79 ASTRONOMY AND ASTROPHYSICS↗

Gene targets for engineering osmotolerance in Caldicellulosiruptor bescii

Abstract Background Caldicellulosiruptor bescii , a promising biocatalyst being developed for use in consolidated bioprocessing of lignocellulosic materials to ethanol, grows poorly and has reduced conversion at elevated medium osmolarities. Increasing tolerance to elevated fermentation osmolarities is desired to enable performance necessary of a consolidated bioprocessing (CBP) biocatalyst. Results Two strains of C. bescii showing growth phenotypes in elevated osmolarity conditions were identified. The first strain, ORCB001, carried a deletion of the FapR fatty acid biosynthesis and malonyl-CoA metabolism repressor and had a severe growth defect when grown in high-osmolarity conditions—introduced as the addition of either ethanol, NaCl, glycerol, or glucose to growth media. The second strain, ORCB002, displayed a growth rate over three times higher than its genetic parent when grown in high-osmolarity medium. Unexpectedly, a genetic complement ORCB002 exhibited improved growth, failing to revert the observed phenotype, and suggesting that mutations other than the deleted transcription factor (the fruR/cra gene) are responsible for the growth phenotype observed in ORCB002. Genome resequencing identified several other genomic alterations (three deleted regions, three substitution mutations, one silent mutation, and one frameshift mutation), which may be responsible for the observed increase in osmolarity tolerance in the fruR / cra -deficient strain, including a substitution mutation in dnaK , a gene previously implicated in osmoresistance in bacteria. Differential expression analysis and transcription factor binding site inference indicates that FapR negatively regulates malonyl-CoA and fatty acid biosynthesis, as it does in many other bacteria. FruR/Cra regulates neighboring fructose metabolism genes, as well as other genes in global manner. Conclusions Two systems able to effect tolerance to elevated osmolarities in C. bescii are identified. The first is fatty acid biosynthesis. The other is likely the result of one or more unintended, secondary mutations present in another transcription factor deletion strain. Though the locus/loci and mechanism(s) responsible remain unknown, candidate mutations are identified, including a mutation in the dnaK chaperone coding sequence. These results illustrate both the promise of targeted regulatory manipulation for osmotolerance (in the case of fapR ) and the challenges (in the case of fruR/cra ).

59 BASIC BIOLOGICAL SCIENCES↗

Design and Characterization of a Lens-Coupled System for Dynamic X-Ray Diffraction

X-ray diffraction (XRD) is a necessary technique for understanding states of materials under static and dynamic loading conditions. The higher-pressure Equation of State (EOS) of many materials can only be explored via shock or ramp compression at temperatures and pressures of interest. While static XRD work has yielded EOS measurements in the 100 - 200 GPa regime, dynamic X-ray diffraction (DXRD) can explore EOS phases in the TPa regime, which closely resembles inner-core planetary conditions. DXRD hinges on the ability to measure the exact phase or phase change of a material while under dynamic loading conditions. Macroscopic diagnostic systems (e.g. velocimetry and pyrometry) can infer a phase change but not identify the specific phase entered by a material. While microscopic (atomic-level) diagnostic systems (e.g. DXRD) have been designed and implemented in Department of Energy’s (DOE) National Laboratories complex, the unique nature of Sandia National Laboratories’ Pulsed Power Facility (Z Machine) prohibits the use of such devices. The destructive nature of Z experiments presents a challenge to data capture and retrieval. Furthermore there are electromagnetic interference, X-ray background, and mechanical constraints to consider. Thus, a multi-part X-ray diagnostic for use on the Z Machine and Z-Beamlet Laser system has been designed and analyzed. Portions of this new DYnamic SCintillator Optic (DYSCO) have been built, tested and fielded. A data analysis software has been written. Finally, the radiance profile of the DYSCO’s scintillator has been characterized through experiments performed at the University of Arizona.

36 MATERIALS SCIENCE↗

Predicting Cell Death and Mutation Frequency for a Wide Spectrum of LET by Assuming DNA Break Clustering Inside Repair Domains

The high relative biological effectiveness (RBE) of high charged and energy (HZE) particles for cell death, DNA mutations and cancer remain based on experimental data. In this work, we propose that the existence of DNA repair domains is sufficient to predict both cell death and mutation frequencies for any LET by only taking into account experimental data from low-LET, offering one mechanism for RBE across LET. We hypothesize that whenever multiple DNA double-strand breaks (DSBs) are generated within the same DNA repair domain, DSBs are actively regrouped for more efficient repair [1]. This hypothesis has been supported by the low-LET sublinear dose response observed at doses greater than ~1Gy for 53BP1 radiation-induced foci (RIF) reflecting increasing DSB/RIF with dose [2]. Previously, we modeled radiation-induced cell death of human breast cells by first inferring the size of these domains from the dose dependence of low-LET RIF, and by associating a lethality factor to the number of pairs of DSBs in each RIF [1]. In this work, we first integrate the new NASA computer models RITCARD (Relativistic Ion Tracks, Chromosome Aberrations, Repair, and Damage) [3] and BDSTracks (Biological Damage by Stochastic Tracks) for a more accurate microdosimetry and a better model of the nuclear organization to predict the location of DSBs. A large array of particles and energy are simulated, covering more than three orders of magnitude for LET (~1-1000 keV/µm). Next, we extend our previous model to predict mutation frequencies by assuming that clustered DSBs increase mutation probability, which is formalized by the mutation frequency being linearly dependent on both the number of DSBs and the number of pairs of DSBs inside individual RIF. Linear coefficients are estimated so that simulations predict accurately mutation frequencies observed in Chinese hamster cells exposed to low-LET. Keeping these coefficients unchanged, we then predict mutation frequencies induced by HZE by simulating DSBs and obtain RBEs for mutations and cell death following the expected experimental bell shape for LET dependence. We also observe an orientation effect that needs to be confirmed, showing different RBE depending on the angle of the HZE beam hitting the main axis of the cell.

Plante, Ianik↗

Getting Warmer: IceCube Nears Freeze Out

IceCube has recently detected a diffuse population of high-energy neutrinos arising from the Milky Way. We use this high-significance detection to place new limits on dark matter (DM) annihilation to neutrinos with two complementary approaches. The first method uses the background-subtracted Galactic longitude distribution of shower events to place a conservative bound on the DM annihilation cross section that does not rely on any assumed Galactic cosmic ray emission model; the resulting limits on the velocity-averaged annihilation cross section improve upon existing bounds by factors of a few. The second method uses the template-dependent neutrino energy spectra from the Inner Galaxy, inferred under different Galactic cosmic ray emission models. This complementary approach shows that the inferred Galactic neutrino intensities are already sensitive to DM contributions near the thermal-relic benchmark for a range of TeV-scale DM masses, though this comparison is more model-dependent. Our results demonstrate that measurements of diffuse Galactic neutrino emission can be used as a powerful probe of DM annihilation into neutrinos. Future observations with IceCube-Gen2 and KM3NeT will substantially extend this sensitivity, potentially allowing a decisive test of the thermal freeze-out mechanism with Galactic neutrino observations.

Mukhopadhyay, Mainak [Fermilab; Chicago U., KICP; ↗

Equation of State Optimization and Uncertainty Quantization: Implementation in the LANL EOS Production Code OpenSesame

We detail the approaches of particle swarm optimization and Bayesian inference through Markov chain Monte Carlo for equation of state development. This work includes formulation of the modeling for the equation of state, numeric optimization of the parametric models via particle swarm optimization, and generation of probability distributions of equations of state from Markov chain Monte Carlo.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗