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 793 records · Page 44

In Situ Inference for Earth System Predictability

Focal Area: Focal Area 3: Insight gleaned from complex simulated data using AI, big data analytics, and other advanced methods, including explainable AI and physics- or knowledge-guided AI. Science Challenge: An understanding of future evolution in precipitation extremes is critical to numerous DOE mission questions. Extreme events are by nature short time-scale events that are difficult to diagnose in available model data. Accurate modeling of extreme events necessarily requires high spatial resolution at the storm scale locally. However, the environment in which storms grow is dependent on global, remote, processes. These complex spatiotemporal relationships are impossible to diagnose at resolutions required to accurately model storms responsible for extreme precipitation. At exascale, climate simulations will produce results at fine enough resolution to investigate these relationships. However, the resulting data from these simulations will be far too large to save for post-simulation analysis. We advocate for fitting statistical models inside the simulations as they run, a context known as in situ, which will facilitate scientific investigations using the full fine-scale data stream. Figure 1 shows an example of the type of model we could consider, a Bayesian hierarchical spatial regression model. Precipitation extremes at each grid cell are modeled using extreme value distributions. Since extremes are rare, fitting models to individual grid cells can result in high variance and poor estimates. Instead, the model can be made more robust by smoothing the parameters of the extreme value model across space. Additionally, the parameters themselves can be functionally linked to other variables elsewhere in the simulation. Thus, we can use the fine-scale data to build more robust models for extremes that link extreme behavior to other climate patterns.

54 ENVIRONMENTAL SCIENCES↗

A direct measurement of the mean occupation function of quasars: Breaking degeneracies between halo occupation distribution models

Recent work on quasar clustering suggests a degeneracy in the halo occupation distribution constrained from two-point correlation functions. To break this degeneracy, we make the first empirical measurement of the mean occupation function (MOF) of quasars at z ∼ 0.2 by matching quasar positions with groups and clusters identified in the MaxBCG sample. We fit two models to the MOF, a power law and a four-parameter model. The number distribution of quasars in host halos is close to Poisson, and the slopes of the MOF obtained from our best-fit models (for the power-law case) favor an MOF that monotonically increases with halo mass. The best-fit slopes are 0.53 ± 0.04 and 1.03 ± 1.12 for the power-law model and the four-parameter model, respectively. We measure the radial distribution of quasars within dark matter halos and find it to be adequately described by a power law with a slope –2.3 ± 0.4. We measure the conditional luminosity function (CLF) of quasars and show that there is no evidence that quasar luminosity depends on host halo mass, similar to the inferences drawn from clustering measurements. We also measure the conditional black hole mass function (CMF) of our quasars. Although the results are consistent with no dependence on halo mass, we observe a slight indication of downsizing of the black hole mass function. The lack of halo mass dependence in the CLF and CMF shows that quasars residing in galaxy clusters have characteristic luminosity and black hole mass scales.

79 ASTRONOMY AND ASTROPHYSICS↗

Precision measurements of EFT parameters and BAO peak shifts for the Lyman- α forest

We present precision measurements of the bias parameters of the one-loop power spectrum model of the Lyman- α (Ly- α ) forest, derived within the effective field theory (EFT) of large-scale structure. We fit our model to the three-dimensional flux power spectrum measured from the ACCEL 2 hydrodynamic simulations. The EFT model fits the data with an accuracy of below 2% up to k = 2 h Mpc − 1 . Further, we analytically derive how nonlinearities in the three-dimensional clustering of the Ly- α forest introduce biases in measurements of the baryon acoustic oscillations (BAOs) scaling parameters in radial and transverse directions. From our EFT parameter measurements, we obtain a theoretical error budget of Δ α ∥ = − 0.2 % ( Δ α ⊥ = − 0.3 % ) for the radial (transverse) parameters at redshift z = 2.0 . This corresponds to a shift of − 0.3 % (0.1%) for the isotropic (anisotropic) distance measurements. We provide an estimate for the shift of the BAO peak for Ly- α -quasar cross-correlation measurements assuming analytical and simulation-based scaling relations for the nonlinear quasar bias parameters resulting in a shift of − 0.2 % ( − 0.1 % ) for the radial (transverse) dilation parameters, respectively. This analysis emphasizes the robustness of Ly- α forest BAO measurements to the theory modeling. We provide informative priors and an error budget for measuring the BAO feature—a key science driver of the currently observing Dark Energy Spectroscopic Instrument (DESI). Our work paves the way for full-shape cosmological analyses of Ly- α forest data from DESI and upcoming surveys such as the Prime Focus Spectrograph, WEAVE-QSO, and 4MOST. Published by the American Physical Society 2025

