Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluid Network”

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 199 records · Page 11

A New Physics-Based Modeling of Multiple Non-Planar Hydraulic Fractures Propagation

Because of the low permeability in shale plays, closely spaced hydraulic fractures and multilateral horizontal wells are generally required to improve production. Therefore, understanding the potential fracture interaction and stress evolution is critical in optimizing fracture/well design and completion strategy in multi-stage horizontal wells. In this paper, a novel fully coupled reservoir flow and geomechanics model based on the dual-lattice system is developed to simulate multiple non-planar fractures propagation. The numerical model from Discrete Element Method (DEM) is used to simulate the mechanics of fracture propagations and interactions, while a conjugate irregular lattice network is generated to represent fluid flow in both fractures and formation. The fluid flow in the formation is controlled by Darcy’s law, but within fractures it is simulated by using cubic law for laminar flow through parallel plates. Initiation, growth and coalescence of the microcracks will lead to the generation of macroscopic fractures, which is explicitly mimicked by failure and removal of bonds between particles from the discrete element network. We investigate the fracture propagation path in both homogeneous and heterogeneous reservoirs using the simulator developed. Stress shadow caused by the transverse fracture will change the orientation of principal stress in the fracture neighborhood, which may inhibit or alter the growth direction of nearby fracture clusters. However, the initial in-situ stress anisotropy often helps overcome this phenomenon. Under large in-situ stress anisotropy, the hydraulic fractures are more likely to propagate in a direction that is perpendicular to the minimum horizontal stress. Under small in-situ stress anisotropy, there is a greater chance for fractures from nearby clusters to merge with each other. Then, we examine the differences in fracture geometry caused by fracturing in cemented or uncemented wellbore. Moreover, the impact of intrinsic reservoir heterogeneity caused by the rock fabric and mineralogy on fracture nucleation and propagation paths is examined through a three-layered reservoir. Finally, we apply the method to a realistic heterogeneous dataset.

02 PETROLEUM↗

Fouling resistance prediction using artificial neural network nonlinear auto-regressive with exogenous input model based on operating conditions and fluid properties correlations

Fouling in a heat exchanger in Crude Preheat Train (CPT) refinery is an unsolved problem that reduces the plant efficiency, increases fuel consumption and CO{sub 2} emission. The fouling resistance behavior is very complex. It is difficult to develop a model using first principle equation to predict the fouling resistance due to different operating conditions and different crude blends. In this paper, Artificial Neural Networks (ANN) MultiLayer Perceptron (MLP) with input structure using Nonlinear Auto-Regressive with eXogenous (NARX) is utilized to build the fouling resistance model in shell and tube heat exchanger (STHX). The input data of the model are flow rates and temperatures of the streams of the heat exchanger, physical properties of product and crude blend data. This model serves as a predicting tool to optimize operating conditions and preventive maintenance of STHX. The results show that the model can capture the complexity of fouling characteristics in heat exchanger due to thermodynamic conditions and variations in crude oil properties (blends). It was found that the Root Mean Square Error (RMSE) are suitable to capture the nonlinearity and complexity of the STHX fouling resistance during phases of training and validation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

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↗

Liquid crystalline cellulose-based nematogels

Physical properties of composite materials can be pre-engineered by controlling their structure and composition at the mesoscale. However, approaches to achieving this are limited and rarely scalable. We introduce a new breed of self-assembled nematogels formed by an orientationally ordered network of thin cellulose nanofibers infiltrated with a thermotropic nematic fluid. The interplay between orientational ordering within the nematic network and that of the small-molecule liquid crystal around it yields a composite with highly tunable optical properties. By means of combining experimental characterization and modeling, we demonstrate submillisecond electric switching of transparency and facile responses of the composite to temperature changes. Finally, we discuss a host of potential technological uses of these self-assembled nematogel composites, ranging from smart and privacy windows to novel flexible displays.

36 MATERIALS SCIENCE↗

Prospects of federated machine learning in fluid dynamics

