Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “empirical validation”

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 451 records · Page 25

A New Empirical Thermospheric Density Model JB2006 using New Solar Indices

A new empirical atmospheric density model is developed using the CIRA72 (Jacchia 71) model as the basis for the diffusion equations. New solar indices based on orbit based sensor data are used for the solar irradiances in the extreme and far ultraviolet wavelengths. New exospheric temperature and semiannual density equations are employed to represent the major thermospheric density variations. Temperature correction equations are also developed for diurnal and latitudinal effects, and finally density correction factors are used for model corrections required at high altitude (1500-4000 km). The new model, Jacchia-Bowman 2006, is validated through comparisons of accurate daily density drag data previously computed for numerous satellite. For 400 km altitude the standard deviation of 16% for the standard Jacchia model is reduced to 10% for the new JB2006 model for periods of low geomagnetic storm activity.

Bowman, Bruce R.↗

Estimating the Soil Temperature Profile from a Single Depth Observation: A Simple Empirical Heatflow Solution

Two data sets of experimental field observations with a range of meteorological conditions are used to investigate the possibility of modeling near-surface soil temperature profiles in a bare soil. It is shown that commonly used heat flow methods that assume a constant ground heat flux can not be used to model the extreme variations in temperature that occur near the surface. This paper proposes a simple approach for modeling the surface soil temperature profiles from a single depth observation. This approach consists of two parts: 1) modeling an instantaneous ground flux profile based on net radiation and the ground heat flux at 5cm depth; 2) using this ground heat flux profile to extrapolate a single temperature observation to a continuous near surface temperature profile. The new model is validated with an independent data set from a different soil and under a range of meteorological conditions.

Holmes, Thomas↗

Machine Learning for the Validation of Expert-Elicited Causal Risk Diagrams

Exposure to spaceflight poses risk to human health in complex ways. To help manage this risk, the Human Systems Risk Board (HSRB) at the National Aeronautics and Space Administration (NASA) maintains a set of causal diagrams that attempt to explain how spaceflight hazards generate health risks and lead to adverse outcomes both in-mission, immediately post-mission, and over the long term. These causal risk diagrams are formulated as directed acyclic graphs (DAGs) and can function as knowledge graphs of connected risks and outcomes. These DAGs have proven useful for communication, and, through network analysis, have allowed for the identification of structurally important factors in the risk network. However, the utility these DAGs provide is directly proportional to their verisimilitude, making assessment of this trait using empirical data – whether from actual human spaceflight or various spaceflight analogue exposures and model organisms – a high priority. In this research we explore the use of machine learning algorithms to learn DAG structure from empirical data as a means of evaluating human-elicited DAG structures. To do so, we test several different graph structure-learning algorithms on data concerning changes in the bones of rats and mice after exposure to either spaceflight or a spaceflight analogue. We explore potential methods for indexing the similarity between each algorithm’s output DAG with all the others and with that of the expert-elicited DAG. We discuss next steps in this ongoing line of research and open science initiatives underway to complete them.

directed acyclic graphs↗

Modeling the Quiet Time Outflow Solution in the Polar Cap

We use the Polar Wind Outflow Model (PWOM) to study the geomagnetically quiet conditions in the polar cap during solar maximum, The PWOM solves the gyrotropic transport equations for O(+), H(+), and He(+) along several magnetic field lines in the polar region in order to reconstruct the full 3D solution. We directly compare our simulation results to the data based empirical model of Kitamura et al. [2011] of electron density, which is based on 63 months of Akebono satellite observations. The modeled ion and electron temperatures are also compared with a statistical compilation of quiet time data obtained by the EISCAT Svalbard Radar (ESR) and Intercosmos Satellites (Kitamura et al. [2011]). The data and model agree reasonably well. This study shows that photoelectrons play an important role in explaining the differences between sunlit and dark results, ion composition, as well as ion and electron temperatures of the quiet time polar wind solution. Moreover, these results provide validation of the PWOM's ability to model the quiet time ((background" solution.

Glocer, Alex↗

