Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical 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,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

Semi-empirical calculation of depletion region width in n+-p silicon solar cells

An extension of the width calculation for impurity distributions of complementary error function and the Gaussian function, in addition to the exponential function is used in combination with experimental width values deduced from capacitance measurements at zero bias to numerically determine the width occurring in both the n+ and p junctions. The total possible error o the semi-empirical calculation is less than plus or minus 3%. Results show considerable variation in the depletion widths occurring in the n+ and p regions for different impurity profiles.

Wang, E. Y.↗
DOE OSTI2012

Stability and error analysis of nodal expansion method for convection-diffusion equation

The development, and stability and error analyses of nodal expansion method (NEM) for one dimensional steady-state convection diffusion equation is presented. Following the traditional procedure to develop NEM, the discrete formulation of the convection-diffusion equation, which is similar to the standard finite difference scheme, is derived. The method of discrete perturbation analysis is applied to this discrete form to study the stability of the NEM. The scheme based on the NEM is found to be stable for local Peclet number less than 4.644. A maximum principle is proved for the NEM scheme, followed by an error analysis carried out by applying the Maximum principle together with a carefully constructed comparison function. The scheme for the convection diffusion equation is of second-order. Numerical experiments are carried and the results agree with the conclusions of the stability and error analyses. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Spatial resolution of different discretizations over long-time for the Dirac equation with small potentials

In this report we compare the long-time error bounds and spatial resolution of finite difference methods with different spatial discretizations for the Dirac equation with small electromagnetic potentials characterized by $\varepsilon \in (0, 1]$ a dimensionless parameter. We begin with the simple and widely used finite difference time domain (FDTD) methods, and establish rigorous error bounds of them, which are valid up to the time at $O(1/\varepsilon)$. In the error estimates, we pay particular attention to how the errors depend explicitly on the mesh size $h$ and time step $\tau$ as well as the small parameter $\varepsilon$. Based on the results, in order to obtain "correct" numerical solutions up to the time at $O(1/\varepsilon)$, the $\varepsilon$-scalability (or meshing strategy requirement) of the FDTD methods should be taken as $h = O(\varepsilon^{1/2})$ and $\tau = O(\varepsilon^{1/2})$. To improve the spatial resolution capacity, we apply the Fourier spectral method to discretize the Dirac equation in space. Error bounds of the resulting finite difference Fourier pseudospectral (FDFP) methods show that they exhibit uniform spatial errors in the long-time regime, which are optimal in space as suggested by the Shannon's sampling theorem. Extensive numerical results are reported to confirm the error bounds and demonstrate that they are sharp.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Bayesian strategies for uncertainty quantification of the thermodynamic properties of materials

Reliable models of the thermodynamic properties of materials are critical for industrially relevant applications that require a good understanding of equilibrium phase diagrams, thermal and chemical transport, and microstructure evolution. The goal of thermodynamic models is to capture data from both experimental and computational studies and then make reliable predictions when extrapolating to new regions of parameter space. These predictions will be impacted by artifacts present in real data sets such as outliers, systematic errors, and unreliable or missing uncertainty bounds. Such issues increase the probability of the thermodynamic model producing erroneous predictions. We present a Bayesian framework for the selection, calibration, and quantification of uncertainty of thermodynamic property models. The modular framework addresses numerous concerns regarding thermodynamic models including thermodynamic consistency, robustness to outliers, and systematic errors by the use of hyperparameter weightings and robust Likelihood and Prior distribution choices. Furthermore, the framework’s inherent transparency (e.g. our choice of probability functions and associated parameters) enables insights into the complex process of thermodynamic assessment. We introduce these concepts through examples where the true property model is known. Additionally, we demonstrate the utility of the framework through the creation of a property model from a large set of experimental specific heat and enthalpy measurements of Hafnium metal from 0 to 4900K.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Mass-preserving time second-order explicit–implicit domain decomposition schemes for solving parabolic equations with variable coefficients

