Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “timing error evaluation”

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 865 records · Page 48

Dynamic signal recovery in distribution grids using compressive lossy measurements

Distribution system state estimation requires reliable aggregation of the measured data. However, the large volume of the measured data imposes a significant stress on the underlying communication infrastructure. With the challenges associated with measurement availability, current distribution systems are typically unobservable. To cope with the unobservability issue, compressive sensing theory allows us to recover system state information from a small number of measurements provided the states of the distribution system exhibit sparsity. In this paper, we evaluate the robustness of an updated Kalman filtered modified compressive sensing (KF-ModCS) technique that dynamically estimates the grid states using a small fraction of measured data. In practice, measurements used for sparsity based state estimation may also be intermittent due to communication network induced losses. Further, to understand the effect of packet losses on KF-ModCS, we provide an upper bound for the expected variances of the state estimation error for a given rate of information loss. This upper bound is further improved if the support set of the sparse signal that characterizes the state dynamics does not change over time and/or the reduced model is observable. Simulations based on two practical data sets collected from actual customers in a distribution grid validate the theoretical results.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Coordination Sphere of Lanthanide Aqua Ions Resolved with Ab Initio Molecular Dynamics and X-ray Absorption Spectroscopy

To resolve the fleeting structures of lanthanide Ln 3+ aqua ions in solution, we (i) performed the first ab initio molecular dynamics (AIMD) simulations of the entire series of Ln 3+ aqua ions in explicit water solvent using pseudopotentials and basis sets recently optimized for lanthanides and (ii) measured the symmetry of the hydrating waters about Ln 3+ ions (Nd 3+ , Dy 3+ , Er 3+ , Lu 3+ ) for the first time with extended X-ray absorption fine structure (EXAFS). EXAFS spectra were measured experimentally and generated from AIMD trajectories to directly compare simulation, which concurrently considers the electronic structure and the atomic dynamics in solution, with experiment. In this study, we performed a comprehensive evaluation of EXAFS multiple-scattering analysis (up to 6.5 Å) to measure Ln–O distances and angular correlations (i.e., symmetry) and elucidate the molecular geometry of the first hydration shell. This evaluation, in combination with symmetry-dependent L 3 - and L 1 -edge spectral analysis, shows that the AIMD simulations remarkably reproduces the experimental EXAFS data. The error in the predicted Ln–O distances is less than 0.07 Å for the later lanthanides, while we observed excellent agreement with predicted distances within experimental uncertainty for the early lanthanides. Our analysis revealed a dynamic, symmetrically disordered first coordination shell, which does not conform to a single molecular geometry for most lanthanides. This work sheds critical light on the highly elusive coordination geometry of the Ln 3+ aqua ions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computational fluid dynamics analysis of char conversion in Sandia’s pressurized entrained flow reactor

Design and analysis of practical reactors utilizing solid feedstocks rely on reaction rate parameters that are typically generated in lab-scale reactors. Evaluation of the reaction rate information often relies on assumptions of uniform temperature, velocity, and species distributions in the reactor, in lieu of detailed measurements that provide local information. This assumption might be a source of substantial error, since reactor designs can impose significant inhomogeneities, leading to data misinterpretation. Spatially resolved reactor simulations help understand the key processes within the reactor and support the identification of severe variations of temperature, velocity, and species distributions. In this work, Sandia’s pressurized entrained flow reactor is modeled to identify inhomogeneities in the reaction zone. Tracer particles are tracked through the reactor to estimate the residence times and burnout ratio of introduced coal char particles in gasifying environments. The results reveal a complex mixing environment for the cool gas and particles entering the reactor along the centerline and the main high-speed hot gas reactor flow. Furthermore, the computational fluid dynamics (CFD) results show that flow asymmetries are introduced through the use of a horizontal gas pre-heating section that connects to the vertical reactor tube. Computed particle temperatures and residence times in the reactor differ substantially from the idealized plug flow conditions typically evoked in interpreting experimental measurements. Furthermore, experimental measurements and CFD analysis of heat flow through porous refractory insulation suggest that for the investigated conditions (1350 °C, <20 atm), the thermal conductivity of the insulation does not increase substantially with increasing pressure.

47 OTHER INSTRUMENTATION↗

Numerical evaluation of the Rayleigh integral for planar radiators using the FFT

