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 1,459 records · Page 81

CD8 + lymphocyte control of SIV infection during antiretroviral therapy

CD8 + lymphocytes play an important role in suppressing in vivo viral replication in HIV infection. However, both the extent to which and the mechanisms by which CD8 + lymphocytes contribute to viral control are not completely understood. A recent experiment depleted CD8 + lymphocytes in simian immunodeficiency virus (SIV)-infected rhesus macaques (RMs) on antiretroviral treatment (ART) to study the role of CD8 + lymphocytes. CD8 + lymphocytes depletion resulted in temporary plasma viremia in all studied RMs. Viral control was restored when CD8 + lymphocytes repopulated. We developed a viral dynamic model to fit the viral load (VL) data from the CD8 depletion experiment. We explicitly modeled the dynamics of the latent reservoir and the SIV-specific effector cell population including their exhaustion and their potential cytolytic and noncytolytic functions. We found that the latent reservoir significantly contributes to the size of the peak VL after CD8 depletion, while drug efficacy plays a lesser role. Our model suggests that the overall CD8 + lymphocyte cytolytic killing rate is dynamically changing depending on the levels of antigen-induced effector cell activation and exhaustion. Based on estimated parameters, our model suggests that before ART or without ART the overall CD8 cytolytic killing rate is small due to exhaustion. However, after the start of ART, the overall CD8 cytolytic killing rate increases due to an expansion of SIV-specific CD8 effector cells. Further, we estimate that the cytolytic killing rate can be significantly larger than the cytopathic death rate in some animals during the second phase of ART-induced viral decay. Lastly, our model provides a new explanation for the puzzling findings by Klatt et al. and Wong et al. that CD8 depletion done immediately before ART has no noticeable effect on the first phase viral decay slope seen after ART initiation Overall, by incorporating effector cells and their exhaustion, our model can explain the effects of CD8 depletion on VL during ART, reveals a detailed dynamic role of CD8 + lymphocytes in controlling viral infection, and provides a unified explanation for CD8 depletion experimental data.

59 BASIC BIOLOGICAL SCIENCES↗

A Model-Based Approach to Scintillator/Photomultiplier System Characterization

Scintillator/photomultiplier systems are prevalent in many instruments to measure a variety on nuclear events. When gamma energy from a nuclear event irradiates the scintillator, event radiation interacts with the scintillator material generating photons which are detected by the photomultiplier tube (PMT). The PMT photoelectrons or photocathode current is amplified in the PMT by a number of dynode stages producing a total charge out for a given radiation flux into scintillator. The resulting current is converted to a voltage when passed through a typical load impedance. This voltage is attenuated and provided as input to a digitizer through a variety of signal conditioners (filters, amplifiers, receivers) before analog-to-digital conversion. In this report, we concentrate on the development of analytic models for the scintillator/photomultiplier system and develop both a (1) statistical simulator; and (2) a model-based signal processor to extract the desired signal information while rejecting the noise and uncertainty. We develop the model analytically, transform it to the model-based framework (state-space), and then use actual experimental data to fit the model and apply it to our processing problem.

42 ENGINEERING↗

Deep galaxy counts in the K band with the Kech telescope

We present deep galaxy counts in the K (lambda 2.2 micrometer) band, obtained at the W. M. Kech 10 m telescope. The data reach limiting magnitudes K approximately 24 mag, about 5 times deeper than the deepest published K-band images to date. The counts are performed in three small (approximately 1 min), widely separated high-latitude fields. Extensive Monte Carlo tests were used to derive the comleteness corrections and minimize photometric biases. The counts continue to rise, with no sign of a turnover, down to the limits of our data, with the logarithmic slope of d log N/dm = 0.315 +/- 0.02 between K = 20 and 24 mag. This implies a cumulative surface density of approximately 5 x 10(exp 5) galaxies/sq deg, or approximately 2 x 10(exp 10) over the entire sky, down to K = 24 mag. Our counts are in good agreement with, although slightly lower than, those from the Hawaii Deep Survey by Cowie and collaborators; the discrepancies may be due to the small differences in the aperture corrections. We compare our counts with some of the available theoretical predictions. The data do not require models with a high value of Omega(sub 0), but can be well fitted by models with no (or little) evolution, and cosmologies with a low value of Omega(sub 0). Given the uncertainties in the models, it may be premature to put useful constrains on the value of Omega(sub 0) from the counts alone. Optical-to-IR colors are computed, using CCD data obtaind previously at Palomar. We find a few red galaxies with (r-K) approximately greater than 5 mag, or (i-K) approximately greater than 5 mag; these may be ellipticals at z approximately 1. While the redshift distribution of galaxies in our counts is still unknown, the flux limits reached would allow us to detect unobscured L(sub *) galaxies out to substantial redshifts (z greater than 3?).

