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,567 records · Page 87

Extragalactic interstellar extinction curves: Indicators of local physical conditions

Normalized interstellar extinction curves (ISECs) in the Milky Way and other galaxies show a variety of shapes. This variety is attributed to differences along different sight lines in the abundances of the several dust and gas components contributing to extinction. In this paper we propose that these abundance differences are not arbitrary but are a specific consequence of the physical conditions on those sight lines. If this proposal is correct, then it implies that ISECs contain information about physical conditions in the regions generating extinction. This may be particularly important for high redshift galaxies where information on the conditions may be difficult to obtain. We adopt a model of extinction carriers in which the solid and gaseous components are not immutable but respond time-dependently to the local physics. We validate this model by fitting extinction curves measured on sight lines in the Magellanic Clouds and obtained for the gamma-ray burst afterglow GRB 080605. We present results for this model as follows: (1) we show that computed ISECs are controlled by a small number of physical parameters, (2) we demonstrate the sensitivity of computed ISECs to these parameters, (3) we compute as examples ISECs for particular galaxy types, and (4) we note that different galaxy types have different shapes of ISEC.

79 ASTRONOMY AND ASTROPHYSICS↗

Confronting the concordance model of cosmology with Planck data

We confront the concordance (standard) model of cosmology, the spatially flat ΛCDM Universe with power-law form of the primordial spectrum with Planck CMB angular power spectrum data searching for possible smooth deviations beyond the flexibility of the standard model. The departure from the concordance cosmology is modeled in the context of Crossing statistic and statistical significance of this deviation is used as a measure to test the consistency of the standard model to the Planck data. Derived Crossing functions suggest the presence of some broad features in angular spectrum beyond the expectations of the concordance model. Our results indicate that the concordance model of cosmology is consistent to the Planck data only at 2 to 3σ confidence level if we allow smooth deviations from the angular power spectrum given by the concordance model. This might be due to random fluctuations or may hint towards smooth features in the primordial spectrum or departure from another aspect of the standard model. Best fit Crossing functions indicate that there are lack of power in the data at both low-ℓ and high-ℓ with respect to the concordance model. This hints that we may need some modifications in the foreground modeling to resolve the significant inconsistency at high-ℓ. However, presence of some systematics at high-ℓ might be another reason for the deviation we found in our analysis.

79 ASTRONOMY AND ASTROPHYSICS↗

Active Vibration Control of a Large Flexible Manipulator by Inertial Force and Joint Torque

The efficiency and positional accuracy of a lightweight flexible manipulator are limited by its flexural vibrations, which last after a gross motion is completed. The vibration delays subsequent operations. In the proposed work, the vibration is suppressed by inertial force of a small arm in addition to the joint actuators and passive damping treatment. The proposed approach is: (1) Dynamic modeling of a combined system, a large flexible manipulator and a small arm, (2) Determination of optimal sensor location and controller algorithm, and (3) Verification of the fitness of model and the performance of controller.

Lee, Soo Han↗

A New Method of Studying the Flow of the Water Along the Bottom of a Model of a Flying-boat Hull

A new method of studying the flow of the water along the bottom of a model of a flying-boat hull is described. In this method, the model is fitted with a transparent bottom and is divided down the center line by a bulkhead. The flow is observed and photographed through one-half of the model by means of the diffused illumination from a battery of lamps contained in the other half of the model. Photographs of the flow, particularly of the changes that occur when the step ventilates, are shown. The results of the present investigation indicate that the method has considerable promise, chiefly in connection with motion-picture studies.

Ward, Kenneth E↗

How microscopic epistasis and clonal interference shape the fitness trajectory in a spin glass model of microbial long-term evolution

The adaptive dynamics of evolving microbial populations takes place on a complex fitness landscape generated by epistatic interactions. The population generically consists of multiple competing strains, a phenomenon known as clonal interference. Microscopic epistasis and clonal interference are central aspects of evolution in microbes, but their combined effects on the functional form of the population’s mean fitness are poorly understood. Here, we develop a computational method that resolves the full microscopic complexity of a simulated evolving population subject to a standard serial dilution protocol. Through extensive numerical experimentation, we find that stronger microscopic epistasis gives rise to fitness trajectories with slower growth independent of the number of competing strains, which we quantify with power-law fits and understand mechanistically via a random walk model that neglects dynamical correlations between genes. We show that increasing the level of clonal interference leads to fitness trajectories with faster growth (in functional form) without microscopic epistasis, but leaves the rate of growth invariant when epistasis is sufficiently strong, indicating that the role of clonal interference depends intimately on the underlying fitness landscape. The simulation package for this work may be found at https://github.com/nmboffi/spin_glass_evodyn .

