Modeling Modulated Tool Path Machining Using Coupled Smoothed Particle Hydrodynamics and Finite Element Method
Explore the source record for details and available documents.
Engineering, Materials science↗
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.
A coupled smoothed particle hydrodynamics-finite volume approach for shock capturing in one-dimension
Not Available
42 ENGINEERING↗
A Coupled Smoothed Particle Hydrodynamics-Finite Volume Approach for Shock Capturing in One-Dimension.
Abstract not provided.
Myers, Connor↗
High-Fidelity Simulations of Objects in Accretion Disks by Coupling REBOUND to a Hydrodynamics Code
Explore the source record for details and available documents.
Dempsey, Adam Mahlon↗
Characteristic momentum of Hydro + and a bound on the enhancement of the speed of sound near the QCD critical point
Near the critical point in the QCD phase diagram, hydrodynamics breaks down at a momentum where the frequency of the fastest hydrodynamic mode becomes comparable with the decay rate of the slowest nonhydrodynamic mode. Hydro+ was developed as a framework which extends the range of validity of hydrodynamics beyond that momentum value. This was achieved through coupling the hydrodynamic modes to the slowest nonhydrodynamic mode. In this work, analyzing the spectrum of linear perturbations in Hydro+, we find that a slow mode falls out of equilibrium if its momentum is greater than a characteristic momentum value. That characteristic momentum turns out to be set by the branch points of the dispersion relations. These branch points occur at the critical momenta of so-called spectral curves and are related to the radius of convergence of the derivative expansion. The existence of such a characteristic momentum scale suggests that a particular class of slow modes has no remarkable effect on the flow of the plasma. Based on these results, we show that there is an enhancement of the speed of sound due to the Hydro+ corrections compared to the known thermodynamic corrections near the critical point. This enhancement can not exceed a temperature-dependent upper bound which we derive near the critical point in the QCD phase diagram.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Impact of in situ nuclear networks and atomic opacities on neutron star merger ejecta dynamics, nucleosynthesis, and kilonovae
Context. Binary neutron star merger (BNSM) ejecta are key sites of rapid neutron capture (r-process) nucleosynthesis and they produce kilonovae powered by the radioactive decay of freshly synthesized nuclei. Modeling their evolution requires multi-physics simulations involving hydrodynamics, nuclear reactions, and radiative processes. The impact of nuclear burning and atomic opacity is poorly understood and often treated with simplified prescriptions. Aims. We systematically investigate different treatments of nuclear heating, particle thermalization, and atomic opacities in radiation-hydrodynamics simulations of BNSM ejecta and kilonova light curves. Methods. Ejecta profiles from long-term numerical-relativity simulations of asymmetric neutron star binaries with a massive neutron star remnant were evolved to ∼30 days using a 2D ray-by-ray approach. We compared simplified heating-rate and thermalization prescriptions with in situ Nuclear reaction Network (NN) calculations that track nuclear energy deposition and include a composition-dependent thermalization scheme. We also contrasted various gray opacity models with a frequency-dependent treatment based on atomic calculations. Results. Coupling NN and hydrodynamics significantly affects nucleosynthesis and kilonova emission. Assuming homologous expansion alters abundance evolution and produces a narrower, less populated second r-process peak and a third peak shifted to higher mass numbers. The back-reaction of nuclear heating affects the temperature evolution enough to delay and redden the early (t∼ hours) kilonova peaks. A constant thermalization efficiency underestimates and reddens the early emission while overestimating the late-time luminosity compared to the composition-dependent treatment. Analytical opacity prescriptions yield a more extended, colder photosphere, resulting in dimmer, redder kilonovae at early times (t≲ hour), while the delayed recession of the photosphere prolongs the red emission at t ≳ 5 days. Conclusions. Coupling hydrodynamics to an in situ NN is crucial for reliable nucleosynthesis and kilonova predictions. Resolving the first several hundred milliseconds of the hydrodynamics is essential for robust nucleosynthesis calculations. Composition-dependent thermalization and frequency-dependent, atomic-physics-based opacities are needed to accurately capture the temperature evolution of the ejecta and the brightness and color evolution of the kilonova. Calibrated analytic nuclear-power fits with simplified thermalization and opacity prescriptions can still reproduce the density and temperature evolution of the ejecta.
74 ATOMIC AND MOLECULAR PHYSICS↗
Discretization Writeup for Grey Flux-Limited Radiation Diffusion
This report documents the time and space discretizations for grey flux-limited diffusion applied to the thermal radiative transfer (TRT) equations. We begin with a description of the physics being solved before moving into the diffusion approximation. Once we have the TRT system, we show a finite-volume-inspired discretization from Jim Morel (Texas A&M University, NUEN 627 class notes, lecture 8). As systems become hotter they emit more photons in the form of blackbody radiation. Because average photon energy of the blackbody source is proportional to the temperature of the system, we call these thermal photons or thermal radiation. As material temperatures increase, increasing fractions of the total energy in the system go into the radiation field. In addition, radiation can deposit energy and momentum non-locally, making it an important phenomenon for heating and impulse. In order to accurately study systems at high temperatures, we wish to add the physics of thermal radiation to our hydrodynamic system. In practice, coupling radiation and hydrodynamics is often done by operator-splitting each timestep into two consecutive, non-overlapping phases: (1) update the hydrodynamics for a fixed radiation state (2) update the radiation and internal energy for an otherwise fixed hydrodynamic state. Because of this clean separation of physics updates, in this report we show only the latter phase, which involves solely the TRT equations.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
A dissolution model of alite coupling surface topography and ions transport under different hydrodynamics conditions at microscale
Portland cement is the most produced material in the world. The hydration process of cement consists of a group of complex chemical reactions. In order to investigate the mechanism of cement hydration, it is vital to study the hydration of each phase separately. An integrated model is proposed in this paper to simulate the dissolution of alite under different hydrodynamic conditions at microscale, coupling Kinetic Monte Carlo model (KMC), Lattice Boltzmann method (LBM) and diffusion boundary layer (DBL). The dissolution of alite is initialised with KMC. Two Multiple-relaxation-time (MRT) LB models are used to simulate the fluid flow and transport of ions, respectively. For solid-liquid interface, DBL is adapted to calculate the concentration gradient and dissolution flux. The model is validated with experiment from literature. The simulation results show good agreements with the results published in the literature.
36 MATERIALS SCIENCE↗
Sensitivity of Simulations of Double-detonation Type Ia Supernovae to Integration Methodology
Abstract We study the coupling of hydrodynamics and reactions in simulations of the double-detonation model for Type Ia supernovae. When assessing the convergence of simulations, the focus is usually on spatial resolution; however, the method of coupling the physics together as well as the tolerances used in integrating a reaction network also play an important role. In this paper, we explore how the choices made in both coupling and integrating the reaction portion of a simulation (operator/Strang splitting versus the simplified spectral deferred corrections method we introduced previously) influences the accuracy, efficiency, and nucleosynthesis of simulations of double detonations. We find no need to limit reaction rates or reduce the simulation time step to the reaction timescale. The entire simulation methodology used here is GPU-accelerated and made freely available as part of the Castro simulation code.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Evaluating the Effectiveness of Soil Profile Rehabilitation for Pluvial Flood Mitigation Through Two-Dimensional Hydrodynamic Modeling
Pluvial flooding, driven by increasingly impervious surfaces and intense storm events, presents a growing challenge for urban areas worldwide. In Baltimore City, MD, USA, climate change, rapid urbanization, and aging stormwater infrastructure are exacerbating flooding impacts, resulting in significant socio-economic consequences. This study evaluated the effectiveness of a soil profile rehabilitation scenario using a 2D hydrodynamic modeling approach for the Tiffany Run watershed, Baltimore City. This study utilized different extreme storm events, a high-resolution (1 m) LiDAR Digital Terrain Model (DTM), building footprints, and hydrological soil data. These datasets were integrated into a fully coupled 2D hydrodynamic model, the City Catchment Analysis Tool (CityCAT), to simulate urban flood dynamics. The pre-soil rehabilitation simulation revealed a maximum water depth of 3.00 m in most areas, with hydrologic soil groups C and D, especially downstream of the study area. The post-soil rehabilitation simulation was targeted at vacant lots and public parcels, accounting for 33.20% of the total area of the watershed. This resulted in a reduced water depth of 2.50 m. Additionally, the baseline runoff coefficient of 0.49 decreased to 0.47 following the rehabilitation, and the model consistently recorded a peak runoff reduction rate of 4.10 across varying rainfall intensities. The validation using a contingency matrix demonstrated true-positive rates of 0.75, 0.50, 0.64, and 0 for the selected events, confirming the model’s capability at capturing real-world flood occurrences.
Baltimore City↗
Pulsed power to control growth of silicon nanoparticles in low temperature flowing plasmas
Low-temperature plasmas have seen increasing use for synthesizing high-quality, mono-disperse nanoparticles (NPs). Recent work has highlighted that an important process in NP growth in plasmas is particle trapping—small, negatively charged nanoparticles become trapped by the positive electrostatic potential in the plasma, even if only momentarily charged. In this article, results are discussed from a computational investigation into how pulsing the power applied to an inductively coupled plasma (ICP) reactor may be used for controlling the size of NPs synthesized in the plasma. The model system is an ICP at 1 Torr to grow silicon NPs from an Ar/SiH 4 gas mixture. This system was simulated using a two-dimensional plasma hydrodynamics model coupled to a three-dimensional kinetic NP growth and trajectory tracking model. The effects of pulse frequency and pulse duty cycle are discussed. We identified separate regimes of pulsing where particles become trapped for one pulsed cycle, a few cycles, and many cycles—each having noticeable effects on particle size distributions. For the same average power, pulsing can produce a stronger trapping potential for particles when compared to continuous wave power, potentially increasing particle mono-dispersity. Pulsing may also offer a larger degree of control over particle size for the same average power. Experimental confirmation of predicted trends is discussed.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Non-local transport in radiation-hydrodynamics codes for ICF by efficient coupling to an external Vlasov–Fokker–Planck code
Accurately incorporating non-local transport into radiation-hydrodynamics codes, and indeed any fluid system, has long been elusive. To date, a simplified and accurate theory that can be easily integrated has not been available. This limitation affects modeling in inertial confinement fusion (ICF) and magnetic confinement fusion systems, among others, where non-local transport is well-known to be present. Here, we present a coupling methodology between a full Vlasov–Fokker–Planck (VFP) electron kinetic code and radiation-hydrodynamics (rad-hydro) codes. The VFP code is used to adjust native electron transport in the rad-hydro code, thus enabling improved transport without the need to integrate a full electron VFP solver into the rad-hydro code. This approach necessitates only occasional invocation of the VFP code, reducing computational intensity compared to following the dynamic evolution entirely with the VFP code on fluid time scales. We illustrate that the methodology is more accurate than other simplified methods in thermal decay systems relevant to ICF and can replicate standard theoretical results with high accuracy.
Electronic transport↗
Extended experimental inferential structure determination method in determining the structural ensembles of disordered protein states
Proteins with intrinsic or unfolded state disorder comprise a new frontier in structural biology, requiring the characterization of diverse and dynamic structural ensembles. Here we introduce a comprehensive Bayesian framework, the Extended Experimental Inferential Structure Determination (X-EISD) method, which calculates the maximum log-likelihood of a disordered protein ensemble. X-EISD accounts for the uncertainties of a range of experimental data and back-calculation models from structures, including NMR chemical shifts, J-couplings, Nuclear Overhauser Effects (NOEs), paramagnetic relaxation enhancements (PREs), residual dipolar couplings (RDCs), hydrodynamic radii (R h ), single molecule fluorescence Förster resonance energy transfer (smFRET) and small angle X-ray scattering (SAXS). We apply X-EISD to the joint optimization against experimental data for the unfolded drkN SH3 domain and find that combining a local data type, such as chemical shifts or J-couplings, paired with long-ranged restraints such as NOEs, PREs or smFRET, yields structural ensembles in good agreement with all other data types if combined with representative IDP conformers.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
HYDRODYNAMIC TESTING OF ZIRCALOY TO STEEL HUB COUPLINGS FOR USE IN ADVANCED TEST REACTOR I-LOOP TUBE
The ATR I-positions' isolation from the flux trap region provides an opportunity for a test facility that could create conditions representative of commercial light water reactors. Such a facility requires piping couplings and joints that withstand extreme thermal and pressure cycling. Also, transitioning from steel to zircaloy piping inside the reactor vessel allows the neutronic properties of zircaloy to be leveraged while using cost-effective stainless steel outside of the core. Grayloc-style hub couplings will act as transition points between stainless steel piping and zircaloy piping inside the reactor vessel. These couplings have not been sufficiently tested for leakage where there is a transition between materials. This transition is a concern because the thermal expansion coefficients of the materials differ. Seal integrity between the two components may be greatly affected by thermal expansion at the connection. Thus, the effect of thermal cycling must be assessed. Thermal and pressure cycling and hydrodynamic testing will be conducted to verify seal integrity of the steel to zircaloy hub coupling when subjected to ATR conditions.
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Recent progress on coarse graining simulations
We focus on coarse graining simulations based on the primary conservation equations, effectively codesigned physics and algorithms, and low-Mach-number corrected (LMC) hydrodynamics. Simulation methods involve LANL’s x-Radiation-Adaptive-Grid-Eulerian Large-Eddy Simulation, Besnard-Harlow-Rauenzahn (BHR) Reynolds-Averaged Navier-Stokes (RANS) approach, and Dynamic BHR – a paradigm bridging RANS and LES. A relevant question addressed relates to whether 3D RANS and RANS/LES hybrids – the industry standards for aerospace and automotive research, are presently relevant for practical variable-density applications involving shocked and accelerated interface instabilities. Furthermore, recent simulations of the GaTECH inclined mixing-layer shock-tube and NIF ICF-capsule experiments are used to demonstrate issues, challenges, and potential for 3D coarse grained LMC simulation strategies for robustly simulating complex transitional and coupled hydrodynamics-multiphysics with coarser resolution. Present LES readiness to provide accurate predictions at scale is demonstrated – whereas 3D RANS and RANS/LES bridging do not appear impactful in this context.
42 ENGINEERING↗
Anisotropic flow of identified hadrons in $\mathrm{X}$e-$\mathrm{X}$e collisions at $\sqrt{s_{NN}}$ = 5.44 $\mathrm{T}$e$\mathrm{V}$
Measurements of elliptic (v 2 ) and triangular (v 3 ) flow coefficients of π ± , K ± , $p+\bar{p}$, $K$$^{0}_{S}$, and $Λ+\bar{Λ}$ obtained with the scalar product method in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV are presented. The results are obtained in the rapidity range |y| < 0.5 and reported as a function of transverse momentum, p T , for several collision centrality classes. The flow coefficients exhibit a particle mass dependence for p T < 3 GeV/c, while a grouping according to particle type (i.e., meson and baryon) is found at intermediate transverse momenta (3 < p T < 8 GeV/c). The magnitude of the baryon v 2 is larger than that of mesons up to p T = 6 GeV/c. The centrality dependence of the shape evolution of the p T -differential v 2 is studied for the various hadron species. The v 2 coefficients of π ± , K ± , and $p+\bar{p}$ are reproduced by MUSIC hydrodynamic calculations coupled to a hadronic cascade model (UrQMD) for p T < 1 GeV/c. A comparison with v n measurements in the corresponding centrality intervals in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV yields an enhanced v 2 in central collisions and diminished value in semicentral collisions.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Heavy and light flavor jet quenching in heavy-ion collisions in a perturbative QCD approach
In this work, we study the flavor dependence of jet quenching in heavy-ion collisions. A next-to-leading-order perturbative QCD framework is applied to account for both quark and gluon contributions to light and heavy flavor hadron productions. A linear Boltzmann transport model coupled to hydrodynamic simulation is utilized to study the nuclear modification of heavy and light flavor jets in the dynamically evolving quark-gluon plasma. Within our jet quenching framework, we obtain a nice description of the nuclear modification factors of charged hadrons, D mesons, B mesons and B–decayed D mesons over a wide range of transverse momentum (8-300 GeV). Our result shows that perturbative QCD is sufficient to explain the color, mass and energy dependence of parton energy loss and jet quenching in heavy-ion collisions.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Enhanced Diffusive Transport in Fluctuating Porous Media
Mass transport within porous structures is a ubiquitous process in biological, geological, and technological systems. Despite the importance of these phenomena, there is no comprehensive theory that describes the complex and diverse transport behavior within porous environments. While the porous matrix itself is generally considered a static and passive participant, many porous environments are in fact dynamic, with fluctuating walls, pores that open and close, and dynamically changing cross-links. While diffusion has been measured in fluctuating structures, notably in model biological systems, it is rarely possible to isolate the effect of fluctuations because of the absence of control experiments involving an identical static counterpart, and it is generally impossible to observe the dynamics of the structure. In this paper, we present a direct comparison of the diffusion of nanoparticles of various sizes within a trackable, fluctuating porous matrix and a geometrically equivalent static matrix, in conditions spanning a range of regimes from obstructed to highly confined. The experimental system comprised a close-packed layer of colloidal spheres that were either immobilized to a planar surface or allowed to fluctuate locally, within the space defined by their nearest neighbors. Interestingly, the effective long-time diffusion coefficient was approximately 35–65% greater in the fluctuating porous matrix than in the static one (depending on the size of the nanoparticle probes), regardless of the geometric regime. This was explained by considering the enhancing effects of matrix fluctuations on the short-time diffusion coefficient and cooperative “gate-opening” motions of matrix particles and nanoparticle probes.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