Djorgovski, S.↗

Analytic Guided-Search Model of Human Performance Accuracy in Target- Localization Search Tasks

Current models of human visual search have extended the traditional serial/parallel search dichotomy. Two successful models for predicting human visual search are the Guided Search model and the Signal Detection Theory model. Although these models are inherently different, it has been difficult to compare them because the Guided Search model is designed to predict response time, while Signal Detection Theory models are designed to predict performance accuracy. Moreover, current implementations of the Guided Search model require the use of Monte-Carlo simulations, a method that makes fitting the model's performance quantitatively to human data more computationally time consuming. We have extended the Guided Search model to predict human accuracy in target-localization search tasks. We have also developed analytic expressions that simplify simulation of the model to the evaluation of a small set of equations using only three free parameters. This new implementation and extension of the Guided Search model will enable direct quantitative comparisons with human performance in target-localization search experiments and with the predictions of Signal Detection Theory and other search accuracy models.

Eckstein, Miguel P.↗

Shock-Stationary Application of Pseudoshock Models During High-Amplitude Combustion-Driven Unsteadiness

The isolator pseudo-shock provides necessary compression within a dual-mode scramjet engine and buffers the engine system against unstart. Quasi-1D flux-conserved models are the state-of-the-art reduced-order model for optimization and online control of dual-mode scramjet engines. The stability and efficacy of this modeling approach is evaluated against data from a combustor-driven direct-connect experiment. The experiment exhibited strong combustor-driven unsteadiness that produced upstream propagating weak shocks into the isolator, interacting with the pseudo-shock. While this configuration resulted in unsteadiness that is atypical of standard operation, the experiment provided an opportunity to evaluate the modeling techniques in highly transient states. Such transients could occur during maneuvering or result from unexpected combustor events. A flexible quasi-1D formulation, the Fievet flux-conserved model is fit using Bayesian inference in the laboratory and shock-stationary reference frames. Model performance is analyzed using the Bayesian posteriors and model evaluations over the measured shock train speed range. It is concluded that to produce consistent isolator pressure profile estimates in this unsteady environment, the model must be implemented in a shock-stationary reference frame. Implementing this conclusion in model-based engine controllers may reduce needed unstart safety margins and increase maximum performance.

Bayesian↗

Measurements of the Sunyaev–Zel’Dovich effect in MACS J0647.7+7015 and MACS J1206.2–0847 at high angular resolution with MUSTANG

We present high resolution (9″) imaging of the Sunyaev–Zel’dovich Effect (SZE) toward two massive galaxy clusters, MACS J0647.7+7015 (z = 0.591) and MACS J1206.2–0847 (z = 0.439). We compare these 90 GHz measurements, taken with the Multiplexed Squid/TES Array at Ninety Gigahertz (MUSTANG ) receiver on the Green Bank Telescope, with generalized Navarro–Frenk–White (gNFW) models derived from Bolocam 140 GHz SZE data as well as maps of the thermal gas derived from Chandra X-ray observations. We adopt a serial-fitting approach, in which gNFW models are first fit to the Bolocam data and then compared to the MUSTANG data to determine an overall best-fit model. For MACS J0647.7+7015, we find a gNFW profile with core slope parameter γ = 0.9 fits the MUSTANG image with χ{sub red}{sup 2}=1.005 and probability to exceed (PTE) = 0.34. For MACS J1206.2–0847, we find γ=0.7, χ{sub red}{sup 2}=0.993, and PTE = 0.70. In addition, we find a significant (>3σ) residual SZE feature in MACS J1206.2–0847 coincident with a group of galaxies identified in Very Large Telescope data and filamentary structure found in a weak-lensing mass reconstruction. We suggest the detected sub-structure may be the SZE decrement from a low mass foreground group or an infalling group. Giant Metrewave Radio Telescope measurements at 610 MHz reveal diffuse extended radio emission to the west, which we posit is either an active galactic nucleus-driven radio lobe, a bubble expanding away from disturbed gas associated with the SZE signal, or a bubble detached and perhaps re-accelerated by sloshing within the cluster. Using the spectroscopic redshifts available, we find evidence for a foreground (z = 0.423) or infalling group, coincident with the residual SZE feature.

