Showing 1–50 of 50 results for author: Kalliadasis, S

Forecasting with an N-dimensional Langevin Equation and a Neural-Ordinary Differential Equation

Authors: Antonio Malpica-Morales, Miguel A. Duran-Olivencia, Serafim Kalliadasis

Abstract: Accurate prediction of electricity day-ahead prices is essential in competitive electricity markets. Although stationary electricity-price forecasting techniques have received considerable attention, research on non-stationary methods is comparatively scarce, despite the common prevalence of non-stationary features in electricity markets. Specifically, existing non-stationary techniques will often… ▽ More Accurate prediction of electricity day-ahead prices is essential in competitive electricity markets. Although stationary electricity-price forecasting techniques have received considerable attention, research on non-stationary methods is comparatively scarce, despite the common prevalence of non-stationary features in electricity markets. Specifically, existing non-stationary techniques will often aim to address individual non-stationary features in isolation, leaving aside the exploration of concurrent multiple non-stationary effects. Our overarching objective here is the formulation of a framework to systematically model and forecast non-stationary electricity-price time series, encompassing the broader scope of non-stationary behavior. For this purpose we develop a data-driven model that combines an N-dimensional Langevin equation (LE) with a neural-ordinary differential equation (NODE). The LE captures fine-grained details of the electricity-price behavior in stationary regimes but is inadequate for non-stationary conditions. To overcome this inherent limitation, we adopt a NODE approach to learn, and at the same time predict, the difference between the actual electricity-price time series and the simulated price trajectories generated by the LE. By learning this difference, the NODE reconstructs the non-stationary components of the time series that the LE is not able to capture. We exemplify the effectiveness of our framework using the Spanish electricity day-ahead market as a prototypical case study. Our findings reveal that the NODE nicely complements the LE, providing a comprehensive strategy to tackle both stationary and non-stationary electricity-price behavior. The framework's dependability and robustness is demonstrated through different non-stationary scenarios by comparing it against a range of basic naive methods. △ Less

Submitted 12 May, 2024; originally announced May 2024.

Comments: 26 pages, 7 figures

Journal ref: Chaos, 34, 043105 (2024)

Physics-informed Bayesian inference of external potentials in classical density-functional theory

Authors: Antonio Malpica-Morales, Peter Yatsyshin, Miguel A. Duran-Olivencia, Serafim Kalliadasis

Abstract: The swift progression of machine learning (ML) has not gone unnoticed in the realm of statistical mechanics. ML techniques have attracted attention by the classical density-functional theory (DFT) community, as they enable discovery of free-energy functionals to determine the equilibrium-density profile of a many-particle system. Within DFT, the external potential accounts for the interaction of t… ▽ More The swift progression of machine learning (ML) has not gone unnoticed in the realm of statistical mechanics. ML techniques have attracted attention by the classical density-functional theory (DFT) community, as they enable discovery of free-energy functionals to determine the equilibrium-density profile of a many-particle system. Within DFT, the external potential accounts for the interaction of the many-particle system with an external field, thus, affecting the density distribution. In this context, we introduce a statistical-learning framework to infer the external potential exerted on a many-particle system. We combine a Bayesian inference approach with the classical DFT apparatus to reconstruct the external potential, yielding a probabilistic description of the external potential functional form with inherent uncertainty quantification. Our framework is exemplified with a grand-canonical one-dimensional particle ensemble with excluded volume interactions in a confined geometry. The required training dataset is generated using a Monte Carlo (MC) simulation where the external potential is applied to the grand-canonical ensemble. The resulting particle coordinates from the MC simulation are fed into the learning framework to uncover the external potential. This eventually allows us to compute the equilibrium density profile of the system by using the tools of DFT. Our approach benchmarks the inferred density against the exact one calculated through the DFT formulation with the true external potential. The proposed Bayesian procedure accurately infers the external potential and the density profile. We also highlight the external-potential uncertainty quantification conditioned on the amount of available simulated data. The seemingly simple case study introduced in this work might serve as a prototype for studying a wide variety of applications, including adsorption and capillarity. △ Less

Submitted 14 September, 2023; v1 submitted 13 September, 2023; originally announced September 2023.

Journal ref: J. Chem. Phys. 159, 104109 (2023)

Unconditional bound-preserving and energy-dissipating finite-volume schemes for the Cahn-Hilliard equation

Authors: Rafael Bailo, José A. Carrillo, Serafim Kalliadasis, Sergio P. Perez

Abstract: We propose finite-volume schemes for the Cahn-Hilliard equation which unconditionally and discretely preserve the boundedness of the phase field and the dissipation of the free energy. Our numerical framework is applicable to a variety of free-energy potentials, including Ginzburg-Landau and Flory-Huggins, to general wetting boundary conditions, and to degenerate mobilities. Its central thrust is… ▽ More We propose finite-volume schemes for the Cahn-Hilliard equation which unconditionally and discretely preserve the boundedness of the phase field and the dissipation of the free energy. Our numerical framework is applicable to a variety of free-energy potentials, including Ginzburg-Landau and Flory-Huggins, to general wetting boundary conditions, and to degenerate mobilities. Its central thrust is the upwind methodology, which we combine with a semi-implicit formulation for the free-energy terms based on the classical convex-splitting approach. The extension of the schemes to an arbitrary number of dimensions is straightforward thanks to their dimensionally split nature, which allows to efficiently solve higher-dimensional problems with a simple parallelisation. The numerical schemes are validated and tested through a variety of examples, in different dimensions, and with various contact angles between droplets and substrates. △ Less

Submitted 17 February, 2023; v1 submitted 11 May, 2021; originally announced May 2021.

MSC Class: 65M08; 35Q92; 35Q35; 35Q70

Physics-constrained Bayesian inference of state functions in classical density-functional theory

Authors: Peter Yatsyshin, Serafim Kalliadasis, Andrew B. Duncan

Abstract: We develop a novel data-driven approach to the inverse problem of classical statistical mechanics: given experimental data on the collective motion of a classical many-body system, how does one characterise the free energy landscape of that system? By combining non-parametric Bayesian inference with physically-motivated constraints, we develop an efficient learning algorithm which automates the co… ▽ More We develop a novel data-driven approach to the inverse problem of classical statistical mechanics: given experimental data on the collective motion of a classical many-body system, how does one characterise the free energy landscape of that system? By combining non-parametric Bayesian inference with physically-motivated constraints, we develop an efficient learning algorithm which automates the construction of approximate free energy functionals. In contrast to optimisation-based machine learning approaches, which seek to minimise a cost function, the central idea of the proposed Bayesian inference is to propagate a set of prior assumptions through the model, derived from physical principles. The experimental data is used to probabilistically weigh the possible model predictions. This naturally leads to humanly interpretable algorithms with full uncertainty quantification of predictions. In our case, the output of the learning algorithm is a probability distribution over a family of free energy functionals, consistent with the observed particle data. We find that surprisingly small data samples contain sufficient information for inferring highly accurate analytic expressions of the underlying free energy functionals, making our algorithm highly data efficient. We consider excluded volume particle interactions, which are ubiquitous in nature, whilst being highly challenging for modelling in terms of free energy. To validate our approach we consider the paradigmatic case of one-dimensional fluid and develop inference algorithms for the canonical and grand-canonical statistical-mechanical ensembles. Extensions to higher-dimensional systems are conceptually straightforward, whilst standard coarse-graining techniques allow one to easily incorporate attractive interactions. △ Less

Submitted 23 June, 2021; v1 submitted 7 October, 2020; originally announced October 2020.

Journal ref: J. Chem. Phys. (2022)

Enhancement of damaged-image prediction through Cahn-Hilliard Image Inpainting

Authors: José A. Carrillo, Serafim Kalliadasis, Fuyue Liang, Sergio P. Perez

Abstract: We assess the benefit of including an image inpainting filter before passing damaged images into a classification neural network. For this we employ a modified Cahn-Hilliard equation as an image inpainting filter, which is solved via a finite volume scheme with reduced computational cost and adequate properties for energy stability and boundedness. The benchmark dataset employed here is MNIST, whi… ▽ More We assess the benefit of including an image inpainting filter before passing damaged images into a classification neural network. For this we employ a modified Cahn-Hilliard equation as an image inpainting filter, which is solved via a finite volume scheme with reduced computational cost and adequate properties for energy stability and boundedness. The benchmark dataset employed here is MNIST, which consists of binary images of handwritten digits and is a standard dataset to validate image-processing methodologies. We train a neural network based of dense layers with the training set of MNIST, and subsequently we contaminate the test set with damage of different types and intensities. We then compare the prediction accuracy of the neural network with and without applying the Cahn-Hilliard filter to the damaged images test. Our results quantify the significant improvement of damaged-image prediction due to applying the Cahn-Hilliard filter, which for specific damages can increase up to 50% and is in general advantageous for low to moderate damage. △ Less

Submitted 15 March, 2021; v1 submitted 21 July, 2020; originally announced July 2020.

Comments: An interactive jupyter notebook with the code of this work is available at https://github.com/sergiopperez/Image_Inpainting. The MNIST dataset employed in this work can be downloaded from http://yann.lecun.com/exdb/mnist/

MSC Class: 68U10; 94A08; 65M22; 76M25; 76M12

High-order well-balanced finite volume schemes for hydrodynamic equations with nonlocal free energy

Authors: José A. Carrillo, Manuel J. Castro, Serafim Kalliadasis, Sergio P. Perez

