Search NASA⌕ Search

SEARCH · Search NASA

Results for “Approximate computing”

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,549 records · Page 86

Approximate Green's function methods for HZE transport in multilayered materials

A nonperturbative analytic solution of the high charge and energy (HZE) Green's function is used to implement a computer code for laboratory ion beam transport in multilayered materials. The code is established to operate on the Langley nuclear fragmentation model used in engineering applications. Computational procedures are established to generate linear energy transfer (LET) distributions for a specified ion beam and target for comparison with experimental measurements. The code was found to be highly efficient and compared well with the perturbation approximation.

Wilson, John W.↗

Observational Signatures of Coronal Heating in Magnetohydrodynamic Simulations without Radiation or a Lower Atmosphere

It is extremely difficult to simulate the details of coronal heating and also make meaningful predictions of the emitted radiation. Thus, testing realistic models with observations is a major challenge. Observational signatures of coronal heating depend crucially on radiation, thermal conduction, and the exchange of mass and energy with the transition region and chromosphere below. Many magnetohydrodynamic simulation studies do not include these effects, opting instead to devote computational resources to the magnetic aspects of the problem. We have developed a simple method of accounting approximately for the missing effects. It is applied to the simulation output ex post facto and therefore may be a valuable tool for many studies. We have used it to predict the emission from a model corona that is driven by vortical boundary motions meant to represent photospheric convection. We find that individual magnetic strands experience short-term brightenings, both scattered throughout the computational volume and in localized clusters. The former may explain the diffuse component of the observed corona, while the latter may explain bright coronal loops. Several observed properties of loops are reproduced reasonably well: width, lifetime, and quasi-circular cross section (aspect ratio not high). Our results lend support to the idea that loops are multistranded structures heated by "storms" of nanoflares.

James A. Klimchuk↗

Computational dispersion properties of horizontal staggered grids for atmospheric and ocean models

The computational dispersion properties of horizontally and time-horizontally staggered grids utilizing corresponding centered-difference techniques for approximation of the adjustment, or gravity wave equations, are examined in terms of their group velocity characteristics. Results are acquired for oceanic and atmospheric models, the former being characterized by a much smaller Rossby radius of deformation. For all grids considered additional filtering is required to control and even eliminate waves with poor computational dispersion characteristics. Computational dispersion properties along with other computational characteristics and requirements give some guidance for an optimal selection of an appropriate grid for an ocean or atmospheric model.

Fox-Rabinovitz, Michael S.↗

Approximations for column effect in airplane wing spars

The significance attaching to "column effect" in airplane wing spars has been increasingly realized with the passage of time, but exact computations of the corrections to bending moment curves resulting from the existence of end loads are frequently omitted because of the additional labor involved in an analysis by rigorously correct methods. The present report represents an attempt to provide for approximate column effect corrections that can be graphically or otherwise expressed so as to be applied with a minimum of labor. Curves are plotted giving approximate values of the correction factors for single and two bay trusses of varying proportions and with various relationships between axial and lateral loads. It is further shown from an analysis of those curves that rough but useful approximations can be obtained from Perry's formula for corrected bending moment, with the assumed distance between points of inflection arbitrarily modified in accordance with rules given in the report. The discussion of general rules of variation of bending stress with axial load is accompanied by a study of the best distribution of the points of support along a spar for various conditions of loading.

Warner, Edward P↗

Sensitivity analysis in multipole-accelerated panel methods for potential flow

In the design of an airframe, the effect of changing the geometry on resulting computations is necessary for design optimization. The geometry is defined in terms of a series of design variables, including design variables to define the wing planform, tail, canard, pylon, and nacelle. Design optimization in this research is based on how these design variable affect the potential flow. The potential flow is computed as a function of the geometry and location of a series of panels describing the airframe, which are in turn a function of the design variables. Multipole accelerated panel methods improve the computational complexity of the problem and thus are an attractive approach. To utilize the methods in design optimization, it was necessary to define the appropriate sensitivity derivatives. The overhead incurred from finding the sensitivity derivatives in conjunction with the original computation should be small. This research developed the background for multipole-accelerated panel methods and the framework for finding sensitivity derivatives in the methods. Potential flow panel codes are commonly used for powered-lift aerodynamic predictions for three dimensional geometries. Given an airframe which has been discretized into a series of panels to define the airframe geometry, potential is computed as a function of the influence of all panels on all other panels. This is a computationally intensive problem for which efficient solutions are desired to improve the computational time and to allow greater resolution by use of more panels. One such solution is the use of hierarchical multipole methods which entail approximations of the effects of far-field terms. Hierarchical multipole methods have become prevalent in molecular dynamics and gravitational physics, and have been introduced into the fields of capacitance calculations, computational fluid dynamics, and electromagnetics. The methods utilize multipole expansions to describe the effect of bodies (i.e. particles, astrophysical bodies, panels, etc.) within a sphere on points distant from the sphere, where the influence diminishes as a function of distance. The expansions are exact with infinite series, however, for practical computations, the series are truncated and accuracy is selected based on the number of terms retained in the expansions. A hierarchical tree structure groups bodies together based on proximity to allow definition of multipole expansions for each group. The multipole expansions are then used to compute the effect of the bodies in a group on distant bodies.

