Search NASA⌕ Search

SEARCH · Search NASA

Results for “hydrodynamic coupling”

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 73 records · Page 4

Sea–Level Rise Impacts on Tidal Marshes and Estuarine Biogeochemical Processes

We used a numerical model to investigate the effects of sea-level rise (SLR) on the biogeochemical processes in the York River Estuary with extensive tidal marshes. The fully-coupled hydrodynamic-water quality-marsh model accounts for the spatial and temporal variations of physical-biogeochemical interactions between the tidal marshes and surrounding waters. This study focuses on an SLR scenario where the vertical accretion of tidal marshes keeps pace with the rising sea levels. Results show that SLR amplifies the tidal range and prolongs flooding duration, which results in enhanced porewater exchanges of materials between the tidal marshes and the surrounding waters. The increased availability of shallow-water habitats and enhanced light utilization in the shallow areas under SLR promote phytoplankton production (PP) in the shallow-water regions of the York River. Consequently, the organic carbon in the open water is fueled by the contributions from shallow waters and the enhanced export of organic carbon from the marshes. The change in the dissolved oxygen (DO) budget in the York River Estuary is attributed to changes in water column respiration, net metabolism of the benthic layer, reaeration, PP, and increased stratification under SLR. The net DO flux out of the York River increases at the York River mouth. Diel DO variation, especially in the marshes in the upper estuary, promotes phosphorus release from the sediment. The changes in dissolved nitrogen under SLR are relatively minimal.

54 ENVIRONMENTAL SCIENCES↗

Self-mixing in microtubule-kinesin active fluid from nonuniform to uniform distribution of activity

Active fluids have applications in micromixing, but little is known about the mixing kinematics of systems with spatiotemporally-varying activity. To investigate, UV-activated caged ATP is used to activate controlled regions of microtubule-kinesin active fluid and the mixing process is observed with fluorescent tracers and molecular dyes. At low Péclet numbers (diffusive transport), the active-inactive interface progresses toward the inactive area in a diffusion-like manner that is described by a simple model combining diffusion with Michaelis-Menten kinetics. At high Péclet numbers (convective transport), the active-inactive interface progresses in a superdiffusion-like manner that is qualitatively captured by an active-fluid hydrodynamic model coupled to ATP transport. Results show that active fluid mixing involves complex coupling between distribution of active stress and active transport of ATP and reduces mixing time for suspended components with decreased impact of initial component distribution. This work will inform application of active fluids to promote micromixing in microfluidic devices.

59 BASIC BIOLOGICAL SCIENCES↗

Correlations between asymmetric compression, burn amplification, and hot-spot velocities in inertial confinement fusion implosions

Here, this manuscript examines the correlations between the hot-spot velocity (an observable signature of residual kinetic energy), low-mode implosion asymmetries, and burn amplification in inertial confinement fusion implosions on the National Ignition Facility (NIF). Using a combination of two-dimensional axis-symmetric and three-dimensional radiation-hydrodynamic simulations coupled to neutronics, we find that for typical NIF implosions, the stagnation asymmetry multiplies the observed hot-spot velocity anywhere from 80% to 120%, while burn amplification always increases it. Additionally, we find stagnation asymmetry typically deflects the observed hot-spot flow. The two mechanisms (low-mode implosion asymmetries and burn amplification) can be decoupled, and application of a simple model to a database of cryogenic implosions on the NIF infers the total hot-spot velocity amplification. This finding modifies the interpretation of data collected from inertial confinement fusion experiments and impacts the magnitude and origin of low-mode asymmetries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Energy-conserving contact dynamics of nonspherical rigid-body particles

Understanding the contact dynamics of nonspherical particles is crucial for accurately modeling colloidal and granular systems where shape anisotropy dictates structural organization and transport properties. We here introduce an energy-conserving contact dynamics framework for arbitrary convex rigid-body particles by implementing vertex–boundary interactions in 2D and vertex–surface and edge–edge detection in 3D. The established formulation enables continuous force evaluation and prevents particle overlap while conserving total energy during translational and rotational motion. We demonstrate the framework’s stability and its utility to capture packing behavior, anisotropic diffusion, and equations of state of polygonal and polyhedral particles as examples. The framework establishes a robust and extensible foundation for investigating nonequilibrium dynamics of complex nonspherical particulate systems, enabling enhanced understanding of phenomena across spatiotemporal scales in self- and directed-assembly, granular flows, and hydrodynamics, potentially coupled with interactions that represent underlying physical mechanisms.

