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 1,747 records · Page 97

Arctic ocean radiative fluxes and cloud forcing estimated from the ISCCP C2 cloud dataset, 1983-1990

Radiative fluxes and cloud forcings for the ocean areas of the Arctic are computed from the monthly cloud product of the International Satellite Cloud Climatology Project (ISCCP) for 1983-90. Spatially averaged short-wave fluxes are compared well with climatological values, while downwelling longwave fluxes are significantly lower. This is probably due to the fact that the ISCCP cloud amounts are underestimates. Top-of-the-atmosphere radiative fluxes are in excellent agreement with measurements from the Earth Radiation Budget Experiment (ERBE). Computed cloud forcings indicate that clouds have a warming effect at the surface and at the top of the atmosphere during winter and a cooling effect during summer. The net radiative effect of clouds is larger at the surface during winter but greater at the top of the atmosphere during summer. Overall the net radiative effect of clouds at the top of the atmosphere is one of cooling. This is in contrast to a previous result from ERBE data showing arctic cloud forcings have a net warming effect. Sensitivities to errors in input parameters are generally greater during winter with cloud amount being the most important paarameter. During summer the surface radiation balance is most sensitive to errors in the measurements of surface reflectance. The results are encouraging, but the estimated error of 20 W/sq m in surface net radiative fluxes is too large, given that estimates of the net radiative warming effect due to a doubling of CO2 are on the order of 4 W/sq m. Because it is difficult to determine the accuracy of results with existing in situ observations, it is recommended that the development of improved algorithms for the retrieval of surface radiative properties be accompanied by the simultaneous assembly of validation datasets.

Schweiger, Axel J.↗

Inverse Mapping of the Collision Kernel and Wall Flux Scaling in a Tall Convection‐Cloud Chamber Using Local Sensors and Knowledge‐Informed Deep Learning

Droplet collision–coalescence is a crucial process in cloud physics, but accurately representing this process under different dynamical conditions remains challenging. A proposed future convective‐cloud chamber aims to investigate this key process, but the method for observing it remains unclear, even though it is theoretically established that collision‐coalescence will occur. This study serves as a proof‐of‐concept demonstration of how knowledge‐informed deep learning, combined with measurement data from local sensors in the chamber, can be used to estimate the collision kernels, which determine how the droplet size distribution evolves during collision‐coalescence. In addition to estimating the collision kernel, we also address wall fluxes, another uncertain but important process that acts as a source of heat and moisture in the chamber. Ensemble runs of large‐eddy simulations are conducted by scaling the wall fluxes and the collision kernel, while the measured flow and cloud properties are used as inputs for a neural network. Results indicate that this approach successfully maps the scaling of wall fluxes and the collision kernel with biases of approximately 1% or less relative to the range of the target data. This proof‐of‐concept lays the groundwork for future applications; when the real measurements are available, real sensor data combined with the trained model presented in this work will enable estimation of the actual wall fluxes and collision kernel.

cloud chamber↗

Revised prediction (estimation) of Cape Kennedy, Florida, wind speed profile

The prediction of the wind profile maximum speed at Cape Kennedy, Florida, is made for any selected calendar data. The prediction is based on a normal probability distribution model with 15 years of smoothed input data and is static in the sense that no dynamic principles of persistence or synoptic features are considered. Comparison with similar predictions based on 6 years of data shows the same general pattern, but the variability decreased with the increase of sample size.

Guttman, N. B.↗

An Estimation of Turbulent Kinetic Energy and Energy Dissipation Rate Based on Atmospheric Boundary Layer Similarity Theory

Algorithms are developed to extract atmospheric boundary layer profiles for turbulence kinetic energy (TKE) and energy dissipation rate (EDR), with data from a meteorological tower as input. The profiles are based on similarity theory and scalings for the atmospheric boundary layer. The calculated profiles of EDR and TKE are required to match the observed values at 5 and 40 m. The algorithms are coded for operational use and yield plausible profiles over the diurnal variation of the atmospheric boundary layer.

Han, Jongil↗

Maximum likelihood tuning of a vehicle motion filter

