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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
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 OSTI2015

High-resolution mass models of dwarf galaxies from LITTLE THINGS

We present high-resolution rotation curves and mass models of 26 dwarf galaxies from “Local Irregulars That Trace Luminosity Extremes, The H i Nearby Galaxy Survey” (LITTLE THINGS). LITTLE THINGS is a high-resolution (∼6″ angular; <2.6 km s{sup −1} velocity resolution) Very Large Array H i survey for nearby dwarf galaxies in the local volume within 11 Mpc. The high-resolution H i observations enable us to derive reliable rotation curves of the sample galaxies in a homogeneous and consistent manner. The rotation curves are then combined with Spitzer archival 3.6 μm and ancillary optical U, B, and V images to construct mass models of the galaxies. This high quality multi-wavelength data set significantly reduces observational uncertainties and thus allows us to examine the mass distribution in the galaxies in detail. We decompose the rotation curves in terms of the dynamical contributions by baryons and dark matter (DM) halos, and compare the latter with those of dwarf galaxies from THINGS as well as ΛCDM Smoothed Particle Hydrodynamic (SPH) simulations in which the effect of baryonic feedback processes is included. Being generally consistent with THINGS and simulated dwarf galaxies, most of the LITTLE THINGS sample galaxies show a linear increase of the rotation curve in their inner regions, which gives shallower logarithmic inner slopes α of their DM density profiles. The mean value of the slopes of the 26 LITTLE THINGS dwarf galaxies is α=−0.32±0.24 which is in accordance with the previous results found for low surface brightness galaxies (α=−0.2±0.2) as well as the seven THINGS dwarf galaxies (α=−0.29±0.07). However, this significantly deviates from the cusp-like DM distribution predicted by DM-only ΛCDM simulations. Instead our results are more in line with the shallower slopes found in the ΛCDM SPH simulations of dwarf galaxies in which the effect of baryonic feedback processes is included. In addition, we discuss the central DM distribution of DDO 210 whose stellar mass is relatively low in our sample to examine the scenario of inefficient supernova feedback in low mass dwarf galaxies predicted from recent ΛCDM SPH simulations of dwarf galaxies where central cusps still remain.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

THE COSMIC EVOLUTION OF MASSIVE BLACK HOLES AND GALAXY SPHEROIDS: GLOBAL CONSTRAINTS AT REDSHIFT z {approx}< 1.2

We study observational constraints on the cosmic evolution of the relationships between massive black hole (MBH) mass (M{sub .}) and stellar mass (M{sub *,sph}; or velocity dispersion {sigma}) of a host galaxy/spheroid. Assuming that the M{sub .}-M{sub *,sph} (or M{sub .}-{sigma}) relation evolves with redshift as {proportional_to}(1 + z){sup {Gamma}}, the MBH mass density can be obtained from either the observationally determined galaxy stellar mass functions or velocity dispersion distribution functions over redshift z {approx} 0-1.2 for any given {Gamma}. The MBH mass density at different redshifts can also be inferred from the luminosity function of QSOs/active galactic nuclei (AGNs) provided the radiative efficiency {epsilon} is known. By matching the MBH density inferred from galaxies to that obtained from QSOs/AGNs, we find that {Gamma} = 0.64{sup +0.27}{sub -0.29} for the M{sub .}-M{sub *,sph} relation and {Gamma} = -0.21{sup +0.28}{sub -0.33} for the M{sub .}-{sigma} relation, and {epsilon} = 0.11{sup +0.04}{sub -0.03}. Our results suggest that MBH mass growth precedes bulge mass growth but that the galaxy velocity dispersion does not increase with the mass growth of the bulge after nuclear activity is quenched, which is roughly consistent with the two-phase galaxy formation scenario proposed by Oser et al. in which a galaxy roughly doubles its mass after z = 1 due to accretion and minor mergers while its velocity dispersion drops slightly.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

ACTIVE GALACTIC NUCLEUS FEEDBACK AND ENTROPY INJECTION IN GALAXY CLUSTER CORES

