Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluid mixing”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Numerical simulation of flow and mixing in fracture intersections

Fluid transport through fractured geological formations is strongly influenced by the redistribution of solutes at fracture intersections. In this study, we perform detailed numerical simulations of flow and scalar transport within the intersection of two smooth, planar fractures. The analysis focuses on the mixing ratio, the proportion of solute flux exiting along the outlet branch aligned with the primary inlet flow direction, relative to the total solute flux at the outlets. We systematically investigate how the mixing ratio varies with four key parameters: Peclet number, Reynolds number, flow rate ratio between outlet branches, and fracture intersection angle. Results show that the mixing ratio decreases with increasing Peclet number and outlet flow rate ratio, consistent with reduced diffusive spreading and enhanced streamline routing. While low Reynolds numbers have minimal impact, inertial effects at higher Reynolds numbers significantly increase the mixing ratio. Additionally, acute and obtuse intersection angles alter flow partitioning and modify the solute distribution at the outlets. These findings provide a quantitative basis for incorporating physically realistic mixing behavior—intermediate between complete mixing and streamline-following assumptions—into network-scale transport models. The results have direct relevance to subsurface energy systems, including geothermal energy production, carbon sequestration, and contaminant remediation.

58 GEOSCIENCES↗
DOE OSTI2023

Mixed imbibition controls the advance of wetting fluid in multiscale geological media

Imbibition plays a central role during multiphase flows in geological media; yet, our understanding of imbibition dynamics across the multiscale features of rocks, including fractures, microcracks, and heterogeneous matrices is limited. Here, we deployed a unique combination of micro- and clinical X-ray CT imaging to probe directly the spatiotemporal evolution of imbibition in a natural, multiscale shale sample. The method enabled unraveling of mixed imbibition behavior, with spontaneous imbibition rates in the matrix pores and microcracks driven by capillarity proceeding as $\propto t^{-1/2}$, whereas pressure-driven flows in the macroscopic fractures yield apparent imbibition rates with $\propto t^{-1}$ scaling. Regardless of the aperture, we showed that spontaneous imbibition proceeds uniformly away from the fractures. The effects of microcracks are further investigated by comparing the dynamics of the imbibition fronts with and without microcracks. Finally, these results highlight the paramount role of microcracks on anomalous front roughening and the need to account for microcracks appropriately in the development of predictive theoretical models.

58 GEOSCIENCES↗
NASA NTRS1993

Mixing and transient interface condensation of a liquid hydrogen tank

Experiments were conducted to investigate the effect of axial jet-induced mixing on the pressure reduction of a thermally stratified liquid hydrogen tank. The tank was nearly cylindrical, having a volume of about 0.144 cu m with 0.559 m in diameter and 0.711 m long. A mixer/pump unit, which had a jet nozzle outlet of 0.0221 m in diameter was located 0.178 m from the tank bottom and was installed inside the tank to generate the axial jet mixing and tank fluid circulation. The liquid fill and jet flow rate ranged from 42 to 85 percent (by volume) and 0.409 to 2.43 cu m/hr, respectively. Mixing tests began with the tank pressure ranging from 187.5 to 238.5 kPa at which the thermal stratification results in 4.9 to 6.2 K liquid sub cooling. The mixing time and transient vapor condensation rate at the liquid-vapor interface are determined. Two mixing time correlations, based on the thermal equilibrium and pressure equilibrium, are developed. Both mixing time correlations are expressed as functions of system and buoyancy parameters and compared well with other experimental data. The steady state condensation rate correlation of Sonin et al. based on steam-water data is modified and expressed as a function of jet subcooling. The limited liquid hydrogen data of the present study shows that the modified steady state condensation rate correlation may be used to predict the transient condensation rate in a mixing process if the instantaneous values of jet sub cooling and turbulence intensity at the interface are employed.

Lin, C. S.↗
NASA NTRS2005

Fault-Tolerant Heat Exchanger

