Search NASA⌕ Search

SEARCH · Search NASA

Results for “FLUID FLOW”

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 1,801 records · Page 100
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

Two-phase choked flow of cryogenic fluids in converging-diverging nozzles

Data are presented for the two phase choked flow of three cryogenic fluids - nitrogen, methane, and hydrogen - in four converging-diverging nozzles. The data cover a range of inlet stagnation conditions, all single phase, from well below to well above the thermodynamic critical conditions. In almost all cases the nozzle throat conditions were two phase. The results indicate that the choked flow rates were not very sensitive to nozzle geometry. However, the axial pressure profiles, especially the throat pressure and the point of vaporization, were very sensitive to both nozzle geometry and operating conditions. A modified Henry-Fauske model correlated all the choked flow rate data to within + or - 10 percent. Neither the equilibrium model nor the Henry-Fauske model predicted throat pressures well over the whole range of data. Above the thermodynamic critical temperature the homogeneous equilibrium model was preferred for both flow rate and pressure ratio. The data of the three fluids could be normalized by the principle of corresponding states.

Simoneau, R. J.↗
DOE OSTI2020

Seismicity and Stress Associated With a Fluid-Driven Fracture: Estimating the Evolving Geometry

A coupled approach, combining the theory of rate- and state-dependent friction and methods from poroelasticity, forms the basis for a quantitative relationship between displacements and fluid leak-off from a growing fracture and changes in the rate of seismic events in the region surrounding the fracture. Poroelastic Green's functions link fracture aperture changes and fluid flow from the fracture to changes in the stress field and pore pressure in the adjacent formation. The theory of rate- and state-dependent friction provides a connection between Coulomb stress changes and variations in the rate of seismic events. Numerical modeling indicates that the Coulomb stress changes can vary significantly between formations with differing properties. The relationship between the seismicity rate changes and the changes in the formation stresses and fluid pressure is nonlinear, but a transformation produces a quantity that is linearly related to the aperture changes and fluid leak-off from the fracture. The methodology provides a means for mapping changes in seismicity into fracture aperture changes and to image an evolving fracture. An application to observed microseismicity associated with a hydrofracture reveals asymmetric fracture propagation within two main zones, with extended propagation in the upper zone. Furthermore, the time-varying volume of the fracture agrees with the injected volume, given by the integration of rate changes at the injection well, providing validation of the estimated aperture changes.

58 GEOSCIENCES↗
DOE OSTI2020

Three-dimensional signatures of self-similarity in a high-energy-density plasma shear-driven mixing layer

A hydrodynamic shear mixing layer experiment at the National Ignition Facility had previously demonstrated Eulerian scaling of integrated, late-time quantities, including turbulent kinetic energy. In this manuscript, the experiment is repeated with new materials. In this work, using the new dataset, we demonstrate that Euler-number scalings hold not just for late time, but dynamically throughout the experiment, for measurements in all three spatial dimensions. Incorporating the dynamic scaling leads to an enhanced calculation that the heavier of the two scaled experiments has approached three generations of mergers of its primary instability's structures and a consistent observation of such a merger in action in the lighter of the two scaled experiments. Furthermore, the improved scrutiny of the time evolution of instability structures leads to sharper estimates of turbulent kinetic energy, including a demonstration of different behaviors correlating with surface roughness (quantitatively consistent with transitions between laminar and turbulent initial states), as predicted by a Reynolds-averaged turbulent model, which evidently correctly handles the differing shock-roughness interactions to drive its internal state of the model into different regimes. Altogether, a picture arises of the analytical improvements in treating these variations (of times, densities, and roughnesses) as a unified whole and of multiple ways by which deviations from the scaling could indicate an onset of non-hydrodynamic behavior. Such deviations were not expected for these experiments (which models correctly indicated would remain hydrodynamic) but could be introduced by, for example, imposing external fields or increasing drive energy to test conditions relevant to inertial confinement fusion or other high-energy-density experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Toward learning Lattice Boltzmann collision operators

