Search NASA⌕ Search

SEARCH · Search NASA

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

A Novel Multiphysics Multiscale Multiporosity Shale Gas Transport Model for Geomechanics/Flow Coupling in Steady and Transient States

Summary A novel multiphysics multiscale multiporosity shale gas transport (M3ST) model was developed to investigate shale gas transport in both transient and steady states. The microscale model component contains a kerogen domain and an inorganic matrix domain, and each domain has its own geomechanical and gas transport properties. Permeabilities of various shale cores were measured in the laboratory using a pulse decay permeameter (PDP) with different pore pressure and confining stress combinations. The PDP-measured apparent permeability as a function of pore pressure under two effective stresses was fitted using the microscale M3ST model component based on nonlinear least squares fitting (NLSF), and the fitted model parameters were able to provide accurate model predictions for another effective stress. The parameters and petrophysical properties determined in the steady state were then used in the transient-state, continuum-scale M3ST model component, which performed history matching of the evolutions of the upstream and downstream gas pressures. In addition, a double-exponential empirical model was developed as a powerful alternative to the M3ST model to fit laboratory-measured apparent permeability under various effective stresses and pore pressures. The developed M3ST model and the research findings in this study provided critical insights into the role of the multiphysics mechanisms, including geomechanics, fluid dynamics and transport, and the Klinkenberg effect on shale gas transport across different spatial scales in both steady and transient states.

Engineering↗

Detection of enhancement in number densities of background galaxies due to magnification by massive galaxy clusters

We present a detection of the enhancement in the number densities of background galaxies induced from lensing magnification and use it to test the Sunyaev-Zel'dovich effect (SZE-) inferred masses in a sample of 19 galaxy clusters with median redshift z similar or equal to 0.42 selected from the South Pole Telescope SPT-SZ survey. These clusters are observed by the Megacam on the Magellan Clay Telescope though gri filters. Two background galaxy populations are selected for this study through their photometric colours; they have median redshifts zmedian similar or equal to 0.9 (low-z background) and z(median) similar or equal to 1.8 (high-z background). Stacking these populations, we detect the magnification bias effect at 3.3 sigma and 1.3 sigma for the low-and high-z backgrounds, respectively. We fit Navarro, Frenk and White models simultaneously to all observed magnification bias profiles to estimate the multiplicative factor. that describes the ratio of the weak lensing mass to the mass inferred from the SZE observable-mass relation. We further quantify systematic uncertainties in. resulting from the photometric noise and bias, the cluster galaxy contamination and the estimations of the background properties. The resulting. for the combined background populations with 1 sigma uncertainties is 0.83 +/- 0.24(stat) +/- 0.074(sys), indicating good consistency between the lensing and the SZE-inferred masses. We use our best-fitting eta to predict the weak lensing shear profiles and compare these predictions with observations, showing agreement between the magnification and shear mass constraints. This work demonstrates the promise of using the magnification as a complementary method to estimate cluster masses in large surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling the Effects of Morphine-Altered Virus Specific Antibody Responses on HIV/SIV Dynamics

Abstract Drugs of abuse, such as opiates, have been widely associated with enhancing HIV replication, accelerating disease progression and diminishing host-immune responses, thereby making it harder to effectively manage HIV infection. It is thus important to study the effects of drugs of abuse on HIV-infection and immune responses. Here, we develop mathematical models that incorporate the effects of morphine-altered antibody responses on HIV/SIV dynamics. Based on fitting the model to experimental data from simian immunodeficiency virus (SIV) infections in control and morphine-addicted macaques, we found that two of the most significant effects of virus specific antibodies are neutralizing viral particles and enhancing viral clearance. Using our model, we quantified how morphine alters virus-specific antibody responses, and how this alteration affects the key components of virus dynamics such as infection rate, virus clearance, viral load, CD4 + T cell count, and CD4 + T cell loss in SIV-infected macaques under conditioning with morphine. We found that in a subpopulation of SIV-infected morphine addicted macaques, the presence of drugs of abuse may cause significantly diminished antibody responses, resulting in more severe infection with increased SIV infectivity, a decreased viral clearance rate, increased viral load, and higher CD4 + T cell loss.

59 BASIC BIOLOGICAL SCIENCES↗

Temperature- and Intensity-Dependent Photovoltaic Measurements to Identify Dominant Recombination Pathways

In novel photovoltaic absorbers, it is often difficult to assess the root causes of low open-circuit voltages, which may be due to bulk recombination or sub-optimal contacts. In the present work, we discuss the role of temperature- and illumination-dependent device electrical measurements in quantifying and distinguishing these performance losses - in particular, for determining bounds on interface recombination velocities, band alignment, and minority carrier lifetime. We assess the accuracy of this approach by direct comparison to photoelectron spectroscopy. Then, we demonstrate how more computationally intensive model parameter fitting approaches can draw more insights from this broad measurement space. We apply this measurement and modeling approach to high-performance III-V and thin-film chalcogenide devices.

14 SOLAR ENERGY↗

OSIF (Open Source Impedance Fitter)

OSIF is a tool that allows the user to easily fit electrochemical impedance spectra of proton exchange membrane fuel cells to a commonly accepted equivalent circuit model. The fitted parameters provide insight into the physical properties of the cell which allow the user to support experimental hypotheses.

Pfeilsticker, Jason↗

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↗

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↗

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↗

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↗

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↗

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↗

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.↗

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↗

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↗

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.↗

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↗

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↗

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↗