Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamics modeling”

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 739 records · Page 41

Anisotropic flow and flow fluctuations of identified hadrons in Pb–Pb collisions at $ \sqrt{s_{\textrm{NN}}}$ = 5.02 TeV

The first measurements of elliptic flow of π ± , K ± , p+$\bar{p}$, K$_{S}^{0}$, Λ+$\bar{Λ}$, $\phi$, Ξ – +$\barΞ$ + , and Ω – +$\barΩ$ + using multiparticle cumulants in Pb–Pb collisions at $\sqrt{s_{NN}}$= 5.02 TeV are resented. Results obtained with two- (v 2 {2}) and four-particle cumulants (v 2 {4}) are shown as a function of transverse momentum, p T , for various collision centrality intervals. Combining the data for both v 2 {2} and v 2 {4} also allows us to report the first measurements of the mean elliptic flow, elliptic flow fluctuations, and relative elliptic flow fluctuations for various hadron species. These observables probe the event-by-event eccentricity fluctuations in the initial state and the contributions from the dynamic evolution of the expanding quark–gluon plasma. The characteristic features observed in previous p T -differential anisotropic flow measurements for identified hadrons with two-particle correlations, namely the mass ordering at low p T and the approximate scaling with the number of constituent quarks at intermediate p T , are similarly present in the four-particle correlations and the combinations of v 2 {2} and v 2 {4}. In addition, a particle species dependence of flow fluctuations is observed that could indicate a significant contribution from final state hadronic interactions. The comparison between experimental measurements and CoLBT model calculations, which combine the various physics processes of hydrodynamics, quark coalescence, and jet fragmentation, illustrates their importance over a wide p T range.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Validation of SPH code Spheral to model interacting solid bodies in a supersonic flow

Contemporary discussions of planetary defense involve analyzing the risks posed by smaller sized, 20 to 200 m diameter, asteroids which are capable of breaking up in the atmosphere and generating a blast wave. Consequence assessments for this size class of asteroids are performed through fast-running analytic or semi-analytic models which are informed by high-fidelity hydrocode simulations of asteroid entry and breakup. However, insufficient historical data necessitates validating the independent physical processes which dominate airburst events. Here, the Fluid Solid Interface Smoothed Particle Hydrodynamics solver was previously used by Pearl et al. in 2023 to model the Chelyabinsk airburst and is used here to perform a series of validation simulations. The first effort involves modeling a cylinder in a hypersonic flow and comparing the bow shock geometry to that predicted by analytic theory. The second effort involves modeling the separation of two spherical bodies in supersonic flow and validating against experimental footage. Combined, these exercises demonstrate the ability of the code to model the flight-path of interacting solid bodies in a hypersonic flow.

Airburst↗

Characterization of an imploding cylindrical plasma for electron transport studies using x-ray emission spectroscopy

We report here on the characterization of the conditions of an imploding cylindrical plasma by time-resolved x-ray emission spectroscopy. Knowledge about this implosion platform can be applied to studies of particle transport for inertial confinement fusion schemes or to astrophysical plasmas. A cylindrical Cl-doped CH foam within a tube of solid CH was irradiated by 36 beams (I total ~5×10 14 W/cm 2 , 1.5ns square pulse, and E total ~ 16.2kJ) of the OMEGA-60 laser to radially compress the CH toward the axis. The analysis of the time-resolved spectra showed that the compression can be described by four distinct phases, each presenting different plasma conditions. First the ablation of the cylinder is dominant; second, the foam is heated and induces a significant jump in emission intensities; third, the temperature and density of the foam reaches a maximum; and finally, the plasma expands. Ranges for the plasma temperature were inferred with the atomic physics code SCRAM (Spectroscopic Collisional-Radiative Atomic Model) and the experimental data have been compared to hydrodynamic simulations performed with the 2D code FLASH, which showed a similar implosion dynamic over time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

Assess and Benchmark Magneto-Inertial Fusion (MIF) Scaling. Final Technical Report for the SNL/Rochester ALPHA Follow-on Project

