544 research outputs found
Virtual Element Methods for hyperbolic problems on polygonal meshes
In the present paper we develop the Virtual Element Method for hyperbolic
problems on polygonal meshes, considering the linear wave equations as our
model problem. After presenting the semi-discrete scheme, we derive the
convergence estimates in H^1 semi-norm and L^2 norm. Moreover we develop a
theoretical analysis on the stability for the fully discrete problem by
comparing the Newmark method and the Bathe method. Finally we show the
practical behaviour of the proposed method through a large array of numerical
tests
Complete Subdivision Algorithms, II: Isotopic Meshing of Singular Algebraic Curves
Given a real valued function f(X,Y), a box region B_0 in R^2 and a positive
epsilon, we want to compute an epsilon-isotopic polygonal approximation to the
restriction of the curve S=f^{-1}(0)={p in R^2: f(p)=0} to B_0. We focus on
subdivision algorithms because of their adaptive complexity and ease of
implementation. Plantinga and Vegter gave a numerical subdivision algorithm
that is exact when the curve S is bounded and non-singular. They used a
computational model that relied only on function evaluation and interval
arithmetic. We generalize their algorithm to any bounded (but possibly
non-simply connected) region that does not contain singularities of S. With
this generalization as a subroutine, we provide a method to detect isolated
algebraic singularities and their branching degree. This appears to be the
first complete purely numerical method to compute isotopic approximations of
algebraic curves with isolated singularities
Mimetic Finite Difference methods for Hamiltonian wave equations in 2D
In this paper we consider the numerical solution of the Hamiltonian wave
equation in two spatial dimension. We use the Mimetic Finite Difference (MFD)
method to approximate the continuous problem combined with a symplectic
integration in time to integrate the semi-discrete Hamiltonian system. The main
characteristic of MFD methods, when applied to stationary problems, is to mimic
important properties of the continuous system. This approach, associated with a
symplectic method for the time integration yields a full numerical procedure
suitable to integrate Hamiltonian problems. A complete theoretical analysis of
the method and some numerical simulations are developed in the paper.Comment: 26 pages, 8 figure
Explicit local time-stepping methods for time-dependent wave propagation
Semi-discrete Galerkin formulations of transient wave equations, either with
conforming or discontinuous Galerkin finite element discretizations, typically
lead to large systems of ordinary differential equations. When explicit time
integration is used, the time-step is constrained by the smallest elements in
the mesh for numerical stability, possibly a high price to pay. To overcome
that overly restrictive stability constraint on the time-step, yet without
resorting to implicit methods, explicit local time-stepping schemes (LTS) are
presented here for transient wave equations either with or without damping. In
the undamped case, leap-frog based LTS methods lead to high-order explicit LTS
schemes, which conserve the energy. In the damped case, when energy is no
longer conserved, Adams-Bashforth based LTS methods also lead to explicit LTS
schemes of arbitrarily high accuracy. When combined with a finite element
discretization in space with an essentially diagonal mass matrix, the resulting
time-marching schemes are fully explicit and thus inherently parallel.
Numerical experiments with continuous and discontinuous Galerkin finite element
discretizations validate the theory and illustrate the usefulness of these
local time-stepping methods.Comment: overview paper, typos added, references updated. arXiv admin note:
substantial text overlap with arXiv:1109.448
Reduction of dynamical biochemical reaction networks in computational biology
Biochemical networks are used in computational biology, to model the static
and dynamical details of systems involved in cell signaling, metabolism, and
regulation of gene expression. Parametric and structural uncertainty, as well
as combinatorial explosion are strong obstacles against analyzing the dynamics
of large models of this type. Multi-scaleness is another property of these
networks, that can be used to get past some of these obstacles. Networks with
many well separated time scales, can be reduced to simpler networks, in a way
that depends only on the orders of magnitude and not on the exact values of the
kinetic parameters. The main idea used for such robust simplifications of
networks is the concept of dominance among model elements, allowing
hierarchical organization of these elements according to their effects on the
network dynamics. This concept finds a natural formulation in tropical
geometry. We revisit, in the light of these new ideas, the main approaches to
model reduction of reaction networks, such as quasi-steady state and
quasi-equilibrium approximations, and provide practical recipes for model
reduction of linear and nonlinear networks. We also discuss the application of
model reduction to backward pruning machine learning techniques
- …