A Change Impact Analysis to Characterize Evolving Program Behaviors

Change impact analysis techniques estimate the potential effects of changes made to software. Directed Incremental Symbolic Execution (DiSE) is an intraprocedural technique for characterizing the impact of software changes on program behaviors. DiSE first estimates the impact of the changes on the source code using program slicing techniques, and then uses the impact sets to guide symbolic execution to generate path conditions that characterize impacted program behaviors. DiSE, however, cannot reason about the flow of impact between methods and will fail to generate path conditions for certain impacted program behaviors. In this work, we present iDiSE, an extension to DiSE that performs an interprocedural analysis. iDiSE combines static and dynamic calling context information to efficiently generate impacted program behaviors across calling contexts. Information about impacted program behaviors is useful for testing, verification, and debugging of evolving programs. We present a case-study of our implementation of the iDiSE algorithm to demonstrate its efficiency at computing impacted program behaviors. Traditional notions of coverage are insufficient for characterizing the testing efforts used to validate evolving program behaviors because they do not take into account the impact of changes to the code. In this work we present novel definitions of impacted coverage metrics that are useful for evaluating the testing effort required to test evolving programs. We then describe how the notions of impacted coverage can be used to configure techniques such as DiSE and iDiSE in order to support regression testing related tasks. We also discuss how DiSE and iDiSE can be configured for debugging finding the root cause of errors introduced by changes made to the code. In our empirical evaluation we demonstrate that the configurations of DiSE and iDiSE can be used to support various software maintenance tasks

Rungta, Neha Shyam↗

Vertical Axis Rotational Motion Cues in Hovering Flight Simulation

A previous study that examined how yaw motion affected a pilot's ability to perform realistic hovering flight tasks indicated that any amount of pure yaw motion had little-to-no effect on pilot performance or opinion. In that experiment, pilots were located at the vehicle's center of rotation; thus lateral or longitudinal accelerations were absent. The purpose of the new study described here was to investigate further these unanticipated results for additional flight tasks, but with the introduction of linear accelerations associated with yaw rotations when the pilot is not at the center of rotation. The question of whether a yaw motion degree-of-freedom is necessary or not is important to government regulators who specify what simulator motions are necessary according to prescribed levels of simulator sophistication. Currently, specifies two levels of motion sophistication for flight simulators: full 6-degree-of-freedom and 3-degree-of-freedom. For the less sophisticated simulator, the assumed three degrees of freedom are pitch, roll, and heave. If other degrees of freedom are selected, which are different f rom these three, they must be qualified on a case-by-case basis. Picking the assumed three axes is reasonable and based upon experience, but little empirical data are available to support the selection of critical axes. Thus, the research described here is aimed at answering this question. The yaw and lateral degrees of freedom were selected to be examined first, and maneuvers were defined to uncouple these motions from changes in the gravity vector with respect to the pilot. This approach simplifies the problem to be examined. For this experiment, the NASA Ames Vertical Motion Simulator was used in a comprehensive investigation. The math model was an AH-64 Apache in hover, which was identified from flight test data and had previously been validated by several AH-64 pilots. The pilot's head was located 4.5 ft in front of the vehicle center of gravity, which is representative of the AH-64 pilot location. Six test pilots flew three tasks that were specifically designed to represent a broad class of situations in which both lateral and yaw motion cues may be useful. For the first task, the pilot controlled only the yaw axis and was required to rapidly acquire a North heading from 15 deg yaw offsets to either the East or West. This task allowed for full, or 1:1, motion to be used in all axes (yaw, lateral, and longitudinal). The second task was a 10 sec., 180 deg. pedal turn over a runway, but with the pilot only controlling the yaw degree-of-freedom. The position of the vehicle's center-of-mass remained fixed. This maneuver was taken from a current U.S. Army rotary wing design standard5 and is representative of a maneuver performed for acceptance of military helicopters; however, it does not allow for full 1:1 motion, since the simulator cab cannot rotate 180 deg. The third task required the pilot to perform a rapid 9 ft climb at a constant heading. This task was challenging, because rapid collective lever movement in the unaugmented AH64 results in a substantial yawing moment (due to engine torque) that must be countered by the pilot. This task also had full motion in all axes, but, in this case, the pilot had two axes to control simultaneously, rather than one as in the previous tasks. Four motion configurations were examined for each task: full motion (except for the 180 deg turn, for which the motion system was configured to provide as much motion as possible), full linear with no yaw motion, full yaw with no linear motion, and no motion. Each configuration was flown four times in a randomized test matrix, and the pilots were not informed of the configuration given. Vehicle state data were recorded for objective performance comparisons, and pilots provided subjective comments and ratings. As part of the pilots' evaluation, they were asked to rate the compensation required, the overall fidelity of the motion as compared to real flight, and whether motion was detected or not in each of the six degrees of freedom. In addition, the pilots provided a numerical level-of confidence rating, between 1 and 7, corresponding to how sure they were whether or not motion was present in each degree-of-freedom. The latter rating allow classical signal detection analysis to be performed.

