Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Generalized Linear 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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

SAIGE-GPU: accelerating genome- and phenome-wide association studies using GPUs

Genome-wide association studies (GWAS) at biobank scale are computationally intensive, especially for admixed populations requiring robust statistical models. SAIGE is a widely used method for generalized linear mixed-model GWAS but is limited by its CPU-based implementation, making phenome-wide association studies impractical for many research groups. We developed SAIGE-GPU, a GPU-accelerated version of SAIGE that replaces CPU-intensive matrix operations with GPU-optimized kernels. The core innovation is distributing genetic relationship matrix calculations across GPUs and communication layers. Applied to 2068 phenotypes from 635 969 participants in the Million Veteran Program, including diverse and admixed populations, SAIGE-GPU achieved a 5-fold speedup in mixed model fitting on supercomputing infrastructure and cloud platforms. We further optimized the variant association testing step through multi-core and multi-trait parallelization. Deployed on Google Cloud Platform and Azure, the method provided substantial cost and time savings. Source code and binaries are available for download at https://github.com/saigegit/SAIGE/tree/SAIGE-GPU-1.3.3. A code snapshot is archived at Zenodo for reproducibility (DOI: [10.5281/zenodo.17642591]). SAIGE-GPU is available in a containerized format for use across HPC and cloud environments and is implemented in R/C++ and runs on Linux systems.

