Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Hydrodynamic theory”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Optimizing Power Generation of a Bottomraised Oscillating Surge Wave Energy Converter Using a Theoretical Model

Preliminary sizing of an oscillating surge wave energy converter (OSWEC) is an iterative process that relies on knowledge of the relevant hydrodynamic coefficients for a given geometry. Often, the complex definition of the device geometry requires coefficients to be obtained through experiments or numerical boundary element solvers such as WAMIT. These techniques demand significant user and computational effort, therefore inhibiting the fine scale parametrization of object dimensions. In this study, a theoretical model, originally presented in Michele et al. (2016), is developed and demonstrated to efficiently optimize the power production for an OSWEC device (subjected to certain environmental conditions) with variations in device widths, heights, and distances from the seabed. Assuming negligible device thickness, the OSWEC motions are modelled as a bottom-raised two-dimensional flap in regular waves using potential flow theory formulated in elliptical coordinates. The solutions to this diffraction-radiation problem are obtained using Mathieu functions with appropriate boundary conditions. The resulting potentials are then used to derive frequency-dependent expressions for the added mass and radiation damping coefficients, along with wave excitation magnitude in surge, pitch, and coupled surge-pitch motions. Good agreement in hydrodynamic coefficient curves is shown between the theoretical model and the numerical results obtained from the boundary element-based program WAMIT. The theoretical model is then employed to maximize the time-averaged output power while maintaining or reducing the hinge reaction force, with variations in device dimensions, wave frequency and amplitude.

ACE↗
DOE OSTI2022

$T\overline{T}$ deformations and the width of fundamental particles

We provide a simple geometric meaning for deformations of so-called $T\overline{T}$ type in relativistic and non-relativistic systems. Deformations by the cross products of energy and momentum currents in integrable quantum field theories are known to modify the thermodynamic Bethe ansatz equations by a “CDD factor”. In turn, CDD factors may be interpreted as additional, fixed shifts incurred in scattering processes: a finite width added to the fundamental particles (or, if negative, to the free space between them). We suggest that this physical effect is a universal way of understanding $T\overline{T}$ deformations, both in classical and quantum systems. We first show this in non-relativistic systems, with particle conservation and translation invariance, using the deformation formed out of the densities and currents of particles and momentum. This holds at the level of the equations of motion, and for any interaction potential, integrable or not. We then argue, and show by similar techniques in free relativistic particle systems, that $T\overline{T}$ deformations of relativistic systems produce the equivalent phenomenon, accounting for length contractions. We also show that, in both the relativistic and non-relativistic cases, the width of particles is equivalent to a state-dependent change of metric, where the distance function discounts the particles’ widths, or counts the additional free space. This generalises and explains the known field-dependent coordinate change describing $T\overline{T}$ deformations. The results connect such deformations with generalised hydrodynamics, where the relations between scattering shifts, widths of particles and state-dependent changes of metric have been established.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Computational study of inertial migration of prolate particles in a straight rectangular channel

Inertial migration of spherical particles has been investigated extensively using experiments, theory, and computational modeling. Yet, a systematic investigation of the effect of particle shape on inertial migration is still lacking. Herein, we numerically mapped the migration dynamics of a prolate particle in a straight rectangular microchannel using smoothed particle hydrodynamics at moderate Reynolds number flows. After validation, we applied our model to 2:1 and 3:1 shape aspect ratio particles at multiple confinement ratios. Their effects on the final focusing position, rotational behavior, and transitional dynamics were studied. In addition to the commonly reported tumbling motion, for the first time, we identified a new logrolling behavior of a prolate ellipsoidal particle in the confined channel. This new behavior occurs when the confinement ratio is above an approximate threshold value of K = 0.72. Our microfluidic experiments using cell aggregates with similar shape aspect ratio and confinement ratio confirmed this new predicted logrolling motion. In this study, we also found that the same particle can undergo different rotational modes, including kayaking behavior, depending on its initial cross-sectional position and orientation. Furthermore, we examined the migration speed, angular velocity, and rotation period as well as their dependence on both particle shape aspect ratio and confinement ratio. Our findings are especially relevant to the applications where particle shape and alignment are used for sorting and analysis, such as the use of barcoded particles for biochemical assays through optical reading, or the shape-based enrichment of microalgae, bacteria, and chromosomes.