Schroeder, Jeffrey A.↗

Improved Pyrotechnic Shock Scaling Method Based on Shock Response Spectrum and Statistical Energy Analysis

It has been more than forty years since the pyrotechnic shock scaling method was introduced. The scaling method estimates the attenuation of the Shock Response Spectrum (SRS) based on the distance from the source, structural configurations, types of structural joints and interfaces, and intervening structure [1]. The method has been successfully used in the spacecraft community countless times and still is frequently used to develop pyrotechnic shock requirements at various levels of assembly of a spacecraft. However, since the method was derived empirically from a limited set of shock test data [2], the aerospace community has been looking for an alternative approach. In this paper, a computational scaling method based on Statistical Energy Analysis (SEA) is re-visited. SEA is traditionally used as a method for investigating the diffusion of acoustic and vibratory energies of a system at steady-state [3]. Previously, an approach using the SRS as an acceleration constraint condition to the SEA model was introduced to estimate the attenuation from a shock source [4,5]. In the current investigation, additional examples are provided to further validate the approach. In addition, responses in the time domain are produced by the Local Modal Phase Reconstruction (LMPR) approach which was recently introduced to the community [6].

Fernandez, Juan P.↗

MISR Calibration Issues in High-Contrast Scenes, and Empirical Corrections

We diagnose the potential causes for the Multi-angle Imaging SpectroRadiometer's (MISR) persistent high aerosol optical depth (AOD) bias at low AOD with the aid of coincident MODerate-resolution Imaging Spectroradiometer (MODIS) imagery from NASA's Terra satellite. Internal reflections within the MISR instrument are responsible for a large portion of the high AOD bias in high-contrast scenes, which are especially common as broken-cloud situations over ocean. Discrepancies between MODIS and MISR nadir-viewing near-infrared (NIR) images are used to optimize nine parameters, along with a background reflectance modulation term (that was modeled separately), to represent the observed features. Independent, surface-based AOD measurements from the AErosol RObotic NETwork (AERONET) and the Marine Aerosol Network (MAN) are compared with MISR Research Algorithm (RA) AOD retrievals for 1118 coincidences to validate the corrections when applied to the nadir and off-nadir cameras. Additionally, the calibration coefficients for the red and NIR channels used for MISR over-water aerosol retrievals were reassessed with the RA to be consistent on a camera-by-camera basis. With these corrections, plus the baseline RA corrections applied (except enhanced cloud screening), the median AOD bias in the mid-visible (green) band decreases from 0.010 to 0.002, the RMSE decreases by approx. 10%, and the slope and correlation of the MISR vs. sun photometer Ångström Exponent improves. For AOD558 nm < 0.10 and with additional cloud screening, the median bias for the RA-retrieved AOD in the green band decreases from 0.011 to 0.003, compared to ~ 0.023 for the Standard Algorithm (SA). RMSE decreases by ~ 20% compared to the baseline (uncorrected) RA and by 17-53% compared to the SA. After all corrections and cloud screening are implemented, for AOD558 nm < 0.10, which includes about half the validation data, 68% absolute AOD errors for the RA have dropped to < 0.02 (approx. 0.018).