Discrete element method↗

Interpreting dynamic-compression experiments to uncover the time dependence of freezing: Application to gallium

Using pulsed-power magnetic field sources to compress gallium to gigapascal pressures on nanosecond timescales, we report here experiments on shockless dynamic compression of a liquid metal. Time-resolved velocimetry data reveal signatures of rapid freezing from a metastable liquid state, and we demonstrate that the kinetics of this nonequilibrium solidification can be accurately simulated with a computational modeling framework we have developed in previous studies, where classical nucleation theory is coupled with hydrodynamics. Notably, velocity traces in some of our experiments show evidence of a phase transition, while others do not, even though other types of evidence suggest that solidification may be occurring in all of them. We explain how predictions made by our models regarding the presence or absence of these phase-transition signatures motivated additional experiments that later confirmed the theoretical predictions. Our analysis shows that due to the rapid, quasi-isentropic nature of the loading path, our experiments were able to compress liquid gallium to metastable states that are undercooled below the equilibrium melt temperature by more than 300 K and exhibit pressures that approach five times the equilibrium melt pressure. The understanding gained in this study should form the basis for future dynamic-compression experiments aimed at interrogating melt curves at high pressures.

gallium↗

Heavy flavor quenching and flow: the roles of initial condition, pre-equilibrium evolution, and in-medium interaction

Within an advanced Langevin-hydrodynamics framework coupled to a hybrid fragmentation-coalescence hadronization model, we study heavy flavor quenching and flow in relativistic heavy-ion collisions. We investigate how the initial heavy quark spectrum, the in-medium energy loss and hadronization mechanisms of heavy quarks, the evolution profile of the pre-equilibrium stage, the medium flow, and the temperature dependence of heavy quark diffusion coefficients influence the suppression and elliptic flow of heavy mesons at the RHIC and the LHC. Our results show that the different modeling of initial conditions, pre-equilibrium evolution, and in-medium interactions can individually yield uncertainties of approximately 10-40% in D meson suppression and flow at a low transverse momentum. Furthermore, we also find that proper combinations of collisional versus radiative energy loss, coalescence versus fragmentation in hadronization, and the inclusion of medium flow are the most important factors for describing the suppression and elliptic flow of heavy mesons.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Reionization time of the Local Group and Local-Group-like halo pairs

ABSTRACT Patchy cosmic reionization resulted in the ionizing UV background asynchronous rise across the Universe. The latter might have left imprints visible in present-day observations. Several numerical simulation-based studies show correlations between the reionization time and overdensities and object masses today. To remove the mass from the study, as it may not be the sole important parameter, this paper focuses solely on the properties of paired haloes within the same mass range as the Milky Way. For this purpose, it uses CoDaII, a fully coupled radiation hydrodynamics reionization simulation of the local Universe. This simulation holds a halo pair representing the Local Group, in addition to other pairs, sharing similar mass, mass ratio, distance separation, and isolation criteria but in other environments, alongside isolated haloes within the same mass range. Investigations of the paired halo reionization histories reveal a wide diversity although always inside-out, given our reionization model. Within this model, haloes in a close pair tend to be reionized at the same time but being in a pair does not bring to an earlier time their mean reionization. The only significant trend is found between the total energy at z = 0 of the pairs and their mean reionization time: Pairs with the smallest total energy (bound) are reionized up to 50 Myr earlier than others (unbound). Above all, this study reveals the variety of reionization histories undergone by halo pairs similar to the Local Group, that of the Local Group being far from an average one. In our model, its reionization time is ∼625 Myr against 660 ± 4 Myr (z ∼ 8.25 against 7.87 ± 0.02) on average.

79 ASTRONOMY AND ASTROPHYSICS↗

Gamma-ray signatures of 𝑟-process radioactivity from the collapse of magnetized white dwarfs

