Journal
SIAM JOURNAL ON SCIENTIFIC COMPUTING
Volume 32, Issue 6, Pages 3564-3583Publisher
SIAM PUBLICATIONS
DOI: 10.1137/100792135
Keywords
implicit methods; finite differencing; monotonicity
Categories
Funding
- NASA [PF8-90054, NAS8-03060]
- Chandra X-ray Center
- U.S. Department of Energy Office of Advanced Scientific Computing Research [DE-AC02-05CH11231]
Ask authors/readers for more resources
The equation governing the streaming of a quantity down its gradient superficially looks similar to the simple constant velocity advection equation. In fact, it is the same as an advection equation if there are no local extrema in the computational domain or at the boundary. However, in general when there are local extrema in the computational domain, it is a nontrivial nonlinear equation. The standard upwind time evolution with a CFL-limited timestep results in spurious oscillations at the grid scale. These oscillations, which originate at the extrema, propagate throughout the computational domain and are undamped even at late times. These oscillations arise because of unphysically large fluxes leaving (entering) the maxima (minima) with the standard CFL-limited explicit methods. Regularization of the equation shows that it is diffusive at the extrema; because of this, an explicit method for the regularized equation with Delta t proportional to Delta x(2) behaves fine. We show that the implicit methods show stable and converging results with Delta t proportional to Delta x; however, surprisingly, even implicit methods are not stable with large enough timesteps. In addition to these subtleties in the numerical implementation, the solutions to the streaming equation are quite novel: nondifferentiable solutions emerge from initially smooth profiles; the solutions show transport over large length scales, e. g., in form of tails. The fluid model for cosmic rays interacting with a thermal plasma (valid at space scales much larger than the cosmic ray Larmor radius) and the equation of saturated conduction in a collisionless plasma are similar to the streaming equation, so our method will find applications in fluid modeling of important processes in plasma astrophysics.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available