Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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,009 records · Page 56

Error estimates to smooth solutions of semi-discrete discontinuous Galerkin methods with quadrature rules for scalar conservation laws

Here, we focus on error estimates to smooth solutions of semi-discrete discontinuous Galerkin (DG) methods with quadrature rules for scalar conservation laws. The main techniques we use are energy estimate and Taylor expansion first introduced by Zhang and Shu in (Zhang and Shu, SIAM J Num Anal 42 (2004), 641–666). We show that, with p K (piecewise polynomials of degree k) finite elements in 1D problems, if the quadrature over elements is exact for polynomials of degree (2k) error estimates of O(h k+1/2 ) are obtained for general monotone fluxes, and optimal estimates of O(h k+1 ) are obtained for upwind fluxes. For multidimensional problems, if in addition quadrature over edges is exact for polynomials of degree (2k+1), error estimates of O(h k ) are obtained for general monotone fluxes, and O(h k+1/2 ) are obtained for monotone and sufficiently smooth numerical fluxes. Numerical results validate our analysis.

97 MATHEMATICS AND COMPUTING↗

Estimation of Errors in Stress Distributions Computed in Finite Element Simulations of Polycrystals

The accuracy of the stresses predicted from crystal plasticity-based finite element formulation depends on estimation and control of the errors associated with the discretization. In the current work, the errors in the stress distribution are estimated in virtual polycrystalline samples of α-phase titanium (hexagonal close-packed phase of Ti–6Al– 4 V). To estimate the error, the stress field, which does not possess inter-element continuity, is smoothed over a grain using an L 2 projection, thereby providing continuous stress distributions with inter-element continuity. The differences between the continuous (smooth) and discontinuous (raw) stress fields are calculated at individual Gauss quadrature points and used to estimate errors for corresponding elements and grains. Error estimations are performed for a Voronoi-tessellated microstructure, an equiaxed microstructure, and two microstructures with varying grain sizes for tensile loading extending into the fully plastic regime (≈ 5% extension). Magnitudes of the errors are found to depend on microstructural characteristics, particularly the shape and size of grains. Samples having variations in grain size or having less spherical grains exhibited higher errors than samples with uniformly sized, equiaxed grains, with the size variations having a more pronounced effect. In conclusion, errors correlate with proximity to grain boundaries at small (elastic) strains and with deformation-induced features (deformation bands) at large strains.

42 ENGINEERING↗

Parameter estimation with maximal updated densities

A recently developed measure-theoretic framework solves a stochastic inverse problem (SIP) for models where uncertainties in model output data are predominantly due to aleatoric (i.e., irreducible) uncertainties in model inputs (i.e., parameters). The subsequent inferential target is a distribution on parameters. Another type of inverse problem is to quantify uncertainties in estimates of “true” parameter values under the assumption that such uncertainties should be reduced as more data are incorporated into the problem, i.e., the uncertainty is considered epistemic. A major contribution of this work is the formulation and solution of such a parameter identification problem (PIP) within the measure-theoretic framework developed for the SIP. The approach is novel in that it utilizes a solution to a stochastic forward problem (SFP) to update an initial density only in the parameter directions informed by the model output data. In other words, this method performs “selective regularization” only in the parameter directions not informed by data. Further, the solution is defined by a maximal updated density (MUD) point where the updated density defines the measure-theoretic solution to the PIP. Another significant contribution of this work is the full theory of existence and uniqueness of MUD points for linear maps with Gaussian distributions. Data-constructed Quantity of Interest (QoI) maps are also presented and analyzed for solving the PIP within this measure-theoretic framework as a means of reducing uncertainties in the MUD estimate. We conclude with a demonstration of the general applicability of the method on two problems involving either spatial or temporal data for estimating uncertain model parameters. The first problem utilizes spatial data from a stationary partial differential equation to produce a MUD estimate of an uncertain boundary condition. The second problem utilizes temporal data obtained from the state-of-the-art ADvanced CIRCulation (ADCIRC) model to obtain a MUD estimate of uncertain wind drag coefficients for a simulated extreme weather event near the Shinnecock Inlet located in the Outer Barrier of Long Island, NY, USA.

42 ENGINEERING↗

Validation of spatiodemographic estimates produced through data fusion of small area census records and household microdata

