Search NASA⌕ Search

SEARCH · Search NASA

Results for “FLUID MECHANICS”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Thermomechanical Stress and Creep-Fatigue Analysis of a High-Temperature Prototype Receiver for Heating Particles

This work presents a three-dimensional (3D) thermomechanical model of a prototype-scale enclosed light trapping solar receiver for heating particles. Results of the thermoelastic model are used to estimate receiver lifetime under maximum flux conditions. A computational fluid dynamics (CFD) model is first developed to predict the temperature fields in a multi-panel assembly under steady operating conditions. Solar flux distributions on the receiver are obtained from the software package SolTrace and applied to the 3D thermal model. The subsequent particle heating is captured through a simplified 1D energy balance. Panel reradiation is considered through a surface-to-surface radiation model and natural convection loss to the surrounding air is captured in a representative fluid domain surrounding the receiver. The resulting temperature fields from the CFD analysis are used as inputs for a thermoelastic mechanical model with representative boundary conditions. With the resultant temperature and stress fields, a creep-fatigue damage and lifetime analysis is performed using the linear damage accumulation (LDA) theory. The Manson-Coffin formula and Larson Miller correlation are used to calculate the fatigue and creep, respectively. A maximum damage (corresponding to a 30-year service life) is defined for design assessment. The model was first developed and verified in detail by comparing with published results in the literature (temperature and stress profiles and distributions, and creep/fatigue damage fractions) for tubular solar receivers with supercritical carbon dioxide as the working fluid. It was then implemented to model a planar-cavity receiver with various design parameters. Specifically, three different design geometries are considered, and the results show that a maximum temperature of approximately 1200 K could be reached for each design with the given incident solar flux, with the main difference being the distribution of these temperatures. Preliminary resulting stresses for the small-scale prototype without design optimization vary from 20 MPa to 250 MPa for each design, with the maximum stresses occurring on the front face and concave geometry on the side of the panel. In future work, the developed methodology shown here will be applied to analyze a full-scale (50-150 MWth) receiver.

concentrated solar power↗
DOE OSTI2022

Carbon nanofibers (CNFs) dispersed in ultra-high performance concrete (UHPC): Mechanical property, workability and permeability investigation

This paper presents experimental data on the effect of carbon nano-fiber (CNF) dispersion on mechanical behavior, workability, and permeability of ultra-high-performance concrete containing CNF (UHPC-CNF). The addition of CNFs in concrete has been shown to provide improved performance. However, an inhomogeneous distribution of CNFs in the cement matrix can result in no impact or decreased mechanical performance. Due to the van der Waal's forces and hydrophobic surface properties of the CNFs, good dispersion in mix fluids and cement paste is challenging. Experimental studies of UHPC-CNF with different dispersion methods are provided. Further mixing with optimized UHPC paste allowed for the investigation of the optimal dosage of CNFs, impact on porosity, compressive strength, flexural strength, hydration rate, and slump. Scanning electron microscopy (SEM) analysis visually revealed the effect of dispersion. Our results indicate that shear mixing and subsequent ultrasonic dispersion with chemical surfactants can provide a well-dispersion liquid admixture of CNFs resulting in high mechanical performance of UHPC-CNFs composites. The water permeability and chloride resistance of optimized UHPC-CNFs composites were further evaluated with the wicking test and ponding test to reveal improved performance for these properties as well.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Computational Assessment of Hydrodynamic Loads on Rockeries for River Bank Protection

