3,666 research outputs found
Multi-revolution composition methods for highly oscillatory differential equations
We introduce a new class of multi-revolution composition methods for the approximation of the N th-iterate of a given near-identity map. When applied to the numerical integration of highly oscillatory systems of differential equations, the technique benefits from the properties of standard composition methods: it is intrinsically geometric and well-suited for Hamiltonian or divergence-free equations for instance. We prove error estimates with error constants that are independent of the oscillatory frequency. Numerical experiments, in particular for the nonlinear Schrödinger equation, illustrate the theoretical results, as well as the efficiency and versatility of the methods
Highly-oscillatory evolution equations with multiple frequencies: averaging and numerics
International audienceIn this paper, we are concerned with the application of the recently introduced multi-revolution composition methods, on the one hand, and two-scale methods, on the other hand, to a class of highly-oscillatory evolution equations with multiple frequencies. The main idea relies on a well-balanced reformulation of the problem as an equivalent mono-frequency equation which allows for the use of the two aforementioned techniques
On-surface radiation condition for multiple scattering of waves
The formulation of the on-surface radiation condition (OSRC) is extended to
handle wave scattering problems in the presence of multiple obstacles. The new
multiple-OSRC simultaneously accounts for the outgoing behavior of the wave
fields, as well as, the multiple wave reflections between the obstacles. Like
boundary integral equations (BIE), this method leads to a reduction in
dimensionality (from volume to surface) of the discretization region. However,
as opposed to BIE, the proposed technique leads to boundary integral equations
with smooth kernels. Hence, these Fredholm integral equations can be handled
accurately and robustly with standard numerical approaches without the need to
remove singularities. Moreover, under weak scattering conditions, this approach
renders a convergent iterative method which bypasses the need to solve single
scattering problems at each iteration.
Inherited from the original OSRC, the proposed multiple-OSRC is generally a
crude approximate method. If accuracy is not satisfactory, this approach may
serve as a good initial guess or as an inexpensive pre-conditioner for Krylov
iterative solutions of BIE
Mechanism, dynamics, and biological existence of multistability in a large class of bursting neurons
Multistability, the coexistence of multiple attractors in a dynamical system,
is explored in bursting nerve cells. A modeling study is performed to show that
a large class of bursting systems, as defined by a shared topology when
represented as dynamical systems, is inherently suited to support
multistability. We derive the bifurcation structure and parametric trends
leading to multistability in these systems. Evidence for the existence of
multirhythmic behavior in neurons of the aquatic mollusc Aplysia californica
that is consistent with our proposed mechanism is presented. Although these
experimental results are preliminary, they indicate that single neurons may be
capable of dynamically storing information for longer time scales than
typically attributed to nonsynaptic mechanisms.Comment: 24 pages, 8 figure
A uniformly accurate scheme for the numerical integration of penalized Langevin dynamics
In molecular dynamics, penalized overdamped Langevin dynamics are used to
model the motion of a set of particles that follow constraints up to a
parameter . The most used schemes for simulating these dynamics
are the Euler integrator in and the constrained Euler
integrator. Both have weak order one of accuracy, but work properly only in
specific regimes depending on the size of the parameter . We
propose in this paper a new consistent method with an accuracy independent of
for solving penalized dynamics on a manifold of any dimension.
Moreover, this method converges to the constrained Euler scheme when
goes to zero. The numerical experiments confirm the theoretical
findings, in the context of weak convergence and for the invariant measure, on
a torus and on the orthogonal group in high dimension and high codimension.Comment: 27 page
Symmetric integrators with improved uniform error bounds and long-time conservations for the nonlinear Klein-Gordon equation in the nonrelativistic limit regime
In this paper, we are concerned with symmetric integrators for the nonlinear
relativistic Klein--Gordon (NRKG) equation with a dimensionless parameter
, which is inversely proportional to the speed of light.
The highly oscillatory property in time of this model corresponds to the
parameter and the equation has strong nonlinearity when \eps is
small. There two aspects bring significantly numerical burdens in designing
numerical methods. We propose and analyze a novel class of symmetric
integrators which is based on some formulation approaches to the problem,
Fourier pseudo-spectral method and exponential integrators. Two practical
integrators up to order four are constructed by using the proposed symmetric
property and stiff order conditions of implicit exponential integrators. The
convergence of the obtained integrators is rigorously studied, and it is shown
that the accuracy in time is improved to be \mathcal{O}(\varepsilon^{3}
\hh^2) and \mathcal{O}(\varepsilon^{4} \hh^4) for the time stepsize \hh.
The near energy conservation over long times is established for the multi-stage
integrators by using modulated Fourier expansions. These theoretical results
are achievable even if large stepsizes are utilized in the schemes. Numerical
results on a NRKG equation show that the proposed integrators have improved
uniform error bounds, excellent long time energy conservation and competitive
efficiency
A uniformly accurate scheme for the numerical integration of penalized Langevin dynamics
In molecular dynamics, penalized overdamped Langevin dynamics are used to model the motion of a set of particles that follow constraints up to a parameter ε. The most used schemes for simulating these dynamics are the Euler integrator in Rd and the constrained Euler integrator. Both have weak order one of accuracy, but work properly only in specific regimes depending on the size of the parameter ε. We propose in this paper a new consistent method with an accuracy independent of ε for solving penalized dynamics on a manifold of any dimension. Moreover, this method converges to the constrained Euler scheme when ε goes to zero. The numerical experiments confirm the theoretical findings, in the context of weak convergence and for the invariant measure, on a torus and on the orthogonal group in high dimension and high codimension.publishedVersio
Research and Education in Computational Science and Engineering
Over the past two decades the field of computational science and engineering
(CSE) has penetrated both basic and applied research in academia, industry, and
laboratories to advance discovery, optimize systems, support decision-makers,
and educate the scientific and engineering workforce. Informed by centuries of
theory and experiment, CSE performs computational experiments to answer
questions that neither theory nor experiment alone is equipped to answer. CSE
provides scientists and engineers of all persuasions with algorithmic
inventions and software systems that transcend disciplines and scales. Carried
on a wave of digital technology, CSE brings the power of parallelism to bear on
troves of data. Mathematics-based advanced computing has become a prevalent
means of discovery and innovation in essentially all areas of science,
engineering, technology, and society; and the CSE community is at the core of
this transformation. However, a combination of disruptive
developments---including the architectural complexity of extreme-scale
computing, the data revolution that engulfs the planet, and the specialization
required to follow the applications to new frontiers---is redefining the scope
and reach of the CSE endeavor. This report describes the rapid expansion of CSE
and the challenges to sustaining its bold advances. The report also presents
strategies and directions for CSE research and education for the next decade.Comment: Major revision, to appear in SIAM Revie
- …