SEARCH · Search NASA

Results for “hydrodynamic interactions”

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 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Revisiting the empirical particle-fluid coupling model used in DEM-CFD by high-resolution DEM-LBM-IMB simulations: A 2D perspective

The work investigates the applicability of the unresolved Computational Fluid Dynamics and Discrete Element Method (CFDDEM) technique based on empirical equations for fluid-particle coupling. We first carry out a series of representative volume element simulations using the high-resolution particle-resolved Lattice Boltzmann method and Discrete Element Method (LBMDEM) coupled by an Immersed Moving Boundary (IMB) scheme. Then, we compare the results obtained by both LBMDEM and empirical equations used in unresolved CFDDEM with analytical solutions. It is found that the existing empirical equations used in solving fluid-particle interactions in 2D CFDDEM fail to accurately calculate the hydrodynamic force applied to solid particles. The underlying reason is that the existing empirical models are obtained based on 3D experimental results and thus are not applicable to 2D problems. Based on the simulation results, a new drag coefficient model is then proposed. The estimated drag forces using the new model are compared favourably with the simulated ones, indicating the good performance of the proposed model.

42 ENGINEERING↗

DOE OSTI2020

Flow-induced oscillation patterns for two tandem cylinders with turbulence stimulation and variable stiffness and damping

Herein, the oscillation patterns of two 1-DOF cylinder oscillators, undergoing VIV and galloping, are investigated in a free surface water channel for 3.2 × 104 ≤ Re ≤ 1.2 × 10 5 . The cylinders are arranged in tandem and supported by springs for a range of different spring stiffness and damping parameters. The efficiency of a current energy converter (CEC), based on flow-induced oscillation (FIO) of multiple cylinders in tandem, is critically related to the cylinder oscillation patterns. Due to the limited number of studies on the FIO for multiple cylindrical or prismatic bodies, their oscillation patterns have neither been identified nor classified. The surfaces of the cylinders are modified by turbulence stimulation to enhance FIO. Three different center-to-center spacing, five stiffness, and six damping ratios, for a total of 90 sets of experiments were conducted. The current velocity range is from 0.34 m/s to 1.32 m/s. From more than 2000 tests, five major patterns, nine sub-patterns are identified and classified. The patterns are defined based on the amplitudes, frequencies and the phase angle differences between the two cylinders. The characteristics and mechanics of each oscillation pattern are explored and explained from the perspective of fluid-structure interaction (FSI). By systematically varying the parameters, the underlying hydrodynamic mechanisms, including the coupling level between vortices and cylinders, the shielding effect, the wake effect, and the stability states are revealed. Few preliminary observations on the connection between oscillation pattern and harnessed power by the tandem cylinders are reported.

42 ENGINEERING↗

DOE OSTI2021

Simulating dust grain-radiation coupling on a moving mesh

ABSTRACT We present a model for the interaction between dust and radiation fields in the radiation hydrodynamic code arepo-rt, which solves the moment-based radiative transfer equations on an unstructured moving mesh. Dust is directly treated using live simulation particles, each of which represent a population of grains that are coupled to hydrodynamic motion through a drag force. We introduce methods to calculate radiation pressure on and photon absorption by dust grains. By including a direct treatment of dust, we are able to calculate dust opacities and update radiation fields self-consistently based on the local dust distribution. This hybrid scheme coupling dust particles to an unstructured mesh for radiation is validated using several test problems with known analytic solutions, including dust driven via spherically symmetric flux from a constant luminosity source and photon absorption from radiation incident on a thin layer of dust. Our methods are compatible with the multifrequency scheme in arepo-rt, which treats UV, optical photons as single scattered and IR photons as multi scattered. At IR wavelengths, we model heating of and thermal emission from dust. Dust and gas are not assumed to be in local thermodynamic equilibrium but transfer energy through collisional exchange. We estimate dust temperatures by balancing these dust-radiation and dust-gas energy exchange rates. This framework for coupling dust and radiation can be applied in future radiation hydrodynamic simulations of galaxy formation.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Simulations of the Collision between a Debris Stream and an Outer Dusty Torus: A Possible Channel for Forming a Fast-rise and Long-delay Radio Outburst in Tidal Disruption Events

