Source author record

K. J. H. Law

K. J. H. Law 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

18works
12topics
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

18 published item(s)

preprint2015arXiv

Controlling Unpredictability with Observations in the Partially Observed Lorenz '96 Model

In the context of filtering chaotic dynamical systems it is well-known that partial observations, if sufficiently informative, can be used to control the inherent uncertainty due to chaos. The purpose of this paper is to investigate, both theoretically and numerically, conditions on the observations of chaotic systems under which they can be accurately filtered. In particular, we highlight the advantage of adaptive observation operators over fixed ones. The Lorenz '96 model is used to exemplify our findings. We consider discrete-time and continuous-time observations in our theoretical developments. We prove that, for fixed observation operator, the 3DVAR filter can recover the system state within a neighbourhood determined by the size of the observational noise. It is required that a sufficiently large proportion of the state vector is observed, and an explicit form for such sufficient fixed observation operator is given. Numerical experiments, where the data is incorporated by use of the 3DVAR and extended Kalman filters, suggest that less informative fixed operators than given by our theory can still lead to accurate signal reconstruction. Adaptive observation operators are then studied numerically; we show that, for carefully chosen adaptive observation operators, the proportion of the state vector that needs to be observed is drastically smaller than with a fixed observation operator. Indeed, we show that the number of state coordinates that need to be observed may even be significantly smaller than the total number of positive Lyapunov exponents of the underlying system.

preprint2014arXiv

Dynamic and Energetic Stabilization of Persistent Currents in Bose-Einstein Condensates

We study conditions under which vortices in a highly oblate harmonically trapped Bose-Einstein condensate (BEC) can be stabilized due to pinning by a blue-detuned Gaussian laser beam, with particular emphasis on the potentially destabilizing effects of laser beam positioning within the BEC. Our approach involves theoretical and numerical exploration of dynamically and energetically stable pinning of vortices with winding number up to $S=6$, in correspondence with experimental observations. Stable pinning is quantified theoretically via Bogoliubov-de Gennes excitation spectrum computations and confirmed via direct numerical simulations for a range of conditions similar to those of experimental observations. The theoretical and numerical results indicate that the pinned winding number, or equivalently the winding number of the superfluid current about the laser beam, decays as a laser beam of fixed intensity moves away from the BEC center. Our theoretical analysis helps explain previous experimental observations, and helps define limits of stable vortex pinning for future experiments involving vortex manipulation by laser beams.

preprint2013arXiv

Accuracy and Stability of Filters for Dissipative PDEs

Data assimilation methodologies are designed to incorporate noisy observations of a physical system into an underlying model in order to infer the properties of the state of the system. Filters refer to a class of data assimilation algorithms designed to update the estimation of the state as data is acquired sequentially. For linear problems subject to Gaussian noise filtering can be performed exactly using the Kalman filter. For nonlinear systems it can be approximated in a systematic way by particle filters. However in high dimensions these particle filtering methods can break down. Hence, for the large nonlinear systems arising in applications such as oceanography and weather forecasting, various ad hoc filters are used, based on Gaussian approximations. In this work, we study the accuracy and stability of these ad hoc filters in the context of the 2D incompressible Navier-Stokes equation. The ideas readily generalize to a range of dissipative partial differential equations (PDEs). By working in this infinite dimensional setting we provide an analysis which is useful for the understanding of high dimensional filtering, and is robust to mesh-refinement. We describe theoretical results showing that, in the small observational noise limit, the filters can be tuned to perform accurately in tracking the signal itself (filter accuracy), provided the system is observed in a sufficiently large low dimensional space; roughly speaking this space should be large enough to contain the unstable modes of the linearized dynamics. The tuning corresponds to what is known as variance inflation in the applied literature. Numerical results are given which illustrate the theory. The positive results herein concerning filter stability complement recent numerical studies which demonstrate that the ad hoc filters can perform poorly in reproducing statistical variation about the true signal.

preprint2013arXiv

Analysis of the 3DVAR Filter for the Partially Observed Lorenz '63 Model

