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,441 records · Page 80

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

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

79 ASTRONOMY AND ASTROPHYSICS↗

A critical assessment of solutions to the galaxy diversity problem

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

79 ASTRONOMY AND ASTROPHYSICS↗

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

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

79 ASTRONOMY AND ASTROPHYSICS↗

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

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

08 HYDROGEN↗

Backus-Gilbert inversion of travel time data

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

Johnson, L. E.↗

NTF model concept for the X-29A. Part 2: Design and fabrication concept

The first concept was done to determine if a small model to fit the NASA Langley sizing criteria could accommodate aerodynamic and mechanical requirements. For a model wing span of 1.7 ft (and a length of 3 ft), the model is 6.25 percent (1/16 scale). The results were encouraging enough for Grumman to request, in early December 1981, that the design concept be refined at NASA Langley, where up-to-date cryogenic model design information was available. This working session was completed by a Grumman designer by mid-December. In its present form, the concept will be the basis for any continued effort.

Toscano, G.↗

Methods to Improve the Maintenance of the Earth Catalog of Satellites During Severe Solar Storms

The objective of this thesis is to investigate methods to improve the ability to maintain the inventory of orbital elements of Earth satellites during periods of atmospheric disturbance brought on by severe solar activity. Existing techniques do not account for such atmospheric dynamics, resulting in tracking errors of several seconds in predicted crossing time. Two techniques are examined to reduce of these tracking errors. First, density predicted from various atmospheric models is fit to the orbital decay rate for a number of satellites. An orbital decay model is then developed that could be used to reduce tracking errors by accounting for atmospheric changes. The second approach utilizes a Kalman filter to estimate the orbital decay rate of a satellite after every observation. The new information is used to predict the next observation. Results from the first approach demonstrated the feasibility of building an orbital decay model based on predicted atmospheric density. Correlation of atmospheric density to orbital decay was as high as 0.88. However, it is clear that contemporary: atmospheric models need further improvement in modeling density perturbations polar region brought on by solar activity. The second approach resulted in a dramatic reduction in tracking errors for certain satellites during severe solar Storms. For example, in the limited cases studied, the reduction in tracking errors ranged from 79 to 25 percent.

Wilkin, Paul G.↗

Limb Darkening and Planetary Transits: Testing Center-to-limb Intensity Variations and Limb-darkening Directly from Model Stellar Atmospheres

The transit method, employed by Microvariability and Oscillation of Stars ( MOST ), Kepler , and various ground-based surveys has enabled the characterization of extrasolar planets to unprecedented precision. These results are precise enough to begin to measure planet atmosphere composition, planetary oblateness, starspots, and other phenomena at the level of a few hundred parts per million. However, these results depend on our understanding of stellar limb darkening, that is, the intensity distribution across the stellar disk that is sequentially blocked as the planet transits. Typically, stellar limb darkening is assumed to be a simple parameterization with two coefficients that are derived from stellar atmosphere models or fit directly. In this work, we revisit this assumption and compute synthetic planetary-transit light curves directly from model stellar atmosphere center-to-limb intensity variations (CLIVs) using the plane-parallel Atlas and spherically symmetric SAtlas codes. We compare these light curves to those constructed using best-fit limb-darkening parameterizations. We find that adopting parametric stellar limb-darkening laws leads to systematic differences from the more geometrically realistic model stellar atmosphere CLIV of about 50–100 ppm at the transit center and up to 300 ppm at ingress/egress. While these errors are small, they are systematic, and they appear to limit the precision necessary to measure secondary effects. Our results may also have a significant impact on transit spectra.

79 ASTRONOMY AND ASTROPHYSICS↗

A small-box approach to the local crystal structure of Y 3 NbO 7