Rayleigh's integral formula is evaluated numerically for planar radiators of any shape, with any specified velocity in the source plane using the fast Fourier transfrom algorithm. The major advantage of this technique is its speed of computation - over 400 times faster than a straightforward two-dimensional numerical integration. The technique is developed for computation of the radiated pressure in the nearfield of the source and can be easily extended to provide, with little computation time, the vector intensity in the nearfield. Computations with the FFT of the nearfield pressure of baffled rectangular plates with clamped and free boundaries are compared with the 'exact' solution to illuminate any errors. The bias errors, introduced by the FFT, are investigated and a technique is developed to significantly reduce them.

Williams, E. G.↗

Physics-based Simulation of Human Posture Using 3D Whole Body Scanning Technology for Astronaut Space Suit Evaluation

Over the past few years high precision three-dimensional (3D) full body laser scanners have been developed to be used as a powerful anthropometry tool for quantification of the morphology of the human body. The full body scanner can quickly extract body characteristics in non-contact fashion. It is required for the Anthropometry and Biomechanics Facility (ABF) to have capabilities for kinematics simulation of a digital human at various postures whereas the laser scanner only allows capturing a single static posture at each time. During this summer fellowship period a theoretical study has been conducted to estimate an arbitrary posture with a series of example postures through finite element (FE) approximation and found that four-point isoparametric FE approximation would result in reasonable maximum position errors less than 5%. Subsequent pilot scan experiments demonstrated that a bead marker with a nominal size of 6 mm could be used as a marker for digitizing 3-D coordinates of anatomical landmarks for further kinematic analysis. Two sessions of human subject testing were conducted for reconstruction of an arbitrary postures from a set of example postures for each joint motion for the forearm/hand complex and the whole upper extremity.

Kim, Kyu-Jung↗

Finite and infinite matrix product states for Gutzwiller projected mean-field wave functions

Matrix product states (MPS) and “dressed” ground states of quadratic mean fields (e.g., Gutzwiller projected Slater determinants) are both important classes of variational wave functions. This latter class has played important roles in understanding superconductivity and quantum spin liquids. We present a method to obtain both the finite and infinite MPS (iMPS) representation of the ground state of an arbitrary fermionic quadratic mean-field Hamiltonian (which in the simplest case is a Slater determinant and in the most general case is a Pfaffian). We also show how to represent products of such states (e.g., determinants times Pfaffians). From this representation one can project to single occupancy and evaluate the entanglement spectra after Gutzwiller projection. We then obtain the MPS and iMPS representation of Gutzwiller projected mean-field states that arise from the variational slave-fermion approach to the S = 1 bilinear-biquadratic quantum spin chain. To accomplish this, we develop an approach to orthogonalize degenerate iMPS to find all the states in the degenerate ground-state manifold. We find the energies of the MPS and iMPS states match the variational energies closely, indicating the method is accurate and there is minimal loss due to truncation error. We then present an exploration of the entanglement spectra of projected slave-fermion states, exploring their qualitative features and finding good qualitative agreement with the respective exact ground-state spectra found from density matrix renormalization group.

1-dimensional spin chains↗

Component pose reconstruction using a single robotic total station for panelized building envelopes

The construction industry can benefit greatly from increased automation in construction processes in terms of precision, time, and labor costs. This is especially true for prefabricated construction for either new buildings or envelope retrofits. Here, prefabricated components are manufactured offsite according to design specifications and installed onsite. Accurate information on the component’s position and orientation (pose) is needed to achieve this. A tool named “Real-Time Evaluator” (RTE), designed to autonomously track prefabricated components as they are being installed and provide real-time installation instructions, is currently under development using a single robotic total station. This is challenging since at least three points are needed to determine the pose of an object in space. To achieve this goal with a single robotic total station, two key algorithms were developed: the “resection” algorithm aligns the digital twin with the physical twin regardless of the location of the total station, and the “transformation” algorithm gives instructions of translations and rotations to installers to achieve the desired installation pose. The algorithms were evaluated experimentally on a lab-scale demonstration. Results show that the resection algorithm achieved an average error of < 3.1 mm, while the transformation algorithm predicted the rotation angle along a single axis with an error of < 0.5◦.

Tang, Mengjia↗

Frequency-wavenumber spectrum for GATE phase I rainfields

The oceanic rainfall frequency-wavenumber spectrum and its associated space-time correlation have been evaluated from subsets of GATE phase I data. The records, of a duration of four days, were sampled at 15 minutes intervals in 4 x 4 km grid boxes over a 400 km diameter hexagon. In the low frequencies-low wavenumber region the results coincide with those obtained by using the stochastic model proposed by North and Nakomoto (1989). From the derived spectrum the inherent time and space scales of the stochastic model were determined to be approximately 13 hours and 36 km. The formalism proposed by North and Nakamoto was taken together with the derived spectrum to compute the mean square sampling error due to intermittent visits of a spaceborne sensor.