Rodriguez, Alex [Argonne National Laboratory (ANL)↗
DOE OSTI2017

Calculation of standard viscoelastic responses with multiple retardation times through analysis of static force spectroscopy AFM data

Here, we explore the physics of an atomic force microscopy (AFM) cantilever tip interacting with a generalized viscoelastic sample containing an arbitrary number of characteristic times, when the cantilever's base is driven with constant velocity toward the sample. This mode of operation, often called static force spectroscopy (SFS), can be harnessed to thoroughly analyze time-dependent viscoelastic information frequently overlooked in experiments. We generalize the solution of previous authors who have studied the standard linear solid model, and offer a solution applicable to any linear viscoelastic model. This generalization is crucial for the prediction of the model's response over wide ranges of time-scale. As a demonstration, successful predictions of harmonic functions (e.g., loss tangent) over a wide frequency range are obtained through analysis of simulated SFS results. In addition, we show that analysis through the generalized solution and previous expressions is no longer valid when the force does not grow linearly in time, so we also deliver an alternate route for extracting the viscoelastic information, which does not rely on the force linearity assumption. Despite the large amount of theoretical content (included for theoretical rigor's sake), the practical user can also benefit from the new procedures offered and the corresponding explanations.

42 ENGINEERING↗
DOE OSTI2010

Quantum effects on Higgs-strahlung events at linear colliders within the general two-Higgs-doublet model

The associated production of neutral Higgs bosons with the Z{sup 0} gauge boson (h{sup 0}Z{sup 0}, H{sup 0}Z{sup 0}) is investigated in the context of the future linear colliders, such as the ILC and CLIC, within the general two-Higgs-doublet model (2HDM). We compute the corresponding cross section for the processes e{sup +}e{sup -{yields}}Z{sup 0}h{sup 0}/Z{sup 0}H{sup 0} at one-loop, including the full set of corrections at O({alpha}{sub ew}{sup 3}) together with the leading O({alpha}{sub ew}{sup 4}) terms, in full consistency with the available theoretical and phenomenological constraints. We find that the wave-function corrections to the external Higgs fields are the dominant source of the quantum effects, which turn out to be large and negative (e.g., {delta}{sigma}/{sigma}{approx}-20%/-60%) in all the {radical}(s) range, and located predominantly in the region around tan{beta}{approx}1 and moderate values of the parameter {lambda}{sub 5} (being {lambda}{sub 5}<0). This behavior can be ultimately traced back to the enhancement potential of the triple Higgs-boson self-couplings, a trademark feature of the 2HDM with no counterpart in the Higgs sector of the minimal supersymmetric standard model. Even under this substantial depletion of the one-loop-corrected signal (which is also highly distinctive with respect to the SM expectation for e{sup +}e{sup -{yields}}Z{sup 0}H), the predicted Higgs-strahlung rates comfortably reach a few tens of femtobarn, which means barely {approx}10{sup 3}-10{sup 4} events per 500 fb{sup -1} of integrated luminosity. Because of their great complementarity, we argue that the combined analysis of the Higgs-strahlung events (h{sup 0}Z{sup 0}, H{sup 0}Z{sup 0}) and the previously computed one-loop Higgs-pair production processes (H{sup 0}A{sup 0}) could be instrumental to probe the structure of the Higgs sector at future linac facilities.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Modelling intergranular and transgranular micro-cracking in polycrystalline materials

Presented is a grain boundary formulation for intergranular and transgranular micro-cracking in three-dimensional polycrystalline aggregates.. The formulation is based on the displacement and stress boundary integral equations of solid mechanics and it has the advantage of expressing the polycrystalline problem in terms of grain boundary variables only. The individual grains within the polycrystalline morphology are modelled as generally anisotropic linear elastic domains with random spatial orientation. Transgranular micro-cracking is assumed to occur along specific cleavage planes, whose orientation in space within the grains depend upon the crystallographic lattice. Both intergranular and transgranular micro-cracking are modelled using suitably defined cohesive laws, whose parameters characterise the behaviour of the two mechanisms. The algorithm developed to track the inter/transgranular micro-cracking history is presented and discussed. Several numerical tests involving pseudo-3D and fully 3D morphologies are performed and analysed. The presented numerical results show that the developed formulation is capable of tracking the initiation and evolution of both intergranular and transgranular cracking as well as their competition, thus providing a useful tool for the study of damage micro-mechanics.

36 MATERIALS SCIENCE↗
DOE OSTI2017

An Analysis of Terrestrial and Aquatic Environmental Controls of Riverine Dissolved Organic Carbon in the Conterminous United States

Riverine carbon cycling is an important, but insufficiently investigated component of the global carbon cycle. Analyses of environmental controls on riverine carbon cycling are critical for improved understanding of mechanisms regulating carbon processing and storage along the terrestrial-aquatic continuum. Here, we compile and analyze riverine dissolved organic carbon (DOC) concentration data from 1402 United States Geological Survey (USGS) gauge stations to examine the spatial variability and environmental controls of DOC concentrations in the United States (U.S.) surface waters. DOC concentrations exhibit high spatial variability, with an average of 6.42 ± 6.47 mg C/ L (Mean ± Standard Deviation). In general, high DOC concentrations occur in the Upper Mississippi River basin and the Southeastern U.S., while low concentrations are mainly distributed in the Western U.S. Single-factor analysis indicates that slope of drainage areas, wetlands, forests, percentage of first-order streams, and instream nutrients (such as nitrogen and phosphorus) pronouncedly influence DOC concentrations, but the explanatory power of each bivariate model is lower than 35%. Analyses based on the general multi-linear regression models suggest DOC concentrations are jointly impacted by multiple factors. Soil properties mainly show positive correlations with DOC concentrations; forest and shrub lands have positive correlations with DOC concentrations, but urban area and croplands demonstrate negative impacts; total instream phosphorus and dam density correlate positively with DOC concentrations. Notably, the relative importance of these environmental controls varies substantially across major U.S. water resource regions. In addition, DOC concentrations and environmental controls also show significant variability from small streams to large rivers, which may be caused by changing carbon sources and removal rates by river orders. In sum, our results reveal that general multi-linear regression analysis of twenty one terrestrial and aquatic environmental factors can partially explain (56%) the DOC concentration variation. In conclusion, this study highlights the complexity of the interactions among these environmental factors in determining DOC concentrations, thus calls for processes-based, non-linear methodologies to constrain uncertainties in riverine DOC cycling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Do Remote Camera Arrangements and Image Capture Settings Improve Individual Identification of Golden Eagles?

Individual identification of animals from camera traps has become an important task in wildlife research, but camera deployment methods often do not facilitate this important undertaking. Identification of individual golden eagles (Aquila chrysaetos) is possible using uniquely marked rectrices, but no studies have explored methods to maximize the rate of individual identification from camera images. Furthermore, our objectives were to assess whether different camera heights (1 m vs. 3 m), image capture settings (one image after a 1-min delay vs. burst of 5 images after a 30 sec delay), and arrangements relative to bait (dorsally vs. ventrally aimed) affected views of rectrices on golden eagles and our ability to identify individuals. We conducted our study from 15 December 2016 to 3 March 2017 on the Savannah River Site, South Carolina. First, we developed a scoring system based on views of rectrices and used a linear mixed-effects model to compare image scores among different camera arrangements and image settings. Next, after identifying individual eagles, we used generalized linear mixed-effects models to compare total individual eagle detections, total days an individual was detected, and probability of obtaining an unknown individual identification among camera arrangements and settings. Overall, we scored a total of 27,499 images, with 8,083 providing views of marked rectrices that allowed identification of 18 individual eagles. Average image scores and proportion of images suitable for individual identification were higher from elevated (3 m) camera arrangements than standard arrangements (1 m) across sites. Regardless of camera height, faster frequency of image capture provided more images that could be used to identify individuals and the most trap days per individual. Researchers and managers should consider deploying elevated cameras traps with faster frequency of image capture to improve data quality and potential for analysis of golden eagle populations and trends across the species’ range.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

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↗
DOE OSTI2017

Conservation planning for freshwater–marine carryover effects on Chinook salmon survival

Experiences of migratory species in one habitat may affect their survival in the next habitat, in what is known as carryover effects. These effects are especially relevant for understanding how freshwater experience affects survival in anadromous fishes. Here, we study the carryover effects of juvenile salmon passage through a hydropower system (Snake and Columbia rivers, northwestern United States). To reduce the direct effect of hydrosystem passage on juveniles, some fishes are transported through the hydrosystem in barges, while the others are allowed to migrate in-river. Although hydrosystem survival of transported fishes is greater than that of their run-of-river counterparts, their relative juvenile-to-adult survival (hereafter survival) can be less. We tested for carryover effects using generalized linear mixed effects models of survival with over 1 million tagged Chinook salmon, Oncorhynchus tshawytscha (Walbaum) (Salmonidae), migrating in 1999–2013. Carryover effects were identified with rear-type (wild vs. hatchery), passage-type (run-of-river vs. transported), and freshwater and marine covariates. Importantly, the Pacific Decadal Oscillation (PDO) index characterizing cool/warm (i.e., productive/nonproductive) ocean phases had a strong influence on the relative survival of rear- and passage-types. Specifically, transportation benefited wild Chinook salmon more in cool PDO years, while hatchery counterparts benefited more in warm PDO years. Transportation was detrimental for wild Chinook salmon migrating early in the season, but beneficial for later season migrants. Hatchery counterparts benefited from transportation throughout the season. Altogether, wild fish could benefit from transportation approximately 2 weeks earlier during cool PDO years, with still a benefit to hatchery counterparts. Furthermore, we found some support for hypotheses related to higher survival with increased river flow, high predation in the estuary and plume areas, and faster migration and development-related increased survival with temperature. Thus, pre- and within-season information on local- and broad-scale conditions across habitats can be useful for planning and implementing real-time conservation programs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Estimating mortality rates among passerines caught for ringing with mist nets using data from previously ringed birds

Abstract Mist netting is the most commonly used method for catching birds for scientific ringing, but despite decades of use, there have been few attempts to quantify the associated potential risks to the individuals caught. Any incidence of mortality through capture and handling, however low, is of potential ethical concern and may also introduce biases into the data. We estimate the mortality rate associated with capture of previously ringed (recaptured) passerines from the British and Irish Ringing Scheme (c. 1.5 million records) caught using mist nets. The importance of species, age, mass, month, time, previous captures, and an index of predator occurrence was tested using generalized linear mixed‐effects models. The average mortality rate was 0.0011, most of which was reported to occur before the individuals had been extracted from the nets (c. 70% of incidents). Juveniles appeared to be at higher risk and the incidence of predation from mist nets was seasonal, with increased risk during the winter. Species differed in their reported mortality rates with the apparent risk being greatest for Chiffchaff Phylloscopus collybita (0.0029) and Bullfinch Pyrrhula pyrrhula (0.0027). To improve our understanding (and hence minimize risk in future), we recommend collecting more complete data on incidences of mortality, and also injuries; exercising increased care when the species we have identified as being at greater risk are likely to be captured, and ensuring there are robust procedures for the checking of nets (as most reported incidents of mortality occur before handling). We also recommend that all Ringing Schemes should collate and make available data on capture‐related mortality. Overall rates of mortality associated with capture, although, were low and support the use of mist netting as a safe capture technique, without undue bias from mortality, when used by appropriately trained individuals.

Clewley, Gary D.↗
DOE OSTI2017

A new approach to estimate fecundity rate from inter‐birth intervals

Abstract Accurate estimates of fecundity rate are key to population assessments and effectively direct conservation efforts. We present a new approach to estimate fecundity rate based on the probability of a female giving birth, conditional on a previous birth t years ago, from which an expected inter‐birth interval ( IBI ) can be estimated. We use generalized linear mixed‐effects models to account for individual and temporal variability and apply the approach to individual reproductive histories of bottlenose dolphins ( Tursiops truncatus ) from the east coast of Scotland. We estimate a fecundity rate of 0.222 (95% CI = 0.218–0.253) and an expected IBI of 4.49 yr (95% CI = 3.94–4.93 yr). We use simulated data samples to show that the approach produces estimates with a minimum bias of <3%. Simulations are also used to investigate the effect of the most common data‐driven biases in the estimates of birth intervals and fecundity rate; we recommend longitudinal studies of at least 10 yr and capture probabilities of at least 0.3 when using this methodology. The approach may be modified to incorporate other parameters of interest and should be applicable to any population with comprehensive data on birth intervals.

Arso Civil, Mònica↗
DOE OSTI2013

Estimation of a Self-Consistent Set of Radiobiological Parameters From Hypofractionated Versus Standard Radiation Therapy of Prostate Cancer

Purpose: To determine a self-consistent set of radiobiological parameters in prostate cancer. Methods and Materials: A method to estimate intrinsic radiosensitivity (α), fractionation sensitivity (α/β), repopulation doubling time, number of clonogens, and kick-off time for accelerated repopulation of prostate cancer has been developed. Based on the generalized linear-quadratic model and without assuming the isoeffective hypothesis, the potential applications of the method were investigated using the clinical outcome of biochemical relapse-free survival recently reviewed in the literature. The strengths and limitations of the method, regarding the fitted parameters and 95% confidence intervals (CIs), are also discussed. Results: Our best estimate of α/β is 2.96 Gy (95% CI 2.41-3.53 Gy). The corresponding α value is 0.16 Gy{sup −1} (95% CI 0.14-0.18 Gy{sup −1}), which is compatible with a realistic number of clonogens: 6.5 × 10{sup 6} (95% CI 1.5 × 10{sup 6}-2.1 × 10{sup 7}). The estimated cell doubling time is 5.1 days (95% CI 4.2-7.2 days), very low if compared with that reported in the literature. This corresponds to the dose required to offset the repopulation occurring in 1 day of 0.52 Gy/d (95% CI 0.32-0.68 Gy/d). However, a long kick-off time of 31 days (95% CI 22-41 days) from the start of radiation therapy was found. Conclusion: The proposed analytic/graphic method has allowed the fitting of clinical data, providing a self-consistent set of radiobiological parameters for prostate cancer. With our analysis we confirm a low value for α/β with a correspondingly high value of intrinsic radiosensitivity, a realistic average number of clonogens, a long kick-off time for accelerated repopulation, and a surprisingly fast repopulation that suggests the involvement of subpopulations of specifically tumorigenic stem cells during continuing radiation therapy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2023

Nonlinear mode coupling and energetics of driven magnetized shear-flow turbulence

To comprehensively understand the saturation of two-dimensional (2D) magnetized Kelvin–Helmholtz-instability-driven turbulence, energy transfer analysis is extended from the traditional interaction between scales to include eigenmode interactions, by using the nonlinear couplings of linear eigenmodes of the ideal instability. While both kinetic and magnetic energies cascade to small scales, a significant fraction of turbulent energy deposited by unstable modes in the fluctuation spectrum is shown to be re-routed to the conjugate-stable modes at the instability scale. They remove energy from the forward cascade at its inception. The remaining cascading energy flux is shown to attenuate exponentially at a small scale, dictated by the large-scale stable modes. Guided by a widely used instability-saturation assumption, a general quasi-linear model of instability is tested by retaining all nonlinear interactions except those that couple to the large-scale stable modes. These complex interactions are analytically removed from the magnetohydrodynamic equations using a novel technique. Observations are an explosive large-scale vortex separation instead of the well-known merger of 2D, a dramatic enhancement in turbulence level and spectral energy fluxes, and a reduced small-scale dissipation length scale. These show the critical role of the stable modes in instability saturation. Possible reduced-order turbulence models are proposed for fusion and astrophysical plasmas, based on eigenmode-expanded energy transfer analyses.

Physics↗
DOE OSTI2021

Nonadiabatic transitions in Landau-Zener grids: Integrability and semiclassical theory

We demonstrate that the general model of a linearly time-dependent crossing of two energy bands is integrable. Namely, the Hamiltonian of this model has a quadratically time-dependent commuting operator. We apply this property to four-state Landau-Zener (LZ) models that have previously been used to describe the Landau-Stückelberg interferometry experiments with an electron shuttling between two semiconductor quantum dots. The integrability then leads to simple but nontrivial exact relations for the transition probabilities. In addition, the integrability leads to a semiclassical theory that provides analytical approximation for the transition probabilities in these models for all parameter values. Lastly, the results predict a dynamic phase transition, and show that similarly looking models belong to different topological classes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2011

Determinants of Organophosphorus Pesticide Urinary Metabolite Levels in Young Children Living in an Agricultural Community

Organophosphorus (OP) pesticides are used in agriculture and several are registered for home use. As young children age they may experience different pesticide exposures due to varying diet, behavior, and other factors. We measured six OP dialkylphosphate (DAP) metabolites (three dimethyl alkylphosphates (DMAP) and three diethyl alkylphosphates (DEAP)) in urine samples collected from ~400 children living in an agricultural community when they were 6, 12, and 24 months old. We examined bivariate associations between DAP metabolite levels and determinants such as age, diet, season, and parent occupation. To evaluate independent impacts, we then used generalized linear mixed multivariable models including interaction terms with age. The final models indicated that DMAP metabolite levels increased with age. DMAP levels were also positively associated with daily servings of produce at 6- and 24-months. Among the 6-month olds, DMAP metabolite levels were higher when samples were collected during the summer/spring versus the winter/fall months. Among the 12-month olds, DMAP and DEAP metabolites were higher when children lived ≤60 meters from an agricultural field. Among the 24-month-olds, DEAP metabolite levels were higher during the summer/spring months. Our findings suggest that there are multiple determinants of OP pesticide exposures, notably dietary intake and temporal and spatial proximity to agricultural use. The impact of these determinants varied by age and class of DAP metabolite.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Data-driven modeling and control of dynamical systems using Koopman and Perron-Frobenius operators

This dissertation studies the data-driven modeling and control problem of nonlinear systems by exploiting the linear operator theoretic framework involving Koopman and Perro-Frobenius operator. A systematic linear-operator based controller design procedure has been established, which can be used to solve a variety of nonlinear control problems, including feedback stabilization using control Lyapunov functions, optimal quadratic regulation using Koopman eigenfunctions and convex optimization formulation of optimal control problem using P-F and Koopman operator approximation. As the core of data-driven modeling, we first propose a new algorithm for the finite-dimensional approximation of the linear transfer Koopman and Perron-Frobenius operator from time-series data. We argue that the existing approach for the finite-dimensional approximation of these transfer operators such as Dynamic Mode Decomposition (DMD) and Extended Dynamic Mode Decomposition (EDMD) do not capture two important properties of these operators, namely positivity and Markov property. The algorithm we propose preserves these two properties. We call the proposed algorithm as naturally structured DMD (NSDMD) since it retains the inherent properties of these operators. Naturally structured DMD algorithm leads to a better approximation of the steady-state dynamics of the system regarding computing Koopman and Perron- Frobenius operator eigenfunctions and eigenvalues. However, preserving positivity property is critical for capturing the real transient dynamics of the system. This positivity property of the transfer operators and it's finite-dimensional approximation play an important role for controller and estimator design of nonlinear systems. To solve the feedback stabilization problem for nonlinear control systems, we tried to take advantage of the Koopman operator framework. The Koopman operator approach provides a linear representation for a nonlinear dynamical system and a bilinear representation for a nonlinear control system. The problem of feedback stabilization of a nonlinear control system is then transformed to the stabilization of a bilinear control system. We propose a control Lyapunov function (CLF)-based approach for the design of stabilizing feedback controllers for the bilinear system. The search for finding a CLF for the bilinear control system is formulated as a convex optimization problem. This leads to a schematic procedure for designing CLF-based stabilizing feedback controllers for the bilinear system and hence the original nonlinear system. Another advantage of the proposed controller design approach outlined in this dissertation is that it does not require explicit knowledge of system dynamics. In particular, the bilinear representation of a nonlinear control system in the Koopman eigenfunction space can be obtained from time-series data. Next, we study the optimal quadratic regulation problem for nonlinear systems. The linear operator theoretic framework involving the Koopman operator is used to lift the dynamics of nonlinear control system to an infinite-dimensional bilinear system. The optimal quadratic regulation problem for nonlinear system is formulated in terms of the finite-dimensional approximation of the bilinear system. A convex optimization-based approach is proposed for solving the quadratic regulator problem for bilinear system. We applied a variety of examples and compared the simulation results between our framework and conventional LQR control using linearized model. For more general optimal control problems, we provide a density-function based convex formulation for the optimal control problem of the nonlinear system. The convex formulation relies on the duality result in the stability theory of a dynamical system involving density function and Perron-Frobenius operator. The optimal control problem is formulated as an infinite-dimensional convex optimization program. The finite-dimensional approximation of the optimization problem relies on the recent advances made in the data-driven computation of the Koopman operator, which is dual to the Perron-Frobenius operator. Simulation results are presented to demonstrate the application of the developed framework.

Huang, Bowen↗
DOE OSTI2012

Sensitivity of posterior parameter correlations and response parameter correlations to prior response uncertainty

The two parameters one response model of the generalized linear least squares parameter adjustment methodology was used in order to analyze the procedure of improving cross sections agreement with integral experimental results. Two extreme adjustment cases were considered. One case is using extremely accurate measured integral responses, such as k{sub eff} of Godiva or Jezebel for instance. The other one is modifying only one parameter, for instance the modification of {nu}-bar in ENDF/B-VII, in order to improve agreement of calculated responses with corresponding integral measurements results. In both cases new posterior parameters correlations are generated. The conclusion is that it is not advised to use integral measurements information in the evaluation process of a general purpose nuclear data file. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Bayesian Monte Carlo Evaluation of Imperfect (n, 233 U) Data and Model

Conventional nuclear data evaluation methods using generalized linear least squares make the following assumptions: prior and posterior probability distribution functions (PDFs) of all model parameters and data are normal (Gaussian); the linear approximation is sufficiently accurate to minimize the cost function (even for nonlinear models); the model (e.g., of neutron cross section) and experimental data (including covariance data) are without defect and prior PDFs of parameters and measured data are known perfectly. Neglect of covariance between model parameters and measured data in conventional evaluations contributes to imperfections. These assumptions are inherent to the generalized linear least squares minimization method commonly used for resolved resonance region neutron cross section evaluations but are often not justified due to the presence of non-normal PDFs, nonlinear models (e.g., R-matrix formalism), and inherent imperfections in data and models (e.g., imperfect covariance data). Here, these assumptions are removed in a mathematical framework of Bayes’ theorem, which is implemented using the Metropolis-Hastings Monte Carlo method. Most importantly, new parameters are introduced to parameterize discrepancies between the theoretical model and measured data to quantify judgement about discrepancies or imperfections in a reproducible manner. An evaluation of 233U in the eV region using the ENDF-B/VIII.0 library and transmission data (Guber et al.) is presented, and posterior parameters are compared to those obtained by conventional evaluation methods. This example illustrates the effects of removing the most harmful assumption: that of model-data perfection.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