Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamics simulations”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Modeling of a spatially resolved ion temperature diagnostic for inertial confinement fusion

Here, the performance of modern laser-driven inertial confinement fusion (ICF) experiments is degraded by contamination of the deuterium–tritium (DT) fuel with high-Z material during compression. Simulations suggest that this mix can be described by the ion temperature distribution of the implosion, given that such contaminants deviate in temperature from the surrounding DT plasma. However, existing neutron time-of-flight (nTOF) diagnostics only measure the spatially integrated ion temperature. This paper describes the techniques and forward modeling used to develop a novel diagnostic imaging system to measure the spatially resolved ion temperature of an ICF implosion for the first time. The technique combines methods in neutron imaging and nTOF diagnostics to measure the ion temperature along one spatial dimension at yields currently achievable on the OMEGA laser. A detailed forward model of the source and imaging system was developed to guide instrument design. The model leverages neutron imaging reconstruction algorithms, radiation hydrodynamics and Monte Carlo simulations, optical ray tracing, and more. The results of the forward model agree with the data collected on OMEGA using the completed diagnostic. The analysis of the experimental data is still ongoing and will be discussed in a separate publication.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Quark flavor equilibration of the quark-gluon plasma

The early stage of a heavy-ion collision is marked by rapid entropy production and the transition from a gluon saturated initial condition to a plasma of quarks and gluons that evolves hydrodynamically. However, during the early times of the hydrodynamic evolution, the chemical composition of the QCD medium is still largely unknown. We present a study of quark chemical equilibration in the (Q)GP using a novel model of viscous hydrodynamic evolution in partial chemical equilibrium. Motivated by the success of gluon saturated initial condition models, we initialize the QCD medium as a completely gluon dominated state. Local quark production during the hydrodynamic phase is then simulated through the evolution of time-dependent fugacities for each independent quark flavor, with the timescales set as free parameters to compare different rates of equilibration. We present the results of complete heavy-ion collision simulations using this partial chemical equilibrium model, and show the effects on hadronic and electromagnetic observables. In particular, we show that the development of flow is sensitive to the equilibration timescale, providing an empirical way to probe the chemical equilibration of the QCD medium.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Second-harmonic generation in nonlinear plasmonic lattices enhanced by quantum emitter gain medium

In this work, we report on a theoretical study of second-harmonic generation (SHG) in plasmonic nanostructures interacting with two-level quantum emitters (QEs) under incoherent energy pump. We generalize the driven-dissipative Tavis–Cummings model by introducing the anharmonic surface plasmon-polariton (SPP) mode coupled to QEs and examine physical properties of corresponding SPP-QE polariton states. Our calculations of the SHG efficiency for strong QE-SPP coupling demonstrate orders of magnitude enhancement facilitated by the polariton gain. We further discuss time-domain numerical simulations of SHG in a square lattice comprising Ag nanopillars coupled to QEs utilizing a fully vectorial nonperturbative nonlinear hydrodynamic model for conduction electrons coupled to Maxwell–Bloch equations for QEs. The simulations support the idea of gain enhanced SHG and show orders of magnitude increase in the SHG efficiency as the QEs are tuned in resonance with the lattice plasmon mode and brought above the population inversion threshold by incoherent pumping. By varying pump frequency and tuning QEs to a localized plasmon mode, we demonstrate further enhancement of the SHG efficiency facilitated by strong local electric fields. The incident light polarization dependence of the SHG is examined and related to the symmetries of participating plasmon modes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Hydrodynamic characterization of the redox chemistry of crown-encapsulated uranyl complexes

