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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Incidence and predictors of tuberculosis among children receiving antiretroviral therapy in the Wolaita Zone: A retrospective cohort study

Tuberculosis is the leading cause of morbidity and mortality among children living with the human immunodeficiency virus (HIV), mainly in sub-Saharan Africa, including Ethiopia. Tuberculosis remains a significant health concern for HIV-positive children in Ethiopia. There is a paucity of data on the incidence and predictors of tuberculosis among children living with HIV on antiretroviral therapy in the Wolaita zone. Hence, this study aimed to assess the incidence and predictors of tuberculosis among children living with HIV on antiretroviral therapy in the Wolaita zone between January 2010 to December 2020. A retrospective cohort study was conducted among 389 children receiving antiretroviral therapy in Wolaita zone health facilities between January 2010 to December 2020. The checklist was adapted from the standardized antiretroviral treatment (ART) follow-up form currently used by the institutions’ ART clinics. The Kaplan-Meier survival function and Log-rank were used to estimate the survival for each categorical variable to compare the survival between different exposure groups. Both bivariable and multivariable parametric survival Gompertz models were fitted to identify predictors of tuberculosis among HIV-positive children. The association was summarized using an adjusted hazard ratio (AHR), and statistical significance was declared at 95% CI and p-value < 0.05. The goodness of the model fit was assessed using a Cox-Snell residual plot. The incidence rate of tuberculosis among children living with HIV was 3.5 (95% CI 2.7–4.5) per 100 child years. World Health Organization clinical stage III or IV (AHR = 2.31, 95% CI [1.26, 4.22]), hemoglobin level <10 g/dL (AHR = 2.87, 95% CI [1.51, 5.45]), fair or poor ART adherence (AHR = 4.4, 95% CI[2.18, 9.05]), underweight (AHR = 2.55, 95% CI [1.45, 4.51]), age >10 years (AHR = 3.62; 95% CI [1.29, 10.0]), and cotrimoxazole preventive therapy (AHR = 0.23; 95% CI [0.08, 0.60]) were among the independent predictors of TB occurrence. The incidence of tuberculosis among children on ART was high. HIV-positive children presenting with advanced disease staging (III and IV), anemia, “fair” and “poor” ART adherence, underweight, age above ten years, and not receiving cotrimoxazole preventive therapy were at higher risk of TB. Therefore, counseling on ART adherence, early diagnosis, and prompt treatment of anemia and malnutrition are recommended to avert tuberculosis.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Comparative modeling reveals the molecular determinants of aneuploidy fitness cost in a wild yeast model

Although implicated as deleterious in many organisms, aneuploidy can underlie rapid phenotypic evolution. However, aneuploidy will be maintained only if the benefit outweighs the cost, which remains incompletely understood. To quantify this cost and the molecular determinants behind it, we generated a panel of chromosome duplications in Saccharomyces cerevisiae and applied comparative modeling and molecular validation to understand aneuploidy toxicity. We show that 74%–94% of the variance in aneuploid strains’ growth rates is explained by the cumulative cost of genes on each chromosome, measured for single-gene duplications using a genomic library, along with the deleterious contribution of small nucleolar RNAs (snoRNAs) and beneficial effects of tRNAs. Machine learning to identify properties of detrimental gene duplicates provided no support for the balance hypothesis of aneuploidy toxicity and instead identified gene length as the best predictor of toxicity. Our results present a generalized framework for the cost of aneuploidy with implications for disease biology and evolution.

genic load↗
DOE OSTI2023

Integrating AlphaFold and deep learning for atomistic interpretation of cryo-EM maps

Abstract Interpretation of cryo-electron microscopy (cryo-EM) maps requires building and fitting 3D atomic models of biological molecules. AlphaFold-predicted models generate initial 3D coordinates; however, model inaccuracy and conformational heterogeneity often necessitate labor-intensive manual model building and fitting into cryo-EM maps. In this work, we designed a protein model-building workflow, which combines a deep-learning cryo-EM map feature enhancement tool, CryoFEM (Cryo-EM Feature Enhancement Model) and AlphaFold. A benchmark test using 36 cryo-EM maps shows that CryoFEM achieves state-of-the-art performance in optimizing the Fourier Shell Correlations between the maps and the ground truth models. Furthermore, in a subset of 17 datasets where the initial AlphaFold predictions are less accurate, the workflow significantly improves their model accuracy. Our work demonstrates that the integration of modern deep learning image enhancement and AlphaFold may lead to automated model building and fitting for the atomistic interpretation of cryo-EM maps.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Theoretical justification for heat flux limiter 0.15

