Search NASA⌕ Search

SEARCH · Search NASA

Results for “SPH”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Reticular chemistry for the rational design of mechanically robust mesoporous merged-net metal-organic frameworks

Access to metal-organic frameworks (MOFs) with enhanced mechanical stability is key to their successful deployment in practical applications. However, the high porosity of the material often affects mechanical stability. In this article, to achieve highly porous MOFs with enhanced mechanical stability, we explored the merged-net approach where two relatively fragile frameworks were merged into a robust MOF structure. We demonstrate the effectiveness of this approach by computationally evaluating mechanical properties of sph-MOFs with varying lengths of linkers. Prominently, we pinpoint the significance of triangular rigidity on the robustness of large-pore MOFs and, subsequently, designed and synthesized a rare earth (RE)-based RE-sph-MOF-5 by the reticulation of hexanuclear RE clusters, tritopic linkers, and unprecedentedly large planar hexatopic linkers containing 19 phenyl rings. The mechanical properties of sph-MOFs were characterized and quantified using amplitude-frequency modulation (AM-FM) bimodal atomic force microscopy (AFM) analyses. Markedly, the mesoporous RE-sph-MOF-5 expresses high mechanical stability despite its large mesoporous cavities.

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

Electronic Structure, Optical Properties, and Photoelectrochemical Activity of Sn-Doped Fe 2 O 3 Thin Films

Hematite (Fe 2 O 3 ) is a well-known oxide semiconductor suitable for photoelectrochemical (PEC) water splitting and industry gas sensing. It is widely known that Sn doping of Fe 2 O 3 can enhance the device performance, yet the underlying mechanism remains elusive. In this work, we determine the relationship between electronic structure, optical properties and PEC activity of Sn doped Fe 2 O 3 by studying highly crystalline, well-controlled thin films prepared by pulsed laser deposition (PLD). We show that Sn doping substantially increases the n-type conductivity of Fe 2 O 3 , and the conduction mechanism is better described by small-polaron hoping (SPH) model. Only 0.2% Sn doping significantly reduce the activation energy barrier for SPH conduction from at least 0.5 eV for undoped Fe 2 O 3 to 0.14 eV for doped ones. A combination of X-ray photoemission, X-ray absorption spectroscopy and DFT calculations reveals the Fermi level gradually shifts toward the conduction band minimum with Sn doping. A localized Fe 2+ like gap state is observed at the top of valence band, accounting for the SPH conduction. Interestingly, in contrast to the literature, we find that only 0.2% Sn doping in Fe 2 O 3 significantly improves the PEC activity, while more Sn decreases it. The improved PEC activity is partially attributed to an increased band bending potential which facilitates the charge separation at space charge region. The reduced activation energy barrier for SPH will facilitate the transport of photo-excited carriers for the enhanced PEC, which is of interest for further carrier dynamics study.

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

Latent variable Gaussian process models: A rank‐based analysis and an alternative approach

Abstract Gaussian process (GP) models have been extended to emulate expensive computer simulations with both qualitative/categorical and quantitative/continuous variables. Latent variable (LV) GP models, which have been recently developed to map each qualitative variable to some underlying numerical LVs, have strong physics‐based justification and have achieved promising performance. Two versions use LVs in Cartesian (LV‐Car) space and hyperspherical (LV‐sph) space, respectively. Despite their success, the effects of these different LV structures are still poorly understood. This article illuminates this issue with two contributions. First, we develop a theorem on the effect of the ranks of the qualitative factor correlation matrices of mixed‐variable GP models, from which we conclude that the LV‐sph model restricts the interactions between the input variables and thus restricts the types of response surface data with which the model can be consistent. Second, following a rank‐based perspective like in the theorem, we propose a new alternative model named LV‐mix that combines the LV‐based correlation structures from both LV‐Car and LV‐sph models to achieve better model flexibility than them. Through extensive case studies, we show that LV‐mix achieves higher average accuracy compared with the existing two.

Tao, Siyu↗
DOE OSTI2021

