scispace - formally typeset
Search or ask a question

Showing papers in "arXiv: Cellular Automata and Lattice Gases in 1996"


Journal ArticleDOI
TL;DR: In this paper, Chen et al. used the half-way wall bounceback boundary condition for the 2-D Poiseuille flow with forcing to obtain second-order accuracy for the 3-D square duct flow.
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 a new method to specify these conditions are proposed. These conditions are constructed in consistency of the wall boundary condition based on an idea of bounceback of non-equilibrium distribution. When these conditions are used together with the improved incompressible LBGK model by Zou et al., the simulation results recover the analytical solution of the plane Poiseuille flow driven by pressure (density) difference with machine accuracy. Since the half-way wall bounceback boundary condition is very easy to implement and was shown theoretically to give second-order accuracy for the 2-D Poiseuille flow with forcing, it is used with pressure (density) inlet/outlet conditions proposed in this paper and in Chen et al. to study the 2-D Poiseuille flow and the 3-D square duct flow. The numerical results are approximately second-order accurate. The magnitude of the error of the half-way wall bounceback is comparable with that using some other published boundary conditions. Besides, the bounceback condition has a much better stability behavior than that of other boundary conditions.

2,001 citations


Journal ArticleDOI
TL;DR: The most common state space reconstruction method in the analysis of chaotic time series is the Method of Delays (MOD), and many techniques have been suggested to estimate the parameters of MOD, i.e., the time delay and the embedding dimension $m$ as mentioned in this paper.
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.

264 citations


Journal ArticleDOI
TL;DR: The authors compare two theoretically distinct approaches to generating artificial data sets for testing hypotheses about a given data set, i.e., fitting a single model to the original data, and then generating surrogate data sets that are ''typical realizations'' of that model.
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.

17 citations


Journal ArticleDOI
TL;DR: In this paper, the authors used the series of mean daily temperature of air recorded over a period of 215 years for analysing the dimensionality and the predictability of the atmospheric system.
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.

15 citations


Posted Content
TL;DR: The growing need for a better understanding of nonlinear processes in plasma physics has in the last decades stimulated the development of new and more advanced data analysis techniques as mentioned in this paper, and some of the basic properties one may wish to infer from a data set and then presents appropriate analysis techniques with some recent applications.
Abstract: The growing need for a better understanding of nonlinear processes in plasma physics has in the last decades stimulated the development of new and more advanced data analysis techniques. This review lists some of the basic properties one may wish to infer from a data set and then presents appropriate analysis techniques with some recent applications. The emphasis is put on the investigation of nonlinear wave phenomena and turbulence in space plasmas.

9 citations


Journal ArticleDOI
TL;DR: In this article, it is shown that using massive parallelism, particle trajectories can be simulated in a time that scales logarithmically in the length of the trajectory.
Abstract: The Lorentz lattice gas is studied from the perspective of computational complexity theory. It is shown that using massive parallelism, particle trajectories can be simulated in a time that scales logarithmically in the length of the trajectory. This result characterizes the ``logical depth" of the Lorentz lattice gas and allows us to compare it to other models in statistical physics.

1 citations


Journal ArticleDOI
TL;DR: The lattice Boltzmann equation (LBE) is a microscopically-inspired method designed to solve macroscopic fluid dynamics problems as discussed by the authors, and it lives at the interface between the microscopic (molecular) and macroscopy (continuum) worlds, hopefully capturing the best of the two.
Abstract: The lattice Boltzmann equation (LBE) is a microscopically-inspired method designed to solve macroscopic fluid dynamics problems. As a such, it lives at the interface between the microscopic (molecular) and macroscopic (continuum) worlds, hopefully capturing the best of the two. In this paper we shall discuss whether or not, after almost a decade since its inception, LBE has lived up to the initial expectations. Open problems and future research developments are also briefly outlined.

Posted Content
TL;DR: In this paper, the authors present numerical simulations of plunging breakers including the splash-up phenomenon, which is governed by the classical Navier-Stokes equation and is based on a piecewise linear version of the volume of fluid method.
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.