Leathrum, James F., Jr.↗

Computationally Guided Development of Components for High Energy Density Solid-State Lithium-Sulfur Batteries

All electric vertical take-off and landing vehicles (eVTOL) for urban air mobility (UAM) concepts face numerous challenging technical barriers before their introduction into the consumer marketplace. The most challenging of these technical barriers to overcome is developing an energy storage system capable of meeting the rigorous aerospace safety and performance criteria1. The performance metrics for eVTOL craft, such as specific energy, specific power, and safety, exceed those of electric automobiles by a factor of two to four. Current state-of-the-art (SOA) lithium-ion batteries are incapable of meeting the key performance criteria of energy and safety for eVTOL. Therefore, next generation advanced chemistries and designs must be developed to meet required performance metrics for electric aviation2. Beyond lithium-ion chemistries, such as lithium-sulfur, show promise in their high energy, while limitations exist in their power and cyclability due to low electrical conductivity and high intermediate solubility in organic liquid electrolytes. Several strategies to overcome the low electrical conductivity involve the use of selenium as a dopant in the active sulfur material, along with the incorporation of 2-dimensional electron-conducting holey-graphene to improve the composite cathodes electronic conductivity. Furthermore, combining this chemistry with a solid-electrolyte avoids the components’ dissolution issues3. However, the development of composite solid-state cathodes is non-trivial as several components must be intimately mixed so that the active component has sufficient access to both electrons and lithium ions to undergo full electrochemical conversion. Mathematical modeling of battery components can assist experimental design through a robust and rigorous combination of computational modeling techniques covering multiple length scales. The objective is to leverage modern computational materials methods combined with battery multiphysics tools to develop radically advanced compatible cathode and electrolyte materials, build and test solid state lithium-sulfur cells and packs. A NASA-based cross-organizational team of high-powered experts combined integrated computational predictive modeling, fundamental chemistry analysis, advanced material science, and battery cell development to tackle this very challenging, multidisciplinary problem. This presentation will show a multiscale computational modeling approach that has produced a novel particle dynamics method called Solid Electrolyte Sphere Approximation Model (SESAM). SESAM modeling targets the 1-10 µm scale structures and provides electromechanical and grain interactions for predictive design guidelines for the manufacturing of solid-state components. Parameters such as particle size and volume fraction of the constituent materials were modeled and experimentally fabricated to optimize electrochemical performance through improved microstructure design. Experimental feedback was provided through ionic and electronic conductivity assessment and structural analysis of developed materials and cell components.

battery↗

A theoretical investigation of forebody shapes designed for natural laminar boundary-layer flow

The design of forebody shapes for natural laminar flow is discussed. For subsonic flow, computed results for three shapes of different fineness ratios indicate that laminar flow can be attained under conditions that approximate those on the forebody of a cruise missile flying at a low altitude at a high subsonic Mach number. For supersonic (Mach 2.00) design, a one-parameter family of hyperbolic arcs was used to generate forebody shapes having a favorable pressure gradient over the forebody length. Computed results for these shapes indicated laminar and transitional flow over the range of Reynolds numbers considered.

Barger, R. L.↗

Improved Depiction Of Computed Unsteady Flows

PLOT4D and STREAKER programs increase ease of production and utility of animated depictions of unsteady flows simulated numerically by advanced computational fluid dynamics. Programs developed to maximize information displayed while minimizing amount stored. PLOT4D, an extended version of PLOT3D, used to depict steady three-dimensional flows, incorporates features making it more computationally efficient than to visualize unsteady flow as succession of steady flows, animating frame by frame. PLOT4D postprocesses numerical output of fluid-dynamical computation into maps of flow quantities on two-coordinate surfaces depicting cut through flow field or conform to body surface. STREAKER generates streaklines, close approximations of patterns traced by streams of dye or smoke released into simulated flows. Streaklines defined via spatial succession of endpoints of pathlines of massless particles released into flow. STREAKER computes endpoints by fourth-order Runge-Kutta integration of velocity vectors.

