Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical error”

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,207 records · Page 67

Optimal application of Morrison's iterative noise removal for deconvolution. Appendices

Morrison's iterative method of noise removal, or Morrison's smoothing, is applied in a simulation to noise-added data sets of various noise levels to determine its optimum use. Morrison's smoothing is applied for noise removal alone, and for noise removal prior to deconvolution. For the latter, an accurate method is analyzed to provide confidence in the optimization. The method consists of convolving the data with an inverse filter calculated by taking the inverse discrete Fourier transform of the reciprocal of the transform of the response of the system. Various length filters are calculated for the narrow and wide Gaussian response functions used. Deconvolution of non-noisy data is performed, and the error in each deconvolution calculated. Plots are produced of error versus filter length; and from these plots the most accurate length filters determined. The statistical methodologies employed in the optimizations of Morrison's method are similar. A typical peak-type input is selected and convolved with the two response functions to produce the data sets to be analyzed. Both constant and ordinate-dependent Gaussian distributed noise is added to the data, where the noise levels of the data are characterized by their signal-to-noise ratios. The error measures employed in the optimizations are the L1 and L2 norms. Results of the optimizations for both Gaussians, both noise types, and both norms include figures of optimum iteration number and error improvement versus signal-to-noise ratio, and tables of results. The statistical variation of all quantities considered is also given.

Ioup, George E.↗

Design and Modeling of a Coded-Interleaving System in the Presence of Fading

Future deep space links are migrating towards higher frequency bands such as Ka band and optical. These links are susceptible to non linear effects (e.g. cloud coverage, atmospheric turbulence, scintillation, antenna mis-pointing, jitter, relative motion, etc.) that cause various degrees of fading loss over various time scales. We consider a channel experiencing block fading where the channel coherence time is a lot longer than the length of a codeword. As long as the fading duration is not excessive, an interleaver can be used in a coded-communication system that randomizes the fades among symbols of multiple codewords. This achieves an “averaging” effect that spreads the signal energy more evenly among the codewords, and improves the likelihood that each of them can be decoded correctly. For simplicity, we consider a block interleaver where adjacent symbols of a codeword have a fixed delay (or separation) with each other, and investigate the error-rate performance of a coded-interleaving system (with fixed and finite interleaving depth) in the presence of block fading. Unlike prior research which assumes theoretical channel coherence models, we aim to build upon our analysis using empirical spacecraft signal measurements obtained with the Deep Space Network, and to use the fundamental theories of probability and statistics (e.g. variations of Law of Large Numbers) to quantify the error-rate performance of the coded-interleaving system. We propose to utilize these results in the design of a two-dimensional coded-interleaving system for fading channels.

Sanchez-Net, Marc↗

Statistically informed upscaling of damage evolution in brittle materials

The presence and growth of micro-cracks degrade the strength of brittle solids, greatly impacting the overall material response. Hence, the evolution of these micro-cracks must be accounted for in models describing the relationship between stress and strain so that accurate predictions of material failure can be made. The evolution of individual cracks and crack networks can be simulated with high-fidelity microscale models utilizing highly resolved meshes that can be computationally expensive to a point in which it limits the simulation scale. Hence, for many engineering applications that require simulations of large components, continuum-scale models, which cannot explicitly resolve individual cracks and thus lose important physical information, are required. In this work, we bridge these two scales by developing and implementing a continuum-scale effective moduli constitutive model that is informed by crack statistics generated from a high-fidelity model resolved using a finite-discrete element method (FDEM) implementation. Using statistical information describing the evolution of crack lengths and orientations, this model can capture the effects of brittle damage evolution without the need to resolve individual cracks. We have successfully captured the stress-strain behavior of the high-fidelity simulations using the statistics-based constitutive model shown through direct comparison of stress-strain curves. The curves match within error bars present in the strain-softening portions of the stress-strain curve of the high-fidelity results due to the statistical variation of the initial pre-existing crack network. In conclusion, the stress-strain curves are also compared to experimental results for similar loading conditions and show good qualitative agreement.

36 MATERIALS SCIENCE↗

Droplet size distributions in turbulent clouds: experimental evaluation of theoretical distributions

