SEARCH · Search NASA

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 505 records · Page 28

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Improving designer productivity

Designer and design team productivity improves with skill, experience, and the tools available. The design process involves numerous trials and errors, analyses, refinements, and addition of details. Computerized tools have greatly speeded the analysis, and now new theories and methods, emerging under the label Artificial Intelligence (AI), are being used to automate skill and experience. These tools improve designer productivity by capturing experience, emulating recognized skillful designers, and making the essence of complex programs easier to grasp. This paper outlines the aircraft design process in today's technology and business climate, presenting some of the challenges ahead and some of the promising AI methods for meeting those challenges.

Hill, Gary C.↗

NASA NTRS1992

Prediction of airfoil stall using Navier-Stokes equations in streamline coordinates

A Navier-Stokes procedure to calculate the flow about an airfoil at incidence was developed. The parabolized equations are solved in the streamline coordinates generated for an arbitrary airfoil shape using conformal mapping. A modified k-epsilon turbulence model is applied in the entire domain, but the eddy viscosity in the laminar region is suppressed artificially to simulate the region correctly. The procedure was applied to airfoils at various angles of attack, and the results are quite satisfactory for both laminar and turbulent flows. It is shown that the present choice of the coordinate system reduces the error due to numerical diffusion, and that the lift is accurately predicted for a wide range of incidence.

Choi, D. H.↗

NASA NTRS1991

Search-based model identification of smart-structure damage

This paper describes the use of a combined model and parameter identification approach, based on modal analysis and artificial intelligence (AI) techniques, for identifying damage or flaws in a rotating truss structure incorporating embedded piezoceramic sensors. This smart structure example is representative of a class of structures commonly found in aerospace systems and next generation space structures. Artificial intelligence techniques of classification, heuristic search, and an object-oriented knowledge base are used in an AI-based model identification approach. A finite model space is classified into a search tree, over which a variant of best-first search is used to identify the model whose stored response most closely matches that of the input. Newly-encountered models can be incorporated into the model space. This adaptativeness demonstrates the potential for learning control. Following this output-error model identification, numerical parameter identification is used to further refine the identified model. Given the rotating truss example in this paper, noisy data corresponding to various damage configurations are input to both this approach and a conventional parameter identification method. The combination of the AI-based model identification with parameter identification is shown to lead to smaller parameter corrections than required by the use of parameter identification alone.

Glass, B. J.↗

NASA NTRS1992

A numerical model including PID control of a multizone crystal growth furnace

This paper presents a 2D axisymmetric combined conduction and radiation model of a multizone crystal growth furnace. The model is based on a programmable multizone furnace (PMZF) designed and built at NASA Lewis Research Center for growing high quality semiconductor crystals. A novel feature of this model is a control algorithm which automatically adjusts the power in any number of independently controlled heaters to establish the desired crystal temperatures in the furnace model. The control algorithm eliminates the need for numerous trial and error runs previously required to obtain the same results. The finite element code, FIDAP, used to develop the furnace model, was modified to directly incorporate the control algorithm. This algorithm, which presently uses PID control, and the associated heat transfer model are briefly discussed. Together, they have been used to predict the heater power distributions for a variety of furnace configurations and desired temperature profiles. Examples are included to demonstrate the effectiveness of the PID controlled model in establishing isothermal, Bridgman, and other complicated temperature profies in the sample. Finally, an example is given to show how the algorithm can be used to change the desired profile with time according to a prescribed temperature-time evolution.

Panzarella, Charles H.↗

NASA NTRS1996

Surface and Atmospheric Parameter Retrieval From AVIRIS Data: The Importance of Non-Linear Effects

