An official website of the United States government
Abstract We propose a stochastic Galerkin method using sparse wavelet bases for the Boltzmann equation with multi-dimensional random inputs. Themethod uses locally supported piecewise polynomials as an orthonormal basis of the random space. By a sparse approach, only a moderate number of basis functions is required to achieve good accuracy in multi-dimensional random spaces. We discover a sparse structure of a set of basis-related coefficients, which allows us to accelerate the computation of the collision operator. Regularity of the solution of the Boltzmann equation in the random space and an accuracy result of the stochastic Galerkin method are proved in multi-dimensional cases. The efficiency of the method is illustrated by numerical examples with uncertainties from the initial data, boundary data and collision kernel.
more »
« less
A stochastic collocation method based on sparse grids for a stochastic Stokes-Darcy model
In this paper, we develop a sparse grid stochastic collocation method to improve the computational efficiency in handling the steady Stokes-Darcy model with random hydraulic conductivity. To represent the random hydraulic conductivity, the truncated Karhunen-Loève expansion is used. For the discrete form in probability space, we adopt the stochastic collocation method and then use the Smolyak sparse grid method to improve the efficiency. For the uncoupled deterministic subproblems at collocation nodes, we apply the general coupled finite element method. Numerical experiment results are presented to illustrate the features of this method, such as the sample size, convergence, and randomness transmission through the interface.
more »
« less
- Award ID(s):
- 1722647
- PAR ID:
- 10304075
- Date Published:
- Journal Name:
- Discrete and Continuous Dynamical Systems - Series S
- Volume:
- 0
- Issue:
- 0
- ISSN:
- 1937-1632
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Summary We propose and analyze an efficient ensemble algorithm with artificial compressibility (AC) for fast decoupled computation of multiple realizations of the stochastic Stokes‐Darcy model with random hydraulic conductivity (including the one in the interface conditions), source terms, and initial conditions. The solutions are found by solving three smaller decoupled subproblems with two common time‐independent coefficient matrices for all realizations, which significantly improves the efficiency for both assembling and solving the matrix systems. The fully coupled Stokes‐Darcy system can be first decoupled into two smaller subphysics problems by the idea of the partitioned time stepping, which reduces the size of the linear systems and allows parallel computing for each subphysics problem. The AC further decouples the velocity and pressure which further reduces storage requirements and improves computational efficiency. We prove the long time stability and the convergence for this new ensemble method. Three numerical examples are presented to support the theoretical results and illustrate the features of the algorithm, including the convergence, stability, efficiency, and applicability.more » « less
-
null (Ed.)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
-
Sparse learning models have shown promising performance in the high dimensional machine learning applications. The main challenge of sparse learning models is how to optimize it efficiently. Most existing methods solve this problem by relaxing it as a convex problem, incurring large estimation bias. Thus, the sparse learning model with nonconvex constraint has attracted much attention due to its better performance. But it is difficult to optimize due to the non-convexity.In this paper, we propose a linearly convergent stochastic second-order method to optimize this nonconvex problem for large-scale datasets. The proposed method incorporates second-order information to improve the convergence speed. Theoretical analysis shows that our proposed method enjoys linear convergence rate and guarantees to converge to the underlying true model parameter. Experimental results have verified the efficiency and correctness of our proposed method.more » « less
-
A numerical inverse method called FlowPaths is presented to solve for the hydraulic conductivity field of an isotropic heterogeneous porous medium from a known specific discharge field (and constant-head boundary conditions). This method makes possible a new approach to reactive transport experiments, aimed at understanding the dynamic spatial and temporal evolution of hydraulic conductivity, which simultaneously record the evolving reaction and the evolving flow geometry. This inverse method assumes steady, two-dimensional flow through a square matrix of grid blocks. A graph-theoretical approach is used to find a set of flow paths through the porous medium using the known components of the specific discharge, where every vertex is traversed by at least one path from the upstream high-head boundary to the downstream low-head boundary. Darcy’s law is used to create an equation for the unknown head drop across each edge. Summation of these edge equations along each path through the network generates a set of linearly independent head-drop equations that is solved directly for the hydraulic conductivity field. FlowPaths is verified by generating 12,740 hydraulic conductivity fields of varying size and heterogeneity, calculating the corresponding specific discharge field for each, and then using that specific discharge field to estimate the underlying hydraulic conductivity field. When estimates from FlowPaths are compared to the simulated hydraulic conductivity fields, the inverse method is demonstrated to be accurate and numerically stable. Accordingly, within certain limitations, FlowPaths can be used in field or laboratory applications to find hydraulic conductivity from a known velocity field.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3934/dcdss.2021104
