Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userConvergence and non-convergence of many-particle evolutions with multiple signs
330
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We address the question of convergence of evolving interacting particle systems as the number of particles tends to infinity. We consider two types of particles, called positive and negative. Same-sign particles repel each other, and opposite-sign particles attract each other. The interaction potential is the same for all particles, up to the sign, and has a logarithmic singularity at zero. The central example of such systems is that of dislocations in crystals. Because of the singularity in the interaction potential, the discrete evolution leads to blow-up in finite time. We remedy this situation by regularising the interaction potential at a length-scale $delta_n>0$, which converges to zero as the number of particles $n$ tends to infinity. We establish two main results. The first one is an evolutionary convergence result showing that the empirical measures of the positive and of the negative particles converge to a solution of a set of coupled PDEs which describe the evolution of their continuum densities. In the setting of dislocations these PDEs are known as the Groma-Balogh equations. In the proof we rely on the theory of $lambda$-convex gradient flows, a priori estimates for the Groma-Balogh equations and Orlicz spaces. The proof require $delta_n$ to converge to zero sufficiently slowly. The second result is a counterexample, demonstrating that if $delta_n$ converges to zero sufficiently fast, then the limits of the empirical measures of the positive and the negative dislocations do not satisfy the Groma-Balogh equations. These results show how the validity of the Groma-Balogh equations as the limit of many-particle systems depends in a subtle way on the scale at which the singularity of the potential is regularised.
rate research
Read More
The starting point of our analysis is a class of one-dimensional interacting particle systems with two species. The particles are confined to an interval and exert a nonlocal, repelling force on each other, resulting in a nontrivial equilibrium configuration. This class of particle systems covers the setting of pile-ups of dislocation walls, which is an idealised setup for studying the microscopic origin of several dislocation density models in the literature. Such density models are used to construct constitutive relations in plasticity models. Our aim is to pass to the many-particle limit. The main challenge is the combination of the nonlocal nature of the interactions, the singularity of the interaction potential between particles of the same type, the non-convexity of the the interaction potential between particles of the opposite type, and the interplay between the length-scale of the domain with the length-scale $ell_n$ of the decay of the interaction potential. Our main results are the $Gamma$-convergence of the energy of the particle positions, the evolutionary convergence of the related gradient flows for $ell_n$ sufficiently large, and the non-convergence of the gradient flows for $ell_n$ sufficiently small.
The Athena MHD code has been extended to integrates the motion of particles coupled with the gas via aerodynamic drag, in order to study the dynamics of gas and solids in protoplanetary disks and the formation of planetesimals. Our particle-gas hybrid scheme is based on a second order predictor-corrector method. Careful treatment of the momentum feedback on the gas guarantees exact conservation. The hybrid scheme is stable and convergent in most regimes relevant to protoplanetary disks. We describe a semi-implicit integrator generalized from the leap-frog approach. In the absence of drag force, it preserves the geometric properties of a particle orbit. We also present a fully-implicit integrator that is unconditionally stable for all regimes of particle-gas coupling. Using our hybrid code, we study the numerical convergence of the non-linear saturated state of the streaming instability. We find that gas flow properties are well converged with modest grid resolution (128 cells per pressure length eta r for dimensionless stopping time tau_s=0.1), and equal number of particles and grid cells. On the other hand, particle clumping properties converge only at higher resolutions, and finer resolution leads to stronger clumping before convergence is reached. Finally, we find that measurement of particle transport properties resulted from the streaming instability may be subject to error of about 20%.
We consider an evolution problem associated to the Kazdan-Warner equation on a closed Riemann surface $(Sigma,g)$ begin{align*} -Delta_{g}u=8pileft(frac{he^{u}}{int_{Sigma}he^{u}{rm d}mu_{g}}-frac{1}{int_{Sigma}{rm d}mu_{g}}right) end{align*} where the prescribed function $hgeq0$ and $max_{Sigma}h>0$. We prove the global existence and convergence under additional assumptions such as begin{align*} Delta_{g}ln h(p_0)+8pi-2K(p_0)>0 end{align*} for any maximum point $p_0$ of the sum of $2ln h$ and the regular part of the Green function, where $K$ is the Gaussian curvature of $Sigma$. In particular, this gives a new proof of the existence result by Yang and Zhu [Proc. Amer. Math. Soc. 145 (2017), no. 9, 3953-3959] which generalizes existence result of Ding, Jost, Li and Wang [Asian J. Math. 1 (1997), no. 2, 230-248] to the non-negative prescribed function case.
We study the asymptotic behaviour of a gradient system in a regime in which the driving energy becomes singular. For this system gradient-system convergence concepts are ineffective. We characterize the limiting behaviour in a different way, by proving $Gamma$-convergence of the so-called energy-dissipation functional, which combines the gradient-system components of energy and dissipation in a single functional. The $Gamma$-limit of these functionals again characterizes a variational evolution, but this limit functional is not the energy-dissipation functional of any gradient system. The system in question describes the diffusion of a particle in a one-dimensional double-well energy landscape, in the limit of small noise. The wells have different depth, and in the small-noise limit the process converges to a Markov process on a two-state system, in which jumps only happen from the higher to the lower well. This transmutation of a gradient system into a variational evolution of non-gradient type is a model for how many one-directional chemical reactions emerge as limit of reversible ones. The $Gamma$-convergence proved in this paper both identifies the `fate of the gradient system for these reactions and the variational structure of the limiting irreversible reactions.
In this paper, we focus on the problem of stochastic optimization where the objective function can be written as an expectation function over a closed convex set. We also consider multiple expectation constraints which restrict the domain of the problem. We extend the cooperative stochastic approximation algorithm from Lan and Zhou [2016] to solve the particular problem. We close the gaps in the previous analysis and provide a novel proof technique to show that our algorithm attains the optimal rate of convergence for both optimality gap and constraint violation when the functions are generally convex. We also compare our algorithm empirically to the state-of-the-art and show improved convergence in many situations.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