A geometrically thick dusty torus structure is believed to exist in the nuclear regions of galaxies (especially in active galactic nuclei). The debris stream from a tidal disruption event (TDE) will possibly collide with the dusty torus and produce a transient flare. We perform three-dimensional hydrodynamic simulations to model the dynamical evolution of the interaction between unbound debris and a dusty torus. During the continuous interaction, shocked material will be spilled out from the interaction region and form an outflow. We calculate the temporal evolution of synchrotron emission by assuming that the shock accelerates a fraction of electrons in the outflow into a nonthermal distribution. We find that radio emission from the debris–torus collision generates a steep-rise and slow-decline radio light curve due to the sharp edge and dense gas of the dusty torus, where the radio outburst is delayed relative to the main optical/X-ray outburst by several years or even several tens of years. We apply our model to a TDE that happened in a narrow-line Seyfert I galaxy (PS16dtm), where both the radio spectrum and the light curve can be roughly reproduced. Future high-sensitivity, wide-field-of-view radio surveys have the opportunity to detect more such radio flares.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Attribute-Aware RBFs: Interactive Visualization of Time Series Particle Volumes Using RT Core Range Queries

Smoothed-particle hydrodynamics (SPH) is a mesh-free method used to simulate volumetric media in fluids, astrophysics, and solid mechanics. Visualizing these simulations is problematic because these datasets often contain millions, if not billions of particles carrying physical attributes and moving over time. Radial basis functions (RBFs) are used to model particles, and overlapping particles are interpolated to reconstruct a high-quality volumetric field; however, this interpolation process is expensive and makes interactive visualization difficult. Existing RBF interpolation schemes do not account for color-mapped attributes and are instead constrained to visualizing just the density field. To address these challenges, we exploit ray tracing cores in modern GPU architectures to accelerate scalar field reconstruction. We use a novel RBF interpolation scheme to integrate per-particle colors and densities, and leverage GPU-parallel tree construction and refitting to quickly update the tree as the simulation animates over time or when the user manipulates particle radii. We also propose a Hilbert reordering scheme to cluster particles together at the leaves of the tree to reduce tree memory consumption. Finally, we reduce the noise of volumetric shadows by adopting a spatially temporal blue noise sampling scheme. Our method can provide a more detailed and interactive view of these large, volumetric, time-series particle datasets than traditional methods, leading to new insights into these physics simulations.

Particle Volumes↗

DOE OSTI2024

Acceleration mechanisms of energetic ion debris in laser-driven tin plasma EUV sources

Laser-driven tin plasmas are driving new-generation nanolithography as sources of extreme ultraviolet (EUV) radiation centered at 13.5 nm. A major challenge facing industrial EUV source development is predicting energetic ion debris produced during the plasma expansion that may damage the sensitive EUV channeling multilayer optics. Gaining a detailed understanding of the plasma dynamics and ion acceleration mechanisms in these sources could provide critical insights for designing debris mitigation strategies in future high-power EUV sources. Here, we develop a fully kinetic model of tin-EUV sources using one-dimensional particle-in-cell simulations to study ion debris acceleration, which will be valuable for cross-validation of radiation-hydrodynamic simulations. An inverse-bremsstrahlung heating operator is used to model the interaction of a tin target with an Nd:YAG laser, and thermal conduction is included through a Monte Carlo Coulomb collision operator. While the large-scale evolution is in reasonable agreement with analogous hydrodynamic simulations, the significant timescale for collisional equilibration between electrons and ions allows for the development of prominent two-temperature features. A collimated flow of energetic ions is produced with a spectrum that is significantly enhanced at high energies compared to fluid simulations. The dominant acceleration mechanism is found to be a large-scale electric field supported mainly by the electron pressure gradient, which is enhanced in the kinetic simulations due to the increased electron temperature. We discuss the implications of these results for future modeling of tin-EUV sources and the development of debris mitigation schemes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Bayesian calibration of viscous anisotropic hydrodynamic simulations of heavy-ion collisions

