Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Data 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 1,711 records · Page 95
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

On extraction of twist-two Compton form factors from DVCS observables through harmonic analysis

We investigate the exercise of locally extracting the real and imaginary parts of the four twist-2 Compton form factors (CFFs) {$\mathcal{H,E,\tilde H,\tilde E}$} which arise in the deeply virtual Compton scattering (DVCS) process e + p → e + p + y. Neglecting dynamical higher-twist contributions, we find that there are a sufficient number of DVCS observables and degrees of freedom to extract all 8 leading quantities model-independently, exploiting the azimuthal dependence of the absolute cross sections across all possible beam and target polarizations at a common kinematical point in {Q, t,x b , y}. As an example, for typical JLab lab-frame kinematics, we simplify the reduced DVCS observables to their dominant terms, providing a sufficient number of equations for local determination of the twist-2 CFFs. We demonstrate the feasibility using harmonic fitting to both cross sections and beam spin asymmetries with both real and pseudo-data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Joint estimation of biogeochemical model parameters from multiple experiments: A bayesian approach applied to mercury methylation

Here, to characterize complex biogeochemical systems, results from multiple experiments, where each targets a specific subprocess, are commonly combined. The resulting datasets are interpreted through the calibration of biogeochemical models for process inference and predictions. Commonly used calibration approaches of fitting datasets from individual experiments to subprocess models one at a time is prone to missing information shared between datasets and incomplete uncertainty propagation. We propose a Bayesian joint-fitting scheme addressing the above-mentioned concerns by jointly fitting all the available datasets, thus calibrating the entire biogeochemical model in one go using Markov Chain Monte Carlo (MCMC). The identification of null spaces in the parameter distributions from MCMC guided the simplification of certain subprocess models. For example, fast kinetic sorption was replaced by equilibrium sorption, and Monod demethylation was replaced by first-order demethylation. Joint fitting of datasets resulted in complete uncertainty propagation with parameter estimates informed by all available data.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

INTERFACES. A Program for Determining the 3D Structures of Surfaces Sites Using NMR Data

Dynamic nuclear polarization surface enhanced NMR spectroscopy has enabled the determination of high-resolution structures from surface-supported molecules, including singlesite heterogeneous catalysts. Structure determinations have largely mimicked the approaches used in biomolecular NMR spectroscopy, namely, using distance measurements to constrain a conformational search. These early demonstrations made use of purpose-built software, which has limited the adoption of the technique. Herein, we describe the open-source program INTERFACES (Interpret NMR to Elucidate or Reconstruct the Full Atomistic Configurations of External Surfaces) which automates the analysis of RE(SP)DOR data as well as the structure determination for surface sites. Distances, angles, dihedral angles, complex orientation, and distance from the support can all be sampled to find all structures that agree with the experimental data. A χ 2 metric is used to define the error ranges of the REDOR fits and produce structures with an arbitrary level of confidence. Structural solutions are then provided as both overlays and ORTEP-like probability ellipsoids.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Theoretical power-law relationship between permeability and formation factor

Permeability k and formation factor F are important petrophysical properties that have broad applications to reservoir engineering. Although well-log analysis has advanced reservoir characterizations, direct measurement of k is still not available. Accordingly, resistivity logs and more particularly formation factor have been frequently used to estimate k. In this study, we present a theoretical framework to relate k to F. We use power-law scaling from continuum percolation theory to develop a theoretic relationship between k and F, which is similar in form to a published power-law equation with empirical coefficients. Our result gives a physical interpretation to this empirical relationship. We also compare power-law scaling laws from percolation theory with data derived from experimental measurements and numerical simulations in the literature. Finally, using the Curve Fitting toolbox of MATLAB, the parameters of the power-law equations are determined, and the significance of the non-universal power laws are discussed.

42 ENGINEERING↗
DOE OSTI2020

Comparing Experimental Phase Behavior of Ion-Doped Block Copolymers with Theoretical Predictions Based on Selective Ion Solvation

The effects of salt-doping on the morphological behavior of block copolymers are well established but remain poorly understood, partially because of the challenge of resolving electrostatics in a heterogeneous medium with low average permittivity. By employing a recently developed field theory, we analyze the phase behavior of polystyrene-b-poly(ethylene oxide) (SEO) copolymers doped with lithium bis(trifluoromethanesulfonyl)imide salts (LiTFSI). Using a single fitting parameter, the ionic solvation radius, we obtain qualitative agreement between our theory and experimental data over a range of polymer molecular weights and copolymer compositions. Such agreement supports and highlights the need of solvation free energy to accurately describe the self-assembly of ion-doped block copolymers and demonstrates that experimentally observed dependence on molecular weight, not present in neutral block copolymers, can be rationalized by solvation effects. Overall, morphological variations are stronger than those predicted by the leading, linear order theory but can be captured by the full model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Quantification and assessment of the atmospheric boundary layer height measured during the AWAKEN experiment by a scanning LiDAR

