title: Parallel Strands: A Preliminary Investigation into Mining the Web for Bilingual Text

abstract: Parallel corpora are a valuable resource for machine translation, but at present their availability and utility is limited by genre- and domain-specificity, licensing restrictions, and the basic difficulty of locating parallel texts in all but the most dominant of the world's languages. A parallel corpus resource not yet explored is the World Wide Web, which hosts an abundance of pages in parallel translation, offering a potential solution to some of these problems and unique opportunities of its own. This paper presents the necessary first step in that exploration: a method for automatically finding parallel translated documents on the Web. The technique is conceptually simple, fully language independent, and scalable, and preliminary evaluation results indicate that the method may be accurate enough to apply without human intervention.


title: Word Length Frequency and Distribution in English: Observations, Theory, and Implications for the Construction of Verse Lines

abstract: Recent observations in the theory of verse and empirical metrics have suggested that constructing a verse line involves a pattern-matching search through a source text, and that the number of found elements (complete words totaling a specified number of syllables) is given by dividing the total number of words by the mean number of syllables per word in the source text. This paper makes this latter point explicit mathematically, and in the course of this demonstration shows that the word length frequency totals in English output are distributed geometrically (previous researchers reported an adjusted Poisson distribution), and that the sequential distribution is random at the global level, with significant non-randomness in the fine structure. Data from a corpus of just under two million words, and a syllable-count lexicon of 71,000 word-forms is reported. The pattern-matching theory is shown to be internally coherent, and it is observed that some of the analytic techniques described here form a satisfactory test for regular (isometric) lineation in a text.


title: Combining Expression and Content in Domains for Dialog Managers

abstract: We present work in progress on abstracting dialog managers from their domain in order to implement a dialog manager development tool which takes (among other data) a domain description as input and delivers a new dialog manager for the described domain as output. Thereby we will focus on two topics; firstly, the construction of domain descriptions with description logics and secondly, the interpretation of utterances in a given domain.


title: Isometric Lineation in English Texts: An Empirical and Mathematical Examination of its Character and Consequences

abstract: In this paper we build on earlier observations and theory regarding word length frequency and sequential distribution to develop a mathematical characterization of some of the language features distinguishing isometrically lineated text from unlineated text, in other words the features distinguishing isometrical verse from prose. It is shown that the frequency of syllables making complete words produces a flat distribution for prose, while that for verse exhibits peaks at the line length position and subsequent multiples of that position. Data from several verse authors is presented, including a detailed mathematical analysis of the dynamics underlying peak creation, and comments are offered on the processes by which authors construct lines. We note that the word-length sequence of prose is random, whereas lineation necessitates non-random word-length sequencing, and that this has the probable consequence of introducing a degree of randomness into the otherwise highly ordered grammatical sequence. In addition we observe that this effect can be ameliorated by a reduction in the mean word length of the text (confirming earlier observations that verse tends to use shorter words) and the use of lines varying from the core isometrical set. The frequency of variant lines is shown to be coincident with the frequency of polysyllables, suggesting that the use of variant lines is motivated by polysyllabic word placement. The restrictive effects of different line lengths, the relationship between metrical restriction and poetic effect, and the general character of metrical rules are also discussed.


title: Some Properties of Preposition and Subordinate Conjunction Attachments

abstract: Determining the attachments of prepositions and subordinate conjunctions is a key problem in parsing natural language. This paper presents a trainable approach to making these attachments through transformation sequences and error-driven learning. Our approach is broad coverage, and accounts for roughly three times the attachment cases that have previously been handled by corpus-based techniques. In addition, our approach is based on a simplified model of syntax that is more consistent with the practice in current state-of-the-art language processing systems. This paper sketches syntactic and algorithmic details, and presents experimental results on data sets derived from the Penn Treebank. We obtain an attachment accuracy of 75.4% for the general case, the first such corpus-based result to be reported. For the restricted cases previously studied with corpus-based methods, our approach yields an accuracy comparable to current work (83.1%).


title: Deriving the Predicate-Argument Structure for a Free Word Order Language

abstract: In relatively free word order languages, grammatical functions are intricately related to case marking. Assuming an ordered representation of the predicate-argument structure, this work proposes a Combinatory Categorial Grammar formulation of relating surface case cues to categories and types for correctly placing the arguments in the predicate-argument structure. This is achieved by assigning case markers GF-encoding categories. Unlike other CG formulations, type shifting does not proliferate or cause spurious ambiguity. Categories of all argument-encoding grammatical functions follow from the same principle of category assignment. Normal order evaluation of the combinatory form reveals the predicate-argument structure. Application of the method to Turkish is shown.


title: How to define a context-free backbone for DGs: Implementing a DG in the LFG formalism

abstract: This paper presents a multidimensional Dependency Grammar (DG), which decouples the dependency tree from word order, such that surface ordering is not determined by traversing the dependency tree. We develop the notion of a \emph{word order domain structure}, which is linked but structurally dissimilar to the syntactic dependency tree. We then discuss the implementation of such a DG using constructs from a unification-based phrase-structure approach, namely Lexical-Functional Grammar (LFG). Particular attention is given to the analysis of discontinuities in DG in terms of LFG's functional uncertainty.


title: Letter to Sound Rules for Accented Lexicon Compression

abstract: This paper presents trainable methods for generating letter to sound rules from a given lexicon for use in pronouncing out-of-vocabulary words and as a method for lexicon compression. As the relationship between a string of letters and a string of phonemes representing its pronunciation for many languages is not trivial, we discuss two alignment procedures, one fully automatic and one hand-seeded which produce reasonable alignments of letters to phones. Top Down Induction Tree models are trained on the aligned entries. We show how combined phoneme/stress prediction is better than separate prediction processes, and still better when including in the model the last phonemes transcribed and part of speech information. For the lexicons we have tested, our models have a word accuracy (including stress) of 78% for OALD, 62% for CMU and 94% for BRULEX. The extremely high scores on the training sets allow substantial size reductions (more than 1/20). WWW site: http://tcts.fpms.ac.be/synthesis/mbrdico


title: Primitive Part-of-Speech Tagging using Word Length and Sentential Structure

abstract: It has been argued that, when learning a first language, babies use a series of small clues to aid recognition and comprehension, and that one of these clues is word length. In this paper we present a statistical part of speech tagger which trains itself solely on the number of letters in each word in a sentence.


title: Separating Surface Order and Syntactic Relations in a Dependency Grammar

abstract: This paper proposes decoupling the dependency tree from word order, such that surface ordering is not determined by traversing the dependency tree. We develop the notion of a \emph{word order domain structure}, which is linked but structurally dissimilar to the syntactic dependency tree. The proposal results in a lexicalized, declarative, and formally precise description of word order; features which lack previous proposals for dependency grammars. Contrary to other lexicalized approaches to word order, our proposal does not require lexical ambiguities for ordering alternatives.


title: Partial Evaluation for Efficient Access to Inheritance Lexicons

abstract: Multiple default inheritance formalisms for lexicons have attracted much interest in recent years. I propose a new efficient method to access such lexicons. After showing two basic strategies for lookup in inheritance lexicons, a compromise is developed which combines to a large degree (from a practical point of view) the advantages of both strategies and avoids their disadvantages. The method is a kind of (off-line) partial evaluation that makes a subset of inherited information explicit before using the lexicon. I identify the parts of a lexicon which should be evaluated, and show how partial evaluation works for inheritance lexicons. Finally, the theoretical results are confirmed by a complete implementation. Speedups by a factor of 10-100 are reached.


title: Spotting Prosodic Boundaries in Continuous Speech in French

abstract: A radio speech corpus of 9mn has been prosodically marked by a phonetician expert, and non expert listeners. this corpus is large enough to train and test an automatic boundary spotting system, namely a time delay neural network fed with F0 values, vowels and pseudo-syllable durations. Results validate both prosodic marking and automatic spotting of prosodic events.


title: Error-Driven Pruning of Treebank Grammars for Base Noun Phrase Identification

abstract: Finding simple, non-recursive, base noun phrases is an important subtask for many natural language processing applications. While previous empirical methods for base NP identification have been rather complex, this paper instead proposes a very simple algorithm that is tailored to the relative simplicity of the task. In particular, we present a corpus-based approach for finding base NPs by matching part-of-speech tag sequences. The training phase of the algorithm is based on two successful techniques: first the base NP grammar is read from a ``treebank'' corpus; then the grammar is improved by selecting rules with high ``benefit'' scores. Using this simple algorithm with a naive heuristic for matching rules, we achieve surprising accuracy in an evaluation on the Penn Treebank Wall Street Journal.


title: Segregatory Coordination and Ellipsis in Text Generation

abstract: In this paper, we provide an account of how to generate sentences with coordination constructions from clause-sized semantic representations. An algorithm is developed to generate sentences with ellipsis, gapping, right-node-raising, and non-constituent coordination constructions. Various examples from linguistic literature will be used to demonstrate that the algorithm does its job well.


title: A Variant of Earley Parsing

abstract: The Earley algorithm is a widely used parsing method in natural language processing applications. We introduce a variant of Earley parsing that is based on a ``delayed'' recognition of constituents. This allows us to start the recognition of a constituent only in cases in which all of its subconstituents have been found within the input string. This is particularly advantageous in several cases in which partial analysis of a constituent cannot be completed and in general in all cases of productions sharing some suffix of their right-hand sides (even for different left-hand side nonterminals). Although the two algorithms result in the same asymptotic time and space complexity, from a practical perspective our algorithm improves the time and space requirements of the original method, as shown by reported experimental results.


title: Towards an implementable dependency grammar

abstract: The aim of this paper is to define a dependency grammar framework which is both linguistically motivated and computationally parsable. See the demo at http://www.conexor.fi/analysers.html#testing


title: Some Ontological Principles for Designing Upper Level Lexical Resources

