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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

DOE OSTI2025

High-Temperature Gas-Cooled Reactors Multiphysics Simulation Demonstration and Code Validation

This study presents a comprehensive benchmarking and verification effort of several thermal-hydraulic and multiphysics capabilities for high-temperature gas-cooled reactor applications. The first part of this effort focuses on the running-in verification of Griffin’s multiphysics capabilities, specifically for simulating the evolution of pebble-bed reactor cores from startup to equilibrium. Since Fiscal Year 2024, improvements and enhancements have been implemented in Griffin, including simplifying the process to specify streamlines and developing the online cross-section generation capability. In the absence of validation data, code-to-code comparisons are conducted with kugelpy, showing good agreement for integral quantities like k-eff predictions and predictions for maximum power density. However, accuracy issues are noted for more detailed quantities like the spatial distribution of fission rate densities which will require further work to address. The second part of this report presents an improved System Analysis Module (SAM) core channel model where the effects of cross flow are considered during the pressurized loss of forced cooling transient, resulting in an improved agreement of the predicted pebble temperature with respect to the predictions from the SAM 2D porous media model. Additionally, the wall channeling effect due to variable porosity at the near wall region of the core is also investigated. Furthermore, to demonstrate Griffin’s online cross-section generation capability, a Multiphysics simulation is performed by coupling Griffin to the SAM core channel model. In the third part of the report, as a part of the Organisation for Economic Co-operation and Development/Nuclear Energy Agency (OECD/NEA) thermal-hydraulic code validation benchmark activity for a high-temperature gas-cooled reactor, the High Temperature Test Facility (HTTF) is investigated first using the NekRS computational fluid dynamics (CFD) code to study the flow mixing phenomenon in the lower plenum of the facility. Then, code-to-code and code-to-data comparisons are performed for Test PG27, which is a pressurized conduction cooldown (PCC) test, using five different codes by six organizations from five countries. The different simulations show good agreements in terms of the general trend but there are differences in some results such as the peak temperatures of different regions and heat removal rate.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Influence of gasoline fuel formulation on lean autoignition in a mixed-mode-combustion (deflagration/autoignition) engine

We report stoichiometric spark-ignition engines suffer efficiency penalties due to throttling losses at low loads, a low specific-heat ratio of the stoichiometric working fluid, and limits on compression ratio due to end-gas autoignition leading to undesirable knocking. Mixed-Mode Combustion (MMC) mitigates these shortcomings by using a lean working fluid where a spark-initiated pilot-stabilized deflagrative flame front is followed by controlled end-gas autoignition. This MMC study investigates the effects of initial conditions (intake air temperature, intake pressure, equivalence ratio, and intake oxygen fraction) on autoignition tendency of four gasoline-range fuels with varying properties and composition. The use of fuels with varying octane sensitivity (S) allowed exploring the importance of low-temperature heat release in triggering autoignition. Fuels with high S were less reactive for conditions that promote low-temperature chemistry (operation at high intake air pressure or without N2 dilution). Conversely, an Alkylate fuel with low S showed a greater autoignition resistance at operating conditions that were unfavorable for low-temperature chemistry. Next, the effect of residual gas composition on autoignition tendency of fuels was examined with a chemical-kinetics model. Among the various molecules in the residual gas, nitric oxide (NO) enhanced the low-temperature chemistry and increased the autoignition tendency most significantly. The fuels’ autoignition response to increasing NO amount corroborates the experimental observations. Next, the sequential autoignition of the end-gas was assessed to be less impacted by thermal stratification because of lean mixtures showing relatively less low-temperature chemistry, when compared to stoichiometric mixtures. Next, the effect of changing equivalence ratio on the autoignition was found to be similar for all fuels, regardless of their S. With changing intake air temperature, the response of fuels’ autoignition tendency depended on the dilution level used. At high dilution (i.e. low intake [O2]), fuels’ reactivity increased with increasing intake air temperature. In contrast, for operation without dilution, the autoignition tendency of the low-S Alkylate fuel decreased with increasing intake air temperature, while that of high-S High Cycloalkane fuel still increased with increasing intake air temperature. In conclusion, conventional octane metrics (RON and MON) have utility in assessing the autoignition tendency under lean MMC operation. Moreover, the fuel requirements for MMC align with that of stoichiometric operation: i.e., high RON and high S fuels are desirable for stable non-knocking operation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Capturing Secondary Kinetic Instabilities in Three‐Dimensional Dayside Reconnection Using an Improved Gradient‐Based Closure