We make the first estimate of non-gravitational energy profiles in galaxy cluster cores (and beyond) based on observational data. Comparing the observed entropy profiles within r{sub 500}, from the Representative XMM-Newton Cluster Structure Survey to simulated base entropy profiles without feedback from both adaptive mesh refinement (AMR) and smoothed particle hydrodynamic (SPH) non-radiative simulations, we estimate the amount of additional non-gravitational energy, E{sub ICM}, contained in the intracluster medium (ICM), as well as the total energy feedback, E{sub Feedback}, from active galactic nuclei (AGNs; the central AGNs in most cases) into the clusters. The total feedback energy scales with the mean spectroscopic temperature as E{sub Feedback}∝T{sub sp}{sup 2.52±0.08} and E{sub Feedback}∝T{sub sp}{sup 2.17±0.11} for the SPH and AMR baseline profiles. The mean non-gravitational energy per particle within r{sub 500} remaining in the ICM after energy lost during cooling is ε{sub ICM} = 2.8 ± 0.8 keV for the SPH theoretical relation and ε{sub ICM} = 1.7 ± 0.9 keV for the AMR theoretical relation. We use the NRAO/VLA Sky Survey source catalog to determine the radio luminosity, L{sub R} , at 1.4 GHz of the central source(s) of our sample. For T{sub sp} > 3 keV, the E{sub Feedback} correlates with L{sub R} , although with different normalization for cool-core and non-cool-core clusters. We show that AGNs could provide a significant portion of the feedback.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Temperature structure of the intracluster medium from smoothed-particle hydrodynamics and adaptive-mesh refinement simulations

Analyses of cosmological hydrodynamic simulations of galaxy clusters suggest that X-ray masses can be underestimated by 10%-30%. The largest bias originates from both violation of hydrostatic equilibrium (HE) and an additional temperature bias caused by inhomogeneities in the X-ray-emitting intracluster medium (ICM). To elucidate this large dispersion among theoretical predictions, we evaluate the degree of temperature structures in cluster sets simulated either with smoothed-particle hydrodynamics (SPH) or adaptive-mesh refinement (AMR) codes. We find that the SPH simulations produce larger temperature variations connected to the persistence of both substructures and their stripped cold gas. This difference is more evident in nonradiative simulations, whereas it is reduced in the presence of radiative cooling. We also find that the temperature variation in radiative cluster simulations is generally in agreement with that observed in the central regions of clusters. Around R {sub 500} the temperature inhomogeneities of the SPH simulations can generate twice the typical HE mass bias of the AMR sample. We emphasize that a detailed understanding of the physical processes responsible for the complex thermal structure in ICM requires improved resolution and high-sensitivity observations in order to extend the analysis to higher temperature systems and larger cluster-centric radii.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

HTTR 3 D Cross Section Generation with Serpent and MAMMOTH

The High Temperature Test Reactor (HTTR) is a graphite moderated thermal reactor with a nominal power output of 30 MWth. The purpose of this work is to prepare a set of neutron cross sections for use in the transient analyses of three Loss of Forced Cooling (LOFC) events planned for this reactor. The cross sections are generated with the Serpent 2 3 D Monte Carlo code to include the axial heterogeneity and the strong axial coupling of this core. The approach yields good agreement of results when compared to previously used 2 D lattice models, which require a fine spatial discretization and, surprisingly, comes at a higher computational cost. In addition, the temperature distribution must be considered to obtain reasonable cross sections and the axial leakage component dominates the spectral effects in the core. Consequently, a temperature model as a function of the fuel and moderator temperatures is adopted, and the core is depleted to 390 effective full-power days (EFPDs). A full tabulation of cross sections is prepared at this burnup point for a variety of fuel and moderator temperatures. Next, the multiphysics reactor application MAMMOTH is used to evaluate the quality of these cross sections using a diffusion solver. As expected, the homogenization error in the cross sections is significant and the Super Homogenization (SPH) correction from MAMMOTH is necessary to preserve key reaction rates. The SPH corrected MAMMOTH results are in excellent agreement with the Serpent results and reproduce both the reference power profile and temperature coefficients at each tabulation point. This fact confirms the accuracy of the SPH corrected cross sections and demonstrates the maturity of using Serpent 2 and MAMMOTH for 3 D cross section generation, even for reactors as complex as the HTTR. Finally, the authors recommend that the cross sections should be parametrized with the local burnup and that a full core burnup calculation with coupled thermal fluids would provide a better estimate of the initial condition for any subsequent transient analysis. Unfortunately, this was not within the scope of this work, and direct analysis with these tabulations should be considered as a first order approximation, at best.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