Smith, Merritt H.↗

Statistical Model of Evaporating Multicomponent Fuel Drops

An improved statistical model has been developed to describe the chemical composition of an evaporating multicomponent- liquid drop and of the mixture of gases surrounding the drop. The model is intended for use in computational simulations of the evaporation and combustion of sprayed liquid fuels, which are typically mixtures of as many as hundreds of different hydrocarbon compounds. The present statistical model is an approximation designed to afford results that are accurate enough to contribute to understanding of the simulated physical and chemical phenomena, without imposing an unduly large computational burden.

Harstad, Kenneth↗

The Monostatic/Bistatic Approximation

Many Radar Cross Section (RCS) prediction codes are limited to one monostatic return per run. However, such codes can calculate multiple bistatic returns per incident angle for a relatively small amount of additional computer resources. This note describes a method of using bistatic returns to generate multiple monostatic predictions for each incident angle computed. Typical results are presented and show the accuracy is initially good and then degrades as the separation angle between the incident and viewing angles becomes large. Introduction Since 1990, the monostatic/bistatic approximation has been used to reduce the number of runs required by finite-volume time-domain (FVTD) codes for making RCS versus azimuth plots. This approximation was spawned by the observation of a range test where the transmit and receive antennas were separated by a few degrees to prevent cross talk between the antennas. The measurements from this range are presented as monostatic RCS rather than bistatic: RCS. The procedure of reporting experimental bistatic RCS as the monostatic RCS at the angle bisecting the transmit and receive antennas was extended to FVTD codes and produces excellent results.

Schuh, Michael J.↗

Atomistic Modeling of Semiconductors: Si, C, and 3C-SiC

An ongoing task of the Computational Materials Group (CMG) at the NASA Glenn Research Center is to enhance the role of atomistic simulations based on quantum-approximate methods in the study of new materials and their properties. One of the main goals of the activity continues to be breaching limitations that arise from the natural balance between accuracy, range of application, and computational simplicity. Whether that balance can be maintained while breaking new ground depends on the methods available with a minimum of constraints and limitations for the study of the energetics of arbitrary systems. The main tool used in CMG research, the Bozzolo- Ferrante-Smith (BFS) method for alloys, has no inherent constraint in its formulation, a feature that has allowed for successful research on various topics. In this article, we report on the latest development of the CMG program, namely, the extension and application of the BFS method to compound semiconductors, a departure from our previous research based primarily on metallic alloys.

Bozzolo, Guillermo↗

TPSAS-NF1676L-11678-DND

Cloud vertical profiles retrieved from CALIPSO and CloudSat data provide additional information that greatly improves the estimate of atmospheric radiative heating rates and potentially improves the estimate of the global surface radiation budget. To utilize CALIPSO and CloudSat cloud vertical profiles, we combined their profiles and generated merged profiles. Weused the merged cloud profiles and MODIS derived cloud profiles by the CERES cloud algorithm in irradiance computations. In this study, the improvement of the irradiance by CALIPSO and CloudSat derived merged cloud vertical profiles is evaluated by comparing modeled irradiances with and without merged cloud profiles. Because of lower cloud base heights provided by the merged profiles compared with cloud base heights derived from an empirical formula using the MODIS derived cloud top height and optical thickness, the global annual mean surface longwave downward irradiance increases approximately by 7 Wm-2. While CALIPSO version 3 data have not been incorporated in the irradiance computation and the effect of the new CALIPSO data needs to be assessed, a comparison with surface observations suggests an improvement of surface downward longwave irradiance computations. The surface downward longwave irradiance computed with CALIPSO and CloudSat merged clouds agrees better with surface observations than the irradiance computed with MODIS only does. The bias difference over ARM SGP sites with 2 years of data decreases from -5 W m-2to -2 Wm-2, and it decreases from -9 Wm-2to no bias over polar validation sites. The result suggests that CALIPSO CloudSat derived merged clouds significantly improve surface longwave irradiance estimates once they are combined with MODIS derived cloud properties and with temperature and humidity profiles from reanalysis.

Seiji Kato↗

Near-optimal output feedback regulation of ill-conditioned linear systems

A two-time scale approximation for the linear quadratic optimal output feedback regulator program is examined. Necessary conditions for optimality, as well as an algorithm for computing locally near-optimal gains are derived. If it is assumed that the slow and fast subsystem initial conditions are uniformly distributed, optimal gains for the two-time-scale problem provide a second-order approximation to optimal closed-loop performance in the unperturbed system. This is verified with a numerical example.

Moerder, Daniel D.↗

Thermal stresses in composite tubes using complementary virtual work

