Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “input estimation”

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 613 records · Page 34

Extending SLUSCHI for Automated Diffusion Calculations

We present an extension of the SLUSCHI package (Solid and Liquid in Ultra Small Coexistence with Hovering Interfaces) to enable automated diffusion calculations from first-principles molecular dynamics. While the original SLUSCHI workflow was designed for melting temperature estimation via solid-liquid coexistence, we adapt its input and output handling to isolate the volume search stage and generate one production trajectory suitable for diffusion analysis. Post-processing tools parse VASP outputs, compute mean-square displacements (MSD), and extract tracer diffusivities using the Einstein relation with robust error estimates through block averaging. Diagnostic plots, including MSD curves, running slopes, and velocity autocorrelations, are produced automatically to help identify diffusive regimes. The method has been validated through representative case studies: self-diffusion in Al-Cu liquid alloys, sublattice melting in Li7La3Zr2O12 and Er2O3, interstitial oxygen transport in bcc and fcc Fe, and oxygen diffusivity in Fe-O liquids with variable Si and Al contents. Viscosity and diffusivity are linked through the Stokes-Einstein relation, with composition dependence assessed via simple linear mixing. This capability broadens SLUSCHI from melting-point predictions to transport property evaluation, enabling high-throughput, fully first-principles datasets of diffusion coefficients and viscosities across metals and oxides.

36 MATERIALS SCIENCE↗

Feasibility of Earthbound Motion in Lunar Gravity

BACKGROUND: Marginal stability of astronaut movement while performing lunar surface tasks has been well documented, and is clearly demonstrated in videos of falls, and near falls, during Apollo Lunar Extravehicular Activities. Referencing mission reports from Apollo 15 and 16 [1, 2], suspected causes for falls include: surface conditions, visibility, and gravitational effects (hypogravity). In this preliminary test, we employ the open-source biomechanical tool OpenSim [3, 4] to analyze the impact of lunar gravity (Lg) on two object-pickup motions performed by a single shirt-sleeved subject. Specifically, we attempt to answer the following questions based on an estimation of the Center of Mass Projection (CoMP) and ground reaction force Center of Pressure (COP) as it relates to the astronaut Base of Support (BOS) for 1g and Lg conditions: 1. Is the task motion, as performed in 1g, dynamically feasible in 1g and Lg? 2. Can we make the motion dynamically feasible in Lg by slowing it down? 3. Is the Lg COP equal to the 1g COP at a theoretically predicted reduction in motion speed? METHODS AND RESULTS: To answer the first question, the gravitational acceleration in the OpenSim model is modified from a nominal 9.81 m/s2 to 1.64 m/s2, and the 1g joint trajectory is input to an OpenSim based method [5] for estimating ground reaction forces and moments. From this method, the position of the CoMP and COP can be estimated and checked to see if they remain within a simulated BOS formed from the footprint of the OpenSim model to determine whether the motion is dynamically feasible. As expected, both of the motions were estimated to be feasible in a 1g environment, however, both motions had periods of infeasibility in Lg. It is well known that crew members make adjustments to motion trajectories in altered gravity fields to maintain balance. As a first step, we considered the simple adjustment of slowing the motion in Lg by a constant factor. This was accomplished by scaling the time stamps in the motion trajectory file by that factor. For the two motions considered, it was found that scale factors of 1.3 and 1.4 kept the COP just within the subject BOS. The CoMP is unchanged by the gravity field. Simple analysis of an inverted pendulum in the Lg environment, which generalizes to a general multibody system, leads to a theoretical prediction that a reduction in speed factor of √1g/Lg, or 2.445, will make the COP trajectory in Lg equivalent to that in 1g. When the above procedure was performed with a factor of 2.445, the estimated COP in Lg, was observed to be very close to that in 1g. In summary, we have developed a method for estimating the CoMP and COP in Lg, for subject motion collected in 1g. We believe this method can prove to be a valuable check and balance for simulated Lg training and testing by exposing potential simulator-induced artifacts that make the simulated task motion seem possible, when in fact, it would violate the above criteria. We also note that a reduction in task speed should tend the task motion towards stability, with a theoretical slowdown factor of √1g/Lg making the motion stability equal to that in 1g according to the CoMP and COP criteria.

