Source author record

Shin-ichi Ito

Shin-ichi Ito appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

cond-mat.stat-mech cond-mat.other cond-mat.soft math.NA Methodology Numerical Analysis physics.comp-ph physics.data-an

Catalog footprint

What is connected

5works

8topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

preprint2020arXiv

A Phase Prediction Method for Pattern Formation in Time-Dependent Ginzburg-Landau Dynamics for Kinetic Ising Model without a priori Assumptions on Domain Patterns

We propose a phase prediction method for the pattern formation in the uniaxial two-dimensional kinetic Ising model with the dipole-dipole interactions under the time-dependent Ginzburg-Landau dynamics. Taking the effects of the material thickness into account by assuming the uniformness along the magnetization axis, the model corresponds to thin magnetic materials with long-range repulsive interactions. We propose a new theoretical basis to understand the effects of the material parameters on the formation of the magnetic domain patterns in terms of the equation of balance governing the balance between the linear- and nonlinear forces in the equilibrium state. Based on this theoretical basis, we propose a new method to predict the phase in the equilibrium state reached after the time-evolution under the dynamics with a given set of parameters, by approximating the third-order term using the restricted phase-space approximation [R. Anzaki, K. Fukushima, Y. Hidaka, and T. Oka, Ann. Phys. 353, 107 (2015)] for the $ϕ^4$-models. Although the proposed method does not have the perfect concordance with the actual numerical results, it has no arbitrary parameters and functions to tune the prediction. In other words, it is a method with no a priori assumptions on domain patterns.

preprint2020arXiv

Adjoint-based exact Hessian computation

We consider a scalar function depending on a numerical solution of an initial value problem, and its second-derivative (Hessian) matrix for the initial value. The need to extract the information of the Hessian or to solve a linear system having the Hessian as a coefficient matrix arises in many research fields such as optimization, Bayesian estimation, and uncertainty quantification. From the perspective of memory efficiency, these tasks often employ a Krylov subspace method that does not need to hold the Hessian matrix explicitly and only requires computing the multiplication of the Hessian and a given vector. One of the ways to obtain an approximation of such Hessian-vector multiplication is to integrate the so-called second-order adjoint system numerically. However, the error in the approximation could be significant even if the numerical integration to the second-order adjoint system is sufficiently accurate. This paper presents a novel algorithm that computes the intended Hessian-vector multiplication exactly and efficiently. For this aim, we give a new concise derivation of the second-order adjoint system and show that the intended multiplication can be computed exactly by applying a particular numerical method to the second-order adjoint system. In the discussion, symplectic partitioned Runge--Kutta methods play an essential role.

preprint2016arXiv

Data assimilation for massive autonomous systems based on second-order adjoint method

Data assimilation (DA) is a fundamental computational technique that integrates numerical simulation models and observation data on the basis of Bayesian statistics. Originally developed for meteorology, especially weather forecasting, DA is now an accepted technique in various scientific fields. One key issue that remains controversial is the implementation of DA in massive simulation models under limited computation time and resources. In this paper, we propose an adjoint-based DA method for massive autonomous models that produces optimum estimates and their uncertainties within practical computation time and resource constraints. The uncertainties are given as several diagonal components of an inverse Hessian matrix, which is the covariance matrix of a normal distribution that approximates the target posterior probability density function in the neighborhood of the optimum. Conventional algorithms for deriving the inverse Hessian matrix require $O(CN^2+N^3)$ computations and $O(N^2)$ memory, where $N$ is the number of degrees of freedom of a given autonomous system and $C$ is the number of computations needed to simulate time series of suitable length. The proposed method using a second-order adjoint method allows us to directly evaluate the diagonal components of the inverse Hessian matrix without computing all of its components. This drastically reduces the number of computations to $O(C)$ and the amount of memory to $O(N)$ for each diagonal component. The proposed method is validated through numerical tests using a massive two-dimensional Kobayashi's phase-field model. We confirm that the proposed method correctly reproduces the parameter and initial state assumed in advance, and successfully evaluates the uncertainty of the parameter. Such information regarding uncertainty is valuable, as it can be used to optimize the design of experiments.

preprint2014arXiv

Stochastic modeling on fragmentation process over lifetime and its dynamical scaling law of fragment distribution

We propose a stochastic model of a fragmentation process, developed by taking into account fragment lifetime as a function of their size based on the Gibrat process. If lifetime is determined by a power function of fragment size, numerical results indicate that size distributions at different times can be collapsed into a single time-invariant curve by scaling size by average fragment size (i.e., the distribution obeys the dynamical scaling law). If lifetime is determined by a logarithmic function of fragment size, the distribution does not obey the scaling law. The necessary and sufficient condition that the scaling law is obeyed is obtained by a scaling analysis of the master equation.

preprint2012arXiv

Dynamical scaling of fragment distribution in drying paste

We reproduce patterns of drying paste by means of smoothed particle hydrodynamics which is the one of methods for solving the equations of continuum in the Lagrangian description. In addition to reproduce a realistic pattern, we find that average size of fragments decays in proportion to inverse time in the case of a linear drying process. Distributions of the size of the fragments are obtained depending on the time. We find a universal scaling distribution by scaling analysis with the average size of the fragment.