This paper describes the use of maximum likelihood parameter estimation unknown parameters appearing in a nonlinear vehicle motion filter. The filter uses the kinematic equations of motion of a rigid body in motion over a spherical earth. The nine states of the filter represent vehicle velocity, attitude, and position. The inputs to the filter are three components of translational acceleration and three components of angular rate. Measurements used to update states include air data, altitude, position, and attitude. Expressions are derived for the elements of filter matrices needed to use air data in a body-fixed frame with filter states expressed in a geographic frame. An expression for the likelihood functions of the data is given, along with accurate approximations for the function's gradient and Hessian with respect to unknown parameters. These are used by a numerical quasi-Newton algorithm for maximizing the likelihood function of the data in order to estimate the unknown parameters. The parameter estimation algorithm is useful for processing data from aircraft flight tests or for tuning inertial navigation systems.

Trankle, Thomas L.↗

Validation of Spherically Symmetric Inversion by Use of a Tomographically Reconstructed Three-Dimensional Electron Density of the Solar Corona

Determining the coronal electron density by the inversion of white-light polarized brightness (pB) measurements by coronagraphs is a classic problem in solar physics. An inversion technique based on the spherically symmetric geometry (spherically symmetric inversion, SSI) was developed in the 1950s and has been widely applied to interpret various observations. However, to date there is no study of the uncertainty estimation of this method. We here present the detailed assessment of this method using a three-dimensional (3D) electron density in the corona from 1.5 to 4 solar radius as a model, which is reconstructed by a tomography method from STEREO/COR1 observations during the solar minimum in February 2008 (Carrington Rotation, CR 2066).We first show in theory and observation that the spherically symmetric polynomial approximation (SSPA) method and the Van de Hulst inversion technique are equivalent. Then we assess the SSPA method using synthesized pB images from the 3D density model, and find that the SSPA density values are close to the model inputs for the streamer core near the plane of the sky (POS) with differences generally smaller than about a factor of two; the former has the lower peak but extends more in both longitudinal and latitudinal directions than the latter. We estimate that the SSPA method may resolve the coronal density structure near the POS with angular resolution in longitude of about 50 deg. Our results confirm the suggestion that the SSI method is applicable to the solar minimum streamer (belt), as stated in some previous studies. In addition, we demonstrate that the SSPA method can be used to reconstruct the 3D coronal density, roughly in agreement with the reconstruction by tomography for a period of low solar activity (CR 2066). We suggest that the SSI method is complementary to the 3D tomographic technique in some cases, given that the development of the latter is still an ongoing research effort.

corona↗

The Microwave Radiative Properties of Falling Snow Derived from Nonspherical Ice Particle Models. Part II: Initial Testing Using Radar, Radiometer and In Situ Observations

In this study, two different particle models describing the structure and electromagnetic properties of snow are developed and evaluated for potential use in satellite combined radar-radiometer precipitation estimation algorithms. In the first model, snow particles are assumed to be homogeneous ice-air spheres with single-scattering properties derived from Mie theory. In the second model, snow particles are created by simulating the self-collection of pristine ice crystals into aggregate particles of different sizes, using different numbers and habits of the collected component crystals. Single-scattering properties of the resulting nonspherical snow particles are determined using the discrete dipole approximation. The size-distribution-integrated scattering properties of the spherical and nonspherical snow particles are incorporated into a dual-wavelength radar profiling algorithm that is applied to 14- and 34-GHz observations of stratiform precipitation from the ER-2 aircraft-borne High-Altitude Imaging Wind and Rain Airborne Profiler (HIWRAP) radar. The retrieved ice precipitation profiles are then input to a forward radiative transfer calculation in an attempt to simulate coincident radiance observations from the Conical Scanning Millimeter-Wave Imaging Radiometer (CoSMIR). Much greater consistency between the simulated and observed CoSMIR radiances is obtained using estimated profiles that are based upon the nonspherical crystal/aggregate snow particle model. Despite this greater consistency, there remain some discrepancies between the higher moments of the HIWRAP-retrieved precipitation size distributions and in situ distributions derived from microphysics probe observations obtained from Citation aircraft underflights of the ER-2. These discrepancies can only be eliminated if a subset of lower-density crystal/aggregate snow particles is assumed in the radar algorithm and in the interpretation of the in situ data.

Ice particles; Radiances; Radiative transfer; Airc↗

Robust time-domain identification of mass stiffness, and damping matrices

