Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “network flow”

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,639 records · Page 91

Variations in Bedrock and Vegetation Cover Modulate Subsurface Water Flow Dynamics of a Mountainous Hillslope

Abstract Predicting the hydrological response of watersheds to climate disturbances requires a detailed understanding of the processes connecting hillslopes and streams. Using a network of soil moisture and temperature sensors, electrical resistivity tomography monitoring, and a weather station we assess the above and below‐ground processes driving the hydrological response of a hillslope during snowmelt and summer monsoon. The transect covers bedrock and vegetation gradients, with a steep upper part characterized by shallow bedrock, and gentle lower part underlain by colluvium. The main vegetation cover is conifers on the upper, and grass and veratrum on the lower part. Combined with a simplified hydrological model, we show that the thin soil layer of the steep slope acts as a preferential flow path, leading to mostly shallow lateral flow, interrupted by vertical flow, mostly at tree locations, and likely facilitated by flow along fractures and roots. Vertical flow and upstream‐driven groundwater dynamics are prevailing at the colluvium, presenting a very different hydrological behavior compared to the upper part. These results show that subsurface structure and features have a strong control on the hydrological response of a hillslope and that those can create considerably varying hydrological dynamics across small spatial scales.

54 ENVIRONMENTAL SCIENCES↗

CRN Modeling of Ammonia RQL Combustion using a Partially-Stirred Reactor Approach

Ammonia is a promising alternative to hydrogen with high energy density and favorable storage and transport characteristics. However low flammability and a propensity for high nitrogen oxide (NOx) emissions make direct utilization challenging. Recently, two-stage rich-quench-lean (RQL) combustion strategies have shown promise in achieving low NOx emissions with ammonia. In this approach, the rich stage serves to oxidize a portion of the fuel, while thermally decomposing as much of the remaining ammonia as possible, generating hydrogen. In the second (lean) stage, air is rapidly introduced, burning out the hydrogen and residual ammonia. Two-stage RQL combustion of ammonia has been investigated in the open literature both experimentally and numerically. In general, idealized chemical reactor network (CRN) models predict NOx concentrations below that of 2D/3D computational fluid dynamics models and experiments. The primary drivers of these discrepancies may be largely attributed to finite rate mixing non-adiabatic operation. The typical CRN model is comprised of a perfectly-stirred-reactor (PSR), followed by a plug-flow-reactor (PFR), meant to represent the flame, and post-flame zones, respectively. In the two-stage RQL approach two PSR-PFR networks are arranged sequentially, corresponding to the rich and lean stages, with secondary air injection in between. In the authors’ past work, this arrangement has demonstrated the significant sensitivity of exit NOx to the rich stage equivalence ratio, while the amount of secondary air injection was shown to be less critical. In this paper, the CRN model is extended to (1) include the impacts of heat loss and (2) utilize a partially-stirred-reactor (PaSR) approach to study the impacts of mixing on emissions performance. Varying amounts of heat loss are applied to the rich relaxation zone to understand emissions performance and changes to optimization of equivalence ratio and residence time. Premixed and non-premixed configurations are considered in the rich stage PaSR, with varying degrees of mixing intensity to study the interaction between mixing, transport, and kinetic timescales. Critically, the impact of mixing between hot products and secondary air injection is studied to understand practical injector needs. Results show unburnt ammonia leaving the rich stage as a primary contributor to NOx emissions – driven both by increased heat loss and reduced mixing rates. Furthermore, heat losses have shown to create conditions which are conducive to increased N2O formation in the lean stage. The results of this study will be considered in the context of developing optimized two-stage RQL combustors for ammonia..

advanced gas turbines↗

Coordinative Real-time Sub-Transmission Volt-VAr Control for Reactive Power Regulation between Transmission and Distribution Systems