42 ENGINEERING↗
DOE OSTI2025

Supercooled Goldstone Bosons at the QCD Chiral Phase Transition

We discuss a universal nonequilibrium enhancement of long-wavelength Goldstone bosons induced by quenches to the broken phase in Model G—the dynamical universality class of an 𝑂⁡(4) antiferromagnet and the chiral phase transition in QCD. Scaling arguments for the coarsening dynamics describing the formation of the chiral condensate predict a parametric enhancement in the infrared spectra of Goldstone bosons, a prediction confirmed by stochastic simulations of the transition. The details of the enhancement are determined by the nonlinear dynamics of a superfluid effective theory, which is a limit of Model G reflecting the broken 𝑂⁡(4) symmetry. Our results translate to a parametric enhancement of low-momentum pions in heavy-ion collisions at the LHC, which are underpredicted in current hydrodynamic models without critical dynamics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Alpha-heating analysis of burning plasma and ignition experiments on the National Ignition Facility

A recent experiment conducted on the National Ignition Facility (NIF) described in the study by Abu-Shawareb et al. achieved a fusion yield output of 1.3 MJ from ~ 220 kJ of x-ray energy absorbed by the capsule, demonstrating remarkable progress in the field of laser driven inertial confinement fusion. In the study by A. R. Christopherson [“Effects of charged particle heating on the hydrodynamics of inertially confined plasmas,” Ph.D. thesis (2020)], the plasma conditions needed to claim the onset of ignition and burn propagation were outlined and multiple criterion were provided to assess progress in inertial fusion experiments. In this work, we modify the metrics from A. R. Christopherson to accurately calculate performance metrics for indirect-drive experiments on the NIF. We also show that performance metric trends observed in NIF data are consistent with theory and simulations. This analysis indicates that all the identified criterion for ignition and burn propagation have been exceeded by experiment 210 808.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Modeling of streamflow in a 30 km long reach spanning 5 years using OpenFOAM 5.x

Abstract. Developing accurate and efficient modeling techniques for streamflow at the tens-of-kilometers spatial scale and multi-year temporal scale is critical for evaluating and predicting the impact of climate- and human-induced discharge variations on river hydrodynamics. However, achieving such a goal is challenging because of limited surveys of streambed hydraulic roughness, uncertain boundary condition specifications, and high computational costs. We demonstrate that accurate and efficient three-dimensional (3-D) hydrodynamic modeling of natural rivers at 30 km and 5-year scales is feasible using the following three techniques within OpenFOAM, an open-source computational fluid dynamics platform: (1) generating a distributed hydraulic roughness field for the streambed by integrating water-stage observation data, a rough wall theory, and a local roughness optimization and adjustment strategy; (2) prescribing the boundary condition for the inflow and outflow by integrating precomputed results of a one-dimensional (1-D) hydraulic model with the 3-D model; and (3) reducing computational time using multiple parallel runs constrained by 1-D inflow and outflow boundary conditions. Streamflow modeling for a 30 km long reach in the Columbia River (CR) over 58 months can be achieved in less than 6 d using 1.1 million CPU hours. The mean error between the modeled and the observed water stages for our simulated CR reach ranges from −16 to 9 cm (equivalent to approximately ±7 % relative to the average water depth) at seven locations during most of the years between 2011 and 2019. We can reproduce the velocity distribution measured by the acoustic Doppler current profiler (ADCP). The correlation coefficients of the depth-averaged velocity between the model and ADCP measurements are in the range between 0.71 and 0.83 at 75 % of the survey cross sections. With the validated model, we further show that the relative importance of dynamic pressure versus hydrostatic pressure varies with discharge variations and topography heterogeneity. Given the model's high accuracy and computational efficiency, the model framework provides a generic approach to evaluate and predict the impacts of climate- and human-induced discharge variations on river hydrodynamics at tens-of-kilometers and decadal scales.

