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,603 records · Page 89

Quantifying local and global mass balance errors in physics-informed neural networks

Physics-informed neural networks (PINN) have recently become attractive for solving partial differential equations (PDEs) that describe physics laws. By including PDE-based loss functions, physics laws such as mass balance are enforced softly in PINN. This paper investigates how mass balance constraints are satisfied when PINN is used to solve the resulting PDEs. We investigate PINN’s ability to solve the 1D saturated groundwater flow equations (diffusion equations) for homogeneous and heterogeneous media and evaluate the local and global mass balance errors. We compare the obtained PINN’s solution and associated mass balance errors against a two-point finite volume numerical method and the corresponding analytical solution. We also evaluate the accuracy of PINN in solving the 1D saturated groundwater flow equation with and without incorporating hydraulic heads as training data. We demonstrate that PINN’s local and global mass balance errors are significant compared to the finite volume approach. Tuning the PINN’s hyperparameters, such as the number of collocation points, training data, hidden layers, nodes, epochs, and learning rate, did not improve the solution accuracy or the mass balance errors compared to the finite volume solution. Mass balance errors could considerably challenge the utility of PINN in applications where ensuring compliance with physical and mathematical properties is crucial.

54 ENVIRONMENTAL SCIENCES↗

Hydro-mechanical model for wetting/drying and fracture development in geomaterials

This study presents a modeling approach for studying hydro-mechanical coupled processes, including fracture development, within geological formations. This is accomplished through the novel linking of two codes: TOUGH2, which is a widely used simulator of subsurface multiphase flow based on the finite volume method; and an implementation of the Rigid-Body-Spring Network (RBSN) method, which provides a discrete (lattice) representation of material elasticity and fracture development. The modeling approach is facilitated by a Voronoi-based discretization technique, capable of representing discrete fracture networks. The TOUGH–RBSN simulator is intended to predict fracture evolution, as well as mass transport through permeable media, under dynamically changing hydrologic and mechanical conditions. Numerical results are compared with those of two independent studies involving hydro-mechanical coupling: (1) numerical modeling of swelling stress development in bentonite; and (2) experimental study of desiccation cracking in a mining waste. The comparisons show good agreement with respect to moisture content, stress development with changes in pore pressure, and time to crack initiation. Finally, the observed relationship between material thickness and crack patterns (e.g., mean spacing of cracks) is captured by the proposed modeling approach.

58 GEOSCIENCES↗

A data-driven method for modelling dissipation rates in stratified turbulence

We present a deep probabilistic convolutional neural network (PCNN) model for predicting local values of small-scale mixing properties in stratified turbulent flows, namely the dissipation rates of turbulent kinetic energy and density variance, $\varepsilon$ and $\chi$ . Inputs to the PCNN are vertical columns of velocity and density gradients, motivated by data typically available from microstructure profilers in the ocean. The architecture is designed to enable the model to capture several characteristic features of stratified turbulence, in particular the dependence of small-scale isotropy on the buoyancy Reynolds number $Re_b:=\varepsilon /(\nu N^2)$ , where $\nu$ is the kinematic viscosity and $N$ is the background buoyancy frequency, the correlation between suitably locally averaged density gradients and turbulence intensity and the importance of capturing the tails of the probability distribution functions of values of dissipation. Empirically modified versions of commonly used isotropic models for $\varepsilon$ and $\chi$ that depend only on vertical derivatives of density and velocity are proposed based on the asymptotic regimes $Re_b\ll 1$ and $Re_b\gg 1$ , and serve as an instructive benchmark for comparison with the data-driven approach. When trained and tested on a simulation of stratified decaying turbulence which accesses a range of turbulent regimes (associated with differing values of $Re_b$ ), the PCNN outperforms assumptions of isotropy significantly as $Re_b$ decreases, and additionally demonstrates improvements over the fitted empirical models. A differential sensitivity analysis of the PCNN facilitates a comparison with the theoretical models and provides a physical interpretation of the features enabling it to make improved predictions.

42 ENGINEERING↗

Spatio–Temporal Machine Learning for Regional to Continental Scale Terrestrial Hydrology