In this paper, the new mass-preserving time second-order explicit–implicit domain decomposition (DD) schemes for solving parabolic equations with various coefficients are proposed. In the schemes, first, the interface values on the interfaces of sub-domains are explicitly calculated, in which two kinds of explicit difference schemes are used. Second, the interior values in sub-domains are determined by the implicit schemes. Third, the interface values are recalculated explicitly to preserve mass. The important feature is that the developed schemes are more than first order in time step and the correction operators at the third step ensure the schemes’ mass conservative. The stability and convergence are analyzed for the two schemes. The optimal error estimates are obtained in the H{sup 1} semi-norm under a stability condition. Numerical experiments confirm the theoretical results.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A Quantized State Integrator With Second Order Errors Over Monotonic Segments

We introduce a novel quantized state integration method that is analogous to the explicit, second order Runge-Kutta technique. A feature of this method is that, for a single differential equation over an interval where the solution is monotonic, it exhibits an error proportional to the square of the integration quantum. We offer a theoretical proof of this behavior and demonstrate it with a numerical example. At the same time, an extension of the method to a system with input causes the errors to become proportional to the integration quantum. The latter observation fits established theory; the former appears to be a new result that might offer new insights into the construction of quantized state integration methods.

Mahawattege, Rasika↗
DOE OSTI2017

DG-IMEX Stochastic Galerkin Schemes for Linear Transport Equation with Random Inputs and Diffusive Scalings

In this paper, we consider the linear transport equation under diffusive scaling and with random inputs. The method is based on the generalized polynomial chaos approach in the stochastic Galerkin framework. Several theoretical aspects will be addressed. Additionally, a uniform numerical stability with respect to the Knudsen number ϵ, and a uniform in ϵ error estimate is given. For temporal and spatial discretizations, we apply the implicit–explicit scheme under the micro–macro decomposition framework and the discontinuous Galerkin method, as proposed in Jang et al. (SIAM J Numer Anal 52:2048–2072, 2014) for deterministic problem. Lastly, we provide a rigorous proof of the stochastic asymptotic-preserving (sAP) property. Extensive numerical experiments that validate the accuracy and sAP of the method are conducted.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1979

Shuttle entry trajectory reconstruction using inflight accelerometer and gyro measurements

An error analysis has been made of a Shuttle postflight entry trajectory reconstruction process to obtain trajectory state estimation errors and to assess the impact of these errors on Shuttle aerodynamic force coefficient extraction. In this analysis, the entry trajectory is assumed to be reconstructed via numerical integration of onboard accelerometer and gyro measurements and constrained to satisfy ground-based radio tracking. The trajectory state estimation errors are calculated using a Kalman-Schmidt sequential filter assuming various measurement error models and combinations of ground-based tracking. The resultant trajectory estimation errors are analyzed in a simplified perturbation process to establish the accuracy to which postflight aerodynamic force coefficients can be determined. Results are presented which show that the principal error sources affecting the trajectory reconstruction and thus the force coefficient extraction, assuming perfect atmospheric density knowledge, are the accelerometer and gyro resolution, acceleration-sensitive gyro drifts, and the alignment uncertainties associated with integration on the Shuttle.

Compton, H. R.↗
DOE OSTI2023

11-th order of accuracy for numerical solution of 3-D Poisson equation with irregular interfaces on unfitted Cartesian meshes

