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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Review of PBX 9502 Pop Plot data

Pop plot data (distance of run and time to detonation as a function of shock pressure) is a simple characterization of shock initiation of a high explosive. It is used, along with other data, to calibrate reactive burn models. For the insensitive TATB based explosive PBX 9502, there are data at 5 initial temperatures (77 K, -55 C, 23 C, 76 C, 130 C) and for several lots. Here we review the dependence of the Pop plot on temperature and the variation with lot. There is a significant increase in shock sensitivity with temperature.Since the Pop plot is important for modeling, fit parameters are interpolated to cover the temperature range –55 < T < 76 C. The scatter in the data and variation with lot is important for assessing the accuracy to which a model can be calibrated.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
NASA NTRS1975

The 1973 solar occultation of the Crab Nebula pulsar

The mean electron density of the solar corona was determined by measuring the dispersion of radiofrequency pulses from pulsar NP 0532 during the June 1973 solar occultation. Trends continued which were noticed in 1971 as solar activity declined. Model fitting results suggest that the corona continued to become even more concentrated toward the equator in 1973 than in 1971. The number density of electrons in most regions decreased. The best model of the distribution of corona electrons is suggested to be one with zero density at the poles. K-corona isophotes and contours of equal path-integrated density are presented for several models. Electron density versus date and position in the corona are tabulated. It is seen that there is no simple relationship between the onset of major solar activity and density or scattering enhancements.

Weisberg, J. M.↗
DOE OSTI2016

Investigation of SiPM physics parameters down to cryogenic temperatures and for a bio-medical application

Silicon PhotoMultiplier (SiPM) detector has become a suitable visible light/photon detector for many applications like high energy physics and neutrino experiments, fluorescence detection, bio-photonics and medical imaging. The first part of my thesis was oriented to the studies of SiPM physics parameters as a function of temperature. Particularly, recent KETEK devices (year 2015) with different technological characteristics like p/n and n/p junctions, with and without trench technology, and different widths of epitaxial layer were studied in the temperature range from 308.15 K (+35°C) down to 238.15 K (-35°C). In addition, the Hamamatsu devices from 2011 production run as well as new devices from 2015 year, with improved technological characteristics inducing a reduced noise, were investigated in a wider temperature range from 318.15 K (+45°C) down to 98.15 K (-175°C). For these purposes, I participated to the design, installation, commissioning and calibration of a cryogenic experimental setup dedicated to electrical, optical and temperature studies of SiPM devices. Also, I have developed an automatic analysis procedure able to handle in a short time an impressive quantity of experimental data (i.e. tens of Gb/device) and to give a precise and fast information on main SiPM parameters and their temperature dependence. I have also developed a physical modeldescribing the DC I-V curves of SiPM detectors at different temperatures. The proposed model fits well the shape of IV curve in a very large currents range from 10⁻¹² A up to 10⁻⁵ A over the full working range of various devices. Consequently, the IV model can be used as a simple and fast method for determination of SiPM parameters like breakdown voltage VBD, the shape of Geiger triggering probability PGeiger as a function of Vbias as well as the Vbias working range. The comparison of these parameters with those calculated from AC measurements and analyzed by the automatic procedure showed a good agreement. The second part of my thesis was oriented to the study of SiPM devices and their physical parameters required to build a prototype of betasensitive intracerebral probe. Such probe is dedicated to measure the local concentration of radiolabeled molecules on awake and freely moving animal and to study new animal models of human disorders (neurodegenerative diseases, tumor growth, and neuropsychiatric disorders). It is composed of small size, low-noise SiPM device coupled to a scintillating fiber and readout by a dedicated miniaturized low-power consumption electronics. Three SiPM devices have been chosen as the most adapted for our application: two small KETEK devices of 0.5×0.5 mm² size (with and without optical trenches, specially developed by KETEK to fulfill our requirements) and a standard Hamamatsu device of 1.3×1.3 mm² size, all devices having 50 × 50 μm² μcell size. For each SiPM the gain G, dark count rate DCR and beta sensitivity were measured as a function of Vbias and temperature. The obtained results showed that the small field of view and newly developed structure of the KETEK devices allow a large decrease of the dark count rate DCR. However, this small field of view also leads to a reduced light collection due to the thickness of the epoxy protection resin on top of the SiPM and the acceptance angle of the fiber. Since the beta sensitivity represents a tradeoff between photon detection efficiency PDE and dark count rate DCR, KETEK SiPMs exhibit similar performances in comparison with the Hamamatsu device. Preliminary results demonstrate that the beta sensitivity of KETEK devices can be significantly improved by using focusing lens between the scintillating fiber and the SiPM or by reducing the thickness of its epoxy protection resin.

