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 1,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Mechanical Resonant Sensing of Spin Texture Dynamics in a 2D Antiferromagnet

Abstract The coupling between the spin degrees of freedom and macroscopic mechanical motions, including striction, shearing, and rotation, has attracted wide interest with applications in actuation, transduction, and information processing. Experiments so far have established the mechanical responses to the long‐range ordered or isolated single spin states. However, it remains elusive whether mechanical motions can couple to a different type of magnetic structure, the non‐collinear spin textures, which exhibit nanoscale spatial variations of spin (domain walls, skyrmions,etc.) and are promising candidates to realize high‐speed computing devices. Here, collective spin texture dynamics is detected with nanoelectromechanical resonators fabricated from 2D antiferromagnetic (AFM) MnPS 3 with 10 −9 strain sensitivity. By examining radio frequency mechanical oscillations under magnetic fields, new magnetic transitions are identified with sharp dips in resonant frequency. They are attributed to collective AFM domain wall motions as supported by the analytical modeling of magnetostriction and large‐scale spin‐dynamics simulations. Additionally, an abnormally large modulation in the mechanical nonlinearity at the transition field infers a fluid‐like response due to ultrafast domain motion. The work establishes a strong coupling between spin texture and mechanical dynamics, laying the foundation for electromechanical manipulation of spin texture and developing quantum hybrid devices.

Chemistry↗
DOE OSTI2018

Theory of ITG turbulent saturation in stellarators: identifying mechanisms to reduce turbulent transport

A three-field fluid model that allows for general three-dimensional equilibrium geometry is developed to describe ion temperature gradient turbulent saturation processes in stellarators. The theory relies on the paradigm of nonlinear transfer of energy from unstable to damped modes at comparable wavelength as the dominant saturation mechanism. The unstable-to-damped mode interaction is enabled by a third mode that for dominant energy transfer channels primarily serves as a regulator of the nonlinear energy transfer rate. The identity of the third wave in the interaction defines different scenarios for turbulent saturation with the dominant scenario depending upon the properties of the 3D geometry. The nonlinear energy transfer physics is quantified by the product of a turbulent correlation lifetime and a geometric coupling coefficient. The turbulent correlation time is determined by a three-wave frequency mismatch, which at long wavelength can be calculated from the sum of the linear eigenfrequencies of the three modes. Larger turbulent correlation times denote larger levels of nonlinear energy transfer and hence smaller turbulent transport. The theory provides an analytic prediction for how 3D shaping can be tuned to lower turbulent transport through saturation processes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1976

Dissipative, forced turbulence in two-dimensional magnetohydrodynamics

The equations of motion for turbulent two-dimensional magnetohydrodynamic flows are solved in the presence of finite viscosity and resistivity, for the case in which external forces (mechanical and/or magnetic) act on the fluid. The goal is to verify the existence of a magnetohydrodynamic dynamo effect which is represented mathematically by a substantial back-transfer of mean square vector potential to the longest allowed Fourier wavelengths. External forces consisting of a random part plus a fraction of the value at the previous time step are employed, after the manner of Lilly for the Navier-Stokes case. The regime explored is that for which the mechanical and magnetic Reynolds numbers are in the region of 100 to 1000. The conclusions are that mechanical forcing terms alone cannot lead to dynamo action, but that dynamo action can result from either magnetic forcing terms or from both mechanical and magnetic forcing terms simultaneously.

Fyfe, D.↗
DOE OSTI2018

Plasma Energization in Colliding Magnetic Flux Ropes

Magnetic flux ropes are commonly observed throughout the heliosphere, and recent studies suggest that interacting flux ropes are associated with some energetic particle events. In this work, we carry out 2D particle-in-cell (PIC) simulations to study the coalescence of two magnetic flux ropes (or magnetic islands), and the subsequent plasma energization processes. The simulations are initialized with two magnetic islands embedded in a reconnecting current sheet. The two islands collide and eventually merge into a single island. Particles are accelerated during this process as the magnetic energy is released and converted to the plasma energy, including bulk kinetic energy increase by the ideal electric field, and thermal energy increase by the fluid compression and the nonideal electric field. We find that contributions from these different energization mechanisms are all important and comparable with each other. Fluid shear and a nongyrotropic pressure tensor also contribute to the energy conversion process. For simulations with different box sizes ranging from $L_x$ ~ 25–100$d_i$ and ion-to-electron mass ratios $m_i /m_e$ = 25,100, and 400, we find that the general evolution is qualitatively the same for all runs, and the energization depends only weakly on either the system size or the mass ratio. The results may help us understand plasma energization in solar and heliospheric environments.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1953