79 ASTRONOMY AND ASTROPHYSICS↗

The Atacama Cosmology Telescope: likelihood for small-scale CMB data

The Atacama Cosmology Telescope has measured the angular power spectra of microwave fluctuations to arcminute scales at frequencies of 148 and 218 GHz, from three seasons of data. At small scales the fluctuations in the primordial Cosmic Microwave Background (CMB) become increasingly obscured by extragalactic foregounds and secondary CMB signals. We present results from a nine-parameter model describing these secondary effects, including the thermal and kinematic Sunyaev-Zel'dovich (tSZ and kSZ) power; the clustered and Poisson-like power from Cosmic Infrared Background (CIB) sources, and their frequency scaling; the tSZ-CIB correlation coefficient; the extragalactic radio source power; and thermal dust emission from Galactic cirrus in two different regions of the sky. In order to extract cosmological parameters, we describe a likelihood function for the ACT data, fitting this model to the multi-frequency spectra in the multipole range 500 < l < 10000. We extend the likelihood to include spectra from the South Pole Telescope at frequencies of 95, 150, and 220 GHz. Accounting for different radio source levels and Galactic cirrus emission, the same model provides an excellent fit to both datasets simultaneously, with χ{sup 2}/dof= 675/697 for ACT, and 96/107 for SPT. We then use the multi-frequency likelihood to estimate the CMB power spectrum from ACT in bandpowers, marginalizing over the secondary parameters. This provides a simplified 'CMB-only' likelihood in the range 500 < l < 3500 for use in cosmological parameter estimation.

79 ASTRONOMY AND ASTROPHYSICS↗

The Atacama Cosmology Telescope: Likelihood for Small-Scale CMB Data

The Atacama Cosmology Telescope has measured the angular power spectra of microwave fluctuations to arcminute scales at frequencies of 148 and 218 GHz, from three seasons of data. At small scales the fluctuations in the primordial Cosmic Microwave Background (CMB) become increasingly obscured by extragalactic foregounds and secondary CMB signals. We present results from a nine-parameter model describing these secondary effects, including the thermal and kinematic Sunyaev-Zel'dovich (tSZ and kSZ) power; the clustered and Poisson-like power from Cosmic Infrared Background (CIB) sources, and their frequency scaling; the tSZ-CIB correlation coefficient; the extragalactic radio source power; and thermal dust emission from Galactic cirrus in two different regions of the sky. In order to extract cosmological parameters, we describe a likelihood function for the ACT data, fitting this model to the multi-frequency spectra in the multipole range 500 < l < 10000. We extend the likelihood to include spectra from the South Pole Telescope at frequencies of 95, 150, and 220 GHz. Accounting for different radio source levels and Galactic cirrus emission, the same model provides an excellent fit to both datasets simultaneously, with χ2/dof= 675/697 for ACT, and 96/107 for SPT. We then use the multi-frequency likelihood to estimate the CMB power spectrum from ACT in bandpowers, marginalizing over the secondary parameters. This provides a simplified 'CMB-only' likelihood in the range 500 < l < 3500 for use in cosmological parameter estimation

Cosmology↗

Does the first chaotic inflation model in supergravity provide the best fit to the Planck data?

I describe the first model of chaotic inflation in supergravity, which was proposed by Goncharov and the present author in 1983. The inflaton potential of this model has a plateau-type behavior V{sub 0}(1−(8/3) e{sup −√6|ϕ|}) at large values of the inflaton field. This model predicts n{sub s}=1−(2/N)≈0.967 and r=(4/(3N{sup 2}))≈4×10{sup −4}, in good agreement with the Planck data. I propose a slight generalization of this model, which allows to describe not only inflation but also dark energy and supersymmetry breaking.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A viscous blast-wave model for high energy heavy-ion collisions

