Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Network simulation”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Reynolds averaged turbulence modelling using deep neural networks with embedded invariance

There exists significant demand for improved Reynolds-averaged Navier–Stokes (RANS) turbulence models that are informed by and can represent a richer set of turbulence physics. This paper presents a method of using deep neural networks to learn a model for the Reynolds stress anisotropy tensor from high-fidelity simulation data. A novel neural network architecture is proposed which uses a multiplicative layer with an invariant tensor basis to embed Galilean invariance into the predicted anisotropy tensor. It is demonstrated that this neural network architecture provides improved prediction accuracy compared with a generic neural network architecture that does not embed this invariance property. Furthermore, the Reynolds stress anisotropy predictions of this invariant neural network are propagated through to the velocity field for two test cases. For both test cases, significant improvement versus baseline RANS linear eddy viscosity and nonlinear eddy viscosity models is demonstrated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Machine learning of hidden variables in multiscale fluid simulation

Solving fluid dynamics equations often requires the use of closure relations that account for missing microphysics. For example, when solving equations related to fluid dynamics for systems with a large Reynolds number, sub-grid effects become important and a turbulence closure is required, and in systems with a large Knudsen number, kinetic effects become important and a kinetic closure is required. By adding an equation governing the growth and transport of the quantity requiring the closure relation, it becomes possible to capture microphysics through the introduction of ‘hidden variables’ that are non-local in space and time. The behavior of the ‘hidden variables’ in response to the fluid conditions can be learned from a higher fidelity or ab-initio model that contains all the microphysics. In our study, a partial differential equation simulator that is end-to-end differentiable is used to train judiciously placed neural networks against ground-truth simulations. We show that this method enables an Euler equation based approach to reproduce non-linear, large Knudsen number plasma physics that can otherwise only be modeled using Boltzmann-like equation simulators such as Vlasov or particle-in-cell modeling.

42 ENGINEERING↗
DOE OSTI2023

Elucidating the water–anatase TiO 2 (101) interface structure using infrared signatures and molecular dynamics

The structure and dynamics of water on solid surfaces critically affect the chemistry of materials in ambient and aqueous environments. Here, we investigate the hydrogen bonding network of water adsorbed on the majority (101) surface of anatase TiO 2 , a widely used photocatalyst, using polarization- and azimuth-resolved infrared spectroscopy combined with neural network potential molecular dynamics simulations. Our results show that one monolayer of water saturates the undercoordinated titanium (Ti 5c ) sites, forming one-dimensional chains of molecule hydrogen bonded to surface undercoordinated bridging oxygen (O 2c ) atoms. As the coverage increases, water adsorption on O 2c sites leads to significant restructuring of the water monolayer and the formation of a two-dimensional hydrogen bond network characterized by tightly bound pairs of water molecules on adjacent Ti 5c and O 2c sites. This structural motif likely persists at ambient conditions, influencing the reactions occurring there. In conclusion, the results reported here provide critical details of the structure of the water–anatase (101) interface that were previously hypothesized but unconfirmed experimentally.

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

Predicting Local Material Thickness from Steady-State Ultrasonic Wavefield Measurements Using a Convolutional Neural Network

Acoustic steady-state excitation spatial spectroscopy (ASSESS) is a full-field, ultrasonic non-destructive evaluation (NDE) technique used to locate and characterize defects in plate-like structures. ASSESS generates a steady-state, single-tone ultrasonic excitation in a structure and a scanning laser Doppler Vibrometer (LDV) measures the resulting full-field surface velocity response. Traditional processing techniques for ASSESS data rely on wavenumber domain analysis. This paper presents the alternative use of a convolutional neural network (CNN), trained using simulated ASSESS data, to predict the local plate thickness at every pixel in the wavefield measurement directly. The defect detection accuracy of CNN-based thickness predictions are shown to improve for defects of greater size, and for defects with higher thickness reductions. The CNN demonstrates the ability to predict thickness accurately in regions where Lamb wave dispersion relations are complex or unknown, such as near the boundaries of a test specimen, so long as the CNN is trained on data that accounts for these regions. The CNN also shows generalizability to ASSESS experimental data, despite an entirely simulated training dataset.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Solving Newton’s equations of motion with large timesteps using recurrent neural networks based operators