This paper introduces an integrated voltage regulation approach to coordinate the Volt-VAr control between transmission and distribution systems. Distributed photovoltaic (PV) systems are aggregated as a virtual power plant (VPP). At the sub-transmission level, a Coordinated Real-time Sub-Transmission Volt-VAr Control Tool (CReST-VCT) is introduced to coordinate the operation of the shunt devices and the reactive power provided by the VPPs. The CReST-VCT uses an enhanced Volt-VAr control algorithm with weighted voltage control objects that include minimizing losses, maintaining the voltage within desirable levels, and minimizing the number of mechanical switching of shunt elements. The algorithm runs every 5 minutes and is solved using AC optimal power flow (OPF) solution technique. The outputs include the settings for each sub-transmission-level Volt-VAr control devices and the reactive power generation requirements for each VPP. Then, each VPP will run a distribution Volt-VAr control algorithm to select the PV inverters that provide the Volt-VAr control services considering the distribution network topology. The goal of the VPP reactive power control is to meet the transmission level service requirements while maintaining the voltage and line flow in the distribution system within limits. Each VPP will update its voltage regulation capacity every 5 minutes to allow the CReST-VCT to formulate the optimization problem for the next dispatch interval. The proposed tool is simulated on Duke Energy Carolina system to demonstrate the capability of providing voltage support service by dispatching of reactive power of distributed PV as VPP.

Ke, Xinda↗

Grid Optimization Competition on Synthetic and Industrial Power Systems

This paper summarizes a grid optimization (GO) competition effort in the United States to find the best solution strategies for up to interconnect-scale power system networks with around 32,000 buses. The optimization problem is a mixedinteger, non-convex non-linear problem, (MINLP) and includes discrete variables such as unit commitment and line switching, control settings (transformer taps and phase shifters with impedance correction tables), and bus shunts. The case study includes six actual industry grids as well as 16 realistic synthetic grids created by three different dataset teams. The winners are selected and ranked based on scoring criteria, which consider the solution quality (such as objective functions) within time limits. Nine winner teams are selected from 26 competitor teams. The results achieved by different teams are described and the performance of different algorithms on synthetic grids and actual industry grids are compared and analyzed.

mixed-integer non-linear programming↗

The Airspace Concepts Evaluation System Architecture and System Plant

The Airspace Concepts Evaluation System is a simulation of the National Airspace System. It includes models of flights, airports, airspaces, air traffic controls, traffic flow managements, and airline operation centers operating throughout the United States. It is used to predict system delays in response to future capacity and demand scenarios and perform benefits assessments of current and future airspace technologies and operational concepts. Facilitation of these studies requires that the simulation architecture supports plug and play of different air traffic control, traffic flow management, and airline operation center models and multi-fidelity modeling of flights, airports, and airspaces. The simulation is divided into two parts that are named, borrowing from classical control theory terminology, control and plant. The control consists of air traffic control, traffic flow management, and airline operation center models, and the plant consists of flight, airport, and airspace models. The plant can run open loop, in the absence of the control. However, undesired affects, such as conflicts and over congestions in the airspaces and airports, can occur. Different controls are applied, "plug and played", to the plant. A particular control is evaluated by analyzing how well it managed conflicts and congestions. Furthermore, the terminal area plants consist of models of airports and terminal airspaces. Each model consists of a set of nodes and links which are connected by the user to form a network. Nodes model runways, fixes, taxi intersections, gates, and/or other points of interest, and links model taxiways, departure paths, and arrival paths. Metering, flow distribution, and sequencing functions can be applied at nodes. Different fidelity model of how a flight transits are can be used by links. The fidelity of the model can be adjusted by the user by either changing the complexity of the node/link network-or the way that the link models how the flights transit from one node to the other.

Windhorst, Robert↗

Authentication Protocol for ICS without Encryption