Integrated hydrologic models can simulate coupled surface and subsurface processes but are computationally expensive to run at high resolutions over large domains. Here we develop a novel deep learning model to emulate subsurface flows simulated by the integrated ParFlow–CLM model across the contiguous US. We compare convolutional neural networks like ResNet and UNet run autoregressively against our novel architecture called the Forced SpatioTemporal RNN (FSTR). The FSTR model incorporates separate encoding of initial conditions, static parameters, and meteorological forcings, which are fused in a recurrent loop to produce spatiotemporal predictions of groundwater. We evaluate the model architectures on their ability to reproduce 4D pressure heads, water table depths, and surface soil moisture over the contiguous US at 1 km resolution and daily time steps over the course of a full water year. The FSTR model shows superior performance to the baseline models, producing stable simulations that capture both seasonal and event–scale dynamics across a wide array of hydroclimatic regimes. The emulators provide over 1,000× speedup compared to the original physical model, which will enable new capabilities like uncertainty quantification and data assimilation for integrated hydrologic modeling that were not previously possible. Our results demonstrate the promise of using specialized deep learning architectures like FSTR for emulating complex process–based models without sacrificing fidelity.

54 ENVIRONMENTAL SCIENCES↗

Fracture Characterization Via AI‐Assisted Analysis of Temperature Logs

Abstract Fractures control fluid flow, mass transport, and heat transfer in a geothermal reservoir. This makes accurate characterization of fracture networks a prerequisite for optimal design and control of a reservoir's exploitation. We develop a deep‐learning procedure to identify fracture locations via interpretation of temporally and spatially continuous downhole temperature measurements. A long short‐term memory fully convolutional network (LSTM‐FCN) is used both to capture long‐term dependencies in sequential temperature data and to distill local features around fractures. A wellbore and fractured‐reservoir thermal model is established to generate temperature data for network training. The trained LSTM‐FCN exhibits a unique ability to detect multiple fractures intersecting a borehole. We use the LSTM‐FCN algorithm to evaluate the effectiveness of different‐stage wellbore temperature measurements on fracture detection in a complex fractured system. Our experiments reveal that the use of various‐stage temperature information as an input feature set improves the robustness of fracture detection to noise interference. This study indicates the practical feasibility of obtaining accurate fracture‐network reconstructions from temperature signals, at reasonable computational cost.

Yang, Xiaoyu↗

Linear shaped-charge jet optimization using machine learning methods

Linear shaped charges are used to focus energy into rapidly creating a deep linear incision. The general design of a shaped charge involves detonating a confined mass of high explosive (HE) with a metal-lined concave cavity on one side to produce a high velocity jet for the purpose of striking and penetrating a given material target. This jetting effect occurs due to the interaction of the detonation wave with the cavity geometry, which produces an unstable fluid phenomenon known as the Richtmyer–Meshkov instability and results in the rapid growth of a long narrow jet. We apply machine learning and optimization methods to hydrodynamics simulations of linear shaped charges to improve the simulated jet characteristics. The designs that we propose and investigate in this work generally involve modifying the behavior of the detonation waves prior to interaction with the liner material. These designs include the placement of multiple detonators and the use of metal inclusions within the HE. In conclusion, we are able to produce a linear shaped-charge design with a higher penetration depth than the baseline case that we consider and accomplish this using the same amount of or less HE.

36 MATERIALS SCIENCE↗

Neural posterior unfolding

