Source author record

Colin Cotter

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

7works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

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

7 published item(s)

preprint2021arXiv

A note on error analysis for a nonconforming discretisation of the tri-Helmholtz equation with singular data

We apply the nonconforming discretisation of Wu and Xu (2019) to the tri-Helmholtz equation on the plane where the source term is a functional evaluating the test function on a one-dimensional mesh-aligned embedded curve. We present error analysis for the convergence of the discretisation and linear convergence as a function of mesh size is recovered almost everywhere away from the embedded curve which aligns with classic regularity theory.

preprint2020arXiv

A Particle Filter for Stochastic Advection by Lie Transport (SALT): A case study for the damped and forced incompressible 2D Euler equation

In this work, we combine a stochastic model reduction with a particle filter augmented with tempering and jittering, and apply the combined algorithm to a damped and forced incompressible 2D Euler dynamics defined on a simply connected bounded domain. We show that using the combined algorithm, we are able to assimilate data from a reference system state (the ``truth") modelled by a highly resolved numerical solution of the flow that has roughly $3.1\times10^6$ degrees of freedom, into a stochastic system having two orders of magnitude less degrees of freedom, which is able to approximate the true state reasonably accurately for $5$ large scale eddy turnover times, using modest computational hardware. The model reduction is performed through the introduction of a stochastic advection by Lie transport (SALT) model as the signal on a coarser resolution. The SALT approach was introduced as a general theory using a geometric mechanics framework from Holm, Proc. Roy. Soc. A (2015). This work follows on the numerical implementation for SALT presented by Cotter et al, SIAM Multiscale Model. Sim. (2019) for the flow in consideration. The model reduction is substantial: The reduced SALT model has $4.9\times 10^4$ degrees of freedom. Results from reliability tests on the assimilated system are also presented.

preprint2020arXiv

Modelling uncertainty using stochastic transport noise in a 2-layer quasi-geostrophic model

The stochastic variational approach for geophysical fluid dynamics was introduced by Holm (Proc Roy Soc A, 2015) as a framework for deriving stochastic parameterisations for unresolved scales. This paper applies the variational stochastic parameterisation in a two-layer quasi-geostrophic model for a beta-plane channel flow configuration. We present a new method for estimating the stochastic forcing (used in the parameterisation) to approximate unresolved components using data from the high resolution deterministic simulation, and describe a procedure for computing physically-consistent initial conditions for the stochastic model. We also quantify uncertainty of coarse grid simulations relative to the fine grid ones in homogeneous (teamed with small-scale vortices) and heterogeneous (featuring horizontally elongated large-scale jets) flows, and analyse how the spread of stochastic solutions depends on different parameters of the model. The parameterisation is tested by comparing it with the true eddy-resolving solution that has reached some statistical equilibrium and the deterministic solution modelled on a low-resolution grid. The results show that the proposed parameterisation significantly depends on the resolution of the stochastic model and gives good ensemble performance for both homogeneous and heterogeneous flows, and the parameterisation lays solid foundations for data assimilation.

preprint2020arXiv

Space-time metamorphosis

We study the problem of registering images. The framework we use is metamorphosis and we construct a variational Eulerian space-time setting and pose the registration problem as an infinite-dimensional optimisation problem. The geodesic equations correspond to a system of advection and continuity equations and are solved analytically. Well-posedness of a primal conforming finite element method is established and its convergence is investigated numerically. This provides a discrete forward operator for the matching parameterized by a space-time velocity field. We propose a gradient descent method on this control variable and show several promising numerical results for this approach.

preprint2019arXiv

Data assimilation for a quasi-geostrophic model with circulation-preserving stochastic transport noise

This paper contains the latest installment of the authors' project on developing ensemble based data assimilation methodology for high dimensional fluid dynamics models. The algorithm presented here is a particle filter that combines model reduction, tempering, jittering, and nudging. The methodology is tested on a two-layer quasi-geostrophic model for a $β$-plane channel flow with $O(10^6)$ degrees of freedom out of which only a minute fraction are noisily observed. The model is reduced by following the stochastic variational approach for geophysical fluid dynamics introduced in Holm (Proc Roy Soc A, 2015) as a framework for deriving stochastic parametrisations for unresolved scales. The reduction is substantial: the computations are done only for $O(10^4)$ degrees of freedom. We introduce a stochastic time-stepping scheme for the two-layer model and prove its consistency in time. Then, we analyze the effect of the different procedures (tempering combined with jittering and nudging) on the performance of the data assimilation procedure using the reduced model, as well as how the dimension of the observational data (the number of "weather stations") and the data assimilation step affect the accuracy and uncertainty of the results.

preprint2016arXiv

Multilevel Ensemble Transform Particle Filtering

This paper extends the Multilevel Monte Carlo variance reduction technique to nonlinear filtering. In particular, Multilevel Monte Carlo is applied to a certain variant of the particle filter, the Ensemble Transform Particle Filter. A key aspect is the use of optimal transport methods to re-establish correlation between coarse and fine ensembles after resampling; this controls the variance of the estimator. Numerical examples present a proof of concept of the effectiveness of the proposed method, demonstrating significant computational cost reductions (relative to the single-level ETPF counterpart) in the propagation of ensembles.

preprint2011arXiv

Waltzing peakons and compacton pairs in a cross-coupled Camassa-Holm equation

We consider singular solutions of a system of two cross-coupled Camassa-Holm (CCCH) equations. This CCCH system admits peakon solutions, but it is not in the two-component CH integrable hierarchy. The system is a pair of coupled Hamiltonian partial differential equations for two types of solutions on the real line, each of which separately possesses exp(-|x|) peakon solutions with a discontinuity in the first derivative at the peak. However, there are no self-interactions, so each of the two types of peakon solutions moves only under the induced velocity of the other type. We analyse the `waltzing' solution behaviour of the cases with a single bound peakon pair (a peakon couple), as well as the over-taking collisions of peakon couples and the antisymmetric case of the head-on collision of a peakon couple and a peakon anti-couple. We then present numerical solutions of these collisions, which are inelastic because the waltzing peakon couples each possess an internal degree of freedom corresponding to their `tempo' -- that is, the period at which the two peakons of opposite type in the couple cycle around each other in phase space. Finally, we discuss compacton couple solutions of the cross-coupled Euler-Poincaré (CCEP) equations and illustrate the same types of collisions as for peakon couples, with triangular and parabolic compacton couples. We finish with a number of outstanding questions and challenges remaining for understanding couple dynamics of the CCCH and CCEP equations.