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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Herbig Ae/Be Stars toward the Dark Cloud LDN 1667

We report the discovery of a new emission-line object, named SPH 4-South = (GAIA EDR3 5616553300192230272), toward the dark cloud LDN 1667. This object came to our attention after inspecting public images that show a faint diffuse nebula a few arcseconds south of SPH 4, an emission-line object previously classified as a T Tauri star. We present high-resolution spectra and analyzed JHK photometry of SPH 4 and SPH 4-South and new narrowband and archival broadband images of these objects. A comparison of the spectra of SPH 4 and SPH 4-South with high-resolution ones of DG Cir and R Mon strongly suggests that SPH 4 and SPH 4-South are Herbig Ae/Be stars. The classification of SPH 4-South is further supported by using a k-NN algorithm to its position in an H–K versus J–H color–color diagram. Both stars are detected in the four WISE bands and the WISE colors allow us to classify SPH 4 as a Class I and SPH 4-South as a Class II source. We also show that the faint nebula is most probably associated with SPH 4-South. Using published results on LDN 1667 and the Gaia Early Data Release 3, we conclude that SPH 4 is a member of LDN 1667. The case of SPH 4-South is not clear because the determination of its distance and proper motion could be affected by the nebulosity around the star, although membership of SPH 4-South to LDN 1667 cannot be ruled out.

47 OTHER INSTRUMENTATION↗
DOE OSTI2016

Assessment of a Smoothed Particle Hydrodynamics Method for Thin Film on Textured Surfaces

A large number of previous studies suggest that the behavior of the near-wall liquid micro-layer plays a key role in boiling heat transfer and boiling crisis; see examples. The thin liquid layer interacts with heater surface and 'external' two-phase flow. The net effects of these interactions define the micro-hydrodynamics that leads to dryout. As a result, a good understanding of the thin-layer hydrodynamics is significant for the accurate prediction of CHF. Recently, with the advent of nanotechnology including micro-scale manufacturing, it becomes possible to perform high heat-flux boiling experiments with controlled surface conditions. Although many reports show rich observations about the effects of surface characteristics on CHF, physical mechanisms remain not fully understood. In order to have better understanding of the above phenomena, several computational models are developed, ranging from the lubrication theory to CFD methods. The former models are limited to one-dimensional formulation and unable to deal with the dry-out phenomena. It turns out that CFD methods have a greater potential for the mechanistic study of boiling. Towards the study objective, a particle-based CFD method called Smoothed Particle Hydrodynamics (SPH) is adopted. SPH is a Lagrangian particle method for solving systems of partial differential equations. Originally SPH is proposed for fluid dynamics in the context of astrophysical problems. Because mass and density are introduced during the particle approximation, SPH is also applied in fluid dynamic problems, where mass and density are important parameters to determine, without mass loss (a limitation of interface capturing schemes in mesh-based methods). Besides, geometrically complex and/or dynamic boundaries and interfaces can be handled without undue difficulty in SPH. Currently, most SPH simulations focus on large-scale scenarios (meter scale), like flooding and solitary wave. Also some works have been done with SPH in pore-scale. As for micro- or nano-scale simulation, people usually apply Dissipative Particle Dynamics (DPD) and Molecular Dynamics (MD) for particle-based simulation. Almost no work has been done on micro-scale with SPH, so assessments are needed to investigate the capability of SPH in application of micro-hydrodynamics. The Phenomena Identification and Ranking Table (PIRT) is applied to align application's technical requirements with associated physics and materials, models, simulation code capabilities and verification, validation and uncertainty quantification activities. The figures of merit are determined as: (1) Timing of micro-layer rewetting dry surface; (2) Effects of surface features. The corresponding key phenomena are determined as: (1) Hydrodynamics; (2) Interface Dynamics; (3) Energy; (4) Surface Effect. This study deals with isothermal scenario and no energy term is considered. Based on the current needs, a new SPH package, Micro-SPH, is compiled based on SPHysics with surface tension and wettability built in. Test cases are performed with Micro-SPH to assess the capability of SPH to each phenomenon identified. NEUTRINO with incompressible SPH solver is employed to investigate micro-fluid dynamics on a surface with uniformly distributed cylindrical pillars for application purposes. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

COUPLING SMOOTHED PARTICLE HYDRODYNAMICS WITH FINITE ELEMENT METHOD TO SIMULATE RESIDUAL STRESSES FROM FRICTION STIR PROCESSING

