Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Embedded 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,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

Computational efficiency achieved by time splitting of finite difference operators.

A technique is presented for computing multidimensional time-dependent flow fields that avoids much of the inefficiency typically found in finite difference calculations. The technique initially divides the flow field into regions, each containing a mesh of general quadrilateral cells chosen to provide spatial resolution of the local features of the flow. A finite difference operator of second order accuracy, consisting of a sequence of one-dimensional operators (each operating at near maximum Courant-Friedrich-Lewy number) is then constructed for each region. Numerical results illustrating the technique for inviscid flows about simple bodies that generate shock waves, embedded shock waves, and expansion fans are presented and compared with exact theory.

Maccormack, R. W.↗
NASA NTRS1992

Design for interaction between humans and intelligent systems during real-time fault management

Initial results are reported to provide guidance and assistance for designers of intelligent systems and their human interfaces. The objective is to achieve more effective human-computer interaction (HCI) for real time fault management support systems. Studies of the development of intelligent fault management systems within NASA have resulted in a new perspective of the user. If the user is viewed as one of the subsystems in a heterogeneous, distributed system, system design becomes the design of a flexible architecture for accomplishing system tasks with both human and computer agents. HCI requirements and design should be distinguished from user interface (displays and controls) requirements and design. Effective HCI design for multi-agent systems requires explicit identification of activities and information that support coordination and communication between agents. The effects are characterized of HCI design on overall system design and approaches are identified to addressing HCI requirements in system design. The results include definition of (1) guidance based on information level requirements analysis of HCI, (2) high level requirements for a design methodology that integrates the HCI perspective into system design, and (3) requirements for embedding HCI design tools into intelligent system development environments.

Malin, Jane T.↗
NASA NTRSDate not supplied

Building and Breaking Carbon Composites with REACTER

Carbon composites have become indispensable for aerospace and other high-performance applications, and a detailed picture of their morphology and failure mechanisms remains difficult to obtain through experiment. REACTER is a versatile computational modeling tool for atomistic molecular dynamics designed to model chemical reactions at the speed and length scales of classical force fields.1 In this work, several recent features of REACTER were applied to the creation and subsequent mechanical testing of two classes of carbon composites, carbon-fiber reinforced polymers (CFRP) and carbon nanotube (CNT) composites. Carbon fiber core morphologies were created by the method of Desai et al.,2 but using the advanced reaction constraints framework of REACTER, their proposed multistep procedure was reduced to a single uninterrupted molecular dynamics simulation. The carbon fiber filler was embedded into a polymer matrix by simulated in situ polymerization of several thermosetting resins, including bismaleimide and polyarylacetylene, to obtain the final CFRP model. To generate the second class of carbon composite, CNT networks were grown dynamically using the new ‘create atoms’ feature of REACTER, and similarly infiltrated with resin to obtain CNT composites. The resulting models were compared directly to experiment using simulated high resolution transmission electron microscopy and x-ray diffraction. Failure mechanisms were elucidated by simulating mechanically induced bond breaking, as characterized by third-order DFT-based tight-binding (DFTB3) simulations, via a reaction constraint on the total potential energy of the involved atoms.

Molecular Dynamics↗
DOE OSTI2024

Riemannian Optimization Applied to AC Optimal Power Flow

The nonlinear, nonconvex AC optimal power flow problem is of growing importance as the nature of the power grid evolves. This problem can be difficult to solve for interior point methods. However, the advent of optimization algorithms over smooth Riemannian manifolds presents an alternative approach. The nonlinear, nonconvex constraints in the AC power flow problem form an embedded submanifold of Euclidean space. In this paper, the authors explore the performance of Riemannian optimization algorithms for the ACOPF problem where the optimization is performed directly on the AC power flow manifold. This is done by using the Julia programming language and the Julia packages PowerModels.jl and Manopt.jl.

AC optimal power flow↗
DOE OSTI2021

On the eigenvector bias of Fourier feature networks: From regression to solving multi-scale PDEs with physics-informed neural networks

Physics-informed neural networks (PINNs) are demonstrating remarkable promise in integrating physical models with gappy and noisy observational data, but they still struggle in cases where the target functions to be approximated exhibit high-frequency or multi-scale features. Here in this work we investigate this limitation through the lens of Neural Tangent Kernel (NTK) theory and elucidate how PINNs are biased towards learning functions along the dominant eigen-directions of their limiting NTK. Using this observation, we construct novel architectures that employ spatio-temporal and multi-scale random Fourier features, and justify how such coordinate embedding layers can lead to robust and accurate PINN models. Numerical examples are presented for several challenging cases where conventional PINN models fail, including wave propagation and reaction–diffusion dynamics, illustrating how the proposed methods can be used to effectively tackle both forward and inverse problems involving partial differential equations with multi-scale behavior.