Employing a viscosity-based survival scale for initial geometrical perturbations formed in relativistic heavy-ion collisions, we model the radial flow velocity at freeze-out. Subsequently, we use the Cooper-Frye freeze-out prescription, with viscous corrections to the distribution function, to extract the transverse momentum dependence of particle yields and flow harmonics. We fit the model parameters for central collisions, by fitting the spectra of identified particles at the Large Hadron Collider (LHC), and estimate them for other centralities using simple hydrodynamic relations. We use the results of Monte Carlo Glauber model for initial eccentricities. We demonstrate that this improved viscous blast-wave model leads to good agreement with transverse momentum distribution of elliptic and triangular flow for all centralities and estimate the shear viscosity to entropy density ratio η/s ≃ 0.24 at the LHC.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Full Modeling and parameter compression methods in configuration space for DESI 2024 and beyond

In the contemporary era of high-precision spectroscopic surveys, led by projects like DESI, there is an increasing demand for optimizing the extraction of cosmological information from clustering data. This work conducts a thorough comparison of various methodologies for modeling the full shape of the two-point statistics in configuration space. We investigate the performance of both direct fits (Full Modeling) and the parameter compression approaches (ShapeFit and Standard). We utilize the ABACUS-SUMMIT simulations, tailored to exceed DESI's precision requirements. Particularly, we fit the two-point statistics of three distinct tracers (LRG, ELG, and QSO), by employing a Gaussian Streaming Model in tandem with Convolution Lagrangian Perturbation Theory and Effective Field Theory. We explore methodological setup variations, including the range of scales, the set of galaxy bias parameters, the inclusion of the hexadecapole, as well as model extensions encompassing varying ns and allowing for w 0 w a CDM dark energy model. Throughout these varied explorations, while precision levels fluctuate and certain configurations exhibit tighter parameter constraints, our pipeline consistently recovers the parameter values of the mocks within 1σ in all cases for a 1-year DESI volume. Additionally, we compare the performance of configuration space analysis with its Fourier space counterpart using three models: PyBird, FOLPS and velocileptors, presented in companion papers. We find good agreement with the results from all these models.

79 ASTRONOMY AND ASTROPHYSICS↗

Observations of ultraviolet-bright stars in globular clusters with the Hopkins Ultraviolet Telescope

Two UV-bright stars in globular clusters, UV5 in NGC 1851 and vZ 1128 in M3, were observed with the Hopkins Ultraviolet Telescope (HUT) during the Astro-1 space shuttle mission in 1990 December. The stars' spectra show weak absorption features and no significant emission features other than well known geocoronal lines. Detailed fitting of Kurucz (1991) stellar atmosphere models using a chi(sup 2) minimization technique indicates T(sub eff) = 16 000 K, log g = 2.5, and abundance (-1.0) (logarithm of abundance of elements heavier than helium relative to solar) for UV5, and T(sub eff) = 35 000 K, log g = 4.0, and abundance (-3.5) for vZ 1128. The Kurucz model which best fits vZ 1128 overpredicts the flux in the region below approximately 1000 A, an effect seen in previous models of O-type stars. Our results are robust (to within approximately 1000 K) with respect to uncertainties in interstellar reddening and atomic and molecular hydrogen column densities. We do not see significant molecular hydrogen absorption, which might have indicated material in a circumstellar shell, in either star's spectrum. We estimate the stellar luminosities to be log L/solar luminosity = 3.33 +/- 0.15 for UV5 and log L/solar luminosity = 3.21 +/- 0.12 for vZ 1128. These atmospheric parameters place both stars on the Schoenberner (post-AGB) tracks, though the stellar masses derived from the best-fitting Kurucz models are somewhat less than those predicted by the Schoenberner models. Examination of individual absorption line strengths reveals no significant abundance anomalies in either star.

Van Dyke Dixon, W.↗

Color-magnitude distribution of face-on nearby galaxies in Sloan digital sky survey DR7

We have analyzed the distributions in the color-magnitude diagram (CMD) of a large sample of face-on galaxies to minimize the effect of dust extinctions on galaxy color. About 300,000 galaxies with log (a/b) < 0.2 and redshift z < 0.2 are selected from the Sloan Digital Sky Survey DR7 catalog. Two methods are employed to investigate the distributions of galaxies in the CMD, including one-dimensional (1D) Gaussian fitting to the distributions in individual magnitude bins and two-dimensional (2D) Gaussian mixture model (GMM) fitting to galaxies as a whole. We find that in the 1D fitting, two Gaussians are not enough to fit galaxies with the excess present between the blue cloud and the red sequence. The fitting to this excess defines the center of the green valley in the local universe to be (u – r){sub 0.1} = –0.121M {sub r,} 0{sub .1} – 0.061. The fraction of blue cloud and red sequence galaxies turns over around M {sub r,} {sub 0.1} ∼ –20.1 mag, corresponding to stellar mass of 3 × 10{sup 10} M {sub ☉}. For the 2D GMM fitting, a total of four Gaussians are required, one for the blue cloud, one for the red sequence, and the additional two for the green valley. The fact that two Gaussians are needed to describe the distributions of galaxies in the green valley is consistent with some models that argue for two different evolutionary paths from the blue cloud to the red sequence.

