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 145 records · Page 8

Recent observations of cosmic ray antiprotons and a critical assessment of the theories of their origin

The models proposed to explain the observed spectrum of cosmic ray antiprotons are reexamined in light of recent results from balloon-borne experiments. It is found that the prediction of modified closed galaxy model fits the observed data very well. Models in which secondary antiprotons are produced in the sources, could be made consistent with the data provided the secondary particles do not suffer considerable adiabatic deceleration. It has been shown that there cannot be any significant contribution to the observed antiprotons, from the evaporation of mini black holes or from the annihilation of dark matter like photinos. The role of extragalactic cosmic rays has been examined critically in the context of the recent data, and they are not the source of cosmic ray antiprotons. However, determination of the energy spectrum of antiprotons at least up to a few tens of GeV would be valuable to provide information on the possible existence of supersymmetric particles and on the modulation of extragalactic cosmic rays while entering the Galaxy.

Stephens, S. A.↗

Stellar populations in edge-on galaxies from deep CCD surface photometry, 1: NGC 5907

We present extremely deep charge coupled device (CCD) surface photometry of the edge-on Sc galaxy NGC 5907. Our data reach reliably to a surface brightness of R=27 mag/sq arcsec, some two magnitudes fainter than any previous work. We obtained this improvement using a 2048X2048 CCD with a wide (approximately 24 min) field, which made it possible to sky subtract directly from the galaxy frame, and by taking many dark sky flatfields. Our analysis of these data, using a full 2D model fitting procedure with a detailed error model, confirm the thin disk parameters of van der Kruit & Searle (1981). In particular, we confirm that the galaxy's disk has a radial cutoff and a constant scale height with radius. We find evidence for a stellar warp in this system, which has the same orientation as the H I warp. Our deep data also confirm that NGC 5907 has no thick disk. This suggests that theories of thick disk formation from star formation in the early stages of disk collapse, or by secular heating mechanisms, are unlikely to be correct, because they would predict that every galaxy would have a thick disk. Thick disk formation from the accretion of satellite galaxies is more likely.

Morrison, Heather L.↗

A Hubble Space Telescope study of the underlying white dwarf in the dwarf nova VW Hydri during quiescence

We have analyzed a far-ultraviolet spectrum of the dwarf nova VW Hyi obtained during quiescence with the Hubble Space Telescope Faint Object Spectrograph (HST FOS). The observation occurred 10 days after the return to optical quiescence from a superoutburst of VW Hyi. The spectrum reveals a very strong Stark-broadened Ly-alpha absorption with narrow geocoronal emission, and a very rich metallic absorption-line spectrum dominated by strong resonance absorption features of singly and doubly ionized silicon and carbon, the first solid identification of metallic absorption features arising in the accreted atmosphere of the white dwarf. We confirm the reported low-resolution IUE detection of the underlying white dwarf photosphere by Mateo & Szkody. A synthetic spectral analysis with hot, high-gravity LTE model atmospheres yields a best-fit model with the following parameters: T(sub eff) = 22,000 +/- 1000 K, log g = 8.0 +/- 0.3, with chemical abundances for oxygen of 0.3 times solar, for nitrogen of 5 times solar, and for all other heavy elements of 0.15 times solar. Based upon our absorption-line measurements of the observations at different orbital phases, we find no conclusive evidence of equivalent width variations with orbital phase. In the absence of any significant reduction of the white dwarf's core mass by past nova explosions, its lower limit cooling age is approximately 50 million years.

Sion, Edward M.↗

Nonlinear Aerodynamic Modeling From Flight Data Using Advanced Piloted Maneuvers and Fuzzy Logic

Results of the Aeronautics Research Mission Directorate Seedling Project Phase I research project entitled "Nonlinear Aerodynamics Modeling using Fuzzy Logic" are presented. Efficient and rapid flight test capabilities were developed for estimating highly nonlinear models of airplane aerodynamics over a large flight envelope. Results showed that the flight maneuvers developed, used in conjunction with the fuzzy-logic system identification algorithms, produced very good model fits of the data, with no model structure inputs required, for flight conditions ranging from cruise to departure and spin conditions.