42 ENGINEERING↗
DOE OSTI2022

Dynamics retrieval from stochastically weighted incomplete data by low-pass spectral analysis

Time-resolved serial femtosecond crystallography (TR-SFX) provides access to protein dynamics on sub-picosecond timescales, and with atomic resolution. Due to the nature of the experiment, these datasets are often highly incomplete and the measured diffracted intensities are affected by partiality. To tackle these issues, one established procedure is that of splitting the data into time bins, and averaging the multiple measurements of equivalent reflections within each bin. This binning and averaging often involve a loss of information. Here, we propose an alternative approach, which we call low-pass spectral analysis (LPSA). In this method, the data are projected onto the subspace defined by a set of trigonometric functions, with frequencies up to a certain cutoff. This approach attenuates undesirable high-frequency features and facilitates retrieving the underlying dynamics. A time-lagged embedding step can be included prior to subspace projection to improve the stability of the results with respect to the parameters involved. Subsequent modal decomposition allows to produce a low-rank description of the system's evolution. Using a synthetic time-evolving model with incomplete and partial observations, we analyze the LPSA results in terms of quality of the retrieved signal, as a function of the parameters involved. We compare the performance of LPSA to that of a range of other sophisticated data analysis techniques. We show that LPSA allows to achieve excellent dynamics reconstruction at modest computational cost. Finally, we demonstrate the superiority of dynamics retrieval by LPSA compared to time binning and merging, which is, to date, the most commonly used method to extract dynamical information from TR-SFX data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1986

Multistep modeling of protein structure: application to bungarotoxin

Modelling of bungarotoxin in atomic details is presented in this article. The model-building procedure utilizes the low-resolution crystal coordinates of the c-alpha atoms of bungarotoxin, sequence homology within the neurotoxin family, as well as high-resolution x-ray diffraction data of cobratoxin and erabutoxin. Our model-building procedure involves: (a) principles of comparative modelling, (b) embedding procedures of distance geometry, and (c) use of molecular mechanics for optimizing packing. The model is not only consistent with the c-alpha coordinates of crystal structure, but also agrees with solution conformational features of the triple-stranded beta sheet as observed by NOE measurements.

NASA Discipline Exobiology↗
NASA NTRS2018

Aerodynamic Design and Optimization of Fan Stage for Boundary Layer Ingestion Propulsion System

The present paper addresses the process of preliminary design of a low-pressure fan and outlet guide vane (OGV) of a boundary layer ingestion (BLI) propulsion system. The tail-cone thruster systems of NASA's STARC_ABL (Single-aisle Turboelectric Aircraft with an Aft Boundary-Layer propulsor) adopts an axi-symmetric BLI type inlet as opposed to other embedded engine systems. Thus, the focus of the present work is placed on maximizing the efficiency of the fan and OGV stages under a significant radial distortion. A parameterization with B-spline function for camber line angles, metal chord, thickness distribution and stacking axis of blades is presented. The flowpath lines are also parameterized by B-spline function and aggregated in the design system of blades. The design optimization with evolutionary algorithm is performed with constraints of fan pressure ratio, OGV exit swirl angle and nozzle exit properties. The inlet conditions for the turbo-machinery CFD (Computational Fluid Dynamics) domain and the design goal of the fan stage are driven by a propulsion airframe integration (PAI) model that uses a 3-D unstructured RANS (Reynolds Average Navier Stokes) solver and actuator disk model. The expected power saving of the BLI propulsor is quantified via PAI analysis and the resulting preliminary design of the fan stages is compared with a clean-inlet flow propulsor.

Turbo-machinery↗
NASA NTRSDate not supplied

Evaluating Lunar Water Processing System Model Configurations for Small Scale Oxygen and Hydrogen Production Within JAXA’s ISRU Technology Demonstration Plant

In-Situ Resource Utilization (ISRU) refers to novel methods of extracting and processing local resources for use in life support and propulsion systems, reducing or eliminating the required consumables to be transferred from Earth. Current estimates of water-ice availability embedded in regolith within the Moon’s permanently shadowed regions (PSR’s)range between 1-5% by weight. However, the composition and characteristics of the “wet” regolith is unknown. Alternate ISRU excavation techniques and Concept of Operations (ConOps) must be explored to optimize surface system operations based on these factors. To assess the feasibility of different ISRU subsystem technologies and compare system architecture configurations, an interchangeable system model was generated to incorporate technologies spanning excavation of raw materials to storage of products and determine optimal arrangement of total system processing needs. Total Mass, Volume, and Power (M/V/P) requirements were computed for 168 design iterations of this water processing plant.

