Search NASA⌕ Search

SEARCH · Search NASA

Results for “hydrodynamic limit”

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 37 records · Page 2

A future of inertial confinement fusion without laser-plasma instabilities

From the beginning of inertial confinement fusion (ICF) research, laser-plasma instabilities excited by narrowband lasers have limited the hydrodynamic design space of all laser-based approaches to inertial fusion energy (IFE). With advances in broadband laser technologies, the next generation of ICF drivers will likely have large bandwidth that is engineered to mitigate laser-plasma instabilities, thereby expanding the hydrodynamic design space to include both robust high yields (>200−MJ) with large-energy laser systems (>4 MJ) and high gains (>10) with moderate-energy laser facilities (<2 MJ). State-of-the-art simulations indicate that laser bandwidths of a few percent are required to mitigate instabilities for IFE-relevant conditions. To test these models and demonstrate that high-bandwidth lasers mitigate laser-plasma instabilities, the Fourth-generation Laser for Ultrabroadband eXperiments will be used with the OMEGA Laser System. The goal is to provide the community with the confidence to invest in a multiple-beam high-bandwidth laser facility that will demonstrate the necessary ablation pressures for robust direct-drive ignition without detrimental levels of hot electrons that degrade the implosion performance. It is important to recognize that the highest performing ICF implosions will likely never be completely LPI free because of the significant advantages to maximizing the laser intensity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Emergent viscous hydrodynamics from a single quantum particle

We investigate an explicit example of how spatial decoherence can lead to hydrodynamic behavior in the late-time, long-wavelength regime of open quantum systems. We focus on the case of a single nonrelativistic quantum particle linearly coupled to a thermal bath of noninteracting harmonic oscillators at temperature T , a la Caldeira and Leggett. Taking advantage of decoherence in the position representation, we expand the reduced density matrix in powers of the off-diagonal spatial components, so that high-order terms are suppressed at late times. Truncating the resulting power series at second order leads to a set of dissipative transient hydrodynamic equations similar to the nonrelativistic limit of equations widely used in simulations of the quark-gluon plasma formed in ultrarelativistic heavy-ion collisions. Transport coefficients are directly determined by the damping constant γ , which quantifies the influence of the environment. The asymptotic limit of our hydrodynamic equations reduces to the celebrated Navier-Stokes equations for a compressible fluid in the presence of a drag force. Furthermore, our results shed new light on the onset of hydrodynamic behavior in open quantum systems where a system with few degrees of freedom is coupled to a large thermal environment.

Hydrodynamics↗

Physics Guided Simulation of Electrostatic Discharge: Technical Report

Triboelectrically-charged objects may create threshold sparks, electrostatic discharge (ESD) events, to equilibrate charge between themselves and other relatively charged objects. ESD events exhibit many complex physical phenomena. They are a nexus of several fields of physics with disparate characteristic scales: plasma physics, chemical kinetics, hydrodynamics, circuit models, etc. These scales can span many orders of magnitude from the varied collisions thermalizing information within a plasma on the $\mathcal{O}(fs/ps)$ to the physical size of the plasma channel on the $\mathcal{O}(100µm)$, to the speed of a nonlinear hydrodynamic wave propagating at $\mathcal{O}(µm, ns)$. These threshold ESD events may occur in situations of programmatic importance, delivering energy and power profiles to a “victim load” generating deleterious consequences. To predict and mitigate these consequences we must answer questions about the spark’s energy budget: how much energy goes into producing the spark channel; how much gets radiated away; how much energy is advected away into the hydrodynamics; and how much energy is delivered to a victim load. An ESD simulation toolset has been created and evolved in order to answer these questions. An appropriate, physics guided implementation for simulation can be done by gaining insight into its constituent physics and leveraging that intuition to choose a suitable numerical operator. We examine in detail the chemical kinetics, circuit discharge, and hydrodynamics to deter mine dominant regimes, values, timescales, and interactions to uncover the underlying physical dynamics. We also examine and propose model reduction schemes for high-dimensional chemical kinetics. We use past and current work with experimentally validated and theoretically-verified hydrodynamics to calculate applicability limits of the non-ionizing strong shock limit. We quantify the energy budget from a hydrodynamic perspective and demonstrate that a significant fraction of the stored energy is “earmarked” for hydrodynamic advection as an energy terminus. Lastly, we combine the constituent physics of an ESD event (chemical kinetics, circuit model, and hydrodynamics) into a cohesive, actionable toolset and obtain promising results from an isothermal test case. We then propose a viable, modular evolution of the ESD toolset based upon the performed examination of the physics uncovering dominant physical scales and the stiffness of the compositional differential system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spin relaxation rate for heavy quarks in weakly coupled QCD plasma