Study of secondary-flow patterns in an annular cascade of turbine nozzle blades with vortex design

In order to increase understanding of the origin of losses in a turbine, the secondary-flow components in the boundary layers and the blade wakes of an annular cascade of turbine nozzle blades (vortex design) was investigated. A detailed study was made of the total-pressure contours and, particularly, of the inner-wall loss cores downstream of the blades. The inner-wall loss core associated with a blade of the turbine-nozzle cascade is largely the accumulation of low-momentum fluids originating elsewhere in the cascade. This accumulation is effected by a secondary-flow mechanism which acts to transport the low-momentum fluids across the channels on the walls and radially in the blade wakes and boundary layers. The patterns of secondary flow were determined by use of hydrogen sulfide traces, paint, flow fences, and total pressure surveys. At one flow condition investigated, the radial transport of low-momentum fluid in the blade wake and on the suction surface near the trailing edge accounted for 65 percent of the loss core; 30 percent resulted from flow in the thickened boundary layer on the suction surface and 35 percent from flow in the blade wake.

FLOW, SUBSONIC↗
DOE OSTI2021

Investigation of theoretical scaling laws using large eddy simulations for airborne spreading of viral contagion from sneezing and coughing

Using a set of large eddy point-particle simulations, we explore the fluid dynamics of an ejected puff resulting from a cough/sneeze. The ejection contains over 61,000 potentially virus-laden droplets at an injection Reynolds number of about 46,000, comparable to an actual cough/ sneeze. We observe that global puff properties, such as centroid, puff volume, momentum, and buoyancy vary little across realizations. Other properties, such as maximum extent, shape, and edge velocity of the puff, may exhibit substantial variation. In many realizations, a portion of the puff splits off and advances along a random direction, while keeping airborne droplet nuclei afloat. This peeled-off portion provides a mechanism for virus-laden droplets to travel over large distances in a short amount of time. We also observe that the vast majority of droplets remain suspended within the puff after all liquid has evaporated. The main objectives of the study are to (i) evaluate assumptions of Balachandar’s et al. theory [Int. J. Multiphase Flow 132, 103439 (2020)], which include buoyancy effects, shape of the puff, and droplet evaporation rate, (ii) obtain values of closure parameters, which include location and time of the virtual origin, and puff entrainment and drag coefficients, and (iii) evaluate the accuracy of the theory in predicting the shape, size, and location of the puff, as well as droplet number density long after ejection. The theory adequately predicts global puff properties including size, velocity, and distance traveled, the largest size of droplets that exit the puff due to settling, and the droplet size distribution within the puff long after ejection.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Phase-field modeling of rate-dependent fluid-driven fracture initiation and propagation

The rate-dependent behavior associated with deformation and fracturing of materials, such as natural rocks, poses significant challenges for modeling. In addition to the complications of the viscoelastic response, the speed of fracture propagation reflects micromechanical mechanisms in the fracture process zone(FPZ). In order to represent these complicated behaviors, a thermodynamically consistent, rate-dependent fracture model is required. Based on rigorous thermodynamic principles, we derive a rate-dependent phase-field mechanical model coupled with single-phase fluid flow in both the matrix and the fracture. The model is guaranteed to satisfy energy conservation during fracture propagation. Here, the system of equations is solved using the introduced solution procedure and a novel preconditioner that accounts for the complex fluid-structure interaction. The proposed phase-field model is tested against several benchmark problems on solid-fluid coupling, fluid-driven fracture propagation and rate-dependent viscoelastic deformation. The model serves as a strong basis for investigating rate-dependent fracturing experiments and for making predictions of material behaviors under new conditions.

42 ENGINEERING↗
DOE OSTI2021

Evolution of anisotropy with saturation and its implications for the elastoplastic responses of clay rocks