Abstract: We propose high-order well-balanced finite-volume schemes for a broad class of hydrodynamic systems with attractive-repulsive interaction forces and linear and nonlinear damping. Our schemes are suitable for free energies containing convolutions of an interaction potential with the density, which are essential for applications such as the Keller-Segel model, more general Euler-Poisson systems, or… ▽ More We propose high-order well-balanced finite-volume schemes for a broad class of hydrodynamic systems with attractive-repulsive interaction forces and linear and nonlinear damping. Our schemes are suitable for free energies containing convolutions of an interaction potential with the density, which are essential for applications such as the Keller-Segel model, more general Euler-Poisson systems, or dynamic-density functional theory. Our schemes are also equipped with a nonnegative-density reconstruction which allows for vacuum regions during the simulation. We provide several prototypical examples from relevant applications highlighting the benefit of our algorithms elucidate also some of our analytical results. △ Less

Submitted 19 November, 2020; v1 submitted 11 April, 2020; originally announced April 2020.

MSC Class: 65XX (Primary); 35Qxx; 35Q35; 35Q82; 76M12 (Secondary)

A Finite-Volume Method for Fluctuating Dynamical Density Functional Theory

Authors: Antonio Russo, Sergio P. Perez, Miguel A. Durán-Olivencia, Peter Yatsyshin, José A. Carrillo, Serafim Kalliadasis

Abstract: We introduce a finite-volume numerical scheme for solving stochastic gradient-flow equations. Such equations are of crucial importance within the framework of fluctuating hydrodynamics and dynamic density functional theory. Our proposed scheme deals with general free-energy functionals, including, for instance, external fields or interaction potentials. This allows us to simulate a range of physic… ▽ More We introduce a finite-volume numerical scheme for solving stochastic gradient-flow equations. Such equations are of crucial importance within the framework of fluctuating hydrodynamics and dynamic density functional theory. Our proposed scheme deals with general free-energy functionals, including, for instance, external fields or interaction potentials. This allows us to simulate a range of physical phenomena where thermal fluctuations play a crucial role, such as nucleation and other energy-barrier crossing transitions. A positivity-preserving algorithm for the density is derived based on a hybrid space discretization of the deterministic and the stochastic terms and different implicit and explicit time integrators. We show through numerous applications that not only our scheme is able to accurately reproduce the statistical properties (structure factor and correlations) of the physical system, but, because of the multiplicative noise, it allows us to simulate energy barrier crossing dynamics, which cannot be captured by mean-field approaches. △ Less

Submitted 31 December, 2020; v1 submitted 11 October, 2019; originally announced October 2019.

Machine learning memory kernels as closure for non-Markovian stochastic processes

Authors: Antonio Russo, Miguel A. Duran-Olivencia, Ioannis G. Kevrekidis, Serafim Kalliadasis

Abstract: Finding the dynamical law of observable quantities lies at the core of physics. Within the particular field of statistical mechanics, the generalized Langevin equation (GLE) comprises a general model for the evolution of observables covering a great deal of physical systems with many degrees of freedom and an inherently stochastic nature. Although formally exact, the GLE brings its own great chall… ▽ More Finding the dynamical law of observable quantities lies at the core of physics. Within the particular field of statistical mechanics, the generalized Langevin equation (GLE) comprises a general model for the evolution of observables covering a great deal of physical systems with many degrees of freedom and an inherently stochastic nature. Although formally exact, the GLE brings its own great challenges. It depends on the complete history of the observables under scrutiny, as well as the microscopic degrees of freedom, all of which are often inaccessible. We show that these drawbacks can be overcome by adopting elements of machine learning from empirical data, in particular coupling a multilayer perceptron (MLP) with the formal structure of the GLE and calibrating the MLP with the data. This yields a powerful computational tool capable of describing noisy complex systems beyond the realms of statistical mechanics. It is exemplified with a number of representative examples from different fields: from a single colloidal particle and particle chains in a thermal bath to climatology and finance, showing in all cases excellent agreement with the actual observable dynamics. The new framework offers an alternative perspective for the study of non-equilibrium processes opening also a new route for stochastic modelling. △ Less

Submitted 21 November, 2022; v1 submitted 22 March, 2019; originally announced March 2019.

Journal ref: IEEE Trans Neural Netw Learn Syst (2022) 1-13

A linear, second-order, energy stable, fully adaptive finite-element method for phase-field modeling of wetting phenomena

Authors: B. Aymard, U. Vaes, M. Pradas, S. Kalliadasis

Abstract: We propose a new numerical method to solve the Cahn-Hilliard equation coupled with non-linear wetting boundary conditions. We show that the method is mass-conservative and that the discrete solution satisfies a discrete energy law similar to the one satisfied by the exact solution. We perform several tests inspired by realistic situations to verify the accuracy and performance of the method: wetti… ▽ More We propose a new numerical method to solve the Cahn-Hilliard equation coupled with non-linear wetting boundary conditions. We show that the method is mass-conservative and that the discrete solution satisfies a discrete energy law similar to the one satisfied by the exact solution. We perform several tests inspired by realistic situations to verify the accuracy and performance of the method: wetting of a chemically heterogeneous substrate in three dimensions, wetting-driven nucleation in a complex two-dimensional domain and three-dimensional diffusion through a porous medium. △ Less

Submitted 18 October, 2019; v1 submitted 18 January, 2019; originally announced January 2019.

Comments: Fix some typos and change DOI

MSC Class: 65M60

Journal ref: Journal of Computational Physics: X, Volume 2, March 2019, 100010

Well-balanced finite volume schemes for hydrodynamic equations with general free energy

Authors: José A. Carrillo, Serafim Kalliadasis, Sergio P. Perez, Chi-Wang Shu

Abstract: Well balanced and free energy dissipative first- and second-order accurate finite volume schemes are proposed for a general class of hydrodynamic systems with linear and nonlinear damping. The natural Liapunov functional of the system, given by its free energy, allows for a characterization of the stationary states by its variation. An analog property at the discrete level enables us to preserve s… ▽ More Well balanced and free energy dissipative first- and second-order accurate finite volume schemes are proposed for a general class of hydrodynamic systems with linear and nonlinear damping. The natural Liapunov functional of the system, given by its free energy, allows for a characterization of the stationary states by its variation. An analog property at the discrete level enables us to preserve stationary states at machine precision while keeping the dissipation of the discrete free energy. These schemes allow for analysing accurately the stability properties of stationary states in challeging problems such as: phase transitions in collective behavior, generalized Euler-Poisson systems in chemotaxis and astrophysics, and models in dynamic density functional theories; having done a careful validation in a battery of relevant test cases. △ Less

Submitted 4 November, 2020; v1 submitted 3 December, 2018; originally announced December 2018.

Comments: Videos from the simulations of this work are available at https://figshare.com/projects/Well-balanced_finite_volume_schemes_for_hydrodynamic_equations_with_general_free_energy/60122

Dynamics of the Desai-Zwanzig model in multi-well and random energy landscapes

Authors: Susana N. Gomes, Serafim Kalliadasis, Grigorios A. Pavliotis, Petr Yatsyshin

Abstract: We analyze a variant of the Desai-Zwanzig model [J. Stat. Phys. {\bf 19}1-24 (1978)]. In particular, we study stationary states of the mean field limit for a system of weakly interacting diffusions moving in a multi-well potential energy landscape, coupled via a Curie-Weiss type (quadratic) interaction potential. The location and depth of the local minima of the potential are either deterministic… ▽ More We analyze a variant of the Desai-Zwanzig model [J. Stat. Phys. {\bf 19}1-24 (1978)]. In particular, we study stationary states of the mean field limit for a system of weakly interacting diffusions moving in a multi-well potential energy landscape, coupled via a Curie-Weiss type (quadratic) interaction potential. The location and depth of the local minima of the potential are either deterministic or random. We characterize the structure and nature of bifurcations and phase transitions for this system, by means of extensive numerical simulations and of analytical calculations for an explicitly solvable model. Our numerical experiments are based on Monte Carlo simulations, the numerical solution of the time-dependent nonlinear Fokker-Planck (McKean-Vlasov equation), the minimization of the free energy functional and a continuation algorithm for the stationary solutions. △ Less

Submitted 6 February, 2019; v1 submitted 12 October, 2018; originally announced October 2018.

Journal ref: Phys. Rev. E 99, 032109 (2019)

Recent advances in the evolution of interfaces: thermodynamics, upscaling, and universality

Authors: M. Schmuck, G. A. Pavliotis, S. Kalliadasis

Abstract: We consider the evolution of interfaces in binary mixtures permeating strongly heterogeneous systems such as porous media. To this end, we first review available thermodynamic formulations for binary mixtures based on \emph{general reversible-irreversible couplings} and the associated mathematical attempts to formulate a \emph{non-equilibrium variational principle} in which these non-equilibrium c… ▽ More We consider the evolution of interfaces in binary mixtures permeating strongly heterogeneous systems such as porous media. To this end, we first review available thermodynamic formulations for binary mixtures based on \emph{general reversible-irreversible couplings} and the associated mathematical attempts to formulate a \emph{non-equilibrium variational principle} in which these non-equilibrium couplings can be identified as minimizers. Based on this, we investigate two microscopic binary mixture formulations fully resolving heterogeneous/perforated domains: (a) a flux-driven immiscible fluid formulation without fluid flow; (b) a momentum-driven formulation for quasi-static and incompressible velocity fields. In both cases we state two novel, reliably upscaled equations for binary mixtures/multiphase fluids in strongly heterogeneous systems by systematically taking thermodynamic features such as free energies into account as well as the system's heterogeneity defined on the microscale such as geometry and materials (e.g. wetting properties). In the context of (a), we unravel a \emph{universality} with respect to the coarsening rate due to its independence of the system's heterogeneity, i.e. the well-known ${\cal O}(t^{1/3})$-behaviour for homogeneous systems holds also for perforated domains. Finally, the versatility of phase field equations and their \emph{thermodynamic foundation} relying on free energies, make the collected recent developments here highly promising for scientific, engineering and industrial applications for which we provide an example for lithium batteries. △ Less