de Belsunce, Roger (ORCID:0000000336604028)↗

Probing the Diffuse Ly α Emission on Cosmological Scales: Ly α Emission Intensity Mapping Using the Complete SDSS-IV eBOSS

Based on Sloan Digital Sky Survey Data Release 16, we have detected the large-scale structure of Lyα emission in the universe at redshifts z = 2–3.5 by cross-correlating quasar positions and Lyα emission imprinted in the residual spectra of luminous red galaxies. We apply an analytical model to fit the corresponding Lyα surface brightness profile and multipoles of the redshift-space quasar–Lyα emission cross-correlation function. The model suggests an average cosmic Lyα luminosity density of $6.6^{+3.3}_{-3.1}$ x 10 40 erg s -1 cMpc -3 , a ~2σ detection with a median value about 8–9 times those estimated from deep narrowband surveys of Lyα emitters at similar redshifts. Although the low signal-to-noise ratio prevents us from a significant detection of the Lyα forest–Lyα emission cross-correlation, the measurement is consistent with the prediction of our best-fit model from quasar–Lyα emission cross-correlation within current uncertainties. We rule out the scenario where the Lyα photons mainly originate from quasars. We find that Lyα emission from star-forming galaxies, including contributions from that concentrated around the galaxy centers and that in diffuse Lyα-emitting halos, is able to explain the bulk of the Lyα luminosity density inferred from our measurements. Ongoing and future surveys can further improve the measurements and advance our understanding of the cosmic Lyα emission field.

79 ASTRONOMY AND ASTROPHYSICS↗

Fitting local, low-dimensional parameterizations of optical turbulence modeled from optimal transport velocity vectors

This work exploits a connection between optimal transport theory and the physics of image propagation to yield a locally low-dimensional model of turbulence-corrupted imagery. Optimal transport produces an invertible, pixel-wise linear trajectories to approximate the globally nonlinear turbulence between a clean and turbulence corrupted image pair. We use the low-dimensional model to fit subsets of the optimal transport vector fields and stitch the local models into a surrogate for the global map to be used for image cleaning. Experiments are performed on laboratory generated data of beam propagation using different values of the Fried parameter (a scale measuring turbulence coherence) as well as a toy data set. The results suggest this is a fruitful direction, and first step, towards using multiple realizations of turbulence corrupted images to learn a blind surrogate for the optimal transport vector field for image cleaning.

47 OTHER INSTRUMENTATION↗

The Disk Physical Conditions of 48 Persei Constrained by Contemporaneous Hα Spectroscopy and Interferometry

We utilize a multi-step modeling process to produce synthetic interferometric and spectroscopic observables, which are then compared to their observed counterparts. Our extensive set of interferometric observations of the Be star 48 Per, totaling 291 data points, were obtained at the Navy Precision Optical Interferometer from 2006 November 07 to 23. Our models were further constrained by comparison with contemporaneous Hα line spectroscopy obtained at the John S. Hall Telescope at the Lowell Observatory on 2006 November 1. Theoretical spectral energy distributions, SEDs, for 48 Per were confirmed by comparison with observations over a wavelength regime of 0.4–60 μm from Touhami et al. and Vieira et al. Our best-fitting combined model from Hα spectroscopy, Hα interferometry, and SED fitting has a power-law density fall off, n, of 2.3 and an initial density at the stellar surface of ρ{sub 0}=1.0×10{sup −11} g cm{sup −3} with an inclination constrained by Hα spectroscopy and interferometry of 45{sup ∘}±5{sup ∘}. The position angle for the system, measured east from north, is 114° ± 18°. Our best-fit model shows that the disk emission originates in a moderately large disk with a radius of 25 R {sub *}, which is consistent with a disk mass of approximately 5 × 10{sup 24} g or 3 × 10{sup −10} M {sub *}. Finally, we compare our results with previous studies of 48 Per by Quirrenbach et al. and find agreement, whereas our disk size does not agree with Delaa et al., based on a much smaller visibility set.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling HCV cure after an ultra-short duration of therapy with direct acting agents