THE FORMATION OF A MILKY WAY-SIZED DISK GALAXY. I. A COMPARISON OF NUMERICAL METHODS

The long-standing challenge of creating a Milky Way- (MW-) like disk galaxy from cosmological simulations has motivated significant developments in both numerical methods and physical models. We investigate these two fundamental aspects in a new comparison project using a set of cosmological hydrodynamic simulations of an MW-sized galaxy. In this study, we focus on the comparison of two particle-based hydrodynamics methods: an improved smoothed particle hydrodynamics (SPH) code Gadget, and a Lagrangian Meshless Finite-Mass (MFM) code Gizmo. All the simulations in this paper use the same initial conditions and physical models, which include star formation, “energy-driven” outflows, metal-dependent cooling, stellar evolution, and metal enrichment. We find that both numerical schemes produce a late-type galaxy with extended gaseous and stellar disks. However, notable differences are present in a wide range of galaxy properties and their evolution, including star-formation history, gas content, disk structure, and kinematics. Compared to Gizmo, the Gadget simulation produced a larger fraction of cold, dense gas at high redshift which fuels rapid star formation and results in a higher stellar mass by 20% and a lower gas fraction by 10% at z = 0, and the resulting gas disk is smoother and more coherent in rotation due to damping of turbulent motion by the numerical viscosity in SPH, in contrast to the Gizmo simulation, which shows a more prominent spiral structure. Given its better convergence properties and lower computational cost, we argue that the MFM method is a promising alternative to SPH in cosmological hydrodynamic simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Convergence of the Critical Cooling Rate for Protoplanetary Disk Fragmentation Achieved: The Key Role of Numerical Dissipation of Angular Momentum

We carry out simulations of gravitationally unstable disks using smoothed particle hydrodynamics (SPH) and the novel Lagrangian meshless finite mass (MFM) scheme in the GIZMO code. Our aim is to understand the cause of the nonconvergence of the cooling boundary for fragmentation reported in the literature. We run SPH simulations with two different artificial viscosity implementations and compare them with MFM, which does not employ any artificial viscosity. With MFM we demonstrate convergence of the critical cooling timescale for fragmentation at β{sub crit}≈3. Nonconvergence persists in SPH codes. We show how the nonconvergence problem is caused by artificial fragmentation triggered by excessive dissipation of angular momentum in domains with large velocity derivatives. With increased resolution, such domains become more prominent. Vorticity lags behind density, due to numerical viscous dissipation in these regions, promoting collapse with longer cooling times. Such effect is shown to be dominant over the competing tendency of artificial viscosity to diminish with increasing resolution. When the initial conditions are first relaxed for several orbits, the flow is more regular, with lower shear and vorticity in nonaxisymmetric regions, aiding convergence. Yet MFM is the only method that converges exactly. Our findings are of general interest, as numerical dissipation via artificial viscosity or advection errors can also occur in grid-based codes. Indeed, for the FARGO code values of β{sub crit} significantly higher than our converged estimate have been reported in the literature. Finally, we discuss implications for giant planet formation via disk instability.

79 ASTRONOMY AND ASTROPHYSICS↗
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↗
NASA NTRS1994

Tidal disruption of Periodic Comet Shoemaker-Levy 9 and a constraint on its mean density