We compute the relaxation rate of the spin density of heavy quarks in a perturbative QCD plasma to leading-log order in the coupling constant g. The spin relaxation rate Γ s in spin hydrodynamics is shown to be Γ s ~ g 4 log(1/g)T(T/M) 2 in the heavy-quark limit T/M $\ll$ 1, which is smaller than the relaxation rate of other non- hydrodynamic modes by additional powers of T/M. We demonstrate three different methods to evaluate the spin relaxation rate: 1) the Green-Kubo formula in the spin hydrodynamic regime, 2) the spin density correlation function in the strict hydrodynamic limit, and 3) quantum kinetic theory of the spin distribution function in momentum space. We highlight the interesting differences between these methods, while they are ultimately connected to each other by the underlying Ward-Takahashi identity for the non-conserved spin density.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Effects of radiation pressure on the evaporative wind of HD 209458b

ABSTRACT The role of radiation pressure in shaping exoplanet photoevaporation remains a topic of contention. Radiation pressure from the exoplanet’s host star has been proposed as a mechanism to drive the escaping atmosphere into a ‘cometary’ tail and explain the high velocities observed in systems where mass-loss is occurring. In this paper, we present results from high-resolution 3D hydrodynamic simulations of a planet similar to HD 209458b. We self-consistently launch a wind flowing outwards from the planet by calculating the ionization and heating resulting from incident high-energy radiation, and account for radiation pressure. We first present a simplified calculation, setting a limit on the Lyman-α flux required to drive the photoevaporated planetary material to larger radii and line-of-sight velocities. We then present the results of our simulations, which confirm the limits determined by our analytic calculation. We thus demonstrate that, within the limits of our hydrodynamic simulation and for the Lyman-α fluxes expected for HD 209458, radiation pressure is unlikely to significantly affect photoevaporative winds or to explain the high velocities at which wind material is observed, though further possibilities remain to be investigated.

79 ASTRONOMY AND ASTROPHYSICS↗

Effects of charge exchange on the evaporative wind of HD 209458b

ABSTRACT The role of charge exchange in shaping exoplanet photoevaporation remains a topic of contention. Exchange of electrons between stellar wind protons from the exoplanet’s host star and neutral hydrogen from the planet’s wind has been proposed as a mechanism to create ‘energetic neutral atoms’ (ENAs), which could explain the high absorption line velocities observed in systems where mass-loss is occurring. In this paper, we present results from three-dimensional hydrodynamic simulations of the mass-loss of a planet similar to HD 209458b. We self-consistently launch a planetary wind by calculating the ionization and heating resulting from incident high-energy radiation, inject a stellar wind into the simulation, and allow electron exchange between the stellar and planetary winds. We predict the potential production of ENAs by the wind–wind interaction analytically, and then present the results of our simulations, which confirm the analytic limits. Within the limits of our hydrodynamic simulation, we find that charge exchange with the stellar wind properties examined here is unable to explain the absorption observed at high Doppler velocities.

79 ASTRONOMY AND ASTROPHYSICS↗

Dynamics in a stellar convective layer and at its boundary: Comparison of five 3D hydrodynamics codes