The problem of effectively combining data with a mathematical model constitutes a major challenge in applied mathematics. It is particular challenging for high-dimensional dynamical systems where data is received sequentially in time and the objective is to estimate the system state in an on-line fashion; this situation arises, for example, in weather forecasting. The sequential particle filter is then impractical and ad hoc filters, which employ some form of Gaussian approximation, are widely used. Prototypical of these ad hoc filters is the 3DVAR method. The goal of this paper is to analyze the 3DVAR method, using the Lorenz '63 model to exemplify the key ideas. The situation where the data is partial and noisy is studied, and both discrete time and continuous time data streams are considered. The theory demonstrates how the widely used technique of variance inflation acts to stabilize the filter, and hence leads to asymptotic accuracy.

preprint2013arXiv

Well-Posedness And Accuracy Of The Ensemble Kalman Filter In Discrete And Continuous Time

The ensemble Kalman filter (EnKF) is a method for combining a dynamical model with data in a sequential fashion. Despite its widespread use, there has been little analysis of its theoretical properties. Many of the algorithmic innovations associated with the filter, which are required to make a useable algorithm in practice, are derived in an ad hoc fashion. The aim of this paper is to initiate the development of a systematic analysis of the EnKF, in particular to do so in the small ensemble size limit. The perspective is to view the method as a state estimator, and not as an algorithm which approximates the true filtering distribution. The perturbed observation version of the algorithm is studied, without and with variance inflation. Without variance inflation well-posedness of the filter is established; with variance inflation accuracy of the filter, with resepct to the true signal underlying the data, is established. The algorithm is considered in discrete time, and also for a continuous time limit arising when observations are frequent and subject to large noise. The underlying dynamical model, and assumptions about it, is sufficiently general to include the Lorenz '63 and '96 models, together with the incompressible Navier-Stokes equation on a two-dimensional torus. The analysis is limited to the case of complete observation of the signal with additive white noise. Numerical results are presented for the Navier-Stokes equation on a two-dimensional torus for both complete and partial observations of the signal with additive white noise.

preprint2012arXiv

Accuracy and Stability of The Continuous-Time 3DVAR Filter for The Navier-Stokes Equation

The 3DVAR filter is prototypical of methods used to combine observed data with a dynamical system, online, in order to improve estimation of the state of the system. Such methods are used for high dimensional data assimilation problems, such as those arising in weather forecasting. To gain understanding of filters in applications such as these, it is hence of interest to study their behaviour when applied to infinite dimensional dynamical systems. This motivates study of the problem of accuracy and stability of 3DVAR filters for the Navier-Stokes equation. We work in the limit of high frequency observations and derive continuous time filters. This leads to a stochastic partial differential equation (SPDE) for state estimation, in the form of a damped-driven Navier-Stokes equation, with mean-reversion to the signal, and spatially-correlated time-white noise. Both forward and pullback accuracy and stability results are proved for this SPDE, showing in particular that when enough low Fourier modes are observed, and when the model uncertainty is larger than the data uncertainty in these modes (variance inflation), then the filter can lock on to a small neighbourhood of the true signal, recovering from order one initial error, if the error in the observations modes is small. Numerical examples are given to illustrate the theory.

preprint2012arXiv

Characteristics of Two-Dimensional Quantum Turbulence in a Compressible Superfluid

Under suitable forcing a fluid exhibits turbulence, with characteristics strongly affected by the fluid's confining geometry. Here we study two-dimensional quantum turbulence in a highly oblate Bose-Einstein condensate in an annular trap. As a compressible quantum fluid, this system affords a rich phenomenology, allowing coupling between vortex and acoustic energy. Small-scale stirring generates an experimentally observed disordered vortex distribution that evolves into large-scale flow in the form of a persistent current. Numerical simulation of the experiment reveals additional characteristics of two-dimensional quantum turbulence: spontaneous clustering of same-circulation vortices, and an incompressible energy spectrum with $k^{-5/3}$ dependence for low wavenumbers $k$ and $k^{-3}$ dependence for high $k$.

preprint2012arXiv

Evaluating Data Assimilation Algorithms