For the first time the optimal local truncation error method (OLTEM) with 125-point stencils and unfitted Cartesian meshes has been developed in the general 3-D case for the Poisson equation for heterogeneous materials with smooth irregular interfaces. The 125-point stencils equations that are similar to those for quadratic finite elements are used for OLTEM. The interface conditions for OLTEM are imposed as constraints at a small number of interface points and do not require the introduction of additional unknowns, i.e., the sparse structure of global discrete equations of OLTEM is the same for homogeneous and heterogeneous materials. The stencils coefficients of OLTEM are calculated by the minimization of the local truncation error of the stencil equations. These derivations include the use of the Poisson equation for the relationship between the different spatial derivatives. Such a procedure provides the maximum possible accuracy of the discrete equations of OLTEM. In contrast to known numerical techniques with quadratic elements and third order of accuracy on conforming and unfitted meshes, OLTEM with the 125-point stencils provides 11-th order of accuracy, i.e., an extremely large increase in accuracy by 8 orders for similar stencils. The numerical results show that OLTEM yields much more accurate results than high-order finite elements with much wider stencils. The increased numerical accuracy of OLTEM leads to an extremely large increase in computational efficiency. Additionally, a new post-processing procedure with the 125-point stencil has been developed for the calculation of the spatial derivatives of the primary function. The post-processing procedure includes the minimization of the local truncation error and the use of the Poisson equation. It is demonstrated that the use of the partial differential equation (PDE) for the 125-point stencils improves the accuracy of the spatial derivatives by 6 orders compared to post-processing without the use of PDE as in existing numerical techniques. At an accuracy of 0.1% for the spatial derivatives, OLTEM reduces the number of degrees of freedom by 900 - 4∙10 6 times compared to quadratic finite elements. The developed post-processing procedure can be easily extended to unstructured meshes and can be independently used with existing post-processing techniques (e.g., with finite elements).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

The pp → W(→ lν) + γ process at next-to-next-to-leading order

We present details of the calculation of the pp → W(→ lν)γ process at next-to-next-to-leading order in QCD, calculated using the jettiness slicing method. The calculation is based entirely on analytic amplitudes. Because of the radiation zero, the NLO QCD contribution from the gq channel is as important as the contribution from the Born $ q\overline{q} $ process, disrupting the normal counting of leading and sub-leading contributions. We also assess the importance of electroweak (EW) corrections, including the EW corrections to both the six-parton channel 0 →$ \overline{u} d\nu {e}^{+}\gamma g $ and the five-parton channel 0 →$ \overline{u} d\nu {e}^{+}\gamma $. Previous experimental results have been shown to agree with theoretical predictions, taking into account the large experimental errors. With the advent of run II data from the LHC, the statistical errors on the data will decrease, and will be competitive with the error on theoretical predictions for the first time. We present numerical results for $ \sqrt{s} $ = 7 and 13 TeV. Analytic results for the one-loop six-parton QCD amplitude and the tree-level seven-parton QCD amplitude are presented in appendices.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2013

Optical tomography reconstruction algorithm with the finite element method: An optimal approach with regularization tools

Highlights: •New strategies to improve the accuracy of the reconstruction through mesh and finite element parameterization. •Use of gradient filtering through an alternative inner product within the adjoint method. •An integral form of the cost function is used to make the reconstruction compatible with all finite element formulations, continuous and discontinuous. •Gradient-based algorithm with the adjoint method is used for the reconstruction. -- Abstract: Optical tomography is mathematically treated as a non-linear inverse problem where the optical properties of the probed medium are recovered through the minimization of the errors between the experimental measurements and their predictions with a numerical model at the locations of the detectors. According to the ill-posed behavior of the inverse problem, some regularization tools must be performed and the Tikhonov penalization type is the most commonly used in optical tomography applications. This paper introduces an optimized approach for optical tomography reconstruction with the finite element method. An integral form of the cost function is used to take into account the surfaces of the detectors and make the reconstruction compatible with all finite element formulations, continuous and discontinuous. Through a gradient-based algorithm where the adjoint method is used to compute the gradient of the cost function, an alternative inner product is employed for preconditioning the reconstruction algorithm. Moreover, appropriate re-parameterization of the optical properties is performed. These regularization strategies are compared with the classical Tikhonov penalization one. It is shown that both the re-parameterization and the use of the Sobolev cost function gradient are efficient for solving such an ill-posed inverse problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

