Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Error 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 505 records · Page 28

Recovery of the Antarctic Ozone Hole

The Antarctic ozone hole develops each year and culminates by early Spring. Antarctic ozone values have been monitored since 1979 using satellite observations from the TOMS and OMI instruments. The severity of the hole has been assessed using the minimum total ozone value from the October monthly mean (depth of the hole), the average size during the September-October period, and the ozone mass deficit. Ozone is mainly destroyed by halogen catalytic cycles, and these losses are modulated by temperature variations in the collar of the polar lower stratospheric vortex. In this presentation, we show the relationships of halogens and temperature to both the size and depth of the hole. Because atmospheric halogen levels are responding to international agreements that limit or phase out production, the amount of halogens in the stratosphere should decrease over the next few decades. We use two methods to estimate ozone hole recovery. First, we use projections of halogen levels combined with age-of-air estimates in a parametric model. Second, we use a coupled chemistry climate model to assess recovery. We find that the ozone hole is recovering at an extremely slow rate and that large ozone holes will regularly recur over the next 2 decades. Furthermore, full recovery to 1980 levels will not occur until approximately 2068. We will also show some error estimates of these dates and the impact of climate change on the recovery.

Newman, Paul A.↗

A $C^1$-Conforming Arbitrary-Order Two-Dimensional Virtual Element Method for the Fourth-Order Phase-Field Equation

We present a two-dimensional conforming virtual element method for the fourth-order phase-field equation. Our proposed numerical approach to the solution of this high-order phase-field (HOPF) equation relies on the design of an arbitrary-order accurate, virtual element space with $C^1$ global regularity. Such regularity is guaranteed by taking the values of the virtual element functions and their full gradient at the mesh vertices as degrees of freedom. Attaining high-order accuracy requires also edge polynomial moments of the trace of the virtual element functions and their normal derivatives. In this work, we detail the scheme construction, and prove its convergence by deriving error estimates in different norms. A set of representative test cases allows us to assess the behavior of the method.

97 MATHEMATICS AND COMPUTING↗

Program Analyzes Errors In STAGS

EAC computer program designed for analysis of errors in results of STAGS computer program (COSMIC Program HQN-10967). Requires input data for geometry of plate, properties of material, and set of boundary conditions. These input data come from STAGS code. (The specific link between input and output data from STAGS and input data for EAC is POSTP, postprocessor program in STAGS processors.) EAC computes continuous solution from discrete results of STAGS in order to estimate error of results of STAGS. Written in FORTRAN 77.

Thurston, Gaylen A.↗

Transient thermal-structural analysis using adaptive unstructured remeshing and mesh movement

An adaptive unstructured remeshing technique is applied to transient thermal-structural analysis. The effectiveness of the technique, together with the finite element method and an error estimation technique, is evaluated by two applications which have exact solutions: (1) the steady-state thermal analysis of a plate subjected to a highly localized surface heating, and (2) the transient thermal-structural analysis of a simulated convectively cooled leading edge subjected to a translating heat source. These applications demonstrate that the remeshing technique significantly reduces the problem size as well as the analysis solution error as compared to the results produced using standard structured meshes.

Dechaumphai, Pramote↗

The Geometric Factor of Electrostatic Plasma Analyzers: A Case Study from the Fast Plasma Investigation for the Magnetospheric Multiscale mission

We report our findings comparing the geometric factor (GF) as determined from simulations and laboratory measurements of the new Dual Electron Spectrometer (DES) being developed at NASA Goddard Space Flight Center as part of the Fast Plasma Investigation on NASA's Magnetospheric Multiscale mission. Particle simulations are increasingly playing an essential role in the design and calibration of electrostatic analyzers, facilitating the identification and mitigation of the many sources of systematic error present in laboratory calibration. While equations for laboratory measurement of the Geometric Factpr (GF) have been described in the literature, these are not directly applicable to simulation since the two are carried out under substantially different assumptions and conditions, making direct comparison very challenging. Starting from first principles, we derive generalized expressions for the determination of the GF in simulation and laboratory, and discuss how we have estimated errors in both cases. Finally, we apply these equations to the new DES instrument and show that the results agree within errors. Thus we show that the techniques presented here will produce consistent results between laboratory and simulation, and present the first description of the performance of the new DES instrument in the literature.

Magnetospheric Multiscale Mission↗

A binary coherent optical receiver for the free-space channel

