Propagation of a Solitary Fission Wave

Reaction-diffusion phenomena are encountered in an astonishing array of natural systems. Under the right conditions, self stabilizing reaction waves can arise that will propagate at constant velocity. Numerical studies have shown that fission waves of this type are also possible and that they exhibit soliton like properties. Here, we derive the conditions required for a solitary fission wave to propagate at constant velocity. The results place strict conditions on the shapes of the flux, diffusive, and reactive profiles that would be required for such a phenomenon to persist, and this condition would apply to other reaction diffusion phenomena as well. Numerical simulations are used to confirm the results and show that solitary fission waves fall into a bistable class of reaction diffusion phenomena. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4729927]United States Nuclear Regulatory Commission NRC-38-08-946Mechanical Engineerin

Flame Enhancement and Quenching in Fluid Flows

We perform direct numerical simulations (DNS) of an advected scalar field which diffuses and reacts according to a nonlinear reaction law. The objective is to study how the bulk burning rate of the reaction is affected by an imposed flow. In particular, we are interested in comparing the numerical results with recently predicted analytical upper and lower bounds. We focus on reaction enhancement and quenching phenomena for two classes of imposed model flows with different geometries: periodic shear flow and cellular flow. We are primarily interested in the fast advection regime. We find that the bulk burning rate v in a shear flow satisfies v ~ a*U+b where U is the typical flow velocity and a is a constant depending on the relationship between the oscillation length scale of the flow and laminar front thickness. For cellular flow, we obtain v ~ U^{1/4}. We also study flame extinction (quenching) for an ignition-type reaction law and compactly supported initial data for the scalar field. We find that in a shear flow the flame of the size W can be typically quenched by a flow with amplitude U ~ alpha*W. The constant alpha depends on the geometry of the flow and tends to infinity if the flow profile has a plateau larger than a critical size. In a cellular flow, we find that the advection strength required for quenching is U ~ W^4 if the cell size is smaller than a critical value.Comment: 14 pages, 20 figures, revtex4, submitted to Combustion Theory and Modellin

Development of singularities for the compressible Euler equations with external force in several dimensions

We consider solutions to the Euler equations in the whole space from a certain class, which can be characterized, in particular, by finiteness of mass, total energy and momentum. We prove that for a large class of right-hand sides, including the viscous term, such solutions, no matter how smooth initially, develop a singularity within a finite time. We find a sufficient condition for the singularity formation, "the best sufficient condition", in the sense that one can explicitly construct a global in time smooth solution for which this condition is not satisfied "arbitrary little". Also compactly supported perturbation of nontrivial constant state is considered. We generalize the known theorem by Sideris on initial data resulting in singularities. Finally, we investigate the influence of frictional damping and rotation on the singularity formation.Comment: 23 page

The Speed of Fronts of the Reaction Diffusion Equation

We study the speed of propagation of fronts for the scalar reaction-diffusion equation $u_t = u_{xx} + f(u)$\, with $f(0) = f(1) = 0$. We give a new integral variational principle for the speed of the fronts joining the state $u=1$ to $u=0$. No assumptions are made on the reaction term $f(u)$ other than those needed to guarantee the existence of the front. Therefore our results apply to the classical case $f > 0$ in $(0,1)$, to the bistable case and to cases in which $f$ has more than one internal zero in $(0,1)$.Comment: 7 pages Revtex, 1 figure not include

Rarefactions and large time behavior for parabolic equations and monotone schemes

We consider the large time behavior of monotone semigroups associated with degenerate parabolic equations and monotone difference schemes. For an appropriate class of initial data the solution is shown to converge to rarefaction waves at a determined asymptotic rate.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/46470/1/220_2005_Article_BF01229452.pd

Fighting Enemies and Noise: Competition of Residents and Invaders in a Stochastically Fluctuating Environment

The possible control of competitive invasion by infection of the invader and multiplicative noise is studied. The basic model is the Lotka-Volterra competition system with emergent carrying capacities. Several stationary solutions of the non-infected and infected system are identified as well as parameter ranges of bistability. The latter are used for the numerical study of invasion phenomena. The diffusivities, the infection but in particular the white and coloured multiplicative noise are the control parameters. It is shown that not only competition, possible infection and mobilities are important drivers of the invasive dynamics but also the noise and especially its color and the functional response of populations to the emergence of noise

Ruelle–Takens–Newhouse scenario in reaction-diffusion-convection system

Direct numerical simulations of the transition process from periodic to chaotic dynamics are presented for two variable Oregonator-diffusion model coupled with convection. Numerical solutions to the corresponding reaction-diffusion-convection system of equations show that natural convection can change in a qualitative way, the evolution of concentration distribution, as compared with convectionless conditions. The numerical experiments reveal distinct bifurcations as the Grashof number is increased. A transition to chaos similar to Ruelle-Takens-Newhouse scenario is observed. Numerical results are in agreement with the experiments

Global in Time Solutions to Kolmogorov-Feller Pseudodifferential Equations with Small Parameter

The goal in this paper is to demonstrate a new method for constructing global-in-time approximate (asymptotic) solutions of (pseudodifferential) parabolic equations with a small parameter. We show that, in the leading term, such a solution can be constructed by using characteristics, more precisely, by using solutions of the corresponding Hamiltonian system and without using any integral representation. For completeness, we also briefly describe the well-known scheme developed by V.P.Maslov for constructing global-in-time solutions.Comment: 27 page