Search NASA⌕ Search

SEARCH · Search NASA

Results for “computational frameworks”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

PhILMs: Collaboratory on Mathematics and Physics-Informed Learning Machines for Multiscale and Multiphysics Problems

The landscape of computational science and engineering is continually evolving, with the challenge of high-dimensional regression problems standing as a significant hurdle in numerous scientific endeavors. Addressing this challenge, our research, funded by this award, has led to the development of an innovative computational framework known as Probabilistic Partition of Unity Networks (PPOU-Nets). This initiative represents a collaborative effort to harness the potential of mathematics and physics-informed machine learning in tackling multiscale and multiphysics problems prevalent in high-dimensional spaces. Through this work, we have proposed a novel methodology that seamlessly integrates adaptive dimensionality reduction and a mixture of experts model, thereby facilitating a more efficient and accurate approximation of complex functions. This research effort has not only advanced the state of computational science but also opened new avenues for exploration in quantum computing and beyond. This report outlines the motivation, methodology, key findings, and implications of our work, underscoring our contributions to the broader scientific community and the potential pathways for future research.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Collision of localized shocks in AdS 5 as a series expansion in transverse gradients

We introduce a computational framework to more efficiently calculate the collision of localized shocks in five dimensional asymptotically Anti-de Sitter space. We expand the Einstein equations in transverse gradients and find that our numerical results agree well with exact solutions already at first order in the expansion. Moreover, the Einstein equations at first order in transverse gradients can be decoupled into two sets of differential equations. Here, the bulk fields of one of these sets has only a negligible contribution to boundary observables, such that the computation on each time slice can be simplified to the solution of several planar shockwave equations plus four further differential equations for each transverse plane ‘pixel’. At the cost of errors of ≲ 10% at the hydrodynamization time and for low to mid rapidities, useful numerical solutions can be sped up by roughly one order of magnitude.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Numerical Investigation of the Spark Discharge Process in a Crossflow

The present study numerically investigates the spark discharge process under crossflow conditions using a thermal equilibrium plasma solver that fully couples the electromagnetic physics and fluid dynamics in a computational framework. Numerical results are validated by the comparison with experimental data. Here, the spark discharge experiment is performed in a constant volume vessel using an inductive coil ignition system for automotive applications, and the evolution of the spark channel is measured using high-speed imaging. The crossflow in the gap between the spark-plug electrodes is generated by a rotating fan with two different fan speeds, and the flow velocity across the gap is characterized by particle image velocimetry (PIV) measurement. A computational fluid dynamics (CFD) solver is employed to simulate the crossflow and provide the flow field variables (velocity, pressure, temperature) to the plasma solver. The crossflow velocity predicted in the flow simulation agrees well with the PIV data in that the non-uniform velocity profiles at monitoring points are reproduced by the CFD code. With the crossflow initialization in the plasma solver, the simulated spark discharge process from the breakdown to spark discharge matches the experimental data, including the voltage and circuit waveforms and the high-speed images of the spark channel evolution. The stretch of spark channel captured by plasma simulations agrees with the measured data. The plasma simulation reveals that the mean temperature of the spark channel is maintained at 5000 K during the discharge phase, and the temperature varies along the spark channel so that the highest value is obtained at the spark root on the center electrode. Overall, the results presented in this paper are meant to provide valuable information about the properties of the plasma generated by the spark discharge.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

GnarlyX: Eulerian multi-material hydrodynamics coupled to equation of state and hyperelastic, plastic constitutive models

GnarlyX is a new hydrocode for direct numerical simulations of the microstructural behavior of high explosives at the mesoscale. We summarize the computational framework for multidimensional, Eulerian multi-material hydrodynamics coupled to EOS and hyperelastic, plastic constitutive models. We present 1D verification tests of multi-material only and combined multimaterial and strength capabilities with comparisons to exact solutions of shock states resulting from an incident shock impacting the material interface of PBX9502 and aluminum. We show that GnarlyX performs well in capturing the resulting shock waves in examining numerical convergence with exact solutions. In later work, we will summarize the thermomechanics and multi-dimensional, parallel computing capabilities in GnarlyX with multi-dimensional verification tests.

