Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical integration”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

High‐order symmetric cubature rules for tetrahedra and pyramids

Abstract In this article, we present an algorithm to construct high‐order fully symmetric cubature rules for tetrahedral and pyramidal elements, with positive weights and integration points that are in the interior of the domain. Cubature rules are fully symmetric if they are invariant to affine transformations of the domain. We divide the integration points into symmetry orbits where each orbit contains all the points generated by the permutation stars. These relations are represented by equality constraints. The construction of symmetric cubature rules require the solution of nonlinear polynomial equations with both inequality and equality constraints. For higher orders, we use an algorithm that consists of five sequential phases to produce the cubature rules. In the literature, symmetric numerical integration rules are available for the tetrahedron for orders p = 1 – 10, 14 , and for the pyramid up to p = 10 . We have obtained fully symmetric cubature rules for both of these elements up to order p = 20 . Numerical tests are presented that verify the polynomial‐precision of the cubature rules. Convergence studies are performed for the integration of exponential, weakly singular, and trigonometric test functions over both elements with flat and curved faces. With increase in p , improvements in accuracy is realized, though nonmonotonic convergence is observed.

Jaśkowiec, Jan↗
DOE OSTI2025

Probabilistic flux limiters

The stable numerical integration of shocks in compressible flow simulations relies on the reduction or elimination of Gibbs phenomena (unstable, spurious oscillations). A popular method to virtually eliminate Gibbs oscillations caused by numerical discretization in under-resolved simulations is to use a flux limiter. A wide range of flux limiters have been studied in the literature, with recent interest in their optimization via machine learning methods trained on high-resolution datasets. The common use of flux limiters in numerical codes as plug-and-play blackbox components makes them key targets for design improvement. Even for deterministic dynamical models, numerical uncertainty is introduced via coarse-graining required by insufficient computational power to solve all scales of motion. Conventional flux limiters are deterministic and lack the capacity to address uncertainties, both aleatoric (inherent randomness) and epistemic (modeling uncertainty due to limited knowledge), which arise in coarse-grained numerical simulations. Here, we introduce a conceptually distinct type of flux limiter that is designed to handle the effects of randomness in the model and uncertainty in model parameters. Unlike traditional single-function flux limiters, these new probabilistic flux limiters incorporate multiple flux limiting functions, each applied with a learned probability drawn from high-resolution data to mitigate the effects of uncertainty in numerical simulations. This approach departs from traditional single-function limiters by explicitly modeling and incorporating uncertainty into the shock capturing process. Using the example of Burgers' equation as a testbed, we show that a machine learned, probabilistic flux limiter may be used in a shock capturing code to more accurately capture shock profiles. In particular, we show that our probabilistic flux limiter outperforms standard limiters and can be successively improved upon (up to a point) by expanding the set of probabilistically chosen flux limiting functions.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1991

From differential to difference equations for first order ODEs

When constructing an algorithm for the numerical integration of a differential equation, one should first convert the known ordinary differential equation (ODE) into an ordinary difference equation. Given this difference equation, one can develop an appropriate numerical algorithm. This technical note describes the derivation of two such ordinary difference equations applicable to a first order ODE. The implicit ordinary difference equation has the same asymptotic expansion as the ODE itself, whereas the explicit ordinary difference equation has an asymptotic that is similar in structure but different in value when compared with that of the ODE.

Freed, Alan D.↗
DOE OSTI2024

Should We Conserve Entropy or Energy when Computing CAPE with Mixed-Phase Precipitation Physics?

Abstract The rapidly increasing resolution of global atmospheric reanalysis and climate model datasets necessitates finding methods for computing convective available potential energy (CAPE) both efficiently and accurately. To this end, this article compares two common methods for computing CAPE which conserve either energy or entropy. Inaccuracies in these computations arise from both physical and numerical errors. For instance, computing CAPE with entropy conserved results in physical errors from nonequilibrium phase transitions but minimizes numerical errors because solutions are analytic at each height. In contrast, computing CAPE with energy conserved avoids these physical errors, but accumulates numerical errors that are grid-resolution-dependent because the numerical integration of a differential equation is required. Analysis of CAPE computed with large databases of soundings from the tropical Amazon and midlatitude storm environments shows that physical errors from the entropy method are typically 1%–3% as large as CAPE, which is comparable to the numerical errors from conserving energy with grid spacing of 25 and 250 m using explicit first-order and second-order integration schemes, respectively. Errors in entropy-based CAPE calculations are also insensitive to vertical grid spacing, in contrast to energy-based calculations whose error strongly scales with the grid spacing. It is shown that entropy-based methods are advantageous when intercomparing datasets with differing vertical resolution because they produce accurate and reasonably fast results that are insensitive to grid resolution, whereas a second-order energy-based method is advantageous when analyzing data with a consistent vertical resolution because of its superior computational efficiency. Significance Statement Convective available potential energy (CAPE) is a measure of instability in the atmosphere that helps forecasters and researchers understand when and where thunderstorms will form. The purpose of this article is to identify the most efficient and accurate methods for computing CAPE. Two methods are considered here, one that relates to the entropy (a measure of thermodynamic disorder) of an air parcel and one that relates to the energy of an air parcel. Results indicate that the entropy method is most accurate and insensitive to the resolution of the data used for the calculation (which can vary considerably), whereas the energy method uses the least computation time.

