Conservation of B, S, and Q charges in relativistic viscous hydrodynamics solved with smoothed particle hydrodynamics
Not Available
Plumberg, Christopher (ORCID:0000000166783966)↗
SEARCH · Engineering Papers
Results for “Smoothed particle hydrodynamics”
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 New Approach to Constructing Area-Weighted Axisymmetric Discretizations for Smoothed Particle Hydrodynamics
Explore the source record for details and available documents.
Astronomy and AstroPhysics↗
Nuclear Mitigation of Hypothetical Asteroid Threats in Smoothed Particle Hydrodynamics
Explore the source record for details and available documents.
Astronomy and AstroPhysics↗
Numerical simulation of asteroid geometry variance on airburst threat
For an atmospheric airburst the primary source of concern when assessing uncertainty is the size and velocity. Determining these properties provides the basis for threat assessment, as the total energy of the asteroid may then be estimated, and the threat investigated thoroughly. Even with clarity as to how much energy an asteroid may deposit, a great deal of uncertainty still exists for the actual energy deposition process. One such source of uncertainty is the geometry of the incoming asteroid. The geometry of an asteroid will alter the stress distribution during entry, which adds uncertainty to when fracture will occur. Here, in this study, we use Smoothed Particle Hydrodynamics to model the atmospheric airburst of Tunguska-scale asteroids with varying geometric profiles, including a sphere, ellipsoid, binary and superellipsoid. Each asteroid is modeled as a homogenous structure with strength. We assess uncertainty through a series of planar 2D simulation cases for each geometry, comparing the source of stochasticity across geometries. A single 3D airburst simulation for each geometry is also analyzed. Additionally, the 3D cases are compared to the highly uncertain Tunguska event, predicting variance in burst height across geometries, but all bounded by theoretical burst heights proposed for Tunguska.
Airburst↗
Meshfree simulation and prediction of recrystallized grain size in friction stir processed 316L stainless steel
Friction stir processing (FSP) is a promising solid-phase microstructural modification technique that can repair and enhance damaged stainless steel surfaces exposed to harsh environments. The quality of the repaired material is closely correlated to the recrystallized grain size in the stir zone (SZ), which is influenced by the thermomechanical conditions dictated by FSP process parameters. Thus, establishing a reliable relationship between these parameters and recrystallized grain size in the SZ is crucial for optimizing repair quality. However, existing experimental approaches often rely on indirect temperatures measured far from the SZ, along with rough strain rate estimations, which are imprecise and time-consuming. Meanwhile, existing mesh-based modeling methods usually face numerical challenges when dealing with the large material deformations inherent in FSP. Here, to address these issues, this study introduces a meshfree process model for FSP based on the smoothed particle hydrodynamics (SPH) method, aimed at predicting process conditions under different parameters. The model is validated using experimental data from 11 combinations of tool traverse and rotation speeds on 316 L stainless steel. Correlations between process parameters, material flow, temperature, strain, strain rate, and recrystallized grain size are revealed through SPH simulations and electron backscatter diffraction (EBSD) imaging. The results show that in situ SZ temperatures range from 1071 to 1322°C, which exceed the tool temperature by over 300°C. Furthermore, SZ temperature, strain rate, and grain size increase monotonically with higher tool temperature and faster traverse speed. A relationship is then established between the model-predicted Zener-Hollomon parameter and the recrystallized grain size based on EBSD data, expressed as ln(d) = -0.364 ln(Z) + 14.673. Finally, this relationship exhibits satisfactory accuracy with errors of less than 26.9% in predicting grain sizes at various SZ locations, which offers valuable insights for optimizing FSP repair processes for 316 L stainless steel.
316L stainless steel↗
Modeling of melting and erosion in high velocity microparticle impacts
High velocity collisions can cause melting, a situation involving elasticity, plasticity, and phase transformations in tandem, which is very challenging for mechanical modeling. Here, this paper advances the use of smoothed particle hydrodynamics (SPH) simulations for this complex situation, investigating stainless steel microspheres impacting a tin substrate. After extensive validation against a range of mechanical and geometrical measurements of experimental impacts, the model provides a detailed, quantitative view of melt volume production over a broad range of impact velocities. It also enables direct decomposition of the impact energy partitioning, including analysis of the energy absorbed by plasticity in the impactor, the system kinetic energy, and the adiabatic heating of the substrate in addition to the melting itself. Since these quantities are not available experimentally, SPH modeling substantially enhances the quantitative understanding of such impacts. Finally, the model speaks to the phase characteristics of the ejecta formed during impact, confirming the experimental inference that the ejected material likely contains a significant fraction of solid material.
Constitutive modeling↗
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↗
X-ray observation of individual Ti-6Al-4V spherical powder particle impact in an in-situ operando laser directed energy deposition system
Understanding the particle capture mechanism in the laser directed energy deposition additive manufacturing process provides a foundation for improving productivity and reducing material loss. Capturing this phenomenon is challenging due to its highly transient time interval. To analyze the impact behavior of spherical Ti-6Al-4V powder, a custom laser directed energy deposition system was designed to control the deposition of individual particles. Synchrotron X-ray imaging at 24 kHz monitored the impact dynamics of individual powder particles. A multi-physics model based on the Smoothed-Particle Hydrodynamics scheme supported the analysis by providing the temperature and velocity fields of the molten pool. Results revealed that in laser directed energy deposition, surface tension forces dominated the powder–molten pool interaction. Hydrophobic Ti-6Al-4V powder particles had an equilibrium contact angle of 121 degrees with their molten state. Furthermore, Ti-6Al-4V powder particles required a high impact velocity (6 to 12 m/s) to transition from oscillation to submergence behavior. Lastly, oscillating powder particles that impacted near the location of the laser beam exhibited faster melting, highlighting the contribution of the laser beam to the melting mechanism. This work lays the foundation for investigations of powder particle impact in laser directed deposition and supports the validation of numerical models of powder–molten pool interactions.
36 MATERIALS SCIENCE↗
Structural uniformity and compositional homogeneity of solid-phase alloyed rod
Solid-phase processes have emerged as an alternative to fusion-based alloying to avoid coarse microstructures, undesirable phase formation, and high energy consumption. However, achieving uniform distribution of alloying elements during friction-based processing remains challenging due to highly heterogeneous thermomechanical conditions. This work evaluates the structural uniformity and compositional homogeneity of Al–Cu–Zn alloyed rods produced by friction extrusion (FE) and establishes the role of the rotational speed to feed rate ratio (N/V) on alloying effectiveness. A systematic matrix of FE experiments was conducted at constant extrusion ratio with N/V values ranging from 3.7 to 300. Compositional uniformity was assessed along the rod length (ICP-OES), in three dimensions (X-ray computed tomography), and at the microscale (SEM–EDS), supported by a gray-level co-occurrence matrix (GLCM)–based homogeneity metric. Smoothed particle hydrodynamics (SPH) simulations were used to reveal material flow and thermomechanical fields. Results show that N/V = 100 produces a high-shear mixing zone that eliminates the unmixed core and enables near-full dissolution and dispersion of Cu and Zn. At lower N/V, a laminar flow region persists at the rod center, causing segregation and large composition gradients. The combined experimental–computational analysis provides mechanistic insight into the transition from fragmented particle dispersion to thermomechanically assisted metallurgical mixing. This study establishes processing–structure relationships for solid-phase alloying and provides guidance for achieving homogenized compositions comparable to wrought alloys via rapid, scalable FE processing.
Aluminum↗
Viscous Gubser flow with conserved charges to benchmark fluid simulations
We present semi-analytical solutions for the evolution of both the temperature and chemical potentials for viscous Gubser flow with conserved charges. Such a solution can be especially useful in testing numerical codes intended to simulate relativistic fluids with large chemical potentials. The freeze-out hypersurface profiles for constant energy density are calculated, along with the corresponding normal vectors, and presented as a new unit test for numerical codes. We also compare the influence of the equation of state on the semi-analytical solutions. We benchmark the newly developed smoothed particle hydrodynamics code ccake that includes both shear viscosity and three conserved charges. Here, the numerical solutions are in excellent agreement with the semi-analytical solution and are also able to accurately reproduce the hypersurface at freeze-out.
Hydrodynamic models↗
High-Fidelity Numerical Wave Tank Verification & Validation Study: Wave Generation Through Paddle Motion: Preprint
This paper presents a numerical benchmark study of wave propagation due to a paddle motion using different high-fidelity numerical models, which are capable of replicating the nearly actual physical wave tank testing. A full time series of the measured wave generation paddle motion which was used to generate wave propagation in the physical wave tank will be utilized in each of the models contributed by IEA OES Task 10's participants, which includes both computational fluid dynamics (CFD) and smooth hydrodynamic particle (SPH). The high-fidelity simulations of the physical wave testcase will allow for the evaluation of the initial transient effects from wave ramp-up and its evolution in the wave tank over time for two representative regular waves with varying levels of nonlinearity. A couple of interesting metrics like the predicted wave surface elevation at select wave probes, wave period, and phase-shift in time will be assessed to evaluate the relative accuracy of numerical models versus experimental data within specified time intervals. These models will serve as a guide for modelers in the wave energy community and provide a base case to allow further and more detailed numerical modeling of the fixed Kramer Sphere Cases under wave excitation force wave tank testing.
HYDRO ENERGY,TIDAL AND WAVE POWER↗
The Effects of Realistic Impact-processed Asteroid Topology on Atmospheric Airbursts
When small asteroids (∼10–200 m) enter Earth’s atmosphere they can break up before hitting the ground, releasing a considerable amount of energy in what is known as an atmospheric airburst. To date, numerous studies have examined how an asteroid’s strength, size, shape, and entry characteristics affect the breakup process; but few have assessed the effects of heterogenous internal topologies. We present a two-step high-fidelity simulation approach, which we use to model the atmospheric breakup of asteroids with realistic impact-processed internal topologies and their associated ground effects. In the first step, we model the asteroid as it is impacted to generate an “impact-processed” topology. We then map the resulting fractures (i.e., damage field) from that simulation to the initial state of an atmospheric breakup simulation. For both phases, we use the smoothed particle hydrodynamics code Spheral++. As a demonstrative example, we apply our new simulation pipeline to a Tunguska-scale impactor (50 m diameter, 20 km s −1 entry velocity, and 45° angle) in a case study comparing impact-processed and fully intact asteroids in both planar 2D and full 3D. For sufficiently strong asteroids, we find that impact processing generally causes the asteroid to break up at higher altitudes. We also observe that impact-processed asteroids break up into smaller fragments than objects with equivalent properties that have not been impact-processed.
Astronomy and AstroPhysics↗
The Role of Disruptive Impacts on Ocean Generation and Longevity in Icy Moons
Among icy moons of the outer Solar System, subsurface water oceans and large collisions both seem to be common, but the impact of the latter on the presence and persistence of the former is unclear and has rarely been investigated. Here we interface a smoothed-particle hydrodynamics model to simulate collisions with a thermal-structural evolution model to simulate the evolution of moons pre-collision, post-collision and without a collision. Overall, even such large-scale collisions affect only the ocean thickness or longevity, and the presence or absence of an ocean is affected for only a part of a moon’s history. In reaccreted moons, the ocean survives the impact and becomes much thicker inside larger moons with radius near 1,000 km, whereas an ocean that would otherwise arise inside moons with radius near 500 km is absent because the collision promotes ice–rock differentiation. Our simulations have not yielded an ocean developed post-impact—whether directly via collisional or reaccretional heating or indirectly through tidal heating due to collision-induced orbital changes—in a moon that would otherwise have remained frozen. The ocean-enhancing effect is pronounced only for late disruptive impacts onto large, 1,000-km-class targets, which are unlikely in recent solar system history.
Marc Neveu↗
Integrated process-structure–property–performance optimization of cold-sprayed zinc coatings on AZ91 magnesium alloy
High-pressure die cast (HPDC) AZ91 magnesium alloy is widely used in automotive components such as transmission housings and brackets for its excellent strength-to-weight ratio. Zinc-based cold spray coatings can be applied selectively to vulnerable areas to enhance corrosion resistance, minimize galvanic coupling with dissimilar metals, and eliminate the need for full-surface oxide coatings, making the process more efficient and targeted. A comprehensive evaluation of 16 combinations of nitrogen carrier gas temperatures and pressures led to the identification of an optimal range of process parameters, yielding Zn coatings with porosity <0.5 % by area, wear rates reduced by a factor of two compared to uncoated AZ91, and adhesion strengths up to 35 MPa. The enhanced mechanical performance of the coating is attributed to the low porosity and the formation of a metallurgical bond at the coating-substrate interface. Corrosion studies using macroscale potentiodynamic polarization (PDP) and electrochemical impedance spectroscopy (EIS) revealed a significant decrease in corrosion rate and a shift to more noble corrosion potentials (ZCP) for coated substrates. Furthermore, the Zn cold-sprayed samples exhibited significantly lower corrosion-induced evolved hydrogen content compared to the base AZ91 substrate and AZ91 coated with industrial coatings, demonstrating that the Zn layer effectively protects the substrate from the corrosive environment. Overall, cold spray Zn coatings significantly improve the mechanical and corrosion performance of AZ91 Mg alloys, addressing key material challenges and enabling their broader use in automotive applications.
Corrosion performance↗
TEAMER Technical Support of Ramboll's for Numerical Modeling of WECs to Support OES Task 10: Cooperative Research and Development (Final Report)
National Technology & Engineering Solutions of Sandia, LLC (NTESS) in collaboration with the National Renewable Energy Laboratory (NREL), and with guidance from Ramboll, will perform fluid dynamics simulations to support of The Ocean Energy Systems (OES) Energy Technology Collaboration Program Task 10 Wave Energy Converters (WEC) Modelling Verification and Validation effort. Specific numerical simulations include the fixed device wave impingement studies to benchmark the performance of simulation techniques against physical testing results.
16 TIDAL AND WAVE POWER↗
Precision Polishing of Spheres Via In-Situ Process Monitoring and Machine-Learning-Based Optimization
In inertial confinement fusion (ICF) experiments seeking output gains of unity and beyond, the quality of the ablator capsule is paramount for minimizing hydrodynamic mix that quenches the central hot spot. Defects in the form of foreign particles or missing mass on the surface and within the wall of the capsule are primary offenders. High density carbon capsules made for ICF experiments on the National Ignition Facility (NIF) are precision polished to achieve the surface smoothness in the order of a few nm as well as to minimize isolated defects in the form of pits. Given the critical role of this process, we are developing smart manufacturing techniques with goal of elevating the efficiency of this process. Our approach is to use MEMS-based sensors to capture the fine vibrational signals generated during the polishing process and combine it with synchronized visual feedback as needed. Beyond using these sensors for process monitoring, we use specific deep learning methods to analyze the data and extract correlations with both the process parameters and the final performance of the polishing run. Here, we describe the multiple fronts that we have explored in this regard and the results we have gotten so far. This approach promises to have the potential to ultimately provide real-time feedback that can be used for ensuring the progress of the run as well as a means for faster optimization.
36 MATERIALS SCIENCE↗