A compact, lightweight heat exchanger has been designed to be fault-tolerant in the sense that a single-point leak would not cause mixing of heat-transfer fluids. This particular heat exchanger is intended to be part of the temperature-regulation system for habitable modules of the International Space Station and to function with water and ammonia as the heat-transfer fluids. The basic fault-tolerant design is adaptable to other heat-transfer fluids and heat exchangers for applications in which mixing of heat-transfer fluids would pose toxic, explosive, or other hazards: Examples could include fuel/air heat exchangers for thermal management on aircraft, process heat exchangers in the cryogenic industry, and heat exchangers used in chemical processing. The reason this heat exchanger can tolerate a single-point leak is that the heat-transfer fluids are everywhere separated by a vented volume and at least two seals. The combination of fault tolerance, compactness, and light weight is implemented in a unique heat-exchanger core configuration: Each fluid passage is entirely surrounded by a vented region bridged by solid structures through which heat is conducted between the fluids. Precise, proprietary fabrication techniques make it possible to manufacture the vented regions and heat-conducting structures with very small dimensions to obtain a very large coefficient of heat transfer between the two fluids. A large heat-transfer coefficient favors compact design by making it possible to use a relatively small core for a given heat-transfer rate. Calculations and experiments have shown that in most respects, the fault-tolerant heat exchanger can be expected to equal or exceed the performance of the non-fault-tolerant heat exchanger that it is intended to supplant (see table). The only significant disadvantages are a slight weight penalty and a small decrease in the mass-specific heat transfer.

Izenson, Michael G.↗
DOE OSTI2021

Development and Field Testing Novel Natural Gas Surface Process Equipment for Replacement of Water as Primary Hydraulic Fracturing Fluid (Final Report)

Southwest Research Institute® (SwRI ® ), Schlumberger Technology Corporation (SLB), and Chevron Corporation ® (Chevron) developed a novel, optimized, and lightweight process for natural gas (NG) to replace water as a low-cost fracturing medium with a low environmental impact. Hydraulic fracturing is used to increase oil and NG production by injecting high-pressure fluid, primarily water, into a rock formation, which fractures the rock and releases trapped oil and NG. This method was developed to increase yield and make feasible production areas that would not otherwise be viable for large-scale oil and NG extraction using traditional drilling technologies. Since the fracturing fluid is composed of approximately 90% water, one of the principal drawbacks to hydraulic fracturing is its excessive water use and associated large environmental footprint. According to recent data collected at fracture sites within the United States, fracturing applications in North America can consume as much as 9 million gallons of water per well. During the fracturing process, some of the fracturing fluid is permanently lost and the portion that is recovered is contaminated by both fracturing chemicals and dissolved solids from the formation. The recovered water or flow-back represents a significant environmental challenge, as it must be treated before it can be reintroduced into the natural water system. Although there is some recycling for future fracturing, the majority of the flow-back water is hauled from the well site to a treatment facility or to an injection well for permanent underground disposal. To mitigate these issues, an optimized, lightweight and modular surface process using NG to replace a majority of the water was developed as a cost-effective and environmentally clean fracturing fluid. Using NG will result in significantly less consumption since the gas that is injected as a fracturing fluid will be mixed with the formation gas and extracted as if it were from the formation itself. This process will minimize the collection, waste, and treatment of large amounts of water and will reduce the environmental impact of transporting and storing the fracturing fluid. There are two major steps involved in utilizing NG as the primary fracturing medium: (1) increasing the supply pressure of NG to wellhead pressures suitable for fracturing and (2) mixing the required chemicals and proppant needed for the fracturing process at these elevated pressures. The second step (NG-proppant mixing at elevated pressures) still requires demonstration; but very similar processes have been demonstrated in the field with other gases such as nitrogen (N 2 ) and carbon dioxide (CO 2 ). However, the first step (a compact, on-site unit for generating high-pressure NG at costs feasible for fracturing) has not been developed and is currently not commercially available. Due to the inherent compressibility of NG, more energy is required to compress the gas than what is required to pump water (or other incompressible liquids) to the very high pressure required for downhole injection. This project aimed to develop a novel, hybrid method to overcome this challenge. The project accomplishments and findings are discussed in this report. Ultimately, the research and development efforts described herein demonstrate that fracturing with NG foam is a feasible alternative to using water.

