Related papers: Hermite Methods for the Scalar Wave Equation
The gradient expansion of the kinetic energy functional, when applied for atoms or finite systems, usually grossly overestimates the energy in the fourth order and generally diverges in the sixth order. We avoid the divergence of the…
We propose an efficient and easy-to-implement gradient-enhanced least squares Monte Carlo method for computing price and Greeks (i.e., derivatives of the price function) of high-dimensional American options. It employs the sparse Hermite…
The discrete-dipole approximation (DDA) is a flexible technique for computing scattering and absorption by targets of arbitrary geometry. In this paper we perform systematic study of various non-stationary iterative (conjugate gradient)…
A particle scheme for scalar conservation laws in one space dimension is presented. Particles representing the solution are moved according to their characteristic velocities. Particle interaction is resolved locally, satisfying exact…
We introduce an efficient numerical method for second order linear ODEs whose solution may vary between highly oscillatory and slowly changing over the solution interval. In oscillatory regions the solution is generated via a nonoscillatory…
In this introductory work I will present the Finite Difference method for hyperbolic equations, focusing on a method which has second order precision both in time and space (the so-called staggered leapfrog method) and applying it to the…
We analyze the asymptotic behavior of solutions to wave equations with strong damping terms. If the initial data belong to suitable weighted $L^1$ spaces, lower bounds for the difference between the solutions and the leading terms in the…
We study the asymptotic behavior of solutions for the semilinear damped wave equation with variable coefficients. We prove that if the damping is effective, and the nonlinearity and other lower order terms can be regarded as perturbations,…
Standard discontinuous Galerkin methods, based on piecewise polynomials of degree $ \qq=0,1$, are considered for temporal semi-discretization for second order hyperbolic equations. The main goal of this paper is to present a simple and…
A fast and accurate algorithm for the computation of Gauss-Hermite and generalized Gauss-Hermite quadrature nodes and weights is presented. The algorithm is based on Newton's method with carefully selected initial guesses for the nodes and…
We consider the Navier-Stokes equations in a channel with a narrowing of varying height. The model is discretized with high-order spectral element ansatz functions, resulting in 6372 degrees of freedom. The steady-state snapshot solutions…
Elastic waves of short wavelength propagating through the upper layer of the Earth appear to move faster at large separations of source and receiver than at short separations. This scale dependent velocity is a manifestation of Fermat's…
We consider time-harmonic electromagnetic wave equations in composites of a dispersive material surrounded by a classical material. In certain frequency ranges this leads to sign-changing permittivity and/or permeability. Previously meshing…
Based on a new approximation method, namely pseudospectral method, a solution for the three order nonlinear ordinary differential laminar boundary layer Falkner-Skan equation has been obtained on the semi-infinite domain. The proposed…
Substantially extending previous results of the authors for smooth solutions in the viscous case, we develop linear damping estimates for periodic roll-wave solutions of the inviscid Saint-Venant equations and related systems of hyperbolic…
This paper describes an efficient algorithm for computing steady two-dimensional surface gravity wave in irrotational motion. The algorithm complexity is O(N log N), N being the number of Fourier modes. The algorithm allows the arbitrary…
In this paper, we consider numerical homogenization of acoustic wave equations with heterogeneous coefficients, namely, when the bulk modulus and the density of the medium are only bounded. We show that under a Cordes type condition the…
In this paper we propose stochastic gradient-free methods and accelerated methods with momentum for solving stochastic optimization problems. All these methods rely on stochastic directions rather than stochastic gradients. We analyze the…
Ordinary Differential Equations are a simple but powerful framework for modeling complex systems. Parameter estimation from times series can be done by Nonlinear Least Squares (or other classical approaches), but this can give…
Quadratization refers to a transformation of an arbitrary system of polynomial ordinary differential equations to a system with at most quadratic right-hand side. Such a transformation unveils new variables and model structures that…