skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Operator-Based Uncertainty Quantification of Stochastic Fractional Partial Differential Equations
Abstract Fractional calculus provides a rigorous mathematical framework to describe anomalous stochastic processes by generalizing the notion of classical differential equations to their fractional-order counterparts. By introducing the fractional orders as uncertain variables, we develop an operator-based uncertainty quantification framework in the context of stochastic fractional partial differential equations (SFPDEs), subject to additive random noise. We characterize different sources of uncertainty and then, propagate their associated randomness to the system response by employing a probabilistic collocation method (PCM). We develop a fast, stable, and convergent Petrov–Galerkin spectral method in the physical domain in order to formulate the forward solver in simulating each realization of random variables in the sampling procedure.  more » « less
Award ID(s):
1923201
PAR ID:
10213731
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Journal of Verification, Validation and Uncertainty Quantification
Volume:
4
Issue:
4
ISSN:
2377-2158
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; (, Biology)
    Inverse problems involving time-fractional differential equations have become increasingly important for modeling systems with memory-dependent dynamics, particularly in biotransport and viscoelastic materials. Despite their potential, these problems remain challenging due to issues of stability, non-uniqueness, and limited data availability. Recent advancements in Physics-Informed Neural Networks (PINNs) offer a data-efficient framework for solving such inverse problems, yet most implementations are restricted to integer-order derivatives. In this work, we develop a PINN-based framework tailored for inverse problems involving time-fractional derivatives. We consider two representative applications: anomalous diffusion and fractional viscoelasticity. Using both synthetic and experimental datasets, we infer key physical parameters including the generalized diffusion coefficient and the fractional derivative order in the diffusion model and the relaxation parameters in a fractional Maxwell model. Our approach incorporates a customized residual loss function scaled by the standard deviation of observed data to enhance robustness. Even under 25% Gaussian noise, our method recovers model parameters with relative errors below 10%. Additionally, the framework accurately predicts relaxation moduli in porcine tissue experiments, achieving similar error margins. These results demonstrate the framework’s effectiveness in learning fractional dynamics from noisy and sparse data, paving the way for broader applications in complex biological and mechanical systems. 
    more » « less
  2. ; ; ; (, BIT numerical mathematics)
    In uncertainty quantification, it is commonly required to solve a forward model consisting of a partial differential equation (PDE) with a spatially varying uncertain coefficient that is represented as an affine function of a set of random variables, or parameters. Discretizing such models using stochastic Galerkin finite element methods (SGFEMs) leads to very high-dimensional discrete problems that can be cast as linear multi-term matrix equations (LMTMEs). We develop efficient computational methods for approximating solutions of such matrix equations in low rank. To do this, we follow an alternating energy minimization (AEM) framework, wherein the solution is represented as a product of two matrices, and approximations to each component are sought by solving certain minimization problems repeatedly. Inspired by proper generalized decomposition methods, the iterative solution algorithms we present are based on a rank-adaptive variant of AEM methods that successively computes a rank-one solution component at each step. We introduce and evaluate new enhancement procedures to improve the accuracy of the approximations these algorithms deliver. The efficiency and accuracy of the enhanced AEM methods is demonstrated through numerical experiments with LMTMEs associated with SGFEM discretizations of parameterized linear elliptic PDEs. 
    more » « less
  3. ; (, Journal of Computational and Nonlinear Dynamics)
    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
  4. ; ; (, ESAIM: Mathematical Modelling and Numerical Analysis)
    null (Ed.)
    In this paper, we propose a local discontinuous Galerkin (LDG) method for nonlinear and possibly degenerate parabolic stochastic partial differential equations, which is a high-order numerical scheme. It extends the discontinuous Galerkin (DG) method for purely hyperbolic equations to parabolic equations and shares with the DG method its advantage and flexibility. We prove the L 2 -stability of the numerical scheme for fully nonlinear equations. Optimal error estimates ( O ( h (k+1) )) for smooth solutions of semi-linear stochastic equations is shown if polynomials of degree k are used. We use an explicit derivative-free order 1.5 time discretization scheme to solve the matrix-valued stochastic ordinary differential equations derived from the spatial discretization. Numerical examples are given to display the performance of the LDG method. 
    more » « less
  5. (, Frontiers in applied mathematics and statistics)
    null (Ed.)
    A plethora of complex dynamical systems from disordered media to biological systems exhibit mathematical characteristics (e.g., long-range dependence, self-similar and power law magnitude increments) that are well-fitted by fractional partial differential equations (PDEs). For instance, some biological systems displaying an anomalous diffusion behavior, which is characterized by a non-linear mean-square displacement relation, can be mathematically described by fractional PDEs. In general, the PDEs represent various physical laws or rules governing complex dynamical systems. Since prior knowledge about the mathematical equations describing complex dynamical systems in biology, healthcare, disaster mitigation, transportation, or environmental sciences may not be available, we aim to provide algorithmic strategies to discover the integer or fractional PDEs and their parameters from system's evolution data. Toward deciphering non-trivial mechanisms driving a complex system, we propose a data-driven approach that estimates the parameters of a fractional PDE model. We study the space-time fractional diffusion model that describes a complex stochastic process, where the magnitude and the time increments are stable processes. Starting from limited time-series data recorded while the system is evolving, we develop a fractional-order moments-based approach to determine the parameters of a generalized fractional PDE. We formulate two optimization problems to allow us to estimate the arguments of the fractional PDE. Employing extensive simulation studies, we show that the proposed approach is effective at retrieving the relevant parameters of the space-time fractional PDE. The presented mathematical approach can be further enhanced and generalized to include additional operators that may help to identify the dominant rule governing the measurements or to determine the degree to which multiple physical laws contribute to the observed dynamics. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email