42 ENGINEERING↗
DOE OSTI2021

A new coupled multiphase flow–finite strain deformation–fault slip framework for induced seismicity

Production of hydrocarbons and water from subsurface reservoirs are known to cause permanent deformation of the reservoir and seismicity along faults both of which are detrimental to sustainable development of natural resources. Most of the prior studies on understanding fluid flow-induced plasticity and seismicity have focused on one or the other phenomenon due to the numerical difficulty associated with simultaneous modeling of the two failure phenomena because stress and deformation evolve non-linearly in both plasticity and seismicity. However, in reservoirs undergoing long-term production, plastic failure can alter the stress paths of points on a fault such that the onset, location, and magnitude of actual seismic events can no longer be predicted by a poroelastic simulation due to inaccurate stress and deformation history. We present a computational framework for coupled multiphase flow, finite strain poroplastic deformation, and dynamic fault slip and use it to understand the impact of plastic deformation on the onset, location, and magnitude of induced fault slip events. We evaluate the impact of plasticity on reservoir pressure, deformation, induced stress, and fault slip by comparing infinitesimal strain elastic and finite strain poro-elastoplastic models. For real-world applications, we consider different scenarios where the reservoir is either mechanically weaker or stronger than the caprock. We analyze the stress evolution as a function of the change in reservoir pressure to understand the role of contrast in reservoir and caprock elastic moduli on geomechanical stability. The results show that the poroplastic reservoir exhibits larger vertical deformation and delayed slip than the poroelastic reservoir after the same amount of oil production. For the same amount of pressure drop, a reservoir with a smaller modulus than the caprock displays a larger vertical displacement and an earlier onset of both plastic failure and fault slip. For a reservoir with a larger modulus than the caprock, vertical displacement is larger on the reservoir top boundary and smaller on the ground surface, and a higher pressure drop is needed to induce plastic failure and fault slip.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Adaptive Conformer Sampling for Property Prediction Using the Conductor-like Screening Model for Real Solvents

The valorization of lignocellulose-derived bioproducts requires effective separation from excessive water. Liquid–liquid extraction is a promising low-energy separation technology, but effective extraction requires solvent selection based on the thermodynamic properties of the bioproduct and solvent components. We propose a computational framework for predicting such properties by developing an adaptive conformer selection approach for use with COSMO-RS (conductor-like screening model for real solvents) calculations. In this framework, molecular dynamics simulations are used to generate many molecular structures (conformers) at representative temperatures in varying solvent environments. Conformers are then clustered based on structural metrics in a low-dimensional space and selected using a mixed-integer quadratic programming problem to iteratively insert a sampled conformer. At each iteration, we determine bioproduct properties using COSMO-RS. Here, we demonstrate the capability of the proposed framework on representative bioproducts to show convergence of the adaptive sampling toward experimentally measured properties with fewer calculations than required by random conformer sampling, enabling the improved screening of solvent systems for liquid-phase separation.

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

Multimodal neuroelectric interface development

We are developing electromyographic and electroencephalographic methods, which draw control signals for human-computer interfaces from the human nervous system. We have made progress in four areas: 1) real-time pattern recognition algorithms for decoding sequences of forearm muscle activity associated with control gestures; 2) signal-processing strategies for computer interfaces using electroencephalogram (EEG) signals; 3) a flexible computation framework for neuroelectric interface research; and d) noncontact sensors, which measure electromyogram or EEG signals without resistive contact to the body.

Evaluation Studies↗
DOE OSTI2024

APACE: AlphaFold2 and advanced computing as a service for accelerated discovery in biophysics