In this work, we explore the possibility of learning from data collision operators for the Lattice Boltzmann Method using a deep learning approach. We compare a hierarchy of designs of the neural network (NN) collision operator and evaluate the performance of the resulting LBM method in reproducing time dynamics of several canonical flows. In the current study, as a first attempt to address the learning problem, the data were generated by a single relaxation time BGK operator. We demonstrate that vanilla NN architecture has very limited accuracy. On the other hand, by embedding physical properties, such as conservation laws and symmetries, it is possible to dramatically increase the accuracy by several orders of magnitude and correctly reproduce the short and long time dynamics of standard fluid flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1991

A finite element model of conduction, convection, and phase change near a solid/melt interface

Detailed understanding of heat transfer and fluid flow is required for many aerospace thermal systems. These systems often include phase change and operate over a range of accelerations or effective gravitational fields. An approach to analyzing such systems is presented which requires the simultaneous solution of the conservation laws of energy, momentum, and mass, as well as an equation of state. The variable property form of the governing equations are developed in two-dimensional Cartesian coordinates for a Newtonian fluid. A numerical procedure for solving the governing equations is presented and implemented in a computer program. The Galerkin form of the finite element method is used to solve the spatial variation of the field variables, along with the implicit Crank-Nicolson time marching algorithm. Quadratic Langrangian elements are used for the internal energy and the two components of velocity. Linear Lagrangian elements are used for the pressure. The location of the solid/liquid interface as well as the temperatures are determined form the calculated internal energy and pressure. This approach is quite general in that it can describe heat transfer without phase change, phase change with a sharp interface, and phase change without an interface. Analytical results from this model are compared to those of other researchers studying transient conduction, convection, and phase change and are found to be in good agreement. The numerical procedure presented requires significant computer resources, but this is not unusual when compared to similar studies by other researchers. Several methods are suggested to reduce the computational times.

Viterna, Larry A.↗
NASA NTRS1997

Pulsatile Fluid Shear in Bone Remodeling

The objective of this investigation was to elucidate the sensitivity to transients in fluid shear stress in bone remodeling. Bone remodeling is clearly a function of the local mechanical environment which includes interstitial fluid flow. Traditionally, load-induced remodeling has been associated with low frequency (1-2 Hz) signals attributed to normal locomotion. McLeod and Rubin, however, demonstrated in vivo remodeling events associated with high frequency (15-30 Hz) loading. Likewise, other in vivo studies demonstrated that slowly applied strains did not trigger remodeling events. We therefore hypothesized that the mechanosensitive pathways which control bone maintenance and remodeling are differentially sensitive to varying rates of applied fluid shear stress.

Frangos, John A.↗
DOE OSTI2022

Probing multiscale dissolution dynamics in natural rocks through microfluidics and compositional analysis

Mineral dissolution significantly impacts many geological systems. Carbon released by diagenesis, carbon sequestration, and acid injection are examples where geochemical reactions, fluid flow, and solute transport are strongly coupled. The complexity in these systems involves interplay between various mechanisms that operate at timescales ranging from microseconds to years. Current experimental techniques characterize dissolution processes using static images that are acquired with long measurement times and/or low spatial resolution. These limitations prevent direct observation of how dissolution reactions progress within an intact rock with spatially heterogeneous mineralogy and morphology. We utilize microfluidic cells embedded with thin rock samples to visualize dissolution with significant temporal resolution (100 ms) in a large observation window (3 × 3 mm). Here we injected acidic fluid into eight shale samples ranging from 8 to 86 wt % carbonate. The pre- and postreaction microstructures are characterized at the scale of pores (0.1 to 1 µm) and fractures (1 to 1,000 µm). We observe that nonreactive particle exposure, fracture morphology, and loss of rock strength are strongly dependent on both the relative volume of reactive grains and their distribution. Time-resolved images of the rock unveil the spatiotemporal dynamics of dissolution, including two-phase flow effects in real time and illustrate the changes in the fracture interface across the range of compositions. Moreover, the dynamical data provide an approach for characterizing reactivity parameters of natural heterogeneous samples when porous media effects are not negligible. The platform and workflow provide real-time characterization of geochemical reactions and inform various subsurface engineering processes.

58 GEOSCIENCES↗
DOE OSTI2014

Carbon Dioxide Sealing Capacity: Textural or Compositional Controls? (Final Scientific/Technical Report)

