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

A. J. Roberts

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

ResearcherUnclaimed source record
math.DSmath.NAmath.APmath.PRNumerical Analysisphysics.flu-dynchao-dynmath-phmath.MPnlin.CDphysics.comp-phBiological Physicscond-mat.softMathematical SoftwareTissues and Organs

Catalog footprint

What is connected

33works
15topics
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

33 published item(s)

preprint2022arXiv

Backwards theory supports modelling via invariant manifolds for non-autonomous dynamical systems

This article establishes the foundation for a new theory of invariant/integral manifolds for non-autonomous dynamical systems. Current rigorous support for dimensional reduction modelling of slow-fast systems is limited by the rare events in stochastic systems that may cause escape, and limited in many applications by the unbounded nature of PDE operators. To circumvent such limitations, we initiate developing a backward theory of invariant/integral manifolds that complements extant forward theory. Here, for deterministic non-autonomous ODE systems, we construct a conjugacy with a normal form system to establish the existence, emergence and exact construction of center manifolds in a finite domain for systems `arbitrarily close' to that specified. A benefit is that the constructed invariant manifolds are known to be exact for systems `close' to the one specified, and hence the only error is in determining how close over the domain of interest for any specific application. Built on the base developed here, planned future research should develop a theory for stochastic and/or PDE systems that is useful in a wide range of modelling applications.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Embed to rigorously and accurately homogenise quasi-periodic multi-scale heterogeneous PDEs, with computer algebra

For microscale heterogeneous PDEs, this article further develops novel theory and methodology for their macroscale mathematical/asymptotic homogenization. This article specifically encompasses the case of quasi-periodic heterogeneity with finite scale separation: no scale separation limit is required. Dynamical systems theory frames the homogenization as a slow manifold of the ensemble of all phase-shifts of the heterogeneity. Depending upon any perceived scale separation within the quasi-periodic heterogeneity, the homogenization may be done in either one step, or two sequential steps: the results are equivalent. The theory not only assures us of the existence and emergence of the homogenization, it also provides a practical systematic method to construct the homogenization to any specified order. For a class of heterogeneities, we show that the macroscale homogenization is potentially valid down to lengths which are just twice that of the microscale heterogeneity! This methodology provides a new rigorous and flexible approach to homogenization that potentially also provides correct initial and boundary conditions, treatment of forcing and control, and analysis of uncertainty.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Learning high-order spatial discretisations of PDEs with symmetry-preserving iterative algorithms

Common techniques for the spatial discretisation of PDEs on a macroscale grid include finite difference, finite elements and finite volume methods. Such methods typically impose assumed microscale structures on the subgrid fields, so without further tailored analysis are not suitable for systems with subgrid-scale heterogeneity or nonlinearities. We provide a new algebraic route to systematically approximate, in principle exactly, the macroscale closure of the spatially-discrete dynamics of a general class of heterogeneous non-autonomous reaction-advection-diffusion PDEs. This holistic discretisation approach, developed through rigorous theory and verified with computer algebra, systematically constructs discrete macroscale models through physics informed by the PDE out-of-equilibrium dynamics, thus relaxing many assumptions regarding the subgrid structure. The construction is analogous to recent gray-box machine learning techniques in that predictions are directed by iterative layers (as in neural networks), but informed by the subgrid physics (or 'data') as expressed in the PDEs. A major development of the holistic methodology, presented herein, is novel inter-element coupling between subgrid fields which preserve self-adjointness of the PDE after macroscale discretisation, thereby maintaining the spectral structure of the original system. This holistic methodology also encompasses homogenisation of microscale heterogeneous systems, as shown here with the canonical examples of heterogeneous 1D waves and diffusion.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Staggered grids for multidimensional multiscale modelling

