Related papers: Higher order operator splitting Fourier spectral m…
In this work, we explore the use of operator splitting algorithms for solving regularized structural topology optimization problems. The context is the classical structural design problems (e.g., compliance minimization and compliant…
We present high order explicit geometric integrators to solve linear-quadratic optimal control problems and $N$-player differential games. These problems are described by a system coupled non-linear differential equations with boundary…
It is well known that second order linear ordinary differential equations with slowly varying coefficients admit slowly varying phase functions. This observation is the basis of the Liouville-Green method and many other techniques for the…
We study Cauchy problems of fractional differential equations in both space and time variables by expressing the solution in terms of ``stochastic composition" of the solutions to two simpler problems. These Cauchy sub-problems respectively…
A new analytical operator method is discussed which solves linear ordinary differential equations with regular singularities. Solutions are obtained in analytic series form and also in Mellin-Barnes-type contour integral form. Exact series…
We demonstrate application of the method of higher-order operators to nonlinear standard optomechanics. It is shown that a symmetry breaking in frequency shifts exists, corresponding to inequivalency of red and blue side-bands. This arises…
In this work, we develop an $\mathcal{O}(N)$ implicit real space method in 1D and 2D for the Cahn Hilliard (CH) and vector Cahn Hilliard (VCH) equations, based on the Method Of Lines Transpose (MOL$^\text{T}$) formulation. This formulation…
We propose a discontinuous Galerkin method for convection-subdiffusion equations with a fractional operator of order $\alpha (1<\alpha<2)$ defined through the fractional Laplacian. The fractional operator of order $\alpha$ is expressed as a…
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…
Motivated by Fredholm theory, we develop a framework to establish the convergence of spectral methods for operator equations $\mathcal L u = f$. The framework posits the existence of a left-Fredholm regulator for $\mathcal L$ and the…
In this paper, we consider a second-order scalar auxiliary variable (SAV) Fourier spectral method to solve the nonlinear fractional generalized wave equation. Unconditional energy conservation or dissipation properties of the fully discrete…
It is common in classical mechanics to encounter systems whose Hamiltonian $H$ is the sum of an often exactly integrable Hamiltonian $H_0$ and a small perturbation $\epsilon H_1$ with $\epsilon\ll1$. Such near-integrability can be exploited…
While recent advances in deep learning have shown promising efficiency gains in solving time-dependent partial differential equations (PDEs), matching the accuracy of conventional numerical solvers still remains a challenge. One strategy to…
The conventional second-order Path Integral Monte Carlo method is plagued with the sign problem in solving many-fermion systems. This is due to the large number of anti-symmetric free fermion propagators that are needed to extract the…
In this paper we consider a class of boundary value problems for third order nonlinear functional differential equation. By the reduction of the problem to operator equation we establish the existence and uniqueness of solution and…
Our paper "Solving Third Order Linear Difference Equations in Terms of Second Order Equations" gave two algorithms for solving difference equations in terms of lower order equations: an algorithm for absolute factorization, and an algorithm…
Hamiltonian systems are known to conserve the Hamiltonian function, which describes the energy evolution over time. Obtaining a numerical spatio-temporal scheme that accurately preserves the discretized Hamiltonian function is often a…
This paper proposes a frequency/time hybrid integral-equation method for the time dependent wave equation in two and three-dimensional spatial domains. Relying on Fourier Transformation in time, the method utilizes a fixed…
We discuss systematic extensions of the standard (St{\"o}rmer-Verlet) splitting method for differential equations of Hamiltonian mechanics, with relative accuracy of order $\tau^2$ for a timestep of length $\tau$, to higher orders in…
Variational quantum algorithms are a promising tool for solving partial differential equations. The standard approach for its numerical solution are finite difference schemes, which can be reduced to the linear algebra problem. We consider…