The redox properties of actinide-containing species strongly influence their reactivity, speciation, and interfacial behavior, but the experimental quantification of the electrochemical characteristics of molecular actinide complexes in nonaqueous media has not received the attention it deserves. Here, results from hydrodynamic methods and electrochemical simulations of U(VI)/U(V) redox are reported, including quantification of heterogeneous electron-transfer kinetics and estimation of chemical reversibility of U(VI)/U(V) interconversion at electrodes in acetonitrile-based electrolyte. The complexes investigated are recently reported U(VI) and U(V) complexes in which the uranyl ion (UO 2 n+ ) is encapsulated in a macrocyclic 18-crown-6-like moiety templated by a Pt(II) center. These complexes feature the most positive value U VI /U V reduction potential yet reported and are thus particularly relevant to study of facile U(V) generation from U(VI) precursors as well as uranium electroanalysis. Rotating disk electrode (RDE) studies have been used to quantify the diffusion coefficients of the U(VI) and U(V) complexes, and standard heterogeneous electron transfer rate constants ( k 0 ) for the redox have been determined using a conventional Koutecký-Levich analysis. Rotating ring-disk electrode (RRDE) studies have been used to directly interrogate the chemical reversibility of U(VI)-U(V) interconversion, confirming that reduction of the U(VI) complex at an Au disk is associated with formation of the U(V) analogue that can be readily re-oxidized at a Pt ring under hydrodynamic (rotating) conditions. Because measurements of the type reported here are generally associated with current flows that are larger than those found in corresponding quiescent (unstirred) conditions, our findings suggest that hydrodynamic methods could be advantageous for design of electroanalytical approaches to detection of actinide species and study of their redox properties.

actinides↗
DOE OSTI2024

Effects of stratification on overshooting and waves atop the convective core of M ⊙ main-sequence stars

As a massive star evolves along the main sequence, its core contracts, leaving behind a stable stratification in helium. We simulate two-dimensional convection in the core at three different stages of evolution of a $5\,\mathrm{ M}_{\odot }$ star, with three different stratifications in helium atop the core. We study the propagation of internal gravity waves in the stably stratified envelope, along with the overshooting length of convective plumes above the convective boundary. We find that the stratification in helium in evolved stars hinders radial motions and effectively shields the radiative envelope against plume penetration. This prevents convective overshooting from being an efficient mixing process in the radiative envelope. In addition, internal gravity waves are less excited in evolved models compared to the zero-age-main-sequence model, and are also more damped in the stratified region above the core. As a result, the wave power is several orders of magnitude lower in mid- and terminal-main-sequence models compared to zero-age-main-sequence stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Modeling cross-beam energy transfer with sector ray tracing

Ray-based cross-beam energy transfer (CBET) models are an essential feature of the radiation-hydrodynamic codes used to simulate inertial confinement fusion implosions, but full 3D ray-based CBET calculations can have a prohibitively high computational cost. Sector ray tracing can be used to reduce the cost by orders of magnitude in cases where the coronal plasma and laser drive can be approximated as spherically symmetric. An extension of sector ray tracing (section ray tracing) can be used to relax the assumption of a spherically symmetric drive while still retaining most of the computational savings of a sector ray trace. We discuss the foundations of sector and section ray tracing and compare them to full ray tracing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A continuum meshfree method for sandbox-style numerical modeling of accretionary and doubly vergent wedges

Here, in this study, we introduce and utilize GEOSPH, a Lagrangian particle-based continuum Smoothed Particle Hydrodynamics (SPH) code to simulate the deformation of geologic materials in sandbox-style experiments using an elastoplastic constitutive model. Due to the meshfree nature of SPH, our proposed framework is capable of handling the large deformations characteristic of sandbox-style experiments. SPH is also able to capture strain localization and discontinuities, allowing us to resolve faulting patterns emerging in the sandbox experiments. We successfully validate our SPH implementation against the benchmark experiments of Buiter and co-workers, and then apply SPH to understand the deformational processes occurring in doubly vergent orogens resultant from the S-point experimental setup and from finitely high wedge-like indenter backstops of different inclinations. We find that the topographic slope of the prowedge results from the competing processes of vertical uplift due to the indenter and forward thrusting, while the topographic slope of the retrowedge is a product of repeated shallow slope failure. Our findings show that SPH is a promising numerical method which should be added to the current palette of computational tools used by structural geologists for solving problems relating to tectonics and crustal deformation.

58 GEOSCIENCES↗
DOE OSTI2023

Effects of Different Closure Choices in Core-collapse Supernova Simulations

The two-moment method is widely used to approximate the full neutrino transport equation in core-collapse supernova (CCSN) simulations, and different closures lead to subtle differences in the simulation results. In this paper, we compare the effects of closure choices on various physical quantities in 1D and 2D time-dependent CCSN simulations with our multigroup radiation hydrodynamics code FORNAX. We find that choices of the third-order closure relations influence the time-dependent simulations only slightly. Choices of the second-order closure relation have larger consequences than choices of the third-order closure, but these are still small compared to the remaining variations due to ambiguities in some physical inputs such as the nuclear equation of state. We also find that deviations in Eddington factors are not monotonically related to deviations in physical quantities, which means that simply comparing the Eddington factors does not inform one concerning which closure is better.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Performance analysis and tank test validation of a hybrid ocean wave-current energy converter with a single power takeoff