The atmospheric boundary layer (ABL) height plays a key role in many atmospheric processes as one of the dominant flow length scales. However, a systematic quantification of the ABL height over the entire range of scales (i.e., with periods ranging from one minute to one year) is still lacking in literature. In this work, the ABL height is quantified based on high-resolution measurements collected by a scanning pulsed Doppler LiDAR during the recent American WAKE experimeNt (AWAKEN) campaign. The high availability of ABL height estimates (≈2200 collected over one year and each of them based on 10-min averaged statistics) allows to robustly assess five different ABL height models, i.e., one for convective thermal conditions and four for stable conditions. Thermal condition is quantified by a stability parameter spanning three orders of magnitude and probed by near-ground 3D sonic anemometry. The free-atmosphere stability, quantified by the Brunt–Väisälä frequency, is both calculated from simultaneous radiosonde measurements and obtained from the best fit of two of the chosen ABL height models. Good agreement is found between the data and three of the chosen models, quantified by mean absolute errors on the ABL height between 281 and 585 m. Furthermore, the seasonal variability of the convective ABL height model parameters (−15% to +23% with respect to the year baseline) agrees with the variability of buoyancy-generated turbulence caused by the variation in solar radiation throughout the year.

17 WIND ENERGY↗
DOE OSTI2025

Emulating the Lyman-Alpha forest 1D power spectrum from cosmological simulations: new models and constraints from the eBOSS measurement

We present the Lyssa suite of high-resolution cosmological simulations of the Lyman-α forest designed for cosmological analyses. These 18 simulations have been run using the Nyx code with 40963 hydrodynamical cells in a 120 Mpc (∼ 81 Mpc/h) comoving box and individually provide sub-percent level convergence of the Lyman-α forest 1d flux power spectrum. We build a Gaussian process emulator for the Lyssa simulations in the lym1d likelihood framework to interpolate the power spectrum at arbitrary parameter values. We validate this emulator based on leave-one-out tests and based on the parameter constraints for simulations outside of the training set. We also perform comparisons with a previous emulator, showing a percent level accuracy and a good recovery of the expected cosmological parameters. Using this emulator we derive constraints on the linear matter power spectrum amplitude and slope parameters A Lyα and n Lyα . While the best-fit Planck ΛCDM model has A Lyα = 8.79 and n Lyα = -2.363, from DR14 eBOSS data we find that A Lyα < 7.6 (95% CI) and n Lyα = -2.369 ± 0.008. The low value of A Lyα , in tension with Planck, is driven by the correlation of this parameter with the mean transmission of the Lyman-α forest. This tension disappears when imposing a well-motivated external prior on this mean transmission, in which case we find A Lyα = 9.8 ± 1.1 in accordance with Planck.

Walther, Michael↗
DOE OSTI2025

A kinetic line-driven radiation operator and its application to Gyrokinetics

A velocity dependent, kinetic model for line radiation is developed for continuum kinetic codes. It has been implemented in the full-f gyrokinetic code Gkeyll. The total radiation for a charge state is modeled as an advection in velocity space with a form of $\nabla_v \cdot(v\nu(v)f(v))$, guaranteeing particle conservation. The velocity dependence (in the form of an effective frequency $\nu(v)$) is found through fitting the energy loss of the operator, i.e. the second velocity moment, to the radiation data in the OpenADAS database. Therefore, each individual transition does not need to be evaluated every time step, significantly reducing the computational cost of including line radiation in a kinetic model. The dependence on velocity instead of the usual, temperature, allows the radiation to be computed from non-Maxwellian electron distribution functions: We benchmark the model against a collisional radiative model using isotropic non-Maxwellian distribution functions. A velocity dependent model of radiation can more accurately describe the radiation in the more kinetic regimes expected in reactor-scale devices. The velocity dependence qualitatively captures the quantum mechanical need for a minimum velocity before any radiation occurs.

kinetic↗
DOE OSTI2025

Thermodynamics of the dipole-octupole pyrochlore magnet Ce 2 Hf 2 O 7 in applied magnetic fields

