Search NASA⌕ Search

SEARCH · Search NASA

Results for “estimation steps”

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 433 records · Page 24

Far-infrared lines from H II regions: Abundance variations in the galaxy

Far-infrared lines of (N III) (57 microns), (O III) (52, 88 microns), (Ne III) (36 microns), and (S III) (19, 33 microns) have been measured in the H II regions G1.13 - 0.11, W31B, G23.95 + 0.15, G25.38 - 0.18, G29.96 - 0.02, W43, W51e, S156, S158, NGC 3576, NGC 3603, and G298.22-0.34. These observations were made with the facility Cryogenic Grating Spectrometer on the Kuiper Airborne Observatory to examine variations in abundances throughout the Galaxy. Previously published observations of G0.095 + 0.012, G333.60 - 0.21, G45.13 + 0.14A, K3-50, and M17 are also discussed. The giant H II region 30 Doradus in the Large Magellanic Cloud (LMC) was observed for comparison. Fluxes for (Ne II) (12.8 microns), (S IV) (10.5 microns), and the radio free-free continuum were collected from the literature for those sources. Electron densities were estimated from FIR line-pair ratios, and ionic abundances were estimated from the FIR line and radio fluxes. The excitation was estimated from the O(2+)/S(2+) ratio. Corrections for unseen ionization stages were calculated with the use of constnat-density H II region models. The validity and range of applicability of such semiempirical ionization correction schemes are discussed. The abundances with respect to hydrogen exhibit gradients with R(sub G) comparable to those previously measured for our Galaxy and for other galaxies. The overall gradients are d (log N/H)/dR = -0.10 +/- 0.02 dex/kpc, d(log Ne/H)/dR = -0.08 +/- 0.02 dex/kpc and d(log S/H)/dR = 0.07 +/- 0.02 dex/kpc. Compared to the Orion Nebula, the intermediate R(sub G) H II regions with 6 is less than R(sub G) is less than 11 kpc have similar or lower S/H and N/O ratios. The N/O ratios in the inner Galaxy are more than twice those observed in the Orion Nebula and intermediate R(sub G) H II regions. In fact, all the abundance ratios are as well or better fitted by a step fit with two levels than by a linear gradient. As has been noted in previous studies, the N/O ratio estimated from infrared observations of the doubly ionized N and O lines in H II regions is larger than the ratio estimated from optical observations of the singly ionized N and O lines. The Ne(2+)/O(2+) ratio is observed to be essentially constant over a wide range of excitation. This contradicts predictions of model H II regions calculated with the use of Local Thermodynamic Equilibrium (LTE) model stellar atmospheres. We conclude that these stellar atmospheres significantly underestimate the actual emergent fluxes for energies greater than 41 eV.

Simpson, Janet P.↗

Particle-based modelling of axisymmetric tandem mirror devices

In this work, we describe the use of a 1D-2V quasi-neutral hybrid electrostatic PIC with Monte-Carlo Coulomb collisions and non-uniform magnetic field to model the parallel transport and confinement in an axisymmetric tandem mirror device. End-plugs, based on simple-mirrors, are positioned at each end of the device and fueled with neutral beams (25 and 100 keV) to produce a sloshing ion population and increase the density of the end-plugs relative to the central cell. Results show the formation of a potential difference barrier between the central cell and the end-plugs. This potential confines a large fraction of the low energy thermal ions in the central cell which would otherwise be lost in a simple mirror, demonstrating the advantage of the beam-driven tandem mirror configuration relative to simple mirrors. In addition, we explore the effect of end-plug electron temperature on the confinement time of the device and compare it with theoretical estimates. Finally, we discuss the limitations of the code in its present form and describe the next logical steps to improve its predictive capability such as a fully nonlinear Fokker–Planck collision operator, multiply nested flux surface solutions and modeling the exhaust region up to the wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Exponentially Reduced Circuit Depths Using Trotter Error Mitigation

