In this paper, we propose a new class of operator factorization methods to discretize the integral fractional Laplacian \begin{document}$$ (- \Delta)^\frac{{ \alpha}}{{2}} $$\end{document} for \begin{document}$$ \alpha \in (0, 2) $$\end{document}. One main advantage is that our method can easily increase numerical accuracy by using high-degree Lagrange basis functions, but remain its scheme structure and computer implementation unchanged. Moreover, it results in a symmetric (multilevel) Toeplitz differentiation matrix, enabling efficient computation via the fast Fourier transforms. If constant or linear basis functions are used, our method has an accuracy of \begin{document}$$ {\mathcal O}(h^2) $$\end{document}, while \begin{document}$$ {\mathcal O}(h^4) $$\end{document} for quadratic basis functions with \begin{document}$ h $$\end{document} a small mesh size. This accuracy can be achieved for any \begin{document}$$ \alpha \in (0, 2) $$\end{document} and can be further increased if higher-degree basis functions are chosen. Numerical experiments are provided to approximate the fractional Laplacian and solve the fractional Poisson problems. It shows that if the solution of fractional Poisson problem satisfies \begin{document}$$ u \in C^{m, l}(\bar{ \Omega}) $$\end{document} for \begin{document}$$ m \in {\mathbb N} $$\end{document} and \begin{document}$$ 0 < l < 1 $$\end{document}, our method has an accuracy of \begin{document}$$ {\mathcal O}(h^{\min\{m+l, \, 2\}}) $$\end{document} for constant and linear basis functions, while \begin{document}$$ {\mathcal O}(h^{\min\{m+l, \, 4\}}) $$\end{document}$ for quadratic basis functions. Additionally, our method can be readily applied to approximate the generalized fractional Laplacians with symmetric kernel function, and numerical study on the tempered fractional Poisson problem demonstrates its efficiency.
more »
« less
A unified framework for optimal control of fractional in time subdiffusive semilinear PDEs
We consider optimal control of fractional in time (subdiffusive, i.e., for \begin{document}$$ 0<\gamma <1 $$\end{document}) semilinear parabolic PDEs associated with various notions of diffusion operators in an unifying fashion. Under general assumptions on the nonlinearity we \begin{document}$$\mathsf{first\;show}$$\end{document} the existence and regularity of solutions to the forward and the associated \begin{document}$$\mathsf{backward\;(adjoint)}$$\end{document} problems. In the second part, we prove existence of optimal \begin{document}$$\mathsf{controls }$$\end{document} and characterize the associated \begin{document}$$\mathsf{first\;order}$$\end{document} optimality conditions. Several examples involving fractional in time (and some fractional in space diffusion) equations are described in detail. The most challenging obstacle we overcome is the failure of the semigroup property for the semilinear problem in any scaling of (frequency-domain) Hilbert spaces.
more »
« less
- PAR ID:
- 10345650
- Date Published:
- Journal Name:
- Discrete and Continuous Dynamical Systems - S
- Volume:
- 15
- Issue:
- 8
- ISSN:
- 1937-1632
- Page Range / eLocation ID:
- 1883
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
This paper investigates the global existence of weak solutions for the incompressible \begin{document}$ p $$\end{document}-Navier-Stokes equations in \begin{document}$$ \mathbb{R}^d $$\end{document} \begin{document}$$ (2\leq d\leq p) $$\end{document}. The \begin{document}$$ p $$\end{document}-Navier-Stokes equations are obtained by adding viscosity term to the \begin{document}$$ p $$\end{document}-Euler equations. The diffusion added is represented by the \begin{document}$$ p $$\end{document}-Laplacian of velocity and the \begin{document}$$ p $$\end{document}-Euler equations are derived as the Euler-Lagrange equations for the action represented by the Benamou-Brenier characterization of Wasserstein-\begin{document}$$ p $$\end{document}$ distances with constraint density to be characteristic functions.more » « less
-
We address the global existence and uniqueness of solutions for the anisotropically reduced 2D Kuramoto-Sivashinsky equations in a periodic domain with initial data \begin{document}$$ u_{01} \in L^2 $$\end{document} and \begin{document}$$ u_{02} \in H^{-1 + \eta} $$\end{document} for \begin{document}$$ \eta > 0 $$\end{document}.more » « less
-
For any finite horizon Sinai billiard map \begin{document}$ T $$\end{document} on the two-torus, we find \begin{document}$$ t_*>1 $$\end{document} such that for each \begin{document}$$ t\in (0,t_*) $$\end{document} there exists a unique equilibrium state \begin{document}$$ \mu_t $$\end{document} for \begin{document}$$ - t\log J^uT $$\end{document}, and \begin{document}$$ \mu_t $$\end{document} is \begin{document}$$ T $$\end{document}-adapted. (In particular, the SRB measure is the unique equilibrium state for \begin{document}$$ - \log J^uT $$\end{document}.) We show that \begin{document}$$ \mu_t $$\end{document} is exponentially mixing for Hölder observables, and the pressure function \begin{document}$$ P(t) = \sup_\mu \{h_\mu -\int t\log J^uT d \mu\} $$\end{document} is analytic on \begin{document}$$ (0,t_*) $$\end{document}. In addition, \begin{document}$$ P(t) $$\end{document} is strictly convex if and only if \begin{document}$$ \log J^uT $$\end{document} is not \begin{document}$$ \mu_t $$\end{document}-a.e. cohomologous to a constant, while, if there exist \begin{document}$$ t_a\ne t_b $$\end{document} with \begin{document}$$ \mu_{t_a} = \mu_{t_b} $$\end{document}, then \begin{document}$$ P(t) $$\end{document} is affine on \begin{document}$$ (0,t_*) $$\end{document}. An additional sparse recurrence condition gives \begin{document}$$ \lim_{t\downarrow 0} P(t) = P(0) $$\end{document}$.more » « less
-
We consider the linear third order (in time) PDE known as the SMGTJ-equation, defined on a bounded domain, under the action of either Dirichlet or Neumann boundary control \begin{document}$ g $$\end{document}. Optimal interior and boundary regularity results were given in [1], after [41], when \begin{document}$$ g \in L^2(0, T;L^2(\Gamma)) \equiv L^2(\Sigma) $$\end{document}, which, moreover, in the canonical case \begin{document}$$ \gamma = 0 $$\end{document}, were expressed by the well-known explicit representation formulae of the wave equation in terms of cosine/sine operators [19], [17], [24,Vol Ⅱ]. The interior or boundary regularity theory is however the same, whether \begin{document}$$ \gamma = 0 $$\end{document} or \begin{document}$$ 0 \neq \gamma \in L^{\infty}(\Omega) $$\end{document}, since \begin{document}$$ \gamma \neq 0 $$\end{document} is responsible only for lower order terms. Here we exploit such cosine operator based-explicit representation formulae to provide optimal interior and boundary regularity results with \begin{document}$$ g $$\end{document} "smoother" than \begin{document}$$ L^2(\Sigma) $$\end{document}, qualitatively by one unit, two units, etc. in the Dirichlet boundary case. To this end, we invoke the corresponding results for wave equations, as in [17]. Similarly for the Neumann boundary case, by invoking the corresponding results for the wave equation as in [22], [23], [37] for control smoother than \begin{document}$$ L^2(0, T;L^2(\Gamma)) $$\end{document}, and [44] for control less regular in space than \begin{document}$$ L^2(\Gamma) $$\end{document}$. In addition, we provide optimal interior and boundary regularity results when the SMGTJ equation is subject to interior point control, by invoking the corresponding wave equations results [42], [24,Section 9.8.2].more » « less
An official website of the United States government

