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,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Quantifying the Turnover of Transcriptional Subclasses of HIV-1-Infected Cells

HIV-1-infected cells in peripheral blood can be grouped into different transcriptional subclasses. Quantifying the turnover of these cellular subclasses can provide important insights into the viral life cycle and the generation and maintenance of latently infected cells. We used previously published data from five patients chronically infected with HIV-1 that initiated combination antiretroviral therapy (cART). Patient-matched PCR for unspliced and multiply spliced viral RNAs combined with limiting dilution analysis provided measurements of transcriptional profiles at the single cell level. Furthermore, measurement of intracellular transcripts and extracellular virion-enclosed HIV-1 RNA allowed us to distinguish productive from non-productive cells. We developed a mathematical model describing the dynamics of plasma virus and the transcriptional subclasses of HIV-1-infected cells. Fitting the model to the data allowed us to better understand the phenotype of different transcriptional subclasses and their contribution to the overall turnover of HIV-1 before and during cART. The average number of virus-producing cells in peripheral blood is small during chronic infection. We find that a substantial fraction of cells can become defectively infected. Assuming that the infection is homogenous throughout the body, we estimate an average in vivo viral burst size on the order of 104 virions per cell. Our study provides novel quantitative insights into the turnover and development of different subclasses of HIV-1-infected cells, and indicates that cells containing solely unspliced viral RNA are a good marker for viral latency. The model illustrates how the pool of latently infected cells becomes rapidly established during the first months of acute infection and continues to increase slowly during the first years of chronic infection. Having a detailed understanding of this process will be useful for the evaluation of viral eradication strategies that aim to deplete the latent reservoir of HIV-1.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Reconstructing the Universe: Testing the Mutual Consistency of the Pantheon and SDSS/eBOSS BAO Data Sets with Gaussian Processes

We test the mutual consistency between the baryon acoustic oscillation measurements from the eBOSS SDSS final release and the Pantheon supernova compilation in a model-independent fashion using Gaussian process regression. We also test their joint consistency with the ΛCDM model in a model-independent fashion. We also use Gaussian process regression to reconstruct the expansion history that is preferred by these two data sets. While this methodology finds no significant preference for model flexibility beyond ΛCDM, we are able to generate a number of reconstructed expansion histories that fit the data better than the best-fit ΛCDM model. These example expansion histories may point the way toward modifications to ΛCDM. We also constrain the parameters Ω{sub k} and H {sub 0} r {sub d} both with ΛCDM and with Gaussian process regression. We find that H {sub 0} r {sub d} = 10,030 ± 130 km s{sup −1} and Ω{sub k} = 0.05 ± 0.10 for ΛCDM and that H {sub 0} r {sub d} = 10,040 ± 140 km s{sup −1} and Ω{sub k} = 0.02 ± 0.20 for the Gaussian process case.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1999

The Impact of SN1987A with its Circumstellar Ring

Following is final report on the study "The Impact of SN1987A with its Circumstellar Ring", which is now complete. In 1994, it was predicted that the blast wave from SN1987A would strike the circumstellar ring in AD1999+/-3, and that the ring would brighten by several hundreds optically as the shock entered the ring. It was also predicted that the emission lines from the shocked ring would have linewidths of a few hundred km/s. Today, we see a "hot spot" on the ring that first appeared in Hubble Space Telescope (HST) images in 1995 and has doubled in brightness between August 1997 and February 1998. Moreover, spectra from STIS show that the emission lines from the hot spot have widths and blue shifts of order 200 km/s, just as we predicted in 1994. Our guess of 1999+/-3 for the impact time was lucky, because we assumed that the hot gas in the bubble between the supernova and the ring had a low density (approximately 10 cm(exp -3)). But it was pointed out that the ROSAT observations of soft X-rays from SN 1987A implied that the intervening gas had a higher density, (approximately 100 cm(exp -3)), which would delay the impact until ca. 2007. Others developed hydrodynamic models to fit the ROSAT X-ray emission spectrum and came to the same conclusion. But these models were oversimplified in that they assumed that the ring was perfectly round. Now we see clearly that the hot spot is a peninsula that protrudes inward from the ring -- the first spot on the ring to be struck by the blast wave. No doubt there are other protrusions on the ring, which we may expect to light up in the next few years until they finally merge to set the entire ring ablaze, probably within 5 - 10 years. It was also predicted that the X-ray emission must be accompanied by optical and ultraviolet emission from atoms in the supernova debris and in the circumstellar gas that cross the reverse shock and the blast wave, respectively, and that the Lyman-alpha and NV emission lines should be bright enough to see with the STIS.

