Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model fitting”

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 253 records · Page 14

Blowout Jets: Hinode X-Ray Jets that Don't Fit the Standard Model

Nearly half of all H-alpha macrospicules in polar coronal holes appear to be miniature filament eruptions. This suggests that there is a large class of X-ray jets in which the jet-base magnetic arcade undergoes a blowout eruption as in a CME, instead of remaining static as in most solar X-ray jets, the standard jets that fit the model advocated by Shibata. Along with a cartoon depicting the standard model, we present a cartoon depicting the signatures expected of blowout jets in coronal X-ray images. From Hinode/XRT movies and STEREO/EUVI snapshots in polar coronal holes, we present examples of (1) X-ray jets that fit the standard model, and (2) X-ray jets that do not fit the standard model but do have features appropriate for blowout jets. These features are (1) a flare arcade inside the jet-base arcade in addition to the small flare arcade (bright point) outside that standard jets have, (2) a filament of cool (T is approximately 80,000K) plasma that erupts from the core of the jetbase arcade, and (3) an extra jet strand that should not be made by the reconnection for standard jets but could be made by reconnection between the ambient unipolar open field and the opposite-polarity leg of the filament-carrying flux-rope core field of the erupting jet-base arcade. We therefore infer that these non-standard jets are blowout jets, jets made by miniature versions of the sheared-core-arcade eruptions that make CMEs

Moore, Ronald L.↗

Measurement uncertainty of rossi-alpha neutron experiments

Rossi-alpha neutron experiments are used to estimate the prompt neutron decay constant of a fissile assembly, a quantity of widespread interest in applications including in nuclear nonproliferation and criticality safety. This work develops a mathematical model to efficiently estimate measurement uncertainty of Rossi-alpha neutron experiments inferred from a two-exponential fit model with histogram binning. The derived uncertainty estimates were validated using repeated Rossi-alpha measurements of a subcritical, 4.5-kg sphere of weapons-grade, alpha-phase plutonium with nickel, copper, tungsten, and polyethylene reflectors. The estimates of uncertainty for the histogram data produced by the model were conservative and agree with the reference uncertainties within noise. The estimates of the prompt neutron decay constant uncertainty agreed with the reference uncertainties within one standard deviation. In conclusion, the proposed model will reduce total measurement times, ultimately reducing operational and procedural costs in application.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A fusion relevant data-driven engineering void swelling model for 9Cr tempered martensitic steels

The UCSB database on cavity evolution in 9-12Cr tempered martensitic steels (TMS), includes the results for both dual heavy and helium ion (DII), and High Flux Isotope Reactor (HFIR) in situ helium injection (ISHI) neutron irradiations at 500°C. These results were combined with literature single ion and fission neutron irradiation data to derive a model for the void volume fraction, f v , as a function of displacements per atom (dpa) and transmutant helium concentrations in atomic parts per million (appm). The scientific foundation for the paper is described in a companion paper entitled “Cavity Evolution and Void Swelling in Dual Ion Irradiated Tempered Martensitic Steels”. Here, in this study, we show that f v (dpa, He/dpa) is described by the incubation dose, dpa i , for the onset of void growth, and the post-incubation growth rate, f v ’(%/dpa). Both dpa i and f v ’ decrease with increasing He/dpa at > ~ 5. The dpa i is also lower for the ISHI neutron irradiations at the same He/dpa. Single heavy ion and fission reactor neutron irradiations, with low He/dpa ratios, have a much larger dpa i . Based on a combined analysis of DII, single ion, ISHI and fission neutron data, we further show that the post-incubation f v data analyzed here have a common empirical curve shape, with f v ’ reaching up to ~ 0.2%/dpa at very high dpa. We also show that f v ’ can be predicted based on a physical model of defect partitioning between evolving sinks. At 500°C and fusion relevant He/dpa ≈ 10, the best-fit model predicts nominal swelling, S = f v /(1-f v ), of ~ 1.1, 4.9 and 16% at 50, 100 and 200 dpa, respectively. The physically motivated, data-driven model includes estimated uncertainties for both dpa i and f v ’.