Nagai, Andrii↗
DOE OSTI2018

An Indication of Anisotropy in Arrival Directions of Ultra-high-energy Cosmic Rays through Comparison to the Flux Pattern of Extragalactic Gamma-Ray Sources

A new analysis of the dataset from the Pierre Auger Observatory provides evidence for anisotropy in the arrival directions of ultra-high-energy cosmic rays on an intermediate angular scale, which is indicative of excess arrivals from strong, nearby sources. The data consist of 5514 events above 20 EeV with zenith angles up to 80 deg recorded before 2017 April 30. Sky models have been created for two distinct populations of extragalactic gamma-ray emitters: active galactic nuclei from the second catalog of hard Fermi-LAT sources (2FHL) and starburst galaxies from a sample that was examined with Fermi-LAT. Flux-limited samples, which include all types of galaxies from the Swift-BAT and 2MASS surveys, have been investigated for comparison. The sky model of cosmic-ray density constructed using each catalog has two free parameters, the fraction of events correlating with astrophysical objects and an angular scale characterizing the clustering of cosmic rays around extragalactic sources. A maximum-likelihood ratio test is used to evaluate the best values of these parameters and to quantify the strength of each model by contrast with isotropy. It is found that the starburst model fits the data better than the hypothesis of isotropy with a statistical significance of 4.0 sigma, the highest value of the test statistic being for energies above 39 EeV. The three alternative models are favored against isotropy with 2.7-3.2 sigma significance. The origin of the indicated deviation from isotropy is examined and prospects for more sensitive future studies are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1972

Estimation of automobile-driver describing function from highway tests using the double steering wheel

The automobile-driver describing function for lateral position control was estimated for three subjects from frequency response analysis of straight road test results. The measurement procedure employed an instrumented full size sedan with known steering response characteristics, and equipped with a lateral lane position measuring device based on video detection of white stripe lane markings. Forcing functions were inserted through a servo driven double steering wheel coupling the driver to the steering system proper. Random appearing, Gaussian, and transient time functions were used. The quasi-linear models fitted to the random appearing input frequency response characterized the driver as compensating for lateral position error in a proportional, derivative, and integral manner. Similar parameters were fitted to the Gabor transformed frequency response of the driver to transient functions. A fourth term corresponding to response to lateral acceleration was determined by matching the time response histories of the model to the experimental results. The time histories show evidence of pulse-like nonlinear behavior during extended response to step transients which appear as high frequency remnant power.

Delp, P.↗
NASA NTRS1981

The Voyager 1 EUV spectrum of the Io plasma torus

The Voyager 1 EUV spectrum of the hot Io plasma torus obtained before encounter in 1979 at a range of 4,800,000 km is analyzed in terms of a collisional model with improved collision strengths. Model calculations are based entirely on steady-state equilibrium between electron excitation-deactivation and radiative decay in the subspecies structure, and plasma temperature and species number densities are derived by fitting model spectra to the observations. Lines of S II, S III, S IV, O II, and O III are found to be dominant in the plasma torus, at torus half width densities of 44, 160, 220, 50, and 340/cu cm, respectively, and lines of S V and K III may also be present. An effective electron temperature of 80,000 K is estimated for the central dense region near 6 Jupiter radii. Model calculations are in substantial agreement with other spectroscopic observations of the hot torus; however, the observations exhibit strong morphological differences with recent ground-based observations of S II emission, and calculated ion densities differ from in situ Voyager 1 plasma science measurements.