R K Huffman↗

Probabilistic Assessment and Uncertainty Analysis of CO2 Storage Capacity of the Morrow B Sandstone—Farnsworth Field Unit

This paper presents probabilistic methods to estimate the quantity of carbon dioxide (CO2) that can be stored in a mature oil reservoir and analyzes the uncertainties associated with the estimation. This work uses data from the Farnsworth Field Unit (FWU), Ochiltree County, Texas, which is currently undergoing a tertiary recovery process. The input parameters are determined from seismic, core, and fluid analyses. The results of the estimation of the CO2 storage capacity of the reservoir are presented with both expectation curve and log probability plot. The expectation curve provides a range of possible outcomes such as the P90, P50, and P10. The deterministic value is calculated as the statistical mean of the storage capacity. The coefficient of variation and the uncertainty index, P10/P90, is used to analyze the overall uncertainty of the estimations. A relative impact plot is developed to analyze the sensitivity of the input parameters towards the total uncertainty and compared with Monte Carlo. In comparison to the Monte Carlo method, the results are practically the same. The probabilistic technique presented in this paper can be applied in different geological settings as well as other engineering applications.

03 NATURAL GAS↗

Twist-three cross-sections in deeply virtual Compton scattering

We study the deeply virtual Compton scattering process with both twist-two and twist-three Compton form factors and present our cross-sections formulas with all polarization configurations. While the twist-three contributions are generally assumed to be negligible in the literature due to the kinematical suppression, we compare them with the twist-two ones at typical JLab 6 GeV and 12 GeV kinematics as well as EIC kinematics and show their kinematical suppression explicitly, justifying the leading-twist approximation made in the literature. In addition, we also estimate the twist-three Compton form factors using Wandzura-Wilczek relations and inputs of twist-two generalized parton distributions based on a reggeized spectator model. With those estimated Compton form factors, we analyze the kinematical behavior of twist-two and twist-three cross-sections in a wide range of kinematics, and discuss the optimal regions for separating the leading-twist effects from the higher-twist ones.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Determination of stability and control parameters of a light airplane from flight data using two estimation methods

Two identification methods, the equation error method and the output error method, are used to estimate stability and control parameter values from flight data for a low-wing, single-engine, general aviation airplane. The estimated parameters from both methods are in very good agreement primarily because of sufficient accuracy of measured data. The estimated static parameters also agree with the results from steady flights. The effect of power different input forms are demonstrated. Examination of all results available gives the best values of estimated parameters and specifies their accuracies.

Klein, V.↗

Robust state estimation with imperfect sensors

Robustness properties of extended Kalman-type filters are investigated for time-varying plants when nonlinearities or linear dynamics are present at the plant input and output. Sufficient conditions for the asymptotic stability (of the state estimation error) in the absence of noise are obtained when the nonlinearities or the linear dynamics are perfectly known. When the nonlinearities are known only within bands of uncertainty, the estimation error is shown to be ultimately bounded under certain conditions. The results obtained highlight the robustness properties of such filters, which are similar to, but not exact duals of those of linear quadratic regulators.

Joshi, S. M.↗

Real-Time Frequency Response Estimation Using Joined-Wing SensorCraft Aeroelastic Wind-Tunnel Data

A new method is presented for estimating frequency responses and their uncertainties from wind-tunnel data in real time. The method uses orthogonal phase-optimized multi- sine excitation inputs and a recursive Fourier transform with a least-squares estimator. The method was first demonstrated with an F-16 nonlinear flight simulation and results showed that accurate short period frequency responses were obtained within 10 seconds. The method was then applied to wind-tunnel data from a previous aeroelastic test of the Joined- Wing SensorCraft. Frequency responses describing bending strains from simultaneous control surface excitations were estimated in a time-efficient manner.

Grauer, Jared A↗

Estimating Consequences of MMOD Penetrations on ISS

