Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Analysis of the weighted shifted boundary method for the Poisson and Stokes problems

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods, and relies on reformulating the original boundary value problem over a surrogate (approximate) computational domain. Accuracy is maintained by properly shifting the location and values of the boundary conditions. This avoids integration over cut cells and the associated implementation issues. Recently, the Weighted SBM (WSBM) was proposed for the Navier-Stokes equations with free surfaces and the Stokes flow with moving boundaries. The attribute “weighted” in the name WSBM stems from the fact that its variational form is weighted with the elemental volume fraction of active fluid. The motivation for the development of the WSBM was the preservation of the volume of active fluid to a higher degree of accuracy, which in turn resulted in improved stability and robustness characteristics in moving-boundary, time-dependent simulations. In this article, we present the numerical analysis of the WSBM formulations for the Poisson and Stokes problems. We give mathematical conditions under which the bilinear forms defining the discrete variational formulations are uniformly coercive (Poisson problem) or inf-sup stable (Stokes problem). By these results, stability and optimal convergence is proven in the natural norm; L2-error estimates can also be derived.

Approximate domain boundaries↗
DOE OSTI2025

Dehydroxylation kinetics of kaolinite and montmorillonite examined using isoconversional methods

The use of calcined clays as supplementary cementitious materials (SCMs) in concrete is a promising strategy towards decarbonizing the cement and concrete industry. This is especially relevant considering the ever-increasing demand for concrete. Comprehensive understanding of the kinetics of calcination is essential towards maximizing the potential reactivity of clay minerals while ensuring energy efficiency. In this study, the kinetics of the dehydroxylation of kaolinite and montmorillonite are investigated under non-isothermal conditions at constant heating rate. Activation energies ( E a ) are determined via Friedman differential and advanced Vyazovkin incremental methods over the isoconversional range; these are devoid of computational approximations, thus allowing kinetic analysis without assuming a specific reaction model. Kinetic equations—in the differential form as well as a combination of differential and integral forms are compared against the experimentally determined reaction models to identify the most probable dehydroxylation mechanism for kaolinite and montmorillonite. A reaction order mechanism is established for dehydroxylation of kaolinite, while montmorillonite is noted to undergo dehydroxylation via a single-step reversible diffusion-controlled process. Kinetic triplet—comprising activation energy, reaction model and pre-exponential factor—is used to predict isothermal calcination conditions, which is further verified using analytical techniques. Heat release rates of clay-portlandite blends from isothermal calorimetry are used within a thermodynamic framework to quantify reactivity of the calcined clays. Here, the study demonstrates a general approach based on isoconversional methods to predict calcination conditions for different clays that can be used in efficient and optimized production of blended cements or SCMs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

GPR_calculator: An on-the-fly surrogate model to accelerate massive nudged elastic band calculations

We present GPR_calculator, a package based on Python and C++ programming languages to build an on-the-fly surrogate model using Gaussian Process Regression (GPR) to approximate computationally expensive electronic structure calculations. The key idea is to dynamically train a GPR model during the simulation that can accurately predict energies and forces with uncertainty quantification. When the uncertainty is high, the costly electronic structure calculation is performed to obtain the ground truth data, which is then used to update the GPR model. To illustrate the effectiveness of GPR_calculator, we demonstrate its application in Nudged Elastic Band (NEB) simulations of surface diffusion and reactions, achieving 3-10 times acceleration compared to pure ab initio calculations. The source code is available at https://github.com/MaterSim/GPR_calculator.

Gaussian process regression↗
DOE OSTI2022

How to see hidden patterns in metamaterials with interpretable machine learning

Machine learning models can assist with metamaterials design by approximating computationally expensive simulators or solving inverse design problems. However, past work has usually relied on black box deep neural networks, whose reasoning processes are opaque and require enormous datasets that are expensive to obtain. Here, in this work, we develop two novel machine learning approaches to metamaterials discovery that have neither of these disadvantages. These approaches, called shape-frequency features and unit-cell templates, can discover 2D metamaterials with user-specified frequency band gaps. Our approaches provide logical rule-based conditions on metamaterial unit-cells that allow for interpretable reasoning processes, and generalize well across design spaces of different resolutions. The templates also provide design flexibility where users can almost freely design the fine resolution features of a unit-cell without affecting the user’s desired band gap.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A New Hybrid Quantum-Classical Algorithm for Solving the Unit Commitment Problem

