Search NASA⌕ Search

SEARCH · Search NASA

Results for “TRUNCATION 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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Bayesian analysis of nucleon-nucleon scattering data in pionless effective field theory

We perform Bayesian model calibration of two-nucleon (NN) low-energy constants (LECs) appearing in an NN interaction based on pionless effective field theory (πEFT). The calibration is carried out for potentials constructed using naive dimensional analysis in NN relative momenta (p) up to next-to-leading order [NLO, O(p 2 )] and next-to-next-to-next-to-leading order [N3LO, O(p 4 )]. We consider two classes of pionless πEFT potential: one that acts in all partial waves and another that is dominated by s-wave physics. The two classes produce broadly similar results for calibrations to NN data up to E lab = 5 MeV. Our analysis accounts for the correlated uncertainties that arise from the truncation of the pionless πEFT. We simultaneously estimate both the πEFT LECs and the parameters that quantify the truncation error. This permits the first quantitative estimates of the pionless πEFT breakdown scale, Λ b : the 95% intervals are Λ b ∈[50.11,63.03] MeV at NLO and Λ b ∈[72.27,88.54] MeV at N3LO. Furthermore, invoking naive dimensional analysis for the NN potential, therefore, does not lead to consistent results across orders in pionless πEFT. This exemplifies the possible use of Bayesian tools to identify inconsistencies in a proposed EFT power counting.

Bayesian methods↗
DOE OSTI2025

From chiral effective field theory to perturbative QCD: A Bayesian model mixing approach to symmetric nuclear matter

Constraining the equation of state (EOS) of strongly interacting, dense matter is the focus of intense experimental, observational, and theoretical effort. Chiral effective field theory (𝜒⁢EFT ) can describe the EOS between the typical densities of nuclei and those in the outer cores of neutron stars, while perturbative QCD (pQCD) can be applied to properties of deconfined quark matter, both with quantified theoretical uncertainties. However, describing the full range of densities in between with a single EOS that has well-quantified uncertainties is a challenging problem. Bayesian multimodel inference from 𝜒⁢EFT and pQCD can help bridge the gap between the two theories. In this work, we introduce a correlated Bayesian model mixing framework that uses a Gaussian process (GP) to assimilate different information into a single QCD EOS for symmetric nuclear matter. The present implementation uses a stationary GP to infer this mixed EOS solely from the EOSs of 𝜒⁢EFT and pQCD while accounting for the truncation errors of each theory. The GP is trained on the pressure as a function of number density in the low- and high-density regions where 𝜒⁢EFT and pQCD are, respectively, valid. We impose priors on the GP kernel hyperparameters to suppress unphysical correlations between these regimes. This, together with the assumption of stationarity, results in smooth 𝜒⁢EFT-to-pQCD curves for both the pressure and the speed of sound. We show that using uncorrelated mixing requires uncontrolled extrapolation of at least one of 𝜒⁢EFT or pQCD into regions where the perturbative series breaks down and leads to an acausal EOS. Here, we also discuss extensions of this framework to nonstationary and less differentiable GP kernels, its future application to neutron-star matter, and the incorporation of additional constraints from nuclear theory, experiment, and multimessenger astronomy.

Bayesian methods↗
NASA NTRS1969

Thermal convection at infinite Prandtl number

Numerical solutions are developed for the full scaled three dimensional thermal convection problem, at infinite Prandtl number, and for rigid boundaries. The procedure is designed to give an approximate account of horizontally homogeneous and isotropic flow situations. Results, which included a maximum of 36 horizontal wave number vectors and 20 vertical wave numbers, appear to adequately describe the flow fill up to R = 100,000; beyond this R the results appear to show horizontal wave number truncation error. This error seems to affect the boundary slope of the mean temperature field more than other mean quantities. Despite some numerical uncertainities, certain of the qualitative features of the flow fill are predicted with reasonable confidence.

Herring, J. R.↗
NASA NTRS1972

Classical eight- and lower-order Runge-Kutta-Nystroem formulas with stepsize control for special second-order differential equations

The formulas include a stepsize control procedure, based on a complete coverage of the leading term of the truncation error in x. The formulas require fewer evaluations per stop than other Runge-Kutta-Nystrom formulas if the latter are operated by using the standard procedure for stepsize control. An example is presented. With results being of the same accuracy, Runge-Kutta-Nystrom formulas discussed save 50 percent or more computer time compared with other Runge-Kutta-Nystrom formulas.

Fehlberg, E.↗
NASA NTRS1972

Power series evaluation of transition and covariance matrices.

Reexamination power series solutions to the matrix covariance differential equation and the transition differential equation. Truncation error bounds are derived which are computationally attractive and which extend previous results. Polynomial approximations are obtained by exploiting the functional equations satisfied by the transition and covariance matrices. The series-functional equation propagation technique represents a fast and accurate alternative to the numerical integration of the time-invariant transition and covariance equations.

Bierman, G. J.↗
NASA NTRS1972

Numerical and Experimental Studies of the Natural Convection Flow Within a Horizontal Cylinder Subjected to a Uniformly Cold Wall Boundary Condition

Numberical solutions are obtained for the quasi-compressible Navier-Stokes equations governing the time dependent natural convection flow within a horizontal cylinder. The early time flow development and wall heat transfer is obtained after imposing a uniformly cold wall boundary condition on the cylinder. Solutions are also obtained for the case of a time varying cold wall boundary condition. Windware explicit differ-encing is used for the numerical solutions. The viscous truncation error associated with this scheme is controlled so that first order accuracy is maintained in time and space. The results encompass a range of Grashof numbers from 8.34 times 10,000 to 7 times 10 to the 7th power which is within the laminar flow regime for gravitationally driven fluid flows. Experiments within a small scale instrumented horizontal cylinder revealed the time development of the temperature distribution across the boundary layer and also the decay of wall heat transfer with time.