Shemansky, D. E.↗
NASA NTRS2004

Progressive Aerodynamic Model Identification From Dynamic Water Tunnel Test of the F-16XL Aircraft

Development of a general aerodynamic model that is adequate for predicting the forces and moments in the nonlinear and unsteady portions of the flight envelope has not been accomplished to a satisfactory degree. Predicting aerodynamic response during arbitrary motion of an aircraft over the complete flight envelope requires further development of the mathematical model and the associated methods for ground-based testing in order to allow identification of the model. In this study, a general nonlinear unsteady aerodynamic model is presented, followed by a summary of a linear modeling methodology that includes test and identification methods, and then a progressive series of steps suggesting a roadmap to develop a general nonlinear methodology that defines modeling, testing, and identification methods. Initial steps of the general methodology were applied to static and oscillatory test data to identify rolling-moment coefficient. Static measurements uncovered complicated dependencies of the aerodynamic coefficient on angle of attack and sideslip in the stall region making it difficult to find a simple analytical expression for the measurement data. In order to assess the effect of sideslip on the damping and unsteady terms, oscillatory tests in roll were conducted at different values of an initial offset in sideslip. Candidate runs for analyses were selected where higher order harmonics were required for the model and where in-phase and out-of-phase components varied with frequency. From these results it was found that only data in the angle-of-attack range of 35 degrees to 37.5 degrees met these requirements. From the limited results it was observed that the identified models fit the data well and both the damping-in-roll and the unsteady term gain are decreasing with increasing sideslip and motion amplitude. Limited similarity between parameter values in the nonlinear model and the linear model suggest that identifiability of parameters in both terms may be a problem. However, the proposed methodology can still be used with careful experiment design and carefully selected values of angle of attack, sideslip, amplitude, and frequency of the oscillatory data.

Murphy, Patrick C.↗
DOE OSTI2017

Quantitative assessment of soil CO 2 concentration and stable carbon isotope for leakage detection at geological carbon sequestration sites

This paper presents a quantitative method to evaluate CO 2 concentration ([CO 2 ]) and stable carbon isotope ratio (δ 13 C) as indicators for leakage detection at a geological carbon sequestration site by combining use of field release tests and a numerical modeling approach. A numerical model was developed to simulate CO 2 dynamics by considering diffusion, dissolution in soil water, and soil respiration. The numerical model fits the background dynamics of [CO 2 ] (360 to 550 ppm) and δ 13 C (-16‰ to -6‰) well and reproduces fairly the overall trend observed during the CO 2 release test. The model was further applied to assess detection probability (DP) of [CO 2 ] and δ 13 C for leakage detection in terms of various factors, such as CO 2 leakage rate, background variations, δ 13 C of the CO 2 leaked, and the threshold value of signal-to-noise ratio. Modeling results suggest that δ 13 C may have a higher DP than [CO 2 ]. This study also shows that DP of δ 13 C for the IEA Weyburn project is close to 0, implying δ 13 C is inappropriate to be an indicator for CO 2 leakage at the site. Finally, the quantitative method developed can also be used to design a monitoring plan or strategy in a near-surface environment for geological carbon sequestration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Barium stars as tracers of s -process nucleosynthesis in AGB stars: II. Using machine learning techniques on 169 stars

