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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

Migration of Accreting Planets in Radiative Discs From Dynamical Torques

We present the results of hydrodynamical simulations of the orbital evolution of planets undergoing runaway gas accretion in radiative discs. We consider accreting disc models with constant mass flux through the disc, and where radiative cooling balances the effect of viscous heating and stellar irradiation. We assume that 20–30 M(sub ⨁) giant planet cores are formed in the region where viscous heating dominates and migrate outward under the action of a strong entropy-related corotation torque. In the case where gas accretion is neglected and for an α viscous stress parameter α = 2 × 10(exp -3), we find evidence for strong dynamical torques in accreting discs with accretion rates M ̇ ≳ 7 × 10(exsp -8)M(sub ⨀) yr(exp -1). Their main effect is to increase outward migration rates by a factor of ∼2 typically. In the presence of gas accretion, however, runaway outward migration is observed with the planet passing through the zero-torque radius and the transition between the viscous heating and stellar heating dominated regimes. The ability for an accreting planet to enter a fast migration regime is found to depend strongly on the planet growth rate but can occur for values of the mass flux through the disc of M ̇ ≳ 5 × 10(exp -8)M(sub ⨀) yr(exp -1).We find that an episode of runaway outward migration can cause an accreting planet formed in the 5–10 au region to temporarily orbit at star–planet separations as large as∼60–70 au. However, increase in the amplitude of the Lindblad torque associated with planet growth plus change in the streamline topology near the planet systematically cause the direction of migration to be reversed. Subsequent evolution corresponds to the planet migrating inward rapidly until it becomes massive enough to open a gap in the disc and migrate in the type II regime. Our results indicate that a planet can reach large orbital distances under the combined effect of dynamical torques and gas accretion, but an alternative mechanism is required to explain the presence of massive planets on wide orbits.

Accretion↗
DOE OSTI2022

Heavy Elements and Electromagnetic Transients from Neutron Star Mergers

Compact binary mergers involving neutron stars can eject a fraction of their mass to space. Being extremely neutron rich, this material undergoes rapid neutron capture nucleosynthesis, and the resulting radioactivity powers fast, short-lived electromagnetic transients known as kilonova or macronova. Such transients are exciting probes of the most extreme physical conditions and their observation signals the enrichment of the Universe with heavy elements. Here the current understanding of the mass ejection mechanisms, the properties of the ejecta, and the resulting radioactive transients are reviewed. The first well-observed event in the aftermath of GW170817 delivered a wealth of insights, but much of today's picture of such events is still based on a patchwork of theoretical studies. Apart from summarizing the current understanding, questions where no consensus has been reached yet are also pointed out, and possible directions for the future research are sketched. In an appendix, a publicly available heating rate library based on the WinNet nuclear reaction network is described, and a simple fit formula to alleviate the implementation in hydrodynamic simulations is provided.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Revisiting the Kinematics of the Cylinder Test

The cylinder expansion experiment is a well-established performance test for condensed explosives that is utilized routinely to determine pressure-energy-volume relationships for detonation products. The modern cylinder test employs optical interferometric techniques to measure velocity, as opposed to older realizations where streak cameras were commonplace. Despite their widespread use, questions sometimes remain as to what kind of data the velocity diagnostics in a cylinder test provide along with how to interpret such measurements and relate them to physical quantities of interest. Here in this study, equations are derived that fully describe the kinematics of the cylinder wall during expansion. These equations are then applied to experimental as well as synthetic data generated via a hydrodynamic simulation in order to verify the mathematical framework developed.

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

Sediment transport and morphodynamic changes due to upstream variation and sea level rise on the Western Mekong Delta, Vietnam

The hydrological and morphological dynamics across the Vietnam Mekong Delta (VMD) coastline are increasingly shaped by human interventions and climate change, affecting the sustainability of socio-economic environments. This study employs the TELEMAC-MASCARET modeling system (TELEMAC-2D, TOMAWAC, and SISYPHE modules) to simulate hydrodynamics, sediment transport, and morphological changes for a 10-year period under baseline (SC1) and future climate (SC2) scenarios, considering wave climatology of Northeast (NEM) and Southwest monsoons (SWM). The results of sedimentation and erosion were analyzed across five key coastal zones from Ca Mau Cape (KV1) to Ha Tien (KV5). Model results highlight distinct seasonal patterns of sediment redistribution shaped by monsoonal flows, wave action, and seabed morphology. During the NEM, sediment deposition is concentrated near estuarine zones, while the SWM favors offshore sediment accumulation due to perpendicular wave actions, and shallow, gently sloping bathymetry that dissipates wave energy. Coastal erosion predominantly occurs at convex shoreline features subjected to stronger coastal currents, particularly in fine-grained sediment areas. Under future sea level rise and land subsidence scenarios, reduced flow velocities lead to lower bottom shear stress, promoting deposition and mitigating erosion problems despite increased wave energy. Besides that, sediment transport analyses across 14 cross-sections (C-S) reveal a southward trend during the NEM and a northward reversal during the SWM, with upstream sediment reductions primarily influencing near-estuarine zones and exerting minimal impact beyond the mid-coastal regions. Furthermore, these findings advance the understanding of sediment redistribution processes under climate change and seasonal change, and provide a scientific basis for adaptive coastal management and sustainable development.

