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

A Climatology and Extreme Value Analysis of Large Hail in China

In this paper, hail size records are analyzed at 2254 stations in China and a hail size climatology is developed based on gridded hail observations for the period 1960–2015. It is found that the annual percentiles of hail size records changed sharply and national-wide after 1980, therefore two periods, 1960–79 and 1980–2015, are studied. There are some similarities between the two periods in terms of the characteristics of hail size such as the spatial distribution patterns of mean annual maximum hail size and occurrence week of annual maximum hail size. The 1980–2015 period had higher observation density than the 1960–79 period, but showed smaller mean annual maximum hail size, especially in northern China. In the majority of grid boxes, the annual maximum hail size experienced a decreasing trend during the 1980–2015 period. A Gumbel extreme value model is fitted to each grid box to estimate the return periods of maximum hail size. The scale and location parameter of the fitted Gumbel distributions are higher in eastern China than in western China, thereby reflecting a greater likelihood of large hail in eastern China. In southern China, the maximum hail size exceeds 127 mm for a 10-yr return period, whereas in northern China maximum hail size exceeds this threshold for a 50-yr return period. Furthermore, the Gumbel model is found to potentially underestimate the maximum hail size for certain return periods, but provides a more informed picture of the spatial distribution of extreme hail size and the regional differences.

54 ENVIRONMENTAL SCIENCES↗

A spectral study of the persistent X-ray flux from 4 U/MXB 1636-53

A spectral analysis of the persistent X-ray flux from the X-ray burst source 4 U/MXB 1636-53, which was continuously observed with EXOSAT for about 25 h in May, 1984, is presented. It is found that no standard one-component spectral model provides an acceptable representation of all the data. However, almost all 'complex one-component' and 'two-component' models do accurately fit the data. A blackbody plus a 'multicolor' disk spectrum does not fit these data. In addition, some 'two-component' spectral models can be fitted to the data with more than one set of different parameter values. Based on these results, it is suggested that the various physical scenarios that have been proposed to interpret 'two-component' spectra of low-mass X-ray binaries may be based on mathematical parameterizations which may not provide unique representations of the data. If this is the case, the reported 'two-component' mathematical fits may not have any true physical significance.

Vacca, W. D.↗

Modeling the Earth's magnetospheric magnetic field confined within a realistic magnetopause

Empirical data-based models of the magnetosphereic magnetic field have been widely used during recent years. However, the existing models (Tsyganenko, 1987, 1989a) have three serious deficiencies: (1) an unstable de facto magnetopause, (2) a crude parametrization by the K(sub p) index, and (3) inaccuracies in the equatorial magnetotail B(sub z) values. This paper describes a new approach to the problem; the essential new features are (1) a realistic shape and size of the magnetopause, based on fits to a large number of observed crossing (allowing a parametrization by the solar wind pressure), (2) fully controlled shielding of the magnetic field produced by all magnetospheric current systems, (3) new flexible representations for the tail and ring currents, and (4) a new directional criterion for fitting the model field to spacecraft data, providing improved accuracy for field line mapping. Results are presented from initial efforts to create models assembled from these modules and calibrated against spacecraft data sets.

Tsyganenko, N. A.↗

Structure of the Bastnäsite (001) Surface by Crystal Truncation Rod X-ray Diffraction and Ab Initio Molecular Dynamics: Implications for Separations of a Rare Earth Ore Mineral

Bastnäsite ((Ce,La)FCO 3 ) is the primary mineral source of light rare earth elements, but its surface structure is not well understood. This presents a major challenge in improving beneficiation strategies. In this work, a synergistic combination of X-ray scattering and ab initio molecular dynamics (AIMD) was used to gain atomistic insight into the interfacial structure of bastnäsite. Surface X-ray scattering was used to measure crystal truncation rods (CTRs) of the bastnäsite (001) surface, a significant crystal face with a previously unknown termination. The best-fit atomic-scale model of the CTR data features a carbonate layer at the surface, which is stabilized by the relaxation of carbonate groups from their bulk structural positions. AIMD simulations predict similar surface relaxations, which are shown to be influenced by the protonation of oxygen atoms at the surface. Evidence of ordered water at the interface is also observed in the best-fit model and AIMD simulations. Further, the presence of a carbonate layer at this dominant crystal surface is significant for improving separation technologies because most commonly used ligands utilize anionic functional groups to chelate metal cations at particle surfaces. Without modification, anionic ligands are expected to have poor affinity for the carbonate-terminated (001) surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Models of the quasar population. II - The effects of dust obscuration