In cases of sustained-virological response (SVR or cure) after an ultra-short duration (≤ 27 days) of direct-acting antiviral (DAA)-based therapy, despite HCV being detected at end of treatment (EOT), have been reported. Established HCV mathematical models that predict the treatment duration required to achieve cure do not take into account the possibility that the infectivity of virus produced during treatment might be reduced. The aim of this study was to develop a new mathematical model that considers the fundamental and critical concept that HCV RNA in serum represents both infectious virus (V i ) and non-infectious virus (V ni ) in order to explain the observation of cure with ultrashort DAA therapy. Established HCV models were compared to the new mathematical model to retrospectively explain cure in 2 patients who achieved cure after 24 or 27 days of paritaprevir, ombitasvir, dasabuvir, ritonavir and ribavirin or sofosbuvir plus ribavirin, respectively. Fitting established models with measured longitudinal HCV viral loads indicated that in both cases, cure would not have been expected without an additional 3–6 weeks of therapy after the actual EOT. In contrast, the new model fits the observed outcome by considering that in addition to blocking V i and V ni production (ε~0.998), these DAA + ribavirin treatments further enhanced the ratio of V ni to V i , thus increasing the log (V ni /V i ) from 1 at pretreatment to 6 by EOT, which led to <1 infectious-virus particle in the extracellular body fluid (i.e., cure) prior to EOT. This new model can explain cure after short duration of DAA + ribavirin therapy by suggesting that a minimum 6-fold increase of log (V ni /V i ) results from drug-induced enhancement of the V ni /V i .

59 BASIC BIOLOGICAL SCIENCES↗

Correcting for telluric absorption: Methods, case studies, and release of the TelFit code

Ground-based astronomical spectra are contaminated by the Earth's atmosphere to varying degrees in all spectral regions. We present a Python code that can accurately fit a model to the telluric absorption spectrum present in astronomical data, with residuals of ∼3%-5% of the continuum for moderately strong lines. We demonstrate the quality of the correction by fitting the telluric spectrum in a nearly featureless A0V star, HIP 20264, as well as to a series of dwarf M star spectra near the 819 nm sodium doublet. We directly compare the results to an empirical telluric correction of HIP 20264 and find that our model-fitting procedure is at least as good and sometimes more accurate. The telluric correction code, which we make freely available to the astronomical community, can be used as a replacement for telluric standard star observations for many purposes.

79 ASTRONOMY AND ASTROPHYSICS↗

Bayesian characterization and forward propagation of the uncertainty in thermodynamic models for redox active materials

Demonstrating the thermodynamic efficiency of hydrogen conversion processes using various materials is a critical step in developing new technologies for storing concentrated solar energy, and is largely accomplished by using a thermodynamic model derived from experimental data. A main goal of this project is to calculate the uncertainty of the thermodynamic efficiency by calculating the uncertainty of the components that feed into the efficiency. Many different models and data sets were used to test the workflow. First, the models were fit to the data using a Bayesian Inference and a method called Markov Chain Monte Carlo (MCMC), which found the maximum a priori parameters, and a posterior probability distribution of the parameters. Next, the different models were compared to each other using model evidence values. It was found that for cleaner data sets, overfitting had not yet been reached, and the most complicated model was ideal, but on the noisier data sets, the less complex models were favored because the more complicated models resulted in overfitting. Next, forward propagation was used to calculate the enthalpy change and its associated uncertainty. A few variations on the models were tried, such as fitting in a different variable, producing negligible or negative effects on the fits of the models. Thus, the original models were used. A sensitivity analysis was performed, and used to calculate the model error. On the cleaner data sets, there was very minimal experimental noise, and thus, all resulting error was from the model. With consideration of the model error, the models fit the data very well, and the simpler model had a high model error, as expected. All these components will then be used to calculate the thermodynamic efficiency of the different materials.

