Related papers: Non-iterative and exact method for constraining pa…
In this paper, we propose an inexact perturbed path-following algorithm in the framework of Lagrangian dual decomposition for solving large-scale structured convex optimization problems. Unlike the exact versions considered in literature,…
Spurious numerical mixing is a frequent phenomenon in ocean models. In this paper, we present an efficient and robust methodology that defines the vertical grid motion so that this mixing is reduced. This motion is defined as the solution…
We develop a rigid multiblob method for numerically solving the mobility problem for suspensions of passive and active rigid particles of complex shape in Stokes flow in unconfined, partially confined, and fully confined geometries. As in a…
Understanding the macroscopic behavior of dynamical systems is an important tool to unravel transport mechanisms in complex flows. A decomposition of the state space into coherent sets is a popular way to reveal this essential macroscopic…
This paper introduces a mathematical framework for explicit structural dynamics, employing approximate dual functionals and rowsum mass lumping. We demonstrate that the approach may be interpreted as a Petrov-Galerkin method that utilizes…
The numerical approximation of incompressible fluid-structure interaction problems with Lagrange multiplier is generally based on strongly coupled schemes. This delivers unconditional stability but at the expense of solving a…
Locking is a common effect in finite element and isogeometric analysis. In the case of plates, transverse shear locking is most prominent, for shells several other types of locking exist. A common cure are mixed methods that introduce…
The Verlet method is still widely used to integrate the equations of motion in ab initio molecular dynamics simulations. We show that the stability limit of the Verlet method may be significantly increased by setting an upper limit on the…
In this work, we present a high-fidelity and efficient point-particle direct numerical simulation framework based on a multi-block overset curvilinear grid system, enabling large-scale Lagrangian particle tracking in complex geometries with…
Implicit particle-in-cell codes offer advantages over their explicit counterparts in that they suffer weaker stability constraints on the need to resolve the higher frequency modes of the system. This feature may prove particularly valuable…
To address the magnetization dynamics in ferromagnetic materials described by the Landau-Lifshitz-Gilbert equation under large damping parameters, a third-order accurate numerical scheme is developed by building upon a second-order method…
We introduce a particle-based simulation method for granular material in interactive frame rates. We divide the simulation into two decoupled steps. In the first step, a relatively small number of particles is accurately simulated with a…
This paper establishes the iteration-complexity of a Jacobi-type non-Euclidean proximal alternating direction method of multipliers (ADMM) for solving multi-block linearly constrained nonconvex programs. The subproblems of this ADMM variant…
We introduce a data-driven method for learning the equations of motion of mechanical systems directly from position measurements, without requiring access to velocity data. This is particularly relevant in system identification tasks where…
In the present chapter we focus on the fundamentals of non-grid-conforming numerical approaches to simulating particulate flows, implementation issues and grid convergence vs. available reference data. The main idea is to avoid adapting the…
A macroscopic model to describe the dynamics of ion transport in ion channels is the Poisson-Nernst-Planck(PNP) equations. In this paper, we develop a finite-difference method for solving PNP equations, which is second-order accurate in…
We extend the divergence preserving cut finite element method presented in [T. Frachon, P. Hansbo, E. Nilsson, S. Zahedi, SIAM J. Sci. Comput., 46 (2024)] for the Darcy interface problem to unfitted outer boundaries. We impose essential…
A new methodology is developed to integrate numerically the equations of motion for classical many-body systems in molecular dynamics simulations. Its distinguishable feature is the possibility to preserve, independently on the size of the…
Relativistic plasmas around compact objects can sometimes be approximated as being force-free. In this limit, the plasma inertia is negligible and the overall dynamics is governed by global electric currents. We present a novel numerical…
We construct a particle integrator for nonrelativistic particles by means of the splitting method based on the exact flow of the equation of motion of particles in the presence of constant electric and magnetic field. This integrator is…