The numerical solutions to partial differential equations using fractals.

Date

2022

Authors

Karve, Prachi

Journal Title

Journal ISSN

Volume Title

Publisher

Abstract

The goal of this work is to construct and use continuous fractal functions to find a numerical solution of partial differential equations in one dimension. We focus on univariate, R-valued, continuous fractal functions defined on interval [0,1]. These fractal functions provide more flexibility for approximation and interpolation. They are used as interpolation functions in the concept of hierarchical basis in the finite element method. Based on a continuous Read-Bajraktarevic operator defined on a space of bounded continuous functions, a fractal bubble is initially constructed on the interval [0,1]. Further, linear nodal basis functions defined on the interval [0,1] are expressed as fixed points of an iterated function system. Consequently, they are established as continuous fractal functions. Thus, a hierarchical fractal basis is constructed on the interval [0,1], consisting of three shape functions: two linear nodal basis functions and a hierarchical element namely, the fractal bubble. The nodal basis functions ensure the continuity at connecting nodes while the fractal bubble yields an additional degree of freedom without increasing the numbers of nodes in a mesh. Subsequently, a finite element mesh is constructed by discretizing the interval [0,1] in 2^n elements. The hierarchical fractal basis on a finite element mesh is then obtained by translating and dilating the three hierarchical basis functions. This basis defined on the finite element mesh spans the space of fractal functions and is used as a finite element space. This space has a direct sum hierarchical decomposition into a space of piecewise linear nodal basis functions and a space of piecewise fractal functions which are zero at the nodes. Subsequently, to set up element matrices in a finite element model, the fractal functions are integrated explicitly using the self-similarity property. However, a trapezoidal rule is used for numerical integration of fractal functions where the explicit integration of fractal functions is unattainable. As an illustration, a numerical solution to the Poisson equation with homogeneous Dirichlet boundary conditions is found. A numerical error is computed in L^2 norm. Bound of an interpolation error in L^2 norm and H^1 seminorm is established theoretically. A python code is used for generating the fractal bubble from an iterated function system and for the computational purpose to find an approximate solution to the Poisson equation in one dimension.

Description

Keywords

Citation

Source

Type

Thesis (MPhil)

Book Title

Entity type

Access Statement

License Rights

Restricted until

Downloads

File
Description