Many clay rocks have distinct bedding planes. Experimental studies have shown that their mechanical properties evolve with the degree of saturation (DOS), often with higher stiffness and strength after drying. For transversely isotropic rocks, the effects of saturation can differ between the bed-normal (BN) and bed-parallel (BP) directions, which gives rise to saturation-dependent stiffness and strength anisotropy. Accurate prediction of the mechanical behavior of clay rocks under partially saturated conditions requires numerical models that can capture the evolving elastic and plastic anisotropy with DOS. In this study, we present an anisotropy framework for coupled solid deformation-fluid flow in unsaturated elastoplastic media. Here we incorporate saturation-dependent strength anisotropy into an anisotropic modified Cam-Clay (MCC) model and consider the evolving anisotropy in both the elastic and plastic responses. The model was calibrated using experimental data from triaxial tests to demonstrate its capability in capturing strength anisotropy at various levels of saturation. Through numerical simulations, we demonstrate the role of evolving stiffness and strength anisotropy in the mechanical behavior of clay rocks. Plane strain simulations of triaxial compression tests were also conducted to demonstrate the impacts of material anisotropy and DOS on the mechanical and fluid flow responses.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Subsurface Characterization and Machine Learning Predictions at Brady Hot Springs Results

Geothermal power plants typically show decreasing heat and power production rates over time. Mitigation strategies include optimizing the management of existing wells - increasing or decreasing the fluid flow rates across the wells - and drilling new wells at appropriate locations. The latter is expensive, time-consuming, and subject to many engineering constraints, but the former is a viable mechanism for periodic adjustment of the available fluid allocations. Data and supporting literature from a study describing a new approach combining reservoir modeling and machine learning to produce models that enable strategies for the mitigation of decreased heat and power production rates over time for geothermal power plants. The computational approach used enables translation of sets of potential flow rates for the active wells into reservoir-wide estimates of produced energy and discovery of optimal flow allocations among the studied sets. In our computational experiments, we utilize collections of simulations for a specific reservoir (which capture subsurface characterization and realize history matching) along with machine learning models that predict temperature and pressure timeseries for production wells. We evaluate this approach using an "open-source" reservoir we have constructed that captures many of the characteristics of Brady Hot Springs, a commercially operational geothermal field in Nevada, USA. Selected results from a reservoir model of Brady Hot Springs itself are presented to show successful application to an existing system. In both cases, energy predictions prove to be highly accurate: all observed prediction errors do not exceed 3.68% for temperatures and 4.75% for pressures. In a cumulative energy estimation, we observe prediction errors that are less than 4.04%. A typical reservoir simulation for Brady Hot Springs completes in approximately 4 hours, whereas our machine learning models yield accurate 20-year predictions for temperatures, pressures, and produced energy in 0.9 seconds. This paper aims to demonstrate how the models and techniques from our study can be applied to achieve rapid exploration of controlled parameters and optimization of other geothermal reservoirs. Includes a synthetic, yet realistic, model of a geothermal reservoir, referred to as open-source reservoir (OSR). OSR is a 10-well (4 injection wells and 6 production wells) system that resembles Brady Hot Springs (a commercially operational geothermal field in Nevada, USA) at a high level but has a number of sufficiently modified characteristics (which renders any possible similarity between specific characteristics like temperatures and pressures as purely random). We study OSR through CMG simulations with a wide range of flow allocation scenarios. Includes a dataset with 101 simulated scenarios that cover the period of time between 2020 and 2040 and a link to the published paper about this project, where we focus on the Machine Learning work for predicting OSR's energy production based on the simulation data, as well as a link to the GitHub repository where we have published the code we have developed (please refer to the repository's readme file to see instructions on how to run the code). Additional links are included to associated work led by the USGS to identify geologic factors associated with well productivity in geothermal fields. Below are the high-level steps for applying the same modeling + ML process to other geothermal reservoirs: 1. Develop a geologic model of the geothermal field. The location of faults, upflow zones, aquifers, etc. need to be accounted for as accurately as possible 2. The geologic model needs to be converted to a reservoir model that can be used in a reservoir simulator, such as, for instance, CMG STARS, TETRAD, or FALCON 3. Using native state modeling, the initial temperature and pressure distributions are evaluated, and they become the initial conditions for dynamic reservoir simulations 4....

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Modeling Subsurface Performance of a Geothermal Reservoir Using Machine Learning