02 PETROLEUM↗
NASA NTRS1971

Fluid phase analyzer Patent

Mixed liquid and vapor phase analyzer design with thermocouples for relative heat transfer measurement

Hays, L. G.↗
NASA NTRS1944

The "limiting line" in mixed subsonic and supersonic flow of compressible fluids

It is well known that the vorticity for any fluid element is constant if the fluid is non-viscous and the change of state of the fluid is isentropic. When a solid body is placed in a uniform stream, the flow far ahead of the body is irrotational. Then if the flow is further assumed to be isentropic, the vorticity will be zero over the whole filed of flow. In other words, the flow is irrotational. For such flow over a solid body, it is shown by Theodorsen that the solid body experiences no resistance. If the fluid has a small viscosity, its effect will be limited in the boundary layer over the solid body and the body will have a drag due to the skin friction. This type of essentially isentropic irrotational flow is generally observed for a streamlined body placed in a uniform stream, if the velocity of the stream is kept below the so-called "critical speed." At the critical speed or rather at a certain value of the ratio of the velocity of the undisturbed flow and the corresponding velocity of sound, shock waves appear. This phenomenon is called the "compressibility bubble." Along a shock wave, the change of state of the fluid is no longer isentropic, although still adiabatic. This results in an increase in entropy of the fluid and generally introduces vorticity in an originally irrotational flow. The increase in entropy of the fluid is, of course, the consequence of changing part of the mechanical energy into heat energy. In other words, the part of fluid affected by the shock wave has a reduced mechanical energy. Therefore, with the appearance of shock waves, the wake of the streamline body is very much widened, and the drag increases drastically. Furthermore, the accompanying change in the pressure distribution over the body changes the aerodynamic moment acting on it and in the case of an airfoil decreases the lift force. All these consequences of the breakdown of isentropic irrotational flow are generally undesirable in applied aerodynamics. Its occurrence should be delayed as much as possible by modifying the shape or contour of the body. However, such endeavor will be very much facilitated if the cause or the criterion for the breakdown can be found first.

MACH NUMBER↗
NASA NTRS2016

Modeling Ullage Dynamics of Tank Pressure Control Experiment during Jet Mixing in Microgravity

A CFD model for simulating the fluid dynamics of the jet induced mixing process is utilized in this paper to model the pressure control portion of the Tank Pressure Control Experiment (TPCE) in microgravity1. The Volume of Fluid (VOF) method is used for modeling the dynamics of the interface during mixing. The simulations were performed at a range of jet Weber numbers from non-penetrating to fully penetrating. Two different initial ullage positions were considered. The computational results for the jet-ullage interaction are compared with still images from the video of the experiment. A qualitative comparison shows that the CFD model was able to capture the main features of the interfacial dynamics, as well as the jet penetration of the ullage.

active control↗
NASA NTRS1994

Computations of drop collision and coalescence

Computations of drops collision and coalescence are presented. The computations are made possible by a recently developed finite difference/front tracking technique that allows direct solutions of the Navier-Stokes equations for a multi-fluid system with complex, unsteady internal boundaries. This method has been used to examine the boundaries between the various collision modes for drops of equal size and two examples, one of a 'reflective' collision and another of a 'grazing' collision is shown. From drops of unequal size, coalescence can result in considerable mixing between the fluid from the small and the large drop. This problem is discussed and one example showed. In many cases it is necessary to account also for heat transfer along with the fluid mechanics. We show two preliminary results where we are using extensions of the method to simulate such a problem. One example shows pattern formation among many drops moving due to thermal migration, the other shows unstable evolution of a solidification front.

