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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Phase I Project: Fiber Optic Distributed Acoustic Sensing for Periodic Hydraulic Tests

The extraction of heat from hot rock requires circulation of fluid through fracture networks. Because the geometry and connectivity of these fractures determines the efficiency of fluid circulation, many tools are used to characterize fractures before and after development of the reservoir. Under this project, a new tool was developed that allows hydraulic connectivity between geothermal boreholes to be identified. Nanostrain in rock fractures is measured using fiber optic distributed acoustic sensing (DAS). This strain is measured in one borehole in response to periodic pressure pulses induced in another borehole. The strain in the fractures represents hydraulic connectivity between wells. DAS is typically used at frequencies of Hz to kHz, but strain at mHz frequencies were measured for this project. The tool was demonstrated in the laboratory and in the field. In the laboratory, strain in fiber optic cables was measured in response to compression due to oscillating fluid pressure. DAS recorded strains as small as 10 picometer/m in response to 1 cm of water level change. At a fractured crystalline rock field site, strain was measured in boreholes. Fiber-optic cable was mechanically coupled borehole walls using pressured flexible liners. In one borehole 30 m from the oscillating pumping source, pressure and strain were measured simultaneously. The DAS system measured fracture displacement at frequencies of less than 1 mHz (18 min periods) and amplitudes of less than 1 nm, in response to fluid pressure changes of less 20 Pa (2 mm of water). The attenuation and phase shift of the monitored strain signal is indicative of the permeability and storage (compliance) of the fracture network that connects the two wells. The strain response as a function of oscillation frequency is characteristic of the hydraulic structure of the formation. This is the first application of DAS to the measurement of low frequency strain in boreholes. It has enormous potential for monitoring geothermal reservoirs for purposes of understanding reservoir compliance and for assuring security of injection fluids. Periodic pressure pulses can be induced by oscillating injection or pumping during operation so the system could provide real time reservoir data. DAS cable may already be installed at a site as it is becoming increasingly used for seismic observation. Simulations conducted for this project indicate that strain should be propagated through borehole cements, so observations can be made outside of well casing. In uncased holes, the cable would need to be mechanically coupled to the borehole wall to provide measurements. One option would be to install fiber into cemented and abandoned boreholes to extend their utility.

15 GEOTHERMAL ENERGY↗
DOE OSTI2016

RELATIONSHIPS BETWEEN FLUID VORTICITY, KINETIC HELICITY, AND MAGNETIC FIELD ON SMALL-SCALES (QUIET-NETWORK) ON THE SUN

We derive horizontal fluid motions on the solar surface over large areas covering the quiet-Sun magnetic network from local correlation tracking of convective granules imaged in continuum intensity and Doppler velocity by the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory . From these we calculate the horizontal divergence, the vertical component of vorticity, and the kinetic helicity of fluid motions. We study the correlations between fluid divergence and vorticity, and between vorticity (kinetic helicity) and the magnetic field. We find that the vorticity (kinetic helicity) around small-scale fields exhibits a hemispherical pattern (in sign) similar to that followed by the magnetic helicity of large-scale active regions (containing sunspots). We identify this pattern to be a result of the Coriolis force acting on supergranular-scale flows (both the outflows and inflows), consistent with earlier studies using local helioseismology. Furthermore, we show that the magnetic fields cause transfer of vorticity from supergranular inflow regions to outflow regions, and that they tend to suppress the vortical motions around them when magnetic flux densities exceed about 300 G (from HMI). We also show that such an action of the magnetic fields leads to marked changes in the correlations between fluid divergence and vorticity. These results are speculated to be of importance to local dynamo action (if present) and to the dynamical evolution of magnetic helicity at the small-scale.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1987

FLUINT - General fluid system analysis with SINDA '85

This paper introduces FLUINT (FLUid INTegrator), a general fluid system simulation program that works in conjunction with SINDA '85 (Systems Improved Numerical Differencing Analyzer, 1985 version). FLUINT solves arbitrary networks representing internal fluid systems while SINDA '85 simultaneously solves traditional thermal networks. This transportable program is intended to provide spacecraft thermal management engineers with the analytic tools needed to simulate single- and two-phase coolant loops. The fundamental concepts for fluid system modeling are described, and methods of modeling complex, hardware-specific components are outlined. Underlying correlations and solution methods are briefly discussed. Example analyses and comparisons to test data are provided.

