The aim of this paper is to develop high-order methods for solving time-fractional partial differential equations. The proposed high-order method is based on high-order finite element method for ...space and finite difference method for time. Optimal convergence rate
O
(
(
Δ
t
)
2
−
α
+
N
−
r
)
is proved for the
(
r
−
1
)
th-order finite element method
(
r
≥
2
)
.
High order accurate weighted essentially nonoscillatory (WENO) schemes are relatively new but have gained rapid popularity in numerical solutions of hyperbolic partial differential equations (PDEs) ...and other convection dominated problems. The main advantage of such schemes is their capability to achieve arbitrarily high order formal accuracy in smooth regions while maintaining stable, nonoscillatory, and sharp discontinuity transitions. The schemes are thus especially suitable for problems containing both strong discontinuities and complex smooth solution features. WENO schemes are robust and do not require the user to tune parameters. At the heart of the WENO schemes is actually an approximation procedure not directly related to PDEs, hence the WENO procedure can also be used in many non-PDE applications. In this paper we review the history and basic formulation of WENO schemes, outline the main ideas in using WENO schemes to solve various hyperbolic PDEs and other convection dominated problems, and present a collection of applications in areas including computational fluid dynamics, computational astronomy and astrophysics, semiconductor device simulation, traffic flow models, computational biology, and some non-PDE applications. Finally, we mention a few topics concerning WENO schemes that are currently under investigation.
This paper is concerned with the existence of mild solutions for a class of impulsive fractional partial semilinear differential equations. Some errors in Mophou (2010)
2 are corrected, and some ...previous results are generalized.
High-dimensional PDEs have been a longstanding computational challenge. We propose to solve high-dimensional PDEs by approximating the solution with a deep neural network which is trained to satisfy ...the differential operator, initial condition, and boundary conditions. Our algorithm is meshfree, which is key since meshes become infeasible in higher dimensions. Instead of forming a mesh, the neural network is trained on batches of randomly sampled time and space points. The algorithm is tested on a class of high-dimensional free boundary PDEs, which we are able to accurately solve in up to 200 dimensions. The algorithm is also tested on a high-dimensional Hamilton–Jacobi–Bellman PDE and Burgers' equation. The deep learning algorithm approximates the general solution to the Burgers' equation for a continuum of different boundary conditions and physical conditions (which can be viewed as a high-dimensional space). We call the algorithm a “Deep Galerkin Method (DGM)” since it is similar in spirit to Galerkin methods, with the solution approximated by a neural network instead of a linear combination of basis functions. In addition, we prove a theorem regarding the approximation power of neural networks for a class of quasilinear parabolic PDEs.
•We develop a deep learning algorithm for solving high-dimensional PDEs.•The algorithm is meshfree, which is key since meshes become infeasible in higher dimensions.•We accurately solve a class of high-dimensional free boundary PDEs in up to 200 dimensions.•We prove a theorem regarding the approximation power of neural networks for quasilinear PDEs.
We present an unconditionally energy stable finite difference scheme for the Modified Phase Field Crystal equation, a generalized damped wave equation for which the usual Phase Field Crystal equation ...is a special degenerate case. The method is based on a convex splitting of a discrete pseudoenergy and is semi-implicit. The equation at the implicit time level is nonlinear but represents the gradient of a strictly convex function and is thus uniquely solvable, regardless of time step-size. We present a local-in-time error estimate that ensures the pointwise convergence of the scheme.
This work proposes and analyzes a stochastic collocation method for solving elliptic partial differential equations with random coefficients and forcing terms. These input data are assumed to depend ...on a finite number of random variables. The method consists of a Galerkin approximation in space and a collocation in the zeros of suitable tensor product orthogonal polynomials (Gauss points) in the probability space, and naturally leads to the solution of uncoupled deterministic problems as in the Monte Carlo approach. It treats easily a wide range of situations, such as input data that depend nonlinearly on the random variables, diffusivity coefficients with unbounded second moments, and random variables that are correlated or even unbounded. We provide a rigorous convergence analysis and demonstrate exponential convergence of the "probability error" with respect to the number of Gauss points in each direction of the probability space, under some regularity assumptions on the random input data. Numerical examples show the effectiveness of the method. Finally, we include a section with developments posterior to the original publication of this work. There we review sparse grid stochastic collocation methods, which are effective collocation strategies for problems that depend on a moderately large number of random variables.
•Deep auto-regressive dense encoder-decoder surrogate for predicting transient PDEs.•Physics-constrained learning enables the model to learn dynamics without training data.•A Bayesian framework is ...proposed for interpretable uncertainty quantification of the models' predictions at each time-step.•The auto-regressive model is tested on several non-linear PDE systems with features including turbulence and shock discontinuities.
In recent years, deep learning has proven to be a viable methodology for surrogate modeling and uncertainty quantification for a vast number of physical systems. However, in their traditional form, such models can require a large amount of training data. This is of particular importance for various engineering and scientific applications where data may be extremely expensive to obtain. To overcome this shortcoming, physics-constrained deep learning provides a promising methodology as it only utilizes the governing equations. In this work, we propose a novel auto-regressive dense encoder-decoder convolutional neural network to solve and model non-linear dynamical systems without training data at a computational cost that is potentially magnitudes lower than standard numerical solvers. This model includes a Bayesian framework that allows for uncertainty quantification of the predicted quantities of interest at each time-step. We rigorously test this model on several non-linear transient partial differential equation systems including the turbulence of the Kuramoto-Sivashinsky equation, multi-shock formation and interaction with 1D Burgers' equation and 2D wave dynamics with coupled Burgers' equations. For each system, the predictive results and uncertainty are presented and discussed together with comparisons to the results obtained from traditional numerical analysis methods.
This paper employs a Takagi-Sugeno (T-S) fuzzy partial differential equation (PDE) model to solve the problem of sampled-data exponential stabilization in the sense of spatial ∥·∥ ∞ for a class of ...nonlinear parabolic distributed parameter systems (DPSs), where only a few actuators and sensors are discretely distributed in space. Initially, a T-S fuzzy PDE model is assumed to be derived by the sector nonlinearity method to accurately describe complex spatiotemporal dynamics of the nonlinear DPSs. Subsequently, a static sampled-data fuzzy local state feedback controller is constructed based on the T-S fuzzy PDE model. By constructing an appropriate Lyapunov-Krasovskii functional candidate and employing vector-valued Wirtinger's inequalities, a variation of vector-valued Poincaré-Wirtinger inequality in one-dimensional spatial domain, as well as a vector-valued Agmon's inequality, it is shown that the suggested sampled-data fuzzy controller exponentially stabilizes the nonlinear DPSs in the sense of ∥·∥ ∞ , if sufficient conditions presented in term of standard linear matrix inequalities (LMIs) are fulfilled. Moreover, an LMI relaxation technique is utilized to enhance exponential stabilization ability of the suggested sampled-data fuzzy controller. Finally, the satisfactory and better performance of the suggested sampled-data fuzzy controller are demonstrated by numerical simulation results of two examples.
Simulation of multiphase flow in porous media is crucial for the effective management of subsurface energy and environment-related activities. The numerical simulators used for modeling such ...processes rely on spatial and temporal discretization of the governing mass and energy balance partial-differential equations (PDEs) into algebraic systems via finite-difference/volume/element methods. These simulators usually require dedicated software development and maintenance, and suffer low efficiency from a runtime and memory standpoint for problems with multi-scale heterogeneity, coupled-physics processes or fluids with complex phase behavior. Therefore, developing cost-effective, data-driven models can become a practical choice, and in this work, we choose deep learning approaches as they can handle high dimensional data and accurately predict state variables with strong nonlinearity. In this paper, we describe a gradient-based deep neural network (GDNN) constrained by the physics related to multiphase flow in porous media. We tackle the nonlinearity of flow in porous media induced by rock heterogeneity, fluid properties, and fluid-rock interactions by decomposing the nonlinear PDEs into a dictionary of elementary differential operators. We use a combination of operators to handle rock spatial heterogeneity and fluid flow by advection. Since the augmented differential operators are inherently related to the physics of fluid flow, we treat them as first principles prior knowledge to regularize the GDNN training. We use the example of pressure management at geologic CO2 storage sites, where CO2 is injected in saline aquifers and brine is produced, and apply GDNN to construct a predictive model that is trained with physics-based simulation data and emulates the physics process. We demonstrate that GDNN can effectively predict the nonlinear patterns of subsurface responses, including the temporal and spatial evolution of the pressure and CO2 saturation plumes. We also successfully extend the GDNN to convolutional neural network (CNN), namely gradient-based CNN (GCNN), and validate its capability to improve the prediction accuracy. GDNN has great potential to tackle challenging problems that are governed by highly nonlinear physics and enable the development of data-driven models with higher fidelity.
•A workflow based on Gradient-based deep neural network (GDNN) emulates multiphase flow in porous media with high fidelity.•Differential operators derived from the governing PDEs are used as first principle prior knowledge to regularize the training of GDNN.•GDNN is further successfully extended to image-based approaches such as convolutional neural networks (CNN).