Owing to large pressure gradients at early times, standard hydrodynamic model simulations of relativistic heavy-ion collisions do not become reliable until O(1) fm/c after the collision. To address this one often introduces a prehydrodynamic stage that models the early evolution microscopically, typically as a conformal, weakly interacting gas. In such an approach the transition from the prehydrodynamic to the hydrodynamic stage is discontinuous, introducing considerable theoretical model ambiguity. Alternatively, fluids with large anisotropic pressure gradients can be handled macroscopically using the recently developed viscous anisotropic hydrodynamics (VAH). In high-energy heavy-ion collisions VAH is applicable already at very early times, and at later times transitions smoothly into conventional second-order viscous hydrodynamics. Here we present a Bayesian calibration of the VAH model with experimental data for Pb-Pb collisions at the LHC at $\sqrt{^{S}NN}$ = 2.76 A TeV. We find that the VAH model has the unique capability of constraining the specific viscosities of the quark-gluon plasma at higher temperatures than other previously used models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Multi frame radiography of supersonic water jets interacting with a foil target

Pulsed-power-driven underwater electrical explosion of cylindrical or conical wire arrays produces supersonic water jets that emerge from a bath, propagating through the air above it. Interaction of these jets with solid targets may represent a new platform for attaining materials at high pressure (>1010 Pa) conditions in a university-scale laboratory. However, measurements of the internal structure of such jets and how they interact with targets are difficult optically due to large densities and density contrasts involved. We utilized multi-frame x-ray radiographic imaging capabilities of the ID19 beamline at the European Synchrotron Radiation Facility to explore the water jet and its interaction with a 50 μm thick copper foil placed a few mm from the surface of water. The jet was generated with a ∼130 kA-amplitude current pulse of ∼450 ns rise time applied to a conical wire array. X-ray imaging revealed a droplet-type structure of the jet with an average density of <400 kg/m3 propagating with a velocity of ∼1400 m/s. Measurements of deformation and subsequent perforation of the target by the jet suggested pressures at the jet–target interface of ∼5 × 109 Pa. The results were compared to hydrodynamic simulations for better understanding of the jet parameters and their interaction with the foil target. These results can be used in future research to optimize the platform, and extend it to larger jet velocities in the case of higher driving currents supplied to the wire array.

Physics↗

DOE OSTI2022

Modeling and Analysis of a Novel Offshore Binary Species Free-Floating Longline Macroalgal Farming System

The investigation of innovative macroalgal cultivation is important and needed to optimize farming operations, increase biomass production, reduce the impact on the ecosystem, and lower system and operational costs. However, most macroalgal farming systems (MFSs) are stationary, which need to occupy a substantial coastal area, require extensive investment in farm infrastructure, and cost high fertilizer and anchoring expenses. This study aims to model, analyze, and support a novel binary species free-floating longline macroalgal cultivation concept. The expected outcomes could provide a basis for the design and application of the novel MFS to improve biomass production, decrease costs, and reduce the impact on the local ecosystem. In this paper, Saccharina latissima and Nereocystis luetkeana were modeled and validated, and coupled with longline to simulate the binary species MFS free float in various growth periods and associated locations along the US west coast. Further, the numerical predictions indicated the possibility of failure on the longline and breakage at the kelp holdfasts is low. However, the large forces due to an instantaneous change in dynamic loads caused by loss of hydrostatic buoyancy when the longline stretches out of the water would damage the kelps. Buoy-longline contact interactions could damage the buoy, resulting in the loss of the system by sinking. Furthermore, the kelp-longline and kelp-kelp entanglements could potentially cause kelp damage.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Collective dynamics of polarized spin-half fermions in relativistic heavy-ion collisions