Cullimore, B. A.↗
DOE OSTI2021

Impacts of Pore Network-Scale Wettability Heterogeneity on Immiscible Fluid Displacement: A Micromodel Study

The mixed-wet nature of reservoir formations imposes a wide range of rock wettability from strong resident-fluid wetting to strong invading-fluid wetting. The characteristics of two-phase flow in porous media composed of mixed-wetting surfaces remain poorly understood. In this study, we investigated the displacement of resident ethylene glycol (EG) by hexane in two mixed-wet micromodels of identical 2.5-D geometry heterogeneity, with uniformly or heterogeneously distributed patches strongly wetting to hexane. These patches are mixed among pores with unaltered EG-wetting surfaces. Along with control tests in the originally EG-wet micromodel, here we show the classic fingering and transitions in flow regimes at $logCa$ (capillary number) from -7.2 to -3.9. Moreover, pore-scale distributions of wettability and their spatial correlation influence displacement efficiency. In the two mixed-wet micromodels, we found (a) an increase of steady-state hexane saturation at the end of experiments by up to 0.12 in the capillary fingering regime and a decrease of at most by 0.06 in the viscous fingering regime, compared to the EG-wet micromodel, and (b) dispersed and fragmented hexane distribution after displacement. Brine drainage during supercritical CO 2 (scCO 2 ) injections in these micromodels occurs with lower wettability contrasts, and under similar viscosity ratios and interfacial tensions resulted in higher displacement efficiency relative to displacement of EG by hexane. While mixed-wettability can enhance displacement efficiency compared to uniform wettability, the dynamics of immiscible fluids in strong mixed-wet reservoirs are expected to be less pronounced in contributing to the efficiency of geological CO 2 sequestration, oil recovery, and remediation of hydrocarbon-contaminated aquifers.

2.5-D micromodel↗
DOE OSTI2023

Gravity-driven controls on fluid and carbonate precipitation distributions in fractures

Abstract Many challenges related to carbon-dioxide ( $$\hbox {CO}_2$$ CO 2 ) sequestration in subsurface rock are linked to the injection of fluids through induced or existing fracture networks and how these fluids are altered through geochemical interactions. Here, we demonstrate that fluid mixing and carbonate mineral distributions in fractures are controlled by gravity-driven chemical dynamics. Using optical imaging and numerical simulations, we show that a density contrast between two miscible fluids causes the formation of a low-density fluid runlet that increases in areal extent as the fracture inclination decreases from 90 $$^\circ$$ ∘ (vertical fracture plane) to 30 $$^\circ$$ ∘ . The runlet is sustained over time and the stability of the runlet is controlled by the gravity-driven formation of 3D vortices that arise in a laminar flow regime. When homogeneous precipitation was induced, calcium carbonate covered the entire surface for horizontal fractures (0 $$^\circ$$ ∘ ). However, for fracture inclinations greater than 10 $$^\circ$$ ∘ , the runlet formation limited the areal extent of the precipitation to less than 15% of the fracture surface. These insights suggest that the ability to sequester $$\hbox {CO}_2$$ CO 2 through mineralization along fractures will depend on the fracture orientation relative to gravity, with horizontal fractures more likely to seal uniformly.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1984

Freezing of Martian streams under climatic conditions

The valley networks of Mars are widely believed to have formed at a time when climatic conditions on the planet were significantly different from those that currently prevail. This view arises from the following observations: (1) the valleys form integrated branching networks which suggests fluid drainage, and water is the most plausible fluid, (2) the present atmosphere contains only minute amounts of water, (3) the networks appear to be more akin to terrestrial valleys that are eroded by streams of modest discharges than features that form by catastrophic floods, and (4) small streams of water will rapidly freeze under present climatic conditions. Climatic conditions at the time of formation of the valleys are studied based on the assumption that they were cut by running water.

Carr, M. H.↗
NASA NTRS2016

Generalized Fluid System Simulation Program, Version 6.0