In this section, we describe a semi-nonlocal kinetic model for the electron heat flux. The model has no free parameters. It takes the more physically motivated approach of placing a limiter on the perturbation to the electron distribution function, rather than the usual ad-hoc limiter on the heat flux. If the kinetic model is fitted by a heat flux limiter model, we find that the higher value of f e = 0.15 is a much closer fit than other models with f e ≃ 0.02.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Training material models using gradient descent algorithms

High temperature design requires accurate constitutive models to describe material inelastic deformation and failure behavior. Oftentimes, calibrating accurate models devolves into the problem of fitting the model parameters against experimental test data. Here, we present the pyopmat package, an open source framework for calibrating constitutive models against experiment data subjected to various loading conditions using machine learning techniques. The package calculates the exact gradient of the model response with respect to the parameters using a combination of automatic differentiation and the adjoint method. Given this exact gradient, we compare the performance of several gradient-based optimization techniques in fitting realistic constitutive models against data. Here, we demonstrate the efficiency and accuracy of our package through example problems using both synthetic data, generated using known parameter sets, under monotonic and cyclic loading conditions and also with an example applying the techniques developed here to actual high temperature creep-fatigue test data.

36 MATERIALS SCIENCE↗
DOE OSTI2022

ML Pipeline (Machine Learning Pipeline) [SWR-22-35]

The Machine Learning Pipeline (ML Pipeline) allows the user to fit a model to predict a dependent variable, y, based on a feature set, X. ML Pipeline gives the ability to automatically generate features with its Feature Engineering function, automatically select the most important features with its Feature Selection function, fit the model with its Fit function, and evaluate the model with its Evaluate function. ML Pipeline's functionality can not only help the user predict values for the desired variable it also gives the user a better understanding of which features are important and how effective the model is. ML Pipeline accomplishes all this, while remaining simple to use and interpret.

Schiek, Andrew↗
DOE OSTI2021

Fourier-like Thermal Relaxation of Nanoscale Explosive Hot Spots

Hot spots are local regions of high temperature that are widely considered to govern explosive initiation. Hot spot dynamics rests on a delicate balance between heat generation due to chemical reactions and heat loss through thermal conduction, making accurate determinations of the conductivity under extreme conditions a key component of predictive explosive models. We develop here an approach to directly determine the thermal transport properties of explosive hot spots with realistic initial structures through a combination of molecular dynamics (MD) and diffusive heat equation (HEq) modeling. Effective thermal conductivity values are determined by fitting HEq models to MD predictions of long timescale hot spot relaxation. The approach is applied to model hot spots in the molecular crystalline explosive 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) for a range of shock strengths and two limiting cases for impact orientation. Isotropic and anisotropic HEq models yield similar results, despite TATB exhibiting some of the largest and most anisotropic thermal conductivity values for explosive near normal conditions. The conductivity is found to be a strong function of density, which parametrically captures dependence on temperature, pressure, and material state. The associated root-mean-square errors of the fitted HEq models are approximately 5% of MD predicted final equilibrium temperatures. Here, the conductivity values determined here for TATB hot spots are considerably larger than those used in a prior hot spot criticality study, which may significantly impact predictions for critical hot spot sizes. The approach provides a convenient foundation for determining the effective thermal conductivity for hot spot problems in other explosives and directly yields information on reasonable approximations that might be taken in higher-level models for those materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Glass Property-Composition Models for Support of Hanford WTP LAW Facility Operation

Current plans for the River Protection Project envision starting to vitrifying low-activity waste (LAW) by 2023 using a Direct Feed Low-Activity Waste (DFLAW) approach and subsequently using a full-pretreatment approach. The Hanford Tank Waste Treatment and Immobilization Plant (WTP) LAW Facility will be operated and controlled using a LAW glass formulation algorithm (GFA), which requires several inputs based on research and development results. LAW glass property-composition models for several product quality and processing properties are key inputs for the LAW GFA. It is envisioned that the preliminary LAW GFA discussed by Kim and Vienna (2012) will be used for commissioning and initial radioactive operations of the WTP LAW Facility under Bechtel National, Inc. using the DFLAW approach. Then, an updated LAW GFA will be developed for implementation by the WTP operating contractor that takes over after WTP LAW Facility commissioning. This report documents the enhanced LAW glass property-composition models developed for use in the updated LAW GFA. The properties for which models were developed include Product Consistency Test (PCT) response, Vapor Hydration Test (VHT) response, viscosity at 1150 °C, electrical conductivity at 1150 °C, melter SO 3 tolerance at 1150 °C, and K-3 refractory corrosion at 1208 °C. Table S.1 lists the tables in this report that contain the recommended models for each of these properties. The model types recommended include partial quadratic mixture (PQM) models for viscosity, electrical conductivity, melter SO 3 tolerance and K-3 corrosion, bias corrected PQM model (bcPQM) for PCT, and logistic PQM model for VHT. The fits of model and validation subsets were found to be well predicted by the recommended models.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Glass Property-Composition Models for Support of Hanford WTP LAW Facility Operation