58 GEOSCIENCES↗
DOE OSTI2020

Horizon constraints on holographic Green’s functions

We explore a new class of general properties of thermal holographic Green’s functions that can be deduced from the near-horizon behaviour of classical perturbations in asymptotically anti-de Sitter spacetimes. We show that at negative imaginary Matsubara frequencies and appropriate complex values of the wavenumber the retarded Green’s functions of generic operators are not uniquely defined, due to the lack of a unique ingoing solution for the bulk perturbations. From a boundary perspective these ‘pole-skipping’ points correspond to locations in the complex frequency and momentum planes at which a line of poles of the retarded Green’s function intersects with a line of zeroes. As a consequence the dispersion relations of collective modes in the boundary theory at energy scales ω ~ T are directly constrained by the bulk dynamics near the black-brane horizon. For the case of conserved U(1) current and energy-momentum tensor operators we give examples where the dispersion relations of hydrodynamic modes pass through a succession of poleskipping points as real wavenumber is increased. We discuss implications of our results for transport, hydrodynamics and quantum chaos in holographic systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Staggered scheme for the compressible fluctuating hydrodynamics of multispecies fluid mixtures

Here, we present a numerical formulation for the solution of nonisothermal, compressible Navier-Stokes equations with thermal fluctuations to describe mesoscale transport phenomena in multispecies fluid mixtures. The novelty of our numerical method is the use of staggered grid momenta along with a finite volume discretization of the thermodynamic variables to solve the resulting stochastic partial differential equations. The key advantages of the numerical scheme are that it significantly simplifies the discretization of diffusive and stochastic momentum fluxes into a more compact form, and it provides an unambiguous prescription of boundary conditions involving pressure. The staggered grid scheme more accurately reproduces the equilibrium static structure factor of hydrodynamic fluctuations in gas mixtures compared to a collocated scheme described previously by Balakrishnan et al. [Phys. Rev. E 89, 013017 (2014)1539-375510.1103/PhysRevE.89.013017]. The numerical method is tested for ideal noble gases mixtures under various nonequilibrium conditions, such as applied thermal and concentration gradients, to assess the role of cross-diffusion effects, such as Soret and Dufour, on the long-ranged correlations of hydrodynamic fluctuations, which are also more accurately reproduced compared to the collocated scheme. We numerically study giant nonequilibrium fluctuations driven by concentration gradients and fluctuation-driven Rayleigh-Taylor instability in gas mixtures. Wherever applicable, excellent agreement is observed with theory and measurements from the direct simulation Monte Carlo method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Multiphase tin equation of state using density functional theory

In this work, we perform density functional theory (DFT) calculations of five solid phases and the liquid phase of tin. The calculations include cold curves of the five solid phases, phonon calculations in the quasiharmonic approximation over a range of volumes for each solid phase, and DFT-based molecular dynamics (DFT-MD) simulations of the liquid phase, including those of the melt curve using the Z method. Using the DFT results, we construct a tabular multiphase sesame equation of state for tin, referred to as sesame 2162. Comparisons to experimental data are made and show a high level of agreement in isobaric data, isothermal data, shock data, and phase boundary measurements, including measurements of the melt curve. The 2162 EOS will be useful for hydrodynamics simulations and has been designed with an eye toward hydrodynamics simulations that incorporate materials strength models and allow for modeling of the kinetics of phase transitions.

3-dimensional systems↗
DOE OSTI2019

Coupling constant corrections in a holographic model of heavy ion collisions with nonzero baryon number density

