Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluid Properties”

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 451 records · Page 25

Fuel Property Effects on Knock Propensity and Thermal Efficiency in a Direct-Injection Spark-Ignition Engine

Engine knock remains one of the major barriers to further improvement in thermal efficiency of Direct-Injection Spark-Ignition (DISI) engines. While Research Octane Number and Motor Octane Number are often used as standard rating methods for knock resistance of fuels, the impacts of other fuel properties on knock propensity in modern engines such as heat of vaporization (HoV) and laminar flame speed (LFS) require better understanding in order to co-optimize fuels and engine designs to achieve higher thermal efficiency and lower CO2 emission. In the present study, computational fluid dynamics (CFD) is used to model a boosted DISI engine with a focus on knock prediction and fuel property effects. A level-set G-equation model is employed to capture turbulent premixed combustion, and is coupled with a transported Livengood-Wu (L-W) integral approach to predict autoignition in the unburnt region. A criterion associated with the L-W integral is developed to accurately predict knock onset and knock-limited spark-advance. This model is then applied to a sensitivity analysis of HoV and LFS on knock tendency and thermal efficiency. The pressure-temperature trajectory framework is applied and extended to study the fuel effects on auto-ignition process in the engine. An existing efficiency-based merit function, which is derived from experiments for boosted SI engines, is evaluated and improved based on the current CFD results.

DISI↗

Modeling Adsorption in Silica Pores via Minkowski Functionals and Molecular Electrostatic Moments

Capillary condensation phenomena are important in various technological and environmental processes. Using molecular simulations, we study the confined phase behavior of fluids relevant to carbon sequestration and shale gas production. As a first step toward translating information from the molecular to the pore scale, we express the thermodynamic potential and excess adsorption of methane, nitrogen, carbon dioxide, and water in terms of the pore’s geometric properties via Minkowski functionals. This mathematical reconstruction agrees very well with molecular simulations data. Our results show that the fluid molecular electrostatic moments are positively correlated with the number of adsorption layers in the pore. Moreover, stronger electrostatic moments lead to adsorption at lower pressures. These findings can be applied to improve pore-scale thermodynamic and transport models.

42 ENGINEERING↗

Convergence of micro-geochemistry and micro-geomechanics towards understanding proppant shale rock interaction: A Caney shale case study in southern Oklahoma, USA

As a direct outcome of economic development coupled with an increase in population, global energy demand will continue to rise in the coming decades. Although renewable energy sources are increasingly investigated for optimal production, the immediate needs require focus on energy sources that are currently available and reliable, with a minimal environmental impact; the efficient exploration and production of unconventional hydrocarbon resources is bridging the energy needs and energy aspirations, during the current energy transition period. The main challenges are related to the accurate quantification of the critical rock properties that influence production, their heterogeneity and the multiscale driven physico-chemical nature of rock–fluid interactions. A key feature of shale reservoirs is their low permeability due to dominating nanoporosity of the clay-rich matrix. As a means of producing these reservoirs in a cost-effective manner, a prerequisite is creation of hydraulic fracture networks capable of the highest level of continued conductivity. Fracturing fluid chemical design, formation brine geochemical composition, and rock mineralogy all contribute to swelling-induced conductivity damage. The Caney Shale is an organic-rich, often calcareous mudrock. Many studies have examined the impact that clay has on different kinds of shale productivity but there is currently no data reported on the Caney Shale in relation to horizontal drilling; all reported data on the Caney Shale is on vertical wells which are shallow, compared to an emerging play that is at double the depth. Here, in this work we develop geochemical–geomechanical integration of rock properties at micro-and nanoscales that can provide insights into the potential proppant embedment and its mitigation. The novel methodology amalgamates the following: computed X-ray tomography, scanning electron microscopy, energy dispersive spectroscopy, micro-indentation, and Raman spectroscopy techniques. Here, our results show that due to the multiscale heterogeneity in the Caney Shale, these geochemical and structural properties translate into a variation in mechanical properties that will impact interaction between the proppant and the host shale rock.

03 NATURAL GAS↗

