Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Modeling electrokinetic flows by consistent implicit incompressible smoothed particle hydrodynamics

We present an efficient implicit incompressible smoothed particle hydrodynamics (I2SPH) discretization of Navier-Stokes, Poisson-Boltzmann, and advection-diffusion equations subject to Dirichlet or Robin boundary conditions. It is applied to model various two and three dimensional electrokinetic flows in simple or complex geometries. The I2SPH's accuracy and convergence are examined via comparison with analytical solutions, grid-based numerical solutions, or empirical models. The new method provides a framework to explore broader applications of SPH in microfluidics and complex fluids with charged objects, such as colloids and biomolecules, in arbitrary complex geometries.

Pan, Wenxiao↗
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 OSTI2012

Mesoscale Simulation of Proton Transport in Proton Exchange Membranes

Previous efforts to model proton transport through fuel cell membranes have largely focused on disparate length scales: molecular dynamics at the atomistic level and fuel cell stack engineering approaches at the macroscale. A new multiscale approach to bridge these extremes is proposed in this work which combines concepts from coarse-grained (CG) modeling with smoothed particle hydrodynamics (SPH) to capture the qualitative morphology and transport behavior of a proton exchange membrane at the length scale of tens of nanometers. This method allows for connection to atomistic simulations via the inclusion of transport coefficients from molecular dynamics and coarse-grained forces derived for the polymer backbone, side chain, proton, and water interactions. Information pertaining to macroscopic conductivity is obtained by volume averaging based on the flux and chemical potential fields within the membrane. Proton transport is effectively coarse-grained via introduction of a composition variable associated with each interacting site which carries the field information. By combining this technique with local electrostatics and coordinate dependent diffusion constants, the effects of double layer formation within the water pores and the influence of proximity to sulfonate groups on transport is recovered. The combined CG-SPH method is validated and subsequently applied to an equilibrated hydrated Nafion structure with a box length of 40 nm. The resulting conductivities calculated for the material agree very well with trends from experiment and provide insight into the complex interplay of morphology, proton distribution, and diffusion coefficients at a length scale that can be expanded beyond feasible atomistic molecular dynamics simulations to capture the effects of mesoscopic morphology on proton conduction.

Jorn, Ryan↗
DOE OSTI2012

Tendency of spherically imploding plasma liners formed by merging plasma jets to evolve toward spherical symmetry

We have performed three-dimensional (3D) simulations using smoothed particle hydrodynamics (SPH) in order to study the effects of discrete plasma jets on the processes of plasma liner formation, implosion on vacuum, and expansion. It was found that the pressure histories of the inner portion of the liner from 3D SPH simulations with a uniform liner and with 30 discrete plasma jets were qualitatively and quantitatively similar from peak compression through the complete stagnation of the liner. The 3D simulations with a uniform liner were first benchmarked against results from one-dimensional radiation-hydrodynamic simulations [T. J. Awe et al., Phys. Plasmas 18, 072705 (2011)]. Two-dimensional plots of the pressure field show that the discrete jet SPH case evolves towards a profile that is almost indistinguishable from the SPH case with a uniform liner, thus indicating that non-uniformities due to discrete jets are smeared out by late stages of the implosion. The processes of plasma liner formation and implosion on vacuum were shown to be robust against Rayleigh-Taylor instability growth. Finally, interparticle mixing for a liner imploding on vacuum was investigated. The mixing rate was found to be very small until after the peak compression for the 30 jet simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

NUMERICAL CONVERGENCE IN SMOOTHED PARTICLE HYDRODYNAMICS

We study the convergence properties of smoothed particle hydrodynamics (SPH) using numerical tests and simple analytic considerations. Our analysis shows that formal numerical convergence is possible in SPH only in the joint limit N → ∞, h → 0, and N{sub nb} → ∞, where N is the total number of particles, h is the smoothing length, and N{sub nb} is the number of neighbor particles within the smoothing volume used to compute smoothed estimates. Previous work has generally assumed that the conditions N → ∞ and h → 0 are sufficient to achieve convergence, while holding N{sub nb} fixed. We demonstrate that if N{sub nb} is held fixed as the resolution is increased, there will be a residual source of error that does not vanish as N → ∞ and h → 0. Formal numerical convergence in SPH is possible only if N{sub nb} is increased systematically as the resolution is improved. Using analytic arguments, we derive an optimal compromise scaling for N{sub nb} by requiring that this source of error balance that present in the smoothing procedure. For typical choices of the smoothing kernel, we find N{sub nb} ∝N {sup 0.5}. This means that if SPH is to be used as a numerically convergent method, the required computational cost does not scale with particle number as O(N), but rather as O(N {sup 1} {sup +} {sup δ}), where δ ≈ 0.5, with a weak dependence on the form of the smoothing kernel.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Numerical calculation of protein-ligand binding rates through solution of the Smoluchowski equation using smooth particle hydrodynamics

