## Numerical Analysis Seminar

Date, time & place Tuesday 16:30 - 18:00 056Room #056 (Graduate School of Math. Sci. Bldg.)

Seminar information archive

### 2021/07/06

16:30-18:00   Online
Ken Hayami (National Institute of Informatics (Professor Emeritus))
Iterative solution methods for least squares problems and their applications
(Japanese)
[ Reference URL ]
https://forms.gle/B5Hwxa7o8F36hZKr7

### 2021/06/22

17:00-18:30   Online
Taiji Suzuki (The University of Tokyo)
On approximation ability and adaptivity of deep neural network (Japanese)
[ Reference URL ]
https://forms.gle/HwetNGXCzbCyMC7B7

### 2021/06/08

16:30-18:00   Online
Kohei Soga (Keio University)
Action minimizing random walks and numerical analysis of Hamilton-Jacobi equations (Japanese)
[ Reference URL ]
https://forms.gle/kjhqne4nV6fqEFWB8

### 2021/05/11

16:30-18:00   Online
Takuya Tsuchiya (Ehime University )
Topics on finite element error analysis on anisotropic meshes (Japanese)
[ Reference URL ]
https://forms.gle/CoaM4vSE1GvDRuDR6

### 2021/04/27

16:30-18:00   Online
Akitoshi Takayasu (University of Tsukuba)
Rigorous numerics for nonlinear heat equations in the complex plane of time (Japanese)
[ Reference URL ]
https://forms.gle/qW5ktphBB6dsh8Np7

### 2021/01/12

16:30-18:00   Online
Takaharu Yaguchi (Kobe University)
DGNet: Deep Energy-Based Modeling of Discrete-Time Physics and Related Topics (Japanese)
[ Reference URL ]
https://forms.gle/DpuhGupZ7NYbot5d7

### 2020/12/15

16:30-18:00   Online
Ming-Cheng Shiue (National Chiao Tung University)
Iterated pressure-correction projection methods for the 2d Navier-Stokes equations based on the scalar auxiliary variable approach (English)
[ Abstract ]
In this talk, the first-order iterated pressure-correction projection methods based on the scalar auxiliary variable approach is proposed and studied for the 2d Navier-Stokes equations and Boussinesq equations.
In the literature, enormous amounts of work have contributed to the study of numerical schemes for computing the Navier-Stokes equations. In general, two of the main numerical difficulties for solving Navier-Stokes equations are the incompressible condition and the nonlinear term. One of the approaches to deal with the incompressible condition is the so-called projection. The typical projection method only needs to solve the Poisson type of equations depending on the nonlinear term's treatment, which is efficient. However, the pressure-correction projection methods suffer from the splitting error, leading to spurious numerical boundary layers and the limitation of accuracy in time. In the literature, an iterated pressure-correction projection method has been proposed to overcome the difficulty.
As for the nonlinear term treatment, it is better to treat the nonlinear term explicitly so that one only requires to solve the corresponding linear system with constant coefficients at each time step. However, such treatment often results in a restricted time step due to the stable issue. Recently, the scalar auxiliary variable approach has been constructed to have an unconditional energy stable numerical scheme.
In this work, a new iterated pressure-correction projection method based on the scalar auxiliary variable's simple choice is proposed. We find that this new scheme can enjoy two properties, including reducing the splitting errors and having unconditional energy stability. The proofs of the energy stability and error convergence are provided and analyzed. Finally, numerical examples are provided to illustrate the theoretical work. This is joint work with Tony Chang.
[ Reference URL ]
https://forms.gle/y7w2nmaYtHNeoDSn8

### 2020/12/01

16:30-18:00   Online
Hiroyuki Sato (Kyoto University)
Conjugate gradient methods for optimization problems on manifolds (Japanese)
[ Reference URL ]
https://forms.gle/Ubeccm8neLkacjbk8

### 2020/10/27

16:30-18:00   Online
Buyang Li (The Hong Kong Polytechnic University)
Convergent evolving finite element algorithms for mean curvature flow and Willmore flow of closed surfaces (English)
[ Abstract ]
We construct evolving surface finite element methods for the mean curvature and Willmore flow through equivalently reformulating the original equations into coupled systems governing the evolution of surface position, velocity, normal vector and mean curvature. Then we prove $H^1$-norm convergence of the proposed evolving surface finite element methods for the reformulated systems, by combining stability estimates and consistency estimates. The stability analysis is based on the matrix–vector formulation of the finite element method and does not use geometric arguments. The geometry enters only into the consistency estimates. Numerical experiments illustrate and complement the theoretical results.
[1] https://doi.org/10.1007/s00211-019-01074-2
[2] https://arxiv.org/abs/2007.15257
[ Reference URL ]
https://forms.gle/HeuUxWLGa696KPvz8