Precipitation efficiency and optical properties of clouds, both central to determining Earth's weather and climate, depend on the size distribution of cloud particles. Here in this paper, theoretical expressions for cloud droplet size distribution shape are evaluated using measurements from controlled experiments in a convective-cloud chamber. The experiments are a unique opportunity to constrain theory because they are in steady-state and because the initial and boundary conditions are well characterized compared to typical atmospheric measurements. Three theoretical distributions obtained from a Langevin drift-diffusion approach to cloud formation via stochastic condensation are tested: (a) stochastic condensation with a constant removal time-scale; (b) stochastic condensation with a size-dependent removal time-scale; (c) droplet growth in a fixed supersaturation condition and with size-dependent removal. In addition, a similar Weibull distribution that can be obtained from the drift-diffusion approach, as well as from mechanism-independent probabilistic arguments (e.g., maximum entropy), is tested as a fourth hypothesis. Statistical techniques such as the χ2 test, sum of squared errors of prediction, and residual analysis are employed to judge relative success or failure of the theoretical distributions to describe the experimental data. An extensive set of cloud droplet size distributions are measured under different aerosol injection rates. Five different aerosol injection rates are run both for size-selected aerosol particles, and six aerosol injection rates are run for broad-distribution, polydisperse aerosol particles. In relative comparison, the most favourable comparison to the measurements is the expression for stochastic condensation with size-dependent droplet removal rate. However, even this optimal distribution breaks down for broad aerosol size distributions, primarily due to deviations from the measured large-droplet tail. Lastly, a possible explanation for the deviation is the Ostwald ripening effect coupled with deactivation/activation in polluted cloud conditions.

54 ENVIRONMENTAL SCIENCES↗

Melt Pool and Heat Treatment Optimization for the Fabrication of High-Strength and High-Toughness Additively Manufactured 4340 Steel

Additively manufactured (AM) components offer superior design flexibility compared to their conventionally manufactured counterparts, and optimizing processing parameters is key to achieving high-quality depositions with desirable and predictable mechanical properties. This study was focused on 4340 steel fabricated using laser powder bed fusion (LPBF), and 42 laser power and scan speed combinations have been systematically investigated to determine an optimized melt pool geometry that would ensure fully-dense parts. The AM material was compared with a wrought 4340 equivalent and studied in two customized heat treated conditions, optimized for strength and toughness, respectively. The microstructures of the as-fabricated and heat treated AM and wrought materials were characterized to assess differences introduced by the layer-by-layer fabrication process and subsequent heat treatment. Tensile properties of both materials were also evaluated and demonstrate that the AM materials offer equal or superior properties compared to the wrought equivalents. Differences in fracture surface morphologies indicate the distinct failure mechanisms associated with the materials’ characteristic microstructures, and the role of inclusions in the failures was studied to elucidate these differences. Complementary to the experimental investigations, the dataset was leveraged to make recommendations for future design of experiments to optimize AM build parameters in other material systems. A statistical Monte Carlo analysis was used to predict the interpolation error produced using reduced datasets and to enable informed processing parameters selection. These findings are discussed to make recommendations for the use of AM materials for high-integrity structural applications.

36 MATERIALS SCIENCE↗

Resolving Mesoscale Convective Systems: Grid Spacing Sensitivity in the Tropics and Midlatitudes

Abstract Mesoscale convective systems (MCSs) are a critical global water cycle component and drive extreme precipitation events in tropical and midlatitude regions. However, simulating deep convection remains challenging for modern numerical weather and climate models due to the complex interactions of processes from microscales to synoptic scales. Recent models with kilometer‐scale horizontal grid spacings offer notable improvements in simulating deep convection compared to coarser‐resolution models. Still, deficiencies in representing key physical processes, such as entrainment, lead to systematic biases. Additionally, evaluating model outputs using process‐oriented observational data remain difficult. This study presents an ensemble of MCS simulations with spanning the deep convective gray zone ( from 12 km to 125 m) in the Southern Great Plains of the U.S. and the Amazon Basin. Comparing these simulations with Atmospheric Radiation Measurement (ARM) wind profiler observations, we find greater sensitivity in the Amazon Basin compared to the Great Plains. Convective drafts converge structurally at sub‐kilometer scales, but some deficiencies remain. In both regions, simulated up and downdrafts are too deep and extreme downdrafts are not strong enough. Furthermore, Amazonian updrafts are too strong. Overall, we observe higher sensitivity in the tropics, including an artificial buildup in vertical kinetic energy at scales of , suggesting a need for 250 m in this region. Nevertheless, bulk convergence—agreement of storm‐average statistics—is achievable with kilometer‐scale simulations within a 10% error margin with 1 km providing a good balance between accuracy and computational cost.