36 MATERIALS SCIENCE↗

Isotherm Modeling and Techno-Economic Analysis of Contactor Technologies for New Tetraamine-Appended MOF for NGCC Applications

Recently, a family of tetraamine-functionalized metal-organic frameworks (MOF) has been reported as promising sorbent materials for capturing CO2 from flue gas conditions relevant to natural gas combined cycle (NGCC) applications. The main advantages of these materials are their two-step cooperative CO2 adsorption, which gives rise to unusual two step-shaped CO2 adsorption profiles and their high thermal stability. This work presents the modelling of the two-transition isotherm of the tetraamine-appended MOF, N,N'-bis(3-aminopropyl)-1,4-diaminobutane (3-4-3)-appended Mg2(dobpdc), and the techno-economic analysis (TEA) of carbon capture processes utilizing this sorbent. Due to the unusual isotherm shapes of the experimental CO2 adsorption data for tetraamine-appended Mg2(dobpdc) and the strong nonlinearity of CO2 loading with respect to temperature and pressure, we tested two different models which use logistic functions for representing the different isotherm behaviors in the different pressure ranges. The first model uses the quadratic isotherm model in the low-pressure region, the Langmuir isotherm model in the middle pressure range, and the dual site Langmuir isotherm model in the high-pressure range. To model the transition between regions we used the arctangent functions independent of temperature and the thermal effect was accounted by using the Clausius-Clapeyron relation]. The second model is an extension of the weighted dual-site Langmuir isotherm model presented by Hughes et al.. In this extended model, the dual-site Langmuir isotherm is employed in the three transition regions, using temperature-dependent logistic functions to activate or deactivate the isotherm model in the low, middle, and high-pressure ranges. Both models fit the experimental data quite well with root mean squared errors (RMSE’s) of 0.41 and 0.17 for model 1 and model 2, respectively. Since model 2 resulted in a lower RMSE, it was leveraged for the development of the gas/solid contactor models used by the TEA. Specifically, two different contactor models, an axial-flow fixed bed and moving bed contactor, were developed as part of this work. These models are dynamic, pressure-driven, and consist of mass, energy, and momentum conservation equations. A kinetic model was also developed by performing parameter estimation using experimental fixed bed breakthrough data. These models are then used to simulate CO2 capture processes from the flue gas generated from a ~600 gross MW NGCC power plant. A cost model was developed which considers the capital cost of the reactors and the significant operating costs such as steam and electricity. Using NETL’s Framework for the Optimization and Quantification of Uncertainty of Uncertainty and Surrogates tool (FOQUS), which has the capability of linking models built using numerous modelling platforms with derivative-free optimization solvers, a techno-economic optimization of the carbon capture processes was performed which minimizes the cost of capture.

Caballero, Daison↗

On the Response of the Special Sensor Microwave/Imager to the Marine Environment: Implications for Atmospheric Parameter Retrievals

A reasonably rigorous basis for understanding and extracting the physical information content of Special Sensor Microwave/Imager (SSM/I) satellite images of the marine environment is provided. To this end, a comprehensive algebraic parameterization is developed for the response of the SSM/I to a set of nine atmospheric and ocean surface parameters. The brightness temperature model includes a closed-form approximation to microwave radiative transfer in a non-scattering atmosphere and fitted models for surface emission and scattering based on geometric optics calculations for the roughened sea surface. The combined model is empirically tuned using suitable sets of SSM/I data and coincident surface observations. The brightness temperature model is then used to examine the sensitivity of the SSM/I to realistic variations in the scene being observed and to evaluate the theoretical maximum precision of global SSM/I retrievals of integrated water vapor, integrated cloud liquid water, and surface wind speed. A general minimum-variance method for optimally retrieving geophysical parameters from multichannel brightness temperature measurements is outlined, and several global statistical constraints of the type required by this method are computed. Finally, a unified set of efficient statistical and semi-physical algorithms is presented for obtaining fields of surface wind speed, integrated water vapor, cloud liquid water, and precipitation from SSM/I brightness temperature data. Features include: a semi-physical method for retrieving integrated cloud liquid water at 15 km resolution and with rms errors as small as approximately 0.02 kg/sq m; a 3-channel statistical algorithm for integrated water vapor which was constructed so as to have improved linear response to water vapor and reduced sensitivity to precipitation; and two complementary indices of precipitation activity (based on 37 GHz attenuation and 85 GHz scattering, respectively), each of which are relatively insensitive to variations in other environmental parameters.

