Related papers: Improved scalar auxiliary variable schemes for ori…
An accurate system to study the stability of pipe flow that ensures regularity is presented. The system produces a spectrum that is as accurate as Meseguer \& Trefethen (2000), while providing flexibility to amend the boundary conditions…
Parametric finite element discretizations of constrained geometric flows must simultaneously address high-order geometric stiffness, mesh degeneration, and nonlinear global constraints. This paper develops a stabilized dual-SAV (scalar…
We prove the stability with respect to the flux of solutions to initial-boundary value problems for scalar non-autonomous conservation laws in one space dimension. Key estimates are obtained through a careful construction of the solutions.
With the purpose of examining biased updates in variance-reduced stochastic gradient methods, we introduce SVAG, a SAG/SAGA-like method with adjustable bias. SVAG is analyzed in a cocoercive root-finding setting, a setting which yields the…
We develop a novel and efficient iterative scheme for solving incompressible steady Navier-Stokes equations. The method is an adaptation of the Incremental Viscosity Splitting approximation for unsteady flows to steady equations. At each…
This paper considers a modular grad-div stabilization method for approximating solutions of the time-dependent Boussinesq model of non-isothermal flows. The proposed method adds a minimally intrusive step to an existing Boussinesq code,…
We construct new higher-order implicit-explicit (IMEX) schemes using the generalized scalar auxiliary variable (GSAV) approach for the Landau-Lifshitz equation. These schemes are linear, length preserving and only require solving one…
This paper presents two modular grad-div algorithms for calculating solutions to the Navier-Stokes equations (NSE). These algorithms add to an NSE code a minimally intrusive module that implements grad-div stabilization. The algorithms do…
The main objective of this paper is to present an efficient structure-preserving scheme, which is based on the idea of the scalar auxiliary variable approach, for solving the space fractional nonlinear Schr\"{o}dinger equation. First, we…
We present a parametric family of semi-implicit second order accurate numerical methods for non-conservative and conservative advection equation for which the numerical solutions can be obtained in a fixed number of forward and backward…
Here we propose the Reweighted Autoencoded Variational Bayes for Enhanced Sampling (RAVE) method, a new iterative scheme that uses the deep learning framework of variational autoencoders to enhance sampling in molecular simulations. RAVE…
The long-term success of climate models that operate on multi-resolution grids will depend on access to subgrid parameterizations that act appropriately across a wide range of spatial and temporal scales. As the first step in a series of…
The increasing integration of distributed energy resources (DERs) calls for new monitoring and operational planning tools to ensure stability and sustainability in distribution grids. One idea is to use existing monitoring tools in…
We construct first- and second-order time discretization schemes for the Cahn-Hilliard-Navier-Stokes system based on the multiple scalar auxiliary variables approach (MSAV) approach for gradient systems and (rotational) pressure-correction…
We consider constrained optimization problems with a nonsmooth objective function in the form of mathematical expectation. The Sample Average Approximation (SAA) is used to estimate the objective function and variable sample size strategy…
A suitable notion of weak solution to infinite-dimensional rate-independent systems, called Inertial Balanced Viscosity (IBV) solution, is introduced. The key feature of such notion is that the energy dissipated at jump discontinuities…
In this paper, a self-adaptive contractive (SAC) algorithm is proposed for enhanced dynamic phasor estimation in the diverse operating conditions of modern power systems. At a high-level, the method is composed of three stages: parameter…
We focus on the numerical approximation of the Cahn-Hilliard type equations, and present a family of second-order unconditionally energy-stable schemes. By reformulating the equation into an equivalent system employing a scalar auxiliary…
This paper considers a type of incremental aggregated gradient (IAG) method for large-scale distributed optimization. The IAG method is well suited for the parameter server architecture as the latter can easily aggregate potentially staled…
Phase retrieval (PR) is an inverse problem about recovering a signal from phaseless linear measurements. This problem can be effectively solved by minimizing a nonconvex amplitude-based loss function. However, this loss function is…