Industrial control systems (ICS) describes the different types of control systems and associated instrumentation, which includes the devices, systems, networks, and controls used to operate and/or automate industrial processes. Depending on the industry, each ICS functions differently and are built to electronically manage tasks efficiently. Historically, industrial plants have included networking, and changes have been made over time without consideration for security. Many ICS devices currently have inadequate security features. Adding encryption in ICS is undesirable because many of the new detection method involves deep packet inspection. Therefore, we propose to develop a method of providing authentication between devices without encrypting the data. This testbed design will allow rapid development and experimentation with Cyber Physical Systems. Modbus connects to a supervisory computer with a Remote Terminal Unit (RTU) in Supervisory Control and Data Acquisition (SCADA) systems. Our project seeks to enhance security in the ICS environment by adding device authentication to the lower levels of the plant floor. This approach was selected as it allows verification of the sender of data from level 1/level 0 devices up to supervisory systems. This addition is desirable as it adds enhanced verification to the integrity of system data without obscuring (via encryption) the data being transmitted across the network. This is important because many ICS systems are real time critical (decryption adds to latency) and many other good security options require deep packet inspection (which encryption completely defeats). To accomplish this, a proof of concept of adding an authentication scheme was designed and will be tested with a common ICS protocol. The setup includes a single serial cable connecting the serial ports on two devices, a Master and a Minion. Authentication will occur in the following methods: - Data layer involving Transport Layer Security/Secure Socket Layer (TLS/SSL); - TCP/IP route implementation; - Key exchange algorithm; - Raspberry Pi programming for authentication. A portion of the setup includes adding Raspberry Pi into the system as the device is handling the authentication for the network traffic. Modbus is the chosen ICS protocol. It is used to transmit signals from instrumentation and control devices back to a main controller or data gathering system. Typically, the device identity is assumed but not verified. For verifying device identity, we constructed a digital signature architecture with TLS. A novel adaptation of this is to use only the digital signature that is provided and not the encryption capability. This allows ICS traffic to flow unobscured while still providing verification for the integrity of the data. Verification of this process will involve adding additional I/O level network devices with spoofed IP addresses and detecting the authenticated versus unauthenticated traffic. Our TLS implementation will have encryption disabled but should still allow for digital signing of data. This will enable our system to distinguish between real end devices and spoofed end devices. This experiment is presently building the test setup. It is focusing on a small setting to identify authenticated and unauthenticated devices on the network. This should allow for scaling to a larger network size after the initial proof of concept. One concern with this scheme is a certificate deployment has a large scale key management. Network size will thus increase the key management task for security. In addition, we believe that further research should be done on exploring the use of this method to defeat man-in-the-middle attacks.

42 ENGINEERING↗

Towards catchment classification in data-scarce regions

Assessing spatial variation in hydrologic processes can help to inform freshwater management and advance ecological understanding, yet many areas lack sufficient flow records on which to base classifications. Seeking to address this challenge, we apply concepts developed in data-rich settings to public, global data in order to demonstrate a broadly replicable approach to characterizing hydrologic variation. The proposed approach groups the basins associated with reaches in a river network according to key environmental drivers of hydrologic conditions. This initial study examines Colorado (USA), where long-term streamflow records permit comparison to previously distinguished flow regime types, and the Republic of Ecuador, where data limitations preclude such analysis. The flow regime types assigned to gages in Colorado corresponded reasonably well to the classes distinguished from environmental features. The divisions in Ecuador reflected major known biophysical gradients while also providing a higher resolution supplement to an existing depiction of freshwater ecoregions. Although freshwater policy and management decisions occur amidst uncertainty and imperfect knowledge, this classification framework offers a rigorous and transferrable means to distinguish catchments in data-scarce regions. The maps and attributes of the resulting ecohydrologic classes offer a departure point for additional study and data collection programs such as the placement of stations in under-monitored classes, and the divisions may serve as a preliminary template with which to structure conservation efforts such as environmental flow assessments.

13 HYDRO ENERGY↗

Characterizing the impact of particle behavior at fracture intersections in three-dimensional discrete fracture networks

We characterize the influence of different intersection mixing rules for particle tracking simulations on transport properties through three-dimensional discrete fracture networks. It is too computationally burdensome to explicitly resolve all fluid dynamics within a large three-dimensional fracture network. In discrete fracture network (DFN) models, mass transport at fracture intersections is modeled as a subgrid scale process based on a local Péclet number. The two most common mass transfer mixing rules are (1) complete mixing, where diffusion dominates mass transfer, and (2) streamline routing, where mass follows pathlines through an intersection. Although it is accepted that mixing rules impact local mass transfer through single intersections, the effect of the mixing rule on transport at the fracture network scale is still unresolved. Through the use of explicit particle tracking simulations, we study transport through a quasi-two-dimensional lattice network and a three-dimensional network whose fracture radii follow a truncated power-law distribution. We find that the impact of the mixing rule is a function of the initial particle injection condition, the heterogeneity of the velocity field, and the geometry of the network. Furthermore, our particle tracking simulations show that the mixing rule can particularly impact concentrations on secondary flow pathways. As a result, we relate these local differences in concentration to reactive transport and show that streamline routing increases the average mixing rate in DFN simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Fracture Sustainability in Enhanced Geothermal Systems: Experimental and Modeling Constraints