MISR↗

Evaluation of Cloud Microphysics Simulated using a Meso-Scale Model Coupled with a Spectral Bin Microphysical Scheme through Comparison with Observation Data by Ship-Borne Doppler and Space-Borne W-Band Radars

Equivalent radar reflectivity factors (Ze) measured by W-band radars are directly compared with the corresponding values calculated from a three-dimensional non-hydrostatic meso-scale model coupled with a spectral-bin-microphysical (SBM) scheme for cloud. Three case studies are the objects of this research: one targets a part of ship-borne observation using 95 GHz Doppler radar over the Pacific Ocean near Japan in May 2001; other two are aimed at two short segments of space-borne observation by the cloud profiling radar on CloudSat in November 2006. The numerical weather prediction (NWP) simulations reproduce general features of vertical structures of Ze and Doppler velocity. A main problem in the reproducibility is an overestimation of Ze in ice cloud layers. A frequency analysis shows a strong correlation between ice water contents (IWC) and Ze in the simulation; this characteristic is similar to those shown in prior on-site studies. From comparing with the empirical correlations by the prior studies, the simulated Ze is overestimated than the corresponding values in the studies at the same IWC. Whereas the comparison of Doppler velocities suggests that large-size snowflakes are necessary for producing large velocities under the freezing level and hence rules out the possibility that an overestimation of snow size causes the overestimation of Ze. Based on the results of several sensitivity tests, we conclude that the source of the overestimation is a bias in the microphysical calculation of Ze or an overestimation of IWC. To identify the source of the problems needs further validation research with other follow-up observations.

Iguchi, T.↗

The MODIS Aerosol Algorithm, Products, Validation and Applications

The MODerate resolution Imaging Spectroradiometer (MODIS) currently aboard both the Terra and Aqua satellites produces a suite of products designed to characterize global aerosol distribution, optical thickness and particle size. Never before has a space-borne instrument been able to provide such detailed information, complementing field and modeling efforts to produce a comprehensive picture of aerosol characteristics. The three years of Terra-MODIS data have been validated by comparing with co-located AERONET observations of aerosol optical thickness and derivations of aerosol size parameters. Some 8000 comparison points located at 133 AERONET sites around the globe show that the MODIS aerosol optical thickness retrievals are accurate to within the pre-launch expectations. MODIS-derived size parameters are also compared with AERONET retrievals and found to agree well for fine-mode dominated aerosol regimes. Aerosol regimes dominated by dust aerosol are less accurate, attributed to what is thought to be nonsphericity. Errors due to nonsphericity will be reduced by introducing a new set of empirical phase functions, derived without any assumptions of particle shape. The major innovation that MODIS bring to the field of remote sensing of aerosol is the measure of particle size and the separation of finemode and coarsemode dominated aerosol regimes. Particle size can separate finemode man-made aerosols created during combustion, from larger natural aerosols originating from salt spray or wind erosion. This separation allows for the calculation of aerosol radiative effect and the estimation of the man-made aerosol radiative forcing. MODIS can also be used in regional studies of aerosol-cloud interaction that affect the global radiative and hydrological cycles.

Remer, L. A.↗

Modeling Solar Zenith Angle Effects on the Polar Wind

We use the Polar Wind Outflow Model (PWOM) to study the geomagnetically quiet conditions in the polar cap during solar maximum. The PWOM solves the gyrotropic transport equations for O+, H+, and He+ along several magnetic field lines in the polar region in order to reconstruct the full 3D solution. We directly compare our simulation results to the data based empirical model of Kitamura et al. (2011) of electron density which is based on 63 months of Akebono satellite observations. The modeled ion and electron temperatures are also compared with a statistical compilation of quiet time data obtained by the EISCAT Svalbard Radar (ESR) and Intercosmos Satellites. The data and model agree reasonably well, albeit with some differences. This study shows that photoelectrons play an important role in explaining the differences between sunlit and dark results of electron density, ion composition, as well as ion and electron temperatures of the quiet time polar wind solution. Moreover, these results provide an initial validation of the PWOM s ability to model the quiet time "background" solution.

