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 217 records · Page 12

Anti-symmetric representation of the extended magnetohydrodynamic equations

We introduce the anti-symmetric representation of the extended magnetohydrodynamic (MHD) equations. In this representation, the use of the anti-symmetric flux operator (∇·ν+ν·∇) results in conservation theorems with discrete analogs. Inherently robust numerical applications are achieved with little effort, and conservation to machine precision is possible with simple numerical schemes. Starting from the two-fluid equations, we construct a single-fluid MHD model based on generalized center-of-mass variables for the mass (√ρ), momentum (√ρν), and pressure (√ p ). This model is shown to possess identical conservation properties to the two-fluid system, with the only restriction being the use of a single temperature. Common approximations to the Braginskii heat fluxes and to the gyroviscous stress tensor are cast into our representation for convenience. The discrete conservation properties are verified using the classic Orszag–Tang vortex problem. In addition to the favorable mass, momentum, and energy conservation properties, the time reversibility of the simulations is demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Magneto-structural and induction heating properties of MFe2O4 (M?=?Co, Mn, Zn) MNPs for magnetic particle hyperthermia application

The thermal decomposition approach was used for the synthesis of MFe2O4 magnetic nanoparticles (MNPs) by substituting M as Co, Mn, and Zn. The obtained MNPs were characterized for magneto-structural properties using X-Ray diffraction patterns, FTIR, Raman and Mossbauer spectroscopy techniques which confirm the formation of phase pure cubic spinel ferrite with space group Fd3m and five Raman active modes. The size, morphology, and compositional analysis was performed using HRTEM and EDX where the size of MNPs was found to be less than 10 nm that attains superparamagnetism. The magnetic hyperthermia performance of obtained MNPs was evaluated by induction heating experiments at magnetic field range 13.3 to 26.7 kAm-1. The specific absorption rate (SAR) and intrinsic loss power (ILP) values were determined at different magnetic fields within the human tolerable range and correlated with magneto-structural properties to evaluate its potential for possible magnetic particle hyperthermia therapy.

Magnetic nanoparticles, complex-decomposition, mag↗

Evolution of Pore Structure in Organic-Lean and Organic-Rich Mudrocks

The pore structure of mudrocks is a key characteristic to evaluate flow behavior through these rocks. Although significant advances have been made to resolve pore characteristics, porosity evolution, or pore connectivity, there is still insufficient knowledge linking porosity evolution to flow and transport in mudrocks. To better understand these links, we conducted very small-angle (VSANS) and small-angle neutron scattering (SANS) experiments on 13 sets of mudrocks from global locations, characterized by differences in composition, maturity, and depositional environment. Our results indicate that a homogeneous pore structure reflects a stable and low-energy depositional environment. Organic-lean mudrocks mainly contain pores <10 nm, which restricts flow and is beneficial for, e.g., caprock seal integrity. Increasing maturity as a result of increasing burial depth, and therefore, temperature leads to porosity evolution. In organic-rich mudrocks (ORMs), this happens through the generation of organic matter pores and pore preservation around diagenetically formed carbonates. ORMs mainly contain macropores (>50 nm), favoring production from, e.g., unconventional reservoirs. SANS-derived pore size distributions will help to better understand fluid imbibition and flow properties in mudrocks through comprehensive quantitative characterization of the nano- to micron scale pore network conduits and their relationships to burial diagenesis.

58 GEOSCIENCES↗

Nanoscale Mg-Depleted Layers Slow Carbonation of Forsterite (Mg 2 SiO 4 ) When Water Is Limited

