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 1,063 records · Page 59

Topography-based Flood Planning and Optimization Capability Development Report

Globally, water-related disasters are among the most frequent and costly natural hazards. Flooding inflicts catastrophic damage on critical infrastructure and population, resulting in substantial economic and social costs. NISAC is developing LeveeSim, a suite of nonlinear and network optimization models, to predict optimal barrier placement to protect critical regions and infrastructure during flood events. LeveeSim currently includes a high-performance flood model to simulate overland flow, as well as a network optimization model to predict optimal barrier placement during a flood event. The LeveeSim suite models the effects of flooding in predefined regions. By manipulating a domain’s underlying topography, developers altered flood propagation to reduce detrimental effects in areas of interest. This numerical altering of a domain’s topography is analogous to building levees, placing sandbags, etc. To induce optimal changes in topography, NISAC used a novel application of an optimization algorithm to minimize flooding effects in regions of interest. To develop LeveeSim, NISAC constructed and coupled hydrodynamic and optimization algorithms. NISAC first implemented its existing flood modeling software to use massively parallel graphics processing units (GPUs), which allowed for the simulation of larger domains and longer timescales. NISAC then implemented a network optimization model to predict optimal barrier placement based on output from flood simulations. As proof of concept, NISAC developed five simple test scenarios, and optimized topographic solutions were compared with intuitive solutions. Finally, as an early validation example, barrier placement was optimized to protect an arbitrary region in a simulation of the historic Taum Sauk dam breach.

42 ENGINEERING↗

The behaviour of basic autocatalytic signalling modules in isolation and embedded in networks

In this paper, we examine the behaviour of basic autocatalytic feedback modules involving a species catalyzing its own production, either directly or indirectly. We first perform a systematic study of the autocatalytic feedback module in isolation, examining the effect of different factors, showing how this module is capable of exhibiting monostable threshold and bistable switch-like behaviour. We then study the behaviour of this module embedded in different kinds of basic networks including (essentially) irreversible cycles, open and closed reversible chains, and networks with additional feedback. We study the behaviour of the networks deterministically and also stochastically, using simulations, analytical work, and bifurcation analysis. We find that (i) there are significant differences between the behaviour of this module in isolation and in a network: thresholds may be altered or destroyed and bistability may be destroyed or even induced, even when the ambient network is simple. The global characteristics and topology of this network and the position of the module in the ambient network can play important and unexpected roles. (ii) There can be important differences between the deterministic and stochastic dynamics of the module embedded in networks, which may be accentuated by the ambient network. This provides new insights into the functioning of such enzymatic modules individually and as part of networks, with relevance to other enzymatic signalling modules as well.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantum Communication Networks for Energy Applications: Review and Perspective

Abstract The energy sector is expected to undergo significant changes in the coming decades with the advent of new technologies, including smart grid development, microgrid expansion, increasing electric vehicle and renewable energy usage, and enhanced measures to minimize greenhouse gas emission, among others. In tandem, these changes are expected to create new opportunities for the deployment of quantum technologies within the energy sector. Building on the authors' previous reviews on the current state of and future opportunities for quantum sensing, quantum computing and quantum simulations for energy sector applications, this work provides an overview of recent progress in quantum networking and communications for the energy industry, with a focus on platforms, devices, and protocols, including quantum teleportation and quantum key distribution. Specific areas of relevance to the energy sector are then analyzed, including the role of quantum networks for greenhouse gas monitoring, secure data collection and transmission in smart grids, nuclear power plants’ safety, facilitating oil and gas exploration, and other energy‐relevant applications. This review concludes with a brief overview of areas for future innovation, including the need for platforms for simulating quantum networks, quantum material and platform design, and computational approaches to accelerate quantum protocol discovery and development.

Paudel, Hari P.↗

SU-E-J-163: A Biomechanical Lung Model for Respiratory Motion Study

