Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical 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 1,189 records · Page 66

FACETS Dynamical Downscaling Simulations over North America by the CAM-MPAS Variable-Resolution Model

Comprehensive assessment of climate datasets created by statistical or dynamical models is important for effectively communicating model projection and associated uncertainty to stakeholders and decision-makers. The Department of Energy FACETS project aims to foster such communication through development of metrics and their demonstration on a hierarchy of downscaled climate datasets to quantify aspects of climate change projections that are credible, particularly for supporting decisions related to the energy-water-land nexus. As a part of this effort, we have produced a regional climate dataset using the Model for Prediction Across Scales coupled to the Community Atmosphere Model (CAM-MPAS). This global modeling framework is configured with variable-resolution meshes featuring higher resolutions over North America, as well as quasi-uniform resolution meshes across the globe. The variable-resolution configurations allow fine-scale features to be better resolved inside the refinement and interact with global-scale circulations. The dataset includes multiple uniform- (240km and 120km) and variable-resolution (200-50km, 100-25km, and 46-12km) simulations that are designed to be compatible with other regional climate simulations that contribute to the hierarchy of downscaled climate datasets of the project. Furthermore, the dataset consists of simulations for both the present-day (1989-2010) and future (2079-2100) climate and post-processing of the model output has been coordinated across the project for consistency to facilitate common analysis across the hierarchy of datasets. Altogether, this CAM-MPAS model dataset provides a unique opportunity to assess the influence of resolutions and modeling framework on model credibility and climate change projection.

Sakaguchi, Koichi↗

$\Sigma (1385)^{\pm }$ resonance production in Pb–Pb collisions at $\sqrt{s_{\textrm{NN}}} = 5.02$ TeV

Hadronic resonances are used to probe the hadron gas produced in the late stage of heavy-ion collisions since they decay on the same timescale, of the order of 1–10 fm/c, as the decoupling time of the system. In the hadron gas, (pseudo)elastic scatterings among the products of resonances that decayed before the kinetic freeze-out and regeneration processes counteract each other, the net effect depending on the resonance lifetime, the duration of the hadronic phase, and the hadronic cross sections at play. In this context, the Σ(1385) ± particle is of particular interest as models predict that regeneration dominates over rescattering despite its relatively short lifetime of about 5.5 fm/c. The first measurement of the Σ(1385) ± resonance production at midrapidity in Pb–Pb collisions at $\sqrt{s_{NN}}$=5.02 TeV with the ALICE detector is presented in this Letter. The resonances are reconstructed via their hadronic decay channel, Λπ, as a function of the transverse momentum (p T ) and the collision centrality. The results are discussed in comparison with the measured yield of pions and with expectations from the statistical hadronization model as well as commonly employed event generators, including PYTHIA8/Angantyr and EPOS3 coupled to the UrQMD hadronic cascade afterburner. None of the models can describe the data. For Σ(1385) ± , a similar behaviour as K*(892) 0 is observed in data unlike the predictions of EPOS3 with afterburner.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Theory of Alike Selectivity in Biological Channels

We introduce a statistical mechanical model of the selectivity filter that accounts for the interaction between ions within the channel and derive Eisenman equation of the filter selectivity directly from the condition of barrier-less conduction.

alike selectivity↗

Evaluation of the Maximum Soluble Water Recovery from Pretreated Urine Brines in Low Earth Orbit

The Urine Processor Assembly (UPA) on the International Space Station recovers water from pretreated urine by thin-film, vacuum compression distillation (VCD). A number of iterative improvements to the distiller assembly (DA) have enabled increases in on-orbit water recovery and improved distillate water quality at on-orbit recoveries up to 87% by volume relative to the volume of pretreated urine. The phosphoric acid-based pretreatment formulation has enabled these higher recoveries on-orbit since May of 2016. Application of simple empirical chemistry models have allowed water recovery of sulfuric-acid pretreated urine from the Russian Orbiting segment to be distilled with phosphoric acid pretreated urine from the U.S. Orbiting Segment. Returned liquid brine from brine filters in the UPA brine loop have enabled empirical measurement of gypsum solubility as a function of recovery, conductivity, or mass density in the real brines of ISS. Calculations of maximum soluble recovery, MSR, are based on preventing the precipitation of gypsum, which is the first mineral to precipitate for both types of pretreated urine. A statistical urine model of urine composition provides an estimate of worst-case urine from a crew of three during one year of operations on ISS. The MSR is defined based on the worst-case concentrations of urinary calcium and urinary sulfate that are pretreated, pooled, and distilled. The potential to increase water recoveries above 87% without gypsum precipitation are evaluated. The potential benefits and costs of the addition of a conductivity sensor or a Coriolis mass flow meter to the brine loop are evaluated.