abstract: The purpose of this paper is to explore some semantic problems related to the use of linguistic ontologies in information systems, and to suggest some organizing principles aimed to solve such problems. The taxonomic structure of current ontologies is unfortunately quite complicated and hard to understand, especially for what concerns the upper levels. I will focus here on the problem of ISA overloading, which I believe is the main responsible of these difficulties. To this purpose, I will carefully analyze the ontological nature of the categories used in current upper-level structures, considering the necessity of splitting them according to more subtle distinctions or the opportunity of excluding them because of their limited organizational role.


title: A Comparison of WordNet and Roget's Taxonomy for Measuring Semantic Similarity

abstract: This paper presents the results of using Roget's International Thesaurus as the taxonomy in a semantic similarity measurement task. Four similarity metrics were taken from the literature and applied to Roget's The experimental evaluation suggests that the traditional edge counting approach does surprisingly well (a correlation of r=0.88 with a benchmark set of human similarity judgements, with an upper bound of r=0.90 for human subjects performing the same task.)


title: Some comments on the correlation dimension of $1/f^\alpha$ noise

abstract: It has recently been observed that a stochastic (infinite degree of freedom) time series with a $1/f^\alpha$ power spectrum can exhibit a finite correlation dimension, even for arbitrarily large data sets. [A.R. Osborne and A.~Provenzale, {\sl Physica D} {\bf 35}, 357 (1989).] I will discuss the relevance of this observation to the practical estimation of dimension from a time series, and in particular I will argue that a good dimension algorithm need not be trapped by this anomalous fractal scaling. Further, I will analytically treat the case of gaussian \onefas noise, with explicit high and low frequency cutoffs, and derive the scaling of the correlation integral $C(N,r)$ in various regimes of the $(N,r)$ plane. Appears in: {\sl Phys. Lett. A} {\bf 155} (1991) 480--493.


title: Statistical error in a chord estimator of correlation dimension: the ``rule of five''

abstract: The statistical precision of a chord method for estimating fractal dimension from a correlation integral is derived. The optimal chord length is determined, and a comparison is made to other estimators. These calculations use the approximation that all pairwise distances between the points are statistically independent; the adequacy of this approximation is assessed numerically. The chord method provides a very quick and easy dimension estimate which is only slightly less precise than the optimal estimator. Keywords: correlation dimension, statistical error


title: Detecting Nonlinearity in Data with Long Coherence Times

abstract: We consider the limitations of two techniques for detecting nonlinearity in time series. The first technique compares the original time series to an ensemble of surrogate time series that are constructed to mimic the linear properties of the original. The second technique compares the forecasting error of linear and nonlinear predictors. Both techniques are found to be problematic when the data has a long coherence time; they tend to indicate nonlinearity even for linear time series. We investigate the causes of these difficulties both analytically and with numerical experiments on ``real'' and computer-generated data. In particular, although we do see some initial evidence for nonlinear structure in the SFI dataset E, we are inclined to dismiss this evidence as an artifact of the long coherence time.


title: A Lattice Boltzmann Model for Multi-phase Fluid Flows

abstract: We develop a lattice Boltzmann equation method for simulating multi-phase immiscible fluid flows with variation of density and viscosity, based on the model proposed by Gunstensen {\em et al} for two-component immiscible fluids. The numerical measurements of surface tension and viscosity agree well with theoretical predictions. Several basic numerical tests, including spinodal decomposition, two-phase fluid flows in two-dimensional channels and two-phase viscous fingering, are shown in agreement of experiments and analytical solutions.


title: Message-Passing Multi-Cell Molecular Dynamics on the Connection Machine 5

abstract: We present a new scalable algorithm for short-range molecular dynamics simulations on distributed memory MIMD multicomputer based on a message-passing multi-cell approach. We have implemented the algorithm on the Connection Machine 5 (CM-5) and demonstrate that meso-scale molecular dynamics with more than $10^8$ particles is now possible on massively parallel MIMD computers. Typical runs show single particle update-times of $0.15 \mu s$ in 2 dimensions (2D) and approximately $1 \mu s$ in 3 dimensions (3D) on a 1024 node CM-5 without vector units, corresponding to more than 1.8 GFlops overall performance. We also present a scaling equation which agrees well with actually observed timings.


title: Comparison of Spectral Method and Lattice Boltzmann Simulations of Two-Dimensional Hydrodynamics

abstract: We present numerical solutions of the two-dimensional Navier-Stokes equations by two methods; spectral and the novel Lattice Boltzmann Equation (LBE) scheme. Very good agreement is found for global quantities as well as energy spectra. The LBE scheme is, indeed, providing reasonably accurate solutions of the Navier-Stokes equations with an isothermal equation of state, in the nearly incompressible limit. Relaxation to a previously reported ``sinh-Poisson'' state is also observed for both runs.


title: Domain Growth, Wetting and Scaling in Porous Media

abstract: The lattice Boltzmann (LB) method is used to study the kinetics of domain growth of a binary fluid in a number of geometries modeling porous media. Unlike the traditional methods which solve the Cahn-Hilliard equation, the LB method correctly simulates fluid properties, phase segregation, interface dynamics and wetting. Our results, based on lattice sizes of up to $4096\times 4096$, do not show evidence to indicate the breakdown of late stage dynamical scaling, and suggest that confinement of the fluid is the key to the slow kinetics observed. Randomness of the pore structure appears unnecessary.


title: Lattice Boltzmann Equation for Quantum Mechanics

abstract: It is shown that the Lattice Boltzmann equation for hydrodynamics can be extended in such a way as to describe non-relativistic quantum mechanics.


title: The Lattice Boltzmann Equation Method for the Simulation of Compressible Fluid Flow

abstract: We systematically derived hydrodynamic equations and transport coefficients for a class of multi-speed lattice Boltzmann models in D dimensions, using the multi-scale technique. The constitutive relation of physical fluid is recovered by a modified equilibrium distribution in Maxwell-Boltzmann type. With the use of the rest particles and the particle reservoir, we were able to add one degree of freedom into the sound speed of the modeled fluid. When the sound speed is tuned small enough, the compressible region of fluid flow can be reached. An example 2-D model is presented, together with the numerical verification for its transport coefficients.


title: Wavelets and Fast Numerical Algorithms

abstract: Wavelet based algorithms in numerical analysis are similar to other transform methods in that vectors and operators are expanded into a basis and the computations take place in this new system of coordinates. However, due to the recursive definition of wavelets, their controllable localization in both space and wave number (time and frequency) domains, and the vanishing moments property, wavelet based algorithms exhibit new and important properties. For example, the multiresolution structure of the wavelet expansions brings about an efficient organization of transformations on a given scale and of interactions between different neighbouring scales. Moreover, wide classes of operators which naively would require a full (dense) matrix for their numerical description, have sparse representations in wavelet bases. For these operators sparse representations lead to fast numerical algorithms, and thus address a critical numerical issue. We note that wavelet based algorithms provide a systematic generalization of the Fast Multipole Method (FMM) and its descendents. These topics will be the subject of the lecture. Starting from the notion of multiresolution analysis, we will consider the so-called non-standard form (which achieves decoupling among the scales) and the associated fast numerical algorithms. Examples of non-standard forms of several basic operators (e.g. derivatives) will be computed explicitly.


title: Hydrodynamic Spinodal Decomposition: Growth Kinetics and Scaling Functions

abstract: We examine the effects of hydrodynamics on the late stage kinetics in spinodal decomposition. From computer simulations of a lattice Boltzmann scheme we observe, for critical quenches, that single phase domains grow asymptotically like $t^{\alpha}$, with $\alpha \approx .66$ in two dimensions and $\alpha \approx 1.0$ in three dimensions, both in excellent agreement with theoretical predictions.


title: Lattice Boltzmann Thermohydrodynamics

abstract: We introduce a lattice Boltzmann computational scheme capable of modeling thermohydrodynamic flows of monatomic gases. The parallel nature of this approach provides a numerically efficient alternative to traditional methods of computational fluid dynamics. The scheme uses a small number of discrete velocity states and a linear, single-time-relaxation collision operator. Numerical simulations in two dimensions agree well with exact solutions for adiabatic sound propagation and Couette flow with heat transfer.


title: Don't bleach chaotic data

abstract: A common first step in time series signal analysis involves digitally filtering the data to remove linear correlations. The residual data is spectrally white (it is ``bleached''), but in principle retains the nonlinear structure of the original time series. It is well known that simple linear autocorrelation can give rise to spurious results in algorithms for estimating nonlinear invariants, such as fractal dimension and Lyapunov exponents. In theory, bleached data avoids these pitfalls. But in practice, bleaching obscures the underlying deterministic structure of a low-dimensional chaotic process. This appears to be a property of the chaos itself, since nonchaotic data are not similarly affected. The adverse effects of bleaching are demonstrated in a series of numerical experiments on known chaotic data. Some theoretical aspects are also discussed.


title: Identification of Continuous-Time Dynamical Systems: Neural Network Based Algorithms and Parallel Implementation

abstract: Time-delay mappings constructed using neural networks have proven successful in performing nonlinear system identification; however, because of their discrete nature, their use in bifurcation analysis of continuous-time systems is limited. This shortcoming can be avoided by embedding the neural networks in a training algorithm that mimics a numerical integrator. Both explicit and implicit integrators can be used. The former case is based on repeated evaluations of the network in a feedforward implementation; the latter relies on a recurrent network implementation. Here the algorithms and their implementation on parallel machines (SIMD and MIMD architectures) are discussed.


title: Convergence of Convective-Diffusive Lattice Boltzmann Methods

abstract: Lattice Boltzmann methods are numerical schemes derived as a kinetic approximation of an underlying lattice gas. A numerical convergence theory for nonlinear convective-diffusive lattice Boltzmann methods is established. Convergence, consistency, and stability are defined through truncated Hilbert expansions. In this setting it is shown that consistency and stability imply convergence. Monotone lattice Boltzmann methods are defined and shown to be stable, hence convergent when consistent. Examples of diffusive and convective-diffusive lattice Boltzmann methods that are both consistent and monotone are presented.


