30,517 research outputs found
Accurate macroscale modelling of spatial dynamics in multiple dimensions
Developments in dynamical systems theory provides new support for the
macroscale modelling of pdes and other microscale systems such as Lattice
Boltzmann, Monte Carlo or Molecular Dynamics simulators. By systematically
resolving subgrid microscale dynamics the dynamical systems approach constructs
accurate closures of macroscale discretisations of the microscale system. Here
we specifically explore reaction-diffusion problems in two spatial dimensions
as a prototype of generic systems in multiple dimensions. Our approach unifies
into one the modelling of systems by a type of finite elements, and the
`equation free' macroscale modelling of microscale simulators efficiently
executing only on small patches of the spatial domain. Centre manifold theory
ensures that a closed model exist on the macroscale grid, is emergent, and is
systematically approximated. Dividing space either into overlapping finite
elements or into spatially separated small patches, the specially crafted
inter-element/patch coupling also ensures that the constructed discretisations
are consistent with the microscale system/PDE to as high an order as desired.
Computer algebra handles the considerable algebraic details as seen in the
specific application to the Ginzburg--Landau PDE. However, higher order models
in multiple dimensions require a mixed numerical and algebraic approach that is
also developed. The modelling here may be straightforwardly adapted to a wide
class of reaction-diffusion PDEs and lattice equations in multiple space
dimensions. When applied to patches of microscopic simulations our coupling
conditions promise efficient macroscale simulation.Comment: some figures with 3D interaction when viewed in Acrobat Reader. arXiv
admin note: substantial text overlap with arXiv:0904.085
ADAM: Analysis of Discrete Models of Biological Systems Using Computer Algebra
Background: Many biological systems are modeled qualitatively with discrete
models, such as probabilistic Boolean networks, logical models, Petri nets, and
agent-based models, with the goal to gain a better understanding of the system.
The computational complexity to analyze the complete dynamics of these models
grows exponentially in the number of variables, which impedes working with
complex models. Although there exist sophisticated algorithms to determine the
dynamics of discrete models, their implementations usually require
labor-intensive formatting of the model formulation, and they are oftentimes
not accessible to users without programming skills. Efficient analysis methods
are needed that are accessible to modelers and easy to use. Method: By
converting discrete models into algebraic models, tools from computational
algebra can be used to analyze their dynamics. Specifically, we propose a
method to identify attractors of a discrete model that is equivalent to solving
a system of polynomial equations, a long-studied problem in computer algebra.
Results: A method for efficiently identifying attractors, and the web-based
tool Analysis of Dynamic Algebraic Models (ADAM), which provides this and other
analysis methods for discrete models. ADAM converts several discrete model
types automatically into polynomial dynamical systems and analyzes their
dynamics using tools from computer algebra. Based on extensive experimentation
with both discrete models arising in systems biology and randomly generated
networks, we found that the algebraic algorithms presented in this manuscript
are fast for systems with the structure maintained by most biological systems,
namely sparseness, i.e., while the number of nodes in a biological network may
be quite large, each node is affected only by a small number of other nodes,
and robustness, i.e., small number of attractors
Computer algebra and operators
The symbolic computation of operator expansions is discussed. Some of the capabilities that prove useful when performing computer algebra computations involving operators are considered. These capabilities may be broadly divided into three areas: the algebraic manipulation of expressions from the algebra generated by operators; the algebraic manipulation of the actions of the operators upon other mathematical objects; and the development of appropriate normal forms and simplification algorithms for operators and their actions. Brief descriptions are given of the computer algebra computations that arise when working with various operators and their actions
Holistic finite differences accurately model the dynamics of the Kuramoto-Sivashinsky equation
We analyse the nonlinear Kuramoto-Sivashinsky equation to develop an accurate
finite difference approximation to its dynamics. The analysis is based upon
centre manifold theory so we are assured that the finite difference model
accurately models the dynamics and may be constructed systematically. The
theory is applied after dividing the physical domain into small elements by
introducing insulating internal boundaries which are later removed. The
Kuramoto-Sivashinsky equation is used as an example to show how holistic finite
differences may be applied to fourth order, nonlinear, spatio-temporal
dynamical systems. This novel centre manifold approach is holistic in the sense
that it treats the dynamical equations as a whole, not just as the sum of
separate terms
- …