A model for dust obscuration is proposed to account for the apparent lack of high-redshift quasars which are assumed to be obscured by dust in intervening galaxies. The distribution function for the optical depth as a function of redshift is constructed and used to determine the expected properties of quasars. The apparent magnitude distribution at three points is used to fit the model. Optical results are given for the comoving number density, redshift distribution, and mean color excess. If the theory is correct, evidence for dust may be seen in quasar spectra and a correlation between quasar continuum slopes and strength of metal line absorption systems should be observed.

Heisler, Julia↗

Chandra Grating Spectroscopy of Three Hot White Dwarfs

High-resolution soft X-ray spectroscopic observations of single hot white dwarfs are scarce. With the Chandra Low-Energy Transmission Grating, we have observed two white dwarfs, one is of spectral type DA (LB1919) and the other is a non-DA of spectral type PG1159 (PG1520+525). The spectra of both stars are analyzed, together with an archival Chandra spectrum of another DA white dwarf (GD246). Aims. The soft X-ray spectra of the two DA white dwarfs are investigated in order to study the effect of gravitational settling and radiative levitation of metals in their photospheres. LB1919 is of interest because it has a significantly lower metallicity than DAs with otherwise similar atmospheric parameters. GD246 is the only white dwarf known that shows identifiable individual iron lines in the soft X-ray range. For the PG1159 star, a precise effective temperature determination is performed in order to confine the position of the blue edge of the GW Vir instability region in the HRD. Methods. The Chandra spectra are analyzed with chemically homogeneous as well as stratified NLTE model atmospheres that assume equilibrium between gravitational settling and radiative acceleration of chemical elements. Archival EUV and UV spectra obtained with EUVE, FUSE, and HST are utilized to support the analysis. Results. No metals could be identified in LB1919. All observations are compatible with a pure hydrogen atmosphere. This is in stark contrast to the vast majority of hot DA white dwarfs that exhibit light and heavy metals and to the stratified models that predict significant metal abundances in the atmosphere. For GD246 we find that neither stratified nor homogeneous models can fit the Chandra spectrum. The Chandra spectrum of PG1520+525 constrains the effective temperature to T(sub eff) = 150 000 +/- 10 000 K. Therefore, this nonpulsating star together with the pulsating prototype of the GWVir class (PG1159-035) defines the location of the blue edge of the GWVir instability region. The result is in accordance with predictions from nonadiabatic stellar pulsation models. Such models are therefore reliable tools to investigate the interior structure of GW Vir variables. Conclusions. Our soft X-ray study reveals that the understanding of metal abundances in hot DA white dwarf atmospheres is still incomplete. On the other hand, model atmospheres of hydrogen-deficient PG1159-type stars are reliable and reproduce well the observed spectra from soft X-ray to optical wavelengths.

Adamczak, J.↗

Validity of the Temperature Reconstruction from Water Isotopes in Ice Cores