Purpose: This work presents a biomechanical model to investigate the complex respiratory motion for the lung tumor tracking in radiosurgery by computer simulation. Methods: The models include networked massspring-dampers to describe the tumor motion, different types of surrogate signals, and the force generated by the diaphragm. Each mass-springdamper has the same mechanical structure and each model can have different numbers of mass-spring-dampers. Both linear and nonlinear stiffness parameters were considered, and the damping ratio was tuned in a range so that the tumor motion was over-damped (no natural tumor oscillation occurs without force from the diaphragm). The simulation was run by using ODE45 (ordinary differential equations by Runge-Kutta method) in MATLAB, and all time courses of motions and inputs (force) were generated and compared. Results: The curvature of the motion time courses around their peaks was sensitive to the damping ratio. Therefore, the damping ratio can be determined based on the clinical data of a high sampling rate. The peak values of different signals and the time the peaks occurred were compared, and it was found that the diaphragm force had a time lead over the tumor motion, and the lead time (0.1–0.4 seconds) depended on the distance between the tumor and the diaphragm. Conclusion: We reported a model based analysis approach for the spatial and temporal relation between the motion of the lung tumor and the surrogate signals. Due to the phase lead of the diaphragm in comparing with the lung tumor motion, the measurement of diaphragm motion (or its electromyography signal) can be used as a beam gating signal in radiosurgery, and it can also be an additional surrogate signal for better tumor motion tracking. The research is funded by the American Cancer Society (ACS) grant. The grant name is: Frameless SRS Based on Robotic Head Motion Cancellation. The grant number is: RSG-13-313-01-CCE.

60 APPLIED LIFE SCIENCES↗

Optical Abelian lattice gauge theories

We discuss a general framework for the realization of a family of Abelian lattice gauge theories, i.e., link models or gauge magnets, in optical lattices. We analyze the properties of these models that make them suitable for quantum simulations. Within this class, we study in detail the phases of a U(1)-invariant lattice gauge theory in 2+1 dimensions, originally proposed by P. Orland. By using exact diagonalization, we extract the low-energy states for small lattices, up to 4 Multiplication-Sign 4. We confirm that the model has two phases, with the confined entangled one characterized by strings wrapping around the whole lattice. We explain how to study larger lattices by using either tensor network techniques or digital quantum simulations with Rydberg atoms loaded in optical lattices, where we discuss in detail a protocol for the preparation of the ground-state. We propose two key experimental tests that can be used as smoking gun of the proper implementation of a gauge theory in optical lattices. These tests consist in verifying the absence of spontaneous (gauge) symmetry breaking of the ground-state and the presence of charge confinement. We also comment on the relation between standard compact U(1) lattice gauge theory and the model considered in this paper. - Highlights: Black-Right-Pointing-Pointer We study the quantum simulation of dynamical gauge theories in optical lattices. Black-Right-Pointing-Pointer We focus on digital simulation of abelian lattice gauge theory. Black-Right-Pointing-Pointer We rediscover and discuss the puzzling phase diagram of gauge magnets. Black-Right-Pointing-Pointer We detail the protocol for time evolution and ground-state preparation in any phase. Black-Right-Pointing-Pointer We provide two experimental tests to validate gauge theory quantum simulators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Post-Fracturing Experiment Simulation of Hydraulic Fracture Propagation in a Deep Mine Using a Fully Coupled 3D Network-Flow and Quasi-static Discrete Element Model

Aimed at supporting the hydraulic fracturing experiments in a vertical borehole in the Poorman phyllite at approximately 1500 m depth at the Sanford Underground Research Facility (SURF), we perform post-hydraulic fracturing experiment simulations and compare the simulated borehole pressure changes and flow rates with actual measurements using a newly developed fully coupled three-dimensional (3D) network flow and quasi-static discrete element model (DEM). The quasi-static DEM model, which is constructed by Delaunay tessellation of the rock volume, considers rock fabric heterogeneities by using the “disordered” DEM mesh and adding random perturbations to the stiffness and tensile/shear strengths of individual DEM elements and the elastic beam between them. A conjugate 3D flow network based on DEM lattice is constructed to calculate the fluid flow in both the fracture and porous matrix. One distinctive advantage of the DEM model is that fracturing is naturally described by the breakage of elastic beams between DEM elements. It is also extremely convenient to introduce mechanical anisotropy into the model by simply assigning orientation-dependent tensile/shear strengths to the elastic beams.

15 GEOTHERMAL ENERGY↗

Two-dimensional isometric tensor networks on an infinite strip

The exact contraction of a generic two-dimensional (2D) tensor network state (TNS) is known to be exponentially hard, making simulation of 2D systems difficult. The recently introduced class of isometric TNS (isoTNS) represents a subset of TNS that allows for efficient simulation of such systems on finite square lattices. The isoTNS ansatz requires the identification of an “orthogonality column” of tensors, within which one-dimensional matrix product state (MPS) methods can be used for calculation of observables and optimization of tensors. Here we extend isoTNS to infinitely long strip geometries and introduce an infinite version of the Moses Move algorithm for moving the orthogonality column around the network. Using this algorithm, we iteratively transform an infinite MPS representation of a 2D quantum state into a strip isoTNS and investigate the entanglement properties of the resulting state. In addition, we demonstrate that the local observables can be evaluated efficiently. Lastly, we introduce an infinite time-evolving block decimation algorithm (iTEBD 2 ) and use it to approximate the ground state of the 2D transverse field Ising model on lattices of infinite strip geometry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Fracture Intensity Impacts on Reaction Front Propagation and Mineral Weathering in Three‐Dimensional Fractured Media

