BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Takuya Tsuchiya

Takuya Tsuchiya contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
gr-qcmath.NANumerical Analysisastro-ph.EPmath-phmath.DSmath.MPquant-ph

Trust snapshot

Quick read

Trust 21 - Emerging
12works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2022arXiv

Anisotropic interpolation error estimates using a new geometric parameter

We present precise anisotropic interpolation error estimates for smooth functions using a new geometric parameter and derive inverse inequalities on anisotropic meshes. In our theory, the interpolation error is bounded in terms of the diameter of a simplex and the geometric parameter. Imposing additional assumptions makes it possible to obtain anisotropic error estimates. This paper also includes corrections to an error in Theorem 2 of our previous paper, "General theory of interpolation error estimates on anisotropic meshes" (Japan Journal of Industrial and Applied Mathematics, 38 (2021) 163-191).

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Lectures on error analysis of interpolation on simplicial triangulations without the shape-regularity assumption Part 1: Lagrange interpolation on triangles

In the error analysis of finite element methods, the shape regularity assumption on triangulations is typically imposed to obtain a priori error estimations. In practical computations, however, very thin or degenerated elements that violate the shape regularity assumption may appear when we use adaptive mesh refinement. In this manuscript, we attempt to establish an error analysis approach without the shape regularity assumption on triangulations. We have presented several papers on the error analysis of finite element methods on non-shape regular triangulations. The main points in these papers are that, in the error estimates of finite element methods, the circumradius of the triangles is one of the most important factors. The purpose of this manuscript is to provide a simple and plain explanation of the results to researchers and, in particular, graduate students who are interested in the subject. Therefore, this manuscript is not intended to be a research paper. We hope that, in the future, it will be merged into a textbook on the mathematical theory of the finite element methods.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Numerical simulations of semilinear Klein--Gordon equation in the de Sitter spacetime with structure-preserving scheme

We perform some simulations of the semilinear Klein--Gordon equation in the de Sitter spacetime. We reported the accurate numerical results of the equation with the structure-preserving scheme (SPS) in an earlier publication (Tsuchiya and Nakamura in J. Comput. Appl. Math. \textbf{361}: 396--412, 2019). To investigate the factors for the stability and accuracy of the numerical results with SPS, we perform some simulations with three discretized formulations. The first formulation is the discretized equations with SPS, the second one is with SPS that replaces the second-order difference as the standard second-order central difference, and the third one is with SPS that replaces the discretized nonlinear term as the standard discretized expression. As a result, the above two replacements in SPS are found to be effective for accurate simulations. On the other hand, the ingenuity of replacing the second-order difference in the first formulation is not effective for maintaining the stability of the simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

On the stability of covariant BSSN formulation

In this study, we investigate the numerical stability of the covariant BSSN (cBSSN) formulation proposed by Brown. We calculate the constraint amplification factor (CAF), which is an eigenvalue of the coefficient matrix of the evolution equations of the constraints on the cBSSN formulation and on some adjusted formulations with constraints added to the evolution equations. The adjusted formulations have a higher numerical stability than the cBSSN formulation from the viewpoint of the CAF.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Stable numerical simulation of Einstein equations in gravitational collapse space--time

We perform simulations in a gravitational collapsing model using the Einstein equations. In this paper, we review the equations for constructing the initial values and the evolution form of the Einstein equations called the BSSN formulation. In addition, since we treat a nonvacuum case, we review the evolution equations of the matter fields of a perfect fluid. To make the simulations stable, we propose a modified system, which decreases numerical errors in analysis, and we actually perform stable simulations with decreased numerical errors.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Constraint-Preserving Scheme for Maxwell's Equations

We derive the discretized Maxwell's equations using the discrete variational derivative method (DVDM), calculate the evolution equation of the constraint, and confirm that the equation is satisfied at the discrete level. Numerical simulations showed that the results obtained by the DVDM are superior to those obtained by the Crank-Nicolson scheme. In addition, we study the two types of the discretized Maxwell's equations by the DVDM and conclude that if the evolution equation of the constraint is not conserved at the discrete level, then the numerical results are also unstable.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Post-Newtonian effects on the stability of the triangular solution in the three-body problem for general masses

Continuing work initiated in earlier publications [Ichita, Yamada and Asada, Phys. Rev. D {\bf 83}, 084026 (2011); Yamada and Asada, Phys. Rev. D {\bf 86}, 124029 (2012)], we examine the post-Newtonian (PN) effects on the stability of the triangular solution in the relativistic three-body problem for general masses. For three finite masses, a condition for stability of the triangular solution is obtained at the first post-Newtonian (1PN) order, and it recovers previous results for the PN restricted three-body problem when one mass goes to zero. The stability regions still exist even at the 1PN order, though the PN triangular configuration for general masses is less stable than the PN restricted three-body case as well as the Newtonian one.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the circumradius condition for piecewise linear triangular elements

We discuss the error analysis of linear interpolation on triangular elements. We claim that the circumradius condition is more essential then the well-known maximum angle condition for convergence of the finite element method. Numerical experiments show that this condition is the best possible. We also point out that the circumradius condition is closely related to the definition of surface area.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Constraint Propagation of $C^2$-adjusted Formulation II --- Another Recipe for Robust Baumgarte-Shapiro-Shibata-Nakamura Evolution System ---

In order to obtain an evolution system which is robust against the violation of constraints, we present a new set of evolution systems based on the so-called Baumgarte-Shapiro-Shibata-Nakamura (BSSN) equations.The idea is to add functional derivatives of the norm of constraints, $C^2$, to the evolution equations, which was proposed by Fiske (2004) and was applied to the ADM formulation in our previous study. We derive the constraint propagation equations, discuss the behavior of constraint damping, and present the results of numerical tests using the gauge-wave and polarized Gowdy wave spacetimes. The construction of the $C^2$-adjusted system is straightforward. However, in BSSN, there are two kinetic constraints and three algebraic constraints; thus, the definition of $C^2$ is a matter of concern. By analyzing constraint propagation equations, we conclude that $C^2$ should include all the constraints, which is also confirmed numerically. By tuning the parameters, the lifetime of the simulations can be increased as 2-10 times as longer than those of the standard BSSN evolutions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Constraint Propagation of $C^2$-adjusted Formulation - Another Recipe for Robust ADM Evolution System

With a purpose of constructing a robust evolution system against numerical instability for integrating the Einstein equations, we propose a new formulation by adjusting the ADM evolution equations with constraints. We apply an adjusting method proposed by Fiske (2004) which uses the norm of the constraints, C2. One of the advantages of this method is that the effective signature of adjusted terms (Lagrange multipliers) for constraint-damping evolution is pre-determined. We demonstrate this fact by showing the eigenvalues of constraint propagation equations. We also perform numerical tests of this adjusted evolution system using polarized Gowdy-wave propagation, which show robust evolutions against the violation of the constraints than that of the standard ADM formulation.

0 citations0 reviewsNo code linkedOpen access