Related papers: A Jacobian-free Newton-Krylov method for high-orde…
This paper addresses the study of a new class of nonsmooth optimization problems, where the objective is represented as a difference of two generally nonconvex functions. We propose and develop a novel Newton-type algorithm to solving such…
We study finite-sum non-convex optimization $\min_{x\in\mathbb{R}^d} F(x) \;=\; \frac{1}{n}\sum_{i=1}^n f_i(x)$ and analyze a variance-reduced cubic Newton method based on EMA-smoothed SARAH estimators for both gradient and Hessian…
This work presents a novel stabilization strategy for the Galerkin formulation of the incompressible Navier-Stokes equations, developed to achieve high accuracy while ensuring convergence and compatibility with high-order elements on…
We present a high-order spatial discretization of a continuum gyrokinetic Vlasov model in axisymmetric tokamak edge plasma geometries. Such models describe the phase space advection of plasma species distribution functions in the absence of…
In this paper, we introduce a new finite element method for solving the Stokes equations in the primary velocity-pressure formulation. This method employs $H(div)$ finite elements to approximate velocity, which leads to two unique…
We propose a fixed-point-based numerical framework for computing stationary states of nonlocal Fokker-Planck-type equations. Instead of discretising the differential operators directly, we reformulate the stationary problem as a nonlinear…
A modification of Newton's method for solving systems of $n$ nonlinear equations is presented. The new matrix-free method relies on a given decomposition of the invertible Jacobian of the residual into invertible sparse local Jacobians…
Recently, a new stabilizer free weak Galerkin method (SFWG) is proposed, which is easier to implement and more efficient. The main idea is that by letting $j\geq j_{0}$ for some $j_{0}$, where $j$ is the degree of the polynomials used to…
A robust finite volume method for viscoelastic flow analysis on general unstructured meshes is developed. It is built upon a general-purpose stabilization framework for high Weissenberg number flows. The numerical framework provides full…
This paper proposes a strong second-order two-step explicit/implicit technique with spectral orthogonal basis Galerkin finite element method for solving a two-dimensional Gray-Scott model subject to appropriate initial and boundary…
In this paper, we propose a high-order extension of the multiscale method introduced by the authors in [SIAM J. Numer. Anal., 63(4) (2025), pp. 1617--1641] for heterogeneous Stokes problems, while also providing several other improvements,…
This paper develops a robust three-level time split high-order Leapfrog/Crank-Nicolson technique for solving the two-dimensional unsteady sobolev and regularized long wave equations arising in fluid mechanics. A deep analysis of the…
The Volume-Averaged Navier-Stokes equations are used to study fluid flow in the presence of fixed or moving solids such as packed or fluidized beds. We develop a high-order finite element solver using both forms A and B of these equations.…
In this work, a simple fourth-order accurate finite volume semi-discrete scheme is introduced to solve astrophysical magnetohydrodynamics (MHD) problems on Cartesian meshes. Hydrodynamic quantities like density, momentum and energy are…
A sequential quadratic optimization algorithm is proposed for solving smooth nonlinear equality constrained optimization problems in which the objective function is defined by an expectation of a stochastic function. The algorithmic…
In this paper, we introduce and analyze a class of numerical schemes that demonstrate remarkable superiority in terms of efficiency, the preservation of positivity, energy stability, and high-order precision to solve the time-dependent…
A high-order numerical method is developed for solving the Cahn-Hilliard-Navier-Stokes equations with the Flory-Huggins potential. The scheme is based on the $Q_k$ finite element with mass lumping on rectangular grids, the second-order…
This article presents a high order conservative flux optimization (CFO) finite element method for the elliptic diffusion equations. The numerical scheme is based on the classical Galerkin finite element method enhanced by a flux…
Classical theory for quasi-Newton schemes has focused on smooth deterministic unconstrained optimization while recent forays into stochastic convex optimization have largely resided in smooth, unconstrained, and strongly convex regimes.…
This paper is concerned with numerical analysis of two fully discrete Chorin-type projection methods for the stochastic Stokes equations with general non-solenoidal multiplicative noise. The first scheme is the standard Chorin scheme and…