14 SOLAR ENERGY↗

Long-time mean-square displacements in proteins

We propose a method for obtaining the intrinsic, long-time mean square displacement (MSD) of atoms and molecules in proteins from finite-time molecular dynamics (MD) simulations. Typical data from simulations are limited to times of 1 to 10 ns, and over this time period the calculated MSD continues to increase without a clear limiting value. The proposed method consists of fitting a model to MD simulation-derived values of the incoherent intermediate neutron scattering function, I inc (Q,t), for finite times. The infinite-time MSD, r 2 , appears as a parameter in the model and is determined by fits of the model to the finite-time I inc (Q,t). Specifically, the r 2 is defined in the usual way in terms of the Debye-Waller factor as I(Q,t = ∞) = exp( –Q 2 r 2 /3). The method is illustrated by obtaining the intrinsic MSD r 2 of hydrated lysozyme powder (h = 0.4 g water/g protein) over a wide temperature range. The intrinsic r 2 obtained from data out to 1 and to 10 ns is found to be the same. In conclusion, the intrinsic r 2 is approximately twice the value of the MSD that is reached in simulations after times of 1 ns which correspond to those observed using neutron instruments that have an energy resolution width of 1 μeV.

59 BASIC BIOLOGICAL SCIENCES↗

HighDimMixedModels.jl: Robust high-dimensional mixed-effects models across omics data

High-dimensional mixed-effects models are an increasingly important form of regression in which the number of covariates rivals or exceeds the number of samples, which are collected in groups or clusters. The penalized likelihood approach to fitting these models relies on a coordinate descent algorithm that lacks guarantees of convergence to a global optimum. Here, we empirically study the behavior of this algorithm on simulated and real examples of three types of data that are common in modern biology: transcriptome, genome-wide association, and microbiome data. Our simulations provide new insights into the algorithm’s behavior in these settings, and, comparing the performance of two popular penalties, we demonstrate that the smoothly clipped absolute deviation (SCAD) penalty consistently outperforms the least absolute shrinkage and selection operator (LASSO) penalty in terms of both variable selection and estimation accuracy across omics data. To empower researchers in biology and other fields to fit models with the SCAD penalty, we implement the algorithm in a Julia package, HighDimMixedModels.jl .

Gorstein, Evan↗

Improving neutrino-nuclei interaction models: Recommendations and case studies on Peelle’s Pertinent Puzzle

Improving the modeling of neutrino-nuclei interactions using data-driven methods is crucial for high-precision neutrino oscillation experiments. This paper investigates Peelle’s Pertinent Puzzle (PPP) in the context of neutrino measurements, a longstanding challenge to fitting theoretical models to experimental data. Inconsistencies in data-model comparisons hinder efforts to enhance the accuracy and reliability of model predictions. We analyze various sources contributing to these inconsistencies and propose strategies to address them, supported by practical case studies. We advocate for incorporating model fitting exercises as a standard practice in cross section publications to enhance the robustness of results. We use a common analysis framework to explore PPP-related challenges with MicroBooNE and T2K data in an unified manner. Our findings offer valuable insights for improving the accuracy and reliability of neutrino-nuclei interaction models, particularly by systematically tuning models using data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

TU-D-207B-02: Delta-Radiomics: The Prognostic Value of Therapy-Induced Changes in Radiomics Features for Stage III Non-Small Cell Lung Cancer Patients

