Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Implicit”

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,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Acceleration of Power System Dynamic Simulations Using a Deep Equilibrium Layer and Neural ODE Surrogate

The dominant paradigm for power system dynamic simulation is to build system-level simulations by combining physics-based models of individual components. The sheer size of the system along with the rapid integration of inverter-based resources exacerbates the computational burden of running time domain simulations. Here, in this paper, we propose a data-driven surrogate model based on implicit machine learningspecifically deep equilibrium layers and neural ordinary differential equationsto learn a reduced order model of a portion of the full underlying system. The data-driven surrogate achieves similar accuracy and reduction in simulation time compared to a physics-based surrogate, without the constraint of requiring detailed knowledge of the underlying dynamic models. This work also establishes key requirements needed to integrate the surrogate into existing simulation workflows; the proposed surrogate is initialized to a steady state operating point that matches the power flow solution by design.

Neural ordinary differential equations↗
DOE OSTI2022

CoordNet: Data Generation and Visualization Generation for Time-Varying Volumes via a Coordinate-Based Neural Network

Although deep learning has demonstrated its capability in solving diverse scientific visualization problems, it still lacks generalization power across different tasks. To address this challenge, we propose CoordNet, a single coordinate-based framework that tackles various tasks relevant to time-varying volumetric data visualization without modifying the network architecture. The core idea of our approach is to decompose diverse task inputs and outputs into a unified representation (i.e., coordinates and values) and learn a function from coordinates to their corresponding values. We achieve this goal using a residual block-based implicit neural representation architecture with periodic activation functions. We evaluate CoordNet on data generation (i.e., temporal super-resolution and spatial super-resolution) and visualization generation (i.e., view synthesis and ambient occlusion prediction) tasks using time-varying volumetric data sets of various characteristics. Here, the experimental results indicate that CoordNet achieves better quantitative and qualitative results than the state-of-the-art approaches across all the evaluated tasks. Source code and pre-trained models are available at https://github.com/stevenhan1991/CoordNet.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Estimating Higher-Order Moments Using Symmetric Tensor Decomposition

In this paper, we consider the problem of decomposing higher-order moment tensors, i.e., the sum of symmetric outer products of data vectors. Such a decomposition can be used to estimate the means in a Gaussian mixture model and for other applications in machine learning. The dth-order empirical moment tensor of a set of p observations of n variables is a symmetric d-way tensor. Our goal is to nd a low-rank tensor approximation comprising r $\ll$ p symmetric outer products. The challenge is that forming the empirical moment tensor costs O(pn d ) operations and O(n d ) storage, which may be prohibitively expensive; additionally, the algorithm to compute the low-rank approximation costs O(n d ) per iteration. Our contribution is avoiding formation of the moment tensor, computing the low-rank tensor approximation of the moment tensor implicitly using O(pnr) operations per iteration and no extra memory. This advance opens the door to more applications of higher-order moments since they can now be efficiently computed. We present numerical evidence of the computational savings and show an example of estimating the means for higher-order moments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

The Transport Equation in Optically Thick Media: Discussion of IMC and its Diffusion Limit

We discuss the limits of validity of the Implicit Monte Carlo (IMC) method for the transport of thermally emitted radiation. The weakened coupling between the radiation and material energy of the IMC method causes defects in handling problems with strong transients. We introduce an approach to asymptotic analysis for the transport equation that emphasizes the fact that the radiation and material temperatures are always different in time-dependent problems, and we use it to show that IMC does not produce the correct diffusion limit. As this is a defect of IMC in the continuous equations, no improvement to its discretization can remedy it.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Demonstration of a new unstructured mesh IMC x-ray transport capability in LAP codes

In this document, the Advanced Simulation and Computing (ASC) Transport Project’s Jayenne team presents evidence that the Los Alamos National Laboratory (LANL) Level 2 milestone statement, Demonstrate a new unstructured mesh IMC x-ray transport capability in LAP codes, due Q4 FY20, has been satisfied with the release and integration of the Jayenne team’s Implicit Monte Carlo (IMC) solver libraries into the lumos multiphysics solver. The Jayenne project’s software includes the algorithms and features prescribed by the milestone description. The lumos software is maintained by the Lagrangian Applications Project (LAP) along with the hydrodynamics code flag and operates under the same ASC program as Jayenne. A set of demonstration problems has been specified, executed, and analyzed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2022

Electron Heat Fluxes Generated by Intense Whistler Waves at the Upper Ionospheric Altitudes