The structure of a coherent binary optical MAP receiver was derived and the effects of background radiation and phase estimation errors on receiver performance were evaluated. It was found that random phase errors do not affect the performance of the coherent receiver significantly, as long as the standard deviation of the phase error remains less than roughly 0.14 radians (or 8 degrees). It was also shown that background radiation is effectively suppressed by the optical homodyne receiver: Significant deterioration in receiver performance is observed only in the presence of extremely intense background fields. The most notable characteristic of the binary optical homodyne receiver therefore appears to be its ability to achieve quantum limited performance in high background environments, provided that sufficiently accurate estimates of the received phase can be obtained.

Vilnrotter, V. A.↗

Cramer-Rao Lower Bound for SoOp-R-Based Root-Zone Soil Moisture Remote Sensing

Signals of opportunity (SoOp) reflectometry (SoOp-R) is a maturing field for geophysical remote sensing as evidenced by the growing number of airborne and spaceborne experiments. As this approach receives more attention, it is worth analyzing SoOp-R’s capabilities to retrieve subsurface soil moisture (SM) by leveraging communication and navigation satellite transmitters. In this research, the CRLB is used to identify the effects of variable SoOp-R parameters on the best achievable estimation error for root-zone soil moisture (RZSM). This study investigates the use of multiple frequency, polarization, and incidence angle measurement configurations on a two-layered dielectric profile. The results also detail the effects of variable SM conditions on the capability of SoOp-R systems to predict subsurface SM. The most prevalent observation is the importance of using at least two frequencies to limit uncertainties from subsurface SM estimates. If at least two frequencies are used, the CRLB of a profile is retrievable within the root-zone depending on the surface SM content as well as the number of independent measurements of the profile. For a depth of 30 cm, it is observed that a CRLB corresponding to 4% RZSM estimation accuracy is achievable with as few as 2 dual-frequency-based SoOp-R measurements. For this depth, increasing number of measurements provided by polarization and incidence angle allow for sensing of increasingly wet SM profile structures. This study, overall, details a methodology by which SoOp-R receiver system can be designed to achieve a desired CRLB using a trade-off study between the available measurements and SM profile.

Signals of Opportunity (SoOp)↗

Coherent Lidar Design and Performance Verification

This final report summarizes the investigative results from the 3 complete years of funding and corresponding publications are listed. The first year saw the verification of beam alignment for coherent Doppler lidar in space by using the surface return. The second year saw the analysis and computerized simulation of using heterodyne efficiency as an absolute measure of performance of coherent Doppler lidar. A new method was proposed to determine the estimation error for Doppler lidar wind measurements without the need for an independent wind measurement. Coherent Doppler lidar signal covariance, including wind shear and turbulence, was derived and calculated for typical atmospheric conditions. The effects of wind turbulence defined by Kolmogorov spatial statistics were investigated theoretically and with simulations. The third year saw the performance of coherent Doppler lidar in the weak signal regime determined by computer simulations using the best velocity estimators. Improved algorithms for extracting the performance of velocity estimators with wind turbulence included were also produced.

Frehlich, Rod↗

Intra-hour Solar Irradiance Forecast in Multiple Locations using Deep Transfer Learning

In recent years, solar power system installation imposes several challenges on the operations of local and regional power grids due to the inherent variability of ground-level solar irradiance. This work proposes a novel real-time solar forecast methodology for intra-hour solar irradiance based on deep transfer learning from ground-based sky imager for time horizons ranging from 5-15 min. There are three unique aspects of the proposed methodology: (1) a Deep Learning based algorithm development which is modeled as a classification approach rather than a traditional regression approach; (2) the use of the Transfer Learning technique to show generalization capability, robustness, and portability of baseline model in the newly deployed location where availability of enough data for training is typically scarce, and (3) redefinition of point-based irradiation forecast error estimation technique with a window-based one that is more intuitive and user-friendly. The system is developed using multiple years of irradiance and sky image recording in New Jersey and one-year data from Colorado, USA. The method is validated against ground telemetry from these two locations of diverse geographic and climatic conditions. Results show that the forecasting method proposed in this work is robust and highly accurate (8% MAPE error) for multiple locations deployment.

Deep Learning, Convolution Neural Networks, transf↗

Kalman filter divergence and aircraft motion estimators

Kalman filters designed for many aerospace systems turn out to be unsatisfactory. The estimate errors become large compared to the errors predicted by the theory ('divergence'). One of the principal causes of this failure is that the system model contains states or modes that are undisturbed by the modeled process noise, and are neutrally stable (NS). One cure for such problems is periodic restarting of a time-varying Kalman filter. Other cures include minimum variance observers with eigenvalue constraints, added noise, pole-shifting, and destabilization. Several examples are given, including effective time-invariant estimators for the longitudinal and lateral motions of an airplane where several NS modes are undisturbed by wind gusts. An interpretation of these estimators as a 'strapdown IMU' without accelerometers, gimbaled gyros, or servos is given.

