The local radial basis function (RBF) method is a promising solver for variable‐order time fractional diffusion equation (TFDE), as it overcomes the computational burden of the traditional global method. Application of the local RBF method is limited to Fickian diffusion, while real‐world diffusion is usually non‐Fickian in multiple dimensions. This article is the first to extend the application of the local RBF method to two‐dimensional, variable‐order, time fractional diffusion equation in complex shaped domains. One of the main advantages of the local RBF method is that only the nodes located in the subdomain, surrounding the local point, need to be considered when calculating the numerical solution at this point. This approach can perform well with large scale problems and can also mitigate otherwise ill‐conditioned problems. The proposed numerical approach is checked against two examples with curved boundaries and known analytical solutions. Shape parameter and subdomain node number are investigated for their influence on the accuracy of the local RBF solution. Furthermore, quantitative analysis, based on root‐mean‐square error, maximum absolute error, and maximum error of the partial derivative indicates that the local RBF method is accurate and effective in approximating the variable‐order TFDE in two‐dimensional irregular domains.
- Award ID(s):
- 2012291
- NSF-PAR ID:
- 10295130
- Date Published:
- Journal Name:
- Journal of Inverse and Ill-posed Problems
- Volume:
- 29
- Issue:
- 2
- ISSN:
- 0928-0219
- Page Range / eLocation ID:
- 219 to 231
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
-
null (Ed.)We prove well-posedness and regularity of solutions to a fractional diffusion porous media equation with a variable fractional order that may depend on the unknown solution. We present a linearly implicit time-stepping method to linearize and discretize the equation in time, and present rigorous analysis for the convergence of numerical solutions based on proved regularity results.more » « less
-
Abstract The modeling of nonlinear dynamical systems subject to strong and evolving nonsmooth nonlinearities is typically approached via integer-order differential equations. In this study, we present the possible application of variable-order (VO) fractional operators to a class of nonlinear lumped parameter models that have great practical relevance in mechanics and dynamics. Fractional operators are intrinsically multiscale operators that can act on both space- and time-dependent variables. Contrarily to their integer-order counterpart, fractional operators can have either fixed or VO. In the latter case, the order can be function of either independent or state variables. We show that when using VO equations to describe the response of dynamical systems, the order can evolve as a function of the response itself; therefore, allowing a natural and seamless transition between widely dissimilar dynamics. Such an intriguing characteristic allows defining governing equations for dynamical systems that are evolutionary in nature. Within this context, we present a physics-driven strategy to define VO operators capable of capturing complex and evolutionary phenomena. Specific examples include hysteresis in discrete oscillators and contact problems. Despite using simplified models to illustrate the applications of VO operators, we show numerical evidence of their unique modeling capabilities as well as their connection to more complex dynamical systems.more » « less
-
more » « less
Abstract It is known that the solutions to space‐fractional diffusion equations exhibit singularities near the boundary. Therefore, numerical methods discretized on the composite mesh, in which the mesh size is refined near the boundary, provide more precise approximations to the solutions. However, the coefficient matrices of the corresponding linear systems usually lose the diagonal dominance and are ill‐conditioned, which in turn affect the convergence behavior of the iteration methods.In this work we study a finite volume method for two‐sided fractional diffusion equations, in which a locally refined composite mesh is applied to capture the boundary singularities of the solutions. The diagonal blocks of the resulting three‐by‐three block linear system are proved to be positive‐definite, based on which we propose an efficient block Gauss–Seidel method by decomposing the whole system into three subsystems with those diagonal blocks as the coefficient matrices. To further accelerate the convergence speed of the iteration, we use T. Chan's circulant preconditioner31as the corresponding preconditioners and analyze the preconditioned matrices' spectra. Numerical experiments are presented to demonstrate the effectiveness and the efficiency of the proposed method and its strong potential in dealing with ill‐conditioned problems. While we have not proved the convergence of the method in theory, the numerical experiments show that the proposed method is convergent.
-
Oxygen isotopes are a well-known geochemical tool with applications to equilibrium thermometry, fluid tracing, and magma and ore petrogenesis. High-precision, high-spatial resolution oxygen isotope analysis by SIMS has also enabled the development of oxygen isotopes as a tool for geospeedometry. Here, we detail Fast Grain Boundary (FGB), an updated computational approach and software tool for determining time-temperature (T-t) histories through modeling of oxygen isotope diffusion. FBG models a rock system, rather than a single phase (cf. thermochronometry based on He, Ar, Pb), and has the potential to constrain continuous thermal histories over a wide range of temperatures, including at high temperature (500-800°C). The new FGB also allows for inversion of the FGB model to extract thermal histories from intragrain oxygen isotope zoning data using the Levenberg-Marquardt (LM) algorithm. Tests with synthetic datasets show that the LM algorithm is able to distinguish between simple linear cooling and more complex thermal histories containing, for instance, reheating events. Inversion of an actual oxygen isotope data set from titanite are consistent with the previously determined T-t path for the sample region, showing a brief period of >700 °C conditions, followed by cooling below 500 °C in <5 m.y.. However, the inversion suffers from a flat-bottomed minimum and does not produce a well-converged T-t path. These results point to analytical precision as a continuing challenge in recovering tightly constrained thermal histories for the real data set and emphasize the need for further development of high-precision microanalytical oxygen isotope standards. In the meantime, we use FGB modeling to explore sampling and analytical approaches that improve the resolution of inversion solutions for current analytical capabilities. For instance, inversion most successfully recovers a well constrained T-t path solution when SIMS analysis targets oxygen isotope gradients developed near grain rims, as opposed to oxygen isotope values in grain centers. Additional tests that probe the sensitivity of the inversion results to modal mineralogy and relative grain sizes suggest that careful targeting of samples in the field can enhance the recovery of unique T-t paths.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1515/jiip-2019-0040