Enhanced geothermal systems (EGS) offer the potential for a much larger energy source than conventional hydrothermal systems. Hot, low-permeability rocks are prevalent at depth around the world, but the challenge of extracting thermal energy depends on the ability to create and sustain open fracture networks. Laboratory experiments were conducted using a suite of selected rock cores (granite, metasediment, rhyolite ash-flow tuff, and silicified rhyolitic tuff) at relevant pressures (uniaxial loading up to 20.7 MPa and fluid pressures up to 10.3 MPa) and temperatures (150–250 °C) to evaluate the potential impacts of circulating fluids through fractured rock by monitoring changes in fracture aperture, mineralogy, permeability, and fluid chemistry. Because a fluid in disequilibrium with the rocks (deionized water) was used for these experiments, there was net dissolution of the rock sample: this increased with increasing temperature and experiment duration. Thermal-hydrological-mechanical-chemical (THMC) modeling simulations were performed for the rhyolite ash-flow tuff experiment to test the ability to predict the observed changes. These simulations were performed in two steps: a thermal-hydrological-mechanical (THM) simulation to evaluate the effects of compression of the fracture, and a thermal-hydrological-chemical (THC) simulation to evaluate the effects of hydrothermal reactions on the fracture mineralogy, porosity, and permeability. Furthermore, these experiments and simulations point out how differences in rock mineralogy, fluid chemistry, and geomechanical properties influence how long asperity-propped fracture apertures may be sustained. Such core-scale experiments and simulations can be used to predict EGS reservoir behavior on the field scale.

15 GEOTHERMAL ENERGY↗

Master State Distributed Estimator (masde)

The MSE technology provides a distributed state estimation scheme for power utilities by making it a real time dynamic system, versus the static snapshot method used today. This will allow utilities to verify power grid readings and identify false data on the communication network. The prevailing method of estimation takes all the data from the system into a single load flow equation at the utilities command center. This becomes a large algorithm that is very time consuming to solve, providing the engineers only snap shots of the system. The MFD distributes an algorithm to devises that are already installed at various locations on the power grid. This allows for much faster computer times because the algorithms contain magnitudes less data. When each distributed estimator completes its calculation, it sends the results back to the master device to check against all other distributed estimators.

Reen, DylanW.↗

Software For Automated Leak Detection Using Infrared Camera

This code can read in videos or images in either a batch or real-time format. Videos are broken up into frames, and the frames are processed using an optical flow algorithm to decipher movement between adjacent frames. This adherent movement is ran through a convolutional neural network that automatically classifies the contents of the video. Additional content inside the code aids with noisy images and removal of nuisance movement.

