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,297 records · Page 72

Increased depth-diameter ratios in the Medusae Fossae Formation deposits of Mars

Depth to diameter ratios for fresh impact craters on Mars are commonly cited as approximately 0.2 for simple craters and 0.1 for complex craters. Recent computation of depth-diameter ratios in the Amazonis-Memnonia region of Mars indicates that craters within the Medusae Fossae Formation deposits found in this region display greater depth-diameter ratios than expected for both simple and complex craters. Photoclinometric and shadow length techniques have been used to obtain depths of craters within the Amazonis-Memnonia region. The 37 craters in the 2 to 29 km diameter range and displaying fresh impact morphologies were identified in the area of study. This region includes the Amazonian aged upper and middle members of the Medusae Fossae Formation and Noachian aged cratered and hilly units. The Medusae Fossae Formation is characterized by extensive, flat to gently undulating deposits of controversial origin. These deposits appear to vary from friable to indurated. Early analysis of crater degradation in the Medusae Fossae region suggested that simple craters excavated to greater depths than expected based on the general depth-diameter relationships derived for Mars. However, too few craters were available in the initial analysis to estimate the actual depth-diameter ratios within this region. Although the analysis is continuing, we are now beginning to see a convergence towards specific values for the depth-diameter ratio depending on geologic unit.

Barlow, N. G.↗

Structured adaptive grid generation using algebraic methods

The accuracy of the numerical algorithm depends not only on the formal order of approximation but also on the distribution of grid points in the computational domain. Grid adaptation is a procedure which allows optimal grid redistribution as the solution progresses. It offers the prospect of accurate flow field simulations without the use of an excessively timely, computationally expensive, grid. Grid adaptive schemes are divided into two basic categories: differential and algebraic. The differential method is based on a variational approach where a function which contains a measure of grid smoothness, orthogonality and volume variation is minimized by using a variational principle. This approach provided a solid mathematical basis for the adaptive method, but the Euler-Lagrange equations must be solved in addition to the original governing equations. On the other hand, the algebraic method requires much less computational effort, but the grid may not be smooth. The algebraic techniques are based on devising an algorithm where the grid movement is governed by estimates of the local error in the numerical solution. This is achieved by requiring the points in the large error regions to attract other points and points in the low error region to repel other points. The development of a fast, efficient, and robust algebraic adaptive algorithm for structured flow simulation applications is presented. This development is accomplished in a three step process. The first step is to define an adaptive weighting mesh (distribution mesh) on the basis of the equidistribution law applied to the flow field solution. The second, and probably the most crucial step, is to redistribute grid points in the computational domain according to the aforementioned weighting mesh. The third and the last step is to reevaluate the flow property by an appropriate search/interpolate scheme at the new grid locations. The adaptive weighting mesh provides the information on the desired concentration of points to the grid redistribution scheme. The evaluation of the weighting mesh is accomplished by utilizing the weight function representing the solution variation and the equidistribution law. The selection of the weight function plays a key role in grid adaptation. A new weight function utilizing a properly weighted boolean sum of various flowfield characteristics is defined. The redistribution scheme is developed utilizing Non-Uniform Rational B-Splines (NURBS) representation. The application of NURBS representation results in a well distributed smooth grid by maintaining the fidelity of the geometry associated with boundary curves. Several algebraic methods are applied to smooth and/or nearly orthogonalize the grid lines. An elliptic solver is utilized to smooth the grid lines if there are grid crossings. Various computational examples of practical interest are presented to demonstrate the success of these methods.

Yang, Jiann-Cherng↗

Orbital coplanarity in solar-type binary systems: Implications for planetary system formation and detection

The equatorial inclinations of solar-type stars within visual binary systems are computed by combining v sin i measurements with rotational period information, or with expected rotational velocities based upon the age of the star in question. These inclinations are then compared with the orbital inclinations of the systems to test the alignment between the equatorial and orbital planes, and how the tendency for or against coplanarity varies as a function of parameters such as spectral type, separation, eccentricity, etc. The results are extended to planetary systems in order to determine the appropriateness of basing planetary search strategies upon a parent star's equatorial inclination, and to address issues in planetary system formation and evolution, including the stability of planetary orbits within binary systems. During the course of this project new or improved v sin i measurements are made for over 30 solar-type stars within binary systems, and (for the purposes of the study) tentative orbits are computed for thirteen long-period systems. The results suggest that approximate coplanarity between the equatorial and orbital planes exists solar-type binary systems with separations less than 30-40 AU. The coplanarity tendency, as well as this 'critical separation,' is not significantly affected by most of the other parameters studied. The one significant exception occurs with hierarchical multiple systems, where noncoplanarity may exist at relatively small separations. If it is assumed that planetary distances in our solar system are typical, the results suggest there is no reason to expect planets to orbit in planes significantly different from that of the parent star's equator, in turn suggesting that planetary formation models and search strategies dependent upon this assumption are valid from this standpoint. The results also suggest that noncoplanarity between the components of a binary system is not a significant issue in addressing the stability of planetary orbits within the system.