Abstract Studying reaction front propagation in heterogeneous natural settings is challenging, but numerical simulations can provide insight into the varying spatial and temporal scales of reaction front propagation. Here, the impact of increasing fracture intensity on mineral dissolution rates, and the extent of reaction front propagation is investigated using reactive transport simulations in upscaled discrete fracture network domains with varied fracture intensity. Domain‐averaged dissolution rates vary less than 0.5 log units regardless of the fracture intensity, but the spatial distribution of reactions is controlled by the location and number of dead‐end fractures and the number of connected flowpaths through the domain. Higher fracture intensities lead to more weathering in the domain because of more available mineral for water‐rock interactions. We find that reaction fronts propagate through the primary flowpaths in the first 10,000 years of the simulation for a 10‐m length domain, then propagate into secondary flowpaths and dead‐end fractures between 10,000 and 100,000 years, and finally into the matrix over timescales of hundreds of thousands of years. The domain‐averaged reaction rates decrease through time corresponding to a transition from dissolution in advection‐dominated, fast‐flowing pathways, to dissolution in transport‐limited zones of disconnected fractures and matrix. Matrix dissolution, or dissolution under transport‐limited conditions, is the dominant process at late times in these simulations. The results of these simulations recreate the observed paradox found in nature where highly fractured hillslopes tend to be more weathered but have slower weathering rates, while hillslopes with fewer fractures, are less weathered but have higher dissolution rates.

58 GEOSCIENCES↗

ASTROCHEMICAL CORRELATIONS IN MOLECULAR CLOUDS

We investigate the spectral correlations between different species used to observe molecular clouds. We use hydrodynamic simulations and a full chemical network to study the abundances of over 150 species in typical Milky Way molecular clouds. We perform synthetic observations in order to produce emission maps of a subset of these tracers. We study the effects of different lines of sight and spatial resolution on the emission distribution and perform a robust quantitative comparison of the species to each other. We use the Spectral Correlation Function (SCF), which quantifies the root mean squared difference between spectra separated by some length scale, to characterize the structure of the simulated cloud in position-position-velocity (PPV) space. We predict the observed SCF for a broad range of observational tracers, and thus identify homologous species. In particular, we show that the pairs C and CO, C{sup +} and CN, and NH{sub 3} and H{sub 2}CS have very similar SCFs. We measure the SCF slope variation as a function of beam size for all species and demonstrate that the beam size has a distinct effect on different species emission. However, for beams of up to 10'', placing the cloud at 1 kpc, the change is not large enough to move the SCF slopes into different regions of parameter space. The results from this study provide observational guidance for choosing the best tracer to probe various cloud length scales.

79 ASTRONOMY AND ASTROPHYSICS↗

Data-scarce surrogate modeling of shock-induced pore collapse process

Understanding the mechanisms of shock-induced pore collapse is of great interest in various disciplines in sciences and engineering, including materials science, biological sciences, and geophysics. However, numerical modeling of the complex pore collapse processes can be costly. To this end, a strong need exists to develop surrogate models for generating economic predictions of pore collapse processes. Here, in this work, we study the use of a data-driven reduced-order model, namely dynamic mode decomposition, and a deep generative model, namely conditional generative adversarial networks, to resemble the numerical simulations of the pore collapse process at representative training shock pressures. Since the simulations are expensive, the training data are scarce, which makes training an accurate surrogate model challenging. To overcome the difficulties posed by the complex physics phenomena, we make several crucial treatments to the plain original form of the methods to increase the capability of approximating and predicting the dynamics. In particular, physics information is used as indicators or conditional inputs to guide the prediction. In realizing these methods, the training of each dynamic mode composition model takes only around 30 s on CPU. In contrast, training a generative adversarial network model takes 8 h on GPU. Moreover, using dynamic mode decomposition, the final-time relative error is around 0.3% in the reproductive cases. We also demonstrate the predictive power of the methods at unseen testing shock pressures, where the error ranges from 1.3 to 5% in the interpolatory cases and 8 to 9% in extrapolatory cases.

97 MATHEMATICS AND COMPUTING↗

Computationally efficient CFD prediction of bubbly flow using physics-guided deep learning