title: The non-linear response of the magnetosphere: 30 October 1978

abstract: Previous efforts to find evidence of deterministic nonlinear dynamics in the global geomagnetic system have treated the geomagnetic system as autonomous. However, the geomagnetic system is strongly driven by the stochastic solar wind. We consider the response of the magnetosphere, as given by the AE index, for one day when the IMF had a nearly constant southward value. Using both a series of non-linear statistics and non-linear prediction of the response to the input signal $v B_s$, we find that there is some evidence for deterministic non-linear response of the Earth's magnetosphere on that day.


title: Global Bifurcations in Rayleigh-Benard Convection: Experiments, Empirical Maps and Numerical Bifurcation Analysis

abstract: We use nonlinear signal processing techniques, based on artificial neural networks, to construct an empirical mapping from experimental Rayleigh-Benard convection data in the quasiperiodic regime. The data, in the form of a one-parameter sequence of Poincare sections in the interior of a mode-locked region (resonance horn), are indicative of a complicated interplay of local and global bifurcations with respect to the experimentally varied Rayleigh number. The dynamic phenomena apparent in the data include period doublings, complex intermittent behavior, secondary Hopf bifurcations, and chaotic dynamics. We use the fitted map to reconstruct the experimental dynamics and to explore the associated local and global bifurcation structures in phase space. Using numerical bifurcation techniques we locate the stable and unstable periodic solutions, calculate eigenvalues, approximate invariant manifolds of saddle type solutions and identify bifurcation points. This approach constitutes a promising data post-processing procedure for investigating phase space and parameter space of real experimental systems; it allows us to infer phase space structures which the experiments can only probe with limited measurement precision and only at a discrete number of operating parameter settings.


title: A New Class of Cellular Automata for Reaction-Diffusion Systems

abstract: We introduce a new class of cellular automata to model reaction-diffusion systems in a quantitatively correct way. The construction of the CA from the reaction-diffusion equation relies on a moving average procedure to implement diffusion, and a probabilistic table-lookup for the reactive part. The applicability of the new CA is demonstrated using the Ginzburg-Landau equation.


title: Stability Analysis of Lattice Boltzmann Methods

abstract: The lattice Boltzmann equation describes the evolution of the velocity distribution function on a lattice in a manner that macroscopic fluid dynamical behavior is recovered. Although the equation is a derivative of lattice gas automata, it may be interpreted as a Lagrangian finite-difference method for the numerical simulation of the discrete-velocity Boltzmann equation that makes use of a BGK collision operator. As a result, it is not surprising that numerical instability of lattice Boltzmann methods have been frequently encountered by researchers. We present an analysis of the stability of perturbations of the particle populations linearized about equilibrium values corresponding to a constant-density uniform mean flow. The linear stability depends on the following parameters: the distribution of the mass at a site between the different discrete speeds, the BGK relaxation time, the mean velocity, and the wavenumber of the perturbations. This parameter space is too large to compute the complete stability characteristics. We report some stability results for a


title: Initial and Boundary Conditions for the Lattice Boltzmann Method

abstract: A new approach of implementing initial and boundary conditions for the lattice Boltzmann method is presented. The new approach is based on an extended collision operator that uses the gradients of the fluid velocity. The numerical performance of the lattice Boltzmann method is tested on several problems with exact solutions and is also compared to an explicit finite difference projection method. The discretization error of the lattice Boltzmann method decreases quadratically with finer resolution both in space and in time. The roundoff error of the lattice Boltzmann method creates problems unless double precision arithmetic is used.


title: Exact Results for the Asymmetric Simple Exclusion Process with a Blockage

abstract: We present new results for the current as a function of transmission rate in the one dimensional totally asymmetric simple exclusion process (TASEP) with a blockage that lowers the jump rate at one site from one to r < 1. Exact finite volume results serve to bound the allowed values for the current in the infinite system. This proves the existence of a gap in allowed density corresponding to a nonequilibrium ``phase transition'' in the infinite system. A series expansion in r, derived from the finite systems, is proven to be asymptotic for all sufficiently large systems. Pade approximants based on this series, which make specific assumptions about the nature of the singularity at r = 1, match numerical data for the ``infinite'' system to a part in 10^4.


title: Numerical Simulations of Particulate Suspensions via a Discretized Boltzmann Equation Part I. Theoretical Foundation

abstract: A new and very general technique for simulating solid-fluid suspensions is described; its most important feature is that the computational cost scales linearly with the number of particles. The method combines Newtonian dynamics of the solid particles with a discretized Boltzmann equation for the fluid phase; the many-body hydrodynamic interactions are fully accounted for, both in the creeping-flow regime and at higher Reynolds numbers. Brownian motion of the solid particles arises spontaneously from stochastic fluctuations in the fluid stress tensor, rather than from random forces or displacements applied directly to the particles. In this paper, the theoretical foundations of the technique are laid out, illustrated by simple analytical and numerical examples; in the companion paper, extensive numerical tests of the method, for stationary flows, time-dependent flows, and finite Reynolds number flows, are reported.


title: Numerical Simulations of Particulate Suspensions via a Discretized Boltzmann Equation Part II. Numerical Results

abstract: A new and very general technique for simulating solid-fluid suspensions has been described in a previous paper (Part I); the most important feature of the new method is that the computational cost scales with the number of particles. In this paper (Part II), extensive numerical tests of the method are described; for creeping flows, both with and without Brownian motion, and at finite Reynolds numbers. Hydrodynamic interactions, transport coefficients, and the short-time dynamics of random dispersions of up to 1024 colloidal particles have been simulated.


title: Numerical Solution of the Schroedinger Equation using a Quantum Lattice Boltzmann Equation}

abstract: The quantum Lattice Boltzmann equation (QLBe), a new variant of the lattice Boltzmann equation, specifically designed to describe non relativistic quantum motion, is validated for the case of a free-particle in (1+1) space-time dimensions. The proper parameter regime under which the method needs to be operated in order to reproduce faithful non-relativistic quantum motion is also discussed.


title: A Vectorized Algorithm for Molecular Dynamics of Short Range Interacting Particles

abstract: We report on a lattice based algorithm, completely vectorized for molecular dynamics simulations. Its algorithmic complexity is of the order $O(N)$, where $N$ is the number of particles. The algorithm works very effectively when the particles have short range interaction, but it is applicable to each kind of interaction. The code was tested on a Cray ymp el in a simulation of flowing granular material.


title: Accuracy of Discrete-Velocity BGK Models for the Simulation of the Incompressible Navier-Stokes Equations

abstract: Two discretizations of a 9-velocity Boltzmann equation with a BGK collision operator are studied. A Chapman-Enskog expansion of the PDE system predicts that the macroscopic behavior corresponds to the incompressible Navier-Stokes equations with additional terms of order Mach number squared. We introduce a fourth-order scheme and compare results with those of the commonly used lattice Boltzmann discretization and with finite-difference schemes applied to the incompressible Navier-Stokes equations in primitive-variable form. We numerically demonstrate convergence of the BGK schemes to the incompressible Navier-Stokes equations and quantify the errors associated with compressibility and discretization effects. When compressibility error is smaller than discretization error, convergence in both grid spacing and time step is shown to be second-order for the LB method and is confirmed to be fourth-order for the fourth-order BGK solver. However, when the compressibility error is simultaneously reduced as the grid is refined, the LB method behaves as a first-order scheme in time.


title: Abstracts on Pattern Formation and Lattice-Gas Automata

abstract: This article contains 39 pages of abstracts of invited and poster presentations for the FIELDS INSTITUTE FOR RESEARCH IN MATHEMATICAL SCIENCES AND NATO ADVANCED RESEARCH WORKSHOP PROGRAM ON PATTERN FORMATION and LATTICE-GAS AUTOMATA, JUNE 07-12, 1993


title: Hydrodynamic behaviour of Lattice Boltzmann and Lattice BGK models

abstract: We present a numerical analysis of the validity of classical and generalized hydrodynamics for Lattice Boltzmann Equation (LBE) and Lattice BGK methods in two and three dimensions, as a function of the collision parameters of these models. Our analysis is based on the wave-number dependence of the evolution operator. Good ranges of validity are found for BGK models as long as the relaxation time is chosen smaller than or equal to unity. The additional freedom in the choice of collision parameters for LBE models does not seem to give significant improvement.


title: Lattice Boltzmann-Langevin Equations

abstract: Intrinsic fluctuations around the solution of the lattice Boltzmann equation are described or modeled by addition of a white Gaussian noise source. For stationary states a fluctuation-dissipation theorem relates the variance of the fluctuations to the linearized Boltzmann collision operator and the pair correlation function.


title: Instabilities and Patterns (minor technical modifications)

abstract: Violation of (semi)-detailed balance conditions in lattice gas automata gives rise to unstable spatial fluctuations that lead to phase separation and pattern formation in spinodal decomposition, unstable propagating modes, driven diffusive systems and unstable uniform flows.


title: Elimination of Nonlinear Deviations in Thermal Lattice BGK Models

abstract: Abstracet: We present a new thermal lattice BGK model in D-dimensional space for the numerical calculation of fluid dynamics. This model uses a higher order expansion of equilibrium distribution in Maxwellian type. In the mean time the lattice symmetry is upgraded to ensure the isotropy of 6th order tensor. These manipulations lead to macroscopic equations free from nonlinear deviations. We demonstrate the improvements by conducting classical Chapman-Enskog analysis and by numerical simulation of shear wave flow. The transport coefficients are measured numerically, too.


title: Particle-Mesh Methods on the Connection Machine

