Exclusion Process with Slow Boundary

We study the hydrodynamic and the hydrostatic behavior of the simple symmetric exclusion process with slow boundary. The term slow boundary means that particles can be born or die at the boundary sites, at a rate proportional to \(N^{-\theta }\), where \(\theta > 0\) and N is the scaling parameter. In the bulk, the particles exchange rate is equal to 1. In the hydrostatic scenario, we obtain three different linear profiles, depending on the value of the parameter \(\theta \); in the hydrodynamic scenario, we obtain that the time evolution of the spatial density of particles, in the diffusive scaling, is given by the weak solution of the heat equation, with boundary conditions that depend on \( \theta \). If \(\theta \in (0,1)\), we get Dirichlet boundary conditions, (which is the same behavior if \(\theta =0\), see Farfán in Hydrostatics, statical and dynamical large deviations of boundary driven gradient symmetric exclusion processes, 2008); if \(\theta =1\), we get Robin boundary conditions; and, if \(\theta \in (1,\infty )\), we get Neumann boundary conditions.     

Dynamics of one- and two-dimensional kinks in bistable reaction-diffusion equations with quasidiscrete sources of reaction

We study the evolution of fronts in a bistable reaction-diffusion system when the nonlinear reaction term is spatially inhomogeneous. This equation has been used to model wave propagation in various biological systems. Extending previous works on homogeneous reaction terms, we derive asymptotically an equation governing the front motion, which is strongly nonlinear and, for the two-dimensional case, generalizes the classical mean curvature flow equation. We study the motion of one- and two-dimensional fronts, finding that the inhomogeneity acts as a "potential function" for the motion of the front; i.e., there is wave propagation failure and the steady state solution depends on the structure of the function describing the inhomogeneity. (c) 2001 American Institute of Physics.     

Calculation of superpropagators in non-linear quantum field theories

A new method of constructing the superpropagators, i.e. the Fourier transforms of the expressions of the form is suggested. The method makes it possible to derive by use of the same technique explicit analytic expressions for the superpropagators for a wide class of field theories — from strictly local up to essentially non-local. The essence of the method is the construction of a differential equation for the superpropagator which in general is of an infinite order. By use of the boundary condition atp ²=0 we find the solution of this equation depending on one arbitrary real parameter. Simple examples are given to illustrate the method.     

Charging of aerosol particles in the near free-molecule regime

The charging of small neutral and charged particles suspended in weakly ionized plasma is investigated under the assumption that the Coulomb + image forces give rise to the ion transport in the carrier plasma and define the rate of charging processes. Our approach is based on a BGK version of the kinetic equation [1,2] describing the ion transport in the presence of force fields created by the particle charge and the image force. A special type of the perturbation theory (with respect to the reciprocal Knudsen number) is used for calculating the rate of ion deposition onto neutral and charged particles. As the starting approximation, the free-molecule ion distribution with a floating ion flux is used for evaluating the collision term in the Boltzmann equation. The value of the ion flux as a function of the particle size is then fixed self-consistently from the solution of the Boltzmann equation with the approximated collision term. The expression for the ion flux J(a) to the spherical particle of radius a is derived in the form , where J_fm is the free-molecule flux (no carrier plasma) and is a correction factor taking into account the ion-molecular collisions. The latter is shown to never exceed unity and to depend weakly on the particle-ion interaction.Received: 29 December 2003, Published online: 15 April 2004PACS: 36.40.Wa Charged clusters - 82.30.Fi Ion-molecule, ion-ion, and charge-transfer reactions - 92.60.Mt Particles and aerosols     

Pomraning-Eddington approximation for radiative transfer in a spherical turbid medium

Abstract