Peters, John M.↗
NASA NTRS1971

Comparison of various methods for mathematical analysis of the Foucault knife edge test pattern to determine optical imperfections

The linearized integral equation for the Foucault test of a solid mirror was solved by various methods: power series, Fourier series, collocation, iteration, and inversion integral. The case of the Cassegrain mirror was solved by a particular power series method, collocation, and inversion integral. The inversion integral method appears to be the best overall method for both the solid and Cassegrain mirrors. Certain particular types of power series and Fourier series are satisfactory for the Cassegrain mirror. Numerical integration of the nonlinear equation for selected surface imperfections showed that results start to deviate from those given by the linearized equation at a surface deviation of about 3 percent of the wavelength of light. Several possible procedures for calibrating and scaling the input data for the integral equation are described.

Gatewood, B. E.↗
NASA NTRS1996

Closed-Form Evaluation of Mutual Coupling in a Planar Array of Circular Apertures

The integral expression for the mutual admittance between circular apertures in a planar array is evaluated in closed form. Very good accuracy is realized when compared with values that were obtained by numerical integration. Utilization of this closed-form expression, for all element pairs that are separated by more than one element spacing, yields extremely accurate results and significantly reduces the computation time that is required to analyze the performance of a large electronically scanning antenna array.

Bailey, M. C.↗
NASA NTRS1992

Numerical analysis of the hot-gas-side and coolant-side heat transfer in liquid rocket engine combustors

The objective of this paper are to develop a multidisciplinary computational methodology to predict the hot-gas-side and coolant-side heat transfer and to use it in parametric studies to recommend optimized design of the coolant channels for a regeneratively cooled liquid rocket engine combustor. An integrated numerical model which incorporates CFD for the hot-gas thermal environment, and thermal analysis for the liner and coolant channels, was developed. This integrated CFD/thermal model was validated by comparing predicted heat fluxes with those of hot-firing test and industrial design methods for a 40 k calorimeter thrust chamber and the Space Shuttle Main Engine Main Combustion Chamber. Parametric studies were performed for the Advanced Main Combustion Chamber to find a strategy for a proposed combustion chamber coolant channel design.

Wang, Ten-See↗
NASA NTRS1993

Symmetric Galerkin boundary formulations employing curved elements

Accounts of the symmetric Galerkin approach to boundary element analysis (BEA) have recently been published. This paper attempts to add to the understanding of this method by addressing a series of fundamental issues associated with its potential computational efficiency. A new symmetric Galerkin theoretical formulation for both the (harmonic) heat conduction and the (biharmonic) elasticity problem that employs regularized singular and hypersingular boundary integral equations (BIEs) is presented. The novel use of regularized BIEs in the Galerkin context is shown to allow straightforward incorporation of curved, isoparametric elements. A symmetric reusable intrinsic sample point (RISP) numerical integration algorithm is shown to produce a Galerkin (i.e., double) integration strategy that is competitive with its counterpart (i.e., singular) integration procedure in the collocation BEA approach when the time saved in the symmetric equation solution phase is also taken into account. This new formulation is shown to be capable of employing hypersingular BIEs while obviating the requirement of C 1 continuity, a fact that allows the employment of the popular continuous element technology. The behavior of the symmetric Galerkin BEA method with regard to both direct and iterative equation solution operations is also addressed. A series of example problems are presented to quantify the performance of this symmetric approach, relative to the more conventional unsymmetric BEA, in terms of both accuracy and efficiency. It is concluded that appropriate implementations of the symmetric Galerkin approach to BEA indeed have the potential to be competitive with, if not superior to, collocation-based BEA, for large-scale problems.

Kane, J. H.↗
NASA NTRS1976

Pines method for circular and elliptical orbits