abstract: We describe an efficient Particle-Mesh algorithm for the Connection Machine CM-5. Our particular method parallelizes well and the computation time per time step decreases as the particles become more clustered. We achieve floating-point computation rates of 4--5 MFlops/sec/processing node and total operations (the sum of floating-point and integer arithmetic plus communications) of 5--10 MOps/sec/processing node. The rates scale almost linearly from 32 to 256 processors. Although some of what we discuss is specific to the CM-5, many aspects (e.g., the computation of the force on a mesh) are generic to all implementations, and other aspects (e.g., the algorithm for assignment of the density to the mesh) are useful on any parallel computer.


title: Lattice Methods and Their Applications To Reacting Systems

abstract: The recent development of the lattice gas automata method and its extension to the lattice Boltzmann method have provided new computational schemes for solving a variety of partial differential equations and modeling chemically reacting systems. The lattice gas method, regarded as the simplest microscopic and kinetic approach which generates meaningful macroscopic dynamics, is fully parallel and can, as a result, be easily programmed on parallel machines. In this paper, we introduce the basic principles of the lattice gas method and the lattice Boltzmann method, their numerical implementations and applications to chemically reacting systems. Comparisons of the lattice Boltzmann method with the lattice gas technique and other traditional numerical schemes, including the finite difference scheme and the pseudo-spectral method, for solving the Navier-Stokes hydrodynamic fluid flows will be discussed. Recent developments of the lattice gas and the lattice Boltzmann method and their applications to pattern formation in chemical reaction-diffusion systems, multiphase fluid flows and polymeric dynamics will be presented.


title: Cellular automaton model of reaction-transport processes

abstract: The transport and chemical reactions of solutes are modelled as a cellular automaton in which molecules of different species perform a random walk on a regular lattice and react according to a local probabilistic rule. The model describes advection and diffusion in a simple way, and as no restriction is placed on the number of particles at a lattice site, it is also able to describe a wide variety of chemical reactions. Assuming molecular chaos and a smooth density function, we obtain the standard reaction-transport equations in the continuum limit. Simulations on one- and two-dimensional lattices show that the discrete model can be used to approximate the solutions of the continuum equations. We discuss discrepancies which arise from correlations between molecules and how these disappear as the continuum limit is approached. Of particular interest are simulations displaying long-time behaviour which depends on long-wavelength statistical fluctuations not accounted for by the standard equations. The model is applied to the reactions $a + b \rightleftharpoons c$ and $a + b \rightarrow c$ with homogeneous and inhomogeneous initial conditions as well as to systems subject to autocatalytic reactions and displaying spontaneous formation of spatial concentration patterns.


title: Simulation of non-ideal gases and liquid-gas phase transitions by lattice Boltzmann equation

abstract: We describe in detail a recently proposed lattice-Boltzmann model for simulating flows with multiple phases and components. In particular, the focus is on the modeling of one-component fluid systems which obey non-ideal gas equations of state and can undergo a liquid-gas type phase transition. The model is shown to be momentum-conserving. From the microscopic mechanical stability condition, the densities in bulk liquid and gas phases are obtained as functions of a temperature-like parameter. Comparisons with the thermodynamic theory of phase transition show that the LBE model can be made to correspond exactly to an isothermal process. The density profile in the liquid-gas interface is also obtained as function of the temperature-like parameter and is shown to be isotropic. The surface tension, which can be changed independently, is calculated. The analytical conclusions are verified by numerical simulations. (To appear in Phys. Rev. E)


title: Lattice Boltzmann Magnetohydrodynamics