A kinetic model for the transport of electrons in a graphene layer

In this article, we propose a new numerical scheme for the computation of the transport of electrons in a graphene device. The underlying quantum model for graphene is a massless Dirac equation, whose eigenvalues display a conical singularity responsible for non-adiabatic transitions between the two modes. We first derive a kinetic model which takes the form of two Boltzmann equations coupled by a collision operator modeling the non-adiabatic transitions. This collision term includes a Landau–Zener transfer term and a jump operator whose presence is essential in order to ensure a good energy conservation during the transitions. We propose an algorithmic realization of the semi-group solving the kinetic model, by a particle method. We give analytic justification of the model and propose a series of numerical experiments studying the influences of the various sources of errors between the quantum and the kinetic models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

WeakIdent: Weak formulation for identifying differential equation using narrow-fit and trimming

Data-driven identification of differential equations is an interesting but challenging problem, especially when the given data are corrupted by noise. When the governing differential equation is a linear combination of various differential terms, the identification problem can be formulated as solving a linear system, with the feature matrix consisting of linear and nonlinear terms multiplied by a coefficient vector. This product is equal to the time derivative term, and thus generates dynamical behaviors. The goal is to identify the correct terms that form the equation to capture the dynamics of the given data. We propose a general and robust framework to recover differential equations using a weak formulation with two new mechanisms, narrow-fit and trimming, for both ordinary and partial differential equations (ODEs and PDEs). The weak formulation facilitates an efficient and robust way to handle noise, and two new mechanisms, narrow-fit and trimming, improve the coefficient support and value recoveries respectively. For each sparsity level, Subspace Pursuit is utilized to find an initial set of support from the large dictionary. Then, we focus on highly dynamic regions (rows of the feature matrix), and error normalize the feature matrix in the narrow-fit step. The support is further updated via trimming the terms that contribute the least. Finally, the support set of features with the smallest Cross-Validation error is chosen as the result. A comprehensive set of numerical experiments are presented for both systems of ODEs and PDEs with various noise levels. The proposed method gives a robust recovery of the coefficients, and a significant denoising effect which can handle up to 100% noise-to-signal ratio for some equations. We compare the proposed method with several state-of-the-art algorithms for the recovery of differential equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Stability boundary and optimal operating parameter identification in milling using Bayesian learning

Here, this paper describes a novel Bayesian learning approach for stability boundary and optimal parameter identification in milling without the knowledge of the underlying tool dynamics or material cutting force coefficients. The paper is divided into two parts. First, a Bayesian learning method for stability lobe identification using test results is described. Each axial depth and spindle speed combination is characterized by a probability of stability which is updated using Bayes’ rule when a test result (stable or unstable) is made available. A novel likelihood function is defined which incorporates knowledge of the stability behavior. Numerical results show convergence to the analytical stability lobe diagram. Second, an adaptive experimental strategy to identify optimal operating parameters that maximize material removal rate is described. Numerical evaluation shows convergence to the optimal operating point with error less than 15 % within ten tests on average. The approach is validated using experimental results. Results show that the proposed method is an efficient and robust learning method to identify the stability lobe diagram and optimal operating parameters with a limited number of tests/data points.

Bayesian machine learning↗
DOE OSTI2016

Spline based least squares integration for two-dimensional shape or wavefront reconstruction

In this paper, we present a novel method to handle two-dimensional shape or wavefront reconstruction from its slopes. The proposed integration method employs splines to fit the measured slope data with piecewise polynomials and uses the analytical polynomial functions to represent the height changes in a lateral spacing with the pre-determined spline coefficients. The linear least squares method is applied to estimate the height or wavefront as a final result. Numerical simulations verify that the proposed method has less algorithm errors than two other existing methods used for comparison. Especially at the boundaries, the proposed method has better performance. The noise influence is studied by adding white Gaussian noise to the slope data. Finally, experimental data from phase measuring deflectometry are tested to demonstrate the feasibility of the new method in a practical measurement.