The recently discovered dipole-octupole pyrochlore magnet Ce 2 Hf 2 O 7 is a promising three-dimensional quantum spin liquid candidate which shows no signs of ordering at low temperature. Here we investigate the thermodynamic response to magnetic fields applied along the global [110] direction using specific heat measurements and fits using numerical methods, and solve the corresponding magnetic structure using neutron diffraction. Specific heat data in moderate fields are reproduced well, however, at high fields the agreement is not satisfactory. We especially observe a two-step release of entropy, a finding that demands a review of both theory and experiment. We address it within the framework of three possible scenarios, including an analysis of the crystal field Hamiltonian not restricted to the two-dimensional single-ion doublet subspace. We conclusively rule out two of these scenarios and find qualitative agreement with a simple model of field misalignment with respect to the crystalline direction. As a result, we discuss the implications of our findings for [111] applied fields and for future experiments on Ce 2 Hf 2 O 7 and its sister compounds.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Pion-nucleon scattering in baryon chiral perturbation theory combined with the 1/𝑁 𝑐 expansion

This work implements the combined baryon chiral perturbation theory (BChPT) and 1/𝑁 𝑐 expansions for pion-nucleon elastic scattering. The effective theory is based on the baryon sector dynamical spin-flavor 𝑆⁢𝑈⁡(4) symmetry emergent in the large 𝑁 𝑐 limit, whose breaking is controlled by the 1/𝑁 𝑐 expansion. The noncommutativity of the chiral and 1/𝑁 𝑐 expansions in unitarity corrections (loops) requires a linking of both expansions. As it was shown in the case of baryon masses and currents, the natural linking is the 𝜉 expansion, in which 𝒪⁡(𝑝) = 𝒪⁡(1/𝑁 𝑐 ) = 𝒪⁡(𝜉). The spin-flavor symmetry requires that the ground state baryons span an 𝑆⁢𝑈⁡(4) symmetric irreducible representation which implies that in particular 𝑁 and Δ are active degrees of freedom in the effective theory. The scattering amplitude is expanded to the next-to-next-to leading order in the 𝜉 expansion, corresponding to the one-loop contributions with the leading-order Lagrangian. The results are given for generic 𝑁 𝑐 in order to demonstrate the consistency of the framework. The spin-flavor symmetry plays a central role in maintaining the consistency of the effective theory with respect to the 1/𝑁 𝑐 expansion. This consistency manifests itself in an improvement in the convergence of the low energy expansion with respect to the case of the ordinary BChPT without an explicit dynamical Δ, which is known to be inconsistent with the constraints of 𝑁 𝑐 scaling. Fits to the 𝜋⁢𝑁 → 𝜋⁢𝑁, 𝑆, 𝑃, and 𝐷 partial wave amplitudes from the SAID data base are finally used to test the framework and to determine the energy range of its applicability.