The prediction of protein 3D structure from amino acid sequence is a computational grand challenge in biophysics and plays a key role in robust protein structure prediction algorithms, from drug discovery to genome interpretation. The advent of AI models, such as AlphaFold, is revolutionizing applications that depend on robust protein structure prediction algorithms. To maximize the impact, and ease the usability, of these AI tools we introduce APACE, AlphaFold2 and advanced computing as a service, a computational framework that effectively handles this AI model and its TB-size database to conduct accelerated protein structure prediction analyses in modern supercomputing environments. We deployed APACE in the Delta and Polaris supercomputers and quantified its performance for accurate protein structure predictions using four exemplar proteins: 6AWO, 6OAN, 7MEZ, and 6D6U. Using up to 300 ensembles, distributed across 200 NVIDIA A100 GPUs, we found that APACE is up to two orders of magnitude faster than off-the-self AlphaFold2 implementations, reducing time-to-solution from weeks to minutes. This computational approach may be readily linked with robotics laboratories to automate and accelerate scientific discovery.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Multi-scale computational analysis of unidirectional carbon fiber reinforced polymer composites under various loading conditions

In this study, a multi-scale computational analysis based on representative volume element (RVE) modeling and molecular dynamics (MD) simulations is developed to investigate the microscopic failure mechanisms of unidirectional (UD) carbon fiber reinforced polymer (CFRP) composites. The average properties of the 200-nm thickness interphase region between fiber and matrix are characterized through MD simulations and an analytical gradient model. The results demonstrate that the interphase region has higher Young’s modulus and strength, compared to the bulk matrix. This stiffened interphase region influences the composite response significantly. Specifically, the traditional two-phase model with zero-thickness interface fails to capture the stress-strain behavior compared to the experimental data. However, by adding the interphase region to a modified RVE model, the accuracy of simulation results will be improved significantly. Furthermore, a coupled experimental-computational micromechanics approach is adopted to calibrate and validate the cohesive parameters of the interface. By including the cohesive interface, our modified RVE model accurately captures the failure strength of the composites. Finally, different failure mechanisms for specimens are investigated using our multi-scale computational framework. In conclusion, the results show that the failure modes of UD CFRP composites are very complex and multiple failure mechanisms co-exist depending on the loading conditions, agreeing well with our experimental analyses.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Hierarchical ensemble Kalman methods with sparsity-promoting generalized gamma hyperpriors

This paper introduces a computational framework to incorporate flexible regularization techniques in ensemble Kalman methods, generalizing the iterative alternating scheme to nonlinear inverse problems. The proposed methodology approximates the maximum a posteriori (MAP) estimate of a hierarchical Bayesian model characterized by a conditionally Gaussian prior and generalized gamma hyperpriors. Suitable choices of hyperparameters yield sparsity-promoting regularization. We propose an iterative algorithm for MAP estimation, which alternates between updating the unknown with an ensemble Kalman method and updating the hyperparameters in the regularization to promote sparsity. Here, the effectiveness of our methodology is demonstrated in several computed examples, including compressed sensing and subsurface flow inverse problems.

Ensemble Kalman methods↗
DOE OSTI2023

Reduced-order autodifferentiable ensemble Kalman filters

This paper introduces a computational framework to reconstruct and forecast a partially observed state that evolves according to an unknown or expensive-to-simulate dynamical system. Our reduced-order autodifferentiable ensemble Kalman filters (ROAD-EnKFs) learn a latent low-dimensional surrogate model for the dynamics and a decoder that maps from the latent space to the state space. The learned dynamics and decoder are then used within an EnKF to reconstruct and forecast the state. Numerical experiments show that if the state dynamics exhibit a hidden low-dimensional structure, ROAD-EnKFs achieve higher accuracy at lower computational cost compared to existing methods. If such structure is not expressed in the latent state dynamics, ROAD-EnKFs achieve similar accuracy at lower cost, making them a promising approach for surrogate state reconstruction and forecasting.

Mathematics↗
DOE OSTI2018

Distributed caching for processing raw arrays