The efficacy of event isotropy as an event shape observable

Event isotropy $\mathcal{I}$ sph , an event shape observable that measures the distance of a final state from a spherically symmetric state, is designed for new physics signals that are far from QCD-like. Using a new technique for producing a wide variety of signals that can range from near-spherical to jetty, we compare event isotropy to other observables. We show that thrust T and the C parameter (and λ max , the largest eigenvalue of the sphericity matrix) are strongly correlated and thus redundant, to a good approximation. By contrast, event isotropy adds considerable information, often serving to break degeneracies between signals that would have almost identical T and C distributions. Signals with broad distributions in T (or λ max ) and in $\mathcal{I}$ sph separately often have much narrower distributions, and are more easily distinguished, in the ($\mathcal{I}$ sph , λ max ) plane. An intuitive, semi-analytic estimation technique clarifies why this is the case and assists with the interpretation of the distributions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Embedded smoothed particle hydrodynamics

In this paper, a hybrid method for embedding a meshfree smoothed particle hydrodynamics into mesh-based ALE hydro schemes is presented. The method relies on the embedded mesh framework presented in Puso et al. (2015, 2012) and utilizes SPH (Owen, 2010) in the solid foreground; when an appropriate damage or failure model is used with SPH, it becomes a natural way to model material fracture and cracking in high-explosive driven experiments. Overall, we provide various examples which validate the coupling between SPH and the embedded mesh method.

42 ENGINEERING↗
DOE OSTI2023

Atomic Three-Dimensional Investigations of Pd Nanocatalysts for Acetylene Semi-hydrogenation

Deciphering the three-dimensional (3D) insight into nanocatalyst surfaces at the atomic level is crucial to understanding catalytic reaction mechanisms and developing high-performance catalysts. Nevertheless, better understanding the inherent insufficiency of a long-range ordered lattice in nanocatalysts is a big challenge. Here, in this work, we report the local structure of Pd nanocatalysts, which is beneficial for demonstrating the shape-structure-adsorption relationship in acetylene hydrogenation. The 5.27 nm spherical Pd catalyst (Pd sph ) shows an ethylene selectivity of 88% at complete acetylene conversion, which is much higher than those of the Pd octahedron and Pd cube and superior to other reported monometallic Pd nanocatalysts so far. By virtue of the local structure revelation combined with the atomic pair distribution function (PDF) and reverse Monte Carlo (RMC) simulation, the atomic surface distribution of the unique compressed strain of Pd-Pd pairs in Pd sph was revealed. Density functional theory calculations verified the obvious weakening of the ethylene adsorption energy on account of the surface strain of Pd sph . It is the main factor to avoid the over-hydrogenation of acetylene. The present work, entailing shape-induced surface strain manipulation and atomic 3D insight, opens a new path to understand and optimize chemical activity and selectivity in the heterogeneous catalysis process.

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

Preservation of kinetics parameters generated by Monte Carlo calculations in two-step deterministic calculations

The generation of accurate kinetic parameters such as mean generation time Λ and effective delayed neutron fraction β eff via Monte Carlo codes is established. Employing these in downstream deterministic codes warrants another step to ensure no additional error is introduced by the low-order transport operator when computing forward and adjoint fluxes for bilinear weighting of these parameters. Another complexity stems from applying superhomogenization (SPH) equivalence in non-fundamental mode approximations, where reference and low-order calculations rely on a 3D full core model. In these cases, SPH factors can optionally be computed for only part of the geometry while preserving reaction rates and K-effective, but the impact of such approximations on kinetics parameters has not been thoroughly studied. This paper aims at studying the preservation of bilinearly-weighted quantities in the Serpent–Griffin calculation procedure. Diffusion and transport evaluations of IPEN/MB-01, Godiva, and Flattop were carried out with the Griffin reactor physics code, testing available modeling options using Serpent-generated multigroup cross sections and equivalence data. Verifying Griffin against Serpent indicates sensitivities to multigroup energy grid selection and regional application of SPH equivalence, introducing significant errors; these were demonstrated to be reduced through the use of a transport method together with a finer energy grid.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Modeling Solids in Nuclear Astrophysics with Smoothed Particle Hydrodynamics

