Dynamic properties of Lennard Jones fluids and liquid metals

The dependence of the dynamic properties of liquid metals and Lennard-Jones fluids on the characteristics of the interaction potentials is analyzed. Molecular-dynamics simulations of liquids in analogous conditions but assuming that their particles interact either through a Lennard-Jones or a liquid-metal potential were carried out. The Lennard-Jones potentials were chosen so that both the effective size of the particles and the depth of the potential well were very close to those of the liquid-metal potentials. In order to investigate the extent to which the dynamic properties of liquids depend on the short-range attractive interactions as well as on the softness of the potential cores, molecular-dynamics simulations of the same systems but assuming purely repulsive interactions with the same potential cores were also performed. The study includes both singleparticle dynamic properties, such as the velocity autocorrelation functions, and collective dynamic properties, such as the intermediate scattering funcfunctions, and collective dynamic properties, such as the intermediate scattering functions, the dynamic structure factors, the longitudinal and transverse current correlations, and the transport coefficients.

Manifolds on the verge of a hyperbolicity breakdown

We study numerically the disappearance of normally hyperbolic invariant tori in quasiperiodic systems and identify a scenario for their breakdown. In this scenario, the breakdown happens because two invariant directions of the transversal dynamics come close to each other, losing their regularity. On the other hand, the Lyapunov multipliers associated with the invariant directions remain more or less constant. We identify notable quantitative regularities in this scenario, namely that the minimum angle between the two invariant directions and the Lyapunov multipliers have power law dependence with the parameters. The exponents of the power laws seem to be universal.

Stochastic generation of homogeneous isotropic turbulence with well-defined-spectra

A precise and simple computational model to generate well-behaved two-dimensional turbulent flows is presented. The whole approach rests on the use of stochastic differential equations and is general enough to reproduce a variety of energy spectra and spatiotemporal correlation functions. Analytical expressions for both the continuous and the discrete versions, together with simulation algorithms, are derived. Results for two relevant spectra, covering distinct ranges of wave numbers, are given.

Front propagation in spatially modulated media

A two-dimensional reaction-diffusion front which propagates in a modulated medium is studied. The modulation consists of a spatial variation of the local front velocity in the transverse direction to that of the front propagation. We study analytically and numerically the final steady-state velocity and shape of the front, resulting from a nontrivial interplay between the local curvature effects and the global competition process between different maxima of the control parameter. The transient dynamics of the process is also studied numerically and analytically by means of singular perturbation techniques.

Nonequilibrium orientational patterns in two-component Langmuir monolayers

A model of a phase-separating two-component Langmuir monolayer in the presence of a photoinduced reaction interconverting two components is formulated. An interplay between phase separation, orientational ordering, and reaction is found to lead to a variety of nonequilibrium self-organized patterns, both stationary and traveling. Examples of the patterns, observed in numerical simulations, include flowing droplets, traveling stripes, wave sources, and vortex defects.

Self-organizing propagation patterns from dynamic self-assembly in monolayers

Propagation of localized orientational waves, as imaged by Brewster angle microscopy, is induced by low intensity linearly polarized light inside axisymmetric smectic-C confined domains in a photosensitive molecular thin film at the air/water interface (Langmuir monolayer). Results from numerical simulations of a model that couples photoreorientational effects and long-range elastic forces are presented. Differences are stressed between our scenario and the paradigmatic wave phenomena in excitable chemical media.

Velocity fluctuations of fracturelike disruptions of associating polymer solutions

Velocity has been measured as a function of time for propagating crack tips as water is injected into solutions of end-capped associating polymers in a rectanguar Hele-Shaw cell. Measurements were performed for flows with different values of cell gap, channel width, polymer molecular weight, and polymer concentration. The condition for the onset of fracturelike behavior is well described by a Deborah number which uses the shear-thinning shear rate of the polymer solution as a characteristic frequency for network relaxation. At low molecular weight, the onset of fracturelike pattern evolution is accompanied by an abrupt jump in tip velocity, followed by a lower and approximately constant acceleration. At high molecular weight, the transition to fracturelike behavior involves passing through a regime that may be understood in terms of stick-slip dynamics. The crack-tip wanders from side to side and fluctuates (in both speed and velocity along the channel) with a characteristic frequency which depends linearly on the invading fluid injection rate.

Chemoconvection patterns in the Methylene-blue-glucose system: Weakly nonlinear analysis