Standard relativistic hydrodynamics has been successful in describing the properties of the strongly interacting matter produced in the heavy-ion collision experiments. Recently, there has been a significant theoretical advancement in this field to explain spin polarization of hadrons emitted in these processes. Although current models have successfully explained some of the experimental data based on the coupling between spin polarization and vorticity of the medium, they still lack a clear understanding of the differential measurements. This is commonly interpreted as an indication that the spin needs to be treated as an independent degree of freedom whose dynamics is not entirely bound to flow circulation. In particular, if the spin is a macroscopic property of the system, in equilibrium its dynamics should follow hydrodynamic laws. Here, we develop a framework of relativistic hydrodynamics which includes spin degrees of freedom from the quantum kinetic theory for Dirac fermions and use it for modeling the dynamics of matter. Following experimental observations, we assume that the polarization effects are small and derive conservation laws for the net baryon current, the energy–momentum tensor and the spin tensor based on the de Groot–van Leeuwen–van Weert definitions of these currents. We present various properties of the spin polarization tensor and its components, analyze the propagation properties of the spin polarization components, and derive the spin-wave velocity for arbitrary statistics. We find that only the transverse spin components propagate, analogously to the electromagnetic waves. Finally, using our framework, we study the space–time evolution of the spin polarization for the systems respecting certain space–time symmetries and calculate the mean spin polarization per particle, which can be compared to the experimental data. We find that, for some observables, our spin polarization results agree qualitatively with the experimental findings and other model calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Emergence of Microphysical Bulk Viscosity in Binary Neutron Star Postmerger Dynamics

Abstract In nuclear matter in isolated neutron stars, the flavor content (e.g., proton fraction) is subject to weak interactions, establishing flavor ( β -)equilibrium. However, there can be deviations from this equilibrium during the merger of two neutron stars. We study the resulting out-of-equilibrium dynamics during the collision by incorporating direct and modified Urca processes (in the neutrino-transparent regime) into general-relativistic hydrodynamics simulations with a simplified neutrino transport scheme. We demonstrate how weak-interaction-driven bulk viscosity in postmerger simulations can emerge and assess the bulk viscous dynamics of the resulting flow. We further place limits on the impact of the postmerger gravitational-wave strain. Our results show that weak-interaction-driven bulk viscosity can potentially lead to a phase shift of the postmerger gravitational-wave spectrum, although the effect is currently on the same level as the numerical errors of our simulation.

Astronomy & Astrophysics↗

DOE OSTI2022

Surface fluctuating hydrodynamics methods for the drift-diffusion dynamics of particles and microstructures within curved fluid interfaces

Here we introduce fluctuating hydrodynamics approaches on surfaces for capturing the drift-diffusion dynamics of particles and microstructures immersed within curved fluid interfaces of spherical shape. We take into account the interfacial hydrodynamic coupling, traction coupling with the surrounding bulk fluid, and thermal fluctuations. For fluid-structure interactions, we introduce Immersed Boundary Methods (IBM) and related Stochastic Eulerian-Lagrangian Methods (SELM) for curved surfaces. We use these approaches to investigate the statistics of surface fluctuating hydrodynamics and microstructures. For velocity autocorrelations, we find characteristic power-law scalings $τ^{-1}, τ^{-2}$, and plateaus can emerge. This depends on the physical regime associated with the geometry, surface viscosity, and bulk viscosity. This differs from the characteristic $τ^{-3/2}$ scaling for bulk three dimensional fluids. We develop theory explaining these observed power-laws associated with time-scales for dissipation within the fluid interface and coupling to the surrounding fluid. We then use our introduced methods to investigate a few example systems and roles of hydrodynamic coupling and thermal fluctuations including for the kinetics of passive particles and active microswimmers in curved fluid interfaces.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Diffusive hydrodynamics from integrability breaking