As applications continue to generate multi-dimensional data at exponentially increasing rates, fast analytics to extract meaningful results is becoming extremely important. The database community has developed array databases that alleviate this problem through a series of techniques. In-situ mechanisms provide direct access to raw data in the original format---without loading and partitioning. Parallel processing scales to the largest datasets. In-memory caching reduces latency when the same data are accessed across a workload of queries. However, we are not aware of any work on distributed caching of multi-dimensional raw arrays. In this paper, we introduce a distributed framework for cost-based caching of multi-dimensional arrays in native format. Given a set of files that contain portions of an array and an online query workload, the framework computes an effective caching plan in two stages. First, the plan identifies the cells to be cached locally from each of the input files by continuously refining an evolving R-tree index. In the second stage, an optimal assignment of cells to nodes that collocates dependent cells in order to minimize the overall data transfer is determined. We design cache eviction and placement heuristic algorithms that consider the historical query workload. A thorough experimental evaluation over two real datasets in three file formats confirms the superiority -- by as much as two orders of magnitude -- of the proposed framework over existing techniques in terms of cache overhead and workload execution time.

Zhao, Weijie↗
DOE OSTI2024

Inverse-Velocity Transformation Wall Model for Reacting Turbulent Hypersonic Boundary Layers

The present study builds on prior work by taking advantage of the novel framework proposed by Griffin et al. (hereafter referred to as the GFM) as a baseline. The model is progressively extended to multicomponent reacting mixtures, accounting for differential diffusion and finite-rate chemistry in a similar fashion to Di Renzo & Urzay and Di Renzo et al. The accuracy of the present approach, as well the prior model of Di Renzo & Urzay, is assessed in an a priori sense for the first time in a turbulent reacting boundary layer, using the boundary-layer data of Williams et al. Five species are included in the present analysis, i.e., N s = 5, namely N 2 , O 2 , NO, N and O, a neutral mixture most representative of dissociation/recombination phenomena for temperatures below 6000 K. A schematic of a flow over a wedge representative of the described configuration is presented in Figure 1. The brief is organized as follows: In Section 2, the wall-model equations and the computational framework are presented. In Section 3, the a priori results of the proposed model are described and compared to the extended EWM. Finally, in Section 4, some conclusions are offered.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Methane utilization in Methylomicrobium alcaliphilum 20Z R : a systems approach

Biological methane utilization, one of the main sinks of the greenhouse gas in nature, represents an attractive platform for production of fuels and value-added chemicals. Despite the progress made in our understanding of the individual parts of methane utilization, our knowledge of how the whole-cell metabolic network is organized and coordinated is limited. Attractive growth and methane-conversion rates, a complete and expert-annotated genome sequence, as well as large enzymatic, 13 C-labeling, and transcriptomic datasets make Methylomicrobium alcaliphilum 20Z R an exceptional model system for investigating methane utilization networks. Here we present a comprehensive metabolic framework of methane and methanol utilization in M. alcaliphilum 20Z R . A set of novel metabolic reactions governing carbon distribution across central pathways in methanotrophic bacteria was predicted by in-silico simulations and confirmed by global non-targeted metabolomics and enzymatic evidences. Our data highlight the importance of substitution of ATP-linked steps with PPi-dependent reactions and support the presence of a carbon shunt from acetyl-CoA to the pentose-phosphate pathway and highly branched TCA cycle. The diverged TCA reactions promote balance between anabolic reactions and redox demands. As a result, the computational framework of C 1 -metabolism in methanotrophic bacteria can represent an efficient tool for metabolic engineering or ecosystem modeling.

09 BIOMASS FUELS↗
DOE OSTI2016

An Interface Debonding Model for Integrity Assessment of Nuclear Fuel Rods

