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 397 records · Page 22

Electrostatic Relativistic Fluid Models of Electron Emission in a Warm Diode

A semi-analytic fluid model has been developed for characterizing relativistic electron emission across a warm diode gap. Here we demonstrate the use of this model in (i) verifying multi-fluid codes in modeling compressible relativistic electron flows (the EMPIRE-Fluid code is used as an example; see also Ref. 1), (ii) elucidating key physics mechanisms characterizing the influence of compressibility and relativistic injection speed of the electron flow, and (iii) characterizing the regimes over which a fluid model recovers physically reasonable solutions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A sticky situation: The heterogeneity and implications of magma mush unlocking

Magma, partially molten rock below Earth's surface, can be made up of both fluid (melt) and particles (crystals). The relative fraction of the fluid and particles matters for the magma's mobility and ability to reach Earth's surface. Here we conduct physical laboratory experiments to examine the mechanisms through which a magma becomes mobile and the speed at which pure melt and infiltrate a mixture of fluids and particles. The results have implications for understanding magma mobility and the mechanisms that move magma to Earth's surface. Poster presented at the 2025 International Association of Volcanology and Chemistry of Earth's Interior Meeting (IAVCEI).

multiphase flows↗

Review of Wave Energy Converter Power Take-Off Systems, Testing Practices, and Evaluation Metrics: Preprint

While the field of wave energy has been the subject of numerical simulation, scale model testing, and precommercial project testing for decades, wave energy technologies remain in the early stages of development and must continuing proving themselves as a promising modern renewable energy field. One of the difficulties that wave energy systems have been struggling to overcome is the design of highly efficient energy conversion systems that can convert the mechanical power, derived from the oscillation of wave activated bodies, into another useful product. Often the power take-off (PTO) is defined as the single unit responsible for converting mechanical power into another usable form such as electricity, pressurized fluid, compressed air, and others. The PTO, and the entire power conversion chain (PCC), is of great importance as it affects not only how efficient wave power is converted into electricity, but also contributes to the mass, size, structural dynamics, and levelized cost of energy (LCOE) of the wave energy converter (WEC). Unlike wind and solar, there is no industrial standard device, or devices, for wave energy conversion and this diversity is transferred to the PTO system. The majority of current WEC PTO systems incorporate a mechanical or hydraulic drive train, power generator, and an electrical control system. The challenge of WEC PTO designs is designing a mechanical-to-electrical component that can efficiently convert irregular, bi-directional, low frequency and low alternating velocity wave motions. While gross average power levels can be predicted in advance, the variable wave elevation input has to be converted into smooth electrical output and hence usually necessitates some type of energy storage system, such as battery storage, accumulator super capacitors, etc., or other means of compensation such as an array of devices. One of the primary challenges for wave energy converter systems is the fluctuating nature of wave resources, which require WEC components to be designed to handle loads (i.e. torques, forces, and powers) that are many times greater than the average load. This approach requires a much greater PTO capacity than the average power output and indicates a higher cost. In addition, supporting mechanical coupling and or gearing can be added to the PCC to help alleviate the difficulties with transmission and control of fluctuating large loads with low frequencies (indicative of wave forcing) into smaller loads at higher frequencies (optimum for conventional electrical machine design) can quickly increase the complexity of the PCC which could result in a greater number of failure modes and increased maintenance costs. All of the previous points demonstrate how the PTO influences WEC dynamics, reliability, performance and cost which are critical design factors. This paper further explores these topics by providing a review of the state-of-the-art PTO systems currently under development, how these novel PTO systems are tested and derisked prior to precommercial deployment, and the evaluation metrics historically used to differentiate between PTO designs and how they can be improved to support control co-design focused development of wave energy systems.

laboratory testing↗

Unveiling the Essential Parameters Driving Mineral Reactions during CO2 Storage in Carbonate Aquifers through Proxy Models