Passivation of silicate surfaces by accumulated reaction products is an obstacle to efficient CO 2 mineralization. In this study, we investigate a unique passivation effect during the carbonation of the basalt mineral forsterite (Mg 2 SiO 4 ) in humid supercritical CO 2 (50 °C, 90 bar). Using in situ high-pressure infrared spectroscopy, we demonstrate that dissolution of forsterite into a thin water film slows significantly after reaction for ˜24 h, even under far-from-equilibrium conditions. Further, 29 Si magic angle spinning nuclear magnetic resonance spectroscopy detects a highly polymerized amorphous silica at this stage. On the basis of transmission electron microscopy and energy dispersive X-ray spectroscopy, we show that the silica is present as a Mg-depleted layer that is just 2–3 nm thick on the reacted forsterite particles. The decrease in the level of forsterite dissolution in the presence of an extraordinarily thin Mg-depleted layer can be strongly linked to properties of the thin fluid film at the surface, highlighting the importance of water during mineral carbonation. This study furthers our understanding of silicate mineral carbonation under select low-water, humidified fluid conditions relevant to basaltic geologic reservoirs, recovery of critical elements by carbonation of mafic ores, and sequestration of atmospheric CO 2 by enhanced rock weathering.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DC Hall coefficient of the strongly correlated Hubbard model

The Hall coefficient is related to the effective carrier density and Fermi surface topology in non-interacting and weakly interacting systems. In strongly correlated systems, the relation between the Hall coefficient and single-particle properties is less clear. Clarifying this relation would give insight into the nature of transport in strongly correlated materials that lack well-formed quasiparticles. In this work, we investigate the DC Hall coefficient of the Hubbard model using determinant quantum Monte Carlo in conjunction with a recently developed expansion of magneto-transport coefficients in terms of thermodynamic susceptibilities. At leading order in the expansion, we observe a change of sign in the Hall coefficient as a function of temperature and interaction strength, which we relate to a change in the topology of the apparent Fermi surface. We also combine our Hall coefficient results with optical conductivity values to evaluate the Hall angle, as well as effective mobility and effective mass based on Drude theory of metals.

36 MATERIALS SCIENCE↗

Low energy electrodynamics of CrI 3 layered ferromagnet

We report on the optical properties from terahertz (THz) to Near-Infrared (NIR) of the layered magnetic compound CrI 3 at various temperatures, both in the paramagnetic and ferromagnetic phase. In the NIR spectral range, we observe an insulating electronic gap around 1.1 eV which strongly hardens with decreasing temperature. The blue shift observed represents a record in insulating materials and it is a fingerprint of a strong electron-phonon interaction. Moreover, a further gap hardening is observed below the Curie temperature, indicating the establishment of an effective interaction between electrons and magnetic degrees of freedom in the ferromagnetic phase. Similar interactions are confirmed by the disappearance of some phonon modes in the same phase, as expected from a spin-lattice interaction theory. Therefore, the optical properties of CrI 3 reveal a complex interaction among electronic, phononic and magnetic degrees of freedom, opening many possibilities for its use in 2-Dimensional heterostructures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Fracture Sustainability in Enhanced Geothermal Systems: Experimental and Modeling Constraints

Enhanced geothermal systems (EGS) offer the potential for a much larger energy source than conventional hydrothermal systems. Hot, low-permeability rocks are prevalent at depth around the world, but the challenge of extracting thermal energy depends on the ability to create and sustain open fracture networks. Laboratory experiments were conducted using a suite of selected rock cores (granite, metasediment, rhyolite ash-flow tuff, and silicified rhyolitic tuff) at relevant pressures (uniaxial loading up to 20.7 MPa and fluid pressures up to 10.3 MPa) and temperatures (150–250 °C) to evaluate the potential impacts of circulating fluids through fractured rock by monitoring changes in fracture aperture, mineralogy, permeability, and fluid chemistry. Because a fluid in disequilibrium with the rocks (deionized water) was used for these experiments, there was net dissolution of the rock sample: this increased with increasing temperature and experiment duration. Thermal-hydrological-mechanical-chemical (THMC) modeling simulations were performed for the rhyolite ash-flow tuff experiment to test the ability to predict the observed changes. These simulations were performed in two steps: a thermal-hydrological-mechanical (THM) simulation to evaluate the effects of compression of the fracture, and a thermal-hydrological-chemical (THC) simulation to evaluate the effects of hydrothermal reactions on the fracture mineralogy, porosity, and permeability. Furthermore, these experiments and simulations point out how differences in rock mineralogy, fluid chemistry, and geomechanical properties influence how long asperity-propped fracture apertures may be sustained. Such core-scale experiments and simulations can be used to predict EGS reservoir behavior on the field scale.

