Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fluid characterization”

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 307 records · Page 17

Modeling the Nonlinear Rheology of Polymer Additive Manufacturing

This report summarizes molecular and continuum simulation studies focused on developing physics - based predictive models for the evolution of polymer molecular order during the nonlinear processing flows of additive manufacturing. Our molecular simulations of polymer elongation flows identified novel mechanisms of fluid dissipation for various polymer architectures that might be harnessed to enhance material processability. In order to predict the complex thermal and flow history of polymer realistic additive manufacturing processes, we have developed and deployed a high - performance mesh - free hydrodynamics module in Sandia's LAMMPS software. This module called RHEO – short for Reproducing Hydrodynamics and Elastic Objects – hybridizes an updated - Lagrange reproducing - kernel method for complex fluids with a bonded particle method (BPM) to capture solidification and solid objects in multiphase flows. In combination, our two methods allow rapid, multiscale characterization of the hydrodynamics and molecular evolution of polymers in realistic processing geometries.

36 MATERIALS SCIENCE↗

Initial use of Nek5000/Cardinal to improve closure models in Pronghorn

Heat transfer coefficient closure models for pebble bed reactors are built using a data-driven approach by leveraging the capabilities of an Evolutionary Algorithm entitled Particle Swarm Optimization (PSO). In the present work, the Computational Fluid Dynamics code nekRS was used in order to collect the high-fidelity flow data for a core with 1,568 pebbles. To characterize the heat transfer, multiple concentric regions were considered to extract the physical quantities of interest, e.g./ the Reynolds number. The PSO algorithm is employed as part of an inverse problem targeting determine what are the coefficients for a Nusselt number correlation to match the collected data. Such correlation should follow any given format that is defined a priori. Finally, two correlations are proposed, one with an implicit dependence on the pebbles’ wall temperatures and another expressed as a fully explicit correlation depending on the flow conditions and the position within the core. Anyway, given the generic nature of the proposed approach, correlations following different formats could be tested. Preliminary results for the high-fidelity simulation of a fast MSR core are presented. The target Reynolds number is currently 20K, with the expectation that this will increase, pending the availability of further computational resources. These simulations will be used to inform lower fidelity models, including a coarse CFD turbulence model in Pronghorn. Additionally, they will serve as a reference for the RANS models in Nek5000/NekRS.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Subsurface Characterization for Evaluating Geothermal Resource Potential from Existing Oil and Gas Wells in Tuttle, Oklahoma

Oil and gas (O&G) wells often encounter co-produced hot water, possibly suitable for geothermal direct-use applications. The City of Tuttle is located on the eastern part of the Anadarko sedimentary basin in Oklahoma with high heat-in-place potential and recovery capability at depth. This study aims at demonstrating the potential of geothermal energy production for direct-use applications in two public schools and 250 nearby houses in Tuttle via repurposing existing O&G wells. In this scope, geochemistry, geology, and borehole log data were collected and incorporated into a 3D conceptual subsurface model. A digital elevation model (DEM) was used to represent the study area topography with four O&G wells. In addition, hydrogeochemical characteristics of the geothermal fluid and scaling potential were analyzed using ternary diagrams and chemical ratios to develop mixing models. The subsurface geology model indicated that the study area primarily consists of Permian to Mississippian Sandstone and Limestone formations, implying a porosity ranging between 12% and 22%, and a permeability up to 3.90E-14 m2 in certain reservoir levels. The reservoir temperature is expected to be ranging between 80 degrees C to 95 degrees C around 3 km depth with an average temperature gradient of 22.8 degrees C/km. Chemical geothermometers also estimated the reservoir temperature as 90 degrees C. Findings of the chemical model demonstrated that the geothermal fluid is Sodium-Potassium-Chloride-Sulfate type and possibly mixed with shallow groundwater resulting in higher Ca and Mg concentrations and lower Na/K ratio implying lower calcite scaling. These results comprehensively characterize the potential of geothermal resources in the study area and imply that geothermal energy production by repurposing existing O&G wells is suitable for low-temperature direct-use applications.