Physics-based models have been mainstream in fluid dynamics for developing predictive models. In recent years, machine learning has offered a renaissance to the fluid community due to the rapid developments in data science, processing units, neural network based technologies, and sensor adaptations. So far in many applications in fluid dynamics, machine learning approaches have been mostly focused on a standard process that requires centralizing the training data on a designated machine or in a data center. In this article, we present a federated machine learning approach that enables localized clients to collaboratively learn an aggregated and shared predictive model while keeping all the training data on each edge device. We demonstrate the feasibility and prospects of such a decentralized learning approach with an effort to forge a deep learning surrogate model for reconstructing spatiotemporal fields. Our results indicate that federated machine learning might be a viable tool for designing highly accurate predictive decentralized digital twins relevant to fluid dynamics.

36 MATERIALS SCIENCE↗

Agent-Based Chemical Plume Tracing Using Fluid Dynamics

This paper presents a rigorous evaluation of a novel, distributed chemical plume tracing algorithm. The algorithm is a combination of the best aspects of the two most popular predecessors for this task. Furthermore, it is based on solid, formal principles from the field of fluid mechanics. The algorithm is applied by a network of mobile sensing agents (e.g., robots or micro-air vehicles) that sense the ambient fluid velocity and chemical concentration, and calculate derivatives. The algorithm drives the robotic network to the source of the toxic plume, where measures can be taken to disable the source emitter. This work is part of a much larger effort in research and development of a physics-based approach to developing networks of mobile sensing agents for monitoring, tracking, reporting and responding to hazardous conditions.

Zarzhitsky, Dimitri↗

Water Hammer Simulations of Monomethylhydrazine Propellant

Fluid Transient analysis is important for the design of spacecraft propulsion system to ensure structural stability of the system in the event of sudden closing or opening of the valve. Generalized Fluid System Simulation Program (GFSSP), a general purpose flow network code developed at NASA/MSFC is capable of simulating pressure surge due to sudden opening or closing of valve when thermodynamic properties of real fluid are available for the entire range of simulation. Specifically GFSSP needs an accurate representation of pressure density relationship in order to predict pressure surge during a fluid transient. Unfortunately, the available thermodynamic property programs such as REFPROP, GASP or GASPAK do not provide the thermodynamic properties of Monomethylhydrazine(MMH). This work illustrates the process used for building a customized table of properties of state variables from available properties and speed of sound that is required by GFSSP for simulation. Good agreement was found between the simulations and measured data. This method can be adopted for modeling flow networks and systems with other fluids whose properties are not known in detail in order to obtain general technical insight.

Burkhardt, Zachary↗

Water Hammer Simulations of MMH Propellant - New Capability Demonstration of the Generalized Fluid Flow Simulation Program

Fluid Transient analysis is important for the design of spacecraft propulsion system to ensure structural stability of the system in the event of sudden closing or opening of the valve. Generalized Fluid System Simulation Program (GFSSP), a general purpose flow network code developed at NASA/MSFC is capable of simulating pressure surge due to sudden opening or closing of valve when thermodynamic properties of real fluid are available for the entire range of simulation. Specifically GFSSP needs an accurate representation of pressure-density relationship in order to predict pressure surge during a fluid transient. Unfortunately, the available thermodynamic property programs such as REFPROP, GASP or GASPAK does not provide the thermodynamic properties of Monomethylhydrazine (MMH). This paper will illustrate the process used for building a customized table of properties of state variables from available properties and speed of sound that is required by GFSSP for simulation. Good agreement was found between the simulations and measured data. This method can be adopted for modeling flow networks and systems with other fluids whose properties are not known in detail in order to obtain general technical insight. Rigorous code validation of this approach will be done and reported at a future date.

Burkhardt, Z.↗

Microfluidic Sensors with Impregnated Fluorophores for Simultaneous Imaging of Spatial Structure and Chemical Oxygen Gradients