Differential cross section measurements are the currency of scientific exchange in particle and nuclear physics. A key challenge for these analyses is the correction for detector distortions, known as deconvolution or unfolding. Binned unfolding of cross section measurements traditionally rely on the regularized inversion of the response matrix that represents the detector response, mapping pre-detector (`particle level') observables to post-detector (`detector level') observables. In this paper we introduce Neural Posterior Unfolding, a modern, Bayesian approach that leverages normalizing flows for unfolding. By using normalizing flows for neural posterior estimation, NPU offers several key advantages including implicit regularization through the neural network architecture, fast amortized inference that eliminates the need for repeated retraining, and direct access to the full uncertainty in the unfolded result. In addition to introducing NPU, we implement a classical Bayesian unfolding method called Fully Bayesian Unfolding (FBU) in modern Python so it can also be studied. These tools are validated on simple Gaussian examples and then tested on simulated jet substructure examples from the Large Hadron Collider (LHC). We find that the Bayesian methods are effective and worth additional development to be analysis ready for cross section measurements at the LHC and beyond.

Analysis and statistical methods↗

Generic Crystalline Disposal Reference Case

A generic reference case for disposal of spent nuclear fuel and high-level radioactive waste in crystalline rock is outlined. The generic cases are intended to support development of disposal system modeling capability by establishing relevant baseline conditions and parameters. Establishment of a generic reference case requires that the emplacement concept, waste inventory, waste form, waste package, backfill/buffer properties, EBS failure scenarios, host rock properties, and biosphere be specified. The focus in this report is on those elements that are unique to crystalline disposal, especially the geosphere representation. Three emplacement concepts are suggested for further analyses: a waste packages containing 4 PWR assemblies emplaced in boreholes in the floors of tunnels (KBS-3 concept), a 12-assembly waste package emplaced in tunnels, and a 32-assembly dual purpose canister emplaced in tunnels. In addition, three failure scenarios were suggested for future use: a nominal scenario involving corrosion of the waste package in the tunnel emplacement concepts, a manufacturing defect scenario applicable to the KBS-3 concept, and a disruptive glaciation scenario applicable to both emplacement concepts. The computational approaches required to analyze EBS failure and transport processes in a crystalline rock repository are similar to those of argillite/shale, with the most significant difference being that the EBS in a crystalline rock repository will likely experience highly heterogeneous flow rates, which should be represented in the model. The computational approaches required to analyze radionuclide transport in the natural system are very different because of the highly channelized nature of fracture flow. Computational workflows tailored to crystalline rock based on discrete transport pathways extracted from discrete fracture network models are recommended.

Environmental Protection↗

Pressure Amplification Off High Impedance Barriers in DDT

The Deflagration-to-Detonation Transition (DDT) in one-dimensional porous explosive, where combustion in an explosive transitions to detonation, can be described by the following model. This simplified model proceeds in five steps, as follows: 1) Ignition of the explosive, surface burning. 2) Convective burning, with the flame front penetrating through the porous network of the explosive. This proceeds until the pressure grows high enough to result in choked flow in the pores restricting the convective burn. 3) The choked flow results in the formation of a high-density compact of explosive. This compact is driven into undisturbed material by the pressure of the burning explosive. See Figure1. 4) The compression of the undisturbed porous explosive by the compact leads to the ignition of a compressive burn. This builds in pressure until a supported shock forms. 5) The shock builds in pressure until detonation occurs. See Figure 2 for an overview streak of the proceeding steps.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

The Deflagration-to-Detonation Transition in Two Dimensions

The Deflagration-to-Detonation Transition (DDT) in one-dimensional porous explosive, where combustion in an explosive transitions to detonation, can be described by the following model. This simplified model proceeds in five steps, as follows: 1) Ignition of the explosive, surface burning. 2) Convective burning, with the flame front penetrating through the porous network of the explosive. This proceeds until the pressure grows high enough to result in choked flow in the pores restricting the convective burn. 3) The choked flow results in the formation of a high-density compact of explosive. This compact is driven into undisturbed material by the pressure of the burning explosive. 4) The compression of the undisturbed porous explosive by the compact leads to the ignition of a compressive burn. This builds in pressure until a supported shock forms. 5) The shock builds in pressure until detonation occurs. See Figures 1 and 2 for an overview of the proceeding steps. It has been assumed in the past that the same mechanism which drives DDT in the one dimensional case will apply to both two and three dimensions. However, this has never been experimentally verified and it is possible that an entirely new mechanism drives DDT higher dimensions. In order to test this a series of two dimensional DDT tests were performed to help provide empirical evidence of the mechanism in higher dimensions.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Fractionation of Lignin for Selective Shape Memory Effects at Elevated Temperatures