79 ASTRONOMY AND ASTROPHYSICS↗

A critical assessment of solutions to the galaxy diversity problem

Galactic rotation curves exhibit a diverse range of inner slopes. Observational data indicates that explaining this diversity may require a mechanism that correlates a galaxy's surface brightness with the central-most region of its dark matter halo. Here, in this work, we compare several concrete models that capture the relevant physics required to explain the galaxy diversity problem. We focus specifically on a Self-Interacting Dark Matter (SIDM) model with an isothermal core and two Cold Dark Matter (CDM) models with/without baryonic feedback. In contrast to the CDM case, the SIDM model can lead to the formation of an isothermal core in the halo, and is also mostly insensitive to baryonic feedback processes, which act on longer time-scales. Using rotation curves from 90 galaxies in the Spitzer Photometry & Accurate Rotation Curves (SPARC) catalog, we perform a comprehensive model comparison that addresses issues of statistical methodology from prior works. The best-fit halo models that we recover are consistent with standard CDM concentration-mass and abundance matching relations. We find that both the SIDM and feedback-affected CDM models are better than a CDM model with no feedback in explaining the rotation curves of low and high surface brightness galaxies in the sample. However, when compared to each other, there is no strong statistical preference for either the SIDM or the feedback-affected CDM halo model as the source of galaxy diversity in the SPARC catalog.

79 ASTRONOMY AND ASTROPHYSICS↗

Robustness of cosmic birefringence measurement against Galactic foreground emission and instrumental systematics

The polarization of the cosmic microwave background (CMB) can be used to search for parity-violating processes like that predicted by a Chern-Simons coupling to a light pseudoscalar field. Such an interaction rotates E modes into E modes in the observed CMB signal through an effect known as cosmic birefringence. Even though isotropic birefringence can be confused with the rotation produced by a miscalibration of the detectors' polarization angles, the degeneracy between both effects is broken when Galactic foreground emission is used as a calibrator. Here, in this work, we use realistic simulations of the High-Frequency Instrument of the Planck mission to test the impact that Galactic foreground emission and instrumental systematics have on the recent birefringence measurements obtained through this technique. Our results demonstrate the robustness of the methodology against the miscalibration of polarization angles and other systematic effects, like intensity-to-polarization leakage, beam leakage, or cross-polarization effects. However, our estimator is sensitive to the EB correlation of polarized foreground emission. Here we propose to correct the bias induced by dust EB by modeling the foreground signal with templates produced in Bayesian component-separation analyses that fit parametric models to CMB data. Acknowledging the limitations of currently available dust templates like that of the Commander sky model, high-precision CMB data and a characterization of dust beyond the modified blackbody paradigm are needed to obtain a definitive measurement of cosmic birefringence in the future.

79 ASTRONOMY AND ASTROPHYSICS↗

Reduced-Order Modeling of Hydrogen Releases from Vent Stacks and with Wind

This report describes computational modeling in the HyRAM+ software for study of hydrogen behavior in two common scenarios. First, models of unignited plumes exiting a vent stack were considered. It was shown that entrainment and vent backpressure were major factors in plume physics. Second, HyRAM+ was extended to include wind effects on plumes by modifying the plume momentum balance and entrainment modeling. Use of the model showed that plume shape and length changed with wind speed and direction; in all cases, wind causes a shortening of the plume along the streamline. While the no-wind case in HyRAM+ has been validated and the newly developed wind model was fitted to very limited experimental data, more controlled experimental configurations would help validate the models and ensure accurate simulation of hydrogen plume behavior for vent stack releases or in wind.

08 HYDROGEN↗

Backus-Gilbert inversion of travel time data

Application of the Backus-Gilbert theory for geophysical inverse problems to the seismic body wave travel-time problem is described. In particular, it is shown how to generate earth models that fit travel-time data to within one standard error and having generated such models how to describe their degree of uniqueness. An example is given to illustrate the process.

Johnson, L. E.↗