18,494 research outputs found
Coupling techniques for nonlinear hyperbolic equations. IV. Multi-component coupling and multidimensional well-balanced schemes
This series of papers is devoted to the formulation and the approximation of
coupling problems for nonlinear hyperbolic equations. The coupling across an
interface in the physical space is formulated in term of an augmented system of
partial differential equations. In an earlier work, this strategy allowed us to
develop a regularization method based on a thick interface model in one space
variable. In the present paper, we significantly extend this framework and, in
addition, encompass equations in several space variables. This new formulation
includes the coupling of several distinct conservation laws and allows for a
possible covering in space. Our main contributions are, on one hand, the design
and analysis of a well-balanced finite volume method on general triangulations
and, on the other hand, a proof of convergence of this method toward entropy
solutions, extending Coquel, Cockburn, and LeFloch's theory (restricted to a
single conservation law without coupling). The core of our analysis is, first,
the derivation of entropy inequalities as well as a discrete entropy
dissipation estimate and, second, a proof of convergence toward the entropy
solution of the coupling problem.Comment: 37 page
A Hybrid Godunov Method for Radiation Hydrodynamics
From a mathematical perspective, radiation hydrodynamics can be thought of as
a system of hyperbolic balance laws with dual multiscale behavior (multiscale
behavior associated with the hyperbolic wave speeds as well as multiscale
behavior associated with source term relaxation). With this outlook in mind,
this paper presents a hybrid Godunov method for one-dimensional radiation
hydrodynamics that is uniformly well behaved from the photon free streaming
(hyperbolic) limit through the weak equilibrium diffusion (parabolic) limit and
to the strong equilibrium diffusion (hyperbolic) limit. Moreover, one finds
that the technique preserves certain asymptotic limits. The method incorporates
a backward Euler upwinding scheme for the radiation energy density and flux as
well as a modified Godunov scheme for the material density, momentum density,
and energy density. The backward Euler upwinding scheme is first-order accurate
and uses an implicit HLLE flux function to temporally advance the radiation
components according to the material flow scale. The modified Godunov scheme is
second-order accurate and directly couples stiff source term effects to the
hyperbolic structure of the system of balance laws. This Godunov technique is
composed of a predictor step that is based on Duhamel's principle and a
corrector step that is based on Picard iteration. The Godunov scheme is
explicit on the material flow scale but is unsplit and fully couples matter and
radiation without invoking a diffusion-type approximation for radiation
hydrodynamics. This technique derives from earlier work by Miniati & Colella
2007. Numerical tests demonstrate that the method is stable, robust, and
accurate across various parameter regimes.Comment: accepted for publication in Journal of Computational Physics; 61
pages, 15 figures, 11 table
The VOLNA code for the numerical modelling of tsunami waves: generation, propagation and inundation
A novel tool for tsunami wave modelling is presented. This tool has the
potential of being used for operational purposes: indeed, the numerical code
\VOLNA is able to handle the complete life-cycle of a tsunami (generation,
propagation and run-up along the coast). The algorithm works on unstructured
triangular meshes and thus can be run in arbitrary complex domains. This paper
contains the detailed description of the finite volume scheme implemented in
the code. The numerical treatment of the wet/dry transition is explained. This
point is crucial for accurate run-up/run-down computations. Most existing
tsunami codes use semi-empirical techniques at this stage, which are not always
sufficient for tsunami hazard mitigation. Indeed the decision to evacuate
inhabitants is based on inundation maps which are produced with this type of
numerical tools. We present several realistic test cases that partially
validate our algorithm. Comparisons with analytical solutions and experimental
data are performed. Finally the main conclusions are outlined and the
perspectives for future research presented.Comment: 47 pages, 27 figures. Other author's papers can be downloaded at
http://www.lama.univ-savoie.fr/~dutykh
A numerical comparison between degenerate parabolic and quasilinear hyperbolic models of cell movements under chemotaxis
We consider two models which were both designed to describe the movement of
eukaryotic cells responding to chemical signals. Besides a common standard
parabolic equation for the diffusion of a chemoattractant, like chemokines or
growth factors, the two models differ for the equations describing the movement
of cells. The first model is based on a quasilinear hyperbolic system with
damping, the other one on a degenerate parabolic equation. The two models have
the same stationary solutions, which may contain some regions with vacuum. We
first explain in details how to discretize the quasilinear hyperbolic system
through an upwinding technique, which uses an adapted reconstruction, which is
able to deal with the transitions to vacuum. Then we concentrate on the
analysis of asymptotic preserving properties of the scheme towards a
discretization of the parabolic equation, obtained in the large time and large
damping limit, in order to present a numerical comparison between the
asymptotic behavior of these two models. Finally we perform an accurate
numerical comparison of the two models in the time asymptotic regime, which
shows that the respective solutions have a quite different behavior for large
times.Comment: One sentence modified at the end of Section 4, p. 1
Conditions for duality between fluxes and concentrations in biochemical networks
Mathematical and computational modelling of biochemical networks is often
done in terms of either the concentrations of molecular species or the fluxes
of biochemical reactions. When is mathematical modelling from either
perspective equivalent to the other? Mathematical duality translates concepts,
theorems or mathematical structures into other concepts, theorems or
structures, in a one-to-one manner. We present a novel stoichiometric condition
that is necessary and sufficient for duality between unidirectional fluxes and
concentrations. Our numerical experiments, with computational models derived
from a range of genome-scale biochemical networks, suggest that this
flux-concentration duality is a pervasive property of biochemical networks. We
also provide a combinatorial characterisation that is sufficient to ensure
flux-concentration duality. That is, for every two disjoint sets of molecular
species, there is at least one reaction complex that involves species from only
one of the two sets. When unidirectional fluxes and molecular species
concentrations are dual vectors, this implies that the behaviour of the
corresponding biochemical network can be described entirely in terms of either
concentrations or unidirectional fluxes
- …