McCray, Richard↗
NASA NTRS2019

Modelling above-ground biomass stock over Norway using national forest inventory data with ArcticDEM and Sentinel-2 data

Boreal forests constitute a large portion of the global forest area, yet they are undersampled through field surveys, and only a few remotely sensed data sources provide structural information wall-to-wall throughout the boreal domain. ArcticDEM is a collection of high-resolution (2 m) space-borne stereogrammetric digital surface models (DSM) covering the entire land area north of 60° of latitude. The free-availability of ArcticDEM data offers new possibilities for aboveground biomass mapping (AGB) across boreal forests, and thus it is necessary to evaluate the potential for these data to map AGB over alternative open-data sources (i.e., Sentinel-2). This study was performed over the entire land area of Norway north of 60° of latitude, and the Norwegian national forest inventory (NFI) was used as a source of field data composed of accurately geolocated field plots (n=7710) systematically distributed across the study area. Separate random forest models were fitted using NFI data, and corresponding remotely sensed data consisting of either: i) a canopy height model (ArcticCHM) obtained by subtracting a high-quality digital terrain model (DTM) from the ArcticDEM DSM height values, ii) Sentinel-2 (S2), or iii) a combination of the two (ArcticCHM+S2). Furthermore, we assessed the effect of the forest- and terrain-specific factors on the models’ predictive accuracy. The best model (,i.e., ArcticCHM+S2) explained nearly 60% of the variance of the training set, which translated in the largest accuracy in terms of root mean square error (RMSE=41.4 t/ha). This result highlights the synergy between 3D and multispectral data in AGB modelling. Furthermore, this study showed that despite the importance of ArcticCHM variables, the S2 model performed slightly better than ArcticCHM model. This finding highlights some of the limitations of ArcticDEM, which, despite the unprecedented spatial resolution, is highly heterogeneous due to the blending of multiple acquisitions across different years and seasons. We found that both forest- and terrain-specific characteristics affected the uncertainty of the ArcticCHM+S2 model and concluded that the combined use of ArcticCHM and Sentinel-2 represents a viable solution for AGB mapping across boreal forests. The synergy between the two data sources allowed for a reduction of the saturation effects typical of multispectral data while ensuring the spatial consistency in the output predictions due to the removal of artifacts and data voids present in ArcticCHM data. While the main contribution of this study is to provide the first evidence of the best-case-scenario (i.e., availability of accurate terrain models) that ArcticDEM data can provide for large-scale AGB modelling, it remains critically important for other studies to investigate how ArcticDEM may be used in areas where no DTMs are available as is the case for large portions of the boreal zone.

space-borne imagery↗
DOE OSTI2017

Optical and ultraviolet spectroscopic analysis of SN 2011fe at late times

This paper presents optical spectra of the nearby Type Ia supernova SN 2011fe at 100, 205, 311, 349 and 578 d post-maximum light, as well as an ultraviolet (UV) spectrum obtained with the Hubble Space Telescope at 360 d post-maximum light. We compare these observations with synthetic spectra produced with the radiative transfer code PHOENIX. The day +100 spectrum can be well fitted with models that neglect collisional and radiative data for forbidden lines. Curiously, including these data and recomputing the fit yields a quite similar spectrum, but with different combinations of lines forming some of the stronger features. At day +205 and later epochs, forbidden lines dominate much of the optical spectrum formation; however, our results indicate that recombination, not collisional excitation, is the most influential physical process driving spectrum formation at these late times. Consequently, our synthetic optical and UV spectra at all epochs presented here are formed almost exclusively through recombinationdriven fluorescence. Furthermore, our models suggest that the UV spectrum even as late as day +360 is optically thick and consists of permitted lines from several iron-peak species. These results indicate that the transition to the 'nebular' phase in Type Ia supernovae is complex and highly wavelength dependent.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2012

Using Random Forest Models to Predict Organizational Violence

We present a methodology to access the proclivity of an organization to commit violence against nongovernment personnel. We fitted a Random Forest model using the Minority at Risk Organizational Behavior (MAROS) dataset. The MAROS data is longitudinal; so, individual observations are not independent. We propose a modification to the standard Random Forest methodology to account for the violation of the independence assumption. We present the results of the model fit, an example of predicting violence for an organization; and finally, we present a summary of the forest in a "meta-tree,"

Levine, Burton↗
NASA NTRS1991

On the reconciliation of simultaneous microwave imaging and hard X-ray observations of a solar flare