Barium (Ba) stars are characterised by an abundance of heavy elements made by the slow neutron capture process (s-process). This peculiar observed signature is due to the mass transfer from a stellar companion, bound in a binary stellar system, to the Ba star observed today. The signature is created when the stellar companion is an asymptotic giant branch (AGB) star. We aim to analyse the abundance pattern of 169 Ba stars using machine learning techniques and the AGB final surface abundances predicted by the FRUITY and Monash stellar models. We developed machine learning algorithms that use the abundance pattern of Ba stars as input to classify the initial mass and metallicity of each Ba star’s companion star using stellar model predictions. We used two algorithms. The first exploits neural networks to recognise patterns, and the second is a nearest-neighbour algorithm that focuses on finding the AGB model that predicts the final surface abundances closest to the observed Ba star values. In the second algorithm, we included the error bars and observational uncertainties in order to find the best-fit model. The classification process was based on the abundances of Fe, Rb, Sr, Zr, Ru, Nd, Ce, Sm, and Eu. We selected these elements by systematically removing s-process elements from our AGB model abundance distributions and identifying the elements whose removal had the biggest positive effect on the classification. We excluded Nb, Y, Mo, and La. Our final classification combined the output of both algorithms to identify an initial mass and metallicity range for each Ba star companion. With our analysis tools, we identified the main properties for 166 of the 169 Ba stars in the stellar sample. The classifications based on both stellar sets of AGB final abundances show similar distributions, with an average initial mass of M = 2.23 M ⊙ and 2.34 M ⊙ and an average [Fe/H] = –0.21 and –0.11, respectively. We investigated why the removal of Nb, Y, Mo, and La improves our classification and identified 43 stars for which the exclusion had the biggest effect. We found that these stars have statistically significant and different abundances for these elements compared to the other Ba stars in our sample. We discuss the possible reasons for these differences in the abundance patterns.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Neptune's Latitudinal Variations as Viewed with ALMA

We present spatially resolved millimeter maps of Neptune between 95 and 242 GHz taken with the Atacama Large Millimeter/submillimeter Array (ALMA) in 2016–2017. The millimeter weighting functions peak between 1 and 10 bar on Neptune, lying in between the altitudes probed at visible/infrared and centimeter wavelengths. Thus, these observations provide important constraints on the atmospheric structure and dynamics of Neptune. We identify seven well-resolved latitudinal bands of discrete brightness temperature variations, on the order of 0.5–3 K in all three observed ALMA spectral bands. We model Neptune’s brightness temperature using the radiative-transfer code Radio-BEAR and compare how various H{sub 2}S, CH{sub 4}, and ortho-/para-H{sub 2} abundance profiles can fit the observed temperature variations across the disk. We find that observed variations in brightness temperature with latitude can be explained by variations in the H{sub 2}S profile that range from sub- to supersaturations at altitudes above the 10 bar pressure level, while variations in CH{sub 4} improve the quality of fit near the equator. At the south polar cap, our best-fit model has a depleted deep atmospheric abundance of H{sub 2}S from 30 to only 1.5 times the protosolar value, while simultaneously depleting the CH{sub 4} abundance. This pattern of enhancement and depletion of condensible species is consistent with a global circulation structure where enriched air rises at the midlatitudes (32°–12°S) and north of the equator (2°–20°N), and dry air descends at the poles (90°–66°S) and just south of the equator (12°S–2°N). Our analysis finds more complex structure near the equator than accounted for in previous circulation models.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2012

Calibration of Two-dimensional Floodplain Modeling in the Atchafalaya River Basin Using SAR Interferometry