The oxidation of solutions of glucose with methylene-blue as a catalyst in basic media can induce hydrodynamic overturning instabilities, termed chemoconvection in recognition of their similarity to convective instabilities. The phenomenon is due to gluconic acid, the marginally dense product of the reaction, which gradually builds an unstable density profile. Experiments indicate that dominant pattern wavenumbers initially increase before gradually decreasing or can even oscillate for long times. Here, we perform a weakly nonlinear analysis for an established model of the system with simple kinetics, and show that the resulting amplitude equation is analogous to that obtained in convection with insulating walls. We show that the amplitude description predicts that dominant pattern wavenumbers should decrease in the long term, but does not reproduce the aforementioned increasing wavenumber behavior in the initial stages of pattern development. We hypothesize that this is due to horizontally homogeneous steady states not being attained before pattern onset. We show that the behavior can be explained using a combination of pseudo-steady-state linear and steady-state weakly nonlinear theories. The results obtained are in qualitative agreement with the analysis of experiments.

Nonlinear chemoconvection in the Methylene-blue-glucose system

Interfacial hydrodynamic instabilities arise in a range of chemical systems. One mechanism for instability is the occurrence of unstable density gradients due to the accumulation of reaction products. In this paper we conduct two-dimensional nonlinear numerical simulations for a member of this class of system: the methylene-blue¿glucose reaction. The result of these reactions is the oxidation of glucose to a relatively, but marginally, dense product, gluconic acid, that accumulates at oxygen permeable interfaces, such as the surface open to the atmosphere. The reaction is catalyzed by methylene-blue. We show that simulations help to disassemble the mechanisms responsible for the onset of instability and evolution of patterns, and we demonstrate that some of the results are remarkably consistent with experiments. We probe the impact of the upper oxygen boundary condition, for fixed flux, fixed concentration, or mixed boundary conditions, and find significant qualitative differences in solution behavior; structures either attract or repel one another depending on the boundary condition imposed. We suggest that measurement of the form of the boundary condition is possible via observation of oxygen penetration, and improved product yields may be obtained via proper control of boundary conditions in an engineering setting. We also investigate the dependence on parameters such as the Rayleigh number and depth. Finally, we find that pseudo-steady linear and weakly nonlinear techniques described elsewhere are useful tools for predicting the behavior of instabilities beyond their formal range of validity, as good agreement is obtained with the simulations.

Spatial patterns mediated by subcritical and supercritical thermal fluctuations in the splay Fréedericksz transition

The magnetically induced splay Fréedericksz transition is reexamined to look for pattern forming phenomena slightly above or below criticality. By using our traditional scheme of stochastic nematodynamic equations, situations are, respectively, found of transient and permanent predominance of transversal periodicities (wave numbers) along the direction perpendicular to the initial orientation within the sample. The relevance of these predictions in relation with recent observations in the electrically driven splay Fréedericksz transition, and in general with other pattern forming phenomena, is stressed.

Wave propagation in a medium with disordered excitability

The effect of quenched disorder on the propagation of autowaves in excitable media is studied both experimentally and numerically in relation to the light-sensitive Belousov-Zhabotinsky reaction. The spatial disorder is introduced through a random distribution with two different levels of transmittance. In one dimension the (time-averaged) wave speed is smaller than the corresponding to a homogeneous medium with the mean excitability. Contrarily, in two dimensions the velocity increases due to the roughening of the front. Results are interpreted using kinematic and scaling arguments. In particular, for d = 2 we verify a theoretical prediction of a power-law dependence for the relative change of the propagation speed on the disorder amplitude.

Regular wave propagation out of noise in chemical active media

A pacemaker, regularly emitting chemical waves, is created out of noise when an excitable photosensitive Belousov-Zhabotinsky medium, strictly unable to autonomously initiate autowaves, is forced with a spatiotemporal patterned random illumination. These experimental observations are also reproduced numerically by using a set of reaction-diffusion equations for an activator-inhibitor model, and further analytically interpreted in terms of genuine coupling effects arising from parametric fluctuations. Within the same framework we also address situations of noise-sustained propagation in subexcitable media.

Travelling-stripe forcing generates hexagonal patterns

We study the response of Turing stripe patterns to a simple spatiotemporal forcing. This forcing has the form of a traveling wave and is spatially resonant with the characteristic Turing wavelength. Experiments conducted with the photosensitive chlorine dioxide-iodine-malonic acid reaction reveal a striking symmetry-breaking phenomenon of the intrinsic striped patterns giving rise to hexagonal lattices for intermediate values of the forcing velocity. The phenomenon is understood in the framework of the corresponding amplitude equations, which unveils a complex scenario of dynamical behaviors.

Spatiotemporal order out of noise

Experimental observations of self-organized behavior arising out of noise are also described, and details on the numerical algorithms needed in the computer simulation of these problems are given.

External Fluctuations in Front Propagation

We study the effects of external noise in a one-dimensional model of front propagation. Noise is introduced through the fluctuations of a control parameter leading to a multiplicative stochastic partial differential equation. Analytical and numerical results for the front shape and velocity are presented. The linear-marginal-stability theory is found to increase its range of validity in the presence of external noise. As a consequence noise can stabilize fronts not allowed by the deterministic equation.