Brandon, Jay M.↗

Dust in the small Magellanic Cloud. 2: Dust models from interstellar polarization and extinction data

We study the dust in the Small Magellanic Cloud using our polarization and extinction data (Paper 1) and existing dust models. The data suggest that the monotonic SMC extinction curve is related to values of lambda(sub max), the wavelength of maximum polarization, which are on the average smaller than the mean for the Galaxy. On the other hand, AZV 456, a star with an extinction similar to that for the Galaxy, shows a value of lambda(sub max) similar to the mean for the Galaxy. We discuss simultaneous dust model fits to extinction and polarization. Fits to the wavelength dependent polarization data are possible for stars with small lambda(sub max). In general, they imply dust size distributions which are narrower and have smaller mean sizes compared to typical size distributions for the Galaxy. However, stars with lambda(sub max) close to the Galactic norm, which also have a narrower polarization curve, cannot be fit adequately. This holds true for all of the dust models considered. The best fits to the extinction curves are obtained with a power law size distribution by assuming that the cylindrical and spherical silicate grains have a volume distribution which is continuous from the smaller spheres to the larger cylinders. The size distribution for the cylinders is taken from the fit to the polarization. The 'typical', monotonic SMC extinction curve can be fit well with graphite and silicate grains if a small fraction of the SMC carbon is locked up in the grain. However, amorphous carbon and silicate grains also fit the data well. AZV456, which has an extinction curve similar to that for the Galaxy, has a UV bump which is too blue to be fit by spherical graphite grains.

Rodrigues, C. V.↗

Local models for scatter estimation and descattering in polyenergetic X-ray tomography

We propose a new modeling approach for scatter estimation and descattering in polyenergetic X-ray computed tomography (CT) based on fitting models to local neighborhoods of a training set. X-ray CT is widely used in medical and industrial applications. X-ray scatter, if not accounted for during reconstruction, creates a loss of contrast in CT reconstructions and introduces severe artifacts including cupping, shading, and streaks. Even when these qualitative artifacts are not apparent, scatter can pose a major obstacle in obtaining quantitatively accurate reconstructions. Our approach to estimating scatter is, first, to generate a training set of 2D radiographs with and without scatter using particle transport simulation software. To estimate scatter for a new radiograph, we adaptively fit a scatter model to a small subset of the training data containing the radiographs most similar to it. We compared local and global (fit on full data sets) versions of several X-ray scatter models, including two from the recent literature, as well as a recent deep learning-based scatter model, in the context of descattering and quantitative density reconstruction of simulated, spherically symmetrical, single-material objects comprising shells of various densities. Our results show that, when applied locally, even simple models provide state-of-the-art descattering, reducing the error in density reconstruction due to scatter by more than half.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Vertical Structure and Optical Properties of Titans Aerosols from Radiance Measurements Made Inside and Outside the Atmosphere