Colorado experienced an extreme flooding event in September of 2013 that ravaged canyon roadways southeast of Estes Park, CO. Central Federal Lands was asked to rebuild SH 34 and CO 43 after the floods. To minimize river encroachments in narrow canyons, roadway fill embankments must be over steepened. Using a rockery to stabilize the roadway embankment rather than traditional vertical, concrete structures was the solution. However, no analytical techniques were available for evaluating the hydraulic stability of such structures. Argonne National Laboratory was employed to estimate the hydrodynamic forces by using advanced, computational fluid dynamics (CFD). An investigation was conducted on a section of a model rockery comprised of large boulders exposed to a variety of hydraulic flow conditions. A repeatable geometric section of a typical rockery structure was modeled and used to create boundaries for the CFD model. The test case matrix covers computer simulations of various rockery and channel geometries. The geometrical parameters taken into account are: rockery height, outside face batter angle of the rockery, rockery base width, river bed width, and the hydraulic parameters are: flow depth, flow velocity, and hydraulic angle of attack. Two failure mechanisms of rockery structures, sliding and overturning, are considered. The existing design equations for both modes are modified by addition of hydrodynamic force terms. The new terms were calculated for a parametric set of geometry and condition variation using the CFD model. Moreover, formulas for dimensionless hydraulic force coefficients are proposed that can be used in the design of new rockeries or in evaluation of existing dry-stack rockeries in river environments.

42 ENGINEERING↗
DOE OSTI2020

Assessing mechanical response of CO 2 storage into a depleted carbonate reef using a site-scale geomechanical model calibrated with field tests and InSAR monitoring data

Geomechanical risks of injection have raised concerns regarding secure CO 2 storage. In this work, a combined monitoring and modeling approach is used to assess the stress changes and surface uplift associated with CO 2 injection into a depleted carbonate reef of the Michigan basin. A site-scale geomechanical model is built by assigning mechanical properties of formations using well-log and experimental data. Gravity load is applied to the model to estimate the vertical component of stress as well as different lateral boundary displacement scenarios to estimate horizontal stresses. We used a poroelastic pressure-dependent model (instead of a linear elastic mechanical earth model) to calibrate initial stresses using hydraulic fracture test data measured at depleted reservoir status. Multi-phase fluid flow-geomechanical simulations are performed to estimate the poroelastic response during (1) primary depletion (2) field-scale CO 2 injection phase (3) a hypothetical forecast scenario in which well bottom hole pressure (BHP) reach 45000 KPa. The predicted surface uplift is less than 1 mm at the end of the field-scale CO 2 injection phase which is in good agreement with Interferometric Synthetic Aperture Radar (InSAR) uplift measurement. Although the InSAR data shows an insignificant uplift, hydromechanical modeling of injection shows that CO 2 injection still causes reservoir deformation emphasizing the role of carbonate overburden and reservoir formation mechanical properties and limited size of reef on diminishing the surface deformation. Modeling indicates poroelastic response of caprock matters to estimate uplift. The lower permeability of the top two layers provides additional barrier to large uplift. Also, history of subsidence due to production should be accounted to predict uplift due to a follow up injection correctly. This report shows the significance of combining a calibrated geomechanical model with field measured stresses and monitoring data to be used as a tool to ensure the safety of CO 2 storage.

42 ENGINEERING↗
DOE OSTI2021

Fluid Dynamics in Curvilinear Coordinates without Fictitious Forces

The use of curvilinear coordinates is sometimes indicated by the inherent geometry of a fluid dynamics problem, but this introduces fictitious forces into the momentum equations that spoil the strict conservative form. If one is willing to work in three dimensions, these fictitious forces can be eliminated by solving for rectangular (Cartesian) momentum components on a curvilinear mesh. A thoroughly geometric approach to fluid dynamics on spacetime demonstrates this transparently, while also giving insight into a greater unity of the relativistic and nonrelativistic cases than is usually appreciated.

pseudo-forces↗
DOE OSTI2021

A coupled thermo-hydro-mechanical model for simulating leakoff-dominated hydraulic fracturing with application to geologic carbon storage