abstract: Lattice gas and lattice Boltzmann methods are recently developed numerical schemes for simulating a variety of physical systems. In this paper a new lattice Boltzmann model for modeling two-dimensional incompressible magnetohydrodynamics (MHD) is presented. The current model fully utilizes the flexibility of the lattice Boltzmann method in comparison with previous lattice gas and lattice Boltzmann MHD models, reducing the number of moving directions from $36$ in other models to $12$ only. To increase computational efficiency, a simple single time relaxation rule is used for collisions, which directly controls the transport coefficients. The bi-directional streaming process of the particle distribution function in this paper is similar to the original model [ H. Chen and W. H. Matthaeus, Phys. Rev. Lett., {\bf 58}, 1845(1987), S.Chen, H.Chen, D.Mart\'{\i}nez and W.H.Matthaeus, Phys. Rev. Lett. {\bf 67},3776 (1991)], but has been greatly simplified, affording simpler implementation of boundary conditions and increasing the feasibility of extension into a workable three-dimensional model. Analytical expressions for the transport coefficients are presented. Also, as example cases, numerical calculation for the Hartmann flow is performed, showing a good agreement between the theoretical


title: Simulation of Cavity Flow by the Lattice Boltzmann Method

abstract: A detailed analysis is presented to demonstrate the capabilities of the lattice Boltzmann method. Thorough comparisons with other numerical solutions for the two-dimensional, driven cavity flow show that the lattice Boltzmann method gives accurate results over a wide range of Reynolds numbers. Studies of errors and convergence rates are carried out. Compressibility effects are quantified for different maximum velocities, and parameter ranges are found for stable simulations. The paper's objective is to stimulate further work using this relatively new approach for applied engineering problems in transport phenomena utilizing parallel computers.


title: A Lattice Boltzmann Subgrid Model for High Reynolds Number Flows

abstract: A subgrid turbulence model for the lattice Boltzmann method is proposed for high Reynolds number fluid flow applications. The method, based on the standard Smagorinsky subgrid model and a single-time relaxation lattice Boltzmann method, incorporates the advantages of the lattice Boltzmann method for handling arbitrary boundaries and is easily implemented on parallel machines. The method is applied to a two-dimensional driven cavity flow for studying dynamics and the Reynolds number dependence of the flow structures. The substitution of other subgrid models, such as the dynamic subgrid model, in the framework of the LB method is discussed.


title: Estimating Functions of Probability Distributions from a Finite Set of Samples, Part 1: Bayes Estimators and the Shannon Entropy

abstract: We present estimators for entropy and other functions of a discrete probability distribution when the data is a finite sample drawn from that probability distribution. In particular, for the case when the probability distribution is a joint distribution, we present finite sample estimators for the mutual information, covariance, and chi-squared functions of that probability distribution.


title: Estimating Functions of Distributions from A Finite Set of Samples, Part 2: Bayes Estimators for Mutual Information, Chi-Squared, Covariance and other Statistics


title: Correlations and Renormalization in Lattice Gases

abstract: A complete formulation is given of an exact kinetic theory for lattice gases. This kinetic theory makes possible the calculation of corrections to the usual Boltzmann / Chapman-Enskog analysis of lattice gases due to the buildup of correlations. It is shown that renormalized transport coefficients can be calculated perturbatively by summing terms in an infinite series. A diagrammatic notation for the terms in this series is given, in analogy with the diagrammatic expansions of continuum kinetic theory and quantum field theory. A closed-form expression for the coefficients associated with the vertices of these diagrams is given. This method is applied to several standard lattice gases, and the results are shown to correctly predict experimentally observed deviations from the Boltzmann analysis.

categories: comp-gas adap-org chao-dyn cond-mat hep-th nlin.AO nlin.CD nlin.CG nlin.PS patt-sol


title: A Lattice-Gas with Long-Range Interactions Coupled to a Heat Bath

abstract: Introduced is a lattice-gas with long-range 2-body interactions. An effective inter-particle force is mediated by momentum exchanges. There exists the possibility of having both attractive and repulsive interactions using finite impact parameter collisions. There also exists an interesting possibility of coupling these long-range interactions to a heat bath. A fixed temperature heat bath induces a permanent net attractive interparticle potential, but at the expense of reversibility. Thus the long-range dynamics is a kind of a Monte Carlo Kawasaki updating scheme. The model has a $P\rho T$ equation of state. Presented are analytical and numerical results for the a lattice-gas fluid governed by a nonideal equation of state. The model's complexity is not much beyond that of the FHP lattice-gas. It is suitable for massively parallel processing and may be used to study critical phenomena in large systems.


title: Lattice-Gas Automata Fluids on Parallel Supercomputers

abstract: A condensed history and theoretical development of lattice-gas automata in the Boltzmann limit is presented. This is provided as background to set up the context for understanding the implementation of the lattice-gas method on two parallel supercomputers: the MIT cellular automata machine CAM-8 and the Connection Machine CM-5. The macroscopic limit of two-dimensional fluids is tested by simulating the Rayleigh-B\'enard convective instability, Kelvin-Helmholtz shear instability, and the Von Karman vortex shedding instability. Performance of the two machines in terms of both site update rate and maximum problem size are comparable. The CAM-8, being a low-cost desktop machine, demonstrates the potential of special-purpose digital hardware.


title: Application of Complex Daubechies' Wavelets to Numerical Simulation of a Nonlinear Signal Propagation Model

abstract: We report the first application of complex symmetric wavelets to the numerical simulation of a nonlinear signal propagation model. This model is the so-called nonlinear Schrodinger equation that describes, for instance, the evolution of the electric field amplitude in nonlinear optical fibers. We propose and study a new way to implement a global space-time adaptive grid, based on interpolation properties of higher-order scaling functions.


title: Transforming Signs to Phase Distributions in Quantum Simulations

abstract: A method is developed which speeds up averaging in quantum simulations where minus signs cause difficulties. A Langevin equation method in conjunction with a replication algorithm is used enabling one to average over a continuously varying complex number. Instead of ensemble averaging this number directly, the phase of the complex number is followed over time. The method is illustrated in some simple cases where the answers obtained can be compared to exact results, and also compared to conventional averaging procedures which converge orders of magnitude slower than this method. Limitations of this method are also described.


title: Generalized Boltzmann Equation for Lattice Gas Automata

abstract: In this paper, for the first time a theory is formulated that predicts velocity and spatial correlations between occupation numbers that occur in lattice gas automata violating semi-detailed balance. Starting from a coupled BBGKY hierarchy for the $n$-particle distribution functions, cluster expansion techniques are used to derive approximate kinetic equations. In zeroth approximation the standard nonlinear Boltzmann equation is obtained; the next approximation yields the ring kinetic equation, similar to that for hard sphere systems, describing the time evolution of pair correlations. As a quantitative test we calculate equal time correlation functions in equilibrium for two models that violate semi-detailed balance. One is a model of interacting random walkers on a line, the other one is a two-dimensional fluid type model on a triangular lattice. The numerical predictions agree very well with computer simulations.


title: Lattice Gas Automata: a discrete model for simulation of continuous transport phenomena in packages of agricultural products

abstract: A Lattice Boltzmann scheme has been applied to the problem of heat and mass transport in packagings of cut flowers. The LB scheme, simulating convection-diffusion processes and heat and mass transfer between flowers and air flow, is described in detail. The objective of the research programme is to optimize the packaging design for improved flower quality. Comparison of experimental results and simulation results show that the LB-scheme is a promising simulation technique for achieving the research objective.


title: Generating surrogate data for time series with several simultaneously measured variables

abstract: We propose an extension to multivariate time series of the phase-randomized Fourier-transform algorithm for generating surrogate data. Such surrogate data sets must mimic not only the autocorrelations of each of the variables in the original data set, they must mimic the cross-correlations {\em between} all the variables as well. The method is applied both to a simulated example (the three components of the Lorenz equations) and to data from a multichannel electroencephalogram.


title: Growth Kinetics in Multicomponent Fluids

abstract: The hydrodynamic effects on the late stage kinetics in spinodal decomposition of multicomponent fluids are examined using a lattice Boltzmann scheme with stochastic fluctuations in the fluid and at the interface. In two dimensions, the three and four component immiscible fluid mixture (with a $1024^2$ lattice) behaves like an off-critical binary fluid with an estimated domain growth of $t^{0.4\pm.03}$ rather than $t^{1/3}$ as previously predicted, showing the significant influence of hydrodynamics. In three dimensions (with a $256^3$ lattice), we estimate the growth as $t^{0.96\pm0.05}$ for both critical and off-critical quenching, in agreement with phenomenological theory.


title: A Class of Fast Methods for Processing Irregularly Sampled or Otherwise Inhomogeneous One-Dimensional Data

abstract: With the ansatz that a data set's correlation matrix has a certain parametrized form (one general enough, however, to allow the arbitrary specification of a slowly-varying decorrelation distance and population variance) the general machinery of Wiener or optimal filtering can be reduced from $O(n^3)$ to $O(n)$ operations, where $n$ is the size of the data set. The implied vast increases in computational speed can allow many common sub-optimal or heuristic data analysis methods to be replaced by fast, relatively sophisticated, statistical algorithms. Three examples are given: data rectification, high- or low- pass filtering, and linear least squares fitting to a model with unaligned data points.


title: The correlation dimension of differenced data

abstract: The effect differencing has on the estimated correlation dimension of a nongaussian time series is discussed. Two different methods for generating surrogate data sets are compared. The results suggest that any filtering should be applied to both the raw and surrogate time series, otherwise spurious results can occur.


title: Generalized redundancies for time series analysis

abstract: Extensions to various information theoretic quantities used for nonlinear time series analysis are discussed, as well as their relationship to the generalized correlation integral. It is shown that calculating redundancies from the correlation integral can be more accurate and more efficient than direct box counting methods. It is also demonstrated that many commonly used nonlinear statistics have information theory based analogues. Furthermore, the relationship between the correlation integral and information theoretic statistics allows us to define ``local'' versions of many information theory based statistics; including a local version of the Kolmogorov-Sinai entropy, which gives an estimate of the local predictability.


title: Abstracts for the June 27-29 Princeton Conference: Discrete Models for Fluid Mechanics

abstract: Submitted by the authors for the June 27-29 Princeton Conference. Questions should be directed to: qian@cfd.princeton.edu


title: Theoretical approach to two-dimensional traffic flow models

abstract: In this paper we present a theoretical analysis of a recently proposed two-dimensional Cellular Automata model for traffic flow in cities with the novel ingredient of turning capability. Numerical simulations of this model show that there is a transition between a freely moving phase with high velocity to a jammed state with low velocity. We study the dynamics of such a model starting with the microscopic evolution equation, which will serve as a basis for further analysis. It is shown that a kinetic approach, based on the Boltzmann assumption, is able to provide a reasonably good description of the jamming transition. We further introduce a space-time continuous phenomenological model leading to a couple of partial differential equations whose preliminary results agree rather well with the numerical simulations.

categories: comp-gas adap-org cond-mat nlin.AO nlin.CG


title: Testing For Nonlinearity Using Redundancies: Quantitative and Qualitative Aspects

abstract: A method for testing nonlinearity in time series is described based on information-theoretic functionals -- redundancies, linear and nonlinear forms of which allow either qualitative, or, after incorporating the surrogate data technique, quantitative evaluation of dynamical properties of scrutinized data. An interplay of quantitative and qualitative testing on both the linear and nonlinear levels is analyzed and robustness of this combined approach against spurious nonlinearity detection is demonstrated. Evaluation of redundancies and redundancy-based statistics as functions of time lag and embedding dimension can further enhance insight into dynamics of a system under study.


title: Flow in the Driven Cavity Calculated by the Lattice Boltzmann Method

abstract: The lattice Boltzmann method with enhanced collisions and rest particles is used to calculate the flow in a two-dimensional lid-driven cavity. The abilitity of this method to compute the velocity and the pressure of an incompressible fluid in a geometry with Dirichlet and Neumann boundary conditions is verified by calculating a test-problem where the analytical solution is known. Different parameter configurations have been tested for Reynolds numbers from $Re=10$ to $Re=2000$. The vortex structure for a more generalized lid-driven cavity problem with a non-uniform top speed has been studied for various acpect ratios.


title: Cellular Automaton for Realistic Modelling of Landslides

abstract: A numerical model is developed for the simulation of debris flow in landslides over a complex three dimensional topography. The model is based on a lattice, in which debris can be transferred among nearest neighbors according to established empirical relationships for granular flows. The model is then validated by comparing a simulation with reported field data. Our model is in fact a realistic elaboration of simpler ``sandpile automata'', which have in recent years been studied as supposedly paradigmatic of ``self-organized criticality''. Statistics and scaling properties of the simulation are examined, and show that the model has an intermittent behavior.


title: Nonlinear Analysis of Experimental Noisy Time Series in Fluidized Bed Systems

abstract: The paper describes the application of some numerical techniques to analyze and to characterize the observed dynamical behaviour of fluidized bed systems. The preliminary results showed clearly that the dynamics of the considered process can be nonrecurrent and governed by chaotic deterministic rules rather than by stochastic ones. This significant conclusion was allowed by the application of a proper filtering procedure which was able to reduce the unwanted influence of significant broadband noise detected in the sampled experimental data.


title: Recognition of Temporal Sequences of Patterns

abstract: Using an asymmetric associative network with synchronous updating, it is possible to recall a sequence of patterns. To obtain a stable sequence generation with a large storage capacity, we introduce a threshold that eliminates the contribution of weakly correlated patterns. For this system we find a set of evolution equations for the overlaps of the states with the patterns to be recognized. We solve these equations in the limit of the stationary cycle, and obtain the critical value of the capacity as a function of the threshold and temperature. Finally, a numerical simulation is made, confirming the theoretical results.


title: Heat Transfer in Lattice BGK Modeled Fluid

abstract: The thermal lattice BGK model is a recently suggested numerical tool which aims to simulate thermohydrodynamic problems. We investigate the quality of the lattice BGK simulation by calculating the temperature profiles in the Couette flow under different Eckert and Mach numbers. A revised lower order model is proposed and the higher order model is proved once again to be advantageous.


title: Deterministic Chaos in Tropical Atmospheric Dynamics

abstract: We examine an 11-year data set from the tropical weather station of Tlaxcala, Mexico. We find that mutual information drops quickly with the delay, to a positive value which relaxes to zero with a time scale of 20 days. We also examine the mutual dependence of the observables and conclude that the data set gives the equivalent of 8 variables per day, known to a precision of $2\%$. We determine the effective dimension of the attractor to be $D_{eff} \approx 11.7$ at the scale $3.5\% < R/R_{max} < 8\%$. We find evidence that the effective dimension increases as $R/R_{max} \to 0$, supporting a conjecture by Lorenz that the climate system may consist of a large number of weakly coupled subsystems, some of which have low-dimensional attractors. We perform a local reconstruction of the dynamics in phase space; the short-term predictability is modest and agrees with theoretical estimates. Useful skill in predictions of 10-day rainfall accumulation anomalies reflects the persistence of weather patterns, which follow the 20-day decay rate of the mutual information.


title: Renormalization of Lattice Gas Transport Coefficients

abstract: A method is described for calculating corrections to the Boltzmann/Chapman-Enskog analysis of lattice gases due to the buildup of correlations. It is shown that renormalized transport coefficients can be calculated perturbatively by summing terms in an infinite series. A diagrammatic notation for the terms in this series is given, in analogy with the Feynman diagrams of quantum field theory. This theory is applied to an example lattice gas and shown to correctly predict experimental deviation from the Boltzmann prediction.


title: Algebraic Spatial Correlations and Non-Gibbsian Equilibrium States

abstract: Non-Gibbsian stationary states occur in dissipative non-equilibrium systems. They are closely connected with the lack of detailed balance and the absence of a fluctuation-dissipation theorem. These states exhibit spatial correlations that are long ranged under generic conditions, even in systems with short range interactions, provided the system has slow modes and some degree of spatial anisotropy. In this paper we present a theory for static pair correlations in lattice gas automata violating detailed balance, and we show that the spatially uniform non-Gibbsian equilibrium state exhibits long range correlations, even in the absence of an external driving field.


title: Renormalized Equilibria of a Schloegl Model Lattice Gas

abstract: A lattice gas model for Schloegl's second chemical reaction is described and analyzed. Because the lattice gas does not obey a semi-detailed-balance condition, the equilibria are non-Gibbsian. In spite of this, a self-consistent set of equations for the exact homogeneous equilibria are described, using a generalized cluster-expansion scheme. These equations are solved in the two-particle BBGKY approximation, and the results are compared to numerical experiment. It is found that this approximation describes the equilibria far more accurately than the Boltzmann approximation. It is also found, however, that spurious solutions to the equilibrium equations appear which can only be removed by including effects due to three-particle correlations.


title: New Results for Diffusion in Lorentz Lattice Gas Cellular Automata

abstract: New calculations to over ten million time steps have revealed a more complex diffusive behavior than previously reported, of a point particle on a square and triangular lattice randomly occupied by mirror or rotator scatterers. For the square lattice fully occupied by mirrors where extended closed particle orbits occur, anomalous diffusion was still found. However, for a not fully occupied lattice the super diffusion, first noticed by Owczarek and Prellberg for a particular concentration, obtains for all concentrations. For the square lattice occupied by rotators and the triangular lattice occupied by mirrors or rotators, an absence of diffusion (trapping) was found for all concentrations, except on critical lines, where anomalous diffusion (extended closed orbits) occurs and hyperscaling holds for all closed orbits with {\em universal} exponents ${\displaystyle{d_f = \frac{7}{4}}}$ and ${\displaystyle{\tau = \frac{15}{7}}}$. Only one point on these critical lines can be related to a corresponding percolation problem. The questions arise therefore whether the other critical points can be mapped onto a new percolation-like problem, and of the dynamical significance of hyperscaling.


title: Analytical solutions of the lattice Boltzmann BGK model

abstract: Analytical solutions of the two dimensional triangular and square lattice Boltzmann BGK models have been obtained for the plain Poiseuille flow and the plain Couette flow. The analytical solutions are written in terms of the characteristic velocity of the flow, the single relaxation time $\tau $ and the lattice spacing. The analytic solutions are the exact representation of these two flows without any approximation.


title: Prandtl number of lattice Bhatnagar-Gross-Krook fluid

abstract: The lattice Bhatnagar-Gross-Krook modeled fluid has an unchangeable unit Prandtl number. A simple method is introduced in this letter to formulate a flexible Prandtl number for the modeled fluid. The effectiveness was demonstrated by numerical simulations of the Couette flow.


title: Diffusion in Lorentz Lattice Gas Cellular Automata: the honeycomb and quasi-lattices compared with the square and triangular lattices

abstract: We study numerically the nature of the diffusion process on a honeycomb and a quasi-lattice, where a point particle, moving along the bonds of the lattice, scatters from randomly placed scatterers on the lattice sites according to strictly deterministic rules. For the honeycomb lattice fully occupied by fixed rotators two (symmetric) isolated critical points appear to be present, with the same hyperscaling relation as for the square and the triangular lattices. No such points appear to exist for the quasi-lattice. A comprehensive comparison is made with the behavior on the previously studied square and triangular lattices. A great variety of diffusive behavior is found ranging from propagation, super-diffusion, normal, quasi-normal, anomalous to absence of diffusion. The influence of the scattering rules as well as of the lattice structure on the diffusive behavior of a point particle moving on the all lattices studied so far, is summarized.


title: A Method for Near-Equilibrium Discrete-Velocity Gas Flows

abstract: We present a simulation scheme for discrete-velocity gases based on {\em local thermodynamic equilibrium}. Exploiting the kinetic nature of discrete-velocity gases, in that context, results in a natural splitting of fluxes, and the resultant scheme strongly resembles the original processes. The kinetic nature of the scheme and the modeling of the {\em infinite collision rate} limit, result in a small value of the coefficient of (numerical)-viscosity, the behavior of which is remarkably physical [18]. A first order method, and two second order methods using the total variation diminishing principle are developed and an example application presented. Given the same computer resources, it is expected that with this approach, much higher Reynold's number will be achievable than presently possible with either lattice gas automata or lattice Boltzmann approaches. The ideas being general, the scheme is applicable to any discrete-velocity model, and to lattice gases as well.


title: Shock structure in a nine-velocity gas

abstract: The exact structure of a shock is computed in a multiple-speed discrete-velocity gas, the nine-velocity gas, wherein the multiplicity of speeds ensures nontrivial thermodynamics. Obtained as a solution of the model Boltzmann equations, the procedure consists of tracking the shock as a trajectory of a three dimensional dynamical system connecting an equilibrium upstream state to an equilibrium downstream state. The two equilibria satisfy the jump conditions obtained from the model Euler equations. Comparison of the shock structure to that in a monatomic perfect gas, as given by the Navier-Stokes equation, shows excellent agreement. The shock in the nine-velocity gas has an overshoot in entropy alone, like in a monatomic gas. The near-equilibrium flow technique for discrete-velocity gases (Nadiga \& Pullin [2]), a kinetic flux-splitting method based on the local thermodynamic equilibrium, is also seen to capture the shock structure remarkably well.


title: Plane Waves in a Multispeed Discrete-Velocity Gas

abstract: A kinetic flux-splitting procedure used in conjunction with local thermodynamic equilibrium in a finite volume allows us to investigate numerically discrete-velocity gas flows. The procedure, outlined for a general discrete-velocity gas, is used to simulate flows of the nine-velocity gas, a simple two dimensional multiple-speed discrete-velocity gas, wherein a multiplicity of speeds ensures nontrivial thermodynamics. After verifying the linear wave limit and the non-linear steepening of wavefronts, the stability and propagation of planar discontinuities in that model gas is studied. The supersonic-subsonic requirement for the stable propagation of a discontinuity, being kinematic in nature is the same in the model gas, as e.g., in a perfect gas. However, the finiteness of the velocity space in the model gas does not allow a translation of the above kinematic condition to the thermodynamic requirement of increasing entropy across a compressive shock: a case of an entropy decreasing compressive shock in the model gas is presented. Finally, the interaction of various types of waves--shock waves, rarefactions and contact surfaces--in the model gas are shown in a simulation of the shock tube problem.


title: An Euler Solver Based on Locally Adaptive Discrete Velocities

abstract: A new discrete-velocity model is presented to solve the three-dimensional Euler equations. The velocities in the model are of an adaptive nature---both the origin of the discrete-velocity space and the magnitudes of the discrete-velocities are dependent on the local flow--- and are used in a finite volume context. The numerical implementation of the model follows the near-equilibrium flow method of Nadiga and Pullin [1] and results in a scheme which is second order in space (in the smooth regions and between first and second order at discontinuities) and second order in time. (The three-dimensional code is included.) For one choice of the scaling between the magnitude of the discrete-velocities and the local internal energy of the flow, the method reduces to a flux-splitting scheme based on characteristics. As a preliminary exercise, the result of the Sod shock-tube simulation is compared to the exact solution.


title: Nonlinear stability analysis of plane Poiseuille flow by normal forms

abstract: In the subcritical interval of the Reynolds number 4320\leq R\leq R_c\equiv 5772, the Navier--Stokes equations of the two--dimensional plane Poiseuille flow are approximated by a 22--dimensional Galerkin representation formed from eigenfunctions of the Orr--Sommerfeld equation. The resulting dynamical system is brought into a generalized normal form which is characterized by a disposable parameter controlling the magnitude of denominators of the normal form transformation. As rigorously proved, the generalized normal form decouples into a low--dimensional dominant and a slaved subsystem. {}From the dominant system the critical amplitude is calculated as a function of the Reynolds number. As compared with the Landau method, which works down to R=5300, the phase velocity of the critical mode agrees within 1 per cent; the critical amplitude is reproduced similarly well except close to the critical point, where the maximal error is about 16 per cent. We also examine boundary conditions which partly differ from the usual ones.


title: Lattice Boltzmann Simulation of Non-Ideal Fluids

abstract: A lattice Boltzmann scheme able to model the hydrodynamics of phase separation and two-phase flow is described. Thermodynamic consistency is ensured by introducing a non-ideal pressure tensor directly into the collision operator. We also show how an external chemical potential can be used to supplement standard boundary conditions in order to investigate the effect of wetting on phase separation and fluid flow in confined geometries. The approach has the additional advantage of reducing many of the unphysical discretisation problems common to previous lattice Boltzmann methods.


title: Lattice Boltzmann Approach to Viscous Flows Between Parallel Plates

abstract: Four different kinds of laminar flows between two parallel plates are investigated using the Lattice Boltzmann Method (LBM). The LBM accuracy is estimated in two cases using numerical fits of the parabolic velocity profiles and the kinetic energy decay curves, respectively. The error relative to the analytical kinematic viscosity values was found to be less than one percent in both cases. The LBM results for the unsteady development of the flow when one plate is brought suddenly at a constant velocity, are found in excellent agreement with the analytical solution. Because the classical Schlichting's approximate solution for the entrance--region flow is not valid for small Reynolds numbers, a Finite Element Method solution was used in order to check the accuracy of the LBM results.


title: Cellular automaton model of precipitation/dissolution coupled with solute transport

abstract: Precipitation/dissolution reactions coupled with solute transport are modelled as a cellular automaton in which solute molecules perform a random walk on a regular lattice and react according to a local probabilistic rule. Stationary solid particles dissolve with a certain probability and, provided solid is already present or the solution is saturated, solute particles have a probability to precipitate. In our simulation of the dissolution of a solid block inside uniformly flowing water we obtain solid precipitation downstream from the original solid edge, in contrast to the standard reaction-transport equations. The observed effect is the result of fluctuations in solute density and diminishes when we average over a larger ensemble. The additional precipitation of solid is accompanied by a substantial reduction in the relatively small solute concentration. The model is appropriate for the study of the r\^ole of intrinsic fluctuations in the presence of reaction thresholds and can be employed to investigate porosity changes associated with the carbonation of cement.


title: Multi-component lattice-Boltzmann model with interparticle interaction

abstract: A previously proposed [X. Shan and H. Chen, Phys. Rev. E {\bf 47}, 1815, (1993)] lattice Boltzmann model for simulating fluids with multiple components and interparticle forces is described in detail. Macroscopic equations governing the motion of each component are derived by using Chapman-Enskog method. The mutual diffusivity in a binary mixture is calculated analytically and confirmed by numerical simulation. The diffusivity is generally a function of the concentrations of the two components but independent of the fluid velocity so that the diffusion is Galilean invariant. The analytically calculated shear kinematic viscosity of this model is also confirmed numerically.


title: Lattice gas with ``interaction potential''

abstract: We present an extension of a simple automaton model to incorporate non-local interactions extending over a spatial range in lattice gases. {}From the viewpoint of Statistical Mechanics, the lattice gas with interaction range may serve as a prototype for non-ideal gas behavior. {}From the density fluctuations correlation function, we obtain a quantity which is identified as a potential of mean force. Equilibrium and transport properties are computed theoretically and by numerical simulations to establish the validity of the model at macroscopic scale.


title: Lattice Boltzmann Model For Magnetic Fluids

abstract: A lattice Boltzmann model with interacting particles was developed in order to simulate the magneto-rheological characteristics of magnetic fluids. In the frame of this model, $6\, +\,1$ species of particles are allowed to move across a $2D$ triangular lattice. Among these species, $6$ of them carry an individual magnetic dipole moment and interact themselves not only as a result oflocal collisions, as in current Lattice Boltzmann models, but also as a result of nearest neighbours magnetic dipole-dipole interaction. The relative distribution of the individual magnetic moments is determined by the intensity of an external static magnetic field acting on the whole system. This model exhibits some relevant characteristics of real magnetic fluids, i.e., anisotropic structure formation as a result of magnetic field induced gas-liquid phase transition and magnetic field dependence of the sound velocity and the attenuation coefficient.


title: On reconstruction of chaotic attractor from time series represented as "clusters"

abstract: The method of reconstruction of an attractor from a set of short time series ({\it clusters}) is proposed and discussed. This method is most useful for correlation dimension estimation of experimental data.


title: Predicting and characterizing data sequences from structure-variable systems

abstract: In principle, all the natural systems such as biological, ecological and economical systems are structure-variable systems (in which some environment parameters are not fixed). In this Letter we show that data sequences from many structure-variable systems are short-term predictable. We also argue how to characterize the data sequences from structure-variable systems.


title: SAMPLING ALMOST PERIODIC FUNCTIONS WITH RANDOM PROBES OF FINITE DENSITY

abstract: We consider the problem of reconstructing a function given its values on a set of points with finite density. We prove that with probability one, the values of an almost periodic function on a random array of points (with finite density) completely determine the function. We also give some properties of the associated Blaschke product.


title: A Lattice-Gas Model of Microemulsions

abstract: We develop a lattice gas model for the nonequilibrium dynamics of microemulsions. Our model is based on the immiscible lattice gas of Rothman and Keller, which we reformulate using a microscopic, particulate description so as to permit generalisation to more complicated interactions, and on the prescription of Chan and Liang for introducing such interparticle interactions into lattice gas dynamics. We present the results of simulations to demonstrate that our model exhibits the correct phenomenology, and we contrast it with both equilibrium lattice models of microemulsions, and to other lattice gas models.

categories: comp-gas chem-ph cond-mat mtrl-th nlin.CG


title: Analysis and boundary condition of the lattice Boltzmann BGK model with two velocity components

abstract: In this paper, we study the two dimensional lattice Boltzmann BGK model (LBGK) by analytically solving a simple flow in a 2~-D channel. The flow is driven by the movement of upper boundary with vertical injection fluid at the porous boundaries. The velocity profile is shown to satisfy a second-order finite-difference form of the simplified incompressible Navier-Stokes equation. With the analysis, different boundary conditions can be studied theoretically. A momentum exchange principle is also revealed at the boundaries. A general boundary condition for any given velocity boundary is proposed based on the analysis.


title: Estimating Predictability: Redundancy and Surrogate Data Method

abstract: A method for estimating theoretical predictability of time series is presented, based on information-theoretic functionals---redundancies and surrogate data technique. The redundancy, designed for a chosen model and a prediction horizon, evaluates amount of information between a model input (e.g., lagged versions of the series) and a model output (i.e., a series lagged by the prediction horizon from the model input) in number of bits. This value, however, is influenced by a method and precision of redundancy estimation and therefore it is a) normalized by maximum possible redundancy (given by the precision used), and b) compared to the redundancies obtained from two types of the surrogate data in order to obtain reliable classification of a series as either unpredictable or predictable. The type of predictability (linear or nonlinear) and its level can be further evaluated. The method is demonstrated using a numerically generated time series as well as high-frequency foreign exchange data and the theoretical predictability is compared to performance of a nonlinear predictor.