Prompted by the detection of stratospheric cloud layers by Cassini's Composite Infrared Spectrometer (CIRS; see Anderson, C.M., Samuelson, R.E. [2011]. Icarus 212, 762-778), we have re-examined the observations made by the Descent Imager/Spectral Radiometer (DISR) in the atmosphere of Titan together with two constraints from measurements made outside the atmosphere. No evidence of thin layers (<1 km) in the DISR image data sets is seen beyond the three previously reported layers at 21 km, 11 km, and 7 km by Karkoschka and Tomasko (Karkoschka, E., Tomasko, M.G. [2009]. Icarus 199, 442-448). On the other hand, there is evidence of a thicker layer centered at about 55 km. A rise in radiance gradients in the Downward-Looking Visible Spectrometer (DLVS) data below 55 km indicates an increase in the volume extinction coefficient near this altitude. To fit the geometric albedo measured from outside the atmosphere the decrease in the single scattering albedo of Titan's aerosols at high altitudes, noted in earlier studies of DISR data, must continue to much higher altitudes. The altitude of Titan's limb as a function of wavelength requires that the scale height of the aerosols decrease with altitude from the 65 km value seen in the DISR observations below 140 km to the 45 km value at higher altitudes. We compared the variation of radiance with nadir angle observed in the DISR images to improve our aerosol model. Our new aerosol model fits the altitude and wavelength variations of the observations at small and intermediate nadir angles but not for large nadir angles, indicating an effect that is not reproduced by our radiative transfer model. The volume extinction profiles are modeled by continuous functions except near the enhancement level near 55 km altitude. The wavelength dependence of the extinction optical depth is similar to earlier results at wavelengths from 500 to 700 nm, but is smaller at shorter wavelengths and larger toward longer wavelengths. A Hapke-like model is used for the ground reflectivity, and the variation of the Hapke single scattering albedo with wavelength is given. Fits to the visible spectrometers looking upward and downward are achieved except in the methane bands longward of 720 nm. This is possibly due to uncertainties in extrapolation of laboratory measurements from 1 km-am paths to much longer paths at lower pressures. It could also be due to changes in the single scattering phase functions at low altitudes, which strongly affect the path length through methane that the photons travel. We demonstrate the effects on the model fits by varying each model parameter individually in order to illustrate the sensitivity of our determination of each model parameter.

Titan↗

INTEGRAL Observations of the Enigmatic Be Stars (gamma) Cassiopeiae and HD 110432

We present the results of a hard X-ray study of the Be stars gamma Cassiopeiae and HD 110432 based on observations made with the INTEGRAL observatory. These stars are known to be moderately strong, X-ray sources (L(sub x) approx. equal to = 10(sup 32)-10(sup 33) erg per second). These values are at the extreme high end of the known luminosity distribution for active coronal systems, but several orders of magnitude below typical X-ray binaries. The hard X-ray spectra for these systems are quite similar. They can be well fitted by either optically thin thermal plasma models with kT = 12.5 - 14 keV or a cutoff powerlaw + gaussian line model with photon indices in the 1.3 - 1.5 range and a line energy of 6.7 keV. The 20-50 keV light curves show no evidence for flaring and no significant evidence for periodic variability. It has been proposed that the X-ray emission is due to either accretion onto a white dwarf companion or magnetic activity near the surface of the Be star. We discuss in detail the pros and cons of each scenario towards explaining our spectral and temporal results. Given that both thermal and nonthermal models fit the data equally well, we cannot use the spectra to delineate between these two scenarios. Recent observations indicate that gamma Cas has a approx. 1 solar mass companion in a 203.59 day orbit. This is consistent with the white dwarf - Be star binary model but the lack of periodic modulation of the flux on this timescale calls this conclusion into question. On the other hand the lack of flaring activity may rule against the magnetic activity model. We discuss advances in observations and theory that need to be made to resolve the origin of these systems.

Sturner, S. J.↗

Response of CR-39 nuclear track detectors to protons with non-normal incidence

This paper presents data from experiments with protons at non-normal incidence to CR-39 nuclear track detectors, analyzing the properties of detection efficiency, proton track diameter, track contrast, and track eccentricity. Understanding the CR-39 response to protons incident at an angle is important for designing charged particle detectors for inertial confinement fusion (ICF) applications. This study considers protons with incident energies less than 3 MeV. In this regime, an incident angle of 10° has no effect on CR-39 detection efficiency, and >85% detection efficiency is preserved up through 25° in the range of 1.0 MeV–2.1 MeV. For ICF applications, incident angles above 30° are deemed impractical for detector design due to significant drops in proton detection at all energies. We observe significant reductions in detection efficiency compared to theoretical predictions, particularly at low energies where proton tracks are etched away. The proton track diameter measured by the scan system is observed to decrease with higher incident angles. The track diameters are analyzed with two fitting models, and it is shown that the diameter–energy relation can be fit with the existing models at angles up to 30° . The optical contrast of the tracks tends to increase with the angle, meaning that the tracks are fainter, and a larger increase is observed for higher energies. Eccentricity, a measure of how elongated proton tracks are, increases with the incident angle and drops after the critical angle. The lowest energy tracks remain nearly circular even at higher angles.

