Search NASA⌕ Search

SEARCH · Search NASA

Results for “Parameter estimation using data”

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 397 records · Page 22

The catalog-to-cosmology framework for weak lensing and galaxy clustering for LSST

We present TXPipe, a modular, automated and reproducible pipeline for ingesting catalog data and performing all the calculations required to obtain quality-assured two-point measurements of lensing and clustering, and their covariances, with the metadata necessary for parameter estimation. The pipeline is developed within the Rubin Observatory Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC), and designed for cosmology analyses using LSST data. In this paper, we present the pipeline for the so-called 3x2pt analysis -- a combination of three two-point functions that measure the auto- and cross-correlation between galaxy density and shapes. We perform the analysis both in real and harmonic space using TXPipe and other LSST-DESC tools. We validate the pipeline using Gaussian simulations and show that it accurately measures data vectors and recovers the input cosmology to the accuracy level required for the first year of LSST data under this simplified scenario. We also apply the pipeline to a realistic mock galaxy sample extracted from the CosmoDC2 simulation suite (Korytov et al. 2019). TXPipe establishes a baseline framework that can be built upon as the LSST survey proceeds. Furthermore, the pipeline is designed to be easily extended to science probes beyond the 3x2pt analysis.

79 ASTRONOMY AND ASTROPHYSICS↗

Dissipation and Bathymetric Sensitivities in an Unstructured Mesh Global Tidal Model

Abstract The mechanisms and geographic distribution of global tidal dissipation in barotropic tidal models are examined using a high resolution unstructured mesh finite element model. Mesh resolution varies between 2 and 25 km and is especially focused on inner shelves and steep bathymetric gradients. Tidal response sensitivities to bathymetric changes are examined to put into context response sensitivities to frictional processes. We confirm that the Ronne Ice Shelf dramatically affects Atlantic tides but also find that bathymetry in the Hudson Bay system is a critical control. We follow a sequential frictional parameter optimization process and use TPXO9 data‐assimilated tidal elevations as a reference solution. From simulated velocities and depths, dissipation within the global model is estimated and allows us to pinpoint dissipation at high resolution. Boundary layer dissipation is extremely focused with 1.4% of the ocean accounting for 90% of the total. Internal tide friction is much more distributed with 16.7% of the ocean accounting for 90% of the total. Often highly regional dissipation can impact basin‐scale and even ocean wide tides. Optimized boundary layer friction parameters correlate very well with the physical characteristics of the locality with high friction factors associated with energetic tidal regions, deep ocean island chains, and ice covered areas. Global complex M 2 tide errors are 1.94 cm in deep waters. Total global boundary layer and internal tide dissipation are estimated, respectively, at 1.83 and 1.49 TW. This continues the trend in the literature toward attributing more dissipation to internal tides.

54 ENVIRONMENTAL SCIENCES↗

Identifying mitigation strategies for $\mathrm{COVID}$-$\mathrm{19}$ superspreading on flights using models that account for passenger movement

Despite commercial airlines mandating masks, there have been multiple documented events of COVID-19 superspreading on flights. Conventional models do not adequately explain superspreading patterns on flights, with infection spread wider than expected from proximity based on passenger seating. An important reason for this is that models typically do not consider the movement of passengers during the flight, boarding, or deplaning. Understanding the risks for each of these aspects could provide insight into effective mitigation measures. We modeled infection risk from seating and fine-grained movement patterns – boarding, deplaning, and inflight movement. We estimated infection model parameters from a prior superspreading event. We validated the model and the impact of interventions using available data from three flights, including cabin layout and seat locations of infected and uninfected passengers, to suggest interventions to mitigate COVID-19 superspreading events during air travel. Specifically, we studied: 1) London to Hanoi with 201 passengers, including 13 secondary infections among passengers; 2) Singapore to Hangzhou with 321 passengers, including 12 to 14 secondary infections; 3) a non-superspreading event on a private jet in Japan with 9 passengers and no secondary infections. Here our results show that the inclusion of passenger movement better explains the infection spread patterns than conventional models do. We also found that FFP2/N95 mask usage would have reduced infection by 95–100%, while cloth masks would have reduced it by only 40–80%. Results indicate that leaving the middle seat vacant is effective in reducing infection, and the effectiveness increases when combined with good quality masks. However, with a good mask, the risk is quite low even without the middle seats being empty. Our results suggest the need for more stringent guidelines to reduce aviation-related superspreading events of COVID-19.