We demonstrate the accuracy and effectiveness of a Lagrangian particle-based method, smoothed particle hydrodynamics (SPH), to study diffusion in biomolecular systems by numerically solving the time-dependent Smoluchowski equation for continuum diffusion. The numerical method is first verified in simple systems and then applied to the calculation of ligand binding to an acetylcholinesterase monomer. Unlike previous studies, a reactive Robin boundary condition (BC), rather than the absolute absorbing (Dirichlet) boundary condition, is considered on the reactive boundaries. This new boundary condition treatment allows for the analysis of enzymes with "imperfect" reaction rates. Rates for inhibitor binding to mAChE are calculated at various ionic strengths and compared with experiment and other numerical methods. We find that imposition of the Robin BC improves agreement between calculated and experimental reaction rates. Although this initial application focuses on a single monomer system, our new method provides a framework to explore broader applications of SPH in larger-scale biomolecular complexes by taking advantage of its Lagrangian particle-based nature.

Pan, Wenxiao↗
DOE OSTI2016

Comparing simulations of AGN feedback

We perform adaptive mesh refinement (AMR) and smoothed particle hydrodynamics (SPH) cosmological zoom simulations of a region around a forming galaxy cluster, comparing the ability of the methods to handle successively more complex baryonic physics. In the simplest, non-radiative case, the two methods are in good agreement with each other, but the SPH simulations generate central cores with slightly lower entropies and virial shocks at slightly larger radii, consistent with what has been seen in previous studies. The inclusion of radiative cooling, star formation, and stellar feedback leads to much larger differences between the two methods. Most dramatically, at z=5, rapid cooling in the AMR case moves the accretion shock to well within the virial radius, while this shock remains near the virial radius in the SPH case, due to excess heating, coupled with poorer capturing of the shock width. On the other hand, the addition of feedback from active galactic nuclei (AGNs) to the simulations results in much better agreement between the methods. For our AGN model, both simulations display halo gas entropies of 100 keV cm{sup 2}, similar decrements in the star formation rate, and a drop in the halo baryon content of roughly 30%. This is consistent with the AGN growth being self-regulated, regardless of the numerical method. However, the simulations with AGN feedback continue to differ in aspects that are not self-regulated, such that in SPH a larger volume of gas is impacted by feedback, and the cluster still has a lower entropy central core.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1993

Smoothed particle hydrodynamics with GRAPE-1A

We describe the implementation of a smoothed particle hydrodynamics (SPH) scheme using GRAPE-1A, a special-purpose processor used for gravitational N-body simulations. The GRAPE-1A calculates the gravitational force exerted on a particle from all other particles in a system, while simultaneously making a list of the nearest neighbors of the particle. It is found that GRAPE-1A accelerates SPH calculations by direct summation by about two orders of magnitudes for a ten thousand-particle simulation. The effective speed is 80 Mflops, which is about 30 percent of the peak speed of GRAPE-1A. Also, in order to investigate the accuracy of GRAPE-SPH, some test simulations were executed. We found that the force and position errors are smaller than those due to representing a fluid by a finite number of particles. The total energy and momentum were conserved within 0.2-0.4 percent and 2-5 x 10 exp -5, respectively, in simulations with several thousand particles. We conclude that GRAPE-SPH is quite effective and sufficiently accurate for self-gravitating hydrodynamics.

Umemura, Masayuki↗
NASA NTRS1994

A comparison of cosmological hydrodynamic codes