gas wells↗

Subsurface Characterization for Evaluating Geothermal Resource Potential from Existing Oil and Gas Wells in Tuttle, Oklahoma: Preprint

Oil and gas (O&G) wells often encounter co-produced hot water, possibly suitable for geothermal direct-use applications. The City of Tuttle is located on the eastern part of the Anadarko sedimentary basin in Oklahoma with high heat-in-place potential and recovery capability at depth. This study aims at demonstrating the potential of geothermal energy production for direct-use applications in two public schools and 250 nearby houses in Tuttle via repurposing existing O&G wells. In this scope, geochemistry, geology, and borehole log data were collected and incorporated into a 3D conceptual subsurface model. A digital elevation model (DEM) was used to represent the study area topography with four O&G wells. In addition, hydrogeochemical characteristics of the geothermal fluid and scaling potential were analyzed using ternary diagrams and chemical ratios to develop mixing models. The subsurface geology model indicated that the study area primarily consists of Permian to Mississippian Sandstone and Limestone formations, implying a porosity ranging between 12% and 22%, and a permeability up to 3.90E-14 m2 in certain reservoir levels. The reservoir temperature is expected to be ranging between 80 degrees C to 95 degrees C around 3 km depth with an average temperature gradient of 22.8 degrees C/km. Chemical geothermometers also estimated the reservoir temperature as 90 degrees C. Findings of the chemical model demonstrated that the geothermal fluid is Sodium-Potassium-Chloride-Sulfate type and possibly mixed with shallow groundwater resulting in higher Ca and Mg concentrations and lower Na/K ratio implying lower calcite scaling. These results comprehensively characterize the potential of geothermal resources in the study area and imply that geothermal energy production by repurposing existing O&G wells is suitable for low-temperature direct-use applications.

GEOTHERMAL ENERGY↗

Multiscale Characterization of the Caney Shale — An Emerging Play in Oklahoma

From a hydrocarbon perspective, the Caney Shale has historically been evaluated as a sealing unit, which resulted in limited studies characterizing the rock properties of the Caney Shale and its suitability for hydraulic fracturing. The objective of our research is to help bridge the current knowledge gap through the integration of multiscale laboratory techniques and to characterize the macro- and microscale rock properties of the Caney Shale. We employed an integrated approach for the characterization of the Caney using 200 ft (61 m) of Caney core from a target well in southern Oklahoma. Core observation and petrographic analysis of thin sections were combined to characterize the general rock types and associated fabrics and textures. Mineralogical composition, pore system architecture, and rock fabric were analyzed using x-ray diffraction (XRD), scanning electron microscopy/energy dispersive x-ray spectroscopy (SEM/EDS), and focused ion beam (FIB)-SEM. In addition, rebound hardness and indentation testing were carried out to determine rock hardness (brittleness) and elasticity, respectively. With the integrated multiscale characterization, three mixed carbonate-siliciclastic rock types were identified — mudstone, calcareous siltstone, and silty carbonate — likely representing a spectrum of deposition from low to relatively high energy environments in the distal portions of a ramp system. Silty carbonate contains mostly interparticle pores. The calcareous siltstones and silty mudstones contain a combination of organic matter pores and interparticle pores. Each of the rock types shows unique mineralogical compositions based on XRD. The mudstone lithofacies has the highest clay content and the least carbonate content. Calcareous siltstones show moderate carbonate and clay content. Silty carbonate indicates the highest carbonate content with the least clay content. In an order of mudstone, calcareous siltstone, and silty carbonate, rebound hardness and Young’s modulus show an increasing trend. As a result of rock-fluid interactions, there are potential scaling reactions during completion and production that could ultimately affect permeability and production rates. Overall, the proposed multiscale integration approach is critical for the geologic characterization of most rocks. However, in shale reservoirs dominated by microporosity and microstructure where engineered fractures are expected to provide permeability at a reservoir scale, successful integration is essential. An optimized, integrated geological characterization of the Caney Shale that is well aligned with the engineering designs in drilling, completing, and producing wellbores will ultimately lead to optimal production while providing safe and environmentally responsible operations.