Splitting tensile strength of shale cores: intact versus fractured and sealed with ureolysis-induced calcium carbonate precipitation (UICP)

Ureolysis-induced calcium carbonate precipitation (UICP) is a biomineral solution where the urease enzyme converts urea and calcium into calcium carbonate. The resulting biomineral can bridge gaps in fractured shale, reduce undesired fluid flow, limit fracture propagation, better store carbon dioxide, and potentially enhance well efficiency. The mechanical properties of shale cores were investigated using a modified Brazilian indirect tensile strength test. An investigation of intact shale using Eagle Ford and Wolfcamp cores was conducted at varying temperatures. Results show no significant difference between shale types (average tensile strength = 6.19 MPa). Eagle Ford displayed higher strength at elevated temperature, but temperature did not influence Wolfcamp. Comparatively, cores with a single, lengthwise heterogeneous fracture were sealed with UICP and further tested for tensile strength. UICP was delivered via a flow-through method which injected 20–30 sequential patterns of ureolytic microorganisms and UICP-promoting fluids into the fracture until permeability reduced by three orders of magnitude or with an immersion method which placed cores treated with guar gum and UICP-promoting fluids into a batch reactor, demonstrating that guar gum is a suitable inclusion and may reduce the number of flow-through injections required. Tensile results for both delivery methods were variable (0.15–8 MPa), and in some cores the biomineralized fracture split apart, possibly due to insufficient sealing and/or heterogeneity in the composite UICP-shale cores. Notably in other cores the biomineralized fracture remained intact, demonstrating more cohesion than the surrounding shale, indicating that UICP may produce a strong seal for subsurface application.

58 GEOSCIENCES↗

Using defects as a ‘fossil record’ to help interpret complex processes during additive manufacturing: as applied to raster-scanned electron beam powder bed additively manufactured Ti–6Al–4V

Abstract Defects in parts produced by additive manufacturing, instead of simply being perceived as deleterious, can act as important sources of information associated with the complex physical processes that occur during materials deposition and subsequent thermal cycles. Indeed, they act as materials-state ‘fossil’ records of the dynamic AM process. The approach of using defects as epoch-like records of prior history has been developed while studying additively manufactured Ti–6Al–4V and has given new insights into processes that may otherwise remain either obscured or unquantified. Analogous to ‘epochs,’ the evolution of these defects often is characterized by physics that span across a temporal length scale. To demonstrate this approach, a broad range of analyses including optical and electron microscopy, X-ray computed tomography, energy-dispersive spectroscopy, and electron backscatter diffraction have been used to characterize a raster-scanned electron beam Ti–6Al–4V sample. These analysis techniques provide key characteristics of defects such as their morphology, location within the part, complex compositional fields interacting with the defects, and structures on the free surfaces of defects. Observed defects have been classified as banding, spherical porosity, and lack of fusion. Banding is directly related to preferential evaporation of Al, which has an influence on mechanical properties. Lack-of-fusion defects can be used to understand columnar grain growth, fluid flow of melt pools, humping, and spattering events. Graphical abstract

36 MATERIALS SCIENCE↗

Transient mixing behavior of a supercritical fluid injected into supercritical and subcritical environments

Transient, near-field mixing of a supercritical fluid injected into supercritical and subcritical environments is investigated experimentally. Fuel (propane) from a 100 μm diameter Bosch diesel injector is injected into a high-pressure test rig with a steady flow of heated, low-velocity nitrogen to control the ambient pressure and temperature. The fuel supply pressure is approximately 190 bar, allowing short injections lasting 3–4 ms to be performed in quick succession to acquire statistically significant datasets to represent the unsteady fuel–air mixing. Two test cases with the same ambient density and injector inlet pressure and temperature were chosen in an attempt to isolate effects of mixture phase behavior. Microscopic rainbow Schlieren deflectometry (RSD) is used for visualization and quantitative data acquisition. In the case of the supercritical ambient, the injector inlet and exit and ambient conditions are all supercritical. The RSD analysis reveals a mixing process similar to the single phase fluid mixing, though condensation is also observed locally, and can be explained by the adiabatic mixing theory presented in the paper. For the subcritical ambient case, the injector inlet is supercritical, but the injector exit and ambient conditions are both subcritical. The results show a rapid expansion in the near-field promoting mixing between propane and nitrogen. At downstream locations, significant condensation is also observed as predicted by the adiabatic mixing theory. As a result, the supercritical ambient case exhibited longer and narrower jet penetration compared to the subcritical ambient case.