AVIRIS data represent a new and important approach for the retrieval of atmospheric and surface parameters from optical remote sensing data. Not only as a test for future space systems, but also as an operational airborne remote sensing system, the development of algorithms to retrieve information from AVIRIS data is an important step to these new approaches and capabilities. Many things have been learned since AVIRIS became operational, and the successive technical improvements in the hardware and the more sophisticated calibration techniques employed have increased the quality of the data to the point of almost meeting optimum user requirements. However, the potential capabilities of imaging spectrometry over the standard multispectral techniques have still not been fully demonstrated. Reasons for this are the technical difficulties in handling the data, the critical aspect of calibration for advanced retrieval methods, and the lack of proper models with which to invert the measured AVIRIS radiances in all the spectral channels. To achieve the potential of imaging spectrometry, these issues must be addressed. In this paper, an algorithm to retrieve information about both atmospheric and surface parameters from AVIRIS data, by using model inversion techniques, is described. Emphasis is put on the derivation of the model itself as well as proper inversion techniques, robust to noise in the data and an inadequate ability of the model to describe natural variability in the data. The problem of non-linear effects is addressed, as it has been demonstrated to be a major source of error in the numerical values retrieved by more simple, linear-based approaches. Non-linear effects are especially critical for the retrieval of surface parameters where both scattering and absorption effects are coupled, as well as in the cases of significant multiple-scattering contributions. However, sophisticated modeling approaches can handle such non-linear effects, which are especially important over vegetated surfaces. All the data used in this study were acquired during the 1991 Multisensor Airborne Campaign (MAC-Europe), as part of the European Field Experiment on a Desertification-threatened Area (EFEDA), carried out in Spain in June-July 1991.

Green Robert O.↗

NASA NTRS1998

Effects of Shear, Defocus, and Wavefront Error on the Theoretical Performance of the Composite Infrared Spectrometer for Cassini

The combined effects on performance of shear between the two arms, defocus of die detector, and difference in wavefront between the two arms of a Fourier transform spectrometer using cube corner retroreflectors were investigated. Performance was characterized by the amplitude of the fringe signals coming from a detector as the path-length difference was scanned. A closed-form expression was found for the combined effects of shear and defocus, and it was found that defocus had no effect in the absence of shear. The effect of wavefront error was modeled numerically and assumed to be independent of shear and defocus. Results were compared with measurements made on the breadboard and engineering model of the Composite Infrared Spectrometer for the Cassini mission to Saturn, and good agreement was found.

Martino, Antonio J.↗

NASA NTRS2001

Analysis of High Contrast Imaging Techniques for Space Based Direct Planetary Imaging

We report on our ongoing investigations of a number of techniques for direct detection and imaging of Earth-like planets around nearby stellar sources. Herein, we give a quantitative analysis of these techniques and compare and contrast them via computer simulations. The techniques we will be reporting on are Bracewell Interferometry, Nisenson Apodized Square Aperture, and Coronagraphic masking techniques. We parameterize our results with respect to wavelength, aperture size, effects of mirror speckle, both mid- and high-spatial frequency, detector and photon noise as well pointing error. The recent numerous detections of Jupiter and Saturn like planets has driven a resurgence in research of space based high contrast imaging techniques for direct planetary imaging. Work is currently ongoing for concepts for NASA's Terrestrial Planet Finder mission and a number of study teams have been funded. The authors are members of one team.

Lyon, Richard G.↗

NASA NTRS2001

Segmented Domain Decomposition Multigrid For 3-D Turbomachinery Flows

A Segmented Domain Decomposition Multigrid (SDDMG) procedure was developed for three-dimensional viscous flow problems as they apply to turbomachinery flows. The procedure divides the computational domain into a coarse mesh comprised of uniformly spaced cells. To resolve smaller length scales such as the viscous layer near a surface, segments of the coarse mesh are subdivided into a finer mesh. This is repeated until adequate resolution of the smallest relevant length scale is obtained. Multigrid is used to communicate information between the different grid levels. To test the procedure, simulation results will be presented for a compressor and turbine cascade. These simulations are intended to show the ability of the present method to generate grid independent solutions. Comparisons with data will also be presented. These comparisons will further demonstrate the usefulness of the present work for they allow an estimate of the accuracy of the flow modeling equations independent of error attributed to numerical discretization.

Celestina, M. L.↗

NASA NTRS2002

High-Order Semi-Discrete Central-Upwind Schemes for Multi-Dimensional Hamilton-Jacobi Equations

We present high-order semi-discrete central-upwind numerical schemes for approximating solutions of multi-dimensional Hamilton-Jacobi (HJ) equations. This scheme is based on the use of fifth-order central interpolants like those developed in [1], in fluxes presented in [3]. These interpolants use the weighted essentially nonoscillatory (WENO) approach to avoid spurious oscillations near singularities, and become "central-upwind" in the semi-discrete limit. This scheme provides numerical approximations whose error is as much as an order of magnitude smaller than those in previous WENO-based fifth-order methods [2, 1]. Thee results are discussed via examples in one, two and three dimensions. We also pregnant explicit N-dimensional formulas for the fluxes, discuss their monotonicity and tl!e connection between this method and that in [2].

Bryson, Steve↗

NASA NTRS2014

