36 research outputs found
Convergence of a Particle Method and Global Weak Solutions for a Family of Evolutionary PDEs
Abstract: In this talk, I will present global existence and uniqueness results for a family of fluid transport equations by establishing convergence results for the particle method applied to these equations. The considered family of PDEs is a collection of strongly nonlinear equations, which yield traveling wave solutions and can be used to model a variety of flows in fluid dynamics. We apply a particle method to the studied evolutionary equations and provide a new self-contained method for proving its convergence. The latter is accomplished by using the concept of space-time bounded variation and the associated compactness properties. From this result, we prove the existence of a unique global weak solution in some special cases and obtain stronger regularity properties of the solution than previously established
An asymptotic preserving scheme for kinetic models with singular limit
We propose a new class of asymptotic preserving schemes to solve kinetic
equations with mono-kinetic singular limit. The main idea to deal with the
singularity is to transform the equations by appropriate scalings in velocity.
In particular, we study two biologically related kinetic systems. We derive the
scaling factors and prove that the rescaled solution does not have a singular
limit, under appropriate spatial non-oscillatory assumptions, which can be
verified numerically by a newly developed asymptotic preserving scheme. We set
up a few numerical experiments to demonstrate the accuracy, stability,
efficiency and asymptotic preserving property of the schemes.Comment: 24 pages, 6 figure
A Finite-Volume Method for Nonlinear Nonlocal Equations with a Gradient Flow Structure
We propose a positivity preserving entropy decreasing finite volume scheme
for nonlinear nonlocal equations with a gradient flow structure. These
properties allow for accurate computations of stationary states and long-time
asymptotics demonstrated by suitably chosen test cases in which these features
of the scheme are essential. The proposed scheme is able to cope with
non-smooth stationary states, different time scales including metastability, as
well as concentrations and self-similar behavior induced by singular nonlocal
kernels. We use the scheme to explore properties of these equations beyond
their present theoretical knowledge
Integration of the EPDiff equation by particle methods
The purpose of this paper is to apply particle methods to the numerical solution of the EPDiff equation. The weak solutions of EPDiff are contact discontinuities that carry momentum so that wavefront interactions represent collisions in which momentum is exchanged. This behavior allows for the description of many rich physical applications, but also introduces difficult numerical challenges. We present a particle method for the EPDiff equation that is well-suited for this class of solutions and for simulating collisions between wavefronts. Discretization by means of the particle method is shown to preserve the basic Hamiltonian, the weak and variational structure of the original problem, and to respect the conservation laws associated with symmetry under the Euclidean group. Numerical results illustrate that the particle method has superior features in both one and two dimensions, and can also be effectively implemented when the initial data of interest lies on a submanifold
A Diffuse-Domain Based Numerical Method for a Chemotaxis-Fluid Model
In this paper, we consider a coupled chemotaxis-fluid system that models
self-organized collective behavior of oxytactic bacteria in a sessile drop.
This model describes the biological chemotaxis phenomenon in the fluid
environment and couples a convective chemotaxis system for the oxygen-consuming
and oxytactic bacteria with the incompressible Navier-Stokes equations subject
to a gravitational force, which is proportional to the relative surplus of the
cell density compared to the water density.
We develop a new positivity preserving and high-resolution method for the
studied chemotaxis-fluid system. Our method is based on the diffuse-domain
approach, which we use to derive a new chemotaxis-fluid diffuse-domain (cf-DD)
model for simulating bioconvection in complex geometries. The drop domain is
imbedded into a larger rectangular domain, and the original boundary is
replaced by a diffuse interface with finite thickness. The original
chemotaxis-fluid system is reformulated on the larger domain with additional
source terms that approximate the boundary conditions on the physical
interface. We show that the cf-DD model converges to the chemotaxis-fluid model
asymptotically as the width of the diffuse interface shrinks to zero. We
numerically solve the resulting cf-DD system by a second-order hybrid
finite-volume finite-difference method and demonstrate the performance of the
proposed approach on a number of numerical experiments that showcase several
interesting chemotactic phenomena in sessile drops of different shapes, where
the bacterial patterns depend on the droplet geometries
Three-layer approximation of two-layer shallow water equations
Two-layer shallow water equations describe flows that consist of two layers of inviscid fluid of different (constant) densities flowing over bottom topography. Unlike the single-layer shallow water system, the two-layer one is only conditionally hyperbolic: the system loses its hyperbolicity because of the momentum exchange terms between the layers and as a result its solutions may develop instabilities. We study a three-layer approximation of the two-layer shallow water equations by introducing an intermediate layer of a small depth. We examine the hyperbolicity range of the three-layer model and demonstrate that while it still may lose hyperbolicity, the three-layer approximation may improve stability properties of the two-layer shallow water system