Techniques such as Iterative Proportional Fitting have been previously suggested as a means to generate new data with the demographic granularity of individual surveys and the spatial granularity of small area tabulations of censuses and surveys. This article explores internal and external validation approaches for synthetic, small area, household- and individual-level microdata using a case study for Bangladesh. Using data from the Bangladesh Census 2011 and the Demographic and Health Survey, we produce estimates of infant mortality rate and other household attributes for small areas using a variation of an iterative proportional fitting method called P-MEDM. We conduct an internal validation to determine: whether the model accurately recreates the spatial variation of the input data, how each of the variables performed overall, and how the estimates compare to the published population totals. We conduct an external validation by comparing the estimates with indicators from the 2009 Multiple Indicator Cluster Survey (MICS) for Bangladesh to benchmark how well the estimates compared to a known dataset which was not used in the original model. The results indicate that the estimation process is viable for regions that are better represented in the microdata sample, but also revealed the possibility of strong overfitting in sparsely sampled sub-populations.

42 ENGINEERING↗

Using Pressure and Volumetric Approaches to Estimate CO2 Storage Capacity in Deep Saline Aquifers

Various approaches are used to evaluate the capacity of saline aquifers to store CO 2 , resulting in a wide range of capacity estimates for a given aquifer. The two approaches most used are the volumetric “open aquifer” and “closed aquifer” approaches. We present four full-scale aquifer cases, where CO 2 storage capacity is evaluated both volumetrically (with “open” and/or “closed” approaches) and through flow modeling. These examples show that the “open aquifer” CO 2 storage capacity estimation can strongly exceed the cumulative CO 2 injection from the flow model, whereas the “closed aquifer” estimates are a closer approximation to the flow-model derived capacity. An analogy to oil recovery mechanisms is presented, where the primary oil recovery mechanism is compared to CO 2 aquifer storage without producing formation water; and the secondary oil recovery mechanism (water flooding) is compared to CO 2 aquifer storage performed simultaneously with extraction of water for pressure maintenance. This analogy supports the finding that the “closed aquifer” approach produces a better estimate of CO 2 storage without water extraction, and highlights the need for any CO 2 storage estimate to specify whether it is intended to represent CO 2 storage capacity with or without water extraction.

58 GEOSCIENCES↗

How to estimate soil organic carbon stocks of agricultural fields? perspectives using ex-ante evaluation

Estimating soil organic carbon (SOC) stocks of agricultural fields has a range of important applications from development of sustainable management practices to monitoring carbon stocks. There are many estimation strategies with the potential for more reliable estimates of SOC stock and more efficient use of soil sampling and analysis resources, especially by leveraging readily available auxiliary information such as remote sensing. However, concrete guidance for strategy selection is lacking. This study narrows this gap with a comparison of strategies for estimating deep SOC stock (0–60 cm) in a prototypical field. Using high density SOC stock measurements and simulation, we built on past studies by 1) ex-ante evaluating a large number of strategy options, 2) using a Bayesian approach to quantify the uncertainty of the comparison, and 3) considering multiple Bayesian models to assess sensitivity to this modeling choice. We found that, using readily available auxiliary information, both balanced and stratified sampling offer substantial improvements over simple random sampling. The auxiliary information most important for this improvement is a Sentinel-2 SOC index = blue / (green × red), followed by the topographic wetness index. We found that these results are robust to the choice of mapping method, but that there is uncertainty in the magnitude of improvement. Here, we recommend future studies implement this Bayesian approach for simulated ex-ante evaluation of SOC stock estimation strategies across more fields to investigate the generalizability of these findings.

54 ENVIRONMENTAL SCIENCES↗

Impact of cyber attacks on distributed compressive sensing based state estimation in power distribution grids

Modern power distribution grids suffer from multiple vulnerabilities due to the tight integration between the physical system and the cyber infrastructure. Sophisticated and malicious cyber attacks continue to adversely impact the grid operation leading to performance degradation, service interruption, and grid failure. State estimation plays an essential role in grid monitoring and advancing cyber-attack situational awareness. In this regard, this paper first proposes a distributed compressive sensing (CS) state estimation approach for an unobservable distribution grid. Further, the proposed distributed CS approach divides the distribution grid into sub-areas to perform local state estimation. Then an alternating direction method of multipliers (ADMM) based iterative information exchange among neighboring areas is employed to complete the estimation process. In this estimation process, the impact of loss of measurement data, false data injection (FDI), replay, and neighborhood cyber-attacks is analyzed. Extensive simulations are performed on the IEEE 37-bus and IEEE 123-bus standard networks to demonstrate the algorithm’s robustness to the aforementioned cyber-attacks. A quantitative analysis of computational complexity and simulation time of the distributed CS based approach is also presented.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Physics-constrained deep neural network method for estimating parameters in a redox flow battery

