Related papers: The Parareal Algorithm Applied to the FESOM 2 Ocea…
We present a numerical method which is able to approximate traveling waves (e.g. viscous profiles) in systems with hyperbolic and parabolic parts by a direct long-time forward simulation. A difficulty with long-time simulations of traveling…
The ParaOpt algorithm was recently introduced as a time-parallel solver for optimal-control problems with a terminal-cost objective, and convergence results have been presented for the linear diffusive case with implicit-Euler time…
In the paper, we present a strategy for accelerating posterior inference for unknown inputs in time fractional diffusion models. In many inference problems, the posterior may be concentrated in a small portion of the entire prior support.…
We introduce two multiscale numerical schemes for the time integration of weakly nonlinear Schr\"odinger equations, built upon the discretization of Picard iterates of the solution. These high-order schemes are designed to achieve high…
Numerical simulations of the air in the atmosphere and water in the oceans are essential for numerical weather prediction. The state-of-the-art for performing these fluid simulations relies on an Eulerian viewpoint, in which the fluid…
We construct a pseudospectral method for the solution of time-dependent, non-linear partial differential equations on a three-dimensional spherical shell. The problem we address is the treatment of tensor fields on the sphere. As a test…
Numerical simulation of wave propagation in an infinite medium is made possible by surrounding a finite region by a perfectly matched layer (PML). Using this approach a generalized three-dimensional (3D) formulation is proposed for…
We present a computational algorithm for computing short range forces between particles. The algorithm has two distinguishing features. First, it is optimized for multi-processor computers, and will use as many processors as are available.…
This paper deals with a time-split explicit/implicit approach for solving a two-dimensional hydrodynamic flow model with appropriate initial and boundary conditions. The time-split technique is employed to upwind the convection term and to…
We adopt the frozen Gaussian approximation (FGA) for modeling seismic waves. The method belongs to the category of ray-based beam methods. It decomposes seismic wavefield into a set of Gaussian functions and propagates these Gaussian…
The paper investigates a non-intrusive parallel time integration with multigrid for space-fractional diffusion equations in two spatial dimensions. We firstly obtain a fully discrete scheme via using the linear finite element method to…
Large molecular dynamics simulations (millions of atoms, tens of microseconds, thousands of processors) hit the strong scalability wall: simulation on twice as many processors does not take half the time. Inspired by large N-body space…
We propose a class of temporally high-order parametric finite element methods for simulating solid-state dewetting of thin films in two dimensions using a sharp-interface model. The process is governed by surface diffusion and contact point…
A novel approach to simulate dynamic, recurrent flows with very large time steps is presented. Data-driven forecasts based on the method of analogues are made employing a set of precomputed time series segments of flow fields. These…
The state of the art for physical hazard prediction from weather and climate requires expensive km-scale numerical simulations driven by coarser resolution global inputs. Here, a generative diffusion architecture is explored for downscaling…
We discuss the scalable parallel solution of the Poisson equation within a Particle-In-Cell (PIC) code for the simulation of electron beams in particle accelerators of irregular shape. The problem is discretized by Finite Differences.…
We have developed an efficacious algorithm for simulation of the beam-beam interaction in synchrotron colliders based on the nonlinear $\delta f$ method, where $\delta f$ is the much smaller deviation of the beam distribution from the…
In this work, a coarse-graining method previously proposed by the authors in a companion paper based on solving diffusion equations is applied to CFD-DEM simulations, where coarse graining is used to obtain solid volume fraction, particle…
Networks of interconnected resistors, springs and beams, or pores are standard models of studying scalar and vector transport processes in heterogeneous materials and media, such as fluid flow in porous media, and conduction, deformations,…
Spectral methods for solving partial differential equations (PDEs) and stochastic partial differential equations (SPDEs) often use Fourier or polynomial spectral expansions on either uniform and non-uniform grids. However, while very widely…