Accurate estimates of the mass, stiffness, and damping characteristics of a structure is necessary for determining the control laws best suited for active control methodologies. There are several modal identification techniques available for determining the frequencies, damping ratios, and mode shapes of a structure. However, modal identification methods in both the frequency and time domains have difficulties for certain circumstances. Frequency domain techniques which utilize the steady-state response from various harmonic inputs often encounter difficulties when the frequencies are closely distributed, the structure exhibits a high degree of damping, or the steady-state condition is hard to establish. Time domain techniques have produced successful results, but lack robustness with respect to measurement noise. In this paper, two identification techniques and an estimation method are combined to form a time-domain technique to accurately identify the mass, stiffness, and damping matrices from noisy measurements.

Roemer, Michael J.↗

Engineering description of the ascent/descent bet product

The Ascent/Descent output product is produced in the OPIP routine from three files which constitute its input. One of these, OPIP.IN, contains mission specific parameters. Meteorological data, such as atmospheric wind velocities, temperatures, and density, are obtained from the second file, the Corrected Meteorological Data File (METDATA). The third file is the TRJATTDATA file which contains the time-tagged state vectors that combine trajectory information from the Best Estimate of Trajectory (BET) filter, LBRET5, and Best Estimate of Attitude (BEA) derived from IMU telemetry. Each term in the two output data files (BETDATA and the Navigation Block, or NAVBLK) are defined. The description of the BETDATA file includes an outline of the algorithm used to calculate each term. To facilitate describing the algorithms, a nomenclature is defined. The description of the nomenclature includes a definition of the coordinate systems used. The NAVBLK file contains navigation input parameters. Each term in NAVBLK is defined and its source is listed. The production of NAVBLK requires only two computational algorithms. These two algorithms, which compute the terms DELTA and RSUBO, are described. Finally, the distribution of data in the NAVBLK records is listed.

Seacord, A. W., II↗

The Application of Differential Scanning Calorimetry to Investigate Precipitation Behavior in Nickel-Base Superalloys Under Continuous Cooling and Heating Conditions

In this study, a suite of experimental tools and fast-acting, numerical-simulation techniques was used to quantify the precipitation behavior of three nickel-base superalloys: IN-100, LSHR, and 718. Experimental methods comprised differential scanning calorimetry (DSC) to establish the specific heat as a function of temperature and selected direct-resistance heating trials (using a Gleeble ® machine) to obtain samples for microstructural analysis. For the DSC experiments, each alloy was cooled at a prescribed constant rate (between 5 and 20 K/min) after an initial soak/equilibration in the high-temperature, single-phase (supersolvus) temperature regime. On-heating DSC trials beginning at ambient temperature were also performed on alloy 718 in three different starting conditions: super- δ -solvus solution treated and water quenched (denoted as ST), solution treated and aged (STA), and solution treated and overaged (STOA). DSC results, revealing the thermal signatures associated with the kinetics of precipitation of γ ' (IN-100, LSHR) or γ ' and γ " (718), were interpreted using a previously-developed fast-acting routine that treats concurrent nucleation, growth, coarsening, and dissolution. For these simulations, special attention was paid to various thermo-kinetic input parameters including equilibrium solvus-approach curves, bulk free energies of transformation, matrix-precipitate interface energies, and effective diffusivities. For the γ - γ ' superalloys (IN-100 and LSHR), estimates of precipitate volume fraction as a function of temperature from the specific-heat data revealed semi-quantitative agreement with simulation predictions. For the γ - γ '- γ " superalloy (718), simulation predictions of precipitate volume fractions were converted to specific heat as a function of temperature and showed semi-quantitative agreement with the direct measurements.

36 MATERIALS SCIENCE↗

Assessing the Effect of a Deep‐Rooted Grass on Belowground Carbon Storage in Cultivated Land: Insights From a Multi‐Site US Study