Purpose: To determine how radiomics features change during radiation therapy and whether those changes (delta-radiomics features) can improve prognostic models built with clinical factors. Methods: 62 radiomics features, including histogram, co-occurrence, run-length, gray-tone difference, and shape features, were calculated from pretreatment and weekly intra-treatment CTs for 107 stage III NSCLC patients (5–9 images per patient). Image preprocessing for each feature was determined using the set of pretreatment images: bit-depth resample and/or a smoothing filter were tested for their impact on volume-correlation and significance of each feature in univariate cox regression models to maximize their information content. Next, the optimized features were calculated from the intratreatment images and tested in linear mixed-effects models to determine which features changed significantly with dose-fraction. The slopes in these significant features were defined as delta-radiomics features. To test their prognostic potential multivariate cox regression models were fitted, first using only clinical features and then clinical+delta-radiomics features for overall-survival, local-recurrence, and distant-metastases. Leave-one-out cross validation was used for model-fitting and patient predictions. Concordance indices(c-index) and p-values for the log-rank test with patients stratified at the median were calculated. Results: Approximately one-half of the 62 optimized features required no preprocessing, one-fourth required smoothing, and one-fourth required smoothing and resampling. From these, 54 changed significantly during treatment. For overall-survival, the c-index improved from 0.52 for clinical factors alone to 0.62 for clinical+delta-radiomics features. For distant-metastases, the c-index improved from 0.53 to 0.58, while for local-recurrence it did not improve. Patient stratification significantly improved (p-value<0.05) for overallsurvival and distant-metastases when delta-radiomics features were included. The delta-radiomics versions of autocorrelation, kurtosis, and compactness were selected most frequently in leave-one-out iterations. Conclusion: Weekly changes in radiomics features can potentially be used to evaluate treatment response and predict patient outcomes. High-risk patients could be recommended for dose escalation or consolidation chemotherapy. This project was funded in part by grants from the National Cancer Institute (NCI) and the Cancer Prevention Research Institute of Texas (CPRIT).

60 APPLIED LIFE SCIENCES↗

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↗

Cultural Artifact Detection in Long Wave Infrared Imagery.

Detection of cultural artifacts from airborne remotely sensed data is an important task in the context of on-site inspections. Airborne artifact detection can reduce the size of the search area the ground based inspection team must visit, thereby improving the efficiency of the inspection process. This report details two algorithms for detection of cultural artifacts in aerial long wave infrared imagery. The first algorithm creates an explicit model for cultural artifacts, and finds data that fits the model. The second algorithm creates a model of the background and finds data that does not fit the model. Both algorithms are applied to orthomosaic imagery generated as part of the MSFE13 data collection campaign under the spectral technology evaluation project.

47 OTHER INSTRUMENTATION↗

WE-AB-207B-10: On Spinal Nerve Toxicity from Single-Session SAbR in Pigs and the Translation of Small Animal NTCP Models

Purpose: The incidence of peripheral neuropathy has risen with increased utilization of SAbR. There is no consensus regarding the dose-tolerance of the peripheral nervous system. In 2015, we commenced an investigation to test the hypotheses that single-session irradiation to the pig spinal nerves exhibit a similar dose-tolerance as that of the spinal cord and that a dose-length effect exists. This work evaluates the direct application of small animal NTCP models to both large animal spinal cord and preliminary peripheral nerve data. Methods: To date, 16 of 25 Yucatan minipigs have received single-session SAbR to a 1.5cm length and 4 of 25 have received irradiation to a 0.5cm length of left-sided C6-C8 spinal nerves. Toxicity related gait change has been observed in 13 animals (9 from the long length group and 4 from the short). This preliminary data is overlaid on several dose-response models which have been fit to rodent spinal cord tolerance experiments. Model parameters define a toxicity profile between a completely serial or parallel behaving organ. Adequacy of model application, including how length effects are handled, to published minipig spinal cord dose-response data and to preliminary peripheral nerve response data was evaluated through residual analysis. Results: No rodent-derived dose-response models were directly applicable to all pig data for the different lengths irradiated. Several models fit the long-length irradiated spinal cord data well, with the more serial-like models fitting best. Preliminary data on the short-length irradiation suggests no length effect exists, disproving our hypothesis. Conclusion: Direct application of small-animal NTCP models to pig data suggests dose-length effect predictions from small animal data may not translate clinically. However, the small animal models used have not considered dose heterogeneity and it is expected that including the low-to-mid dose levels in the penumbral region will improve this match. This work was funded by the Cancer Prevention Research Institute of Texas (CPRIT).

60 APPLIED LIFE SCIENCES↗

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↗