skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Integral approximations to classical diffusion and smoothed particle hydrodynamics.

Abstract

Abstract not provided.

Authors:
;  [1];  [2]
  1. (Columbia)
  2. (UCFL)
Publication Date:
Research Org.:
Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
Sponsoring Org.:
unfunded
OSTI Identifier:
1183107
Report Number(s):
SAND2014-17136J
537086
DOE Contract Number:
AC04-94AL85000
Resource Type:
Journal Article
Resource Relation:
Journal Name: CMAME journal
Country of Publication:
United States
Language:
English

Citation Formats

Lehoucq, Richard B., Du, Qiang, and Tartakowsky, Dan. Integral approximations to classical diffusion and smoothed particle hydrodynamics.. United States: N. p., 2014. Web.
Lehoucq, Richard B., Du, Qiang, & Tartakowsky, Dan. Integral approximations to classical diffusion and smoothed particle hydrodynamics.. United States.
Lehoucq, Richard B., Du, Qiang, and Tartakowsky, Dan. Fri . "Integral approximations to classical diffusion and smoothed particle hydrodynamics.". United States. doi:.
@article{osti_1183107,
title = {Integral approximations to classical diffusion and smoothed particle hydrodynamics.},
author = {Lehoucq, Richard B. and Du, Qiang and Tartakowsky, Dan},
abstractNote = {Abstract not provided.},
doi = {},
journal = {CMAME journal},
number = ,
volume = ,
place = {United States},
year = {Fri Aug 01 00:00:00 EDT 2014},
month = {Fri Aug 01 00:00:00 EDT 2014}
}
  • The contribution of the paper is the approximation of a classical diffusion operator by an integral equation with a volume constraint. A particular focus is on classical diffusion problems associated with Neumann boundary conditions. By exploiting this approximation, we can also approximate other quantities such as the flux out of a domain. Our analysis of the model equation on the continuum level is closely related to the recent work on nonlocal diffusion and peridynamic mechanics. In particular, we elucidate the role of a volumetric constraint as an approximation to a classical Neumann boundary condition in the presence of physical boundary.more » The volume-constrained integral equation then provides the basis for accurate and robust discretization methods. An immediate application is to the understanding and improvement of the Smoothed Particle Hydrodynamics (SPH) method.« less
  • The contribution of the paper is the approximation of a classical diffusion operator by an integral equation with a volume constraint. A particular focus is on classical diffusion problems associated with Neumann boundary conditions. By exploiting this approximation, we can also approximate other quantities such as the flux out of a domain. Our analysis of the model equation on the continuum level is closely related to the recent work on nonlocal diffusion and peridynamic mechanics. In particular, we elucidate the role of a volumetric constraint as an approximation to a classical Neumann boundary condition in the presence of physical boundary.more » The volume-constrained integral equation then provides the basis for accurate and robust discretization methods. As a result, an immediate application is to the understanding and improvement of the Smoothed Particle Hydrodynamics (SPH) method.« less
  • We propose a novel Smoothed Particle Hydrodynamics (SPH) discretization of the fully-coupled Landau-Lifshitz-Navier-Stokes (LLNS) and advection-diffusion equations. The accuracy of the SPH solution of the LLNS equations is demonstrated by comparing the scaling of velocity variance and 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 are found to be in a good agreement with theoretical models. To validate the accuracy of the SPH method for the coupled LLNS and advection-diffusion equations, we simulate the interface between two miscible fluids. We study the formation ofmore » the so-called giant fluctuations of the front between light and heavy fluids with and without gravity, where the light fluid lays on the top of the heavy fluid. We find that the power spectra of the simulated concentration field is in good agreement with the experiments and analytical solutions. In the absence of gravity the the power spectra decays as the power -4 of the wave number except for small wave numbers which diverge from this power law behavior due to the effect of finite domain size. Gravity suppresses the fluctuations resulting in the much weaker dependence of the power spectra on the wave number. 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 overlying a light fluid. The front dynamics is shown to agree well with the analytical solutions.« less
  • 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 formationmore » 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.« less
  • We perform a series of smoothed particle hydrodynamics simulations of isolated dwarf galaxies to compare different metal mixing models. In particular, we examine the role of diffusion in the production of enriched outflows and in determining the metallicity distributions of gas and stars. We investigate different diffusion strengths by changing the pre-factor of the diffusion coefficient, by varying how the diffusion coefficient is calculated from the local velocity distribution, and by varying whether the speed of sound is included as a velocity term. Stronger diffusion produces a tighter [O/Fe]–[Fe/H] distribution in the gas and cuts off the gas metallicity distributionmore » function at lower metallicities. Diffusion suppresses the formation of low-metallicity stars, even with weak diffusion, and also strips metals from enriched outflows. This produces a remarkably tight correlation between “metal mass-loading” (mean metal outflow rate divided by mean metal production rate) and the strength of diffusion, even when the diffusion coefficient is calculated in different ways. The effectiveness of outflows at removing metals from dwarf galaxies and the metal distribution of the gas is thus dependent on the strength of diffusion. By contrast, we show that the metallicities of stars are not strongly dependent on the strength of diffusion, provided that some diffusion is present.« less