Solving problems related to planning and operations of large-scale power systems is challenging on classical computers due to their inherent nature as mixed-integer and nonlinear problems. Quantum computing provides new avenues to approach these problems. We develop a hybrid quantum-classical algorithm for the Unit Commitment (UC) problem in power systems which aims at minimizing the total cost while optimally allocating generating units to meet the hourly demand of the power loads. The hybrid algorithm combines a variational quantum algorithm (VQA) with a classical Benders-type heuristic. The resulting algorithm computes approximate solutions to UC in three stages: i) a collection of UC vectors capable meeting the power demand with lowest possible operating costs is generated based on VQA; ii) a classical sequential least squares programming (SLSQP) routine is leveraged to find the optimal power level corresponding to a predetermined number of candidate vectors; iii) in the last stage, the approximate solution of UC along with generating units power level combination is given. To demonstrate the effectiveness of the presented method, three different systems with 3 generating units, 10 generating units, and 26 generating units were tested for different time periods. In addition, convergence of the hybrid quantum-classical algorithm for select time periods is proven out on IonQ's Forte system.

Aboumrad, Willie [IonQ, Inc]↗
DOE OSTI2022

Fusion RF Modeling Machine Learning (FusionML_RF) v1.0

FusionML_RF consists of multiple codes and trained machine learning (ML) models that perform low-cost output modeling from the Genray-CQL3D. Three machine learning techniques (multilayer perceptron, random forest, and Gaussian process) provide fast surrogate models for lower hybrid current drive (LHCD) simulations. For example, completing a single GENRAY/CQL3D simulation without radial diffusion of fast electrons requires several minutes of wall-clock time. On the other hand, these ML models achieve ~ms of inference time with high accuracy across the input parameter space. This software collection consists of multiple components. (1) codes that use ML methods and precomputed Genray-CQL3D simulation output to build regression models that enable approximate computations of Genray-CLQ3D outputs from arbitrary but physically meaningful input parameters (surrogate modeling); (2) three trained models created by the team, using a database of 16,000+ GENRAY/CQL3D simulations, to study the performance of ML models for surrogate modeling; (3) codes that load the trained models and simulation data, and then compute mean squared error between the models' predictions and the ground truth of simulation output data. This collection is being made available in conjunction with a scientific publication about the work to promote reusability and provide an artifact of the scientific work.

Bai, Zhe↗
DOE OSTI2025

Tough Errors Are no Match (TEAM): Optimizing the Quantum Compiler for Noise Resilience

This report summarizes research performed under the Tough Errors Are no Match (TEAM) project. The primary focus of TEAM has been to research and develop a compilation toolbox leveraging techniques from quantum characterization and control, probabilistic programming, and approximate computing. Our goal was to develop robust protocols that can be integrated into quantum compilers to optimize and enhance the robustness of noisy computation. Here, we provide a summary of TEAM work focused on characterization and control of quantum systems.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1976

A block iterative finite element algorithm for numerical solution of the steady-state, compressible Navier-Stokes equations

An iterative method for numerically solving the time independent Navier-Stokes equations for viscous compressible flows is presented. The method is based upon partial application of the Gauss-Seidel principle in block form to the systems of nonlinear algebraic equations which arise in construction of finite element (Galerkin) models approximating solutions of fluid dynamic problems. The C deg-cubic element on triangles is employed for function approximation. Computational results for a free shear flow at Re = 1,000 indicate significant achievement of economy in iterative convergence rate over finite element and finite difference models which employ the customary time dependent equations and asymptotic time marching procedure to steady solution. Numerical results are in excellent agreement with those obtained for the same test problem employing time marching finite element and finite difference solution techniques.

Cooke, C. H.↗
NASA NTRS1977

A block iterative finite element algorithm for numerical solution of the steady-state, compressible Navier-Stokes equations

An iterative method for numerically solving the time independent Navier-Stokes equations for viscous compressible flows is presented. The method is based upon partial application of the Gauss-Seidel principle in block form to the systems of the nonlinear algebraic equations which arise in construction of finite element (Galerkin) models approximating solutions of fluid dynamic problems. The continuous cubic element on triangles is employed for function approximation. Computational results for a free shear flow at Re = 1000 indicate significant achievement of economy in iterative convergence rate over finite element and finite difference models which employ the customary time dependent equations and symptotic time marching procedure to steady solution. Numerical results are in excellent agreement with those obtained for the same test problem employing time marching finite element and finite difference solution techniques.