The Pomraning-Eddington approximation is used to solve the radiative transfer problem for anisotropic scattering in a spherical homogeneous turbid medium with diffuse and specular reflecting boundaries. This approximation replaces the radiative transfer integro-differential equation by a second-order differential equation which has an analytical solution in terms of the modified Bessel function. Here, we calculate the partial heat flux at the boundary of anisotropic scattering on a homogeneous solid sphere. The calculations are carried out for spherical media of radii 0.1, 1.0 and 10 mfp and for scattering albedos between 0.1 and 1.0. In addition, the calculations are given for media with transparent, diffuse reflecting and diffuse and specular reflecting boundaries. Two different weight functions are used to verify the boundary conditions. Our results are compared with those given by the Galerkin technique and show greater accuracy for thick and highly scattering media.     

Exact and asymptotic solutions of the Cauchy–Poisson problem with localized initial conditions and a constant function of the bottom

In the paper, the asymptotic solutions for a problem of Cauchy–Poisson type with localized initial conditions are constructed. The bottom of the basin under consideration which was constant before the perturbation, is instantly perturbed at the initial time moment by a spatially localized function. Simplifications of the corresponding formulas are presented inside and outside the vicinity of the leading front, as well as in the case of a special choice of the initial condition. It is shown that, in the vicinity of the leading front, the asymptotic solution coincides with the asymptotic solution of the linear Boussinesq equation.     

Narrow Escape, Part III: Non-Smooth Domains and Riemann Surfaces

We consider the narrow escape problem in two-dimensional Riemannian manifolds (with a metric g) with corners and cusps, in an annulus, and on a sphere. Specifically, we calculate the mean time it takes a Brownian particle diffusing in a domain Ω to reach an absorbing window when the ratio between the absorbing window and the otherwise reflecting boundary is small. If the boundary is smooth, as in the cases of the annulus and the sphere, the leading term in the expansion is the same as that given in part I of the present series of papers, however, when it is not smooth, the leading order term is different. If the absorbing window is located at a corner of angle α, then if near a cusp, then grows algebraically, rather than logarithmically. Thus, in the domain bounded between two tangent circles, the expected lifetime is , where is the ratio of the radii. For the smooth boundary case, we calculate the next term of the expansion for the annulus and the sphere. It can also be evaluated for domains that can be mapped conformally onto an annulus. This term is needed in real life applications, such as trafficking of receptors on neuronal spines, because is not necessarily large, even when is small. In these two problems there are additional parameters that can be small, such as the ratio δ of the radii of the annulus. The contributions of these parameters to the expansion of the mean escape time are also logarithmic. In the case of the annulus the mean escape time is .     

Ricci curvature non-minimal derivative coupling cosmology with field re-scaling