54 ENVIRONMENTAL SCIENCES↗

Empirical study of simulated two-planet microlensing events

We undertake the first study of two-planet microlensing models recovered from simulations of microlensing events generated by realistic multiplanet systems in which 292 planetary events, including 16 two-planet events, were detected from 6690 simulated light curves. We find that when two planets are recovered, their parameters are usually close to those of the two planets in the system most responsible for the perturbations. However, in 1 of the 16 examples, the apparent mass of both detected planets was more than doubled by the unmodeled influence of a third, massive planet. This fraction is larger than but statistically consistent with the roughly 1.5% rate of serious mass errors due to unmodeled planetary companions for the 274 cases from the same simulation in which a single planet is recovered. We conjecture that an analogous effect due to unmodeled stellar companions may occur more frequently. For 7 out of 23 cases in which two planets in the system would have been detected separately, only one planet was recovered because the perturbations due to the two planets had similar forms. This is a small fraction (7/274) of all recovered single-planet models, but almost a third of all events that might plausibly have led to two-planet models. Still, in these cases, the recovered planet tends to have parameters similar to one of the two real planets most responsible for the anomaly.

79 ASTRONOMY AND ASTROPHYSICS↗

Black hole and galaxy coevolution from continuity equation and abundance matching

We investigate the coevolution of galaxies and hosted supermassive black holes (BHs) throughout the history of the universe by a statistical approach based on the continuity equation and the abundance matching technique. Specifically, we present analytical solutions of the continuity equation without source terms to reconstruct the supermassive BH mass function from the active galactic nucleus (AGN) luminosity functions. Such an approach includes physically motivated AGN light curves tested on independent data sets, which describe the evolution of the Eddington ratio and radiative efficiency from slim- to thin-disk conditions. We nicely reproduce the local estimates of the BH mass function, the AGN duty cycle as a function of mass and redshift, along with the Eddington ratio function and the fraction of galaxies with given stellar mass hosting an AGN with given Eddington ratio. We exploit the same approach to reconstruct the observed stellar mass function at different redshift from the ultraviolet and far-IR luminosity functions associated with star formation in galaxies. These results imply that the build-up of stars and BHs in galaxies occurs via in situ processes, with dry mergers playing a marginal role at least for stellar masses ≲3×10{sup 11} M{sub ⊙} and BH masses ≲10{sup 9} M{sub ⊙}, where the statistical data are more secure and less biased by systematic errors. In addition, we develop an improved abundance matching technique to link the stellar and BH content of galaxies to the gravitationally dominant dark matter (DM) component. The resulting relationships constitute a testbed for galaxy evolution models, highlighting the complementary role of stellar and AGN feedback in the star formation process. In addition, they may be operationally implemented in numerical simulations to populate DM halos or to gauge subgrid physics. Moreover, they may be exploited to investigate the galaxy/AGN clustering as a function of redshift, mass, and/or luminosity. In fact, the clustering properties of BHs and galaxies are found to be in full agreement with current observations, thus further validating our results from the continuity equation. Finally, our analysis highlights that (i) the fraction of AGNs observed in the slim-disk regime, where most of the BH mass is accreted, increases with redshift; and (ii) already at z≳6 a substantial amount of dust must have formed over timescales ≲10{sup 8} yr in strongly star-forming galaxies, making these sources well within the reach of ALMA surveys in (sub)millimeter bands.

79 ASTRONOMY AND ASTROPHYSICS↗

Dark energy survey year 3 results: High-precision measurement and modeling of galaxy-galaxy lensing