Product formulas are a popular class of digital quantum simulation algorithms due to their conceptual simplicity, low overhead, and performance, which often exceeds theoretical expectations. Recently, Richardson extrapolation and polynomial interpolation have been proposed to mitigate the Trotter error incurred by the use of these formulas. This work provides a rigorous, general analysis of these techniques for computing time-evolved observables, simplifying the interpolation algorithm in the process, and shows that extrapolation generically improves the performance of product formulas for this task. We demonstrate that, to achieve error 𝜖 in a simulation of time 𝑇 using a 𝑝 ⁢th-order product formula with extrapolation, circuit depths of 𝑂⁡(𝑇 1+1/𝑝 ⁢polylog (1/𝜖)) are sufficient—an exponential improvement in the precision over product formulas alone. Furthermore, we prove that these algorithms achieve commutator scaling, and improve the 𝑇 complexity for the interpolation algorithm. By relaxing the requirement of performing exact Chebyshev interpolation, our simplified algorithm eliminates the need for fractional implementations of Trotter steps, reducing computational overhead. Finally, we show these techniques can be combined with the classical shadows method to estimate many time-evolved local observables. Taken together, our findings provide the strongest evidence yet for the utility of Trotter error-mitigation techniques in algorithmic applications.

quantum algorithms & computation↗

S U ( 2 ) principal chiral model with tensor renormalization group on a cubic lattice

We study the continuous phase transition and thermodynamic observables in the three-dimensional Euclidean S U ( 2 ) principal chiral field model with the triad tensor renormalization group and the anisotropic tensor renormalization group methods. Using these methods, we find results that are consistent with previous Monte Carlo estimates and the predicted renormalization group scaling of the magnetization close to criticality. These results bring us one step closer to studying finite-density quantum chromodynamics in four dimensions using tensor network methods. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

AFUE Analysis Software Tool

The Annual Fuel Utilization Efficiency (AFUE) analysis of residential and light commercial furnaces follows ANSI/ASHRAE Standard 103-2017 (i.e., Method of Testing for Annual Fuel Utilization Efficiency of Residential Central Furnaces and Boilers) . The analysis is a complex, comprehensive method based on furnace configuration and specific components equipped, requiring detailed furnace testing and measurement data. For this reason, an AFUE Analysis Tool using Microsoft Excel enabled with Visual Basic for Applications (VBA) was developed with a user-friendly interface and comprehensive coverage. The tool consists of three worksheets: unit and configuration selection, geometry and measurement data input, and AFUE plus key results. This tool can be used to estimate the AFUE of both condensing and noncondensing furnaces with single-stage, two-stage, and step-modulating functions. The tool was validated with experimental data from Oak Ridge National Laboratory’s natural gas furnace projects that are commercially available. The results indicate the tool is reasonably accurate in the evaluation of a new R&D modified furnace unit.