Submitted 24 April, 2018; originally announced April 2018.

The pressure tensor across a liquid-vapour interface

Authors: Carlos Braga, Edward Smith, Andreas Nold, David N. Sibley, Serafim Kalliadasis

Abstract: Inhomogeneous fluids exhibit physical properties that are neither uniform nor isotropic. The pressure tensor is a case in point, key to the mechanical description of the interfacial region. Kirkwood and Buff, and later Irving and Kirkwood, obtained a formal treatment based on the analysis of the pressure across a planar surface [J.G. Kirkwood and F.P. Buff, J. Chem. Phys., 17(3), (1949), J.H. Irvi… ▽ More Inhomogeneous fluids exhibit physical properties that are neither uniform nor isotropic. The pressure tensor is a case in point, key to the mechanical description of the interfacial region. Kirkwood and Buff, and later Irving and Kirkwood, obtained a formal treatment based on the analysis of the pressure across a planar surface [J.G. Kirkwood and F.P. Buff, J. Chem. Phys., 17(3), (1949), J.H. Irving and J.G. Kirkwood, J. Chem. Phys. 18, 817 (1950)]. We propose a generalisation of Irving and Kirkwood's argument to fluctuating, non-planar surfaces and obtain an expression for the pressure tensor that is not smeared by thermal fluctuations at the molecular scale and corresponding capillary waves [F.P. Buff, R.A. Lovett, and F.H. Stillinger, Jr., Phys. Rev. Lett. 15, 621 (1965)]. We observe the emergence of surface tension, defined as an excess tangential stress, acting exactly across the dividing surface at the sharpest molecular resolution. The new statistical mechanical expressions extend current treatments to fluctuating inhomogeneous systems far from equilibrium. △ Less

Submitted 27 April, 2018; v1 submitted 16 November, 2017; originally announced November 2017.

Comments: 9 pages, 4 figure. Contact: chcbraga@gmail.com

Instability, rupture and fluctuations in thin liquid films: Theory and computations

Authors: Miguel A. Durán-Olivencia, Rishabh S. Gvalani, Serafim Kalliadasis, Grigorios A. Pavliotis

Abstract: Thin liquid films are ubiquitous in natural phenomena and technological applications. They have been extensively studied via deterministic hydrodynamic equations, but thermal fluctuations often play a crucial role that needs to be understood. An example of this is dewetting, which involves the rupture of a thin liquid film and the formation of droplets. Such a process is thermally activated and re… ▽ More Thin liquid films are ubiquitous in natural phenomena and technological applications. They have been extensively studied via deterministic hydrodynamic equations, but thermal fluctuations often play a crucial role that needs to be understood. An example of this is dewetting, which involves the rupture of a thin liquid film and the formation of droplets. Such a process is thermally activated and requires fluctuations to be taken into account self-consistently. In this work we present an analytical and numerical study of a stochastic thin-film equation derived from first principles. Following a brief review of the derivation, we scrutinise the behaviour of the equation in the limit of perfectly correlated noise along the wall-normal direction. The stochastic thin-film equation is also simulated by adopting a numerical scheme based on a spectral collocation method. The scheme allows us to explore the fluctuating dynamics of the thin film and the behaviour of its free energy in the vicinity of rupture. Finally, we also study the effect of the noise intensity on the rupture time, which is in agreement with previous works. △ Less

Submitted 4 February, 2019; v1 submitted 27 July, 2017; originally announced July 2017.

Comments: 23 pages, 8 figures, Published version J Stat Phys (2019)

Two-dimensional pulse dynamics and the formation of bound states on electrified falling films

Authors: M. G. Blyth, D. Tseluiko, T. -S. Lin, S. Kalliadasis

Abstract: The flow of an electrified liquid film down an inclined plane wall is investigated with the focus on coherent structures in the form of travelling waves on the film surface, in particular, single-hump solitary pulses and their interactions. The flow structures are analysed first using a long-wave model, which is valid in the presence of weak inertia, and second using the Stokes equations. For obtu… ▽ More The flow of an electrified liquid film down an inclined plane wall is investigated with the focus on coherent structures in the form of travelling waves on the film surface, in particular, single-hump solitary pulses and their interactions. The flow structures are analysed first using a long-wave model, which is valid in the presence of weak inertia, and second using the Stokes equations. For obtuse angles, gravity is destablising and solitary pulses exist even in the absence of an electric field. For acute angles, spatially non-uniform solutions exist only beyond a critical value of the electric field strength; moreover, solitary-pulse solutions are present only at sufficiently high supercritical elec- tric field strengths. The electric field increases the amplitude of the pulses, can generate recirculation zones in the humps, and alters the far-field decay of the pulse tails from ex- ponential to algebraic with a significant impact on pulse interactions. A weak-interaction theory predicts an infinite sequence of bound-state solutions for non-electrified flow, and a finite set for electrified flow. The existence of single-hump pulse solutions and two-pulse bound states is confirmed for the Stokes equations via boundary-element computations. In addition, the electric field is shown to trigger a switch from absolute instability to convective instability, thereby regularising the dynamics, and this is confirmed by time- dependent simulations of the long-wave model. △ Less

Submitted 19 July, 2018; v1 submitted 13 June, 2017; originally announced June 2017.

Journal ref: J. Fluid Mech. (2018), vol. 855, 210-235

Rate of Convergence of General Phase Field Equations towards their Homogenized Limit

Authors: M. Schmuck, S. Kalliadasis

Abstract: Over the last few decades, phase-field equations have found increasing applicability in a wide range of mathematical-scientific fields (e.g. geometric PDEs and mean curvature flow, materials science for the study of phase transitions) but also engineering ones (e.g. as a computational tool in chemical engineering for interfacial flow studies). Here, we focus on phase-field equations in strongly he… ▽ More Over the last few decades, phase-field equations have found increasing applicability in a wide range of mathematical-scientific fields (e.g. geometric PDEs and mean curvature flow, materials science for the study of phase transitions) but also engineering ones (e.g. as a computational tool in chemical engineering for interfacial flow studies). Here, we focus on phase-field equations in strongly heterogeneous materials with perforations such as porous media. To the best of our knowledge, we provide the first derivation of error estimates for fourth order, homogenized, and nonlinear evolution equations. Our fourth order problem induces a slightly lower convergence rate, i.e., $ε^{1/4}$, where $ε$ denotes the material's specific heterogeneity, than established for second-order elliptic problems (e.g. \cite{Zhikov2006}) for the error between the effective macroscopic solution of the (new) upscaled formulation and the solution of the microscopic phase field problem. We hope that our study will motivate new modelling, analytic, and computational perspectives for interfacial transport and phase transformations in strongly heterogeneous environments. △ Less

Submitted 27 February, 2017; originally announced February 2017.

Comments: arXiv admin note: substantial text overlap with arXiv:1310.1499

Pseudospectral methods for density functional theory in bounded and unbounded domains

Authors: Andreas Nold, Benjamin D. Goddard, Peter Yatsyshin, Nikos Savva, Serafim Kalliadasis

Abstract: Classical Density Functional Theory (DFT) is a statistical-mechanical framework to analyze fluids, which accounts for nanoscale fluid inhomogeneities and non-local intermolecular interactions. DFT can be applied to a wide range of interfacial phenomena, as well as problems in adsorption, colloidal science and phase transitions in fluids. Typical DFT equations are highly non-linear, stiff and conta… ▽ More Classical Density Functional Theory (DFT) is a statistical-mechanical framework to analyze fluids, which accounts for nanoscale fluid inhomogeneities and non-local intermolecular interactions. DFT can be applied to a wide range of interfacial phenomena, as well as problems in adsorption, colloidal science and phase transitions in fluids. Typical DFT equations are highly non-linear, stiff and contain several convolution terms. We propose a novel, efficient pseudo-spectral collocation scheme for computing the non-local terms in real space with the help of a specialized Gauss quadrature. Due to the exponential accuracy of the quadrature and a convenient choice of collocation points near interfaces, we can use grids with a significantly lower number of nodes than most other reported methods. We demonstrate the capabilities of our numerical methodology by studying equilibrium and dynamic two-dimensional test cases with single- and multispecies hard-sphere and hard-disc particles modelled with fundamental measure theory, with and without van der Waals attractive forces, in bounded and unbounded physical domains. We show that our results satisfy statistical mechanical sum rules. △ Less

Submitted 22 January, 2017; originally announced January 2017.

Discrete self-similarity in interfacial hydrodynamics and the formation of iterated structures

Authors: Michael C. Dallaston, Marco A. Fontelos, Dmitri Tseluiko, Serafim Kalliadasis

Abstract: The formation of iterated structures, such as satellite and sub-satellite drops, filaments and bubbles, is a common feature in interfacial hydrodynamics. Here we undertake a computational and theoretical study of their origin in the case of thin films of viscous fluids that are destabilized by long-range molecular or other forces. We demonstrate that iterated structures appear as a consequence of… ▽ More The formation of iterated structures, such as satellite and sub-satellite drops, filaments and bubbles, is a common feature in interfacial hydrodynamics. Here we undertake a computational and theoretical study of their origin in the case of thin films of viscous fluids that are destabilized by long-range molecular or other forces. We demonstrate that iterated structures appear as a consequence of discrete self-similarity, where certain patterns repeat themselves, subject to rescaling, periodically in a logarithmic time scale. The result is an infinite sequence of ridges and filaments with similarity properties. The character of these discretely self-similar solutions as the result of a Hopf bifurcation from ordinarily self-similar solutions is also described. △ Less

Submitted 21 November, 2017; v1 submitted 19 September, 2016; originally announced September 2016.

Comments: LaTeX, 5 pages, replaced with minor changes, accepted for publication in Physical Review Letters

Journal ref: Phys. Rev. Lett. 120, 034505 (2018)

Self-similar finite-time singularity formation in degenerate parabolic equations arising in thin-film flows

Authors: Michael C Dallaston, Dmitri Tseluiko, Zhong Zheng, Marco A Fontelos, Serafim Kalliadasis

Abstract: A thin liquid film coating a planar horizontal substrate may be unstable to perturbations in the film thickness due to unfavourable intermolecular interactions between the liquid and the substrate, which may lead to finite-time rupture. The self-similar nature of the rupture has been studied before by utilizing the standard lubrication approximation along with the Derjaguin (or disjoining) pressur… ▽ More A thin liquid film coating a planar horizontal substrate may be unstable to perturbations in the film thickness due to unfavourable intermolecular interactions between the liquid and the substrate, which may lead to finite-time rupture. The self-similar nature of the rupture has been studied before by utilizing the standard lubrication approximation along with the Derjaguin (or disjoining) pressure formalism used to account for the intermolecular interactions, and a particular form of the disjoining pressure with exponent $n=3$ has been used, namely, $Π(h)\propto -1/h^{3}$, where $h$ is the film thickness. In the present study, we use a numerical continuation method to compute discrete solutions to self-similar rupture for a general disjoining pressure exponent $n$. We focus on axisymmetric point-rupture solutions and show that pairs of solution branches merge as $n$ decreases, leading to a critical value $n_c \approx 1.485$ below which stable similarity solutions do not appear to exist. We verify that this observation also holds true for plane-symmetric line-rupture solutions for which the critical value turns out to be slightly larger than for the axisymmetric case, $n_c^{\mathrm{plane}}\approx 1.5$. Computation of the full time-dependent problem also demonstrates the loss of stable similarity solutions and the subsequent onset of cascading oscillatory structures. △ Less

Submitted 13 September, 2016; originally announced September 2016.

Dynamical density functional theory for orientable colloids including inertia and hydrodynamic interactions

Authors: Miguel A. Durán-Olivencia, Benjamin D. Goddard, Serafim Kalliadasis

Abstract: Over the last few decades, classical density-functional theory (DFT) and its dynamic extensions (DDFTs) have become powerful tools in the study of colloidal fluids. Recently, previous DDFTs for spherically-symmetric particles have been generalised to take into account both inertia and hydrodynamic interactions, two effects which strongly influence non-equilibrium properties. The present work furth… ▽ More Over the last few decades, classical density-functional theory (DFT) and its dynamic extensions (DDFTs) have become powerful tools in the study of colloidal fluids. Recently, previous DDFTs for spherically-symmetric particles have been generalised to take into account both inertia and hydrodynamic interactions, two effects which strongly influence non-equilibrium properties. The present work further generalises this framework to systems of anisotropic particles. Starting from the Liouville equation and utilising Zwanzig's projection-operator techniques, we derive the kinetic equation for the Brownian particle distribution function, and by averaging over all but one particle, a DDFT equation is obtained. Whilst this equation has some similarities with DDFTs for spherically-symmetric colloids, it involves a translational-rotational coupling which affects the diffusivity of the (asymmetric) particles. We further show that, in the overdamped (high friction) limit, the DDFT is considerably simplified and is in agreement with a previous DDFT for colloids with arbitrary shape particles. △ Less

Submitted 8 July, 2016; originally announced July 2016.

Comments: dynamical density functional theory ; colloidal fluids ; arbitrary-shape particles ; orientable colloids

Journal ref: J Stat Phys (2016) 164: 785

Noise-induced transitions in rugged energy landscapes

Authors: A. B. Duncan, S. Kalliadasis, G. A. Pavliotis, M. Pradas

Abstract: We consider the problem of an overdamped Brownian particle moving in multiscale potential with N + 1 characteristic length scales: the macroscale and N separated microscales. We show that the coarse-grained dynamics is given by an overdamped Langevin equation with respect to the free energy and with a space dependent diffusion tensor, the calculation of which requires the solution of N fully coupl… ▽ More We consider the problem of an overdamped Brownian particle moving in multiscale potential with N + 1 characteristic length scales: the macroscale and N separated microscales. We show that the coarse-grained dynamics is given by an overdamped Langevin equation with respect to the free energy and with a space dependent diffusion tensor, the calculation of which requires the solution of N fully coupled Poisson equations. We study in detail the structure of the bifurcation diagram for one-dimensional problems and we show that the multiscale structure in the potential leads to hysteresis effects and to noise-induced transitions. Furthermore, we obtain an explicit formula for the effective diffusion coefficient for a self-similar separable potential and we investigate the limit of infinitely many small scales. △ Less

Submitted 14 April, 2016; originally announced April 2016.

Comments: 26 pages

Journal ref: Phys. Rev. E 94, 032107 (2016)

Controlling roughening processes in the stochastic Kuramoto-Sivashinsky equation

Authors: Susana N. Gomes, Serafim Kalliadasis, Demetrios T. Papageorgiou, Grigorios A. Pavliotis, Marc Pradas

Abstract: We present a novel control methodology to control the roughening processes of semilinear parabolic stochastic partial differential equations in one dimension, which we exemplify with the stochastic Kuramoto-Sivashinsky equation. The original equation is split into a linear stochastic and a nonlinear deterministic equation so that we can apply linear feedback control methods. Our control strategy i… ▽ More We present a novel control methodology to control the roughening processes of semilinear parabolic stochastic partial differential equations in one dimension, which we exemplify with the stochastic Kuramoto-Sivashinsky equation. The original equation is split into a linear stochastic and a nonlinear deterministic equation so that we can apply linear feedback control methods. Our control strategy is then based on two steps: first, stabilize the zero solution of the deterministic part and, second, control the roughness of the stochastic linear equation. We consider both periodic controls and point actuated ones, observing in all cases that the second moment of the solution evolves in time according to a power-law until it saturates at the desired controlled value. △ Less

Submitted 4 January, 2017; v1 submitted 7 April, 2016; originally announced April 2016.

Nanoscale fluid structure of liquid-solid-vapour contact lines for a wide range of contact angles

Authors: Andreas Nold, David N. Sibley, Benjamin D. Goddard, Serafim Kalliadasis

Abstract: We study the nanoscale behaviour of the density of a simple fluid in the vicinity of an equilibrium contact line for a wide range of Young contact angles between 40 and 135 degrees. Cuts of the density profile at various positions along the contact line are presented, unravelling the apparent step-wise increase of the film height profile observed in contour plots of the density. The density profil… ▽ More We study the nanoscale behaviour of the density of a simple fluid in the vicinity of an equilibrium contact line for a wide range of Young contact angles between 40 and 135 degrees. Cuts of the density profile at various positions along the contact line are presented, unravelling the apparent step-wise increase of the film height profile observed in contour plots of the density. The density profile is employed to compute the normal pressure acting on the substrate along the contact line. We observe that for the full range of contact angles, the maximal normal pressure cannot solely be predicted by the curvature of the adsorption film height, but is instead softened -- likely by the width of the liquid-vapour interface. Somewhat surprisingly however, the adsorption film height profile can be predicted to a very good accuracy by the Derjaguin-Frumkin disjoining pressure obtained from planar computations, as was first shown in [Nold et al., Phys. Fluids, 26, 072001, 2014] for contact angles less than 90 degrees, a result which here we show to be valid for the full range of contact angles. This suggests that while two-dimensional effects cannot be neglected for the computation of the normal pressure distribution along the substrate, one-dimensional planar computations of the Derjaguin-Frumkin disjoining pressure are sufficient to accurately predict the adsorption height profile. △ Less

Submitted 5 August, 2017; v1 submitted 11 March, 2015; originally announced March 2015.

Comments: 19 pages, 6 figures

Journal ref: Mathematical Modelling of Natural Phenomena, 10(4), 111-125 (2015)

Nonlinear forecasting of the generalised Kuramoto-Sivashinsky equation

Authors: H. Gotoda, M. Pradas, S. Kalliadasis

Abstract: We study the emergence of pattern formation and chaotic dynamics in the one-dimensional (1D) generalized Kuramoto-Sivashinsky (gKS) equation by means of a time-series analysis, in particular a nonlinear forecasting method which is based on concepts from chaos theory and appropriate statistical methods. We analyze two types of temporal signals, a local one and a global one, finding in both cases th… ▽ More We study the emergence of pattern formation and chaotic dynamics in the one-dimensional (1D) generalized Kuramoto-Sivashinsky (gKS) equation by means of a time-series analysis, in particular a nonlinear forecasting method which is based on concepts from chaos theory and appropriate statistical methods. We analyze two types of temporal signals, a local one and a global one, finding in both cases that the dynamical state of the gKS solution undergoes a transition from high-dimensional chaos to periodic pulsed oscillations through low-dimensional deterministic chaos with increasing the control parameter of the system. Our results demonstrate that the proposed nonlinear forecasting methodology allows to elucidate the dynamics of the system in terms of its predictability properties. △ Less

Submitted 18 February, 2015; originally announced February 2015.

Comments: 20 pages, 15 figures

A new framework for extracting coarse-grained models from time series with multiscale structure

Authors: Serafim Kalliadasis, Sebastian Krumscheid, Grigorios A. Pavliotis

Abstract: In many applications it is desirable to infer coarse-grained models from observational data. The observed process often corresponds only to a few selected degrees of freedom of a high-dimensional dynamical system with multiple time scales. In this work we consider the inference problem of identifying an appropriate coarse-grained model from a single time series of a multiscale system. It is known… ▽ More In many applications it is desirable to infer coarse-grained models from observational data. The observed process often corresponds only to a few selected degrees of freedom of a high-dimensional dynamical system with multiple time scales. In this work we consider the inference problem of identifying an appropriate coarse-grained model from a single time series of a multiscale system. It is known that estimators such as the maximum likelihood estimator or the quadratic variation of the path estimator can be strongly biased in this setting. Here we present a novel parametric inference methodology for problems with linear parameter dependency that does not suffer from this drawback. Furthermore, we demonstrate through a wide spectrum of examples that our methodology can be used to derive appropriate coarse-grained models from time series of partial observations of a multiscale system in an effective and systematic fashion. △ Less

Submitted 5 May, 2015; v1 submitted 5 September, 2014; originally announced September 2014.

Fluid structure in the immediate vicinity of an equilibrium three-phase contact line and assessment of disjoining pressure models using density functional theory

Authors: Andreas Nold, David N. Sibley, Benjamin D. Goddard, Serafim Kalliadasis

Abstract: We examine the nanoscale behavior of an equilibrium three-phase contact line in the presence of long-ranged intermolecular forces by employing a statistical mechanics of fluids approach, namely density functional theory (DFT) together with fundamental measure theory (FMT). This enables us to evaluate the predictive quality of effective Hamiltonian models in the vicinity of the contact line. In par… ▽ More We examine the nanoscale behavior of an equilibrium three-phase contact line in the presence of long-ranged intermolecular forces by employing a statistical mechanics of fluids approach, namely density functional theory (DFT) together with fundamental measure theory (FMT). This enables us to evaluate the predictive quality of effective Hamiltonian models in the vicinity of the contact line. In particular, we compare the results for mean field effective Hamiltonians with disjoining pressures defined through (I) the adsorption isotherm for a planar liquid film, and (II) the normal force balance at the contact line. We find that the height profile obtained using (I) shows good agreement with the adsorption film thickness of the DFT-FMT equilibrium density profile in terms of maximal curvature and the behavior at large film heights. In contrast, we observe that while the height profile obtained by using (II) satisfies basic sum rules, it shows little agreement with the adsorption film thickness of the DFT results. The results are verified for contact angles of 20, 40 and 60 degrees. △ Less

Submitted 4 August, 2014; v1 submitted 20 June, 2014; originally announced June 2014.

Journal ref: Phys. Fluids 26, 072001 (2014)

Rate of Convergence of Phase Field Equations in Strongly Heterogeneous Media towards their Homogenized Limit

Authors: Markus Schmuck, Grigorios A. Pavliotis, Serafim Kalliadasis

Abstract: We study phase field equations based on the diffuse-interface approximation of general homogeneous free energy densities showing different local minima of possible equilibrium configurations in perforated/porous domains. The study of such free energies in homogeneous environments found a broad interest over the last decades and hence is now widely accepted and applied in both science and engineeri… ▽ More We study phase field equations based on the diffuse-interface approximation of general homogeneous free energy densities showing different local minima of possible equilibrium configurations in perforated/porous domains. The study of such free energies in homogeneous environments found a broad interest over the last decades and hence is now widely accepted and applied in both science and engineering. Here, we focus on strongly heterogeneous materials with perforations such as porous media. To the best of our knowledge, we present a general formal derivation of upscaled phase field equations for arbitrary free energy densities and give a rigorous justification by error estimates for a broad class of polynomial free energies. The error between the effective macroscopic solution of the new upscaled formulation and the solution of the microscopic phase field problem is of order $ε^1/2$ for a material given characteristic heterogeneity $ε$. Our new, effective, and reliable macroscopic porous media formulation of general phase field equations opens new modelling directions and computational perspectives for interfacial transport in strongly heterogeneous environments. △ Less

Submitted 5 October, 2013; originally announced October 2013.

The contact line behaviour of solid-liquid-gas diffuse-interface models

Authors: David N. Sibley, Andreas Nold, Nikos Savva, Serafim Kalliadasis

Abstract: A solid-liquid-gas moving contact line is considered through a diffuse-interface model with the classical boundary condition of no-slip at the solid surface. Examination of the asymptotic behaviour as the contact line is approached shows that the relaxation of the classical model of a sharp liquid-gas interface, whilst retaining the no-slip condition, resolves the stress and pressure singularities… ▽ More A solid-liquid-gas moving contact line is considered through a diffuse-interface model with the classical boundary condition of no-slip at the solid surface. Examination of the asymptotic behaviour as the contact line is approached shows that the relaxation of the classical model of a sharp liquid-gas interface, whilst retaining the no-slip condition, resolves the stress and pressure singularities associated with the moving contact line problem while the fluid velocity is well defined (not multi-valued). The moving contact line behaviour is analysed for a general problem relevant for any density dependent dynamic viscosity and volume viscosity, and for general microscopic contact angle and double well free-energy forms. Away from the contact line, analysis of the diffuse-interface model shows that the Navier--Stokes equations and classical interfacial boundary conditions are obtained at leading order in the sharp-interface limit, justifying the creeping flow problem imposed in an intermediate region in the seminal work of Seppecher [Int. J. Eng. Sci. 34, 977--992 (1996)]. Corrections to Seppecher's work are given, as an incorrect solution form was originally used. △ Less

Submitted 4 October, 2013; originally announced October 2013.

Comments: 33 pages, 3 figures

Journal ref: Phys. Fluids 25, 092111 (2013)

Wetting in a two-dimensional capped capillary. Part II: Three-phase coexistence

Authors: P. Yatsyshin, N. Savva, S. Kalliadasis

Abstract: In Part II of this study we consider two cases of three-phase coexistence. First, the capped capillary may allow for vapour, drop-like, and slab-like phases to coexist at the same values of temperature and chemical potential. Second, the slit pore forming the bulk of the capped capillary may allow for the coexistence between vapour, planar prewetting film and capillary-liquid. While the considerat… ▽ More In Part II of this study we consider two cases of three-phase coexistence. First, the capped capillary may allow for vapour, drop-like, and slab-like phases to coexist at the same values of temperature and chemical potential. Second, the slit pore forming the bulk of the capped capillary may allow for the coexistence between vapour, planar prewetting film and capillary-liquid. While the consideration of the former case allows us to summarise the phenomenology presented in Part I and to show that the transition line of wedge prewetting is shifted in capillary-like geometries by a constant value, depending on the capillary width, the careful examination of the latter case allows us to uncover a new phase transition in confined fluids, a continuous planar prewetting transition. A planar prewetting transition is known to be a distinctly first-order phenomenon, and typically taking place on the scale of several atomic diameters. A continuous prewetting transition, on the other hand, is scale invariant. Thus, apart from being of fundamental significance, this finding has potential for facilitating experimental detection as well as measurements of planar prewetting. Further, we provide proof for the existence of a tri-critical point of the three-phase coexistence line of the capped capillary while by considering a dynamic model of wetting we show how the relaxation of the system can be pinned by a metastable state. We present a full parametric study of our model system and support our findings with exhaustive examples of density profiles, adsorption and free energy isotherms, and full phase diagrams. △ Less

Submitted 20 September, 2013; originally announced September 2013.

Comments: 27 pages, 9 figures

Wetting in a two-dimensional capped capillary. Part I: Wetting temperature and capillary prewetting

Authors: P. Yatsyshin, N. Savva, S. Kalliadasis

Abstract: In this two-part study we investigate the phase behaviour of a fluid spatially confined in a semi-infinite rectangular pore formed by three orthogonal walls and connected to a reservoir maintaining constant values of pressure and temperature in the fluid. Far from the capping wall this prototypical two-dimensional system reduces to a one-dimensional slit pore. However, the broken translational sym… ▽ More In this two-part study we investigate the phase behaviour of a fluid spatially confined in a semi-infinite rectangular pore formed by three orthogonal walls and connected to a reservoir maintaining constant values of pressure and temperature in the fluid. Far from the capping wall this prototypical two-dimensional system reduces to a one-dimensional slit pore. However, the broken translational symmetry leads to a wetting behavior strikingly different from that of a slit pore. Using a realistic model of an atomic fluid with long-ranged Lennard-Jones fluid-fluid and fluid-substrate interactions, we present for the first time detailed computations of full phase diagrams of two-dimensional capped capillaries. Our analysis is based on the statistical mechanics of fluids, in particular density functional theory. We show the existence of capillary wetting temperature, which is a property of the pore, and relatively to the fluid temperature determines whether capillary condensation is a first-order or a continuous phase transition. We also report for the first time a first-order capillary wetting transition, which can be preceded by a first-order capillary prewetting. A full parametric study is undertaken and we support our findings with exhaustive examples of density profiles, adsorption and free energy isotherms, as well as full phase diagrams. △ Less

Submitted 20 September, 2013; originally announced September 2013.

Comments: 34 pages, 10 figures

New and general framework for adsorption processes on dynamic interfaces

Authors: Markus Schmuck, Serafim Kalliadasis

Abstract: We introduce a new and general continuum thermodynamic framework for the mathematical analysis and computation of adsorption on dynamic interfaces. To the best of our knowledge, there is no formulation available that accounts for the coupled dynamics of interfaces and densities of adsorbants. Our framework leads to analytic adsorption isotherms which also take the interfacial geometry fully into a… ▽ More We introduce a new and general continuum thermodynamic framework for the mathematical analysis and computation of adsorption on dynamic interfaces. To the best of our knowledge, there is no formulation available that accounts for the coupled dynamics of interfaces and densities of adsorbants. Our framework leads to analytic adsorption isotherms which also take the interfacial geometry fully into account. We demonstrate the utility and physical consistency of our framework with a new computational multi-level discretization strategy. In the computations, we recover the experimentally observed feature that the adsorption of particles minimizes the interfacial tension. △ Less

Submitted 30 July, 2013; originally announced July 2013.

A comparison of slip, disjoining pressure, and interface formation models for contact line motion through asymptotic analysis of thin two-dimensional droplet spreading

Authors: David N. Sibley, Andreas Nold, Nikos Savva, Serafim Kalliadasis

Abstract: The motion of a contact line is examined, and comparisons drawn, for a variety of models proposed in the literature. Pressure and stress behaviours at the contact line are examined in the prototype system of quasistatic spreading of a thin two-dimensional droplet on a planar substrate. The models analysed include three disjoining pressure models based on van der Waals interactions, a model introdu… ▽ More The motion of a contact line is examined, and comparisons drawn, for a variety of models proposed in the literature. Pressure and stress behaviours at the contact line are examined in the prototype system of quasistatic spreading of a thin two-dimensional droplet on a planar substrate. The models analysed include three disjoining pressure models based on van der Waals interactions, a model introduced for polar fluids, and a liquid-gas diffuse-interface model; Navier-slip and two non-linear slip models are investigated, with three microscopic contact angle boundary conditions imposed (two of these contact angle conditions having a contact line velocity dependence); and the interface formation model is also considered. In certain parameter regimes it is shown that all of the models predict the same quasistatic droplet spreading behaviour. △ Less

Submitted 1 September, 2014; v1 submitted 30 July, 2013; originally announced July 2013.

Comments: 29 pages, 3 figures, J. Eng. Math. 2014

A new stochastic mode reduction strategy for dissipative systems

Authors: M. Schmuck, M. Pradas, S. Kalliadasis, G. A. Pavliotis

Abstract: We present a new methodology for studying non-Hamiltonian nonlinear systems based on an information theoretic extension of a renormalization group technique using a modified maximum entropy principle. We obtain a rigorous dimensionally reduced description for such systems. The neglected degrees of freedom by this reduction are replaced by a systematically deefined stochastic process under a constr… ▽ More We present a new methodology for studying non-Hamiltonian nonlinear systems based on an information theoretic extension of a renormalization group technique using a modified maximum entropy principle. We obtain a rigorous dimensionally reduced description for such systems. The neglected degrees of freedom by this reduction are replaced by a systematically deefined stochastic process under a constraint on the second moment. This then forms the basis of a computationally efficient method. Numerical computations for the generalized Kuramoto-Sivashinsky equation sup- port our method and reveal that the long-time underlying stochastic process of the fast (unresolved) modes obeys a universal distribution which does not depend on the initial conditions and which we rigorously derive by the maximum entropy principle. △ Less

Submitted 17 May, 2013; originally announced May 2013.

Journal ref: Phys. Rev. Lett. 110, 244101 (2013) [5 pages]

Derivation of effective macroscopic Stokes-Cahn-Hilliard equations for periodic immiscible flows in porous media

Authors: Markus Schmuck, Marc Pradas, Gregorios A. Pavliotis, Serafim Kalliadasis

Abstract: Using thermodynamic and variational principles we examine a basic phase field model for a mixture of two incompressible fluids in strongly perforated domains. With the help of the multiple scale method with drift and our recently introduced splitting strategy for Ginzburg-Landau/Cahn-Hilliard-type equations [Schmuck et al., Proc. R. Soc. A 468:3705-3724, 2012.], we rigorously derive an effective m… ▽ More Using thermodynamic and variational principles we examine a basic phase field model for a mixture of two incompressible fluids in strongly perforated domains. With the help of the multiple scale method with drift and our recently introduced splitting strategy for Ginzburg-Landau/Cahn-Hilliard-type equations [Schmuck et al., Proc. R. Soc. A 468:3705-3724, 2012.], we rigorously derive an effective macroscopic phase field formulation under the assumption of periodic flow and a sufficiently large Péclet number. As for classical convection-diffusion problems, we obtain systematically diffusion-dispersion relations (including Taylor-Aris-dispersion). Our results also provide a convenient analytical and computational framework to macroscopically track interfaces in porous media. In view of the well-known versatility of phase field models, our study proposes a promising model for many engineering and scientific applications such as multiphase flows in porous media, microfluidics, and fuel cells. △ Less

Submitted 19 October, 2013; v1 submitted 23 October, 2012; originally announced October 2012.

Comments: 25 pages, 5 figures

Journal ref: Nonlinearity 26 (2013) 3259-3277

On the moving contact line singularity: Asymptotics of a diffuse-interface model

Authors: David N. Sibley, Andreas Nold, Nikos Savva, Serafim Kalliadasis

Abstract: The behaviour of a solid-liquid-gas system near the three-phase contact line is considered using a diffuse-interface model with no-slip at the solid and where the fluid phase is specified by a continuous density field. Relaxation of the classical approach of a sharp liquid-gas interface and careful examination of the asymptotic behaviour as the contact line is approached is shown to resolve the st… ▽ More The behaviour of a solid-liquid-gas system near the three-phase contact line is considered using a diffuse-interface model with no-slip at the solid and where the fluid phase is specified by a continuous density field. Relaxation of the classical approach of a sharp liquid-gas interface and careful examination of the asymptotic behaviour as the contact line is approached is shown to resolve the stress and pressure singularities associated with the moving contact line problem. Various features of the model are scrutinised, alongside extensions to incorporate slip, finite-time relaxation of the chemical potential, or a precursor film at the wall. △ Less

Submitted 1 October, 2013; v1 submitted 5 October, 2012; originally announced October 2012.

Comments: 14 pages, 3 figures

Journal ref: Eur. Phys. J. E (2013) 36: 26

Geometry-induced phase transition in fluids: capillary prewetting

Authors: Petr Yatsyshin, Nikos Savva, Serafim Kalliadasis

Abstract: We report a new first-order phase transition preceding capillary condensation and corresponding to the discontinuous formation of a curved liquid meniscus. Using a mean-field microscopic approach based on the density functional theory we compute the complete phase diagram of a prototypical two-dimensional system exhibiting capillary condensation, namely that of a fluid with long-ranged dispersion… ▽ More We report a new first-order phase transition preceding capillary condensation and corresponding to the discontinuous formation of a curved liquid meniscus. Using a mean-field microscopic approach based on the density functional theory we compute the complete phase diagram of a prototypical two-dimensional system exhibiting capillary condensation, namely that of a fluid with long-ranged dispersion intermolecular forces which is spatially confined by a substrate forming a semi-infinite rectangular pore exerting long-ranged dispersion forces on the fluid. In the T-mu plane the phase line of the new transition is tangential to the capillary condensation line at the capillary wetting temperature, Tcw. The surface phase behavior of the system maps to planar wetting with the phase line of the new transition, termed capillary prewetting, mapping to the planar prewetting line. If capillary condensation is approached isothermally with T>Tcw, the meniscus forms at the capping wall and unbinds continuously, making capillary condensation a second-order phenomenon. We compute the corresponding critical exponent for the divergence of adsorption. △ Less

Submitted 1 March, 2013; v1 submitted 4 September, 2012; originally announced September 2012.

Comments: 5 pages, 4 figures, 5 movies

Journal ref: Phys. Rev. E (Rapid Communications), 87, 020402(R) (2012)

Slip or not slip? A methodical examination of the interface formation model using two-dimensional droplet spreading on a horizontal planar substrate as a prototype system

Authors: David N. Sibley, Nikos Savva, Serafim Kalliadasis

Abstract: We consider the spreading of a thin two-dimensional droplet on a planar substrate as a prototype system to compare the contemporary model for contact line motion based on interface formation of Shikhmurzaev [Int. J. Multiphas. Flow 19, 589 (1993)], to the more commonly used continuum fluid dynamical equations augmented with the Navier-slip condition. Considering quasistatic droplet evolution and u… ▽ More We consider the spreading of a thin two-dimensional droplet on a planar substrate as a prototype system to compare the contemporary model for contact line motion based on interface formation of Shikhmurzaev [Int. J. Multiphas. Flow 19, 589 (1993)], to the more commonly used continuum fluid dynamical equations augmented with the Navier-slip condition. Considering quasistatic droplet evolution and using the method of matched asymptotics, we find that the evolution of the droplet radius using the interface formation model reduces to an equivalent expression for a slip model, where the prescribed microscopic dynamic contact angle has a velocity dependent correction to its static value. This result is found for both the original interface formation model formulation and for a more recent version, where mass transfer from bulk to surface layers is accounted for through the boundary conditions. Various features of the model, such as the pressure behaviour and rolling motion at the contact line, and their relevance, are also considered in the prototype system we adopt. △ Less

Submitted 9 July, 2012; originally announced July 2012.

Comments: 45 pages, 18 figures

Journal ref: Phys. Fluids 24, 082105 (2012)

Upscaled phase-field models for interfacial dynamics in strongly heterogeneous domains

Authors: Markus Schmuck, Marc Pradas, Greg A. Pavliotis, Serafim Kalliadasis

Abstract: We derive a new effective macroscopic Cahn-Hilliard equation whose homogeneous free energy is represented by 4-th order polynomials, which form the frequently applied double-well potential. This upscaling is done for perforated/strongly het- erogeneous domains. To the best knowledge of the authors, this seems to be the first attempt of upscaling the Cahn-Hilliard equation in such domains. The new… ▽ More We derive a new effective macroscopic Cahn-Hilliard equation whose homogeneous free energy is represented by 4-th order polynomials, which form the frequently applied double-well potential. This upscaling is done for perforated/strongly het- erogeneous domains. To the best knowledge of the authors, this seems to be the first attempt of upscaling the Cahn-Hilliard equation in such domains. The new homog- enized equation should have a broad range of applicability due to the well-known versatility of phase-field models. The additionally introduced feature of systemati- cally and reliably accounting for confined geometries by homogenization allows for new modeling and numerical perspectives in both, science and engineering. Our results are applied to wetting dynamics in porous media and to a single channel with strongly heterogeneous walls. △ Less

Submitted 25 May, 2012; v1 submitted 21 April, 2012; originally announced April 2012.

Journal ref: Proc. R. Soc. A 468(2147):3705-3724 2012

Generalized dynamical density functional theory for classical fluids and the significance of inertia and hydrodynamic interactions

Authors: Benjamin D. Goddard, Andreas Nold, Nikos Savva, Grigorios A. Pavliotis, Serafim Kalliadasis

Abstract: We study the dynamics of a colloidal fluid including inertia and hydrodynamic interactions, two effects which strongly influence the non-equilibrium properties of the system. We derive a general dynamical density functional theory (DDFT) which shows very good agreement with full Langevin dynamics. In suitable limits, we recover existing DDFTs and a Navier-Stokes-like equation with additional non-l… ▽ More We study the dynamics of a colloidal fluid including inertia and hydrodynamic interactions, two effects which strongly influence the non-equilibrium properties of the system. We derive a general dynamical density functional theory (DDFT) which shows very good agreement with full Langevin dynamics. In suitable limits, we recover existing DDFTs and a Navier-Stokes-like equation with additional non-local terms. △ Less

Submitted 8 August, 2012; v1 submitted 15 February, 2012; originally announced February 2012.

Comments: 5 pages, 4 figures, 4 supplementary movie files, I supplementary pdf

Journal ref: Phys. Rev. Lett. 109, 120603 (2012)

A New Mode Reduction Strategy for the Generalized Kuramoto-Sivashinsky Equation

Authors: M. Schmuck, M. Pradas, G. A. Pavliotis, S. Kalliadasis

Abstract: Consider the generalized Kuramoto-Sivashinsky (gKS) equation. It is a model prototype for a wide variety of physical systems, from flame-front propagation, and more general front propagation in reaction-diffusion systems, to interface motion of viscous film flows. Our aim is to develop a systematic and rigorous low-dimensional representation of the gKS equation. For this purpose, we approximate it… ▽ More Consider the generalized Kuramoto-Sivashinsky (gKS) equation. It is a model prototype for a wide variety of physical systems, from flame-front propagation, and more general front propagation in reaction-diffusion systems, to interface motion of viscous film flows. Our aim is to develop a systematic and rigorous low-dimensional representation of the gKS equation. For this purpose, we approximate it by a renormalization group (RG) equation which is qualitatively characterized by rigorous error bounds. This formulation allows for a new stochastic mode reduction guaranteeing optimality in the sense of maximal information entropy. Herewith, noise is systematically added to the reduced gKS equation and gives a rigorous and analytical explanation for its origin. These new results would allow to reliably perform low-dimensional numerical computations by accounting for the neglected degrees of freedom in a systematic way. Moreover, the presented reduction strategy might also be useful in other applications where classical mode reduction approaches fail or are too complicated to be implemented. △ Less

Submitted 25 September, 2013; v1 submitted 8 November, 2011; originally announced November 2011.

Comments: 28 pages

Semi-Parametric Drift and Diffusion Estimation for Multiscale Diffusions

Authors: Sebastian Krumscheid, Grigorios A. Pavliotis, Serafim Kalliadasis

Abstract: We consider the problem of statistical inference for the effective dynamics of multiscale diffusion processes with (at least) two widely separated characteristic time scales. More precisely, we seek to determine parameters in the effective equation describing the dynamics on the longer diffusive time scale, i.e. in a homogenization framework. We examine the case where both the drift and the diffus… ▽ More We consider the problem of statistical inference for the effective dynamics of multiscale diffusion processes with (at least) two widely separated characteristic time scales. More precisely, we seek to determine parameters in the effective equation describing the dynamics on the longer diffusive time scale, i.e. in a homogenization framework. We examine the case where both the drift and the diffusion coefficients in the effective dynamics are space-dependent and depend on multiple unknown parameters. It is known that classical estimators, such as Maximum Likelihood and Quadratic Variation of the Path Estimators, fail to obtain reasonable estimates for parameters in the effective dynamics when based on observations of the underlying multiscale diffusion. We propose a novel algorithm for estimating both the drift and diffusion coefficients in the effective dynamics based on a semi-parametric framework. We demonstrate by means of extensive numerical simulations of a number of selected examples that the algorithm performs well when applied to data from a multiscale diffusion. These examples also illustrate that the algorithm can be used effectively to obtain accurate and unbiased estimates. △ Less

Submitted 20 March, 2013; v1 submitted 8 November, 2011; originally announced November 2011.

Comments: 32 pages, 10 figures

Journal ref: Multiscale Model. Simul., 11(2):442-473, 2013

Wavy regimes of film flow down a fibre

Authors: Christian Ruyer-Quil, Serafim Kalliadasis

Abstract: We consider axisymmetric traveling waves propagating on the gravity-driven flow of a liquid down a vertical fibre. Our starting point is the two-equation model for the flow derived in the study by Ruyer-Quil \emph{et al.} [\emph{J. Fluid Mech.} {\bf 603}, 431 (2008)]. The speed, amplitude and shape of the traveling waves are obtained for a wide range of parameters by using asymptotic analysis and… ▽ More We consider axisymmetric traveling waves propagating on the gravity-driven flow of a liquid down a vertical fibre. Our starting point is the two-equation model for the flow derived in the study by Ruyer-Quil \emph{et al.} [\emph{J. Fluid Mech.} {\bf 603}, 431 (2008)]. The speed, amplitude and shape of the traveling waves are obtained for a wide range of parameters by using asymptotic analysis and elements from dynamical systems theory. Four different regimes are identified corresponding to the predominance of four different physical effects: Advection by the flow, azimuthal curvature, inertia and viscous dispersion. Construction of the traveling-wave branches of solutions reveals complex transitions from one regime to another. A phase diagram of the different regimes in the parameter space is oferred. △ Less

Submitted 7 October, 2011; originally announced October 2011.

Comments: Paper submitted to Phys. Rev. E the 5th of October 2011

Journal ref: Phys. Rev. E 85, 046302 (2012)

The overdamped limit of dynamic density functional theory: Rigorous results

Authors: Benjamin D. Goddard, Grigorios A. Pavliotis, Serafim Kalliadasis

Abstract: Consider the overdamped limit for a system of interacting particles in the presence of hydrodynamic interactions. For two-body hydrodynamic interactions and one- and two-body potentials, a Smoluchowski-type evolution equation is rigorously derived for the one-particle distribution function. This new equation includes a novel definition of the diffusion tensor. A comparison with existing formulatio… ▽ More Consider the overdamped limit for a system of interacting particles in the presence of hydrodynamic interactions. For two-body hydrodynamic interactions and one- and two-body potentials, a Smoluchowski-type evolution equation is rigorously derived for the one-particle distribution function. This new equation includes a novel definition of the diffusion tensor. A comparison with existing formulations of dynamic density functional theory is also made. △ Less

Submitted 8 August, 2012; v1 submitted 16 August, 2011; originally announced August 2011.

Comments: 31 pages, 1 figure

MSC Class: 82C22; 76M45; 76M50; 35B40

Journal ref: Mulitscale Model. Simul. 10 (2) 633-663, 2012

Additive noise effects in active nonlinear spatially extended systems

Authors: M. Pradas, G. A. Pavliotis, S. Kalliadasis, D. T. Papageorgiou, D. Tseluiko

Abstract: We examine the effects of pure additive noise on spatially extended systems with quadratic nonlinearities. We develop a general multiscale theory for such systems and apply it to the Kuramoto-Sivashinsky equation as a case study. We first focus on a regime close to the instability onset (primary bifurcation), where the system can be described by a single dominant mode. We show analytically that th… ▽ More We examine the effects of pure additive noise on spatially extended systems with quadratic nonlinearities. We develop a general multiscale theory for such systems and apply it to the Kuramoto-Sivashinsky equation as a case study. We first focus on a regime close to the instability onset (primary bifurcation), where the system can be described by a single dominant mode. We show analytically that the resulting noise in the equation describing the amplitude of the dominant mode largely depends on the nature of the stochastic forcing. For a highly degenerate noise, in the sense that it is acting on the first stable mode only, the amplitude equation is dominated by a pure multiplicative noise, which in turn induces the dominant mode to undergo several critical state transitions and complex phenomena, including intermittency and stabilisation, as the noise strength is increased. The intermittent behaviour is characterised by a power-law probability density and the corresponding critical exponent is calculated rigorously by making use of the first-passage properties of the amplitude equation. On the other hand, when the noise is acting on the whole subspace of stable modes, the multiplicative noise is corrected by an additive-like term, with the eventual loss of any stabilised state. We also show that the stochastic forcing has no effect on the dominant mode dynamics when it is acting on the second stable mode. Finally, in a regime which is relatively far from the instability onset, so that there are two unstable modes, we observe numerically that when the noise is acting on the first stable mode, both dominant modes show noise-induced complex phenomena similar to the single-mode case. △ Less

Submitted 19 September, 2012; v1 submitted 3 August, 2011; originally announced August 2011.

Journal ref: European Journal of Applied Mathematics, 23, pp 563-591 (2012)

Wetting on a spherical wall: influence of liquid-gas interfacial properties

Authors: Andreas Nold, Alexandr Malijevský, Serafim Kalliadasis

Abstract: We study the equilibrium of a liquid film on an attractive spherical substrate for an intermolecular interaction model exhibiting both fluid-fluid and fluid-wall long-range forces. We first reexamine the wetting properties of the model in the zero-curvature limit, i.e., for a planar wall, using an effective interfacial Hamiltonian approach in the framework of the well known sharp-kink approximatio… ▽ More We study the equilibrium of a liquid film on an attractive spherical substrate for an intermolecular interaction model exhibiting both fluid-fluid and fluid-wall long-range forces. We first reexamine the wetting properties of the model in the zero-curvature limit, i.e., for a planar wall, using an effective interfacial Hamiltonian approach in the framework of the well known sharp-kink approximation (SKA). We obtain very good agreement with a mean-field density functional theory (DFT), fully justifying the use of SKA in this limit. We then turn our attention to substrates of finite curvature and appropriately modify the so-called soft-interface approximation (SIA) originally formulated by Napiórkowski and Dietrich [Phys. Rev. B 34, 6469 (1986)] for critical wetting on a planar wall. A detailed asymptotic analysis of SIA confirms the SKA functional form for the film growth. However, it turns out that the agreement between SKA and our DFT is only qualitative. We then show that the quantitative discrepancy between the two is due to the overestimation of the liquid-gas surface tension within SKA. On the other hand, by relaxing the assumption of a sharp interface, with, e.g., a simple smoothing of the density profile there, markedly improves the predictive capability of the theory, making it quantitative and showing that the liquid-gas surface tension plays a crucial role when describing wetting on a curved substrate. In addition, we show that in contrast to SKA, SIA predicts the expected mean-field critical exponent of the liquid-gas surface tension. △ Less

Submitted 2 October, 2013; v1 submitted 31 March, 2011; originally announced March 2011.

Journal ref: Phys. Rev. E 84, 021603 (2011)

Noise induced state transitions, intermittency and universality in the noisy Kuramoto-Sivashinsky equation

Authors: M. Pradas, D. Tseluiko, S. Kalliadasis, D. T. Papageorgiou, G. A. Pavliotis

Abstract: We analyze the effect of pure additive noise on the long-time dynamics of the noisy Kuramoto-Sivashinsky (KS) equation in a regime close to the instability onset. We show that when the noise is highly degenerate, in the sense that it acts only on the first stable mode, the solution of the KS equation undergoes several transitions between different states, including a critical on-off intermittent s… ▽ More We analyze the effect of pure additive noise on the long-time dynamics of the noisy Kuramoto-Sivashinsky (KS) equation in a regime close to the instability onset. We show that when the noise is highly degenerate, in the sense that it acts only on the first stable mode, the solution of the KS equation undergoes several transitions between different states, including a critical on-off intermittent state that is eventually stabilized as the noise strength is increased. Such noise-induced transitions can be completely characterized through critical exponents, obtaining that both the KS and the noisy Burgers equation belong to the same universality class. The results of our numerical investigations are explained rigorously using multiscale techniques. △ Less

Submitted 19 September, 2012; v1 submitted 5 November, 2010; originally announced November 2010.

Comments: 4 pages, 4 figures

Journal ref: Physical Review Letters 106, 060602 (2011)

Rigorous coherent-structure theory for falling liquid films: Viscous dispersion effects on bound-state formation and self-organization

Authors: Marc Pradas, Dmitri Tseluiko, Serafim Kalliadasis

Abstract: We examine the interaction of two-dimensional solitary pulses on falling liquid films. We make use of the second-order model derived by Ruyer-Quil and Manneville [Eur. Phys. J. B 6, 277 (1998); Eur. Phys. J. B 15, 357 (2000); Phys. Fluids 14, 170 (2002)] by combining the long-wave approximation with a weighted residuals technique. The model includes (second-order) viscous dispersion effects which… ▽ More We examine the interaction of two-dimensional solitary pulses on falling liquid films. We make use of the second-order model derived by Ruyer-Quil and Manneville [Eur. Phys. J. B 6, 277 (1998); Eur. Phys. J. B 15, 357 (2000); Phys. Fluids 14, 170 (2002)] by combining the long-wave approximation with a weighted residuals technique. The model includes (second-order) viscous dispersion effects which originate from the streamwise momentum equation and tangential stress balance. These effects play a dispersive role that primarily influences the shape of the capillary ripples in front of the solitary pulses. We show that different physical parameters, such as surface tension and viscosity, play a crucial role in the interaction between solitary pulses giving rise eventually to the formation of bound states consisting of two or more pulses separated by well-defined distances and travelling at the same velocity. By developing a rigorous coherent-structure theory, we are able to theoretically predict the pulse-separation distances for which bound states are formed. Viscous dispersion affects the distances at which bound states are observed. We show that the theory is in very good agreement with computations of the second-order model. We also demonstrate that the presence of bound states allows the film free surface to reach a self-organized state that can be statistically described in terms of a gas of solitary waves separated by a typical mean distance and characterized by a typical density. △ Less

Submitted 1 October, 2013; v1 submitted 4 November, 2010; originally announced November 2010.

Journal ref: Phys. Fluids 23, 044104 (2011)

Equilibrium gas-liquid-solid contact angle from density-functional theory

Authors: Antonio Pereira, Serafim Kalliadasis

Abstract: We investigate the equilibrium of a fluid in contact with a solid boundary through a density-functional theory. Depending on the conditions, the fluid can be in one phase, gas or liquid, or two phases, while the wall induces an external field acting on the fluid particles. We first examine the case of a liquid film in contact with the wall. We construct bifurcation diagrams for the film thickness… ▽ More We investigate the equilibrium of a fluid in contact with a solid boundary through a density-functional theory. Depending on the conditions, the fluid can be in one phase, gas or liquid, or two phases, while the wall induces an external field acting on the fluid particles. We first examine the case of a liquid film in contact with the wall. We construct bifurcation diagrams for the film thickness as a function of the chemical potential. At a specific value of the chemical potential, two equally stable films, a thin one and a thick one, can coexist. As saturation is approached, the thickness of the thick film tends to infinity. This allows the construction of a liquid-gas interface that forms a well defined contact angle with the wall. △ Less

Submitted 11 March, 2010; originally announced March 2010.

Comments: 35 pages, 20 figures

Journal ref: J. Fluid Mech. vol. 692, pgs 53-77 (2012)

Two-dimensional droplet spreading over random topographical substrates

Authors: Nikos Savva, Serafim Kalliadasis, Grigorios A. Pavliotis

Abstract: We examine theoretically the effects of random topographical substrates on the motion of two-dimensional droplets via appropriate statistical approaches. Different random substrate families are represented as stationary random functions. The variance of the droplet shift at both early times and in the long-time limit is deduced and the droplet footprint is found to be a normal random variable at… ▽ More We examine theoretically the effects of random topographical substrates on the motion of two-dimensional droplets via appropriate statistical approaches. Different random substrate families are represented as stationary random functions. The variance of the droplet shift at both early times and in the long-time limit is deduced and the droplet footprint is found to be a normal random variable at all times. It is shown that substrate roughness decreases droplet wetting, illustrating also the tendency of the droplet to slide without spreading as equilibrium is approached. Our theoretical predictions are verified by numerical experiments. △ Less

Submitted 29 December, 2009; originally announced December 2009.

Comments: 12 pages, 5 figures

The structure of flame filaments in chaotic flows

Authors: Istvan Zoltan Kiss, J. H. Merkin, S. K. Scott, P. L. Simon, S. Kalliadasis, Zoltan Neufeld

Abstract: The structure of flame filaments resulting from chaotic mixing within a combustion reaction is considered. The transverse profile of the filaments is investigated numerically and analytically based on a one-dimensional model that represents the effect of stirring as a convergent flow. The dependence of the steady solutions on the Damkohler number and Lewis number is treated in detail. It is foun… ▽ More The structure of flame filaments resulting from chaotic mixing within a combustion reaction is considered. The transverse profile of the filaments is investigated numerically and analytically based on a one-dimensional model that represents the effect of stirring as a convergent flow. The dependence of the steady solutions on the Damkohler number and Lewis number is treated in detail. It is found that, below a critical Damkohler number Da(crit), the flame is quenched by the flow. The quenching transition appears as a result of a saddle-node bifurcation where the stable steady filament solution collides with an unstable one. The shape of the steady solutions for the concentration and temperature profiles changes with the Lewis number and the value of Da(crit) increases monotonically with the Lewis number. Properties of the solutions are studied analytically in the limit of large Damkohler number and for small and large Lewis number. △ Less

Submitted 31 October, 2002; originally announced October 2002.

Comments: 17 pages, 13 figures, to be published in Physica D