Agriculture depletes soil organic carbon (SOC), partly due to the exclusion of deep-rooted perennials. Reintroducing deep-rooted perennials to cultivated land may help to mitigate SOC loss. We quantified the effect of deep roots on SOC by comparing 8 to 30 year-old stands of switchgrass ( Panicum virgatum L .) with paired annual row crop fields at 12 sites across the central and eastern USA. We hypothesized that switchgrass would store more root C and SOC than neighboring shallow-rooted annual crops, and that these effects would extend deeper than 30 cm. We also evaluated whether switchgrass stimulates decomposition of SOC at depth using radiocarbon ( 14 C). Finally, we explored whether the effect of switchgrass on SOC is moderated by soil chemical and physical properties. While the effect of switchgrass on SOC in the surface 100 cm was positive at most sites, the average effect was not statistically significant (difference in SOC = 0.6 kg C m −2 [95% CI −0.8 to +1.9 kg C m −2 ]). By contrast, we found that root C was consistently more abundant under switchgrass, yielding an estimated additional 0.6 kg C m −2 in the surface 100 cm of soil [95% CI +0.5 to +0.7 kg C m −2 ]. 14 C measurements suggested that root C inputs were adding to existing SOC without stimulating decomposition. The effect of switchgrass on belowground C was not strongly related to any of the soil properties that we evaluated. Our observations show that root C can contribute substantially to belowground C stocks when deep-rooted perennials replace shallow-rooted crops.

60 APPLIED LIFE SCIENCES↗

Dealing with Measurement Uncertainties as Nuisance Parameters in Bayesian Model Calibration

In the presence of model discrepancy, the calibration of physics-based models for physical parameter inference is a challenging problem. Lack of identifiability between calibration parameters and model discrepancy requires additional identifiability constraints to be placed on the model discrepancy to obtain unique physical parameter estimates. If these assumptions are violated, the inference for the calibration parameters can be systematically biased. In many applications, such as in dynamic material property experiments, many of the calibration inputs refer to measurement uncertainties. Here, we develop a metric for identifying overfitting of these measurement uncertainties, propose a prior capable of reducing this overfitting, and show how this leads to a diagnostic tool for validation of physical parameter inference. The approach is demonstrated for a benchmark example and applied for a material property application to perform inference on the equation of state parameters of tantalum.

97 MATHEMATICS AND COMPUTING↗

Aerosol and Cloud Optical Properties from the ARM Raman Lidars: The Feature Detection and Extinction (RLPROF-FEX) Value-Added Product

Aerosols and their interactions and influence on clouds are among the main sources of uncertainties in radiative direct and indirect forcing (IPCC 2013). Continuous height-resolved measurements of cloud and aerosol optical properties are needed to reduce these uncertainties. Here we describe the Raman Lidar Profiles – Feature detection and Extinction (RLPROF-FEX) Value-Added Product (VAP) derived using Raman lidar data at multiple U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) user facility sites. RLPROF-FEX provides estimates of extinction, backscatter, and depolarization using the algorithm described by Thorsen et al. 2015 and Thorsen and Fu 2015. This document provides a description of the FEX algorithm, its input and output data, and related details about the Raman lidar (RL) system.

54 ENVIRONMENTAL SCIENCES↗

Release of Evaluated 235 U(n,f) Average Prompt Fission Neutron Multiplicities Including the CGMF Model

This report documents an evaluation of the average prompt fission neutron multiplicity, $\overline{v}_p$, of 235 U from 200 keV to 15 MeV that is a potential release candidate for the upcoming U.S. nuclear data library, ENDF/B-VIII.1. This evaluation had to be re-done from "scratch", as the input to the $\overline{v}_p$ evaluation of the previous library, ENDF/B-VIII.0, was lost. That means that all available experimental data were re-analyzed and uncertainties were re-estimated. Another major difference to ENDF/B-VIII.0 is that this evaluation includes model information from the Hauser-Feshbach fission fragment decay code CGMF, while ENDF/B-VIII.0 is based purely on experimental data. CGMF links several fission quantities with each other; $\overline{v}_p$ is predicted by assumptions made on, e.g., pre-neutron emission yields as a function of mass, the total kinetic energy, or spin and parity of fission fragments. This allows to perform two types of validation for the new 235 U $\overline{v}_p$: On the one hand, one can employ evaluated CGMF parameters obtained from fitting to experimental 235 U $\overline{v}_p$ to predict yields as a function of mass, the average total kinetic energy, or the mean energy of the prompt fission neutron spectrum. These model-predicted values can then be compared to experimental and evaluated data. The model-predicted fission-observable values using evaluated parameters obtained here are reasonably close to experimental data indicating the evaluated 235 U(n,f) $\overline{v}_p$ are physical. On the other hand, one can validate 235 U $\overline{v}_p$ with respect to integral responses such as fast ICSBEP critical assemblies or LLNL pulsed spheres. LLNL pulsed-sphere neutron-leakage spectra are minimally impacted by the new 235 U $\overline{v}_p$ as these experimental data are shape data and the $\overline{v}_p$ would mostly lead to a change in normalization of the data as the spheres are relatively thin (0.7 and 1.5 mean-free path) and, thus, mostly depend on 235 U $\overline{v}_p$ from 12-15 MeV. The change in the predicted effective neutron multiplication factor, k eff , of selected ICSBEP critical assemblies, however, is large compared to values using ENDF/B-VIII.0 and experimental k eff : The average bias is 108 pcm across all studied k eff values versus 12 pcm for ENDF/B-VIII.0. A reasonable performance in simulating keff (mean bias of 14 pcm) can be retained by tweaking 235 U $\overline{v}_p$ from 3-5 MeV, and combining it with a recent 235 U PFNS evaluation that is also a ENDF/B-VIII.1 release candidate.