Here, in this paper, we present a physics-constrained deep neural network (PCDNN) method for parameter estimation in the zero-dimensional (0D) model of the vanadium redox flow battery (VRFB). In this approach, we use deep neural networks to approximate the model parameters as functions of the operating conditions. This method allows the integration of VRFB computational models as the physical constraints in the parameter learning process, leading to enhanced accuracy of parameter estimation and cell voltage prediction. Using an experimental dataset, we demonstrate that the PCDNN method can estimate model parameters for a range of operating conditions and improve the 0D model prediction of voltage compared to the 0D model prediction with constant operation-condition-independent parameters estimated with traditional inverse methods. We also demonstrate that the PCDNN approach has an improved generalization ability for estimating parameter values for operating conditions not used in the training process.

25 ENERGY STORAGE↗

srlife : a software tool for estimating the life of high temperature concentrating solar receivers. Part I – metallic receivers

Here, this paper introduces srlife, a tool for estimating the structural service life of concentrating solar power (CSP) receivers operating at high temperatures. Supporting both metallic and ceramic receiver designs, srlife is available as open-source software at https://github.com/applied-material-modeling/srlife and can be installed via the PyPi package manager (https://pypi.org). Given basic receiver geometry and incident heat flux, the tool performs thermohydraulic and structural analysis and estimates the life of a receiver. Designed for easy integration into a software stack, including solar field and levelized cost analysis, the tool can be utilized for optimizing receiver designs to meet service life and economic targets. This paper is Part I in a two-part series. Part I discusses the analysis process used to estimate the life of metallic receivers, along with a description of the required input data. Additionally, several heuristics applied within srlife can reduce analysis time significantly while maintaining accurate life estimations for metallic receivers when compared to full analyses. Several examples demonstrating the utility of srlife in receiver design are also discussed. Part II focuses on the life estimation of ceramic receivers, using time-dependent reliability analysis and various ceramic failure models implemented in srlife.

Creep-fatigue analysis↗

Cost estimation for solid waste management in industrialising regions - Precedents, problems and prospects

Highlights: Black-Right-Pointing-Pointer We review cost estimation approaches for solid waste management. Black-Right-Pointing-Pointer Unit cost method and benchmarking techniques used in industrialising regions (IR). Black-Right-Pointing-Pointer Variety in scope, quality and stakeholders makes cost estimation challenging in IR. Black-Right-Pointing-Pointer Integrate waste flow and cost models using cost functions to improve cost planning. - Abstract: The importance of cost planning for solid waste management (SWM) in industrialising regions (IR) is not well recognised. The approaches used to estimate costs of SWM can broadly be classified into three categories - the unit cost method, benchmarking techniques and developing cost models using sub-approaches such as cost and production function analysis. These methods have been developed into computer programmes with varying functionality and utility. IR mostly use the unit cost and benchmarking approach to estimate their SWM costs. The models for cost estimation, on the other hand, are used at times in industrialised countries, but not in IR. Taken together, these approaches could be viewed as precedents that can be modified appropriately to suit waste management systems in IR. The main challenges (or problems) one might face while attempting to do so are a lack of cost data, and a lack of quality for what data do exist. There are practical benefits to planners in IR where solid waste problems are critical and budgets are limited.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Application of Recursive Estimation to Heat Tracing for Groundwater/Surface‐Water Exchange