Data assimilation leads naturally to a Bayesian formulation in which the posterior probability distribution of the system state, given the observations, plays a central conceptual role. The aim of this paper is to use this Bayesian posterior probability distribution as a gold standard against which to evaluate various commonly used data assimilation algorithms. A key aspect of geophysical data assimilation is the high dimensionality and low predictability of the computational model. With this in mind, yet with the goal of allowing an explicit and accurate computation of the posterior distribution, we study the 2D Navier-Stokes equations in a periodic geometry. We compute the posterior probability distribution by state-of-the-art statistical sampling techniques. The commonly used algorithms that we evaluate against this accurate gold standard, as quantified by comparing the relative error in reproducing its moments, are 4DVAR and a variety of sequential filtering approximations based on 3DVAR and on extended and ensemble Kalman filters. The primary conclusions are that: (i) with appropriate parameter choices, approximate filters can perform well in reproducing the mean of the desired probability distribution; (ii) however they typically perform poorly when attempting to reproduce the covariance; (iii) this poor performance is compounded by the need to modify the covariance, in order to induce stability. Thus, whilst filters can be a useful tool in predicting mean behavior, they should be viewed with caution as predictors of uncertainty. These conclusions are intrinsic to the algorithms and will not change if the model complexity is increased, for example by employing a smaller viscosity, or by using a detailed NWP model.

preprint2011arXiv

Stability of Filters for the Navier-Stokes Equation

Data assimilation methodologies are designed to incorporate noisy observations of a physical system into an underlying model in order to infer the properties of the state of the system. Filters refer to a class of data assimilation algorithms designed to update the estimation of the state in a on-line fashion, as data is acquired sequentially. For linear problems subject to Gaussian noise filtering can be performed exactly using the Kalman filter. For nonlinear systems it can be approximated in a systematic way by particle filters. However in high dimensions these particle filtering methods can break down. Hence, for the large nonlinear systems arising in applications such as weather forecasting, various ad hoc filters are used, mostly based on making Gaussian approximations. The purpose of this work is to study the properties of these ad hoc filters, working in the context of the 2D incompressible Navier-Stokes equation. By working in this infinite dimensional setting we provide an analysis which is useful for understanding high dimensional filtering, and is robust to mesh-refinement. We describe theoretical results showing that, in the small observational noise limit, the filters can be tuned to accurately track the signal itself (filter stability), provided the system is observed in a sufficiently large low dimensional space; roughly speaking this space should be large enough to contain the unstable modes of the linearized dynamics. Numerical results are given which illustrate the theory. In a simplified scenario we also derive, and study numerically, a stochastic PDE which determines filter stability in the limit of frequent observations, subject to large observational noise. The positive results herein concerning filter stability complement recent numerical studies which demonstrate that the ad hoc filters perform poorly in reproducing statistical variation about the true signal.

preprint2010arXiv

Stable structures with high topological charge in nonlinear photonic quasicrystals

Stable vortices with topological charge of 3 and 4 are examined numerically and analytically in photonic quasicrystals created by interference of 5 as well as 8 beams, in the cases of cubic as well as saturable nonlinearities. These structures are experimentally realizable, including a prototypical example of a stable charge 4 vortex. Direct numerical simulations corroborate the analytical and numerical linear stability analysis predictions.

preprint2010arXiv

Stable Vortex-Bright Soliton Structures in Two-Component Bose Einstein Condensates

We report the numerical realization and demonstration of robustness of certain 2-component structures in Bose-Einstein Condensates in 2 and 3 spatial dimensions with non-trivial topological charge in one of the components. In particular, we identify a stable symbiotic state in which a higher-dimensional bright soliton exists even in a homogeneous setting with defocusing interactions, as a result of the effective potential created by a stable vortex in the other component. The resulting vortex-bright solitary waves, which naturally generalize the recently experimentally observed dark-bright solitons, are examined both in the homogeneous medium and in the presence of parabolic and periodic external confinement and are found to be very robust.

preprint2009arXiv

A two-dimensional paradigm for symmetry breaking: the nonlinear Schrodinger equation with a four-well potential

We study the existence and stability of localized states in the two-dimensional (2D) nonlinear Schrodinger (NLS)/Gross-Pitaevskii equation with a symmetric four-well potential. Using a fourmode approximation, we are able to trace the parametric evolution of the trapped stationary modes, starting from the corresponding linear limits, and thus derive the complete bifurcation diagram for the families of these stationary modes. The predictions based on the four-mode decomposition are found to be in good agreement with the numerical results obtained from the NLS equation. Actually, the stability properties coincide with those suggested by the corresponding discrete model in the large-amplitude limit. The dynamics of the unstable modes is explored by means of direct simulations. Finally, while we present the full analysis for the case of the focusing nonlinearity, the bifurcation diagram for the defocusing case is briefly considered too.

