Related papers: OBMeshfree: An optimization-based meshfree solver …
In this paper, we propose a meshfree method based on the Gaussian radial basis function (RBF) to solve both classical and fractional PDEs. The proposed method takes advantage of the analytical Laplacian of Gaussian functions so as to…
We present a Finite Volume formulation for determining discontinuous distributions of magnetic fields within non-orthogonal and non-uniform meshes. The numerical approach is based on the discretization of the vector potential variant of the…
Nonlocal diffusion model provides an appropriate description of the diffusion process of solute in the complex medium, which cannot be described properly by classical theory of PDE. However, the operators in the nonlocal diffusion models…
We develop a novel Hybrid High-Order method for the simulation of Darcy flows in fractured porous media. The discretization hinges on a mixed formulation in the bulk region and on a primal formulation inside the fracture. Salient features…
In this paper, we present numerical simulations with local and nonlocal models under dynamic loading conditions. We show that for finite element (FE) computations of high-velocity, impact problems with softening material models will result…
This work develops a multiscale solution decomposition (MSD) method for nonlocal-in-time problems to separate a series of known terms with multiscale singularity from the original singular solution such that the remaining unknown part…
This paper presents a framework for analyzing and designing cylindrical omega-bianisotropic metasurfaces, inspired by mode matching and digital signal processing techniques. Using the discrete Fourier transform, we decompose the the…
The discrete fracture model (DFM) has been widely used in the simulation of fluid flow in fractured porous media. Traditional DFM uses the so-called hybrid-dimensional approach to treat fractures explicitly as low-dimensional entries (e.g.…
A numerical framework is developed to solve various types of PDEs on complicated domains, including steady and time-dependent, non-linear and non-local PDEs, with different boundary conditions that can also include non-linear and non-local…
We consider a strongly heterogeneous medium saturated by an incompressible viscous fluid as it appears in geomechanical modeling. This poroelasticity problem suffers from rapidly oscillating material parameters, which calls for a thorough…
Peridynamic (PD) theory is significant and promising in engineering and materials science; however, it imposes challenges owing to the enormous computational cost caused by its nonlocality. Our main contribution, which overcomes the…
We propose a method for efficiently coupling the finite element method with atomistic simulations, while using molecular dynamics or kinetic Monte Carlo techniques. Our method can dynamically build an optimized unstructured mesh that…
In this work, we introduce a novel algorithm for the Biot problem based on a Hybrid High-Order discretization of the mechanics and a Symmetric Weighted Interior Penalty discretization of the flow. The method has several assets, including,…
This paper proposes a novel approach that combines variational integration with the bonded discrete element method (BDEM) to achieve faster and more accurate fracture simulations. The approach leverages the efficiency of implicit…
We present a method for generating orthogonal quadrilateral meshes subject to user-defined feature alignment and sizing constraints. The approach relies on computing integrable orthogonal frame fields, whose symmetries are implicitly…
We demonstrate that the recently developed Optimal Uncertainty Quantification (OUQ) theory, combined with recent software enabling fast global solutions of constrained non-convex optimization problems, provides a methodology for rigorous…
In this paper, we compute stiffness matrix of the nonlocal Laplacian discretized by the piecewise linear finite element on nonuniform meshes, and implement the FEM in the Fourier transformed domain. We derive useful integral expressions of…
We present a computational framework for efficient optimization-based "inverse design" of large-area "metasurfaces" (subwavelength-patterned surfaces) for applications such as multi-wavelength and multi-angle optimizations, and…
Optimization-based (OB) alternatives to traditional flux limiters couch preservation of properties such as local bounds and maximum principles into optimization problems, which impose these properties through inequality constraints. In this…
A method for model reduction in nonlinear ODE systems is demonstrated through computational examples. The method does not require an implicit separation of time-scales in the fine dynamics to be effective. From the computational standpoint,…