In this letter, cosmology of a simple NMDC gravity with \(\xi R \phi _{,\mu }\phi ^{,\mu }\) term and a free kinetic term is considered in flat geometry and in presence of dust matter. A logarithm field transformation \(\phi ' = \mu \ln \phi \) is proposed phenomenologically. Assuming slow-roll approximation, equation of motion, scalar field solution and potential are derived as function of kinematic variables. The field solution and potential are found straightforwardly for power-law, de-Sitter and super-acceleration expansions. Slow-roll parameters and slow-roll condition are found to depend on more than one variable. At large field the re-scaling effect can enhance the acceleration. For slow-rolling field, the negative coupling \(\xi \) could enhance the effect of acceleration.     

Uniform Regularity Close to Cross Singularities in an Unstable Free Boundary Problem

We introduce a new method for the analysis of singularities in the unstable problem $$ \Delta u = -\chi_{\{u >0 \}}, $$ which arises in solid combustion as well as in the composite membrane problem. Our study is confined to points of “supercharacteristic” growth of the solution, i.e. points at which the solution grows faster than the characteristic/invariant scaling of the equation would suggest. At such points the classical theory is doomed to fail, due to incompatibility of the invariant scaling of the equation and the scaling of the solution. In the case of two dimensions our result shows that in a neighborhood of the set at which the second derivatives of u are unbounded, the level set {u = 0} consists of two C ¹-curves meeting at right angles. It is important that our result is not confined to the minimal solution of the equation but holds for all solutions.     

Interface Fluctuations for the D = 1 Stochastic Ginzburg–Landau Equation with Nonsymmetric Reaction Term

We consider a Ginzburg–Landau equation in the interval [?ε^?1, ε^?1], ε>0, with Neumann boundary conditions, perturbed by an additive white noise of strength $\sqrt {\varepsilon } $ and reaction term being the derivative of a function which has two equal–depth wells at ±1, but is not symmetric. When ε=0, the equation has equilibrium solutions that are increasing, and connect ?1 with +1. We call them instantons, and we study the evolution of the solutions of the perturbed equation in the limit ε→0⁺, when the initial datum is close to an instanton. We prove that, for times that may be of the order of ε^?1, the solution stays close to some instanton whose center, suitably normalized, converges to a Brownian motion plus a drift. This drift is known to be zero in the symmetric case, and, using a perturbative analysis, we show that if the nonsymmetric part of the reaction term is sufficiently small, it determines the sign of the drift.     

Metastability for the Exclusion Process with Mean-Field Interaction

We consider an exclusion particle system with long-range, mean-field-type interactions at temperature 1/β. The hydrodynamic limit of such a system is given by an integrodifferential equation with one conservation law on the circle $C$ : it is the gradient flux of the Kac free energy functional F _β. For β≤1, any constant function with value m ∈ [?1, +1] is the global minimizer of F _β in the space $\{ u:\int_C {u(x)} \,dx = m\} $ . For β>1, F _β restricted to $\{ u:\int_C {u(x)} \,dx = m\} $ may have several local minima: in particular, the constant solution may not be the absolute minimizer of F _β. We therefore study the long-time behavior of the particle system when the initial condition is close to a homogeneous stable state, giving results on the time of exit from (suitable) subsets of its domain of attraction. We follow the Freidlin–Wentzell approach: first, we study in detail F _β together with the time asymptotics of the solution of the hydrodynamic equation; then we study the probability of rare events for the particle system, i.e., large deviations from the hydrodynamic limit.     

The Enskog Equation for Confined Elastic Hard Spheres

A kinetic equation for a system of elastic hard spheres or disks confined by a hard wall of arbitrary shape is derived. It is a generalization of the modified Enskog equation in which the effects of the confinement are taken into account and it is supposed to be valid up to moderate densities. From the equation, balance equations for the hydrodynamic fields are derived, identifying the collisional transfer contributions to the pressure tensor and heat flux. A Lyapunov functional, \({\mathcal {H}}[f]\), is identified. For any solution of the kinetic equation, \({\mathcal {H}}\) decays monotonically in time until the system reaches the inhomogeneous equilibrium distribution, that is a Maxwellian distribution with a density field consistent with equilibrium statistical mechanics.     

Pomraning-Eddington approximation for radiative transfer in a spherical turbid medium

The Pomraning-Eddington approximation is used to solve the radiative transfer problem for anisotropic scattering in a spherical homogeneous turbid medium with diffuse and specular reflecting boundaries. This approximation replaces the radiative transfer integro-differential equation by a second-order differential equation which has an analytical solution in terms of the modified Bessel function. Here, we calculate the partial heat flux at the boundary of anisotropic scattering on a homogeneous solid sphere. The calculations are carried out for spherical media of radii 0.1, 1.0 and 10 mfp and for scattering albedos between 0.1 and 1.0. In addition, the calculations are given for media with transparent, diffuse reflecting and diffuse and specular reflecting boundaries. Two different weight functions are used to verify the boundary conditions. Our results are compared with those given by the Galerkin technique and show greater accuracy for thick and highly scattering media.     

The Largest Fragment of a Homogeneous Fragmentation Process

We show that in homogeneous fragmentation processes the largest fragment at time t has size

$$\begin{aligned} e^{-t \Phi '(\overline{p})}t^{-\frac{3}{2} (\log \Phi )'(\overline{p})+o(1)}, \end{aligned}$$

where \(\Phi \) is the Lévy exponent of the fragmentation process, and \(\overline{p}\) is the unique solution of the equation \((\log \Phi )'(\bar{p})=\frac{1}{1+\bar{p}}\). We argue that this result is in line with predictions arising from the classification of homogeneous fragmentation processes as logarithmically correlated random fields.

     

Scalar conservation laws with discontinuous flux function: I. The viscous profile condition

The equation , whereH is Heaviside's step function, appears for example in continuous sedimentation of solid particles in a liquid, in two-phase flow, in traffic-flow analysis and in ion etching. The discontinuity of the flux function atx=0 causes a discontinuity of a solution, which is not uniquely determined by the initial data. The equation can be written as a triangular 2×2 non-strictly hyperbolic system. This augmentation is non-unique and a natural definition is given by means of viscous profiles. By a viscous profile we mean a stationary solution ofu _t +(F ) _x =u _xx , whereF is a smooth approximation of the discontinuous flux, i.e.,H is smoothed. In terms of the 2×2 system, the discontinuity atx=0 is either a regular Lax, an under-or overcompressive, a marginal under- or overcompressive or a degenerate shock wave. In some cases, depending onf andg, there is a unique viscous profile (e.g. undercompressive and regular Lax waves) and in some cases there are infinitely many (e.g. overcompressive waves). The main purpose of the paper is to show the equivalence between a previously introduced uniqueness condition for the discontinuity of the solution atx=0 and the viscous profile condition.     

Vacuum Boundary Effects

The Casimir effect is highly dependent on the shape and structure of space boundaries. This dependence is encoded in the variation of vacuum energy with the different types of boundary conditions. We analyze from a global perspective the properties of the Casimir energy as a function on the largest space of the consistent boundary conditions M_F\mathcal{M}_{F} for a massless scalar field confined between to homogeneous parallel plates. In particular, we analyze the analytic properties of this function and point out the existence of a third order phase transition at periodic boundary conditions. We also characterize the boundary conditions which give rise to attractive or repulsive Casimir forces. In the interface between both regimes we find a very interesting family of boundary conditions without Casimir effect, and fully characterize the boundary conditions which do not induce any type of Casimir force.     

An elastodynamic solution of finite long orthotropic hollow cylinder under torsion impact

The paper presents an analytical method to solve the elastodynamic problem of a finite-length orthotropic hollow cylinder subjected to a torsion impact often occurring in engineering fields. The elastodynamic solution is composed of a quasi-static solution of homogeneous equation satisfied with the non-homogeneous boundary condition and a dynamic solution of non-homogeneous equation satisfied with homogeneous boundary condition. The quasi-static solution can be obtained by directly solving the quasi-static equation satisfied with the non-homogeneous boundary condition. The solution of a non-homogeneous dynamic equation is obtained by means of a finite Hankel transform to a radial variable r, Laplace transform to a time variable t and finite Fourier transform to an axial variable z. Thus, the elastodynamic solution of the finite length of an orthotropic hollow cylinder subjected to a torsion impact is obtained. On the other hand, a dynamic finite element for the same problem is also carried out by applying the ANSYS finite-element analysis system. Comparing the theoretical solution with finite-element solution, it can be found that two kinds of results obtained by making use of two different solving methods are suitably approached. Therefore, it is further concluded that the methods and computing processes of the theoretical solution are effective and accurate.     

Beam vibrations with time-dependent boundary conditions

A procedure is outlined for the solution of the vibration problem of a Bernoulli-Euler beam with time-dependent boundary conditions. The solution is greatly simplified if the dependent variable in the original partial differential equation can be changed to produce homogeneous boundary conditions and at the same time maintain a homogeneous differential equation. A method for making such a change is given and illustrated by solving a cantilever beam problem with a time-dependent tip displacement.     

Field of a point source in a semi-infinite elastic medium

An exact solution for the tensor Green's function of a harmonic field for a semi-infinite elastic medium is presented in an easy-to-use form in the theory of wave scattering. The solution is derived in the form of a sum of the Green's functions for an infinite medium and the term satisfying the homogeneous wave equation for a semi-infinite elastic medium. The results reproduce the known far-field asymptotics containing longitudinal, transversal and surface Rayleigh-type wave modes. The near-field asymptotic is essentially different for the regions far and near the boundary.