Numerical simulation is a commonly employed technique for studying carbon dioxide (CO2) storage processes in porous media, particularly saline aquifers. It enables the representation of diverse trapping mechanisms and the assessment of CO2 retention capacity within the subsurface. The intricate physicochemical phenomena involved necessitate the incorporation of multiphase flow, accurate depiction of fluid and rock properties, and their interactions. Among these factors, geochemical reaction rates and mechanisms are pivotal for successful CO2 trapping in carbonate reactive rocks. However, research on kinetic parameters and the influence of lithology on CO2 storage remains limited. This limitation is partly due to the challenges faced in laboratory experiments, where the time scale of the reactions and the lack of in situ conditions hinder accurate measurement of mineral reaction rates. This study employs proxy models constructed using response surfaces calibrated with simulation results to address uncertainties associated with geochemical reactions. Monte Carlo simulation is utilized to explore a broader range of parameters and identify influential factors affecting CO2 mineralization. The findings indicate that an open database containing kinetic parameters can support uncertainty assessment. Additionally, the proxy models effectively represent objective functions related to CO2 injectivity and mineralization, with calcite dissolution playing a predominant role. pH, calcite concentration, and CO2 injection rate significantly impact dolomite precipitation, while quartz content remains unaffected.

58 GEOSCIENCES↗

Simulation of Multiphase Flow and Poromechanical Effects Around Injection Wells in CO 2 Storage Sites

In geological CO 2 storage operations, wellbore deformations and leakage pathways formations can occur around injection and abandoned wells subjected to high rates and long-term CO 2 injection. To guide engineering design and prevent CO 2 leakage risks, a full understanding of the underlying physics and robust numerical models is necessary to evaluate the response of underground formations in the near wellbore region and in the reservoir. In this study, a multi-scale and multi-physics open-source simulator (GEOS) is used to simulate multiphase flow and poromechanical deformations over time in three dimensions. The governing equations for mechanical deformations of the rock body and multiphase compositional fluid flow within the rock matrix are solved with a fully coupled finite element and finite volume approach. The Drucker–Prager model with friction hardening is applied to simulate elastoplastic deformation and a multiphase fluid model with power-law correlations for relative permeability is used to model the migration of CO 2 plume, which are coupled with numerical implicit scheme. Simulation results are verified against multiple analytical solutions for multiphase flow and wellbore problems, thus demonstrating the accuracy of this advanced simulator. In two engineering applications, here we highlight the impact of elastoplastic deformation and coupled modeling for assessing induced displacements and stress perturbations, which are more pronounced in the near wellbore regions. This work focuses on short-term processes in the vicinity of injection wells where stress evolutions, rock deformations and multiphase compositional flow and transport are simulated jointly to ensure wellbore stability and prevent damage. This fully coupled geomechanical model can simulate multiphase flow and any associated poromechanical effects within the CO 2 storage site and in the surrounding formations. Such a large-scale, long-term, multi-physics simulation model is useful in many ways: it can guide operational decisions for CO 2 injection, assess the containment potential and risks of a site, and analyze the wellbore stability and integrity during and after CO 2 injection.

58 GEOSCIENCES↗

Mitigating keyhole pore formation by nanoparticles during laser powder bed fusion additive manufacturing

Keyhole pore formation is one of the most detrimental subsurface defects in the laser metal additive manufacturing process. However, effective ways to mitigate keyhole pore formation beyond tuning laser processing conditions during keyhole mode laser melting are still lacking. Here we report a novel approach to mitigate keyhole pore formation during laser powder bed fusion (LPBF) process by using stable nanoparticles. The critical keyhole depth for keyhole pore generation (i.e., the largest keyhole depth without keyhole pore formation) during LPBF of Al6061 increases from 246 µm to 454 µm (85% increase) after adding TiC nanoparticles. In-depth x-ray imaging studies and thermo-fluid dynamics simulation enable us to identify that two mechanisms work together to mitigate keyhole pore generation: (1) adding nanoparticles prevents the keyhole from collapsing by increasing the liquid viscosity to impede the protrusion development; (2) adding nanoparticles slows down the keyhole pore movement by increasing the liquid viscosity, resulting in the recapturing of the pore by the keyhole. We further demonstrate that adding TiC nanoparticles can also eliminate the keyhole fluctuation induced keyhole pore during LPBF of Al6061. Our research provides a potential way to mitigate keyhole pore formation for defect lean metal additive manufacturing.