We present a detailed comparison of the simulation results of various hydrodynamic codes. Starting with identical initial conditions based on the cold dark matter scenario for the growth of structure, with parameters h = 0.5 Omega = Omega(sub b) = 1, and sigma(sub 8) = 1, we integrate from redshift z = 20 to z = O to determine the physical state within a representative volume of size L(exp 3) where L = 64 h(exp -1) Mpc. Five indenpendent codes are compared: three of them Eulerian mesh-based and two variants of the smooth particle hydrodynamics 'SPH' Lagrangian approach. The Eulerian codes were run at N(exp 3) = (32(exp 3), 64(exp 3), 128(exp 3), and 256(exp 3)) cells, the SPH codes at N(exp 3) = 32(exp 3) and 64(exp 3) particles. Results were then rebinned to a 16(exp 3) grid with the exception that the rebinned data should converge, by all techniques, to a common and correct result as N approaches infinity. We find that global averages of various physical quantities do, as expected, tend to converge in the rebinned model, but that uncertainites in even primitive quantities such as (T), (rho(exp 2))(exp 1/2) persists at the 3%-17% level achieve comparable and satisfactory accuracy for comparable computer time in their treatment of the high-density, high-temeprature regions as measured in the rebinned data; the variance among the five codes (at highest resolution) for the mean temperature (as weighted by rho(exp 2) is only 4.5%. Examined at high resolution we suspect that the density resolution is better in the SPH codes and the thermal accuracy in low-density regions better in the Eulerian codes. In the low-density, low-temperature regions the SPH codes have poor accuracy due to statiscal effects, and the Jameson code gives the temperatures which are too high, due to overuse of artificial viscosity in these high Mach number regions. Overall the comparison allows us to better estimate errors; it points to ways of improving this current generation ofhydrodynamic codes and of suiting their use to problems which exploit their best individual features.

Kang, Hyesung↗
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 OSTI2016

Splenic Arterial Embolization in the Treatment of Severe Portal Hypertension Due to Pancreatic Diseases: The Primary Experience in 14 Patients

ObjectiveThis retrospective study reports our experience using splenic arterial particle embolization and coil embolization for the treatment of sinistral portal hypertension (SPH) in patients with and without gastric bleeding.MethodsFrom August 2009 to May 2012, 14 patients with SPH due to pancreatic disease were diagnosed and treated with splenic arterial embolization. Two different embolization strategies were applied; either combined distal splenic bed particle embolization and proximal splenic artery coil embolization in the same procedure for acute hemorrhage (1-step) or interval staged distal embolization and proximal embolization in the stable patient (2-step). The patients were clinically followed.ResultsIn 14 patients, splenic arterial embolization was successful. The one-step method was performed in three patients suffering from massive gastric bleeding, and the bleeding was relieved after embolization. The two-step method was used in 11 patients, who had chronic gastric variceal bleeding or gastric varices only. The gastric varices disappeared in the enhanced CT scan and the patients had no gastric bleeding during follow-up.ConclusionsSplenic arterial embolization, particularly the two-step method, proved feasible and effective for the treatment of SPH patients with gastric varices or gastric variceal bleeding.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2019

A Consistent Spatially Adaptive Smoothed Particle Hydrodynamics Method for fluid-structure interactions

A new consistent, spatially adaptive, smoothed particle hydrodynamics (SPH) method for Fluid-Structure Interactions (FSI) is presented. The method combines several attributes that have not been simultaneously satisfied by other SPH methods. Specifically, it is second-order convergent; it allows for resolutions spatially adapted with moving (translating and rotating) boundaries of arbitrary geometries; and, it accelerates the FSI solution as the adaptive approach leads to fewer degrees of freedom without sacrificing accuracy. The key ingredients in the method are a consistent discretization of differential operators, a posteriori error estimator/distance-based criterion of adaptivity, and a particle-shifting technique. The method is applied in simulating six different flows or FSI problems. The new method's convergence, accuracy, and efficiency attributes are assessed by comparing the results it produces with analytical, finite element, and consistent SPH uniform high-resolution solutions as well as experimental data.

Hu, Wei↗
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 OSTI2014

Smoothed Particle Hydrodynamics Continuous Boundary Force method for Navier-Stokes equations subject to Robin boundary condition

Robin boundary condition for the Navier-Stokes equations is used to model slip conditions at the fluid-solid boundaries. A novel Continuous Boundary Force (CBF) method is proposed for solving the Navier-Stokes equations subject to Robin boundary condition. In the CBF method, the Robin boundary condition at boundary is replaced by the homogeneous Neumann boundary condition at the boundary and a volumetric force term added to the momentum conservation equation. Smoothed Particle Hydrodynamics (SPH) method is used to solve the resulting Navier-Stokes equations. We present solutions for two-dimensional and three-dimensional flows in domains bounded by flat and curved boundaries subject to various forms of the Robin boundary condition. The numerical accuracy and convergence are examined through comparison of the SPH-CBF results with the solutions of finite difference or finite element method. Taken the no-slip boundary condition as a special case of slip boundary condition, we demonstrate that the SPH-CBF method describes accurately both no-slip and slip conditions.

Pan, Wenxiao↗
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↗