Inside surfaces of polystyrene microfluidic structures were impregnated with the oxygen sensing dye (PtTFPP) using a solvent-induced fluorophore impregnation (SIFI) method. By this technique, microfluidic oxygen sensors are obtained that enable simultaneous imaging of both chemical oxygen gradients and the physical structure of the microfluidic interior. A gentle method of fluorophore impregnation using acetonitrile solutions of PtTFPP at 50oC was developed leading to 10 micron deep region containing fluorophore. This region is localized at the surface to sense oxygen in the interior fluid during use. Regions of the device that do not contact the interior fluid pathways lack fluorophores and are dark in fluorescent imaging. The technique was demonstrated on straight microchannel and pore network devices, the latter having pillars of 300 µm diameter spaced center to center at 340 µm providing pore throats of 40 µm. Sensing within channels or pores, and imaging across the pore network devices were performed using a Lambert LIFA-P frequency domain fluorescence lifetime imaging system on a Leical microscope platform. Calibrations of different devices prepared by the SIFI method were indistinguishable. Gradient imaging showed fluorescent regions corresponding to the fluid pore network, dark pillars, and fluorescent lifetime varying across the gradient, thus providing both physical and chemical imaging. More generally, the SIFI technique can impregnate the interior surfaces of other polystyrene containers, such as cuvettes or cell and tissue culture containers, to enable sensing interior conditions.

microfluidic, sensor, imaging, oxygen↗

An analytical thermohydraulic model for discretely fractured geothermal reservoirs

In discretely fractured reservoirs such as those found in Enhanced/Engineered Geothermal Systems (EGS), knowledge of the fracture network is important in understanding the thermal hydraulics, i.e., how the fluid flows and the resulting temporal evolution of the subsurface temperature. Here, the purpose of this study was to develop an analytical model of the fluid flow and heat transport in a discretely fractured network that can be used for a wide range of modeling applications and serve as an alternative analysis tool to more computationally intensive numerical codes. Given the connectivity and structure of a fracture network, the flow in the system was solved using a linear system of algebraic equations for the pressure at the nodes of the network. With the flow determined, the temperature in the fracture was solved by coupling convective heat transport in the fracture with one-dimensional heat conduction perpendicular to the fracture, employing the Green’s function derived solution for a single discrete fracture. The predicted temperatures along the fracture surfaces from the analytical solution were compared to numerical simulations using the TOUGH2 reservoir code. Through two case studies, we showed the capabilities of the analytical model and explored the effect of uncertainty in the fracture apertures and network structure on thermal performance. While both sources of uncertainty independently produce large variations in production temperature, uncertainty in the network structure, whenever present, had a predominant influence on thermal performance.

58 GEOSCIENCES↗

Production and Characterization of Additively Manufactured Radiator Panels with Integral Branching Heat Pipes for High-Temperature Heat Rejection

Emerging concepts for fission surface power and nuclear electric propulsion necessitate lightweight, mechanically robust, and thermally efficient heat rejection radiators. State-of-the-art intermediate-temperature (~400 K) composite radiator assemblies have been developed based on titanium-water heat pipes bonded to metal, graphite, and carbon-fiber-based panels. NASA has identified a need for new radiator concepts that can operate at even higher temperatures (500 – 600 K), minimize thermal resistances and thermal stress failures at bond interfaces, and approach areal densities of 2 – 3 kg m -2 . To meet these needs, our team is developing additively manufactured (AM) radiator panels with integral branching wicking heat pipe networks. Water is selected as the working fluid for this temperature range. Based on simulations and thermal vacuum experiments, these branching embedded heat pipe networks can efficiently distribute heat over panels for finned surface efficiencies of η f >70% at T H = 500 K input heat. This paper first presents laser powder-bed fusion AM strategies to produce embedded porous structures for wicking heat pipes in Inconel 718 and titanium alloys (commercially pure and Ti-6Al-4V alloys). Post-build chemical and thermal treatments are described that yield hydrophilic wicking surfaces for operation with water. Transient rate-of-rise experiments with water and acetone are reported that yield estimates for AM wick porosity (ϵ), permeability (K), and effective pore radius (r pore ). Based on the wick characterization results, small prototype radiator panels (75 × 125 mm) with integrated heat pipe networks were manufactured. Heat rejection performance data are presented from cold thermal vacuum testing, with heat input temperatures up to ~510 K. Future efforts will focus on improving heat pipe performance, optimizing radiator mass, and evaluating larger panels to assess scalability.

thermal management↗

Production and Characterization of Additively Manufactured Radiator Panels With Integral Branching Heat Pipes for High-Temperature Heat Rejection