Numerical schemes for wave-like systems with small dissipation are often inaccurate and unstable due to truncation errors and numerical roundoff errors. Hence, numerical simulations of wave-like systems lacking proper handling of these numerical issues often fail to represent the physical characteristics of wave phenomena. This challenge gets even more intricate for multiscale modelling, especially in multiple dimensions. When using the usual collocated grid, about two-thirds of the resolved wave modes are incorrect with significant dispersion. But, numerical schemes on staggered grids (with alternating variable arrangement) are significantly less dispersive and preserve much of the wave characteristics. Also, the group velocity of the energy propagation in the numerical waves on a staggered grid is in the correct direction, in contrast to the collocated grid. For high accuracy and to preserve much of the wave characteristics, this article extends the concept of staggered grids in full-domain modelling to multidimensional multiscale modelling. Specifically, this article develops 120 multiscale staggered grids and demonstrates their stability, accuracy, and wave-preserving characteristic for equation-free multiscale modelling of weakly damped linear waves. But most characteristics of the developed multiscale staggered grids must also hold in general for multiscale modelling of many complex spatio-temporal physical phenomena such as the general computational fluid dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A multiscale scheme accurately simulates macroscale shocks in an equation-free framework

Scientists and engineers often create accurate, trustworthy, computational simulation schemes - but all too often these are too computationally expensive to execute over the time or spatial domain of interest. The equation-free approach is to marry such trusted simulations to a framework for numerical macroscale reduction - the patch dynamics scheme. This article extends the patch scheme to scenarios in which the trusted simulation resolves abrupt state changes on the microscale that appear as shocks on the macroscale. Accurate simulation for problems in these scenarios requires extending the patch scheme by capturing the shock within a novel patch, and also modifying the patch coupling rules in the vicinity in order to maintain accuracy. With these two extensions to the patch scheme, straightforward arguments derive consistency conditions that match the usual order of accuracy for patch schemes. The new scheme is successfully tested on four archetypal problems. This technique will empower scientists and engineers to accurately and efficiently simulate, over large spatial domains, multiscale multiphysics systems that have rapid transition layers on the microscale.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A toolbox of Equation-Free functions in Matlab\Octave for efficient system level simulation

The `equation-free toolbox' empowers the computer-assisted analysis of complex, multiscale systems. Its aim is to enable you to immediately use microscopic simulators to perform macro-scale system level tasks and analysis, because micro-scale simulations are often the best available description of a system. The methodology bypasses the derivation of macroscopic evolution equations by computing the micro-scale simulator only over short bursts in time on small patches in space, with bursts and patches well-separated in time and space respectively. We introduce the suite of coded equation-free functions in an accessible way, link to more detailed descriptions, discuss their mathematical support, and introduce a novel and efficient algorithm for Projective Integration. Some facets of toolbox development of equation-free functions are then detailed. Download the toolbox functions (https://github.com/uoa1184615/EquationFreeGit) and use to empower efficient and accurate simulation in a wide range of your science and engineering problems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Equation-free patch scheme for efficient computational homogenisation via self-adjoint coupling

Equation-free macroscale modelling is a systematic and rigorous computational methodology for efficiently predicting the dynamics of a microscale system at a desired macroscale system level. In this scheme, the given microscale model is computed in small patches spread across the space-time domain, with patch coupling conditions bridging the unsimulated space. For accurate simulations, care must be taken in designing the patch coupling conditions. Here we construct novel coupling conditions which preserve translational invariance, rotational invariance, and self-adjoint symmetry, thus guaranteeing that conservation laws associated with these symmetries are preserved in the macroscale simulation. Spectral and algebraic analyses of the proposed scheme in both one and two dimensions reveal mechanisms for further improving the accuracy of the simulations. Consistency of the patch scheme's macroscale dynamics with the original microscale model is proved. This new self-adjoint patch scheme provides an efficient, flexible, and accurate computational homogenisation in a wide range of multiscale scenarios of interest to scientists and engineers.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Rigorous modelling of nonlocal interactions determines a macroscale advection-diffusion PDE

A slowly-varying or thin-layer multiscale assumption empowers macroscale understanding of many physical scenarios from dispersion in pipes and rivers, including beams, shells, and the modulation of nonlinear waves, to homogenisation of micro-structures. Here we begin a new exploration of the scenario where the given physics has non-local microscale interactions. We rigorously analyse the dynamics of a basic example of shear dispersion. Near each cross-section, the dynamics is expressed in the local moments of the microscale non-local effects. Centre manifold theory then supports the local modelling of the system's dynamics with coupling to neighbouring cross-sections as a non-autonomous forcing. The union over all cross-sections then provides powerful new support for the existence and emergence of a macroscale model advection-diffusion PDE global in the large, finite-sized, domain. The approach quantifies the accuracy of macroscale advection-diffusion approximations, and has the potential to open previously intractable multiscale issues to new insights.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Better buffers for patches in macroscale simulation of systems with microscale randomness

We consider one dimensional lattice diffusion model on a microscale grid with many discrete diffusivity values which repeat periodicially. Computer algebra explores how the dynamics of small coupled `patches' predict the slow emergent macroscale dynamics. We optimise the geometry and coupling of patches by comparing the macroscale predictions of the patch solutions with the macroscale solution on the infinite domain, which is derived for a general diffusivity period. The results indicate that patch dynamics is a viable method for numerical macroscale modelling of microscale systems with fine scale roughness. Moreover, the minimal error on the macroscale is generally obtained by coupling patches via `buffers' that are as large as half of each patch.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Slowly varying, macroscale models emerge from microscale dynamics over multiscale domains