Tryggvason, Gretar↗
NASA NTRS1997

Flowfield-Dependent Mixed Explicit-Implicit (FDMEL) Algorithm for Computational Fluid Dynamics

Despite significant achievements in computational fluid dynamics, there still remain many fluid flow phenomena not well understood. For example, the prediction of temperature distributions is inaccurate when temperature gradients are high, particularly in shock wave turbulent boundary layer interactions close to the wall. Complexities of fluid flow phenomena include transition to turbulence, relaminarization separated flows, transition between viscous and inviscid incompressible and compressible flows, among others, in all speed regimes. The purpose of this paper is to introduce a new approach, called the Flowfield-Dependent Mixed Explicit-Implicit (FDMEI) method, in an attempt to resolve these difficult issues in Computational Fluid Dynamics (CFD). In this process, a total of six implicitness parameters characteristic of the current flowfield are introduced. They are calculated from the current flowfield or changes of Mach numbers, Reynolds numbers, Peclet numbers, and Damkoehler numbers (if reacting) at each nodal point and time step. This implies that every nodal point or element is provided with different or unique numerical scheme according to their current flowfield situations, whether compressible, incompressible, viscous, inviscid, laminar, turbulent, reacting, or nonreacting. In this procedure, discontinuities or fluctuations of an variables between adjacent nodal points are determined accurately. If these implicitness parameters are fixed to certain numbers instead of being calculated from the flowfield information, then practically all currently available schemes of finite differences or finite elements arise as special cases. Some benchmark problems to be presented in this paper will show the validity, accuracy, and efficiency of the proposed methodology.

Garcia, S. M.↗
DOE OSTI2021

Rapid charging made practical in graphite-based lithium batteries: surface-acoustic wave turbulent electrolyte mixing to overcome diffusion limited charging rates [Final Report]

One of the key limits in rapidly recharging a lithium-ion battery is the depletion of lithium ions within the electrolyte adjacent the anode during charging and long diffusion time to overcome this depletion. It also causes dendrite formation, inefficient use of the lithium, and battery degradation over many charge-discharge cycles. Because the liquid electrolyte remains quiescent and unmixed, this depletion layer’s depth rapidly grows to match the anode-cathode separation distance at even modest charge rates. The solution proposed by PI Prof. James Friend and Co-PI Prof. Ping Liu, both from the University of California, San Diego, is to mix the electrolyte and minimize the Li+ ion concentration gradient during charging. Even in the presence of the separator, the charging rate could then be significantly increased. They propose to accomplish this using surface acoustic wave (SAW)-driven acoustic streaming, a technique employing 10-mW fingernail-sized solid-state devices from the telecommunications industry to drive turbulent mixing to submicron length scales in a manner completely compatible with the typical 20700 and 18650 cells used in electric vehicles. 0.1. Turbulent acoustic streaming mixes the electrolyte during charging—even with separator Our proposed 100-MHz SAW device used to recirculate the electrolyte is compatible with lithium- ion battery electrochemistry, as it is made in our lab of single-crystal lithium niobate. Their solution is straightforward to drive from a DC power source alongside the signal provided during battery charging. The device is only needed during charging and does not consume power dur- ing battery discharge. Uniquely, SAW generates extreme accelerations of over 1 billion meters per second squared in the fluid, driving turbulent mixing from centimeter to submicron length scales, even through the porous separator materials that tend to be used in batteries, all while avoiding interfering with the anode’s solid electrolyte interphase layer that forms during use. By employing novel fluid mechanics, the investigators are proposing a new direction for battery research away from direct use of materials science and electrochemistry. The chemistry agnostic solution may be employed in any battery chemistry that makes use of liquid electrolytes, providing a broader transformative benefit to the battery research community. They are to produce a series of prismatic and 20700-cell sized 2 Ah batteries capable of being charged and discharged at least 500 times without more than 20% loss in battery capacity, and to provide analysis tools useful to the battery research and development community for adopting this approach for other battery chemistries and configurations.