A potential risk of injecting CO2 into storage reservoirs with marginal permeability (≲ 10 mD (1 mD = 10 –15 m 2 )) is that commercial injection rates could induce fracturing of the reservoir and/or the caprock. Such fracturing is essentially fluid-driven fracturing in the leakoff-dominated regime. Recent studies suggested that fracturing, if contained within the lower portion of the caprock complex, could substantially improve the injectivity without compromising the overall seal integrity. Modeling this phenomenon entails complex coupled interactions among the fluids, the fracture, the reservoir, and the caprock. Here, we develop a simple method to capture all these interplays in high fidelity by sequentially coupling a hydraulic fracturing module with a coupled thermal-hydrological-mechanical (THM) model for nonisothermal multiphase flow. The model was made numerically tractable by taking advantage of self-stabilizing features of leakoff-dominated fracturing. The model is validated against the PKN solution in the leakoff-dominated regime. Moreover, we employ the model to study thermo-poromechanical responses of a fluid-driven fracture in a field-scale carbon storage reservoir that is loosely based on the In Salah project's Krechba reservoir. The model reveals complex yet intriguing behaviors of the reservoir-caprock-fluid system with fracturing induced by cold CO 2 injection. We also study the effects of the in situ stress contrast between the reservoir and caprock and thermal contraction on the vertical containment of the fracture. The proposed model proves effective in simulating practical problems on length and time scales relevant to geological carbon storage.

58 GEOSCIENCES↗
DOE OSTI2019

Induced seismicity in geologic carbon storage

Geologic carbon storage, as well as other geo-energy applications, such as geothermal energy, seasonal natural gas storage and subsurface energy storage imply fluid injection and/or extraction that causes changes in rock stress field and may induce (micro)seismicity. If felt, seismicity has a negative effect on public perception and may jeopardize wellbore stability and damage infrastructure. Thus, induced earthquakes should be minimized to successfully deploy geo-energies. However, numerous processes may trigger induced seismicity, which contribute to making it complex and translates into a limited forecast ability of current predictive models. We review the triggering mechanisms of induced seismicity. Specifically, we analyze (1) the impact of pore pressure evolution and the effect that properties of the injected fluid have on fracture and/or fault stability; (2) non-isothermal effects caused by the fact that the injected fluid usually reaches the injection formation at a lower temperature than that of the rock, inducing rock contraction, thermal stress reduction and stress redistribution around the cooled region; (3) local stress changes induced when low-permeability faults cross the injection formation, which may reduce their stability and eventually cause fault reactivation; (4) stress transfer caused by seismic or aseismic slip; and (5) geochemical effects, which may be especially relevant in carbonate-containing formations. We also review characterization techniques developed by the authors to reduce the uncertainty in rock properties and subsurface heterogeneity both for the screening of injection sites and for the operation of projects. Based on the review, we propose a methodology based on proper site characterization, monitoring and pressure management to minimize induced seismicity.

58 GEOSCIENCES↗
DOE OSTI2022

Preliminary Primary System Thermal Fluids Analysis of a Horizontal Compact HTGR

The Horizontal Compact High Temperature Gas Reactor (HC-HTGR) is being designed by a multi-disciplinary team of nuclear, mechanical, and structural engineers, under the support of a DOE-NE Advanced Reactor Demonstration Program’s Advanced Reactor Concepts-20 (ARC-20) award. The objective of this ARC-20 project is to deliver a conceptual design for the proposed HC-HTGR in 3 years and support its commercialization as a safe, low-cost HTGR. Argonne National Laboratory (Argonne) is collaborating on the thermal hydraulic design and analysis of the HC-HTGR reactor pressure vessel internals. The scope of this work includes ensuring the reactor is able to maintain maximum core temperatures below designated safety thresholds during normal operation, shutdown, and accident conditions. This report documents the preliminary thermal hydraulic analysis of the HC-HTGR core design performed with a 1D fluid-3D solid coupled model built using the System Analysis Module (SAM). This assembly level model was utilized to inform the core assembly design, predict the temperature distribution in the peak power assembly including the peak fuel temperature, coolant channel outlet temperatures, and graphite temperature gradients. A key result of this analysis was the determination that the peak fuel temperature remains below the safety threshold of 1250°C. Additionally, this model was used to assess the assembly coolant channel and bypass flow mass flow rate distribution. Following the assembly level analysis, attention turned to the development of a full core reduced order model to be used to predict the core wide coolant flow distribution and to model certain operational and accidental transients.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2018