Current plans for the River Protection Project envision starting to vitrifying low-activity waste (LAW) by 2023 using a Direct Feed Low-Activity Waste (DFLAW) approach and subsequently using a full-pretreatment approach. The Hanford Tank Waste Treatment and Immobilization Plant (WTP) LAW Facility will be operated and controlled using a LAW glass formulation algorithm (GFA), which requires several inputs based on research and development results. LAW glass property-composition models for several product quality and processing properties are key inputs for the LAW GFA. It is envisioned that the preliminary LAW GFA discussed by Kim and Vienna (2012) will be used for commissioning and initial radioactive operations of the WTP LAW Facility under Bechtel National, Inc. using the DFLAW approach. Then, an updated LAW GFA will be developed for implementation by the WTP operating contractor that takes over after WTP LAW Facility commissioning. This report documents the enhanced LAW glass property-composition models developed for use in the updated LAW GFA. The properties for which models were developed include Product Consistency Test (PCT) response, Vapor Hydration Test (VHT) response, viscosity at 1150 °C, electrical conductivity at 1150 °C, melter SO3 tolerance at 1150 °C, and K-3 refractory corrosion at 1208 °C. Table S.1 lists the tables in this report that contain the recommended models for each of these properties. The model types recommended include partial quadratic mixture (PQM) models for viscosity, electrical conductivity, melter SO 3 tolerance and K-3 corrosion, bias corrected PQM model (bcPQM) for PCT, and logistic PQM model for VHT. The fits of model and validation subsets were found to be well predicted by the recommended models.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

First Sagittarius A* Event Horizon Telescope Results. IV. Variability, Morphology, and Black Hole Mass

In this paper we quantify the temporal variability and image morphology of the horizon-scale emission from Sgr A*, as observed by the EHT in 2017 April at a wavelength of 1.3 mm. We find that the Sgr A* data exhibit variability that exceeds what can be explained by the uncertainties in the data or by the effects of interstellar scattering. The magnitude of this variability can be a substantial fraction of the correlated flux density, reaching ~100% on some baselines. Through an exploration of simple geometric source models, we demonstrate that ring-like morphologies provide better fits to the Sgr A* data than do other morphologies with comparable complexity. We develop two strategies for fitting static geometric ring models to the time-variable Sgr A* data; one strategy fits models to short segments of data over which the source is static and averages these independent fits, while the other fits models to the full data set using a parametric model for the structural variability power spectrum around the average source structure. Both geometric modeling and image-domain feature extraction techniques determine the ring diameter to be 51.8 ± 2.3 μas (68% credible intervals), with the ring thickness constrained to have an FWHM between ~30% and 50% of the ring diameter. To bring the diameter measurements to a common physical scale, we calibrate them using synthetic data generated from GRMHD simulations. This calibration constrains the angular size of the gravitational radius to be ${4.8}_{-0.7}^{+1.4}$ μas, which we combine with an independent distance measurement from maser parallaxes to determine the mass of Sgr A* to be ${4.0}_{-0.6}^{+1.1}\times {10}^{6}$ M⊙.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

eXtremeMAT: Uncertainty Quantification of the LApx Model

Presentation on the uncertainty quantification efforts to parameterize and fit the LApx model for ferritic-martensitic steels and austenitic stainless steels. These workflows support ease of adoption of the code to new materials and enable rapid model fitting and understanding of the level of confidence in predictions.

mechanical properties↗
DOE OSTI2025

An analysis of parameter compression and Full-Modeling techniques with Velocileptors for DESI 2024 and beyond

In anticipation of forthcoming data releases of current and future spectroscopic surveys, we present the validation tests and analysis of systematic effects within velocileptors modeling pipeline when fitting mock data from the AbacusSummit N-body simulations. We compare the constraints obtained from parameter compression methods to the direct fitting (Full-Modeling) approaches of modeling the galaxy power spectra, and show that the ShapeFit extension to the traditional template method is consistent with the Full-Modeling method within the standard ΛCDM parameter space. We show the dependence on scale cuts when fitting the different redshift bins using the ShapeFit and Full-Modeling methods. We test the ability to jointly fit data from multiple redshift bins as well as joint analysis of the pre-reconstruction power spectrum with the post-reconstruction BAO correlation function signal. We further demonstrate the behavior of the model when opening up the parameter space beyond ΛCDM and also when combining likelihoods with external datasets, namely the Planck CMB priors. Finally, we describe different parametrization options for the galaxy bias, counterterm, and stochastic parameters, and employ the halo model in order to physically motivate suitable priors that are necessary to ensure the stability of the perturbation theory.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Sensitivity of modelled passive margin stratigraphy to variations in sea level, sediment supply and subsidence

