Skip to content
DraconiticResearch
Search Sources API
186,631 studies31 sourcesUpdated 40 min ago
PreprintCC BY 4.0arXiv

A Split-Flux Method for Dust-Gas Two-Fluid Simulations from Strong to Weak Coupling

Abstract

We develop a new method for simulating pressureless fluids of neutral dust coupled to magnetohydrodynamic gas through dust-gas drag. It has been pointed out that, in the strong-coupling regime, dust-gas multifluid methods suffer from artificial variations in the dust-to-gas mass ratio when the dust and gas are evolved with independent numerical fluxes. This problem arises from inconsistencies between the dust and gas mass fluxes. In the strong-coupling limit, the dust-to-gas mass ratio is advected with the local gas velocity. Recently, a new method based on the Harten-Lax-van Leer (HLL) solver was proposed. This method constructs the dust mass flux so that it is consistent with the gas mass flux, and successfully reduces artificial variations in the dust-to-gas mass ratio. Here, we develop a Split-Flux method that can be coupled to an HLLD gas solver. Through various numerical experiments, we show that the Split-Flux method suppresses artificial variations in the dust-to-gas mass ratio, while preserving the low numerical diffusivity of the HLLD solver.

Related

  1. A comparison of particle and fluid approaches for modeling the streaming instability in global patch simulations
    Dominik Hetzel Ostertag, Mario Flock, Paul Rosendahl +3
    Astronomy & AstrophysicsPeer-reviewedPublisher termsSolar & stellar physicsCrossref
  2. Coupling between stellar and HI lopsidedness in Milky Way--type galaxies from the Auriga Superstars cosmological simulations
    Arianna Dolfi, Facundo A. Gómez, Rebekka Bieri +5
    Astronomy & AstrophysicsPeer-reviewedPublisher termsGalaxiesCrossref
  3. Horndeski Gravity as an Effective Dark Energy Fluid in Cosmological N-body Simulations
    Vicente Pedreros Aballay
    arXivPreprintarXiv licenceFull textCosmologyarXiv
  4. What can strong lensing recover in edge-on disc galaxies? A mock study based on TNG50 simulations
    Zheng Fang, Xiaoyue Cao, Ran Li +1
    arXivPreprintarXiv licenceFull textGalaxiesarXiv
  5. A comparison between particle and fluid approaches for streaming instability in global patch simulations
    Dominik Hetzel Ostertag, Mario Flock, Paul Rosendahl +3
    arXivPreprintarXiv licenceFull textInstrumentation & methodsPlanetary science & exoplanetsarXiv
BibTeX
@misc{iwasaki2026splitflux,
  title = {{A Split-Flux Method for Dust-Gas Two-Fluid Simulations from Strong to Weak Coupling}},
  author = {Kazunari Iwasaki},
  howpublished = {arXiv},
  year = {2026},
  month = sep,
  doi = {10.48550/arxiv.2607.24418},
  eprint = {2607.24418},
  archivePrefix = {arXiv},
  url = {https://arxiv.org/abs/2607.24418},
  note = {CC BY 4.0},
}
RIS
TY  - UNPB
TI  - A Split-Flux Method for Dust-Gas Two-Fluid Simulations from Strong to Weak Coupling
AU  - Kazunari Iwasaki
PY  - 2026
DA  - 2026/09/23
PB  - arXiv
DO  - 10.48550/arxiv.2607.24418
UR  - https://arxiv.org/abs/2607.24418
AB  - We develop a new method for simulating pressureless fluids of neutral dust coupled to magnetohydrodynamic gas through dust-gas drag. It has been pointed out that, in the strong-coupling regime, dust-gas multifluid methods suffer from artificial variations in the dust-to-gas mass ratio when the dust and gas are evolved with independent numerical fluxes. This problem arises from inconsistencies between the dust and gas mass fluxes. In the strong-coupling limit, the dust-to-gas mass ratio is advected with the local gas velocity. Recently, a new method based on the Harten-Lax-van Leer (HLL) solver was proposed. This method constructs the dust mass flux so that it is consistent with the gas mass flux, and successfully reduces artificial variations in the dust-to-gas mass ratio. Here, we develop a Split-Flux method that can be coupled to an HLLD gas solver. Through various numerical experiments, we show that the Split-Flux method suppresses artificial variations in the dust-to-gas mass ratio, while preserving the low numerical diffusivity of the HLLD solver.
ER  -