Related papers: A new formulation for coupled magma/mantle dynamic…
Targeting simulations on parallel hardware architectures, this paper presents computational kernels for efficient computations in mortar finite element methods. Mortar methods enable a variationally consistent imposition of coupling…
In order to meet the requirements of practical applications, a model of deforming manifold in the embedded space is proposed. The deforming vector and deforming field are presented to precisely describe the deforming process, which have…
Dislocation dynamic is a typically gradient flow problem, and most of work solves it just as ODE, which means that the interacting energy of dislocations is ignored. We take the interaction energy into account and use it to introduce new…
We describe a novel splitting approach to numerical relativistic magnetohydrodynamics (RMHD) designed to expand its applicability to the domain of ultra-high magnetisation (high-$\sigma$). In this approach, the electromagnetic field is…
This study focuses on solving the numerical challenges of imposing absorbing boundary conditions for dynamic simulations in the material point method (MPM). To attenuate elastic waves leaving the computational domain, the current work…
Contraction analysis establishes exponential incremental convergence of a nonlinear system by solving a linear matrix inequality for a contraction metric, and has become a standard resource for solving problems in nonlinear control and…
While 3D Gaussian Splatting (3DGS) excels in static scene modeling, its extension to dynamic scenes introduces significant challenges. Existing dynamic 3DGS methods suffer from either over-smoothing due to low-rank decomposition or feature…
This study presents an innovative direct numerical simulation approach for complex particle systems with irregular shapes and large numbers. Using partially saturated methods, it accurately models arbitrary shapes, albeit at considerable…
Granular impact -- the dynamic intrusion of solid objects into granular media -- is widespread across scientific and engineering applications including geotechnics. Existing approaches for simulating granular impact dynamics have relied on…
We study the fully mixed formulation of the Biot equations, which is characterized by a symmetric coupling between flow and deformation. This structure enables the use of stable mixed finite elements for each subproblem without a strong…
Accurate and robust modelling of large deformation three dimensional contact interaction is an important area of engineering, but it is also challenging from a computational mechanics perspective. This is particularly the case when there is…
The Oxygen Depletion problem is an implicit free boundary value problem. The dynamics allow topological changes in the free boundary. We show several mathematical formulations of this model from the literature and give a new formulation…
A new algorithm for efficient and fully time-reversible integration of first-principles molecular dynamics based on orbital-free density functional theory (OFDFT) is presented. The algorithm adapts to this nontrivial case the recently…
The mixing of magmas is a fundamental process in the Earth system causing extreme compositional variations in igneous rocks. This process can develop with different intensities both in space and time, making the interpretation of…
An analysis of the dynamics is performed, of exactly solvable models for fragile and strong glasses, exploiting the partitioning of the free energy landscape in inherent structures. The results are compared with the exact solution of the…
We consider the scientifically challenging and policy-relevant task of understanding the past and projecting the future dynamics of the Antarctic ice sheet. The Antarctic ice sheet has shown a highly nonlinear threshold response to past…
We introduce a new overlapping Domain Decomposition Method (DDM) to solve the fully nonlinear Monge-Amp\`ere equation. While DDMs have been extensively studied for linear problems, their application to fully nonlinear partial differential…
Atmosphere modelling applications become increasingly memory-bound due to the inconsistent development rates between processor speeds and memory bandwidth. In this study, we mitigate memory bottlenecks and reduce the computational load of…
We employ adaptive mesh refinement, implicit time stepping, a nonlinear multigrid solver and parallel computation, to solve a multi-scale, time dependent, three dimensional, nonlinear set of coupled partial differential equations for three…
Quantum embedding approaches involve the self-consistent optimization of a local fragment of a strongly correlated system, entangled with the wider environment. The `energy-weighted' density matrix embedding theory (EwDMET) was established…