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,639 records · Page 91

Summary of Ice Shape Geometric Fidelity Studies on an Iced Swept Wing

Understanding the aerodynamic impact of swept-wing ice accretions is a crucial component of the design of modern aircraft. Computer-simulation tools are commonly used to approximate ice shapes, so the necessary level of detail or fidelity of those simulated ice shapes must be understood relative to high-fidelity representations of the ice. Previous tests were performed in the NASA Icing Research Tunnel to acquire high-fidelity ice shapes. Some of those ice shapes are based on aircraft certification requirements. From this database, full-span artificial ice shapes were designed and manufactured for both an 8.9%-scale and 13.3%-scale semispan wing model of the CRM65 which has been established as the full-scale baseline for this swept-wing project. These models were tested in the Walter H. Beech wind tunnel at Wichita State University and at the ONERA (Office national d'etudes et de recherches aerospatiales) F1 facility, respectively. The data collected in the Wichita State University wind tunnel provided a low-Reynolds number baseline study while the pressurized F1 facility produced data over a wide range of Reynolds and Mach numbers with the highest Reynolds number studied being approximately Re = 11.9 by 10 (sup 6). Three different fidelity representations were created based on three different icing conditions. Lower-fidelity ice shapes were created by lofting a smooth ice shape between cross-section cuts of the high-fidelity ice shape. Grit roughness was attached to this smooth ice shape as another fidelity variant. The data indicates that the geometric fidelity of the ice shapes resulted in significant differences in lift and drag. These results were similar at both facilities over the wide range of test conditions utilized.

Aerodynamics↗

Autonomous Spacecraft Attitude Control Using Deep Reinforcement Learning

While machine learning and spacecraft autonomy continue to gain research interest, significant work remains to be done in efficiently applying modern machine learning techniques to problems in space ight. This study presents a framework for deriving a discrete neural spacecraft attitude controller using reinforcement learning, a paradigm of machine learning, without the need for high-performance computing. The developed attitude controller is an approximately time-optimal solution to a highly constrained control problem, able to achieve well above industry-standard pointing accuracies. Control examples are also presented of the agent performing large-angle spacecraft slews in the developed simulation environment and future extensions of this work are discussed.

ATAP↗

Reduced-Order Aerodynamic Modeling Based on CFD Frequency Responses from Multisine Inputs

A system identification analysis was performed to determine a reduced-order model (ROM) of a computational fluid dynamics (CFD) solver in support of linear aeroservoelastic model development and feedback control design. The approach was applied to the FUN3D code for the half-span wind tunnel test article used in the NASA-Boeing collaboration called the Integrated Adaptive Wing Technology Maturation (IAWTM) project. In a transonic flow condition, multiple inputs (11 structural mode displacements and 3 control surface deflections) were simultaneously excited with orthogonal phase-optimized multisines while multiple outputs (the corresponding 14 generalized aerodynamic forces) were recorded. From these recorded times series, the matrix of frequency responses was computed and subsequently fit using rational function approximations (RFAs). It was found that the entire (14 x 14) matrix of frequency responses could be determined from a single CFD run and that results generally followed trends predicted using other methods. Differences were attributed to the modeling fidelity and nonlinearities from structural mode and control surface interactions at higher reduced frequencies. More accurate fits of the RFAs to the frequency response data were obtained by making two CFD runs, one with only structural mode excitations and one with only control surface excitations, which reduced the degree of nonlinearity in the modeling data.

Aeroservoelasticity↗

On Characterization of Flow Disturbances in Arc-Jet Testing

This paper reports computational simulations and analysis of flow characterization tests in high enthalpy arc-jet facilities at NASA Ames Research Center. Flow disturbances and their characterization are presented through case studies, and their implications for testing are discussed. The case studies include free-jet test configurations from three different arc heaters and nozzles: the 60-MW IHF 30-inch conical nozzle, the 10-MW TP3 15-inch conical nozzle, and the 20-MW AHF 12-inch conical nozzle. For all cases, test articles are placed in the jet exiting the conical nozzle, and the existence of a flow disturbance is confirmed through flow survey data using pitot pressure and heat flux probes and accompanying analysis. The paper focuses on cases where the effects of the disturbances on test article surface quantities are important. The present analysis comprises computational fluid dynamics simulations of the nonequilibrium flowfield in the facility nozzles and test box, including the models tested. Comparisons of computations with the experimental measurements are presented. The computations that reproduce the probe sweep data approximately are essential to interpret the arc-jet test data accurately, while providing insights into several observed calorimeter anomalies caused by flow disturbances.

Arc-jet testing↗

On Characterization of Flow Disturbances in Arc-Jet Testing

This paper reports computational simulations and analysis of flow characterization tests in high enthalpy arc-jet facilities at NASA Ames Research Center. Flow disturbances and their characterization are presented through case studies, and their implications for testing are discussed. The case studies include free-jet test configurations from three different arc heaters and nozzles: the 60-MW IHF 30-inch conical nozzle, the 10-MW TP3 15-inch conical nozzle, and the 20-MW AHF 12-inch conical nozzle. For all cases, test articles are placed in the jet exiting the conical nozzle, and the existence of a flow disturbance is confirmed through flow survey data using pitot pressure and heat flux probes and accompanying analysis. The paper focuses on cases where the effects of the disturbances on test article surface quantities are important. The present analysis comprises computational fluid dynamics simulations of the nonequilibrium flowfield in the facility nozzles and test box, including the models tested. Comparisons of computations with the experimental measurements are presented. The computations that reproduce the probe sweep data approximately are essential to interpret the arc-jet test data accurately, while providing insights into several observed calorimeter anomalies caused by flow disturbances.

Arc-jet testing↗

Kriging-Based Inverse Uncertainty Quantification of BISON Fission Gas Release Model

Nuclear reactor fuel performance analysis studies the thermo-mechanical behavior of fuel rods and verify their compliance with safety criteria under both normal operation and accidental conditions. One of the major concerns in fuel performance analysis is the behavior of the fission gases xenon and krypton in uranium dioxide fuel, which significantly affect the thermo-mechanical performance of the nuclear fuel rods employed in current LWRs. Previously, the modeling of fission gas release (FGR) in nuclear fuel performance codes have been relying on empirical models which are not applicable beyond their range of calibration. Recently a more efficient and flexible physics-based FGR and swelling model that can describe a wider range of reactor operation conditions have been implemented in fuel performance code BISON. In this work, we performed Bayesian calibration to inversely quantify the uncertainties associated with the 5 tuning parameters in BISON FGR model, based on the time-dependent FGR measurement data from Riso-AN3 benchmark. It takes around 1 hour and 40 minutes for BISON to simulate the Riso-AN3 benchmark with a moderate mesh size using 32 processors. As hundreds of thousands of BISON simulations are needed during Markov Chain Monte Carlo (MCMC) sampling, running BISON directly is not a viable solution. Therefore we propose to use Kriging-based metamodels to replace BISON during MCMC sampling. Surrogate model is an approximation of the input/output relation of a computer code/model. It is also called metamodel, response surface, emulator. Surrogate models usually take much less computational time than the original model (like BISON) while maintaining the input/output relation of the original model to a desirable accuracy. Constructing the metamodels normally requires small number of executions of the original model. Once validated, metamodels can be used to perform uncertainty and sensitivity analysis, Bayesian calibration, optimization, etc.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Sparse chronology strategy for integrating seasonal energy storage in capacity expansion models

Here, this study develops the sparse chronology method to enhance the representative period framework in capacity expansion models, enabling the effective integration of long-duration energy storage modeling. Traditional representative period methods cannot capture the state of charge of seasonal energy storage systems because they do not establish effective inter-day linkages to connect the state of charge between periods. The sparse chronology approach addresses this limitation by establishing inter-day linkages that allow state of charge to shift inter-seasonally. At the same time, it groups identical representative days into partitions, applying constraints sparsely and implicitly to reduce computational load further. Validation results demonstrate that this method successfully simulates long-duration energy storage patterns, achieving close alignment with a continuous yearly benchmark model, with seasonal trends and state of charge cycles clearly represented. The computational load analysis reveals that the sparse chronology method efficiently applies constraints on maximum and minimum state of charge limits within the representative day framework, eliminating the need for detailed constraints on each individual day. By partitioning representative days and constraining only the start and end of each partition, the method significantly decreases computational requirements. Simulation results show that sparse chronology closely approximates the continuous yearly method's accuracy, even with as few as 20 representative days, achieving correlation values with the benchmark of nearly 0.9 in state of charge plots. Furthermore, it maintains computational efficiency, requiring only 4 % of the solver time compared to the continuous yearly method with 20 representative days. This approach allows capacity expansion models to incorporate long-duration energy storage with high temporal, spatial, and technological resolution, enabling more detailed modeling for large-scale power systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Design and Computational/Experimental Analysis of Low Sonic Boom Configurations

Recent studies have shown that inviscid CFD codes combined with a planar extrapolation method give accurate sonic boom pressure signatures at distances greater than one body length from supersonic configurations if either adapted grids swept at the approximate Mach angle or very dense non-adapted grids are used. The validation of CFD for computing sonic boom pressure signatures provided the confidence needed to undertake the design of new supersonic transport configurations with low sonic boom characteristics. An aircraft synthesis code in combination with CFD and an extrapolation method were used to close the design. The principal configuration of this study is designated LBWT (Low Boom Wing Tail) and has a highly swept cranked arrow wing with conventional tails, and was designed to accommodate either 3 or 4 engines. The complete configuration including nacelles and boundary layer diverters was evaluated using the AIRPLANE code. This computer program solves the Euler equations on an unstructured tetrahedral mesh. Computations and wind tunnel data for the LBWT and two other low boom configurations designed at NASA Ames Research Center are presented. The two additional configurations are included to provide a basis for comparing the performance and sonic boom level of the LBWT with contemporary low boom designs and to give a broader experiment/CFD correlation study. The computational pressure signatures for the three configurations are contrasted with on-ground-track near-field experimental data from the NASA Ames 9x7 Foot Supersonic Wind Tunnel. Computed pressure signatures for the LBWT are also compared with experiment at approximately 15 degrees off ground track.

Cliff, Susan E.↗

Low-Rank Solution Methods for Stochastic Eigenvalue Problems

In this report, we study efficient solution methods for stochastic eigenvalue problems arising from discretization of self-adjoint partial differential equations with random data, where the underlying operators depend linearly on the random parameters. With the stochastic Galerkin approach, the solutions are represented as generalized polynomial chaos expansions. When these solutions can be approximated well by low-rank objects, we introduce a low-rank variant of the inverse subspace iteration algorithm for computing one or several minimal eigenvalues and corresponding eigenvectors of parameter-dependent matrices. In the algorithm, the iterates are approximated by low-rank matrices, which leads to significant cost savings. Furthermore, the algorithm is tested on two benchmark problems, a stochastic diffusion problem with some poorly separated eigenvalues, and an operator derived from a discrete stochastic Stokes problem whose minimal eigenvalue is related to the inf-sup stability constant. Numerical experiments show that the low-rank algorithm produces accurate solutions compared to the Monte Carlo method, and it uses much less computational time than the original algorithm without low-rank approximation.

97 MATHEMATICS AND COMPUTING↗

Non-linear analysis of symmetric structures with unsymmetric boundary conditions

Two computational techniques are developed analytically to generate global approximation vectors for use in the nonlinear large-deflection analysis of symmetric structures with asymmetric boundary conditions. One method is based on linear combinations of FEM-generated symmetric and antisymmetric vector components, while the other employs a preconditioned-conjugate-gradient technique. The derivation of the procedures is explained in detail, and their effectiveness is evaluated in sample computations for an elliptic toroid and a cylindrical panel; the results are presented in tables and graphs and characterized.

Noor, Ahmed K.↗

A Wave Energy Converter Design Load Case Study

This article presents an example by which design loads for a wave energy converter (WEC) might be estimated through the various stages of the WEC design process. Unlike previous studies, this study considers structural loads, for which, an accurate assessment is crucial to the optimization and survival of a WEC. Three levels of computational fidelity are considered. The first set of design load approximations are made using a potential flow frequency-domain boundary-element method with generalized body modes. The second set of design load approximations are made using a modified version of the linear-based time-domain code WEC-Sim. The final set of design load simulations are realized using computational fluid dynamics coupled with finite element analysis to evaluate the WEC’s loads in response to both regular and focused waves. This study demonstrates an efficient framework for evaluating loads through each of the design stages. In comparison with experimental and high-fidelity simulation results, the linear-based methods can roughly approximate the design loads and the sea states at which they occur. The high-fidelity simulations for regular wave responses correspond well with experimental data and appear to provide reliable design load data. The high-fidelity simulations of focused waves, however, result in highly nonlinear interactions that are not predicted by the linear-based most-likely extreme response design load method.

16 TIDAL AND WAVE POWER↗

Stationary Bragg reflection of laser light in inhomogeneous absorbing plasmas inside inertial confinement fusion Hohlraums

Laser-produced plasma in inertial confinement fusion (ICF) Hohlraums are marked with density non-uniformity whose length scale can go down to micrometers. This scale is of the order of the laser wavelength. The WKB approximation, which is classically used in radiation-hydrodynamic codes to compute the laser trajectory, cannot correctly take into account such small-scale inhomogeneity of the plasma. Going beyond this approximation, we predict a novel mechanism for the laser reflection. Here we show that an electromagnetic plane wave with wave number k resonates with the k B =2 k Fourier component of a multimode perturbation of the background density and generates a reflected wave. It is the first time that this reflection is considered for stationary inhomogeneous ICF plasmas, and the energy absorption is taken into account. This mechanism, which is a form of Bragg reflection, can occur away from the critical surface and generate a drift of the location of the laser absorption. Furthermore, this absorption will be periodically modulated with a k B wave number. The stationary Bragg reflection can explain ongoing discrepancies between experimental and numerical data about laser trajectory and absorption in ICF Hohlraums.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of equilibrium air computer programs suitable for numerical computation using time-dependent or shock-capturing methods

Computer programs were developed which compute the thermodynamic properties of equilibrium air for use in either the time-dependent or shock-capturing computational methods. For the time-dependent method, tne NASA-ARC RGAS computer program was modified to allow internal energy and density to be used as the independent variables. In addition, simplified-curve fits for p = p(e,rho), a = a(e,rho), and T = T(p,rho) were devised to reduce computer time. For the shock-capturing method a simplified curve fit for h = h(p,rho) was made. These approximate curve fits may be particularly useful when employed on advanced computers such as the Illiac 4 or the CDC Star since they avoid the cumbersome table-lookup feature of the RGAS program.

Tannehill, J. C.↗

Uncertainty Quantification using Deep Ensembles for Decision Making in Cyber-Physical-Human Systems

In this paper and its companion, Differential Equation Approximation Using Gradient-Boosted Quantile Regression, Robison et al., we examine an approach to quantifying model uncertainty with the aim of increasing the trustworthiness of computational models in human-machine interactions. In Differential Equation Approximation Using Gradient-Boosted Quantile Regression, we focus on gradient-boosted decision trees, while in this one, we give more details about deep ensembles. Uncertainty quantification is crucial for building trustworthy autonomous decision-making agents in human-machine teams. There are two types of uncertainties: aleatoric and epistemic. The former is related to the inherent stochasticity (noise) of the process, whereas the latter is associated with the lack of knowledge or representation capability of models, such as neural networks. By lack of knowledge, we mean the model’s inability to accurately predict outputs for all possible inputs. The aleatory uncertainty can be estimated fairly easily with, for example, filters, whereas epistemic uncertainty is challenging to compute. This paper uses deep ensembles to quantify both aleatory and epistemic uncertainty. It can act as an uncertainty-aware surrogate transition model for decision-making frameworks. "Uncertainty-aware" means that the surrogate transition model should make predictions along with confidence in those predictions. In the context of decision-making, the transition models are ordinary differential equations (ODEs). Since ODEs can be simulated to make one-step or multi-step predictions, a good surrogate model for them should perform reasonably well in both modes. In a multi-step approach, the trajectory sampling method TS∞ was used to propagate uncertainty over multiple steps. The cartpole dynamical system was selected to demonstrate the ability of deep ensembles as good surrogate transition models for decision-making frameworks. The deep ensembles modeled the dynamics of cartpole ODEs and made uncertainty-aware predictions in single-step and multi-step transition modes.

CPH systems↗

Efficient Monte Carlo event generation for neutrino-nucleus exclusive cross sections

Modern neutrino-nucleus cross section computations need to incorporate sophisticated nuclear models to achieve greater predictive precision. However, the computational complexity of these advanced models often limits their practicality for experimental analyses. To address this challenge, we introduce a new Monte Carlo method utilizing normalizing flows to generate surrogate cross sections that closely approximate those of the original model while significantly reducing computational overhead. As a case study, we built a Monte Carlo event generator for the neutrino-nucleus cross section model developed by the Ghent group. This model employs a Hartree-Fock procedure to establish a quantum mechanical framework in which both the bound and scattering nucleon states are solutions to the mean-field nuclear potential. The surrogate cross sections generated by our method demonstrate excellent accuracy with a relative effective sample size of more than 98.4%, providing a computationally efficient alternative to traditional Monte Carlo sampling methods for differential cross sections.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

An X-Ray Analysis Database of Photoionization Cross Sections Including Variable Ionization

Results of research efforts in the following areas are discussed: review of the major theoretical and experimental data of subshell photoionization cross sections and ionization edges of atomic ions to assess the accuracy of the data, and to compile the most reliable of these data in our own database; detailed atomic physics calculations to complement the database for all ions of 17 cosmically abundant elements; reconciling the data from various sources and our own calculations; and fitting cross sections with functional approximations and incorporating these functions into a compact computer code.Also, efforts included adapting an ionization equilibrium code, tabulating results, and incorporating them into the overall program and testing the code (both ionization equilibrium and opacity codes) with existing observational data. The background and scientific applications of this work are discussed. Atomic physics cross section models and calculations are described. Calculation results are compared with available experimental data and other theoretical data. The functional approximations used for fitting cross sections are outlined and applications of the database are discussed.

Wang, Ping↗

Some numerical and physical aspects of unsteady Navier-Stokes computations over airfoils using dynamic meshes

An upwind-biased implicit approximate factorization algorithm is applied to several unsteady flows on dynamic meshes. The thin-layer form of the compressible Navier-Stokes equations is used to solve both laminar and turbulent flows over airfoils pitching about the quarter chord. Numerical aspects of the solutions are investigated, including grid and time step effects. Two methods for determining fluxes - flux-vector splitting and flux-difference splitting - are compared. Flux-difference splitting predicts results more accurately than flux-vector splitting on a coarse mesh, but both methods agree on a fine mesh. Physical aspects of the computations are also examined. An equilibrium turbulent boundary layer model computes generally better unsteady results in comparison with experiment than a nonequilibrium model for the transonic case analyzed. Also, the size and location of the primary shed vortex for an airfoil pitching up at a constant rate is calculated in good agreement with experiment for two pitch rates.

Rumsey, Christopher L.↗