The apparent tidal disruption of Periodic Comet Shoemaker-Levy 9 (1993e) during a close encounter within approximately 1.62 planetary radii of Jupiter can be used along with theoretical models of tidal disruption to place an upper bound on the density of the predisruption body. Depending on the theoretical model used, these upper bounds range from rho(sub c) less than 0.702 +/- 0.080 g/cu cm for a simple analytical model calibrated by numerical smoothed particle hydrodynamics (SPH) simulations to rho(sub c) less than 1.50 +/- 0.17 g/cu cm for a detailed semianalytical model. The quoted uncertainties stem from an assumed uncertainty in the perijove radius. However, the uncertainty introduced by the different theoretical models is the major source of error; this uncertainty could be eliminated by future SPH simulations specialized to cometary disruptions, including the effects of initially prolate, spinning comets. If the SPH-based upper bound turns out to be most appropriate, it would be consistent with the predisruption body being a comet with a relatively low density and porous structure, as has been asserted previously based on observations of cometary outgassing. Regardless of which upper bound is preferable, the models all agree that the predisruption body could not have been a relatively high-density body, such as an asteroid with rho approximately = 2 g/cu cm.

Boss, Alan P.↗
NASA NTRS1995

Hydrodynamics of binary coalescence. 2: Polytropes with gamma = 5/3

We present a new numerical study of the equilibrium and stability properties of close binary systems. We use the smoothed-particle hydrodynamics (SPH) technique both to construct accurate equilibrium configurations in three dimensions and to follow their hydrodynamic evolution. We adopt a simple polytropic equation of state p = K(sub rho)(exp gamma) with gamma = 5/3 and K = constant within each star, applicable to low-mass degenerate dwarfs as well as low-mass main-sequence stars. For degenerate configurations, we set the two polytropic constants equal, K = K prime, independent of the mass ratio. For main-sequence stars, we adjust K and K prime so as to obtain a simple mass-radius relation of the form R/R prime = M/M prime, where R prime and M prime are the radius and mass of the secondary. Along a sequence of binary equilibrium configurations for two identical stars, we demonstrate the existence of both secular and dynamical instabilities, confirming directly the results of recent analytic work. We use the SPH method to calculate the nonlinear development of the dynamical instability and to determine the final fate of the system. We find that the two stars merge together into a single, rapidly rotating object in just a few orbital periods. Equilibrium sequences are also constructed for systems containing two nonidentical stars. These sequences terminate at a Roche limit, which we can determine very accurately using SPH. For two low-mass main-sequence stars with mass ratio q approximately less than 0.4 we find that the (synchronized) Roche limit configuration is secularly unstable. We discuss the implications of our results for the evolution of double white-dwarf systems and W Ursae Majoris binaries.

Rasio, Frederic A.↗
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 OSTI2018

Advanced Validation Risk-Informed Approach for the Flooding Tool Based Upon Smooth Particle Hydrodynamics: Validation and Development Status of NEUTRINO

Risk-Informed validation is to focus advanced validation on high-priority and critical elements of the RISMC methods and tools that can best be addressed by this targeted approach to validation. Initial validation focus will be on the RISMC flooding analysis computer software based upon Smoothed Particle Hydrodynamics (SPH), NEUTRINO. The NEUTRINO code is a SPH-based CFD software adopted as one of the RISMC modeling tools for external hazards analysis. This document discusses the validation and development activities of NEUTRINO, especially in relation to its simulation capability for flooding hazards analysis. In order to draw meaningful conclusions from the discussion, we first analyzed the potential risks to nuclear power plants in flooding scenarios (e.g., hazard modes, industry/regulatory concerns, etc.). Then, an initial high-level PIRT was proposed to define the key phenomena relevant to the analysis of flooding hazards. Finally, based on the framework of this PIRT analysis, the validation and development efforts that have been conducted to date for NEUTRINO were discussed. This analysis helped to determine the sufficiency, efficiency, and adequacy of the current code development and validation efforts for NEUTRINO.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2017

Effect of Unsaturated Flow Modes on Partitioning Dynamics of Gravity-Driven Flow at a Simple Fracture Intersection: Laboratory Study and Three-Dimensional Smoothed Particle Hydrodynamics Simulations