15 GEOTHERMAL ENERGY↗

Model-predictive optimal control of ferrofluidic microrobots in three-dimensional space

Ferrofluid microrobots have emerged as promising tools for minimally invasive medical procedures. Their unique properties to navigate complex fluids and reach otherwise inaccessible regions of the human body have enabled new applications in targeted drug delivery, tissue engineering, and diagnostics. Here, this paper proposes a model-predictive controller for the external magnetic manipulation of ferrofluid microrobots in three dimensions (3D). The internal optimization routine of the controller determines appropriate changes in the applied electromagnetic field to minimize the deviation between the actual and desired trajectories of the microrobot. A linear system governing locomotion is derived and used as the equality constraints of the optimization problems associated with the feedback index. In addition to ferrofluid droplets, the controller presented in this work may be applied to other magnetically-pulled microrobots. Several experiments are performed to validate the controller and showcase its ability to adapt to changes in system parameters such as the desired tracking trajectory and the size, orientation, deformation, and velocity of the microrobot. The accuracy of the controller is analyzed for each experiment, and the average error is found to be within 0.25 mm for small velocities. An additional experiment is performed to demonstrate significant improvement over a PID controller that is optimally tuned using Bayesian optimization. The results presented in this paper suggest that the proposed control algorithm could enable new microrobotic capabilities in minimally invasive medical procedures, lab-on-a-chip applications, and microfluidics.

60 APPLIED LIFE SCIENCES↗

FY20 SAM Code Developments and Validations for Transient Safety Analysis of Advanced non-LWRs

The System Analysis Module (SAM) is under development at Argonne National Laboratory as a modern system-level modeling and simulation tool for advanced non-light water reactor safety analyses. It utilizes the object-oriented application framework MOOSE to leverage the modern software environment and advanced numerical methods. The capabilities of SAM are being extended to enable the transient modeling, analysis, and design of various advanced nuclear reactor systems. This report summarizes major progress in SAM code development, capability enhancements, demonstration, and validation to support transient safety analysis of advanced non-LWRs. Rapid developments continued in fiscal year 2020 (FY20) to support various needs of the advanced reactor community, especially the NRC and industry on the licensing safety analysis of advanced reactor designs. Significant code changes were made to provide various capability enhancements, bug fixes, and user friendliness improvements. Major code updates are summarized in Section 1, while four important enhancements are detailed in Sections 2-5, including a multi-dimension flow model; reactivity feedback and decay heat models; control and trip system modeling, and additional fluid and solid thermophysical property models. Code validation activities in FY20 include using test data from the Fast Flux Test Facility (FFTF), the High Temperature Test Facility (HTTF), and several separate effects test facilities for pebble-bed modeling.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

X-ray Computed Tomography on UNESE Core: FY2020 Data Report to Support Fracture and Multiphase Fluid Flow Studies

Natural and induced fractures are potential preferential pathways for migration of radioactive gases to earths surface from underground nuclear explosions (UNEs). This report documents X-ray computed tomography (XRCT) imaging on 26 samples of rock core that was collected to support the Underground Nuclear Explosion Signatures Experiment (UNESE) program. The XRCT datasets are intended to help fill a data gap on the three-dimensional (3D) characteristics of natural and/or induced fractures at the centimeter and smaller scale, which may strongly influence multiphase fluid flow and transport properties of preferential flow paths and interaction with the matrix of the surrounding host rock. Pre- and post-UNE rock samples were carefully chosen to enable comparison of fractures as a function of lithologic and petrophysical properties, as well as distance to the past UNEs. This report serves as documentation for the data, including an introduction with the research motivation, a methods and materials section, descriptions of the XRCT datasets without post-processing, and recommendations for 3D quantification via image analysis and digital rock physics.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Numerical Investigation of Wall-Cooling Effect on Aero-Optical Distortions for Hypersonic Boundary Layer