59 BASIC BIOLOGICAL SCIENCES↗

How microscopic epistasis and clonal interference shape the fitness trajectory in a spin glass model of microbial long-term evolution

The adaptive dynamics of evolving microbial populations takes place on a complex fitness landscape generated by epistatic interactions. The population generically consists of multiple competing strains, a phenomenon known as clonal interference. Microscopic epistasis and clonal interference are central aspects of evolution in microbes, but their combined effects on the functional form of the population’s mean fitness are poorly understood. Here, we develop a computational method that resolves the full microscopic complexity of a simulated evolving population subject to a standard serial dilution protocol. Through extensive numerical experimentation, we find that stronger microscopic epistasis gives rise to fitness trajectories with slower growth independent of the number of competing strains, which we quantify with power-law fits and understand mechanistically via a random walk model that neglects dynamical correlations between genes. We show that increasing the level of clonal interference leads to fitness trajectories with faster growth (in functional form) without microscopic epistasis, but leaves the rate of growth invariant when epistasis is sufficiently strong, indicating that the role of clonal interference depends intimately on the underlying fitness landscape. The simulation package for this work may be found at https://github.com/nmboffi/spin_glass_evodyn .

59 BASIC BIOLOGICAL SCIENCES↗

Earth as an Extrasolar Planet: Earth Model Validation Using EPOXI Earth Observations

The EPOXI Discovery Mission of Opportunity reused the Deep Impact flyby spacecraft to obtain spatially and temporally resolved visible photometric and moderate resolution near-infrared (NIR) spectroscopic observations of Earth. These remote observations provide a rigorous validation of whole disk Earth model simulations used to better under- stand remotely detectable extrasolar planet characteristics. We have used these data to upgrade, correct, and validate the NASA Astrobiology Institute s Virtual Planetary Laboratory three-dimensional line-by-line, multiple-scattering spectral Earth model (Tinetti et al., 2006a,b). This comprehensive model now includes specular reflectance from the ocean and explicitly includes atmospheric effects such as Rayleigh scattering, gas absorption, and temperature structure. We have used this model to generate spatially and temporally resolved synthetic spectra and images of Earth for the dates of EPOXI observation. Model parameters were varied to yield an optimum fit to the data. We found that a minimum spatial resolution of approx.100 pixels on the visible disk, and four categories of water clouds, which were defined using observed cloud positions and optical thicknesses, were needed to yield acceptable fits. The validated model provides a simultaneous fit to the Earth s lightcurve, absolute brightness, and spectral data, with a root-mean-square error of typically less than 3% for the multiwavelength lightcurves, and residuals of approx.10% for the absolute brightness throughout the visible and NIR spectral range. We extend our validation into the mid-infrared by comparing the model to high spectral resolution observations of Earth from the Atmospheric Infrared Sounder, obtaining a fit with residuals of approx.7%, and temperature errors of less than 1K in the atmospheric window. For the purpose of understanding the observable characteristics of the distant Earth at arbitrary viewing geometry and observing cadence, our validated forward model can be used to simulate Earth s time dependent brightness and spectral properties for wavelengths from the far ultraviolet to the far infrared.brightness

Robinson, Tyler D.↗

Impact of Pulsar and Fallback Sources on Multifrequency Kilonova Models

We explore the impact of pulsar electromagnetic dipole and fallback accretion emission on the luminosity of a suite of kilonova models. The pulsar models are varied over pulsar magnetic field strength, pulsar lifetime, ejecta mass, and elemental abundances; the fallback models are varied over fallback accretion rate and ejecta mass. For the abundances, we use Fe and Nd as representatives of the wind and dynamical ejecta, respectively. We simulate radiative transfer in the ejecta in either 1D spherical or 2D cylindrical spatial geometry. For the grid of 1D simulations, the mass fraction of Nd is 0, 10 –4 , or 10 –3 , and the rest is Fe. Our models that fit the bolometric luminosity of AT 2017gfo (the kilonova associated with the first neutron star merger discovered in gravitational waves, GW170817) do not simultaneously fit the B , V , and I time evolution. Here, we find that the trends of the evolution in B and V magnitudes are better matched by the fallback model relative to the pulsar model, implying the time-dependence of the remnant source influences the color evolution. Further exploration of the parameter space and model deficiencies is needed before we can describe AT 2017gfo with a remnant source.