ISRU↗
NASA NTRSDate not supplied

Direct Reconstruction of Ablative Thermal Protection System Aeroheating Using A Green's Function Approach

A Green’s function inverse heat transfer (IHT) approach is used to reconstruct the surface heating conditions on ablative thermal protection system (TPS) materials from embedded heat flux sensor and temperature probe measurements. The approach models the temperature time-history at the measurement location as a discrete linear system, allowing for the heat flux boundary condition to be recovered directly without the need for time-marching schemes. The effects of material decomposition and pyrolysis gas transport are modeled using an energy source/sink analogue. The performance of the reconstruction approach is analyzed on a 1D test case representative of an atmospheric entry heating scenario. The approach can recover the TPS surface heat flux to within 4% of the input heating condition with a computation time of 2-3 seconds (>3 orders of magnitude faster than current time-marching IHT methods). As a byproduct of the surface heating reconstruction, the algorithm also captures the surface pyrolysis gas mass flux and solid decomposition at multiple through-thickness locations within the TPS.

Kenneth McAfee↗
NASA NTRSDate not supplied

Evaluating Lunar Water Processing System Model Configurations for Small Scale Oxygen and Hydrogen Production Within JAXA'S ISRU Technology

Introduction: In-Situ Resource Utilization (ISRU) refers to novel methods of extracting and processing local resources for use in life support and propulsion systems, reducing or eliminating the required consumables to be transferred from Earth. Current estimates of water-ice availability embedded in regolith within the Moon’s permanently shadowed regions (PSR’s) range between 1-5% by weight. However, the composition and characteristics of the “wet” regolith is unknown. Alternate ISRU excavation techniques and Concept of Operations (ConOps) must be explored to optimize surface system operations based on these factors. To assess the feasibility of different ISRU subsystem technologies and compare system architecture configurations, an interchangeable system model was generated to incorporate technologies spanning excavation of raw materials to storage of products and determine optimal arrangement of total system processing needs. Total Mass, Volume, and Power (M/V/P) requirements were computed for 168 design iterations of this water processing plant. System Model: In FY24, the System Engineering and Integration (SE&I) ISRU Modeling and Analysis (SIMA) team developed a lunar water processing system model using the Mission Analysis and Integration Tool (MAIT) to estimate the M/V/P for ISRU subsystems operating under a wide range of Hydrogen (H2) and Oxygen (O2) production targets for the Space Technology Mission Directorate (STMD) [1]. Based on Japan Aerospace Exploration Agency’s (JAXA) surface operational requirements, this system architecture was modified to include the ability to excavate consolidated icy regolith (versus granular ice excavation using Kennedy Space Center’s (KSC) ISRU Pilot Excavator, IPEx) and explore the feasibility of processing the lunar water both inside and outside of the PSR. For the consolidated icy regolith case study, excavation was performed via a mobility transport chassis outfitted with The Regolith Ice Drill for Exploring New Terrain (TRIDENT) for drilling [2] and the Cold Operable Lunar Deployable Arm (COLDArm) [3] for regolith transfer. The system model determines the required rover and payload. M/V/P to handle the required regolith processing rates. The regolith is then sorted and heated to sublimate the ice (via an auger dryer). The exiting high temperature, low pressure vapor is cleaned of volatiles (via cold trap) and electrolyzed to produce H2 and O2. These products are then dried, liquified with 20 K and 90 K cryocoolers (for H2 and O2, respectively), and stored in cylindrical tanks. Study Goals: Due to the different ConOps options of regolith transport to the ridge for processing versus processing it directly inside the PSR, as well as the unknown regolith/water-ice composition, new excavation techniques and their power configurations are being evaluated within a ISRU system architecture for production targets less than NASA’s pilot plant (1 mT). This analysis investigates the feasibility of numerous excavation techniques, power architectures, and logistical operations and determines an optimal system configuration with regards to M/V/P. It aims to investigate which parameters, both locally and globally, have the greatest effect on each subsystem within the plant. This can be used to identify the most critical components of the plant, and guide future decisions on allocating funding for research and development. The results from this study may provide subsystem developers with appropriate interfaces with excavation subsystems and downstream processes, and assessing the overall feasibility of each excavation technique, power architecture, and logistical timeframe. References: [1] Carlson, A. et. al. (2024) ICES. [2] Zacny, K., et. al. (2024) “ASCE Earth and Space”. [3] McCormick, R., et. Al. (2024) IEEE Xplore.

ISRU↗
DOE OSTI2023

