• Search Research Projects
  • Search Researchers
  • How to Use
  1. Back to previous page

微分代数方程式に対する高速な構造保存数値解法の構築

Research Project

Project/Area Number 22K13955
Research Category

Grant-in-Aid for Early-Career Scientists

Allocation TypeMulti-year Fund
Review Section Basic Section 12040:Applied mathematics and statistics-related
Research InstitutionThe University of Tokyo

Principal Investigator

佐藤 峻  東京大学, 大学院情報理工学系研究科, 助教 (40849072)

Project Period (FY) 2022-04-01 – 2027-03-31
Project Status Granted (Fiscal Year 2023)
Budget Amount *help
¥4,550,000 (Direct Cost: ¥3,500,000、Indirect Cost: ¥1,050,000)
Fiscal Year 2026: ¥910,000 (Direct Cost: ¥700,000、Indirect Cost: ¥210,000)
Fiscal Year 2025: ¥910,000 (Direct Cost: ¥700,000、Indirect Cost: ¥210,000)
Fiscal Year 2024: ¥910,000 (Direct Cost: ¥700,000、Indirect Cost: ¥210,000)
Fiscal Year 2023: ¥910,000 (Direct Cost: ¥700,000、Indirect Cost: ¥210,000)
Fiscal Year 2022: ¥910,000 (Direct Cost: ¥700,000、Indirect Cost: ¥210,000)
Keywords構造保存数値解法 / exponential integrator / SAV法 / 2次保存量 / 陰的線形スキーム / 常微分方程式 / 微分代数方程式 / 計算量
Outline of Research at the Start

微分方程式の数値解法は現代科学の様々な分野において重要な役割を担っている.中でも,数値的に解きづらい問題に対しては,微分方程式の構造 (保存量や対称性など) を尊重した構造保存数値解法が有効である.常微分方程式に対する構造保存数値解法の理論は概ね良く整備されており,近年は計算量を小さくするための工夫が盛んに研究されている.一方で,常微分方程式の一般化である微分代数方程式に対する構造保存数値解法の研究は未だ限定的かつ散発的である.
本研究では,微分代数方程式に対する構造保存数値解法の枠組を整備し,特に近年発展している計算量の小さい構造保存数値解法を微分代数方程式にまで拡張することを目指す.

Outline of Annual Research Achievements

微分方程式の数値解法は現代科学の様々な分野において重要な役割を担っている.中でも,数値的に解きづらい問題に対しては,微分方程式の構造 (保存量や対称性など) を尊重した構造保存数値解法が有効である.常微分方程式に対する構造保存数値解法の理論は概ね良く整備されており,近年は計算量を小さくするための工夫が盛んに研究されている.一方で,常微分方程式の一般化である微分代数方程式に対する構造保存数値解法の研究は未だ限定的かつ散発的である.本研究では,微分代数方程式に対する構造保存数値解法の枠組を整備し,特に近年発展している計算量の小さい構造保存数値解法を微分代数方程式にまで拡張することを目指す.

本研究では,まず,常微分方程式に対して,計算量の小さい構造保存数値解法の研究を行なった.前年度に検討した,2次の保存量を保存する高精度かつ計算量の小さい構造保存数値解法,Scalar Auxiliary Variable (SAV) 法,保存的exponential Runge--Kutta法について,さらに研究を進めた.SAV法は,主に偏微分方程式に対して近年盛んに研究されている手法で,補助変数を導入することで計算量の小さい構造保存数値解法を構成する手法である.また,exponential Runge--Kutta法は,線形項をもつ発展方程式に対して,行列指数関数を用いて構成するRunge--Kutta法の変種で,線形項が支配的な場合に有効であることが知られている.本研究では,これらを巧妙に組み合わせることで,高精度でありながら計算量の小さい構造保存数値解法が構成できることを確認した.今年度は,この手法を非線形Klein--Gordon方程式に適用した際に,他の保存的PDEの場合よりも計算効率がよくなることを発見した.また,この結果を論文にまとめ,採択された.

Current Status of Research Progress
Current Status of Research Progress

2: Research has progressed on the whole more than it was originally planned.

Reason