Many physical systems are well described on domains which are relatively large in some directions but relatively thin in other directions. In this scenario we typically expect the system to have emergent structures that vary slowly over the large dimensions. For practical mathematical modelling of such systems we require efficient and accurate methodologies for reducing the dimension of the original system and extracting the emergent dynamics. Common mathematical approximations for determining the emergent dynamics often rely on self-consistency arguments or limits as the aspect ratio of the 'large' and 'thin' dimensions becomes unphysically infinite. Here we build on a new approach, previously establish for systems which are large in only one dimension, which analyses the dynamics at each cross-section of the domain with a rigorous multivariate Taylor series. Then centre manifold theory supports the local modelling of the system's emergent dynamics with coupling to neighbouring cross-sections treated as a non-autonomous forcing. The union over all cross-sections then provides powerful support for the existence and emergence of a centre manifold model global in the large finite domain. Quantitative error estimates are determined from the interactions between the cross-section coupling and both fast and slow dynamics. Two examples provide practical details of our methodology. The approach developed here may be used to quantify the accuracy of known approximations, to extend such approximations to mixed order modelling, and to open previously intractable modelling issues to new tools and insights.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Smooth subgrid fields underpin rigorous closure in spatial discretisation of reaction-advection-diffusion PDEs

Finite difference/element/volume methods of discretising PDEs impose a subgrid scale interpolation on the dynamics. In contrast, the holistic discretisation approach developed herein constructs a natural subgrid scale field adapted to the whole system out-of-equilibrium dynamics. Consequently, the macroscale discretisation is fully informed by the underlying microscale dynamics. We establish a new proof that in principle there exists an exact closure of the dynamics of a general class of reaction-advection-diffusion PDEs, and show how our approach constructs new systematic approximations to the in-principle closure starting from a simple, piecewise-linear, continuous approximation. Under inter-element coupling conditions that guarantee continuity of several field properties, the holistic discretisation possesses desirable properties such as a natural cubic spline first-order approximation to the field, and the self-adjointness of the diffusion operator under periodic, Dirichlet and Neumann macroscale boundary conditions. As a concrete example, we demonstrate the holistic discretisation procedure on the well-known Burgers' PDE, and compare the theoretical and numerical stability of the resulting discretisation to other approximations. The approach developed here promises to be able to systematically construct automatically good, macroscale discretisations to a wide range of PDEs, including wave PDEs.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Accuracy of patch dynamics with mesoscale temporal coupling for efficient exascale simulation