47 OTHER INSTRUMENTATION↗

Response of CR-39 nuclear track detectors to protons with non-normal incidence

This paper presents data from experiments with protons at non-normal incidence to CR-39 nuclear track detectors, analyzing the properties of detection efficiency, proton track diameter, track contrast, and track eccentricity. Understanding the CR-39 response to protons incident at an angle is important for designing charged particle detectors for inertial confinement fusion (ICF) applications. This study considers protons with incident energies less than 3 MeV. In this regime, an incident angle of 10° has no effect on CR-39 detection efficiency, and >85% detection efficiency is preserved up through 25° in the range of 1.0 MeV–2.1 MeV. For ICF applications, incident angles above 30° are deemed impractical for detector design due to significant drops in proton detection at all energies. We observe significant reductions in detection efficiency compared to theoretical predictions, particularly at low energies where proton tracks are etched away. The proton track diameter measured by the scan system is observed to decrease with higher incident angles. The track diameters are analyzed with two fitting models, and it is shown that the diameter–energy relation can be fit with the existing models at angles up to 30°. The optical contrast of the tracks tends to increase with the angle, meaning that the tracks are fainter, and a larger increase is observed for higher energies. Eccentricity, a measure of how elongated proton tracks are, increases with the incident angle and drops after the critical angle. The lowest energy tracks remain nearly circular even at higher angles.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Missing components in ΛCDM from DESI Y1 baryonic acoustic oscillation measurements: Insights from redshift remapping

We explore transformations of the Friedman-Lemaître-Robertson-Walker (FLRW) metric and cosmological parameters that align with observational data while aiming to gain insights into potential extensions of standard cosmological models. We modified the FLRW metric by introducing a scaling factor, e 2Θ(a) –the cosmological scaling function (CSF), which alters the standard relationship between cosmological redshift and the cosmic scale factor without affecting angular measurements or cosmic microwave background (CMB) anisotropies. Using data from DESI Year 1, Pantheon+ supernovae, and the Planck CMB temperature power spectrum, we constrained both the CSF and cosmological parameters through a Markov chain Monte Carlo approach. Our results indicate that the CSF model fits observational data with a lower Hubble constant (although it is compatible with the value given by Planck 2018 within 1σ) and is predominantly dark matter dominated. Additionally, the CSF model produces temperature and lensing power spectra similar to those predicted by the standard model, though with lower values in the CSF model at large scales. We also checked that when fitting a CSF model without dark energy to the data, we obtain a more negative conformal function. This suggests that the CSF model may offer hints about missing elements and opens up a new avenue for exploring physical interpretations of cosmic acceleration.

79 ASTRONOMY AND ASTROPHYSICS↗

Similarity Metric for Data Optimization and Efficient Training of Reactive Machine Learning Force Fields for Hydrocarbon Radiolysis