Gao, Zhiming [Oak Ridge National Laboratory (ORNL)↗

Improving Self-Driving Labs: Quantifying System-Level Experiment Repeatability and Broadening Instrument-Level Compatibility

Modular Autonomous Research System (MARS) is a self-driving laboratory (SDL) which performs wet-lab science with peptide-lanthanide combinations in an automated and, ultimately, an autonomous manner to aid in soil analysis for domestic lithium mining. Autonomous experimentation involves automated experimentation, experiment planning, and active learning. MARS consists of a 6-axis robotic arm (UR5e) on a linear rail, pipette robots (Opentrons 2), and microplate readers. These components transport, operate on, and collect data with chemical solutions in standard labware. For effective autonomy, MARS must perform system-level labware operations repeatably, plan experiments autonomously, and be portable between research-domains. Repeatability is evaluated by labware placement precision, such that future operations can properly locate labware, as well as the elapsed time, so that low variance mean estimates of experiment duration can inform high-level researcher decision making. Autonomous experiment planning is the next step to decouple experimentation from human management; however, there is a conflict between the ideal system-level experiment goals and the constraints imposed by instruments’ limitations. Sub-domain portability is a long-term goal to extend MARS’ research beyond the chemistry of peptide-lanthanide binding to other sub-domains without having to invest significant overhead to system retrofitting. To address these goals, we manually trained the robotic arm labware placement and modelled statistical failurerate and uncertainty Additionally, we benchmarked the duration and variance of each experiment sub-operation as a heuristic for research decision making. Next, we use a parameterized geometric program (PGP) approach to design experiments that optimize system-level objectives and satisfy instrument-level constraints. Lastly, we proposed a Python framework to maximize MARS’ extensibility to other scientific sub-domains through a JSON-based experiment specification.

36 MATERIALS SCIENCE↗

Modeling human decision making behavior in supervisory control

An optimal decision control model was developed, which is based primarily on a dynamic programming algorithm which looks at all the available task possibilities, charts an optimal trajectory, and commits itself to do the first step (i.e., follow the optimal trajectory during the next time period), and then iterates the calculation. A Bayesian estimator was included which estimates the tasks which might occur in the immediate future and provides this information to the dynamic programming routine. Preliminary trials comparing the human subject's performance to that of the optimal model show a great similarity, but indicate that the human skips certain movements which require quick change in strategy.

Tulga, M. K.↗

Mechanism of the atmospheric oxidation of sulfur dioxide - Catalysis by hydroxyl radicals

A flash photolysis/resonance fluorescence technique was used to investigate the decay of OH due to the reaction OH + SO2 (+M) to HOSO2 (+M). In the presence of small amounts of NO (10 to the 14th per cu cm), the decays deviated from the normal semilogarithmic linearity due to reformation of OH. On the basis of computer simulations of the decay curves, it is suggested that the reactions HOSO2 + O2 to HO2 + SO3 (k3), and HO2 + HO to OH + NO2 are the likely subsequent steps in SO2 oxidation. The upper limit for the binding energy of HOSO2 relative to OH + SO2 is estimated to be 32 kcal/mol. The atmospheric implications of a catalytic oxidation mechanism are briefly discussed.

Margitan, J. J.↗

Computations of atmospheric angular momentum, with emphasis on the merit period

Three-dimensional wind and pressure fields used to calculate effective angular momentum functions (EAMFs) for earth rotation/polar motion are produced by a complex 3-step system at the U.S. National Meteorological Center. The resulting series of functions have errors which can be estimated by comparisons with similar series produced by the European Center for Medium Range Weather Forecasts. For the axial wind component, related to length of day (lod) changes, the small differences between the two series of atmospheric angular momentum (AAM) are confined to the tropics and the higher latitudes of the Southern Hemisphere. In addition, differences between the seasonal components of the AAM series and those of lod are found to be significantly reduced by including the entire stratosphere's momentum. Comparisons of the EAMFs for 8 calendar years indicate substantial interannual variability exists in the annual amplitudes of these forcing terms. In particular, the axial AAM showed quite different behavior from normal in 1983 and 1984, the years of the MERIT campaign. These two years showed very high and low annual amplitudes for AAM, respectively, which can be attributed to the effect of the tropical Pacific sea surface warming (El Nino) and its aftermath.

Salstein, David A.↗

A simple hydrologically based model of land surface water and energy fluxes for general circulation models

A generalization of the single soil layer variable infiltration capacity (VIC) land surface hydrological model previously implemented in the Geophysical Fluid Dynamics Laboratory (GFDL) general circulation model (GCM) is described. The new model is comprised of a two-layer characterization of the soil column, and uses an aerodynamic representation of the latent and sensible heat fluxes at the land surface. The infiltration algorithm for the upper layer is essentially the same as for the single layer VIC model, while the lower layer drainage formulation is of the form previously implemented in the Max-Planck-Institut GCM. The model partitions the area of interest (e.g., grid cell) into multiple land surface cover types; for each land cover type the fraction of roots in the upper and lower zone is specified. Evapotranspiration consists of three components: canopy evaporation, evaporation from bare soils, and transpiration, which is represented using a canopy and architectural resistance formulation. Once the latent heat flux has been computed, the surface energy balance is iterated to solve for the land surface temperature at each time step. The model was tested using long-term hydrologic and climatological data for Kings Creek, Kansas to estimate and validate the hydrological parameters, and surface flux data from three First International Satellite Land Surface Climatology Project Field Experiment (FIFE) intensive field campaigns in the summer-fall of 1987 to validate the surface energy fluxes.

Liang, XU↗

Creep of a Silicon Nitride Under Various Specimen/Loading Configurations

Extensive creep testing of a hot-pressed silicon nitride (NC132) was performed at 1300 C in air using five different specimen/loading configurations, including pure tension, pure compression, four-point uniaxial flexure, ball-on-ring biaxial flexure, and ring-on-ring biaxial flexure. Nominal creep strain and its rate for a given nominal applied stress were greatest in tension, least in compression, and intermediate in uniaxial and biaxial flexure. Except for the case of compressive loading, nominal creep strain generally decreased with time, resulting in less-defined steady-state condition. Of the four different creep formulations - power-law, hyperbolic sine, step, redistribution models - the conventional power-law model still provides the most convenient and reasonable means to estimate simple, quantitative creep parameters of the material. Predictions of creep deformation for the case of multiaxial stress state (biaxial flexure) were made based on pure tension and compression creep data by using the design code CARES/Creep.

Choi, Sung R.↗

Error Analysis for Discontinuous Galerkin Method for Parabolic Problems

In the proposal, the following three objectives are stated: (1) A p-version of the discontinuous Galerkin method for a one dimensional parabolic problem will be established. It should be recalled that the h-version in space was used for the discontinuous Galerkin method. An a priori error estimate as well as a posteriori estimate of this p-finite element discontinuous Galerkin method will be given. (2) The parameter alpha that describes the behavior double vertical line u(sub t)(t) double vertical line 2 was computed exactly. This was made feasible because of the explicitly specified initial condition. For practical heat transfer problems, the initial condition may have to be approximated. Also, if the parabolic problem is proposed on a multi-dimensional region, the parameter alpha, for most cases, would be difficult to compute exactly even in the case that the initial condition is known exactly. The second objective of this proposed research is to establish a method to estimate this parameter. This will be done by computing two discontinuous Galerkin approximate solutions at two different time steps starting from the initial time and use them to derive alpha. (3) The third objective is to consider the heat transfer problem over a two dimensional thin plate. The technique developed by Vogelius and Babuska will be used to establish a discontinuous Galerkin method in which the p-element will be used for through thickness approximation. This h-p finite element approach, that results in a dimensional reduction method, was used for elliptic problems, but the application appears new for the parabolic problem. The dimension reduction method will be discussed together with the time discretization method.

Kaneko, Hideaki↗

Practical Application of Model-based Programming and State-based Architecture to Space Missions

A viewgraph presentation to develop models from systems engineers that accomplish mission objectives and manage the health of the system is shown. The topics include: 1) Overview; 2) Motivation; 3) Objective/Vision; 4) Approach; 5) Background: The Mission Data System; 6) Background: State-based Control Architecture System; 7) Background: State Analysis; 8) Overview of State Analysis; 9) Background: MDS Software Frameworks; 10) Background: Model-based Programming; 10) Background: Titan Model-based Executive; 11) Model-based Execution Architecture; 12) Compatibility Analysis of MDS and Titan Architectures; 13) Integrating Model-based Programming and Execution into the Architecture; 14) State Analysis and Modeling; 15) IMU Subsystem State Effects Diagram; 16) Titan Subsystem Model: IMU Health; 17) Integrating Model-based Programming and Execution into the Software IMU; 18) Testing Program; 19) Computationally Tractable State Estimation & Fault Diagnosis; 20) Diagnostic Algorithm Performance; 21) Integration and Test Issues; 22) Demonstrated Benefits; and 23) Next Steps