This research project is aiming to assess the carbon dioxide sealing capacity of most common seal-rocks, such as shales and non-fractured limestones, by analyzing the role of textural and compositional parameters of those rocks. We hypothesize that sealing capacity is controlled by textural and/or compositional parameters of caprocks. In this research, we seek to evaluate the importance of textural and compositional parameters affecting the sealing capacity of caprocks. The conceptual framework involves two testable end-member hypotheses concerning the sealing capacity of carbon dioxide reservoir caprocks. Better understanding of the elements controlling sealing quality will advance our knowledge regarding the sealing capacity of shales and carbonates. Due to relatively low permeability, shale and non-fractured carbonate units are considered relatively impermeable formations which can retard reservoir fluid flow by forming high capillary pressure. Similarly, these unites can constitute reliable seals for carbon dioxide capture and sequestration purposes. This project is a part of the comprehensive project with the final aim of studying the caprock sealing properties and the relationship between microscopic and macroscopic characteristics of seal rocks in depleted gas fields of Oklahoma Panhandle. Through this study we examined various seal rock characteristics to infer about their respective effects on sealing capacity in special case of replacing reservoir fluid with super critical carbon dioxide (scCO 2 ). To assess the effect of textural and compositional properties on scCO 2 maximum retention column height we collected 30 representative core samples in caprock formations in three counties (Cimarron, Texas, Beaver) in Oklahoma Panhandle. Core samples were collected from various seal formations (e.g., Cherokee, Keys, Morrowan) at different depths. We studied the compositional and textural properties of the core samples using several techniques. Mercury Injection Porosimetry (MIP), Scanning Electron Microscopy SEM, and Sedigraph measurements are used to assess the pore-throat-size distribution, sorting, texture, and grain size of the samples. Also, displacement pressure at 10% mercury saturation (Pd) and graphically derived threshold pressure (P c ) were determined by MIP technique. SEM images were used for qualitative study of the minerals and pores texture of the core samples. Moreover, EDS (Energy Dispersive X-Ray Spectrometer), BET specific surface area, and Total Organic Carbon (TOC) measurements were performed to study various parameters and their possible effects on sealing capacity of the samples. We found that shales have the relatively higher average sealing threshold pressure (P c ) than carbonate and sandstone samples. Based on these observations, shale formations could be considered as a promising caprock in terms of retarding scCO 2 flow and leakage into above formations. We hypothesized that certain characteristics of shales (e.g., 3 fine pore size, pore size distribution, high specific surface area, and strong physical chemical interaction between wetting phase and mineral surface) make them an efficient caprock for sealing super critical CO 2 . We found that the displacement pressure at 10% mercury saturation could not be the ultimate representative of the sealing capacity of the rock sample. On the other hand, we believe that graphical method, introduced by Cranganu (2004) is a better indicator of the true sealing capacity. Based on statistical analysis of our samples from Oklahoma Panhandle we assessed the effects of each group of properties (textural and compositional) on maximum supercritical CO 2 height that can be hold by the caprock. We conclude that there is a relatively strong positive relationship (+.40 to +.69) between supercritical CO 2 column height based on P c and hard/ soft mineral content index (ratio of minerals with Mohs hardness more than 5 over minerals with Mohs hardness less than 5) in both shales and limestone samples. Average median pore radius and porosity display a strong negative correlation with supercritical CO 2 retention column height. Also, increasing bulk density is positively correlated with the supercritical CO 2 retention column height. One of the most important factors affecting sealing capacity and consequently the height of supercritical CO 2 column is sorting of the pore throats. We observed a strong positive correlation between pore throat sorting and height of CO 2 retention column, especially in shales. This correlation could not be observed in limestone samples. It suggests that the pore throat sorting is more controlling the sealing capacity in shales and shales with well sorted pore throats are the most reliable lithology as seal. We observed that Brunauer–Emmett–Teller (BET) surface area shows a very strong correlation with CO 2 retention column height in limestone samples while BET surface area did not display significant correlation in shales. Pore structure based on SEM micrographs exhibits strong correlation with CO 2 retention column height in limestones. Both intercrystalline and vuggy structures have negative correlations while intergranular texture has positive correlation in limestone with respect to CO 2 retention column height. Textural effects observed on SEM micrographs did not show statistically significant correlation with supercritical CO 2 retention column height in shale samples. Finally, we showed that increasing hard/soft mineral index is strongly correlated with the displacement pressure in limestone samples. Vuggy texture displays a relatively strong and negative correlation with displacement pressure values at 10% mercury saturation in shale samples.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Interface models for three-dimensional Rayleigh–Taylor instability