02 PETROLEUM↗

Shivalik Plasma Device-I, a glow discharge device to study the collective dynamics of dusty plasma

Using a uniquely configured glow discharge-based Shivalik Plasma Device-I, we demonstrate a variety of collective phenomena in dusty plasma away from the glow discharge region. The cylindrical glass device produces plasma using parallel disc-shaped electrodes with a smaller anode size than the cathode. The dust microparticles are initially sprinkled over the grounded cathode. These particles acquire a significant negative charge upon plasma formation, resulting in their levitation due to the balance between the Coulomb force and gravity. The new device supports the levitation of a big-sized (10 × 8 × 5 cm 3 ) three-dimensional dust cloud over the glass surface. It contrasts the dusty plasma formations in-between electrodes reported earlier. As the discharge voltage varies from high to low, the dust cloud travels from over the glass surface to between the electrodes. A complex interplay of dust void over the cathode, a sharp density gradient, and gravity lead to self-excitation of collective dust phenomena. It includes dust density waves (phase velocity, v ph ~4 cm/s), dust cloud oscillation (frequency, f = 5 Hz), sheared flow (flow velocity, v f ~ 1 cm/s), and multiple-sized dust vortices. These dust vortices provided an excellent platform for studying turbulent mixing phenomena. The power spectrum analysis agreed with two-dimensional Kolmogorov power-law scaling. This is an ideal dusty plasma apparatus where we can create or move the dust cloud to a location of choice from the glass surface to in-between the electrodes and excite one among many collective dust dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A novel means to generate high pressure

Diamond anvil cells are the most popular means of generating pressures above 2 GPa. However, in many experiments, such as nuclear magnetic resonance and x-ray absorption, the metallic pressurizing gasket (which confines much of the sample) represents an occluding barrier that requires a low Z gasket material (e.g., Be), a split gasket, or other means to enable better coupling of the sample to electromagnetic radiation. In this paper, we demonstrate a novel method for generating high pressures that confines the sample just above the plane of the gasket by using a diamond with a laser hole drilled into the center of the tip. The sample is then confined by the hole, which is sealed by a flat gasket that fits over the hole. When load is applied to the diamonds, metal flows from the deformed gasket into the hole thereby pressurizing the sample similarly to how a piston pressurizes gas inside a cylinder. The pressurized sample is above the metallic gasket plane just inside the tip of the diamond, and thus easily accessible via x rays or visible light that skims just above the plane of the gasket providing an enhanced aperture of radiation collection. Furthermore, we have demonstrated the utility of this method by obtaining Raman spectra of SnC 2 O 4 and x-ray diffraction spectra of seleno-DL-cystine, all at high pressures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Deuterium and tritium anomalous transport in the generalized Hasegawa–Wakatani resistive drift wave turbulence model with finite ion Larmor radius

Anomalous transport of mixed deuterium–tritium plasma in the edge of magnetic fusion reactors is investigated using numerical solutions of resistive drift wave turbulence model equations, including finite Larmor radius effects, that are derived within the generalized Hasegawa–Wakatani framework. The anomalous cross field diffusivities of deuterium and tritium are compared in turbulence regimes with different values of the electron adiabaticity parameter controlling the existence of zonal flow. The dependence of the tritium-to-deuterium diffusivity ratio on the deuterium and tritium densities and the logarithmic density gradients is analyzed, and a scaling relation is obtained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantum kinetic modeling of KEEN waves in a warm-dense regime