This paper introduces a hybrid ocean wave-current energy converter (HWCEC) that harvests energy from ocean waves and current simultaneously with a single power takeoff (PTO). Specifically, the wave energy is extracted through relative heaving motion between a floating buoy and a submerged second body, while the current energy is extracted using a marine current turbine (MCT). Energy from both sources are integrated by a hybrid PTO whose concept is based on a mechanical motion rectifier (MMR). The hybrid PTO with two one-way clutches converts bidirectional, up-and-down motion from the waves into unidirectional rotation of the generator. Meanwhile a third one-way clutch couples the MCT with the same generator. The wave and current can simultaneously or separately drive the same generator through different engagement and disengagement statuses of the one-way clutches. Time-domain simulation is conducted with hydrodynamic coefficients obtained from computational fluid dynamics analysis and boundary element method. Tank tests were conducted for a HWCEC under co-existing wave and current inputs. For comparison, separate baseline tests of a turbine and a two-body point absorber, each acting in isolation, are conducted. Experimental results validate the dynamic modeling and show that a HWCEC can increase the output power with a range between 29-87% over either current turbine and wave energy converter acting individually, and it can reduce by up to 70% the peak-to-average power ratio compared with the wave energy converter on the tested conditions.

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

Evidence of a toroidal magnetic field in the core of 3C 84

The spatial scales of relativistic radio jets, probed by relativistic magneto-hydrodynamic (RMHD) jet launching simulations and by most very long baseline interferometry (VLBI) observations differ by an order of magnitude. Bridging the gap between these RMHD simulations and VLBI observations requires selecting nearby active galactic nuclei (AGN), the parsec-scale region of which can be resolved. The radio source 3C 84 is a nearby bright AGN fulfilling the necessary requirements: it is launching a powerful, relativistic jet powered by a central supermassive black hole, while also being very bright. Using 22 GHz globe-spanning VLBI measurements of 3C 84 we studied its sub-parsec region in both total intensity and linear polarisation to explore the properties of this jet, with a linear resolution of ~0.1 parsec. We tested different simulation set-ups by altering the bulk Lorentz factor Γ of the jet, as well as the magnetic field configuration (toroidal, poloidal, helical). We confirm the persistence of a limb brightened structure, which reaches deep into the sub-parsec region. The corresponding electric vector position angles (EVPAs) follow the bulk jet flow inside but tend to be orthogonal to it near the edges. Our state-of-the-art RMHD simulations show that this geometry is consistent with a spine-sheath model, associated with a mildly relativistic flow and a toroidal magnetic field configuration.

3C 84 (NGC 1275)↗
DOE OSTI2020

Dark matter cores and cusps in spiral galaxies and their explanations

We compare proposed solutions to the core vs cusp issue of spiral galaxies, which has also been framed as a diversity problem, and demonstrate that the cuspiness of dark matter halos is correlated with the stellar surface brightness. We compare the rotation curve fits to the SPARC sample from a self-interacting dark matter (SIDM) model, which self-consistently includes the impact of baryons on the halo profile, and hydrodynamical N-body simulations with cold dark matter (CDM). The SIDM model predicts a strong correlation between the core size and the stellar surface density, and it provides the best global fit to the data. The CDM simulations without strong baryonic feedback effects fail to explain the large dark matter cores seen in low surface brightness galaxies. On the other hand, with strong feedback, CDM simulations do not produce galaxy analogs with high stellar and dark matter densities, and therefore they have trouble in explaining the rotation curves of high surface brightness galaxies. This implies that current feedback implementations need to be modified. We also explicitly show how the concentration-mass and stellar-to-halo mass relations together lead to a radial acceleration relation (RAR) in an averaged sense, and reiterate the point that the RAR does not capture the diversity of galaxy rotation curves in the inner regions. These results make a strong case for SIDM as the explanation for the cores and cusps of field galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Baryonic effects in the Effective Field Theory of Large-Scale Structure and an analytic recipe for lensing in CMB-S4