The Pines method for trajectory integration is derived from the equations of two-body motion. The concept of the perturbation derivative is explained and a formal definition is presented. An algorithm for applying the Pines method to numerical integration of trajectories is included.

Semar, C.↗
NASA NTRS1980

A calculation procedure for viscous flow in turbomachines, volume 3

A method for analyzing the nonadiabatic viscous flow through turbomachine blade passages was developed. The field analysis is based upon the numerical integration of the full incompressible Navier-Stokes equations, together with the energy equation on the blade-to-blade surface. A FORTRAN IV computer program was written based on this method. The numerical code used to solve the governing equations employs a nonorthogonal boundary fitted coordinate system. The flow may be axial, radial or mixed and there may be a change in stream channel thickness in the through-flow direction. The inputs required for two FORTRAN IV programs are presented. The first program considers laminar flows and the second can handle turbulent flows. Numerical examples are included to illustrate the use of the program, and to show the results that are obtained.

Khalil, I.↗
NASA NTRS2021

Solving Two-Equation Turbulence Models With a Perspective on Solving Transport Equations

There are three principal objectives of this report. The first objective is to investigate frequently used two-equation models when solving the RANS equations. In this study, we consider the 2006Wilcox and the 2003 Menter Shear Stress Transport (SST) models. Also, a simple change of the Wilcox model is introduced to improve computational predictions for transonic flows. Computations for flows over two different airfoils are examined to compare these turbulence models. Effects on modeling the physics of each flow due to variations in the models, such as using either strainrate or vorticity in the turbulence production term, are considered and discussed. The second objective of this report is to not only explore but also evaluate the performance of the solution algorithm for the mean flow and the transport equations. The RANS and turbulence modeling equations aresolved in a weakly coupled manner with a diagonal implicit Runge-Kutta (DIRK) solution algorithm. Throughout this report, emphasis is given to reducing the residuals to machine zero in all flow calculations, so as to eliminate the error due to numerical integration of the discrete governing equations. The final objective is to provide a perspective on solving transport equations for turbulence modeling. Aspects of solving such stiff systems of equations, as well as various numerical difficulties and possible techniques to overcome them, are considered. Discussion is also provided concerning what is called ’numerical compatibility’, which is an essential requirement when designing a solution algorithm for solving transport equations.

turbulence modeling↗
DOE OSTI2026

Magnus method for electronic structure calculations at extreme conditions

We present the application of Magnus based methods to the solution of first order coupled ordinary differential equations in High Energy Density (HED) physics applications. Our focus is on the application to quantum mechanical methods, specifically on the solution of the radial Dirac equation for real and complex energies. HED applications require accurate solutions across a wide range of spatial and energy domains, including regimes where the solutions exhibit pronounced oscillatory behavior. Such cases pose significant computational challenges. We demonstrate that Magnus-based integrators can efficiently and accurately address these challenges. We discuss the implementation of the Magnus method for the solution of the radial Dirac equation, including practical considerations such as the evaluation of matrix exponentials, numerical integration, error estimation, and adaptive step size control. We also discuss the application of these methods to complex energy Green’s function techniques and the efficient approximation of integrals of the solutions relevant to HED electronic structure calculations. Here, we demonstrate the accuracy and robustness of the resulting method in applications to the free-particle case, for which analytic solutions are available for comparison, as well as the challenging case of gold at HED conditions.

general physics↗
NASA NTRS1993

The obliquity of Mars and 'climate friction'

A mathematical theory is presented which explains the increase of the mean obliquity of Mars over geologic time due to the so called 'climate friction' (i.e., the climatic changes associated with obliquity oscillations of Mars). The theory is compared with a 10 m.y. numerical integration of the equations performed for a hypothetically large amount of climate friction for two cases of the obliquity oscillations: (1) a single sinusoid and (2) a sum of three sinusoids. The theory and numerics agree for both cases within about 12 percent on the size of the secular increase in obliquity. One possible mechanism of climate friction investigated is 'postglacial rebound' on Mars. According to this theory, giant polar caps form when the obliquity is low, and slowly squeeze out an equatorial bulge. When the obliquity is high, the caps disappear, but the bulge takes some time to collapse, due to mantle viscosity, causing it to oscillate out of phase with the obliquity oscillations. This causes a secular increase in the average obliquity.

Rubincam, David P.↗
NASA NTRS1979

A semi-analytic theory for the motion of a close-earth artificial satellite with drag