Glocer, A↗

Li-ion Battery Aging with Hybrid Physics-Informed Neural Networks and Fleet-wide Data

In this work, we propose a hybrid model for Li-ion battery discharge and aging prediction that leverages fleet-wide data to predict future capacity drops.The model is built upon an hybrid approach merging physics-based and empirical equations, as well as neural network models in a recurrent neural network cell. The hybrid physics-informed neural network can predict voltage discharge cycles given the loading profile, and estimate the used capacity of the battery under random-loading conditions by tracking aging parameters connected to the residual capacity of the battery. By merging information on the battery aging parameters with existing fleet-wide aging data, the model can predict the future residual capacity of the battery that is being monitored, and therefore enable predictions of voltage discharge curves far ahead in the battery life cycle. We validated the approach using the NASA Prognostics Data Repository Battery data-set, which contains experimental data on Li-ion batteries discharged at random loading conditions in a controlled environment. The approach also allows the identification of discrepancies between the battery aging trend and the trend observed at the fleet level, so that batteries behaving differently from the rest of the fleet can be subject to closer monitoring and further testing to refine predictions.

PINN↗

Convergence properties of finite-difference hydrodynamics schemes in the presence of shocks

We investigate the asymptotic convergence of finite-difference schemes for the Euler equations when the limiting solution contains shocks. The Lax-Wendroff theorem guarantees that certain conservative schemes converge to correct, physically valid solutions. We focus on two one-dimensional operator-split schemes with explicit artificial-viscosity terms. One, an internal-energy scheme, does not satisfy the assumptions of Lax-Wendroff; the other, a conservative total-energy scheme, does. With viscous lengths chosen proportional to the grid size, we find that both schemes converge to their zero-grid-size limits at the theoretically expected rate, but only the conversative scheme converges toward correct solutions of the inviscid fluid equations. We show that the difference in their behaviors results directly from the presence of shocks in the limiting solution. Empirically, we find that when the viscous lenghts tend toward zero more slowly than the grid size, however the nonconservative scheme also converges toward correct solutions. We characterize the asymptotic behavior of the total-energy scheme in a particular problem in which a shock forms. As the grid is refined, a Cauchy error approaches the expected rate of change slowly. We show that the changes in the artificial viscosity alter the diffusion of small-amplitude waves. The differences associated with such waves make the dominant contribution to the Cauchy error. We formulate an analytic model to relate the rate of approach to the effect of varying diffusion in waves and find quantitative agreement with our numerical results.

Kimoto, Paul A.↗

Chlorophyll-a Algorithms for Oligotrophic Oceans: A Novel Approach Based on Three-Band Reflectance Difference