Upcoming Large-Scale Structure surveys will likely become the next leading sources of cosmological information, making it crucial to have a precise understanding of the influence of baryons on cosmological observables. The Effective Field Theory of Large-Scale Structure (EFTofLSS) provides a consistent way to predict the clustering of dark matter and baryons on large scales, where their leading corrections in perturbation theory are given by a simple and calculable functional form even after the onset of baryonic processes. In this paper, we extend the two-fluid-like system up to two-loop order in perturbation theory. Along the way, we show that a new linear counterterm proportional to the relative velocity of the fluids could generically be present, but we show that its effects are expected to be small in our universe. Regardless, we show how to consistently perform perturbation theory in the presence of this new term. We find that the EFTofLSS at two-loop order can accurately account for the details of baryonic processes on large scales. Additionally, we compare our results to a hydrodynamical N-body simulation at many redshifts and find that the counterterms associated with the leading corrections to dark matter and baryons start to differ between redshifts z ≈ 3 and z ≈ 2, signaling the onset of star-formation physics. We then use these fits to compute the lensing power spectrum, show that the understanding of baryonic processes will be important for analyzing CMB-S4 data, and show that the two-loop EFTofLSS accurately captures these effects for ℓ ≲ 2000. Our results are also potentially of interest for current and future weak lensing surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Collapse of magnetized white dwarfs as site of heavy-element formation and kilonova signal

We present the first end-to-end calculation connecting the accretion-induced collapse (AIC) of a magnetized, rapidly rotating white dwarf to observable kilonova signatures, combining two-dimensional (2D) general-relativistic neutrino-magnetohydrodynamic simulations, followed by radiation hydrodynamics with in-situ nuclear network and 2D Monte Carlo radiative transfer with spatially resolved heating rates. Unlike all previous unmagnetized AIC models – which predicted proton-rich, $^{56}$Ni-dominated ejecta – strong magnetic fields eject ${\approx }\, 0.2\, \mathrm{ M}_\odot$ of neutron-rich material ($\langle Y_e \rangle \sim 0.24$) on dynamical time-scales, before neutrino irradiation can raise the electron fraction, enabling strong r-process nucleosynthesis up to and beyond the third peak. The resulting kilonova is lanthanide-rich ($X_{\rm lan} \approx 8~{{\ \rm per\ cent}}$) and dominated by near-infrared emission. We compute synthetic light curves in the Large Synoptic Survey Telescope and J ames Webb Space Telescope bands and find striking agreement, without parameter tuning, between the observations of AT 2023vfi/GRB 230307A and our broadband light curves for polar viewing angles. These results establish magnetized AIC as a viable channel for heavy r-process element production and a compelling progenitor candidate for long-duration gamma-ray bursts with kilonova signatures.

MHD↗
DOE OSTI2022

The local wavenumber model for computation of turbulent mixing

We present an overview of the current status in the development of a two-point spectral closure model for turbulent flows, known as the local wavenumber (LWN) model. The model is envisioned as a practical option for applications requiring multi-physics simulations in which statistical hydrodynamics quantities such as Reynolds stresses, turbulent kinetic energy, and measures of mixing such as density-correlations and mix-width evolution, need to be captured with relatively high fidelity. In this review, we present the capabilities of the LWN model since it was first formulated in the early 1990s, for computations of increasing levels of complexity ranging from homogeneous isotropic turbulence, inhomogeneous and anisotropic single-fluid turbulence, to two-species mixing driven by buoyancy forces. The review concludes with a discussion of some of the more theoretical considerations that remain in the development of this model.

42 ENGINEERING↗
DOE OSTI2022

Critical lensing and kurtosis near a critical point in the QCD phase diagram in and out of equilibrium

In this work, we study the lensing effect of the QCD critical point on hydrodynamic trajectories and its consequences on the net-proton kurtosis $\kappa$ 4 . Including critical behavior by means of the BEST Collaboration equation of state (EoS), we first consider a scenario in equilibrium, then compare with hydrodynamic 0+1⁢D simulations with Bjorken expansion, including both shear and bulk viscous terms. Here, we find that, both in and out of equilibrium, the size and shape of the critical region directly affect if the signal will survive through the dynamical evolution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Future climate doubles the risk of hydraulic failure in a wet tropical forest

