oalib
Search Results: 1 - 10 of 100 matches for " "
All listed articles are free for downloading (OA Articles)
Page 1 /100
Display every page Item
An Asymptotic Preserving Scheme for the Diffusive Limit of Kinetic systems for Chemotaxis  [PDF]
Jose A. Carrillo,Bokai Yan
Mathematics , 2011,
Abstract: In this work we numerically study the diffusive limit of run & tumble kinetic models for cell motion due to chemotaxis by means of asymptotic preserving schemes. It is well-known that the diffusive limit of these models leads to the classical Patlak-Keller-Segel macroscopic model for chemotaxis. We will show that the proposed scheme is able to accurately approximate the solutions before blow-up time for small parameter. Moreover, the numerical results indicate that the global solutions of the kinetic models stabilize for long times to steady states for all the analyzed parameter range. We also generalize these asymptotic preserving schemes to two dimensional kinetic models in the radial case. The blow-up of solutions is numerically investigated in all these cases.
On the Asymptotic Preserving property of the Unified Gas Kinetic Scheme for the diffusion limit of linear kinetic models  [PDF]
Luc Mieussens
Mathematics , 2013, DOI: 10.1016/j.jcp.2013.07.002
Abstract: The unified gas kinetic scheme (UGKS) of K. Xu et al. [K. Xu and J.-C. Huang, J. Comput. Phys., 229, pp. 7747--7764, 2010], originally developed for multiscale gas dynamics problems, is applied in this paper to a linear kinetic model of radiative transfer theory. While such problems exhibit purely diffusive behavior in the optically thick (or small Knudsen) regime, we prove that UGKS is still asymptotic preserving (AP) in this regime, but for the free transport regime as well. Moreover, this scheme is modified to include a time implicit discretization of the limit diffusion equation, and to correctly capture the solution in case of boundary layers. Contrary to many AP schemes, this method is based on a standard finite volume approach, it does neither use any decomposition of the solution, nor staggered grids. Several numerical tests demonstrate the properties of the scheme.
A Comparative Study of an Asymptotic Preserving Scheme and Unified Gas-kinetic Scheme in Continuum Flow Limit  [PDF]
Songze Chen,Kun Xu
Physics , 2013, DOI: 10.1016/j.jcp.2015.02.014
Abstract: Asymptotic preserving (AP) schemes are targeting to simulate both continuum and rarefied flows. Many AP schemes have been developed and are capable of capturing the Euler limit in the continuum regime. However, to get accurate Navier-Stokes solutions is still challenging for many AP schemes. In order to distinguish the numerical effects of different AP schemes on the simulation results in the continuum flow limit, an implicit-explicit (IMEX) AP scheme and the unified gas kinetic scheme (UGKS) based on Bhatnagar-Gross-Krook (BGk) kinetic equation will be applied in the flow simulation in both transition and continuum flow regimes. As a benchmark test case, the lid-driven cavity flow is used for the comparison of these two AP schemes. The numerical results show that the UGKS captures the viscous solution accurately. The velocity profiles are very close to the classical benchmark solutions. However, the IMEX AP scheme seems have difficulty to get these solutions. Based on the analysis and the numerical experiments, it is realized that the dissipation of AP schemes in continuum limit is closely related to the numerical treatment of collision and transport of the kinetic equation. Numerically it becomes necessary to couple the convection and collision terms in both flux evaluation at a cell interface and the collision source term treatment inside each control volume.
Asymptotic-preserving projective integration schemes for kinetic equations in the diffusion limit  [PDF]
Pauline Lafitte,Giovanni Samaey
Mathematics , 2010, DOI: 10.1137/100795954
Abstract: We investigate a projective integration scheme for a kinetic equation in the limit of vanishing mean free path, in which the kinetic description approaches a diffusion phenomenon. The scheme first takes a few small steps with a simple, explicit method, such as a spatial centered flux/forward Euler time integration, and subsequently projects the results forward in time over a large time step on the diffusion time scale. We show that, with an appropriate choice of the inner step size, the time-step restriction on the outer time step is similar to the stability condition for the diffusion equation, whereas the required number of inner steps does not depend on the mean free path. We also provide a consistency result. The presented method is asymptotic-preserving, in the sense that the method converges to a standard finite volume scheme for the diffusion equation in the limit of vanishing mean free path. The analysis is illustrated with numerical results, and we present an application to the Su-Olson test.
An asymptotic preserving scheme for kinetic models for chemotaxis phenomena  [PDF]
Abdelghani Bellouquid,Jacques Tagoudjeu
Mathematics , 2015,
Abstract: In this paper, we propose a numerical scheme to solve the kinetic model for chemotaxis phenomena. Formally, this scheme is shown to be uniformly stable with respect to the small parameter, consistent with the fluid-diffusion limit (Keller-Segel model). Our approach is based on the micro-macro decomposition which leads to an equivalent formulation of the kinetic model that couples a kinetic equation with macroscopic ones. This method is validated with various test cases and compared to other standard methods.
An asymptotic preserving scheme based on a new formulation for NLS in the semiclassical limit  [PDF]
Christophe Besse,Rémi Carles,Florian Méhats
Mathematics , 2012, DOI: 10.1137/120899017
Abstract: We consider the semiclassical limit for the nonlinear Schrodinger equation. We introduce a phase/amplitude representation given by a system similar to the hydrodynamical formulation, whose novelty consists in including some asymptotically vanishing viscosity. We prove that the system is always locally well-posed in a class of Sobolev spaces, and globally well-posed for a fixed positive Planck constant in the one-dimensional case. We propose a second order numerical scheme which is asymptotic preserving. Before singularities appear in the limiting Euler equation, we recover the quadratic physical observables as well as the wave function with mesh size and time step independent of the Planck constant. This approach is also well suited to the linear Schrodinger equation.
A high-order asymptotic-preserving scheme for kinetic equations using projective integration  [PDF]
Pauline Lafitte,Annelies Lejon,Giovanni Samaey
Mathematics , 2014,
Abstract: We investigate a high-order, fully explicit, asymptotic-preserving scheme for a kinetic equation with linear relaxation, both in the hydrodynamic and diffusive scalings in which a hyperbolic, resp. parabolic, limiting equation exists. The scheme first takes a few small (inner) steps with a simple, explicit method (such as direct forward Euler) to damp out the stiff components of the solution and estimate the time derivative of the slow components. These estimated time derivatives are then used in an (outer) Runge-Kutta method of arbitrary order. We show that, with an appropriate choice of inner step size, the time-step restriction on the outer time step is similar to the stability condition for the limiting macroscopic equation. Moreover, the number of inner time steps is also independent of the scaling parameter. We analyse stability and consistency, and illustrate with numerical results.
Hybrid model for the Coupling of an Asymptotic Preserving scheme with the Asymptotic Limit model: The One Dimensional Case  [cached]
Degond Pierre,Deluzet Fabrice,Maldarella Dario,Narski Jacek
ESAIM : Proceedings , 2011, DOI: 10.1051/proc/2011010
Abstract: In this paper a strategy is investigated for the spatial coupling of an asymptotic preserving scheme with the asymptotic limit model, associated to a singularly perturbed, highly anisotropic, elliptic problem. This coupling strategy appears to be very advantageous as compared with the numerical discretization of the initial singular perturbation model or the purely asymptotic preserving scheme introduced in previous works [3, 5]. The model problem addressed in this paper is well suited for the simulation of a plasma in the presence of a magnetic field, whose intensity may vary considerably within the simulation domain.
An asymptotic-preserving scheme for the semiconductor Boltzmann equation toward the energy-transport limit  [PDF]
Jingwei Hu,Li Wang
Mathematics , 2013,
Abstract: We design an asymptotic-preserving scheme for the semiconductor Boltzmann equation which leads to an energy-transport system for electron mass and internal energy as mean free path goes to zero. To overcome the stiffness induced by the convection terms, we adopt an even-odd decomposition to formulate the equation into a diffusive relaxation system. New difficulty arises in the two-scale stiff collision terms, whereas the simple BGK penalization does not work well to drive the solution to the correct limit. We propose a clever variant of it by introducing a threshold on the stiffer collision term such that the evolution of the solution resembles a Hilbert expansion at the continuous level. Formal asymptotic analysis and numerical results are presented to illustrate the efficiency and accuracy of the new scheme.
High Order Asymptotic Preserving DG-IMEX Schemes for Discrete-Velocity Kinetic Equations in a Diffusive Scaling  [PDF]
Juhi Jang,Fengyan Li,Jing-Mei Qiu,Tao Xiong
Mathematics , 2013,
Abstract: In this paper, we develop a family of high order asymptotic preserving schemes for some discrete-velocity kinetic equations under a diffusive scaling, that in the asymptotic limit lead to macroscopic models such as the heat equation, the porous media equation, the advection-diffusion equation, and the viscous Burgers equation. Our approach is based on the micro-macro reformulation of the kinetic equation which involves a natural decomposition of the equation to the equilibrium and non-equilibrium parts. To achieve high order accuracy and uniform stability as well as to capture the correct asymptotic limit, two new ingredients are employed in the proposed methods: discontinuous Galerkin spatial discretization of arbitrary order of accuracy with suitable numerical fluxes; high order globally stiffly accurate implicit-explicit Runge-Kutta scheme in time equipped with a properly chosen implicit-explicit strategy. Formal asymptotic analysis shows that the proposed scheme in the limit of epsilon -> 0 is an explicit, consistent and high order discretization for the limiting equation. Numerical results are presented to demonstrate the stability and high order accuracy of the proposed schemes together with their performance in the limit.
Page 1 /100
Display every page Item


Home
Copyright © 2008-2017 Open Access Library. All rights reserved.