79 ASTRONOMY AND ASTROPHYSICS↗

Limb-darkening and the structure of the Jovian atmosphere

By observing the transit of various cloud features across the Jovian disk, limb-darkening curves were constructed for three regions in the 4.6 to 5.1 mu cm band. Several models currently employed in describing the radiative or dynamical properties of planetary atmospheres are here examined to understand their implications for limb-darkening. The statistical problem of fitting these models to the observed data is reviewed and methods for applying multiple regression analysis are discussed. Analysis of variance techniques are introduced to test the viability of a given physical process as a cause of the observed limb-darkening.

Newman, W. I.↗

Use of collateral information to improve LANDSAT classification accuracies

Methods to improve LANDSAT classification accuracies were investigated including: (1) the use of prior probabilities in maximum likelihood classification as a methodology to integrate discrete collateral data with continuously measured image density variables; (2) the use of the logit classifier as an alternative to multivariate normal classification that permits mixing both continuous and categorical variables in a single model and fits empirical distributions of observations more closely than the multivariate normal density function; and (3) the use of collateral data in a geographic information system as exercised to model a desired output information layer as a function of input layers of raster format collateral and image data base layers.

Strahler, A. H.↗

Modeling of G333.6-0.2 as a spherical H II region

The radio and IR observations of the H II region G333.6-0.2 are matched with a detailed spherical model with a density distribution which has a uniform-density core of radius 0.05 pc, a power-law intermediate zone, and a uniform-density halo. A stellar radiation field is required that is somewhat different from those predicted by available model atmospheres. Of the stellar models of Kurucz, the Teff = 34,000 K and log g = 3.5 model best fits the observed ratio of helium to hydrogen recombination lines. A good fit to all the observations in obtained with an S(2+) ionizing flux which is a factor of 20 times less than predicted by the Kurucz atmosphere. Current model atmospheres may not be appropriate because a single stable star with Teff = 34,000 K fails by at least an order of magnitude to produce the ionizing luminosity. The small radius of the constant-density core implies a short dynamical lifetime of approximately 5000 years for this spherical model; this may indicate that a nonspherical blister geometry is more appropriate for this H II region.

Rubin, R. H.↗

On the Cycling of 231 Pa and 230 Th in Benthic Nepheloid Layers

The naturally-occurring radionuclides protactinium-231 ( 231 Pa) and thorium-230 ( 230 Th) are produced at approximately uniform rates in the ocean and thought to be removed from the water column through a reversible exchange with settling particles. Recent measurements along the U.S. GEOTRACES North Atlantic transect (GA03) revealed two features which are at odds with current understanding about 231 Pa and 230 Th cycling in the ocean: (i) a sharp decrease in dissolved 231 Pa ( 231 Pa d ) and 230 Th ( 230 Th d ) activities with depth below 2000-4000 m and (ii) very high particulate 231 Pa ( 231 Pa p ) and 230 Th ( 230 Th p ) activities near the bottom, at a number of stations between the New England continental shelf and Bermuda. Concomitant measurements of light attenuation from beam transmissometry showed that both features occur in benthic nepheloid layers (BNLs), which suggests that these features may stem, at least partly, from the presence of resuspended sediment in the deep water column. Here we explore the behaviour of 231 Pa and 230 Th in BNLs by using (i) radionuclide, optical, and hydrographic data from the western segment of GA03 (west of Bermuda) and (ii) a simplified model of particle and radionuclide cycling that includes a lateral particle source. First, the BNLs observed at GA03 stations are characterized from measurements of the beam attenuation coefficient converted to particle concentrations. At all stations, particle concentrations below the clear water minimum were the highest in the bottom mixed layer, whose thickness ranged from 95 to 320 m, and decreased generally with height above the bottom. The thickness of strong BNLs varied from 482 to 1358 m and the vertical integral of particle concentration in excess to that at the clear water minimum varied from 1 x 10 4 to 2 x 10 6 mg m −2 , among different stations. Second, the particle-radionuclide model is fitted to data from stations GT11-04 (New England continental rise) and GT11-08 (Hatteras abyssal plain), where samples for radionuclide analyses were collected in the BNL. The model can reproduce simultaneously the increase of particle concentration with depth, the low 231 Pa d and 230 Th d in the BNLs, and the high 231 Pa p and 230 Th p near the bottom. According to the model, at heights less than about 300 m above the seafloor, the dissolved phase was set primarily by a balance between adsorption and desorption, with vertical turbulent mixing playing a secondary role, whilst the particulate phase behaved largely as a non-reactive constituent supplied laterally and transported vertically by particle settling and turbulent mixing. Sensitivity tests with the model suggest that lateral particle sources near continental slopes and similar reliefs can produce significant biases both in the 230 Th normalization method and in the interpretation of sediment 231 Pa/ 230 Th records. Our findings yield insights into the influence of sediment resuspension and transport on 231 Pa and 230 Th in the deep ocean and highlight the need for considering these processes in paleoceanographic applications.

