Norikazu Saito

Tomoya Kemmochi, Norikazu Saito

Maximal regularity is a fundamental concept in the theory of partial differential equations. In this paper, we establish a fully discrete version of maximal regularity for a parabolic equation. We derive various stability results in $L^p(0,T;L^q(Ω))$ norm, $p,q\in (1,\infty)$ for the finite element approximation with the mass-lumping to the linear heat equation. Our method of analysis is an operator theoretical one using pure imaginary powers of operators and might be a discrete version of G.~Dore and A.~Venni (On the closedness of the sum of two closed operators. \emph{Math.\ Z.}, 196(2):189--201, 1987). As an application, optimal order error estimates in that norm are proved. Furthermore, we study the finite element approximation for semilinear heat equations with locally Lipschitz continuous nonlinearity and offer a new method for deriving optimal order error estimates. Some interesting auxiliary results including discrete Gagliardo-Nirenberg and Sobolev inequalities are also presented.

Norikazu Saito

This note was prepared for a lecture given at Kyoto University (RIMS Workshop: "The State of the Art in Numerical Analysis: Theory, Methods, and Applications", November 8-10, 2017). That lecture described the variational analysis of the discontinuous Galerkin time-stepping method for parabolic equations based on an earlier paper by the author (arXiv:1710.10543). I also presented the Banach-Necas-Babuska (BNB) Theorem or the Babuska-Lax-Milgram (BLM) Theorem as the key theorem of our analysis. For proof of the BNB theorem, it is useful to introduce the minimum modulus of operators by T. Kato. This note presents a review of the proofs of Closed Range Theorem and BNB Theorem following the idea of Kato. Moreover, I present an application to BNB theorem to parabolic equations. The well-posedness is proved by BNB theorem. This note is not an original research paper. It includes no new results. This is a revised manuscript and several incorrect descriptions in the original version are fixed.

Masasru Miyashita, Norikazu Saito

New hybridized discontinuous Galerkin (HDG) methods for the interface problem for elliptic equations are proposed. Unknown functions of our schemes are $u_h$ in elements and $\hat{u}_h$ on inter-element edges. That is, we formulate our schemes without introducing the flux variable. Our schemes naturally satisfy the Galerkin orthogonality. The solution $u$ of the interface problem under consideration may not have a sufficient regularity, say $u|_{Ω_1}\in H^2(Ω_1)$ and $u|_{Ω_2}\in H^2(Ω_2)$, where $Ω_1$ and $Ω_2$ are subdomains of the whole domain $Ω$ and $Γ=\partialΩ_1\cap\partialΩ_2$ implies the interface. We study the convergence, assuming $u|_{Ω_1}\in H^{1+s}(Ω_1)$ and $u|_{Ω_2}\in H^{1+s}(Ω_2)$ for some $s\in (1/2,1]$, where $H^{1+s}$ denotes the fractional order Sobolev space. Consequently, we succeed in deriving optimal order error estimates in an HDG norm and the $L^2$ norm. Numerical examples to validate our results are also presented.

Norikazu Saito, Yoshiki Sugitani

Convergence results for the immersed boundary method applied to a model Stokes problem with the homogeneous Dirichlet boundary condition are presented. As a discretization method, we deal with the finite element method. First, the immersed force field is approximated using a regularized delta function and its error in the $W^{-1,p}$ norm is examined for $1\le p<n/(n-1)$, $n$ being the space dimension. Then, we consider the immersed boundary discretization of the Stokes problem and study the regularization and discretization errors separately. Consequently, error estimate of order $h^{1-α}$ in the $W^{1,1}\times L^1$ norm for the velocity and pressure is derived, where $α$ is an arbitrarily small positive number. Error estimate of order $h^{1-α}$ in the $L^r$ norm for the velocity is also derived with $r=n/(n-1-α)$. The validity of those theoretical results are confirmed by numerical examples.

Yuki Ueda, Norikazu Saito

The Nitsche method is a method of "weak imposition" of the inhomogeneous Dirichlet boundary conditions for partial differential equations. This paper explains stability and convergence study of the Nitsche method applied to evolutionary diffusion-advection-reaction equations. We mainly discuss a general space semidiscrete scheme including not only the standard finite element method but also Isogeometric Analysis. Our method of analysis is a variational one that is a popular method for studying elliptic problems. The variational method enables us to obtain the best approximation property directly. Actually, results show that the scheme satisfies the inf-sup condition and Galerkin orthogonality. Consequently, the optimal order error estimates in some appropriate norms are proven under some regularity assumptions on the exact solution. We also consider a fully discretized scheme using the backward Euler method. Numerical example demonstrate the validity of those theoretical results.

Yuki Chiba, Norikazu Saito

We derive several $L^\infty$ error estimates for the symmetric interior penalty (SIP) discontinuous Galerkin (DG) method applied to the Poisson equation in a two-dimensional polygonal domain. Both local and global estimates are examined. The weak maximum principle (WMP) for the discrete harmonic function is also established. We prove our $L^\infty$ estimates using this WMP and several $W^{2,p}$ and $W^{1,1}$ estimates for the Poisson equation. Numerical examples to validate our results are also presented.

Daisuke Inoue, Yuji Ito, Takahito Kashiwabara et al.

When controlling multi-agent systems, the trade-off between performance and scalability is a major challenge. Here, we address this difficulty by using mean field games (MFGs), which is a framework that deduces the macroscopic dynamics describing the density profile of agents from their microscopic dynamics. To effectively use the MFG, we propose a model predictive MFG (MP-MFG), which estimates the agent population density profile with using kernel density estimation and manages the input generation with model predictive control. The proposed MP-MFG generates control inputs by monitoring the agent population at each time step, and thus achieves higher robustness than the conventional MFG. Numerical results show that the MP-MFG outperforms the MFG when the agent model has modeling errors or the number of agents in the system is small.