54 ENVIRONMENTAL SCIENCES↗

110th Anniversary: Microkinetic Modeling of the Vapor Phase Upgrading of Biomass-Derived Oxygenates

Bio-oil produced from fast pyrolysis of biomass is a complex mixture of more than 200 compounds, including oxygenates and acids. As these species are highly undesirable in fuels, catalytic upgrading of biomass pyrolysis product vapors, also known as catalytic fast pyrolysis, is performed to upgrade the vapors to valuable fuels and chemicals. This work presents a detailed microkinetic model, composed of elementary steps, of the catalytic upgrading of acetic acid and acetone, two common oxygenates present in bio-oil. An automated network generator was utilized to construct a reaction network composed of 580 unique species and 2160 unique reactions. The kinetic parameters for each reaction in the network were estimated using transition state theory, the Evans-Polanyi relationship, and thermodynamic data. The resulting mechanistic model is able to describe experimental data presented in the literature for the transformation of acetic acid and acetone on HZSM-5 in a fixed-bed reactor, which is modeled as a plug-flow reactor. Additionally, the model solutions reveal vital information regarding the mechanism by which acetic acid and acetone are upgraded to valuable fuels and chemicals. In the first phase of the mechanism, acetic acid is converted to acetone via acylium ion addition to acetic acid; this is followed by decarboxylation of acetoacetic acid. The second phase is dominated by the self-aldol condensation of acetone, which is shown to occur predominantly through the keto form of acetone rather than the enol form, and subsequent deoxygenation reactions leading to olefins and aromatics. Finally, net rate analysis shows that aromatics are primarily formed via a pathway including aldol condensation of mesityl oxide, whereas olefins are produced from the addition of isobutene and subsequent cracking.

09 BIOMASS FUELS↗

Hybrid em wave - polar semiconductor interaction: A polaronic study

Present paper considers incidence of a most realistic hybrid pump wave on a weakly polar semiconductor having a very small coupling constant. Possibility of optical parametric interaction has been explored in the presence of an external transverse magnetic field. The effect of doping concentrations and transverse magnetostatic field on threshold characteristics of optical parametric interaction in polar semiconductor plasma has been studied, using hydrodynamic model of semiconductors, in the far infrared regime. Numerical estimations have been carried out by using data of weakly polar III-V GaAs semiconductor and influence of control parameters on electron-LO phonon interaction has been analyzed. A particular range of physical parameters is found to be suitable for minimum threshold. The choice of nonlinear medium and favorable range of operating parameters are crucial aspects in design and fabrication of parametric amplifiers and oscillators. The hybrid mode of the pump is found to be favorable for the onset of the said process and realization of a low cost amplifier.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Preliminary design overview and performance assessment of the ITER low-field side reflectometer

The design of the ITER low-field side reflectometer (LFSR) has matured to a complete end-to-end preliminary design. LFSR will supply three key plasma measurements: (1) electron density profile, (2) electron density fluctuations, and (3) poloidal rotation. Simultaneous measurements of the three quantities are enabled by an array of six monostatic antennas which inject from an equatorial port on the outboard side of the ITER vessel. Low-loss transmission lines, consisting of corrugated, overmoded waveguide and miter bends, transmit the 30–165 GHz, O- and X-mode signals to and from the ITER plasma. Integrated transmission-line components serve a range of purposes, such as protection from high-power stray radiofrequency radiation, accommodation of transmission-line displacement, and simultaneous measurement of reference and plasma phases during the discharge. Broadband transmission signals are realized by full-band microwave transceivers combined with quasi-optical multiplexing. A field-programmable gate array (FPGA) processor demodulates the profile reflectometer signals, enabling real-time density profile measurements for plasma control system feedback. A full-scale transmission line test facility provides an integrated environment to assess the performance of critical LFSR components. Theoretical modeling together with insertion loss measurements provide the basis for a comprehensive power budget, which accounts for transmitted output power, transmission-line losses, antenna coupling, and plasma effects. Results indicate that high signal-to-noise ratios are achievable with the current design. A synthetic reflectometer model, using real design parameters and baseline ITER profiles, has been developed to estimate the return signal. With evolving microwave and data acquisition technologies, full-band, ultrafast sweeps (<1 μs) will be realizable for ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Variable stars in Local Group galaxies – V. The fast and early evolution of the low-mass Eridanus II dSph galaxy