We present and characterize the galaxy-galaxy lensing signal measured using the first three years of data from the Dark Energy Survey (DES Y3) covering 4132 deg 2 . These galaxy-galaxy measurements are used in the DES Y3 3 × 2 pt cosmological analysis, which combines weak lensing and galaxy clustering information. We use two lens samples: a magnitude-limited sample and the redmagic sample, which span the redshift range ∼ 0.2 – 1 with 10.7 and 2.6 M galaxies, respectively. For the source catalog, we use the metacalibration shape sample, consisting of ≃ 100 M galaxies separated into four tomographic bins. Our galaxy-galaxy lensing estimator is the mean tangential shear, for which we obtain a total SNR of ∼ 148 for maglim ( ∼ 120 for redmagic), and ∼ 67 ( ∼ 55 ) after applying the scale cuts of 6 Mpc / h . Thus we reach percent-level statistical precision, which requires that our modeling and systematic-error control be of comparable accuracy. The tangential shear model used in the 3 × 2 pt cosmological analysis includes lens magnification, a five-parameter intrinsic alignment model, marginalization over a point mass to remove information from small scales and a linear galaxy bias model validated with higher-order terms. We explore the impact of these choices on the tangential shear observable and study the significance of effects not included in our model, such as reduced shear, source magnification, and source clustering. We also test the robustness of our measurements to various observational and systematics effects, such as the impact of observing conditions, lens-source clustering, random-point subtraction, scale-dependent metacalibration responses, point spread function residuals, and B modes.

79 ASTRONOMY AND ASTROPHYSICS↗

Lagrangian particle tracking applied to high-speed tomographic particle imaging velocimetry

We are investigating the accuracy of Lagrangian analysis of time-resolved, 3D velocity measurements in turbulent flames obtained using high-speed tomographic particle image velocimetry (TPIV). A particle-tracking algorithm was implemented in MATLAB to extract Lagrangian statistics from the velocity fields. The accuracy of Lagrangian statistics depends on the accuracy of the estimated particle trajectories. Errors associated with the spatial and temporal resolutions of high-speed TPIV experiments were evaluated using a direct numerical simulation (DNS) of an isotropic turbulent flow. LPT was then implemented onto TPIV measurements of a laminar counterflow flame. Results calculated with LPT were compared to known results of quantities measured with an Eulerian approach.

47 OTHER INSTRUMENTATION↗

The Role of Deep Convection and Large-scale Circulation in Driving Model Spread in Low Cloud Feedback and Equilibrium Climate Sensitivity

This project aims to advance the understanding of the processes that drive the large uncertainties in climate change projections, use observations to constrain model physics and reduce the inter-model spread in equilibrium climate sensitivity (ECS). There are three major goals: 1) Characterize the representation of the physical pathways that link deep convection, large-scale circulation and low cloud feedback in CMIP6 model simulations and determine the relative contribution of each pathway to the CMIP6 model spread in low cloud feedback and ECS; 2) Use process-oriented diagnostics and multiple observations to evaluate CMIP6 model performance in capturing the observed cloud-circulation relation and deep convection characteristics including convective transition statistics and the bulk properties of mesoscale convective systems (MCSs). Error decomposition in CMIP6 models will be performed. 3) Conduct E3SM short-range hindcasts following the DOE Cloud-Associated Parameterizations Testbed (CAPT) protocol to pinpoint specific model parameters/processes that are crucial to the representation of deep convection, circulation, clouds and the pathways that connect them. We will modify convective parameters in E3SM and analyze the perturbed physics experiments (PPEs) to isolate model parameters that are critical to the uncertainty of ECS.

54 ENVIRONMENTAL SCIENCES↗

On approaches to analyze the sensitivity of simulated hydrologic fluxes to model parameters in the community land model

Here, effective sensitivity analysis approaches are needed to identify important parameters or factors and their uncertainties in complex Earth system models composed of multi-phase multi-component phenomena and multiple biogeophysical-biogeochemical processes. In this study, the impacts of 10 hydrologic parameters in the Community Land Model on simulations of runoff and latent heat flux are evaluated using data from a watershed. Different metrics, including residual statistics, the Nash-Sutcliffe coefficient, and log mean square error, are used as alternative measures of the deviations between the simulated and field observed values. Four sensitivity analysis (SA) approaches, including analysis of variance based on the generalized linear model, generalized cross validation based on the multivariate adaptive regression splines model, standardized regression coefficients based on a linear regression model, and analysis of variance based on support vector machine, are investigated. Results suggest that these approaches show consistent measurement of the impacts of major hydrologic parameters on response variables, but with differences in the relative contributions, particularly for the secondary parameters. The convergence behaviors of the SA with respect to the number of sampling points are also examined with different combinations of input parameter sets and output response variables and their alternative metrics. This study helps identify the optimal SA approach, provides guidance for the calibration of the Community Land Model parameters to improve the model simulations of land surface fluxes, and approximates the magnitudes to be adjusted in the parameter values during parametric model optimization.

58 GEOSCIENCES↗