Our ability to predict the structure and evolution of stars is in part limited by complex, 3D hydrodynamic processes such as convective boundary mixing. Hydrodynamic simulations help us understand the dynamics of stellar convection and convective boundaries. However, the codes used to compute such simulations are usually tested on extremely simple problems and the reliability and reproducibility of their predictions for turbulent flows is unclear. We define a test problem involving turbulent convection in a plane-parallel box, which leads to mass entrainment from, and internal-wave generation in, a stably stratified layer. We compare the outputs from the codes FLASH, MUSIC, PPMSTAR, PROMPI, and SLH, which have been widely employed to study hydrodynamic problems in stellar interiors. The convection is dominated by the largest scales that fit into the simulation box. All time-averaged profiles of velocity components, fluctuation amplitudes, and fluxes of enthalpy and kinetic energy are within ≲3σ of the mean of all simulations on a given grid (128 3 and 256 3 grid cells), where σ describes the statistical variation due to the flow’s time dependence. They also agree well with a 512 3 reference run. The 128 3 and 256 3 simulations agree within 9% and 4%, respectively, on the total mass entrained into the convective layer. The entrainment rate appears to be set by the amount of energy that can be converted to work in our setup and details of the small-scale flows in the boundary layer seem to be largely irrelevant. Our results lend credence to hydrodynamic simulations of flows in stellar interiors. We provide in electronic form all outputs of our simulations as well as all information needed to reproduce or extend our study.

79 ASTRONOMY AND ASTROPHYSICS↗

Insights into the mechanism of electrochemical chloride oxidation in ethanol from X-ray photoelectron spectroscopy, quiescent solution voltammetry, and rotating ring-disk electrodes

The wide availability of bio-derived alcohols provides the impetus to develop processes that convert them to valuable chemicals. The chloride ion is a redox mediator for electrocatalytic ethanol oxidation to 1,1-diethoxyethane (1,1-DEE) through an ethyl hypochlorite (EtOCl) intermediate, and this paper describes the chloride oxidation reaction (COR) to EtOCl on a glassy carbon (GC) electrode. Voltammetry measurements on a GC electrode in inert acetonitrile solvent combined with ex situ X-ray photoelectron spectroscopy (XPS) establish a Volmer step, where chloride ion from solution chemisorbs and is oxidized. In reactive ethanol solvent, ethanol adsorbs, and analyzing the current response in an LSV experiment supports a two-electron-transfer to form EtOCl, with chemisorption of the regenerated chloride. Koutecký–Levich (K–L) analysis on a rotating ring disk electrode (RRDE) shows that the kinetic rate constant of the COR in ethanol is on the order of 10 −8 cm s −1 , which is five orders of magnitude faster than the direct alcohol oxidation reaction in a kinetically limited regime. This hydrodynamic approach in understanding the electrochemistry of this non-aqueous system extends the possibilities for mediated electrocatalysis in neat alcohol solvents.

