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 703 records · Page 39

10 Yr Spatial and Temporal Trends of PM2.5 Concentrations in the Southeastern US Estimated Using High-resolution Satellite Data

Long-term PM2.5 exposure has been reported to be associated with various adverse health outcomes. However, most ground monitors are located in urban areas, leading to a potentially biased representation of the true regional PM2.5 levels. To facilitate epidemiological studies, accurate estimates of spatiotemporally continuous distribution of PM2.5 concentrations are essential. Satellite-retrieved aerosol optical depth (AOD) has been widely used for PM2.5 concentration estimation due to its comprehensive spatial coverage. Nevertheless, an inherent disadvantage of current AOD products is their coarse spatial resolutions. For instance, the spatial resolutions of the Moderate Resolution Imaging Spectroradiometer (MODIS) and the Multiangle Imaging SpectroRadiometer (MISR) are 10 km and 17.6 km, respectively. In this paper, a new AOD product with 1 km spatial resolution retrieved by the multi-angle implementation of atmospheric correction (MAIAC) algorithm was used. A two-stage model was developed to account for both spatial and temporal variability in the PM2.5-AOD relationship by incorporating the MAIAC AOD, meteorological fields, and land use variables as predictors. Our study area is in the southeastern US, centered at the Atlanta Metro area, and data from 2001 to 2010 were collected from various sources. The model was fitted for each year individually, and we obtained model fitting R2 ranging from 0.71 to 0.85, MPE from 1.73 to 2.50 g m3, and RMSPE from 2.75 to 4.10 g m3. In addition, we found cross validation R2 ranging from 0.62 to 0.78, MPE from 2.00 to 3.01 g m3, and RMSPE from 3.12 to 5.00 g m3, indicating a good agreement between the estimated and observed values. Spatial trends show that high PM2.5 levels occurred in urban areas and along major highways, while low concentrations appeared in rural or mountainous areas. A time series analysis was conducted to examine temporal trends of PM2.5 concentrations in the study area from 2001 to 2010. The results showed that the PM2.5 levels in the study area followed a generally declining trend from 2001 to 2010 and decreased about 20 during the period. However, there was an exception of an increase in year 2005, which is attributed to elevated sulfate concentrations in the study area in warm months of 2005. An investigation of the impact of wild and prescribed fires on PM2.5 levels in 2007 suggests a positive relationship between them.

particulate matter (PM2.5)↗

Effects of Inventory Bias on Landslide Susceptibility Calculations

Many landslide inventories are known to be biased, especially inventories for large regions such as Oregon's SLIDO or NASA's Global Landslide Catalog. These biases must affect the results of empirically derived susceptibility models to some degree. We evaluated the strength of the susceptibility model distortion from postulated biases by truncating an unbiased inventory. We generated a synthetic inventory from an existing landslide susceptibility map of Oregon, then removed landslides from this inventory to simulate the effects of reporting biases likely to affect inventories in this region, namely population and infrastructure effects. Logistic regression models were fitted to the modified inventories. Then the process of biasing a susceptibility model was repeated with SLIDO data. We evaluated each susceptibility model with qualitative and quantitative methods. Results suggest that the effects of landslide inventory bias on empirical models should not be ignored, even if those models are, in some cases, useful. We suggest fitting models in well-documented areas and extrapolating across the study region as a possible approach to modeling landslide susceptibility with heavily biased inventories.

Oregon↗

Effects of Inventory Bias on Landslide Susceptibility Calculations

Many landslide inventories are known to be biased, especially inventories for large regions such as Oregons SLIDO or NASAs Global Landslide Catalog. These biases must affect the results of empirically derived susceptibility models to some degree. We evaluated the strength of the susceptibility model distortion from postulated biases by truncating an unbiased inventory. We generated a synthetic inventory from an existing landslide susceptibility map of Oregon, then removed landslides from this inventory to simulate the effects of reporting biases likely to affect inventories in this region, namely population and infrastructure effects. Logistic regression models were fitted to the modified inventories. Then the process of biasing a susceptibility model was repeated with SLIDO data. We evaluated each susceptibility model with qualitative and quantitative methods. Results suggest that the effects of landslide inventory bias on empirical models should not be ignored, even if those models are, in some cases, useful. We suggest fitting models in well-documented areas and extrapolating across the study region as a possible approach to modelling landslide susceptibility with heavily biased inventories.