We predict the gamma-ray line emission from 𝑟-process nuclei synthesized in the ejecta of the accretion-induced collapse (AIC) of a magnetized, rapidly rotating white dwarf. Using ejecta from a two-dimensional general-relativistic neutrino-magnetohydrodynamic simulation, further evolved with a radiation-hydrodynamics code coupled to an in situ nuclear reaction network, we construct angle-dependent gamma-ray spectra in the 0.01–10 MeV band via composition-dependent ray tracing through the ejecta. The emission between ∼1 and 10 d is dominated by 132 I (𝑡 1/2 = 2.3 h), continuously replenished by the decay of its parent 132 Te (𝑡 1/2 = 3.2 d), with additional contributions from 131 I, 133 Xe, and 132 Te . At 𝑡 ≳ 20 d, 56 Co (from 56 Ni decay) becomes the primary emitter. The simultaneous presence of 𝑟 process and iron-peak gamma-ray lines is distinctive of AIC ejecta and absent in binary neutron star mergers, where iron-peak nuclei are generally not synthesized. Comparing with the 3⁢𝜎 continuum sensitivities of planned MeV gamma-ray telescopes (COSI, AMEGO-X, e-ASTROGAM, GRAMS, GammaTPC), we find the brightest 𝑟-process lines detectable to ∼10 Mpc by GammaTPC and GRAMS, with the signal approaching their sensitivity threshold at 30 Mpc. As a result, the 𝑟-process spectral features survive time integration over ∼30 d exposures, demonstrating robustness against the long observation times required by gamma-ray detectors.

Nuclear reactions↗

Fast explicit solutions for neutrino-electron scattering: Explicit asymptotic methods

Here, we present results of explicit asymptotic approximations applied to neutrino-electron scattering in a representative model of neutrino population evolution under conditions characteristic of core-collapse supernova explosions or binary neutron star mergers. It is shown that this approach provides stable solutions of these stiff systems of equations, with accuracy and time stepping comparable to that for standard implicit treatments such as backward Euler, fixed point iteration, and Anderson-accelerated fixed point iteration. Because each time step can be computed more rapidly with the explicit asymptotic approximation than with implicit methods, this suggests that algebraically stabilized explicit integration methods could be used to compute neutrino evolution coupled to hydrodynamics more efficiently in stellar explosions and mergers than the methods currently in use.

79 ASTRONOMY AND ASTROPHYSICS↗

Quantifying Nitrogen Bioextraction by Seaweed Farms – A Real-time Modeling-Monitoring Case Study in Hood Canal, WA

Anthropogenic nitrogen input is the leading cause of increasing marine eutrophication. Cultivating seaweeds at large scales may provide a viable solution for mitigating eutrophication especially when nitrogen extraction is achieved as an added environmental benefit of biofuel and food production. As the first step to promoting and implementing this solution, it is critical to quantify and assess the effectiveness of nitrogen removal by existing seaweed farms. Pacific Northwest National Laboratory and Hood Canal Mariculture Inc. conducted a case study by applying a coupled modeling-monitoring approach to quantify nitrogen bioextraction by a one-hectare commercial sugar kelp farm in Hood Canal, WA. Through this two-year effort, the project team conducted intensive field monitoring of water quality parameters at the farm site and subsequent laboratory analysis. Meanwhile, the team applied a coupled hydrodynamic-macroalgal growth model to simulate kelp growth and nitrogen bioextraction over the two kelp growing seasons. The model results showed an overall good comparison with field monitoring data and showed the effectiveness of nitrogen removal by kelp farms. The established modeling-monitoring approach can also be transferable for supporting broader research activities in macroalgal farming and eutrophication mitigation.

09 BIOMASS FUELS↗

Assessing shellfish water exposure to fecal bacteria pollution in Salish Sea: three-dimensional modeling and implications for monitoring

Fecal bacteria (FB) contamination poses significant risks to shellfish safety and management in coastal and estuarine waters. Despite extensive pollution identification and correction efforts, FB contamination in shellfish-growing areas persists in the Salish Sea, highlighting the need to identify overlooked sources and better understand FB transport from riverine and shoreline inputs to shellfish beds. To address this, a high-resolution three-dimensional hydrodynamic model coupled with FB kinetics was developed and applied to a case study site in Salish Sea—Portage Bay—to simulate freshwater plume circulation, flushing dynamics, and bacterial transport. Daily FB loading from the major freshwater inflow—Nooksack River was generated by both linear interpolation and integrating a machine learning approach (XGBoost), trained on historical hydrological and meteorological data. The model successfully reproduced both the magnitude and seasonal variation of FB concentrations in Portage Bay for the year of 2021, demonstrating that simplified FB kinetics with first-order decay due to mortality was effective in this dynamic coastal environment with short flushing time. Model results identified the Nooksack River as the dominant far-field FB source, while scenario simulations showed that near-field coastal stormwater outfalls elevated local FB levels following rainfall, particularly under low-flow conditions. The XGBoost prediction provided comparable or superior accuracy to linear interpolation, particularly during periods of missing observational data, by capturing short-term variability and event-driven loading more effectively. Integrating data-driven riverine FB inputs with mechanistic coastal numerical modeling provides a robust framework for operational forecasting of shellfish bed exposure risk and supports adaptive monitoring and management of shellfish growing areas in the Salish Sea and similar coastal systems.