Climate change↗
DOE OSTI2023

Cladding and Butt-Joining Dissimilar Aluminum Alloys Simultaneously via Shear Assisted Processing and Extrusion

Shear assisted processing and extrusion has been used to fabricate multi-metallic tubing that incorporates regions with clad and butt-joined structures. In this work, aluminum alloys 6061, 7075, and 1100 are combined to create a monolithic extruded tubing with weldable ends (6061), a high strength center (7075), and a corrosion-resistant cladding (1100). Optical microscopy shows a gradual, defect-free, transition between alloys at the butt-joint interfaces and a sharp transition between alloys at the clad interface. Tensile testing proved the mechanical coherency of the joint, with failure occurring outside the tri-metallic structure. In conclusion, smoothed particle hydrodynamic simulations identified a spiral pattern of the cladding alloy that leads to inconsistent cladding thickness, and identified a simplified 2-alloy combination of aluminum alloys 6082 and 2014 that accomplished the design objectives of the tri-metallic structure.

Coatings↗
DOE OSTI2025

3D mesh regularization within an ALE code using a weighted line sweeping method

The Lagrangian formalism is widely used to simulate hydrodynamic responses in complex engineering applications, particularly those involving strong shock waves. However, as the mesh moves with the fluid, it can become highly distorted, requiring a regularization step. This involves constructing a new grid and remapping conservative quantities onto it to restore mesh quality. This work introduces a regularization method for block-structured meshes within a 3D ALE (Arbitrary Lagrangian-Eulerian) code. The proposed approach prevents mesh tangling while preserving the anisotropic features of the initial Lagrangian mesh. This regularization technique incorporates aspect ratio-based weights to control mesh smoothing. Unlike uniform rezoning techniques, this weighted approach maintains proximity to the Lagrangian mesh while improving mesh quality. Here, the method effectively handles concave geometries by mitigating the grid attraction phenomenon, which typically leads to mesh concentration along concave edges. Numerical experiments demonstrate its efficiency in regularizing severely deformed meshes, and its integration within the ALE framework is validated on challenging hydrodynamic test cases, including the triple point problem.

42 ENGINEERING↗
DOE OSTI2021

Integrated modeling analysis of estuarine responses to extreme hydrological events and sea-level rise

Estuaries are vulnerable to climate change and are particularly sensitive to extreme river discharge (RD) and sea-level rise (SLR). To understand the effects of extreme RD and SLR on estuaries, an integrated modeling approach is adopted to simulate hydrodynamic changes in the Yangtze River Estuary. Model results showed that SLR will intensify the estuarine hydrodynamics in water level, currents, and salinity. Water level, current speed, and salinity increased linearly with the magnitude of SLR. The hydrological model suggested a lower RD in the future. This will increase the saltwater intrusion, which will be more severe during the lower RD and drought extreme hydrological event (EHE) conditions. Model results also suggested that the saltwater intrusion is more sensitive to RD than SLR. In addition, the interplay between the RD and SLR will exacerbate the impacts of coastal oceans on estuarine hydrodynamics, especially during extremely low RD conditions. Moreover, the decrease of RD and deepening of the channel due to human activities will also amplify the effects of SLR on the Estuary, increasing the risk of coastal disasters.

3-D hydrodynamic model, extreme climate, draught, ↗
DOE OSTI2020

Efficient and robust stability analysis in the internal energy, volume, and moles (UVN) space

Hydrodynamic simulation of fluids in vessels is often based on the internal energy (U), volume (V), and the number of moles of all components ( re N = N 1 , N 2 , ..., N n ) as the working variables to determine the thermodynamic state by finding the global maximum of entropy. The procedure of finding the equilibrium state may be divided into two steps: stability analysis and phase-split computations. Stability analysis is performed first, and the phase-split computation proceeds if the system is unstable. Most of the past work on UVN space has been based on phase-split computations, while the stability analysis is barely investigated. Furthermore, previous studies for the stability analysis use only Newton's method to solve the non-linear algebraic equations. We present an efficient and robust approach for stability analysis where UVN is specified. The successive substitution iteration (SSI) is used to provide good initial guesses for Newton's method. The proposed approach results in a reduced number of unknowns and does not require a large number of iterations to achieve convergence in Newton's method. Overall, our proposed formulation is compatible with different equations of state and is applicable to both pure component and multi-components. The robustness and efficiency of the algorithm for stability analysis are demonstrated in various examples.