The Generalized Fluid System Simulation Program (GFSSP) is a general purpose computer program for analyzing steady state and time-dependent flow rates, pressures, temperatures, and concentrations in a complex flow network. The program is capable of modeling real fluids with phase changes, compressibility, mixture thermodynamics, conjugate heat transfer between solid and fluid, fluid transients, pumps, compressors, and external body forces such as gravity and centrifugal. The thermofluid system to be analyzed is discretized into nodes, branches, and conductors. The scalar properties such as pressure, temperature, and concentrations are calculated at nodes. Mass flow rates and heat transfer rates are computed in branches and conductors. The graphical user interface allows users to build their models using the 'point, drag, and click' method; the users can also run their models and post-process the results in the same environment. Two thermodynamic property programs (GASP/WASP and GASPAK) provide required thermodynamic and thermophysical properties for 36 fluids: helium, methane, neon, nitrogen, carbon monoxide, oxygen, argon, carbon dioxide, fluorine, hydrogen, parahydrogen, water, kerosene (RP-1), isobutene, butane, deuterium, ethane, ethylene, hydrogen sulfide, krypton, propane, xenon, R-11, R-12, R-22, R-32, R-123, R-124, R-125, R-134A, R-152A, nitrogen trifluoride, ammonia, hydrogen peroxide, and air. The program also provides the options of using any incompressible fluid with constant density and viscosity or ideal gas. The users can also supply property tables for fluids that are not in the library. Twenty-four different resistance/source options are provided for modeling momentum sources or sinks in the branches. These options include pipe flow, flow through a restriction, noncircular duct, pipe flow with entrance and/or exit losses, thin sharp orifice, thick orifice, square edge reduction, square edge expansion, rotating annular duct, rotating radial duct, labyrinth seal, parallel plates, common fittings and valves, pump characteristics, pump power, valve with a given loss coefficient, Joule-Thompson device, control valve, heat exchanger core, parallel tube, and compressible orifice. The program has the provision of including additional resistance options through User Subroutines. GFSSP employs a finite volume formulation of mass, momentum, and energy conservation equations in conjunction with the thermodynamic equations of state for real fluids as well as energy conservation equations for the solid. The system of equations describing the fluid network is solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods. The application and verification of the code has been demonstrated through 30 example problems.

Majumdar, A. K.↗
DOE OSTI2019

Multivascular networks and functional intravascular topologies within biocompatible hydrogels

Solid organs transport fluids through distinct vascular networks that are biophysically and biochemically entangled, creating complex three-dimensional (3D) transport regimes that have remained difficult to produce and study. We establish intravascular and multivascular design freedoms with photopolymerizable hydrogels by using food dye additives as biocompatible yet potent photoabsorbers for projection stereolithography. We demonstrate monolithic transparent hydrogels, produced in minutes, comprising efficient intravascular 3D fluid mixers and functional bicuspid valves. We also elaborate entangled vascular networks from space-filling mathematical topologies and explore the oxygenation and flow of human red blood cells during tidal ventilation and distension of a proximate airway. Moreover, we deploy structured biodegradable hydrogel carriers in a rodent model of chronic liver injury to highlight the potential translational utility of this materials innovation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Demand Response for Thermal Fairness in District Heating Networks

District Heating (DH) networks are complex thermal grids wherein a centrally heated fluid is circulated through a network of pipes and heat exchangers to meet the heating needs of residential and commercial buildings. Several factors can lead to inefficient and unfair energy distribution among consumers in these networks. These include varying levels of building insulation, distance of individual buildings from the central energy source, and thermal losses in network pipes. Moreover, shortage of energy at the central energy source and extreme weather conditions can exacerbate these issues, leading to differing levels of thermal comfort and customer disgruntlement in the long run. In this paper, we propose and study a demand response scheme that attempts to ensure thermal fairness among energy consumers in modern thermal grids. We develop optimization formulations based on thermodynamic models of DH networks, which determine optimal energy flow for individual buildings in order to achieve thermal fairness across the network. Our experimental results using physics based models for DH networks show that it is possible to achieve network level thermal fairness based social welfare objectives by controlling network parameters such as mass flow rates of water to the consumer premises and the supply water temperature, even when some network parameters like building thermal characteristics remain unknown.

Bhattacharya, Saptarshi↗
NASA NTRS2005

Numerical Modeling of Flow Distribution in Micro-Fluidics Systems

This paper describes an application of a general purpose computer program, GFSSP (Generalized Fluid System Simulation Program) for calculating flow distribution in a network of micro-channels. GFSSP employs a finite volume formulation of mass and momentum conservation equations in a network consisting of nodes and branches. Mass conservation equation is solved for pressures at the nodes while the momentum conservation equation is solved at the branches to calculate flowrate. The system of equations describing the fluid network is solved by a numerical method that is a combination of the Newton-Raphson and successive substitution methods. The numerical results have been compared with test data and detailed CFD (computational Fluid Dynamics) calculations. The agreement between test data and predictions is satisfactory. The discrepancies between the predictions and test data can be attributed to the frictional correlation which does not include the effect of surface tension or electro-kinetic effect.

