In this work, a new formulation for central schemes based on staggered grids is proposed. It is based on a novel approach, in which first a time discretization is carried out, followed by the space discretization. The schemes obtained in this fashion have a simpler structure than previous central schemes. For high order schemes, this simplification results in higher computational efficiency. In...

In this work a new class of numerical methods for the BGK model of kinetic equations is presented. In principle, schemes of any order of accuracy in both space and time can be constructed with this technique. The methods proposed are based on an explicit-implicit time discretization. In particular the convective terms are treated explicitly, while the source terms are implicit. In this fashion ...

Inertial measurement units (IMUs) are small wearable sensors that have tremendous potential to be applied to clinical gait analysis. They allow objective evaluation of gait and movement disorders outside the clinic and research laboratory, and permit evaluation on large numbers of steps. However, repeatability and validity data of these systems are sparse for gait metrics. The purpose of this s...

In combination with simultaneous approximation terms, summation-by-parts (SBP) operators provide a flexible and efficient methodology that leads to consistent, conservative, and provably stable high-order discretizations. Traditional diagonal-norm SBP operators with a repeating interior point operator lead to solutions that have a global order of accuracy lower than the order of the interior po...

Gasdynamic features of slowly rotating gas axisymmetric accretion onto a gravitating center are investigated. The process of flow restructuring is studied as the angular velocity of accreting matter approaches the Keplerian angular velocity. For spherically-symmetric accretion, the conditions are found of the existence of a steady-state supersonic solution for various external boundary conditio...

In this paper, we give sixth order compact finite difference formula for second order integro-differential equations (IDE) with different boundary conditions, and both of error estimates and numerical experiments confirm our compact finite difference method can get fifth order of accuracy. We also adjust compact finite difference method for first order IDE and a system of IDE and give numerical...

Abstract: The use of the conventional advection diffusion equation in many physical situations has been questioned by many investigators in recent years and alternative diffusion models have been proposed. Fractional space derivatives are used to model anomalous diffusion or dispersion, where a particle plume spreads at a rate inconsistent with the classical Brownian motion model. When a fracti...

Knowledge of a truncated Fourier series expansion for a 2nperiodic function of finite regularity, which is assumed to be piecewise smooth in each period, is used to accurately reconstruct the corresponding function. An algebraic equation of degree M is constructed for the M singularity locations in each period for the function in question. The M coefficients in this algebraic equation are obtai...

In this paper, we propose a functional fitting s-stage Runge-Kutta method which is based on the exact integration of the set of the linearly independent functions u i (t), (i = 1,. .. , s). The method is exact when the solution of the ODE can be expressed as the linear combination of u i (t), although the method has an error for general ODE. In this work we investigate the order of accuracy of ...

We develop high order accurate source discretizations for hyperbolic wave propagation problems in first order formulation that are discretized by finite di↵erence schemes. By studying the Fourier series expansions of the source discretization and the finite di↵erence operator, we derive su cient conditions for achieving design accuracy in the numerical solution. Only half of the conditions in F...