Massive parallelisation has lead to a dramatic increase in available computational power. However, data transfer speeds have failed to keep pace and are the major limiting factor in the development of exascale computing. New algorithms must be developed which minimise the transfer of data. Patch dynamics is a computational macroscale modelling scheme which provides a coarse macroscale solution of a problem defined on a fine microscale by dividing the domain into many nonoverlapping, coupled patches. Patch dynamics is readily adaptable to massive parallelisation as each processor can evaluate the dynamics on one, or a few, patches. However, patch coupling conditions interpolate across the unevaluated parts of the domain between patches, and are typically reevaluated at every microscale time step, thus requiring almost continuous data transfer. We propose a modified patch dynamics scheme which minimises data transfer by only reevaluating the patch coupling conditions at `mesoscale' time scales which are significantly larger than the microscale time of the microscale problem. We analyse the error arising from patch dynamics with mesoscale temporal coupling as a function of the mesoscale time interval, patch size, and ratio between the microscale and macroscale.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Ion mediated crosslink driven mucous swelling kinetics

We present an experimentally guided, multi-phasic, multi-species ionic gel model to compare and make qualitative predictions on the rheology of mucus of healthy individuals (Wild Type) versus those infected with Cystic Fibrosis. The mixture theory consists of the mucus (polymer phase) and water (solvent phase) as well as several different ions: H+, Na+ and Ca++. The model is linearized to study the hydration of spherically symmetric mucus gels and calibrated against the experimental data of mucus diffusivities. Near equilibrium, the linearized form of the equation describing the radial size of the gel, reduces to the well-known expression used in the kinetic theory of swelling hydrogels. Numerical studies reveal that the Donnan potential is the dominating mechanism driving the mucus swelling/deswelling transition. However, the altered swelling kinetics of the Cystic Fibrosis infected mucus is not merely governed by the hydroelectric composition of the swelling media, but also due to the altered movement of electrolytes as well as due to the defective properties of the mucin polymer network.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Macroscale boundary conditions for a non-linear heat exchanger

Multiscale modelling methodologies build macroscale models of materials with complicated fine microscale structure. We propose a methodology to derive boundary conditions for the macroscale model of a prototypical non-linear heat exchanger. The derived macroscale boundary conditions improve the accuracy of macroscale model. We verify the new boundary conditions by numerical methods. The techniques developed here can be adapted to a wide range of multiscale reaction-diffusion-advection systems.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Modelling suspended sediment in environmental turbulent fluids

Modelling sediment transport in environmental turbulent fluids is a challenge. This article develops a sound model of the lateral transport of suspended sediment in environmental fluid flows such as floods and tsunamis. The model is systematically derived from a 3D turbulence model based on the Smagorinski large eddy closure. Embedding the physical dynamics into a family of problems and analysing linear dynamics of the system, centre manifold theory indicates the existence of slow manifold parametrised by macroscale variables. Computer algebra then constructs the slow manifold in terms of fluid depth, depth-averaged lateral velocities, and suspended sediment concentration. The model includes the effects of sediment erosion, advection, dispersion, and also the interactions between the sediment and turbulent fluid flow. Vertical distributions of the velocity and concentration in steady flow agree with the established experimental data. Numerical simulations of the suspended sediment under large waves show that the developed model predicts physically reasonable phenomena.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Multiscale modelling couples patches of nonlinear wave-like simulations

The multiscale gap-tooth scheme is built from given microscale simulations of complicated physical processes to empower macroscale simulations. By coupling small patches of simulations over unsimulated physical gaps, large savings in computational time are possible. So far the gap-tooth scheme has been developed for dissipative systems, but wave systems are also of great interest. This article develops the gap-tooth scheme to the case of nonlinear microscale simulations of wave-like systems. Classic macroscale interpolation provides a generic coupling between patches that achieves arbitrarily high order consistency between the multiscale scheme and the underlying microscale dynamics. Eigen-analysis indicates that the resultant gap-tooth scheme empowers feasible computation of large scale simulations of wave-like dynamics with complicated underlying physics. As an pilot study, we implement numerical simulations of dam-breaking waves by the gap-tooth scheme. Comparison between a gap-tooth simulation, a microscale simulation over the whole domain, and some published experimental data on dam breaking, demonstrates that the gap-tooth scheme feasibly computes large scale wave-like dynamics with computational savings.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Multiscale modelling couples patches of two-layer thin fluid flow

The multiscale gap-tooth scheme uses a given microscale simulator of complicated physical processes to enable macroscale simulations by computing only only small sparse patches. This article develops the gap-tooth scheme to the case of nonlinear microscale simulations of thin fluid flow. The microscale simulator is derived by artificially assuming the fluid film flow having two artificial layers but no distinguishing physical feature. Centre manifold theory assures that there exists a slow manifold in the two-layer fluid film flow. Eigenvalue analysis confirms the stability of the microscale simulator. This article uses the gap-tooth scheme to simulate the two-layer fluid film flow. Coupling conditions are developed by approximating the values at the edges of patches by neighbouring macroscale values. Numerical eigenvalue analysis suggests that the gap-tooth scheme with the developed two-layer microscale simulator empowers feasible computation of large scale simulations of fluid film flows. We also implement numerical simulations of the fluid film flow by the gap-tooth scheme. Comparison between a gap-tooth simulation and a microscale simulation over the whole domain demonstrates that the gap-tooth scheme feasibly computes fluid film flow dynamics with computational savings.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Numerical integration of ordinary differential equations with rapidly oscillatory factors

We present a methodology for numerically integrating ordinary differential equations containing rapidly oscillatory terms. This challenge is distinct from that for differential equations which have rapidly oscillatory solutions: here the differential equation itself has the oscillatory terms. Our method generalises Filon quadrature for integrals, and is analogous to integral techniques designed to solve stochastic differential equations and, as such, is applicable to a wide variety of ordinary differential equations with rapidly oscillating factors. The proposed method flexibly achieves varying levels of accuracy depending upon the truncation of the expansion of certain integrals. Users will choose the level of truncation to suit the parameter regime of interest in their numerical integration.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Resolution of subgrid microscale interactions enhances the discretisation of nonautonomous partial differential equations

Coarse grained, macroscale, spatial discretisations of nonlinear nonautonomous partial differential\difference equations are given novel support by centre manifold theory. Dividing the physical domain into overlapping macroscale elements empowers the approach to resolve significant subgrid microscale structures and interactions between neighbouring elements. The crucial aspect of this approach is that centre manifold theory organises the resolution of the detailed subgrid microscale structure interacting via the nonlinear dynamics within and between neighbouring elements. The techniques and theory developed here may be applied to soundly discretise on a macroscale many dissipative nonautonomous partial differential\difference equations, such as the forced Burgers' equation, adopted here as an illustrative example.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Center manifolds for stochastic evolution equations

Stochastic invariant manifolds are crucial in modelling the dynamical behavior of dynamical systems under uncertainty. Under the assumption of exponential trichotomy, existence and smoothness of center manifolds for a class of stochastic evolution equations with linearly multiplicative noise are proved. The exponential attraction and approximation to center manifolds are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Computer algebra compares the stochastic superslow manifold of an averaged SPDE with that of the original slow-fast SPDE

The computer algebra routines documented here empower you to reproduce and check many of the details described by an article on large deviations for slow-fast stochastic systems [abs:1001.4826]. We consider a 'small' spatial domain with two coupled concentration fields, one governed by a 'slow' reaction-diffusion equation and one governed by a stochastic 'fast' linear equation. In the regime of a stochastic bifurcation, we derive two superslow models of the dynamics: the first is of the averaged model of the slow dynamics derived via large deviation principles; and the second is of the original fast-slow dynamics. Comparing the two superslow models validates the averaging in the large deviation principle in this parameter regime.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Modify the Improved Euler scheme to integrate stochastic differential equations

A practical and new Runge--Kutta numerical scheme for stochastic differential equations is explored. Numerical examples demonstrate the strong convergence of the method. The first order strong convergence is then proved using Ito integrals for both Ito and Stratonovich interpretations. As a straightforward modification of the deterministic Improved Euler/Heun method, the method is a good entry level scheme for stochastic differential equations, especially in conjunction with Higham's introduction [SIAM Review, 43:525--546, 2001].

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Averaging approximation to singularly perturbed nonlinear stochastic wave equations

An averaging method is applied to derive effective approximation to the following singularly perturbed nonlinear stochastic damped wave equation νu_{tt}+u_t=\D u+f(u)+ν^α\dot{W} on an open bounded domain $D\subset\R^n$\,, $1\leq n\leq 3$\,. Here $ν>0$ is a small parameter characterising the singular perturbation, and $ν^α$\,, $0\leq α\leq 1/2$\,, parametrises the strength of the noise. Some scaling transformations and the martingale representation theorem yield the following effective approximation for small $ν$, u_t=\D u+f(u)+ν^α\dot{W} to an error of $\ord{ν^α}$\,.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Computer algebra derives the slow manifold of patch or element dynamics on lattices in two dimensions

Developments in dynamical systems theory provides new support for the discretisation of \pde{}s and other microscale systems. Here we explore the methodology applied to the gap-tooth scheme in the equation-free approach of Kevrekidis in two spatial dimensions. The algebraic detail is enormous so we detail computer algebra procedures to handle the enormity. However, modelling the dynamics on 2D spatial patches appears to require a mixed numerical and algebraic approach that is detailed in this report. Being based upon the computation of residuals, the procedures here may be simply adapted to a wide class of reaction-diffusion equations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On the approximation for singularly perturbed stochastic wave equations

We explore the relation between fast waves, damping and imposed noise for different scalings by considering the singularly perturbed stochastic nonlinear wave equations νu_{tt}+u_t=\D u+f(u)+ν^α\dot{W} on a bounded spatial domain. An asymptotic approximation to the stochastic wave equation is constructed by a special transformation and splitting of $νu_{t}$. This splitting gives a clear description of the structure of $u$. The approximating model, for small $ν>0$\,, is a stochastic nonlinear heat equation for exponent $0\leqα<1$\,, and is a deterministic nonlinear wave equation for exponent $α>1$\,.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Self similarity and attraction in stochastic nonlinear reaction-diffusion systems

Similarity solutions play an important role in many fields of science: we consider here similarity in stochastic dynamics. Important issues are not only the existence of stochastic similarity, but also whether a similarity solution is dynamically attractive, and if it is, to what particular solution does the system evolve. By recasting a class of stochastic PDEs in a form to which stochastic centre manifold theory may be applied we resolve these issues in this class. For definiteness, a first example of self-similarity of the Burgers' equation driven by some stochastic forced is studied. Under suitable assumptions, a stationary solution is constructed which yields the existence of a stochastic self-similar solution for the stochastic Burgers' equation. Furthermore, the asymptotic convergence to the self-similar solution is proved. Second, in more general stochastic reaction-diffusion systems stochastic centre manifold theory provides a framework to construct the similarity solution, confirm its relevance, and determines the correct solution for any compact initial condition. Third, we argue that dynamically moving the spatial origin and dynamically stretching time improves the description of the stochastic similarity. Lastly, an application to an extremely simple model of turbulent mixing shows how anomalous fluctuations may arise in eddy diffusivities. The techniques and results we discuss should be applicable to a wide range of stochastic similarity problems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Large deviations for slow-fast stochastic partial differential equations

A large deviation principle is derived for stochastic partial differential equations with slow-fast components. The result shows that the rate function is exactly that of the averaged equation plus the fluctuating deviation which is a stochastic partial differential equation with small Gaussian perturbation. This also confirms the effectiveness of the approximation of the averaged equation plus the fluctuating deviation to the slow-fast stochastic partial differential equations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Macroscopic discrete modelling of stochastic reaction-diffusion equations on a periodic domain

Dynamical systems theory provides powerful methods to extract effective macroscopic dynamics from complex systems with slow modes and fast modes. Here we derive and theoretically support a macroscopic, spatially discrete, model for a class of stochastic reaction-diffusion partial differential equations with cubic nonlinearity. Dividing space into overlapping finite elements, a special coupling condition between neighbouring elements preserves the self-adjoint dynamics and controls interelement interactions. When the interelement coupling parameter is small, an averaging method and an asymptotic expansion of the slow modes show that the macroscopic discrete model will be a family of coupled stochastic ordinary differential equations which describe the evolution of the grid values. This modelling shows the importance of subgrid scale interaction between noise and spatial diffusion and provides a new rigourous approach to constructing semi-discrete approximations to stochastic reaction-diffusion partial differential equations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Resolve subgrid microscale interactions to discretise stochastic partial differential equations

Constructing discrete models of stochastic partial differential equations is very delicate. Stochastic centre manifold theory provides novel support for coarse grained, macroscale, spatial discretisations of nonlinear stochastic partial differential or difference equations such as the example of the stochastically forced Burgers' equation. Dividing the physical domain into finite length overlapping elements empowers the approach to resolve fully coupled dynamical interactions between neighbouring elements. The crucial aspect of this approach is that the underlying theory organises the resolution of the vast multitude of subgrid microscale noise processes interacting via the nonlinear dynamics within and between neighbouring elements. Noise processes with coarse structure across a finite element are the most significant noises for the discrete model. Their influence also diffuses away to weakly correlate the noise in the spatial discretisation. Nonlinear interactions have two further consequences: additive forcing generates multiplicative noise in the discretisation; and effectively new noise processes appear in the macroscale discretisation. The techniques and theory developed here may be applied to soundly discretise many dissipative stochastic partial differential and difference equations.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

An Accurate Lubrication Model of Contaminated Coating Flows

The levelling of short-wave irregularities on a thin film of fluid is primarily due to the action of surface tension. Surface tension gradients are often created by a number of different factor including evaporation, thermal gradients or deposition of surfactants. Lubrication theory, which ignores inertia terms in favour of viscous terms, produces a system of two nonlinear pde's for the unsteady flow of a thin viscous Newtonian fluid containing an insoluble surfactant. A complex model, which systematically includes all relevant effects to these pde's, is developed using centre manifold techniques. The benefits of using these techniques to develop accurate models are in their application. Subtle variations in parameters or assumptions are able to be catered for by including or deleting the relevant terms rather than having to redevelop these models.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Low-Dimensional Modelling of Dynamical Systems

Consider briefly the equations of fluid dynamics-they describe the enormous wealth of detail in all the interacting physical elements of a fluid flow-whereas in applications we want to deal with a description of just that which is interesting. In a wide variety of situations, simple approximate models are needed to perform practical simulations and make forecasts. I review the derivation, from a mathematical description of the detailed dynamics, of accurate, complete and useful low-dimensional models of the interesting dynamics in a system. The development of centre manifold theory and associated techniques puts this modelling process on a firm basis. As in Guckenheinmer & Holmes (1983,S2.5): "... these new methods will really be conventional perturbation style analyses interpreted geometrically..." But the geometric viewpoint of dynamical systems theory greatly enriches our approach by providing a rationale for also deriving correct initial conditions, forcing and boundary conditions for the models-all essential elements of a model.

0 citations0 reviewsNo code linkedOpen access
preprint1994arXiv

Initial Conditions for Models of Dynamical Systems

The long-time behaviour of many dynamical systems may be effectively predicted by a low-dimensional model that describes the evolution of a reduced set of variables. We consider the question of how to equip such a low-dimensional model with appropriate initial conditions, so that it faithfully reproduces the long-term behaviour of the original high-dimensional dynamical system. Our method involves putting the dynamical system into normal form, which not only generates the low-dimensional model, but also provides the correct initial conditions for the model. We illustrate the method with several examples. Keywords: normal form, isochrons, initialisation, centre manifold

0 citations0 reviewsNo code linkedOpen access