-
Inverse norm estimation of perturbed Laplace operators and corresponding eigenvalue problems
Authors:
Kouta Sekine,
Kazuaki Tanaka,
Shin'ichi Oishi
Abstract:
In numerical existence proofs for solutions of the semi-linear elliptic system, evaluating the norm of the inverse of a perturbed Laplace operator plays an important role. We reveal an eigenvalue problem to design a method for verifying the invertibility of the operator and evaluating the norm of its inverse based on Liu's method and the Temple-Lehman-Goerisch method. We apply the inverse-norm's e…
▽ More
In numerical existence proofs for solutions of the semi-linear elliptic system, evaluating the norm of the inverse of a perturbed Laplace operator plays an important role. We reveal an eigenvalue problem to design a method for verifying the invertibility of the operator and evaluating the norm of its inverse based on Liu's method and the Temple-Lehman-Goerisch method. We apply the inverse-norm's estimation to the Dirichlet boundary value problem of the Lotka-Volterra system with diffusion terms and confirm the efficacy of our method.
△ Less
Submitted 30 March, 2021; v1 submitted 4 October, 2019;
originally announced October 2019.
-
A new formulation for the numerical proof of the existence of solutions to elliptic problems
Authors:
Kouta Sekine,
Mitsuhiro T. Nakao,
Shin'ichi Oishi
Abstract:
Infinite-dimensional Newton methods can be effectively used to derive numerical proofs of the existence of solutions to partial differential equations (PDEs). In computer-assisted proofs of PDEs, the original problem is transformed into the infinite Newton-type fixed point equation $w = - {\mathcal L}^{-1} {\mathcal F}(\hat{u}) + {\mathcal L}^{-1} {\mathcal G}(w)$, where ${\mathcal L}$ is a linear…
▽ More
Infinite-dimensional Newton methods can be effectively used to derive numerical proofs of the existence of solutions to partial differential equations (PDEs). In computer-assisted proofs of PDEs, the original problem is transformed into the infinite Newton-type fixed point equation $w = - {\mathcal L}^{-1} {\mathcal F}(\hat{u}) + {\mathcal L}^{-1} {\mathcal G}(w)$, where ${\mathcal L}$ is a linearized operator, ${\mathcal F}(\hat{u})$ is a residual, and ${\mathcal G}(w)$ is a local Lipschitz term. Therefore, the estimations of $\| {\mathcal L}^{-1} {\mathcal F}(\hat{u}) \|$ and $\| {\mathcal L}^{-1}{\mathcal G}(w) \|$ play major roles in the verification procedures. In this paper, using a similar concept as the `Schur complement' for matrix problems, we represent the inverse operator ${\mathcal L}^{-1}$ as an infinite-dimensional operator matrix that can be decomposed into two parts, one finite dimensional and one infinite dimensional. This operator matrix yields a new effective realization of the infinite-dimensional Newton method, enabling a more efficient verification procedure compared with existing methods for the solution of elliptic PDEs. We present some numerical examples that confirm the usefulness of the proposed method. Related results obtained from the representation of the operator matrix as ${\mathcal L}^{-1}$ are presented in the appendix.
△ Less
Submitted 1 October, 2019;
originally announced October 2019.
-
Estimation of Sobolev embedding constant on a domain dividable into bounded convex domains
Authors:
Makoto Mizuguchi,
Kazuaki Tanaka,
Kouta Sekine,
Shin'ichi Oishi
Abstract:
This paper is concerned with an explicit value of the embedding constant from $W^{1,q}(Ω)$ to $L^{p}(Ω)$ for a bounded domain $Ω\subset\mathbb{R}^N~(N\in\mathbb{N})$, where $1\leq q\leq p\leq \infty$. To obtain this value, we previously proposed a formula for estimating the embedding constant on bounded and unbounded Lipschitz domains by estimating the norm of Stein's extension operator, in the ar…
▽ More
This paper is concerned with an explicit value of the embedding constant from $W^{1,q}(Ω)$ to $L^{p}(Ω)$ for a bounded domain $Ω\subset\mathbb{R}^N~(N\in\mathbb{N})$, where $1\leq q\leq p\leq \infty$. To obtain this value, we previously proposed a formula for estimating the embedding constant on bounded and unbounded Lipschitz domains by estimating the norm of Stein's extension operator, in the article (K. Tanaka, K. Sekine, M. Mizuguchi, and S. Oishi, Estimation of Sobolev-type embedding constant on domains with minimally smooth boundary using extension operator, Journal of Inequalities and Applications, Vol. 389, pp. 1-23, 2015). This formula is also applicable to a domain that can be divided into Lipschitz domains. However, the values computed by the previous formula are very large. In this paper, we propose several sharper estimations of the embedding constant on a bounded domain that can be divided into convex domains.
△ Less
Submitted 16 September, 2016;
originally announced September 2016.
-
Rigorous numerical enclosures for positive solutions of Lane-Emden's equation with sub-square exponents
Authors:
Kazuaki Tanaka,
Michael Plum,
Kouta Sekine,
Masahide Kashiwagi,
Shin'ichi Oishi
Abstract:
The purpose of this paper is to obtain rigorous numerical enclosures for solutions of Lane-Emden's equation $-Δu=|u|^{p-1} u$ with homogeneous Dirichlet boundary conditions. We prove the existence of a nondegenerate solution $u$ nearby a numerically computed approximation $\hat{u}$ together with an explicit error bound, i.e., a bound for the difference between $ u $ and $\hat{u}$. In particular, w…
▽ More
The purpose of this paper is to obtain rigorous numerical enclosures for solutions of Lane-Emden's equation $-Δu=|u|^{p-1} u$ with homogeneous Dirichlet boundary conditions. We prove the existence of a nondegenerate solution $u$ nearby a numerically computed approximation $\hat{u}$ together with an explicit error bound, i.e., a bound for the difference between $ u $ and $\hat{u}$. In particular, we focus on the sub-square case in which $1<p<2$ so that the derivative $p|u|^{p-1}$ of the nonlinearity $|u|^{p-1} u$ is not Lipschitz continuous. In this case, it is problematic to apply the classical Newton-Kantorovich theorem for obtaining the existence proof, and moreover several difficulties arise in the procedures to obtain numerical integrations rigorously. We design a method for enclosing the required integrations explicitly, proving the existence of a desired solution based on a generalized Newton-Kantorovich theorem. A numerical example is presented where an explicit solution-enclosure is obtained for $ p=3/2 $ on the unit square domain $Ω=(0,1)^2$.
△ Less
Submitted 7 September, 2021; v1 submitted 15 July, 2016;
originally announced July 2016.
-
Numerical verification method for positiveness of solutions to elliptic equations
Authors:
Kazuaki Tanaka,
Kouta Sekine,
Shin'ichi Oishi
Abstract:
In this paper, we propose a numerical method for verifying the positiveness of solutions to semilinear elliptic equations. We provide a sufficient condition for a solution to an elliptic equation to be positive in the domain of the equation, which can be checked numerically without requiring a complicated computation. We present some numerical examples.
In this paper, we propose a numerical method for verifying the positiveness of solutions to semilinear elliptic equations. We provide a sufficient condition for a solution to an elliptic equation to be positive in the domain of the equation, which can be checked numerically without requiring a complicated computation. We present some numerical examples.
△ Less
Submitted 3 July, 2016; v1 submitted 13 June, 2016;
originally announced June 2016.
-
Sharp numerical inclusion of the best constant for embedding $H_{0}^{1}(Ω) \hookrightarrow L^{p}(Ω)$ on bounded convex domain
Authors:
Kazuaki Tanaka,
Kouta Sekine,
Makoto Mizuguchi,
Shin'ichi Oishi
Abstract:
In this paper, we propose a verified numerical method for obtaining a sharp inclusion of the best constant for the embedding $H_{0}^{1}(Ω) \hookrightarrow L^{p}(Ω)$ on bounded convex domain in $\mathbb{R}^{2}$. We estimate the best constant by computing the corresponding extremal function using a verified numerical computation. Verified numerical inclusions of the best constant on a square domain…
▽ More
In this paper, we propose a verified numerical method for obtaining a sharp inclusion of the best constant for the embedding $H_{0}^{1}(Ω) \hookrightarrow L^{p}(Ω)$ on bounded convex domain in $\mathbb{R}^{2}$. We estimate the best constant by computing the corresponding extremal function using a verified numerical computation. Verified numerical inclusions of the best constant on a square domain are presented.
△ Less
Submitted 3 November, 2020; v1 submitted 18 March, 2015;
originally announced March 2015.
-
Estimation of the Sobolev embedding constant on domains with minimally smooth boundary
Authors:
Kazuaki Tanaka,
Kouta Sekine,
Makoto Mizuguchi,
Shin'ichi Oishi
Abstract:
In this paper, we propose a method for estimating the Sobolev type embedding constant on a domain with minimally smooth boundary. We estimate the embedding constant by constructing an extension operator and computing its operator norm. We also present some examples of estimating the embedding constant for certain domains.
In this paper, we propose a method for estimating the Sobolev type embedding constant on a domain with minimally smooth boundary. We estimate the embedding constant by constructing an extension operator and computing its operator norm. We also present some examples of estimating the embedding constant for certain domains.
△ Less
Submitted 9 June, 2015; v1 submitted 22 November, 2014;
originally announced November 2014.
-
Gate depletion of an InSb two-dimensional electron gas
Authors:
M. M. Uddin,
H. W. Liu,
K. F. Yang,
K. Nagase,
K. Sekine,
C. K. Gaspe,
T. D. Mishima,
M. B. Santos,
Y. Hirayama
Abstract:
We investigated the gate control of a two-dimensional electron gas (2DEG) confined to InSb quantum wells with an Al2O3 gate dielectric formed by atomic layer deposition on a surface layer of Al0.1In0.9Sb or InSb. The wider bandgap of Al0.1In0.9Sb compared to InSb resulted in a linear, sharp, and non-hysteretic response of the 2DEG density to gate bias in the structure with an Al0.1In0.9Sb surface…
▽ More
We investigated the gate control of a two-dimensional electron gas (2DEG) confined to InSb quantum wells with an Al2O3 gate dielectric formed by atomic layer deposition on a surface layer of Al0.1In0.9Sb or InSb. The wider bandgap of Al0.1In0.9Sb compared to InSb resulted in a linear, sharp, and non-hysteretic response of the 2DEG density to gate bias in the structure with an Al0.1In0.9Sb surface layer. In contrast, a nonlinear, slow, and hysteretic (nonvolatile-memory-like) response was observed in the structure with an InSb surface layer. The 2DEG with the Al0.1In0.9Sb surface layer was completely depleted by application of a small gate voltage (-0.9 V).
△ Less
Submitted 18 September, 2013;
originally announced September 2013.