Two-dimensional (2D) satellite imagery has been increasingly employed to improve prediction of floodplain inundation models. However, most focus has been on validation of inundation extent, with little attention on the 2D spatial variations of water elevation and slope. The availability of high resolution Interferometric Synthetic Aperture Radar (InSAR) imagery offers unprecedented opportunity for quantitative validation of surface water heights and slopes derived from 2D hydrodynamic models. In this study, the LISFLOOD-ACC hydrodynamic model is applied to the central Atchafalaya River Basin, Louisiana, during high flows typical of spring floods in the Mississippi Delta region, for the purpose of demonstrating the utility of InSAR in coupled 1D/2D model calibration. Two calibration schemes focusing on Manning s roughness are compared. First, the model is calibrated in terms of water elevations at a single in situ gage during a 62 day simulation period from 1 April 2008 to 1 June 2008. Second, the model is calibrated in terms of water elevation changes calculated from ALOS PALSAR interferometry during 46 days of the image acquisition interval from 16 April 2008 to 1 June 2009. The best-fit models show that the mean absolute errors are 3.8 cm for a single in situ gage calibration and 5.7 cm/46 days for InSAR water level calibration. The optimum values of Manning's roughness coefficients are 0.024/0.10 for the channel/floodplain, respectively, using a single in situ gage, and 0.028/0.10 for channel/floodplain the using SAR. Based on the calibrated water elevation changes, daily storage changes within the size of approx 230 sq km of the model area are also calculated to be of the order of 107 cubic m/day during high water of the modeled period. This study demonstrates the feasibility of SAR interferometry to support 2D hydrodynamic model calibration and as a tool for improved understanding of complex floodplain hydrodynamics

Hahn Chul Jung↗
DOE OSTI2015

The luminosity and stellar mass functions of GRB host galaxies: Insight into the metallicity bias

Long-duration gamma-ray bursts (GRBs) are powerful probes of the star formation history of the universe, but the correlation between the two depends on the highly debated presence and strength of a metallicity bias. To investigate this correlation, we use a phenomenological model that successfully describes star formation rates, luminosities, and stellar masses of star-forming galaxies and apply it to GRB production. We predict the luminosities, stellar masses, and metallicities of host galaxies depending on the presence (or absence) of a metallicity bias. Our best-fitting model includes a moderate metallicity bias, broadly consistent with the large majority of the long-duration GRBs in metal-poor environments originating from a collapsar (probability ∼83%, with [0.74;0.91] range at 90% confidence level), but with a secondary contribution (∼17%) from a metal-independent production channel, such as binary evolution. Because of the mass–metallicity relation of galaxies, the maximum likelihood model predicts that the metal-independent channel becomes dominant at z≲2, where hosts have higher metallicities and collapsars are suppressed. This possibly explains why some studies find no clear evidence of a metal bias based on low-z samples. However, while metallicity predictions match observations well at high redshift (z≳2), there is tension with low-redshift observations, since a significant fraction of GRB hosts are predicted to have (near) solar metallicity. This is in contrast to observations, unless obscured, metal-rich hosts are preferentially missed in current data sets, and suggests that lower efficiencies of the metal-independent GRB channel might be preferred following a comprehensive fit that includes metallicity of GRB hosts from complete samples. Overall, we are able to clearly establish the presence of a metallicity bias for GRB production, but continued characterization of GRB host galaxies is needed to quantify its strength. Tabulated model predictions are available in electronic format.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Toward Accurate Modeling of Galaxy Clustering on Small Scales: Constraining the Galaxy-halo Connection with Optimal Statistics

Applying halo models to analyze the small-scale clustering of galaxies is a proven method for characterizing the connection between galaxies and their host halos. Such works are often plagued by systematic errors or limited to clustering statistics that can be predicted analytically. In this work, we employ a numerical mock-based modeling procedure to examine the clustering of Sloan Digital Sky Survey DR7 galaxies. We apply a standard halo occupation distribution (HOD) model to dark matter only simulations with a ΛCDM cosmology. To constrain the theoreStical models, we utilize a combination of galaxy number density and selected scales of the projected correlation function, redshift-space correlation function, group multiplicity function, average group velocity dispersion, mark correlation function, and counts-in-cells statistics. We design an algorithm to choose an optimal combination of measurements that yields tight and accurate constraints on our model parameters. Compared to previous work using fewer clustering statistics, we find a significant improvement in the constraints on all parameters of our halo model for two different luminosity-threshold galaxy samples. Most interestingly, we obtain unprecedented high-precision constraints on the scatter in the relationship between galaxy luminosity and halo mass. However, our best-fit model results in significant tension (>4σ) for both samples, indicating the need to add second-order features to the standard HOD model. To guarantee the robustness of these results, we perform an extensive analysis of the systematic and statistical errors in our modeling procedure, including a first of its kind study of the sensitivity of our constraints to changes in the halo mass function due to baryonic physics.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2007