This project was a follow-on to the Sandia National Laboratories (SNL) and the Laboratory for Laser Energetics (LLE) ARPA-E ALPHA project entitled “Demonstrating Fuel Magnetization and Laser Heating Tools for Low-Cost Fusion Energy”. The primary purpose of this follow-on project was to obtain additional data at the OMEGA facility to help better understand how MagLIF, a platform that has already demonstrated the scientific viability of magneto-inertial fusion, scales across a factor of 1000 in driver energy. A secondary aspect of this project was to extend simulations and analysis at SNL to cover a wider magneto-inertial fusion (MIF) parameter space and test scaling of those models across this wide range of input energies and conditions of the target. This work was successful in improving understanding of how key physics elements of MIF scales and improves confidence in setting requirements for fusion gain with larger drivers. The OMEGA experiments at the smaller scale verified the hypothesis that preheating the fuel plays a significant role in introducing wall contaminants that mix into the fuel and significantly degrade fusion performance. This contamination not only impacts target performance but the optimal input conditions for the target. However, analysis at the Z-scale showed that target performance at high preheat levels is limited by the Nernst effect, which advects magnetic flux from the hot spot, reducing magnetic insulation and consequently reduces the temperature of the fuel. The combination of MagLIF experiments at the disparate scales of OMEGA and Z along with a multiscale 3D simulation analysis has led to new insight into the physical mechanisms responsible for limiting target performance and provides important benchmarks to assess target scaling more generally for MIF schemes. Finally, in addition to the MagLIF related work, a semi-analytic model of liner driven Field Reversed Configuration (FRC) was developed that predicts the fusion gain for such systems. This model was also validated with 2D radiation magneto-hydrodynamic simulations and predicts that fusion gains of near unity could be driven by the Z machine.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Taylor wave solution for a general equation of state

This document describes a solution procedure for calculating the Taylor wave behind an unsupported Chapman–Jouguet (CJ) detonation in planar, cylindrical, and spherical geometries given a general equation of state. The resulting semi-analytic solution can be utilized to examine new equation of state models for detonation products and during the verification of hydrodynamic codes. The governing partial differential equations are reduced to ordinary differential equations in both characteristic and self-similar forms. The first-order systems corresponding to each geometry are amenable to solution numerically using commonly available methods. A difficulty arises at the CJ point in radial coordinates where the similarity equations become singular. Two separate strategies are proposed to integrate the first-order system. The first one uses an asymptotic approximation near the CJ point that can be used to perturb the boundary conditions. The second one applies a change of variables which removes the singularity at the expense of an additional equation to be integrated. A test problem is provided for the Davis products equation of state to illustrate the qualitative features of the Taylor wave in each geometric configuration and compared with a Lagrangian hydrodynamics research code. A Python code listing gives an implementation using the SciPy library to assists users in generating the results.

97 MATHEMATICS AND COMPUTING↗

Model-based analysis of solute transport and potential carbon mineralization in the active layer of a hillslope underlain by permafrost with seasonal variability and climate change

Permafrost carbon, stored in frozen organic matter across vast Arctic and sub-Arctic regions, represents a substantial and increasingly vulnerable carbon reservoir. As global temperatures rise, the accelerated thawing of permafrost releases greenhouse gases, exacerbating climate change. However, freshly thawed permafrost carbon may also experience lateral transport by groundwater flow to surface water recipients such as rivers and lakes, increasing the terrestrial-to-aquatic transfer of permafrost carbon. Mobilization and subsurface transport mechanisms are poorly understood and not accounted for in global climate models, leading to high uncertainties in the predictions of the permafrost carbon feedback. Here, we focus on a hillslope in Endalen Valley, Svalbard, as a representative example of a high-Arctic hillslope underlain by continuous permafrost. We analyze solute transport in the form of a non-reactive tracer representing dissolved organic carbon (DOC) using a physics-based numerical model with the objective to study governing cryotic and hydrodynamic transport mechanisms relevant for warming permafrost regions. We first analyze transport times for DOC pools at different locations within the active layer under present-day climatic conditions and proceed to study susceptibility for deeper ancient carbon release in the upper permafrost due to thaw under different warming scenarios. Results suggest that DOC in the active layer near the permafrost table experiences rapid lateral transport upon thaw due to saturated conditions and lateral flow, while DOC close to the ground surface experiences slower transport due to flow in unsaturated soil. Deeper permafrost carbon release exhibits vastly different transport behaviors depending on warming and thaw rate. Gradual warming leads to small fractions of DOC being mobilized every year, while the majority moves vertically through percolation and cryosuction. Abrupt thaw resulting from a single very warm year leads to faster lateral transport times, similar to active layer DOC released in saturated conditions. Lastly, we analyze the potential susceptibility of DOC to mineralization to CO 2 prior to export due to soil moisture and temperature conditions. We find that high liquid saturation during transport coincides with very low mineralization rates and potentially inhibits mineralization into CO 2 before export. Overall, the results highlight the importance of subsurface hydrologic and thermal conditions for the retention and lateral export of permafrost carbon by subsurface flow.