Mission Data System (MDS)↗

LightForce Photon-Pressure Collision Avoidance: Updated Efficiency Analysis Utilizing a Highly Parallel Simulation Approach

This paper provides an updated efficiency analysis of the LightForce space debris collision avoidance scheme. LightForce aims to prevent collisions on warning by utilizing photon pressure from ground based, commercial off the shelf lasers. Past research has shown that a few ground-based systems consisting of 10 kilowatt class lasers directed by 1.5 meter telescopes with adaptive optics could lower the expected number of collisions in Low Earth Orbit (LEO) by an order of magnitude. Our simulation approach utilizes the entire Two Line Element (TLE) catalogue in LEO for a given day as initial input. Least-squares fitting of a TLE time series is used for an improved orbit estimate. We then calculate the probability of collision for all LEO objects in the catalogue for a time step of the simulation. The conjunctions that exceed a threshold probability of collision are then engaged by a simulated network of laser ground stations. After those engagements, the perturbed orbits are used to re-assess the probability of collision and evaluate the efficiency of the system. This paper describes new simulations with three updated aspects: 1) By utilizing a highly parallel simulation approach employing hundreds of processors, we have extended our analysis to a much broader dataset. The simulation time is extended to one year. 2) We analyze not only the efficiency of LightForce on conjunctions that naturally occur, but also take into account conjunctions caused by orbit perturbations due to LightForce engagements. 3) We use a new simulation approach that is regularly updating the LightForce engagement strategy, as it would be during actual operations. In this paper we present our simulation approach to parallelize the efficiency analysis, its computational performance and the resulting expected efficiency of the LightForce collision avoidance system. Results indicate that utilizing a network of four LightForce stations with 20 kilowatt lasers, 85% of all conjunctions with a probability of collision Pc > 10 (sup -6) can be mitigated.

conjunction analysis↗

BOLT II Forebody Surface Measurement Analysis