Emerging concepts for fission surface power and nuclear electric propulsion necessitate lightweight, mechanically robust, and thermally efficient heat rejection radiators. State-of-the-art intermediate-temperature (~400 K) composite radiator assemblies have been developed based on titanium-water heat pipes bonded to metal, graphite, and carbon-fiber-based panels. NASA has identified a need for new radiator concepts that can operate at even higher temperatures (500 – 600 K), minimize thermal resistances and thermal stress failures at bond interfaces, and approach areal densities of 2 – 3 kg m -2 . To meet these needs, our team is developing additively manufactured (AM) radiator panels with integral branching wicking heat pipe networks. Water is selected as the working fluid for this temperature range. Based on simulations and thermal vacuum experiments, these branching embedded heat pipe networks can efficiently distribute heat over panels for finned surface efficiencies of η f >70% at T H = 500 K input heat. This paper first presents laser powder-bed fusion AM strategies to produce embedded porous structures for wicking heat pipes in Inconel 718 and titanium alloys (commercially pure and Ti-6Al-4V alloys). Post-build chemical and thermal treatments are described that yield hydrophilic wicking surfaces for operation with water. Transient rate-of-rise experiments with water and acetone are reported that yield estimates for AM wick porosity (ϵ), permeability (K), and effective pore radius ( r pore ). Based on the wick characterization results, small prototype radiator panels (75 × 125 mm) with integrated heat pipe networks were manufactured. Heat rejection performance data are presented from cold thermal vacuum testing, with heat input temperatures up to ~510 K. Future efforts will focus on improving heat pipe performance, optimizing radiator mass, and evaluating larger panels to assess scalability.

nuclear electric propulsion↗

Silane Modification of Glass and Silica Surfaces to Obtain Equally Oil-Wet Surfaces in Glass-Covered Silicon Micromodel Applications

The wettability of silicon and glass surfaces can be modified by silanization. However, similar treatments of glass and silica surfaces using the same silane do not necessarily yield the same wettability as determined by the oil-water contact angle. In this technical note, surface cleaning pretreatments were investigated to determine conditions that would yield oil-wet surfaces on glass with similar wettability to silica surfaces treated with the same silane, and both air-water and oil-water contact angles were determined. Air-water contact angles were less sensitive to differences between silanized silica and glass surfaces, often yielding similar values while the oil-water contact angles were quite different. Borosilicate glass surfaces cleaned with standard cleaning solution 1 (SC1) yield intermediate-wet surfaces when silanized with hexamethyldisilazane, while the same cleaning and silanization yields oil-wet surfaces on silica. However, cleaning glass in boiling concentrated nitric acid creates a surface that can be silanized to obtain oil-wet surfaces using HDMS. Moreover, this method is effective on glass with prior thermal treatment at an elevated temperature of 400oC. In this way, silica and glass can be silanized to obtain equally oil-wet surfaces using HMDS. It is demonstrated that pretreatment and silanization is feasible in silicon-silica/glass micromodels previously assembled by anodic bonding, and that the change in wettability has a significant observable effect on immiscisble fluid displacements in the pore network.

contact angle↗

Cavitation in a soft porous material

Abstract We study the collapse and expansion of a cavitation bubble in a deformable porous medium. We develop a continuum-scale model that couples compressible fluid flow in the pore network with the elastic response of a solid skeleton. Under the assumption of spherical symmetry, our model can be reduced to an ordinary differential equation that extends the Rayleigh–Plesset equation to bubbles in soft porous media. The extended Rayleigh–Plesset equation reveals that finite-size effects lead to the breakdown of the universal scaling relation between bubble radius and time that holds in the infinite-size limit. Our data indicate that the deformability of the porous medium slows down the collapse and expansion processes, a result with important consequences for wide-ranging phenomena, from drug delivery to spore dispersion.

Leng, Yu↗

Fundamental Study of Disposition and Release of Methane in a Shale Gas Reservoir