Compressible wall-modeled large-eddy simulations of Mach 8 turbulent boundary-layer flows over a flat plate were carried out for the conditions of the hypersonic wind tunnel at Sandia National Laboratories. The simulations provide new insight into the effect of wall cooling on the aero-optical path distortions for hypersonic turbulent boundary-layer flows. Four different wall-to-recovery temperature ratios, 0.3, 0.48, 0.71, and 0.89, are considered. Despite the much lower grid resolution, the mean velocity, temperature, and resolved Reynolds stress profiles from the simulation for a temperature ratio of 0.48 are in good agreement with those from a reference direct numerical simulation. The normalized root-mean-square optical path difference obtained from the present simulations is compared with that from reference direct numerical simulations, Sandia experiments, as well as predictions obtained with a semi-analytical model by Notre Dame University. Here the present analysis focuses on the effect of wall cooling on the wall-normal density correlations, on key underlying assumptions of the aforementioned model such as the strong Reynolds analogy, and on the elevation angle effect on the optical path difference. Wall cooling is found to increase the velocity fluctuations and decrease the density fluctuations, resulting in an overall reduction of the normalized optical path distortion. Compared to the simulations, the basic strong Reynolds analogy overpredicts the temperature fluctuations for cooled walls. Also different from the strong Reynolds analogy, the velocity and temperature fluctuations are not perfectly anticorrelated. Finally, as the wall temperature is raised, the density correlation length, away from the wall but inside the boundary layer, increases significantly for beam paths tilted in the downstream direction.

42 ENGINEERING↗

Temperature Measurements in Hypersonic Wind Tunnels via Femtosecond Coherent Anti-Stokes Raman Scattering

A femtosecond coherent anti-Stokes Raman scattering (fs CARS) instrument is developed to perform gas-phase thermometry in cold-flow hypersonic wind tunnels. Measurements are reported for Mach 8 and 14 pure-nitrogen flows. The fs CARS instrument includes a 100 fs pump/Stokes pulse and a spectrally narrow probe pulse from a second harmonic bandwidth compressor. Important experimental considerations such as limits on the pump/Stokes pulse energy are discussed. The fs CARS focusing and collimating optics are mounted on a two-axis translation stage system to scan the measurement location during a 30 second wind tunnel run. Single-laser-shot rotational CARS spectra are recorded at the laser repetition rate of 1 kHz in the wind tunnel freestream and near simple cone models. Spectral fitting is used to determine quantitative gas temperatures. Freestream temperatures at Mach 8 and 14 spanned ranges of 40–75 and 35–50 K, respectively, depending on tunnel operating conditions. Temperature variations across the central 100 mm span of the wind tunnel were quantified. Measured temperature jumps across conical bow shocks from various models varied by less than 1% from predicted values. Hypersonic boundary layer measurements were demonstrated. In conclusion, these measurements illustrate the utility and robustness of this instrument for the study of complex fluid flow phenomena in challenging ground test facilities.

Aerodynamics↗

Blade and Rim Seal Design of a First Stage High Pressure Turbine for a 300 MWe Supercritical CO2 Power Cycle

A first stage high-pressure turbine (HPT) blade is optimized for a 300 MWe supercritical CO2 (sCO2) power cycle using the surrogate-assisted genetic algorithm optimizer in Numeca FINE/Design 3D with objectives of increasing efficiency and decreasing heat load to the blade. The National Institute of Standards and Technology Reference Fluid Thermodynamic and Transport Properties Database (NIST REFPROP) [1] data for supercritical CO2 is formatted into tables of bicubic polynomial coefficients for use in condensable gas simulations in FINE/Turbo. Nearly 3000 unique shapes are evaluated via three-dimensional Reynolds Averaged Navier Stokes simulations, yielding increases in efficiency of up to 0.85 percentage points and decreases in heat load of 14%. A final blade, deemed the advanced blade, is chosen for future experimental analysis. Following this, a squealer tip optimization is performed on both the baseline and advanced blade designs. This optimization resulted in a performance gain of 1.25 points in efficiency and 15% reduction in tip heat load compared to the baseline flat tip design at the same clearance. In tandem, an optimization of the rotor-stator platform rim seal is performed using a parametrized geometry allowing for straight, meandering, and knife seal cavities. This multi-objective optimization focuses on decreasing the cooling mass flow and increasing the heat flux from the rotor and stator disk. The optimization resulted in cooling mass flow decreases of up to 26% while maintaining the average heat flux on the rim seal.