Dean L Muirhead↗

Constraining the hadronic spectrum and repulsive interactions in a hadron resonance gas via fluctuations of conserved charges

We simultaneously incorporate two common extensions of the hadron resonance gas model, namely the addition of extra, unconfirmed resonances to the particle list and the excluded volume repulsive interactions. We emphasize the complementary nature of these two extensions and identify combinations of conserved charge susceptibilities that allow to constrain them separately. In particular, ratios of second-order susceptibilities like $\chi_{11}^{BQ}/\chi_2^B$ and $\chi_{11}^{BS}/\chi_2^B$ are sensitive only to the baryon spectrum, while fourth-to-second order ratios like $\chi_4^B/\chi_2^B$, $\chi_{31}^{BS}/\chi_{11}^{BS}$, or $\chi_{31}^{BQ}/\chi_{11}^{BQ}$ are mainly determined by repulsive interactions. Analysis of the available lattice results suggests the presence of both the extra states in the baryon-strangeness sector and the repulsive baryonic interaction, with indications that hyperons have a smaller repulsive core than non-strange baryons. Furthermore, the modified hadron resonance gas model presented here significantly improves the description of lattice QCD susceptibilities at chemical freeze-out and can be used for the analysis of event-by-event fluctuations in heavy-ion collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Towards Inner Magnetosphere Particle and Field Models

We report on the initial stages of an effort to construct comprehensive empirical models of the plasma and fields in the inner magnetosphere. The models are based not only on archival data (e.g., SCATHA) but also on current spacecraft mission data (e.g., POLAR). In this component of the effort, we incorporate the recently archived data provided by the SCATRA (Spacecraft Charging AT High Altitude) satellite. The SCAT14A satellite was in a near-geostationary orbit and was in operation for more than a decade. In this paper, we focus on the SCATRA plasma and magnetic field data from approximately the first two years of operation. The time-series data are binned according to spatial location and geomagnetic activity. Examples of statistical, empirical models from this initial effort are presented; even the simply-constructed preliminary models reveal such subtle features as the previously reported warm (greater than 100 eV) plasma density enhancement near the magnetic equator.

Spence, H. E.↗

A statistical equilibrium analysis of interstellar CO toward zeta Ophiuchi as recorded by the Goddard High Resolution Spectrograph

High-resolution and high signal-to-noise ratio echelle observations made by the Goddard High Resolution Spectrograph (GHRS), aboard the Hubble Space Telescope, of molecular (12)C(16)O and (13)C(16)O in the absorbing gas toward zeta Oph are reported. A statistical equilibrium model, incorporating a different way of dealing with the radiation field in the interstellar clouds, is presented. Based upon previous observations, two velocity components of the gas at V(sub LSR) = - 0.79 and 0.39 km/s, corresponding to solar velocity = -14.6 and -13.4 km/s, together with three different sets of velocity dispersions are adopted in our models. The ratios of the computed emergent intensities (I(sub em)) resulting from the emission lines in the lowest transitions (J = 1 goes to 0, 2 goes to 1, 3 goes to 2) of ground state of CO are presented. The ratios for the two velocity components, (I(sub em)(1 goes to 0)(sub 1))/(I(sub em)(1 goes to 0)(sub 2)), and (I(sub em)(2 goes to 1)(sub 1))/I(sub em)(1 goes to 0)(sub 1)) are comparable to the observed corrected antenna temperature ratios. From the solutions to the statistical equilibrium equations including the six lowest rotational transitions of CO and from multiple Voigt profile fits, the observed (12)CO and (13)CO UV abroption features are successfully reproduced. Profile fitting and double cloud curve-of-growth analysis using the GHRS and Copernicus data indicate that the abundance of CO is (1.8 +/- 0.2) x 10(exp 15)/sq cm. The most recently measured wavelength of the (6 goes to 0) (13)CO band by Haridass & Huber (R(0) = 1370.616 A) shows that the tentative (13)CO features in the GHRS data are unquestionably identified. If there is no blending from any other unknown interstellar species, the isotope ratio of (12)CO/(13)CO is approximately 82 +/- 25. This result is close to the terrestrial value of 90 and implies that the CO is not substantially fractionated, and (12)C(16)O experiences little preferential self-shielding through the selective isotopic photodissociation process. This conclusion differs strongly with that of Sheffer et al. Other implications of the results are discussed.