To realize efficient computational fluid dynamics (CFD) prediction of two-phase flow, a multi-scale framework was proposed in this paper by applying a physics-guided data-driven approach. Instrumental to this framework, Feature Similarity Measurement (FSM) technique was developed for error estimation in two-phase flow simulation using coarse-mesh CFD, to achieve a comparable accuracy as fine-mesh simulations with fast-running feature. In this work, by defining physics-guided parameters and variable gradients as physical features, FSM has the capability to capture the underlying local patterns in the coarse-mesh CFD simulation. Massive low-fidelity data and respective high-fidelity data are used to explore the underlying information relevant to the main simulation errors and the effects of phenomenological scaling. By learning from previous simulation data, a surrogate model using deep feedforward neural network (DFNN) can be developed and trained to estimate the simulation error of coarse-mesh CFD. In a demonstration case of two-phase bubbly flow, the DFNN model well captured and corrected the unphysical “peaks” in the velocity and void fraction profiles near the wall in the coarse-mesh configuration, even for extrapolative predictions. The research documented supports the feasibility of the physics-guided deep learning methods for coarse mesh CFD simulations which has a potential for the efficient industrial design.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Grid-Forming Frequency Shaping Control for Low-Inertia Power Systems

As power systems transit to a state of high renewable penetration, little or no presence of synchronous generators makes the prerequisite of well-regulated frequency for grid-following inverters unrealistic. Thus, there is a trend to resort to grid-forming inverters which set frequency directly. We propose a novel grid-forming frequency shaping control that is able to shape the aggregate system frequency dynamics into a first-order one with the desired steady-state frequency deviation and Rate of Change of Frequency (RoCoF) after a sudden power imbalance. The no overshoot property resulting from the first-order dynamics allows the system frequency to monotonically move towards its new steady-state without experiencing frequency Nadir, which largely improves frequency security. We prove that our grid-forming frequency-shaping control renders the system internally stable under mild assumptions. The performance of the proposed control is verified via numerical simulations on a modified Icelandic Power Network test case.

control systems↗

DEMONSTRATION OF A DATA-DRIVEN PHYSICS-BASED APPROACH FOR COMPUTATIONALLY EFFICIENT CFD PREDICTION OF TWO-PHASE BUBBLY FLOW

To realize efficient computational fluid dynamics (CFD) prediction of two-phase flow, a multi-scale physics-guided data-driven approach, Feature Similarity Measurement (FSM) technique was developed for error estimation in two-phase flow simulation using coarse-mesh CFD, to achieve a comparable accuracy as fine-mesh simulations with fast-running feature. By defining physics-guided parameters and variable gradients as physical features, FSM has the capability to capture the underlying local patterns in the coarse-mesh CFD simulation. Massive low-fidelity data and respective high-fidelity data are used to explore the underlying information relevant to the main simulation errors and the effects of phenomenological scaling. By learning from previous simulation data, a surrogate model using deep feedforward neural network (DFNN) can be developed and trained to estimate the simulation error of coarse-mesh CFD. In a demonstration case of two-phase bubbly flow, the DFNN model well captured and corrected the unphysical “peaks” in the velocity and void fraction profiles near the wall in the coarse-mesh configuration, even for extrapolative predictions. The research documented supports the feasibility of the physics-guided deep learning methods for coarse mesh CFD simulations which has a potential for the efficient industrial design.

42 ENGINEERING↗

On the Network Topology of Cross-Linked Acrylate Photopolymers: A Molecular Dynamics Case Study

A reactive molecular dynamics approach is used to simulate cross-linking of acrylate polymer networks. By employing the same force field and reactive scheme and studying three representative multifunctional acrylate monomers, we isolate the importance of the nonreactive moieties within these model monomers. Analyses of reactive trajectories benchmark the estimated gel points, cyclomatic character, and spatially resolved cross-linking tendencies of the acrylates as a function of conversion. The insights into the similarities and differences of the polymerization and resulting networks presented here suggest molecular mechanics as a useful tool in the rational design of photopolymerization resins.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Collective neural network behavior in a dynamically driven disordered system of superconducting loops

Collective properties of complex systems composed of many interacting components such as neurons in our brain can be modeled by artificial networks based on disordered systems. We show that a disordered neural network of superconducting loops with Josephson junctions can exhibit computational properties like categorization and associative memory in the time evolution of its state in response to information from external excitations. Superconducting loops can trap multiples of fluxons in many discrete memory configurations defined by the local free energy minima in the configuration space of all possible states. A memory state can be updated by exciting the Josephson junctions to fire or allow the movement of fluxons through the network as the current through them surpasses their critical current thresholds. Simulations performed with a lumped element circuit model of a 4-loop network show that information written through excitations is translated into stable states of trapped flux and their time evolution. Experimental implementation on a high-Tc superconductor YBCO-based 4-loop network shows dynamically stable flux flow in each pathway characterized by the correlations between junction firing statistics. Neural network behavior is observed as energy barriers separating state categories in simulations in response to multiple excitations, and experimentally as junction responses characterizing different flux flow patterns in the network. The state categories that produce these patterns have different temporal stabilities relative to each other and the excitations. This provides strong evidence for time-dependent (short-to-long-term) memories, that are dependent on the geometrical and junction parameters of the loops, as described with a network model.