Petty, Grant W.↗

Transport and Growth Kinetics in Microgravity Protein Crystal Growth

The dynamic coupling between mass transport and incorporation of growth units into the surface of a crystal growing from solution in microgravity is used to derive quantitative information on the crystal growth kinetics. To this end, new procedures for experiment preparation, interferometric data processing and model fitting have been developed. The use of experimental data from the bulk diffusive maw transport together with a model for steady state stagnant crystal growth allows the detailed quantitative understanding of the kinetics of both the concentration depletion zone around the crystal and the growth of the crystal interface. The protein crystal used in the experiment is shown to be growing in the mixed kinetic regime (0.2 x 10(exp -6) centimeters per second less than beta R/D less than 0.9 x 10(exp -6) centimeters per second).

Otalora, F.↗

Coulomb Stress Accumulation along the San Andreas Fault System

Stress accumulation rates along the primary segments of the San Andreas Fault system are computed using a three-dimensional (3-D) elastic half-space model with realistic fault geometry. The model is developed in the Fourier domain by solving for the response of an elastic half-space due to a point vector body force and analytically integrating the force from a locking depth to infinite depth. This approach is then applied to the San Andreas Fault system using published slip rates along 18 major fault strands of the fault zone. GPS-derived horizontal velocity measurements spanning the entire 1700 x 200 km region are then used to solve for apparent locking depth along each primary fault segment. This simple model fits remarkably well (2.43 mm/yr RMS misfit), although some discrepancies occur in the Eastern California Shear Zone. The model also predicts vertical uplift and subsidence rates that are in agreement with independent geologic and geodetic estimates. In addition, shear and normal stresses along the major fault strands are used to compute Coulomb stress accumulation rate. As a result, we find earthquake recurrence intervals along the San Andreas Fault system to be inversely proportional to Coulomb stress accumulation rate, in agreement with typical coseismic stress drops of 1 - 10 MPa. This 3-D deformation model can ultimately be extended to include both time-dependent forcing and viscoelastic response.

Smith, Bridget↗

Regression Analysis of Long-term Profile Ozone Data Set from BUV Instruments

We have produced a profile merged ozone data set (MOD) based on the SBUV/SBUV2 series of nadir-viewing satellite backscatter instruments, covering the period from November 1978 - December 2003. In 2004, data from the Nimbus 7 SBUV and NOAA 9,11, and 16 SBUV/2 instruments were reprocessed using the Version 8 (V8) algorithm and most recent calibrations. More recently, data from the Nimbus 4 BUV instrument, which operated from 1970 - 1977, were also reprocessed using the V8 algorithm. As part of the V8 profile calibration, the Nimbus 7 and NOAA 9 (1993-1997 only) instrument calibrations have been adjusted to match the NOAA 11 calibration, which was established from comparisons with SSBUV shuttle flight data. Given the level of agreement between the data sets, we simply average the ozone values during periods of instrument overlap to produce the MOD profile data set. We use statistical time-series analysis of the MOD profile data set (1978-2003) to estimate the change in profile ozone due to changing stratospheric chlorine levels. The Nimbus 4 BUV data offer an opportunity to test the physical properties of our statistical model. We extrapolate our statistical model fit backwards in time and compare to the Nimbus 4 data. We compare the statistics of the residuals from the fit for the Nimbus 4 period to those obtained from the 1978-2003 period over which the statistical model coefficients were estimated.