Wang, Yulun↗

Radiative Heat Transport and Optical Characterization of High Temperature Molten Salts

The goal of this project was to advance the field of molten fluoride radiative heat transfer (RHT), through optical property measurements, heat transfer loop construction and operation, and computational fluid dynamics simulations. The following outcomes have been accomplished. A comprehensive review has been conducted to compile all available RHT-relevant optical absorption data for molten fluorides. The various measurement methods have also been compiled and critiqued. Thermal emission measurement methods have been refined and documented. A unique thermal emissivity measurement facility has been constructed and used to measure the emissivity of flowing FLiBe-exposed SS316 samples. A slight increase has been observed in the emissivity compared to unexposed SS316. A dual reflectance-emission setup was partially demonstrated for the purpose of molten fluoride absorption measurement. Preliminary infrared edge data for FLiNaK was obtained. A visible near-infrared (NIR) measurement setup has been constructed and will be used to measure absorption of corrosion products in FLiBe. A new multi-band thermal radiation method using blackbody-weighted spectral bands for the absorption-coefficient was deployed to investigate the thermal performance of the UW Natural Circulation FLiBe Loop under various power levels and operating conditions. The results showed an increasing and decreasing trend in the overall heat transfer with optical thickness with a maximum at optical thickness τ D ≈Ο(1). A forced convection FLiBe loop was constructed for the purpose of studying corrosion and thermal hydraulic phenomena. A vertical cantilever pump is used to achieve up to 5 gpm and 5 m/s flow velocity in hot and cold test sections. A methodology has been established for modeling radiative heat transfer (RHT) in molten salt coolants for the purpose of quantifying the contribution of thermal radiation in conditions relevant to FHR or MSR systems. The comparison between the results obtained for forced and mixed convection showed that in the studied ranges of optical thickness, the maximum effect of RHT on overall heat transfer was 60% improvement observed at z/D = 250 for forced convection; whereas for mixed convection, that value was higher at 170% observed at z/D = 30. The highest effect of RHT on the overall energy distribution was observed at τ D ≈3.4. Negligible effect (<20%) was observed above τ D =18 and below τ D ≈0.002. A CFD code-to-code comparison was conducted using COMSOL Multiphysics at UW-Madison and STAR-CCM at MIT, for a forced and mixed convection laminar flow for a grey medium in a vertical heated using both COMSOL Multiphysics at UW-Madison and STAR-CCM at MIT. The relative error was the largest at τ D =40 for the radiative heat flux, bulk temperature, and Nusselt number. On the other hand, the relative error was the least at τ D =0.002 for the conductive heat flux, bulk temperature, and Nusselt number.

Anderson, Mark↗

Incorporating the molecular-scale into a hydrodynamic description of confined aqueous systems

Hydrodynamics provides a continuum-level description of fluid motion, but its applicability at the nanoscale becomes uncertain due to the emerging importance of molecular-level effects such as spatial heterogeneity. Hydrodynamic boundary conditions that incorporate molecular details allow us to partition the system into a near-wall region and a bulk fluid region. We identify a hydrodynamic wall located inside the fluid that determines where slip begins. By extending the hydrodynamic wall with the slip length, the position of the extrapolated wall is established. This offers a unified description of both slip and stagnant flow behaviors, with wall hydrophobicity characterized by the relative location of the extrapolated wall with respect to the physical wall. Employing this concept in analyses of equilibrium molecular dynamics (MD) and non-equilibrium MD simulations of Couette and Poiseuille flows, our results demonstrate consistency between equilibrium and non-equilibrium approaches across different flow types and confinement levels. This demonstrates the robust nature of linear response theory. We then explore the effects of fluid-wall and bulk fluid interactions on the hydrodynamic properties. These findings enhance the effectiveness of molecular-based simulations for investigating complex confined systems in nanofluidics, biology, and colloidal science, offering a complementary molecular-scale perspective to traditional continuum approaches.

