Related papers: Large-scale earthquake sequence simulations on 3D …
This study presents a generalized multiscale multimesh finite element method ($\text{M}^2$-FEM) that addresses several long-standing challenges in the numerical simulation of integral structural theories, often used to model multiscale and…
The trend towards highly parallel multi-processing is ubiquitous in all modern computer architectures, ranging from handheld devices to large-scale HPC systems; yet many applications are struggling to fully utilise the multiple levels of…
Computational chemistry allows researchers to experiment in sillico: by running a computer simulations of a biological or chemical processes of interest. Molecular dynamics with molecular mechanics model of interactions simulates N-body…
In this paper, a highly efficient fast boundary element method (BEM) for solving large-scale engineering acoustic problems in a broad frequency range is developed and implemented. The acoustic problems are modeled by the Burton-Miller…
The interpolated bounce-back scheme and the immersed boundary method are the two most popular algorithms in treating a no-slip boundary on curved surfaces in the lattice Boltzmann method. While those algorithms are frequently implemented in…
Cascaded or central-moment-based lattice Boltzmann method (CLBM) proposed in [Geier \textit{et al.}, Phys. Rev. E \textbf{63}, 066705 (2006)] possesses very good numerical stability. However, two constraints exist in three-dimensional (3D)…
Natural earthquake fault systems are highly non-homogeneous. The inhomogeneities occur be- cause the earth is made of a variety of materials which hold and dissipate stress differently. In this work, we study scaling in earthquake fault…
A major goal in earthquake physics is to derive a constitutive framework for fault slip that captures the dependence of shear strength on fault rheology, sliding velocity, and pore-fluid pressure. In this study, we present H-MEC…
This paper presents an octree construction method, called Cornerstone, that facilitates global domain decomposition and interactions between particles in mesh-free numerical simulations. Our method is based on algorithms developed for 3D…
Calculating portions of eigenvalues and eigenvectors of matrices or matrix pencils has many applications. An approach to this calculation for Hermitian problems based on a density matrix has been proposed in 2009 and a software package…
The quantum circuit model is the de-facto way of designing quantum algorithms. Yet any level of abstraction away from the underlying hardware incurs overhead. In the era of near-term, noisy, intermediate-scale quantum (NISQ) hardware with…
Despite decades of research, creating accurate, robust, and efficient lattice Boltzmann methods (LBM) on non-uniform grids with seamless GPU acceleration remains challenging. This work introduces a novel strategy to address this challenge…
Intensity estimation is a common problem in statistical analysis of spatial point pattern data. This paper proposes a nonparametric Bayesian method for estimating the spatial point process intensity based on mixture of finite mixture (MFM)…
This work outlines a Lattice Boltzmann Method (LBM) for geometrically and constitutively nonlinear solid mechanics to simulate large deformations under dynamic loading conditions. The method utilizes the moment chain approach, where the…
Fracture coalescence is a critical phenomenon for creating large fractures from smaller flaws, affecting fracture network flow and seismic energy release potential. In this paper, simulations of fracture coalescence processes in granite…
The immersed isogeometric Boundary Element Method is presented and applied to the simulation of underground excavations. Nonuniform rational B-splines (NURBS) are used for the accurate definition of complex geometries with few parameters.…
This article presents the development and validation of a hybrid multi-sphere discrete element framework - Rigid3D, for the simulation of granular systems with arbitrarily shaped particles in 3D space. In this DEM framework, a non-spherical…
Strongly coupled immersed boundary (IB) methods solve the nonlinear fluid and structural equations of motion simultaneously for strongly enforcing the no-slip constraint on the body. Handling this constraint requires solving several large…
Several researchers have developed a rich toolbox of matrix compression techniques that exploit structure and redundancy in large matrices. Classical methods such as the block low-rank format and the Fast Multipole Method make it possible…
Combining sum factorization, weighted quadrature, and row-based assembly enables efficient higher-order computations for tensor product splines. We aim to transfer these concepts to immersed boundary methods, which perform simulations on a…