1,322 research outputs found
Finite volume methods for unidirectional dispersive wave model
We extend the framework of the finite volume method to dispersive unidirectional water wave propagation in one space dimension. In particular, we consider a KdV–BBM-type equation. Explicit and implicit–explicit Runge–Kutta-type methods are used for time discretizations. The fully discrete schemes are validated by direct comparisons to analytic solutions. Invariants’ conservation properties are also studied. Main applications include important nonlinear phenomena such as dispersive shock wave formation, solitary waves, and their various interaction
Finite volume methods for unidirectional dispersive wave models
We extend the framework of the finite volume method to dispersive
unidirectional water wave propagation in one space dimension. In particular we
consider a KdV-BBM type equation. Explicit and IMEX Runge-Kutta type methods
are used for time discretizations. The fully discrete schemes are validated by
direct comparisons to analytic solutions. Invariants conservation properties
are also studied. Main applications include important nonlinear phenomena such
as dispersive shock wave formation, solitary waves and their various
interactions.Comment: 25 pages, 12 figures, 51 references. Other authors papers can be
downloaded at http://www.lama.univ-savoie.fr/~dutykh
A direct primitive variable recovery scheme for hyperbolic conservative equations: the case of relativistic hydrodynamics
In this article we develop a Primitive Variable Recovery Scheme (PVRS) to
solve any system of coupled differential conservative equations. This method
obtains directly the primitive variables applying the chain rule to the time
term of the conservative equations. With this, a traditional finite volume
method for the flux is applied in order avoid violation of both, the entropy
and "Rankine-Hugoniot" jump conditions. The time evolution is then computed
using a forward finite difference scheme. This numerical technique evades the
recovery of the primitive vector by solving an algebraic system of equations as
it is often used and so, it generalises standard techniques to solve these kind
of coupled systems. The article is presented bearing in mind special
relativistic hydrodynamic numerical schemes with an added pedagogical view in
the appendix section in order to easily comprehend the PVRS. We present the
convergence of the method for standard shock-tube problems of special
relativistic hydrodynamics and a graphical visualisation of the errors using
the fluctuations of the numerical values with respect to exact analytic
solutions. The PVRS circumvents the sometimes arduous computation that arises
from standard numerical methods techniques, which obtain the desired primitive
vector solution through an algebraic polynomial of the charges.Comment: 19 pages, 6 figures, 2 tables. Accepted for publication in PLOS ON
HARM: A Numerical Scheme for General Relativistic Magnetohydrodynamics
We describe a conservative, shock-capturing scheme for evolving the equations
of general relativistic magnetohydrodynamics. The fluxes are calculated using
the Harten, Lax, and van Leer scheme. A variant of constrained transport,
proposed earlier by T\'oth, is used to maintain a divergence free magnetic
field. Only the covariant form of the metric in a coordinate basis is required
to specify the geometry. We describe code performance on a full suite of test
problems in both special and general relativity. On smooth flows we show that
it converges at second order. We conclude by showing some results from the
evolution of a magnetized torus near a rotating black hole.Comment: 38 pages, 18 figures, submitted to Ap
Uncertainty Quantification for Linear Hyperbolic Equations with Stochastic Process or Random Field Coefficients
In this paper hyperbolic partial differential equations with random
coefficients are discussed. Such random partial differential equations appear
for instance in traffic flow problems as well as in many physical processes in
random media. Two types of models are presented: The first has a time-dependent
coefficient modeled by the Ornstein--Uhlenbeck process. The second has a random
field coefficient with a given covariance in space. For the former a formula
for the exact solution in terms of moments is derived. In both cases stable
numerical schemes are introduced to solve these random partial differential
equations. Simulation results including convergence studies conclude the
theoretical findings
A New Spherical Harmonics Scheme for Multi-Dimensional Radiation Transport I: Static Matter Configurations
Recent work by McClarren & Hauck [29] suggests that the filtered spherical
harmonics method represents an efficient, robust, and accurate method for
radiation transport, at least in the two-dimensional (2D) case. We extend their
work to the three-dimensional (3D) case and find that all of the advantages of
the filtering approach identified in 2D are present also in the 3D case. We
reformulate the filter operation in a way that is independent of the timestep
and of the spatial discretization. We also explore different second- and
fourth-order filters and find that the second-order ones yield significantly
better results. Overall, our findings suggest that the filtered spherical
harmonics approach represents a very promising method for 3D radiation
transport calculations.Comment: 29 pages, 13 figures. Version matching the one in Journal of
Computational Physic
- …