flow simulations↗

Stochastic finite volume method for uncertainty quantification of transient flow in gas pipeline networks

We develop a weakly intrusive framework to simulate the propagation of uncertainty in solutions of generic hyperbolic partial differential equation systems on graph-connected domains with nodal coupling and boundary conditions. The method is based on the Stochastic Finite Volume (SFV) approach and can be applied for uncertainty quantification (UQ) of the dynamical state of fluid flow over actuated transport networks. The numerical scheme has specific advantages for modeling intertemporal uncertainty in time-varying boundary parameters, which cannot be characterized by strict upper and lower (interval) bounds. We describe the scheme for a single pipe, and then formulate the controlled junction Riemann problem (JRP) that enables the extension to general network structures. In conclusion, we demonstrate the method's capabilities and performance characteristics using a standard benchmark test network.

97 MATHEMATICS AND COMPUTING↗

The Local Electronic Structure of Supercritical CO 2 from X-ray Raman Spectroscopy and Atomistic-Scale Modeling

Supercritical CO 2 is encountered in several technical and natural systems related to biology, geophysics, and engineering. While the structure of gaseous CO 2 has been studied extensively, the properties of supercritical CO 2 , particularly close to the critical point, are not well-known. Here, in this work, we combine X-ray Raman spectroscopy, molecular dynamics simulations, and first-principles density functional theory (DFT) calculations to characterize the local electronic structure of supercritical CO 2 at conditions around the critical point. The X-ray Raman oxygen K-edge spectra manifest systematic trends associated with the phase change of CO 2 and the intermolecular distance. Extensive first-principles DFT calculations rationalize these observations on the basis of the 4sσ Rydberg state hybridization. X-ray Raman spectroscopy is found to be a sensitive tool for characterizing electronic properties of CO 2 under challenging experimental conditions and is demonstrated to be a unique probe for studying the electronic structure of supercritical fluids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Wettability Effects on Falling Film Flow and Heat Transfer Over Horizontal Tubes in Jet Flow Mode

The performance of a falling-film heat exchanger is strongly linked to the surface characteristics and the heat transfer processes that take place over the tubes. The primary aim of this numerical study is to characterize the influence of surface wettability on the film flow behavior and its associated surface heat transfer in the jet-flow mode. Volume of fluid (VOF) based simulations are carried out for horizontal tubes with different surface wettabilities. The wettability of the tube surfaces is represented using the Kistler's dynamic contact angle model. Surface wettability effects ranging from superhydrophilic to superhydrophobic are studied by varying the equilibrium contact angle from 2 deg to 175 deg. Two different liquid mass flow rates of 0.06 and 0.18 kg/m-s corresponding to the inline and staggered jet flow modes are studied. Results are presented in terms of the liquid film thickness, the contact areas between the different phases (solid–liquid and liquid–air), and the heat transfer coefficient or Nusselt number. The resistance imposed by the increasing contact angles inhibits the extent of the liquid spreading over the tube surface, and this, in turn, influences the liquid film thickness, and the wetted area of the tube surface. A significant decrement in the heat transfer rate from the tube surfaces was observed as the equilibrium contact angle increased from 2 deg to 175 deg. The local distributions of the Nusselt number over the tube surface are strongly influenced by the flow recirculation in the liquid bulk.

Engineering↗

Utah FORGE: Documentation on Discrete Fracture Network and Fracture Propagation Modelling

This dataset includes reports and a slide presentation on discrete fracture network (DFN) generation and hydraulic fracture modeling at the Utah FORGE site. It details the characterization of natural fractures using well log and core data, as well as stochastic modeling techniques. The reports describe simulations of hydraulic fracture propagation, fluid-mechanical interactions, and induced microseismicity. The dataset also includes history-matching of net pressure and analyses of fracture growth in naturally fractured geothermal reservoirs. The slides summarize key findings and future research directions.

15 GEOTHERMAL ENERGY↗