Sufficiently energetic collisions of heavy ions result in the formation of a droplet of a strongly coupled liquid state of QCD matter known as quark-gluon plasma. By using gauge-gravity duality (holography), a model of a rapidly hydrodynamizing and thermal- izing process like this can be constructed by colliding sheets of energy density moving at the speed of light and tracking the subsequent evolution. In this work, we consider the dual gravitational description of such collisions in the most general bulk theory with a four-derivative gravitational action containing a dynamical metric and a gauge field in five dimensions. Introducing the bulk gauge field enables the analysis of collisions of sheets which carry nonzero “baryon” number density in addition to energy density. Introducing the four-derivative terms enables consideration of such collisions in a gauge theory with finite gauge coupling, working perturbatively in the inverse coupling. While the dynamics of energy and momentum in the presence of perturbative inverse-coupling corrections has been analyzed previously, here we are able to determine the effect of such finite coupling corrections on the dynamics of the density of a conserved global charge, which we take as a model for the dynamics of nonzero baryon number density. In accordance with expectations, as the coupling is reduced we observe that after the collisions less baryon density ends up stopped at mid-rapidity and more of it ends up moving near the lightcone.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Numerical and Analytical Modeling of Flow Partitioning in Partially Saturated Fracture Networks

Infiltration processes in fractured-porous media remain a crucial, yet not very well understood component of recharge and vulnerability assessment. Under partially-saturated conditions flows in fractures, percolating fracture networks, and fault zones contribute to the fastest spectrum of infiltration velocities via preferential pathways. Specifically, the partitioning dynamics at fracture intersections determine the magnitude of flow fragmentation into vertical and horizontal components and hence the bulk flow velocity and dispersion of fracture networks. In this work, we derive an approximate analytical solution for the partitioning process and validate it using smoothed particle hydrodynamics simulations. The developed transfer function is conceptually based on simulation results and laboratory experiments carried out in previous works. It allows to efficiently simulate flow through fracture networks with simple cubic structure and arbitrary number of fractures and aperture sizes via linear response theory and convolution of a given input signal. We derive a non-dimensional bulk flow velocity ($\widetilde{v}$) and dispersion coefficient ($\widetilde{D}$) to characterize fracture networks in terms of dimensionless horizontal and vertical time scales $\tau_m$ and $\tau_0$. The dispersion coefficient is shown to strongly depend on the horizontal time scale and converges towards a constant value of $0.08$ within reasonable fluid and geometrical parameter ranges, while the non-dimensional velocity exhibits a characteristic $\widetilde{v} \sim \tau_m^{-1/2}$ scaling. Given that hydraulic information is often only available at limited places within (fractured-porous) aquifer system, such as boreholes or springs, our study intends to provide a rudimentary analytical concept to potentially reconstruct internal fracture network geometries from external boundary information, e.g., the dispersive properties of discharge (groundwater level fluctuations).

42 ENGINEERING↗
DOE OSTI2025

A fluctuating hydrodynamics model for nanoscale surfactant-laden interfaces

A multispecies diffuse interface model is formulated in a fluctuating hydrodynamics framework for the purpose of simulating surfactant interfaces at the nanoscale. The model generalizes previous work to ternary mixtures, employing a Cahn-Hilliard free energy density combined with incompressible, isothermal fluctuating hydrodynamics where dissipative fluxes include both deterministic and stochastic terms. The intermolecular parameters in the free energy are chosen such that one species acts as a partially miscible surfactant. From Laplace pressure measurements, we show that in this model the surface tension decreases linearly with surfactant concentration, leading to Marangoni convection for interfaces with concentration gradients. In the capillary wave spectrum for interfaces with and without surfactant, we find that for the former, the spectrum deviates significantly from classical capillary wave theory, presumably due to Gibbs elasticity. In non-equilibrium simulations of the Rayleigh-Plateau instability, deterministic simulations showed that the surfactant delays pinching of a fluid cylinder into droplets. However, stochastic simulations indicate that thermal fluctuations disrupt the surfactant's stabilizing effect. Similarly, the spreading of a patch of surfactant, driven by Marangoni convection, was found to be partially suppressed by thermal fluctuations.