We report a facile approach to control the shape memory effects and thermomechanical characteristics of a lignin-based multiphase polymer. Solvent fractionation of a syringylpropane-rich technical organosolv lignin resulted in selective lignin structures having excellent thermal stability coupled with high stiffness and melt-flow resistance. The fractionated lignins were reacted with rubber in melt-phase to form partially networked elastomer enabling selective programmability of the material shape either at 70 °C, a temperature that is high enough for rubbery matrix materials, or at an extremely high temperature, 150 °C. Utilizing appropriate functionalities in fractionated lignins, tunable shape fixity with high strain and stress recovery, particularly high-stress tolerance were maintained. Detailed studies of lignin structures and chemistries were correlated to molecular rigidity, morphology, and stress relaxation, as well as shape memory effects of the materials. The fractionation of lignin enabled enrichment of specific lignin properties for efficient shape memory effects that broaden the materials’ application window. Electron microscopy, melt-rheology, dynamic mechanical analysis and ultra-small angle neutron scattering were conducted to establish morphology of acrylonitrile butadiene rubber (NBR)-lignin elastomers from solvent fractionated lignins.

36 MATERIALS SCIENCE↗

Gauge Measures Hot Erosion Of Insulation

Compact ladder networks containing thermocouples and resistors serve as gauges to measure erosion of insulating materials by hot, flowing gases. Designed for use in blast tubes during tests of rocket motors. Concept readily specialized to other applications in which insulation subject to erosion.

Gould, Reginald J.↗

MEANS User Documentation

MEANS is an acronym for the MIT Extensible Air Network Simulation. MIT is the Massachusetts Institute of Technology, where MEANS was developed. MEANS models the flow of aircraft in a simulation of the US National Airspace System. It can optionally follow the flow of passengers as well. The flight model takes the form of several queues through which the aircraft pass, as well as delays at certain times; the passengers can be transported via flights. MEANS is an event-based model. This means that events are scheduled at some point in the future, and once an event has executed, the simulation t h e advances to that of the next event. This is a different approach than some simulations (particularly those of physical systems), where time is advanced in fixed increments and the state of the system is updated at each step.This documentation is a "user guide" for MEANS - it describes how to USE means, not write new modules. For programming instructions, see the programming documentation, available as a separate document.This document assumse no programming knowledge. It will describe how to run MEANS, understand the data formats, and use such utility programs as already exist. For more detailed information on the internal workings of MEANS, and how to programmatically add to it, refer to the programming documentation.

Melconian, Terran↗

Air Traffic Management Technology Demonstration - 3 (ATD-3) Dynamic Weather Routes Domestic En Route Concept of Operations Synopsis Version 1.0

NASA conducted a survey in 2015 of stakeholders in industry, flight operators, and others to collect data on the most pressing air traffic management (ATM) issues impacting en route operations. Weather related delays and interruptions were identified as the most significant and costly problem impacting en route operations in today's system. Additionally, future traffic demands will place greater demands on the FAA's air traffic control system, necessitate a more efficient and integrated traffic management system and provide controller tools to leverage human productivity. Consequently, NASA is assembling a suite of integrated ground and aircraft-based technologies and decision-making aids under the Airspace Technology Demonstrations 3 (ATD-3) project to continuously search for more efficient weather-avoidance routes and to rapidly respond to weather changes and associated traffic management initiatives. These route efficiency enhancements are targeted at en route airspace and will enable continued use of en route and arrival metering in the presence of weather. NASA's Dynamic Weather Routes (DWR) concept represents a subset of those technologies and focuses primarily on the cruise phase of flight. It leverages existing air traffic weather, airspace, and traffic data, as well as improvements in navigation, surveillance, communication and digital network technologies, and builds upon existing ATM automation to address the shortcomings associated with strategic traffic flow management initiatives and weather forecasting uncertainties. The potential benefits in the form of time, fuel, and cost savings are significant. This concept of operations synopsis provides a detailed description of DWR, its potential benefits and notional implementation paths.

DWR↗

Basaltic Fissure Types on Earth: Suitable Analogs to Evaluate the Origins of Volcanic Terrains on the Moon and Mars?