Magnetospheric whistler waves, chorus and hiss, can't provide the resonance heating of the core electron plasma population. However, these whistler-mode branches, can implicitly participate in the heating processes of the core plasma thermal electron population by triggering the electron precipitation over a broad energy range from the magnetosphere and subsequent atmospheric ionization processes leading to the production of superthermal electron population. These superthermal electrons play a large role in the magnetosphere-ionosphere-atmosphere energy interplay with participation of both magnetically conjugate hemispheres, and their Coulomb interaction with background magnetospheric thermal electrons. Using strong hiss and chorus wave events measured by the Van Allen Probes and SuperThermal Electron Transport code, we evaluate the formation of electron heat fluxes at the upper ionospheric altitudes and discuss their consequences on the formation of electron temperature. It is found that chorus and hiss waves that initiate the precipitation of magnetospheric electrons with energies below 30 keV and the follow-up production of secondary electrons play an important role in the energy balance of ionosphere-magnetosphere system.

Magnetospheric whistler waves as an implicit sourc↗
DOE OSTI2012

Benchmark On Sensitivity Calculation (Phase III)

The sensitivities of the keff eigenvalue to neutron cross sections have become commonly used in similarity studies and as part of the validation algorithm for criticality safety assessments. To test calculations of the sensitivity coefficients, a benchmark study (Phase III) has been established by the OECD-NEA/WPNCS/EG UACSA (Expert Group on Uncertainty Analysis for Criticality Safety Assessment). This paper presents some sensitivity results generated by the benchmark participants using various computational tools based upon different computational methods: SCALE/TSUNAMI-3D and -1D, MONK, APOLLO2-MORET 5, DRAGON-SUSD3D and MMKKENO. The study demonstrates the performance of the tools. It also illustrates how model simplifications impact the sensitivity results and demonstrates the importance of 'implicit' (self-shielding) sensitivities. This work has been a useful step towards verification of the existing and developed sensitivity analysis methods.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

Extended magnetohydrodynamics with embedded particle-in-cell simulation of Ganymede's magnetosphere

Here, we have recently developed a new modeling capability to embed the implicit particle–in–cell (PIC) model iPIC3D into the Block–Adaptive–Tree–Solarwind–Roe–Upwind–Scheme magnetohydrodynamic (MHD) model. The MHD with embedded PIC domains (MHD–EPIC) algorithm is a two–way coupled kinetic–fluid model. As one of the very first applications of the MHD–EPIC algorithm, we simulate the interaction between Jupiter's magnetospheric plasma and Ganymede's magnetosphere. We compare the MHD–EPIC simulations with pure Hall MHD simulations and compare both model results with Galileo observations to assess the importance of kinetic effects in controlling the configuration and dynamics of Ganymede's magnetosphere. We find that the Hall MHD and MHD–EPIC solutions are qualitatively similar, but there are significant quantitative differences. In particular, the density and pressure inside the magnetosphere show different distributions. For our baseline grid resolution the PIC solution is more dynamic than the Hall MHD simulation and it compares significantly better with the Galileo magnetic measurements than the Hall MHD solution. The power spectra of the observed and simulated magnetic field fluctuations agree extremely well for the MHD–EPIC model. The MHD–EPIC simulation also produced a few flux transfer events (FTEs) that have magnetic signatures very similar to an observed event. The simulation shows that the FTEs often exhibit complex 3–D structures with their orientations changing substantially between the equatorial plane and the Galileo trajectory, which explains the magnetic signatures observed during the magnetopause crossings. The computational cost of the MHD–EPIC simulation was only about 4 times more than that of the Hall MHD simulation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Complex functionality with minimal computation: Promise and pitfalls of reduced‐tracer ocean biogeochemistry models

