Workshop with the Kanazawa University group 2024

The main goal of this workshop is to deepen the collaboration between the mathematicians from the Institute of Mathematics, CAS, and the group of mathematicians from the Kanazawa university in Japan in the field of differential equations. The emphasis is laid upon both analysis of differential equations as well as their numerical solutions.

Time and Venue

The workshop takes place on Tuesday, 26th November since 9 AM at IM CAS, library, Žitná 25, Prague.

Schedule

  • 09:00 - 09:30 Hirofumi Notsu (Kanazawa University): An energy estimate and a stabilized Lagrange-Galerkin scheme for a multiphase flow model
    Abstract: Multiphase flow models are commonly employed for understanding complex fluid flows. In this talk, we prove an energy decay property for a multiphase flow model in [Gidaspow, Multiphase Flow and Fluidization, 1994] and propose a stabilized Lagrange– Galerkin scheme. Here, a hyperbolic tangent transformation is employed to preserve the boundedness of the volume fraction, and a novel artificial term is added to obtain the stability property. Two-dimensional numerical examples are presented to see the experimental order of convergence and applicability in modeling sedimentation phenomena. This talk is based on the paper [A. Rudiawan, A. Zak, M. Benes, M. Kimura, and H. Notsu. An energy estimate and a stabilized Lagrange-Galerkin scheme for a multiphase flow model. Applied Mathematics Letters, 153:109059, 2024].
  • 09:30 - 10:00 Masato Kimura (Kanazawa University): Well-posedness of a moving boundary problem related to an inverse shape determination problem
    Abstract: The purpose of this study is twofold. First, we revisit a shape optimization reformulation of a shape inverse problem and propose a simple and efficient numerical approach for realizing the minimization problem. Second, we examine the short-time existence and uniqueness of a classical solution to a Hele-Shaw-like system derived from the shape optimization formulation of the shape inverse problem. This is a joint work with Dr. Julius Fergy T. Rabago, Kanazawa University.
  • 10:00 - 10:30 Coffee break
  • 10:30 - 11:00 Koya Sakakibara (Kanazawa University): Advanced numerical analysis of the Plateau problem using the method of fundamental solutions
    Abstract: This presentation introduces a novel high-speed, high-accuracy computational approach for numerically solving the Plateau problem to find minimal surfaces bounded by a given contour. The proposed method employs the Method of Fundamental Solutions (MFS), a mesh-free approach that constructs smooth approximations to minimal surfaces without requiring traditional meshing techniques. Unlike conventional methods, our MFS-based framework enables the generation of continuous surface approximations, facilitating precise convergence analysis for Dirichlet energy and $L^\infty$ error evaluation in mean curvature. This methodology not only assures reliable surface approximation but also enables effective exploration of multiple minimal surfaces sharing the same boundary by employing varied initial conditions. Through this approach, we examine the dependence of solution multiplicity on boundary geometry, presenting a comprehensive analysis of the relationship between boundary contours and corresponding minimal surfaces.
  • 11:00 - 11:30 Petr Pauš (Czech Technical University): Numerical simulation of dislocation multiple cross-slip
    Abstract: Our contribution deals with the phenomenon in material science called multiple cross-slip of dislocations in slip planes. The numerical model is based on a mean curvature flow equation with additional forcing terms included. The curve motion in 3D space is treated using our tilting method, i.e., mapping of a 3D situation onto a single plane where the curve motion is computed. The physical forces acting on a dislocation curve are evaluated in the 3D setting. This is a joint work with Michal Beneš and Miroslav Kolář.
  • 11:30 - 12:00 John Sebastian Simon (Johann Radon Institute for Computation and Applied Mathematics): Robust shape optimization for missing Dirichlet data
    Abstract: A shape optimization problem subject to an elliptic equation that takes into account a missing data on the Dirichlet boundary condition is considered. The shape optimization problem is formulated by looking for the optimal deformation field that varies the domain where the Poisson equation is posed. The optimization problem is then recasted into the so- called low-regret and no-regret problems to obtain deformation fields which are robust in terms of the Dirichlet data. The analysis of the optimization problem were carried out, this includes providing appropriate topology for establishing the continuity of the deformation-to- state map and providing a proof establishing the existence of optimal deformation fields for both low-regret and no-regret problems. Due to the fact that the low-regret problem is a regularized version of the no-regret problem, convergence of the solutions to the former problem towards a solution of the latter was also provided. The first-order Gâteaux derivative of the reformulated objective function for the low-regret problem was also computed providing us a first-order necessary condition and a gradient-method based on Barzilai– Borwein to solve the problem numerically. Numerical examples are provided to illustrate the convergence previously mentioned, and to illustrate the robustness of the low-regret deformation fields with respect to the missing data.
  • 12:00 - 14:00 Lunch break
  • 14:00 - 14:30 Matteo Caggio (Institute of Mathematics, CAS): Vanishing viscosity limit for compressible fluids
    Abstract: We will present recent results and work in progress that concern the vanishing viscosity limit for compressible fluids. A density dependent viscosity model for a barotropic fluid and an heat-conductive model for a quasi-perfect gas will be considered.
  • 14:30 - 15:00 Ivan Gudoshvnikov (Institute of Mathematcs, CAS): Regularity lost: why functions are not enough in continuum plasticity (so that we need measures).
    Abstract: The connection between plasticity phenomenon in mechanics and the sweeping process was noticed since the very discovery of the sweeping process as a mathematical problem. Using the ideas of J.-J. Moreau, we will discuss how the sweeping process can model the evolution of stresses in an elastoplastic medium. However, with the seemingly simplest constitutive law of ``perfect plasticty'' it is, in generall, impossible to find the strain rate as an $L_2$-function. We will demostrate this ill-posedness using Moreau's framework and also show that various constraint qualifications distinguish this case from the well-posed cases (spatially discrete models, plasticty with hardening).
  • 15:00 - 15:30 Coffee break
  • 15:30 - 15:40 Yuki Karasawa (Kanazawa University): A finite element scheme for a viscoelastic fluid flow
    Abstract: Resin products, which are numerous in our surroundings, are increasingly used in automobiles and airplanes due to their low density and lightweight. It is known that when melted, resins as raw materials behave as viscoelastic fluids, which are different from ordinary viscous fluids. Therefore, in numerical simulations related to polymer processing, it is desired to treat resins as viscoelastic fluids. However, numerical simulation methods for viscoelastic fluids are still under development. Viscoelastic fluid models exist called Oldroyd- B model, PTT, Giesekus, and FENE-P. All of these models include a characteristic term called the upper convection derivative, and the discretization of this derivative is important. A second-order accuracy approximation of the upper convection derivative in time based on the Lie derivative has been proposed and implemented in the finite difference framework [D.D. Medeiros, H. Notsu and C.M. Oishi. SIAM J. Numer. Anal., 59, 2955-2988, 2021]. In this talk, we propose a finite element scheme of second-order in time.
  • 15:40 - 15:50 Ryuhei Wakida (Kanazawa University): A phase field model for fault rupture with dynamic elasticity equation
    Abstract: In Japan, earthquakes often occur, and they sometimes cause heavy damage. An earthquake can be roughly divided into two types, fault-type and trench-type. The purpose of our research is to reproduce a fault-type earthquake numerically. This study is based on a phase field crack growth model, and we expand the model to a wave equation. Furthermore, we will try to artificially introduce a friction term that acts on the fault plane into the phase field model. Adding a friction term to the fracture phase field model is challenging. For the finite element simulation, we use FreeFEM software. In addition, we introduce the parallel computing technique on FreeFEM, applying the overlapping Schwarz domain decomposition method and RAS preconditioner.
  • 15:50 - 16:00 Hisanori Miyata (Kanazawa University): Analysis of mass conservation and gelation in Smoluchowski coagulation equation
    Abstract: This is a presentation of my research on the Smoluchowski coagulation equation. The equation is sometimes observed to show a gelation phenomenon, where mass is no longer conserved at a point in time, and the conditions for this will be discussed in this research. The first of the main theorems shows in mathematical terms, using the method of moments, the conditions under which mass is conserved in weak solutions of the equation. The other is conversely, showing the conditions for gelation.
  • 16:00 - 16:10 Takeshi Nishikawa (Kanazawa University): Eigenvalue problems of the biharmonic equation with mixed boundary conditions
    Abstract: It is known that serious accidents can occur due to the eigen modes of an object, such as the collapse of the Tacoma Bridge in the United States in 1940. Therefore, the purpose of this study is to investigate the eigen modes of an object. For simplicity, we investigated the eigen vibrations of a plate in a two dimensional domain. This problem can be reduced to an eigenvalue problem of the biharmonic equation with a mixed boundary condition. To achieve this goal, we proposed a stable numerical computation method using the finite element method for the biharmonic equation. As a result, we were able to investigate how the plate vibrates through numerical calculation.
  • 16:10 - 16:30 Final discussion
  • 18:00 Dinner

Organizer

Šárka Nečasová