Friction stir processing (FSP) is a solid-state material processing technique that locally modifies the microstructure but also induces undesirable residual stresses. A robust numerical model for the FSP can help in mitigating these residual stresses. Heat source models within a finite element method (FEM) framework suffer from inaccuracies. In contrast, smoothed particle hydrodynamics (SPH) model that explicitly captures the material flow near the tool and the associated heat generation are accurate. However, the computational expense of SPH simulations can be prohibitive. In this work, we propose a coupled SPH-FEM framework. SPH is used to model the heat generation accurately near the tool and which is then inserted into to the FEM model as a heat source. To verify this proposed coupling approach, a test case is set up with typical FSP conditions and it is modeled in both SPH and SPH-FEM. The temperatures profiles were compared after the simulations have reached steady-state temperatures. The similarity of the temperature profiles from SPH-FEA and SPH validated the proposed coupling approach. This proposed approach achieves the accuracy of the SPH method while potentially retaining the low computational expense of FEM.

smoothed particle hydrodynamics, Finite Element Me↗
DOE OSTI2014

Smoothed particle hydrodynamics model for Landau-Lifshitz-Navier-Stokes and advection-diffusion equations

We propose a novel smoothed particle hydrodynamics (SPH) discretization of the fully coupled Landau-Lifshitz-Navier-Stokes (LLNS) and stochastic advection-diffusion equations. The accuracy of the SPH solution of the LLNS equations is demonstrated by comparing the scaling of velocity variance and the self-diffusion coefficient with kinetic temperature and particle mass obtained from the SPH simulations and analytical solutions. The spatial covariance of pressure and velocity fluctuations is found to be in a good agreement with theoretical models. To validate the accuracy of the SPH method for coupled LLNS and advection-diffusion equations, we simulate the interface between two miscible fluids. We study formation of the so-called “giant fluctuations” of the front between light and heavy fluids with and without gravity, where the light fluid lies on the top of the heavy fluid. We find that the power spectra of the simulated concentration field are in good agreement with the experiments and analytical solutions. In the absence of gravity, the power spectra decay as the power −4 of the wavenumber—except for small wavenumbers that diverge from this power law behavior due to the effect of finite domain size. Gravity suppresses the fluctuations, resulting in much weaker dependence of the power spectra on the wavenumber. Finally, the model is used to study the effect of thermal fluctuation on the Rayleigh-Taylor instability, an unstable dynamics of the front between a heavy fluid overlaying a light fluid. The front dynamics is shown to agree well with the analytical solutions.

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

A hybrid finite volume method and smoothed particle hydrodynamics approach for efficient and accurate blast simulations

Modeling strong shock waves in fluids remains a persistent challenge in computational physics. Essential to research efforts in industry and defense, numerous methods have been devised to improve the accuracy and efficiency of shock simulations. A novel, hybrid Finite Volume Method (FVM)-Smoothed Particle Hydrodynamics (SPH) approach is capable of further improving efficiency and retaining accuracy by exploiting the favorable characteristics of each respective method. This hybrid approach is presented for shock capturing in compressible fluids. The Python framework Pyro2 is employed to simulate a coarse FVM mesh, while the Python framework PySPH is utilized to model the fluid in regions with high gradients through SPH particles. The performance of the hybrid FVM-SPH scheme, compared to the individual FVM and SPH methods, is assessed in 1 kt and 10 kt blast simulations. Our results indicate that the hybrid approach offers higher computational efficiency than SPH while preserving its accuracy and characteristics. The hybrid approach had a relative speedup of 11.3x and 22.3x over the FVM and SPH approaches for the 1 kt simulation and a relative speedup of 14.7x and 20.9x over the FVM and SPH approaches for the 10 kt simulation. The hybrid SPH algorithm enables future compressible fluid simulations with more extensive capabilities than grid-based methods alone, presenting potential applications in modeling fluid-structure interactions and solid deformation and fracturing in blast simulations.

Myers, Conner↗
DOE OSTI2016

Smooth Particle Hydrodynamics-based Wind Representation