A semi-analytic method is used to estimate the decay history/lifetime and to generate orbital ephemeris for close-earth satellites perturbed by the atmospheric drag and earth oblateness due to the spherical harmonics J2, J3, and J4. The theory maintains efficiency through the application of the theory of a method of averaging and employs sufficient numerical emphasis to include a rather sophisticated atmospheric density model. The averaged drag effects with respect to mean anomaly are evaluated by a Gauss-Legendre quadrature while the averaged variational equations of motion are integrated numerically with automatic step size and error control.

Liu, J. J. F.↗
NASA NTRS1985

Advanced three-dimensional dynamic analysis by boundary element methods

Advanced formulations of boundary element method for periodic, transient transform domain and transient time domain solution of three-dimensional solids have been implemented using a family of isoparametric boundary elements. The necessary numerical integration techniques as well as the various solution algorithms are described. The developed analysis has been incorporated in a fully general purpose computer program BEST3D which can handle up to 10 subregions. A number of numerical examples are presented to demonstrate the accuracy of the dynamic analyses.

Banerjee, P. K.↗
NASA NTRS1980

Computer prediction of large reflector antenna radiation properties

A FORTRAN program for calculating reflector antenna radiation patterns was rewritten and extended to include reflectors composed of a number of panels. These individual panels must be analytic surfaces. The theoretical foundation for the program is as follows: Geometrical optics techniques are used to trace rays from a feed antenna to the reflector surface and back to a mathematical plane just in front of the reflector. The resulting tangential electric field distribution, assumed to be the only source of forward radiation, is integrated numerically to calculate the radiation pattern for a desired set of angles. When the reflector is composed of more than one panel, each panel is treated as a separated antenna, the ray-tracing procedure and integration being repeated for each panel. The results of the individual aperture plane integrations are stored and summed to yield the relative electric field strength over the angles of interest. An example and several test cases are included to demonstrate the use of the program and verify the new method of computation.

Botula, A.↗
NASA NTRS2003

Updated Integrated Mission Program

Integrated Mission Program (IMP) is a computer program for simulating spacecraft missions around the Earth, Moon, Mars, and/or other large bodies. IMP solves the differential equations of motion by use of a Runge-Kutta numerical-integration algorithm. Users control missions through selection from a large menu of events and maneuvers. Mission profiles, time lines, propellant requirements, feasibility analyses, and perturbation analyses can be computed quickly and accurately. A prior version of IMP, written in FORTRAN 77, was reported in Program Simulates Spacecraft Missions (MFS-28606), NASA Tech Briefs, Vol. 17, No. 4 (April 1993), page 60. The present version, written in double-precision Lahey FORTRAN 90, incorporates a number of improvements over the prior version. Some of the improvements modernize the code to take advantage of today's greater central-processing-unit speeds. Other improvements render the code more modular; provide additional input, output, and debugging capabilities; and add to the variety of maneuvers, events, and means of propulsion that can be simulated. The IMP user manuals (of which there are now ten, each addressing a different aspect of the code and its use) have been updated accordingly.

Dauro, Vincent A., Sr.↗
DOE OSTI2021

Evaluation of gamma-ray transmission through rectangular collimator slits for application in nuclear fuel spectrometry

Gamma-ray spectrometry is widely applied in several science fields, and in particular in non-destructive gamma scanning and gamma emission tomography of irradiated nuclear fuel. Usually, a collimator is used in the experimental setup, to selectively interrogate a region of interest in the fuel. For the optimization of instrument design, as well as for planning measurement campaigns, predictive models for the transmitted gamma-ray intensity through the collimator are needed. Commonly, Monte Carlo Radiation Transport tools are used for accurate prediction of gamma-ray transport, however, the long computation time requirements when used in low-efficiency experimental setups present challenges. Here, the full-energy peak intensity transmitted through a rectangular collimator slit was examined. A uniform planar surface source emitting isotropically was considered, and the rate of photons reaching an ideal counter plane on the opposite side of the collimator was evaluated by analytical integration. To find a closed-form primitive function, some idealizations were required, and thereby parametric models were obtained for the optical field of view, dependent on slit dimensions (length, height and width) and source-to-collimator distance. For contributions from outside the optical field of view, where a closed-form expression cannot be found, fast numerical integral methods were instead used. The results were validated using the Monte Carlo code MCNP6 and show an agreement within three percent for the numerical method. For the analytical method, deviations up to tens of percent were obtained, which is deemed to still be sufficient for instrument design and measurement planning, where often the order of magnitude of the count rate is not a priori known. The method is planned for use in iterative optimization routines in the design of Gamma Emission Tomography devices, as well as for the prediction of gamma spectra obtained in the planning of fuel inspections. An application of the proposed method was demonstrated in spectrum prediction for a short cooling-time fuel rod test from the Halden reactor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