Magnetic reconnection is a highly dynamic process that excites a wide variety of kinetic waves and instabilities. Transverse current sheet instabilities such as the lower-hybrid drift and secondary drift-kink instabilities in particular have been shown by kinetic simulations to modify the reconnection and introduce significant turbulence and mixing to the reconnection layer. Past studies using the ten-moment fluid model to capture important kinetic physics such as the electron inertia and full representation of the pressure tensor proved advantageous to a two-fluid representation of reconnection, but the model struggled when using a local relaxation closure for the heat flux to replicate the current sheet instabilities and subsequent mixing seen in kinetic simulations. This work uses the Gkeyll software framework to perform simulations of asymmetric reconnection based on the 16 October 2015 MMS crossing of a diffusion region, the Burch event. An improved gradient-based heat flux closure is implemented, showing significant improvement in secondary kinetic instabilities that grow in the current sheet. These instabilities generate turbulence which leads to growth of secondary magnetic islands and flux ropes.

Bradshaw, K. [Princeton University, NJ (United Sta↗

DOE OSTI2024

Hydrodynamical Simulations of Proton Ingestion Flashes in Type I X-Ray Bursts

We perform the first multidimensional fluid simulations of thermonuclear helium ignition underneath a hydrogen-rich shell. This situation is relevant to Type I X-ray bursts on neutron stars that accrete from a hydrogen-rich companion. Using the low-Mach number fluid code MAESTROeX, we investigate the growth of the convection zone due to nuclear burning, and the evolution of the chemical abundances in the atmosphere of the star. We also examine the convective boundary mixing processes that cause the evolution to differ significantly from previous one-dimensional simulations that rely on mixing-length theory. We find that the convection zone grows outward as penetrating fluid elements cool the overlying radiative layer, rather than directly from the increasing entropy of the convection zone itself. Simultaneously, these flows efficiently mix composition, carrying carbon out of and protons into the convection zone even before contact with the hydrogen shell. We discuss the implications of these effects for future modeling of these events and observations.

astrophysical fluid dynamics↗

DOE OSTI2022

Analysis of Compressibility Effects and Nonlinear Property Variations in a Supercritical CO2 Mixing Layer

Calculations are performed using the large eddy simulation technique to qualitatively assess the effects of fluid compressibility and thermodynamic nonlinearities on the dynamics of a supercritical flow field. Resulting data are also used to interrogate instantaneous subfilter-scale fields, as modeled by the mixed-dynamic Smagorinsky closure and gradient diffusion model. A three-dimensional wall-resolved calculation of a spatially evolving mixing layer composed solely of carbon dioxide is performed. This is followed by a coarse grid calculation to examine model performance as a function of resolution. Theoretical analysis indicates that the partial derivatives of density with respect to pressure and temperature are important modulators of the temperature and pressure field evolution, respectively. Results pertaining to subfilter velocity-velocity and velocity-temperature fields indicate that the models are most stressed in the regions where compressibility and thermo-fluid nonlinearities dominate the physics.

42 ENGINEERING↗

DOE OSTI2026

Characteristics of Fluid‐Solid Interaction Constitutive Models Within Poroelastodynamics at Higher Strain‐Rates and Large Deformations Implemented in 1D

The large deformation, mixed formulation, finite element (FE) modeling approach presented in Irwin et al. 2024 is extended herein to include improved constitutive models for representing dynamic solid-fluid interactions at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger overpressure magnitudes (𝒪⁡(1⁢0 2 )⁢kPa) within a biphasic soft porous material using Theory of Porous Media (TPM) at finite strain. Specifically, these constitutive modeling improvements are the following: (i) a more physically robust constitutive model for pore fluid seepage velocity via inclusion of pore fluid viscous stress, and (ii) a modified deformation-dependent-permeability model and updated hyperelastic constitutive model better suited for handling larger volumetric compressions and extensions. The novelty of the present work is mainly the contribution (i): inclusion of pore fluid viscous stress at higher strain-rate and large deformations, which requires 𝐶 1 continuity in the weak formulation, accomplished by employing Hermite cubic interpolation functions within a mixed nonlinear poromechanical finite element formulation. In (ii), the model is updated to weakly enforce solid phase incompressibility, such that this assumption is not violated numerically, which provides improved numerical stability for achieving larger overpressure magnitudes on 𝒪⁡(1⁢0 2 ) kPa, which were not achievable with the previous Kozeny–Carman model in Irwin et al. 2024. Also in (ii), the volumetric part of the solid skeleton free energy function is modified to ensure proper bounds on the solid skeleton Jacobian of deformation 𝐽 s related to incompressibility of the solid phase. Uniaxial strain, unidirectional flow examples at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger deformations (up to 0.2 (or 20%) nominal axial strain) demonstrate the improved physical representation—and numerical stability—of these constitutive model improvements.

42 ENGINEERING↗

DOE OSTI2020

The αs and θs in Rayleigh–Taylor and Richtmyer–Meshkov instabilities

The alphas (α b , α s ) and thetas (θ b , θ s ) are the essential parameters describing the asymptotic growth behavior of the mixing layer driven by Rayleigh–Taylor (RT) and Richtmyer–Meshkov (RM) instabilities at the interface between the fluids. Extensive simulations and experiments have been conducted in the determination of these parameters. The results depend on numerical algorithms, resolutions, and the accuracy of experimental measurements. In this article, we present a unique relationship between the alphas and thetas by applying a simple buoyancy-drag mixing model. Our first main result is that, for the same fluids, these parameters are not independent because they are fundamentally connected by physics principles and they all are functions of the drag coefficients of the fluids. As a consequence, the alphas and thetas are intrinsically coupled. All are functions of any single one of these. Additionally, our second main result is that the RT bubble growth rate (α b ) is determined by the bubble competition process – a merger model – at the edge of the mixing layer. Together, these two main results determine all the alphas and thetas by strictly theoretical reasoning. Applications to both hydro and laser driven experiments at OMEGA are presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Advanced CO 2 Capture Solvent Systems for Dynamic Power Generation: Quarterly Research Performance Progress Report, QR4 (Q4FY24)

We developed an integrated Computational Fluid Dynamics (CFD) model to simulate the multi-physics coupled cooling process of mixed gas by cold water within a Direct Contact Cooler (DCC) equipped with a rotating packing bed (RPB). The model captures the interactions between fluid dynamics, heat transfer, mass transport, and phase transitions, while accounting for key operational variables such as RPB rotational speed and the mass flow rates of both liquid and gas. The CFD model has been validated using experimental data, specifically by comparing predicted outflow gas and liquid temperatures to measured results. Our findings demonstrate the significant effects of RPB rotational speed and mass flow rates on cooling performance, providing valuable insights for optimizing DCC efficiency in industrial applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2025

Real-time gross actinide monitoring in extreme (10 6 :1) beta-gamma background fluid streams using mini-ATMFDs

Real-time, trace alpha activity monitoring in extreme mixed radiation fields is a well-known challenge in fission product-rich streams containing mixtures of beta-gamma-alpha radiation-emitting isotopes, arising due to the extreme 10 6 :1 ratio of background beta-gamma vs. alpha radiation activity levels. This paper presents results from a novel approach to addressing this challenge, based on Tensioned Metastable Fluid Detector (TMFD) sensing technology. An acoustically driven technique was developed and tested at Purdue University (using dissolved Rn-Po and U-chain alpha emitters), and validated (with dissolved Pu) at Savannah River National Laboratory for real-time monitoring of alpha activity in the 100–10,000 dpm/mL range.

Environmental Monitoring↗

DOE OSTI2021

Film Cooling Performance Predictions For Air And Supercritical CO 2

The current study explores the possibility of cooling the vanes and blades of a direct-fired sCO 2 turbine using film cooling geometries. The operating conditions of a direct-fired sCO 2 cycle and thermophysical properties of the fluid at those conditions can alter the flow field characteristics of the coolant jet and its mixing with the mainstream. Very little information is present in the literature regarding the performance of film cooling geometries employing a supercritical fluid. The objective of this study is to estimate the resulting film cooling effectiveness while also capturing the effects of the crossflow-to-mainstream velocity ratio on the coolant jet. A computational fluid dynamic model is used to study the coolant jet exiting a cylindrical hole located on a flat plate, with the coolant fed by an internal channel. Steady-state RANS equations were solved along with the SST kω model to provide the turbulence closure. The operating conditions for the direct-fired sCO 2 turbine are obtained using an in-house Cooled Turbine Model. Numerical predictions revealed that the crossflow effects and jet lift-off were more pronounced in the case of sCO 2 when compared to air. Spatial distribution of flow field and cooling effectiveness will be presented at different operating conditions. Recommendations for future studies will be given wherever applicable.

42 ENGINEERING↗

DOE OSTI2023

Subsurface Characterization for Evaluating Geothermal Resource Potential from Existing Oil and Gas Wells in Tuttle, Oklahoma: Preprint

Oil and gas (O&G) wells often encounter co-produced hot water, possibly suitable for geothermal direct-use applications. The City of Tuttle is located on the eastern part of the Anadarko sedimentary basin in Oklahoma with high heat-in-place potential and recovery capability at depth. This study aims at demonstrating the potential of geothermal energy production for direct-use applications in two public schools and 250 nearby houses in Tuttle via repurposing existing O&G wells. In this scope, geochemistry, geology, and borehole log data were collected and incorporated into a 3D conceptual subsurface model. A digital elevation model (DEM) was used to represent the study area topography with four O&G wells. In addition, hydrogeochemical characteristics of the geothermal fluid and scaling potential were analyzed using ternary diagrams and chemical ratios to develop mixing models. The subsurface geology model indicated that the study area primarily consists of Permian to Mississippian Sandstone and Limestone formations, implying a porosity ranging between 12% and 22%, and a permeability up to 3.90E-14 m2 in certain reservoir levels. The reservoir temperature is expected to be ranging between 80 degrees C to 95 degrees C around 3 km depth with an average temperature gradient of 22.8 degrees C/km. Chemical geothermometers also estimated the reservoir temperature as 90 degrees C. Findings of the chemical model demonstrated that the geothermal fluid is Sodium-Potassium-Chloride-Sulfate type and possibly mixed with shallow groundwater resulting in higher Ca and Mg concentrations and lower Na/K ratio implying lower calcite scaling. These results comprehensively characterize the potential of geothermal resources in the study area and imply that geothermal energy production by repurposing existing O&G wells is suitable for low-temperature direct-use applications.

GEOTHERMAL ENERGY↗

DOE OSTI2023

Subsurface Characterization for Evaluating Geothermal Resource Potential from Existing Oil and Gas Wells in Tuttle, Oklahoma

gas wells↗

DOE OSTI2019

Implications of real-gas behavior on refractive index calculations for optical diagnostics of fuel–air mixing at high pressures

Three models to compute the refractive index of gaseous mixtures at real gas conditions are presented with the purpose to improve the accuracy of state relationships between refractive index and thermodynamic properties. Models are compared with experimental data to determine one that is applicable to high pressure mixtures with non-ideal thermodynamic behavior near or above a fluids’ critical point. The optimal model is applied to analyze adiabatic thermal mixing of fuel and air at typical diesel engine conditions. Our results show that the ideal gas mixture law is appropriate in the vapor region, assuming it is known, for example, from experimental measurements. Finally, the model is applied for binary fuel-air mixing at supercritical conditions to demonstrate its full potential.

02 PETROLEUM↗

DOE OSTI2022

An Open-Source Framework for the Computational Analysis and Design of Autothermal Chemical Processes

The main project object was to develop software tools for simulating non-equilibrium autothermal processes, improving the prospects for identifying and designing such systems. The project demonstrates the use of these tools to simulate autothermal pyrolysis, a process recently developed at the pilot-plant scale at Iowa State University. In such process, instead of externally heating a reactor to pro-vide the enthalpy of pyrolysis, sufficient oxygen in the form of air is introduced into the reactor to support partial oxidation of reactants and products with the exothermic energy released supporting endothermic pyrolysis reactions. A fluidized bed is used to assure good mixing of biomass and oxidant and provide an isothermal reaction environment. The amount of oxygen required depends upon the kind of biomass being pyrolyzed and parasitic heat losses from the reactor. For example, for woody biomass pyrolyzed under conditions that simulate adiabatic operation, equivalence ratios can be as low as 0.06, compared to 0.20 or higher for autothermal gasifiers. By removing the heat transfer bottleneck of conventional pyrolysis, operation in autothermal mode allowed a significant increase in reactor throughput process, approaching five times the throughput of the conventionally operated pyrolyzer. Different simulation strategies were considered and developed: a zero-dimensional chemistry model was used to verify the applicability of kinetic schemes to predict biomass fast pyrolysis in autothermal conditions. Conventional chemical reactor models such as the plug flow reactor and the partially stirred reactor were used to investigate the role of mixing in the fluidized bed pyrolyzer and to establish the impact of mixing time on the gas-phase reactions. A comprehensive multiphase computational fluid dynamics (mCFD) framework, including polydisperse granular phase modeling and detailed chemical kinetics was formulated and used to model the experimental setup for autothermal biomass fast pyrolysis at ISU. Multiphase CFD was also used to investigate the role of biomass feed positioning on the mixing of biomass in the pyrolizer. Finally, a reduced order model (ROM), suitable to be implemented in process simulators was obtained. Both the mCFD and the ROM were validated against experiments.

09 BIOMASS FUELS↗

DOE OSTI2020

Scaling relationships for the elastic moduli and viscosity of mixed lipid membranes

The elastic and viscous properties of biological membranes play a vital role in controlling cell functions that require local reorganization of the membrane components as well as dramatic shape changes such as endocytosis, vesicular trafficking, and cell division. These properties are widely acknowledged to depend on the unique composition of lipids within the membrane, yet the effects of lipid mixing on the membrane biophysical properties remain poorly understood. Here, we present a comprehensive characterization of the structural, elastic, and viscous properties of fluid membranes composed of binary mixtures of lipids with different tail lengths. We show that the mixed lipid membrane properties are not simply additive quantities of the single-component analogs. Instead, the mixed membranes are more dynamic than either of their constituents, quantified as a decrease in their bending modulus, area compressibility modulus, and viscosity. While the enhanced dynamics are seemingly unexpected, we show that the measured moduli and viscosity for both the mixed and single-component bilayers all scale with the area per lipid and collapse onto respective master curves. This scaling links the increase in dynamics to mixing-induced changes in the lipid packing and membrane structure. More importantly, the results show that the membrane properties can be manipulated through lipid composition the same way bimodal blends of surfactants, liquid crystals, and polymers are used to engineer the mechanical properties of soft materials, with broad implications for understanding how lipid diversity relates to biomembrane function.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Detailed Simulations of a Three-Stage Supercritical Carbon Dioxide Axial Compressor with a Focus on the Shrouded Stator Cavities Flow

This paper describes the findings of detailed simulations performed to investigate the impact of seal teeth cavity leakage flow on the aerodynamic and thermal performance of a three-stage supercritical CO 2 axial compressor. The study compares a shrouded stator configuration (with cavities) to a cantilevered stator configuration (without cavities) to highlight their differences. High-fidelity computational fluid dynamics simulations were performed using non-linear harmonic (NLH) and mixing plane assumptions, considering various possible rotor/stator interface configurations for mixing plane calculations. The key performance parameters for each case were compared, and the best-performing configuration selected for further analysis. The individual stage performance parameters are also examined and compared between the cantilevered and shrouded configurations. It was observed that in the shrouded case, the leakage flow enters the cavity downstream of the stator trailing edge and gets entrained into the primary flow upstream of the stator, leading to boundary layer changes at the hub and degradation of stator and downstream rotor performance. Vortical flow structures were also observed in the stator wells, which tended to change the flow angles around the region, thereby affecting mixing and velocity distribution, resulting in a slight deterioration of compressor performance. Additionally, the study examines windage heating due to shear work from rotating walls, including the seal teeth surface. The amount of shear work done on the leakage flow and the corresponding rise in fluid temperature were quantified, tabulated, and further compared with a simple analytical model, showing good agreement between them and, hence, validating the numerical approach used.

42 ENGINEERING↗

DOE OSTI2021

Entropy generation from hydrodynamic mixing in inertial confinement fusion indirect-drive targets

The increase in entropy from the physical mixing of two adjacent materials in inertial confinement fusion (ICF) implosions and gas-filled hohlraums is analytically assessed. An idealized model of entropy generation from the mixing of identical ideal-gas particles across a material interface in the presence of pressure and temperature gradients is applied. Physically, mix-driven entropy generation refers to the work done by the gases in expanding into a larger common volume from atomic mixing under the condition of no internal energy change, or work needed to restore the initial unmixed state. Furthermore, the effect of a mix-generated entropy increase is analytically shown to lead to less compression of the composite ICF fluid under adiabatic conditions. The amount of entropy generation is estimated to be ~10 J for a Rayleigh–Taylor-induced micrometer-scale annular mixing layer between the solid deuterium–tritium fuel and (undoped) high-density carbon pusher of an imploding capsule at the National Ignition Facility (NIF). This level of entropy generation is consistent with lower-than-expected fuel compressions measured on the NIF [Hurricane et al., Phys. Plasmas 26, 052704 (2019)]. The degree of entropy increase from mixing of high-Z hohlraum wall material and low-Z, moderate- to high-density gas fills is estimated to lead to ~100 kJ of heat generation for NIF-scale experiments [Moody et al., Phys. Plasmas 21, 056317 (2014)]. This value represents a significant fraction of the inferred missing x-ray drive energy based on observed delays in capsule implosion times compared with mainline simulations [Jones et al., Phys. Plasmas 19, 056315 (2012)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Numerical investigation of mixing and heat transfer in a 7.9 m JP-8 pool fire

The response of objects engulfed in, and adjacent to, large-scale pool fires is of interest in accident and safety assessments. In this study, Fuego, a low-Mach turbulent reacting flow code was used to study conjugate heat transfer in a 7.9 m diameter JP-8 pool fire. Simulations were designed to replicate past experimental measurements (Blanchat et al., 2006) of incident heat flux to three cylindrical calorimeters in and around the pool fire. Two turbulent combustion models were compared directly - the eddy dissipation concept and a more recently developed unsteady flamelet model. First and second order spatial and temporal discretization schemes were also compared to assess the performance of low-dissipation numerical operators. Heat flux predictions to the transportation size calorimeter outside the fire were within experimental uncertainties. Inside the fire, experimental measurements were higher than predicted values and may have been a consequence of soot deposition and augmented participating media radiation from soot and fuel vapor. Simulation predictions improved in cases where turbulent kinetic energy and mixing were more resolved. This work, and others referenced herein, suggest that spatial resolution on the order of 0.5–1.0 cm may be required to fully resolve fluid instabilities, vortex production between fire plumes and crosswind, soot production, and fuel-air mixing. This presents a substantial computational challenge for safety assessments of engulfed objects in fully turbulent pool fires.

Conjugate heat transfer↗