Related papers: Taylor States in Stellarators: A Fast High-order B…
In this paper, we propose a trigonometric-interpolation approach for solutions of second order nonlinear ODEs with mixed boundary conditions. The method interpolates secondary derivative $y''$ of a target solution $y$ by a trigonometric…
This work develops a numerical solver based on the combination of isogeometric analysis (IGA) and the tensor train (TT) decomposition for the approximation of partial differential equations (PDEs) on parameter-dependent geometries. First,…
We present a higher-order boundary condition for atomistic simulations of dislocations that address the slow convergence of standard supercell methods. The method is based on a multipole expansion of the equilibrium displacement, combining…
The present work is devoted to the problem of boundary stabilization of the semilinear 1-D heat equation with nonlocal boundary conditions. The stabilizing controller is finite-dimensional, linear, given in an explicit form, involving only…
We present two (a decoupled and a coupled) integral-equation-based methods for the Morse-Ingard equations subject to Neumann boundary conditions on the exterior domain. Both methods are based on second-kind integral equation (SKIE)…
The state of art of time domain integral equation (TDIE) solvers has grown by leaps and bounds over the past decade. During this time, advances have been made in (i) the development of accelerators that can be retrofitted with these solvers…
We consider quadrature formulas of high order in time based on Radau-type, L-stable implicit Runge-Kutta schemes to solve time dependent stiff PDEs. Instead of solving a large nonlinear system of equations, we develop a method that performs…
A direct construction of equilibrium magnetic fields with toroidal topology at arbitrary order in the distance from the magnetic axis is carried out, yielding an analytical framework able to explore the landscape of possible magnetic flux…
In this paper, we develop a regularized higher-order Taylor based method for solving composite (e.g., nonlinear least-squares) problems. At each iteration, we replace each smooth component of the objective function by a higher-order Taylor…
We consider stationary axisymmetric solutions of the Euler-Poisson equations, which govern the internal structure of barotropic gaseous stars. We take the general form of the equation of states which cover polytropic gaseous stars indexed…
The model under consideration is based on approximate analytical solution of two dimensional stationary Navier-Stokes and Fourier-Kirchhoff equations. Approximations are based on the typical for natural convection assumptions: the fluid…
Embedded, or immersed, approaches have the goal of reducing to the minimum the computational costs associated with the generation of body-fitted meshes by only employing fixed, possibly Cartesian, meshes over which complex boundaries can…
Here we present the mathematical and numerical details of a general hybrid model for plasmas. All grid quantities are stored at cell centers on the grid. The most common discretization of the fields in PIC solvers is to have the electric…
For 2D topological insulators with strong electron-hole hybridization, such as HgTe/CdTe quantum wells, the widely used 4 x 4 k.p Hamiltonian based on the first electron and heavy hole sub-bands yields an equal number of physical and…
We apply topological methods to better understand how the magnetic field in the stellarator edge can be diverted away from the confined region. Our primary method is calculating the winding numbers of closed contours, which gives…
We introduce a high-order spline geometric approach for the initial boundary value problem for Maxwell's equations. The method is geometric in the sense that it discretizes in structure preserving fashion the two de Rham sequences of…
We apply results of Malliavin-Thalmaier-Watanabe for strong and weak Taylor expansions of solutions of perturbed stochastic differential equations (SDEs). In particular, we work out weight expressions for the Taylor coefficients of the…
We apply the near-axis expansion method for quasisymmetric magnetic fields with anisotropic pressure (developed in the companion paper, Part I [E. Rodriguez, A. Bhattacharjee, arXiv:2008.04715 (2020)]) to construct numerical solutions to…
In this paper, we introduce an algorithm that provides approximate solutions to semi-linear ordinary differential equations with highly oscillatory solutions, which, after an appropriate change of variables, can be rewritten as…
In several studies it has been observed that, when using stabilised $\mathbb{P}_k^{}\times\mathbb{P}_k^{}$ elements for both velocity and pressure, the error for the pressure is smaller, or even of a higher order in some cases, than the one…