Status on Iterative Transform Phase Retrieval Applied to the GBT Data

This slide presentation reviews the use of iterative transform phase retrieval in the analysis of the Green Bank Radio Telescope (GBT) Data. It reviews the NASA projects that have used phase retrieval, and the testbed for the algorithm to be used for the James Webb Space Telescope. It shows the comparison of phase retrieval with an interferometer, and reviews the two approaches used for phase retrieval, iterative transform (ITA) or parametric (non-linear least squares model fitting). The concept of ITA Phase Retrieval is reviewed, and the application to Radio Antennas is reviewed. The presentation also examines the National Radio Astronomy Observatory (NRAO) data from the GBT, and the Fourier model that NRAO uses to analyze the data. The challenge for ITA phase retrieval is reviewed, and the coherent approximation for incoherent data is shown. The validity of the approximation is good for a large tilt. There is a review of the proof of concept of the Phase Review simulation using the input wavefront, and the initial sampling parameters estimate from the focused GBT data.

Dean, Bruce↗
NASA NTRS2016

Global Comparison of CFD and Wind-Tunnel-Derived Force and Moment Databases for the Space Launch System

Recently a very large (739 runs) collection of high-fidelity RANS CFD solutions was obtained for Space Launch System ascent aerodynamics for the vehicle to be used for the first exploratory (unmanned) mission (EM-1). The extensive computations, at full-scale conditions, were originally developed to obtain detailed line and protuberance loads and surface pressures for venting analyses. The line loads were eventually integrated for comparison of the resulting forces and moments to the database that was derived from wind tunnel tests conducted at sub-scale conditions. The comparisons presented herein cover the ranges 0.5 < or = M(infinity) < or = 5, −6deg < or = alpha < or = 6deg, and −6deg < or = beta < or = 6deg. For detailed comparisons, slender-body-theory-based component build-up aero models from missile aerodynamics are used. The differences in the model fit coefficients are shown to be relatively small except for the low supersonic Mach number range, 1.1 < or = M(infinity) < or = 2.0. The analysis is intended to support process improvement and development of uncertainty models.

Hemsch, Michael J.↗
DOE OSTI2016

Predictive models of radiative neutrino masses

We discuss two models of radiative neutrino mass generation. The first model features one–loop Zee model with Z{sub 4} symmetry. The second model is the two–loop neutrino mass model with singly- and doubly-charged scalars. These two models fit neutrino oscillation data well and predict some interesting rates for lepton flavor violation processes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

CME PROPAGATION: WHERE DOES AERODYNAMIC DRAG “TAKE OVER”?

We investigate the Sun–Earth dynamics of a set of eight well observed solar coronal mass ejections (CMEs) using data from the Solar Terrestrial Relations Observatory spacecraft. We seek to quantify the extent to which momentum coupling between these CMEs and the ambient solar wind (i.e., the aerodynamic drag) influences their dynamics. To this end, we use results from a 3D flux rope model fit to the CME data. We find that solar wind aerodynamic drag adequately accounts for the dynamics of the fastest CME in our sample. For the relatively slower CMEs, we find that drag-based models initiated below heliocentric distances ranging from 15 to 50 R{sub ⊙} cannot account for the observed CME trajectories. This is at variance with the general perception that the dynamics of slow CMEs are influenced primarily by solar wind drag from a few R{sub ⊙} onwards. Several slow CMEs propagate at roughly constant speeds above 15–50 R{sub ⊙}. Drag-based models initiated above these heights therefore require negligible aerodynamic drag to explain their observed trajectories.

79 ASTRONOMY AND ASTROPHYSICS↗