Nakamoto, Shoichiro↗

Mission-Maps For Outbound Cislunar Transfer Trajectories

This study quantifies the robustness and sensitivity of an outbound cislunar trajectory for a lunar lander in the form of mission-maps, or topological maps that allows either a computer program or mission designer to intuitively optimize the placement of critical outbound correction burns from the derived sensitivity data. The non-linear multi-body dynamics are applied to generate an outbound cislunar reference profile used by a linear covariance analysis (LinCov) tool to compute the expected Δv and trajectory dispersions due to the initial state uncertainty, sensor errors, maneuver execution errors, and disturbance accelerations along the outbound cislunar profile. The rapid performance analysis capabilities of LinCov are complimented with parallel processing techniques to evaluate hundreds and thousands of different translational burn locations, placements, and targeting constraints to identify the combination that minimizes the total Δv usage (nominal plus 3σ Δv) and trajectory dispersions at lunar orbit insertion. This study utilizes a generalized reference targeting algorithm to quickly assess the integrated closed-loop GN&C system performance due to different targeting configurations and constraints. The resulting mission maps provide an intuitive insight to ascertain each trajectory correction maneuver’s (TCM) sensitivity to different burn times along an outbound cislunar trajectory and quickly identify desirable engineering tradeoffs when performing analysis on the number and placement of these burns that nominally zero. Multiple mission maps are generated for a variety of different performance parameters that allow engineers to visually identify optimal solutions for trajectory correction maneuver placements, the number of correction burns, and the targeting constraints for each burn.

GN&C↗

Evaluating Differences Among Crop Models in Simulating Soybean in-Season Growth

Crop models are useful tools for simulating agricultural systems that require continued model development and testing to increase their robustness and improve how they describe our current understanding of processes. Coordinated and “blind” evaluation of multiple models using same protocols and experimental datasets provides unique opportunities to further improve models and enhance their reliability. For soybean [Glycine max (L.) Merr.], there has been limited coordinated multi-model evaluations for the simulation of in-season plant growth dynamics. We evaluated ten dynamic soybean crop models for their simulation of in-season plant growth using data from five experiments conducted in Argentina, Brazil, France, and USA. We evaluated models after a Blind (using only phenology data) and a Full calibration (with in-season and end-of-season variables). Calibration reduced model uncertainty by reducing standard bias for the simulation of in-season variables (biomass, leaf, pod, and stem weights, and leaf area index, LAI). However, we found that most models had difficulty in reproducing leaf growth dynamics, with normalized root mean squared error (nRMSE) of 56% for leaf weight and 43% for LAI (across locations and models after Full calibration). Models with different levels of complexity and experience were capable of simulating final seed yield at maturity with reasonable accuracy (nRMSE of 8–31% after Full calibration). However, the nRMSE for pod weight (of 17–64% after Full calibration) was two-fold larger than that of seed yield. Moreover, the models differed in how they simulated timing from sowing to beginning seed growth (47–93 days) and effective seed filling period (18–54 days), owing to model structural differences in defining the reproductive developmental stages. Overall, we identified the following processes that can benefit from further model improvement: leaf expansion and senescence, reproductive phenology, and partitioning to reproductive growth. Simulation of pod wall tissue and individual seed cohorts is another aspect that many models currently lack. Model improvement can benefit from high-temporal resolution experimental datasets that concurrently account for phenology, plant growth, and partitioning. Further, we recommend collecting reproductive phenology in the field consistent with actual dry matter allocation to organs in the models and collecting multiple observations of seed and pod weight to aid model improvement for simulation of seed growth and yield formation.

Agricultural Model Intercomparison and Improvement↗

Calibration and evaluation of Skylab altimetry for geodetic determination of the geoid

The author has identified the following significant results. The analysis was based on a time series intrinsic relationship between the satellite ephemeris, altimeter measured ranges, and the corresponding a priori values of subsatellite geoidal heights. Using, least squares processing with parameter weighting, the objective was to recover: (1) the absolute geoidal heights of the subsatellite points; and (2) the associated altimeter calibration constants. Preliminary results from Skylab mission SL-2 are given, using various combinations from two sets of orbit ephemeris and altimeter ranges. It is shown that correctly scaled geoidal heights cannot be deduced by merely subtracting the altimeter range from the geodetic height of the satellite unless the satellite ephemeris and the altimeter have no unknown significant systematic errors or biases and drifts. It is emphasized that the primary objective of the Skylab altimeter is to determine the instrument feasibility. Any additional applications of the data such as for geodesy, geophysics, and oceanography are desirable. Although accurate orbit is required for such applications, it is not a prerequisite for determining the instrument feasibility.