Classical molecular dynamics simulations are based on solving Newton’s equations of motion. Using a small timestep, numerical integrators such as Verlet generate trajectories of particles as solutions to Newton’s equations. We introduce operators derived using recurrent neural networks that accurately solve Newton’s equations utilizing sequences of past trajectory data, and produce energy-conserving dynamics of particles using timesteps up to 4000 times larger compared to the Verlet timestep. We demonstrate significant speedup in many example problems including 3D systems of up to 16 particles.

Newton’s equations↗
DOE OSTI2018

Topological vulnerability of power grids to disasters: Bounds, adversarial attacks and reinforcement

Natural disasters like hurricanes, floods or earthquakes can damage power grid devices and create cascading blackouts and islands. The nature of failure propagation and extent of damage, among other factors, is dependent on the structural features of the grid, that are distinct from that of random networks. This paper analyzes the structural vulnerability of real power grids to impending disasters and presents intuitive graphical metrics to quantify the extent of topological damage. We develop two improved graph eigen-value based bounds on the grid vulnerability. Further we study adversarial attacks aimed at weakening the grid’s structural robustness and present three combinatorial algorithms to determine the optimal topological attack. Simulations on power grid networks and comparison with existing work show the improvements of the proposed measures and attack schemes.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2014

Interstellar simulations using a unified microscopic-macroscopic Monte Carlo model with a full gas-grain network including bulk diffusion in ice mantles

We have designed an improved algorithm that enables us to simulate the chemistry of cold dense interstellar clouds with a full gas-grain reaction network. The chemistry is treated by a unified microscopic-macroscopic Monte Carlo approach that includes photon penetration and bulk diffusion. To determine the significance of these two processes, we simulate the chemistry with three different models. In Model 1, we use an exponential treatment to follow how photons penetrate and photodissociate ice species throughout the grain mantle. Moreover, the products of photodissociation are allowed to diffuse via bulk diffusion and react within the ice mantle. Model 2 is similar to Model 1 but with a slower bulk diffusion rate. A reference Model 0, which only allows photodissociation reactions to occur on the top two layers, is also simulated. Photodesorption is assumed to occur from the top two layers in all three models. We found that the abundances of major stable species in grain mantles do not differ much among these three models, and the results of our simulation for the abundances of these species agree well with observations. Likewise, the abundances of gas-phase species in the three models do not vary. However, the abundances of radicals in grain mantles can differ by up to two orders of magnitude depending upon the degree of photon penetration and the bulk diffusion of photodissociation products. We also found that complex molecules can be formed at temperatures as low as 10 K in all three models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Sprites and State Channels: Payment Networks that Go Faster than Lightning

Bitcoin, Ethereum and other blockchain-based cryptocurrencies, as deployed today, cannot scale for wide-spread use. A leading approach for cryptocurrency scaling is a smart contract mechanism called a payment channel which enables two mutually distrustful parties to transact efficiently (and only requires a single transaction in the blockchain to set-up). Payment channels can be linked together to form a payment network, such that payments between any two parties can (usually) be routed through the network along a path that connects them. Crucially, both parties can transact without trusting hops along the route. In this paper, we propose a novel variant of payment channels, called Sprites, that reduces the worst-case “collateral cost” that each hop along the route may incur. The benefits of Sprites are two-fold. 1) In Lightning Network, a payment across a path of ℓ channels requires locking up collateral for Θ(ℓΔ) time, where Δ is the time to commit an on-chain transaction. Sprites reduces this cost to Θ(ℓ+Δ). 2) Unlike prior work, Sprites supports partial withdrawals and deposits, during which the channel can continue to operate without interruption. In evaluating Sprites we make several additional contributions. First, our simulation-based security model is the first formalism to model timing guarantees in payment channels. Our construction is also modular, making use of a generic abstraction from folklore, called the “state channel,” which we are the first to formalize. We also provide a simulation framework for payment network protocols, which we use to confirm that the Sprites construction mitigates against throughput-reducing attacks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Resonant interaction between phonons and PbTe/PbSe (001) misfit dislocation networks