Geothermal power plants typically show decreasing heat and power production rates over time. Mitigation strategies include optimizing the management of existing wells—increasing or decreasing the fluid flow rates across the wells—and drilling new wells at appropriate locations. The latter is expensive, time-consuming, and subject to many engineering constraints, but the former is a viable mechanism for periodic adjustment of the available fluid allocations. In this study, we describe a new approach combining reservoir modeling and machine learning to produce models that enable such a strategy. Our computational approach allows us, first, to translate sets of potential flow rates for the active wells into reservoir-wide estimates of produced energy, and second, to find optimal flow allocations among the studied sets. In our computational experiments, we utilize collections of simulations for a specific reservoir (which capture subsurface characterization and realize history matching) along with machine learning models that predict temperature and pressure timeseries for production wells. We evaluate this approach using an “open-source” reservoir we have constructed that captures many of the characteristics of Brady Hot Springs, a commercially operational geothermal field in Nevada, USA. Selected results from a reservoir model of Brady Hot Springs itself are presented to show successful application to an existing system. In both cases, energy predictions prove to be highly accurate: all observed prediction errors do not exceed 3.68% for temperatures and 4.75% for pressures. In a cumulative energy estimation, we observe prediction errors that are less than 4.04%. A typical reservoir simulation for Brady Hot Springs completes in approximately 4 h, whereas our machine learning models yield accurate 20-year predictions for temperatures, pressures, and produced energy in 0.9 s. This paper aims to demonstrate how the models and techniques from our study can be applied to achieve rapid exploration of controlled parameters and optimization of other geothermal reservoirs.

15 GEOTHERMAL ENERGY↗
DOE OSTI2018

On the Relationship Between Fault Permeability Increases, Induced Stress Perturbation, and the Growth of Aseismic Slip During Fluid Injection

Fluid injections into the deep subsurface can, at times, generate earthquakes, but often, they only produce aseismic deformations. Here we analyze the influence of fault hydromechanical properties on the growth of injection-induced aseismic slip. Using hydromechanical modeling, we show how permeability enhancement in addition to the background stress and frictional weakening has an important effect on the pressure diffusion and slip growth during injection. We find that the more pronounced the fault permeability enhancement, the stronger is the growth of the aseismic slip zone. The effect of enhanced permeability is more pronounced when the fault is initially close to failure. Our results show that aseismic slip grows beyond the pressurized zone when the fault permeability increases, while slip remains behind the pressurized zone when permeability does not vary from its initial preslip value. Furthermore, fault permeability increases should be considered as complementary mechanism to current models of fluid-induced aseismic slip.

58 GEOSCIENCES↗
DOE OSTI2018

Low-GWP HVAC System with Ultra-small Centrifugal Compression

Mechanical Solutions, Incorporated (MSI) was awarded a contract from the United States Department of Energy under Contract DE-EE0007042 to advance the potential of vapor compression utilizing an advanced high speed centrifugal compressor with an advanced heat exchanger design. The design was mandated to maintain high efficiency while utilizing a low GWP refrigerant. The specific goal of this program was to develop an advanced 5 ton air conditioning system. The target of the air conditioner is residential and commercial. To accomplish these goals, MSI teamed with Lennox International, a leading air conditioning manufacturer. MSI would design, oversee manufacturing, assemble, and test the compressor, while Lennox designed the heat exchanger and tested the overall air conditioner. Lennox initiated the project with a detailed cycle analysis that examined systems and refrigerants. At various break points, MSI performed a preliminary compressor design to identify realizable targets. Initially, the refrigerant that was targeted was R1234ze(E), but this was subsequently rejected over concerns of flammability. R515a was ultimately accepted as a qualified fluid based upon its GWP and flammability rating. MSI performed a complete design of the centrifugal compressor based upon the cycle system design and analysis. The design was based upon advanced aerodynamics utilizing transient Computational Fluid Dynamics (CFD), as well as mechanical and rotordynamic analysis. A detailed materials investigation was conducted to ensure mechanical integrity and refrigerant compatibility. To achieve the high rotational speed required by the compressor (105,300 rpm), foil bearings were designed by MSI, as was a high speed motor to provide the required torque for the compressor. All elements of the design complied with Design for Manufacturing and Assembly (DFMA) principles. The compressor was manufactured in components, utilizing specializing subcontractors throughout the country. The components were assembled into an operational compressor in MSI’s Albany facility, where it was tested in the MSI refrigeration test loop. The compressor achieved 90% of its maximum speed during testing at MSI. The interaction between the inverter and the high speed motor was the reason that full speed could not be attained. The compressor operated smoothly and quietly at all speeds, with the bearings performing as expected. The compressor was then shipped to Lennox’s facilities in Texas, where it was tested in the entire air conditioning system. Lennox was able to achieve 95% of the maximum rotational speed. They employed the same inverter as MSI. Lennox compiled the results of the testing and based upon the results, a go forward plan has been developed to correct the issues with the inverter, as well as redesign some parts that proved to be a challenge during the program. The consensus of the team is that a fully functional air conditioning system can be developed to achieve all of the goals identified in this program.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2017