Mourad, A. G.↗

Mitigation of Attentional Tunneling in the Flight Deck Using a Spatial Auditory Display

The role of attentional bias to a specific form of information, e.g., visual versus aural instructions, was evaluated to address the issue of attentional fixation (i.e., "cognitive tunneling") in the flight deck. A method of mitigating attentional tunneling using a spatial auditory display was evaluated in the presence and absence of a biasing factor. Data were gathered for participant recognition, accuracy and response times to aural and visual instructions under different conditions of distractor stimuli. Results showed a significant main effect for target acquisition with the use of spatial audio techniques. No significant main effects or interactions were found for target acquisition timing. Target acquisition was not affected by a "gaming score" biasing factor associated with the distractor tasks.

readback errors↗

Implementation and validation of optimal start control strategy for air conditioners and heat pumps

Commercial buildings are responsible for approximately 20 % of the total energy consumption and greenhouse gas emissions in the United States. Over 85 % of these buildings lack building automation systems, and many are small (<50,000 square feet), underserved, and use rooftop units (RTUs) for heating, ventilation, and air-conditioning needs. Because these buildings lack proper energy management systems, several operational deficiencies lead to excess energy consumption. Studies have shown that managing the RTUs’ heating and cooling set points, schedules, setbacks, and optimal start can result in a 20 % to 25 % reduction in electricity consumption in small commercial buildings. These buildings typically use fixed schedules to start the RTUs 60 to 120 min before occupancy begins, which results in excess energy consumption. This paper presents research that demonstrates and evaluates the performance of four optimal start methods, which utilize data-based modeling as a key element in facilitating adaptive control in response to time-varying inputs while requiring minimal sensor inputs. The evaluation found energy savings in two commercial buildings equipped with RTUs by periodically alternating four different optimal start models during the cooling and heating season. The resulting energy savings are positive for all models and range from 2 to 5 kWh/day/unit. The units on the east side of the building showed higher savings, while interior units showed greater variability in savings due to the differences in capacities and room sizes. Savings were considerably greater during the heating season compared to the cooling season. The performance of all four models on Mondays was poor; models suggested a shorter optimal start time, which resulted in relatively larger errors. Finally, the future work will look at using a different model for the days after weekends and holidays.

42 ENGINEERING↗

Measurement of the 230 Th( p ,2n)Pa229 and 230 Th( p ,3n)Pa228 reaction cross sections from 14.1 to 16.9 MeV

Actinium-225 is of interest for medical isotope production and there is on-going research into methods of producing Ac 225 , either directly or via the decay of its parent isotopes ( Th 229 , Pa 229 , and Ra 225 ). One method that has been suggested is the Th 230 ( p , 2 n ) Pa 229 reaction. However, there is no available cross-section data for this reaction in the literature. Purpose: Measure the Th 230 ( p , 2 n ) and Th 230 ( p , 3 n ) reaction cross sections in the energy range where the ( p , 2 n ) reaction is predicted to peak to determine the feasibility of Ac 225 production via the Th 230 ( p , 2 n ) reaction. Methods: Targets naturally enriched in Th 230 were irradiated at the Center for Accelerator Mass Spectrometry at Lawrence Livermore National Laboratory with energies ranging from 14.1 to 16.9 MeV. Furthermore, chemical processing was used to separate the protactinium activation products, followed by γ -ray spectroscopy to measure the activities of Pa 228 , 229 , 230 , 232 produced in the irradiation. Results: We find that excitation functions are reported for the first time in the literature for the Th 230 ( p , 2 n ) and Th 230 ( p , 3 n ) reactions in this energy range. The peak measured value of the Th 230 ( p , 2 n ) reaction was found to be 182 ± 12 mb at 14.4 ± 0.1 MeV. The Th 232 ( p , n ) Pa 232 reaction was used to verify the experimental conditions, the measured values are reported and are comparable to the existing literature values. From the γ -ray spectrometry data, the half-life of Pa 229 was measured as 1.5 ± 0.1 days, which is within the error of the half-life reported in the evaluated nuclear data as well as in the recent measurements, and the half-life of Pa 228 was measured as 19.5 ± 0.4 hours. Conclusions: Overall, the Th 230 ( p , 2 n ) Pa 229 reaction could reasonably be used for Ac 225 isotope production, although significant amounts of relatively isotopically pure Th 230 would be needed for significant production because the low alpha-decay branching ratio of Pa 229 and long half-life of Th 229 inhibit the in-growth of significant amounts of Ac 225 .

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