The recent boom in shale gas production through hydrofracturing has reshaped the energy production landscape in the United States. Wellbore production rates vary greatly among the wells within a single field and decline rapidly with time, thus bring up a serious concern with the sustainability of shale gas production. Shale gas production starts with creating a fracture network by injecting a pressurized fluid in a wellbore. The induced fractures are then held open by proppant particles. During production, gas releases from the mudstone matrix, migrates to nearby fractures, and ultimately reaches a production wellbore. Given the relatively high permeability of the induced fractures, gas release and migration in low-permeability shale matrix is likely to be a limiting step for long-term wellbore production. Therefore, a clear understanding of the underlying mechanisms of methane disposition and release in shale matrix is crucial for the development of new technologies to maximize gas production and recovery. Shale is a natural nanocomposite material with distinct characteristics of nanometer-scale pore sizes, extremely low permeability, high clay contents, significant amounts of organic carbon, and large spatial heterogeneities. Our work has shown that nanopore confinement plays an important role in methane disposition and release in shale matrix. Using molecular simulations, we show that methane release in nanoporous kerogen matrix is characterized by fast release of pressurized free gas (accounting for ~ 30 - 47% recovery) followed by slow release of adsorbed gas as the gas pressure decreases. The first stage is driven by the gas pressure gradient while the second stage is controlled by gas desorption and diffusion. The long-term production decline appears controlled by the second stage of gas release. We further show that diffusion of all methane in nanoporous kerogen behaves differently from the bulk phase, with much smaller diffusion coefficients. The MD simulations also indicate that a significant fraction (3 - 35%) of methane deposited in kerogen can potentially become trapped in isolated nanopores and thus not recoverable. We have successfully established experimental capabilities for measuring gas sorption and desorption on shale and model materials under a wide range of physical and chemical conditions. Both low and high pressure measurements show significant sorption of CH 4 and CO 2 onto clays, implying that methane adsorbed on clay minerals could contribute a significant portion of gas-in-place in an unconventional reservoir. We have also studied the potential impact of the interaction of shale with hydrofracking fluid on gas sorption. We have found that the CH 4 -CO 2 sorption capacity for the reacted sample is systematically lower (by a factor of ~2) than that for the unreacted (raw) sample. This difference in sorption capacity may result from a mineralogical or surface chemistry change of the shale sample induced by fluid-rock interaction. Our results shed a new light on mechanistic understanding gas release and production decline in unconventional reservoirs.

04 OIL SHALES AND TAR SANDS↗

PATHFINDER: Probing Atmospheric Flows in an Integrated and Distributed Environment

PATHFINDER is a software effort to create a flexible, modular, collaborative, and distributed environment for studying atmospheric, astrophysical, and other fluid flows in the evolving networked metacomputer environment of the 1990s. It uses existing software, such as HDF (Hierarchical Data Format), DTM (Data Transfer Mechanism), GEMPAK (General Meteorological Package), AVS, SGI Explorer, and Inventor to provide the researcher with the ability to harness the latest in desktop to teraflop computing. Software modules developed during the project are available in the public domain via anonymous FTP from the National Center for Supercomputing Applications (NCSA). The address is ftp.ncsa.uiuc.edu, and the directory is /SGI/PATHFINDER.

Wilhelmson, R. B.↗

Two-Dimensional Network Flow Modeling of No-Vent Tank Filling of A Cryogenic Tank With Thermo-Dynamic Vent System Assisted Injector

A finite volume-based network flow modeling tool, Generalized Fluid System Simulation Program (GFSSP), was used to simulate the no-vent filling of a cryogenic tank with a thermodynamic vent system assisted injector. The tank was discretized into nodes in the radial and axial directions to calculate two-dimensional axisymmetric flow with heat transfer between the solid wall and fluid. The nodes are connected by branches where momentum equations are solved to calculate the flowrates. The mass and energy conservation equations and the equation of state are solved to calculate pressure, temperature, and resident mass at the nodes. The energy conservation equations at the solid nodes were solved to calculate the temperature of the solid wall. The system of equations was solved by a combination of successive substitution and simultaneous Newton-Raphson method. The numerical model accounts for a) different regimes of pool boiling heat transfer, b) condensation of vapor around the spray droplets, c) condensation of vapor at the interface of the cooled injector and ullage, and d) condensation at the liquid-vapor interphase. The numerical predictions of tank pressure, filling rate, and wall temperatures were compared with the test data. The paper also presents the effect of nodal discretization by comparing the results of the 1D and 2D models and examines the effect of the droplet diameter on tank filling.

Nodal Model↗