Capillary waves↗
DOE OSTI2020

Fully compressible simulations of waves and core convection in main-sequence stars

Context: Recent, nonlinear simulations of wave generation and propagation in full-star models have been carried out in the anelastic approximation using spectral methods. Although it makes long time steps possible, this approach excludes the physics of sound waves completely and requires rather high artificial viscosity and thermal diffusivity for numerical stability. A direct comparison with observations is thus limited. Aims: In this study, we explore the capabilities of our compressible multidimensional Seven-League Hydro (SLH) code to simulate stellar oscillations. Methods: We compare some fundamental properties of internal gravity and pressure waves in 2D SLH simulations to linear wave theory using two test cases: (1) an interval gravity wave packet in the Boussinesq limit and (2) a realistic 3 M ⊙ stellar model with a convective core and a radiative envelope. Oscillation properties of the stellar model are also discussed in the context of observations. Results: Our tests show that specialized low-Mach techniques are necessary when simulating oscillations in stellar interiors. Basic properties of internal gravity and pressure waves in our simulations are in good agreement with linear wave theory. As compared to anelastic simulations of the same stellar model, we can follow internal gravity waves of much lower frequencies. The temporal frequency spectra of velocity and temperature are flat and compatible with the observed spectra of massive stars. Conclusion: The low-Mach compressible approach to hydrodynamical simulations of stellar oscillations is promising. Our simulations are less dissipative and require less luminosity boosting than comparable spectral simulations. The fully-compressible approach allows for the coupling of gravity and pressure waves in the outer convective envelopes of evolved stars to be studied in the future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Physics of nanoscale immiscible fluid displacement

In this work, we investigate immiscible fluid displacement at small scales where slip lengths are on the order of characteristic system sizes, whereby Cox's law is not expected to be valid. Molecular dynamics simulations show that in this limit hydrodynamic bending becomes small and interfaces remain approximately spherical. In this case the only relevant angle for describing the interface shape is the dynamic microscopic angle at the fluid-solid interface. In our simulations, this angle is found to be described well by the molecular-kinetic theory originally proposed by Blake and Haynes. In general, this implies a different functional dependence between the contact angle (and related quantities) and the flow speed (or capillary number); this is demonstrated for the case of the force on the boundary for immiscible fluid displacement in a two-dimensional channel.

42 ENGINEERING↗
DOE OSTI2025

ForestFlow: predicting the Lyman-α forest clustering from linear to nonlinear scales

On large scales, the Lyman-α forest provides insights into the expansion history of the Universe, while on small scales, it imposes strict constraints on the growth history, the nature of dark matter, and the sum of neutrino masses. This work introduces ForestFlow, a novel framework that bridges the gap between large- and small-scale analyses, which have traditionally relied on distinct modeling approaches. Using conditional normalizing flows, ForestFlow predicts the two Lyman-α linear biases (b δ and b η ) and six parameters describing small-scale deviations of the three-dimensional flux power spectrum (P 3D ) from linear theory as a function of cosmology and intergalactic medium physics. These are then combined with a Boltzmann solver to make consistent predictions, from arbitrarily large scales down to the nonlinear regime, for P 3D and any other statistics derived from it. Trained on a suite of 30 fixed-and-paired cosmological hydrodynamical simulations spanning redshifts from z = 2 to 4.5, ForestFlow achieves 3 and 1.5% precision in describing P 3D and the one-dimensional flux power spectrum (P 1D ) from linear scales to k = 5 Mpc −1 and k ∥ = 4 Mpc −1 , respectively. Thanks to its conditional parameterization, ForestFlow shows similar performance for ionization histories and two ΛCDM model extensions – massive neutrinos and curvature – even though none of these are included in the training set. This framework will enable full-scale cosmological analyses of Lyman-α forest measurements from the DESI survey.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Validation of SPH code Spheral to model interacting solid bodies in a supersonic flow