36 MATERIALS SCIENCE↗

In situ quantification of fracture slip induced by hydraulic injections in a deep borehole: A comparison of two different borehole techniques

In situ measurements of fracture deformation during fluid injection are rare, yet essential for understanding the mechanical response of fractured rock. In this study, we evaluate the reliability of two methods by comparing their slip vector estimates: high-resolution borehole acoustic televiewer images captured before and after injection tests, and displacement data from a three-component borehole deformation probe recorded during the injections. Acoustic televiewer images capture only final in-plane displacement, whereas three-component borehole deformation measurements provide full 3D, transient fracture movement. Four injection tests in a fractured granitic rock mass along an inclined borehole at the Bedretto Underground Laboratory (Switzerland), beneath more than 1100 m of overburden, were analyzed. The two methods yielded consistent kinematics and comparable slip magnitudes, typically in the range of 0.2 – 0.6 mm. Angular differences between estimated slip directions ranged from 0.2° to 28°. For two test intervals with pre-existing fractures optimally oriented for slip, the measured slip vectors agreed with forward-modeled slip vectors based on the regional in situ stress field, supporting the Wallace–Bott hypothesis. Measurable shear displacement was identified in two test intervals containing fractures oriented nearly perpendicular to the least principal stress, including a newly formed hydraulic fracture. This demonstrates that even fractures unfavorably oriented for slip can undergo shear due to local stress perturbations around an inclined borehole.

Acoustic televiewer↗

Transition state stabilization depends on solvent identity, pore size, and hydrophilicity for epoxidations in zeolites

Ti-silicates activate H 2 O 2 to form Ti-hydroperoxo and Ti-peroxo intermediates that can react with alkenes to form epoxide products. Comparisons of kinetics for 1-octene epoxidation with H 2 O 2 on Ti-BEA and Ti-MFI catalysts with different hydrophilicities in methanol (CH 3 OH) or acetonitrile (CH 3 CN) solvents show the significance of the solvent for stabilizing catalytically-relevant species and the complex interdependencies between solvent, catalyst topology, and hydrophilicity. Epoxidation turnover rates are higher in CH 3 CN than CH 3 OH for Ti-BEA, but the opposite trend is observed for Ti-MFI. Ti-silicates with greater silanol densities, however, give greater epoxidation turnover rates than their hydrophobic counterparts in both solvents. Kinetic, spectroscopic, and thermodynamic analyses show that differences in turnover rates mainly arise from changes in the stabilization of reactive surface species by solvent mediated interactions, because the mechanism of the reaction and stability of the fluid-phase reactants remain similar in CH 3 CN and CH 3 OH. Specifically, apparent activation free energy values ($ΔG^{‡}_{App}$) indicate that surface intermediates responsible for alkene epoxidation are stabilized to a greater extent in CH 3 CN on Ti-BEA and in CH 3 OH on Ti-MFI. Hydrophilic Ti-silicates present lower ($ΔG^{‡}_{App}$) values regardless of solvent identity, which suggests that these differences correspond to the number of hydrogen-bonding solvent molecules found near reactive species bound to Ti active sites. Finally, taken together, these findings demonstrate the role of solvent molecules in allowing reactive intermediates to recognize the properties of active sites beyond the length-scale of covalent bonds, which carry implications for epoxidation but also other reactions within solvent-filled pores of microporous materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A comparative study of machine learning models for predicting the state of reactive mixing