Hale, Alan↗

Use of Density Functional Method to Study Molecular Vibrations

Use of a scale factor reduces the errors in computed vibrational frequencies (errors due to anharmonicities, basis set deficiencies, and approximate treatment of electron correlation). Additional errors come from modelling solids with small molecules. Nevertheless, we hope that the present calculations add to the understanding of surface IR measurements of the heat tiles.

Chong, Delano P.↗

Photoionization and Recombination

Theoretically self-consistent calculations for photoionization and (e + ion) recombination are described. The same eigenfunction expansion for the ion is employed in coupled channel calculations for both processes, thus ensuring consistency between cross sections and rates. The theoretical treatment of (e + ion) recombination subsumes both the non-resonant recombination ("radiative recombination"), and the resonant recombination ("di-electronic recombination") processes in a unified scheme. In addition to the total, unified recombination rates, level-specific recombination rates and photoionization cross sections are obtained for a large number of atomic levels. Both relativistic Breit-Pauli, and non-relativistic LS coupling, calculations are carried out in the close coupling approximation using the R-matrix method. Although the calculations are computationally intensive, they yield nearly all photoionization and recombination parameters needed for astrophysical photoionization models with higher precision than hitherto possible, estimated at about 10-20% from comparison with experimentally available data (including experimentally derived DR rates). Results are electronically available for over 40 atoms and ions. Photoionization and recombination of He-, and Li-like C and Fe are described for X-ray modeling. The unified method yields total and complete (e+ion) recombination rate coefficients, that can not otherwise be obtained theoretically or experimentally.

Nahar, Sultana N.↗

Employing Sensitivity Derivatives to Estimate Uncertainty Propagation in CFD

Two methods that exploit the availability of sensitivity derivatives are successfully employed to predict uncertainty propagation through Computational Fluid Dynamics (CFD) code for an inviscid airfoil problem. An approximate statistical second-moment method and a Sensitivity Derivative Enhanced Monte Carlo (SDEMC) method are successfully demonstrated on a two-dimensional problem. First- and second-order sensitivity derivatives of code output with respect to code input are obtained through an efficient incremental iterative approach. Given uncertainties in statistically independent, random, normally distributed flow parameters (input variables); these sensitivity derivatives enable one to formulate first- and second-order Taylor Series approximations for the mean and variance of CFD output quantities. Additionally, incorporation of the first-order sensitivity derivatives into the data reduction phase of a conventional Monte Carlo (MC) simulation allows for improved accuracy in determining the first moment of the CFD output. Both methods are compared to results generated using a conventional MC method. The methods that exploit the availability of sensitivity derivatives are found to be valid when considering small deviations from input mean values.

Putko, Michele M.↗

A First and Second Order Moment Approach to Probabilistic Control Synthesis

This paper presents a robust control design methodology based on the estimation of the first two order moments of the random variables and processes that describe the controlled response. Synthesis is performed by solving an multi-objective optimization problem where stability and performance requirements in time- and frequency domains are integrated. The use of the first two order moments allows for the efficient estimation of the cost function thus for a faster synthesis algorithm. While reliability requirements are taken into account by using bounds to failure probabilities, requirements related to undesirable variability are implemented by quantifying the concentration of the random outcome about a deterministic target. The Hammersley Sequence Sampling and the First- and Second-Moment- Second-Order approximations are used to estimate the moments, whose accuracy and associated computational complexity are compared numerically. Examples using output-feedback and full-state feedback with state estimation are used to demonstrate the ideas proposed.

Crespo, Luis G.↗

Probabilistic Reasoning for Plan Robustness