Geochemical Impact of Acid Spearhead and Slickwater Stimulation on Wolfcamp Shale from the Hydraulic Fracturing Test Site

The Hydraulic Fracturing Test Site 1 (HFTS-1) was a field study performed in the Wolfcamp Formation in the West Texas Permian (Midland) Basin, USA, with a focus on improving the efficiency of hydraulic fracturing. Investigating site-specific rock-fluid geochemical interactions during hydraulic fracturing is an important step to understanding the impact on formation shale porosity, permeability, and long-term shale gas production. During field operations in this region, hydraulic fracturing fluid (HFF) injection usually starts with a concentrated acid spearhead for rapid rock dissolution, followed by the injection of near-neutral pH slickwater containing chemicals and proppants. A multistep sequential injection approach was used to investigate different stages of rock-fluid interactions. The carbonate content in the host rock is important when acid spearhead is considered, as carbonate mineral dissolution is rapid and can result in porosity and permeability changes in the shale matrix. Here, in this study, we designed flow-through experiments using fractured carbonate-rich and clay-rich Wolfcamp shale cores with (1) a short-time acid soaking step and (2) a long-term slickwater flow-through step to simulate the injection method used at HFTS-1. The fluid chemistry was analyzed. A thorough mineralogical progression [e.g., Calcium (Ca) dissolution and iron (Fe) redox progression] in the cores during HFF injection was also characterized and imaged by synchrotron microprobe. Reactive transport modeling was performed based on the experimental setup. The results showed that the acid spearhead is a crucial step in creating a reaction front by mineral dissolution, especially in carbonate-rich shales. A slight layer of ferrihydrite precipitated during the slickwater flow-through period. This study provides insights into potential geochemical impact due to hydraulic fracturing operations in the Permian Basin.

clastic rock↗

Computing the Central Charge of the 3D Ising CFT Using Quantum Finite Elements

The 3D Ising conformal field theory (CFT) describes different physical systems, such as uniaxial magnets or fluids, at their critical points. In absence of an analytical solution for the 3D Ising model, the scaling dimensions and operator product expansion (OPE) coefficients characterizing this CFT must be determined numerically. The currently most-cited values for these quantities have been obtained from the conformal bootstrap, while lattice calculations have so far only produced reliable results for the scaling dimensions involved in calculating the critical exponents. Using Quantum Finite Elements to investigate critical \(\phi^4\)-theory on \(\mathbb{R}\times\mathbb{S}^2\), we show in this work that it is possible to extract scaling dimensions and OPE coefficients of the 3D Ising CFT by fitting the lattice four-point function with expectations from the operator product expansion for the radially quantized CFT and extrapolating to the continuum limit. This way, we have for the first time been able to use Monte Carlo simulations to compute the central charge of the theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The Thermal Response of a Packed Bed Thermal Energy Storage System upon Saturated Steam Injection Using Distributed Temperature Sensing

The effectiveness of a thermal energy storage (TES) system is typically characterized with the help of thermal stratification or temperature gradients along the direction of heat injection, which is typically the flow direction of heat transfer fluid. The steepness of temperature gradients are a direct indicator of the effectiveness or efficiency of the heat storage or dispatch process. The temperature gradient evolution along the packed bed of ceramic particles upon saturated steam injection is presented in this work. Distributed temperature sensing based on optical frequency domain reflectometry was deployed in a packed bed of ceramic particles to capture the thermal front evolution in the axial direction. The physical processes accompanying steam injection in packed beds are complex due to phase change, transitioning two-phase flow, and changes in condensate accumulation. Therefore, the variation of thermal response of the TES system for various steam injection flow rates was experimentally studied using a high-resolution distributed temperature sensing system in a chemically inert alumina particle-packed bed. Distinct zones of different heat transfer modes were observed during the steam injection experiments. A distinct conduction zone, evident from diffuse thermal fronts, was observed at low flow rates, and these thermal gradients became sharper as the flow rate increased. The diffuse thermal fronts in the heat storage media suggest a low exergy efficiency of the TES system, as energy losses started initiating before a significant fraction of the bed was saturated with steam.