We derive interface models for three-dimensional Rayleigh–Taylor instability (RTI), making use of a novel asymptotic expansion in the non-locality of the fluid flow. These interface models are derived for the purpose of studying universal features associated with RTI such as the Froude number in single-mode RTI, the predicted quadratic growth of the interface amplitude under multi-mode random perturbations, the optimal (viscous) mixing rates induced by the RTI and the self-similarity of horizontally averaged density profiles and the remarkable stabilization of the mixing layer growth rate which arises for the three-fluid two-interface heavy–light–heavy configuration, in which the addition of a third fluid bulk slows the growth of the mixing layer to a linear rate. Our interface models can capture the formation of small-scale structures induced by severe interface roll-up, reproduce experimental data in a number of different regimes and study the effects of multiple interface interactions even as the interface separation distance becomes exceedingly small. Compared with traditional numerical schemes used to study such phenomena, our models provide a computational speed-up of at least two orders of magnitude.

Mechanics↗
DOE OSTI2025

Beam Shaping and Oscillation for Metals Additive Manufacturing

Additive Manufacturing (AM) is pivotal for next-generation industrial manufacturing but is hampered by process control challenges leading to unwanted defects and sub-optimal material properties. Our research builds on prior work that utilized beam oscillation techniques to mitigate these issues (Wu et al., 2023). We focus on investigating the influence of beam oscillation frequency on melt pool dynamics, utilizing in-situ synchrotron material characterization from the ESRF (Chen et al., 2026). Preliminary findings indicated that beam oscillation results in unique fluid flow patterns and enhanced melting stability. These findings suggested the possibility of optimized AM process control, including minimized defects and tailored thermal histories. We proposed further research to explore the effects of beam oscillation frequencies comprehensively, using the in-situ characterization capabilities at the ESRF. The goal was to quantify this relationship, offering a pathway to more precise AM process control and optimization. Our research is expected to offer new scientific insights into the fluid mechanics and thermal behaviour in AM, contributing to its broader applicability and commercial viability.

36 MATERIALS SCIENCE↗
NASA NTRS2000

Convection Induced by Traveling Magnetic Fields in Semiconductor Melts

Axisymmetric traveling magnetic fields (TMF) can be beneficial for crystal growth applications. such as the vertical Bridgman, float zone or traveling heater methods. TMF induces a basic flow in the form of a single roll. This type of flow can enhance mass and heat transfer to the growing crystal. More importantly, the TMF Lorentz body force induced in the system can counterbalance the buoyancy forces, so the resulting convection can be much smaller and even the direction of it can be changed. In this presentation, we display basic features of this novel technique. In particular, numerical calculations of the Lorentz force for arbitrary frequencies will be presented along with induced steady-state fluid flow profiles. Also, numerical modeling of the TMF counter-balancing natural convection in vertical Bridgman systems will be demonstrated.

Konstantin, Mazuruk↗
NASA NTRS2016

Cryovolcanic Emplacement of Domes on Europa

Here we explore the hypothesis that certain domes on Europa may have been produced by the extrusion of viscous cryolavas. A new mathematical method for the emplacement and relaxation of viscous lava domes is presented and applied to putative cryovolcanic domes on Europa. A similarity solution approach is applied to the governing equation for fluid flow in a cylindrical geometry, and dome relaxation is explored assuming a volume of cryolava has been rapidly emplaced onto the surface. Nonphysical sin- gularities inherent in previous models for dome relaxation have been eliminated, and cryolava cooling is represented by a time-variable viscosity. We find that at the onset of relaxation, bulk kinematic viscosities may lie in the range between 10(exp 3) and 10(exp 6) sq m/s, while the actual fluid lava viscosity may be much lower. Plausible relaxation times to form the domes, which are linked to bulk cryolava rheology, are found to range from 3.6 days to 7.5 years. We find that cooling of the cryolava, while dominated by conduction through an icy skin, should not prevent fluids from advancing and relaxing to form domes within the timescales considered. Determining the range of emplacement conditions for putative cryolava domes will shed light on Europa's resurfacing history. In addition, the rheologies and compositions of erupted cryolavas have implications for subsurface cryomagma ascent and local surface stress conditions on Europa.