Frith, Stacey↗

Accelerating Markov Chain Monte Carlo sampling with diffusion models

Global fits of physics models require efficient methods for exploring high-dimensional and/or multimodal posterior functions. We introduce a novel method for accelerating Markov Chain Monte Carlo (MCMC) sampling by pairing a Metropolis-Hastings algorithm with a diffusion model that can draw global samples with the aim of approximating the posterior. We briefly review diffusion models in the context of image synthesis before providing a streamlined diffusion model tailored towards low-dimensional data arrays. We then present our adapted Metropolis-Hastings algorithm which combines local proposals with global proposals taken from a diffusion model that is regularly trained on the samples produced during the MCMC run. Our approach leads to a significant reduction in the number of likelihood evaluations required to obtain an accurate representation of the Bayesian posterior across several analytic functions, as well as for a physical example based on a global fit of parton distribution functions. Our method is extensible to other MCMC techniques, and we briefly compare our method to similar approaches based on normalising flows. A code implementation can be found at https://github.com/NickHunt-Smith/MCMC-diffusion.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Constraints on the origin and interior structure of the major planets

Fitting models to the external gravity field of the major planets (Uranus, Neptune, Jupiter, and Saturn), it is found that certain interior characteristics may be common to all four. For Uranus and Neptune, a model with a central iron-silicate core of approximately four earth masses, an 'ice' layer of H2O, CH4, and NH3 in solar proportions of about 10 earth masses, and an H2-He atmosphere of about 1-2 earth masses gives a good fit to available constraints, including heat flow measurements. Models of Jupiter and Saturn have cores very similar to those of Uranus and Neptune; the H2-He layer, however, is markedly more extensive. Modes of origin that are consistent with these features are discussed. Models of this type predict a considerable enrichment of deuterium relative to primordial solar abundances in Uranus and Neptune. Such enrichment is not observed in Uranus; implications are discussed for interior structure and origin.

Hubbard, W. B.↗

A line-blanketed model stellar atmosphere of Sirius

The method of artificial absorption edges is applied to a grid of model atmospheres to include the effects of hydrogen and metal line blanketing on the atmospheric structure. A best-fit model of Sirius is obtained which is free of the discrepancies inherent in unblanketed or partially blanketed models. All major spectral features calculated from the model agree with the corresponding observed features to within the observational uncertainty. Some additional properties of the models are discussed.

Fowler, J. W.↗

IUE spectra and optical imaging of the oxygen-rich supernova remnant N132D

We present new optical Charge Coupled Devices (CCD) interference filter imagery and International Ultraviolet Explorer (IUE) spectroscopy for the oxygen-rich supernova remnant N132D in the Large Magellanic Cloud. The optical images show a wealth of structure, and comparison with an archival Einstein High Resolution Imager (HRI) X-ray image shows that a few optical features have X-ray counter-parts, but in general there is little correlation between X-ray and optical features. The IUE spectra at two positions show strong lines of carbon and oxygen, with lines of neon, magnesium, silicon, and helium also present and variable in relative intensities. We use optical data for N132D from Dopita & Tuohy (1984) with our UV observations to compare with shock models (both with and without thermal conduction) and X-ray photoionization model calculations. While none of the model fits is entirely satisfactory, the generally weak UV emission relative to optical disagrees with the general character of shock model predictions and indicates that photoionization is the dominant excitation mechanism for the UV/optical emission. This conclusion is similar to what was found for E0102 - 7219, the oxygen-rich remnant in the Small Magellanic Cloud. We derive rough abundances for the emitting material in N132D, compare to stellar nucleosynthesis models, and discuss the implications for its precursor. A precursor near 20 solar mass is consistent with the data.

Blair, William P.↗

SEE Rate Estimation: Model Complexity and Data Requirements

