NAJan 2, 2018
Strong convergence of a fully discrete finite element approximation of the stochastic Cahn-Hilliard equationDaisuke Furihata, Mihály Kovács, Stig Larsson et al.
We consider the stochastic Cahn-Hilliard equation driven by additive Gaussian noise in a convex domain with polygonal boundary in dimension $d\le 3$. We discretize the equation using a standard finite element method in space and a fully implicit backward Euler method in time. By proving optimal error estimates on subsets of the probability space with arbitrarily large probability and uniform-in-time moment bounds we show that the numerical solution converges strongly to the solution as the discretization parameters tend to zero.
NAApr 23, 2017
Geometric numerical integrators for Hunter-Saxton-like equationsYuto Miyatake, David Cohen, Daisuke Furihata et al.
We present novel geometric numerical integrators for Hunter--Saxton-like equations by means of new multi-symplectic formulations and known Hamiltonian structures of the problems. We consider the Hunter--Saxton equation, the modified Hunter--Saxton equation, and the two-component Hunter--Saxton equation. Multi-symplectic discretisations based on these new formulations of the problems are exemplified by means of the explicit Euler box scheme, and Hamiltonian-preserving discretisations are exemplified by means of the discrete variational derivative method. We explain and justify the correct treatment of boundary conditions in a unified manner. This is necessary for a proper numerical implementation of these equations and was never explicitly clarified in the literature before, to the best of our knowledge. Finally, numerical experiments demonstrate the favourable behaviour of the proposed numerical integrators.
NAJul 16, 2016
A novel discrete variational derivative method using "average-difference methods"Daisuke Furihata, Shun Sato, Takayasu Matsuo
We consider structure-preserving methods for conservative systems, which rigorously replicate the conservation property yielding better numerical solutions. There, corresponding to the skew-symmetry of the differential operator, that of difference operators is essential to the discrete conservation law. Unfortunately, however, when we employ the standard central difference operator, the simplest one, the numerical solutions often suffer from undesirable spatial oscillations. In this letter, we propose a novel "average-difference method," which is tougher against such oscillations, and combine it with an existing conservative method. Theoretical and numerical analysis in the linear case show the superiority of the proposed method.
CLJun 8, 2023
A modified model for topic detection from a corpus and a new metric evaluating the understandability of topicsTomoya Kitano, Yuto Miyatake, Daisuke Furihata
This paper presents a modified neural model for topic detection from a corpus and proposes a new metric to evaluate the detected topics. The new model builds upon the embedded topic model incorporating some modifications such as document clustering. Numerical experiments suggest that the new model performs favourably regardless of the document's length. The new metric, which can be computed more efficiently than widely-used metrics such as topic coherence, provides variable information regarding the understandability of the detected topics.
NAJan 27, 2016
Pathwise existence of solutions to the Implicit Euler method for the stochastic Cahn-Hilliard EquationDaisuke Furihata, Fredrik Lindgren, Shuji Yoshikawa
We consider the implicit Euler approximation of the stochastic Cahn-Hilliard equation driven by additive Gaussian noise in a spatial domain with smooth boundary in dimension $d\le 3$. We show pathwise existence and uniqueness of solutions for the method under a restriction on the step size that is independent of the size of the initial value and of the increments of the Wiener process. This result also relaxes the imposed assumption on the time step for the deterministic Cahn-Hilliard equation assumed in earlier existence proofs.
LGJan 10, 2024
Structure-Preserving Physics-Informed Neural Networks With Energy or Lyapunov StructureHaoyu Chu, Yuto Miyatake, Wenjun Cui et al.
Recently, there has been growing interest in using physics-informed neural networks (PINNs) to solve differential equations. However, the preservation of structure, such as energy and stability, in a suitable manner has yet to be established. This limitation could be a potential reason why the learning process for PINNs is not always efficient and the numerical results may suggest nonphysical behavior. Besides, there is little research on their applications on downstream tasks. To address these issues, we propose structure-preserving PINNs to improve their performance and broaden their applications for downstream tasks. Firstly, by leveraging prior knowledge about the physical system, a structure-preserving loss function is designed to assist the PINN in learning the underlying structure. Secondly, a framework that utilizes structure-preserving PINN for robust image recognition is proposed. Here, preserving the Lyapunov structure of the underlying system ensures the stability of the system. Experimental results demonstrate that the proposed method improves the numerical accuracy of PINNs for partial differential equations. Furthermore, the robustness of the model against adversarial perturbations in image data is enhanced.
LGJan 29, 2022
Composing a surrogate observation operator for sequential data assimilationKosuke Akita, Yuto Miyatake, Daisuke Furihata
In data assimilation, state estimation is not straightforward when the observation operator is unknown. This study proposes a method for composing a surrogate operator when the true operator is unknown. A neural network is used to improve the surrogate model iteratively to decrease the difference between the observations and the results of the surrogate model. A twin experiment suggests that the proposed method outperforms approaches that tentatively use a specific operator throughout the data assimilation process.