Nepheloid layer↗

Unveiling Structural Heterogeneity and Imbalance of Gold Decahedral Nanoparticles using Four-dimensional Scanning Transmission Electron Microscopy

Multi-twinned structures have been observed in technologically important crystal systems, for example diamond cubic and face-centered cubic (FCC) lattices, that include materials such as diamond, silicon, a wide range of noble metals, and their nanoscale counterparts. Beyond atomic building blocks, the special arrangements also occur in the self-assembly of nanoparticles (NP) and μm-sized colloidal particles and occupy parts of their phase diagrams. Spanning a wide range of length scales, the universality of the structures arises when the systems attempt to achieve multitwinned structures by overcoming geometric misfits during minimizing surface energies with entirely {111} or close-packing facets. While it is fundamental to understand how strain is sustained upon twinned structures and symmetry breaking, the knowledge will be paramount in practical aspects such as guiding and controlling the thin film growth, anisotropic NP growth, and self-assembly of NPs. Au decahedral (Dh) NP, as the most prevalent multi-twinned model system, fits five tetrahedral motifs into a circle by sharing an axis resulting in a geometric misfit angle of 7.35°, or a disclination with power of -7.35°. Postulating how Au FCC lattice adopts the misfit, theoretical models have been developed to address the underlying lattice symmetry and inhomogeneous strain distribution separately. Yet, experimental reports regarding the former have been limited due to the relatively large X-ray beam sizes that do not fit the sizes of NPs. On the other hand, though the latter has been widely adapted in the thermodynamics of small (<10 nm) multi-twinned nanoparticles, previous literature has shown that, at edge length of 17 nm, the theory’s is invalidated by shear strain that is observed in a defect-free Au Dh NP by high-resolution transmission electron microscopy (HR-TEM) imaging. Though the advancement of aberrationcorrected scanning transmission electron microscopy (AC-STEM) imaging and ab initio calculation techniques brings new opportunities, along with challenges in complicated image analysis and limitation in particle size (usually below 10 nm), the gap between nanoscale and mesoscale has never been extended to gain insight from atomic system with straightforward interaction potentials.

4D-STEM↗

Method and models for R-curve instability calculations

This paper presents a simple method for performing elastic R-curve instability calculations. For a single material-structure combination, the calculations can be done on some pocket calculators. On microcomputers and larger, it permits the development of a comprehensive program having libraries of driving force equations for different configurations and R-curve model equations for different materials. The paper also presents several model equations for fitting to experimental R-curve data, both linear elastic and elastoplastic. The models are fit to data from the literature to demonstrate their viability.

Orange, Thomas W.↗

Method and models for R-curve instability calculations

This paper presents a simple method for performing elastic R-curve instability calculations. For a single material-structure combination, the calculations can be done on some pocket calculators. On microcomputers and larger, it permits the development of a comprehensive program having libraries of driving force equations for different configurations and R-curve model equations for different materials. The paper also presents several model equations for fitting to experimental R-curve data, both linear elastic and elastoplastic. The models are fit to data from the literature to demonstrate their viability.

Orange, Thomas W.↗