Mixing phenomena are important mechanisms controlling flow, species transport, and reaction processes in fluids and porous media. Accurate predictions of reactive mixing are critical for many Earth and environmental science problems such as contaminant fate and remediation, macroalgae growth, and plankton biomass evolution. Here, to investigate the evolution of mixing dynamics under different scenarios (e.g., anisotropy, fluctuating velocity fields), a finite-element-based numerical model was built to solve the fast, irreversible bimolecular reaction-diffusion equations to simulate a range of reactive-mixing scenarios. A total of 2,315 simulations were performed using different sets of model input parameters comprising various spatial scales of vortex structures in the velocity field, time-scales associated with velocity oscillations, the perturbation parameter for the vortex-based velocity, anisotropic dispersion contrast (i.e., ratio of longitudinal-to-transverse dispersion), and molecular diffusion. The outputs comprised concentration profiles of reactants and products. The inputs to and outputs from these simulations were concatenated into feature and label matrices, respectively, to train 20 different machine learning (ML) models intended to emulate system behavior. These 20 ML emulators, based on linear methods, Bayesian methods, ensemble learning methods, and multilayer perceptrons (MLPs), were trained to classify the state of mixing and predict three quantities of interest (QoIs) characterizing species production, decay (i.e., average concentration, square of average concentration), and degree of mixing (i.e., variances of species concentration). Unsurprisingly, linear classifiers and regressors failed to reproduce the QoIs; however, ensemble methods (classifiers and regressors) and the MLP model accurately classified the state of reactive mixing and the QoIs. Among ensemble methods, random forest and decision-tree-based AdaBoost faithfully predicted the QoIs. At run time, trained ML emulators produced results times faster than the finite-element simulations. Due to their low computational expense and high accuracy, ensemble and MLP models are excellent emulators for these numerical simulations and great utilities in uncertainty quantification exercises, which can require 1,000s of forward model runs.

97 MATHEMATICS AND COMPUTING↗

OVERVIEW OF PARTICLE DEPOSITION MODELS FOR SPENT NUCLEAR FUEL STORAGE SYSTEMS

Deposition models were built to evaluate contaminant deposition on spent nuclear fuel (SNF) canisters. The primary contaminant of concern is chloride, which is dispersed in the atmosphere and then deposits onto the canisters. During dry storage, the primary degradation process is likely to be Chloride Induced Stress Corrosion Cracking (CISCC) at the heat-affected zones of the canister welds. It is known that stainless steel canisters are susceptible to CISCC; however, the rate of chloride deposition onto the canisters is poorly known, based on sparse field data from a small number of sites. The models presented in this study could be useful for determining the rate of deposition on the canisters and the likelihood of CISCC to help with SNF canister ageing management. The deposition models were developed with the commercial computational fluid dynamics (CFD) code STAR-CCM+. Various deposition mechanisms were considered and incorporated into the models, and a sensitivity study was conducted to determine the most important mechanisms for deposition within a SNF storage system. The models included both a vertical and horizontal configuration storage system: NAC International’s Modular, Advanced Generation, Nuclear All-purpose STORage System (MAGNASTOR®) and a NUHOMS® horizontal storage module respectively. The resulting canister deposition on the horizontal canister is visually compared with inspection data taken onsite at the Calvert Cliffs Nuclear Power Plant. These models are preliminary, and development of the models will continue. Future validation exercises are currently being planned, including the Canister Deposition Field Demonstration (CDFD) effort funded by U.S. Department of Energy office of Nuclear Energy office of Spent Fuel Waste Science and Technology. The goal of the modeling presented is to demonstrate a potential modeling technique that could be used to plan and inform SNF canister ageing management programs with predictive models for the timing and occurrence of canister CISCC.

Suffield, Sarah R.↗

Finite Element Modeling of Electrochemical Polishing of Niobium in Hydrofluoric-Sulfuric Acid Electrolyte

Niobium (Nb) used in superconducting radio-frequency cavities requires smooth surface to achieve optimal performance. In this work, a finite element model that coupled electrochemistry, heat transfer, and fluid dynamics was developed to investigate the electrochemical polishing mechanisms of Nb, using experimentally measured polarization results of coupon samples as validations. The current and potential distribution, oxide growth kinetics of Nb in a complex cavity geometry was investigated as a function of temperature and coolant flow. A low temperature coolant with intermediate flow rate was found to reduce surface current and ensure oxide uniformity. These results could shed light on the design of future particle accelerators.

SRF accelerator cavity↗