ABSTRACT We present a detailed study of the variable star population of Eridanus II (Eri II), an ultra-faint dwarf galaxy that lies close to the Milky Way virial radius. We analyse multi-epoch g, r, i ground-based data from Goodman and the Dark Energy Camera, plus F475W, F606W, F814W space data from the Advanced Camera for Surveys. We report the detection of 67 RR Lyrae (RRL) stars and 2 Anomalous Cepheids, most of them new discoveries. With the RRL stars, we measure the distance modulus of Eri II, μ0 = 22.84 ± 0.05 mag (D⊙ = 370 ± 9 kpc), and derive a metallicity spread of 0.3 dex (0.2 dex intrinsic). The colour distribution of the horizontal branch (HB) and the period distribution of the RRL stars can be nicely reproduced by a combination of two stellar models of [Fe/H] = (−2.62, −2.14). The overall low metallicity is consistent with the red giant branch bump location, 0.65 mag brighter than the HB. These results are in agreement with previous spectroscopic studies. The more metal-rich RRL and the RRab stars have greater central concentration than the more metal-poor RRL and the RRc stars that are mainly located outside ∼1 rh. This is similar to what is found in larger dwarf galaxies such as Sculptor, and in agreement with an outside-in galaxy formation scenario. This is remarkable in such a faint dwarf galaxy with an apparently single and extremely short (<1 Gyr) star formation burst. Finally, we have derived new and independent structural parameters for Eri II and its star cluster using our new data that are in very good agreement with previous estimates.

79 ASTRONOMY AND ASTROPHYSICS↗

Market-Based Coordination of Thermostatically Controlled Loads—Part II: Unknown Parameters and Case Studies

This two-part paper considers the coordination of a population of Thermostatically Controlled Loads (TCLs) with unknown parameters to achieve group objectives. The problem involves designing the bidding and market clearing strategy to motivate self-interested users to realize efficient energy allocation subject to a peak power constraint. The companion paper (Part I) formulates the problem and proposes a load coordination framework using the mechanism design approach. To address the unknown parameters, Part II of this paper presents a joint state and parameter estimation framework based on the expectation maximization algorithm. The overall framework is then validated using real-world weather data and price data, and is compared with other approaches in terms of aggregated power response. Simulation results indicate that our coordination framework can effectively improve the efficiency of the power grid operations and reduce power congestion at key times.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

SU-F-T-132: Variable RBE Models Predict Possible Underestimation of Vaginal Dose for Anal Cancer Patients Treated Using Single-Field Proton Treatments

Purpose: Anal cancer patients treated using a posterior proton beam may be at risk of vaginal wall injury due to the increased linear energy transfer (LET) and relative biological effectiveness (RBE) at the beam distal edge. We investigate the vaginal dose received. Methods: Five patients treated for anal cancer with proton pencil beam scanning were considered, all treated to a prescription dose of 54 Gy(RBE) over 28–30 fractions. Dose and LET distributions were calculated using the Monte Carlo simulation toolkit TOPAS. In addition to the standard assumption of a fixed RBE of 1.1, variable RBE was considered via the application of published models. Dose volume histograms (DVHs) were extracted for the planning treatment volume (PTV) and vagina, the latter being used to calculate the vaginal normal tissue complication probability (NTCP). Results: Compared to the assumption of a fixed RBE of 1.1, the variable RBE model predicts a dose increase of approximately 3.3 ± 1.7 Gy at the end of beam range. NTCP parameters for the vagina are incomplete in the current literature, however, inferring value ranges from the existing data we use D{sub 50} = 50 Gy and LKB model parameters a=1–2 and m=0.2–0.4. We estimate the NTCP for the vagina to be 37–48% and 42–47% for the fixed and variable RBE cases, respectively. Additionally, a difference in the dose distribution was observed between the analytical calculation and Monte Carlo methods. We find that the target dose is overestimated on average by approximately 1–2%. Conclusion: For patients treated with posterior beams, the vaginal wall may coincide with the distal end of the proton beam and may receive a substantial increase in dose if variable RBE models are applied compared to using the current clinical standard of RBE equal to 1.1. This could potentially lead to underestimating toxicities when treating with protons.