Tuite, Logan↗

Asymptotic preserving methods for fluid electron-fluid models in the large magnetic field limit with mathematically guaranteed properties (Final Report)

The current manuscript is a final report on the activities carried out under the Project LDRD-CIS #226834. In scientific terms, the work reported in this manuscript is a continuation of the efforts started with Project LDRD-express #223796 with final report of activities SAND2021-11481, see [83]. In this section we briefly explain what pre-existing developments motivated the current body of work and provide an overview of the activities developed with the funds provided. The overarching goal of the current project LDRD-CIS #226834 and the previous project LDRD-express #223796 is the development of numerical methods with mathematically guaranteed properties in order to solve the Euler-Maxwell system of plasma physics and generalizations thereof. Even though Project #223796 laid out general foundations of space and time discretization of Euler-Maxwell system, overall, it was focused on the development of numerical schemes for purely electrostatic fluid-plasma models. In particular, the project developed a family of schemes with mathematically guaranteed robustness in order to solve the Euler-Poisson model. This model is an asymptotic limit where only electrostatic response of the plasma is considered. Its primary feature is the presence of a non-local force, the electrostatic force, which introduces effects with infinite speed propagation into the problem. Even though instantaneous propagation of perturbations may be considered nonphysical, there are plenty of physical regimes of technical interest where such an approximation is perfectly valid.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fracture-matrix fluid exchange in oil-bearing unconventional mudstones

The poromechanical properties of unconventional reservoir materials are in large part dictated by their mineralogy. Since these properties govern the response to stress experienced during hydraulic fracturing, fluid production, and fluid injection, they play a central role in the formation of microcracks or bedding delaminations which ultimately dominate mass transport. In this work we study access to the porosity of end member unconventional reservoir materials, where the end members are predominantly dictated by carbonate content. Access to the porosity is quantified using state of the art 3D x-ray computed tomography coupled with physics informed data analytics. Xenon gas, which attenuates x-rays, provides a spatiotemporal map of access to the porosity. The accessible porosity is quantified over a range of net confining stress relevant to the manmade disturbances listed above. These experiments demonstrate that heavily carbonated mudstones are nearly impermeable at the core (~ cm) scale, while carbonate free analogues afford better access to the microstructure. Consistent with previous qualitative 2D radiographs, access to the interior of the clastic mudstones is first observed along planar microcracks, followed by slow penetration into the surrounding matrix. Physics informed data analytics of the 3D tomography measurements presented here show that these microcracks do not permit uniform access to the adjacent rock matrix. In addition, variation of the effective pressure elucidates the mechanisms that govern fracture/matrix fluid exchange. Under conditions consistent with hydrocarbon production fluid accumulates in the immediate vicinity of the nearest microcrack. While there is clear evidence that, as intended, part of this accumulation is from the more distant matrix, fluid is also squeezed out of the microcrack. The fluid build-up at the microcrack indicates that migration out of the rock is hindered by the coupled poroelastic response of the microcrack and adjacent rock matrix. We show that these mechanisms ultimately account for the meager oil recovery factors realized in practice. These insights have implications for making reservoir scale predictions based on core scale observations, and provide a basis for devising new asset development techniques to access more porosity, and enhance fluid extraction. Finally, these findings shed light on key features and mechanisms that govern shale storage capacity, with relevance to other important industrial processes, such as geologic CO 2 storage.