Bryson, A. E., Jr.↗

Consistency and convergence for numerical radiation conditions

The problem of imposing radiation conditions at artificial boundaries for the numerical simulation of wave propagation is considered. Emphasis is on the behavior and analysis of the error which results from the restriction of the domain. The theory of error estimation is briefly outlined for boundary conditions. Use is made of the asymptotic analysis of propagating wave groups to derive and analyze boundary operators. For dissipative problems this leads to local, accurate conditions, but falls short in the hyperbolic case. A numerical experiment on the solution of the wave equation with cylindrical symmetry is described. A unified presentation of a number of conditions which have been proposed in the literature is given and the time dependence of the error which results from their use is displayed. The results are in qualitative agreement with theoretical considerations. It was found, however, that for this model problem it is particularly difficult to force the error to decay rapidly in time.

Hagstrom, Thomas↗

Multilevel Monte Carlo methods for the Grad-Shafranov free boundary problem

The equilibrium configuration of a plasma in an axially symmetric reactor is described mathematically by a free boundary problem associated with the celebrated Grad-Shafranov equation. The presence of uncertainty in the model parameters introduces the need to quantify the variability in the predictions. This is often done by computing a large number of model solutions on a computational grid for an ensemble of parameter values and then obtaining estimates for the statistical properties of solutions. In this study, we explore the savings that can be obtained using multilevel Monte Carlo methods, which reduce costs by performing the bulk of the computations on a sequence of spatial grids that are coarser than the one that would typically be used for a simple Monte Carlo simulation. We examine this approach using both a set of uniformly refined grids and a set of adaptively refined grids guided by a discrete error estimator. Numerical experiments show that multilevel methods dramatically reduce the cost of simulation, with cost reductions typically on the order of 60 or more and possibly as large as 200. Furthermore, adaptive griding results in more accurate computation of geometric quantities such as x-points associated with the model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Numerical Simulations For the F-16XL Aircraft Configuration

Numerical simulations of flow around the F-16XL are presented as a contribution to the Cranked Arrow Wing Aerodynamic Project International II (CAWAPI-II). The NASA Tetrahedral Unstructured Software System (TetrUSS) is used to perform numerical simulations. This CFD suite, developed and maintained by NASA Langley Research Center, includes an unstructured grid generation program called VGRID, a postprocessor named POSTGRID, and the flow solver USM3D. The CRISP CFD package is utilized to provide error estimates and grid adaption for verification of USM3D results. A subsonic high angle-of-attack case flight condition (FC) 25 is computed and analyzed. Three turbulence models are used in the calculations: the one-equation Spalart-Allmaras (SA), the two-equation shear stress transport (SST) and the kε turbulence models. Computational results, and surface static pressure profiles are presented and compared with flight data. Solution verification is performed using formal grid refinement studies, the solution of Error Transport Equations, and adaptive mesh refinement. The current study shows that the USM3D solver coupled with CRISP CFD can be used in an engineering environment in predicting vortex-flow physics on a complex configuration at flight Reynolds numbers.

Elmiligui, Alaa A.↗

On the Training and Generalization of Deep Operator Networks

Here, we present a novel training method for deep operator networks (DeepONets), one of the most popular neural network models for operators. DeepONets are constructed by two subnetworks, namely the branch and trunk networks. Typically, the two subnetworks are trained simultaneously, which amounts to solving a complex optimization problem in a high dimensional space. In addition, the nonconvex and nonlinear nature makes training very challenging. To tackle such a challenge, we propose a two-step training method that trains the trunk network first and then sequentially trains the branch network. The core mechanism is motivated by the divide-and-conquer paradigm and is the decomposition of the entire complex training task into two subtasks with reduced complexity. Therein the Gram–Schmidt orthonormalization process is introduced which significantly improves stability and generalization ability. On the theoretical side, we establish a generalization error estimate in terms of the number of training data, the width of DeepONets, and the number of input and output sensors. Numerical examples are presented to demonstrate the effectiveness of the two-step training method, including Darcy flow in heterogeneous porous media.

deep operator networks↗

A Fast Solver for the Fractional Helmholtz Equation

The purpose of this paper is to study a Helmholtz problem with a spectral fractional Laplacian, instead of the standard Laplacian. Recently, it has been established that such a fractional Helmholtz problem better captures the underlying behavior in Geophysical Electromagnetics. In this work, we establish the well-posedness and regularity of this problem. We introduce a hybrid spectral-finite element approach to discretize it and show well-posedness of the discrete system. In addition, we derive a priori discretization error estimates. Finally, we introduce an efficient solver that scales aswell as the best possible solver for the classical integer-order Helmholtz equation. We conclude withseveral illustrative examples that confirm our theoretical findings.