Salish Sea↗

thornado+FLASH-X: A Hybrid Discontinuous Galerkin–Implicit-explicit and Finite-volume Framework for Neutrino-radiation Hydrodynamics in Core-collapse Supernovae

We present neutrino-transport algorithms implemented in the toolkit for high-order neutrino-radiation hydrodynamics (thornado) and their coupling to self-gravitating hydrodynamics within the adaptive mesh refinement–based multiphysics simulation framework FLASH-X. thornado, developed primarily for simulations of core-collapse supernovae (CCSNe), employs a spectral, six-species two-moment formulation with algebraic closure and special-relativistic observer corrections accurate to $\mathcal{O}(v/c)$, and uses discontinuous Galerkin (DG) methods for phase-space discretization combined with implicit-explicit time stepping. A key development is a nonlinear neutrino–matter coupling algorithm based on nested fixed-point iteration with Anderson acceleration, enabling fully implicit treatment of collisional processes, including energy-coupling interactions such as neutrino–electron scattering and pair production. Coupling to finite-volume (FV) hydrodynamics is achieved through a hybrid DG-FV representation of the fluid variables and operator-split evolution within FLASH-X. The implementation is verified using basic transport tests with idealized opacities and relaxation and deleptonization problems with tabulated microphysics. Spherically symmetric CCSN simulations demonstrate accuracy and robustness of the coupled scheme, including close agreement with the CCSN simulation code Chimera. An axisymmetric CCSN simulation further demonstrates the viability of DG-based neutrino transport for multidimensional supernova modeling within FLASH-X. thornado’s neutrino-transport solver is GPU-enabled using OpenMP offloading or OpenACC, and all CCSN applications included in this work use the GPU implementation. Together, these results establish a foundation for future enhancements in physics fidelity, numerical algorithms, and computational performance, for increasingly realistic large-scale CCSN simulations.