Lyu, Cheng-Hsuan↗

Sensitivity of Simulated Global Ocean Carbon Flux Estimates to Forcing by Reanalysis Products

Reanalysis products from MERRA, NCEP2, NCEP1, and ECMWF were used to force an established ocean biogeochemical model to estimate air-sea carbon fluxes (FCO2) and partial pressure of carbon dioxide (pCO2) in the global oceans. Global air-sea carbon fluxes and pCO2 were relatively insensitive to the choice of forcing reanalysis. All global FCO2 estimates from the model forced by the four different reanalyses were within 20% of in situ estimates (MERRA and NCEP1 were within 7%), and all models exhibited statistically significant positive correlations with in situ estimates across the 12 major oceanographic basins. Global pCO2 estimates were within 1% of in situ estimates with ECMWF being the outlier at 0.6%. Basin correlations were similar to FCO2. There were, however, substantial departures among basin estimates from the different reanalysis forcings. The high latitudes and tropics had the largest ranges in estimated fluxes among the reanalyses. Regional pCO2 differences among the reanalysis forcings were muted relative to the FCO2 results. No individual reanalysis was uniformly better or worse in the major oceanographic basins. The results provide information on the characterization of uncertainty in ocean carbon models due to choice of reanalysis forcing.

MERRA↗

Constraining the destruction rate of K 40 in stellar nucleosynthesis through the study of the Ar 40 ( p , n ) K 40 reaction

Background: We present that K 40 plays a significant role in the radiogenic heating of Earth-like exoplanets, which can affect the development of a habitable environment on their surfaces. The initial amount of K 40 in the interior of these planets depends on the composition of the interstellar clouds from which they formed. Within this context, nuclear reactions that regulate the production of K 40 during stellar evolution can play a critical role. Purpose: In this study, we constrain for the first time the astrophysical reaction rate of K 40 ( n , p ) Ar 40 , which is responsible for the destruction of K 40 during stellar nucleosynthesis. We provide to the nuclear physics community high-resolution data on the cross section and angular distribution of the Ar 40 ( p , n ) K 40 reaction. These are important to various applications involving Ar 40 . The associated reaction rate of the Ar 40 ( p , n ) K 40 process addresses a reaction rate gap in the Joint Institute for Nuclear Astrophysics REACLIB database in the region of intermediate-mass isotopes. Methods: We performed differential cross-section measurements on the Ar 40 ( p , n ) K 40 reaction, for six energies in the center-of-mass system between 3.2 and 4.0 MeV and various angles between 0 ° and 135 ° . The experiment took place at the Edwards Accelerator Laboratory at Ohio University using the beam swinger target location and a standard neutron time-of-flight technique. We extracted total and partial cross sections by integrating the double differential cross sections we measured. Results: The total and partial cross sections varied with energy due to the contribution from isobaric analog states and Ericson type fluctuations. The energy-averaged neutron angular distributions were symmetrical relative to 90 ° . Based on the experimental data, local transmission coefficients were extracted and were used to calculate the astrophysical reaction rates of Ar 40 ( p , n ) K 40 and K 40 ( n , p ) Ar 40 reactions. The new rates were found to vary significantly from the theoretical rates in the REACLIB library. We implemented the new rates in network calculations to study nucleosynthesis via the slow neutron capture process, and we found that the produced abundance of K 40 is reduced by up to 10% compared to calculations with the library rates. At the same time, the above result removes a significant portion of the previous theoretical uncertainty on the K 40 yields from stellar evolution calculations. Conclusions: Our results support a destruction rate of K 40 in massive stars via the K 40 ( n , p ) Ar 40 reaction that is larger compared to previous estimates. The rate of K 40 destruction via the K 40 ( n , p ) Ar 40 reaction now has a dramatically reduced uncertainty based on our measurement. Lastly, this result directly affects the predicted stellar yields of K 40 from nucleosynthesis, which is a critical input parameter for the galactic chemical evolution models that are currently employed for the study of significant properties of exoplanets.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Machine learning prediction of electron density and temperature from He I line ratios