Structure optimization with stochastic density functional theory

Linear-scaling techniques for Kohn–Sham density functional theory are essential to describe the ground state properties of extended systems. Still, these techniques often rely on the localization of the density matrix or accurate embedding approaches, limiting their applicability. In contrast, stochastic density functional theory (sDFT) achieves linear- and sub-linear scaling by statistically sampling the ground state density without relying on embedding or imposing localization. In return, ground state observables, such as the forces on the nuclei, fluctuate in sDFT, making optimizing the nuclear structure a highly non-trivial problem. In this work, we combine the most recent noise-reduction schemes for sDFT with stochastic optimization algorithms to perform structure optimization within sDFT. We compare the performance of the stochastic gradient descent approach and its variations (stochastic gradient descent with momentum) with stochastic optimization techniques that rely on the Hessian, such as the stochastic Broyden–Fletcher–Goldfarb–Shanno algorithm. In conclusion, we further provide a detailed assessment of the computational efficiency and its dependence on the optimization parameters of each method for determining the ground state structure of bulk silicon with varying supercell dimensions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Engineering Privacy at the Edge: A Practical Guide to Differential Privacy in System Architectures

The rapid expansion of distributed and edge computing platforms—spanning autonomous vehicles, IoT sensors, and healthcare monitors—has heightened concerns about data privacy. Differential Privacy (DP) offers a rigorous mathematical framework to protect sensitive information while retaining analytical utility. This tutorial introduces the foundations of DP for both numerical and categorical datasets and extends the discussion to correlation-aware techniques tailored for structured and high-dimensional data. Hands-on demonstrations will begin with the PETINA (Privacy prEservaTIoN Algorithms) package for numerical data and continue with MIC-DP (Maximum Information Correlated Differential Privacy) for tabular data. Designed for researchers and practitioners in secure systems, embedded architectures, and AI accelerators, the tutorial emphasizes practical and scalable methods for integrating DP into real-world system designs.

Kotevska, Olivera [ORNL] (ORCID:0000000316772243)↗
DOE OSTI2024

Exponential concentration in quantum kernel methods

Kernel methods in Quantum Machine Learning (QML) have recently gained significant attention as a potential candidate for achieving a quantum advantage in data analysis. Among other attractive properties, when training a kernel-based model one is guaranteed to find the optimal model’s parameters due to the convexity of the training landscape. However, this is based on the assumption that the quantum kernel can be efficiently obtained from quantum hardware. In this work we study the performance of quantum kernel models from the perspective of the resources needed to accurately estimate kernel values. We show that, under certain conditions, values of quantum kernels over different input data can be exponentially concentrated (in the number of qubits) towards some fixed value. Thus on training with a polynomial number of measurements, one ends up with a trivial model where the predictions on unseen inputs are independent of the input data. We identify four sources that can lead to concentration including expressivity of data embedding, global measurements, entanglement and noise. For each source, an associated concentration bound of quantum kernels is analytically derived. Lastly, we show that when dealing with classical data, training a parametrized data embedding with a kernel alignment method is also susceptible to exponential concentration. Our results are verified through numerical simulations for several QML tasks. Altogether, we provide guidelines indicating that certain features should be avoided to ensure the efficient evaluation of quantum kernels and so the performance of quantum kernel methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Localizing Clinical Patterns of Blast Traumatic Brain Injury Through Computational Modeling and Simulation

Blast traumatic brain injury is ubiquitous in modern military conflict with significant morbidity and mortality. Yet the mechanism by which blast overpressure waves cause specific intracranial injury in humans remains unclear. Reviewing of both the clinical experience of neurointensivists and neurosurgeons who treated service members exposed to blast have revealed a pattern of injury to cerebral blood vessels, manifested as subarachnoid hemorrhage, pseudoaneurysm, and early diffuse cerebral edema. Additionally, a seminal neuropathologic case series of victims of blast traumatic brain injury (TBI) showed unique astroglial scarring patterns at the following tissue interfaces: subpial glial plate, perivascular, periventricular, and cerebral gray-white interface. The uniting feature of both the clinical and neuropathologic findings in blast TBI is the co-location of injury to material interfaces, be it solid-fluid or solid-solid interface. This motivates the hypothesis that blast TBI is an injury at the intracranial mechanical interfaces. In order to investigate the intracranial interface dynamics, we performed a novel set of computational simulations using a model human head simplified but containing models of gyri, sulci, cerebrospinal fluid (CSF), ventricles, and vasculature with high spatial resolution of the mechanical interfaces. Simulations were performed within a hybrid Eulerian—Lagrangian simulation suite (CTH coupled via Zapotec to Sierra Mechanics). Because of the large computational meshes, simulations required high performance computing resources. Twenty simulations were performed across multiple exposure scenarios—overpressures of 150, 250, and 500 kPa with 1 ms overpressure durations—for multiple blast exposures (front blast, side blast, and wall blast) across large variations in material model parameters (brain shear properties, skull elastic moduli). All simulations predict fluid cavitation within CSF (where intracerebral vasculature reside) with cavitation occurring deep and diffusely into cerebral sulci. These cavitation events are adjacent to high interface strain rates at the subpial glial plate. Larger overpressure simulations (250 and 500kPa) demonstrated intraventricular cavitation—also associated with adjacent high periventricular strain rates. Additionally, models of embedded intraparenchymal vascular structures—with diameters as small as 0.6 mm—predicted intravascular cavitation with adjacent high perivascular strain rates. The co-location of local maxima of strain rates near several of the regions that appear to be preferentially damaged in blast TBI (vascular structures, subpial glial plate, perivascular regions, and periventricular regions) suggest that intracranial interface dynamics may be important in understanding how blast overpressures leads to intracranial injury.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Denoising diffusion algorithm for inverse design of microstructures with fine-tuned nonlinear material properties