Josephson junctions↗

Social vulnerability and power loss mitigation: A case study of Puerto Rico

The increasing occurrence of extreme weather events urges us to reevaluate the resiliency and vulnerability aspects of our most critical infrastructures — such as power grids — as their failures result in both economic loss and severe human hardship. Seen through the lens of alleviating human suffering, it is crucial to be able to identify critical system components of the infrastructure for targeted hardening given resource constraints. This effort is of particular importance in islanded areas such as Puerto Rico where hurricanes are frequent and resources are limited, and where the spatially diverse effects of power loss on human suffering are all the more severe. Recent studies on evaluating infrastructure networks during extreme weather events have taken a simulation based approach that incorporates a variety of component models, such as weather realizations, topological network models, fragility models, and power flow models to estimate expected loss of service. Here, in this work, we expand such a Component Based Event Simulation (CBES) methodology proposed in the literature and integrate it with a social vulnerability modeling component. This paradigm-advancing approach of synthesizing the cutting edge capability of power network modeling and the social impacts of the power transmission network failure is demonstrated for the island of Puerto Rico. Our work exemplifies the efficacy of this integrated modeling framework in developing a decision metric for targeted transmission line hardening.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Enabling Grid-Forming Control with Fault Ride-Through in Unbalanced Distribution Networks

Distribution networks are often unbalanced, causing oscillatory responses in inverter control designed for balanced conditions. Here, to address this problem, this paper proposes a novel time-domain transformation appropriate for inverter control and enables the decomposition of three-phase unbalanced signals into constant positive and negative components. Relations useful for calculating unbalanced active and reactive power are derived from first principle, providing insight into vector products of unbalanced three-phase signals. Furthermore, a grid-forming control effective under unbalanced conditions is developed, which delivers superior performance while meeting UNIFI1 specifications for grid-forming control under unbalanced conditions. specifications applicable to category 4 inverter-based resource, like setting and regulating frequency/voltage, providing voltage support, sharing active power, injecting negative sequence current, and riding through faults. A current limiter is proposed for safe fault ride-through and integrates with the grid-forming control featuring frequency/voltage droop controllers and current and voltage control loops. The transformation of interconnected inverters is formulated and stability of the proposed control analyzed to support robust parameter selections. The effectiveness of the proposed transformation and grid-forming control is demonstrated through analytical results and real-time simulation of a IEEE 123 distribution network on the Real-Time Digital Simulator. Comparison with existing methods shows that the proposed strategy satisfies the UNIFI specifications with a much better performance.

24 - POWER TRANSMISSION AND DISTRIBUTION↗

Snapshot multispectral imaging using a diffractive optical network

Abstract Multispectral imaging has been used for numerous applications in e.g., environmental monitoring, aerospace, defense, and biomedicine. Here, we present a diffractive optical network-based multispectral imaging system trained using deep learning to create a virtual spectral filter array at the output image field-of-view. This diffractive multispectral imager performs spatially-coherent imaging over a large spectrum, and at the same time, routes a pre-determined set of spectral channels onto an array of pixels at the output plane, converting a monochrome focal-plane array or image sensor into a multispectral imaging device without any spectral filters or image recovery algorithms. Furthermore, the spectral responsivity of this diffractive multispectral imager is not sensitive to input polarization states. Through numerical simulations, we present different diffractive network designs that achieve snapshot multispectral imaging with 4, 9 and 16 unique spectral bands within the visible spectrum, based on passive spatially-structured diffractive surfaces, with a compact design that axially spans ~72 λ m , where λ m is the mean wavelength of the spectral band of interest. Moreover, we experimentally demonstrate a diffractive multispectral imager based on a 3D-printed diffractive network that creates at its output image plane a spatially repeating virtual spectral filter array with 2 × 2 = 4 unique bands at terahertz spectrum. Due to their compact form factor and computation-free, power-efficient and polarization-insensitive forward operation, diffractive multispectral imagers can be transformative for various imaging and sensing applications and be used at different parts of the electromagnetic spectrum where high-density and wide-area multispectral pixel arrays are not widely available.

36 MATERIALS SCIENCE↗