Here this work aims at a quantitative and mechanistic understanding of the dynamic process of the phonon-dislocation interaction in PbTe/PbSe (001) heterostructures using the Concurrent Atomistic-Continuum (CAC) method as the simulation tool. The misfit dislocation network and the atomic-scale dislocation core structure obtained in the simulations are found to agree reasonably well with the experimental observations of the PbTe/PbSe (001) interface. Through visualizing the dynamic interaction between phonons and dislocations, as well as quantifying the dislocation vibration amplitude, the phonon energy transmission, and the thermal resistance of the misfit interfaces, this work has illustrated and quantified two mechanisms for phonon-dislocation interaction: (1) phonon scattering by the strain field of dislocations, and (2) phonon scattering by dislocations that vibrate via the local modes of a dislocation network; the latter, leads to resonant phonon-dislocation interaction, which is manifested as local maxima of out-of-phase vibration of the atoms on the two sides of the slip plane, leading to local minima of the energy transmission in the heterostructure that contains one interface. The local vibrational modes are found to be excited only by shear stress induced by transverse phonons. Among various resonant modes, the one with the lowest frequency has the strongest effect. This work has also demonstrated the collective motion of dislocations under ultrafast phonon pulses. In addition, the dynamic properties of the misfit dislocation network localized within one interface are found to be significantly altered by the presence of misfit dislocations at other interfaces, thus further confirming the cooperative dynamic nature of the motion of dislocations and phonons.

36 MATERIALS SCIENCE↗
DOE OSTI2017

SQERTSS: Dynamic rank based throttling of transition probabilities in kinetic Monte Carlo simulations

Lattice based Kinetic Monte Carlo (KMC) simulations offer a powerful simulation technique for investigating large reaction networks while retaining spatial configuration information, unlike ordinary differential equations. However, large chemical reaction networks can contain reaction processes with rates spanning multiple orders of magnitude. This can lead to the problem of “KMC stiffness” (similar to stiffness in differential equations), where the computational expense has the potential to be overwhelmed by very short time-steps during KMC simulations, with the simulation spending an inordinate amount of KMC steps / cpu-time simulating fast frivolous processes (FFPs) without progressing the system (reaction network). In order to achieve simulation times that are experimentally relevant or desired for predictions, a dynamic throttling algorithm involving separation of the processes into speed-ranks based on event frequencies has been designed and implemented with the intent of decreasing the probability of FFP events, and increasing the probability of slow process events -- allowing rate limiting events to become more likely to be observed in KMC simulations. This Staggered Quasi-Equilibrium Rank-based Throttling for Steady-state (SQERTSS) algorithm designed for use in achieving and simulating steady-state conditions in KMC simulations. Lastly, as shown in this work, the SQERTSS algorithm also works for transient conditions: the correct configuration space and final state will still be achieved if the required assumptions are not violated, with the caveat that the sizes of the time-steps may be distorted during the transient period.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Tailoring the mechanical properties of nanoparticle networks that encompass biomimetic catch bonds