97 MATHEMATICS AND COMPUTING↗

Flood basalt eruptions, comet showers, and mass extinction events

A chronology of initiation dates of the major continental flood basalt episodes has been established from compilation of published K-Ar and Ar-Ar ages of basaltic flows and related basic intrusions. The dating is therefore independent of the biostratigraphic and paleomagnetic time scales, and the estimated errors of the inititation dates are approximately + or - 4 pct. There are 11 distinct episodes of continental flood basalts known during the past 250 Myr. The data show that flood basalt episodes are generally relatively brief geologic events, with intermittent eruptions during peak output periods lasting ony 2 to 3 Myr or less. Statistical analyses suggest that these episodes may have occurred quasi-periodically with a mean cycle time of 32 + or - 1 Myr. The initiation dates of the flood basalts are close to the estimated dates of marine mass extinctions and impact-crater clusters. Although a purely internal forcing might be argued for the flood basalt volcanism, quasi-periodic comet impacts may be the trigger for both the flood basalts and the extinctions. Impact cratering models suggest that large-body impactors lead to deep initial cratering, and therefore may cause mantle disturbances and initiate mantle plume activity. The flood basalt episodes commonly mark the initiation or jump of a mantle hotspot, and are often followed by continental rifting and separation. Evidence from dynamical studies of impacts, occurrences of craters and hotspots, and the geochemistry of boundary layers is synthesized to provide a possible model of impact-generated volcanism. Flood basalt eruptions may themselves have severe effects on climate, and possibly on life. Impacts might, as a result, have led to mass extinctions through direct atmospheric disturbances, and/or indirectly through prolonged flood basalt volcanism.

Rampino, Michael R.↗

Error and Correction Analysis for the FFA@CEBAF Energy Upgrade

An energy upgrade design for the Continuous Electron Beam Accelerator Facility (CEBAF) is under development, using fixed field alternating gradient (FFA) return arcs to recirculate electron beam up to an additional five times through the accelerating structures at CEBAF. A necessary component of any large accelerator is a beam steering and optical correction system. Small environmental changes and system errors can lower beam quality or even shut down the machine; and in pursuit of the scientific mission of JLab, high quality electron beams must be delivered to the experimental halls on a predictable schedule. Correction in the novel FFA arcs of the current upgrade design is complicated by several factors. These complexities inform the choice of correction algorithm structure and parameter values. A baseline algorithm in addition to diagnostic and correction hardware configuration is presented. The effect of this correction protocol is shown with respect to estimated errors, and several possible extensions of the algorithm are discussed. This work presents an important proof of concept for the FFA@CEBAF design effort, and provides a functional correction strategy which may be simply adjusted and optimized for future design changes.

Coxe, Alex [Old Dominion Univ., Norfolk, VA (Unite↗

Trend Detection of Atmospheric Time Series: Incorporating Appropriate Uncertainty Estimates and Handling Extreme Events

This paper is aimed at atmospheric scientists without formal training in statistical theory. Its goal is to, 1) provide a critical review of the rationale for trend analysis of the time series typically encountered in the field of atmospheric chemistry; 2) describe a range of trend-detection methods; and 3) demonstrate effective means of conveying the results to a general audience. Trend detections in atmospheric chemical composition data are often challenged by a variety of sources of uncertainty, which often behave differently to other environmental phenomena such as temperature, precipitation rate, or stream flow, and may require specific methods depending on the science questions to be addressed. Some sources of uncertainty can be explicitly included in the model specification, such as autocorrelation and seasonality, but some inherent uncertainties are difficult to quantify, such as data heterogeneity and measurement uncertainty due to the combined effect of short- and long-term natural variability, instrumental stability, and aggregation of data from sparse sampling frequency. Failure to account for these uncertainties might result in an inappropriate inference of the trends and their estimation errors. On the other hand, the variation in extreme events might be interesting for different scientific questions, for example, the frequency of extremely high surface ozone events and their relevance to human health. In this study we aim to, 1) review trend detection methods for addressing different levels of data complexity in different chemical species; 2) demonstrate that the incorporation of scientifically interpretable covariates can outperform pure numerical curve fitting techniques in terms of uncertainty reduction and improved predictability; 3) illustrate the study of trends based on extreme quantiles that can provide insight beyond standard mean or median based trend estimates; and 4) present an advanced method of quantifying regional trends based on the inter-site correlations of multi-site data. All demonstrations are based on time series of observed trace gases relevant to atmospheric chemistry, but the methods can be applied to other environmental data sets.

Trace gas↗