25 ENERGY STORAGE↗

Bioderived heat transfer fluids and methods of making the same

The present disclosure relates to a composition that includes a lignin-derived mixture that includes at least one of a dimer, a trimer, and/or a tetramer, where the composition is characterized by a thermal stability up to a maximum temperature between about 260° C. and about 300° C.

Thornburg, Nicholas Earl↗

Microstructure and water absorption of ancient concrete from Pompeii: An integrated synchrotron microtomography and neutron radiography characterization

There is renewed interest in using advanced techniques to characterize ancient Roman concrete due to its exceptional durability and low-carbon footprint. In the present work, samples were drilled from the “Hospitium” in Pompeii and were analyzed by synchrotron microtomography (μCT) and neutron radiography to study how the microstructure, including the presence of induced cracks, affects their water adsorption. The water distribution and absorptivity were quantified by neutron radiography. The 3D crack propagation, pore size distribution and orientation, tortuosity, and connectivity were analyzed from μCT results using advanced imaging methods. Porosity was also measured by mercury intrusion porosimetry (MIP) as a reference. Ductile fracture patterns were observed once cracks were introduced. Compared to Portland cement mortar/concrete, the Pompeii samples had relatively high porosity, low connectivity, and a similar coefficient of capillary penetration. In addition, permeability was predicted from models based on percolation theory and pore structure data to evaluate the fluid transport properties. Understanding the microstructure of ancient Pompeii concrete is important because it could inspire the development of modern concrete with high durability.

36 MATERIALS SCIENCE↗

Characterization of Oil Well Cement Performance during Early Hydration under Simulated Borehole Conditions

Summary Experiments on oil well cement (OWC) slurries were performed using the newly developed laboratory-scale wellbore simulation chamber (WSC). The WSC can simulate hydrostatic pressure reduction in the cemented annulus and possible gas migration under representative conditions. Forensic analysis shows that pressurized fluids can result in porous cement and gas channeling during cement slurry gelation. The effects of different factors on slurry pore pressure were also studied, including formation permeability, initial overburden pressure (OBP) representing the depth of interest, wellbore temperature, water/cement (w/c) ratio, cement composition, and the use of a calcium chloride (CaCl2)-based accelerator. By analyzing the temperature history of hydrating cement using degree of hydration, the evolution of cement hydration was characterized for slurry designs cured at different hydration rates. This provides the opportunity to parameterize the slurry designs and other important factors associated with wellbore conditions.

Engineering↗

Using the Shallow Strain Tensor to Characterize Deep Geologic Reservoirs

Abstract Storing and recovering water, carbon, and heat from geologic reservoirs is central to managing resources in a changing climate. We tested the hypothesis that the strain tensor caused by injecting or producing fluids can be measured at shallow depths and interpreted to advance understanding of underlying deep aquifers or reservoirs. Geodetic‐grade strainmeters were deployed at 30 m depth overlying the Bartlesville Formation, a 500‐m‐deep sandstone near Tulsa, OK. The strainmeters are 220 m east of injection well 9A completed in a permeable lens at the base of the Bartlesville Formation. Water was injected into well 9A at approximately 1.0 L/s during four tests that ranged in duration from a few hours to a few weeks. The horizontal strain increased (tension) and the circumferential strain was a few times larger than the radial strain. The vertical strain decreased (compression) during injection. Strain rates were approximately 100 nε/day during the first few hours, but the rates decreased and were approximately 10 nε/day during most of the tests. Four independent methods of poroelastic simulation and inversion predict reservoir properties and geometries that are similar to each other and consistent with independent information about the reservoir. All strain interpretations predict that a boundary to the permeable lens occurs beneath the vicinity of the strainmeters, which is consistent with core data from the site. The boundary of the permeable lens is located by matching the vertical, radial and circumferential strains, which demonstrates the value of measuring the strain tensor.

Murdoch, Lawrence C.↗