58 GEOSCIENCES↗

Uncoupling Electrokinetic Flow Solutions

The continuum-scale electrokinetic porous-media flow and excess charge redistribution equations are uncoupled using eigenvalue decomposition. The uncoupling results in a pair of independent diffusion equations for “intermediate” potentials subject to modified material properties and boundary conditions. The fluid pressure and electrostatic potential are then found by recombining the solutions to the two intermediate uncoupled problems in a matrix-vector multiplication. Expressions for the material properties or source terms in the intermediate uncoupled problem may require extended precision or careful rewriting to avoid numerical cancellation, but the solutions themselves can typically be computed in double precision. The approach works with analytical or gridded numerical solutions and is illustrated through two examples. The solution for flow to a pumping well is manipulated to predict streaming potential and electroosmosis, and a periodic one-dimensional analytical solution is derived and used to predict electroosmosis and streaming potential in a laboratory flow cell subjected to low frequency alternating current and pressure excitation. The examples illustrate the utility of the eigenvalue decoupling approach, repurposing existing analytical solutions or numerical models and leveraging solutions that are simpler to derive for coupled physics.

58 GEOSCIENCES↗

Fluid-driven fractures in granular media: Insights from numerical investigations

We investigate the mechanisms of opening-mode fracture initiation in granular media. The study is based on a simulation of grain-scale fluid-grain interactions through a coupled numerical approach in which the discrete element method is used to solve for the mechanics of a solid granular medium, and computational fluid dynamics is used to model fluid flow and drag forces. In this work, we present benchmark problems with analytical solutions and validate this numerical model against experiments on a viscous-drag-driven cavity in the literature. Additional simulation results show fracture initiation mechanisms in a random granular packing subjected to constant boundary stresses and to fluid injection with a localized source. The dimensionless variable $F_s/F_{sk}$ (ratio of seepage force $F_s$ and skeletal force $F_{sk}$) incorporates the impacts of physical properties and injection parameters including fluid viscosity, injection velocity, grain size, and effective stresses, and it has been used as a criterion separating regimes of fluid invasion and drag-driven fracture opening. Our simulation results show that $F_s/F_{sk}$ in combination with $τ_1$ (ratio of diffusion time from hydromechanical coupling and injection time) serves as a prediction of fracture opening within granular packing. We suggest a simple criterion ($F_s/F_{sk}$ > 1 or $τ_1$ > 0.17) that is valid for various types of granular media and injection conditions to determine if fracture opening will occur. Among other applications, this study is useful to predict the initiation and propagation of fractures in natural sediments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Bubble Point Measurements of Three Binary Mixtures of Refrigerants: R-32/1234yf, R-32/1234ze(E), and R-1132a/1234yf

In this article, bubble point pressures of three binary mixtures at two compositions each have been measured utilizing a static method. The mixtures studied were difluoromethane + 2,3,3,3-tetrafluoroprop-1-ene (R-32/1234yf), difluoromethane + (E)-1,3,3,3-tetrafluoroprop-1-ene (R-32/1234ze(E)), and 1,1-difluoroethene + 2,3,3,3-tetrafluoroprop-1-ene (R-1132a/1234yf). For all of the mixtures, the minimum temperature measured was 270 K. The maximum temperatures measured for the R-32/1234ze(E), R-32/1234yf, and R-1132a/1234yf mixtures were 360, 345, and 325 K, respectively. A total of 100 bubble point pressures are reported and compared to available literature data relative to the REFPROP database. Lemmon, E. W. [ NIST Standard Reference Database 23: Reference Fluid Thermodynamic and Transport Properties-REFPROP, 10.0; National Institute of Standards and Technology, Standard Reference Data Program: Gaithersburg, 2018]. Additionally, the bubble point pressures obtained in this study are modeled using the perturbed-chain statistical associating fluid theory (PC-SAFT) equation of state. The PC-SAFT model parameters, as well as deviations of this work and literature data from the model, are reported.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