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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,063 records · Page 59

All sources NASA NTRS DOE OSTI

NASA NTRS1999

Capella: Structure and Abundances

This grant covers the analysis of EUVE spectra of the cool star binary system Capella. This project has also required the analysis of simultaneous Advanced Satellite for Cosmology and Astrophysics (ASCA) data. The ASCA spectrum of Capella could not be fit with standard models; by imposing models based on strong lines observed with EUVE, a problem wavelength region was identified. Correcting the problem required calculations of atomic collision strengths of higher principal quantum number than had ever been calculated. With these new models applied to the ASCA spectrum, better fits were obtained. Findings are that: (1) ASCA and EUVE spectra are both dominated by a region at 6 x 10(exp 6) K. (2) The high energy cut-off of the ASCA spectrum is consistent with emission from the highest ionization stages of EUVE, namely Fe XXIV. (3) EUVE requires a continuous emission measure distribution with more than two temperatures. (4) The ASCA spectra are of such high statistical significance that systematic uncertainties dominate, including atomic physics issues and calibration issues. (5) While the ASCA spectral fits achieve lower Chi(exp 2 with two-temperature fits, the EUVE-derived emission measure distribution models are also consistent with the spectra. (6) The Fe/H ratio obtained from the ASCA fit is within 20 % of the Fe/H abundance obtained from the summed spectra of Capella over 5 EUVE pointings, as well as the 1996 EUVE data. This result confirms our claims that quasi-continua composed of weak emission lines in the short wavelength spectrometer of EUVE are not major contributors to the measured Capella continuum. Other abundance ratios are also determined from the ASCA data, using models derived with EUVE. Si, Si, and Mg appear to be close to solar photospheric values, while the ratio of Ne/Fe is three to four times lower than solar photospheric values. Whether there is a general First Ionization Potential (FIP) effect or a specific neon anomaly cannot be determined from these data. (8) EUVE line ratios indicating low optical depth place limits on the extent to which the X-ray lines are optically, thick. While the brightest lines may be marginally thick, the ASCA spectrum is not of sufficient resolution to determine the optical depth in these lines.

Brickhouse, Nancy S.↗

NASA NTRS1992

The remarkable twisted disk of NGC 4753 and the shapes of galactic halos

The complex dust lanes in the S0 galaxy NGC 4753 are shown to be consistent with a disk that is strongly twisted by differential precession. Yet another peculiar S0 can therefore be explained as the result of an accretion event. An evolving disk model is fitted to the observed distribution. This disk is inclined by 15 deg relative to the galaxy's equatorial plane and twisted such that its line of nodes changes smoothly by 3.8 pi over a factor of seven in radius. The model indicates that most of the galaxy's mass is unseen, is nearly spherically distributed, and has a nearly scale-free spatial distribution. The ellipticity of the total galactic mass distribution must be constant to within 20 percent over the radial extent of the twisted disk - a conclusion which may significantly constrain galaxy formation mechanisms.

Steiman-Cameron, Thomas Y.↗

NASA NTRS1995

Observations and physical interpretations of the solar wind flow properties as obtained from white light coronagraph aboard SPARTAN 201-01

The solar corona was observed with an externally occulted White Light Coronagraph (WLC) carried on the SPARTAN 201-1 spacecraft on 11-12 Apr. 1993. With observations from WLC and the ground based Mauna Loa White Light Coronagraph, a large number of polar plumes both in the north and south polar holes were traced from 1.16 to 5.5 Rs. Flow properties of the solar wind in coronal holes have been determined (Habbal et al., 1995) by using a two fluid model constrained by density profiles and scale height temperatures from the white light observations, and interplanetary measurements of the flow speed and proton mass flux from Ulysses' south polar passage. Provisions for acceleration by Alfven waves, as well as electron and proton heating, are included in the momentum and the energy equations respectively. The model computations fit remarkably well the empirical constraints of the two different density structures (plumes and coronal holes) for a range of input parameters. In this study we investigate the physical nature of the heating function used in the two-fluid model. Alfven waves have been suggested as the possible source of heating that accelerates the solar wind (Ofman and Davila, 1995). We utilize the density contrast observed in WLC data in the plume and ambient coronal hole region to estimate the Alfven wave frequencies responsible for heating these structures. The source heating function utilized in the two fluid model of the solar wind acceleration will be compared with the resonant Alfven wave heating function.

Guhathakurta, Madhulika↗

NASA NTRS2014

Dust and Gas in the Magellanic Clouds from the Heritage Herschel Key Project. I. Dust Properties and Insights into the Origin of the Submm (Submillimeter) Excess Emission

The dust properties in the Large and Small Magellanic Clouds are studied using the HERITAGE Herschel Key Project photometric data in five bands from 100 to 500 micromillimeters. Three simple models of dust emission were fit to the observations: a single temperature blackbody modified by a powerlaw emissivity (SMBB), a single temperature blackbody modified by a broken power-law emissivity (BEMBB), and two blackbodies with different temperatures, both modified by the same power-law emissivity (TTMBB). Using these models we investigate the origin of the submillimeter excess; defined as the submillimeter (submm) emission above that expected from SMBB models fit to observations < 200 micromillimeters. We find that the BEMBB model produces the lowest fit residuals with pixel-averaged 500 micromillimeters submillimeter excesses of 27% and 43% for the Large and Small Magellanic Clouds, respectively. Adopting gas masses from previous works, the gas-to-dust ratios calculated from our fitting results show that the TTMBB fits require significantly more dust than are available even if all the metals present in the interstellar medium (ISM) were condensed into dust. This indicates that the submillimeter excess is more likely to be due to emissivity variations than a second population of colder dust. We derive integrated dust masses of (7.3 plus or minus 1.7) x 10 (sup 5) and (8.3 plus or minus 2.1) x 10 (sup 4) solar masses for the Large and Small Magellanic Clouds, respectively. We find significant correlations between the submillimeter excess and other dust properties; further work is needed to determine the relative contributions of fitting noise and ISM physics to the correlations.

Magellanic↗

NASA NTRS2020

Bayesian Statistical Models for Community Annoyance Survey Data

This paper demonstrates the use of two Bayesian statistical models to analyze single-event sonic boom exposure and human annoyance data from community response surveys. Each model is fit to data from a NASA pilot study.Unlike many community noise surveys, this study used a panel sample to collect multiple observations per participant instead of a single observation. Thus, a multilevel (also known as hierarchical or mixed-effects) model is used to account for the within-subject correlation in the panel sample data. This paper describes a multilevel logistic regression model and a multilevel ordinal regression model. The paper also proposes a method for calculating a summary dose-response curve from the multilevel models that represents the population. The two models’ summary dose-response curves are visually similar. However, their estimates differ when calculating the noise dose at a fixed percent highly annoyed.

Musical instruments↗

DOE OSTI2020

Predicting molecular dipole moments by combining atomic partial charges and atomic dipoles

The molecular dipole moment ( μ ) is a central quantity in chemistry. It is essential in predicting infrared and sum-frequency generation spectra as well as induction and long-range electrostatic interactions. Furthermore, it can be extracted directly—via the ground state electron density—from high-level quantum mechanical calculations, making it an ideal target for machine learning (ML). Here, we choose to represent this quantity with a physically inspired ML model that captures two distinct physical effects: local atomic polarization is captured within the symmetry-adapted Gaussian process regression framework which assigns a (vector) dipole moment to each atom, while the movement of charge across the entire molecule is captured by assigning a partial (scalar) charge to each atom. The resulting “MuML” models are fitted together to reproduce molecular μ computed using high-level coupled-cluster theory and density functional theory (DFT) on the QM7b dataset, achieving more accurate results due to the physics-based combination of these complementary terms. The combined model shows excellent transferability when applied to a showcase dataset of larger and more complex molecules, approaching the accuracy of DFT at a small fraction of the computational cost. We also demonstrate that the uncertainty in the predictions can be estimated reliably using a calibrated committee model. The ultimate performance of the models—and the optimal weighting of their combination—depends, however, on the details of the system at hand, with the scalar model being clearly superior when describing large molecules whose dipole is almost entirely generated by charge separation. These observations point to the importance of simultaneously accounting for the local and non-local effects that contribute to μ ; furthermore, they define a challenging task to benchmark future models, particularly those aimed at the description of condensed phases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1990

Universal adsorption at the vapor-liquid interface near the consolute point

The ellipticity of the vapor-liquid interface above mixtures of methylcyclohexane (C7H14) and perfluoromethylcyclohexane (C7F14) has been measured near the consolute point T(c) = 318.6 K. The data are consistent with a model of the interface that combines a short-ranged density-vs height profile in the vapor phase with a much longer-ranged composition-versus-height profile in the liquid. The value of the free parameter produced by fitting the model to the data is consistent with results from two other simple mixtures and a mixture of a polymer and solvent. This experiment combines precision ellipsometry of the vapor-liquid interface with in situ measurements of refractive indices of the liquid phases, and it precisely locates the consolute point.

Schmidt, James W.↗

NASA NTRS1994

Ion transport and loss in the Earth's quiet ring current. 2: Diffusion and magnetosphere-ionosphere coupling

We have studied the transport and loss of H(+), He(+), and He(++) ions in the Earth's quiet time ring current (1 to 300 keV/e, 3 to 7 R(sub E), Kp less than 2+, absolute value of Dst less than 11, 70 to 110 degs pitchangles, all LT) comparing the standard radial diffusion model developed for the higher-energy radiation belt particles with measurements of the lower energy ring current ions in a previous paper. Large deviations of that model, which fit only 50% of the data to within a factor of 10, suggested that another transport mechanism is operating in the ring current. Here we derive a modified diffusion coefficient corrected for electric field effects on ring current energy ions that fit nearly 80% of the data to within a factor of 2. Thus we infer that electric field fluctuations from the low-latitude to midlatitude ionosphere (ionospheric dynamo) dominated the ring current transport, rather than high-latitude or solar wind fluctuations. Much of the remaining deviation may arise from convective electric field transport of the E less than 30 keV particles. Since convection effects cannot be correctly treated with this azimuthally symmetric model, we defer treatment of the lowest-energy ions to a another paper. We give chi(exp 2) contours for the best fit, showing the dependence of the fit upon the internal/external spectral power of the predicted electric and magnetic field fluctuations.

Sheldon, R. B.↗

NASA NTRS2002

A Hall Thruster Performance Model Incorporating the Effects of a Multiply-Charged Plasma

A Hall thruster performance model that predicts anode specific impulse, anode efficiency, and thrust is discussed. The model is derived as a function of a voltage loss parameter, an electron loss parameter, and the charge state of the plasma. Experimental data from SPT and TAL type thrusters up to discharge powers of 21.6 kW are used to determine the best fit for model parameters. General values for the model parameters are found, applicable to high power thrusters and irrespective of thruster type. Performance of a 50 kW thruster is calculated for an anode specific impulse of 2500 seconds or a discharge current of 100 A.

Hofer, Richard R.↗

NASA NTRS1982

IUE spectroscopy, visible-band polarimetry, and radiometry of V641 Mon

This hot, double line, ellipsoidal variable member of NGC 2264 has been shown previously to be either a semi-detached or contact close binary. Low-resolution IUE spectra are best fitted to a Kurucz model atmosphere for very small (approximately 0.08 mag) E(B-V). The familiar interstellar absorption dip near lambda 2200 is apparently absent. A suitable model atmosphere can be fitted to the IUE fluxes, but flux excesses (compared to the model) appear for all the published U through L magnitudes. The spectrum of the B through L excess appears to follow a .0001 lambda dependence. It is shown that this cannot be interpreted as arising from another star fortuitously observed in the visible band or IR. Ground based polarization measures indicate V641 Mon to be a polarization variable. Previous and new V light curves show the amplitude of light variability itself to be variable by about a factor of 2. It is suggested that all these observed characteristics are best explained by postulating "third light" and identifying part of it with Rayleigh scattered starlight very near to the stars. From this same region there arise circumstellar absorptions which give rise to nontheoretical strengths for Si II and Si III lines.

Koch, R. H.↗

DOE OSTI2026

Reformulating a Gurson-based dynamic damage model and demonstrating improved predictive power and numerical robustness

The original Tepla (TEnsile PLAsticity) ductile damage model, based on the Gurson yield surface, has long been used to model damage evolution and material failure under dynamic loading. Unfortunately, Tepla suffered from mesh sensitivity, numerical instability, and limited predictive capability. Here, we theoretically reformulate Tepla to address these issues. We especially focus on the prediction of porosity, which is the key state variable used for modeling ductile damage, as compared to more easily measured surface velocities, which are at best an indirect measure of damage. Key model changes include separating the viscosity during volumetric void growth from underlying shear strength behavior and switching to an iterative bisection solver. The new Tepla is then calibrated on incipient spall experiments on half-hard copper and tantalum, which demonstrate its ability to simultaneously fit the model to recovered porosity distributions and measured surface velocities, a stringent test. Improved numerical behavior, such as greatly reduced mesh sensitivity, is also shown in those simulations. Finally, the new Tepla model is applied to several high-explosive loaded, sweeping wave experiments, showing the ability of the model to predict behavior on tests with significantly different loading conditions and histories than the calibration data.

36 MATERIALS SCIENCE↗

NASA NTRS1985

Dynamics of liquid sloshing in upright and inverted bladdered tanks

A system identification methodology was used to examine the dynamics of liquid sloshing in the upright and inverted bladdered hydrazine tanks of the Tracking and Data Relay Satellite, (TDRS) and to evaluate the effects of bladder stiffness on the sloshing parameters. Mechanical models of the two systems were developed using the numerical values derived from static stability tests and from slosh frequency response tests of a full-size model tank fitted with a prototype bladder. For the upright tank (liquid below the bladder) a modified conventional pendulum was used. In the inverted tank (liquid above the bladder) where sloshing is unconventional due to the highly nonsymmetrical orientation of the liquid held by the bladder, a mechanical model using an inverted pendulum which is able to undergo small oscillations as well as large reorientations was necessary. Both thrusting and low-gravity conditions are considered.

Dodge, F. T.↗

DOE OSTI2026

Dynamic crushing of metal lattice metamaterials: Shock mode diagrams and transition to topology-independent compaction regime

Additively manufactured lattice metamaterials offer design versatility in strength and energy absorption and provide an additional degree of freedom through the selection of the lattice topology. Under quasistatic loading, the unit cell structure can strongly affect the stiffness, yield, and post-yield behavior, but whether and to what degree the effect of lattice topology persists into dynamic loading scenarios, up to the compaction shock regime, has not been established. LLNL’ s ALE3D hydrocode was used to perform a computational investigation of dynamic loading in multiple lattice types, including the gyroid, octet, Schwarz D, and rhombic dodecahedron, under impact velocities from 0.25 to 2.25 km/s. Shock Hugoniots for each lattice topology are generated and compared, suggesting that above a critical velocity, distinctions between architectures may not persevere and compacted lattices behave similarly. Here, to investigate the transition between topology-dependent quasistatic compression and the topology-independent regime above the critical velocity, a one-dimensional elastic-linear hardening plasticity-densified solid (E-LHP-DS) shock model for lattice materials was developed that relies upon confined compression to link the quasistatic and shock mechanics. Unlike similar works, the model does not assume rigid behavior prior to yield or locking behavior at densification, allowing a richer exploration of lattice mechanics. With only six parameters, the analytical model simultaneously fit quasistatic confined compression simulations for relative densities 0.1 $≤ \bar{ρ} ≤$ 0.9 and predicted dynamic compaction behavior to traverse several distinct shock modes, each defined by a critical impact speed (equivalently, critical stresses). Comparing the numerical results to the one-dimensional E-LHP-DS shock model predictions suggests that the topology-independence under strong shocks is linked to the onset of densification, which can be predicted based on quasistatic confined compression results.

Cellular material↗

DOE OSTI2021

The Modified Void Nucleation and Growth Model (MNAG) for Damage Evolution in BCC Ta

A void coalescence term was proposed as an addition to the original void nucleation and growth (NAG) model to accurately describe void evolution under dynamic loading. The new model, termed as modified void nucleation and growth model (MNAG model), incorporated analytic equations to explicitly account for the evolution of the void number density and the void volume fraction (damage) during void nucleation, growth, as well as the coalescence stage. The parameters in the MNAG model were fitted to molecular dynamics (MD) shock data for single-crystal and nanocrystalline Ta, and the corresponding nucleation, growth, and coalescence rates were extracted. The results suggested that void nucleation, growth, and coalescence rates were dependent on the orientation as well as grain size. Compared to other models, such as NAG, Cocks–Ashby, Tepla, and Tonks, which were only able to reproduce early or later stage damage evolution, the MNAG model was able to reproduce all stages associated with nucleation, growth, and coalescence. The MNAG model could provide the basis for hydrodynamic simulations to improve the fidelity of the damage nucleation and evolution in 3-D microstructures.

36 MATERIALS SCIENCE↗

NASA NTRS2008

A Chapman-Layers Ionspheric Model for Mars

A numerical model (CLIMM) is developed that adopts functions of two Chapman layers to compute Mars ionospheric electron densities at given local solar zenith angle and height. Electron density profiles derived from Mars Global Survey (MGS)-to-Earth radio occultation measurements collected during 1998 through 2005 are used to fit the model. The present model does not include variations with solar extreme ultraviolet (EUV) radiation cycles and seasons, and may have increased errors at lower latitudes. A more sophisticated model taking into account these variations is being developed and will be available in the future.

solar extreme ultraviolet (EUV)↗

NASA NTRS2016

Characterization of the High-Albedo NEA 3691 Bede

Characterization of NEAs provides important inputs to models for atmospheric entry, risk assessment and mitigation. Diameter is a key parameter because diameter translates to kinetic energy in atmospheric entry. Diameters can be derived from the absolute magnitude, H(PA=0deg), and from thermal modeling of observed IR fluxes. For both methods, the albedo (pv) is important - high pv surfaces have cooler temperatures, larger diameters for a given Hmag, and shallower phase curves (larger slope parameter G). Thermal model parameters are coupled, however, so that a higher thermal inertia also results in a cooler surface temperature. Multiple parameters contribute to constraining the diameter. Observations made at multiple observing geometries can contribute to understanding the relationships between and potentially breaking some of the degeneracies between parameters. We present data and analyses on NEA 3691 Bede with the aim of best constraining the diameter and pv from a combination of thermal modeling and light curve analyses. We employ our UKIRT+Michelle mid-IR photometric observations of 3691 Bede's thermal emission at 2 phase angles (27&43 deg 2015-03-19 & 04-13), in addition to WISE data (33deg 2010-05-27, Mainzer+2011). Observing geometries differ by solar phase angles and by moderate changes in heliocentric distance (e.g., further distances produce somewhat cooler surface temperatures). With the NEATM model and for a constant IR beaming parameter (eta=constant), there is a family of solutions for (diameter, pv, G, eta) where G is the slope parameter from the H-G Relation. NEATM models employing Pravec+2012's choice of G=0.43, produce D=1.8 km and pv≈0.4, given that G=0.43 is assumed from studies of main belt asteroids (Warner+2009). We present an analysis of the light curve of 3691 Bede to constrain G from observations. We also investigate fitting thermophysical models (TPM, Rozitis+11) to constrain the coupled parameters of thermal inertia (Gamma) and surface roughness, which in turn affect diameter and pv. Surface composition can be related to pv. This study focuses on understanding and characterizing the dependency of parameters with the aim of constraining diameter, pv and thermal inertia for 3691 Bede.

Wooden, Diane H.↗

NASA NTRS2016

Characterization of the High-Albedo NEA 3691 Bede

magnitude↗

NASA NTRS1978

An energy-dependent map of the X-ray emission of the Perseus cluster

Results are discussed for X-ray scans of the Perseus cluster with the SAS 3 slat detectors. Counting rates in two broadband channels (1.5 to 5.5 keV and 5.5 to 15.5 keV) are compared as a function of the instrument scan angle measured from the center of emission, and it is found that the emission in the higher energy range is broader in extent. Tests for a systematic energy-dependent error in the background demonstrate that the discovered energy dependence of the Perseus cluster emission is valid. Iron line emission is determined by two methods, various models are fitted to the data, and the total cluster emission in the two energy ranges is computed. The iron line emission is determined to be 0.003 photon/sq cm per sec. It is concluded that the diffuse emission increases in the range from 1.5 to 15.5 keV within the bounds of a model involving a point source with a radius not exceeding about 3 arcmin plus a diffuse source.

Ulmer, M. P.↗