All-Digital Baseband 65nm PLL/FPLL Clock Multiplier using 10-cell Library

PLLs for clock generation are essential for modern circuits, to generate specialized frequencies for many interfaces and high frequencies for chip internal operation. These circuits depend on analog circuits and careful tailoring for each new process, and making them fault tolerant is an incompletely solved problem. Until now, all digital PLLs have been restricted to sampled data DSP techniques and not available for the highest frequency baseband applications. This paper presents the design and preliminary evaluation of an all-digital baseband technique built entirely with an easily portable 10-cell digital library. The library is also described, as it aids in research and low volume design porting to new processes. The advantages of the digital approach are the wide variety of techniques available to give varying degrees of fault tolerance, and the simplicity of porting the design to new processes, even to exotic processes that may not have analog capability. The only tuning parameter is digital gate delay. An all-digital approach presents unique problems and standard analog loop stability design criteria cannot be directly used. Because of the quantization of frequency, there is effectively infinite gain for very small loop error feedback. The numerically controlled oscillator (NCO) based on a tapped delay line cannot be reliably updated while a pulse is active in the delay line, and ordinarily does not have enough frequency resolution for a low-jitter output.

Shuler, Robert L., Jr.↗

NASA NTRS2014

ALL-Digital Baseband 65nm PLL/FPLL Clock Multiplier Using 10-Cell Library

Schuler, Robert L., Jr.↗

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 OSTI2023

Efficient, direct compilation of SU(N) operations into SNAP & Displacement gates

We present a function which connects the parameter of a previously published short sequence of selective number-dependent arbitrary phase (SNAP) and displacement gates acting on a qudit encoded into the Fock states of a superconducting cavity, $V_k(\alpha)=D(\alpha)R_\pi(k)D(-2\alpha)R_\pi(k)D(\alpha)$ to the angle of the Givens rotation $G(\theta)$ on levels $|k\rangle,|k+1\rangle$ that sequence approximates, namely $\alpha=\Phi(\theta) = \frac{\theta}{4\sqrt{k+1}}$. Previous publications left the determination of an appropriate $\alpha$ to numerical optimization at compile time. The map $\Phi$ gives us the ability to compile directly any $d$-dimensional unitary into a sequence of SNAP and displacement gates in $O(d^3)$ complex floating point operations with low constant prefactor, avoiding the need for numerical optimization. Numerical studies demonstrate that the infidelity of the generated gate sequence $V_k$ per Givens rotation $G$ scales as approximately $O(\theta^6)$. We find numerically that the error on compiled circuits can be made arbitrarily small by breaking each rotation into $m$$\theta/m$ rotations, with the full $d\times d$ unitary infidelity scaling as approximately $O(m^{-4})$. This represents a significant reduction in the computational effort to compile qudit unitaries either to SNAP and displacement gates or to generate them via direct low-level pulse optimization via optimal control.

Job, Joshua↗

DOE OSTI2026

Error quantification and comparison of binary neutron star gravitational waveforms from numerical relativity codes

Future gravitational wave detections of merging binary neutron star systems have the possibility to tightly constrain the equation of state of dense nuclear matter. In order to extract such constraints, gravitational waveform models need to be calibrated to accurate numerical relativity simulations of the late inspiral and merger. In this work, we take an essential step toward classifying the error and potential systematics in current generation numerical relativity simulations of merging binary neutron stars. To this end, we perform a direct comparison of two codes (fil, spec), which differ in many aspects, including the numerical methods and discretizations used and equations solved. We find that despite these different approaches, the codes are—within current numerical resolution bounds—fully consistent, and broadly comparable in cost for a given accuracy level. Our results indicate that the error in the waveforms is primarily dominated by the hydrodynamic evolution, consistent with earlier findings in the literature. Finally, we also discuss current limitations and cost estimates for numerical relativity simulations to reach the accuracies required in the era of next-generation gravitational detectors.

gravitational waves↗

DOE OSTI2017

Analysis of errors introduced by geographic coordinate systems on weather numeric prediction modeling

Most atmospheric models, including the Weather Research and Forecasting (WRF) model, use a spherical geographic coordinate system to internally represent input data and perform computations. However, most geographic information system (GIS) input data used by the models are based on a spheroid datum because it better represents the actual geometry of the earth. WRF and other atmospheric models use these GIS input layers as if they were in a spherical coordinate system without accounting for the difference in datum. When GIS layers are not properly reprojected, latitudinal errors of up to 21 km in the midlatitudes are introduced. Recent studies have suggested that for very high-resolution applications, the difference in datum in the GIS input data (e.g., terrain land use, orography) should be taken into account. However, the magnitude of errors introduced by the difference in coordinate systems remains unclear. This research quantifies the effect of using a spherical vs. a spheroid datum for the input GIS layers used by WRF to study greenhouse gas transport and dispersion in northeast Pennsylvania.