Radiolysis is a common approach to sterilize polymers, chemically modify them for upcycling, and accelerate their decomposition for recycling purposes. Reactive molecular dynamics (MD) simulations provide a powerful tool to generate atomic-level trajectories of the reactive processes and quantify radiolytic chemical degradation pathways. For this, machine learning (ML) surrogate models for reactive force fields with quantum mechanical accuracy are now widely used, which require ML training data sets that can provide information on atomic environments for target chemical systems. However, radiolysis chemistry can be highly complex and diverse, which poses significant challenges for generating training data to parametrize ML models. In this regard, we developed a method for optimizing the training data set using a cosine similarity metric to help guide training set selection for radiolysis of polyethylene, a model hydrocarbon polymer, as well as to enhance the transferability of our reactive ML force field (MLFF) to a variety of molecular and polymeric systems. Our approach performs atom-by-atom comparisons between local atomic environments to pinpoint important data points associated with rare and localized events, such as radiolysis damage within structures. We apply this approach to train the Chebyshev Interaction Model for Efficient Simulation (ChIMES) MLFF model, which expresses the atomic interaction potentials in terms of linear combinations of many-body Chebyshev polynomials. We first show that our method can reduce our training set size by ∼70% while improving overall accuracy compared to more standard MD model fitting approaches. We then validate our optimum model against diverse hydrocarbon simulation data, including simple alkanes and systems with unsaturated carbon bonds, over a wide range of thermodynamic conditions. Finally, we use our ChIMES model to perform MD simulations of radiolytic damage with large-scale systems that help avoid system size effects. Overall, our approach yields an MD force field that retains most of the accuracy of the underlying quantum method while yielding many orders of improvement in computational efficiency. In conclusion, our efforts will have impact on future hydrocarbon polymer radiolysis studies, where the chemical details of the polymer–radiation interactions can have a strong effect on the resulting products observed in experiments.

Hydrocarbons↗

Verification Testing of OLI Systems Mixed Solvent Electrolyte Model for the Na-K-Mg-Ca-H-Cl-SO 4 -OH-HCO 3 -CO 3 -CO 2 -H 2 ) System to High Ionic Strength at 25°C

This technical report summarizes model verification results and summary statistics for 41 evaporite mineral solubility cases evaluated by Savannah River National Laboratory using OLI Systems’ aqueous electrolyte thermodynamic modeling software. The 41 verification cases containing a total of 60 solubility curves comprise mineral solubility data from low to high ionic strength at 25°C for the eight-component system Na-K-Mg-Ca-H-Cl-SO 4 -OH-HCO 3 -CO 3 -CO 2 -H 2 O as reported by Harvie et al. (1984). Thermodynamic calculations were executed using OLI Systems’ Stream Analyzer computation module within the OLI Studio software platform (Ver. 11.0, Rev. 11.0.1.9). The Mixed Solvent Electrolyte (MSE) thermodynamic framework was chosen for this investigation because of its superiority in modeling high ionic-strength inorganic salt solutions and actinide redox chemistry and solubility, both of which are relevant to the geological repository conditions at the Waste Isolation Pilot Plant in Carlsbad, New Mexico. Mineral solubility data in various inorganic salt solutions were digitized and extracted from figures generated by Harvie et al. (1984). For each of the 60 solubility curves, a case-specific chemistry model and input file were generated in OLI Studio using OLI Stream Analyzer and the MSE (H 3 O + ion) public databank provided by OLI Systems. Model simulation results were exported to Microsoft Excel to calculate summary statistics and to generate graphs comparing the OLI model predictions to the solubility data. Summary statistics include residuals (model – data) and concordance (accuracy × precision, where precision is indicated by the Pearson correlation coefficient and accuracy accounts for bias and scale differential). Private databanks were not developed, and activity coefficient model regressions were not performed to improve OLI model fits to the data. Of the 41 model verification plots, 83% have a mean of the percent residuals less than or equal to 25%. Similarly, 75% display a concordance greater than or equal to 0.75. Only seven of the 41 verification plots fail to show good agreement between the model and data. Of these seven, three are relevant to the WIPP repository because they involve the Mg-OH-Cl-SO 4 -CO 3 aqueous system. The remaining four address salt solubilities at the pH extremes (strong acid and strong base). It should be noted that in two of the three Mg-OH-Cl-SO 4 -CO 3 system cases, the regressed Harvie et al. (1984) solubility curve also deviated from the data. Lack of agreement between the OLI model-predicted solubility curves and the data is attributable to one or more of the following: specific solid species are not included in the OLI MSE databank; there is significant variation among the different solubility datasets chosen by Harvie et al. (1984); the OLI MSE model’s thermodynamic parameters were determined using different solubility datasets; and the activity coefficient parameters for certain relevant ion-ion and ion-molecule pairs have not been optimized via data regression. Two recommendations for future work are to (1) evaluate solubility data for the Mg-OH-Cl-SO 4 -CO 3 system at high ionic strength and, if necessary, develop a private OLI MSE database that includes missing species and, where necessary, regressed standard state properties and interaction parameters; (2) perform similar verification testing of the OLI model for actinide solubility data.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Glass Property-Composition Models Update for use in Direct Feed High-Level Waste Flowsheet Development

