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,099 records · Page 61

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↗

X-ray perspective of the twisted magnetospheres of magnetars

Anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs) are recognized as the most promising magnetar candidates, as indicated by their energetic bursts and rapid spin-downs. It is expected that the strong magnetic field leaves distinctive imprints on the emergent radiation both by affecting the radiative processes in atmospheres of magnetars and by scattering in the upper magnetospheres. We construct a self-consistent physical model that incorporates emission from the magnetar surface and its reprocessing in the three-dimensional (3D) twisted magnetosphere using a Monte Carlo technique. The synthetic spectra are characterized by four parameters: surface temperature kT, surface magnetic field strength B, magnetospheric twist angle Δϕ, and the normalized electron velocity β. We also create a tabular model (STEMS3D) and apply it to a large sample of XMM-Newton spectra of magnetars. The model successfully fits nearly all spectra, and the obtained magnetic field for 7 out of the 11 sources are consistent with the values inferred from the spin-down rates. We conclude that the continuum-fitting by our model is a robust method to measure the magnetic field strength and magnetospheric configuration of AXPs and SGRs. Investigating the multiple observations of variable sources, we also study the mechanism of their spectral evolution. Our results suggest that the magnetospheres in these sources are highly twisted (Δϕ>1), and the behavior of magnetospheric twisting and untwisting is revealed in the 2002 outburst of 1E 2259+586.

79 ASTRONOMY AND ASTROPHYSICS↗

Wiggly whipped inflation

Motivated by BICEP2 results on the CMB polarization B-mode which imply primordial gravitational waves are produced when the Universe has the expansion rate of about H ≈ 10{sup 14} GeV, and by deviations from a smooth power-law behavior for multipoles ℓ <50 in the CMB temperature anisotropy power spectrum found in the WMAP and Planck experiments, we have expanded our class of large field inflationary models that fit both the BICEP2 and Planck CMB observations consistently. These best-fitted large field models are found to have a transition from a faster roll to the slow roll V(φ)=m{sup 2} φ{sup 2}/2 inflation at a field value around 14.6 M{sub Pl} and thus a potential energy of V(φ) ∼ (10{sup 16} GeV){sup 4}. In general this transition with sharp features in the inflaton potential produces not only suppression of scalars relative to tensor modes at small k but also introduces wiggles in the primordial perturbation spectrum. These wiggles are shown to be useful to explain some localized features in the CMB angular power spectrum and can also have other observational consequences. Thus, primordial GW can be used now to make a tomography of inflation determining its fine structure. The resulting Wiggly Whipped Inflation scenario is described in details and the anticipated perturbation power spectra, CMB power spectra, non-Gaussianity and other observational consequences are calculated and compared to existing and forthcoming observations.

79 ASTRONOMY AND ASTROPHYSICS↗

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↗

SU-E-T-751: Three-Component Kinetic Model of Tumor Growth and Radiation Response for Stereotactic Radiosurgery

Purpose: To develop and study a kinetic model of tumor growth and its response to stereotactic radiosurgery (SRS) by assuming that the cells in irradiated tumor volume were made of three types. Methods: A set of ordinary differential equations (ODEs) were derived for three types of cells and a tumor growth rate. It is assumed that the cells were composed of actively proliferating cells, lethally damaged-dividing cells, and non-dividing cells. We modeled the tumor volume growth with a time-dependent growth rate to simulate the saturation of growth. After SRS, the proliferating cells were permanently damaged and converted to the lethally damaged cells. The amount of damaged cells were estimated by the LQ-model. The damaged cells gradually stopped dividing/proliferating and died with a constant rate. The dead cells were cleared from their original location with a constant rate. The total tumor volume was the sum of the three components. The ODEs were numerically solved with appropriate initial conditions for a given dosage. The proposed model was used to model an animal experiment, for which the temporal change of a rhabdomyosarcoma tumor volume grown in a rat was measured with time resolution sufficient to test the model. Results: To fit the model to the experimental data, the following characteristics were needed with the model parameters. The α-value in the LQ-model was smaller than the commonly used value; furthermore, it decreased with increasing dose. At the same time, the tumor growth rate after SRS had to increase. Conclusions: The new 3-component model of tumor could simulate the experimental data very well. The current study suggested that the radiation sensitivity and the growth rate of the proliferating tumor cells may change after irradiation and it depended on the dosage used for SRS. These preliminary observations must be confirmed by future animal experiments.