As a result of the 2011 accident at the Fukushima Dai-ichi NPP and other operational NPP experience, there is an identified need to better characterize and evaluate the potential impacts of externally generated hazards on NPP safety. Due to the ubiquitous occurrence of high winds around the world and the possible extreme magnitude of the hazard that has been observed, the assessment of the impact of the high-winds hazard has been identified as an important activity by both NPP owner-operators and regulatory authorities. However, recent experience obtained from the conduct of high-winds risk assessments indicates that such activities have been both labor-intensive and expensive to perform. Additionally, the existing suite of methods and tools to conduct such assessments (which were developed decades ago) do not make use of modern computational architectures (e.g., parallel processing, object-oriented programming techniques, or simple user interfaces) or methods (e.g., efficient and robust numerical-solution schemes). As a result, the current suite of methods and tools will rapidly become obsolete. Physics-based 3D simulation methods can provide information to assist in the RISMC PRA methodology. This research is intended to determine what benefits SPH methods could bring to high-winds simulations for the purposes of assessing their potential impact on NPP safety. The initial investigation has determined that SPH can simulate key areas of high-wind events with reasonable accuracy, compared to other methods. Some problems, such as simulation voids, need to be addressed, but possible solutions have been identified and will be tested with continued work. This work also demonstrated that SPH simulations can provide a means for simulating debris movement; however, further investigations into the capability to determine the impact of high winds and the impacts of wind-driven debris that lead to SSC failures need to be done. SPH simulations alone would be limited in size and computation time. An advanced method of combing results from grid-based methods with SPH through a data-driven model is proposed. This method could allow for more accurate simulation of particle movement near rigid bodies even with larger SPH particle sizes. If successful, the data-driven model would eliminate the need for a SPH turbulence model and increase the simulation domain size. Continued research beyond the scope of this project will be needed in order to determine the viability of a data-driven model.

17 WIND ENERGY↗
DOE OSTI2017

Ligand-Enhanced Optical Response of Gold Nanomolecules and Its Fragment Projection Analysis: The Case of Au 30 (SR) 18

We investigate via first-principles simulations the optical absorption spectra of three different Au 30 (SR) 18 monolayer-protected clusters (MPC): Au 30 (StBu) 18 , Au 30 (SPh) 18 , and Au 30 (SPh-pNO 2 ) 18 . Au 30 (StBu) 18 is known in the literature, and its crystal structure is available. In contrast, Au 30 (SPh) 18 and Au 30 (SPh-pNO 2 ) 18 are two species that have been designed by replacing the tert-butyl organic residues of Au 30 (StBu) 18 with aromatic ones so as to investigate the effects of ligand replacement on the optical response of Au nanomolecules. By analogy to a previously studied Au 23 (SR)16– anionic species, despite distinct differences in charge and chemical composition, a substantial ligand enhancement of the absorption intensity in the optical region is also obtained for the Au 30 (SPh-pNO 2 ) 18 MPC. Furthermore, the use of conjugated aromatic ligands with properly chosen electron-withdrawing substituents and exhibiting steric hindrance so as to also achieve charge decompression at the surface is therefore demonstrated as a general approach to enhancing the MPC photoabsorption intensity in the optical region. In addition, we here subject the ligand-enhancement phenomenon to a detailed analysis based on the fragment projection of electronic excited states and on induced transition densities, leading to a better understanding of the physical origin of this phenomenon, thus opening avenues to its more precise control and exploitation.

36 MATERIALS SCIENCE↗
DOE OSTI2013

THE SUPERMASSIVE BLACK HOLE MASS-SPHEROID STELLAR MASS RELATION FOR SERSIC AND CORE-SERSIC GALAXIES

We have examined the relationship between supermassive black hole mass (M{sub BH}) and the stellar mass of the host spheroid (M{sub sph,*}) for a sample of 75 nearby galaxies. To derive the spheroid stellar masses we used improved Two Micron All Sky Survey K{sub s}-band photometry from the ARCHANGEL photometry pipeline. Dividing our sample into core-Sersic and Sersic galaxies, we find that they are described by very different M{sub BH}-M{sub sph,*} relations. For core-Sersic galaxies-which are typically massive and luminous, with M{sub BH} {approx}> 2 Multiplication-Sign 10{sup 8} M{sub Sun }-we find M{sub BH}{proportional_to} M{sub sph,*}{sup 0.97{+-}0.14}, consistent with other literature relations. However, for the Sersic galaxies-with typically lower masses, M{sub sph,*} {approx}< 3 Multiplication-Sign 10{sup 10} M{sub Sun }-we find M{sub BH}{proportional_to}M{sub sph,*}{sup 2.22{+-}0.58}, a dramatically steeper slope that differs by more than 2 standard deviations. This relation confirms that, for Sersic galaxies, M{sub BH} is not a constant fraction of M{sub sph,*}. Sersic galaxies can grow via the accretion of gas which fuels both star formation and the central black hole, as well as through merging. Their black hole grows significantly more rapidly than their host spheroid, prior to growth by dry merging events that produce core-Sersic galaxies, where the black hole and spheroid grow in lockstep. We have additionally compared our Sersic M{sub BH}-M{sub sph,*} relation with the corresponding relation for nuclear star clusters, confirming that the two classes of central massive object follow significantly different scaling relations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