Abstract Smoothed particle hydrodynamics (SPH) is a frequently applied tool in computational astrophysics to solve the fluid dynamics equations governing the systems under study. For some problems, for example when involving asteroids and asteroid impacts, the additional inclusion of material strength is necessary in order to accurately describe the dynamics. In compact stars, that is white dwarfs and neutron stars, solid components are also present. Neutron stars have a solid crust, which is the strongest material known in nature. However, their dynamical evolution, when modeled via SPH or other computational fluid dynamics codes, is usually described as a purely fluid dynamics problem. Here, we present the first 3D simulations of neutron star crustal toroidal oscillations including material strength with the Los Alamos National Laboratory SPH code FleCSPH. In the first half of the paper, we present the numerical implementation of solid material modeling together with standard tests. The second half is on the simulation of crustal oscillations in the fundamental toroidal mode. Here, we dedicate a large fraction of the paper to approaches that can suppress numerical noise in the solid. If not minimized, the latter can dominate the crustal motion in the simulations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Meshfree Process Modeling and Experimental Validation of Friction Riveting of Aluminum 5052 to Aluminum 6061

Friction riveting (Fric-riveting) is an innovative, fast, and energy-efficient process for spot-joining metal-metal structures. Although fric-riveting has been studied experimentally in recent years, its process modeling is rarely found in the literature primarily because of the associated large material deformation, extreme thermomechanical conditions, and complex contact conditions. In this work, a mesh-free smoothed particle hydrodynamics (SPH) framework that can well handle the abovementioned numerical challenges is used to simulate the fric-riveting of AA5052 to AA6061. Predicted material morphology, multi-point temperatures, and plunge force are thoroughly validated by experimental observations. The material severe plastic deformation zone in the vicinity of the riveting zone is further predicted by the SPH model, which indicates the material mixing and potential grain refinement zone. Based on the validated model, process parameters can be optimized which yields better performance over the baseline case.

Friction riveting (fric-riveting), Smoothed partic↗
DOE OSTI2023

Gaussian process hydrodynamics

Abstract We present a Gaussian process (GP) approach, called Gaussian process hydrodynamics (GPH) for approximating the solution to the Euler and Navier-Stokes (NS) equations. Similar to smoothed particle hydrodynamics (SPH), GPH is a Lagrangian particle-based approach that involves the tracking of a finite number of particles transported by a flow. However, these particles do not represent mollified particles of matter but carry discrete/partial information about the continuous flow. Closure is achieved by placing a divergence-free GP prior ξ on the velocity field and conditioning it on the vorticity at the particle locations. Known physics (e.g., the Richardson cascade and velocity increment power laws) is incorporated into the GP prior by using physics-informed additive kernels. This is equivalent to expressing ξ as a sum of independent GPs ξ l , which we call modes, acting at different scales (each mode ξ l self-activates to represent the formation of eddies at the corresponding scales). This approach enables a quantitative analysis of the Richardson cascade through the analysis of the activation of these modes, and enables us to analyze coarse-grain turbulence statistically rather than deterministically. Because GPH is formulated by using the vorticity equations, it does not require solving a pressure equation. By enforcing incompressibility and fluid-structure boundary conditions through the selection of a kernel, GPH requires significantly fewer particles than SPH. Because GPH has a natural probabilistic interpretation, the numerical results come with uncertainty estimates, enabling their incorporation into an uncertainty quantification (UQ) pipeline and adding/removing particles (quanta of information) in an adapted manner. The proposed approach is suitable for analysis because it inherits the complexity of state-of-the-art solvers for dense kernel matrices and results in a natural definition of turbulence as information loss. Numerical experiments support the importance of selecting physics-informed kernels and illustrate the major impact of such kernels on the accuracy and stability. Because the proposed approach uses a Bayesian interpretation, it naturally enables data assimilation and predictions and estimations by mixing simulation data and experimental data.