We propose to utilize machine learning to predict the electron density, ne, and temperature, T e , from He I line intensity ratios. In this approach, training data consist of measured He I line ratios as input and ne and T e measured using other diagnostic(s) as desired output, which is a Langmuir probe in our study. Support vector machine regression analysis is, then, performed with the training data to develop a predictive model for n e and T e , separately. It is confirmed that n e and T e predicted using the developed models agree well with those from the Langmuir probe in the ranges of 0.28 × 10 18 ≤ n e (m -3 ) ≤ 3.8 × 10 18 and 3.2 ≤ T e (eV) ≤ 7.5. The developed models are, further, examined with an evaluation data, which are not included in the training data, and are found to well reproduce absolute values and radial profiles of probe-measured n e and T e .

47 OTHER INSTRUMENTATION↗

Generation of Fragment Angular Momentum in Fission

A recent analysis of experimental data [J. Wilson et al., Nature (London) 590, 566 (2021)] found that the angular momenta of nuclear fission fragments are uncorrelated. Based on this finding, the authors concluded that the spins are therefore determined only after scission has occurred. We show here that the nucleon-exchange mechanism, as implemented in the well-established event-by-event fission model freya, while agitating collective rotational modes in which the two spins are highly correlated, nevertheless leads to fragment spins that are largely uncorrelated. We find that this counterexample invalidates the conclusion in [J. Wilson et al.] that uncorrelated spins must necessarily have been generated after scission (a potentious conclusion that would rule out all models that generate the fragment spins prior to scission). Furthermore, it was reported [J. Wilson et al.] that the mass dependence of the average fragment spin has a sawtooth structure. We demonstrate that such a behavior naturally emerges when shell and deformation effects are included in the moments of inertia of the fragments at scission.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Modified Brink-Axel hypothesis for astrophysical Gamow-Teller transitions

Weak interaction charged current transition strengths from highly excited nuclear states are fundamental ingredients for accurate modeling of compact object composition and dynamics, but are difficult to obtain either from experiment or theory. For lack of alternatives, calculations have often fallen back upon a generalized Brink-Axel hypothesis, that is, assuming the strength function (transition probability) is independent of the initial nuclear state but depends only upon the transition energy and the weak interaction properties of the parent nucleus ground state. Here we present numerical evidence for a modified “local” Brink-Axel hypothesis for Gamow-Teller transitions for pf-shell nuclei relevant to astrophysical applications. Specifically, while the original Brink-Axel hypothesis does not hold globally, strength functions from initial states nearby in energy are similar within statistical fluctuations. This agrees with previous work on strength function moments. Using this modified hypothesis, we can tackle strength functions at previously intractable initial energies, using semiconverged initial states at arbitrary excitation energy. Our work provides a well-founded method for computing accurate thermal weak transition rates for medium-mass nuclei at temperatures occurring in stellar cores near collapse. We finish by comparing results to previous calculations of astrophysical rates.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Dynamic density functional theory of polymers with salt in electric fields

Here we present a dynamic density functional theory for modeling the effects of applied electric fields on the local structure of polymers with added salt (polymer electrolytes). Time-dependent equations for the local electrostatic potential and volume fractions of polymer, cation, and anion of added salt are developed using the principles of linear irreversible thermodynamics. For such a development, a field theoretic description of the free energy of polymer melts doped with salts is used, which captures the effects of local variations in the dielectric function. Connections of the dynamic density functional theory with experiments are established by relating the three phenomenological Onsager’s transport coefficients of the theory to the mutual diffusion of electrolyte, ionic conductivity, and transference number of one of the ions. The theory is connected with a statistical mechanical model developed by Bearman and Kirkwood [J. Chem. Phys. 28, 136 (1958)] after relating the three transport coefficients to friction coefficients. The steady-state limit of the dynamic density functional theory is used to understand the effects of dielectric inhomogeneity on the phase separation in polymer electrolytes. The theory developed here provides not only a way to connect with experiments but also to develop multi-scale models for studying connections between local structure and ion transport in polymer electrolytes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Remote sensing of earth terrain

Two monographs and 85 journal and conference papers on remote sensing of earth terrain have been published, sponsored by NASA Contract NAG5-270. A multivariate K-distribution is proposed to model the statistics of fully polarimetric data from earth terrain with polarizations HH, HV, VH, and VV. In this approach, correlated polarizations of radar signals, as characterized by a covariance matrix, are treated as the sum of N n-dimensional random vectors; N obeys the negative binomial distribution with a parameter alpha and mean bar N. Subsequently, and n-dimensional K-distribution, with either zero or non-zero mean, is developed in the limit of infinite bar N or illuminated area. The probability density function (PDF) of the K-distributed vector normalized by its Euclidean norm is independent of the parameter alpha and is the same as that derived from a zero-mean Gaussian-distributed random vector. The above model is well supported by experimental data provided by MIT Lincoln Laboratory and the Jet Propulsion Laboratory in the form of polarimetric measurements.