Here we introduce a denoising diffusion algorithm to discover microstructures with nonlinear fine-tuned properties. Denoising diffusion probabilistic models are generative models that use diffusion-based dynamics to gradually denoise images and generate realistic synthetic samples. By learning the reverse of a Markov diffusion process, we design an artificial intelligence to efficiently manipulate the topology of microstructures to generate a massive number of prototypes that exhibit constitutive responses sufficiently close to designated nonlinear constitutive behaviors. To identify the subset of microcstructures with sufficiently precise fine-tuned properties, a convolutional neural network surrogate is trained to replace high-fidelity finite element simulations to filter out prototypes outside the admissible range. Results of this study indicate that the denoising diffusion process is capable of creating microstructures of fine-tuned nonlinear material properties within the latent space of the training data. More importantly, this denoising diffusion algorithm can be easily extended to incorporate additional topological and geometric modifications by introducing high-dimensional structures embedded in the latent space. Numerical experiments are conducted on the open-source mechanical MNIST data set (Lejeune, 2020). Consequently, this algorithm is not only capable of performing inverse design of nonlinear effective media, but also learns the nonlinear structure–property map to quantitatively understand the multiscale interplay among the geometry, topology, and their effective macroscopic properties.

42 ENGINEERING↗
DOE OSTI2023

Multifidelity Neural Network Formulations for Prediction of Reactive Molecular Potential Energy Surfaces

Here, this paper focuses on the development of multifidelity modeling approaches using neural network surrogates, where training data arising from multiple model forms and resolutions are integrated to predict high-fidelity response quantities of interest at lower cost. We focus on the context of quantum chemistry and the integration of information from multiple levels of theory. Important foundations include the use of symmetry function-based atomic energy vector constructions as feature vectors for representing structures across families of molecules and single-fidelity neural network training capabilities that learn the relationships needed to map feature vectors to potential energy predictions. These foundations are embedded within several multifidelity topologies that decompose the high-fidelity mapping into model-based components, including sequential formulations that admit a general nonlinear mapping across fidelities and discrepancy-based formulations that presume an additive decomposition. Methodologies are first explored and demonstrated on a pair of simple analytical test problems and then deployed for potential energy prediction for C 5 H 5 using B2PLYP-D3/6-311++G(d,p) for high-fidelity simulation data and Hartree–Fock 6-31G for low-fidelity data. For the common case of limited access to high-fidelity data, our computational results demonstrate that multifidelity neural network potential energy surface constructions achieve roughly an order of magnitude improvement, either in terms of test error reduction for equivalent total simulation cost or reduction in total cost for equivalent error.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1992

Development of an improved model for runback water on aircraft surfaces

A computer simulation for 'running wet' and evaporative aircraft anti-icing systems is developed. The model is based on the analysis of the liquid water film which forms in the regions of direct impingement and, then, breaks up near the impingement limits into rivulets. The wetness factor distribution resulting from the film breakup and the rivulet configuration on the surface are predicted using a stability analysis theory and the laws of mass energy conservation. The solid structure is modeled as a multiple layer wall. The anti-icing system modeled is of the thermal type utilizing hot air and/or electrical heating elements embedded within the wall layers. Experimental observations revealing some of the basic physics of the water flow on the surface are presented. Detailed qualitative documentation of the tests are given. Several numerical examples are considered, and the effect of some of the involved parameters on the system performance are investigated.

Al-Khalil, Kamel M.↗