Microwave imaging data for a small flare with simultaneous hard X-ray spectral observations are compared. The X-ray data suggest that the power-law index delta of the energy distribution of the radiating electrons is 5.3 (thick-target) which differs significantly from the estimate (delta = 1.4) from a homogeneous optically-thin gyrosynchrotron model which fits the radio observations well. In order to reconcile these results, a double power-law energy spectrum is investigated for the energetic electrons in the flare, as assumed by other authors: the power law is steep at low energies and much flatter at the higher energies which produce the bulk of the microwaves. The emission of soft photons by the flat tail strongly contributes to the observed hard X-ray range and would flatten the spectrum there. A thin-target model for the X-ray emission is also inconsistent with radio data. An inhomogenous gyrosychrotron model with a number of free parameters and containing an electron distribution given by the thick-target X-ray model could be made to fit the radio data.

Nitta, N.↗
DOE OSTI2021

Gemini Planet Imager Spectroscopy of the Dusty Substellar Companion HD 206893 B

We present new near-infrared Gemini Planet Imager (GPI) spectroscopy of HD 206893 B, a substellar companion orbiting within the debris disk of its F5V star. The J, H, K1, and K2 spectra from GPI demonstrate the extraordinarily red colors of the object, confirming it as the reddest substellar object observed to date. The significant flux increase throughout the infrared presents a challenging atmosphere to model with existing grids. Best-fit values vary from 1200 to 1800 K for effective temperature and from 3.0 to 5.0 for log(g), depending on which individual wavelength band is fit and which model suite is applied. The extreme redness of the companion can be partially reconciled by invoking a high-altitude layer of submicron dust particles, similar to dereddening approaches applied to the peculiar red field L dwarf population. However, reconciling the HD 206893 B spectra with even those of the reddest low-gravity L dwarf spectra still requires the contribution of additional atmospheric dust, potentially due to the debris disk environment in which the companion resides. Orbit fitting from 4 yr of astrometric monitoring is consistent with a ∼30 yr period, an orbital inclination of 147°, and a semimajor axis of 10 au, well within the estimated disk inner radius of ∼50 au. As one of a very few substellar companions imaged interior to a circumstellar disk, the properties of this system offer important dynamical constraints on companion–disk interaction and provide a benchmark for substellar and planetary atmospheric study.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Kepler and TESS Observations of PG 1159-035

Abstract PG 1159-035 is the prototype of the PG 1159 hot (pre-)white dwarf pulsators. This important object was observed during the Kepler satellite K2 mission for 69 days in 59 s cadence mode and by the TESS satellite for 25 days in 20 s cadence mode. We present a detailed asteroseismic analysis of those data. We identify a total of 107 frequencies representing 32 ℓ = 1 modes, 27 frequencies representing 12 ℓ = 2 modes, and eight combination frequencies. The combination frequencies and the modes with very high k values represent new detections. The multiplet structure reveals an average splitting of 4.0 ± 0.4 μ Hz for ℓ = 1 and 6.8 ± 0.2 μ Hz for ℓ = 2, indicating a rotation period of 1.4 ± 0.1 days in the region of period formation. In the Fourier transform of the light curve, we find a significant peak at 8.904 ± 0.003 μ Hz suggesting a surface rotation period of 1.299 ± 0.002 days. We also present evidence that the observed periods change on timescales shorter than those predicted by current evolutionary models. Our asteroseismic analysis finds an average period spacing for ℓ = 1 of 21.28 ± 0.02 s. The ℓ = 2 modes have a mean spacing of 12.97 ± 0.4 s. We performed a detailed asteroseismic fit by comparing the observed periods with those of evolutionary models. The best-fit model has T eff = 129, 600 ± 11 100 K, M * = 0.565 ± 0.024 M ⊙ , and log g = 7.41 − 0.54 + 0.38 , within the uncertainties of the spectroscopic determinations. We argue for future improvements in the current models, e.g., on the overshooting in the He-burning stage, as the best-fit model does not predict excitation for all of the pulsations detected in PG 1159-035.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1991

Dust scattering in the high-state eclipse of Cen X-3

Exosat data from two high-state eclipses of Cen X-3 are analyzed for the effects of dust scattering. The light curves show bright (2 per cent of the out-of-eclipse flux), smoothly varying emission which reaches its minimum after mideclipse. Of this emission, the soft flux falls monotonically with orbital phase. This soft component with dust scattering is identified and spectral models are fitted to the asymmetric emission, that is, to spectra obtained by subtracting from a spectrum at a given orbital phase before phase zero the corresponding spectrum from the equal and opposite phase after phase zero. It is found that the model of Molnar and Mauche, used to fit dust scattering in Cyg X-3, can explain the observed spectra and yields values for the distance to Cen X-3 and for the size of the dust grains of 5.4 kpc and 0.33 micron, consistent with previous observations in other wavebands.