Exploration of thermal counterflow in He II using particle tracking velocimetry

Flow visualization using particle image velocimetry (PIV) and particularly particle tracking velocimetry (PTV) has been applied to thermal counterflow in He II for nearly two decades now, but the results remain difficult to interpret because tracer particle motion can be influenced by both the normal fluid and superfluid components of He II as well as the quantized vortex tangle. For instance, in one early experiment it was observed (using PTV) that tracer particles move at the normal fluid velocity v n , while in another it was observed (using PIV) that particles move at v n /2. Besides the different visualization methods, the range of applied heat flux investigated by these experiments differed by an order of magnitude. To resolve this apparent discrepancy and explore the statistics of particle motion in thermal counterflow, we apply the PTV method to a wide range of heat flux at a number of different fluid temperatures. In our analysis, we introduce a scheme for analyzing the velocity of particles presumably moving with the normal fluid separately from those presumably influenced by the quantized vortex tangle. Our results show that for lower heat flux there are two distinct peaks in the streamwise particle velocity probability density function (PDF), with one centered at the normal fluid velocity v n (named G2 for convenience) while the other is centered near v n /2 (G1). For higher heat flux there is a single peak centered near v n /2 (G3). Using our separation scheme, we show quantitatively that there is no size difference between the particles contributing to G1 and G2. We also show that nonclassical features of the transverse particle velocity PDF arise entirely from G1, while the corresponding PDF for G2 exhibits the classical Gaussian form. The G2 transverse velocity fluctuation, backed up by second sound attenuation in decaying counterflow, suggests that large-scale turbulence in the normal fluid is absent from the two-peak region. We offer a brief discussion of the physical mechanisms that may be responsible for our observations, revealing that G1 velocity fluctuations may be linked to fluctuations of quantized vortex line velocity, and suggest a number of numerical simulations that may reveal the underlying physics in detail.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

A Deterministic and Probabilistic Framework for Forecasting of Time-Series Damage States and Associated End-of-lIfe of a Pressurizer Water Reactor Surge Line under Design-Basis and Grid-Load-Following Loading Conditions

Argonne National Laboratory in support of the Department of Energy’s Light Water Reactor Sustainability (LWRS) program. The present semi-annual report covers the period between May 2018 and September 2018. In this report, we present the work performed to further improve the capability for structural integrity prediction. Part of this work focused on validating a computational fluid dynamics model for modeling thermal stratification. We also present finite element (FE) model based thermal-mechanical stress analysis of a pressurized water reactor (PWR) surge line under grid-load-following loading condition. In addition, we present work related to end-of-life prediction of a PWR surge line based on ASME code and NUREG-6909 deterministic approaches under design-basis and grid-load-following conditions. Our research in this report period also majorly involved with the development of a strategy for probabilistic fatigue life estimation. First, we applied a Weibull probabilistic modeling approach based on end-of-life fatigue data (traditional S~N data under in-air and PWR water conditions) to estimate the probabilistic life for a given strain amplitude. Second, we developed a data analytics based Markov-Chain-Monte-Carlo (MCMC) model for probabilistic modeling of time-series fatigue damage states and the associated probabilistic end-of-life of a component. Unlike the Weibull type approach, the MCMC type may not require a large number of fatigue tests for estimating the end-of-life for given loading and environmental conditions. We demonstrated the use of this framework through symmetric fatigue tests (stress- or strain-ontrolled tests either under in-air or PWR-water conditions) obtained from our earlier experiments. The MCMC model can use base time-series data either obtained from a model or fatigue experiment. Since experimental data are typically more accurate than model-calculated data, we used experiment-based time-series data to demonstrate the MCMC framework. However, when performing a complex experiment is not possible (due to time and costs involved), a model-based approach can be used to generate the time-series data. Similar model-based data are presented in our earlier report (e.g. for cyclic hardening damage states of a PWR surge line under constant amplitude and design-basis loading [1]). However, for the MCMC work in this report, we used experiment-based time-series damage states obtained from fatigue experiments under design-basis and grid-load-following loading conditions. The experimental loading inputs were based on thermal-mechanical FE model results of a typical PWR surge line. Then, we used the experimental time-series damage state data for forecasting the probabilistic time-series damage states and end-of-life of the PWR surge line. Based on this MCMC modeling strategy, and with the assumed loading, environment and thermal-mechanical boundary conditions, we have estimated that a 316SS-pure base metal PWR surge line would survive at least 159 years with zero failure probability.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