A set of preliminary glass property models and constraints were developed and augmented by models from literature for use in design of direct-feed high-level waste (DFHLW) glasses for flowsheet evaluation, testing, and design of the Tank Waste Treatment and Immobilization Plant (WTP) high-level waste (HLW) Facility. These models and constraints are meant to be used as a place-holder while glass property-composition data gaps are filled and final plant operating models are developed. This report describes the motivation and intended use of the models, the compilation of data, model fitting and selection, methods to apply the models and constraints in glass design and offers example calculations demonstrating their intended use.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Glass Property-Composition Models Update for use in Direct Feed High-Level Waste Flowsheet Development: EWG2.6

A set of preliminary glass property models and constraints were developed and augmented by models from literature for use in design of direct-feed high-level waste (DFHLW) glasses for flowsheet evaluation, testing, and design of the Waste Treatment and Immobilization Plant (WTP) high-level waste (HLW) Facility. These models and constraints are meant to be used as a place-holder while glass property-composition data gaps are filled and final plant operating models are developed. This report describes the motivation and intended use of the models, the compilation of data, model fitting and selection, methods to apply the models and constraints in glass design and offers example calculations demonstrating their intended use.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Galilean satellites - Observations of mutual occultations and eclipses in 1973

Seven Galilean satellite mutual events, two occultations and two eclipses of Europa and three eclipses of Io, were observed at three wavelengths (0.35, 0.50, 0.91 micrometers) with a time resolution of 0.1 sec. Preliminary model fits to the light curves are presented. Model satellites with different albedo distributions (uniform disk, bright solar caps, bright quadrant) are used in generating model occultation and eclipse curves to demonstrate the sensitivity of observed light curves to the brightness distribution on the surface of the occulted or eclipsed satellite. At the present the observations yield no conclusive information of the limb-darkening of Io. The best data for Europa indicate that the satellite is limb-darkened at both 0.50 and 0.91 micrometers.

Wasserman, L. H.↗

The phase of the ten-hour modulation in the Jovian magnetosphere /Pioneers 10 and 11/

The paper describes the study of the phase of the 10-hour modulation of energetic electrons seen by Pioneers 10 and 11 in the Jovian magnetosphere. Attention is given to the peaks rather than the valleys of each cycle because the peaks are where physically interesting features occur, such as particle acceleration, current sheets, etc. To identify the peaks, it is required that the instantaneous intensity be higher than the 5-hour running average and the 5-hour running average be greater than the 10-hour running average. These criteria select an interval rather than a point and it is determined that this interval is an appropriate estimate of the experimental uncertainty. When the phases of the peaks are plotted together, they create patterns which are discussed in terms of disk-like, clock-like, and rotating anomaly models of the magnetosphere. Each model fits some of the data, but no model explains all of the data convincingly. It is concluded that there is still no understanding of the configuration of the outer Jovian magnetosphere.

Fillius, W.↗

Interacting binary galaxies. V - NGC 4782/4783 (3C 278): Unbound colliders, not a supermassive pair

CCD imaging data of the elliptical galaxies NGC 4782 and NGC 4783 have been analyzed along with published spectroscopic measurements in order to develop a viable interaction model for this system. The best-fit model suggests that the binary orbit was hyperbolic before the encounter, but that the energy loss during the collision has led to a loosely bound system. The total system mass in the region occupied by the luminous matter is shown to be 1.4 x 10 to the 12th solar masses. The results support the previous theory that the internal velocity dispersions of colliding galaxies tend to increase just past closest approach.

Borne, Kirk D.↗