preprint2009arXiv

Existence and stability of multisite breathers in honeycomb and hexagonal lattices

We study the existence and stability of multisite discrete breathers in two prototypical non-square Klein-Gordon lattices, namely a honeycomb and a hexagonal one. In the honeycomb case we consider six-site configurations and find that for soft potential and positive coupling the out-of-phase breather configuration and the charge-two vortex breather are linearly stable, while the in-phase and charge-one vortex states are unstable. In the hexagonal lattice, we first consider three-site configurations. In the case of soft potential and positive coupling, the in-phase configuration is unstable and the charge-one vortex is linearly stable. The out-of-phase configuration here is found to always be linearly unstable. We then turn to six-site configurations in the hexagonal lattice. The stability results in this case are the same as in the six-site configurations in the honeycomb lattice. For all configurations in both lattices, the stability results are reversed in the setting of either hard potential or negative coupling. The study is complemented by numerical simulations which are in very good agreement with the theoretical predictions. Since neither the form of the on-site potential nor the sign of the coupling parameter involved have been prescribed, this description can accommodate inverse-dispersive systems (e.g., supporting backward waves) such as transverse dust-lattice oscillations in dusty plasma (Debye) crystals or analogous modes in molecular chains.

preprint2009arXiv

Geometric stabilization of extended S=2 vortices in two-dimensional photonic lattices: theoretical analysis, numerical computation and experimental results

In this work, we focus our studies on the subject of nonlinear discrete self-trapping of S=2 (doubly-charged) vortices in two-dimensional photonic lattices, including theoretical analysis, numerical computation and experimental demonstration. We revisit earlier findings about S=2 vortices with a discrete model, and find that S=2 vortices extended over eight lattice sites can indeed be stable (or only weakly unstable) under certain conditions, not only for the cubic nonlinearity previously used, but also for a saturable nonlinearity more relevant to our experiment with a biased photorefractive nonlinear crystal. We then use the discrete analysis as a guide towards numerically identifying stable (and unstable) vortex solutions in a more realistic continuum model with a periodic potential. Finally, we present our experimental observation of such geometrically extended S=2 vortex solitons in optically induced lattices under both self-focusing and self-defocusing nonlinearities, and show clearly that the S=2 vortex singularities are preserved during nonlinear propagation.

preprint2009arXiv

Localized structures in Kagome lattices

We investigate the existence and stability of gap vortices and multi-pole gap solitons in a Kagome lattice with a defocusing nonlinearity both in a discrete case and in a continuum one with periodic external modulation. In particular, predictions are made based on expansion around a simple and analytically tractable anti-continuum (zero coupling) limit. These predictions are then confirmed for a continuum model of an optically-induced Kagome lattice in a photorefractive crystal obtained by a continuous transformation of a honeycomb lattice.

preprint2009arXiv

Manipulation of Vortices by Localized Impurities in Bose-Einstein Condensates

We consider the manipulation of Bose-Einstein condensate vortices by optical potentials generated by focused laser beams. It is shown that for appropriate choices of the laser strength and width it is possible to successfully transport vortices to various positions inside the trap confining the condensate atoms. Furthermore, the full bifurcation structure of possible stationary single-charge vortex solutions in a harmonic potential with this type of impurity is elucidated. The case when a moving vortex is captured by a stationary laser beam is also studied, as well as the possibility of dragging the vortex by means of periodic optical lattices.

preprint2008arXiv

Structure and stability of two-dimensional Bose-Einstein condensates under both harmonic and lattice confinement

In this work, we study pancake-shaped Bose-Einstein condensates confined by both a cylindrically symmetric harmonic potential and an optical lattice with equal periodicity in two orthogonal directions. We first identify the spectrum of the underlying two-dimensional linear problem through multiple-scale techniques. Then, we use the results obtained in the linear limit as a starting point for a nonlinear existence and stability analysis of the lowest energy states, emanating from the linear ones, in the nonlinear problem. Two-parameter continuations of these states are performed for increasing nonlinearity and optical lattice strengths, and their instabilities and temporal evolution are investigated. It is found that the ground state as well as one of the excited states are either stable or weakly unstable for both attractive and repulsive interatomic interactions.