Van Daele, Ryan D. [University of Michigan, Ann Ar↗

Mesoscopic finite-size effects of unconventional electron transport in PdCoO 2

A wide range of unconventional transport phenomena has recently been observed in single-crystal delafossite metals. Here, we present a theoretical framework to elucidate electron transport using a combination of first-principles calculations and numerical modeling of the anisotropic Boltzmann transport equation. Using PdCoO 2 as a model system, we study different microscopic electron and phonon scattering mechanisms and establish the mean free path hierarchy of quasiparticles at different temperatures. We treat the anisotropic Fermi surface explicitly to numerically obtain experimentally-accessible transport observables, which bridge between the “diffusive,” “ballistic,” and “hydrodynamic” transport regime limits. Further, we illustrate that the distinction between the “quasiballistic” and “quasihydrodynamic” regimes is challenging and often needs to be quantitative in nature. From first-principles calculations, we populate the resulting transport regime plots and demonstrate how the Fermi surface orientation adds complexity to the observed transport signatures in micrometer-scale devices. Our work provides key insights into microscopic interaction mechanisms on open hexagonal Fermi surfaces and establishes their connection to the macroscopic electron transport in finite-size channels.

36 MATERIALS SCIENCE↗

Performance enhancement of direct-drive shock-augmented ignition inertial fusion implosions through shock timing optimization

Shock-augmented ignition (SAI) [R. H. H. Scott et al., Phys. Rev. Lett. 129, 195001 (2022)] is an alternative inertial confinement fusion concept that is designed to achieve high energy gain by combining improved resilience to instabilities with enhanced fuel compression. In SAI, lower implosion velocities can improve hydrodynamic stability and limited laser intensities reduce the excitation of detrimental laser-plasma instabilities, enabling greater areal density accumulation. Here we report an experimental investigation of SAI using a series of warm D 2 implosions. By timing the augmenting shock to arrive at the implosion center immediately prior to peak compression, the areal density and hot-spot pressure increase by 50% and 112%, respectively, relative to unoptimized designs. Ignition-scale simulations reproduce the observed timing dependence and demonstrate the potential of SAI for inertial fusion energy.

Direct drive↗

Validation of Local Structural Loads Computed by OpenFAST Against Measurements From the FOCAL Experimental Campaign

This work presents the validation of the local structural load modeling capability in OpenFAST for floating substructures based on data from the FOCAL experimental campaign. Previously, OpenFAST could only represent the floating substructure as a rigid body, and though this approach can model the global response of the floater in most cases, it is not able to capture the structural loads within the floater's individual members. Consideration of local substructure loads is important for some floating designs, because the pursuit of cost reduction often results in lighter and more flexible structures. To address this limitation, the HydroDyn (hydrodynamics) and SubDyn (substructure dynamics) modules of OpenFAST have been recently extended to account for the flexibility of floating substructures. To validate this new capability, we compare the results obtained by OpenFAST with data measured during the FOCAL experimental campaign, which analyzed a 1:70 scale performance-matched model of the IEA 15-MW reference turbine atop a modified University of Maine VolturnUS-S semisubmersible in a wave basin under the action of both wind and waves. For the purposes of the present work, the most important feature of the experiment is the presence of load cells at the root of each pontoon, and our objective is to assess how well those loads are reproduced by OpenFAST. To model the distributed hydrodynamic and hydrostatic loads along the floating substructure, we adopt a strip-theory approach based on the Morison equation, and we discuss the impact of different hydrodynamic modeling options (wave stretching, MacCamy-Fuchs correction, and second-order wave kinematics) on both motions and loads. For simplicity, we focus on wave-only conditions, both regular and irregular. The results demonstrate good overall agreement for the loads at the root of the pontoons for the waves analyzed in this work, especially given the assumptions and simplifications inherent to a simple strip-theory model.

floating offshore wind turbine↗

Validation of Local Structural Loads Computed by OpenFAST Against Measurements From the Focal Experimental Campaign: Preprint

This work presents the validation of the local structural load modeling capability in OpenFAST for floating substructures based on data from the FOCAL experimental campaign. Previously, OpenFAST could only represent the floating substructure as a rigid body, and though this approach can model the global response of the floater in most cases, it is not able to capture the structural loads within its individual members. Consideration of local substructure loads is important for some floating designs, as the pursuit of cost reduction often results in lighter and more flexible structures. To address this limitation, the HydroDyn (hydrodynamics) and SubDyn (substructure dynamics) modules of OpenFAST have been recently extended to account for the flexibility of floating substructures. To validate this new capability, we compare the results obtained by OpenFAST with data measured during the FOCAL experimental campaign, which analyzed a 1:70 scale performance-matched model of the IEA 15-MW reference turbine atop a modified University of Maine VolturnUS-S semisubmersible in a wave basin under the action of both wind and waves. For the purposes of the present work, the most important feature of the experiment is the presence of load cells at the root of each pontoon, and our objective is to assess how well those loads are reproduced by OpenFAST. To model the distributed hydrodynamic and hydrostatic loads along the floating substructure, we adopt a strip-theory approach based on the Morison equation, and we discuss the impact of different hydrodynamic modeling options (wave stretching, MacCamy-Fuchs correction, and second-order wave kinematics) on both motions and loads. For simplicity, we focus on wave-only conditions, both regular and irregular. The results demonstrate good overall agreement for the loads at the root of the pontoons for the waves analyzed in this work, especially given the assumptions and simplifications inherent to a simple strip-theory model.

floating offshore wind turbine↗

A multishock model for the density variance of anisotropic, highly magnetized, supersonic turbulence

ABSTRACT Shocks form the basis of our understanding for the density and velocity statistics of supersonic turbulent flows, such as those found in the cool interstellar medium (ISM). The variance of the density field, $\sigma ^2_{\rho /\rho _0}$, is of particular interest for molecular clouds (MCs), the birthplaces of stars in the Universe. The density variance may be used to infer underlying physical processes in an MC, and parametrizes the star formation (SF) rate of a cloud. However, models for $\sigma ^2_{\rho /\rho _0}$ all share a common feature – the variance is assumed to be isotropic. This assumption does not hold when a trans-/sub-Alfvénic mean magnetic field, ${B}_0$, is present in the cloud, which observations suggest is relevant for some MCs. We develop an anisotropic model for $\sigma _{\rho /\rho _0}^2$, using contributions from hydrodynamical and fast magnetosonic shocks that propagate orthogonal to each other. Our model predicts an upper bound for $\sigma _{\rho /\rho _0}^2$ in the high Mach number $(\mathcal {M})$ limit as small-scale density fluctuations become suppressed by the strong ${B}_0$. The model reduces to the isotropic $\sigma _{\rho /\rho _0}^2\!-\!\mathcal {M}$ relation in the hydrodynamical limit. To validate our model, we calculate $\sigma _{\rho /\rho _0}^2$ from 12 high-resolution, three-dimensional, supersonic, sub-Alfvénic magnetohydrodynamical (MHD) turbulence simulations and find good agreement with our theory. We discuss how the two MHD shocks may be the bimodally oriented overdensities observed in some MCs and the implications for SF theory in the presence of a sub-Alfvénic ${B}_0$. By creating an anisotropic, supersonic density fluctuation model, this study paves the way for SF theory in the highly anisotropic regime of interstellar turbulence.

Beattie, James R.↗

Disentangling the Effects of Vapor Pressure Deficit and Soil Water Availability on Canopy Conductance in a Seasonal Tropical Forest During the 2015 El Niño Drought

In land models used in Earth system models, vegetation water stress is heuristically related to soil water matric potential, and root water uptake in soil layers is simply proportional to root density to meet the transpiration demand. The efficiency of roots utilizing deeper soil water under water stress cannot be well represented, but it is of critical importance for transpiration regulation and carbon flux. To improve modeling of tropical forest response to drought, a plant hydrodynamic scheme (HYDRO), adapted from a trait-based model in the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), has been implemented in the Energy Exascale Earth System Model (E3SM) Land Model (ELM) to represent the hydraulic dynamics and root water uptake strategy. The adapted model conforms to the big-leaf photosynthesis model in ELM. Vegetation water stress in the model is a function of prognostic leaf water potential and is used to regulate stomatal conductance. Point simulations were performed at two tropical forest sites in the Amazon and Panama. Model sensitivity to the numerical solver, plant hydraulic traits, soil properties, and atmospheric forcing has been evaluated. Our results showed that, even though water stress factors calculated at our study sites using HYDRO are comparable to that simulated by the heuristic water stress model, soil moisture prediction is improved with soil water redistributions through roots. Simulated results using HYDRO is not as sensitive to plant hydraulic traits measured for different tree species as to the soil properties and atmospheric forcing. Simplified numerical solver may cause inaccuracy as to where roots uptake water. This study suggests that plant hydrodynamic representation is critical for improved simulations of soil moisture under water stresses.