Biological “catch” bonds display the distinctive attribute that the bond lifetime can increase under an applied force. Insertion of biomimetic catch bonds into hybrid materials could lead to composites that exhibit remarkable mechanical properties. We model the tensile behavior of polymer-grafted nanoparticle (PGN) networks interconnected by a mixture of catch bonds and conventional “slip” bonds, whose lifetimes decrease with force. We formulate a kinetic master equation that provides the complete probabilistic description of a system of two PGNs. The master equation is used to analyze the parameter space that determines the rupture behavior of the catch bonds in the two-particle system. We then utilize two exemplary sets of catch bond parameters in three-dimensional computer simulations of larger PGN networks under strain-controlled tensile deformation. We demonstrate that the strain at break and toughness of the networks can be altered by “tuning” the attributes of the catch bonds and varying the fraction of catch bonds in the network. We show that networks encompassing the catch bonds could exhibit a several-fold increase in the strain-at-break and toughness relative to those interconnected solely by the slip bonds. Our studies can provide valuable guidelines for tailoring the mechanical properties of novel, bio-inspired nanocomposites.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Oscillatory Spreading and Inertia in Power Grids

The increase in variable renewable generators (VRGs) in power systems has altered the dynamics from a historical experience. VRGs introduce new sources of power oscillations, and the stabilizing response provided by synchronous generators (SGs, e.g., natural gas, coal, etc.), which help avoid some power fluctuations, will lessen as VRGs replace SGs. These changes have led to the need for new methods and metrics to quickly assess the likely oscillatory behavior for a particular network without performing computationally expensive simulations. This work studies the impact of a critical dynamical parameter - the inertia value - on the rest of a power system's oscillatory response to representative VRG perturbations. We use a known localization metric in a novel way to quantify the number of nodes responding to a perturbation and the magnitude of those responses. This metric allows us to relate the spread and severity of a system's power oscillations with inertia. We find that as inertia increases, the system response to node perturbations transitions from localized (only a few close nodes respond) to delocalized (many nodes across the network respond). We introduce a heuristic computed from the network Laplacian to relate this oscillatory transition to the network structure. We show that our heuristic accurately describes the spread of oscillations for a realistic power-system test case. Using a heuristic to determine the likely oscillatory behavior of a system given a set of parameters has wide applicability in power systems, and it could decrease the computational workload of planning and operation.

dynamical systems↗
DOE OSTI2025

Three-dimensional modeling of hyphal fusion, branching, and nutrient transport in filamentous fungi

Fungi exhibit behaviors distinct from other microbes. Filamentous fungi grow by extending complex networks of branched filaments collectively referred to as the mycelium. These networks can expand over large distances and traverse low-nutrient areas by translocating nutrients through the filament network. This spatial characteristic makes filamentous fungi crucial for soil ecosystems, supporting stable microbial communities and promoting plant growth. However, simulating these behaviors is complex. The elongated nature of fungal compartments results in different mechanical interactions compared to the commonly modeled spherical bacteria. These detailed hyphal mechanics require specialized consideration and are often excluded from conventional fungal simulation packages. Additionally, the extensive fungal networks in nature demand computationally intensive simulations, necessitating high-performance algorithms. Therefore, realistic fungi simulations require specialized software. Here, we introduce a fungal modeling expansion to the high-performance biological modelling and interface exchange (bmx) software suite. bmx leverages adaptive mesh refinement in AMReX for chemical diffusion and incorporates a full mechanical model for bacterial cells, accelerated by GPUs. By extending bmx to model filamentous particles, we demonstrate the formation of complex filament networks through interactions like hyphal branching and fusion (anastomosis). We show that the networks produced match real-world fungal structures through various metrics. This work supports computational studies of fungal growth dynamics and can be adapted to investigate the growth of other filamentous structures in biology or materials science. The expanded-BMX package is open-sourced and is available online.

Cell mechanics↗
DOE OSTI2024

Electromagnetic shower separation in LArTPCs: study of the calorimetric method at LArIAT