02 PETROLEUM↗
DOE OSTI2020

Review of the first charged-particle transport coefficient comparison workshop

We present the results of the first Charged-Particle Transport Coefficient Code Comparison Workshop, which was held in Albuquerque, NM October 4–6, 2016. In this first workshop, scientists from eight institutions and four countries gathered to compare calculations of transport coefficients including thermal and electrical conduction, electron–ion coupling, inter-ion diffusion, ion viscosity, and charged particle stopping powers. In this paper, we give general background on Coulomb coupling and computational expense, review where some transport coefficients appear in hydrodynamic equations, and present the submitted data. Large variations are found when either the relevant Coulomb coupling parameter is large or computational expense causes difficulties. Understanding the general accuracy and uncertainty associated with such transport coefficients is important for quantifying errors in hydrodynamic simulations of inertial confinement fusion and high-energy density experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Conservative high-order data transfer method on generalized polygonal meshes

A conservative data transfer (remap) between two meshes is an important step of arbitrary Lagrangian-Eulerian (ALE) hydrodynamics simulations. High-order numerical methods for ALE simulations require both high-order (curvilinear) meshes and high-order remap algorithms. Here we develop a conservative and bounds-preserving method for accurate remapping of discrete fields on generalized polygonal meshes with curvilinear edges. The properties of the proposed method are studied theoretically and numerically for various (smooth and non-smooth) mesh deformations and discrete fields that represent smooth and discontinuous functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Search for CME in U+U and Au+Au collisions in STAR with different approaches of handling backgrounds

We report the chiral magnetic effect (CME) refers to charge separation along a strong magnetic field between left- and right-handed quarks, caused by interactions with topological gluon fields from QCD vacuum fluctuations. We present two approaches to handle the dominant elliptic flow (v 2 ) background in the three-particle correlator (Δγ 112 ), sensitive to CME. In the first approach, we present the Δγ 112 and Δγ 123 measurements in U+U and Au+Au collisions. While hydrodynamic simulations including resonance decays and local charge conservation predict that Δγ 112 scaled by N part /ν 2 will be similar in U+U and Au+Au collisions, the projected B-field exhibits a distinct difference between the two systems and with varying N part . Therefore, U+U and Au+Au collisions provide configurations with different expectations for both CME signal and background. Moreover, the three-particle observable Δγ 123 scaled by N part /ν 3 provide baseline measurement for only the background. In the second approach, we handle the ν 2 background by measuring Δγ 112 with respect to the planes of spectators measured by Zero Degree Calorimeters and participants measured by Time Projection Chamber. These measurements contain different amounts of contributions from CME signal (along B-field, due to spectators) and ν 2 background (determined by the participant geometry). With the two Δγ 112 measurements, the possible CME signal and the background contribution can be determined. We report such a measurement in Au+Au collisions at $\sqrt{s_{NN}}$ = 27GeV with the newly installed event plane detector, and report the new findings in U+U system where the spectator-participant plane correlations are expected to differ from those in Au+Au collisions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

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↗
DOE OSTI2020

Evaluating Exposures of Bay Mussels (Mytilus trossulus) to Contaminants of Emerging Concern through Environmental Sampling and Hydrodynamic Modeling.

Bay mussels (Mytilus trossulus) were transplanted to 18 locations representing a range of potential exposures throughout Puget Sound, WA. Tissues were analyzed for over 200 organic contaminants. Results indicated the widespread exposure of marine organisms to trace levels of organic contaminants including the synthetic opioid oxycodone, present at three urban sites, and the chemotherapy drug melphalan, present at nine locations, at levels that may be of biological concern. Land-use and wastewater outfalls were evaluated as potential sources of CECs to the nearshore. Exposure to alkylphenol ethoxylates was associated with increased impervious surfaces in upland watersheds. A hydrodynamic simulation was performed using the Salish Sea Model to integrate inputs from 99 wastewater sources to Puget Sound. Predictions were consistent with concentrations of several wastewater-associated contaminants and d15N enrichment. These results support the notion that Puget Sound nearshore biota suffer chronic exposures to a suite of contaminants from multiple sources and provide critical to focus future monitoring and management.

toxic contamination, Salish Sea, Sish Sea Model, m↗
DOE OSTI2020