A New Modeling Framework for Multi-Scale Simulation of Hydraulic Fracturing and Production from Unconventional Reservoirs

This paper describes a new modeling framework for microscopic to reservoir-scale simulations of hydraulic fracturing and production. The approach builds upon a fusion of two existing high-performance simulators for reservoir-scale behavior: the GEOS code for hydromechanical evolution during stimulation and the TOUGH+ code for multi-phase flow during production. The reservoir-scale simulations are informed by experimental and modeling studies at the laboratory scale to incorporate important micro-scale mechanical processes and chemical reactions occurring within the fractures, the shale matrix, and at the fracture-fluid interfaces. These processes include, among others, changes in stimulated fracture permeability as a result of proppant behavior rearrangement or embedment, or mineral scale precipitation within pores and microfractures, at µm to cm scales. In our new modeling framework, such micro-scale testing and modeling provides upscaled hydromechanical parameters for the reservoir scale models. We are currently testing the new modeling framework using field data and core samples from the Hydraulic Fracturing Field Test (HFTS), a recent field-based joint research experiment with intense monitoring of hydraulic fracturing and shale production in the Wolfcamp Formation in the Permian Basin (USA). Below, we present our approach coupling the reservoir simulators GEOS and TOUGH+ informed by upscaled parameters from micro-scale experiments and modeling. We provide a brief overview of the HFTS and the available field data, and then discuss the ongoing application of our new workflow to the HFTS data set.

58 GEOSCIENCES↗
DOE OSTI2024

Fluid Mixing during Phase Separation in Crystallizing White Dwarfs

Abstract Accurate models of cooling white dwarfs must treat the energy released as their cores crystallize. This phase transition slows the cooling by releasing latent heat and also gravitational energy, which results from phase separation: liquid C is released from the solid C/O core, driving an outward carbon flux. The Gaia color–magnitude diagram provides striking confirmation of this theory by revealing a mass-dependent overdensity of white dwarfs, indicating slowed cooling at the expected location. However, the observed overdensity is enhanced relative to the models. Additionally, it is associated with increased magnetism, suggesting a link between crystallization and magnetic field generation. Recent works aimed at explaining an enhanced cooling delay and magnetic field generation employ a uniform mixing prescription that assumes large-scale turbulent motions; we show here that these calculations are not self-consistent. We also show that thermohaline mixing is most likely efficient enough to provide the required chemical redistribution during C/O phase separation, and that the resulting velocities and mixing lengths are much smaller than previous estimates. These reduced fluid motions cannot generate measurable magnetic fields, suggesting any link with crystallization needs to invoke a separate mechanism. Finally, this mixing alters the chemical profiles, which in turn affects the frequencies of the pulsation modes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Insights into the hydro-mechanical behavior of a decimeter-scale fracture using the mini-SIMFIP probe