This thesis reports a study of the prospects of using calorimetry for separating photon- from electron-initiated electromagnetic (EM) showers in liquid argon time projection chambers (LArTPCs), a prevalent need for neutrino oscillation experiments based on the LArTPC technology as $\gamma$-initiated showers are a significant background to electron neutrino ($\nu_e$) events, whose signature is electron showers. From events simulated in the LArIAT (Liquid Argon In A Testbeam) experiment's detector, we quantitatively evaluated the dependence of the calorimetric separation on the initial momentum and position (along the beamline) of the primary particle. We could conclude that the separation is more efficient if the primary particle is created inside the chamber, whereas the impact of the momentum is not as significant. Furthermore, a first step was taken towards applying the calorimetric method to the showers measured in LArIAT. A convolutional neural network (CNN) trained on simulation was deployed to find EM showers among the events collected by the experiment. The charge deposition at the beginning of the showers found is centered at a plausible value, although the width of the distribution demands further inquiry.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Blending Pipeline Analysis Tool for Hydrogen (BlendPATH) Documentation and User Manual

The Blending Pipeline Analysis Tool for Hydrogen (BlendPATH) is a flexible, open-source Python tool designed to provide users with case-by-case analysis capabilities to identify the necessary modifications to repurpose existing natural gas transmission pipeline networks to transport hydrogen as a blend of a user-specified volume fraction of hydrogen or as a pure stream and to estimate the associated capital and operating expenditures resulting from those modifications. This tool is intended to be applied during the initial screening stage of a prospective project when pipeline developers compile transmission pipeline technical documentation and history but prior to performing detailed pipeline inspections. Performing analysis with BlendPATH during this initial screening stage can provide the user with an understanding of promising opportunities and probable economic outcomes of repurposing their pipeline network for hydrogen before proceeding with detailed pipeline materials testing and pipeline inspections. BlendPATH consists of multiple modules to simulate, assess, and modify existing natural gas transmission pipeline network designs to be compatible with hydrogen as a blend or pure stream. The tool employs an open-source gas network hydraulic model to simulate existing, user-specified transmission pipeline networks, and it applies ASME B31.12 to assess the pipe segments within these networks for compatibility with hydrogen and to modify the networks to achieve compatibility where the existing infrastructure is inadequate. Users can specify ASME B31.12 design options for pipeline assessment and can select from multiple methods for pipeline modification. This report details the functionalities of BlendPATH Version 2.0.2 and demonstrates an example of applying BlendPATH to a case study. Potential and intended users of this framework include natural gas pipeline developers and operators and researchers at both public and private institutions. BlendPATH is publicly available at https://github.com/NREL/BlendPATH.

08 HYDROGEN↗
DOE OSTI2018

Literature Review on Modeling Cyber Networks and Evaluating Cyber Risks

The National Infrastructure Simulations and Analysis Center (NISAC) conducted a literature review on modeling cyber networks and evaluating cyber risks. The literature review explores where modeling is used in the cyber regime and ways that consequence and risk are evaluated. The relevant literature clusters in three different spaces: network security, cyber-physical, and mission assurance. In all approaches, some form of modeling is utilized at varying levels of detail, while the ability to understand consequence varies, as do interpretations of risk. This document summarizes the different literature viewpoints and explores their applicability to securing enterprise networks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Seasonal and Morphological Controls on Nitrate Retention in Arctic Deltas

Abstract Estimates of nitrate loading to the Arctic Ocean are limited by the lack of field observations within deltas partly due to logistical constraints. To overcome this limitation, we use a remote sensing framework to estimate retention of nitrate in Arctic deltas. We achieve this by coupling hydrological and biogeochemical process models at the network scale for five major Arctic deltas. Binary masks of delta channels were used to simulate flow direction and magnitude through networks. Models were parameterized using historical and seasonal observations. Simulated nitrate retention ranged from 2.9% to 15% of the incoming load. Retention rates were largest during winter but smallest during spring conditions when increased discharges export large nitrate masses to the coast. Under future climate scenarios, retention rates fall by ∼1%–10%. Arctic deltas have an important effect on the magnitude of nitrate entering Arctic seas and the inclusion of processing in deltas can improve flux estimates.

54 ENVIRONMENTAL SCIENCES↗