60 APPLIED LIFE SCIENCES↗

XCal: model-based approach to X-ray CT spectral calibration

Transmission X-ray computed tomography (CT) is widely used to quantitatively reconstruct 3D objects composed of multiple materials. However, accurate CT reconstruction requires the system to be calibrated to account for the effective X-ray spectrum. Unfortunately, measurement of the effective spectrum is ill-posed, and existing calibration methods require that the system be recalibrated when the system parameters are changed. In this paper, we propose XCal, a multi-energy model-based spectral calibration approach for X-ray CT. The XCal approach models the effective spectrum using a separable physics-based model of the CT system. The model parameters are then estimated by fitting calibration data with known objects at multiple energies. An important advantage of XCal is that it allows the user to change scanner settings, such as the source voltage or X-ray filters, without the need for recalibration. Evaluations on simulated and measured datasets demonstrate that XCal significantly improves the accuracy of the estimated spectrum as compared to existing calibration methods.

Li, Wenrui [Purdue Univ., West Lafayette, IN (Unit↗

Validation of the Public Radiation Exposure Calculation for the Incident at the National Institute of Standards and Technology Center for Neutron Research on February 3, 2021

The Department of Energy/National Nuclear Security Administration Consequence Management Program was contacted by the Health Physics Chief of the National Institute of Standards and Technology Center for Neutron Research (NCNR) to review public radiation exposure calculations for an event that occurred on its Gaithersburg, Maryland, campus on February 3, 2021. Subject matter experts from the Nuclear Emergency Support Team (NEST) assets, specifically the Consequence Management Home Team (CMHT) and the National Atmospheric Release Advisory Center (NARAC), were selected to provide support. CMHT used three separate modeling codes to validate the results the scientists at NCNR calculated using the HotSpot model. The analyses were performed using NARAC’s in-house Lagrangian dispersion codes known as LODI and Aeolus, as well as the Turbo FRMAC software from Sandia National Laboratories. The team used parameters provided by the NCNR scientists regarding the site, applicable observable meteorological data, and environmental survey and sampling data to estimate public exposure. Each model estimated public dose at much less than 0.5 mrem. CMHT concurs with the NCNR public radiation exposure calculations which state that members of the public at the 400-meter boundary would have received a radiological dose of less than 0.5 mrem.

61 RADIATION PROTECTION AND DOSIMETRY↗

Large-scale Nonlinear Approaches for Inference of Reporting Dynamics and Unobserved SARS-CoV-2 Infections

This work focuses on estimation of unknown states and parameters in a discrete-time, stochastic, SEIR model using reported case counts and mortality data. An SEIR model is based on classifying individuals with respect to their status in regards to the progression of the disease, where S is the number individuals who remain susceptible to the disease, E is the number of individuals who have been exposed to the disease but not yet infectious, I is the number of individuals who are currently infectious, and R is the number of recovered individuals. For convenience, we include in our notation the number of infections or transmissions, T, that represents the number of individuals transitioning from compartment S to compartment E over a particular interval. Similarly, we use C to represent the number of reported cases.

59 BASIC BIOLOGICAL SCIENCES↗

Use of Satellite, Surface Observations and Numerical Weather Prediction Model Data to Improve Cloud Base Height and Cloud Base Vertical Velocity Estimation

Cloud base height (CBH) and cloud base vertical velocity (CBVV) are important variables that impact the overall climate in a region as they influence the formulation, longevity, and evolution of clouds. Retrieval of both parameters have long used ground instrumentation (e.g., Doppler lidar (DL), ground base radar); however, retrieving CBH from satellites is particularly challenging given that space-based instruments only observe cloud tops. In this manuscript, CBH is retrieved using a multi-linear regression equation, while CBVV used a random forests model. Both retrievals combine satellite and numerical weather prediction data. The satellite data used are the Visible Infrared Imaging Radiometer Suite imagery, while measurements of CBH and CBVV include DL and radiosonde data at the Southern Great Plains (SGP) Atmospheric Radiation Measurement observatory. Data from 83 summer days (May-August) in 2018–2021 featuring cumulus clouds forced by solar heating were examined and used to train the models, with years 2022–2023 used for validation. Various spatial domains were defined with one large (2.4° longitude by 2.0° latitude) SGP domain being split into smaller sections (smallest being 0.99° and 0.61° longitude and latitude respectably). CBH and CBVV values obtained from the DL as compared to the models show root mean square errors between 150 and 200 m, with CBVV values between 0.45 and 1 ms -1 . Finally, it was found that the CBH formulation performs well over all domains, while the CBVV retrievals become less accurate due to more turbulence being introduced into the observations as the number of DL stations decreases in the smaller domains.

54 ENVIRONMENTAL SCIENCES↗

Comparison of source-location algorithms for atmospheric samplers

Numerous algorithms have been developed to determine the source characteristics for an atmospheric release, e.g., (Bieringer et al., 2017). This study compares three models that have been applied to the data collected by the International Monitoring System operated by the Comprehensive Nuclear-Test-Ban Treaty Preparatory Commission to estimate source event parameters. Each model uses a different approach to estimate the parameters. A deterministic model uses a possible source region (PSR) approach (Ringbom et al., 2014) that is based on the correlation between predicted and measured sample values. An unnamed model developed at Pacific Northwest National Laboratory (herein called BAYEST) uses a Bayesian formulation (Eslinger et al., 2019; Eslinger et al., 2020; Eslinger and Schrom, 2016). The FREAR model uses a different Bayesian formulation (De Meutter and Hoffman, 2020; De Meutter et al., 2021; Pieter De Meutter, 2021). The performance of the three source-location models is evaluated with 100 synthetic release cases for the single xenon isotope, 133 Xe. The release cases resulted in detections in a fictitious network with 120 noble gas samplers. All three source-location models use the same sampling data. The two Bayesian models yield more accurate location estimates than the deterministic PSR model, with FREAR having slightly better location performance than BAYEST. Samplers with collection periods of 3, 6, 8, 12, and 24-h were used. Results from BAYEST show that location accuracy improves with each reduction in sample collection length. The BAYEST model is slightly better for estimating the start time of the release. The PSR model has about the same spread in start times as the FREAR model, but the PSR results have a better average start time. The Bayesian source-location algorithms give more accurate results than the PSR approach, and provide release magnitude estimates, while the base PSR model does not estimate the release magnitude. This investigation demonstrates that a reasonably dense sampling grid will sometimes yield poor location and time estimates regardless of the model. The poor estimates generally coincide with cases where there is a much larger distance between the release point and the first detecting sampler than the average sampler spacing.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

LAMOST Time-Domain survey: first results of four K2 plates

From Oct. 2019 to Apr. 2020, LAMOST performed a time-domain (TD) spectroscopic survey of four K2 plates with both low- and medium-resolution observations. The low-resolution spectroscopic survey acquired 282 exposures ( ≈ 46.6 h) over 25 nights, yielding a total of about 767 000 spectra, and the medium-resolution survey took 177 exposures ( ≈ 49.1 h) over 27 nights, collecting about 478 000 spectra. More than 70%/50% of low-resolution/medium-resolution spectra have signal-to-noise ratio higher than 10. We determine stellar parameters (e.g., T eff , log g, [Fe/H]) and radial velocity (RV) with different methods, including LASP, DD-Payne and SLAM. In general, these parameter estimations from different methods show good agreement, and the stellar parameter values are consistent with those of APOGEE. We use the Gaia DR2 RV values to calculate a median RV zero point (RVZP) for each spectrograph exposure by exposure, and the RVZP-corrected RVs agree well with the APOGEE data. The stellar evolutionary and spectroscopic masses are estimated based on the stellar parameters, multi-band magnitudes, distances and extinction values. Finally, we construct a binary catalog including about 2700 candidates by analyzing their light curves, fitting the RV data, calculating the binarity parameters from medium-resolution spectra and cross-matching the spatially resolved binary catalog from Gaia EDR3. In conclusion, the LAMOST TD survey is expected to represent a breakthrough in various scientific topics, such as binary systems, stellar activity, stellar pulsation, etc.

79 ASTRONOMY AND ASTROPHYSICS↗

Estimating a Service-Life Distribution Based on Production Counts and a Failure Database

A manufacturer wanted to compare the service-life distributions of two similar products. These concern product lifetimes after installation (not manufacture). For each product, there were available production counts and an imperfect database providing information on failing units. In the real case, these units were expensive repairable units warrantied against repairs. Failure (of interest here) was relatively rare and driven by a different mode/mechanism than ordinary repair events (not of interest here). Approach: Data models for the service life based on a standard parametric lifetime distribution and a related limited failure population were developed. These models were used to develop expressions for the likelihood of the available data that properly accounts for information missing in the failure database. Results: A Bayesian approach was employed to obtain estimates of model parameters (with associated uncertainty) in order to investigate characteristics of the service-life distribution. Custom software was developed and is included as Supplemental Material to this case study. One part of a responsible approach to the original case was a simulation experiment used to validate the correctness of the software and the behavior of the statistical methodology before using its results in the application, and an example of such an experiment is included here. Because of confidentiality issues that prevent use of the original data, simulated data with characteristics like the manufacturer’s proprietary data are used to illustrate some aspects of our real analyses. Lastly, we also note that, although this case focuses on rare and complete product failure, the statistical methodology provided is directly applicable to more standard warranty data problems involving typically much larger warranty databases where entries are warranty claims (often for repairs) rather than reports of complete failures.

97 MATHEMATICS AND COMPUTING↗

Multi-area parameter error identification for large power systems

Power grid model parameters may contain errors due to various reasons. Detecting and correcting parameter errors typically requires significant computational effort due to the size and complexity of the parameter database. While the normalized Lagrange multiplier (NLM) method can effectively detect, identify and correct parameter errors, its computational burden could rapidly grow with increasing system size. This paper addresses this issue by proposing a multi-area parameter error identification method. Each area has its own outlier detection tool for detecting the incorrect parameters and measurements within the area. On the other hand, due to the reduced redundancy at area boundaries, parameter errors on branches incident to boundary buses may not be detected. Such errors are subsequently detected by a coordination level estimator completing the system-wide parameter detection procedure. In conclusion, performance of the developed method is demonstrated using the IEEE 118-bus and 2000-bus Texas synthetic systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A Robust State Estimation Framework Considering Measurement Correlations and Imperfect Synchronization

This paper develops a robust power system state estimation framework with the consideration of measurement correlations and imperfect synchronization. In the framework, correlations of SCADA and PMU measurements are calculated separately through unscented transformation and a vector autoregression (VAR) model. In particular, PMU measurements during the waiting period of two SCADA measurement scans are buffered to develop the VAR model. The latter takes into account the temporal and spatial correlations of PMU measurements and provides redundant measurements to suppress bad data and mitigate imperfect synchronization. In case of synchronization issues between SCADA and PMU measurements, either forecasted PMU measurements or prior SCADA measurements from the last estimation run are leveraged to restore system observability. After that, robust generalized maximum-likelihood (GM)-estimator is developed to integrate full measurement correlations and handle bad data in SCADA and PMU measurements. Thanks to the robustness and flexibility of the proposed robust estimator, bad data as well as imperfect measurement synchronization are mitigated, yielding high statistical efficiency. Comprehensive comparison results with other alternatives under various conditions demonstrate the benefits of the proposed framework.

Zhao, Junbo (ORCID:0000000284989666)↗