Kong, J. A.↗

Modeling Buried Object Brightness and Visibility for Ground Penetrating Radar

Comparing the observed brightness of various buried objects is a straightforward way to characterize the performance of a ground penetrating radar (GPR) system. However, a limitation arises. A simple comparison of buried object brightness values does not disentangle the effects of the GPR system itself from the system's operating environment and the objects being observed. Therefore, with brightness values exhibiting an unknown synthesis of systemic, environmental, and object factors, GPR system analysis becomes a convoluted affair. In this work, we use an experimentally collected dataset of over 25,000 object observations from five different multistatic radar arrays to develop models of buried object brightness and control for these various effects. Additionally, our modeling efforts provide a means for quantifying the relative brightness of GPR systems, the objects they detect, and the physical properties of those objects that influence observed brightness. To evaluate the models' performance on new object observations, we repeatedly simulate fitting them to half the dataset and predicting the observed brightness values of the unseen half. In addition, we introduce a method for estimating the probability that individual observations constitute a visible object, which aids in failure analysis, performance characterization, and dataset cleaning.

97 MATHEMATICS AND COMPUTING↗

A stochastic damage model for the rupture prediction of a multi-phase solid. I - Parametric studies. II - Statistical approach

A stochastic damage model for predicting the rupture of a brittle multiphase material is developed, based on the microcrack-macrocrack interaction. The model, which incorporates uncertainties in locations, orientations, and numbers of microcracks, characterizes damage by microcracking and fracture by macrocracking. A parametric study is carried out to investigate the change of the stress intensity at the macrocrack tip by the configuration of microcracks. The inherent statistical distribution of the fracture toughness arising from the intrinsic random nature of microcracks is explored using a statistical approach. For this purpose, a computer simulation model is introduced, which incorporates a statistical characterization of geometrical parameters of a random microcrack array.

Lua, Yuan J.↗

Performance of Two Cloud-Radiation Parameterization Schemes in the Finite Volume General Circulation Model for Anomalously Wet May and June 2003 Over the Continental United States and Amazonia

An objective assessment of the impact of a new cloud scheme, called Microphysics of Clouds with Relaxed Arakawa-Schubert Scheme (McRAS) (together with its radiation modules), on the finite volume general circulation model (fvGCM) was made with a set of ensemble forecasts that invoke performance evaluation over both weather and climate timescales. The performance of McRAS (and its radiation modules) was compared with that of the National Center for Atmospheric Research Community Climate Model (NCAR CCM3) cloud scheme (with its NCAR physics radiation). We specifically chose the boreal summer months of May and June 2003, which were characterized by an anomalously wet eastern half of the continental United States as well as northern regions of Amazonia. The evaluation employed an ensemble of 70 daily 10-day forecasts covering the 61 days of the study period. Each forecast was started from the analyzed initial state of the atmosphere and spun-up soil moisture from the first-day forecasts with the model. Monthly statistics of these forecasts with up to 10-day lead time provided a robust estimate of the behavior of the simulated monthly rainfall anomalies. Patterns of simulated versus observed rainfall, 500-hPa heights, and top-of-the-atmosphere net radiation were recast into regional anomaly correlations. The correlations were compared among the simulations with each of the schemes. The results show that fvGCM with McRAS and its radiation package performed discernibly better than the original fvGCM with CCM3 cloud physics plus its radiation package. The McRAS cloud scheme also showed a reasonably positive response to the observed sea surface temperature on mean monthly rainfall fields at different time leads. This analysis represents a method for helpful systematic evaluation prior to selection of a new scheme in a global model.

Sud, Y. C.↗

Full Counting Statistics of Charge in Chaotic Many-Body Quantum Systems

We investigate the full counting statistics of charge transport in U(1)-symmetric random unitary circuits. We consider an initial mixed state prepared with a chemical potential imbalance between the left and right halves of the system and study the fluctuations of the charge transferred across the central bond in typical circuits. Using an effective replica statistical mechanics model and a mapping onto an emergent classical stochastic process valid at large on-site Hilbert space dimension, we show that charge transfer fluctuations approach those of the symmetric exclusion process at long times, with subleading t –1/2 quantum corrections. Here, we discuss our results in the context of fluctuating hydrodynamics and macroscopic fluctuation theory of classical nonequilibrium systems and check our predictions against direct matrix-product state calculations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