今年度は,常微分方程式に関する研究が進展した.現段階では微分代数方程式に対して適用はしていないが,その素地は整ったと言える.

Strategy for Future Research Activity

今年度に行った常微分方程式に対する研究のさらなる発展を目指す.また,それとともに,微分代数方程式への適用も目指す.

Report

(2 results)
  • 2023 Research-status Report
  • 2022 Research-status Report
  • Research Products

    (22 results)

All 2024 2023 2022

All Journal Article (5 results) (of which Int'l Joint Research: 1 results,  Peer Reviewed: 5 results,  Open Access: 3 results) Presentation (17 results) (of which Int'l Joint Research: 3 results,  Invited: 3 results)

  • [Journal Article] Existence results on Lagrange multiplier approach for gradient flows and application to optimization2024

    • Author(s)
      Kenya Onuma, Shun Sato
    • Journal Title

      Japan Journal of Industrial and Applied Mathematics

      Volume: 41 Issue: 1 Pages: 165-189

    • DOI

      10.1007/s13160-023-00595-6

    • Related Report
      2023 Research-status Report
    • Peer Reviewed / Open Access
  • [Journal Article] A unified discretization framework for differential equation approach with Lyapunov arguments for convex optimization2023

    • Author(s)
      Kansei Ushiyama, Shun Sato, Takayasu Matsuo
    • Journal Title

      Advances in Neural Information Processing Systems

      Volume: 36

    • Related Report
      2023 Research-status Report
    • Peer Reviewed / Open Access
  • [Journal Article] High-order linearly implicit schemes conserving quadratic invariants2023

    • Author(s)
      Sato Shun、Miyatake Yuto、Butcher John C.
    • Journal Title

      Applied Numerical Mathematics

      Volume: 187 Pages: 71-88

    • DOI

      10.1016/j.apnum.2023.02.005

    • Related Report
      2022 Research-status Report
    • Peer Reviewed / Open Access / Int'l Joint Research
  • [Journal Article] Mathematical analysis of a conservative numerical scheme for the Ostrovsky equation2022

    • Author(s)
      Shuto Kawai, Shun Sato, Takayasu Matsuo
    • Journal Title

      JSIAM Letters

      Volume: 14 Issue: 0 Pages: 53-56

    • DOI

      10.14495/jsiaml.14.53

    • ISSN
      1883-0609, 1883-0617
    • Related Report
      2022 Research-status Report
    • Peer Reviewed
  • [Journal Article] Essential convergence rate of ordinary differential equations appearing in optimization2022

    • Author(s)
      Kansei Ushiyama, Shun Sato, Takayasu Matsuo
    • Journal Title

      JSIAM Letters

      Volume: 14 Issue: 0 Pages: 119-122

    • DOI

      10.14495/jsiaml.14.119

    • ISSN
      1883-0609, 1883-0617
    • Related Report
      2022 Research-status Report
    • Peer Reviewed
  • [Presentation] 最適化手法の連続時間モデルに対する新しい収束率解析法2024

    • Author(s)
      牛山 寛生, 佐藤 峻, 松尾 宇泰
    • Organizer
      日本応用数理学会 第20回 研究部会連合発表会
    • Related Report
      2023 Research-status Report
  • [Presentation] 光錐座標上のsine-Gordon方程式の構造保存的2段スキームの数学的解析2024

    • Author(s)
      川合 秀人, 佐藤 峻
    • Organizer
      日本応用数理学会 第20回 研究部会連合発表会
    • Related Report
      2023 Research-status Report
  • [Presentation] 部分和分作用素を用いた離散変分導関数法2024

    • Author(s)
      梅津 光汰, 佐藤 峻, 松尾 宇泰
    • Organizer
      日本応用数理学会 第20回 研究部会連合発表会
    • Related Report
      2023 Research-status Report
  • [Presentation] Convergence rates of optimization methods in continuous and discrete time2024

    • Author(s)
      Shun Sato
    • Organizer
      International Conference on Scientific Computing and Machine Learning 2024
    • Related Report
      2023 Research-status Report
    • Int'l Joint Research / Invited
  • [Presentation] 最適化のための Runge--Kutta--Chebyshev 法における曲線探索法2023

    • Author(s)
      牛山 寛生, 佐藤 峻, 松尾 宇泰
    • Organizer
      第49回 数値解析シンポジウム
    • Related Report
      2023 Research-status Report
  • [Presentation] KdV 型発展方程式に対するエネルギー保存解法の数学的解析手法について2023

    • Author(s)
      川合 秀人, 松尾 宇泰, 佐藤 峻
    • Organizer
      第49回 数値解析シンポジウム
    • Related Report
      2023 Research-status Report
  • [Presentation] アフィン座標変換に対して不変な射影を用いた常微分方程式の保存的数値解法2023

    • Author(s)
      石井 直樹, 佐藤 峻, 松尾 宇泰
    • Organizer
      第49回 数値解析シンポジウム
    • Related Report
      2023 Research-status Report
  • [Presentation] 発展方程式に対する乱択化分解解法2023

    • Author(s)
      麻生 豊大, 佐藤 峻, 松尾 宇泰
    • Organizer
      第49回 数値解析シンポジウム
    • Related Report
      2023 Research-status Report
  • [Presentation] 連続最適化手法とその連続極限における収束レートの対応について2023

    • Author(s)
      佐藤 峻, 牛山 寛生, 野沢 諒太, 松尾 宇泰
    • Organizer
      RIMS共同研究(公開型)「新時代における高性能科学技術計算法の探究」
    • Related Report
      2023 Research-status Report
    • Invited
  • [Presentation] A Unified Discretization Framework for Differential Equation Approach with Lyapunov Arguments for Convex Optimization2023

    • Author(s)
      牛山 寛生, 佐藤 峻, 松尾 宇泰
    • Organizer
      第26回情報論的学習理論ワークショップ
    • Related Report
      2023 Research-status Report
  • [Presentation] High-order linearly implicit schemes conserving quadratic invariants2023

    • Author(s)
      S. Sato
    • Organizer
      ANODE2023
    • Related Report
      2022 Research-status Report
    • Int'l Joint Research
  • [Presentation] 最適化手法由来のヘッセ行列を伴う連続力学系モデルに対する数値解析学的アプローチ2023

    • Author(s)
      上島 智哉, 佐藤 峻, 牛山 寛生, 松尾 宇泰, 田中 健一郎
    • Organizer
      日本応用数理学会2023年研究部会連合発表会
    • Related Report
      2022 Research-status Report
  • [Presentation] Nesterovの加速勾配法の可変刻み線形多段法としての解釈とその応用について2022

    • Author(s)
      野沢 諒太, 松尾 宇泰, 佐藤 峻
    • Organizer
      日本応用数理学会 2022年度 年会
    • Related Report
      2022 Research-status Report
  • [Presentation] 勾配流に対する離散勾配を用いた最適化手法の統一的記述について2022

    • Author(s)
      牛山 寛生, 佐藤 峻, 松尾 宇泰
    • Organizer
      日本応用数理学会 2022年度 年会
    • Related Report
      2022 Research-status Report
  • [Presentation] 連続最適化に対する数値解析学的アプローチ2022

    • Author(s)
      佐藤 峻, 牛山 寛生, 松尾 宇泰
    • Organizer
      RIMS共同研究 (公開型)「数値解析が拓く次世代情報社会~エッジから富岳まで~」
    • Related Report
      2022 Research-status Report
  • [Presentation] High-order linearly implicit schemes conserving quadratic invariants2022

    • Author(s)
      S. Sato
    • Organizer
      JSPS Seminar: Topics in computational methods for stochastic and deterministic differential equations
    • Related Report
      2022 Research-status Report
    • Int'l Joint Research / Invited
  • [Presentation] 最適化手法記述のための弱い離散勾配について2022

    • Author(s)
      牛山 寛生, 佐藤 峻, 松尾 宇泰
    • Organizer
      2022年度応用数学合同研究集会
    • Related Report
      2022 Research-status Report

URL: 

Published: 2022-04-19   Modified: 2024-12-25  

Information User Guide FAQ News Terms of Use Attribution of KAKENHI

Powered by NII kakenhi