Statistical Methods outlined in [Ladbury, TNS20071 can be generalized for Monte Carlo Rate Calculation Methods Two Monte Carlo Approaches: a) Rate based on vendor-supplied (or reverse-engineered) model SEE testing and statistical analysis performed to validate model; b) Rate calculated based on model fit to SEE data Statistical analysis very similar to case for CREME96. Information Theory allows simultaneous consideration of multiple models with different complexities: a) Model with lowest AIC usually has greatest predictive power; b) Model averaging using AIC weights may give better performance if several models have similar good performance; and c) Rates can be bounded for a given confidence level over multiple models, as well as over the parameter space of a model.

Ladbury, Ray↗

Age-specific case data reveal varying dengue transmission intensity in US states and territories

Dengue viruses (DENV) are endemic in the US territories of Puerto Rico, American Samoa, and the US Virgin Islands, with focal outbreaks also reported in the states of Florida and Hawaii. However, little is known about the intensity of dengue virus transmission over time and how dengue viruses have shaped the level of immunity in these populations, despite the importance of understanding how and why levels of immunity against dengue may change over time. These changes need to be considered when responding to future outbreaks and enacting dengue management strategies, such as guiding vaccine deployment. We used catalytic models fitted to case surveillance data stratified by age from the ArboNET national arboviral surveillance system to reconstruct the history of recent dengue virus transmission in Puerto Rico, American Samoa, US Virgin Islands, Florida, Hawaii, and Guam. We estimated average annual transmission intensity (i.e., force of infection) of DENV between 2010 and 2019 and the level of seroprevalence by age group in each population. We compared models and found that assuming all reported cases are secondary infections generally fit the surveillance data better than assuming all cases are primary infections. Using the secondary case model, we found that force of infection was highly heterogeneous between jurisdictions and over time within jurisdictions, ranging from 0.00008 (95% CrI: 0.00002–0.0004) in Florida to 0.08 (95% CrI: 0.044–0.14) in American Samoa during the 2010–2019 period. For early 2020, we estimated that seropositivity in 10 year-olds ranged from 0.09% (0.02%–0.54%) in Florida to 56.3% (43.7%–69.3%) in American Samoa. In the absence of serological data, age-specific case notification data collected through routine surveillance combined with mathematical modeling are powerful tools to monitor arbovirus circulation, estimate the level of population immunity, and design dengue management strategies.

60 APPLIED LIFE SCIENCES↗

Nine-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Final Maps and Results

We present the final nine-year maps and basic results from the Wilkinson Microwave Anisotropy Probe (WMAP) mission. The full nine-year analysis of the time-ordered data provides updated characterizations and calibrations of the experiment. We also provide new nine-year full sky temperature maps that were processed to reduce the asymmetry of the effective beams. Temperature and polarization sky maps are examined to separate cosmic microwave background (CMB) anisotropy from foreground emission, and both types of signals are analyzed in detail.We provide new point source catalogs as well as new diffuse and point source foreground masks. An updated template-removal process is used for cosmological analysis; new foreground fits are performed, and new foreground reduced are presented.We nowimplement an optimal C(exp -1)1 weighting to compute the temperature angular power spectrum. The WMAP mission has resulted in a highly constrained Lambda-CDM cosmological model with precise and accurate parameters in agreement with a host of other cosmological measurements. When WMAP data are combined with finer scale CMB, baryon acoustic oscillation, and Hubble constant measurements, we find that big bang nucleosynthesis is well supported and there is no compelling evidence for a non-standard number of neutrino species (N(sub eff) = 3.84 +/- 0.40). The model fit also implies that the age of the universe is (sub 0) = 13.772 +/- 0.059 Gyr, and the fit Hubble constant is H(sub 0) = 69.32 +/- 0.80 km/s/ Mpc. Inflation is also supported: the fluctuations are adiabatic, with Gaussian random phases; the detection of a deviation of the scalar spectral index from unity, reported earlier by the WMAP team, now has high statistical significance (n(sub s) = 0.9608+/-0.0080); and the universe is close to flat/Euclidean (Omega = −0.0027+0.0039/−0.0038). Overall, the WMAP mission has resulted in a reduction of the cosmological parameter volume by a factor of 68,000 for the standard six-parameter Lambda-CDM model, based on CMB data alone. For a model including tensors, the allowed seven-parameter volume has been reduced by a factor 117,000. Other cosmological observations are in accord with the CMB predictions, and the combined data reduces the cosmological parameter volume even further.With no significant anomalies and an adequate goodness of fit, the inflationary flat Lambda-CDM model and its precise and accurate parameters rooted in WMAP data stands as the standard model of cosmology.

