7,936 research outputs found
Efficient computation of high index Sturm-Liouville eigenvalues for problems in physics
Finding the eigenvalues of a Sturm-Liouville problem can be a computationally
challenging task, especially when a large set of eigenvalues is computed, or
just when particularly large eigenvalues are sought. This is a consequence of
the highly oscillatory behaviour of the solutions corresponding to high
eigenvalues, which forces a naive integrator to take increasingly smaller
steps. We will discuss some techniques that yield uniform approximation over
the whole eigenvalue spectrum and can take large steps even for high
eigenvalues. In particular, we will focus on methods based on coefficient
approximation which replace the coefficient functions of the Sturm-Liouville
problem by simpler approximations and then solve the approximating problem. The
use of (modified) Magnus or Neumann integrators allows to extend the
coefficient approximation idea to higher order methods
Recommended from our members
Comparison of Current Gravity Estimation and Determination Models
This paper will discuss the history of gravity estimation and determination models while analyzing methods that are in development. Some fundamental methods for calculating the gravity field include spherical harmonics solutions, local weighted interpolation, and global point mascon modeling (PMC). Recently, high accuracy measurements have become more accessible, and the requirements for high order geopotential modeling have become more stringent. Interest in irregular bodies, accurate models of the hydrological system, and on-board processing has demanded a comprehensive model that can quickly and accurately compute the geopotential with low memory costs. This trade study of current geopotential modeling techniques will reveal that each modeling technique has a unique use case. It is notable that the spherical harmonics model is relatively accurate but poses a cumbersome inversion problem. PMC and interpolation models, on the other hand, are computationally efficient, but require more research to become robust models with high levels of accuracy. Considerations of the trade study will suggest further research for the point mascon model. The PMC model should be improved through mascon refinement, direct solutions that stem from geodetic measurements, and further validation of the gravity gradient. Finally, the potential for each model to be implemented with parallel computation will be shown to lead to large improvements in computing time while reducing the memory cost for each technique.Aerospace Engineering and Engineering Mechanic
Boundary integral methods in high frequency scattering
In this article we review recent progress on the design, analysis and implementation of numerical-asymptotic boundary integral methods for the computation of frequency-domain acoustic scattering in a homogeneous unbounded medium by a bounded obstacle. The main aim of the methods is to allow computation of scattering at arbitrarily high frequency with finite computational resources
Analysis of parametric biological models with non-linear dynamics
In this paper we present recent results on parametric analysis of biological
models. The underlying method is based on the algorithms for computing
trajectory sets of hybrid systems with polynomial dynamics. The method is then
applied to two case studies of biological systems: one is a cardiac cell model
for studying the conditions for cardiac abnormalities, and the second is a
model of insect nest-site choice.Comment: In Proceedings HSB 2012, arXiv:1208.315
A Study of Different Modeling Choices For Simulating Platelets Within the Immersed Boundary Method
The Immersed Boundary (IB) method is a widely-used numerical methodology for
the simulation of fluid-structure interaction problems. The IB method utilizes
an Eulerian discretization for the fluid equations of motion while maintaining
a Lagrangian representation of structural objects. Operators are defined for
transmitting information (forces and velocities) between these two
representations. Most IB simulations represent their structures with
piecewise-linear approximations and utilize Hookean spring models to
approximate structural forces. Our specific motivation is the modeling of
platelets in hemodynamic flows. In this paper, we study two alternative
representations - radial basis functions (RBFs) and Fourier-based
(trigonometric polynomials and spherical harmonics) representations - for the
modeling of platelets in two and three dimensions within the IB framework, and
compare our results with the traditional piecewise-linear approximation
methodology. For different representative shapes, we examine the geometric
modeling errors (position and normal vectors), force computation errors, and
computational cost and provide an engineering trade-off strategy for when and
why one might select to employ these different representations.Comment: 33 pages, 17 figures, Accepted (in press) by APNU
- …