Abstract We produced a 10 Myr synthetic stratigraphic section using a forward stratigraphic model that generates marine deltaic stratigraphy over geological timescales. We recursively fit the model using a Bayesian inversion algorithm to test: (1) if it could be accurately reconstructed; (2) if the parameters used to create it could be recovered; and (3) the sensitivity of the model output to given model parameters and the attendant physical processes. The original synthetic stratigraphic section was produced with cyclical sea‐level variations of 40 and 30 m with 2.4 and 10 Myr periods respectively. Sediment was also supplied cyclically, in 2.4 and 10 Myr cycles with amplitudes of 30 and 80 tons/100 kyr, respectively, varying from a mean of 232 tons/100 kyr. Parameter values were sampled to fit the model using a Markov chain Monte Carlo algorithm, resulting in a ±5 m (1σ) variation between the experimental output and the original. Sea level varied by ±7 m (1σ) within the posterior distribution of parameters. As a result, both the 10 Myr and 2.4 Myr sea‐level cycles could be extracted from the original output. The variation in sediment supply was approximately ±38 tons/100 kyr (1σ) and, as a result, only the larger long‐term supply variations could be accurately recovered in refitting the model. The variation in thermal, flexural and total subsidence across those parameter sets is less than ±10 m (1σ). The original section experienced 150 m of total subsidence at the depocentre. Our results demonstrate the distinct and interpretable imprint of sea level and subsidence on continental margin stratigraphy can be quantified. Moreover, we conclude that sea‐level change produces a defined effect on the geometries of stratigraphic architecture, and that techniques applied for the purpose of delineating sea‐level variation from continental margin strata have a well‐founded conceptual basis.

Schmelz, William J.↗
DOE OSTI2024

Constrained curve fitting for semi-parametric models with radial basis function networks

Common to many analysis pipelines in lattice gauge theory and the broader scientific discipline is the need to fit a semi-parametric model to data. We propose a fit method that utilizes a radial basis function network to approximate the non-parametric component of such models. The approximate parametric model is fit to data using the basin hopping global optimization algorithm. Parameter constraints are enforced through Gaussian priors. The viability of our method is tested by examining its use in a finite-size scaling analysis of the $q$-state Potts model and $p$-state clock model with $q=2,3$ and $p=4,\infty$.

Peterson, Curtis T.↗
DOE OSTI2025

Assessing methods in fusion and fitting for time series construction in remote sensing-based earth observations

This study evaluates the comparative performance of spatiotemporal fusion and time-series fitting methods for constructing high-spatiotemporal-resolution remote sensing time-series data. Due to in-class similarity of fusion methods and fitting methods, we employ the Fit-FC (Fitting, spatial Filtering, and residual Compensation) model as a representative fusion method and the linear harmonic fitting model as a representative fitting method. Both Fit-FC and the linear harmonic fitting are widely used for high-spatiotemporal-resolution time-series data construction, and we modify the original Fit-FC model to enable automatic time-series fusion. To ensure data representativeness, we use 3 years (2019–2021) of Harmonized Landsat and Sentinel-2 surface reflectance datasets and Terra MCD43A4 products. Eight experimental regions are selected worldwide to guarantee generalization of the comparative performance between fusion and fitting methods, covering diverse land-use types (cropland, developed land, forest, and grassland) and varying climatological conditions. Time-series of NDVI and surface reflectance are analyzed under both actual observations and simulated data-missing scenarios. The constructed time-series data reveals that (1) the modified Fit-FC and linear harmonic fitting model achieve excellent performance in constructing high-resolution time-series images; (2) the fusion method outperforms the fitting method in constructing time-series of NDVI and surface reflectance images in cropland-, forest-, and grassland-dominated regions; (3) both methods achieve comparable performance in developed-dominated regions; (4) the fusion method is more robust to missing data, and better captures abrupt phenological transitions under conditions of continuous missing data; (5) the fitting method is computationally more efficient, making it suitable for large-scale time-series image reconstruction. This study provides valuable insights for selecting optimal strategies to generate high-resolution time-series images across diverse application scenarios and lays a foundation for extensions to other vegetation indices or land surface variables.

54 ENVIRONMENTAL SCIENCES↗