SANTA BARBARA CLUSTER COMPARISON TEST WITH DISPH

The Santa Barbara cluster comparison project revealed that there is a systematic difference between entropy profiles of clusters of galaxies obtained by Eulerian mesh and Lagrangian smoothed particle hydrodynamics (SPH) codes: mesh codes gave a core with a constant entropy, whereas SPH codes did not. One possible reason for this difference is that mesh codes are not Galilean invariant. Another possible reason is the problem of the SPH method, which might give too much “protection” to cold clumps because of the unphysical surface tension induced at contact discontinuities. In this paper, we apply the density-independent formulation of SPH (DISPH), which can handle contact discontinuities accurately, to simulations of a cluster of galaxies and compare the results with those with the standard SPH. We obtained the entropy core when we adopt DISPH. The size of the core is, however, significantly smaller than those obtained with mesh simulations and is comparable to those obtained with quasi-Lagrangian schemes such as “moving mesh” and “mesh free” schemes. We conclude that both the standard SPH without artificial conductivity and Eulerian mesh codes have serious problems even with such an idealized simulation, while DISPH, SPH with artificial conductivity, and quasi-Lagrangian schemes have sufficient capability to deal with it.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Simulating bubble dynamics in a buoyant system

Summary Multiphase flows are critical components of many physical systems; however, numerical models of multiphase flows with large parameter gradients can be challenging. Here, two different numerical methods, volume of fluid (VOF) and smoothed particle hydrodynamics (SPH), are used to model the buoyant rise of isolated gas bubbles through quiescent fluids for a range of Bond and Reynolds numbers. The VOF is an Eulerian grid–based method, whereas the SPH is Lagrangian and mesh free. Each method has unique strengths and weaknesses, and a comparison of the two approaches as applied to multiphase phenomena has not previously been performed. The VOF and SPH simulations are compared, verified, and validated. Results using two‐dimensional VOF and SPH simulations are similar to each other and are able to reproduce numerical benchmarks and experimental results for sufficiently large Morton and Reynolds numbers. It is also shown that at low Reynolds numbers, the two methods, SPH and VOF, diverge in the transient regime of the bubble rise. Regimes that require simulations capable of representing three‐dimensional drag are identified as well as regimes in which results from VOF and SPH diverge.

Arai, Erin↗
DOE OSTI2021

A continuum meshfree method for sandbox-style numerical modeling of accretionary and doubly vergent wedges

Here, in this study, we introduce and utilize GEOSPH, a Lagrangian particle-based continuum Smoothed Particle Hydrodynamics (SPH) code to simulate the deformation of geologic materials in sandbox-style experiments using an elastoplastic constitutive model. Due to the meshfree nature of SPH, our proposed framework is capable of handling the large deformations characteristic of sandbox-style experiments. SPH is also able to capture strain localization and discontinuities, allowing us to resolve faulting patterns emerging in the sandbox experiments. We successfully validate our SPH implementation against the benchmark experiments of Buiter and co-workers, and then apply SPH to understand the deformational processes occurring in doubly vergent orogens resultant from the S-point experimental setup and from finitely high wedge-like indenter backstops of different inclinations. We find that the topographic slope of the prowedge results from the competing processes of vertical uplift due to the indenter and forward thrusting, while the topographic slope of the retrowedge is a product of repeated shallow slope failure. Our findings show that SPH is a promising numerical method which should be added to the current palette of computational tools used by structural geologists for solving problems relating to tectonics and crustal deformation.

58 GEOSCIENCES↗
DOE OSTI2012

LoCuSS: THE SUNYAEV-ZEL'DOVICH EFFECT AND WEAK-LENSING MASS SCALING RELATION