The local crystal structure of Y 3 NbO 7 was probed using neutron total scattering. Five different structural models were fit to the experimental pair distribution function of this material up to 15 Å. Fits of defect fluorite (Fm$\overline{3}$m), pyrochlore (Fd$\overline{3}$m), and orthorhombic C weberite models (C222 1 ) were attempted based on prior structural investigations of Y 3 NbO 7 and related fluorite-based materials. Orthorhombic P weberite (P222 1 ) and monoclinic weberite models (P112 1 ) were derived via symmetry reduction of the orthorhombic C weberite structure. All models except monoclinic weberite failed to reproduce the experimental pair distribution function. Misfits were clearly visible for peaks arising from atom pairs within the first coordination shell (defect fluorite and pyrochlore) and beyond (orthorhombic C and P weberite). By contrast, monoclinic weberite (γ ≈ 91.8°) provided an excellent fit, with no major misfits observed. This unit cell features distorted YO 8 dodecahedra, YO 7 monocapped octahedra, and NbO 6 octahedra as building blocks. Comparison between the topologies of the cationic and anionic substructures in monoclinic weberite and its parent orthorhombic C structure revealed displacements of the metal atoms along the [100] direction and a less compact spatial distribution of oxygen atoms along the a axis. In conclusion, the monoclinic weberite model presented herein may serve as a starting point to revisit structural modeling of fluorite-based materials with challenging local structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Parameterized yields of semivolatile products from isoprene oxidation under different NOx levels: impacts of chemical aging and wall-loss of reactive gases

We developed a new parameterization for the yields of semi-volatile/intermediate volatility products from isoprene oxidation that took into account the chemical aging under different NOx levels, including functionalization and fragmentation reactions within the gas phase, oligomerization and photolysis within the condensed phase, as well as the loss of particles and gases to chamber walls. We used a box model to fit the dynamic evolution of measured secondary organic aerosol (SOA) concentrations from chamber experiments. The explicit inclusion of multigenerational gas phase chemistry and SOA oligomerization/photolysis during fitting itself had large impacts on the fitted yields that relate to the initial volatility distribution of isoprene oxidation products. During further multigenerational aging of these initial organics, we varied the branching ratio between functionalization and fragmentation. The best fits suggested increasing amounts of fragmentation with increasing initial NOx levels for mixed experiments. Using a box model to simulate evolution of SOA in the atmosphere at constant dilution rate (in absence of wall losses), we found that the yields determined by including different chemical scenarios during fitting changed predicted SOA mass concentrations substantially over aging timescale of 1-2 days. Our results demonstrated the need for better constraining smog-chamber-derived volatility distributions of both gas- and particle-phase SOA constituents with measurements.

Xing, Li↗

Gaussian quadrature exponential sum modeling of near infrared methane laboratory spectra obtained at temperatures from 106 to 297 K

Transmission measurements made on near-infrared laboratory methane spectra have previously been fit using a Malkmus band model. The laboratory spectra were obtained in three groups at temperatures averaging 112, 188, and 295 K; band model fitting was done separately for each temperature group. These band model parameters cannot be used directly in scattering atmosphere model computations, so an exponential sum model is being developed which includes pressure and temperature fitting parameters. The goal is to obtain model parameters by least square fits at 10/cm intervals from 3800 to 9100/cm. These results will be useful in the interpretation of current planetary spectra and also NIMS spectra of Jupiter anticipated from the Galileo mission.

Giver, Lawrence P.↗

Estimating the reproductive number, total outbreak size, and reporting rates for Zika epidemics in South and Central America

As South and Central American countries prepare for increased birth defects from Zika virus outbreaks and plan for mitigation strategies to minimize ongoing and future outbreaks, understanding important characteristics of Zika outbreaks and how they vary across regions is a challenging and important problem. We developed a mathematical model for the 2015/2016 Zika virus outbreak dynamics in Colombia, El Salvador, and Suriname. We fit the model to publicly available data provided by the Pan American Health Organization, using Approximate Bayesian Computation to estimate parameter distributions and provide uncertainty quantification. The model indicated that a country-level analysis was not appropriate for Colombia. We then estimated the basic reproduction number to range between 4 and 6 for El Salvador and Suriname with a median of 4.3 and 5.3, respectively. We estimated the reporting rate to be around 16% in El Salvador and 18% in Suriname with estimated total outbreak sizes of 73,395 and 21,647 people, respectively. The uncertainty in parameter estimates highlights a need for research and data collection that will better constrain parameter ranges.