Walker, CodyM. [Idaho National Laboratory (INL), I↗

LCLS-II Cryomodule and Cryogenic Distribution Control

LCLS-II is a superconducting upgrade to the existing Linear Coherent Light Source at the SLAC National Accelerator Laboratory. Construction is underway with a planned continuous wave beam rate of up to 1 MHz. Two cryogenic plants provide helium to a distribution system, and 37 cryomodules with superconducting cavities will operate with Liquid helium at 2.2K. The cryomodules and distribution system is controlled with networked PLC's and EPICS as an integrated system that work in concert for controlling valves, pressure, flow, and temperature. Interlocks and critical process information is communicated with the Low Level Radio Frequency, vacuum, and magnet systems. Engaging the controls community proved vital in advancing the controls architecture from a conventional design to a centralized, reliable, and cost-effective distributed platform.

43 PARTICLE ACCELERATORS↗

Cybersecurity Anomaly Detection in SCADA-Assisted OT Networks Using Ensemble-Based State Prediction Model

The cybersecurity threats of power system gradually grow due to the increased sophisticated interactions between Information Technology (IT) and Operational Technology (OT) networks. False data injection attack (FDIA) that aims to compromise the Supervisory Control and Data Acquisition (SCADA) measurement and disturb the system operation is one of such cyber threats. Such attacks can potentially lead to significant operational issues at the control centers and substations, and hence, result in severe physical consequences. To avoid catastrophic failure across the power grid resulting from these attacks, it is essential to arm the OT network with real-time vulnerability assessment tools. To this end, this paper outlines various drawbacks of the Purdue architecture model to defend against cyberattacks in the OT network. Furthermore, a novel ensemble-based state prediction model is proposed to detect cybersecurity anomalies in SCADA assisted OT networks. The proposed model uses control center level generation and load forecasts, scheduled, and forced outages, power flow solutions, and the substation level historical data. The hypothesis of the proposed scheme relies on the fact that additional control center and substation data can hardly be accessed and compromised by attackers. One of the vital features of the proposed scheme is an hour-ahead prediction of the operational feasibility of the SCADA measurement range at the control center and substation in real time helps in detecting anomalies in measurements across both substation and the control center.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Geology of selected lava tubes in the Bend Area, Oregon

Longitudinal profiles representing 5872.5 m of mapped lava tubes and a photogeologic map relating lava tubes to surface geology, regional structure and topography are presented. Three sets of lava tubes were examined: (1) Arnold Lava Tube System (7km long) composed of collapsed and uncollapsed tube segments and lava ponds, (2) Horse Lava Tube System (11 km long) composed of parallel and anastomosing lava tube segments, and (3) miscellaneous lava tubes. Results of this study tend to confirm the layered lava hypothesis of Ollier and Brown (1965) for lava tube formation; however, there are probably several modes of formation for lava tubes in general. Arnold System is a single series of tubes apparently formed in a single basalt flow on a relatively steep gradient. The advancing flow in which the tubes formed was apparently temporarily halted, resulting in the formation of lava ponds which were inflated and later drained by the lava tube system. Horse System probably formed in multiple, interconnected flows. Pre-flow gradient appears to have been less than for Arnold System, and resulted in meandrous, multiple tube networks.

Greely, R.↗

Advanced panel-type influence coefficient methods applied to subsonic flows

An advanced technique for solving the linear integral equations of three-dimensional subsonic potential flows (steady, inviscid, irrotational and incompressible) about arbitrary configurations is presented. It involves assembling select, logically consistent networks whose construction comprises four tasks, which are described in detail: surface geometry definition; singularity strength definition; control point and boundary condition specification; and calculation of induced potential or velocity. The technique is applied to seven wing examples approached by four network types: source/analysis, doublet/analysis, source/design, and doublet/design. The results demonstrate the forgiveness of the model to irregular paneling and the practicality of combined analysis/design boundary conditions. The appearance of doublet strength mismatch is a valuable indicator of locally inadequate paneling.

Johnson, F. T.↗

Evolution of the global water cycle on Mars: The geological evidence

The geological evidence for active water cycling early in the history of Mars (Noachian geological system or heavy bombardment) consists almost exclusively of fluvial valley networks in the heavily cratered uplands of the planet. It is commonly assumed that these landforms required explanation by atmospheric processes operating above the freezing point of water and at high pressure to allow rainfall and liquid surface runoff. However, it has also been documented that nearly all valley networks probably formed by subsurface outflow and sapping erosion involving groundwater outflow prior to surface-water flow. The prolonged ground-water flow also requires extensive water cycling to maintain hydraulic gradients, but is this done via rainfall recharge, as in terrestrial environments?

Baker, V. R.↗

The sensitivity of various mode sets for probing differential rotation shear zones

The potential of future data sets from global oscillations network group (GONG) and solar oscillations investigation (SOI) for resolving long-lived azimuthal jets and shearing flows, is investigated. Various artificial data sets are constructed, containing noise resembling that of a one-year observation run. These are inverted using a two dimensional regularized least squares inversion. The ability of this method to form well localized averages of the rotation rate, as measured by the averaging kernels, is investigated using an extensive mode set and subsets. It is shown that it is possible to keep the noise in the solution down to a few nHz in much of the solar interior, while obtaining a reasonable resolution for a GONG-like data set. At low latitudes in the middle of the convection zone, an angular resolution of less than 10 deg and a radial resolution of about 0.04 solar radii, are obtained. The averaging kernels depend on the mode set, and a reduction in the number of modes tends to introduce small-scale near surface structures into the averaging kernels which would adversely affect the inferred rotation rate.

Christensen-Dalsgaard, Joergen↗