54 ENVIRONMENTAL SCIENCES↗

Thermal transfer enhancement by hydrodynamic plasmons in electron bilayers

Here, we develop a theory of heat transfer induced by thermal charge fluctuations in two-dimensional electron double layers. We consider pristine systems comprised of identical layers and focus on the regime of sufficiently high temperatures and interlayer distances 𝑑, where the relevant charge fluctuations may be described using the hydrodynamic approach. In this limit heat transfer is dominated by the plasmon resonances. For systems with Galilean-invariant electron dispersion the interlayer thermal conductance 𝜘 is proportional to the kinematic viscosity of the electron liquid and decreases as 1/𝑑 4 . In the absence of Galilean invariance 𝜘∝𝜎/𝑑 3 , where 𝜎 is the intrinsic conductivity of the liquid. This strong enhancement can be traced to a drastically different broadening of plasmon resonances in systems with and without Galilean invariance.

2-dimensional systems↗

Testing and Expertise for Marine Energy (TEAMER) Program Support - Numerical Modeling Assistance for iProTech's "PIP" WEC Device: Cooperative Research and Development Final Report, CRADA Number CRD-20-17303

Develop a time-domain, 3 degree of freedom numerical model of IProTech's PIP device, including power take-off (PTO), that can be used to develop control algorithms and establish a baseline technology performance level (TPL). This is important to IProTech because having a good numerical model allows IProTech to optimize design and control algorithms on an ongoing basis. A baseline TPL enables the PIP technology to be benchmarked against other wave energy converter (WEC) concepts, providing (i) justification (or not) for investment in wave tank model testing, (ii) identifying optimum design parameters and control algorithms for the test model, (iii) specifying test conditions and (iv) interpreting wave tank test results. An accredited WEC-Sim model of the PIP device is an essential tool for further development of the technology. 1) The development of numerical models of the PIP concept will enable iProTech to evaluate the potential performance of this concept, and to investigate how sensitive the device is to different variables within the system (e.g. hydraulic components, geometry dimensions, etc.) 2) The numerical models were developed with WEC-Sim and PTO-Sim -- the de-facto industry standard WEC numerical modeling tools. A detailed model of the system's hydraulic circuit was developed with PTO-Sim. However, WEC-Sim is based on linear hydrodynamics and has several limitations; validation against high fidelity models/physical scale models would help to build confidence in the results.