Stewart, R. B.↗
NASA NTRS1973

A linear shift-invariant image preprocessing technique for multispectral scanner systems

A linear shift-invariant image preprocessing technique is examined which requires no specific knowledge of any parameter of the original image and which is sufficiently general to allow the effective radius of the composite imaging system to be arbitrarily shaped and reduced, subject primarily to the noise power constraint. In addition, the size of the point-spread function of the preprocessing filter can be arbitrarily controlled, thus minimizing truncation errors.

Mcgillem, C. D.↗
NASA NTRS1972

Generation and application of the equations of condition for high order Runge-Kutta methods

This thesis develops the equations of condition necessary for determining the coefficients for Runge-Kutta methods used in the solution of ordinary differential equations. The equations of condition are developed for Runge-Kutta methods of order four through order nine. Once developed, these equations are used in a comparison of the local truncation errors for several sets of Runge-Kutta coefficients for methods of order three up through methods of order eight.

Haley, D. C.↗
NASA NTRS1973

Improving the accuracy of angular-momentum projection

A connection is established between Ullah's new method of angular-momentum projection and the conventional Hill-Wheeler method. They are studied for the case where series truncation of some sort is required. It is shown that for a particular choice of angles, the analysis simplifies greatly and at the same time leads to reduced truncation error.

Ford, W. F.↗
NASA NTRS1973

A dynamic transformation method for modal synthesis.

This paper presents a condensation method for large discrete parameter vibration analysis of complex structures that greatly reduces truncation errors and provides accurate definition of modes in a selected frequency range. A dynamic transformation is obtained from the partitioned equations of motion that relates modes not explicity in the condensed solution to the retained modes at a selected system frequency. The generalized mass and stiffness matrices, obtained with existing modal synthesis methods, are reduced using this transformation and solved. Revised solutions are then obtained using new transformations at the calculated eigenvalues and are also used to assess the accuracy of the results. If all the modes of interest have not been obtained, the results are used to select a new set of retained coordinates and a new transformation frequency, and the procedure is repeated for another group of modes.

Kuhar, E. J.↗
NASA NTRS1973

Constrained optimization of image restoration filters.

A preprocessing method to correct for image degradation is proposed which can be thought of as a generalization and extension of previous work by Smith (1966) and Stuller (1972). This method accomodates the problem of noncircularly symmetric imaging system point-spread functions, provides for controlled extent of the preprocessing filter to minimize distortion due to transients resulting from truncation errors and edge effects, can be used with various kinds of system noise, and can be readily extended to provide constraint of other system parameters. The analysis relates to a line-scanner system, although it is applicable in principle to many other system configurations.

Riemer, T. E.↗
NASA NTRS1974

Classical seventh-, sixth-, and fifth-order Runge-Kutta-Nystrom formulas with stepsize control for general second-order differential equations

Runge-Kutta-Nystrom formulas of the seventh, sixth, and fifth order were derived for the general second order (vector) differential equation written as the second derivative of x = f(t, x, the first derivative of x). The formulas include a stepsize control procedure, based on a complete coverage of the leading term of the local truncation error in x, and they require no more evaluations per step than the earlier Runge-Kutta formulas for the first derivative of x = f(t, x). The developed formulas are expected to be time saving in comparison to the Runge-Kutta formulas for first-order differential equations, since it is not necessary to convert the second-order differential equations into twice as many first-order differential equations. The examples shown saved from 25 percent to 60 percent more computer time than the earlier formulas for first-order differential equations, and are comparable in accuracy.

Fehlberg, E.↗
NASA NTRS1974

Time elements

Time elements are introduced for use with Sundman time transformations of the type dt = r(alpha)ds for satellite equations of motion. Two time elements are given, one providing maximum accuracy for alpha = 1, the other for alpha = 2. Time elements and time transformations reduce local truncation error and Liapunov (in track) instability, and provide analytical step size control. Numerical results show accuracy improvements of more than one order of magnitude when time elements are employed with time transformations in the numerical integration of the satellite equations, compared with using time transformations alone.

Nacozy, P. E.↗
NASA NTRS1974

Numerical and experimental studies of the natural convection within a horizontal cylinder

Numerical solutions are obtained for the quasi-compressible Navier-Stokes equations governing the time-dependent natural convection within a horizontal cylinder. The early flow development and wall heat transfer are obtained after a uniformly cold wall is imposed as a boundary condition on the cylinder. Results are also obtained for a time-varying cold wall as a boundary condition with windward explicit differencing used for the numerical solutions. The viscous truncation error associated with this scheme is controlled so that first-order accuracy is maintained in time and space. Experiments within a small-scale instrumented horizontal cylinder revealed the time development of the temperature distribution across the boundary layer as well as the decay of wall heat transfer with time. Agreement between temperature distributions obtained experimentally and numerically was generally good. The time decay of the dimensionless ratio of the Nusselt number to the one-fourth power of the Grashof number is found both numerically and experimentally, and good agreement is obtained between these two results over most of the cylinder wall.

Stewart, R. B.↗
NASA NTRS1975

Downward continuation of gravity information from satellite to satellite tracking or satellite gradiometry in local areas

Integral formulas in the parameter domain are used instead of a representation by spherical harmonics. The neglected regions will cause a truncation error. The application of the discrete form of the integral equations connecting the satellite observations with surface gravity anomalies is discussed in comparison with the least squares prediction method. One critical point of downward continuation is the proper choice of the boundary surface. Practical feasibilities are in conflict with theoretical considerations. The properties of different approaches for this question are analyzed.

Rummel, R.↗