Related papers: A Jacobian-free Newton-Krylov method for high-orde…
The focus in this work is on interior-point methods for inequality-constrained quadratic programs, and particularly on the system of nonlinear equations to be solved for each value of the barrier parameter. Newton iterations give high…
We introduce a quadratically convergent semismooth Newton method for nonlinear semidefinite programming that eliminates the need for the generalized Jacobian regularity, a common yet stringent requirement in existing approaches. Our…
In this paper, we develop zeroth-order algorithms with provably (nearly) optimal sample complexity for stochastic bilevel optimization, where only noisy function evaluations are available. We propose two distinct algorithms: the first is…
Direct methods to obtain global stability modes are restricted by the daunting sizes and complexity of Jacobians encountered in general three-dimensional flows. Jacobian-free iterative approaches such as Arnoldi methods have greatly…
For the simulations of unsteady flow, the global time step becomes really small with a large variation of local cell size. In this paper, an implicit high-order gas-kinetic scheme (HGKS) is developed to remove the restrictions on the time…
In this short exposition, we describe equilibria and periodic orbits in terms of the flow map, {\Phi}, and discuss the essentials of the Jacobian-free Newton-Krylov (JFNK) method that can be used to find them. This method requires little…
A novel numerical approach to solving the shallow-water equations on the sphere using high-order numerical discretizations in both space and time is proposed. A space-time tensor formalism is used to express the equations of motion…
We discuss the order, efficiency, stability and positivity of several meshless schemes for linear scalar hyperbolic equations. Meshless schemes are Generalised Finite Difference Methods (GFDMs) for arbitrary irregular grids in which there…
Efficient and energy stable high order time marching schemes are very important but not easy to construct for the study of nonlinear phase dynamics. In this paper, we propose and study two linearly stabilized second order semi-implicit…
We present a class of hybrid FD-FV (finite difference and finite volume) methods for solving general hyperbolic conservation laws written in first-order form. The presentation focuses on one- and two-dimensional Cartesian grids; however,…
The modeling of electric machines and power transformers typically involves systems of nonlinear magnetostatics or -quasistatics, and their efficient and accurate simulation is required for the reliable design, control, and optimization of…
The magnetohydrodynamics (MHD) equations are continuum models used in the study of a wide range of plasma physics systems, including the evolution of complex plasma dynamics in tokamak disruptions. However, efficient numerical solution…
We present a new approach that allows large time steps in dynamic simulations. Our approach, ConJac, is based on condensation, a technique for eliminating many degrees of freedom (DOFs) by expressing them in terms of the remaining degrees…
Based on the Jacobi polynomial expansion, an arbitrary high-order Discontinuous Galerkin solver for compressible flows on unstructured meshes is proposed in the present work. First, we construct orthogonal polynomials for 2D and 3D…
An effective numerical method is presented for optimizing model parameters that can be applied to any type of system of non-linear equations and any number of data-points, which does not require explicit formulation of the objective…
Recent efforts in applying implicit networks to solve inverse problems in imaging have achieved competitive or even superior results when compared to feedforward networks. These implicit networks only require constant memory during…
We propose a class of weighted compact central (WCC) schemes for solving hyperbolic conservation laws. The linear version can be considered as a high-order extension of the central Lax-Friedrichs (LxF) scheme and the central conservation…
We present a comparison of different multigrid approaches for the solution of systems arising from high-order continuous finite element discretizations of elliptic partial differential equations on complex geometries. We consider the…
In this paper we propose the first better than second order accurate method in space and time for the numerical solution of the resistive relativistic magnetohydrodynamics (RRMHD) equations on unstructured meshes in multiple space…
Logarithmic conformation reformulations for viscoelastic constitutive laws have alleviated the high Weissenberg number problem, and the exploration of highly elastic flows became possible. However, stabilized formulations for logarithmic…