Surface impression molds of the nosetip/forebody interface, discrete tripping elements, and other surface imperfections were formed on the BOLT II flight vehicle prior to its launch in March 2022. The silicon-rubber molds cast on the interface between the nickel nosetip and the stainless steel forebody were measured using a 3D optical measurement system and a stylus profilometer and were evaluated for step direction and height. The step around the nosetip/forebody interface was found to vary in both magnitude and direction, but values were generally modest with most samples below 60 μm across the span and only two samples, both cast in the gutter, with step values larger than 100 μm. Impressions of the secondary-side tripping elements were strictly scanned with the optical device and similar methods were used to calculate trip height as step height. The measured centerline trip height was 2623 ±12 μm, while the leading-edge, quadrant 2 and 3, trip heights were 897 ±5 μm and 891 ±18 μm, respectively. Surface imperfections on the nose and forebody were also evaluated. Surface roughness values and spectra were estimated, and profiles were extracted across the imperfections to quantify the magnitude of those features.

Trip height↗

Adaptive state estimation for control of flexible structures

This paper proposes a new approach of obtaining adaptive state estimation of a system in the presence of unknown system disturbances and measurement noise. In the beginning, a non-optimal Kalman filter with arbitrary initial guess for the process and measurement noises is implemented. At the same time, an adaptive transversal predictor (ATP) based on the recursive least-squares (RLS) algorithm is used to yield optimal one- to p- step-ahead output predictions using the previous input/output data. Referring to these optimal predictions the Kalman filter gain is updated and the performance of the state estimation is thus improved. If forgetting factor is implemented in the recursive least-squares algorithm, this method is also capable of dealing with the situation when the noise statistics are slowly time-varying. This feature makes this new approach especially suitable for the control of flexible structures. A numerical example demonstrates the feasibility of this real time adaptive state estimation method.

Chen, Chung-Wen↗

Novel Estimation of Pilot Performance Characteristics

Two mechanisms internal to the pilot that affect performance during a tracking task are: 1) Pilot equalization (i.e. lead/lag); and 2) Pilot gain (i.e. sensitivity to the error signal). For some applications McRuer's Crossover Model can be used to anticipate what equalization will be employed to control a vehicle's dynamics. McRuer also established approximate time delays associated with different types of equalization - the more cognitive processing that is required due to equalization difficulty, the larger the time delay. However, the Crossover Model does not predict what the pilot gain will be. A nonlinear pilot control technique, observed and coined by the authors as 'amplitude clipping', is shown to improve stability, performance, and reduce workload when employed with vehicle dynamics that require high lead compensation by the pilot. Combining linear and nonlinear methods a novel approach is used to measure the pilot control parameters when amplitude clipping is present, allowing precise measurement in real time of key pilot control parameters. Based on the results of an experiment which was designed to probe workload primary drivers, a method is developed that estimates pilot spare capacity from readily observable measures and is tested for generality using multi-axis flight data. This paper documents the initial steps to developing a novel, simple objective metric for assessing pilot workload and its variation over time across a wide variety of tasks. Additionally, it offers a tangible, easily implementable methodology for anticipating a pilot's operating parameters and workload, and an effective design tool. The model shows promise in being able to precisely predict the actual pilot settings and workload, and observed tolerance of pilot parameter variation over the course of operation. Finally, an approach is proposed for generating Cooper-Harper ratings based on the workload and parameter estimation methodology.

Human pilot models↗

Real-Time Wave Energy Converter Control Using Instantaneous Frequency

Wave Energy Converters (WECs) rely on effective Power Take-Off (PTO) control strategies to maximize energy absorption under dynamic sea conditions. Traditional hydrodynamic modeling techniques may require computationally intensive convolution calculations, making real-time control implementation challenging. This paper presents an alternative approach by leveraging instantaneous frequency estimation to dynamically adjust PTO damping in response to varying wave frequencies. Two real-time frequency estimation methods are explored: the Hilbert Transform (HT) and Phase-Locked Loop (PLL). The Hilbert Transform method provides accurate frequency tracking but introduces a delayed response due to its dependence on causal data. Conversely, the PLL approach demonstrates strong potential in frequency tracking but requires careful gain tuning, particularly in complex sea states. Comparative evaluations across multiple test cases—including sinusoidal variations, amplitude steps, frequency step changes, and real-world JONSWAP spectrum waves—highlight the strengths and limitations of each method. The two different PTO control techniques across the various frequency estimation methods were tested under real-sea states using a state-space model of a point-absorbing Wave Energy Converter. The Capture Width Ratio (CWR) is used as a performance metric, with results showing that the HT achieves a 10.6% improvement, while the PLL estimation yields a 0.9% improvement relative to the fixed parameter control baseline. These results highlight the effectiveness of real-time frequency estimation in improving energy absorption compared to static control parameters.

WEC control↗