A Deep Learning Approach for In-Network Synchrophasor Missing Data Recovery Using Programmable Network Switches

Phasor measurement unit (PMU) networks deliver accurate and timely measurements, which is essential for managing today’s electric power systems. To ensure data quality and enhance the cyber-resilience of PMU networks against malicious attacks and data errors, this study presents an online PMU missing data recovery scheme by leveraging P4 programmable switches. The data plane incorporates a customized PMU protocol parser that abstracts the necessary payload data for recovery. Recovery processes are executed in the control plane using a pre-trained machine learning model. Both traditional and advanced ML models, such as transformer and TimeGPT, are explicitly employed for data prediction. This approach ensures rapid and precise data recovery. Performance evaluations focus on recovery speed and accuracy, using a real dataset from a campus microgrid. With 20% missing PMU data, the mean absolute percentage error for voltage magnitude is 0.0384%, and the phase angle error discrepancy is approximately 0.4064%.

Phasor Measurement Unit, Machine Learning, Program↗

A direct evaluation of the Geosat altimeter wet atmospheric range delay using very long baseline interferometry observations

The overall accuracy of the U.S. Navy Geosat altimeter wet atmospheric range delay caused by refraction through the atmosphere is directly assessed by comparing the estimates made from the DMSP Special Sensor Microwave/Imager and the U.S. Navy Fleet Numerical Ocean Center forecast model for Geosat with measurements of total zenith columnar water vapor content from four VLBI sites. The assessment is made by comparing time series of range delay from various methods at each location. To determine the importance of diurnal variation in water vapor content in noncoincident estimates, the VLBI measurements were made at 15-min intervals over a few days. The VLBI measurements showed strong diurnal variations in columnar water vapor at several sites, causing errors of the order 3 cm rms in any noncoincident measurement of the wet troposphere range delay. These errors have an effect on studies of annual and interannual changes in sea level with Geosat data.

Koblinsky, C. J.↗

Estimating Energy Market Schedules using Historical Price Data

The global climate crisis is expected to reshape the energy generation landscape in the coming decades. Increasing integration of non-dispatchable renewable energy resources into energy infrastructures and markets creates uncertainty as well as new opportunities for flexible energy systems. To conduct proper economic evaluation of flexible energy systems, such as integrated energy systems (IES), advancements in modelling of market interactions, such as bidding, is crucial. This work presents a shortcut algorithm which uses two mixed integer linear programs to compute dispatch schedules (e.g., hourly power production targets) that are constrained by the resource's bid information and characteristics (e.g., minimum up and down times) based on historical locational marginal price (LMP) data. The proposed algorithm is approximately 100 times faster and uses orders of magnitude less data than a full production cost model (PCM). We find the shortcut simulator recapitulates generator dispatch signals for the Prescient PCM with approximately 4% error for the RTS-GMLC test system.

electricity generation↗

The Grassmann interpolation method for spin-unrestricted open-shell systems

The recently reported Grassmann interpolation (G-Int) method [J. A. Tan and K. U. Lao, J. Chem. Phys. 158, 051101 (2023)] has been extended to spin-unrestricted open-shell systems. In contrast to closed-shell systems, where G-Int has to be performed only once since the α and β density matrices are the same, spin-unrestricted open-shell systems require G-Int to be performed twice—one for the α spin and another for the β spin density matrix. In this work, we tested the performance of G-Int to the carbon monoxide radical cation CO●+ and nickelocene complex, which have the doublet and triple ground states, respectively. We found that the Frobenius norm errors associated with the interpolations for the α and β spin density matrices are comparable for a given molecular geometry. These G-Int density matrices, when used as an initial guess for a self-consistent field (SCF) calculation, outperform the conventional SCF guess schemes, such as the superposition of atomic densities, purified superposition of atomic densities, core Hamiltonian, and generalized Wolfsberg–Helmholtz approximation. Depending on the desired accuracy, these G-Int density matrices can be used to directly evaluate the SCF energy without performing SCF iterations. In addition, the spin-unrestricted G-Int density matrices have been used for the first time to directly calculate the atomic charges using the Mulliken and ChElPG population analysis.

Chemistry↗