Abstract Earth System Models increasingly include ocean biogeochemistry models in order to predict changes in ocean carbon storage, hypoxia, and biological productivity under climate change. However, state‐of‐the‐art ocean biogeochemical models include many advected tracers, that significantly increase the computational resources required, forcing a trade‐off with spatial resolution. Here, we compare a state‐of‐the art model with 30 prognostic tracers (TOPAZ) with two reduced‐tracer models, one with 6 tracers (BLING), and the other with 3 tracers (miniBLING). The reduced‐tracer models employ parameterized, implicit biological functions, which nonetheless capture many of the most important processes resolved by TOPAZ. All three are embedded in the same coupled climate model. Despite the large difference in tracer number, the absence of tracers for living organic matter is shown to have a minimal impact on the transport of nutrient elements, and the three models produce similar mean annual preindustrial distributions of macronutrients, oxygen, and carbon. Significant differences do exist among the models, in particular the seasonal cycle of biomass and export production, but it does not appear that these are necessary consequences of the reduced tracer number. With increasing CO 2 , changes in dissolved oxygen and anthropogenic carbon uptake are very similar across the different models. Thus, while the reduced‐tracer models do not explicitly resolve the diversity and internal dynamics of marine ecosystems, we demonstrate that such models are applicable to a broad suite of major biogeochemical concerns, including anthropogenic change. These results are very promising for the further development and application of reduced‐tracer biogeochemical models that incorporate “sub‐ecosystem‐scale” parameterizations.

58 GEOSCIENCES↗
DOE OSTI2016

Comparison of equilibrium ohmic and nonequilibrium swarm models for monitoring conduction electron evolution in high-altitude EMP calculations

Here, atmospheric electromagnetic pulse (EMP) events are important physical phenomena that occur through both man-made and natural processes. Radiation-induced currents and voltages in EMP can couple with electrical systems, such as those found in satellites, and cause significant damage. Due to the disruptive nature of EMP, it is important to accurately predict EMP evolution and propagation with computational models. CHAP-LA (Compton High Altitude Pulse-Los Alamos) is a state-of-the-art EMP code that solves Maxwell inline images equations for gamma source-induced electromagnetic fields in the atmosphere. In EMP, low-energy, conduction electrons constitute a conduction current that limits the EMP by opposing the Compton current. CHAP-LA calculates the conduction current using an equilibrium ohmic model. The equilibrium model works well at low altitudes, where the electron energy equilibration time is short compared to the rise time or duration of the EMP. At high altitudes, the equilibration time increases beyond the EMP rise time and the predicted equilibrium ionization rate becomes very large. The ohmic model predicts an unphysically large production of conduction electrons which prematurely and abruptly shorts the EMP in the simulation code. An electron swarm model, which implicitly accounts for the time evolution of the conduction electron energy distribution, can be used to overcome the limitations exhibited by the equilibrium ohmic model. We have developed and validated an electron swarm model previously in Pusateri et al. (2015). Here we demonstrate EMP damping behavior caused by the ohmic model at high altitudes and show improvements on high-altitude, upward EMP modeling obtained by integrating a swarm model into CHAP-LA.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Mapping of Texture and Phase Fractions in Heterogeneous Stress States during Multiaxial Loading of Biomedical Superelastic NiTi

Thermoelastic deformation mechanisms in polycrystalline biomedical-grade superelastic NiTi are spatially mapped using in situ neutron diffraction during multiaxial loading and heating. The trigonal R-phase is formed from the cubic phase during cooling to room temperature and subsequently deforms in compression, tension, and torsion. The resulting R-phase variant microstructure from the variant reorientation and detwinning processes are equivalent for the corresponding strain in tension and compression, and the variant microstructure is reversible by isothermal loading. The R-phase variant microstructure is consistent between uniaxial and torsional loading when the principal stress directions of the stress state are considered (for the crystallographic directions observed here). The variant microstructure evolution is tracked and the similarity in general behavior between uniaxial and torsional loading, in spite of the implicit heterogeneous stress state associated with torsional loading, pointed to the ability of the reversible thermoelastic transformation in NiTi to accommodate stress and strain mismatch with deformation. This ability of the R-phase, despite its limited variants, to accommodate stress and strain and satisfy strain incompatibility in addition to the existing internal stresses has significance for reducing irrecoverable deformation mechanisms during loading and cycling through the phase transformation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Comment on “Introduction of a Novel Figure of Merit for the Assessment of Transparent Conductive Electrodes in Photovoltaics: Exact and Approximate Form”

Anand et al. introduce an implicit equation for transparent conductive electrode (TCE) merit in photovoltaics (PV) by lumping series resistance in with the diode. However, an explicit equation has been previously derived by lumping series resistance external to the diode. For TCEs of practical interest for PV, the choice of where to lump series resistance does not affect merit, so external lumping is preferred. On the other hand, monolith deadspace, metal grid shadowing, and grid resistance shift TCE merit away from the figure of merit behavior published by Anand et al.

14 SOLAR ENERGY↗
DOE OSTI2018

Prediction of Carbon Dioxide Adsorption via Deep Learning