235U↗

Deep Image Prior Enabled Full Waveform Inversion (Final Technical Report)

MS Student Naveen Gupta worked on the problem of full waveform inversion (FWI) using neural networks as shown in Figure 1. Our goal was to learn a neural network to represent the subsurface velocity model, which when fed into the FWI module (implemented using a numerical forward model of wave equations) produces amplitude estimates that match with ground-truth observations of amplitude. We used neural networks to solve the inverse problem of estimating velocity distributions for a given seismic amplitude data such that, once trained, our neural network model can generate a distribution of velocity profiles for different random vectors fed as inputs to the neural network model.

97 MATHEMATICS AND COMPUTING↗

Response of winter and spring wheat grain yields to meteorological variation

Mathematical models which quantify the relation of wheat yield to selected weather-related variables are presented. Other sources of variation (amount of applied nitrogen, improved varieties, cultural practices) have been incorporated in the models to explain yield variation both singly and in combination with weather-related variables. Separate models were developed for fall-planted (winter) and spring-planted (spring) wheats. Meteorological variation is observed, basically, by daily measurements of minimum and maximum temperatures, precipitation, and tabled values of solar radiation at the edge of the atmosphere and daylength. Two different soil moisture budgets are suggested to compute simulated values of evapotranspiration; one uses the above-mentioned inputs, the other uses the measured temperatures and precipitation but replaces the tabled values (solar radiation and daylength) by measured solar radiation and satellite-derived multispectral scanner data to estimate leaf area index. Weather-related variables are defined by phenological stages, rather than calendar periods, to make the models more universally applicable.

Feyerherm, A. M.↗

Large-scale wind-tunnel investigation of a close-coupled canard-delta-wing fighter model through high angles of attack

Data obtained through very high angles of attack from a large-scale, subsonic wind-tunnel test of a close-coupled canard-delta-wing fighter model are analyzed. The canard delays wing leading-edge vortex breakdown, even for angles of attack at which the canard is completely stalled. A vortex-lattice method was applied which gave good predictions of lift and pitching moment up to an angle of attack of about 20 deg, where vortex-breakdown effects on performance become significant. Pitch-control inputs generally retain full effectiveness up to the angle of attack of maximum lift, beyond which, effectiveness drops off rapidly. A high-angle-of-attack prediction method gives good estimates of lift and drag for the completely stalled aircraft. Roll asymmetry observed at zero sideslip is apparently caused by an asymmetry in the model support structure.

Stoll, F.↗

Bearing detection in the presence of two sources of varying coherence using the complex cepstrum

The effect of the presence of two acoustic sources (one, the primary, whose location is to be detected) of varying coherence on a cepstral bearing finding procedure is experimentally studied. The coherence between the acoustic sources was altered by adding random noise of various SNR (signal-to-noise ratio) to the input signal of the primary source; the same base signal being fed to both sources. The results demonstrate that, when block liftering is used, the primary source bearing is reliably estimated for coherences as low as gamma sup 2 greater than or approx equal to 0.5. The results also imply that background noise (unreflected) of SNR greater than or approx equal to 10 dB will not markedly affect the accuracy of the bearing estimation algorithm.

Fuller, C. R.↗