Cooke, C. H.↗
NASA NTRS1980

An alternative approach to the numerical simulation of steady inviscid flow

A numerical procedure for the efficient simulation of steady inviscid flow is described and its utility demonstrated. Application of the surrogate equation technique allows the formulation of stable, fully conservative, type dependent finite difference equations for use in obtaining numerical solutions to systems of first order partial differential equations, such as the steady state Euler equations or their various approximations. Computational results are presented for the full Euler equations and for the transonic disturbance equations. For the latter case, a computational efficiency greater than that obtained by means of the standard perturbation potential approach is indicated.

Johnson, G. M.↗
NASA NTRS1980

Accuracies of three computationally efficient algorithms for computing atmospheric transmittances

Three algorithms for calculating polychromatic atmospheric transmittance functions have been tested using a set of eleven distinct temperature profiles in order to compare transmittance accuracies achievable by the three methods. The comparison of rms errors demonstrates that the iterative method of McMillin and Fleming (1976) is the most accurate of the efficient algorithms currently available for gases with constant mixing ratios; its accuracy approaches that of the spectroscopic parameters and the computational approximations used in the ground-truth line-by-line calculations. The method of Arking et al. (1974), while less accurate, has the advantage of being perfectly general and easily adapted to cases where spectral bandwidths are varied

Mcmillin, L. M.↗
NASA NTRS1980

An alternative approach to the numerical simulation of steady inviscid flow

A numerical procedure for the efficient simulation of steady inviscid flow is described and its utility is demonstrated. The method is uniformly valid for application in the subsonic, transonic and supersonic flow regimes. It does not rely on the introduction of additional assumptions beyond those necessary to obtain the Euler equations from the Navier-Stokes equations, nor does it make use of a time-asymptotic solution of the unsteady equations of motion. Application of the herein-defined surrogate equation technique allows the formulation of stable, fully-conservative, type-dependent finite difference equations for use in obtaining numerical solutions to systems of first-order partial differential equations, such as the steady-state Euler equations or their various approximations. Computational results are presented for the full Euler equations used to simulate rotational subsonic flow and for the transonic small disturbance equations. For the latter case, a computational efficiency greater than that obtained by means of the standard perturbation potential approach is indicated.

Johnson, G. M.↗
NASA NTRS1985

Spatial variation of cosmic rays near the heliospheric current sheet

A quantitative comparison between theoretical predictions and observations of the intensity of galactic cosmic rays near the interplanetary current sheet is reported. Comparison of model calculations is made with a statistical analysis of observations of galactic cosmic rays at Earth and the simultaneous position of the current sheet. An ensemble of different current sheet inclinations is used, in order to make the analysis of the computations approximate the method used to analyses the data.

Jokipii, J. R.↗
NASA NTRS1986

Cosmic rays near the heliospheric current sheet. II - An ensemble approach to comparing theory and observation

A quantitative comparison is carried out between theoretical predictions and observations of the intensity of galactic cosmic rays near the interplanetary current sheet. Model calculations are compared with a statistical analysis of observations of galactic cosmic rays at the earth and the simultaneous position of the current sheet. Since the observations were made over a period of several years, the inclination of the current sheet varied considerably, and comparison with any one model calculation referring to a given inclination would not be appropriate. An ensemble of different current sheet inclinations are used to generate expected values with the model in order to make the analysis of the computations approximate the method used to analyze the data. Agreement is found between theory and observation at energies of the order of a few GeV.

Jokipii, J. R.↗
NASA NTRS1987

On computing eigensolution sensitivity data using free vibration solutions

A simplified method of computing eigensolution sensitivity derivatives in structural dynamics is developed. It is shown that if the elements of stiffness and mass matrices associated with a design variable are homogeneous functions of that design variable, then eigenvalue derivatives can be computed from element strain and kinetic energies. Furthermore, if cross-mode energies are known, eigensolution derivatives of modified systems can be computed approximately using assume mode reanalysis formulation. A ten bar truss example is used to illustrate the present formulations.

Wang, B. P.↗