In this work, we study unsaturated gravity-driven free-surface flows on a synthetic surface intersected by a horizontal fracture and the importance of flow modes (droplet; rivulet) on the partitioning behavior at the fracture intersection. We present laboratory experiments, three-dimensional smoothed particle hydrodynamics (SPH) simulations using a heavily parallelized code, and an analytical solution. The flow-rate-dependent mode switching from droplets to rivulets is reproduced by the SPH model, and the transition ranges agree with the laboratory experiments. We show that flow modes heavily influence the bypass behavior. Flows favoring the formation of droplets exhibit a much stronger bypass capacity compared to rivulet flows, where nearly the whole fluid mass is initially stored within the horizontal fracture. This behavior is demonstrated for a multi-inlet laboratory setup where the inlet- specific flow rate is chosen so that either a droplet or rivulet flow persists. The effect of fluid buffering within the horizontal fracture is presented in terms of dimensionless fracture in- flow so that characteristic scaling regimes can be recovered. For both cases (rivulets and droplets), the flow within the horizontal fracture transitions into a Washburn regime until a critical threshold is reached and the bypass efficiency increases. For rivulet flows, the initial filling of the horizontal fracture is described by classical plug flow. Meanwhile, for droplet flows, a size-dependent partitioning behavior is observed, and the filling of the fracture takes longer. For the case of rivulet flow, we provide an analytical solution that demonstrates the existence of classical Washburn flow within the horizontal fracture.

unsaturated flow, fractured media, smoothed partic↗
DOE OSTI2013

Pore-Scale Simulation of Intragranular Diffusion: Effects of Incomplete Mixing on Macroscopic Manifestations

Diffusive mass transfer into and out of intragranular micropores ("intragranular diffusion") plays an important role in the large-scale transport of some groundwater contaminants. We are interested in understanding the combined effect of pore-scale advection and intragranular diffusion on solute transport at the effective porous medium scale. To study this problem, we have developed a three-dimensional pore-scale model of fluid flow and solute transport that incorporates diffusion into and out of intragranular pore spaces. Our model is based on the Smoothed Particle Hydrodynamics (SPH) simulation method, which represents fluid and solid phases by a mesh-free particle discretization. In the pore spaces, fluid flow is simulated by discretizing the Navier-Stokes equations using the SPH approach. Solute transport is represented by advection, diffusion within the fluid phase, and diffusion between the fluid and solid phases. Our model is implemented on large-scale parallel computing hardware, allowing us to simulate millions of computational particles and represent fully three-dimensional systems of pores and grains with arbitrarily complex physical geometry. We have used this model system to perform numerical experiments using various model porous media systems, which allows us to draw comparisons between macroscopic measures computed from the pore-scale simulations (such as breakthrough curves) and those predicted by macroscopic formulations that assume complete mixing over the representative volume. In this paper we present results of 3D simulations of pore-scale flow and transport, including cases with and without intragranular diffusion, in two model porous media, one with randomly-packed uniform spherical grains and a second with randomly-packed spheres drawn from a binary grain size distribution. Breakthrough curves were computed from the 3D simulations at various transport distances. Comparable breakthrough curves were computed using 1D macroscopic models with parameters determined independently to the degree possible. Based on comparisons of the pore-scale and macroscopic model results, we draw two primary conclusions. First, non-Fickian behavior is persistent and ubiquitous at the scales considered, and most cases are better represented by a multi-rate mass transfer model even when there is no distinct secondary porosity (i.e., no intragranular diffusion). This suggests that diffusive mass transfer processes between preferential flow paths and relatively immobile zones within the primary porosity may have significant impact on transport, particular in low-concentration tails. Second, the application of mass transfer rate parameters based on an assumption of well-mixed concentrations at the pore-scale tend to overestimate the amount of mass transfer that occurs in heterogeneous pore geometries in which preferential flow leads to incomplete pore-scale lateral mixing.

Environmental Molecular Sciences Laboratory↗
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 OSTI2019

Modeling and Computation of Casting Process by Particle Method

Since particle-based computational methods, such as smoothed particle hydrodynamics (SPH) method, are suitable to treat free surfaces and moving boundaries, it is expected that the simulation with particle methods can represent more accurate molten metal flow. Hence, we apply the SPH method to some casting processes and validate the numerical results by comparing with experiments. While fluid behaviors of the transportation and the mold filling processes are expressed sufficiently, ladle pouring process of aluminum alloy is not reproduced by the existing scheme. So we developed a new numerical model of surface oxide film, and our numerical results show good agreement with the experimental data.

36 MATERIALS SCIENCE↗