19 research outputs found

    Perturbation theory for the effective diffusion constant in a medium of random scatterer

    Full text link
    We develop perturbation theory and physically motivated resummations of the perturbation theory for the problem of a tracer particle diffusing in a random media. The random media contains point scatterers of density ρ\rho uniformly distributed through out the material. The tracer is a Langevin particle subjected to the quenched random force generated by the scatterers. Via our perturbative analysis we determine when the random potential can be approximated by a Gaussian random potential. We also develop a self-similar renormalisation group approach based on thinning out the scatterers, this scheme is similar to that used with success for diffusion in Gaussian random potentials and agrees with known exact results. To assess the accuracy of this approximation scheme its predictions are confronted with results obtained by numerical simulation.Comment: 22 pages, 6 figures, IOP (J. Phys. A. style

    Solving Nonlinear Parabolic Equations by a Strongly Implicit Finite-Difference Scheme

    Full text link
    We discuss the numerical solution of nonlinear parabolic partial differential equations, exhibiting finite speed of propagation, via a strongly implicit finite-difference scheme with formal truncation error O[(Δx)2+(Δt)2]\mathcal{O}\left[(\Delta x)^2 + (\Delta t)^2 \right]. Our application of interest is the spreading of viscous gravity currents in the study of which these type of differential equations arise. Viscous gravity currents are low Reynolds number (viscous forces dominate inertial forces) flow phenomena in which a dense, viscous fluid displaces a lighter (usually immiscible) fluid. The fluids may be confined by the sidewalls of a channel or propagate in an unconfined two-dimensional (or axisymmetric three-dimensional) geometry. Under the lubrication approximation, the mathematical description of the spreading of these fluids reduces to solving the so-called thin-film equation for the current's shape h(x,t)h(x,t). To solve such nonlinear parabolic equations we propose a finite-difference scheme based on the Crank--Nicolson idea. We implement the scheme for problems involving a single spatial coordinate (i.e., two-dimensional, axisymmetric or spherically-symmetric three-dimensional currents) on an equispaced but staggered grid. We benchmark the scheme against analytical solutions and highlight its strong numerical stability by specifically considering the spreading of non-Newtonian power-law fluids in a variable-width confined channel-like geometry (a "Hele-Shaw cell") subject to a given mass conservation/balance constraint. We show that this constraint can be implemented by re-expressing it as nonlinear flux boundary conditions on the domain's endpoints. Then, we show numerically that the scheme achieves its full second-order accuracy in space and time. We also highlight through numerical simulations how the proposed scheme accurately respects the mass conservation/balance constraint.Comment: 36 pages, 9 figures, Springer book class; v2 includes improvements and corrections; to appear as a contribution in "Applied Wave Mathematics II
    corecore