Majumdar, Alok↗
DOE OSTI2014

Failure of cap-rock seals as determined from mechanical stratigraphy, stress history, and tensile-failure analysis of exhumed analogs

In this study, the sedimentologic and tectonic histories of clastic cap rocks and their inherent mechanical properties control the nature of permeable fractures within them. The migration of fluid through mm- to cm-scale fracture networks can result in focused fluid flow allowing hydrocarbon production from unconventional reservoirs or compromising the seal integrity of fluid traps. To understand the nature and distribution of subsurface fluid-flow pathways through fracture networks in cap-rock seals we examine four exhumed Paleozoic and Mesozoic seal analogs in Utah. We combine these outcrop analyses with subsidence analysis, paleoloading histories, and rock-strength testing data in modified Mohr–Coulomb–Griffith analyses to evaluate the effects of differential stress and rock type on fracture mode.

58 GEOSCIENCES↗
DOE OSTI2022

Flow Reduction in Pore Networks of Packed Silica Nanoparticles: Insights from Mesoscopic Fluid Models

A modified many-body dissipative particle dynamics (mDPD) model is rigorously calibrated to achieve realistic fluid–fluid/solid interphase properties and applied for mesoscale flow simulations to elucidate the transport mechanisms of heptane liquid and water, respectively, through pore networks formed by packed silica nanoparticles with a uniform diameter of 30 nm. Two million CPU core hours were used to complete the simulation studies. Results show reduction of permeability by 54–64% in heptane flow and by 88–91% in water flow, respectively, compared to the Kozeny–Carman equation. In these nanopores, a large portion of the fluids are in the near-wall regions and thus not mobile due to the confinement effect, resulting in reduced hydraulic conductivity. Moreover, intense oscillations in the calculated flow velocities also indicate the confinement effect that contests the external driven force to flow. Here, the generic form of Darcy’s law is considered valid for flow through homogeneous nanopore networks, while permeability depends collectively on pore size and surface wettability. This fluid-permeability dependency is unique to flow in nanopores. In addition, potential dependence of permeability on pore connectivity is observed when the porosity remains the same in different core specimens.

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

Performance analysis and acceleration of explicit integration for large kinetic networks using batched GPU computations

Abstract—We demonstrate the systematic implementation of recently-developed fast explicit kinetic integration algorithms that solve efficiently N coupled ordinary differential equations (subject to initial conditions) on modern GPUs. We take representative test cases (Type Ia supernova explosions) and demonstrate two or more orders of magnitude increase in efficiency for solving such systems (of realistic thermonuclear networks coupled to fluid dynamics). This implies that important coupled, multiphysics problems in various scientific and technical disciplines that were intractable, or could be simulated only with highly schematic kinetic networks, are now computationally feasible. As examples of such applications we present the computational techniques developed for our ongoing deployment of these new methods on modern GPU accelerators. We show that similarly to many other scientific applications, ranging from national security to medical advances, the computation can be split into many independent computational tasks, each of relatively small-size. As the size of each individual task does not provide sufficient parallelism for the underlying hardware, especially for accelerators, these tasks must be computed concurrently as a single routine, that we call batched routine, in order to saturate the hardware with enough work.

Haidar Ahmad, Azzam↗
DOE OSTI2022

Assessments of epistemic uncertainty using Gaussian stochastic weight averaging for fluid-flow regression