A planning system must reason about the uncertainty of continuous variables in order to accurately project the possible system state over time. A method is devised for directly reasoning about the uncertainty in continuous activity duration and resource usage for planning problems. By representing random variables as parametric distributions, computing projected system state can be simplified in some cases. Common approximation and novel methods are compared for over-constrained and lightly constrained domains. The system compares a few common approximation methods for an iterative repair planner. Results show improvements in robustness over the conventional non-probabilistic representation by reducing the number of constraint violations witnessed by execution. The improvement is more significant for larger problems and problems with higher resource subscription levels but diminishes as the system is allowed to accept higher risk levels.

autonomous planning↗

Modeling Optical Properties of Mineral Aerosol Particles by Using Nonsymmetric Hexahedra

We explore the use of nonsymmetric geometries to simulate the single-scattering properties of airborne dust particles with complicated morphologies. Specifically, the shapes of irregular dust particles are assumed to be nonsymmetric hexahedra defined by using the Monte Carlo method. A combination of the discrete dipole approximation method and an improved geometric optics method is employed to compute the single-scattering properties of dust particles for size parameters ranging from 0.5 to 3000. The primary optical effect of eliminating the geometric symmetry of regular hexahedra is to smooth the scattering features in the phase function and to decrease the backscatter. The optical properties of the nonsymmetric hexahedra are used to mimic the laboratory measurements. It is demonstrated that a relatively close agreement can be achieved by using only one shape of nonsymmetric hexahedra. The agreement between the theoretical results and their measurement counterparts can be further improved by using a mixture of nonsymmetric hexahedra. It is also shown that the hexahedron model is much more appropriate than the "equivalent sphere" model for simulating the optical properties of dust particles, particularly, in the case of the elements of the phase matrix that associated with the polarization state of scattered light.

Bi, Lei↗

LES of a Jet Excited by the Localized Arc Filament Plasma Actuators

The fluid dynamics of a high-speed jet are governed by the instability waves that form in the free-shear boundary layer of the jet. Jet excitation manipulates the growth and saturation of particular instability waves to control the unsteady flow structures that characterize the energy cascade in the jet.The results may include jet noise mitigation or a reduction in the infrared signature of the jet. The Localized Arc Filament Plasma Actuators (LAFPA) have demonstrated the ability to excite a high-speed jets in laboratory experiments. Extending and optimizing this excitation technology, however, is a complex process that will require many tests and trials. Computational simulations can play an important role in understanding and optimizing this actuator technology for real-world applications. Previous research has focused on developing a suitable actuator model and coupling it with the appropriate computational fluid dynamics (CFD) methods using two-dimensional spatial flow approximations. This work is now extended to three-dimensions (3-D) in space. The actuator model is adapted to a series of discrete actuators and a 3-D LES simulation of an excited jet is run. The results are used to study the fluid dynamics near the actuator and in the jet plume.

Brown, Clifford A.↗

Elastic Model Transitions Using Quadratic Inequality Constrained Least Squares

A technique is presented for initializing multiple discrete finite element model (FEM) mode sets for certain types of flight dynamics formulations that rely on superposition of orthogonal modes for modeling the elastic response. Such approaches are commonly used for modeling launch vehicle dynamics, and challenges arise due to the rapidly time-varying nature of the rigid-body and elastic characteristics. By way of an energy argument, a quadratic inequality constrained least squares (LSQI) algorithm is employed to e ect a smooth transition from one set of FEM eigenvectors to another with no requirement that the models be of similar dimension or that the eigenvectors be correlated in any particular way. The physically unrealistic and controversial method of eigenvector interpolation is completely avoided, and the discrete solution approximates that of the continuously varying system. The real-time computational burden is shown to be negligible due to convenient features of the solution method. Simulation results are presented, and applications to staging and other discontinuous mass changes are discussed

Orr, Jeb S.↗

Interpolation Method Needed for Numerical Uncertainty

Using Computational Fluid Dynamics (CFD) to predict a flow field is an approximation to the exact problem and uncertainties exist. There is a method to approximate the errors in CFD via Richardson's Extrapolation. This method is based off of progressive grid refinement. To estimate the errors, the analyst must interpolate between at least three grids. This paper describes a study to find an appropriate interpolation scheme that can be used in Richardson's extrapolation or other uncertainty method to approximate errors.

Groves, Curtis E.↗

Passive Rocket Diffuser Theory: A Re-Examination of Minimum Second Throat Size