60 APPLIED LIFE 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↗

Improved evidence-based genome-scale metabolic models for maize leaf, embryo, and endosperm

There is a growing demand for genome-scale metabolic reconstructions for plants, fueled by the need to understand the metabolic basis of crop yield and by progress in genome and transcriptome sequencing. Methods are also required to enable the interpretation of plant transcriptome data to study how cellular metabolic activity varies under different growth conditions or even within different organs, tissues, and developmental stages. Such methods depend extensively on the accuracy with which genes have been mapped to the biochemical reactions in the plant metabolic pathways. Errors in these mappings lead to metabolic reconstructions with an inflated number of reactions and possible generation of unreliable metabolic phenotype predictions. Here we introduce a new evidence-based genome-scale metabolic reconstruction of maize, with significant improvements in the quality of the gene-reaction associations included within our model. We also present a new approach for applying our model to predict active metabolic genes based on transcriptome data. This method includes a minimal set of reactions associated with low expression genes to enable activity of a maximum number of reactions associated with high expression genes. We apply this method to construct an organ-specific model for the maize leaf, and tissue specific models for maize embryo and endosperm cells. We validate our models using fluxomics data for the endosperm and embryo, demonstrating an improved capacity of our models to fit the available fluxomics data. All models are publicly available via the DOE Systems Biology Knowledgebase and PlantSEED, and our new method is generally applicable for analysis transcript profiles from any plant, paving the way for further in silico studies with a wide variety of plant genomes.

flux balance analysis↗

Improved Dynamical Constraints on the Mass of the Central Black Hole in NGC 404

We explore the nucleus of the nearby 10{sup 9} M{sub ⊙} early-type galaxy, NGC 404, using Hubble Space Telescope (HST)/STIS spectroscopy and WFC3 imaging. We first present evidence for nuclear variability in UV, optical, and infrared filters over a time period of 15 years. This variability adds to the already substantial evidence for an accreting black hole at the center of NGC 404. We then redetermine the dynamical black hole mass in NGC 404 including modeling of the nuclear stellar populations. We combine HST/STIS spectroscopy with WFC3 images to create a local color–M/L relation derived from stellar population modeling of the STIS data. We then use this to create a mass model for the nuclear region. We use Jeans modeling to fit this mass model to adaptive optics stellar kinematic observations from Gemini/NIFS. From our stellar dynamical modeling, we find a 3σ upper limit on the black hole mass of 1.5×10{sup 5} M{sub ⊙}. Given the accretion evidence for a black hole, this upper limit makes NGC 404 the lowest mass central black hole with dynamical mass constraints. We find that the kinematics of H{sub 2} emission line gas show evidence for non-gravitational motions preventing the use of gas dynamical modeling to constrain the black hole mass. Our stellar population modeling also reveals that the central, counter-rotating region of the nuclear cluster is dominated by ∼1 Gyr old populations.

79 ASTRONOMY AND ASTROPHYSICS↗

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

Eruption of a Multiple-Turn Helical Magnetic Flux Tube in a Large Flare: Evidence for External and Internal Reconnection that Fits the Breakout Model of Solar Magnetic Eruptions

We present observations and an interpretation of a unique multiple-turn spiral flux tube eruption from AR10030 on 2002 July 15. The TRACE CIV observations clearly show a flux tube that is helical and that is erupting from within a sheared magnetic field. These observations are interpreted in the context of the breakout model for magnetic field explosions. The initiation of the helix eruption starts 25 seconds after the peak of the flare s strongest impulsive spike of microwave gryosynchrotron radiation early in the flare s explosive phase, implying that the sheared core field is not the site of the initial reconnection. Within the quadrupolar configuration of the active region, the external and internal reconnection sites are identified in each of two consecutive eruptive flares that produce a double CME. The first external breakout reconnection apparently releases an underlying sheared core field and allows it to erupt, leading to internal reconnection in the wake of the erupting helix. This internal reconnection heats the two-ribbon flare and might or might not produce the helix. These events lead to the first CME and are followed by a second breakout that initiates a second and larger halo CME. The strong magnetic shear in the region is associated with rapid proper motion and evolution of the active region. The multiple-turn helix originates from above a sheared-field magnetic inversion line within a filament channel, and starts to erupt only after fast breakout reconnection has started. These observations are counter to the standard flare model and support the breakout model for eruptive flare initiation. However, the observations are compatible with internal reconnection in a sheared magnetic arcade in the formation and eruption of the helix.