We describe here the crossover from generalized to conventional hydrodynamics in nearly integrable systems. Integrable systems have infinitely many conserved quantities, which spread ballistically, in general. When integrability is broken, only a few of these conserved quantities survive. The remaining conserved quantities are generically transported diffusively; we further derive a compact and general diffusion equation for these. The diffusion constant depends on the matrix elements of the integrability-breaking perturbation; for a certain class of integrability-breaking perturbations, including long-range interactions, the diffusion constant can be expressed entirely in terms of generalized hydrodynamic data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Modeling and simulations of hydrodynamic shocks in a plasma flowing across randomized ICF scale laser beams

High-energy laser beams interacting with flowing plasmas can produce a plasma response that leads to deflection of the beam, beam bending. Such beams have usually a speckle structure generated by optical smoothing techniques that reduce the spatial and temporal coherence in the laser field pattern. The cumulative plasma response from laser speckles slows down the velocity of the incoming flow by momentum conservation. For slightly super-sonic flow the cumulative plasma response to the ponderomotive force exerted by the beam speckle ensemble is the strongest, such that slowing down the flow to subsonic velocities leads eventually to the generation of a shock around the cross section of the beam. This scenario has been predicted theoretically and is confirmed here by our hydrodynamic simulations in two dimensions with speckled beams and in one dimension with a reduced model. The conditions of shock generation are given in terms of the ponderomotive pressure, speckle size and the flow velocity. The nonlinear properties of the shocks are analyzed using Rankine–Hugoniot relations. According to linear theory, temporally smoothed laser beams exhibit a higher threshold for shock generation. Numerical simulations with beams that are smoothed by spectral dispersion compare well with the linear theory results, diverging from those produced by beams with only a random phase plates in the nonlinear regime. The conditions necessary for shock generation and their effects on the laser plasma coupling in inertial confinement fusion (ICF) experiments are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

A conservative approach to scaling magneto-inertial fusion concepts to larger pulsed-power drivers

The Magnetized Liner Inertial Fusion (MagLIF) experimental platform [M. R. Gomez et al., Phys. Rev. Lett. 113, 155003 (2014)] represents the most successful demonstration of magneto-inertial fusion (MIF) techniques to date in pursuit of ignition and significant fusion yields. The pressing question remains regarding how to scale MIF concepts like MagLIF to more powerful pulsed-power drivers while avoiding significant changes in physical regimes that could adversely impact performance. In this work, we propose a conservative approach for scaling general MIF implosions, including MagLIF. Underpinning our scaling approach is a theoretical framework describing the evolution of the trajectory and thickness of a thin-walled, cylindrical, current-driven shell imploding on preheated, adiabatic fuel. By imposing that scaled implosions remain self-similar, we obtain a set of scaling rules expressing key target design parameters and performance metrics as functions of the maximum driver current I max . We identify several scaling paths offering unique, complementary benefits and trade-offs in terms of physics risks and driver requirements. Remarkably, when scaling present-day experiments to higher coupled energies, these paths are predicted to preserve or reduce the majority of known performance-degrading effects, including hydrodynamic instabilities, impurity mix, fuel energy losses, and laser-plasma interactions, with notable exceptions clearly delineated. In the absence of α heating, our scaling paths exhibit neutron yield per-unit-length scaling as $\tilde{Y}$∝[I$3\atop{max}$,I${4.14}\atop{max}$] and ignition parameter scaling as χ∝[I max ,I${2.14}\atop{max}$]. By considering the specific physics risks unique to each scaling path, we provide a roadmap for future investigations to evaluate different scaling options through detailed numerical studies and scaling-focused experiments on present-day facilities. Overall, these results highlight the potential of MIF as a key component of the national ignition effort.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Universal dynamics of coarsening during polymer-polymer thin-film spinodal dewetting kinetics