Day, C. S. R.↗
DOE OSTI2016

A new correlation with lower kilohertz quasi-periodic oscillation frequency in the ensemble of low-mass X-ray binaries

We study the dependence of kilohertz quasi-periodic oscillation (kHz QPO) frequency on accretion-related parameters in the ensemble of neutron-star low-mass X-ray binaries. Based on the mass accretion rate, M-dot , and the magnetic field strength, B, on the surface of the neutron star, we find a correlation between the lower kHz QPO frequency and M-dot /B{sup 2}. The correlation holds in the current ensemble of Z and atoll sources and therefore can explain the lack of correlation between the kHz QPO frequency and X-ray luminosity in the same ensemble. The average run of lower kHz QPO frequencies throughout the correlation can be described by a power-law fit to source data. The simple power law, however, cannot describe the frequency distribution in an individual source. The model function fit to frequency data, on the other hand, can account for the observed distribution of lower kHz QPO frequencies in the case of individual sources as well as the ensemble of sources. The model function depends on the basic length scales, such as the magnetospheric radius and the radial width of the boundary region, both of which are expected to vary with M-dot to determine the QPO frequencies. In addition to modifying the length scales, and hence the QPO frequencies, the variation in M-dot , being sufficiently large, may also lead to distinct accretion regimes, which would be characterized by Z and atoll phases.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Evaluation of neutron light output response functions in EJ-309 organic scintillators

An accurate model of the nonlinear detector response of organic scintillators to neutrons is required to correctly simulate fast neutron detection, as well as interpret measured pulse height data. Several empirical and semi-empirical models are available to fit measured scintillator light output data. In this work, EJ-309 light output data from neutrons depositing 1.15–5.15 MeV on hydrogen were analyzed using empirical models as well as semi-empirical models based on the work of Birks and Voltz. Although all tested models fit the experimental light output data well in the measured range, the models were observed to diverge in low-energy extrapolation. The models were then tested by comparing a measurement and MCNPX-PoliMi simulation of an EJ-309 detector response to fast neutrons from a 252 Cf spontaneous fission source. The agreement between the measured and simulated pulse height distributions varied significantly depending on the light output model used. The best agreement between simulated and measured neutron pulse height distributions was achieved by using the Birks model. The bin-by-bin agreement was better than 5% over the range 0.08–2.18 MeVee, and better than 10% from 2.18 to 3.13 MeVee. Furthermore, the integral count rate over the range 0.08–3.14 MeVee differed by less than 1% in absolute units.

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

CF + excitation in the interstellar medium

The detection of CF + in interstellar clouds potentially allows astronomers to infer the elemental fluorine abundance and the ionization fraction in ultraviolet-illuminated molecular gas. Because local thermodynamic equilibrium (LTE) conditions are hardly fulfilled in the interstellar medium (ISM), the accurate determination of the CF + abundance requires one to model its non-LTE excitation via both radiative and collisional processes. Here, we report quantum calculations of rate coefficients for the rotational excitation of CF + in collisions with para- and ortho-H 2 (for temperatures up to 150 K). As an application, we present non-LTE excitation models that reveal population inversion in physical conditions typical of ISM photodissociation regions (PDRs). We successfully applied these models to fit the CF + emission lines previously observed toward the Orion Bar and Horsehead PDRs. The radiative transfer models achieved with these new rate coefficients allow the use of CF + as a powerful probe to study molecular clouds exposed to strong stellar radiation fields.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

The kinetics of SARS-CoV-2 infection based on a human challenge study

Studying the early events that occur after viral infection in humans is difficult unless one intentionally infects volunteers in a human challenge study. Here, we use data about severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in such a study in combination with mathematical modeling to gain insights into the relationship between the amount of virus in the upper respiratory tract and the immune response it generates. We propose a set of dynamic models of increasing complexity to dissect the roles of target cell limitation, innate immunity, and adaptive immunity in determining the observed viral kinetics. We introduce an approach for modeling the effect of humoral immunity that describes a decline in infectious virus after immune activation. We fit our models to viral load and infectious titer data from all the untreated infected participants in the study simultaneously. We found that a power-law with a power h < 1 describes the relationship between infectious virus and viral load. Viral replication at the early stage of infection is rapid, with a doubling time of ~2 h for viral RNA and ~3 h for infectious virus. We estimate that adaptive immunity is initiated ~7 to 10 d postinfection and appears to contribute to a multiphasic viral decline experienced by some participants; the viral rebound experienced by other participants is consistent with a decline in the interferon response. Altogether, we quantified the kinetics of SARS-CoV-2 infection, shedding light on the early dynamics of the virus and the potential role of innate and adaptive immunity in promoting viral decline during infection.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2014