landslide↗

The Data Synergy Effects of Time-Series Deep Learning Models in Hydrology

When fitting statistical models to variables in geoscientific disciplines such as hydrology, it is a customary practice to stratify a large domain into multiple regions (or regimes) and study each region separately. Traditional wisdom suggests that models built for each region separately will have higher performance because of homogeneity within each region. However, each stratified model has access to fewer and less diverse data points. Here, through two hydrologic examples (soil moisture and streamflow), we show that conventional wisdom may no longer hold in the era of big data and deep learning (DL). We systematically examined an effect we call data synergy, where the results of the DL models improved when data were pooled together from characteristically different regions. The performance of the DL models benefited from modest diversity in the training data compared to a homogeneous training set, even with similar data quantity. Moreover, allowing heterogeneous training data makes eligible much larger training datasets, which is an inherent advantage of DL. A large, diverse data set is advantageous in terms of representing extreme events and future scenarios, which has strong implications for climate change impact assessment. The results here suggest the research community should place greater emphasis on data sharing.

54 ENVIRONMENTAL SCIENCES↗

Fitting a 3-D analytic model of the coronal mass ejection to observations

The application of an analytic magnetohydrodynamic model is presented to observations of the time-dependent explusion of 3D coronal mass ejections (CMEs) out of the solar corona. This model relates the white-light appearance of the CME to its internal magnetic field, which takes the form of a closed bubble, filled with a partly anchored, twisted magnetic flux rope and embedded in an otherwise open background field. The density distribution frozen into the expanding CME expanding field is fully 3D, and can be integrated along the line of sight to reproduce observations of scattered white light. The model is able to reproduce the three conspicuous features often associated with CMEs as observed with white-light coronagraphs: a surrounding high-density region, an internal low-density cavity, and a high-density core. The model also describes the self-similar radial expansion of these structures. By varying the model parameters, the model can be fitted directly to observations of CMEs. It is shown how the model can quantitatively match the polarized brightness contrast of a dark cavity emerging through the lower corona as observed by the HAO Mauna Loa K-coronameter to within the noise level of the data.

Gibson, S. E.↗

Mathematical modeling of hepatitis C RNA replication, exosome secretion and virus release

Hepatitis C virus (HCV) causes acute hepatitis C and can lead to life-threatening complications if it becomes chronic. The HCV genome is a single plus strand of RNA. Its intracellular replication is a spatiotemporally coordinated process of RNA translation upon cell infection, RNA synthesis within a replication compartment, and virus particle production. While HCV is mainly transmitted via mature infectious virus particles, it has also been suggested that HCV-infected cells can secrete HCV RNA carrying exosomes that can infect cells in a receptor independent manner. In order to gain insight into these two routes of transmission, we developed a series of intracellular HCV replication models that include HCV RNA secretion and/or virus assembly and release. Fitting our models to in vitro data, in which cells were infected with HCV, suggests that initially most secreted HCV RNA derives from intracellular cytosolic plus-strand RNA, but subsequently secreted HCV RNA derives equally from the cytoplasm and the replication compartments. Furthermore, our model fits to the data suggest that the rate of virus assembly and release is limited by host cell resources. Including the effects of direct acting antivirals in our models, we found that in spite of decreasing intracellular HCV RNA and extracellular virus concentration, low level HCV RNA secretion may continue as long as intracellular RNA is available. This may possibly explain the presence of detectable levels of plasma HCV RNA at the end of treatment even in patients that ultimately attain a sustained virologic response.

59 BASIC BIOLOGICAL SCIENCES↗

Soft X-ray spectral fits of Geminga with model neutron star atmospheres