Basaltic eruptive fissures of the Great Rift and surroundings on the eastern Snake River Plain of Idaho, USA, and selected volcanic features in Hawai`i, Iceland and northern Africa were surveyed for their relevancy as planetary analogs. Evaluated during field investigations and in satellite imagery for structures, physiography, and geologic setting, fissures were categorized into four broad types: (1) simple, monogenetic fissures with obvious volcanic constructs or deposits, (2) monogenetic fissures now obscured by low shields or relatively large cones, (3) polygenetic volcanic rift zones with multiple vents and deposits, and (4) compound regional fissure systems or dike swarms that comprise major rift zones or large volcanic terrains. Using this classification as an initial base, we surveyed imagery of volcanic features for likely fissure vents in two major geologic settings on the Moon: floor-fractured craters (FFCs) and mare and cryptomare provinces. Two major regions on Mars, the volcanic plains around Alba Mons and the greater Tharsis region, were also surveyed for fissure types and volcanic associations of fissure-like features. The planetary surveys suggest that the proposed classification provides a suitable analog starting point to interpret structures associated with fissure systems on the Moon and Mars. With few exceptions, our survey indicates that each of the studied terrains exhibits a dominant fissure type. Type 1 fissures, most with pyroclastic deposits, prevail in lunar FFCs and mare-like regions; whereas type 2 fissures are ubiquitous in the Tharsis region of Mars and a few exist on the Moon as low shields. Type 3 volcanic rift zones are not common on either the Moon or Mars, although they might become evident in future work on chemically evolved terrains. Type 4 fissures are inferred in mare terrains, often represented as the extensions of major linear rille networks or rimae, with possibly complex dike swarms that were buried beneath voluminous mare basalt lava flows. Likewise, numerous flood lavas on Mars are possibly associated with now-obscured or difficult to define type 4 fissure systems.

Scott S. Hughes↗

On Generalizably Improving RANS Predictions of Flow Separation and Reattachment

This work presents progress in an ongoing data-driven turbulence modeling effort at NASA Langley Research Center to improve predictions for separated flows. Two augmentation strategies are presented – one to improve predictions for the point of reattachment after smooth-body separation, and another to improve predictions for flow separation due to prolonged adverse pressure gradients. The reattachment augmentation was inspired from lessons learned during attempts at using field inversion and machine learning (FIML); however, it is prescribed as a conventional analytic closed-form expression. Results show that, when introduced into the SST k-w model, it shortens the otherwise overpredicted separation bubble length for three cases characterized by different Reynolds numbers and flow geometries. The second augmentation, which addresses flow separation under prolonged adverse pressure gradients, is inferred via FIML and applied as a neural network within the Spalart-Allmaras model. It uses features that are chosen such that the presented results show an improved predictive accuracy for the training case without affecting other cases where the augmentation is not needed. Further testing is required to establish the generalizability of this augmentation. The rationale behind the chosen features that serve as inputs to both these augmentations is also presented for the benefit of the reader.

machine learning↗

Pore-network model of evaporation-induced salt precipitation in porous media: The effect of correlations and heterogeneity

Salt transport and precipitation in porous media constitute a set of complex and fascinating phenomena that are of considerable interest to several important problems, ranging from storage of CO 2 in geological formations, to soil fertility, and protection of pavements and roads, as well as historical monuments. The phenomena occur at the pore scale and are greatly influenced by the heterogeneity of the pore space morphology. In this paper, we present a pore-network (PN) model to study the phenomena. Vapor diffusion, capillary effect at the brine-vapor interface, flow of brine, and transport of salt and its precipitation in the pores that plug the pores partially or completely are all accounted for. The drying process is modeled by the invasion percolation, while transport of salt in brine is accounted for by the convective-diffusion equation. We demonstrate that the drying patterns, the clustering and connectivity of the pore throats in which salt precipitation occurs, the saturation distribution, and the drying rate are all strongly dependent upon the pore-size distribution, the correlations among the pore sizes, and the anisotropy of the pore space caused by stratification that most natural porous media contain. In particular, if the strata are more or less parallel to the direction of injection of the gas that dries out the pore space (air, for example) and/or causes salt precipitation (CO 2 , for example), the drying rate increases significantly. Moreover, salt tends to precipitate in clusters of neighboring pores that are parallel to the open surface of the porous medium.

42 ENGINEERING↗