Thermal scale modeling

Complex system study data indicate that factors associated with multilayer insulation pose major problem in scale modeling, that numerical analysis aids correction for known compromises of scaling criteria, and that probable errors in scale modeling experiments fall within range predicted by statistical analysis.

Mac Gregor, R. K.↗

Retrieval of atmospheric temperature profiles by a scanning microwave spectrometer

The Nimbus-6 satellite carries a scanning microwave spectrometer (SCAMS) experiment. The five frequency bands observed are near 22.2, 31.6, 52.8, 53.8, and 55.4 GHz. The calibration system permitted preflight calibration to an accuracy of about 1 K. In orbit, small empirical corrections were made to the calibration constants to obtain agreement in the mean of SCAMS measurements with computations based on conventional data analyzed by the National Meteorological Center (NMC). Global maps of temperature profiles were retrieved from the SCAMS measurements by a statistical method. Using the NMC analysis as the verification, RMS errors in level temperatures range of about 2-4 K, depending on altitude. Errors for layers of octave extent in pressure are uniformly about 2 K. Theoretical computations show that additional spectrometer channels would improve temperature sensing performance

Rosenkranz, P. W.↗

Maximum theoretical angular accuracy of planar and linear arrays of sensors

The maximum theoretical accuracy in the angular location of a radiating object that can be obtained by using a planar or linear array is studied. The elements are assumed to have identical radiation patterns and the complex voltages observed at their ports are assumed to be subject to phase measurement errors, having normal probability density. An optimum scheme for the statistical extraction of the parameters defining the direction is established noting that the presence of thermal noise does not affect the structure of the estimator. Comparisons with the conventional multiple interferometric techniques are made.

Borgiotti, G. V.↗

On the dynamic estimation of relative weights for observation and forecast in numerical weather prediction

The problem of merging direct and remotely sensed (indirect) data with forecast data to get an estimate of the present state of the atmosphere for the purpose of numerical weather prediction is examined. To carry out this merging optimally, it is necessary to provide an estimate of the relative weights to be given to the observations and forecast. It is possible to do this dynamically from the information to be merged, if the correlation structure of the errors from the various sources is sufficiently different. Some new statistical approaches to doing this are described, and conditions quantified in which such estimates are likely to be good.

Wahba, Grace↗

Improved Modeling in a Matlab-Based Navigation System

An innovative approach to autonomous navigation is available for low earth orbit satellites. The system is developed in Matlab and utilizes an Extended Kalman Filter (EKF) to estimate the attitude and trajectory based on spacecraft magnetometer and gyro data. Preliminary tests of the system with real spacecraft data from the Rossi X-Ray Timing Explorer Satellite (RXTE) indicate the existence of unmodeled errors in the magnetometer data. Incorporating into the EKF a statistical model that describes the colored component of the effective measurement of the magnetic field vector could improve the accuracy of the trajectory and attitude estimates and also improve the convergence time. This model is identified as a first order Markov process. With the addition of the model, the EKF attempts to identify the non-white components of the noise allowing for more accurate estimation of the original state vector, i.e. the orbital elements and the attitude. Working in Matlab allows for easy incorporation of new models into the EKF and the resulting navigation system is generic and can easily be applied to future missions resulting in an alternative in onboard or ground-based navigation.

Deutschmann, Julie↗

The GEOS Ozone Data Assimilation System: Design and Validation

An ozone data assimilation system has been developed at the Data Assimilation Office of the NASA/Goddard Space Flight Center to provide global three-dimensional analyzed ozone mixing ratio and total column ozone. The Total Ozone Mapping Spectrometer (TOMS) total column ozone and the Solar Backscatter Ultraviolet (SBUV) or SBUV/2 partial ozone profile data are assimilated. The analyzed winds from the Goddard Earth Observing System Data Assimilation System (GEOS-DAS) drive the ozone transport. Following every transport model timestep, the model prediction is combined with the observations using a global, physical-space based, statistical analysis scheme. Due to the smaller size of the ozone system than that of a global meteorological data assimilation system, new statistical analysis methodology, including anisotropic and flow-dependent forecast error correlation models, can be implemented and tested in the ozone system more easily. Sample results from the winter 1992 validation period are presented. There is a close agreement between the analyzed fields and the independent observations from ozone sondes and the Halogen Occultation Experiment (HALOE).

Stajner, I.↗