Hamm, Alexandra [Stockholm Univ. (Sweden)] (ORCID:↗

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↗

Simple and Sophisticated Models of Taylor's Cylinder Impact

We present a study based on the simplest possible models for Taylor’s cylinder impact problem, in addition to examination of using convolutional neural networks (CNNs) to map cylinder profiles to strength calibrations. We find that the approximate treatments of Taylor and Hawkyard compare well with hydrodynamic simulations using an equivalent assumption of constant flow stress. The CNN models prove to be well suited to successfully infer parameterizations of the Preston-Tonks-Wallace model of plastic deformation based on the deformed profile of an impacted cylinder

36 MATERIALS SCIENCE↗

Quantitative interpretation of time-lapse seismic data at Farnsworth field unit: Rock physics modeling, and calibration of simulated time-lapse velocity responses

Here, this study investigates the contribution of fluid saturation variation to the time-lapse velocity response by performing fluid substitution modeling. The methodology is exemplified by the time-lapse seismic monitoring of carbon dioxide at Farnsworth field unit (FWU). In order to evaluate the fluid distribution in a matured oil reservoir, the Southwest Regional Partnership (SWP) acquired multiple vertical seismic profile (VSP) surveys at different times during the CO 2 –water alternatinggas (WAG) injection period. In this work, we present a thorough methodology for computing the elastic response of the saturated rock for different fluid saturations using a site-specific petro-elastic model (PEM). The output from the PEM was combined with results from a fluid compositional model to compute the seismic velocities at times corresponding to each VSP survey. To produce a calibrated simulated response, the measured time-lapse seismic velocities were integrated into the numerical simulation model. The mismatches between the predicted and measured time-lapse velocities were minimized through an iterative calibration process using a trained artificial neural network proxy (ANN) coupled with a particle swarm optimizer (PSO). Our study indicates that the hybrid optimization workflow can effectively perform the history matching. With an accurate prediction of the hydrodynamic properties, the migration of CO 2 within the subsurface was modeled by predicting the spatial velocity distribution for a radius of 305 m around the injection well. The technology demonstrated and the expertise gained from this study can guide similar CO 2 -WAG projects.

58 GEOSCIENCES↗

A New Method for Ice-Ice Aggregation in the Adaptive Habit Model

In this paper, a novel methodology for modeling ice-ice aggregation is presented. This methodology combines a modified hydrodynamic collection algorithm with bulk aggregate characteristic information from an offline simulator that collects ice particles, namely, the Ice Particle and Aggregate Simulator, and has been implemented into the Adaptive Habit Microphysics scheme in the Weather Research and Forecasting Model. Aggregates, or snow, are formed via collection of cloud ice particles, where initial ice characteristics and the resulting geometry determine aggregate characteristics. Upon implementation, idealized squall-line simulations are performed to examine the new methodology in comparison with commonly used bulk microphysics schemes. It is found that the adaptive habit aggregation parameterization develops snow and reduces ice mass and number concentrations compared to other schemes. The development of aggregates through the new methodology cascades into other interesting effects, including enhancements in ice and snow growth, as well as homogeneous freezing. Further microphysical analyses reveal varying sensitivities, where snow processes are most sensitive to the new parameterization, followed by ice, then cloud, rain, and graupel processes. Further, the new scheme results in enhancements in surface precipitation due to the persistence of snow at lower altitudes. This persistence is a result of shape-dependent melting and sublimation, increasing the residence time. Moreover, these low-level enhancements are reflected in increases in radar reflectivity at the surface and its spatial distribution. Finally, the ability to predict snow shape and density allows for the simulation of polarimetric radar quantities, resulting in signature enhancements compared to schemes that do not consider spatial and temporal variations in snow shape and density.

54 ENVIRONMENTAL SCIENCES↗

Design of experiments to spectroscopically characterize radiation flow in stochastic media

Precise characterization of experimental radiation flow is required to validate the high energy density physics models, numerical methods, and codes that are used to simulate radiation-hydrodynamics phenomena such as thermal radiation transport in stochastic media. The Cassio code is used to simulate thermal radiation flow through inhomogeneous, stochastic-media-foam configurations containing optically thick clumps dispersed within an optically thin background aerogel. Cassio can model small inhomogeneous problems directly, but most problems require approximations to meet computer limitations on run-times and memory usage. Various examples of these approximations are methods that produce, in one calculation, an ensemble-averaged solution and associated standard deviation; reduced spatial dimensionality with approximate geometries; and full material homogenization with no geometric detail. Cassio simulations are used to design experiments at the OMEGA-60 Laser Facility that can measure the radiation flow using the spatially resolved COAX absorption spectroscopy diagnostic. The experimental platforms flow radiation through foam targets ranging from a background-only aerogel, to a single configuration of a specified stochastic medium, to a fully homogenized foam of the background and clump materials. Under constant total clump mass, larger clumps (here, larger than 10 μm diameter) will mix more slowly with the background such that the bulk radiation flow is faster than it would be in a fully homogenized material. The COAX platform can be used to infer temperature and density profiles in both the background material and clumps, simultaneously, and therefore to differentiate radiation flow in a range of stochastic and homogeneous media.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Reduced-Order CFD Modeling to Support Waste Loading Optimization in Hanford WTP Vitrification

The U.S. DOE Hanford Site stores over 56 million gallons of radioactive liquid tank waste that must be treated and immobilized for long-term disposal The Waste Treatment and Immobilization Plant (WTP) will vitrify this waste by feeding it into Joule-heated melters, where it is incorporated into a stable borosilicate glass Computational fluid dynamics (CFD) simulations of glass melters can provide insight into the maximum achievable waste loading under varying melter operating conditions Fully resolved VOF multiphase simulations were used as the reference model to capture bubble-driven convection in the melter, including bubble formation, rise behavior, and induced glass melt circulation Effective bubble column diameter and rise velocity were extracted from the resolved simulations, compared with empirical correlations, and refit across relevant viscosity and gas flow rate conditions Explicit gas–liquid interface tracking was replaced with a single-phase momentum source term model, enabling faster steady-state CFD simulations while preserving the dominant hydrodynamic effects of bubbling New empirical correlations were developed for effective bubble column diameter and bubble rise velocity by fitting resolved simulation data across expected melter viscosity and gas flow rate ranges, providing improved inputs for the momentum source term model compared with existing literature correlations The momentum source term model reduced fluid-domain mesh size by 89% and achieved an 8.4× computational speedup relative to resolved bubbling simulations The validated momentum source term approach enables prediction of process-relevant heat transfer behavior in the integrated melter model, including heat transfer from the molten glass to the cold cap, plenum, refractory walls, and surrounding structural regions under varying melter operating conditions

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗

Validation of ablation model for polyethylene using pulsed x-ray and proton exposures

The surface erosion of polyethylene is interrogated using pulsed x rays at the Z Machine (Sandia National Laboratories) and with proton beams at the Gamble II generator (Naval Research Laboratory) to validate a coupled model for volumetric thermal ablation, photoionization, finite-rate decomposition, and molecular recombination of radicals. The intense radiation pulses (up to ∼1014W/m2 over tens of nanoseconds) are used to generate one-dimensional vapor flows with low ionization fractions and a simplified geometry compared to typical laser ablation, allowing for evaluation of the model under local thermal equilibrium conditions. Areal momentum carried by the ensuing uniaxial hydrodynamic shock is used to indicate the extent of ablation. The threshold fluence for ablation is found to be in close correspondence with the bulk melt transition, and reasonable agreement with the model is obtained for peak temperatures in polyethylene up to 5500 K and heating rates up to 1011K/s where thermal decomposition reactions are also active.

Moore, N. W. (ORCID:0000000272259674)↗

Uncertainty quantification of bank vegetation impacts on the flood flow field in the American River, California, using large‐eddy simulations

Bank vegetation plays a key role in both hydrodynamics and morphodynamics of natural rivers; however, these effects are often unaccounted for in the computational flow dynamics of natural waterways. Recent studies using the large‐eddy simulation (LES), however, have attempted to gain insights into the impacts of bank vegetation on the mean flow field of the natural rivers using a vegetation model, which applies a sink term to the momentum equations of motion. This approach accounts for the effects of the vegetation and provides a practical approach to account for the complex patches of bank vegetation in large‐scale rivers. To implement the vegetation model, a drag coefficient reflecting the overall resistance of vegetal structures to the flow is needed, but due to the lack of calibrated data and range of size, density and type of vegetation, this parameter can be a significant source of uncertainty in the model results. Here, in this study, we use uncertainty quantification (UQ) to investigate the hydrodynamics and bed shear results when a bank vegetation is incorporated in an LES model. To this end, we used the polynomial chaos expansion and Monte Carlo sampling techniques to determine the uncertainties associated with the drag coefficient in the vegetation model and from uncertainties in the bed roughness and inflow discharge. The UQ analysis provided spatially varying confidence levels for the spanwise and vertical distribution of velocity magnitude and for the bed shear stress distributions. In addition, Sobol indices were computed to indicate the relative influence that each parameter had on the overall uncertainty. In general, it was found that uncertainty in flow discharge was the dominant source of uncertainty; however, the drag coefficient in the vegetation model and the bed roughness parameter also made significant contribution to the uncertainty near the banks and bed, respectively.

54 ENVIRONMENTAL SCIENCES↗

Analytical model for the motion and interaction of two-dimensional active nematic defects

Here, we develop an approximate, analytical model for the velocity of defects in active nematics by combining recent results for the velocity of topological defects in nematic liquid crystals with the flow field generated from individual defects in active nematics. Importantly, our model takes into account the long-range interactions between defects that result from the flows they produce as well as the orientational coupling between defects inherent in nematics. Our work complements previous studies of active nematic defect motion by introducing a linear approximation that allows us to treat defect interactions as two-body interactions and incorporates the hydrodynamic screening length as a tuning parameter. We show that the model can analytically predict bound states between two +1/2 winding number defects, effective attraction between two –1/2 defects, and the scaling of a critical unbinding length between ±1/2 defects with activity. The model also gives predictions for the trajectories of defects, such as the scattering of +1/2 defects by –1/2 defects at a critical impact parameter that depends on activity. In the presence of circular confinement, the model predicts a braiding motion for three +1/2 defects that was recently seen in experiments, as well as stable and ergodic trajectories for four or more defects.

36 MATERIALS SCIENCE↗

Dynamics of a single polyampholyte chain

Polymers that feature both positive and negative charges along chains, known as polyampholytes, represent a class of materials that hold promise for a new generation of energy storage devices, the design of which will require knowledge of the underlying structure and dynamics. Here, we develop a theory based on the Rouse model for the dynamic structure factor of a single polyampholyte chain in the weak coupling regime (negligible intramolecular electrostatics) or subjected to weak external electric fields (governed by linear response). Neglecting effects of small ions, we find deviations in scaling from the classic Rouse theory and make predictions for scattering experiments performed on polyampholytes. We find that, under weak coupling with arbitrarily strong fields, the dynamics are highly dependent on the charge distribution and consequently look at two representative examples—random charge densities and periodic charge densities—with different scaling properties. Additionally, under weak fields, the dynamics are largely independent of charge distribution. Finally, we investigate the influence of hydrodynamic effects and the implications of including inertial effects in the model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Solid-to-plasma transition of polystyrene induced by a nanosecond laser pulse within the context of inertial confinement fusion

Laser Direct Drive (LDD) Inertial Confinement Fusion (ICF) involves irradiating a spherical target of thermonuclear fuel coated with an ablator, usually made of polystyrene. Laser energy absorption near the target surface leads to matter ablation, hydrodynamic shocks, and ultimately capsule implosion. The conservation of spherical symmetry is crucial for implosion efficiency, yet spatial modulations in laser intensity can induce nonuniformities, causing the laser imprint phenomenon. Understanding laser imprint, especially considering the initial solid state, is essential for advancing LDD ICF. A first microscopic model of solid-to-plasma transition was built in 2019, accounting for laser absorption in the solid state with a band structure based ionization model. This model has been improved to include chemical fragmentation and a more accurate description of electron collision frequency in various matter states. The latest development involves assessing the model reliability by comparing theoretical predictions with experimental observations. Despite the success of this approach, questions remain, leading to further investigations and observations under different irradiation conditions. This work presents an experiment with a nanosecond pulse, taking into account hydrodynamic effects, and measures transmission dynamics over the entire laser beam area to observe two-dimensional effects. As a result, the objective is to adapt the theoretical model, couple it with a hydrodynamic code, and observe additional effects related to the initial solid state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