Jayakodige, D. [Hampton Univ., Hampton, VA (United↗
DOE OSTI2020

Model Exploration of an Information-Based Healthcare Intervention Using Parallelization and Active Learning

This paper describes the application of a large-scale active learning method to characterize the parameter space of a computational agent-based model developed to investigate the impact of CommunityRx, a clinical information-based health intervention that provides patients with personalized information about local community resources to meet basic and self-care needs. Additionally, the diffusion of information about community resources and their use is modeled via networked interactions and their subsequent effect on agents' use of community resources across an urban population. A random forest model is iteratively fitted to model evaluations to characterize the model parameter space with respect to observed empirical data. We demonstrate the feasibility of using high-performance computing and active learning model exploration techniques to characterize large parameter spaces; by partitioning the parameter space into potentially viable and non-viable regions, we rule out regions of space where simulation output is implausible to observed empirical data. We argue that such methods are necessary to enable model exploration in complex computational models that incorporate increasingly available micro-level behavior data. We provide public access to the model and high-performance computing experimentation code.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Neutrino physics beyond the current paradigm

An exciting program of international neutrino experiments aims to complete our understanding of three-neutrino mixing in the Standard Model (SM), which has had great success in describing world neutrino data. However, there exist a number of short-baseline neutrino anomalies that do not fit into this simple picture and could be pointing us instead to phenomena beyond the Standard Model (BSM). Neutrinos and neutrino experiments are excellent laboratories to search for BSM physics for a number a reasons. First, neutrinos mass itself is evidence of BSM physics as it implies the addition of new fields to the SM. In addition, neutrino properties are not well-constrained and allow sizable, new interactions. Finally, accelerator- and reactor-based short-baseline neutrino experiments bring high luminosity sources of photons, nuclear, and meson decays near large, sensitive detectors and can therefore probe a variety of dark sector models. In this talk, I will give an overview of worldwide accelerator- and reactor-based short-baseline neutrino experiments and their efforts to search for physics beyond the current paradigm.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

An Applied Strategy for Using Empirical and Hybrid Models in Online Monitoring

The monitoring of plant equipment for failure prediction is one of the key contributors to operation and maintenance (O&M) costs for a nuclear power plant (NPP) because O&M monitoring depends on labor-intensive activities that are required to meet high equipment reliability standards. These activities rely primarily on humans for information gathering, condition diagnosis, and predictive analysis. Online monitoring aims to automate these activities by relying on sensors to replace human information gathering and machine learning to replace human analysis and decision making. To facilitate automated monitoring, a systematic strategy for anomaly detection is needed to optimally use the available sensor data, empirical models, and physics-supported models. This strategy is essential to provide credible reasoning on why and when an empirical (i.e., purely data-driven) versus hybrid (i.e., physics-supported) approach should be used and to determine the ideal mix of these two approaches for a defined anomaly detection scope. The extant methods usually adopt an ad hoc trial-and-error approach that, in addition to being time-consuming and costly, is also highly subjective; it is impacted by the background and the skill set of the personnel making the decisions. Thus, such an approach cannot guarantee an optimum outcome. This represents the motivation of the current research effort, which is focused on devising a scientifically supported strategy for the optimum selection of anomaly detection methods. This report presents a detailed assessment of the main anomaly detection techniques within the empirical or hybrid method streams. Empirical methods include pattern, statistical, and causal inference. Hybrid methods include the use of physics models to train and test data methods, reduce data dimensionality, reduce data-model complexity, augment data, and reduce empirical uncertainty; hybrid methods also include the use of data to tune physics models. The listed techniques within these two streams represent the vast majority of techniques performed for anomaly detection. Using the techniques as outcomes, a strategy was developed to enable a systematic decision-making process to lead to one of these techniques. The strategy is driven by key decision points related to data relevance, simple modeling feasibility, data inference, physics-modeling value, data dimensionality, physics knowledge, method of validation, performance, data availability and suitability for training and testing, cause-effect, entropy inference, and model fitting. Each of these decision points in the strategy is explained in detail in this report with examples, along with the scientific basis behind the decisions and outcomes in common and simplified terminology. The strategy is developed for use by any NPP staff with basic engineering or science knowledge. A user-friendly graphical state flow diagram was also developed as a visual presentation of the strategy. The strategy was tested and demonstrated through two pilot projects for the application of anomaly detection at an NPP. Each pilot had two use cases: an initial case in which certain decisions were made that resulted in one or more empirical techniques and a revised use case where one or more key decisions were modified resulting in using a set of hybrid methods.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

GADRAS-DRF Enhancements for Safeguards – Custom Peak Fitting to Enhance Model Fitting and Isotopics

One major software update was accomplished within the Gamma Detector Response and Analysis Software-Detector Response Function (GADRAS-DRF) package. This update allows users to adjust individual peak fits for use in subsequent analysis processes within GADRAS-DRF. A graphical user interface (GUI) was implemented so users can see the effect of their fit adjustments on the spectrum. This new feature will enhance the capability of the previously funded auto-enrichment and peak-based 1D model fit feature that was implemented in FY23 (Fiscal Year). Isotopic analysis was performed using M400 data obtained from uranium standards and the isotopic assessment is given with and without manual peak adjustments.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2022

A Qualitative Strategy for Fusion of Physics into Empirical Models for Process Anomaly Detection

To facilitate the automated online monitoring of power plants, a systematic and qualitative strategy for anomaly detection is presented. This strategy is essential to provide credible reasoning on why and when an empirical versus hybrid (i.e., physics-supported) approach should be used and to determine the ideal mix of these two approaches for a defined anomaly detection scope. Empirical methods are usually based on pattern, statistical, and causal inference. Hybrid methods include the use of physics models to train and test data methods, reduce data dimensionality, reduce data-model complexity, augment data, and reduce empirical uncertainty; hybrid methods also include the use of data to tune physics models. The presented strategy is driven by key decision points related to data relevance, simple modeling feasibility, data inference, physics-modeling value, data dimensionality, physics knowledge, method of validation, performance, data availability, and suitability for training and testing, cause-effect, entropy inference, and model fitting. The strategy is demonstrated through a pilot use case for the application of anomaly detection to capture a valve packing leak at the high-pressure coolant injection system of a nuclear power plant.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Determination of Constraint-Independent Crack Tip Opening Angle for Stable Crack Growth in High-Strength Ductile Steels

The crack tip opening angle (CTOA) is one of fracture toughness parameters that has been used for decades in describing large stable crack growth in thin-walled aerospace structures under the low-constraint conditions. Recently, the pipeline industry has developed a growing interest in using the CTOA parameter to serve as the minimum required fracture toughness to arrest dynamic crack propagation in modern gas transmission pipelines made of high-strength ductile steel. To meet this industrial need, the CTOA test standard ASTM E3039 was therefore developed for measuring the constant critical CTOA. ASTM E3039 recommends a drop weight tearing test (DWTT) specimen with a shallow crack for standard CTOA testing, but its CTOA may depend on the low constraint condition of the DWTT specimen at the crack tip. Verifying the constraint independence of the DWTT-measured CTOA thus becomes indispensable for applying CTOA toughness to the running fracture control in the pipeline design. For this purpose, the present paper evaluates critical CTOA values in a set of fracture toughness tests on single-edge notched bend (SENB) specimens with shallow and deep cracks, based on four CTOA estimation models. Among these, the Ln(P)-LLD linear fit model is similar to that recommended by ASTM E3039 for CTOA calculation. Fracture test data for X80 pipeline steel and HY80 structural steel were considered in the CTOA evaluation. The results showed that the four CTOA models were able to determine a constraint independent CTOA value for stable crack growth in the SENB specimens. As a result, a single, reliable, constant CTOA value could be determined regardless of the specimen geometry or the crack-tip constraint conditions. Therefore, the CTOA measured using ASTM E3039 is constraint-independent and transferable to use in cases of actual cracks propagating in gas transmission pipelines.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Photometry as a Proxy for Stellar Activity in Radial Velocity Analyses

Stellar activity remains a limiting factor in measuring precise planet parameters from radial velocity spectroscopy, not least in the search for Earth-mass planets orbiting in the habitable zones of Sun-like stars. One approach to mitigate stellar activity is to use combined analyses of both radial velocity and time-series photometry. We present an analysis of simultaneous disk-integrated photometry and radial velocity data of the Sun in order to determine the useful limits of a combined analysis. We find that simple periodogram or autocorrelation analysis of solar photometry give the correct rotation period <50% of the time. We therefore use a Gaussian process to investigate the time variability of solar photometry and to directly compare simultaneous photometry with radial velocity data. We find that the hyperparameter posteriors are relatively stable over 70 yr of solar photometry and the amplitude tracks the solar cycle. We observe good agreement between the hyperparameter posteriors for the simultaneous photometry and radial velocity data. Our primary conclusion is a recommendation to include an additional prior in Gaussian process fits to constrain the evolutionary timescale to be greater than the recurrence timescale (i.e., the rotation period) to recover more physically plausible and useful results. Our results indicate that such simultaneous monitoring may be a useful tool in enhancing the precision of radial velocity surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A New Infrared Criterion for Selecting Active Galactic Nuclei to Lower Luminosities

We present a spectroscopic and photometric analysis of a sample of 416,288 galaxies from the Sloan Digital Sky Survey (SDSS) matched to mid-infrared (mid-IR) data from the Wide-field Infrared Survey Explorer (WISE). By using a new spectroscopic fitting package, GELATO (Galaxy/AGN Emission Line Analysis TOol), we are able to retrieve emission line fluxes and uncertainties for SDSS spectra and robustly determine the presence of broad lines and outflowing components, enabling us to investigate WISE color space as a function of optical spectroscopic properties. In addition, we pursue spectral energy distribution template fitting to assess the relative active galactic nucleus (AGN) contribution and nuclear obscuration to compare to existing mid-IR selection criteria with WISE. We present a selection criterion in mid-IR color space to select AGNs with an ~80% accuracy and a completeness of ~16%. This is the first mid-IR color selection defined by solely using the distribution of Type I and Type II optical spectroscopic AGNs in WISE mid-IR color space. Our selection is an improvement of ~50% in the completeness of targeting spectroscopic AGNs with WISE down to an SDSS r < 17.77 mag. In addition, our new criterion targets a less-luminous population of AGNs, with on average lower [O ιιι] luminosities by~30% ( > 0.1 dex) compared to typical WISE color–color selections. With upcoming large photometric surveys without corresponding spectroscopy, our method presents a way to select larger populations of AGNs at lower AGN luminosities and higher nuclear obscuration levels than traditional mid-IR color selections.

79 ASTRONOMY AND ASTROPHYSICS↗