DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: CRKSPH – A Conservative Reproducing Kernel Smoothed Particle Hydrodynamics Scheme

Abstract

We present a formulation of smoothed particle hydrodynamics (SPH) that utilizes a first-order consistent reproducing kernel, a smoothing function that exactly interpolates linear fields with particle tracers. Previous formulations using reproducing kernel (RK) interpolation have had difficulties maintaining conservation of momentum due to the fact the RK kernels are not, in general, spatially symmetric. Here in this paper, we utilize a reformulation of the fluid equations such that mass, linear momentum, and energy are all rigorously conserved without any assumption about kernel symmetries, while additionally maintaining approximate angular momentum conservation. Our approach starts from a rigorously consistent interpolation theory, where we derive the evolution equations to enforce the appropriate conservation properties, at the sacrifice of full consistency in the momentum equation. Additionally, by exploiting the increased accuracy of the RK method's gradient, we formulate a simple limiter for the artificial viscosity that reduces the excess diffusion normally incurred by the ordinary SPH artificial viscosity. Collectively, we call our suite of modifications to the traditional SPH scheme Conservative Reproducing Kernel SPH, or CRKSPH. CRKSPH retains many benefits of traditional SPH methods (such as preserving Galilean invariance and manifest conservation of mass, momentum, and energy) while improving on many of the shortcomingsmore » of SPH, particularly the overly aggressive artificial viscosity and zeroth-order inaccuracy. We compare CRKSPH to two different modern SPH formulations (pressure based SPH and compatibly differenced SPH), demonstrating the advantages of our new formulation when modeling fluid mixing, strong shock, and adiabatic phenomena.« less

Authors:
; ;
Publication Date:
Research Org.:
Lawrence Livermore National Laboratory (LLNL), Livermore, CA (United States)
Sponsoring Org.:
USDOE National Nuclear Security Administration (NNSA)
OSTI Identifier:
1337531
Alternate Identifier(s):
OSTI ID: 1468909
Report Number(s):
LLNL-JRNL-690757
Journal ID: ISSN 0021-9991; S0021999116306453; PII: S0021999116306453
Grant/Contract Number:  
AC02-06CH11357; AC52-07NA27344
Resource Type:
Published Article
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Name: Journal of Computational Physics Journal Volume: 332 Journal Issue: C; Journal ID: ISSN 0021-9991
Publisher:
Elsevier
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING; 72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS; Hydrodynamics; Meshfree

Citation Formats

Frontiere, Nicholas, Raskin, Cody D., and Owen, J. Michael. CRKSPH – A Conservative Reproducing Kernel Smoothed Particle Hydrodynamics Scheme. United States: N. p., 2017. Web. doi:10.1016/j.jcp.2016.12.004.
Frontiere, Nicholas, Raskin, Cody D., & Owen, J. Michael. CRKSPH – A Conservative Reproducing Kernel Smoothed Particle Hydrodynamics Scheme. United States. https://doi.org/10.1016/j.jcp.2016.12.004
Frontiere, Nicholas, Raskin, Cody D., and Owen, J. Michael. Wed . "CRKSPH – A Conservative Reproducing Kernel Smoothed Particle Hydrodynamics Scheme". United States. https://doi.org/10.1016/j.jcp.2016.12.004.
@article{osti_1337531,
title = {CRKSPH – A Conservative Reproducing Kernel Smoothed Particle Hydrodynamics Scheme},
author = {Frontiere, Nicholas and Raskin, Cody D. and Owen, J. Michael},
abstractNote = {We present a formulation of smoothed particle hydrodynamics (SPH) that utilizes a first-order consistent reproducing kernel, a smoothing function that exactly interpolates linear fields with particle tracers. Previous formulations using reproducing kernel (RK) interpolation have had difficulties maintaining conservation of momentum due to the fact the RK kernels are not, in general, spatially symmetric. Here in this paper, we utilize a reformulation of the fluid equations such that mass, linear momentum, and energy are all rigorously conserved without any assumption about kernel symmetries, while additionally maintaining approximate angular momentum conservation. Our approach starts from a rigorously consistent interpolation theory, where we derive the evolution equations to enforce the appropriate conservation properties, at the sacrifice of full consistency in the momentum equation. Additionally, by exploiting the increased accuracy of the RK method's gradient, we formulate a simple limiter for the artificial viscosity that reduces the excess diffusion normally incurred by the ordinary SPH artificial viscosity. Collectively, we call our suite of modifications to the traditional SPH scheme Conservative Reproducing Kernel SPH, or CRKSPH. CRKSPH retains many benefits of traditional SPH methods (such as preserving Galilean invariance and manifest conservation of mass, momentum, and energy) while improving on many of the shortcomings of SPH, particularly the overly aggressive artificial viscosity and zeroth-order inaccuracy. We compare CRKSPH to two different modern SPH formulations (pressure based SPH and compatibly differenced SPH), demonstrating the advantages of our new formulation when modeling fluid mixing, strong shock, and adiabatic phenomena.},
doi = {10.1016/j.jcp.2016.12.004},
journal = {Journal of Computational Physics},
number = C,
volume = 332,
place = {United States},
year = {Wed Mar 01 00:00:00 EST 2017},
month = {Wed Mar 01 00:00:00 EST 2017}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record
https://doi.org/10.1016/j.jcp.2016.12.004

Citation Metrics:
Cited by: 46 works
Citation information provided by
Web of Science

Save / Share:

Works referencing / citing this record:

Modification of the LSMPS method for the conservation of the thermal energy in laser irradiation processes
journal, September 2018

  • Tanaka, Masayuki; Cardoso, Rui; Bahai, Hamid
  • International Journal for Numerical Methods in Engineering, Vol. 117, Issue 2
  • DOI: 10.1002/nme.5952

A Two-Phase SPH Model for Dynamic Contact Angles Including Fluid–Solid Interactions at the Contact Line
journal, January 2018


Effects of propagation distance and half angle on the merging of hypervelocity plasma jets
journal, May 2019

  • Thompson, Seth; Cassibry, Jason
  • Physics of Plasmas, Vol. 26, Issue 5
  • DOI: 10.1063/1.5074100

The Kelvin–Helmholtz instability and smoothed particle hydrodynamics
journal, July 2019

  • Tricco, Terrence S.
  • Monthly Notices of the Royal Astronomical Society, Vol. 488, Issue 4
  • DOI: 10.1093/mnras/stz2042

The Borg Cube Simulation: Cosmological Hydrodynamics with CRK-SPH
journal, May 2019

  • Emberson, J. D.; Frontiere, Nicholas; Habib, Salman
  • The Astrophysical Journal, Vol. 877, Issue 2
  • DOI: 10.3847/1538-4357/ab1b31

Primordial Earth Mantle Heterogeneity Caused by the Moon-forming Giant Impact?
journal, December 2019

  • Deng, Hongping; Ballmer, Maxim D.; Reinhardt, Christian
  • The Astrophysical Journal, Vol. 887, Issue 2
  • DOI: 10.3847/1538-4357/ab50b9