The Alamo multiphysics solver for phase field simulations with strong-form mechanics and block structured adaptive mesh refinement

Alamo is a high-performance scientific code that uses block-structured adaptive mesh refinement to solve such problems as: the ignition and burn of solid rocket propellant, plasticity, damage and fracture in materials undergoing loading, and the interaction of compressible flow with eroding solid materials. Alamo is powered by AMReX, and provides a set of unique methods, models, and algorithms that enable it to solve solid-mechanics problems (coupled to other physical behavior such as fluid flow or thermal diffusion) using the power of block-structured adaptive mesh refinement.

36 MATERIALS SCIENCE↗

Computational Fluid Dynamics Modeling to Simulate a Combined Reforming Process for Syngas and Hydrogen Production

An Oxygen Transport Membrane (OTM) combined reforming technology for producing syngas and hydrogen integrates the advantages of multiple processes—steam methane reforming (SMR), autothermal reforming (ATR), an air separation unit (ASU)—into a single integrated technology. The OTM consists of a primary reforming tube, in which desulfurized natural gas is partially reformed by steam at high pressure in the presence of a metal catalyst. This process is followed in series by a ceramic OTM with a secondary reformer, in which residual methane reforms and O 2 - ions react with a portion of the CO and H 2 fuel to provide the heat to support both primary and secondary reforming. Although the OTM combined reformer technology for syngas and H 2 production has been substantially developed in the last decade, several challenges that affect the overall production efficiency and reliability are yet to be fully understood, addressed, and resolved. Therefore, developing Computational Fluid Dynamics (CFD) models that incorporate fluid dynamics, mass transport, kinetics, heat transport, and structural mechanics is critical to understanding and minimizing the probability of tube failures during the startup and operation. In this report, an exhaustive literature review was performed to survey the current state of technology for producing syngas and H 2 using either conventional or renewable energy sources. The feedstocks reviewed include natural gas and coal for the conventional technologies, whereas biomass, solar, wind, and nuclear energy for the renewable technologies. The existing industry-grade COMSOL multiphysics models of OTM were upgraded for the latest software release. In addition, they were improved to help achieve grid and solver independence and were successfully ported on the ORNL high-performance computing clusters to speed up their run times. A 42% reduction in the simulation run time was achieved. A new higher-fidelity CFD model of an OTM tube was developed in the StarCCM+ simulation platform. This new model was designed to simulate various physics using first principles, e.g., turbulent flow, heat transfer, and chemical reactions while avoiding unnecessary simplifications. The resulting predictions were qualitatively assessed and provided useful insights into the multiphysics complexity of an OTM tube.

08 HYDROGEN↗

Self-Healing Cements with Improved Toughness at Casing and Formation Interfaces for Subsurface Applications (CRADA 530)

Over the three-year CRADA, Pacific Northwest National Laboratory (PNNL) advanced a new class of Molecular Velcro™ self-healing wellbore cements designed to improve the long-term integrity of wells used for geothermal energy production and carbon storage. Conventional wellbore cement can crack due to repeated mechanical loading, temperature swings, and exposure to aggressive subsurface fluids. Once cracks form, they can create leakage pathways, drive expensive remediation (“workovers”), and in extreme cases shorten well life. This project addressed that problem by developing polymer-modified cement systems that can autonomously re-seal cracks and better tolerate the extreme conditions experienced in subsurface wells, while remaining compatible with practical cementing operations.

15 GEOTHERMAL ENERGY↗

A gas dynamic perspective on particle lifting in electrostatic discharge-like devices using multiphase particle-in-cell approach

