22 research outputs found
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
Flow visualization using momentum and energy transport tubes and applications to turbulent flow in wind farms
As a generalization of the mass-flux based classical stream-tube, the concept
of momentum and energy transport tubes is discussed as a flow visualization
tool. These transport tubes have the property, respectively, that no fluxes of
momentum or energy exist over their respective tube mantles. As an example
application using data from large-eddy simulation, such tubes are visualized
for the mean-flow structure of turbulent flow in large wind farms, in fully
developed wind-turbine-array boundary layers. The three-dimensional
organization of energy transport tubes changes considerably when turbine
spacings are varied, enabling the visualization of the path taken by the
kinetic energy flux that is ultimately available at any given turbine within
the array.Comment: Accepted for publication in Journal of Fluid Mechanic
A particle method for the homogeneous Landau equation
We propose a novel deterministic particle method to numerically approximate
the Landau equation for plasmas. Based on a new variational formulation in
terms of gradient flows of the Landau equation, we regularize the collision
operator to make sense of the particle solutions. These particle solutions
solve a large coupled ODE system that retains all the important properties of
the Landau operator, namely the conservation of mass, momentum and energy, and
the decay of entropy. We illustrate our new method by showing its performance
in several test cases including the physically relevant case of the Coulomb
interaction. The comparison to the exact solution and the spectral method is
strikingly good maintaining 2nd order accuracy. Moreover, an efficient
implementation of the method via the treecode is explored. This gives a proof
of concept for the practical use of our method when coupled with the classical
PIC method for the Vlasov equation.Comment: 27 pages, 14 figures, debloated some figures, improved explanations
in sections 2, 3, and
Numerical Study of a Particle Method for Gradient Flows
We study the numerical behaviour of a particle method for gradient flows
involving linear and nonlinear diffusion. This method relies on the
discretisation of the energy via non-overlapping balls centred at the
particles. The resulting scheme preserves the gradient flow structure at the
particle level, and enables us to obtain a gradient descent formulation after
time discretisation. We give several simulations to illustrate the validity of
this method, as well as a detailed study of one-dimensional
aggregation-diffusion equations.Comment: 27 pages, 21 figure
A deterministic particle method for one-dimensional reaction-diffusion equations
We derive a deterministic particle method for the solution of nonlinear reaction-diffusion equations in one spatial dimension. This deterministic method is an analog of a Monte Carlo method for the solution of these problems that has been previously investigated by the author. The deterministic method leads to the consideration of a system of ordinary differential equations for the positions of suitably defined particles. We then consider the time explicit and implicit methods for this system of ordinary differential equations and we study a Picard and Newton iteration for the solution of the implicit system. Next we solve numerically this system and study the discretization error both analytically and numerically. Numerical computation shows that this deterministic method is automatically adaptive to large gradients in the solution
Simulation particulaire de l’écoulement de l’eau dans un sol partiellement saturé
L'équation de Richard décrit l'écoulement de
l'eau dans un sol partiellement saturé. En plus de la non-linéarité, des difficultés
numériques apparaissent suivant les applications : domaine infini et fronts abrupts. La
résolution numérique est encore un challenge et nécessite des techniques spécifiques. La
méthode particulaire semble être intéressante à cause de sa capacité à transporter des
discontinuités et à résoudre naturellement les conditions aux limites externes