Abstract Porous carbons with different textural properties exhibit great differences in CO 2 adsorption capacity. It is generally known that narrow micropores contribute to higher CO 2 adsorption capacity. However, it is still unclear what role each variable in the textural properties plays in CO 2 adsorption. Herein, a deep neural network is trained as a generative model to direct the relationship between CO 2 adsorption of porous carbons and corresponding textural properties. The trained neural network is further employed as an implicit model to estimate its ability to predict the CO 2 adsorption capacity of unknown porous carbons. Interestingly, the practical CO 2 adsorption amounts are in good agreement with predicted values using surface area, micropore and mesopore volumes as the input values simultaneously. This unprecedented deep learning neural network (DNN) approach, a type of machine learning algorithm, exhibits great potential to predict gas adsorption and guide the development of next‐generation carbons.

Zhang, Zihao↗
DOE OSTI2018

Prediction of Carbon Dioxide Adsorption via Deep Learning

Porous carbons with different textural properties exhibit great differences in CO 2 adsorption capacity. It is generally known that narrow micropores contribute to higher CO 2 adsorption capacity. However, it is still unclear what role each variable in the textural properties plays in CO 2 adsorption. Herein, a deep neural network is trained as a generative model to direct the relationship between CO 2 adsorption of porous carbons and corresponding textural properties. The trained neural network is further employed as an implicit model to estimate its ability to predict the CO 2 adsorption capacity of unknown porous carbons. Interestingly, the practical CO 2 adsorption amounts are in good agreement with predicted values using surface area, micropore and mesopore volumes as the input values simultaneously. This unprecedented deep learning neural network (DNN) approach, a type of machine learning algorithm, exhibits great potential to predict gas adsorption and guide the development of next-generation carbons.

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

A Route to Design Novel Functional Peptides by Applying a Denoising Diffusional Model to mRNA Display Libraries

In vitro directed evolution techniques, such as mRNA display, enable peptide ligand discovery and optimization. However, physical libraries that rely on a genetic code can only search a small fraction of sequence space due to inherent biases in the genetic code and experimental limitations. To address this challenge, denoising diffusion implicit models (DDIMs) are applied to generate novel peptide ligands against B‐cell lymphoma extra‐large (Bcl‐x L ), a key cancer target. Starting with high‐throughput sequencing data from previous selections, a DDIM is trained to produce novel sequences with high affinity binding. Experimental validation confirms that most generated sequences are functionally equivalent to the original library members for Bcl‐x L binding and demonstrated comparable binding kinetics and affinity relative to the wildtype and nearest original neighbors. Importantly, this approach generated rare sequences not easily accessible via mutation and directed evolution. These results indicate that DDIMs can complement and expand directed evolution data, efficiently exploring underrepresented regions of sequence space. This approach provides a broadly applicable framework for accelerating ligand discovery and optimizing molecular properties across diverse targets.

Qi, Pearl [Mork Family Department of Chemical Engi↗
DOE OSTI2019

Simulation of unsteady blood flows in a patient-specific compliant pulmonary artery with a highly parallel monolithically coupled fluid-structure interaction algorithm

Computational fluid dynamics (CFD) is increasingly used to study blood flows in patient-specific arteries for understanding certain cardiovascular diseases. The techniques work quite well for relatively simple problems, but need improvements when the problems become harder in the case when (1) the geometry becomes complex (from a few branches to a full pulmonary artery), (2) the model becomes more complex (from fluid-only calculation to coupled fluid-structure interaction calculation), (3) both the fluid and wall models become highly nonlinear, and (4) the computer on which we run the simulation is a supercomputer with tens of thousands of processor cores. To push the limit of CFD in all four fronts, in this paper, we develop and study a highly parallel algorithm for solving a monolithically coupled fluid-structure system for the modeling of the interaction of the blood flow and the arterial wall. As a case study, we consider a patient-specific, full size pulmonary artery obtained from CT (Computed Tomography) images, with an artificially added layer of wall with a fixed thickness. The fluid is modeled with a system of incompressible Navier-Stokes equations and the wall is modeled by a geometrically nonlinear elasticity equation. As far as we know this is the first time the unsteady blood flow in a full pulmonary artery is simulated without assuming a rigid wall. We report the proposed numerical algorithm and software scale well beyond 10,000 processor cores on a supercomputer for solving the fluid-structure interaction problem discretized with a stabilized finite element method in space and an implicit scheme in time involving hundreds of millions of unknowns.

97 MATHEMATICS AND COMPUTING↗