Mathematics↗
DOE OSTI2025

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

Transient Multiphysics Simulations with Pin Power Reconstruction in the Griffin Reactor Physics Code

This work introduces the pin power reconstruction capability available in the Griffin reactor physics code. This capability is implemented in an unstructured mesh framework, and the methods introduced are applied to the 2D SIMBA reactor core, which has assemblies and pins arranged in a hexagonal lattice. Since this reactor has a non-Cartesian geometry and also operates in the thermal spectrum, a general approach to pin power reconstruction is adopted, where SPH-based equivalence is leveraged to preserve assembly-wise reaction rates, while computing full-core form functions to preserve pin-wise fission production rates within the fuel pins of the reactor core. In a 2D microreactor benchmark problem, this pin power reconstruction approach was shown to reproduce pin powers compared to the Serpent2 Monte Carlo code for fixed temperature conditions and control drum rotation angles, yielding a core-wide RMS error level of 0.6\% and a maximum absolute pin error of 2.3\%. In addition, a tabulated library of multigroup cross sections, SPH factors, and form functions was generated to demonstrate the applicability of pin power reconstruction to a thermal feedback problem. Finally, a control drum transient was successfully simulated, showcasing the application of pin power reconstruction in a transient multiphysics feedback problem.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Modelling and analysis of salt-convection effect on oxide reduction process for uranium oxides using smoothed particle hydrodynamics

Oxygen ion transport via molten salt continuum is a key kinetic feature of electrochemical reduction process for uranium oxides in molten salt. During the electrolysis, the evolution of oxygen ion concentration field in molten salt continuum appears from porous metallic layer to open bulk salt. Despite the important role affecting process throughput, its effect coupled with controlled salt movement was not adequately investigated. To capture the reactive oxygen ion transport phenomena, a multidimensional reduction model based on smoothed particle hydrodynamics (SPH) was combined with the advection-diffusion model. A pelletized oxide fuel bed capturing interconnected salt continuum is considered as a reactant material geometry, which allows forced convective delivery of molten salt through the reactant materials with a relatively low pressure head. Here, in this study, reduction time of a pelletized-fuel bed was quantified with respect to key design parameters of the reactant material configuration. The simulation results demonstrated the advantage of salt convection which effectively removes oxygen ion products from the reactant materials and consequently accelerates the reduction process. Moreover, the reduction-time correlation of a single pelletized fuel was developed as a function of relevant dimensionless parameters and was extended to the application for a pelletized-fuel bed. Finally, several insights and implications in the use of a pelletized oxide fuel were presented. The developed SPH-based framework enabled multi-physics analysis of governing dynamics in appreciating and screening candidate reactant material configurations in a cost-effective manner.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

An innovative and alternative approach toward gear fabrication

In this paper, we present the first demonstration of single-step gear manufacturing via an innovative in-situ friction stir forging approach. This novel process is a natural extension to the friction stir processing technique and relies on shear and normal stresses to form the complex shapes. Friction stir welding pin tool with shoulder is plunged at high revolutions per minute (rpm) into the sample to heat it to the desired temperature and enable plastic flow. Then, the tool is forged at high plunge rate to form the desired structurally-sound and defect free complex shapes. To better visualize the process the experiments are complemented by a 3D thermomechanically coupled, smoothed particle hydrodynamics (SPH) model to understand the material flow pattern both radially and through thickness in the gear teeth, plastic strain distribution, and temperature profile. We have demonstrated application of I-FSF to inherently poor formable AZ31 Mg, AA 7075-T6, and AA5083–10 vol% TiB 2 composite and fabricated 14-tooth spur gears from these materials with minimal farther machining. In conclusion, microstructure of the forged gears exhibited a directional material flow pattern due to the normal forging action and refinement of grains to 6–13 μm through dynamically recrystallization across the radial, through-thickness direction, which is in good agreement with the corresponding SPH simulation results.

42 ENGINEERING↗
DOE OSTI2025

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