Modeling active region transient brightenings observed with X-ray telescope as multi-stranded loops

Strong evidence exists that coronal loops as observed in extreme ultraviolet and soft X-rays may not be monolithic isotropic structures, but can often be more accurately modeled as bundles of independent strands. Modeling the observed active region transient brightenings (ARTBs) within this framework allows for the exploration of the energetic ramifications and characteristics of these stratified structures. Here we present a simple method of detecting and modeling ARTBs observed with the Hinode X-Ray Telescope (XRT) as groups of zero-dimensional strands, which allows us to probe parameter space to better understand the spatial and temporal dependence of strand heating in impulsively heated loops. This partially automated method can be used to analyze a large number of observations to gain a statistical insight into the parameters of coronal structures, including the number of heating events required in a given model to fit the observations. In this article, we present the methodology and demonstrate its use in detecting and modeling ARTBs in a sample data set from Hinode/XRT. These initial results show that, in general, multiple heating events are necessary to reproduce observed ARTBs, but the spatial dependence of these heating events cannot yet be established.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Fitting FFT-derived spectra: Theory, tool, and application to solar radio spike decomposition

Spectra derived from fast Fourier transform (FFT) analysis of time-domain data intrinsically contain statistical fluctuations whose distribution depends on the number of accumulated spectra contributing to a measurement. The tail of this distribution, which is essential for separating the true signal from the statistical fluctuations, deviates noticeably from the normal distribution for a finite number of accumulations. In this paper, we develop a theory to properly account for the statistical fluctuations when fitting a model to a given accumulated spectrum. The method is implemented in software for the purpose of automatically fitting a large body of such FFT-derived spectra. We apply this tool to analyze a portion of a dense cluster of spikes recorded by our FASR Subsystem Testbed instrument during a record-breaking event that occurred on 2006 December 6. The outcome of this analysis is briefly discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Battery Performance, Thermal, and Life Modeling for Southern California Edison (Cooperative Research and Development Final Report)

The objective of this project is to develop battery performance models, thermal models and life models for SCE to use in simulation and analysis of grid services involving energy storage. NREL will fit the models to three separate chemistries using test data provided by SCE. The dataset for each chemistry will include electrical, thermal and aging response to different temperatures and cycling conditions, measured under a variety of cell-level and module-level experiments described in SCE test protocol documents. Once tuned to a test dataset, the battery system model software will provide predictions of battery energy and power loss, and thus cycle and calendar life, for any energy storage grid service power profile of interest to SCE. The automated battery life modeling and simulation tool can also be applied to any other grid battery systems as long as necessary input data is available to accelerate battery model development and battery lifetime analysis.

25 ENERGY STORAGE↗
DOE OSTI2019

Using the Tilted flat-ΛCDM and the Untilted Non-flat ΛCDM Inflation Models to Measure Cosmological Parameters from a Compilation of Observational Data

We use the physically consistent tilted spatially flat and untilted non-flat ΛCDM inflation models to constrain cosmological parameter values with the Planck 2015 cosmic microwave background (CMB) anisotropy data and recent SNe Ia measurements, baryonic acoustic oscillations (BAO) data, growth rate observations, and Hubble parameter measurements. The most dramatic consequence of including the four non-CMB data sets is the significant strengthening of the evidence for non-flatness in the non-flat ΛCDM model, from 1.8σ for the CMB data alone to 5.1σ for the full data combination. The BAO data is the most powerful of the non-CMB data sets in more tightly constraining model-parameter values and in favoring a spatially closed universe in which spatial curvature contributes about a percent to the current cosmological energy budget. The untilted non-flat ΛCDM model better fits the large-angle CMB temperature anisotropy angular spectrum and is more consistent with the Dark Energy Survey constraints on the current value of the rms amplitude of mass fluctuations (σ 8 ) as a function of the current value of the nonrelativistic matter-density parameter (Ω m ) but does not provide as good a fit to the smaller-angle CMB temperature anisotropy data, as does the tilted flat-ΛCDM model. Some measured cosmological parameter values differ significantly between the two models, including the reionization optical depth and the baryonic matter density parameter, both of whose 2σ ranges (in the two models) are disjointed or almost so.

79 ASTRONOMY AND ASTROPHYSICS↗