BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Brynjulf Owren

Brynjulf Owren appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NANumerical AnalysisMachine Learningnlin.SImath.DSphysics.comp-phphysics.flu-dyn

Catalog footprint

What is connected

16works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2022arXiv

Detecting and determining preserved measures and integrals of birational maps

In this paper we use the method of discrete Darboux polynomials to calculate preserved measures and integrals of rational maps. The approach is based on the use of cofactors and Darboux polynomials and relies on the use of symbolic algebra tools. Given sufficient computing power, most, if not all, rational preserved integrals can be found (and even some non-rational ones). We show, in a number of examples, how it is possible to use this method to both determine and detect preserved measures and integrals of the considered rational maps. Many of the examples arise from the Kahan-Hirota-Kimura discretization of completely integrable systems of ordinary differential equations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Learning Hamiltonians of constrained mechanical systems

Recently, there has been an increasing interest in modelling and computation of physical systems with neural networks. Hamiltonian systems are an elegant and compact formalism in classical mechanics, where the dynamics is fully determined by one scalar function, the Hamiltonian. The solution trajectories are often constrained to evolve on a submanifold of a linear vector space. In this work, we propose new approaches for the accurate approximation of the Hamiltonian function of constrained mechanical systems given sample data information of their solutions. We focus on the importance of the preservation of the constraints in the learning strategy by using both explicit Lie group integrators and other classical schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Computational geometric methods for preferential clustering of particle suspensions

A geometric numerical method for simulating suspensions of spherical and non-spherical particles with Stokes drag is proposed. The method combines divergence-free matrix-valued radial basis function interpolation of the fluid velocity field with a splitting method integrator that preserves the sum of the Lyapunov spectrum while mimicking the centrifuge effect of the exact solution. We discuss how breaking the divergence-free condition in the interpolation step can erroneously affect how the volume of the particulate phase evolves under numerical methods. The methods are tested on suspensions of $10^4$ particles evolving in discrete cellular flow field. The results are that the proposed geometric methods generate more accurate and cost-effective particle distributions compared to conventional methods.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Structure preserving deep learning

Over the past few years, deep learning has risen to the foreground as a topic of massive interest, mainly as a result of successes obtained in solving large-scale image processing tasks. There are multiple challenging mathematical problems involved in applying deep learning: most deep learning methods require the solution of hard optimisation problems, and a good understanding of the tradeoff between computational effort, amount of data and model complexity is required to successfully design a deep learning approach for a given problem. A large amount of progress made in deep learning has been based on heuristic explorations, but there is a growing effort to mathematically understand the structure in existing deep learning methods and to systematically design new deep learning methods to preserve certain types of structure in deep learning. In this article, we review a number of these directions: some deep neural networks can be understood as discretisations of dynamical systems, neural networks can be designed to have desirable properties such as invertibility or group equivariance, and new algorithmic frameworks based on conformal Hamiltonian systems and Riemannian manifolds to solve the optimisation problems have been proposed. We conclude our review of each of these topics by discussing some open problems that we consider to be interesting directions for future research.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Lie group integrators

In this survey we discuss a wide variety of aspects related to Lie group integrators. These numerical integration schemes for differential equations on manifolds have been studied in a general and systematic manner since the 1990s and the activity has since then branched out in several different subareas, focussing both on theoretical and practical issues. From two alternative setups, using either frames or Lie group actions on a manifold, we here introduce the most important classes of schemes used to integrate nonlinear ordinary differential equations on Lie groups and manifolds. We describe a number of different applications where there is a natural action by a Lie group on a manifold such that our integrators can be implemented. An issue which is not well understood is the role of isotropy and how it affects the behaviour of the numerical methods. The order theory of numerical Lie group integrators has become an advanced subtopic in its own right, and here we give a brief introduction on a somewhat elementary level. Finally, we shall discuss Lie group integrators having the property that they preserve a symplectic structure or a first integral.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Discretization of polynomial vector fields by polarization

A novel integration method for quadratic vector fields was introduced by Kahan in 1993. Subsequently, it was shown that Kahan's method preserves a (modified) measure and energy when applied to quadratic Hamiltonian vector fields. Here we generalize Kahan's method to cubic resp. higher degree polynomial vector fields and show that the resulting discretization also preserves modified versions of the measure and energy when applied to cubic resp. higher degree polynomial Hamiltonian vector fields.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Geometric integration of non-autonomous Hamiltonian problems

Symplectic integration of autonomous Hamiltonian systems is a well-known field of study in geometric numerical integration, but for non-autonomous systems the situation is less clear, since symplectic structure requires an even number of dimensions. We show that one possible extension of symplectic methods in the autonomous setting to the non-autonomous setting is obtained by using canonical transformations. Many existing methods fit into this framework. We also perform experiments which indicate that for exponential integrators, the canonical and symmetric properties are important for good long time behaviour. In particular, the theoretical and numerical results support the well documented fact from the literature that exponential integrators for non-autonomous linear problems have superior accuracy compared to general ODE schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Geometric properties of Kahan's method

We show that Kahan's discretization of quadratic vector fields is equivalent to a Runge--Kutta method. When the vector field is Hamiltonian on either a symplectic vector space or a Poisson vector space with constant Poisson structure, the map determined by this discretization has a conserved modified Hamiltonian and an invariant measure, a combination previously unknown amongst Runge--Kutta methods applied to nonlinear vector fields. This produces large classes of integrable rational mappings in two and three dimensions, explaining some of the integrable cases that were previously known.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The minimal stage, energy preserving Runge-Kutta method for polynomial Hamiltonian systems is the Averaged Vector Field method

No Runge-Kutta method can be energy preserving for all Hamiltonian systems. But for problems in which the Hamiltonian is a polynomial, the Averaged Vector Field (AVF) method can be interpreted as a Runge-Kutta method whose weights $b_i$ and abscissae $c_i$ represent a quadrature rule of degree at least that of the Hamiltonian. We prove that when the number of stages is minimal, the Runge-Kutta scheme must in fact be identical to the AVF scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A general framework for deriving integral preserving numerical methods for PDEs

A general procedure for constructing conservative numerical integrators for time dependent partial differential equations is presented. In particular, linearly implicit methods preserving a time discretised version of the invariant is developed for systems of partial differential equations with polynomial nonlinearities. The framework is rather general and allows for an arbitrary number of dependent and independent variables with derivatives of any order. It is proved formally that second order convergence is obtained. The procedure is applied to a test case and numerical experiments are provided.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Preserving multiple first integrals by discrete gradients

We consider systems of ordinary differential equations with known first integrals. The notion of a discrete tangent space is introduced as the orthogonal complement of an arbitrary set of discrete gradients. Integrators which exactly conserve all the first integrals simultaneously are then defined. In both cases we start from an arbitrary method of a prescribed order (say, a Runge-Kutta scheme) and modify it using two approaches: one based on projection and one based one local coordinates. The methods are tested on the Kepler problem.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Plane wave stability of some conservative schemes for the cubic Schrödinger equation

The plane wave stability properties of the conservative schemes of Besse and Fei et al. for the cubic Schrödinger equation are analysed. Although the two methods possess many of the same conservation properties, we show that their stability behaviour is very different. An energy preserving generalisation of the Fei method with improved stability is presented.

0 citations0 reviewsNo code linkedOpen access