Analysis of contact conditions and microstructure evolution in shear assisted processing and extrusion using smoothed particle hydrodynamics method

Shear assisted processing and extrusion (ShAPE) is a solid-phase processing technique that adds an additional shear force as compared with a conventional extrusion approach. Recently, ShAPE has demonstrated the capability of extruding high-performance aluminum alloy 7075 (AA7075) tubes at speeds up to 12.2 m/min without surface tearing. However, the relationship among the ShAPE processing parameters, thermomechanical conditions, contact conditions, heat generation, and microstructure evolution remains primarily empirical because an insightful understanding of the associated physics is still lacking. To help elucidate these relationships, this work proposes a thermomechanical meshfree model for the first time for ShAPE processing of AA7075 using the smoothed particle hydrodynamics (SPH) method. The meshfree model is first validated thoroughly by experimental data in terms of material flow, die face temperature, and extrusion force with various processing parameters. The validated model is then used to analyze the steady-state contact conditions and heat generation rates during ShAPE processing. Distributions of the average grain size of AA7075 being extruded are calculated using the SPH model output. The meshfree model results reveal that extrusions conducted at lower temperatures and higher strain rates yield more refined grains and possibly higher material strength, which is also consistent with the experimental observations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Bifurcation of equilibrium positions for ellipsoidal particles in inertial shear flows between two walls

We conducted a systematic numerical investigation of spherical, prolate and oblate particles in an inertial shear flow between two parallel walls, using smoothed particle hydrodynamics (SPH). It was previously shown that above a critical Reynolds number, spherical particles experience a supercritical pitchfork bifurcation of the equilibrium position in shear flow between two parallel walls, namely that the central equilibrium position becomes unstable, leading to the emergence of two new off-centre stable positions (Foxet al.,J. Fluid Mech., vol. 915, 2021). This phenomenon was unexpected given the symmetry of the system. In addition to confirming this finding, we found, surprisingly, that ellipsoidal particles can also return to the centre position from the off-centre positions when the particle Reynolds number is further increased, while spherical particles become unstable under this increased Reynolds number. By utilizing both SPH and the finite element method for flow visualization, we explained the underlining mechanism of this reverse of bifurcation by altered streamwise vorticity and symmetry breaking of pressure. Furthermore, we expanded our investigation to include asymmetric particles, a novel aspect that had not been previously modelled, and we observed similar trends in particle dynamics for both symmetric and asymmetric ellipsoidal particles. While further validation through laboratory experiments is necessary, our research paves the road for development of new focusing and separation methods for shaped particles.

Mechanics↗
DOE OSTI2026

Plutonium(III) versus uranium(III) and samarium(III) in small molecule activation chemistry

We report the PuIII complex, [PuIII(CpMe4)3] (1-Pu), and demonstrate its differences in small molecule reactivity compared to the UIII and SmIII analogs, [UIII(CpMe4)3] (1-U) and [SmIII(CpMe4)3] (1-Sm), respectively. 1-Pu reductively cleaves the small molecule (PhS)2, affording a PuIII complex, [{PuIII(CpMe4)2}2(μ-SPh)2] (2-Pu), while retaining the PuIII center and eliminating (CpMe4)2 as a by-product, a fingerprint of a sterically induced reduction (SIR) reaction. Sm is often used as a surrogate for Pu, but the analogous [SmIII(CpMe4)3], (1-Sm), is unreactive. The (PhS)2 cleavage by 1-U proceeds solely via a metal-based oxidation (i.e., UIII → UIV), to form [UIV(CpMe4)3(SPh)] (3-U). Only 1-U reacts with (PhHN)2, affording the reductive cleavage product, [UIV(CpMe4)3(NHPh)] (4-U). The difference in reactivity of 1-Pu compared to complexes 1-Sm and 1-U was unexpected, and since SIR chemistry can enable complexes to participate in otherwise impossible reductive transformation of substrates, this reinforces the importance of studying small molecule reactivity with the transuranic elements.

Keener, Megan↗