A new empirical algorithm is proposed to estimate surface chlorophyll-a concentrations (Chl) in the global ocean for Chl less than or equal to 0.25 milligrams per cubic meters (approximately 77% of the global ocean area). The algorithm is based on a color index (CI), defined as the difference between remote sensing reflectance (R(sub rs), sr(sup -1) in the green and a reference formed linearly between R(sub rs) in the blue and red. For low Chl waters, in situ data showed a tighter (and therefore better) relationship between CI and Chl than between traditional band-ratios and Chl, which was further validated using global data collected concurrently by ship-borne and SeaWiFS satellite instruments. Model simulations showed that for low Chl waters, compared with the band-ratio algorithm, the CI-based algorithm (CIA) was more tolerant to changes in chlorophyll-specific backscattering coefficient, and performed similarly for different relative contributions of non-phytoplankton absorption. Simulations using existing atmospheric correction approaches further demonstrated that the CIA was much less sensitive than band-ratio algorithms to various errors induced by instrument noise and imperfect atmospheric correction (including sun glint and whitecap corrections). Image and time-series analyses of SeaWiFS and MODIS/Aqua data also showed improved performance in terms of reduced image noise, more coherent spatial and temporal patterns, and consistency between the two sensors. The reduction in noise and other errors is particularly useful to improve the detection of various ocean features such as eddies. Preliminary tests over MERIS and CZCS data indicate that the new approach should be generally applicable to all existing and future ocean color instruments.

Hu, Chuanmin↗

A Machine Learning Approach to Jet-Surface Interaction Noise Modeling

This paper investigates using machine learning to rapidly develop empirical models suitable for system-level aircraft noise studies. In particular, machine learning is used to train a neural network to predict the noise spectra produced by a round jet near a surface over a range of surface lengths, surface standoff distances, jet Mach numbers, and observer angles. These spectra include two sources, jet-mixing noise and jet-surface interaction (JSI) noise, with different scale factors as well as surface shielding and reflection effects to create a multi- dimensional problem. A second model is then trained using data from three rectangular nozzles to include nozzle aspect ratio in the spectral prediction. The training and validation data are from an extensive jet-surface interaction noise database acquired at the NASA Glenn Research Center's Aero-Acoustic Propulsion Laboratory. Although the number of training and validation points is small compared a typical machine learning application, the results of this investigation show that this approach is viable if the underlying data are well behaved.

Brown, Cliff↗

Hypervelocity impact on shielded plates

A ballistic limit equation for hypervelocity impact on thin plates is derived analytically. This equation applies to cases of impulsive impact on a plate that is protected by a multi-shock shield, and it is valid in the range of velocity above 6 km/s. Experimental tests were conducted at the NASA Johnson Space Center on square aluminum plates. Comparing the center deflections of these plates with the theoretical deflections of a rigid-plastic plate subjected to a blast load, one determines the dynamic yield strength of the plate material. The analysis is based on a theory for the expansion of the fragmented projectile and on a simple failure criterion. Curves are presented for the critical projectile radius versus the projectile velocity, and for the critical plate thickness versus the velocity. These curves are in good agreement with curves that have been generated empirically.

Smith, James P.↗

Modeling Magnetic Properties in EZTB

A software module that calculates magnetic properties of a semiconducting material has been written for incorporation into, and execution within, the Easy (Modular) Tight-Binding (EZTB) software infrastructure. [EZTB is designed to model the electronic structures of semiconductor devices ranging from bulk semiconductors, to quantum wells, quantum wires, and quantum dots. EZTB implements an empirical tight-binding mathematical model of the underlying physics.] This module can model the effect of a magnetic field applied along any direction and does not require any adjustment of model parameters. The module has thus far been applied to study the performances of silicon-based quantum computers in the presence of magnetic fields and of miscut angles in quantum wells. The module is expected to assist experimentalists in fabricating a spin qubit in a Si/SiGe quantum dot. This software can be executed in almost any Unix operating system, utilizes parallel computing, can be run as a Web-portal application program. The module has been validated by comparison of its predictions with experimental data available in the literature.

Lee, Seungwon↗

On the design of contractions and settling chambers for optimal turbulence manipulations in wind tunnels

The sensitivity of a moderate axisymmetric contraction of turbulence to upstream conditions is parametrically investigated in terms of turbulence intensities and integral length scales. Semi-empirical correlations are derived to characterize the development of turbulence intensities and integral length scales downstream of typical wind-tunnel turbulence manipulators. A new approach is described, permitting wind-tunnel designers to select the characteristic mesh of the turbulence manipulators, a well as the distance from them to the contraction, in order to obtain certain required turbulence characteristics in the test section at the exit of the contraction. The design charts strongly suggest that it may not be possible to achieve small scales and low intensities simultaneously, when using a contraction ratio of nine or larger. The results also confirm that the rapid distortion theory is only valid, with a proper viscous dissipation correction, up to the 'critical point' of a contraction. At that point, the scales of the various components of the velocity are equal and intercomponent transfer from the lateral to the streamwise velocities commences.

Nagib, H. M.↗