Summary Future climate presents conflicting implications for forest biomass. We evaluate how plant hydraulic traits, elevated CO 2 levels, warming, and changes in precipitation affect forest primary productivity, evapotranspiration, and the risk of hydraulic failure. We used a dynamic vegetation model with plant hydrodynamics (FATES‐HYDRO) to simulate the stand‐level responses to future climate changes in a wet tropical forest in Barro Colorado Island, Panama. We calibrated the model by selecting plant trait assemblages that performed well against observations. These assemblages were run with temperature and precipitation changes for two greenhouse gas emission scenarios (2086–2100: SSP2‐45, SSP5‐85) and two CO 2 levels (contemporary, anticipated). The risk of hydraulic failure is projected to increase from a contemporary rate of 5.7% to 10.1–11.3% under future climate scenarios, and, crucially, elevated CO 2 provided only slight amelioration. By contrast, elevated CO 2 mitigated GPP reductions. We attribute a greater variation in hydraulic failure risk to trait assemblages than to either CO 2 or climate. Our results project forests with both faster growth (through productivity increases) and higher mortality rates (through increasing rates of hydraulic failure) in the neo‐tropics accompanied by certain trait plant assemblages becoming nonviable.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Efficient Probabilistic Prediction and Uncertainty Quantification of Tropical Cyclone–Driven Storm Tides and Inundation

Abstract This study proposes and assesses a methodology to obtain high-quality probabilistic predictions and uncertainty information of near-landfall tropical cyclone–driven (TC-driven) storm tide and inundation with limited time and resources. Forecasts of TC track, intensity, and size are perturbed according to quasi-random Korobov sequences of historical forecast errors with assumed Gaussian and uniform statistical distributions. These perturbations are run in an ensemble of hydrodynamic storm tide model simulations. The resulting set of maximum water surface elevations are dimensionality reduced using Karhunen–Loève expansions and then used as a training set to develop a polynomial chaos (PC) surrogate model from which global sensitivities and probabilistic predictions can be extracted. The maximum water surface elevation is extrapolated over dry points incorporating energy head loss with distance to properly train the surrogate for predicting inundation. We find that the surrogate constructed with third-order PCs using elastic net penalized regression with leave-one-out cross validation provides the most robust fit across training and test sets. Probabilistic predictions of maximum water surface elevation and inundation area by the surrogate model at 48-h lead time for three past U.S. landfalling hurricanes (Irma in 2017, Florence in 2018, and Laura in 2020) are found to be reliable when compared to best track hindcast simulation results, even when trained with as few as 19 samples. The maximum water surface elevation is most sensitive to perpendicular track-offset errors for all three storms. Laura is also highly sensitive to storm size and has the least reliable prediction. Significance Statement The purpose of this study is to develop and evaluate a methodology that can be used to provide high-quality probabilistic predictions of hurricane-induced storm tide and inundation with limited time and resources. This is important for emergency management purposes during or after the landfall of hurricanes. Our results show that sampling forecast errors using quasi-random sequences combined with machine learning techniques that fit polynomial functions to the data are well suited to this task. The polynomial functions also have the benefit of producing exact sensitivity indices of storm tide and inundation to the forecasted hurricane properties such as path, intensity, and size, which can be used for uncertainty estimation. The code implementing the presented methodology is publicly available on GitHub.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Interplay between gain and loss in arrays of nonlinear plasmonic nanoparticles: toward parametric downconversion and amplification

With the help of a theoretical model and finite-difference time-domain (FDTD) simulations based on the hydrodynamic–Maxwell model, we examine the effect of difference-frequency generation (DFG) in an array of L-shaped metal nanoparticles (MNPs) characterized by intrinsic plasmonic nonlinearity. The outcomes of the calculations reveal the spectral interplay between gain and loss in the vicinity of the fundamental frequency of the localized surface plasmon resonances. Subsequently, we identify different array thicknesses and pumping regimes facilitating parametric amplification and spontaneous parametric downconversion. Our results suggest that the parametric amplification regime becomes feasible on a scale of hundreds of nanometers and spontaneous parametric downconversion on the scale of tens of nanometers, opening up new exciting opportunities for developing building blocks of photonic metasurfaces.

Shah, Syed A.↗