volcanism↗
NASA NTRS2019

Verification of Unstructured Grid Adaptation Components

Adaptive unstructured grid techniques have made limited impact on production analysis workflows where the control of discretization error is critical to obtaining reliable simulation results. Recent progress has matured a number of independent implementations of flow solvers, error estimation methods, and anisotropic grid adaptation mechanics. Known differences and previously unknown differences in grid adaptation components and their integrated processes are identified here for study. Unstructured grid adaptation tools are verified using analytic functions and the Code Comparison Principle. Three analytic functions with different smoothness properties are adapted to show the impact of smoothness on implementation differences. A scalar advection-diffusion problem with an analytic solution that models a boundary layer is adapted to test individual grid adaptation components. Laminar flow over a delta wing and turbulent flow over an ONERA M6 wing are verified with multiple, independent grid adaptation procedures to show consistent convergence to fine-grid forces and a moment. The scalar problems illustrate known differences in a grid adaptation component implementation and a previously unknown interaction between components. The wing adaptation cases in the current study document a clear improvement to existing grid adaptation procedures. The stage is set for the infusion of verified grid adaptation into production fluid flow simulations.

Park, Michael A.↗
DOE OSTI2024

High Power Density Thermal Energy Storage With Phase Change Material in Enhanced Compact Heat Exchangers

Abstract Performance of a novel ultracompact thermal energy storage (TES) heat exchanger, designed as a microchannel finned-tube exchanger is presented. With water as the heating–cooling fluid in the microchannels, a salt hydrate phase change material (PCM), lithium nitrate trihydrate (LiNO3 · 3H2O), was encased on the fin side. To establish the hypothesis that small-length-scale encasement (<3 mm) of PCM substantially enhances heat transfer to yield very high power-density energy storage, heat exchanger designs with 10 and 24 fins/inch were considered. They were subjected to thermal cycling, or repeated heating (melting) and cooling (freezing), with inlet fluid flow mimicking diurnal variation between 42 °C and 25 °C (representing typical arid-region conditions) over an accelerated time period. By employing salt self-seeding to obviate subcooling during cooling or recrystallization, the TES was found to exhibit stable long-term (100 heating–cooling cycles) operation with very high PCM-side heat transfer coefficients (∼100–500 W/m2 K) and storage power density (∼160–175 kW/m3). In fact, with optimization of heating–cooling fluid flowrate for given charging–discharging time period and exchanger size, power density >300 kW/m3 can be achieved. The results clearly establish that highly compact heat exchangers used as TES units can provide very high-performance alternatives to conventional ones.

Engineering↗
DOE OSTI2022

Hydromechanical impact of basement rock on injection-induced seismicity in Illinois Basin

Abstract The common explanation of observed injection-induced microseismicity is based on the measured stress state at the injection interval and the assumption that it remains the same in the vicinity. We argue here that representing the stress state in different geologic formations over the injection site with the single Mohr’s circle is insufficient due to local stratigraphic features and contrast in compressibilities of the involved formations. The role of hydromechanical coupling in the microseismic response is also crucial for the proper assessment of the problem. Thoroughly monitored Illinois Basin Decatur Project revealed the majority of CO 2 injection-associated microseismic events being originated in the crystalline basement. Even though basement faults can serve as the conduits for fluid flow—the predicted pressure increase seems to be insufficient to trigger seismicity. To address this issue, accurate laboratory measurements of rock properties from the involved formations are conducted. The pre-injection stress state and its evolution are evaluated with the hydromechanically coupled numerical model. It appears that the presence of an offset in a stiff competent layer affects the stress state in its vicinity. Therefore, both the pre-injection stress state and its evolution during the fluid injection should be addressed during the induced seismicity assessment.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2001