Contemporary discussions of planetary defense involve analyzing the risks posed by smaller sized, 20 to 200 m diameter, asteroids which are capable of breaking up in the atmosphere and generating a blast wave. Consequence assessments for this size class of asteroids are performed through fast-running analytic or semi-analytic models which are informed by high-fidelity hydrocode simulations of asteroid entry and breakup. However, insufficient historical data necessitates validating the independent physical processes which dominate airburst events. Here, the Fluid Solid Interface Smoothed Particle Hydrodynamics solver was previously used by Pearl et al. in 2023 to model the Chelyabinsk airburst and is used here to perform a series of validation simulations. The first effort involves modeling a cylinder in a hypersonic flow and comparing the bow shock geometry to that predicted by analytic theory. The second effort involves modeling the separation of two spherical bodies in supersonic flow and validating against experimental footage. Combined, these exercises demonstrate the ability of the code to model the flight-path of interacting solid bodies in a hypersonic flow.

Airburst↗
DOE OSTI2022

Far-from-equilibrium slow modes and momentum anisotropy in expanding plasma

The momentum distribution of particle production in heavy-ion collisions encodes information about thermalization processes in the early-stage quark-gluon plasma. We use kinetic theory to study the far-from-equilibrium evolution of an expanding plasma with an anisotropic momentumspace distribution. We identify slow and fast degrees of freedom in the far-from-equilibrium plasma from the evolution of moments of this distribution. At late times, the slow modes correspond to hydrodynamic degrees of freedom and are naturally gapped from the fast modes by the inverse of the relaxation time, $τ$$^{–1}_{R}$. At early times, however, there are an infinite number of slow modes with a gap inversely proportional to time, τ –1 . From the evolution of the slow modes we generalize the paradigm of the far-from-equilibrium attractor to vector and tensor components of the energy momentum tensor, and even to higher moments of the distribution function that are not part of the hydrodynamic evolution. We predict that initial-state momentum anisotropy decays slowly in the far-from-equilibrium phase and may persist until the relaxation time.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Alignments between Galaxies and the Cosmic Web at z ~ 1–2 in the IllustrisTNG Simulations

Galaxy formation theories predict that galaxy shapes and angular momenta have nonrandom alignments with the cosmic web. This leads to so-called intrinsic alignment between pairs of galaxies, which is important to quantify as a nuisance parameter for weak lensing. We study galaxy–cosmic web alignment in the IllustrisTNG suite of hydrodynamical simulations at redshifts 1 and 2, finding that alignment trends are consistent with previous studies. However, we find that the magnitude of the spin alignment signal is ~2.4× weaker than seen in previous studies of the Horizon-AGN simulation, suggesting that this signal may have a significant dependence on subgrid physics. Based on IllustrisTNG, we then construct mock observational spectroscopic surveys that can probe shape–cosmic web alignment at z ~ 1–2, modeled on the low-z galaxy redshift and IGM tomography surveys on the upcoming Subaru Prime Focus Spectrograph Galaxy Evolution (PFS GE) survey. However, even over box sizes of L = 205 h -1 Mpc, we find that global anisotropies induce a sample variance in the 2D projected alignment signal that depend on the projected direction; this induces significant errors in the observed alignment. We predict a 5.3σ detection of IllustrisTNG's shape alignment signal at z ~ 1 from Subaru PFS GE, although a detection would be challenging at z ~ 2. However, a rough rescaling of the relative alignment signal strengths between the TNG and Horizon-AGN simulations suggests that PFS GE should be able to more easily constrain the latter's stronger signal.

79 ASTRONOMY AND ASTROPHYSICS↗