The threat from micrometeoroid and orbital debris (MMOD) impacts on space vehicles is often quantified in terms of the probability of no penetration (PNP). However, for large spacecraft, especially those with multiple compartments, a penetration may have a number of possible outcomes. The extent of the damage (diameter of hole, crack length or penetration depth), the location of the damage relative to critical equipment or crew, crew response, and even the time of day of the penetration are among the many factors that can affect the outcome. For the International Space Station (ISS), a Monte-Carlo style software code called Manned Spacecraft Crew Survivability (MSCSurv) is used to predict the probability of several outcomes of an MMOD penetration-broadly classified as loss of crew (LOC), crew evacuation (Evac), loss of escape vehicle (LEV), and nominal end of mission (NEOM). By generating large numbers of MMOD impacts (typically in the billions) and tracking the consequences, MSCSurv allows for the inclusion of a large number of parameters and models as well as enabling the consideration of uncertainties in the models and parameters. MSCSurv builds upon the results from NASA's Bumper software (which provides the probability of penetration and critical input data to MSCSurv) to allow analysts to estimate the probability of LOC, Evac, LEV, and NEOM. This paper briefly describes the overall methodology used by NASA to quantify LOC, Evac, LEV, and NEOM with particular emphasis on describing in broad terms how MSCSurv works and its capabilities and most significant models.

Evans, H.↗

Best estimate of the planetary boundary layer height from multiple remote sensing measurements

Remote sensing measurements have been widely used to estimate the planetary boundary layer height (PBLHT). Each remote sensing approach offers unique strengths and faces different limitations. In this study, we use machine learning (ML) methods to produce a best-estimate PBLHT (PBLHT-BE-ML) by integrating four PBLHT estimates derived from remote sensing measurements at the Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) observatory. Three ML models – random forest (RF) classifier, RF regressor, and light gradient-boosting machine (LightGBM) – were trained on a dataset from 2017 to 2023 that included radiosonde, various remote sensing PBLHT estimates, and atmospheric meteorological conditions. Evaluations indicated that PBLHT-BE-ML from all three models improved alignment with the PBLHT derived from radiosonde data (PBLHT-SONDE), with LightGBM demonstrating the highest accuracy under both stable and unstable boundary layer conditions. Feature analysis revealed that the most influential input features at the SGP site were the PBLHT estimates derived from (a) potential temperature profiles retrieved using Raman lidar (RL) and atmospheric emitted radiance interferometer (AERI) measurements (PBLHT-THERMO), (b) vertical velocity variance profiles from Doppler lidar (PBLHT-DL), and (c) aerosol backscatter profiles from micropulse lidar (PBLHT-MPL). The trained models were then used to predict PBLHT-BE-ML at a temporal resolution of 10 min, effectively capturing the diurnal evolution of PBLHT and its significant seasonal variations, with the largest diurnal variation observed over summer at the SGP site. We applied these trained models to data from the ARM Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE) field campaign (EPC), where the PBLHT-BE-ML, particularly with the LightGBM model, demonstrated improved accuracy against PBLHT-SONDE. Analyses of model performance at both the SGP and EPC sites suggest that expanding the training dataset to include various surface types, such as ocean and ice-covered areas, could further enhance ML model performance for PBLHT estimation across varied geographic regions.