25 ENERGY STORAGE↗
NASA NTRS1974

Heater-mixer for stored fluids

A fluid storage vessel for containing cryogenic fluids is described. The storage vessel contains an auxiliary chamber which is connected to the main container by a jet nozzle. The wall of the auxiliary vessel is heat cycled to produce a corresponding expansion and contraction of the fluid within the auxiliary chamber. This action causes heating and mixing of the stored fluid by means of jetting the expanded fluid to and from relative to the stored fluid contents of the vessel.

Canning, T. N.↗
NASA NTRS2007

Runtime and Pressurization Analyses of Propellant Tanks

Multi-element unstructured CFD has been utilized at NASA SSC to carry out analyses of propellant tank systems in different modes of operation. The three regimes of interest at SSC include (a) tank chill down (b) tank pressurization and (c) runtime propellant draw-down and purge. While tank chill down is an important event that is best addressed with long time-scale heat transfer calculations, CFD can play a critical role in the tank pressurization and runtime modes of operation. In these situations, problems with contamination of the propellant by inclusion of the pressurant gas from the ullage causes a deterioration of the quality of the propellant delivered to the test article. CFD can be used to help quantify the mixing and propellant degradation. During tank pressurization under some circumstances, rapid mixing of relatively warm pressurant gas with cryogenic propellant can lead to rapid densification of the gas and loss of pressure in the tank. This phenomenon can cause serious problems during testing because of the resulting decrease in propellant flow rate. With proper physical models implemented, CFD can model the coupling between the propellant and pressurant including heat transfer and phase change effects and accurately capture the complex physics in the evolving flowfields. This holds the promise of allowing the specification of operational conditions and procedures that could minimize the undesirable mixing and heat transfer inherent in propellant tank operation. It should be noted that traditional CFD modeling is inadequate for such simulations because the fluids in the tank are in a range of different sub-critical and supercritical states and elaborate phase change and mixing rules have to be developed to accurately model the interaction between the ullage gas and the propellant. We show a typical run-time simulation of a spherical propellant tank, containing RP-1 in this case, being pressurized with room-temperature nitrogen at 540 R. Nitrogen, shown in blue on the right-hand side of the figures, enters the tank from the diffuser at the top of the figures and impinges on the RP-1, shown in red, while the propellant is being continuously drained at the rate of 1050 lbs/sec through a pipe at the bottom of the tank. The sequence of frames in Figure 1 shows the resultant velocity fields and mixing between nitrogen and RP-1 in a cross-section of the tank at different times. A vortex is seen to form in the incoming nitrogen stream that tends to entrain propellant, mixing it with the pressurant gas. The RP-1 mass fraction contours in Figure 1 are also indicative of the level of mixing and contamination of the propellant. The simulation is used to track the propagation of the pure propellant front as it is drawn toward the exit with the evolution of the mixing processes in the tank. The CFD simulation modeled a total of 10 seconds of run time. As is seen from Figure 1d, after 5.65 seconds the propellant front is nearing the drain pipe, especially near the center of the tank. Behind this pure propellant front is a mixed fluid of compromised quality that would require the test to end when it reaches the exit pipe. Such unsteady simulations provide an estimate of the time that a high-quality propellant supply to the test article can be guaranteed at the modeled mass flow rate. In the final paper, we will discuss simulations of the LOX and propellant tanks at NASA SSC being pressurized by an inert ullage. Detailed comparisons will be made between the CFD simulations and lower order models as well as with test data. Conditions leading to cryo collapse in the tank will also be identified.

Field, Robert E.↗
DOE OSTI2024

Impulse-driven release of gas-encapsulated drops