Gary, G. Allen↗

Eruption of a Multiple-Turn Helical Magnetic Flux Tube in a Large Flare: Evidence for External and Internal Reconnection that Fits the Breakout Model of Solar Magnetic Eruptions

We present observations and an interpretation of a unique multiple-turn spiral flux tube eruption from active region 10030 on 2002 July 15. The TRACE C IV observations clearly show a flux tube that is helical and erupting from within a sheared magnetic field. These observations are interpreted in the context of the breakout model for magnetic field explosions. The initiation of the helix eruption. as determined by a linear backward extrapolation, starts 25 s after the peak of the flare's strongest impulsive spike of microwave gyrosynchrotron radiation early in the flare s explosive phase, implying that the sheared core field is not the site of the initial reconnection. Within the quadrupolar configuration of the active region, the external and internal reconnection sites are identified in each of two consecutive eruptive flares that produce a double coronal mass ejection (CME). The first external breakout reconnection apparently releases an underlying sheared core field and allows it to erupt, leading to internal reconnection in the wake of the erupting helix. This internal reconnection releases the helix and heats the two-ribbon flare. These events lead to the first CME and are followed by a second breakout that initiates a second and larger halo CME. The strong magnetic shear in the region is compatible with the observed rapid proper motion and evolution of the active region. The multiple-turn helix originates from above a sheared-field magnetic inversion line within a filament channel. and starts to erupt only after fast breakout reconnection has started. These observations are counter to the standard flare model and support the breakout model for eruptive flare initiation.

Gary, G. Allen↗

International Space Station (ISS) Emergency Mask (EM) Development

The Emergency Mask (EM) is considered a secondary response emergency Personal Protective Equipment (PPE) designed to provide respiratory protection to the International Space Station (ISS) crewmembers in response to a post-fire event or ammonia leak. The EM is planned to be delivered to ISS in 2012 to replace the current air purifying respirator (APR) onboard ISS called the Ammonia Respirator (AR). The EM is a one ]size ]fits ]all model designed to fit any size crewmember, unlike the APR on ISS, and uses either two Fire Cartridges (FCs) or two Commercial Off-the-Shelf (COTS) 3M(Trademark). Ammonia Cartridges (ACs) to provide the crew with a minimum of 8 hours of respiratory protection with appropriate cartridge swap ]out. The EM is designed for a single exposure event, for either post ]fire or ammonia, and is a passive device that cannot help crewmembers who cannot breathe on their own. The EM fs primary and only seal is around the wearer fs neck to prevent a crewmember from inhaling contaminants. During the development of the ISS Emergency Mask, several design challenges were faced that focused around manufacturing a leak free mask. The description of those challenges are broadly discussed but focuses on one key design challenge area: bonding EPDM gasket material to Gore(Registered Trademark) fabric hood.

Toon, Katherine P.↗

Computational Bayesian Methods Applied to Complex Problems in Bio and Astro Statistics

In this dissertation we apply computational Bayesian methods to three distinct problems. In the first chapter, we address the issue of unrealistic covariance matrices used to estimate collision probabilities. We model covariance matrices with a Bayesian Normal-Inverse-Wishart model, which we fit with Gibbs sampling. In the second chapter, we are interested in determining the sample sizes necessary to achieve a particular interval width and establish non-inferiority in the analysis of prevalences using two fallible tests. To this end, we use a third order asymptotic approximation. In the third chapter, we wish to synthesize evidence across multiple domains in measurements taken longitudinally across time, featuring a substantial amount of structurally missing data, and fit the model with Hamiltonian Monte Carlo in a simulation to analyze how estimates of a parameter of interest change across sample sizes.

Elrod, Chris↗