Suspect and Nontarget Screening for Contaminants of Emerging Concern in an Urban Estuary

This study used suspect and nontarget screening with high-resolution mass spectrometry to characterize the occurrence of contaminants of emerging concern (CECs) in the nearshore marine environment of Puget Sound (WA). In total, 87 non-polymeric CECs were identified; those confirmed with reference standards (45) included pharmaceuticals, herbicides, vehicle-related compounds, plasticizers, and flame retardants. Eight polyfluoroalkyl substances were detected; perfluorooctanesulfonic acid (PFOS) concentrations were as high as 72-140 ng/L at one location. Low levels of methamphetamine were detected in 41% of the samples. Transformation products of pesticides were tentatively identified, including two novel transformation products of tebuthiuron. While a hydrodynamic simulation, analytical results, and dilution calculations demonstrated the prevalence of wastewater effluent to nearshore marine environments, the identity and abundance of selected CECs revealed the additional contributions from stormwater and localized urban and industrial sources. For the confirmed CECs, risk quotients were calculated based on concentrations and predicted toxicities, and eight CECs had risk quotients >1. Dilution in the marine estuarine environment lowered the risks of most wastewater-derived CECs, but dilution alone is insufficient to mitigate risks of localized inputs. These findings highlighted the necessity of suspect and nontarget screening and revealed the importance of localized contamination sources in urban marine environments.

Tian, Zhenyu↗
DOE OSTI2021

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

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

42 ENGINEERING↗
DOE OSTI2024

The impact of low-mode symmetry on inertial fusion energy output in the burning plasma state

Indirect Drive Inertial Confinement Fusion Experiments on the National Ignition Facility (NIF) have achieved a burning plasma state with neutron yields exceeding 170 kJ, roughly 3 times the prior record and a necessary stage for igniting plasmas. The results are achieved despite multiple sources of degradations that lead to high variability in performance. Results shown here, for the first time, include an empirical correction factor for mode-2 asymmetry in the burning plasma regime in addition to previously determined corrections for radiative mix and mode-1. Analysis shows that including these three corrections alone accounts for the measured fusion performance variability in the two highest performing experimental campaigns on the NIF to within error. Here we quantify the performance sensitivity to mode-2 symmetry in the burning plasma regime and apply the results, in the form of an empirical correction to a 1D performance model. Furthermore, we find the sensitivity to mode-2 determined through a series of integrated 2D radiation hydrodynamic simulations to be consistent with the experimentally determined sensitivity only when including alpha-heating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Stability of the fcc phase in shocked nickel up to 332 GPa

Despite making up 5-20 wt.% of Earth’s predominantly iron core, the melting properties of elemental nickel at core conditions remain poorly understood, due largely to a dearth of experimental data. We present here an in situ X-ray diffraction study performed on laser shock-compressed samples of bulk nickel, reaching pressures up to ~ 500 GPa. Hugoniot states of nickel were targeted using a flat-top laser drive, with in situ X-ray diffraction data collected using the Linac Coherent Light Source. Rietveld methods were used to determine the densities of the shocked states from the measured diffraction data, while peak pressures were determined using a combination of measured particle velocities, shock transit times, hydrodynamic simulations, and laser intensity calibrations. We observed solid compressed face-centered cubic (fcc) Ni up to at least 332 ± 30 GPa along the Hugoniot—significantly higher than expected from the majority of melt lines that have been proposed for nickel. We also bracket the partial melting onset to between 377 ± 38 GPa and 486 ± 35 GPa.

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

Design of inertial fusion implosions reaching the burning plasma regime

Abstract In a burning plasma state 1–7 , alpha particles from deuterium–tritium fusion reactions redeposit their energy and are the dominant source of heating. This state has recently been achieved at the US National Ignition Facility 8 using indirect-drive inertial-confinement fusion. Our experiments use a laser-generated radiation-filled cavity (a hohlraum) to spherically implode capsules containing deuterium and tritium fuel in a central hot spot where the fusion reactions occur. We have developed more efficient hohlraums to implode larger fusion targets compared with previous experiments 9,10 . This delivered more energy to the hot spot, whereas other parameters were optimized to maintain the high pressures required for inertial-confinement fusion. We also report improvements in implosion symmetry control by moving energy between the laser beams 11–16 and designing advanced hohlraum geometry 17 that allows for these larger implosions to be driven at the present laser energy and power capability of the National Ignition Facility. These design changes resulted in fusion powers of 1.5 petawatts, greater than the input power of the laser, and 170 kJ of fusion energy 18,19 . Radiation hydrodynamics simulations 20,21 show energy deposition by alpha particles as the dominant term in the hot-spot energy balance, indicative of a burning plasma state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