### 2020/07/21

16:30-18:00   Online
Tomoya Kemmochi (Nagoya University)
Structure-preserving numerical schemes for constrained gradient flows of planar curves (Japanese)
[ Reference URL ]
https://forms.gle/3JiNEjWnrWLW8cFA9

### 2020/06/30

16:30-18:00   Online
Koya Sakakibara (Okayama University of Science)
Structure-preserving numerical methods for interface problems (Japanese)
[ Reference URL ]
https://forms.gle/ztK741vNdBT7hfGSA

### 2020/06/23

16:30-18:00   Online
Shun Sato (The University of Tokyo)
Linearly implicit and high-order conservative schemes for ordinary differential equations with a quadratic invariant (Japanese)
[ Reference URL ]
https://forms.gle/hvvvFLAhH1314UQK8

### 2020/01/20

16:50-18:20   Room #056 (Graduate School of Math. Sci. Bldg.)
Yves A. B. C. Barbosa (Politecnico di Milano)
Isogeometric Hierarchical Model Reduction: from analysis to patient-specific simulations (English)
[ Abstract ]
In the field of hemodynamics, numerical models have evolved to account for the demands in speed and accuracy of modern diagnostic medicine. In this context, we studied in detail Hierarchical Model Reduction technique combined with Isogeometric Analysis (HigaMOD), a technique recently developed in [Perotto, Reali, Rusconi and Veneziani (2017)]. HigaMod is a reduction procedure used to downscale models when the phenomenon at hand presents a preferential direction of flow, e.g., when modelling the blood flow in arteries or the water flow in a channel network. The method showed a significant improvement in reducing the computational power and simulation time, while giving enough information to analyze the problem at hand.

Recently, we focused our work in solving the ADR problem and the Stokes problem in a patient-specific framework. Specifically, we evaluate the computational efficiency of HigaMod in simulating the blood flow in coronary arteries and cerebral arteries. The main goal is to assess the
mprovement that 1D enriched models can provide, with respect to traditional full models, when dealing with demanding 3D CFD simulations. The results obtained, even though preliminary, are promising [Brandes, Barbosa and Perotto (2019); Brandes, Barbosa, Perotto and Suito (2020)].

### 2019/12/16

16:50-18:20   Room #117 (Graduate School of Math. Sci. Bldg.)
Yuki Ueda (The Hong Kong Polytechnic University)
A second-order stabilization method for linearizing and decoupling nonlinear parabolic systems (Japanese)
[ Abstract ]
We present a new time discretization method for strongly nonlinear parabolic systems. Our method is based on backward finite difference for the first derivative with second-order accuracy and the first-order linear discrete-time scheme for nonlinear systems which has been introduced by H. Murakawa. We propose a second-order stabilization method by combining these schemes.
Our error estimate requires testing the error equation by two test functions and showing $W^{1,\infty}$-boundedness which is proved by ($H^2$ or) $H^3$ energy estimate. We overcome the difficulty for establishing energy estimate by using the generating function technique which is popular in studying ordinary differential equations. Several numerical examples are provided to support the theoretical result.

### 2019/12/09

16:50-18:20   Room #117 (Graduate School of Math. Sci. Bldg.)
Xuefeng Liu (Niigata University)
Point-wise error estimation for the finite element solution to Poisson's equation --- new approach based on Kato-Fujita's method (Japanese)
[ Abstract ]
In 1950s, H. Fujita proposed a method to provide the upper and lower bounds in boundary value problems, which is based on the T*T theory of T. Kato about differential equations. Such a method can be regarded a different formulation of the hypercircle method from Prage-Synge's theorem.
Recently, the speaker extended Kato-Fujita's method to the case of the finite element solution of Poisson's equation and proposed a guaranteed point-wise error estimation. The newly proposed error estimation can be applied to problems defined over domains of general shapes along with general boundary conditions.

### 2019/09/25

16:30-18:00   Room #056 (Graduate School of Math. Sci. Bldg.)
Guanyu Zhou (Tokyo University of Science)
Finite volume method for the Keller-Segel system (Japanese)

### 2019/08/19