ABSTRACT: Understanding hydro-mechanical couplings in fractured rocks is essential for predicting the rock mass response during high-pressure fluid injection, including the stimulation of enhanced geothermal systems. However, fluid-driven fracture dislocations are challenging to measure due to the need for local displacement data at high fluid pressures. In this study, fluid-driven displacement across a decimeter-scale laboratory fracture was investigated using the mini-SIMFIP (step rate injection method for fracture in-situ properties) probe, which is a smaller version of the SIMFIP tool (Guglielmi et al., 2014). The mini-SIMFIP probe is able to resolve the full 3D displacement vector of a fracture. The probe was installed in one of two boreholes across a decimeter-scale saw-cut granite fracture. Two pressure step injection tests were conducted under the same isotropic stress conditions. By varying injection between the boreholes, we estimated the aperture profile across the fracture. The results show a consistent pressure-dependent opening as long as steady-state flow was maintained. At a certain pressure step, steady-state conditions were no longer achievable and the fracture opened rapidly. The pressure-opening relationship diverged between the two tests and indicated a homogenization of the aperture profile and an increase in the non-linearity of the flow regime.

Osten, J↗
DOE OSTI2022

Control rod drive mechanism with heat pipe cooling

A representative cooling system for a nuclear reactor control rod drive mechanism (CRDM) includes an evaporation section located within or next to the CRDM and a condensation section fluidly coupled to the evaporation section. The cooling system includes a set of heat fins coupled to drive coils in the CRDM and heat pipes that extend through the drive coils and heat fins. A fluid evaporates while in the evaporation section of the heat pipes from heat generated by the CRDM and moves out of the evaporation section into the condensation section in the heat fins. The fluid cools and condensates while in the condensation section, recirculating back into the evaporation section. This passive natural circulation cooling system reduces or eliminates the number of water hoses, piping, and other water pumping equipment typically used for cooling a CRDM thereby increasing nuclear reactor reliability and simplifying nuclear reactor operation and maintenance.

Lobscheid, Christian↗
DOE OSTI2021

A weighted Shifted Boundary Method for free surface flow problems

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods and was recently introduced for the Poisson, linear advection/diffusion, Stokes, Navier-Stokes, acoustics, and shallow-water equations. By reformulating the original boundary value problem over a surrogate (approximate) computational domain, the SBM avoids integration over cut cells and the associated problematic issues regarding numerical stability and matrix conditioning. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions. Hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we extend the SBM to the simulation of incompressible Navier-Stokes flows with moving free-surfaces, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach prevents spurious pressure oscillations in time, which would otherwise be produced if the total active fluid volume were to change abruptly over a time step. In fact, the proposed weighted SBM method induces small mass (i.e., volume) conservation errors, which converge quadratically in the case of piecewise-linear finite element interpolations, as the grid is refined. Finally, we present an extensive set of two- and three-dimensional tests to demonstrate the robustness and accuracy of the method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Statistics of Point Vortex Turbulence in Non-neutral Flows and in Flows with Translational and Rotational Symmetries

A theory (Esler and Ashbee in J Fluid Mech 779:275–308, 2015) describing the statistics of N freely-evolving point vortices in a bounded two-dimensional domain is extended. First, the case of a non-neutral vortex gas is addressed, and it is shown that the density of states function can be identified with the probability density function of an infinite sum of independent non-central chi-squared random variables, the details of which depend only on the shape of the domain. Equations for the equilibrium energy spectrum and other statistical quantities follow, the validity of which are verified against direct numerical simulations of the equations of motion. Second, domains with additional conserved quantities associated with a symmetry (e.g., circle, periodic channel) are investigated, and it is shown that the treatment of the non-neutral case can be modified to account for the additional constraint.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Waste heat recovery system

A waste heat recovery system includes a Rankine cycle (RC) circuit having a pump, a boiler, an energy converter, and a condenser fluidly coupled via conduits in that order, to provide additional work. The additional work is fed to an input of a gearbox assembly including a capacity for oil by mechanically coupling to the energy converter to a gear assembly. An interface is positioned between the RC circuit and the gearbox assembly to partially restrict movement of oil present in the gear assembly into the RC circuit and partially restrict movement of working fluid present in the RC circuit into the gear assembly. An oil return line is fluidly connected to at least one of the conduits fluidly coupling the RC components to one another and is operable to return to the gear assembly oil that has moved across the interface from the gear assembly to the RC circuit.

30 DIRECT ENERGY CONVERSION↗