Although particle–laden electrostatic discharges are widely used in laboratories as well as in industrial applications, the mechanism of particle lifting for particles initially at rest in such highly unsteady systems is not well understood. A multiphase gas–particle solver is developed using the multiphase particle-in-cell (MP-PIC) approach to emulate the interaction of a compressible shock-dominated gas phase with the dense particle phase. First, the two-way coupled solver is initially used to study the interaction of a planar traveling shock with a vertical curtain of particulates. The gas and particle phase evolution was found to be in good agreement with a similar experimental study in Ling et al. [Phys. Fluids 24, 113301 (2012)]. Second, the MP-PIC code is used to study the interaction of an expanding blast wave with a thick bed of particles. The simulation considered forces such as quasi-steady drag, pressure-gradient, added-mass, Saffman, and Magnus forces. We observe that the vertical liftoff particles close to the shock impingement point in this configuration are associated with the quasi-steady drag, pressure gradient, and added-mass forces. Also, the Saffman lift and Magnus forces contribute to lifting particles located radially farther away from the shock impingement point. In addition, the study finds a decrease in particle lifting efficiency with decreasing plasma kernel length and shock strength.

42 ENGINEERING↗

Higher-order nonequilibrium term: Effective power density quantifying evolution towards or away from local thermodynamic equilibrium

A common approach to assess the nature of energy conversion in a classical fluid or plasma is to compare power densities of the various possible energy conversion mechanisms. A leading research area is quantifying energy conversion for systems that are not in local thermodynamic equilibrium (LTE), as is common in a number of fluid and plasma systems. Here we introduce the “higher-order nonequilibrium term” (HORNET) effective power density, which quantifies the rate of change of departure of a phase space density from LTE. It has dimensions of power density, which allows for quantitative comparisons with standard power densities. We employ particle-in-cell simulations to calculate HORNET during two processes, magnetic reconnection and decaying kinetic turbulence in collisionless magnetized plasmas, that inherently produce non-LTE effects. We investigate the spatial variation of HORNET and the time evolution of its spatial average. By comparing HORNET with power densities describing changes to the internal energy (pressure dilatation, Pi − D , and divergence of the vector heat flux density), we find that HORNET can be a significant fraction of these other measures (8% and 35% for electrons and ions, respectively, for reconnection; up to 67% for both electrons and ions for turbulence), meaning evolution of the system towards or away from LTE can be dynamically important. Applications to numerous plasma phenomena are discussed. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonadiabatic Energy Fluctuations of Scale-Invariant Quantum Systems in a Time-Dependent Trap

We consider the nonadiabatic energy fluctuations of a many-body system in a time-dependent harmonic trap. In the presence of scale-invariance, the dynamics becomes self-similar and the nondiabatic energy fluctuations can be found in terms of the initial expectation values of the second moments of the Hamiltonian, square position, and squeezing operators. Nonadiabatic features are expressed in terms of the scaling factor governing the size of the atomic cloud, which can be extracted from time-of-flight images. We apply this exact relation to a number of examples: the single-particle harmonic oscillator, the one-dimensional Calogero-Sutherland model, describing bosons with inverse-square interactions that includes the non-interacting Bose gas and the Tonks-Girdardeau gas as limiting cases, and the unitary Fermi gas. We illustrate these results for various expansion protocols involving sudden quenches of the trap frequency, linear ramps and shortcuts to adiabaticity. Our results pave the way to the experimental study of nonadiabatic energy fluctuations in driven quantum fluids.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Planar turbulent wakes under pressure gradient: Integral and self-similarity analyses

By using a combination of integral and self-similarity analyses, the generalized analytical solutions for the mean transverse velocity and Reynolds shear stress are rigorously derived for the first time for the far field of planar turbulent wakes under arbitrary pressure gradients. Specifically, by assuming self-similarity for the mean axial velocity, the analytical formulation for the mean transverse velocity is obtained from the integral of the mean continuity equation, and the analytical formulation for the Reynolds shear stress is obtained from the integral of the momentum equation. The generalized analytical formulations for the mean transverse velocity and Reynolds shear stress consist of multiple components, each with its unique scale and physical mechanism. In the zero pressure gradient limit, the generalized formulations recover the single-scale equations reported by Wei, Liu, and Livescu. Furthermore, simpler approximate formulations for the mean transverse velocity and Reynolds shear stress are also obtained, and show excellent agreement with the experimental measurements. Here, the findings provide new insights into the properties of planar turbulent wakes under pressure gradients, filling some long-standing gaps in the existing literature.

42 ENGINEERING↗