The dewetting dynamics of a supported bilayer polymer thin film on a solid substrate is investigated using grazing incidence x-ray photon correlation spectroscopy. We find that the top layer dewets via the spinodal mechanism. The kinetics of the dewetting is studied by monitoring the time evolution of the surface diffuse x-ray scattering intensity. We study the time evolution of fluctuations about the average surface structure by measuring the two-time x-ray intensity fluctuation correlation functions. Using these two-time correlation functions we quantify the crossover from early-time diffusive dynamics to hydrodynamics. The early diffusive regime satisfies dynamic universality. The two-time correlation functions also quantify the onset of hydrodynamic effects. The hydrodynamic regime is observed during the spinodal dewetting process as these interactions are not screened.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Influence of Lake Ice Biases in Reanalysis Data on Downscaled Climate Simulations over the Great Lakes Region

This data package contains observation-based and model-simulated datasets (all provided in NetCDF format) for evaluating how wintertime lake-ice representation affects regional weather and climate over the Laurentian Great Lakes (freshwater lake ecosystem) during the high–ice-cover winter of 2009. The observational component includes: (1) Stage IV gridded precipitation at 4 km, hourly resolution for January–February 2009 over the Great Lakes region (radar–gauge multisensor precipitation analyses); (2) Great Lakes Surface Environmental Analysis (GLSEA) satellite-derived lake-ice coverage at 1.3 km, daily resolution for the 2009 winter months, providing ice coverage over Lakes Superior, Michigan, Huron, Erie, and Ontario; and (3) in situ measurements at the Standard Rock site on Lake Superior from the Great Lakes Evaporation Network (GLEN) at hourly resolution, including near-surface atmospheric variables and sensible and latent heat fluxes (air–lake exchange) at a fixed point location. The modeling component provides corresponding fields from two simulations, both archived at 4 km, hourly resolution: a standalone Weather Research Forecasting model (WRF) run driven by the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis 5 (ERA5), and a two-way coupled model using WRF and the Finite Volume Community Ocean Model (WRF-FVCOM, a 3-D hydrodynamic lake model). These outputs include variables relevant to air–lake interaction and lake-effect processes (e.g., near-surface temperature, humidity, wind, precipitation, and surface turbulent fluxes), enabling direct comparison with the observational datasets. Users can analyze and visualize these NetCDF files with common tools such as Python (e.g., xarray, netCDF4, numpy, pandas), NCO/CDO, Panoply, or ncview; NetCDF variables can also be converted to other formats (e.g., CSV, GeoTIFF) using these utilities.

EARTH SCIENCE > ATMOSPHERE > ATMOSPHERIC TEMPERATU↗

DOE OSTI2026

Beam-water Interaction System for PFAS Treatment

Per- and polyfluoroalkyl substances (PFAS) are persistent environmental contaminants that resist conventional water treatment due to the strength of their carbon fluorine bonds. Electron-beam irradiation has shown potential for PFAS degradation; however, its effectiveness is limited by shallow energy penetration, restricting treatment to the upper region of the water surface. This project addresses this limitation through the design of a compact mixing chamber and recirculation system to enhance surface renewal and improve treatment efficiency. The system, constrained to approximately 10 cm 20 cm 20 cm, was developed using computer-aided design, computational fluid dynamics (CFD), and MATLAB-based hydraulic analysis to evaluate flow behavior, head loss, and pump compatibility. A perforated-baffle geometry and controlled flow range of 1 3 GPM were implemented to promote uniform mixing while maintaining a stable free surface for electron-beam exposure. Modeling results indicate improved flow distribution and reduced dead zones, supporting more consistent interaction between contaminated water and the beam. The final design enhances hydrodynamic performance within system constraints, contributing to scalable electron-beam water treatment technologies for improved PFAS degradation.

Ramos, Gerardo [Northern Illinois U.]↗