The spectrum of the soft X-ray pulsar Geminga consists of two components, a softer one which can be interpreted as thermal-like radiation from the surface of the neutron star, and a harder one interpreted as radiation from a polar cap heated by relativistic particles. We have fitted the soft spectrum using a detailed magnetized hydrogen atmosphere model. The fitting parameters are the hydrogen column density, the effective temperature T(sub eff), the gravitational redshift z, and the distance to radius ratio, for different values of the magnetic field B. The best fits for this model are obtained when B less than or approximately 1 x 10(exp 12) G and z lies on the upper boundary of the explored range (z = 0.45). The values of T(sub eff) approximately = (2-3) x 10(exp 5) K are a factor of 2-3 times lower than the value of T(sub eff) obtained for blackbody fits with the same z. The lower T(sub eff) increases the compatibility with some proposed schemes for fast neutrino cooling of neutron stars (NSs) by the direct Urca process or by exotic matter, but conventional cooling cannot be excluded. The hydrogen atmosphere fits also imply a smaller distance to Geminga than that inferred from a blackbody fit. An accurate evaluation of the distance would require a better knowledge of the ROSAT Position Sensitive Proportional Counter (PSPC) response to the low-energy region of the incident spectrum. Our modeling of the soft component with a cooler magnetized atmosphere also implies that the hard-component fit requires a characteristic temperature which is higher (by a factor of approximately 2-3) and a surface area which is smaller (by a factor of 10(exp 3), compared to previous blackbody fits.

Meyer, R. D.↗

Combining PCA and nonlinear fitting of peak models to re-evaluate C 1s XPS spectrum of cellulose

Cellulose is an example of a material that responds to XPS by the creation of new chemistry not present in the as-received sample. While improvements in instrumentation may be seen in general as beneficial to surface science, recent studies have shown that the consequences for some materials are detrimental. Here, in this work, these problems are illustrated through an analysis of cellulose spectra obtained during a degradation study. C 1s spectra are decomposed into two well-formed component curves that are open to chemical interpretation. In particular, a component-curve representative of pure cellulose is obtained as well as a second component curve that implies cellulose is degraded through the creation of carbon chemistry involving C—O, C$=$O and O—C$=$O. Since cellulose is a crystalline material, formed through the alignment of molecules under the influence of hydrogen bonds, the analysis and findings presented in this paper are relevant to any material analyzed by XPS whose properties are dependent on hydrogen bonds. The analysis techniques are based on an informed vectorial approach, which extracts directly from data spectral shapes that are used to monitor sample degradation via linear least squares optimization. Related mathematics of Principal Component Analysis and linear analysis are presented.

36 MATERIALS SCIENCE↗

Strongly temperature dependent ferroelectric switching in AlN, Al 1-x Sc x N, and Al 1-x B x N thin films

This manuscript reports the temperature dependence of ferroelectric switching in Al 0.84 Sc 0.16 N, Al 0.93 B 0.07 N, and AlN thin films. Polarization reversal is demonstrated in all compositions and is strongly temperature dependent. Between room temperature and 300 °C, the coercive field drops by almost 50% in all samples, while there was very small temperature dependence of the remanent polarization value. Furthermore, over this same temperature range, the relative permittivity increased between 5% and 10%. Polarization reversal was confirmed by piezoelectric coefficient analysis and chemical etching. Applying intrinsic/homogeneous switching models produces nonphysical fits, while models based on thermal activation suggest that switching is regulated by a distribution of pinning sites or nucleation barriers with an average activation energy near 28 meV.

36 MATERIALS SCIENCE↗

Case study: Accounting for response measurement error in fitting a regression model

This article presents a case study motivated by a plot of data that suggested an emerging trend that the authors were faced with explaining. When the measurement error of the data is accounted for, it turns out there was no real trend. Furthermore, this article shows how to use a Bayesian modeling approach to account for the measurement error.

42 ENGINEERING↗

Characterizing the gene–environment interaction underlying natural morphological variation in Neurospora crassa conidiophores using high-throughput phenomics and transcriptomics

Abstract Neurospora crassa propagates through dissemination of conidia, which develop through specialized structures called conidiophores. Recent work has identified striking variation in conidiophore morphology, using a wild population collection from Louisiana, United States of America to classify 3 distinct phenotypes: Wild-Type, Wrap, and Bulky. Little is known about the impact of these phenotypes on sporulation or germination later in the N. crassa life cycle, or about the genetic variation that underlies them. In this study, we show that conidiophore morphology likely affects colonization capacity of wild N. crassa isolates through both sporulation distance and germination on different carbon sources. We generated and crossed homokaryotic strains belonging to each phenotypic group to more robustly fit a model for and estimate heritability of the complex trait, conidiophore architecture. Our fitted model suggests at least 3 genes and 2 epistatic interactions contribute to conidiophore phenotype, which has an estimated heritability of 0.47. To uncover genes contributing to these phenotypes, we performed RNA-sequencing on mycelia and conidiophores of strains representing each of the 3 phenotypes. Our results show that the Bulky strain had a distinct transcriptional profile from that of Wild-Type and Wrap, exhibiting differential expression patterns in clock-controlled genes (ccgs), the conidiation-specific gene con-6, and genes implicated in metabolism and communication. Combined, these results present novel ecological impacts of and differential gene expression underlying natural conidiophore morphological variation, a complex trait that has not yet been thoroughly explored.

59 BASIC BIOLOGICAL SCIENCES↗

Model Exploration of an Information-Based Healthcare Intervention Using Parallelization and Active Learning

This paper describes the application of a large-scale active learning method to characterize the parameter space of a computational agent-based model developed to investigate the impact of CommunityRx, a clinical information-based health intervention that provides patients with personalized information about local community resources to meet basic and self-care needs. Additionally, the diffusion of information about community resources and their use is modeled via networked interactions and their subsequent effect on agents' use of community resources across an urban population. A random forest model is iteratively fitted to model evaluations to characterize the model parameter space with respect to observed empirical data. We demonstrate the feasibility of using high-performance computing and active learning model exploration techniques to characterize large parameter spaces; by partitioning the parameter space into potentially viable and non-viable regions, we rule out regions of space where simulation output is implausible to observed empirical data. We argue that such methods are necessary to enable model exploration in complex computational models that incorporate increasingly available micro-level behavior data. We provide public access to the model and high-performance computing experimentation code.

97 MATHEMATICS AND COMPUTING↗

Stochastic Models, Indices & Optimization Algorithms for Pricing & Hedging Reliability Risks in Modern Power Grids: Data Plan - Princeton

We collected and cleaned the synthetic grid data produced by NREL for the Texas and New York synthetic grids. We developed a high dimensional joint stochastic model for load at the zone level, and solar and wind power productions at the asset level, capturing the spatial and temporal dependencies between all the variables, and demonstrated how such a model could be fitted to historical data. We designed and implemented a simulation engine which can produce Monte Carlo scenarios for the hourly day-ahead values of load, and solar and wind power productions at the spatial and temporal resolutions of the historical data used to fit the model. Finally we developed an open-source Python package which can, from an input grid model, efficiently use forecasts and large numbers of Monte Carlo scenarios to provide unit commitment and economic dispatch for each of these scenarios. The high dimensional stochastic model and the subsequent Monte Carlo simulation engine were implemented in the package PGscen and the corresponding UC and ED optimization programs in the package Vatic.

14 SOLAR ENERGY↗

Spectral fitting, shock layer modeling, and production of nitrogen oxides and excited nitrogen

An analysis was made of N2 emission from 8.72 MJ/kg shock layer at 2.54, 1.91, and 1.27 cm positions and vibrational state distributions, temperatures, and relative electronic state populations was obtained from data sets. Other recorded arc jet N2 and air spectral data were reviewed and NO emission characteristics were studied. A review of operational procedures of the DSMC code was made. Information on other appropriate codes and modifications, including ionization, were made as well as a determination of the applicability of codes reviewed to task requirement. A review was also made of computational procedures used in CFD codes of Li and other codes on JSC computers. An analysis was made of problems associated with integration of specific chemical kinetics applicable to task into CFD codes.

Blackwell, H. E.↗