Verification Problem


Astronomy research
  Software Infrastructure:
     My instruments
  White dwarf supernova:
     Stable nickel production
     Remnant metallicities
     Colliding white dwarfs
     Merging white dwarfs
     Ignition conditions
     Metallicity effects
     Central density effects
     Detonation density effects
     Tracer particle burning
     Subsonic burning fronts
     Supersonic burning fronts
     W7 profiles
  Massive star supernova:
     Rotating progenitors
     3D evolution
     26Al & 60Fe
     44Ti, 60Co & 56Ni
     Yields of radionuclides
     Effects of 12C +12C
     SN 1987A light curve
     Constraints on Ni/Fe ratios
     An r-process
  Neutron Stars and Black Holes:
     Black Hole mass spectrum
     Compact object IMF
     He burn on neutron stars
     Variable white dwarfs
     Pop III with JWST
     Neutrino HR diagram
     Monte Carlo massive stars
     Pre-supernova neutrinos
     Pre-supernova variations
     Monte Carlo white dwarfs
     SAGB stars
     Nugrid Yields I
     Classical novae
     He shell convection
     Presolar grains
     BBFH at 40 years
  Chemical Evolution:
     Iron Pseudocarbynes
     Radionuclides in the 2020s
     Hypatia catalog
     Zone models H to Zn
     Mixing ejecta
     γ-rays within 100 Mpc
  Thermodynamics & Networks
     Stellar EOS
     12C(α,γ)16O Rate
     Proton-rich NSE
     Reaction networks
     Bayesian reaction rates
  Verification Problems:
     Validating an astro code
Software instruments
cococubed YouTube
Bicycle adventures
Public Outreach
Education materials
2022 ASU Solar Systems Astronomy
2022 ASU Energy in Everyday Life

AAS Journals
AAS YouTube
2022 Earendel, A Highly Magnified Star
2022 TV Columbae, Micronova
2022 White Dwarfs and 12C(α,γ)16O
2022 MESA VI
2022 MESA in Don't Look Up
2022 MESA Marketplace
2012-2023 MESA Schools
2022 MESA Classroom
2021 Bill Paxton, Tinsley Prize

Contact: F.X.Timmes
my one page vitae,
full vitae,
research statement, and
teaching statement.
Noh's (1987) test case is a standard verification problem. A sphere of gas with a gamma-law equation of state is uniformly compressed, testing the ability to transform kinetic energy into internal energy, and the ability to follow supersonic flows. In the standard Noh problem, a cold gas is initialized with a uniform, radially inward speed of 1 cm s$^{-1}$. A shock forms at the origin and propagates outward as the gas stagnates. For an initial gas density of $\rho_0$ = 1 g cm$^{-3}$, the analytic solution in spherical geometry for $\gamma$ = 5/3 predicts a density in the stagnated gas, i.e., after passage of the outward moving shock, of 64 g cm$^{-3}$ .

Most hydrocode implementations produce anomalous "wall-heating" near the origin. As the shock forms at the origin the momentum equation tries to establish the correct pressure level. However, numerical dissipation generates entropy. The density near the origin drops below the correct value to compensate for the excess internal energy (e.g., Rider 2000). Thus, the density profile is altered near the origin while the pressure profile remains at the correct constant level in the post-shock region. See Gehmeyr, Cheng, & Mihalas (1997) for a remarkable exception. This article, this article, this article, and discuss analytic and numerical solutions for the Noh test case.

The tool in noh.tbz provide solutions as a function of time and position for the RMTV verification test case.

analytical and numerical
convergence study