58 GEOSCIENCES↗

NASA NTRS1976

Measurements of solar magnetic fields by Fourier transform techniques. II - Saturated and blended lines

Fourier techniques have been exhaustively calibrated using Unno's (1956) results for the absorption profile of a simple Zeeman triplet. If a simple transformation is applied to the normalized line depths, then magnetic-field strengths and inclination angles can be measured very accurately from noisy saturated line profiles. Systematic errors caused by saturation effects can be estimated and reduced by varying one parameter. When a significant fraction of the line profile is unsplit and unpolarized, large errors may be made in measurements of low fields, unless the line is sufficiently weak. For a weak line, a vertical field of 1600 gauss can be measured to 10% accuracy even when 70% of the line profile is stray light. These stray-light errors are troublesome in measuring fields of gaps and pores but not sunspots. Numerical results of the error analysis are presented graphically.

Tarbell, T. D.↗

DOE OSTI2017

Thermal Modeling of a Dual-Stage Sodium Thermal-Electrochemical Converter (Na-TEC)

The sodium thermal electrochemical converter (Na-TEC) is a heat engine that generates electricity through the isothermal expansion of sodium ions through a beta”-alumina solid-electrolyte. The Na-TEC is a closed system that can theoretically achieve conversion efficiencies above 45% when operating between thermal reservoirs at 1150 K and 550 K. However, thermal designs have confined previous single-stage devices to thermal efficiencies below 20%. To mitigate some of these limitations, we consider dividing the isothermal expansion into two stages; one at the evaporator temperature (1150 K) and another at an intermediate temperature (550 K – 950 K). This dual-stage Na-TEC takes advantage of regeneration and reheating, and could be amenable to better thermal management through a reduction of parasitic losses. These losses are comprised of (i) the heat input that does not participate in power production but simply bypasses the thermo-electrochemical process and is transferred to the condenser, and (ii) heat that is lost from the device directly to the surroundings due to poor insulation. Heat loss to the surroundings occurs from a lower average hot-side temperature given that a fraction of the power conversion occurs at an intermediate stage. Also, the condenser now directly interacts with the lower temperature of the intermediate stage, thereby reducing radiation loss. A simplified geometry is developed for this device to help guide its design from a heat transfer perspective. All thermal contact and boundary conditions are discussed, including the evaporation and condensation of sodium vapor, thermal leakage through the housing, and the thermal contact between the beta”-alumina and the housing. A ray tracing method assuming gray, diffuse surfaces is used to account for the internal radiation heat transfer within the device. Liquid sodium flows in a capillary wick between the condenser and the evaporator with Pe << 1, so only conduction heat transfer is considered within this liquid return path. For simplicity, natural convection between the sodium vapor and the device housing is neglected. A numerical linear resistance model (<5% error when compared to a finite element model) is used to determine the expected conversion temperatures of each stage for a range of candidate materials. Furthermore, a numerical model (<2% error) is developed to find the factor of increase in the thermal resistance for several designs of corrugated features for the device housing. By using an iterative scheme, efficiency vs. power curves are created for a number of material and size combinations. Finally, the fundamental tradeoff between the power and efficiency is explored in conjunction with the total system size.

14 SOLAR ENERGY↗

DOE OSTI2023

Transformer Leakage Inductance Design Methodology

The leakage inductance exhibited by a transformer depends on its winding geometry, which generally involves the selection of several key design parameters in addition to the winding structure and the interleaving configuration. With few resources explaining the effects of these design choices on the observed leakage inductance, numerous trial-and-error iterations become necessary to realize the desired leakage inductance. This paper explores more than a hundred winding geometries feasible in a 2-winding transformer comprising the same magnetic core, number of turns, and wire gauge, and finds the leakage inductance for each unique design using 2-D finite element method (FEM) simulations in association with the semi-analytical double- 2-D model. These leakage inductances are plotted and further analyzed to understand the effects of different design parameters on the effective leakage inductance. The results presented herein and the conclusions drawn from this research can serve as a valuable resource for future design practitioners from both industry and academia.

24 POWER TRANSMISSION AND DISTRIBUTION↗