42 ENGINEERING↗
DOE OSTI2012

The flavonoid eupatorin inactivates the mitotic checkpoint leading to polyploidy and apoptosis

The spindle assembly checkpoint (SAC) is a conserved mechanism that ensures the fidelity of chromosome distribution in mitosis by preventing anaphase onset until the correct bipolar microtubule-kinetochore attachments are formed. Errors in SAC function may contribute to tumorigenesis by inducing numerical chromosome anomalies (aneuploidy). On the other hand, total disruption of SAC can lead to massive genomic imbalance followed by cell death, a phenomena that has therapeutic potency. We performed a cell-based high-throughput screen with a compound library of 2000 bioactives for novel SAC inhibitors and discovered a plant-derived phenolic compound eupatorin (3 Prime ,5-dihydroxy-4 Prime ,6,7-trimethoxyflavone) as an anti-mitotic flavonoid. The premature override of the microtubule drug-imposed mitotic arrest by eupatorin is dependent on microtubule-kinetochore attachments but not interkinetochore tension. Aurora B kinase activity, which is essential for maintenance of normal SAC signaling, is diminished by eupatorin in cells and in vitro providing a mechanistic explanation for the observed forced mitotic exit. Eupatorin likely has additional targets since eupatorin treatment of pre-mitotic cells causes spindle anomalies triggering a transient M phase delay followed by impaired cytokinesis and polyploidy. Finally, eupatorin potently induces apoptosis in multiple cancer cell lines and suppresses cancer cell proliferation in organotypic 3D cell culture model.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

malbacR: A Package for Standardized Implementation of Batch Correction Methods for Omics Data

Mass spectrometry is a powerful tool for identifying and analyzing small molecules, such as metabolites and lipids, in com-plex biological samples. Liquid chromatography and gas chromatography mass spectrometry studies quite commonly in-volve large numbers of samples, which can require significant time for sample preparation and analyses. To accommodate such studies, the samples are commonly split into batches. Inevitably, variations in sample handling, temperature fluctua-tion, imprecise timing, column degradation and other factors result in systematic errors or biases of the measured abundances between the batches. Numerous methods are available via R packages to assist with batch correction for small molecule om-ics data; however, since these methods were developed by different research teams, the algorithms are available in separate R packages, each with different data input and output formats. We introduce the malbacR package which consolidates eleven common batch effect correction methods for small molecule omics data into one place so users can easily implement and compare: pareto scaling, power scaling, range scaling, ComBat, EigenMS, NOMIS, RUV-random, QC-RLSC, WaveI-CA2.0, TIGER, and SERRF. The malbacR package standardizes data input and output formats across these batch correction methods. The package works in conjunction with the pmartR package, allowing users to seamlessly include batch effect cor-rection in a pmartR workflow without needing any additional data manipulation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Modeling Electron Injection at Semiconductor–Molecule Interfaces using First-Principles Dynamics Simulation: Effects of Nonadiabatic Coupling, Self-energy, and Surface Models

Excited electron transfer across semiconductor–molecule heterogeneous interfaces is central to various future electronic and optoelectronic devices. At the same time, first-principles modeling of such dynamical processes remains as a great challenge in theoretical chemistry and condensed matter physics for developing better understanding at the molecular scale. Excited electron transfer from a molecule to semiconductor surface is a particularly difficult case to model accurately, because the initial state of such an electron injection process often lies deep within the dense manifold of the conduction band states in the semiconductor. Nonadiabatic couplings and energy level alignments at such interfaces as well as the finite size error of the surface model all play important roles in numerical modeling of electron injection via first-principles theory. Here, using representative interfaces between a well-defined hydrogen-terminated Si(111) surface and series of covalently adsorbed conjugated molecules, we investigate the extent to which these theoretical and numerical considerations influence the description of electron injection at the semiconductor–molecule interface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