Modeling the influence of nozzle-generated turbulence on diesel sprays

The physical mechanisms governing spray breakup in direct injection engines, such as aerodynamic induced instabilities and nozzle-generated cavitation and turbulence, are not well understood due to the experimental and computational limitations in resolving these processes. Recent x-ray and visible extinction measurements have been con-ducted with a targeted interest in the spray formation region in order to characterize the distribution of droplet sizes throughout the spray. Detailed analysis of these measurements shows promise of yielding insight into likely mechanisms governing atomization, which can inform the improvement of spray models for engine computational fluid dynamic (CFD) codes. In order to investigate potential atomization mechanisms, we employ a joint experimental and computational approach to characterize the structure of the spray formation region using the Engine Combustion Network Spray D injector. X-ray tomography, radiography and ultra-small angle x-ray scattering measurements conducted at the Advanced Photon Source at Argonne National Laboratory quantify the injector geometry, liquid fuel mass and Sauter mean diameter (SMD) distributions under non-vaporizing conditions. Diffused back-illumination imaging measurements, conducted at the Georgia Institute of Technology, characterize the asymmetry of the spray structure. The selected range of injection pressures (50 – 150 MPa) and ambient densities (1.2 – 22.8 kg/m3) allow for the influence of aerodynamic forces on the spray to be studied in a controlled and systematic manner, while isolating the atomization process from the effects of vaporization. In comparison to high ambient density conditions, the spray is observed to be more asymmetric at low ambient density conditions. Although several mechanisms may cause asymmetries in the nozzle exit flow conditions and ultimately the spray distribution, irregularities in the internal nozzle geometry were identified, suggesting an increased sensitivity of the spray structure to internal nozzle surface finish imperfections at such conditions. The presence of these asymmetries may influence the ability to interpret line-of-sight measurements and their derived SMD values and trends from a single viewing angle of the spray. With this consideration in mind, the measured local sensitivities to ambient density suggest that for ambient densities less than 2.4 kg/m3, aerodynamic effects are likely suppressed, allowing the influence of turbulent-induced breakup to be isolated. In concert with the experimental measurements, we utilize three-dimensional, CFD Lagrangian-Eulerian spray simulations in CONVERGE to evaluate the details of the predicted spray structure. In particular, we compare measured and predicted sensitivities of the SMD distribution to changes in injection and ambient conditions from three different atomization models, namely Kelvin Helmholtz (KH), KH Aerodynamics Cavitation Turbulence (KH-ACT), and the newly developed KH-Faeth hybrid model. While none of the existing hybrid spray models were able to replicate the experimentally observed sensitivities, it was found that the scales characterizing the KH-Faeth model show promise of capturing the experimentally observed trends if the effects of secondary droplet breakup are neglected. These results inform recommendations for future experiments and computational studies that can guide the development of an improved spray breakup model.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Embedding condition describing shearing, dissipative collapse

Highlights: • Perturbative approach to the time-dependent Karmarkar condition in shearing collapse. • Analysis of the (in)stability of the collapsing fluid using the TOV equation. • Study of the temperature profiles using a causal heat transport equation. In this work we employ the time-dependent Karmarkar condition to model dissipative gravitational collapse in the presence of shear which leads to a final static configuration. We believe that this is a first attempt at using the time-dependent Karmarkar condition describing a radiating star with shear in which the collapsing core dissipates energy in the form of a radial heat flux. The final static configuration which is described by an embedding class I spacetime was obtained by Hansraj and Moodly who showed that this solution satisfies all the physical requirements for a static compact object. As the collapse proceeds and the heat flux tends to zero the dynamical core evolves into a static compact core with the exterior spacetime being described by the vacuum Schwarzschild solution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