dark matter↗

Characterizing the Roles of Biogeochemical Cycling and Ocean Circulation in Regulating Marine Copper Distributions

Copper (Cu) is a key micronutrient for marine phytoplankton. Its oceanic biogeochemical cycle has elicited considerable attention due to dissolved Cu exhibiting a unique linear profile with depth. Several processes have been proposed for explaining this behavior. Here, in this study, we characterize the relationships between the observed Cu, PO 4 3– , and Si on a global scale. We find that the depth profiles of Cu resemble those of Si more than of PO 4 3– in the global ocean. To understand their relationships, we couple the biogeochemical and internal circulation processes in a model, fitting optimal Cu:PO 4 3– uptake ratios and remineralization length-scales to replicate the marine Cu distributions. The modeling results suggest that Cu uptake needs in the Southern Ocean are substantially higher than those in other oceanic regions. In addition, our modeling results indicate a deep Cu remineralization in the global ocean. We offer an alternative mechanism that relies on biogeochemical cycling and internal circulation to produce the linear depth profiles of dissolved Cu. Our results suggest that diatoms are likely the major phytoplankton dominating oceanic Cu cycling.

54 ENVIRONMENTAL SCIENCES↗

Physics-informed Estimation of the Covariance Matrix for Various Neutron Spectra

A method for estimating covariance matrices which capture the uncertainties in calculated reactor spectra has been developed. This method is based on perturbing the parameters of a physics-based analytic model fitted to a calculated spectrum. The covariance of the perturbed analytic spectra imposes energy-dependent correlations due to the physics of the neutron processes in the reactor, i.e., a fission component, a 1/E down-scatting component, and a thermal Maxwellian component. An analytic model is developed which is shown to produce good fits to several reactor environments. The covariance matrices produced via this method are then used as the prior spectrum in STAYSL least squares spectrum adjustment where it is combined with integral metrics, such as activation measurements, to produce a high-fidelity neutron spectrum characterization. It was concluded that the methodology showed agreeable results for the ACRR free-field spectrum adjustment in STAYSL resulting in a 𝜒 2 value of 2.21 (per degree of freedom), but further work is needed to describe scattering and interface regions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Primary and Low-Strain Creep Models for 9Cr Tempered Martensitic Steels Including the Effects of Irradiation Softening and High-Helium Re-Hardening

Primary and low-strain creep represents a very important integrity challenge to large, complex structures, like fusion reactors. Here, we develop a predictive empirical primary creep model for 9Cr tempered martensitic steels (TMS), relating the applied stress (σ) to strain (ε), time (t) and temperature (T). The most accurate model is based on the applied σ normalized by the steel’s T-dependent ultimate tensile stress (σo), σ/σo(T). The model, fit to 17 heats of 9Cr TMS, yielded a σ root mean square error (RMSE) of ≈±11 MPa. Notably, the model also provides robust predictions for all the other TMS, when calibrated only by the fusion candidate Eurofer97 database. The model was extended to explore two possible effects of neutron irradiation, which produces both displacements per atom (dpa) and helium (He in atomic parts per million, appm) damage. These effects, which have not been previously considered, include: (a) softening, as a function of dpa, at T > ≈400–450 °C, in low-He fission environments (<1 He/dpa); and (b) subsequent re-hardening in high-He (≥10 He/dpa) fusion first-wall environments. The irradiation effect models predict (a) accelerated primary creep due to irradiation softening; and (b) fully arrested creep due to high-He re-hardening.

Alam, Md Ershadul (ORCID:0000000345968026)↗