We report the first fully kinetic, quantum study of kinetic electrostatic electron nonlinear (KEEN) waves, showing that quantum diffraction systematically erodes the classical trapping mechanism, narrows harmonic locking to the fundamental, and hastens post-drive decay. Electrons are evolved with a second-order Strang-split 1D1V Wigner–Poisson solver that couples conservative semi-Lagrangian WENO advection to an analytic Fourier space update for the non-local Wigner term, while ions remain classical. We focus on collisionless dynamics in a weakly coupled regime, providing a controlled baseline before collisional extensions. Short, frequency-tuned ponderomotive pulses drive KEEN formation in a uniform Maxwellian plasma; as the dimensionless quantum parameter H rises from the classical limit to values relevant to warm-dense matter, doped semiconductors, and 2D electron systems, the drive threshold increases, higher harmonics are damped, trapped electron vortices diffuse, and the subplasma electrostatic energy relaxes to a lower stationary level, as confirmed by continuous wavelet analysis. These microscopic changes carry macroscopic weight. Ignition-scale capsules now compress matter to regimes where the electron de Broglie wavelength rivals the Debye length, making classical kinetic descriptions insufficient. By extending KEEN physics into this quantum domain, our results offer a potential diagnostic of non-equilibrium electron dynamics for next-generation inertial-confinement designs and high-energy-density platforms, indicating that predictive fusion modeling may benefit from the integration of kinetic fidelity with quantum effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Statistical properties of homogeneous and isotropic turbulence in He II measured via particle tracking velocimetry

Despite being a quantum two-fluid system, superfluid helium-4 (He II) is observed to behave similarly to classical fluids when a flow is generated by mechanical forcing. This similarity has brought up the feasibility of utilizing He II for high Reynolds number classical turbulence research, considering the small kinematic viscosity of He II. However, it has been suggested that the nonclassical dissipation mechanism in He II at small scales may alter its turbulent statistics and intermittency. In this work, we report our study of a nearly homogeneous and isotropic turbulence (HIT) generated by a towed grid in He II. We measure the velocity field using particle tracking velocimetry with solidified deuterium particles as the tracers. By correlating the velocities measured simultaneously on different particle trajectories or at different times along the same particle trajectory, we are able to conduct both Eulerian and Lagrangian flow analyses. Spatial velocity structure functions obtained through the Eulerian analysis show scaling behaviors in the inertial subrange similar to that for classical HIT but with enhanced intermittency. The Lagrangian analysis allows us to examine the flow statistics down to below the dissipation length scale. Interestingly, abnormal deviations from the classical scaling behaviors are observed in this regime. Lastly, we discuss how these deviations may relate to the motion of quantized vortices in the superfluid component in He II.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Bubble Mass Transport Measurement in the Large-Scale Test Loops at Oak Ridge National Laboratory

Molten salts are complex fluids that incorporate multi-phase behavior depending on the chemistry and the physical properties of the entrained components. These components include the carrier salt, the actinide fuel, and fission, activation, and corrosion products, the concentrations of which depends on the burnup history of the salt. Radionuclide transport from the salt into the cover gas / off-gas system depends on volatility as predicted by thermochemistry, but data from the Molten Salt Reactor Experiment (MSRE) conducted in the late 1960s suggest that bubble formation and transport are also important. Anomalously high fractions of noble metals were found in the off-gas system and were attributed to transportation with parent salt aerosols and their association with rising noble gas bubbles. The prediction of such phenomena requires coupled neutronic, thermal hydraulic, and chemical equilibrium calculations, the framework of which is being developed within the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program supported by the US Department of Energy (DOE) Office of Nuclear Energy (NE). However, separate effects tests and models of experiments using tools such as SAM and Thermochimica can guide model developers through important processes. Conversely, model development guides the choice of experiments and systems to provide data that are relevant for validation. This report describes several experiments that tracked gas transport in molten salts, ranging from small-scale systems to large-scale loops. Gas transportation in a molten salt, LiCl-KCl, has been studied using the shadowgraph technique. Sensors such as residual gas analysis, Raman spectroscopy, and laser-induced breakdown spectroscopy have been tested for off-gas measurements. These data were used to interpret how the gases move through the upstream salt / cover gas and the interface between them. Differential pressure measurements were able to detect individual gas bubbles as they popped at the liquid–gas interface. Salt aerosols were collected on a cascade impactor. Their formation was also observed directly via the shadowgraph method, and most of these aerosols were launched ballistically into the cover gas. Fine mists could also be observed. Convection currents through the salt were visualized and can be used to calculate the thermophysical properties of the salt itself. The apparatuses described in this report and in a previous work (McFarlane et al. 2023) have been commissioned and are available for use in making further measurements of salt/surrogate fission product behavior. Novel approaches using neutron imaging are planned for the study of fluoride salts, which cannot be contained in quartz, so shadowgraph visualization is not available. Bubble transport in convective flow and in a slow-moving salt column are planned. The mobile laser-induced breakdown spectroscopy (LIBS) system is available for several applications, including iodine capture in a molten hydroxide scrubber, H2 transport though molten salts to complement Raman analysis, and online tracking of salt aerosol generation, transport, and deposition. This report summarizes the findings from FY24 and the plans for FY25. The work completes milestone M2AT-24OR0702013 of the DOE-NE Advanced Reactor Technology, Molten Salt Reactor Campaign, DOE-NE-5.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Experimental measurement of the effective contact angle for solvent/packing interactions in a structured packed column for CO 2 capture