Abstract We present and demonstrate a recursive‐estimation framework to infer groundwater/surface‐water exchange based on temperature time series collected at different vertical depths below the sediment/water interface. We formulate the heat‐transport problem as a state‐space model (SSM), in which the spatial derivatives in the convection/conduction equation are approximated using finite differences. The SSM is calibrated to estimate time‐varying specific discharge using the Extended Kalman Filter (EKF) and Extended Rauch‐Tung‐Striebel Smoother (ERTSS). Whereas the EKF is suited to real‐time (“online”) applications and uses only the past and current measurements for estimation (filtering), the ERTSS is intended for near‐real time or batch‐processing (“offline”) applications and uses a window of data for batch estimation (smoothing). The two algorithms are demonstrated with synthetic and field‐experimental data and are shown to be efficient and rapid for the estimation of time‐varying flux over seasonal periods; further, the recursive approaches are effective in the presence of rapidly changing flux and (or) nonperiodic thermal boundary conditions, both of which are problematic for existing approaches to heat tracing of time‐varying groundwater/surface‐water exchange.

54 ENVIRONMENTAL SCIENCES↗

Quantum eigenvalue estimation via time series analysis

We present an efficient method for estimating the eigenvalues of a Hamiltonian H from the expectation values of the evolution operator for various times. For a given quantum state ρ, our method outputs a list of eigenvalue estimates and approximate probabilities. Each probability depends on the support of ρ in those eigenstates of H associated with eigenvalues within an arbitrarily small range. The complexity of our method is polynomial in the inverse of a given precision parameter ε, which is the gap between eigenvalue estimates. Unlike the well-known quantum phase estimation algorithm that uses the quantum Fourier transform, our method does not require large ancillary systems, large sequences of controlled operations, or preserving coherence between experiments, and is therefore more attractive for near-term applications. The output of our method can be used to estimate spectral properties of H and other expectation values efficiently, within additive error proportional to ε.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A digital twin based on OpenFAST linearizations for real-time load and fatigue estimation of land-based turbines

Monitoring a wind turbine requires intensive instrumentation, which would be too cost-prohibitive to deploy to an entire wind plant. This work presents a technique that uses readily available measurements to estimate signals that would otherwise require additional instrumentation. This study presents a digital twin concept with a focus on estimating wind speed, thrust, torque, tower-top position, and loads in the tower using supervisory control and data acquisition (SCADA) measurements. The model combines a linear state-space model obtained using OpenFAST linearizations, a wind speed estimator, and a Kalman filter algorithm that integrates measurements with the state model to perform state estimations. The measurements are: top acceleration, generator torque, pitch, and rotational speed. The article extends previous work that derived the linear state-space model using a different method. The new implementation, based on OpenFAST linearization capability, allows for a systematic extension of the method to more states, inputs, outputs, and to the offshore environment. Results from the two methods are compared, and the validation is made with additional measurements using the GE 1.5-MW turbine located at the National Renewable Energy Laboratory test site. Real-time damage equivalent loads of the tower bottom moment are estimated with an average accuracy of approximately 10%. Overall, the results from this proof of concept are encouraging, and further application of the model will be considered.

17 WIND ENERGY↗

A statistical method for estimating wood thermal diffusivity and probe geometry using in situ heat response curves from sap flow measurements

The heat pulse method is widely used to measure water flux through plants; it works by inferring the velocity of water through a porous medium from the speed at which a heat pulse is propagated through the system. No systematic, non-destructive calibration procedure exists to determine the site-specific parameters necessary for calculating sap velocity, e.g., wood thermal diffusivity and probe spacing. Such parameter calibration is crucial to obtain the correct transpiration flux density from the sap flow measurements at the plant scale; and consequently, to up-scale tree-level water fluxes to canopy and landscape scales. The purpose of this study is to present a statistical framework for estimating the wood thermal diffusivity and probe spacing simutaneously from in-situ heat response curves collected by the implanted probes of a heat ratio apparatus. Conditioned on the time traces of wood temperature following a heat pulse, the parameters are inferred using a Bayesian inversion technique, based on the Markov chain Monte Carlo sampling method. The primary advantage of the proposed methodology is that it does not require known probe spacing or any further intrusive sampling of sapwood. The Bayesian framework also enables direct quantification of uncertainty in estimated sap flow velocity. Experiments using synthetic data show that repeated tests using the same apparatus are essential to obtain reliable and accurate solutions. When applied to field conditions, these tests are conducted during different seasons and automated using the existing data logging system. The seasonality of wood thermal diffusivity is obtained as a by-product of the parameter estimation process, and it is shown to be affected by both moisture content and temperature. Empirical factors are often introduced to account for the influence of non-ideal probe geometry on the estimation of heat pulse velocity, and they are estimated in this study as well. The proposed methodology can be applied for the calibration of existing heat ratio sap flow systems at other sites. It is especially useful when an alternative transpiration calibration device, such as a lysimeter, is not available.