60 APPLIED LIFE SCIENCES↗

Tape-in-tape-out (TITO): a new approach for in-situ contact resistance measurements in high temperature superconducting CORC® cables

One of the ongoing development challenges with ReBCO high-temperature superconducting (HTS) cables is normal zone initiation and local heating, which is associated with over-critical current flowing through and around local performance reductions in individual tapes. Although inter-tape contact resistances are well-reported for individual tapes and HTS cables, these measurements are a challenge in CORC® cables as current percolates through much of the cable. In this work, developments in tape-in-tape-out (TITO) automated individual tape powering experiments are presented, and a simplified modeling approach for current percolation in CORC® cables is developed. An optimization is formulated to fit the model parameters to a large set of TITO experiments on a single cable, allowing the layer-dependent inter-tape contact resistance to be extracted. Measurements are presented and discussed for a straight CORC® cable and a cable bent to a 152 mm and 76 mm radius. The approach provides both quantitative insight and facilitates better understanding of current sharing phenomena in HTS cables.

Teyber, Reed↗

GPS measurements of crustal deformation associated with the 22 April 1991, Valle de la Estrella , Costa Rica earthquake

Results are presented for analysis of Global Positioning System (GPS) measurements made at sites in Costa Rica in February and July 1991. Significant horizontal and vertical displacements relative to February positions were observed. Differences were found in GPS derived vertical and horizontal displacements as compared to other types of geodetic measurements of uplift in the coastal regions. A slip dislocation model which fits the GPS measured displacement was computed. Differences between these data sets and their associated models were investigated in light of unmodeled slip heterogeneity on the fault and post-seismic displacements.

Lundgren, Paul R.↗

System Identification for eVTOL Aircraft Using Simulated Flight Data

This paper describes a system identification method for electric vertical takeoff and landing (eVTOL) aircraft. The approach merges fixed-wing and rotary-wing modeling techniques with new strategies to develop a modeling method for eVTOL vehicles using flight test data. The eVTOL aircraft system identification approach is demonstrated through application to the NASA LA-8 tandem tilt-wing, distributed electric propulsion aircraft using a high-fidelity flight dynamics simulation. Orthogonal phase-optimized multisine inputs are applied to each control surface and propulsor at numerous flight conditions throughout the flight envelope to collect informative flight data. An aero-propulsive model is identified at each flight condition using the equation-error method in the frequency domain. The local model parameters are then blended to create a global model across the nominal flight envelope. Parameter estimation results are shown to provide a good fit to modeling data and have good prediction capability. The methodology is developed with a discussion of unique eVTOL vehicle aerodynamic characteristics and practical strategies intended to inform future flight-based system identification efforts for eVTOL aircraft.

system identification↗

Development and Validation of an Initial Electric Motor Rotor Vibration Model

An overview of the development of a low fidelity prediction tool for electric motor rotor resonant modes and frequencies. These predictions are useful for the prediction of UAM motor noise during the design phase. Finite Element Analysis of motor rotors has led to the development of a parameterized model to predict these modes shapes. This has in turn allowed to a curve-fit based model from parametric sweeps. These curves can predict resonant frequencies and shapes for 1-5kW class motors without running a simulation.

acoustics↗

Rate-equation modelling and ensemble approach to extraction of parameters for viral infection-induced cell apoptosis and necrosis

We develop a theoretical approach that uses physiochemical kinetics modelling to describe cell population dynamics upon progression of viral infection in cell culture, which results in cell apoptosis (programmed cell death) and necrosis (direct cell death). Several model parameters necessary for computer simulation were determined by reviewing and analyzing available published experimental data. By comparing experimental data to computer modelling results, we identify the parameters that are the most sensitive to the measured system properties and allow for the best data fitting. Our model allows extraction of parameters from experimental data and also has predictive power. Using the model we describe interesting time-dependent quantities that were not directly measured in the experiment and identify correlations among the fitted parameter values. Numerical simulation of viral infection progression is done by a rate-equation approach resulting in a system of “stiff” equations, which are solved by using a novel variant of the stochastic ensemble modelling approach. The latter was originally developed for coupled chemical reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