We present the first weak-lensing-based scaling relation between galaxy cluster mass, M{sub WL}, and integrated Compton parameter Y{sub sph}. Observations of 18 galaxy clusters at z {approx_equal} 0.2 were obtained with the Subaru 8.2 m telescope and the Sunyaev-Zel'dovich Array. The M{sub WL}-Y{sub sph} scaling relations, measured at {Delta} = 500, 1000, and 2500 {rho}{sub c}, are consistent in slope and normalization with previous results derived under the assumption of hydrostatic equilibrium (HSE). We find an intrinsic scatter in M{sub WL} at fixed Y{sub sph} of 20%, larger than both previous measurements of M{sub HSE}-Y{sub sph} scatter as well as the scatter in true mass at fixed Y{sub sph} found in simulations. Moreover, the scatter in our lensing-based scaling relations is morphology dependent, with 30%-40% larger M{sub WL} for undisturbed compared to disturbed clusters at the same Y{sub sph} at r{sub 500}. Further examination suggests that the segregation may be explained by the inability of our spherical lens models to faithfully describe the three-dimensional structure of the clusters, in particular, the structure along the line of sight. We find that the ellipticity of the brightest cluster galaxy, a proxy for halo orientation, correlates well with the offset in mass from the mean scaling relation, which supports this picture. This provides empirical evidence that line-of-sight projection effects are an important systematic uncertainty in lensing-based scaling relations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

THE (BLACK HOLE)-BULGE MASS SCALING RELATION AT LOW MASSES

Several recent papers have reported on the occurrence of active galactic nuclei (AGNs) containing undermassive black holes relative to a linear scaling relation between black hole mass (M {sub bh}) and host spheroid stellar mass (M {sub sph,} {sub *}). However, dramatic revisions to the M {sub bh}-M {sub sph,} {sub *} and M {sub bh}-L {sub sph} relations, based on samples containing predominantly inactive galaxies, have recently identified a new steeper relation at M {sub bh} ≲ (2-10) × 10{sup 8} M {sub ☉}, roughly corresponding to M {sub sph,} {sub *} ≲ (0.3-1) × 10{sup 11} M {sub ☉}. We show that this steeper, quadratic-like M {sub bh}-M {sub sph,} {sub *} relation defined by the Sérsic galaxies, i.e., galaxies without partially depleted cores, roughly tracks the apparent offset of the AGN having 10{sup 5} ≲ M {sub bh}/M {sub ☉} ≲ 0.5 × 10{sup 8}. That is, these AGNs are not randomly offset with low black hole masses, but also follow a steeper (nonlinear) relation. As noted by Busch et al., confirmation or rejection of a possible AGN offset from the steeper M {sub bh}-M {sub sph,} {sub *} relation defined by the Sérsic galaxies will benefit from improved stellar mass-to-light ratios for the spheroids hosting these AGNs. Several implications for formation theories are noted. Furthermore, reasons for possible under- and overmassive black holes, the potential existence of intermediate mass black holes (<10{sup 5} M {sub ☉}), and the new steep (black hole)-(nuclear star cluster) relation, M{sub bh}∝M{sub nc}{sup 2.7±0.7}, are also discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Meshfree simulation and experimental validation of extreme thermomechanical conditions in friction stir extrusion

Abstract Friction stir extrusion (FSE) is a novel solid-phase processing technique that consolidates and extrudes metal powders, flakes, chips, or billets into high-performance parts by plastic deformation, which has the potential to save substantial processing time and energy. Currently, most studies on FSE are experimental and only a few numerical models have been developed to explain and predict the complex physics of the process. In this work, a meshfree simulation framework based on smoothed particle hydrodynamics (SPH) was developed for FSE. Unlike traditional grid-based methods, SPH is a Lagrangian particle-based method that can handle severe material deformations, capture moving interfaces and surfaces, and monitor the field variable histories explicitly without complicated tracking schemes. These aspects of SPH make it attractive for the FSE process, where in situ evolution of field variables is difficult to observe experimentally. To this end, a 3-D, fully thermomechanically coupled SPH model was developed to simulate the FSE of aluminum wires. The developed model was thoroughly validated by comparing the numerically predicted material flow, strain, temperature history, and extrusion force with experimental results for a certain set of process parameters. The validated SPH model can serve as an effective tool to predict and better understand the extreme thermomechanical conditions during the FSE process.

99 GENERAL AND MISCELLANEOUS↗