title: Detecting nonlinearity in multivariate time series

abstract: We propose an extension to time series with several simultaneously measured variables of the nonlinearity test, which combines the redundancy -- linear redundancy approach with the surrogate data technique. For several variables various types of the redundancies can be defined, in order to test specific dependence structures between/among (groups of) variables. The null hypothesis of a multivariate linear stochastic process is tested using the multivariate surrogate data. The linear redundancies are used in order to avoid spurious results due to imperfect surrogates. The method is demonstrated using two types of numerically generated multivariate series (linear and nonlinear) and experimental multivariate data from meteorology and physiology.


title: On pressure and velocity flow boundary conditions for the lattice Boltzmann BGK model

abstract: Pressure (density) and velocity boundary conditions inside a flow domain are studied for 2-D and 3-D lattice Boltzmann BGK models (LBGK) and new method to specify these conditions are proposed. These conditions are consistent with the boundary condition we proposed in a previous paper using an idea of bounce-back of non-equilibrium distribution. These conditions give excellent results for the regular LBGK models, and were shown to be second-order accurate by numerical examples. When they are used together with the improved incompressible LBGK model proposed by zou et al. the simulation results recover the analytical solution of the plane Poiseuille flow driven by pressure (density) difference with machine accuracy.


title: A non-slip boundary condition for lattice Boltzmann simulations