This paper addresses the computation of thermally induced stresses in layered, fiber-reinforced composite tubes subjected to a circumferential gradient. The paper focuses on using the principle of complementary virtual work, in conjunction with a Ritz approximation to the stress field, to study the influence on the predicted stresses of including temperature-dependent material properties. Results indicate that the computed values of stress are sensitive to the temperature dependence of the matrix-direction compliance and matrix-direction thermal expansion in the plane of the lamina. There is less sensitivity to the temperature dependence of the other material properties.

Hyer, M. W.↗

Boundary conditions and the simulation of low Mach number flows

The problem of accurately computing low Mach number flows, with the specific intent of studying the interaction of sound waves with incompressible flow structures, such as concentrations of vorticity is considered. This is a multiple time (and/or space) scales problem, leading to various difficulties in the design of numerical methods. Concentration is on one of these difficulties - the development of boundary conditions at artificial boundaries which allow sound waves and vortices to radiate to the far field. Nonlinear model equations are derived based on assumptions about the scaling of the variables. Then these are linearized about a uniform flow and exact boundary conditions are systematically derived using transform methods. Finally, useful approximations to the exact conditions which are valid for small Mach number and small viscosity are computed.

Hagstrom, Thomas↗

Aeroacoustic computation of gust-blade interaction

To better understand and address the challenges faced in computing the acoustics of flow fields, test problems must be considered. In the present study, the sound radiated by the interaction of a flat plate with an oncoming gust containing a two component, mean velocity is computed. The gust has a uniform mean flow in x with Mach number M(infinity) equal to 0.5. The gust's mean velocity in y is of smaller amplitude and is given by: v = 0.1 sin(pi/8(x/M(sub infinity) - t)). This problem has been posed for an upcoming ICASE/LaRC workshop on benchmark problems in computational aeroacoustics. A plate with a length of 30 units in x is used. The plate is assumed to be infinitesimally thin and is centered at the origin. All variables are made dimensionless using the scales specified. Acoustic quantities are obtained by numerically integrating the linearized Euler equations. Integration is performed on the computational domain -100.0 less than or equal to x less than or equal to 100.0, -100.0 less than or equal to y less than or equal to 100.0, using unit length grid spacing in x and in y. An integration scheme is sought which will provide accurate solution to the small quantities of interest at a minimal computational expense. Results indicate that with the given discretization a scheme of minimal fourth order accuracy might be adequate to approximate the waves within the given flow. Thus, a variation of the MacCormack scheme with fourth order accuracy in space and second order accuracy in time was chosen. A scheme with sixth order accuracy in space has also been implemented and results compared with those of the fourth order accurate scheme. To ensure no mass flux, zero normal velocity is assigned at the plate. This condition will induce a discontinuity in the pressure across the plate location. Values for the perturbation pressure p' along the surface of the plate are obtained using a one-sided, third order Taylor expansion, such that p'(sub y) = O. In accordance with the Kutta condition, perturbation pressure at the trailing edge is assigned to zero. In the far field, radiation boundary conditions have been implemented. The effectiveness of the far field conditions are validated by computing in a larger computational domain and comparing the results. Early in time, sound waves begin to radiate from the plate. A Doppler effect is observed. After the initial transients disappear, the strongest waves leave the trailing edge at an approximate 45 degree angle. The intensity pattern of pressure fluctuations shows five lobes (of increasing magnitude with increasing downstream direction) emerging from the plate. Undesirable short wave contaminants are observed in the computed pressure distribution along the plate surface. For a more accurate solution at small scales, a more refined discretization will be required.

Martin, James E.↗

An Interactive MATLAB Program for Fitting Transfer Functions to Frequency Responses

A computer program called FRFit (Frequency Response Fitting) is described for fitting single-input single-output transfer function models to empirical frequency response data. The program is interactive in that the user specifies ``elementary factors'' (gain, delay, pure differentiators and integrators, and first- and second-order zeros and poles) by entering numerical values or moving sliders in a graphical user interface. A nonlinear optimization can then be performed to obtain maximum likelihood estimates of transfer function parameters and uncertainties to provide feedback on the modeling and refine estimates. Several examples are discussed, including the identification of aircraft pitch dynamics from simulation data and data reported in the literature, approximating Theodorsen's function of unsteady aerodynamics, and obtaining a reduced-order model of a computational fluid dynamics code describing the unsteady aerodynamics around an aeroelastic wing. The program has some usefulness as a teaching aid, and can be applied to model structure determination, reduced-order modeling, preliminary analysis, and simple system identification problems. The program was written in MATLAB and is planned for public release through the NASA Software Catalog.

System identification↗