The contact angle is a critical factor for determining the effective mass transfer area for carbon dioxide (CO 2 ) capture via the chemical absorption process in a packed column, and thus the overall capture efficiency of the packed column. Many widely used commercial packings involve microscale features (perforation, corrugation, etc.) that may also affect the wetting behavior. This study proposes a systematic method of using a modified Wilhelmy plate to measure the effective contact angle to characterize the solvent and featured packing interaction. In lieu of computational efforts relying on assumptions guided by semi-informed correlation, the proposed method directly measures the effective contact angle as a function of the solvent and packing thermophysical and hydrodynamic properties. The characterization of the effective contact angle is then integrated in the computational fluid dynamics modeling to reduce uncertainty in the prediction of effective mas transfer area. Experiments were conducted for stainless steel coupons using water and aqueous sodium hydroxide (NaOH) solvent for verification of the proposed experimental protocol. The surface tension of the aqueous NaOH solvent was altered using surfactant and antifoam. The effective contact angle increases with the increased value of the surface tension for flat stainless steel sheets. On the other hand, effective contact angles do not vary in Mellapak coupons for aqueous monoethanolamine (MEA) and NaOH solvent. In this case, surface textures and sheet design play a dominant role in the surface tension of solvents. Furthermore, CO 2 loading has a significant effect on the contact angle for Mellapak coupons. As expected, the contact angle decreases with the increasing temperature of CO 2 capture solvents (MEA, EEMPA). The effective contact angle measurement using the Wilhelmy plate method can provide more accurate and efficient solvents for characterizing the solvent-packing surface interactions. Subsequently, the present method enhances the accuracy of the prediction of the effective mass transfer area in carbon capture by solvent absorption.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Experimental and Computational Heat Transfer Study of sCO2 Single Jet Impingement

The present study experimentally investigates the heat transfer capability of supercritical carbon dioxide (sCO2) single-jet impingement. The evaluated jet Reynolds number range is between 80,000 and 1,000,000, with a non-dimensional jet-to-target surface spacing of 2.8. CO2- impinging jet stagnation conditions were maintained at approximately 20 MPa and a temperature of 673 K for most experiments. The goal is to understand how changes in the aforementioned parameters influence heat transfer between the working fluid and the heated surface. Additionally, due to the elevated Reynolds numbers and difference in thermodynamic properties between air and CO2, air-derived impingement correlations may not be appropriate for CO2 impingement; these correlations will be evaluated against experimental sCO2 impingement data. At the time of this study, no sCO2 impingement data was available relevant to sCO2 power cycles. The target surface is a 1.5-inch diameter copper block centered on the 3mm jet orifice. A mica heating element bolted to the bottom of the copper block provides a uniform heat flux. Thermocouples embedded in the copper block are used to determine the surface temperature. Nusselt numbers obtained from experimental sCO2 data are compared to area-averaged Nusselt numbers from air-derived correlations. The comparisons showed that air correlations drastically underpredict the heat transfer when sCO2 is used as the working fluid. A modified sCO2 correlation using experimental data at discussed conditions is derived based on an existing air correlation.

