Related papers: Non-iterative and exact method for constraining pa…
A new scheme for numerical integration of motion for classical systems composed of rigid polyatomic molecules is proposed. The scheme is based on a matrix representation of the rotational degrees of freedom. The equations of motion are…
In this paper, non-variational systems of differential equations containing small terms are considered, and a consistent approach for deriving approximate conservation laws through the introduction of approximate Lagrange multipliers is…
A novel particle merging algorithm for rarefied gas dynamics simulations is proposed that can conserve arbitrary velocity and spatial moments of the particle distribution via solving a non-negative least squares problem. An extension that…
We describe and summarize a class of minimal numerical models emerged from recent development of simulation methods for dense particle suspensions in overdamped linear flows. The main ingredients include (i) a frame-invariant, short-range…
We study a finite element computational model for solving the coupled problem arising in the interaction between a free fluid and a fluid in a poroelastic medium. The free fluid is governed by the Stokes equations, while the flow in the…
We consider the parallel-in-time solution of scalar nonlinear conservation laws in one spatial dimension. The equations are discretized in space with a conservative finite-volume method using weighted essentially non-oscillatory (WENO)…
We present a new particle-based (discrete element) numerical method for the simulation of granular dynamics, with application to motions of particles on small solar system body and planetary surfaces. The method employs the parallel N-body…
In this article we propose a novel strategy for choosing the Lagrange multipliers in the Levenberg-Marquardt method for solving ill-posed problems modeled by nonlinear operators acting between Hilbert spaces. Convergence analysis results…
The simulation of large nonlinear dynamical systems, including systems generated by discretization of hyperbolic partial differential equations, can be computationally demanding. Such systems are important in both fluid and kinetic…
We introduce a new formulation for the finite element immersed boundary method which makes use of a distributed Lagrange multiplier. We prove that a full discretization of our model, based on a semi-implicit time advancing scheme, is…
Problems involving rolling without slipping or no sideways skidding, to name a few, introduce velocity-dependent constraints that can be efficiently treated by the method of Lagrange multipliers in the Lagrangian formulation of the…
Two are the main objectives of this article: first, we introduce a method for determining and analyzing constrained local extrema that provides a different alternative to all previous works on the topic, by eliminating Lagrange multipliers…
Systems subjected to holonomic constraints follow quite complicated dynamics that could not be described easily with Hamiltonian or Lagrangian dynamics. The influence of holonomic constraints in equations of motions is taken into account by…
A geometrically exact membrane formulation is presented that is based on curvilinear coordinates and isogeometric finite elements, and is suitable for both solid and liquid membranes. The curvilinear coordinate system is used to describe…
We present a novel deep learning approach to approximate the solution of large, sparse, symmetric, positive-definite linear systems of equations. These systems arise from many problems in applied science, e.g., in numerical methods for…
Variational integrators for Lagrangian dynamical systems provide a systematic way to derive geometric numerical methods. These methods preserve a discrete multisymplectic form as well as momenta associated to symmetries of the Lagrangian…
A numerical algorithm for solving mantle convection problems with strongly variable viscosity is presented. Equations for conservation of mass and momentum for highly viscous and incompressible fluids are solved iteratively by a multigrid…
Variational integrators applied to degenerate Lagrangians that are linear in the velocities are two-step methods. The system of modified equations for a two-step method consists of the principal modified equation and one additional equation…
We present an algorithm for the numerical solution of the equations governing combustion in porous inert media. The discretization of the flow problem is performed by the mixed finite element method, the transport problems are discretized…
We investigate a spatially discrete surrogate model for the dynamics of a slender, elastic, inextensible fiber in turbulent flows. Deduced from a continuous space-time beam model for which no solution theory is available, it consists of a…