Gas-encapsulated drops, much like antibubbles, are drops enclosed in a bubble within a liquid. They show potential as payload carriers in fluid transport and mixing techniques where sound waves can be leveraged to induce the collapse of the gas core and the subsequent release of the drop. Here, the interaction of millimetre-sized gas-encapsulated drops with impulsive laser-induced shock waves is investigated to gain fundamental insights into the release process. Experimental synchrotron X-ray phase contrast imaging, which allows the drop dynamics to be visualised inside the encapsulating bubble, is complemented by numerical simulations to study the intricate physics at play. Three drop dynamical release regimes are discovered, namely the drop impact , partial deposition and jet impact regimes. The regime type is mainly dependent on the shape of the bubble interface impacting the drop and the associated Weber and Reynolds numbers. The drop dynamics of the drop impact and partial deposition regimes show similarities with the canonical configuration of drops impacting flat liquid surfaces, whereas the jet impact regime resembles binary drop collisions, which allows existing scaling laws to be applied to describe the underlying processes. The release of the drop is investigated numerically. The time evolution of the drop dissemination within the surrounding liquid discloses enhanced mixing for dynamics involving high Weber and Reynolds numbers such as the drop impact and jet impact regimes.

42 ENGINEERING↗
DOE OSTI2023

Biological carbon pump estimate based on multidecadal hydrographic data

The transfer of photosynthetically produced organic carbon from surface to mesopelagic waters draws carbon dioxide from the atmosphere. However, current observation-based estimates disagree on the strength of this biological carbon pump (BCP). Earth system models (ESMs) also exhibit a large spread of BCP estimates, indicating limited representations of the known carbon export pathways. Here we use several decades of hydrographic observations to produce a top-down estimate of the strength of the BCP with an inverse biogeochemical model that implicitly accounts for all known export pathways. Our estimate of total organic carbon (TOC) export at 73.4 m (model euphotic zone depth) is 15.00 ± 1.12 Pg C year –1 , with only two-thirds reaching 100 m depth owing to rapid remineralization of organic matter in the upper water column. Partitioned by sequestration time below the euphotic zone, τ, the globally integrated organic carbon production rate with τ > 3 months is 11.09 ± 1.02 Pg C year –1 , dropping to 8.25 ± 0.30 Pg C year –1 for τ > 1 year, with 81% contributed by the non-advective-diffusive vertical flux owing to sinking particles and vertically migrating zooplankton. Nevertheless, export of organic carbon by mixing and other fluid transport of dissolved matter and suspended particles remains regionally important for meeting the respiratory carbon demand. Here, the temperature dependence of the sequestration efficiency inferred from our inversion suggests that future global warming may intensify the recycling of organic matter in the upper ocean, potentially weakening the BCP.

58 GEOSCIENCES↗
NASA NTRS1972

Effects of circular geometry in simulation of convection in rotating spacecraft tanks

For computational simulation of the convection and mixing of stratified fluid in a rotating tank (such as used in Apollo flights) with time-dependent rotation, the Navier-Stokes convection problem was formulated for a circular tank configuration. The final equations results from a general approximate theory for combined forced and contained natural convection in a time-dependent rotating system. The equations are cast in terms of vorticity and stream function in a form convenient for computation, with a transformed coordinate system, and appropriate boundary conditions are derived. Accurate representations for the cryogenic supercritical oxygen thermodynamic properties are used in the computations, and an efficient numerical finite difference scheme and computational method are employed.

Martin, E. D.↗
NASA NTRS1972

Bacteria detection instrument and method

A method and apparatus for screening a sample fluid for bacterial presence are disclosed wherein the fluid sample is mixed with culture media of sufficient quantity to permit bacterial growth in order to obtain a test solution. The concentration of oxygen dissolved in the test solution is then monitored using the potential difference between a reference electrode and a noble metal electrode which are in contact with the test solution. The change in oxygen concentration which occurs during a period of time as indicated by the electrode potential difference is compared with a detection criterion which exceeds the change which would occur absent bacteria.

Renner, W.↗