16 TIDAL AND WAVE POWER↗

On the design of stable, consistent, and conservative high-order methods for multi-material hydrodynamics

Obtaining stable and high-order numerical solutions for multi-material hydrodynamics is an open challenge. Although slope limiters are widely used to maintain monotonicity near discontinuities, typical limiting procedures violate closure laws at the discrete level when applied to multi-material hydrodynamics equations. Due to this, the high-order expansions of quantities related by the closure laws are no longer consistent. The commonly observed symptom of this consistency-violation is that the numerical method fails to maintain constant pressure and velocity across material interfaces. This leads to sub-optimal convergence rates for smooth multi-material problems as well. Specialized limiting procedures that satisfy consistency while maintaining conservation need to be developed for such equations. A novel procedure that re-instates consistency into slope-limited high-order discretizations applied to the multi-material hydrodynamics equations is presented here. Using simple examples, it is demonstrated that the presented method satisfies closure laws at the discrete level, while maintaining conservative properties of the high-order method. Furthermore, this procedure involves a projection step which relies on the compact basis of the underlying spatial discretization, i.e. for discontinuous schemes (viz. DG and FV) the projection is local, and does not involve global matrix solves. Comparisons with conventional approaches emphasizes the necessity of the consistent closure-law preserving limiting approach, in order to maintain design order of accuracy for smooth multi-material problems.

36 MATERIALS SCIENCE↗

Nonequilibrium effects and transverse spherocity in ultrarelativistic proton-nucleus collisions

We investigate the effects of nonequilibrium dynamics in small colliding systems by comparing a nonequilibrium transport approach, Parton-Hadron-String Dynamics (PHSD), with a (2+1)-dimensional viscous hydrodynamic model, VISHNew. Focusing on p + Pb collisions at energies available at the Large Hadron Collider, we extract the initial conditions for the hydrodynamic model from PHSD, in order to reduce the impact of the early out-of-equilibrium dynamics and focus on the traces of nonequilibrium in the ensuing medium evolution. We find that in the transport approach quantities like energy density and bulk viscous pressure are highly inhomogeneous on the transverse plane during the whole evolution, whereas the hydrodynamic simulations dissolve more efficiently the initial spatial irregularities, still keeping a high degree of inhomogeneity due to the smaller size and lifetime of the medium produced in p + Pb collisions with respect to heavy-ion reactions. As a first step that will help to identify the impact of these nonequilibrium effects on final observables in proton-nucleus collisions, we perform an analysis of the transverse spherocity, an event-shape observable able to distinguish between jetty and isotropic configurations of transverse momenta. We found that the spherocity distribution in PHSD is slightly shifted towards the isotropic limit with respect to the hydrodynamic result. Even though this dissimilarity is partially due to the difference in the final charged particle production, it mainly comes from the different description within the two frameworks of the medium produced in small colliding systems. This finding supports the idea that multidifferential measurements, such as those based on event categorization according to multiplicity and spherocity, are useful to study final-state observables in ultrarelativistic proton-nucleus collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