Richardson, John↗

Characterizing Pore-Scale Geochemical Alterations in Eagle Ford and Barnett Shale from Exposure to Hydraulic Fracturing Fluid and CO 2 /H 2 O

As demand increases for an affordable energy source that is tied to an environmental obligation to reduce greenhouse gas emissions and water usage, there is a growing consideration in shale production utilizing processes such as 1) enhancing hydrocarbon recovery via carbon dioxide (CO 2 ) flooding, 2) using CO 2 as a fracturing agent to minimize water use, and 3) storing CO 2 in depleted shale formations to mitigate emissions to the atmosphere. Understanding the geochemical reactions and alterations that occur as shale is exposed to fluids and CO 2 is necessary to develop and optimize each of these processes for field applications. While the majority of shale formations are stimulated using traditional fracturing fluid, some may be fractured using CO 2 or other non-traditional means. We examine the effect fracturing fluid has on shale and how it behaves with secondary exposure to dry CO 2 or CO 2 -saturated water using in situ Fourier Transform infrared spectroscopy (FTIR), feature relocation scanning electron microscopy (SEM), and surface area and pore size analysis using volumetric gas sorption. These techniques were performed on Eagle Ford and Barnett shale samples that were exposed to fracturing fluid and unexposed (as received). Shales that have been exposed to traditional fracturing fluid experienced two reaction fronts. The first reaction front was formed during exposure to the fracturing fluid (pH of ~1.4). A secondary reaction front was formed as a result of CO 2 -saturated fluid exposure in the form of carbonic acid (pH ~5.6). These two different reaction mechanisms drove multiple dissolution and precipitation cycles which altered petrophysical properties of the shale and could lead to a significant impact on flow pathways. FTIR showed that equilibration of carbonate dissolution and precipitation cycles could take as long as 35 days. Samples exposed to fracturing fluid showed significantly less carbonate reactivity compared to those exposed to water. Pore size analysis results indicate exposure to fracturing fluid blocked small nanopores (10-15 nm) reducing BET surface area and total pore volume. SEM results show barite precipitated heavily during exposure to fracturing fluid. It appeared that carbonic acid was able to extract sulfur from organic matter to form gypsum evaporites. The mineralogical (barite precipitation and calcite dissolution/precipitation) and pore-scale alterations observed in these samples may lead to enhancement of flow pathways for injected CO 2 or produced hydrocarbons.

04 OIL SHALES AND TAR SANDS↗

Experimental and Computational Heat Transfer Study of sCO 2 Single-Jet Impingement

This study experimentally and computationally investigates the heat transfer capability of supercritical carbon dioxide (sCO 2 ) single jet impingement. The evaluated jet Reynolds number range is between 80,000 and 600,000, with a nondimensional jet-to-target surface spacing of 2.8. CO 2 -impinging jet stagnation conditions were maintained at approximately 200 bar and a temperature of ∼400 °C for most experiments. The goal is to understand how changes in the aforementioned parameters influence heat transfer between the working fluid and the heated surface. Additionally, due to the elevated Reynolds numbers and difference in thermodynamic properties between air and CO 2 , air-derived impingement correlations may not be appropriate for CO 2 impingement; these correlations will be evaluated against experimental sCO 2 impingement data. At the time of this study, no sCO 2 impingement data was available relevant to sCO 2 power cycles. The target surface is a 1.5-in. diameter copper block centered on the 3 mm jet orifice. A mica heating element bolted to the bottom of the copper block provides a uniform heat flux. Thermocouples embedded in the copper block are used to determine the surface temperature. Nusselt numbers obtained from experimental sCO 2 test data are compared to area-averaged Nusselt numbers from air-derived correlations. The comparisons showed that air correlations drastically underpredict the heat transfer when sCO 2 is used as the working fluid. A modified sCO 2 correlation using experimental data at discussed conditions is derived based on an existing air correlation. A CFD study is also performed to further investigate sCO 2 heat transfer characteristics, and assess the numerical model applicability to this problem type.

42 ENGINEERING↗