Continued fraction solution for the radiative transfer equation in three dimensions

Starting from the radiative transfer equation, we obtain an analytical solution for both the free propagator along one of the axes and an arbitrary phase function in the Fourier-Laplace domain. We also find the effective absorption parameter, which turns out to be very different from the one provided by the diffusion approximation. We finally present an analytical approximation procedure and obtain a differential equation that accurately reproduces the transport process. We test our approximations by means of simulations that use the Henyey-Greenstein phase function with very satisfactory results.

Noise and dynamics of self-organized critical phenomena

Different microscopic models exhibiting self-organized criticality are studied numerically and analytically. Numerical simulations are performed to compute critical exponents, mainly the dynamical exponent, and to check universality classes. We find that various models lead to the same exponent, but this universality class is sensitive to disorder. From the dynamic microscopic rules we obtain continuum equations with different sources of noise, which we call internal and external. Different correlations of the noise give rise to different critical behavior. A model for external noise is proposed that makes the upper critical dimensionality equal to 4 and leads to the possible existence of a phase transition above d=4. Limitations of the approach of these models by a simple nonlinear equation are discussed.

Solution to telegrapher's equation in the presence of reflecting and partly reflecting broundaries

We show that the reflecting boundary condition for a one-dimensional telegraphers equation is the same as that for the diffusion equation, in contrast to what is found for the absorbing boundary condition. The radiation boundary condition is found to have a quite complicated form. We also obtain exact solutions of the telegraphers equation in the presence of these boundaries.

Telegrapher's equation with variable propagation speeds

All derivations of the one-dimensional telegraphers equation, based on the persistent random walk model, assume a constant speed of signal propagation. We generalize here the model to allow for a variable propagation speed and study several limiting cases in detail. We also show the connections of this model with anomalous diffusion behavior and with inertial dichotomous processes.

Extreme times in financial markets.

We apply the theory of continuous time random walks (CTRWs) to study some aspects involving extreme events in financial time series. We focus our attention on the mean exit time (MET). We derive a general equation for this average and compare it with empirical results coming from high-frequency data of the U.S. dollar and Deutsche mark futures market. The empirical MET follows a quadratic law in the return length interval which is consistent with the CTRW formalism.

Reply to "Comment on 'Solutions of the telegrapher's equation in the presence of traps'"

This reply adds a number of details to remarks by Foong and Kanno [preceding Comment, Phys. Rev. A 46, 5296 (1992)] on our paper [Phys. Rev. A 45, 2222 (1992)] regarding the discontinuities observed in the curves generated in that paper.

When telegrapher's equation furnishes a better approximation to transport equation than the difussion approximation

It has been suggested that a solution to the transport equation which includes anisotropic scattering can be approximated by the solution to a telegrapher's equation [A.J. Ishimaru, Appl. Opt. 28, 2210 (1989)]. We show that in one dimension the telegrapher's equation furnishes an exact solution to the transport equation. In two dimensions, we show that, since the solution can become negative, the telegrapher's equation will not furnish a usable approximation. A comparison between simulated data in three dimensions indicates that the solution to the telegrapher's equation is a good approximation to that of the full transport equation at the times at which the diffusion equation furnishes an equally good approximation.

Generalization of the persistent random walk to dimensions greater than 1

We propose a generalization of the persistent random walk for dimensions greater than 1. Based on a cubic lattice, the model is suitable for an arbitrary dimension d. We study the continuum limit and obtain the equation satisfied by the probability density function for the position of the random walker. An exact solution is obtained for the projected motion along an axis. This solution, which is written in terms of the free-space solution of the one-dimensional telegraphers equation, may open a new way to address the problem of light propagation through thin slabs.

Exact temporal evolution for some non-linear diffusion process

Exact solutions to FokkerPlanck equations with nonlinear drift are considered. Applications of these exact solutions for concrete models are studied. We arrive at the conclusion that for certain drifts we obtain divergent moments (and infinite relaxation time) if the diffusion process can be extended without any obstacle to the whole space. But if we introduce a potential barrier that limits the diffusion process, moments converge with a finite relaxation time.

Strong coupling solution to Schrödinger Equation: the mixing of states

In this paper we give some ideas that can be useful to solve Schrödinger equations in the case when the Hamiltonian contains a large term. We obtain an expansion of the solution in reciprocal powers of the large coupling constant. The procedure followed consists in considering that the small part of the Hamiltonian engenders a motion adiabatic to the motion generated by the large part of the same.

Diffusion in spatially and temporarily inhomogeneous media

In this paper we consider diffusion of a passive substance C in a temporarily and spatially inhomogeneous two-dimensional medium. As a realization for the latter we choose a phase-separating medium consisting of two substances A and B, whose dynamics is determined by the Cahn-Hilliard equation. Assuming different diffusion coefficients of C in A and B, we find that the variance of the distribution function of the said substance grows less than linearly in time. We derive a simple identity for the variance using a probabilistic ansatz and are then able to identify the interface between A and B as the main cause for this nonlinear dependence. We argue that, finally, for very large times the here temporarily dependent diffusion "constant" goes like t-1/3 to a constant asymptotic value D¿. The latter is calculated approximately by employing the effective-medium approximation and by fitting the simulation data to the said time dependence.

Front propagation in the A+B->2A reaction under subdiffusion

We consider an irreversible autocatalytic conversion reaction A+B->2A under subdiffusion described by continuous-time random walks. The reactants transformations take place independently of their motion and are described by constant rates. The analog of this reaction in the case of normal diffusion is described by the Fisher-Kolmogorov-Petrovskii-Piskunov equation leading to the existence of a nonzero minimal front propagation velocity, which is really attained by the front in its stable motion. We show that for subdiffusion, this minimal propagation velocity is zero, which suggests propagation failure.

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.

Nonlinear effects in the dynamics of transient pattern formation in nematics

A nonlinear calculation of the dynamics of transient pattern formation in the Fréedericksz transition is presented. A Gaussian decoupling is used to calculate the time dependence of the structure factor. The calculation confirms the range of validity of linear calculations argued in earlier work. In addition, it describes the decay of the transient pattern.

Integration of local-scale hydrological and regional-scale geophysical based on a nonlinear Bayesian sequential simulation approach

Significant progress has been made with regard to the quantitative integration of geophysical and hydrological data at the local scale. However, extending the corresponding approaches to the regional scale represents a major, and as-of-yet largely unresolved, challenge. To address this problem, we have developed a downscaling procedure based on a non-linear Bayesian sequential simulation approach. The basic objective of this algorithm is to estimate the value of the sparsely sampled hydraulic conductivity at non-sampled locations based on its relation to the electrical conductivity, which is available throughout the model space. The in situ relationship between the hydraulic and electrical conductivities is described through a non-parametric multivariate kernel density function. This method is then applied to the stochastic integration of low-resolution, re- gional-scale electrical resistivity tomography (ERT) data in combination with high-resolution, local-scale downhole measurements of the hydraulic and electrical conductivities. Finally, the overall viability of this downscaling approach is tested and verified by performing and comparing flow and transport simulation through the original and the downscaled hydraulic conductivity fields. Our results indicate that the proposed procedure does indeed allow for obtaining remarkably faithful estimates of the regional-scale hydraulic conductivity structure and correspondingly reliable predictions of the transport characteristics over relatively long distances.