13:00-17:00   Room #122 (Graduate School of Math. Sci. Bldg.)
Eric Chung (The Chinese University of Hong Kong) 13:00-14:00
Staggered hybridisation for discontinuous Galerkin methods (英語)
[ Abstract ]
In this talk, we present a new staggered hybridization technique for discontinuous Galerkin methods to discretize linear elastodynamic equations and nonlinear Stokes equations. The idea of hybridization is used extensively in many discontinuous Galerkin methods, but the idea of staggered hybridization is new. Our new approach offers several advantages, namely energy conservation, high-order optimal convergence, preservation of symmetry for the stress tensor, block diagonal mass matrices as well as low dispersion error. The key idea is to use two staggered hybrid variables to enforce the continuity of the velocity and the continuity of the normal component of the stress tensor on a staggered mesh. We prove the stability and the convergence of the proposed scheme in both the semi-discrete and the fully-discrete settings. Numerical results confirm the optimal rate of convergence and show that the method has a superconvergent property for dispersion.
Feifei Jing (Northwestern Polytechnical University) 14:30-15:30
DG and HDG methods for the variational inequality problems (英語)
[ Abstract ]
There exist many numerical methods for solving the fluid dynamics equations, the main difference between them lies in the partitions of geometric domain and the discrete forms of governing equations. Due to the discontinuous piecewise polynomial subspaces, DG and HDG methods can be easily implemented on highly unstructured meshes, e.g. general polygonal mesh, and volume integrals could be calculated on physical elements, without reference elements and mappings between physical and reference elements. In this talk, DG and HDG methods employed to a class of variational inequality problems arising in hydrodynamics are studied. Some theoretical results will be shown, as well as the implementations of these methods are also put into practice.
Issei Oikawa (Hitotsubashi University) 16:00-16:30
A new HDG method using a hybridized flux (英語)
[ Abstract ]
We propose a new hybridizable discontinuous Galerkin (HDG) method for steady-state diffusion problems. In our method, both the trace and flux of the exact solution are hybridized. The Lehrenfeld-Schöberl stabilization is implicitly included in the method, so that the orders of convergence in all variables are optimal without postprocessing and computation of any projection. Numerical results are present to show the validation of our method.
Takahito Kashiwabara (The University of Tokyo) 16:30-17:00
Numerical approximation of the Stokes–Darcy problem using discontinuous linear elements (英語)
[ Abstract ]
We consider the Stokes–Darcy interface problem supplemented with the Beavers– Joseph–Saffman condition on the interface separating two domains. This condition allows for discontinuity in the tangential velocities and in the pressures along the interface. To effectively express it, we propose to use discontinuous linear finite elements to approximate all of the velocities/pressures in the Stokes/Darcy regions. The continuity of velocity in the normal direction is weakly enforced by adopting either the penalty method or Nitsche’s method. We present stability and error estimates for the proposed scheme, taking into account the situation where a curved interface is approximated by a polygonal curve or polyhedral surface.

### 2019/07/08

16:50-18:20   Room #056 (Graduate School of Math. Sci. Bldg.)
Takeru Matsuda (University of Tokyo)
Parameter estimation and discretization errors for ordinary differential models (Japanese)

### 2019/07/01

16:50-18:20   Room #117 (Graduate School of Math. Sci. Bldg.)
The effect of preventing scale formation by ceramic balls and its effect on the human body (Japanese)

### 2019/05/13

16:50-18:20   Room #056 (Graduate School of Math. Sci. Bldg.)
Kensuke Aishima (Hosei University)
Iterative refinement for symmetric eigenvalue problems (Japanese)

### 2019/04/22

16:50-18:20   Room #056 (Graduate School of Math. Sci. Bldg.)
Issei Oikawa (Hitotsubashi University )
Superconvergence of the HDG method (Japanese)

### 2019/04/08

16:50-18:20   Room #056 (Graduate School of Math. Sci. Bldg.)
Takeshi Takaishi (Musashino University)
Crack growth model of viscoelastic material with the phase field approach (Japanese)

### 2018/11/05

16:50-18:20   Room #002 (Graduate School of Math. Sci. Bldg.)
Hisashi Okamoto (Gakushuin University)
Tosio Kato as an applied mathematician (Japanese)
[ Abstract ]
Tosio Kato (1917-1999) is nowadays considered to be a rigorous analyst or theorist. Many people consider his contributions in quantum mechanics to be epoch-making, his work on nonlinear partial differential equations elegant and inspiring. However, around the time when he visited USA for the first time in 1954, he was studying problems of applied mathematics, too, notably numerical computation of eigenvalues. I wish to shed light on the historical background of his study of applied mathematics. This is a joint work with Prof. Hiroshi Fujita.

### 2018/10/22

16:50-18:20   Room #002 (Graduate School of Math. Sci. Bldg.)
Kensuke Aihara (Tokyo City University)
Residual smoothing technique for short-recurrence Krylov subspace methods (Japanese)

### 2018/10/15

16:50-18:20   Room #002 (Graduate School of Math. Sci. Bldg.)
Takeyuki Nagasawa (Saitama University)
Möbius invariant discretizations and decomposition of the Möbius energy (Japanese)