Related papers: Preconditioned fully implicit PDE solvers for monu…
In this paper, we consider a nonlinear PDE system governed by a parabolic heat equation coupled in a nonlinear way with a hyperbolic momentum equation describing the behavior of a displacement field coupled with a nonlinear elliptic…
In this work, we are interested in solving large linear systems stemming from the Extra-Membrane-Intra (EMI) model, which is employed for simulating excitable tissues at a cellular scale. After setting the related systems of partial…
A numerical scheme is developed for systems of conservation laws on manifolds which arise in high speed aerodynamics and magneto-aerodynamics. The systems are presented in an arbitrary coordinate system on the manifold and involve source…
The porous medium equation (PME) is a typical nonlinear degenerate parabolic equation. We have studied numerical methods for PME by an energetic variational approach in [C. Duan et al, J. Comput. Phys., 385 (2019) 13-32], where the…
This paper presents a general theory and isogeometric finite element implementation for studying mass conserving phase transitions on deforming surfaces. The mathematical problem is governed by two coupled fourth-order nonlinear partial…
We propose a predictor-corrector adaptive method for the study of hyperbolic partial differential equations (PDEs) under uncertainty. Constructed around the framework of stochastic finite volume (SFV) methods, our approach circumvents…
Parameterization (closure) schemes in numerical weather and climate prediction models account for the effects of physical processes that cannot be resolved explicitly by these models. Methods for finding physical parameterization schemes…
A preconditioned, multipole-accelerated, Krylov-subspace iterative algorithm for the electromagnetic scattering analysis of three dimensional (3D), arbitrary shaped dielectric structures composed of single and multi-layered dielectric…
In this paper, we consider effective discretization strategies and iterative solvers for nonlinear PDE-constrained optimization models for pattern evolution within biological processes. Upon a Sequential Quadratic Programming linearization…
In this work, a systematic protocol is proposed to automatically parametrize implicit solvent models with polar and nonpolar components. The proposed protocol utilizes the classical Poisson model or the Kohn-Sham density functional theory…
This article presents a method for solving large-scale linear inverse problems regular- ized with a nonlinear, edge-preserving penalty term such as the total variation or Perona-Malik. In the proposed scheme, the nonlinearity is handled…
Simulations of the dynamics generated by partial differential equations (PDEs) provide approximate, numerical solutions to initial value problems. Such simulations are ubiquitous in scientific computing, but the correctness of the results…
We present a novel implicit numerical implementation of the parabolic-hyperbolic formulation of the constraints of general relativity. The proposed method is unconditionally stable, has the advantage of not requiring the imposition of any…
A semi-implicit two-phase double-point Material Point Method (MPM) formulation, based on the incremental fractional-step method to model large deformation geotechnical problems has been derived. The semi-implicit formulation has two…
Bayesian statistical inverse problems are often solved with Markov chain Monte Carlo (MCMC)-type schemes. When the problems are governed by large-scale discrete nonlinear partial differential equations (PDEs), they are computationally…
Many applications involving porous media--notably reservoir engineering and geologic applications--involve tight coupling between multiphase fluid flow, transport, and poromechanical deformation. While numerical models for these processes…
Neural PDE solvers offer a powerful tool for modeling complex dynamical systems, but often struggle with error accumulation over long time horizons and maintaining stability and physical consistency. We introduce a multiscale implicit…
Recently, a family of models that couple multifluid systems to the full Maxwell equations draw a lot of attention in laboratory, space, and astrophysical plasma modeling. These models are more complete descriptions of the plasma than…
The solutions to a large class of semi-linear parabolic PDEs are given in terms of expectations of suitable functionals of a tree of branching particles. A sufficient, and in some cases necessary, condition is given for the integrability of…
We consider an efficient preconditioner for boundary integral equation (BIE) formulations of the two-dimensional Stokes equations in porous media. While BIEs are well-suited for resolving the complex porous geometry, they lead to a dense…