abstract: A non-slip boundary condition at a wall for the lattice Boltzmann method is presented. In the present method unknown distribution functions at the wall are assumed to be an equilibrium distribution function with a counter slip velocity which is determined so that fluid velocity at the wall is equal to the wall velocity. Poiseuille flow and Couette flow are calculated with the nine-velocity model to demonstrate the accuracy of the present boundary condition.


title: CAM-8: A Computer Architecture Based on Cellular Automata

abstract: Conventional computers are ill suited to run CA models, and so discourage their development. By creating a hardware platform that makes a broad range of new CA algorithms practical for real applications, we hope to whet the appetite of researchers for the astronomical computing power that can be harnessed in microphysics in a CA format.


title: Detecting Nonlinearity in Pressure Data Inside the Draft Tube of a Real Francis Turbine

abstract: A general method for testing nonlinearity in time series is described and applied to measurements of different pressure data inside the draft tube surge of a real Francis turbine. Comparing the current original time series to an ensemble of surrogates time series, suitably constructed to mimic the linear properties of the original one, we was able to distinguish a linear stochastic from a nonlinear deterministic behaviour and, moreover, to quantify the degree of nonlinearity present in the related dynamics. The problem of detecting nonlinear structure in real data is quite complicated by the influence of various contaminations, like broadband noise and/or long coherence times. These difficulties have been overcame using the combination of a suitable nonlinear filtering technique and a qualitative redundancy statistic analysis. The above investigations allow a quantitative characterization of different dynamical regimes of motion of gas cavities inside real turbines and, moreover, allow to support the reliability of some related mathematical modelizations.


title: On Surrogate Data Testing for Linearity based on the Periodogram

abstract: The method of surrogate data is a tool to test whether data were generated by some class of model. Tests based on the periodogram have been proposed to decide if linear systems driven by Gaussian noise could have generated a sample time series. We show that this procedure based on the periodogram, in general, misspecifies the test statistic. Based on the theory of linear systems we suggest an alternative procedure to obtain the correct distribution of the test statistic and discuss problems of this approach.


title: The developing structure of dynamical systems

abstract: In order to investigate the evolutionary process of many deterministic Dynamical systems with unfixed parameter, a set of dynamical models with parameter changing continuously and the accumulation of this change might be large is introduced and discussed. The boundary crises and the period-doubling bifurcations are found suppressed and scaling properties for these phenomena are exploited. Due to this suppression, the period-doubling bifurcations seem no longer continuous. We further discuss the possible applications of these models.


title: Parallel Algorithms on the ASTRA SIMD Machine

