889 research outputs found
Optimal Collocation Nodes for Fractional Derivative Operators
Spectral discretizations of fractional derivative operators are examined,
where the approximation basis is related to the set of Jacobi polynomials. The
pseudo-spectral method is implemented by assuming that the grid, used to
represent the function to be differentiated, may not be coincident with the
collocation grid. The new option opens the way to the analysis of alternative
techniques and the search of optimal distributions of collocation nodes, based
on the operator to be approximated. Once the initial representation grid has
been chosen, indications on how to recover the collocation grid are provided,
with the aim of enlarging the dimension of the approximation space. As a
results of this process, performances are improved. Applications to fractional
type advection-diffusion equations, and comparisons in terms of accuracy and
efficiency are made. As shown in the analysis, special choices of the nodes can
also suggest tricks to speed up computations
Diffusive approximation of a time-fractional Burger's equation in nonlinear acoustics
A fractional time derivative is introduced into the Burger's equation to
model losses of nonlinear waves. This term amounts to a time convolution
product, which greatly penalizes the numerical modeling. A diffusive
representation of the fractional derivative is adopted here, replacing this
nonlocal operator by a continuum of memory variables that satisfy local-in-time
ordinary differential equations. Then a quadrature formula yields a system of
local partial differential equations, well-suited to numerical integration. The
determination of the quadrature coefficients is crucial to ensure both the
well-posedness of the system and the computational efficiency of the diffusive
approximation. For this purpose, optimization with constraint is shown to be a
very efficient strategy. Strang splitting is used to solve successively the
hyperbolic part by a shock-capturing scheme, and the diffusive part exactly.
Numerical experiments are proposed to assess the efficiency of the numerical
modeling, and to illustrate the effect of the fractional attenuation on the
wave propagation.Comment: submitted to Siam SIA
Regularity theory and high order numerical methods for the (1D)-fractional Laplacian
This paper presents regularity results and associated high-order numerical methods for one-dimensional Fractional-Laplacian boundary-value problems. On the basis of a factorization of solutions as a product of a certain edge-singular weight times a ``regular´´ unknown, a characterization of the regularity of solutions is obtained in terms of the smoothness of the corresponding right-hand sides. In particular, for right-hand sides which are analytic in a Bernstein Ellipse, analyticity in the same Bernstein Ellipse is obtained for the ``regular´´ unknown. Moreover, a sharp Sobolev regularity result is presented which completely characterizes the co-domain of the Fractional-Laplacian operator in terms of certain weighted Sobolev spaces introduced in (Babu{s}ka and Guo, SIAM J. Numer. Anal. 2002). The present theoretical treatment relies on a full eigendecomposition for a certain weighted integral operator in terms of the Gegenbauer polynomial basis. The proposed Gegenbauer-based Nystr"om numerical method for the Fractional-Laplacian Dirichlet problem, further, is significantly more accurate and efficient than other algorithms considered previously. The sharp error estimates presented in this paper indicate that the proposed algorithm is spectrally accurate, with convergence rates that only depend on the smoothness of the right-hand side. In particular, convergence is exponentially fast (resp. faster than any power of the mesh-size) for analytic (resp. infinitely smooth) right-hand sides. The properties of the algorithm are illustrated with a variety of numerical results.Fil: Acosta, Gabriel. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Departamento de Matemática; ArgentinaFil: Borthagaray, Juan Pablo. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Departamento de Matemática; ArgentinaFil: Bruno, Oscar Ricardo. California Institute Of Technology; Estados UnidosFil: Maas, MartÃn Daniel. Consejo Nacional de Investigaciónes CientÃficas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de AstronomÃa y FÃsica del Espacio. - Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de AstronomÃa y FÃsica del Espacio; Argentin
- …