Rigorous derivation of a nonlinear diffusion equation as fast-reaction limit of a continuous coagulation-fragmentation model with diffusion

Weak solutions of the spatially inhomogeneous (diffusive) Aizenmann-Bak model of coagulation-breakup within a bounded domain with homogeneous Neumann boundary conditions are shown to converge, in the fast reaction limit, towards local equilibria determined by their mass. Moreover, this mass is the solution of a nonlinear diffusion equation whose nonlinearity depends on the (size-dependent) diffusion coefficient. Initial data are assumed to have integrable zero order moment and square integrable first order moment in size, and finite entropy. In contrast to our previous result [CDF2], we are able to show the convergence without assuming uniform bounds from above and below on the number density of clusters.

Stabilized Petrov-Galerkin methods for the convection-diffusion-reaction and the Helmholtz equations

We present two new stabilized high-resolution numerical methods for the convection–diffusion–reaction (CDR) and the Helmholtz equations respectively. The work embarks upon a priori analysis of some consistency recovery procedures for some stabilization methods belonging to the Petrov–Galerkin framework. It was found that the use of some standard practices (e.g. M-Matrices theory) for the design of essentially non-oscillatory numerical methods is not feasible when consistency recovery methods are employed. Hence, with respect to convective stabilization, such recovery methods are not preferred. Next, we present the design of a high-resolution Petrov–Galerkin (HRPG) method for the 1D CDR problem. The problem is studied from a fresh point of view, including practical implications on the formulation of the maximum principle, M-Matrices theory, monotonicity and total variation diminishing (TVD) finite volume schemes. The current method is next in line to earlier methods that may be viewed as an upwinding plus a discontinuity-capturing operator. Finally, some remarks are made on the extension of the HRPG method to multidimensions. Next, we present a new numerical scheme for the Helmholtz equation resulting in quasi-exact solutions. The focus is on the approximation of the solution to the Helmholtz equation in the interior of the domain using compact stencils. Piecewise linear/bilinear polynomial interpolation are considered on a structured mesh/grid. The only a priori requirement is to provide a mesh/grid resolution of at least eight elements per wavelength. No stabilization parameters are involved in the definition of the scheme. The scheme consists of taking the average of the equation stencils obtained by the standard Galerkin finite element method and the classical finite difference method. Dispersion analysis in 1D and 2D illustrate the quasi-exact properties of this scheme. Finally, some remarks are made on the extension of the scheme to unstructured meshes by designing a method within the Petrov–Galerkin framework.

On the analysis of travelling waves to a nonlinear flux limited reaction-diffusion equation

In this paper we study the existence and qualitative properties of travelling waves associated to a nonlinear flux limited partial differential equation coupled to a Fisher-Kolmogorov-Petrovskii-Piskunov type reaction term. We prove the existence and uniqueness of finite speed moving fronts of C2 classical regularity, but also the existence of discontinuous entropy travelling wave solutions.

Morphological changes in convection-diffusion-limited deposition

The effect of hydrodynamic flow upon diffusion-limited deposition on a line is investigated using a Monte Carlo model. The growth process is governed by the convection and diffusion field. The convective diffusion field is simulated by the biased-random walker resulting from a superimposed drift that represents the convective flow. The development of distinct morphologies is found with varying direction and strength of drift. By introducing a horizontal drift parallel to the deposition plate, the diffusion-limited deposit changes into a single needle inclined to the plate. The width of the needle decreases with increasing strength of drift. The angle between the needle and the plate is about 45° at high flow rate. In the presence of an inclined drift to the plate, the convection-diffusion-limited deposit leads to the formation of a characteristic columnar morphology. In the limiting case where the convection dominates, the deposition process is equivalent to ballistic deposition onto an inclined surface.

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.

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.

Analytical expressions for the flame front speed in the downward combustion of thin solid fuels and comparison to experiments

We derive analytical expressions for the propagation speed of downward combustion fronts of thin solid fuels with a background flow initially at rest. The classical combustion model for thin solid fuels that consists of five coupled reaction-convection-diffusion equations is here reduced into a single equation with the gas temperature as the single variable. For doing so we apply a two-zone combustion model that divides the system into a preheating region and a pyrolyzing region. The speed of the combustion front is obtained after matching the temperature and its derivative at the location that separates both regions.We also derive a simplified version of this analytical expression expected to be valid for a wide range of cases. Flame front velocities predicted by our analyticalexpressions agree well with experimental data found in the literature for a large variety of cases and substantially improve the results obtained from a previous well-known analytical expression

Diffusion in spatially and temporarily inhomogeneous media: Effects of turbulent mixing