Well-documented present-day distributions of stable water isotopes (HDO and others) show the existence, in middle and high latitudes, of a linear relationship between the mean annual isotope content of precipitation (SD and 51"0) and the mean annual temperature at the precipitation site. Paleoclimatologists have used this relationship, which is particularly well obeyed over Greenland and Antarctica, to infer paleotemperatures from ice core data. There is, however, growing evidence that spatial and temporal isotope/ surface temperature slopes differ, thus complicating the use of stable water isotopes as paleothermometers. In this paper we review empirical estimates of temporal slopes in polar regions and relevant information that can be inferred from isotope models: simple, Rayleigh-type distillation models and (particularly over Greenland) general circulation models (GCMS) fitted with isotope tracer diagnostics. Empirical estimates of temporal slopes appear consistently lower than present-day spatial slopes and are dependent on the timescale considered. This difference is most probably due to changes in the evaporative origins of moisture, changes in the seasonality of the precipitation, changes in the strength of the inversion layer, or some combination of these changes. Isotope models have not yet been used to evaluate the relative influences of these different factors. The apparent disagreement in the temporal and spatial slopes clearly makes calibrating the isotope paleothermometer difficult. Nevertheless, the use of a (calibrated) isotope paleothermometer appears justified; empirical estimates and most (though not all) GCM results support the practice of interpreting ice core isotope records in terms of local temperature changes.

Jouzel, J.↗

Modeling light and temperature effects on leaf emergence in wheat and barley

Phenological development affects canopy structure, radiation interception, and dry matter production; most crop simulation models therefore incorporate leaf emergence rate as a basic parameter. A recent study examined leaf emergence rate as a function of temperature and daylength among wheat (Triticum aestivum L.) and barley (Hordeum vulgare L.) cultivars. Leaf emergence rate and phyllochron were modeled as functions of temperature alone, daylength alone, and the interaction between temperature and daylength. The resulting equations contained an unwieldy number of constants. Here we simplify by reducing the constants by > 70%, and show leaf emergence rate as a single response surface with temperature and daylength. In addition, we incorporate the effect of photosynthetic photon flux into the model. Generic fits for wheat and barley show cultivar differences less than +/- 5% for wheat and less than +/- 10% for barley. Barley is more sensitive to daylength changes than wheat for common environmental values of daylength, which may be related to the difference in sensitivity to daylength between spring and winter cultivars. Differences in leaf emergence rate between cultivars can be incorporated into the model by means of a single, nondimensional factor for each cultivar.

NASA Discipline Number 61-10↗

Characterization of shale using Helium and Argon at high pressures

In order to estimate the shale gas in place and the eventual recovery during shale gas operations, one of the key requirements is to accurately characterize the shale’s petrophysical and transport properties such as porosity, permeability, diffusivity, and storage capacity. Despite the many efforts reported in the technical literature aiming to provide an improved understanding of the complex pore structures and the associated fluid flow in gas shales, complete characterization of organic-rich shale samples still poses a big challenge. Here, we have characterized mass transfer and sorption in shale at different length scales using Helium (He) and Argon (Ar) as probe gases. Thermogravimetric analysis (TGA) with a shale cube of ~1 cm 3 in size and gas expansion experiments with a full-diameter core (3.5” in diameter) were used to measure sorption kinetics/isotherms and mass transfer, respectively. Both samples are from the same depth/location in the Marcellus shale formation. The TGA steady-state technique was initially used to generate excess sorption isotherms for Ar, while dynamic TGA experiments were used to study its sorption kinetics. The TGA experiments demonstrate that Ar, which has a similar sorption potential as Methane, but is generally assumed to be inert, adsorbs onto the surfaces of the mesoporous and microporous regions of the shale samples according to a Langmuir-type behavior. Helium expansion experiments, on the full-diameter core, were used to measure the overall porosity, on the basis that He is a non-sorbing and inert gas as compared to Ar. The He expansion experiments, furthermore, allow us to delineate the mass transfer of gas across the inherent hierarchy of pore sizes, including macropores (macro- and microcracks), mesopores and micropores. Similar expansion experiments were also performed with Ar to study the combined impact of mass transfer and sorption. A triple-porosity model (TPM) was utilized to interpret the He expansion experiments with the shale core and to extract (estimate) relevant transport parameters. We report and compare here the diffusivities and permeabilities of the whole core for both He and Ar, as calculated from the modeling and fitting of the experimental data. On the premise that the shale cube is representative of the matrix region of the core, the Ar sorption kinetics from the cube experiments were subsequently combined with the transport parameters extracted from the He experiments to predict the behavior of the Ar expansion test with the full-diameter core. An excellent agreement is observed between the model predictions and the experimental data. The experimental observations and their interpretation indicate that one must be cautious when using Ar to estimate the true porosity and permeability of shales. In addition, we demonstrate that He and Ar probe gases, when used in tandem, can be employed effectively as a tool to characterize shales in terms of mass transfer and sorption dynamics across scales.

04 OIL SHALES AND TAR SANDS↗

Enhanced indirect optical absorption in AlAs and GaP.

Measurement and analysis of indirect interband optical absorption in GaP. Deviations from a simple indirect absorption law are shown to be consistent with those given by a model which takes into account both the variation of the energy denominators with photon energy and indirect absorption to higher-energy conduction-band valleys. The analysis of the GaP data indicates that such an analysis on data of limited absorption-constant range can provide reasonably accurate values of the direct band gap. The value of the direct band gap best fitting this model for AlAs at 6 K is 3.13 eV. In GaP and AlAs there is evidence of a higher set of extrema giving rise to an additional indirect absorption component 0.34 and 0.20 eV, respectively, above that due to the lowest-energy conduction-band minima. The strength of the absorption due to these higher-energy valleys indicates that they contain states of X sub 3 symmetry.

Dumke, W. P.↗

Permissible Exposure Level for Lunar Dusts: Gaps are Closing

Space faring nations plan to return human explorers to the moon within the next decade. Experience during the Apollo flights suggests that lunar dust will invariably get into the habitat where the finest portion (less than 5 micrometers) could be inhaled by the crew before it is cleared from the atmosphere. NASA is developing a database from which a 6-month, episodic exposure standard for lunar dust can be set. Three kinds of moon dust were prepared from a parent sample of Apollo 14 regolith #14003,96. Our goal was to prepare each type of dust sample with a mean diameter less than 2 m, which is suitable for instillation into the lungs of rats. The three samples were prepared as follows: separation from the parent sample using a fluidized bed, grinding using a jet mill grinder, or grinding with a ball-mill grinder. Grinding simulated restoration of surface activation of dust expected to occur at the surface of the moon on native lunar dust. We used two grinding methods because they seemed to produce different modes of activation. The effects of grinding were preserved by maintaining the dust in ultra-pure nitrogen until immediately before it was placed in suspension for administration to rats. The dust was suspended in physiological saline with 10% Survanta, a lung surfactant. Rats were given intratrachael instillations of the dust suspension at three doses. In addition to the three moon dusts (A, C and E), we instilled the same amount of a negative control (TiO2, B) and a highly-toxic, positive control (quartz, D). These additional mineral dusts were selected because they have well-established and very different permissible exposure levels (PELs). Our goal was to determine where lunar dusts fit between these extremes, and then estimate a PEL for each lunar dust. We evaluated many indices of toxicity to the lung. The figure shows the changes in lactate dehydrogenase (LDH), a marker of cell death, for the five dusts. Benchmark dose software (Version 2.1.2) from the Environmental Protection Agency was used to estimate the 10% effect levels (BMD(sub 10)) using five models. The best-fitting model was used to estimate the optimal BMD(sub 10) (table)

James, John T.↗

Relativistically Skewed Iron Emission and Disk Reflection in Galactic Microquasar XTE J1748-288

We report evidence for an Fe K-alpha fluorescence line feature in the Very High, High, and Low state X-ray spectra of the galactic microquasar XTE JI748-288 during its June 1998 outburst. Spectral analyses were made on observations spread across the outburst, gathered with the Rossi X-ray Timing Explorer. Gaussian line. disk emission line, relativistic disk emission line, and disk reflection models are fit to the data. In the Very High State, the line profile is strongly redshifted and consistent with emission from the innermost radius of a maximally rotating Kerr black hole, 1.235 R(sub g). The line profile is less redshifted in the High State, but increasingly prominent. In the Low State, the line profile is very strong and centered af approx. 6.7 keV; disk line emission models constrain the inner edge of the disk to fluctuate between approx.20 and approx.59 R(sub g). We trace the disk reflection fraction across the full outburst of this source, and find well-constrained fractions below those observed in AGN in the Very High and High States, but consistent with other galactic sources in the Low State. We discuss the possible implications for black hole X-ray binary system dynamics and accretion flow geometry.

Miller, J. M.↗

Efficient, Predictive Tomography of Multi-Qubit Quantum Processors

After decades of R&D, quantum computers comprising more than 2 qubits are appearing. If this progress is to continue, the research community requires a capability for precise characterization (“tomography”) of these enlarged devices, which will enable benchmarking, improvement, and finally certification as mission-ready. As world leaders in characterization -- our gate set tomography (GST) method is the current state of the art – the project team is keenly aware that every existing protocol is either (1) catastrophically inefficient for more than 2 qubits, or (2) not rich enough to predict device behavior. GST scales poorly, while the popular randomized benchmarking technique only measures a single aggregated error probability. This project explored a new insight: that the combinatorial explosion plaguing standard GST could be avoided by using an ansatz of few-qubit interactions to build a complete, efficient model for multi-qubit errors. We developed this approach, prototyped it, and tested it on a cutting-edge quantum processor developed by Rigetti Quantum Computing (RQC), a US-based startup. We implemented our new models within Sandia’s PyGSTi open-source code, and tested them experimentally on the RQC device by probing crosstalk. We found two major results: first, our schema worked and is viable for further development; second, while the Rigetti device is indeed a “real” 8-qubit quantum processor, its behavior fluctuated significantly over time while we were experimenting with it and this drift made it difficult to fit our models of crosstalk to the data.

97 MATHEMATICS AND COMPUTING↗

A canopy reflectance model based on an analytical solution to the multiple scattering equation

An approximation to the radiative transfer equation for solar radiation in relatively full, homogeneous plant canopies is presented and solved analytically for solar zenith angles less than 60 deg. The model predicts reflectance at any depth in the canopy and in any direction and may be inverted with bidirectional reflectance measurements. The model was fit to data at two sun angles and two wavebands (visible and NIR) to within the assumed errors on the reflectance data. The calculated albedos are insensitive to achievable measurement errors. Some of the parameter values themselves found by the inversion agree reasonably well with independent measurements, but the uncertainties introduced by the data noise are rather large. However, the agreement is good enough to demonstrate that the model is physically realistic.

Camillo, Peter↗

Effect of Pressurization on Solid Oxide Cell Oxygen Electrodes: the Role of PrO x Nanoparticle Infiltration

Pressurization of solid oxide cells improves performance by reducing electrode polarization resistance (R P ) and facilitates system integration with balance of plant components such as pressurized storage tanks. However, there are few reports on pressurization effects for electrodes designed for low-temperature operation and utilizing infiltrated catalysts. Here we report an electrochemical impedance spectroscopy study of high performing oxygen electrode materials, SrTi 0.3 Fe 0.63 Co 0.07 O 3-∂ (STFC) and PrO x infiltrated STFC, for oxygen partial pressures ($p$ $o$ 2 ) from 0.1 to 8 atm and temperatures from 550 to 650 °C. decreases more with pressurization for STFC:PrO x , fitting well to with an exponent n~0.3, compared to n~0.25 for STFC. The combination of PrO x infiltration and pressurization yields a substantial R P decrease, e.g., at 600 °C by ~7 times from 0.36 Ω cm 2 at $p$ $o$ 2 = 0.2 atm for STFC to 0.055 Ω cm 2 at $p$ $o$ 2 = 4 atm for STFC:PrO x . Here, a transmission-line-based circuit model impedance fit reveals that the significant oxygen surface reaction (R surf ) resistance contribution decreases substantially with PrO x infiltration; and its $p$ $o$ 2 dependence become more pronounced, with n increasing from ~0.25 to ~0.5. R surf for STFC:PrO x decreases so much at elevated $p$ $o$ 2 that the electrode/electrolyte interface resistance dominates.

25 ENERGY STORAGE↗

Gamma-ray observations of Ophiuchus with EGRET: The diffuse emission and point sources

Observations of the Ophiuchus region made with the Energetic Gamma-Ray Experiment Telescope (EGRET) on the Compton Gamma Ray Observatory (CGRO) during the first 2 1/2 years of operation show the diffuse emission from the interstellar gas in Ophiuchus as well as variable emission from two point sources. The gamma-ray emission is modeled in terms of cosmic-ray interactions with atomic and molecular hydrogen in Ophiuchus and with low-energy photons along the line of sight. The model also includes the flux from the two point sources and an isotropic diffuse contribution. The cosmic-ray density is assumed to be uniform. The derived ratio of molecular hydrogen column density to integrated CO intensity is (1.1 +/- 0.2) x 10(exp 20) H-mols/sq cm (K km/s)(exp -1). At the sensitivity and resolution of the gamma-ray data, no variation of this ratio over the modeled region is discernible, nor are any regions of enhanced cosmic-ray density apparent. The model was fitted to seven narrow energy bands to obtain the energy depedence of the gamma-ray production function and the spectra of the point sources. The derived production function is in good agreement with theoretical calculations and the local cosmic-ray electron and proton spectra. The positions of the point sources were determined from maximum likelihood analysis of the gamma-ray emission observed in excess of the diffuse model. We identify one point source with the quasar PKS 1622-253, which has an average flux, E greater than 100 MeV, of (2.5 +/- 0.5) x 10(exp -7) photons/sq cm/s and photon spectral index -1.9 +/- 0.3. The other source, denoted GRO J1631-27, has not yet been identified at other wavelengths. Its average flux, E greater than 100 MeV, is (1.1 +/- 0.4) x 10(exp -7) photons/sq cm/s; however, its spectral index is poorly determined. The spectral index and intensity of the isotropic contribution to the model agree well with the extragalactic diffuse emission derived from the SAS 2 data.

Hunter, S. D.↗

Design Considerations for Heavily-Doped Cryogenic Schottky Diode Varactor Multipliers

Diode modeling for Schottky varactor frequency multipliers above 500 GHz is presented with special emphasis placed on simple models and fitted equations for rapid circuit design. Temperature- and doping-dependent mobility, resistivity, and avalanche current multiplication and breakdown are presented. Next is a discussion of static junction current, including the effects of tunneling as well as thermionic emission. These results have been compared to detailed measurements made down to 80 K on diodes fabricated at JPL, followed by a discussion of the effect on multiplier efficiency. Finally, a simple model of current saturation in the undepleted active layer suitable for inclusion in harmonic balance simulators is derived.

Schlecht, E.↗

Relating ground truth collection to model sensitivity

The importance of collecting high quality ground truth before a SAR mission over a forested area is two fold. First, the ground truth is used in the analysis and interpretation of the measured backscattering properties; second, it helps to justify the use of a scattering model to fit the measurements. Unfortunately, ground truth is often collected based on visual assessment of what is perceived to be important without regard to the mission itself. Sites are selected based on brief surveys of large areas, and the ground truth is collected by a process of selecting and grouping different scatterers. After the fact, it may turn out that some of the relevant parameters are missing. A three-layer canopy model based on the radiative transfer equations is used to determine, before hand, the relevant parameters to be collected. Detailed analysis of the contribution to scattering and attenuation of various forest components is carried out. The goal is to identify the forest parameters which most influence the backscattering as a function of frequency (P-, L-, and C-bands) and incident angle. The influence on backscattering and attenuation of branch diameters, lengths, angular distribution, and permittivity; trunk diameters, lengths, and permittivity; and needle sizes, their angular distribution, and permittivity are studied in order to maximize the efficiency of the ground truth collection efforts. Preliminary results indicate that while a scatterer may not contribute to the total backscattering, its contribution to attenuation may be significant depending on the frequency.

Amar, Faouzi↗