The integrity degradation of nuclear fuel rods due to material embrittlement resulting from hydride formation has been extensively investigated. Finite element based fracture mechanics models have been proposed to assess the reduction of load capacity of nuclear fuel rod structures with several distribution scenarios of hydrides. In these finite element models, individual hydrides are modeled by opened cracks distributed in bulk matrix. These models are indeed valid for cases where hydrides fractures occur at low load levels. It is important for finite element models to take into account the cases where hydride fractures occur at high load levels, at which bulk matrix material may undergo reasonable high stress or yield. This requires hydride/bulk matrix interface models that can accommodate both bonded and debonded (fractured or crack initiation) behaviors of hydride/matrix interfaces. However, direct modeling debonding behaviors between hydride and bulk matrix is a great challenge. This is because for ductile structural failure simulations ligaments of bulk matrix often undergo large plastic deformation. Complete solutions often unavailable for crack initiation or debonding in media with finite plastic strain using popular implicit finite element technologies. This paper proposes a new computational framework based on discontinuous Galerkin (DG) finite element formulations for easily handling hydride/matrix interface debonding. We present structural failure simulation results for a quarter plane strain rod model with zircaloy materials. It is concluded that the proposed DG interface debonding model is powerful in simulating structural failure due to hydride failure in bulk media with finite plastic deformation. (author)

36 MATERIALS SCIENCE↗
DOE OSTI2019

Shape optimization under uncertainty for rotor blades of horizontal axis wind turbines

We present a computational framework for the shape optimization of a Horizontal-Axis Wind Turbine (HAWT) rotor blade under uncertainty. Our framework integrates aerodynamic simulations based on the blade element method which utilizes reduced order models of the blade structure and wind load with design sensitivity analysis and nonlinear programming. The wind velocity is modeled as a stochastic process to account for variations in time and space. An additional stochastic process accounts for uncertainties in the structural material properties. The uncertainty propagation is based on a non-intrusive polynomial chaos expansion that allows accurate estimation of stochastic performance metrics such as generated power and structural compliance. Sensitivities of cost and constraint functions with respect to shape parameters namely twist angles are computed with an efficient scheme to enable gradient based optimization. To demonstrate the effect of uncertainty, designs obtained from optimization under uncertainty are compared to those obtained from deterministic optimization.

42 ENGINEERING↗
DOE OSTI2016

Comprehensive numerical methodology for direct numerical simulations of compressible Rayleigh-Taylor instability

A comprehensive numerical methodology has been developed that handles the challenges introduced by considering the compressive nature of Rayleigh-Taylor instability (RTI) systems, which include sharp interfacial density gradients on strongly stratified background states, acoustic wave generation and removal at computational boundaries, and stratification-dependent vorticity production. The computational framework is used to simulate two-dimensional single-mode RTI to extreme late-times for a wide range of flow compressibility and variable density effects. The results show that flow compressibility acts to reduce the growth of RTI for low Atwood numbers, as predicted from linear stability analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

MODIFICATION OF THE MOOG SPECTRAL SYNTHESIS CODES TO ACCOUNT FOR ZEEMAN BROADENING OF SPECTRAL LINES

In an attempt to widen access to the study of magnetic fields in stellar astronomy, I present MOOGStokes, a version of the MOOG one-dimensional local thermodynamic equilibrium radiative transfer code, overhauled to incorporate a Stokes vector treatment of polarized radiation through a magnetic medium. MOOGStokes is a suite of three complementary programs, which together can synthesize the disk-averaged emergent spectrum of a star with a magnetic field. The first element (a pre-processing script called CounterPoint) calculates for a given magnetic field strength, wavelength shifts, and polarizations for the components of Zeeman-sensitive lines. The second element (a MOOG driver called SynStokes derived from the existing MOOG driver Synth) uses the list of Zeeman-shifted absorption lines together with the existing machinery of MOOG to synthesize the emergent spectrum at numerous locations across the stellar disk, accounting for stellar and magnetic field geometry. The third and final element (a post-processing script called DiskoBall) calculates the disk-averaged spectrum by weighting the individual emergent spectra by limb darkening and projected area, and applying the effects of Doppler broadening. All together, the MOOGStokes package allows users to synthesize emergent spectra of stars with magnetic fields in a familiar computational framework. MOOGStokes produces disk-averaged spectra for all Stokes vectors ( I, Q, U, V ), normalized by the continuum. MOOGStokes agrees well with the predictions of INVERS10 a polarized radiative transfer code with a long history of use in the study of stellar magnetic fields. In the non-magnetic limit, MOOGStokes also agrees with the predictions of the scalar version of MOOG.

79 ASTRONOMY AND ASTROPHYSICS↗