Related papers: Vertex Block Descent
We propose a forward-backward splitting dynamical system for solving inclusion problems of the form $0\in A(x)+B(x)$ in Hilbert spaces, where $A$ is a maximal operator and $B$ is a single-valued operator. Involved operators are assumed to…
Accurate simulations of ice sheet dynamics, mantle convection, lava flow, and other highly viscous free-surface flows involve solving the coupled Stokes/free-surface equations. In this paper, we theoretically analyze the stability and…
This paper interprets the stabilized finite element method via residual minimization as a variational multiscale method. We approximate the solution to the partial differential equations using two discrete spaces that we build on a…
We demonstrate theory and computations for finite-energy line defect solutions in an improvement of Ericksen-Leslie liquid crystal theory. Planar director fields are considered in two and three space dimensions, and we demonstrate straight…
The well-posedness of a non-local advection-selection-mutation problem deriving from adaptive dynamics models is shown for a wide family of initial data. A particle method is then developed, in order to approximate the solution of such…
The immersed boundary method is a numerical and mathematical formulation for solving fluid-structure interaction problems. It relies on solving fluid equations on an Eulerian fluid grid and interpolating the resulting velocity back onto…
Multilevel Stein variational gradient descent is a method for particle-based variational inference that leverages hierarchies of surrogate target distributions with varying costs and fidelity to computationally speed up inference. The…
We prove a definitive theorem on the asymptotic stability of point vortex solutions to the full Euler equation in 2 dimensions. More precisely, we show that a small, Gevrey smooth, and compactly supported perturbation of a point vortex…
A novel implicit integration scheme for the Discrete Element Method (DEM) based on the variational integrator approach is presented. The numerical solver provides a fully dynamical description that, notably, reduces to an energy…
Numerical simulation of compressible fluid flows is performed using the Euler equations. They include the scalar advection equation for the density, the vector advection equation for the velocity and a given pressure dependence on the…
Intrusive Uncertainty Quantification methods such as stochastic Galerkin are gaining popularity, whereas the classical stochastic Galerkin approach is not ensured to preserve hyperbolicity of the underlying hyperbolic system. We apply a…
The 2D Euler equations are a simple but rich set of non-linear PDEs that describe the evolution of an ideal inviscid fluid, for which one dimension is negligible. Solving numerically these equations can be extremely demanding. Several…
Gradient descent algorithms perform well in convex optimization but can get tied for finding local minima in non-convex optimization. A robust method that combines a spectral approach with nonmonotone line search strategy for solving…
We describe a computational framework for simulating suspensions of rigid particles in Newtonian Stokes flow. One central building block is a collision-resolution algorithm that overcomes the numerical constraints arising from particle…
In this paper we present a new steepest-descent type algorithm for convex optimization problems. Our algorithm pieces the unknown into sub-blocs of unknowns and considers a partial optimization over each sub-bloc. In quadratic optimization,…
Block coordinate descent is a powerful algorithmic template suitable for big data optimization. This template admits a lot of variants including block gradient descent (BGD), which performs gradient descent on a selected block of variables,…
We present a Virtual Element Method (VEM) for the solution of Dirichlet problems for the quasilinear equation $-\text{div} (k(u)\text{grad} u)=f$ with essential boundary conditions. Within the VEM the nonlinear coefficient is evaluated with…
We develop a lowest-order nonconforming virtual element method for planar linear elasticity, which can be viewed as an extension of the idea in Falk (1991) to the virtual element method (VEM), with the family of polygonal meshes satisfying…
We construct a many-body Gaussian variational approach for the two-dimensional trapped Bose gas in the condensate phase. Interaction between particles is modelized by a generalized pseudo-potential of zero range that allows recovering…
Computational modeling and simulation of fluid-structure interactions constitute a fundamental cornerstone for advancing aerospace engineering endeavors. This paper addresses the notion and implementation of the immersed boundary method for…