We consider diffusion of a passive substance C in a phase-separating nonmiscible binary alloy under turbulent mixing. The substance is assumed to have different diffusion coefficients in the pure phases A and B, leading to a spatially and temporarily dependent diffusion ¿coefficient¿ in the diffusion equation plus convective term. In this paper we consider especially the effects of a turbulent flow field coupled to both the Cahn-Hilliard type evolution equation of the medium and the diffusion equation (both, therefore, supplemented by a convective term). It is shown that the formerly observed prolonged anomalous diffusion [H. Lehr, F. Sagués, and J.M. Sancho, Phys. Rev. E 54, 5028 (1996)] is no longer seen if a flow of sufficient intensity is supplied.

Diffusion in spatially and temporarily inhomogeneous media: Effects of turbulent mixing

We consider diffusion of a passive substance C in a phase-separating nonmiscible binary alloy under turbulent mixing. The substance is assumed to have different diffusion coefficients in the pure phases A and B, leading to a spatially and temporarily dependent diffusion ¿coefficient¿ in the diffusion equation plus convective term. In this paper we consider especially the effects of a turbulent flow field coupled to both the Cahn-Hilliard type evolution equation of the medium and the diffusion equation (both, therefore, supplemented by a convective term). It is shown that the formerly observed prolonged anomalous diffusion [H. Lehr, F. Sagués, and J.M. Sancho, Phys. Rev. E 54, 5028 (1996)] is no longer seen if a flow of sufficient intensity is supplied.

Time-Delayed Theory of the Neolithic Transition in Europe

The classical wave-of-advance model of the neolithic transition (i.e., the shift from hunter-gatherer to agricultural economies) is based on Fisher's reaction-diffusion equation. Here we present an extension of Einstein's approach to Fickian diffusion, incorporating reaction terms. On this basis we show that second-order terms in the reaction-diffusion equation, which have been neglected up to now, are not in fact negligible but can lead to important corrections. The resulting time-delayed model agrees quite well with observations

Front dynamics in the presence of spatio-temporal structured noises

Front dynamics modeled by a reaction-diffusion equation are studied under the influence of spatiotemporal structured noises. An effective deterministic model is analytical derived where the noise parameters, intensity, correlation time, and correlation length appear explicitly. The different effects of these parameters are discussed for the Ginzburg-Landau and Schlögl models. We obtain an analytical expression for the front velocity as a function of the noise parameters. Numerical simulation results are in a good agreement with the theoretical predictions.

Front dynamics in turbulent media

A study of a stable front propagating in a turbulent medium is presented. The front is generated through a reaction-diffusion equation, and the turbulent medium is statistically modeled using a Langevin equation. Numerical simulations indicate the presence of two different dynamical regimes. These regimes appear when the turbulent flow either wrinkles a still rather sharp propagating interfase or broadens it. Specific dependences of the propagating velocities on stirring intensities appropriate to each case are found and fitted when possible according to theoretically predicted laws. Different turbulent spectra are considered.

Front dynamics in turbulent media

A study of a stable front propagating in a turbulent medium is presented. The front is generated through a reaction-diffusion equation, and the turbulent medium is statistically modeled using a Langevin equation. Numerical simulations indicate the presence of two different dynamical regimes. These regimes appear when the turbulent flow either wrinkles a still rather sharp propagating interfase or broadens it. Specific dependences of the propagating velocities on stirring intensities appropriate to each case are found and fitted when possible according to theoretically predicted laws. Different turbulent spectra are considered.

Fronts from integrodifference equations and persistence effects on the Neolithic transition

We extend a previous model of the Neolithic transition in Europe [J. Fort and V. Méndez, Phys. Rev. Lett. 82, 867 (1999)] by taking two effects into account: (i) we do not use the diffusion approximation (which corresponds to second-order Taylor expansions), and (ii) we take proper care of the fact that parents do not migrate away from their children (we refer to this as a time-order effect, in the sense that it implies that children grow up with their parents, before they become adults and can survive and migrate). We also derive a time-ordered, second-order equation, which we call the sequential reaction-diffusion equation, and use it to show that effect (ii) is the most important one, and that both of them should in general be taken into account to derive accurate results. As an example, we consider the Neolithic transition: the model predictions agree with the observed front speed, and the corrections relative to previous models are important (up to 70%)

Description of diffusive and propagative behavior on fractals

The known properties of diffusion on fractals are reviewed in order to give a general outlook of these dynamic processes. After that, we propose a description developed in the context of the intrinsic metric of fractals, which leads us to a differential equation able to describe diffusion in real fractals in the asymptotic regime. We show that our approach has a stronger physical justification than previous works on this field. The most important result we present is the introduction of a dependence on time and space for the conductivity in fractals, which is deduced by scaling arguments and supported by computer simulations. Finally, the diffusion equation is used to introduce the possibility of reaction-diffusion processes on fractals and analyze their properties. Specifically, an analytic expression for the speed of the corresponding travelling fronts, which can be of great interest for application purposes, is derived