Recommended DSMC Collision Model Parameters for Planetary Entry

Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity simulations of environments consisting of both continuum and rarified flow. DSMC is a stochastic, particle-based method which solves the fundamental Boltzmann equation and is therefore necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally computationally inefficient to model the continuum regime. In a CFD/DSMC hybrid approach, obtaining an accurate, high-fidelity solution hinges on the consistent treatment of transport properties and the thermo-chemical models employed within the two solvers. In principle, in regions where CFD and DSMC are both employed, the same gas mixture under the same conditions should have the same properties, regardless of simulation type. Observed differences should be due to non-equilibrium processes, rather than differences in physical models. The goal of this work is to provide a comprehensive DSMC transport database encompassing the vast majority of known gas species encountered during Earth or other planetary atmospheric entry. This goal is largely possible due to recently performed ab initio quantum chemistry calculations. Combined with other high-fidelity data, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 20000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute collision-specific Variable Soft Sphere (VSS) collision model parameters, fitted from 300 K to 20000 K. Initial comparisons of transport properties of relevant equilibrium gas mixtures show great agreement between CFD and DSMC-derived results. The completed database can be readily applied to model binary collisions of any gas mixture containing the included species over the specified temperature range, making it a valuable tool for future planetary probe modeling efforts. An example is shown below in Fig. 1. Equilibrium mixture transport properties for a 35-species mixture composed originally of 10% air and 90% pyrolysis species of a carbon-phenolic ablator material [4] are computed using both fitted VSS Parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches are less than 5% for the entire temperature range.

M. R. Gosma↗

Understanding Shock Dynamics in the Inner Heliosphere with Modeling and Type II Radio Data: the 2010-04-03 Event

The 2010 April 03 solar event was studied using observations from STEREO SECCHI, SOHO LASCO, and Wind kilometric Type II data (kmTII) combined with WSA-Cone-ENLIL model simulations performed at the Community Coordinated Modeling Center (CCMC). In particular, we identified the origin of the coronal mass ejection (CME) using STEREO EUVI and SOHO EIT images. A flux-rope model was fit to the SECCHI A and B, and LASCO images to determine the CMEs direction, size, and actual speed. J-maps from STEREO COR2HI-1HI-2 and simulations fromCCMC were used to study the formation and evolution of the shock in the inner heliosphere. In addition, we also studied the time-distance profile of the shock propagation from kmTII radio burst observations. The J-maps together with in-situ datafrom the Wind spacecraft provided an opportunity to validate the simulation results andthe kmTII prediction. Here we report on a comparison of two methods of predictinginterplanetary shock arrival time: the ENLIL model and the kmTII method; andinvestigate whether or not using the ENLIL model density improves the kmTIIprediction. We found that the ENLIL model predicted the kinematics of shock evolutionwell. The shock arrival times (SAT) and linear-fit shock velocities in the ENLILmodel agreed well with those measurements in the J-maps along both the CME leading edge and the Sun-Earth line. The ENLIL model also reproduced most of the largescale structures of the shock propagation and gave the SAT prediction at Earth with an error of 17 hours. The kmTII method predicted the SAT at Earth with an error of 15 hours when using n0 4.16 cm3, the ENLIL model plasma density near Earth; but itimproved to 2 hours when using n0 6.64 cm3, the model density near the CMEleading edge at 1 AU.

Space weather forecast↗

Metrics as tools for bridging climate science and applications

In climate science and applications, the term “metric” is used to describe the distillation of complex, multifaceted evaluations to summarize the overall quality of a model simulation, or other data product, and/or as a means to quantify some response to climate change. Metrics provide insights into the fidelity of processes and outcomes from climate models and can assist with both differentiating models' representation of variables or processes and informing whether models are “fit for purpose.” Metrics can also provide a valuable reference point for co-production of knowledge between climate scientists and climate impact practitioners. Although continued metric developments enable model developers to better understand the impacts of decisions made in the model design process, metrics also have implications for the characterization of uncertainty and facilitating analyses of underlying physical processes. As a result, comprehensive evaluation with multiple metrics enhances usability of climate information by both scientific and stakeholder communities. Here, this paper presents examples of insights gained from the development and appropriate use of metrics, and provides examples of how metrics can be used to engage with stakeholders and inform decision-making.

54 ENVIRONMENTAL SCIENCES↗