Zhang, Damao [Pacific Northwest National Laborator↗

Simulation Evaluation of Pilot Inputs for Real Time Modeling During Commercial Flight Operations

Aircraft dynamics characteristics can only be identified from flight data when the aircraft dynamics are excited sufficiently. A preliminary study was conducted into what types and levels of manual piloted control excitation would be required for accurate Real-Time Parameter IDentification (RTPID) results by commercial airline pilots. This includes assessing the practicality for the pilot to provide this excitation when cued, and to further understand if pilot inputs during various phases of flight provide sufficient excitation naturally. An operationally representative task was evaluated by 5 commercial airline pilots using the NASA Ice Contamination Effects Flight Training Device (ICEFTD). Results showed that it is practical to use manual pilot inputs only as a means of achieving good RTPID in all phases of flight and in flight turbulence conditions. All pilots were effective in satisfying excitation requirements when cued. Much of the time, cueing was not even necessary, as just performing the required task provided enough excitation for accurate RTPID estimation. Pilot opinion surveys reported that the additional control inputs required when prompted by the excitation cueing were easy to make, quickly mastered, and required minimal training.

Martos, Borja↗

Methods for local gravity field approximation

The most widely known modern method for estimating gravity field values from observed data is least-squares collocation. Its advantages are that it can make estimates at arbitrary locations based on irregularly spaced observations, and that it makes use of statistical information about errors in the input data while providing corresponding information about the quality of the output estimates. Disadvantages of collocation include the necessity of inverting square matrices of dimension equal to the number of data values and the need to assume covariance models for the gravity field and the data errors. Fourier methods are an important alternative to collocation; having the advantage of greater computational efficiency, but requiring data estimates to be on a regular grid and not using or providing statistical accuracy information. The GEOFAST algorithm is an implementation of collocation that achieves high computational efficiency by transforming the estimation equations into the frequency domain where an accurate approximation may be made to reduce the workload. The forward and inverse Fast Fourier Transforms (FFTs) are utilized. The accuracy and computational efficiency of the GEOFAST algorithm is demonstrated using two sets of synthetic gravity data: marine gravity for an ocean trench region including wavelengths longer than 200 km; and local land gravity containing wavelengths as short as 5 km. These results are discussed along with issues such as the advantages of first removing reference field models before carrying out the estimation algorithm.

Sailor, R. V.↗

Updates to the Hanford Soil Inventory Model (SIM Version 2) for FY 2023

The purpose of this environmental calculation file (ECF) is to document the updates made to the Soil Inventory Model version 2.1 (SIM-v2.1) (ECF-HANFORD-21-0073, Updates to the Hanford Soil Inventory Model (SIM Version 2) for FY 2021) in FY 2023. These updates to SIM-v2.1 are referred to as the SIM-v2.2 version. The following updates have been made: • Include inventory estimate for a new analyte • hexavalent chromium (Cr(VI)), separate from total chromium inventory • Revise the inventory discharged to the 216-U-10 and 216-T-4 Pond systems based on partitioning of waste streams discharges over time and space among influent ditches and ponds • Enhance the preprocessing and postprocessing of the input and output files Note that the methodology for estimating Cr(VI) inventories is based on ECF-200W-23-0040, Recommendations for Updating Liquid Discharged Inventory and Transport Modeling Parameters for Cumulative Impacts Evaluation of Hexavalent Chromium in the 200 West Area while the revision of inventory discharged to the 216-U-10 and 216-T-4 Pond systems is based on ECF-HANFORD-19-0032, Distribution of Infiltration in the 216-U-10, 216-B-3 Pond, and 216-T-4 Pond Systems 1944-1997.

54 ENVIRONMENTAL SCIENCES↗

Probabilistic Modeling of a Three-Stage Human Landing System Architecture

Unmitigated uncertainties are known to have previously led to failed development programs; in order to combat these uncertainties, risks and their impacts must be understood and handled to ensure program success. In this paper, a probabilistic methodology to handle uncertainties is demonstrated on a three-element Human Landing System (HLS) concept, which allows tracking of current best estimates of the vehicle’s performance and assessment of its robustness against uncertainties. This methodology has two key parts: first, the creation of a dynamic architecture model of a three-element HLS concept; and second, its use with surrogate modeling and range estimating techniques to capture and propagate uncertainties. The DYnamic Rocket EQuation Tool (DYREQT), a space systems synthesis and sizing framework used by NASA, was used as to model the HLS architecture. For the probabilistic analysis, uncertainties of interest within the HLS concept were enumerated and represented as parameters within the DYREQT model as inputs for vehicle stages or mission profile events. Range estimating — a probabilistic method that combines Monte Carlo sampling, focus on critical parameters, and heuristics to assess risk and opportunities — is then adapted with operational parameters as well as vehicle parameters in the DYREQT model to capture mission uncertainty alongside vehicle uncertainty. To perform the range estimation portion of this methodology, the DYREQT model was sampled using a Design of Experiments (DoE) to efficiently explore the architecture design space with respect to the set of uncertainty parameters. Then, the results were used to create surrogate models, multivariate regressions that can visualize hypercube trends in the design space, of the architecture with respect to the uncertainty parameters. Using a correlation matrix constructed for the uncertainty parameters, previously independent samples were transformed to perform a Correlated Monte Carlo on the surrogate models. This probabilistic methodology was proved to provide insight into the underlying uncertainties of the three-element HLS architecture.

Stephanie Y Zhu↗

Probabilistic Modeling of a Three-Stage Human Landing System Architecture

Unmitigated uncertainties are known to have previously led to failed development programs; in order to combat these uncertainties, risks and their impacts must be understood and handled to ensure program success. In this paper, a probabilistic methodology to handle uncertainties is demonstrated on a three-element Human Landing System (HLS) concept, which allows tracking of current best estimates of the vehicle’s performance and assessment of its robustness against uncertainties. This methodology has two key parts: first, the creation of a dynamic architecture model of a three-element HLS concept; and second, its use with surrogate modeling and range estimating techniques to capture and propagate uncertainties. The DYnamic Rocket EQuation Tool (DYREQT), a space systems synthesis and sizing framework used by NASA, was used as to model the HLS architecture. For the probabilistic analysis, uncertainties of interest within the HLS concept were enumerated and represented as parameters within the DYREQT model as inputs for vehicle stages or mission profile events. Range estimating — a probabilistic method that combines Monte Carlo sampling, focus on critical parameters, and heuristics to assess risk and opportunities — is then adapted with operational parameters as well as vehicle parameters in the DYREQT model to capture mission uncertainty alongside vehicle uncertainty. To perform the range estimation portion of this methodology, the DYREQT model was sampled using a Design of Experiments (DoE) to efficiently explore the architecture design space with respect to the set of uncertainty parameters. Then, the results were used to create surrogate models, multivariate regressions that can visualize hypercube trends in the design space, of the architecture with respect to the uncertainty parameters. Using a correlation matrix constructed for the uncertainty parameters, previously independent samples were transformed to perform a Correlated Monte Carlo on the surrogate models. This probabilistic methodology was proved to provide insight into the underlying uncertainties of the three-element HLS architecture.

Stephanie Y. Zhu↗

U.S. IMERG Status

The U.S. Science Team’s Integrated Multi-satellitE Retrievals for the Global Precipitation Measurement (GPM) mission (IMERG) product provides estimates of surface precipitation rate and related information on a global 0.1° half-hour grid. It is run three times at increasingly longer latency and higher information input to serve different communities. This presentation will briefly summarize the major upgrades included in Version 07. The Goddard Profiling (GPROF) algorithm (applied to passive microwave sensor data from the GPM virtual constellation) and the Combined Radar-Radiometer Algorithm (CORRA) that provide input to IMERG V07 are both improved, and the new Precipitation Estimation from Remotely Sensed Information Using Artificial Neural Networks–Dynamic Infrared Rain Rate (PDIR) algorithm is being applied to infrared data, which employs an additional, warmer infrared temperature threshold than the previous infrared scheme and a deep neural network. V07 has improved bias performance as a result of calibrations that now employ the entire swath widths of CORRA and GPROF GPM Microwave Imager precipitation estimates. As well, time continuity in precipitation features is improved as a result of changes to the Kalman filter that approximately preserve the local histogram of precipitation rates (Scheme for Histogram Adjustment with Ranked Precipitation Estimates in the Neighborhood, or SHARPEN), and that better account for differences in sensor performance. A long-standing bug in the geolocation that shifted grid values 0.1° to the east in the latitude band 70°N-S has been corrected. Other changes in V07 include a hierarchical selection among motion vector sources to address deficiencies in the precipitation propagation near orography, an update to the precipitation phase specification for improved consistency with current inputs, and a renewed effort to implement climatological gauge adjustment to the near-real-time Runs. Early evaluations of V07 IMERG show interesting behavior across the TRMM orbit boost that is highly relevant to the planned orbit boost for the GPM Core observatory.

GPM↗

Evaluation of an Image-Derived Input Function for Kinetic Modeling of Nicotinic Acetylcholine Receptor-Binding PET Ligands in Mice

Positron emission tomography (PET) radioligands that bind with high-affinity to α4β2-type nicotinic receptors (α4β2Rs) allow for in vivo investigations of the mechanisms underlying nicotine addiction and smoking cessation. Here, we investigate the use of an image-derived arterial input function and the cerebellum for kinetic analysis of radioligand binding in mice. Two radioligands were explored: 2-[ 18 F]FA85380 (2-FA), displaying similar pKa and binding affinity to the smoking cessation drug varenicline (Chantix), and [ 18 F]Nifene, displaying similar pKa and binding affinity to nicotine. Time–activity curves of the left ventricle of the heart displayed similar distribution across wild type mice, mice lacking the β2-subunit for ligand binding, and acute nicotine-treated mice, whereas reference tissue binding displayed high variation between groups. Binding potential estimated from a two-tissue compartment model fit of the data with the image-derived input function were higher than estimates from reference tissue-based estimations. Rate constants of radioligand dissociation were very slow for 2-FA and very fast for Nifene. We conclude that using an image-derived input function for kinetic modeling of nicotinic PET ligands provides suitable results compared to reference tissue-based methods and that the chemical properties of 2-FA and Nifene are suitable to study receptor response to nicotine addiction and smoking cessation therapies.

2-FA85380↗

Evaluation of an Outer Loop Retrofit Architecture for Intelligent Turbofan Engine Thrust Control

The thrust control capability of a retrofit architecture for intelligent turbofan engine control and diagnostics is evaluated. The focus of the study is on the portion of the hierarchical architecture that performs thrust estimation and outer loop thrust control. The inner loop controls fan speed so the outer loop automatically adjusts the engine's fan speed command to maintain thrust at the desired level, based on pilot input, even as the engine deteriorates with use. The thrust estimation accuracy is assessed under nominal and deteriorated conditions at multiple operating points, and the closed loop thrust control performance is studied, all in a complex real-time nonlinear turbofan engine simulation test bed. The estimation capability, thrust response, and robustness to uncertainty in the form of engine degradation are evaluated.

Litt, Jonathan S.↗

Sensitivity of the Tropical Pacific Ocean to Precipitation Induced Freshwater Flux

We have performed a series of experiments using an ocean model to study the sensitivity of tropical Pacific Ocean to variations in precipitation induced freshwater fluxes. Variations in these fluxes arise from natural causes on all time scales. In addition, estimates of these fluxes are uncertain because of differences among measurement techniques. The model used is a quasi-isopycnal model, covering the Pacific from 40 S to 40 N. The surface forcing is constructed from observed wind stress, evaporation, precipitation, and surface temperature (SST) fields. The heat flux is produced with an iterative technique so as to maintain the model close to the observed climatology, but with only a weak damping to that climatology. Climatological estimates of evaporation are combined with various estimates of precipitation to determine the net surface freshwater flux. Results indicate that increased freshwater input decreases salinity as expected, but increases temperatures in the upper ocean. Using the freshwater flux estimated from the Microwave Sounding Unit leads to a warming of up to 0.6 C in the western Pacific over a case with zero net freshwater flux. SST is sensitive to the discrepancies among different precipitation observations, with root-mean-square differences in SST on the order of 0.2-0.3 C. The change in SST is more pronounced in the eastern Pacific, with differences of over 1 C found among the various precipitation products. Interannual variation in precipitation during El Nino events leads to increased warming. During the winter of 1982-83, freshwater flux accounts for about 0.4 C (approximately 10-15% of the maximum warming) of the surface warming in the central-eastern Pacific. Thus, the error of SST caused by the discrepancies in precipitation products is more than half of the SST anomaly produced by the interannual variability of observed precipitation. Further experiments, in which freshwater flux anomalies are imposed in the western, central, and eastern Pacific, show that the influence of net freshwater flux is also spatially dependent. The imposition of freshwater flux in the far western Pacific leads to a trapping of salinity anomaly to the surface layers near the equator. An identical flux imposed in the central Pacific produces deeper and off-equatorial salinity anomalies. The contrast between these two simulations is consistent with other simulations of the western Pacific barrier layer information.

Yang, Song↗