Endeve, Eirik [Oak Ridge National Laboratory (ORNL↗

Coupling laser physics to radiation-hydrodynamics

In order to accurately model implosion hydrodynamics in a radiation-hydrodynamics code, it is essential to include accurate accounting for energy deposition physics. In inertial confinement fusion (ICF), where capsules are driven by lasers or laser-driven x-rays, energy deposition profiles and energy transport have a strong impact on the development and evolution of capsule dynamics and hydrodynamic instabilities. Nevertheless, accurately modeling laser beam propagation in radiation-hydrodynamics codes presents unique challenges associated with disparate resolution requirements, the potential to seed spurious noise in highly unstable systems, and computational expense. We discuss a new method for coupling laser ray-tracing physics to a radiation hydrodynamics code, developed in the process of implementing the Mazinisin laser ray-trace into the xRAGE radiation hydrodynamics code. In contrast to previous approaches, in which laser ray-tracing is performed on the radiation-hydrodynamics mesh, our method involves a mesh generation and evolution strategy that addresses the unique requirements of the laser ray-trace in a separate mesh, enabling performance enhancements and strategies to reduce noise seeded by the discretization of beams into computational rays. In addition, we have employed several methods to ensure that spurious mesh imprinting is minimized. These involved optimizing the laser and radiation-hydrodynamics meshes as well as interpolation between them and requires the use of an exact initialization method for the radiation-hydrodynamics mesh. Furthermore, these techniques have enabled efficient computation of laser-driven implosions and other experiments with minimal introduction of spurious noise.

47 OTHER INSTRUMENTATION↗

Coupled wave-current modeling for hydrodynamic load analysis of macroalgae cultivation farms

Large-scale cultivation of macroalgae is one of the most promising biofuel sources that could reduce our consumption of fossil fuels and would be particularly economically competitive when grown for the co-production of additional goods such as food and textiles. Current production of biofuel from algal biomass is limited by labor costs and lack of data on potential impacts of cultivation and harvesting on marine and coastal environments as well as a comprehensive characterization of the impact of environmental conditions on macroalgal feedstock and their bioproduct and biofuel yield and quality. Awarded by the United States (U.S.) Department of Energy Advanced Research Projects Agency-Energy (ARPA-E) Macroalgae Research Inspiring Novel Energy Resources (MARINER) program, which seeks to enable the U.S. as a global leader in the production of marine biomass, the Marine Biological Laboratory (MBL) has been developing a test system for tropical seaweed cultivation in the Gulf of Mexico and the Caribbean. An integral step in designing macroalgae cultivation farms is to adequately characterize the hydrodynamic climate at potential growth sites. In regions like the Gulf of Mexico and the Caribbean, which are highly exposed to tropical cyclones, it is critical to not only consider the day-to-day hydrodynamic conditions but to also assess the risk of extreme sea states to ensure the survivability of future macroalgae farms. Under these considerations, the Pacific Northwest National Laboratory (PNNL) is providing modeling support for MBL to inform the design and siting of the farm systems that have been proposed for their ARPA-E MARINER project. This report describes the development of a high-resolution coupled storm surge and wave model to simulate the hydrodynamics and wave climate at proposed macroalgae cultivation sites selected by MBL in Florida and Puerto Rico. Model results, including model validation, water level, current distributions, and sea states, are discussed for both selected sites in Florida and Puerto Rico coast. These simulations provide an accurate insight into the hydrodynamic conditions that a macroalgae farm is likely to experience during its operational lifetime, including current information to support fine-scale hydrodynamic load modeling, risk analysis and system design.

09 BIOMASS FUELS↗

From Hydrodynamics to Jet Quenching, Coalescence, and Hadron Cascade: A Coupled Approach to Solving the $R_{AA}Ⓧv_{2}$ Puzzle

Hydrodynamics and jet quenching are responsible for the elliptic flow $v_{2}$ and suppression of large transverse momentum ($p_{T}$) hadrons, respectively, two of the most important phenomena leading to the discovery of a strongly coupled quark-gluon plasma in high-energy heavy-ion collisions. A consistent description of the hadron suppression factor $R_{AA}$ and $v_{2}$, especially at intermediate $p_{T}$, however, remains a challenge. We solve this long-standing $R_{AA}Ⓧv_{2}$ puzzle by including quark coalescence for hadronization and final state hadron cascade in the coupled linear Boltzmann transport-hydro model that combines concurrent jet transport and hydrodynamic evolution of the bulk medium. We illustrate that quark coalescence and hadron cascade, two keys to solving the puzzle, also lead to a splitting of $v_{2}$ for pions, kaons, and protons in the intermediate $p_{T}$ region. We demonstrate for the first time that experimental data on $R_{AA}, v_{2}$, and their hadron flavor dependence from low to intermediate and high $p_{T}$ in high-energy heavy-ion collisions can be understood within this coupled framework.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

ALEGRA/Sceptre Code Coupling [Brief]

Researchers at Sandia have developed an advanced radiation hydrodynamic simulation capability by coupling the ALEGRA and Sceptre codes.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Exploring Flood Predictability in Taiwan through Coupled Atmospheric–Hydrological and High-Performance Hydrodynamic Models

Effective flood simulation capabilities can tremendously support early warning and disaster prevention. To examine the applicability of a fully physics-based and high-performance flood simulation and forecasting modeling framework for a flood-prone region in Taiwan, we conduct a numerical experiment that couples the Weather Research and Forecasting (WRF) Model, WRF-Hydrological modeling system (WRF-Hydro), and the Two-Dimensional Runoff Inundation Toolkit for Operational Needs (TRITON) to perform integrated rainfall, streamflow, and flood simulations. Furthermore, we first use the coupled WRF and WRF-Hydro (WWH) to predict rainfall and streamflow and then drive TRITON with the predicted streamflow hydrographs to simulate flood depth and inundation area. With the refined spatial resolution and parameterization, this framework can better predict rainfall with reasonable spatial patterns. Although WWH could overestimate the amount of rainfall in some areas, the uncertain rainfall–streamflow predictions produce reasonable flood maps able to pinpoint regions at risk of flooding. In terms of model efficiency, the graphics processing unit–based computation can yield a speed-up factor as high as ∼13 compared to the central processing unit–based computation, promoting the efficacy of the coupled modeling framework in practical real-time flood forecasting.

Coupled models↗