abstract: In view of the tremendous computing power jump of modern RISC processors the interest in parallel computing seems to be thinning out. Why use a complicated system of parallel processors, if the problem can be solved by a single powerful micro-chip? It is a general law, however, that exponential growth will always end by some kind of a saturation, and then parallelism will again become a hot topic. We try to prepare ourselves for this eventuality. The MPPC project started in 1990 in the heydeys of parallelism and produced four ASTRA machines (presented at CHEP'92) with 4k processors (which are expandable to 16k) based on `yesterday's chip-technology' (chip presented at CHEP'91). These machines now provide excellent test-beds for algorithmic developments in a complete, real environment. We are developing for example fast-pattern recognition algorithms which could be used in high-energy physics experiments at the LHC (planned to be operational after 2004 at CERN) for triggering and data reduction. The basic feature of our ASP (Associative String Processor) approach is to use extremely simple (thus very cheap) processor elements but in huge quantities (up to millions of processors) connected together by a very simple string-like communication chain. In this paper we present powerful algorithms based on this architecture indicating the performance perspectives if the hardware (i.e. chip fabrication) quality reaches `present or even future technology levels'.


title: A Lattice Boltzmann Model of Binary Fluid Mixture

abstract: We introduce a lattice Boltzmann for simulating an immiscible binary fluid mixture. Our collision rules are derived from a macroscopic thermodynamic description of the fluid in a way motivated by the Cahn-Hilliard approach to non-equilibrium dynamics. This ensures that a thermodynamically consistent state is reached in equilibrium. The non-equilibrium dynamics is investigated numerically and found to agree with simple analytic predictions in both the one-phase and the two-phase region of the phase diagram.


title: Mutual Information as a Bayesian Measure of Independence

abstract: The problem of hypothesis testing is examined from both the historical and Bayesian points of view in the case that sampling is from an underlying joint probability distribution and the hypotheses tested for are those of independence and dependence of the underlying distribution. Exact results for the Bayesian method are provided. Asymptotic Bayesian results and historical method test quantities are compared, and historical method quantities are interpreted in terms of clearly defined Bayesian quantities. The asymptotic Bayesian test relies upon a statistic that is primarily mutual information.


title: Investigations of a Two-Phase Fluid Model

abstract: We study an interface-capturing two-phase fluid model in which the interfacial tension is modelled as a volumetric stress. Since these stresses are obtainable from a Van der Waals-Cahn-Hilliard free energy, the model is, to a certain degree, thermodynamically realistic. Thermal fluctuations are not considered presently for reasons of simplicity. The utility of the model lies in its momentum-conservative representation of surface tension and the simplicity of its numerical implementation resulting from the volumetric modelling of the interfacial dynamics. After validation of the model in two spatial dimensions, two prototypical applications---instability of an initially high-Reynolds-number liquid jet in the gaseous phase and spinodal decomposition in a liquid-gas system--- are presented.


title: Lattice Gas Automata for Reactive Systems

abstract: Reactive lattice gas automata provide a microscopic approachto the dynamics of spatially-distributed reacting systems. After introducing the subject within the wider framework of lattice gas automata (LGA) as a microscopic approach to the phenomenology of macroscopic systems, we describe the reactive LGA in terms of a simple physical picture to show how an automaton can be constructed to capture the essentials of a reactive molecular dynamics scheme. The statistical mechanical theory of the automaton is then developed for diffusive transport and for reactive processes, and a general algorithm is presented for reactive LGA. The method is illustrated by considering applications to bistable and excitable media, oscillatory behavior in reactive systems, chemical chaos and pattern formation triggered by Turing bifurcations. The reactive lattice gas scheme is contrasted with related cellular automaton methods and the paper concludes with a discussion of future perspectives.


title: Coarse-grained entropy rates for characterization of complex time series

abstract: A method for classification of complex time series using coarse-grained entropy rates (CER's) is presented. The CER's, which are computed from information-theoretic functionals -- redundancies, are relative measures of regularity and predictability, and for data generated by dynamical systems they are related to Kolmogorov-Sinai entropy. A deterministic dynamical origin of the data under study, however, is not a necessary condition for the use of the CER's, since the entropy rates can be defined for stochastic processes as well. Sensitivity of the CER's to changes in data dynamics and their robustness with respect to noise are tested by using numerically generated time series resulted from both deterministic -- chaotic and stochastic processes. Potential application of the CER's in analysis of physiological signals or other complex time series is demonstrated by using examples from pharmaco-EEG and tremor classification.


title: Generating high-dimensional chaotic signals by the sum

abstract: The method of generating of high-dimensional oscillations on the basis of summing of low-dimensional chaotic signals of noncoupled dynamical systems is investigated. It is shown that the correlation dimension of attractor reconstructed from such oscillations is equal to the sum of dimensions of the original signals. The possibilities for obtaining of homogeneous high-dimensional chaotic attractor are numerically demonstrated with the original signals of R{\"{o}}ssler system.


title: Integer Lattice Gases

abstract: We generalize the hydrodynamic lattice gas model to include arbitrary numbers of particles moving in each lattice direction. For this generalization we derive the equilibrium distribution function and the hydrodynamic equations, including the equation of state and the prefactor of the inertial term that arises from the breaking of galilean invariance in these models. We show that this prefactor can be set to unity in the generalized model, therby effectively restoring galilean invariance. Moreover, we derive an expression for the kinematic viscosity, and show that it tends to decrease with the maximum number of particles allowed in each direction, so that higher Reynolds numbers may be achieved. Finally, we derive expressions for the statistical noise and the Boltzmann entropy of these models.


title: State Space Reconstruction Parameters in the Analysis of Chaotic Time Series - the Role of the Time Window Length

abstract: The most common state space reconstruction method in the analysis of chaotic time series is the Method of Delays (MOD). Many techniques have been suggested to estimate the parameters of MOD, i.e. the time delay $\tau$ and the embedding dimension $m$. We discuss the applicability of these techniques with a critical view as to their validity, and point out the necessity of determining the overall time window length, $\tau_w$, for successful embedding. Emphasis is put on the relation between $\tau_w$ and the dynamics of the underlying chaotic system, and we suggest to set $\tau_w \geq \tau_p$, the mean orbital period; $\tau_p$ is approximated from the oscillations of the time series. The procedure is assessed using the correlation dimension for both synthetic and real data. For clean synthetic data, values of $\tau_w$ larger than $\tau_p$ always give good results given enough data and thus $\tau_p$ can be considered as a lower limit ($\tau_w \geq \tau_p$). For noisy synthetic data and real data, an upper limit is reached for $\tau_w$ which approaches $\tau_p$ for increasing noise amplitude.


title: Constrained-Realization Monte-Carlo Method for Hypothesis Testing

abstract: We compare two theoretically distinct approaches to generating artificial (or ``surrogate'') data for testing hypotheses about a given data set. The first and more straightforward approach is to fit a single ``best'' model to the original data, and then to generate surrogate data sets that are ``typical realizations'' of that model. The second approach concentrates not on the model but directly on the original data; it attempts to constrain the surrogate data sets so that they exactly agree with the original data for a specified set of sample statistics. Examples of these two approaches are provided for two simple cases: a test for deviations from a gaussian distribution, and a test for serial dependence in a time series. Additionally, we consider tests for nonlinearity in time series based on a Fourier transform (FT) method and on more conventional autoregressive moving-average (ARMA) fits to the data. The comparative performance of hypothesis testing schemes based on these two approaches is found to depend on whether or not the discriminating statistic is pivotal. A statistic is ``pivotal'' if its distribution is the same for all processes consistent with the null hypothesis. The typical-realization method requires that the discriminating statistic satisfy this property. The constrained-realization approach, on the other hand, does not share this requirement, and can provide an accurate and powerful test without having to sacrifice flexibility in the choice of discriminating statistic.


title: Lattice-gas simulations of Domain Growth, Saturation and Self-Assembly in Immiscible Fluids and Microemulsions

abstract: We investigate the dynamical behavior of both binary fluid and ternary microemulsion systems in two dimensions using a recently introduced hydrodynamic lattice-gas model of microemulsions. We find that the presence of amphiphile in our simulations reduces the usual oil-water interfacial tension in accord with experiment and consequently affects the non-equilibrium growth of oil and water domains. As the density of surfactant is increased we observe a crossover from the usual two-dimensional binary fluid scaling laws to a growth that is {\it slow}, and we find that this slow growth can be characterized by a logarithmic time scale. With sufficient surfactant in the system we observe that the domains cease to grow beyond a certain point and we find that this final characteristic domain size is inversely proportional to the interfacial surfactant concentration in the system.


title: Estimating the Fractal Dimension, K_2-entropy, and the Predictability of the Atmosphere

abstract: The series of mean daily temperature of air recorded over a period of 215 years is used for analysing the dimensionality and the predictability of the atmospheric system. The total number of data points of the series is 78527. Other 37 versions of the original series are generated, including ``seasonally adjusted'' data, a smoothed series, series without annual course, etc. Modified methods of Grassberger and Procaccia are applied. A procedure for selection of the ``meaningful'' scaling region is proposed. Several scaling regions are revealed in the ln C(r) versus ln r diagram. The first one in the range of larger ln r has a gradual slope and the second one in the range of intermediate ln r has a fast slope. Other two regions are settled in the range of small ln r. The results lead us to claim that the series arises from the activity of at least two subsystems. The first subsystem is low-dimensional (d_f=1.6) and it possesses the potential predictability of several weeks. We suggest that this subsystem is connected with seasonal variability of weather. The second subsystem is high-dimensional (d_f>17) and its error-doubling time is about 4-7 days. It is found that the predictability differs in dependence on season. The predictability time for summer, winter and the entire year (T_2 approx. 4.7 days) is longer than for transition-seasons (T_2 approx. 4.0 days for spring, T_2 approx. 3.6 days for autumn). The role of random noise and the number of data points are discussed. It is shown that a 15-year-long daily temperature series is not sufficient for reliable estimations based on Grassberger and Procaccia algorithms.

categories: comp-gas ao-sci chao-dyn nlin.CD nlin.CG


title: Two-Dimensional Navier-Stokes Simulation of Breaking Waves

abstract: Numerical simulations describing plunging breakers including the splash-up phenomenon are presented. The motion is governed by the classical, incompressible, two-dimensional Navier-Stokes equation. The numerical modelling of this two-phase flow is based on a piecewise linear version of the volume of fluid method. Preliminary results concerning the time evolution of the liquid-gas interface and the vorticity field are given for short waves, showing how an initial steep wave undergoes breaking and successive splash-up cycles.


title: Diffusion in a multi-component Lattice Boltzmann Equation model

abstract: Diffusion phenomena in a multiple component lattice Boltzmann Equation (LBE) model are discussed in detail. The mass fluxes associated with different mechanical driving forces are obtained using a Chapman-Enskog analysis. This model is found to have correct diffusion behavior and the multiple diffusion coefficients are obtained analytically. The analytical results are further confirmed by numerical simulations in a few solvable limiting cases. The LBE model is established as a useful computational tool for the simulation of mass transfer in fluid systems with external forces.