Numerical Study of g-Jitter Induced Double-Diffusive Convection

A finite element study is presented of double-diffusive convection driven by g-jitter in a microgravity environment. Mathematical formulations are presented and extensive simulations are carried out for g-jitter induced fluid flow, temperature distribution, and solutal transport in an alloy system under consideration for space flights. Computations include the use of idealized single-frequency and multi-frequency g-jitter as well as the real g-jitter data taken during an actual Space Shuttle fight. Little correlation is seen between these velocity components for the g-jitter components studied. The temperature field is basically undisturbed by convection because of a small Pr number for the fluid. The disturbance of the concentration field, however, is pronounced, and the local variation of the concentration follows the velocity oscillation in time. It is found that although the concentration field varies in both position and time, the local concentration gradient remains approximately constant in time. Numerical study further indicates that with an increase in g-jitter force (or amplitude), the nonlinear convective effects become much more obvious, which in turn drastically change the concentration fields. The simulated results computed using the g-jitter data taken during space flights show that both the velocity and concentration become random, following approximately the same pattern as the g-jitter perturbations.

Shu, Y.↗
DOE OSTI2022

Impact of particle size and particle-flow-wall coupling on pressurized oxy-combustion in the down-fired burner

Concerns over climate change have led to numerous efforts in developing low-carbon energy technologies. Pressurized oxy-combustion (POC) is a promising candidate to reduce carbon emissions in power generation. Designing an effective burner plays a vital role in developing new coal combustion technologies. Because of the high pressure, the volume fraction of coal particles at the fuel inlet of the burner of a pressurized oxy-combustor is close or even higher than the maximum limit that commercial CFD codes (e.g., ANSYS FLUENT) can handle. At this high particle volume fraction, the interactions among particles, fluid flow, and wall need to be re-evaluated. The present computational work is the first step of a systematic analysis of the particle influence, like releasing method, releasing locating, and particle size, in a pilot-scale POC combustor, developed at Washington University in St. Louis (WUSTL). In a POC process, pulverized coal is burned under elevated pressure and O2-CO2 environment. Specifically, a 15-bar, 100 kWth, POC combustor is modeled employing ANSYS FLUENT, using Reynolds-averaged Navier-Stokes (RANS) modeling. It is revealed that for this pilot-scale, pressurized burner, velocity profiles of the near-wall region in this POC facility exhibit some discrepancy against the near-wall turbulent flow velocity profile. In order to investigate the particle influence in the near-wall region, particle releasing location will be investigated. Numerical simulation results exhibit the coupling effect of turbulence flow and particles in this case. The particle size also demonstrates a great effect on particle trajectory, then further has an impact on flame stability and temperature.

Li, Lei↗
NASA NTRS1998

Numerical Investigation of Hot Gas Ingestion by STOVL Aircraft

This report compiles the various research activities conducted under the auspices of the NASA Grant NAG3-1026, "Numerical Investigation of Hot Gas Ingestion by STOVL Aircraft" during the period of April 1989 to April 1994. The effort involved the development of multigrid based algorithms and computer programs for the calculation of the flow and temperature fields generated by Short Take-off and Vertical Landing (STOVL) aircraft, while hovering in ground proximity. Of particular importance has been the interaction of the exhaust jets with the head wind which gives rise to the hot gas ingestion process. The objective of new STOVL designs to reduce the temperature of the gases ingested into the engine. The present work describes a solution algorithm for the multi-dimensional elliptic partial-differential equations governing fluid flow and heat transfer in general curvilinear coordinates. The solution algorithm is based on the multigrid technique which obtains rapid convergence of the iterative numerical procedure for the discrete equations. Initial efforts were concerned with the solution of the Cartesian form of the equations. This algorithm was applied to a simulated STOVL configuration in rectangular coordinates. In the next phase of the work, a computer code for general curvilinear coordinates was constructed. This was applied to model STOVL geometries on curvilinear grids. The code was also validated in model problems. In all these efforts, the standard k-Epsilon model was used.

Vanka, S. P.↗

Refine your search to explore more results.