SEARCH · Engineering Papers
Results for “numerical approximations”
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.
MAGMA: Enabling exascale performance with accelerated BLAS and LAPACK for diverse GPU architectures
MAGMA (Matrix Algebra for GPU and Multicore Architectures) is a pivotal open-source library in the landscape of GPU-enabled dense and sparse linear algebra computations. With a repertoire of approximately 750 numerical routines across four precisions, MAGMA is deeply ingrained in the DOE software stack, playing a crucial role in high-performance computing. Notable projects such as ExaConstit, HiOP, MARBL, and STRUMPACK, among others, directly harness the capabilities of MAGMA. In addition, the MAGMA development team has been acknowledged multiple times for contributing to the vendors’ numerical software stacks. Looking back over the time of the Exascale Computing Project (ECP), we highlight how MAGMA has adapted to recent changes in modern HPC systems, especially the growing gap between CPU and GPU compute capabilities, as well as the introduction of low precision arithmetic in modern GPUs. We also describe MAGMA’s direct impact on several ECP projects. Maintaining portable performance across NVIDIA and AMD GPUs, and with current efforts toward supporting Intel GPUs, MAGMA ensures its adaptability and relevance in the ever-evolving landscape of GPU architectures.
Closed Loop Geothermal Working Group: GeoCLUSTER App, Subsurface Simulation Results, and Publications
To better understand the heat production, electricity generation performance, and economic viability of closed-loop geothermal systems in hot-dry rock, the Closed-Loop Geothermal Working Group -- a consortium of several national labs and academic institutions has tabulated time-dependent numerical solutions and levelized cost results of two popular closed-loop heat exchanger designs (u-tube and co-axial). The heat exchanger designs were evaluated for two working fluids (water and supercritical CO2) while varying seven continuous independent parameters of interest (mass flow rate, vertical depth, horizontal extent, borehole diameter, formation gradient, formation conductivity, and injection temperature). The corresponding numerical solutions (approximately 1.2 million per heat exchanger design) are stored as multi-dimensional HDF5 datasets and can be queried at off-grid points using multi-dimensional linear interpolation. A Python script was developed to query this database and estimate time-dependent electricity generation using an organic Rankine cycle (for water) or direct turbine expansion cycle (for CO2) and perform a cost assessment. This document aims to give an overview of the HDF5 database file and highlights how to read, visualize, and query quantities of interest (e.g., levelized cost of electricity, levelized cost of heat) using the accompanying Python scripts. Details regarding the capital, operation, and maintenance and levelized cost calculation using the techno-economic analysis script are provided. This data submission will contain results from the Closed Loop Geothermal Working Group study that are within the public domain, including publications, simulation results, databases, and computer codes. GeoCLUSTER is a Python-based web application created using Dash, an open-source framework built on top of Flask that streamlines the building of data dashboards. GeoCLUSTER provides users with a collection of interactive methods for streamlining the exploration and visualization of an HDF5 dataset. The GeoCluster app and database are contained in the compressed file geocluster_vx.zip, where the "x" refers to the version number. For example, geocluster_v1.zip is Version 1 of the app. This zip file also contains installation instructions. **To use the GeoCLUSTER app in the cloud, click the link to "GeoCLUSTER on AWS" in the Resources section below. To use the GeoCLUSTER app locally, download the geocluster_vx.zip to your computer and uncompress this file. When uncompressed this file comprises two directories and the geocluster_installation.pdf file. The geo-data app contains the HDF5 database in condensed format, and the GeoCLUSTER directory contains the GeoCLUSTER app in the subdirectory dash_app, as app.py. The geocluster_installation.pdf file provides instructions on installing Python, the needed Python modules, and then executing the app.
Stability and Convergence of Solutions to Stochastic Inverse Problems Using Approximate Probability Densities
Data-consistent inversion is designed to solve a class of stochastic inverse problems where the solution is a pullback of a probability measure specified on the outputs of a quantities of interest (QoI) map. Here, this work presents stability and convergence results for the case where finite QoI data result in an approximation of the solution as a density. Given their popularity in the literature, separate results are proven for three different approaches to measuring discrepancies between probability measures: f-divergences, integral probability metrics, and L p metrics. In the context of integral probability metrics, we also introduce a pullback probability metric that is well-suited for data-consistent inversion. This fills a theoretical gap in the convergence and stability results for data-consistent inversion that have mostly focused on convergence of solutions associated with approximate maps. Numerical results are included to illustrate key theoretical results with intuitive and reproducible test problems that include a demonstration of convergence in the measure-theoretic "almost" sense.
Inward pumping in mechanical face seals.
Hydrodynamic characteristics producing inward pumping in mechanical face seals, solving Reynolds equation for hydrodynamic lubrication with short-bearing approximation and numerical iterative methods
Numerical coefficient of the first exchange inhomogeneity term
First exchange energy inhomogeneity term numerical coefficient approximate derivation
Polymer degradation rate control of hybrid rocket combustion
Polymer degradation to small fragments is treated as a rate controlling step in hybrid rocket combustion. Both numerical and approximate analytical solutions of the complete energy and polymer chain bond conservation equations for the condensed phase are obtained. Comparison with inert atmosphere data is very good. It is found that the intersect of curves of pyrolysis rate versus interface temperature for hybrid combustors, with the thermal degradation theory, falls at a pyrolysis rate very close to that for which a pressure dependence begins to be observable. Since simple thermal degradation cannot give sufficient depolymerization at higher pyrolysis rates, it is suggested that oxidative catalysis of the process occurs at the surface, giving a first order dependence on reactive species concentration at the wall. Estimates of the ratio of this activation energy and interface temperature are in agreement with best fit procedures for hybrid combustion data. Requisite active species concentrations and flux are shown to be compatible with turbulent transport. Pressure dependence of hybrid rocket fuel regression rate is thus shown to be describable in a consistent manner in terms of reactive species catalysis of polymer degradation.
A technique for designing active control systems for astronomical telescope mirrors
This paper considers the problem of designing a control system to achieve and maintain the required surface accuracy of the primary mirror of a large space telescope. Control over the mirror surface is obtained through the application of a corrective force distribution by actuators located on the rear surface of the mirror. The design procedure is an extension of a model control technique developed for distributed parameter plants with known eigenfunctions to include plants whose eigenfunctions must be approximated by numerical techniques. Instructions are given for constructing the mathematical model of the system, and a design procedure is developed for use with. typical numerical data in selecting the number and location of the actuators. Examples of actuator patterns and their effect on various errors are given.
On electromagnetic field problems in inhomogeneous media
Analysis of electromagnetic fields in inhomogeneous media is of practical interest in general scattering and propagation problems and in the study of lenses. For certain types of inhomogeneities, the fields may be represented in terms of two scalars. In a general orthogonal coordinate system, these potentials satisfy second order differential equations. Exact solutions of these equations are known only for a few particular cases and in general, an approximate or numerical technique must be employed. The present work reviews and generalizes some of the main methods of attack of the problem. The results are presented in a form appropriate for numerical computation.
Aircraft maneuver optimization of reduced-order approximation
Review of recent work in recasting the energy type of approximation to aircraft flight in terms of singular perturbation theory and in extending this theory to three-dimensional maneuvers. Singular perturbations for differential equations arising in optimal control are first examined for a system of fairly general form but low order. The attitude dynamics for optimal flight of a rocket in vacuum is then studied as an introductory, and fairly transparent, example. The question of the choice of variables is then discussed. Finally, optimal aircraft flight in various reduced-order approximations is investigated. In particular, the problem of three-dimensional aircraft flight is formulated for singular perturbation treatment, and possibilities for decoupling into several lower-order problems are illustrated. The use of reduced-order approximation facilitates numerical computations by reducing the number of multiplier initial values that must be determined simultaneously and by improving the conditioning of the differential equations.
Neural pulse frequency modulation of an exponentially correlated Gaussian process
The effect of NPFM (Neural Pulse Frequency Modulation) on a stationary Gaussian input, namely an exponentially correlated Gaussian input, is investigated with special emphasis on the determination of the average number of pulses in unit time, known also as the average frequency of pulse occurrence. For some classes of stationary input processes where the formulation of the appropriate multidimensional Markov diffusion model of the input-plus-NPFM system is possible, the average impulse frequency may be obtained by a generalization of the approach adopted. The results are approximate and numerical, but are in close agreement with Monte Carlo computer simulation results.
On approximating hereditary dynamics by systems of ordinary differential equations
The paper deals with methods of obtaining approximate solutions to linear retarded functional differential equations (hereditary systems). The basic notion is to project the infinite dimensional space of initial functions for the hereditary system onto a finite dimensional subspace. Within this framework, two particular schemes are discussed. The first uses well-known piecewise constant approximations, while the second is a new method based on piecewise linear approximating functions. Numerical results are given.
Higher-Order Panel Method for Aerodynamic Flow Analysis
PANAIR uses high-order panel method to predict inviscid subsonic or supersonic flows about arbitrary configuration. Panel method solves linear partial differential equation numerically by approximating configuration surface with panels on which unknown "singularity strengths" are defined. PANAIR includes advanced software technology as well as advanced aerodynamic technology.
Expansion series of integral functions occurring in unsteady aerodynamics
Attention is given to two real integral functions which occur in the kernel of singular integral equations for subsonic unsteady lifting surfaces. The arguments k, r, and X of the functions correspond to the reduced frequency, spanwise distance, and modified coordinate in the flow direction, respectively. The value of the parameter nu in the functions depends on geometrical conditions. The considered investigation has the objective to present a series for general values of the nonnegative integer nu in order to compute efficiently the integral functions. The approach makes it possible to avoid any approximation or numerical quadrature.
The galactic spheroid
Analytic approximations and numerical simulations are used to derive the characteristic behavior of star counts in a galactic spheroidal population, whose visible stars' stellar luminosity function is obtainable for galactocentric distances between about 4 and 12 kpc from star count observations above 30 deg galactic latitude. The total densities of stars and of mass in the spheroid at the solar position are evaluated using different assumed luminosity functions, in order to extrapolate the measured values to a wider range of absolute magnitudes. The upper limits to the frequency of intermediate population stars are derived for the absolute visual magnitude range of 5-8. If such stars occur in either a flat disk with a scale height of 3 kpc or a spheroid with an ellipticity of 0.5, their local surface density is less than 1.8 times that of the spheroid.
Tearing instability in an anisotropic neutral sheet
A study is made of the collisionless tearing-mode stability properties of a field-reversed plasma layer whose temperature distribution is anisotropic. The plasma is confined by its self magnetic field with no external field. A kinetic description is used for both ions and electrons. The effects of the axis-crossing and nonaxis-crossing orbits are discussed. It is found that the conventional technique of matching the inner and outer asymptotic solutions at the electron inner-region is inadequate for the anisotropic case. An intermediate region in which the axis-crossing ion orbits are important is identified. The eigenvalue equation is solved using both analytic approximations and numerical methods to obtain the eigenmode structure and the linear dispersion relation. Previously announced in STAR as N84-14920
Unimolecular diffusion-mediated reactions with a nonrandom time-modulated absorbing barrier
A diffusion-reaction model with time-dependent reactivity is formulated and applied to unimolecular reactions. The model is solved exactly numerically and approximately analytically for the unreacted fraction as a function of time. It is shown that the approximate analytical solution is valid even when the system is far from equilibrium, and when the reactivity probability is more complicated than a square-wave function of time. A discussion is also given of an approach to problems of this type using a stochastically fluctuating reactivity, and the first-passage time for a particular example is derived.
Hard X-ray delays
High time resolution hard X-ray rates with good counting statistics over 5 energy intervals were obtained using a large area balloon-borne scintillation detector during the 27 June 1980 solar flare. The impulsive phase of the flare was comprised of a series of major bursts of several to several tens of seconds long. Superimposed on these longer bursts are numerous smaller approximately 0.5 to 1.0 second spikes. The time profiles for different energies were cross-correlated for the major bursts. The rapid burst decay rates and the simultaneous peaks below 120 keV both indicate a rapid electron energy loss process. Thus, the flux profiles reflect the electron acceleration/injection process. The fast rate data was obtained by a burst memory in 8 and 32 msec resolution over the entire main impulsive phase. These rates will be cross-correlated to look for short time delays and to find rapid fluctuations. However, a cursory examination shows that almost all fluctuations, down to the 5% level, were resolved with 256 msec bins.