Second-throat diffusers serve to isolate rocket engines from the effects of ambient back pressure during testing without using active control systems. Among the most critical design parameters is the relative area of the diffuser throat to that of the nozzle throat. A smaller second throat is generally desirable because it decreases the stagnation-to-ambient pressure ratio the diffuser requires for nominal operation. There is a limit, however. Below a certain size, the second throat can cause pressure buildup within the diffuser and prevent it from reaching the start condition that protects the nozzle from side-load damage. This paper presents a method for improved estimation of the minimum second throat area which enables diffuser start. The new 3-zone model uses traditional quasi-one-dimensional compressible flow theory to approximate the structure of two distinct diffuser flow fields observed in Computational Fluid Dynamics (CFD) simulations and combines them to provide a less-conservative estimate of the second throat size limit. It is unique among second throat sizing methods in that it accounts for all major conical nozzle and second throat diffuser design parameters within its limits of application. The performance of the 3-zone method is compared to the historical normal shock and force balance methods, and verified against a large number of CFD simulations at specific heat ratios of 1.4 and 1.25. Validation is left as future work, and the model is currently intended to function only as a first-order design tool.

Jones, Daniel R.↗

High-Spatial-Resolution OH and CH2O PLIF Visualization in a Dual-Mode Scramjet Combustor

A high-spatial-resolution planar laser-induced fluorescence (PLIF) imaging system was constructed and used to image a cavity-stabilized, premixed ethylene-air flame. The flame was created within a continuous flow, electrically-heated supersonic combustion facility consisting of a Mach 2 nozzle, an isolator with flush-wall fuel injectors, a combustor with a cavity flameholder of height 9 mm and optical access, and an extender. Tests were conducted at total temperature 1200 K, total pressure 300 kPa, equivalence ratio near 0.4 in the combustor, and Mach number near 0.6 in the combustor. A frequency-doubled Nd:YAG laser pumped a dye laser, which produced light at 283.55 nm. The beam was shaped into a light sheet with full width half-maximum 25 microns, which illuminated a streamwise plane that bisected the cavity. An intensified camera system imaged OH in this plane with a square 6.67 mm field of view and in-plane resolution 39 microns. Images were taken between the backward-facing step and 120 mm downstream of the step. OH structures as small as 110 microns were observed. CH2O was excited using 352.48 nm light; the smallest observed CH2O structures were approximately 200 microns wide. Approximately 15,000 images per species were processed and used to compute composite images.

Geipel, Clayton M.↗

Comparison of EKF and UKF for Robotic Missions to Mars

The batched Kalman lter works well for a wide variety of orbit determination scenarios but requires linearized approximations to the underlying dynamical system. Furthermore, implementation depends on the computation of partial derivatives for all measurements and parameters used in the lter, a non-trivial task even when analytical derivations are possible. To address these limitations, many extensions and modi - cations to the linear Kalman lter have been proposed over the years since Kalman's initial publication. Our investigation examines two popular variants of the Kalman lter, namely the Extended Kalman Filter (EKF) and Unscented Kalman Filter (UKF), within the context of spacecraft missions to Mars. The EKF relies on the same mathematical basis as the classic version of the lter, however it sequentially updates the linearization as the ltering process is conducted. In contrast, the UKF is built upon the unscented transform which samples the nearby solution space and numerically propagates perturbed trajectories as part of the update / estimation process. We assess the performance of the ltering algorithms under a variety of mission scenarios with the goal of establishing decision criteria for which implementation to use under di ering circumstances.

Ely, Todd↗

The Calculation of Diffusion Properties for Open Shell Molecules

- Navier-Stokes requires as input various collision integrals for gases - Diffusion, - Viscosity - Thermal conductivity - Thermal diffusion ... - Levin, Stallcop and Partridge state-of-the-art - Non-Relativistic Born-Oppenheimer Approximation (freeze nuclei first, let them move second) - Solve electron problem: Compute realistic Potential Curves (PC) from first principles, - Solve nuclei problem: Compute scattering properties for each PC using sophisticated semi-classical method - Compute cross sections - Compute appropriate thermal collision integrals

STMD↗

Report on a Workshop for Heliophysics Infrastructure

We report observations and findings from a three-day virtual workshop held May 17-19, 2021 which examined the current heliophysics research infrastructure to determine which elements were most utilized, what gaps exist in these elements between current utility and desired capability and, from a user standpoint, what a future state for the infrastructure might look like. Approximately 40 subject matter experts (SMEs) with backgrounds in heliophysics research, computer science and research infrastructure were gathered to consider this topic.

Brian A Thomas↗