Bayesian inversion↗

Fast Gaussian Process Estimation for Large-Scale In Situ Inference using Convolutional Neural Networks

Exascale computing will bring with it significant I/O limitations. One foreseeable consequence of such restrictions is that the user can save only a small fraction of complex simulation data to disk for subsequent analysis. An alternative is to fit statistical models to data in situ, that is, inside the simulation as it runs. This option requires extremely fast statistical estimation to avoid slowing down the simulation. Gaussian processes (GPs) have state-of-the-art predictive performance for modeling spatial data. However, standard estimation methods for GPs scale quite poorly to large data sets as parameter estimation requires inverting a covariance matrix to the size of the data set. In the presented work, we use a convolutional neural network (CNN) to predict the GP parameters for a spatial data set, from a simulation or otherwise, rather than optimize the parameters directly. Here, our presented case study models spatial data from E3SM, the Department of Energy’s Exascale climate model. The CNN is trained on synthetic data simulated from GP models with known parameters and then applied to data from the climate simulation. In the presented examples, the neural network scheme produces parameter estimates that compare well with standard methods such as maximum likelihood estimation in predictive performance but is obtained four orders of magnitude faster.

big data↗

Extended Frequency Divider for Bus Frequency Estimation Considering Virtual Inertia from DFIGs

Accurate estimation of local bus frequency is important for effectively controlling both synchronous and nonsynchronous generators. As the power grid evolves to accommodate essential reliability services such as virtual inertia from nonsynchronous generators, conventional techniques to estimate frequency face challenges. This paper proposes a new frequency estimation method that can effectively include the inertia contributions from double-fed wind generators (DFIGs). This is achieved through the proposed extended frequency divider formula (FDF) to include the contributions of DFIGs via Thevenin equivalents. The proposed extended FDF does not suffer from numerical issues as compared to existing phasor angle derivative-based approaches. Moreover, the knowledge of the rotor speeds of the synchronous machine and DFIGs as well as of the network admittance matrix allows estimating the frequencies of all buses in the grid, thereby significantly improving system situational awareness with a limited number of measurements. Numerical results on the IEEE 39-bus power system with DFIGs show that the proposed method achieves more accurate bus frequency estimations than the original FDF formula and other approaches based on the numerical derivation of the bus voltage phase angles.

double-fed induction generator↗

Optimal Power Flow With State Estimation in the Loop for Distribution Networks

Here in this article, we propose a framework for running optimal control-estimation synthesis in distribution networks. Our approach combines a primal-dual gradient-based optimal power flow solver with a state estimation feedback loop based on a limited set of sensors for system monitoring, instead of assuming exact knowledge of all states. The estimation algorithm reduces uncertainty on unmeasured grid states based on certain online state measurements and noisy "pseudomeasurements." We analyze the convergence of the proposed algorithm and quantify the statistical estimation errors based on a weighted least-squares estimator. The numerical results on a 4521-node network demonstrate that this approach can scale to extremely large networks and provide robustness to both large pseudomeasurement variability and inherent sensor measurement noise.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Unscented Kalman Filter-based Unbiased Minimum-Variance Estimation for Nonlinear Systems with Unknown Inputs

This paper proposes an unscented Kalman filter (UKF)-based unbiased minimum-variance estimation (UMV) method for the nonlinear system with unknown inputs. By utilizing the statistical linerization, the nonlinear system and measurement functions are transformed into a “linear-like” regression form. The latter preserves the nonlinearity of the system and the measurement models. To this end, the unknown inputs can be estimated by the weighted least-squares. This “linear-like” regression form also allows us to resort to the UMV state estimation framework for the development of new nonlinear filter to handle unknown inputs. Specifically, two approaches have been developed: 1) given the estimated inputs, we derive a filter by minimizing the trace of the state error covariance matrix; 2) without input estimation, we derive the filter by minimizing the trace of the state error covariance matrix subject to a constraint imposed on the gain matrix. We prove that these two approaches provide the same results. Numerical results validate the effectiveness of the proposed method.

Unscented Kalman filter, unbiased minimumvariance ↗