Here, we use Gaussian stochastic weight averaging (SWAG) to assess the epistemic uncertainty associated with neural-network-based function approximation relevant to fluid flows. SWAG approximates a posterior Gaussian distribution of each weight, given training data, and a constant learning rate. Having access to this distribution, it is able to create multiple models with various combinations of sampled weights, which can be used to obtain ensemble predictions. The average of such an ensemble can be regarded as the 'mean estimation', whereas its standard deviation can be used to construct 'confidence intervals', which enable us to perform uncertainty quantification (UQ) with regard to the training process of neural networks. We utilize representative neural-network-based function approximation tasks for the following cases: (i) a two-dimensional circular-cylinder wake; (ii) the DayMET dataset (maximum daily temperature in North America); (iii) a three-dimensional square-cylinder wake; and (iv) urban flow, to assess the generalizability of the present idea for a wide range of complex datasets. SWAG-based UQ can be applied regardless of the network architecture, and therefore, we demonstrate the applicability of the method for two types of neural networks: (i) global field reconstruction from sparse sensors by combining convolutional neural network (CNN) and multi-layer perceptron (MLP); and (ii) far-field state estimation from sectional data with two-dimensional CNN. We find that SWAG can obtain physically-interpretable confidence-interval estimates from the perspective of epistemic uncertainty. This capability supports its use for a wide range of problems in science and engineering.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Development of a general method for obtaining the geometry of microfluidic networks

In the present study, a general method for geometry of fluidic networks is developed with emphasis on pressure-driven flows in the microfluidic applications. The design method is based on general features of network's geometry such as cross-sectional area and length of channels. Also, the method is applicable to various cross-sectional shapes such as circular, rectangular, triangular, and trapezoidal cross sections. Using constructal theory, the flow resistance, energy loss and performance of the network are optimized. Also, by this method, practical design strategies for the fabrication of microfluidic networks can be improved. The design method enables rapid prediction of fluid flow in the complex network of channels and is very useful for improving proper miniaturization and integration of microfluidic networks. Minimization of flow resistance of the network of channels leads to universal constants for consecutive cross-sectional areas and lengths. For a Y-shaped network, the optimal ratios of consecutive cross-section areas (A{sub i+1}/A{sub i}) and lengths (L{sub i+1}/L{sub i}) are obtained as A{sub i+1}/A{sub i} = 2{sup −2/3} and L{sub i+1}/L{sub i} = 2{sup −1/3}, respectively. It is shown that energy loss in the network is proportional to the volume of network. It is also seen when the number of channels is increased both the hydraulic resistance and the volume occupied by the network are increased in a similar manner. Furthermore, the method offers that fabrication of multi-depth and multi-width microchannels should be considered as an integral part of designing procedures. Finally, numerical simulations for the fluid flow in the network have been performed and results show very good agreement with analytic results.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Nanopore networks in colloidal silica assemblies characterized by XCT for confined fluid flow modeling

Fluid flow through nanopore structures has exhibited different behavior than that described by Darcy's law, derived from pore networks of micrometer scale. Here, pressure assemblies of colloidal silica spheres of specified diameters from 94 to 620 nanometers were prepared as model nano-porous media. The smaller silica spheres tended to form random packing, but some substantial ordering was found in the packing of the larger silica spheres. X-ray computed tomography (XCT), with a voxel resolution of 16 nm, was used to characterize the nanopore networks. The porosity of the Nano-XCT pore network (34%) was less than the total porosity of the assemblies as determined by X-ray attenuation (52%), possibly due to the limits of voxel resolution during segmentation. Fluid flow in the nanopore networks was simulated using the single-phase Lattice Boltzmann Method, and the simulated permeability was compared with the empirical and experimental values. Based on our characterization, although a well-ordered packing of silica spheres was not achieved, it was found that the nanopore networks in the colloidal silica assemblies had pore size distributions corresponding to the particle sizes. The simulated permeability was less than the experimental measurement for water flow, but the complex packing of silica spheres and surface chemistry issues need to be considered in future research.

02 PETROLEUM↗
DOE OSTI2015

Reactive fluid flow in CO 2 storage reservoirs: A 2‐D pore network model study

Abstract Water acidifies in the presence of CO 2 and prompts mineral dissolution. A 2‐D pore network model scheme is developed to investigate reactive fluid flow in CO 2 storage reservoirs during injection when advective transport prevails. Mineral dissolution satisfies kinetic rate laws and continues until thermodynamic equilibrium is reached. In advection‐dominant regimes, network simulation results show that species concentration, tube enlargement and flow rate can be summarized in terms of the dimensionless Damköhler number Da which is the ratio between advection time along a pore and the reaction time. Reservoirs will tend to experience localized enlargement near injection wells (before water drying) and compact dissolution in the far‐field. The Damköhler number couples with initial pore‐size variability to distort the relationship between mean tube diameter and either local or network‐average flow rates. Both the Damköhler number and pore‐size variability should be considered in field‐scale numerical simulators.

Kim, Seunghee↗