• 検索結果がありません。

発表資料2 Mathematica研究会 seminar2011

N/A
N/A
Protected

Academic year: 2018

シェア "発表資料2 Mathematica研究会 seminar2011"

Copied!
18
0
0

読み込み中.... (全文を見る)

全文

(1)

12/07/11 14:11 Untitled

1. NDSolve で解く微分方程式の例

1.1 エルニーニョのシンプルモデル:Battisti-Hirst 方程式

 エルニーニョは、気候の経年変動(周期が1年より長い変動)の代表的な例として気象学でよく知られています。ペルー沖の海洋の温度が平年よ

り高い現象をエルニーニョ(El Nino)、低い現象をラニーニャ(La Nina)と呼びます。

(2)
(3)

12/07/11 14:11 Untitled

Fortran,DKLAG6

(4)

Fortran,DKLAG6

(5)

12/07/11 14:11 Untitled

参考文献:M. Umeki (2011) , Chaos in the Battisti-Hirst original model for El Nino southern oscillation, to appear in Theoretical and Applied Mechanics Japan

2011.

1.2 ハミルトン勾配系

 ハミルトン系を考える場合、その時間発展においてハミルトニアンは時間的に一定です。

(6)

ハミルトン系の数値積分がハミルトニアンを保存するかどうかは、自明ではなく、粗い計算では保存しません。変数分離系の場合には、シンプレク

ティック積分法が有用であることが知られています。一般の場合に、次のようなハミルトン勾配系を考えることで、保存量を発散しないように保つ

事ができます。

(7)

12/07/11 14:11 Untitled

単振動を例にします。以下は単振動のハミルトン勾配系をオイラー法で解いたものです。

(8)
(9)

12/07/11 14:11 Untitled

保存量が2個(HとI)がある場合は、Gを次のように取って両方を同時に収束させることができます。学会発表のようになってしまうので、詳細を知

りたい人は以下の文献を読んでください。

参考文献:M. Umeki (2010) , A Hamiltonian-Gradient System for Multiple Conservation Laws, Theoretical and Applied Mechanics Japan (2011) 59, 195-201.

(10)
(11)

12/07/11 14:11 Untitled

(12)
(13)

12/07/11 14:11 Untitled

(14)
(15)

12/07/11 14:11 Untitled

1.3 点渦系に対するHG系

 点渦系とは、2次元平面内の点の運動です。それぞれの点渦が、他の点渦からの距離の逆数に比例する速度の和で誘導される運動をします。複素

数z=x+iyで表すのが便利です。無限平面での同一の点渦の場合、以下の式で表されます。

dz_j/dt = !_k i/ (z^*_j-z^*_k) , j=1,...,N, z = x + I y, (x_j,y_j)->(x_{2j-1},x_2j)

(16)
(17)

12/07/11 14:11 Untitled

(18)

Created with Wolfram Mathematica 8.0

参照

関連したドキュメント

In 2003, Agiza and Elsadany 7 studied the duopoly game model based on heterogeneous expectations, that is, one player applied naive expectation rule and the other used

Begehr, Zhang Zhongxiang, Du Jinyuan, On Cauchy-Pompeiu for- mula for functions with values in a universal Clifford algebra, to appear in Mathematica Acta Scientia, 2003, 23B(1)..

A NOTE ON SUMS OF POWERS WHICH HAVE A FIXED NUMBER OF PRIME FACTORS.. RAFAEL JAKIMCZUK D EPARTMENT OF

For i= 1, 2 or 3, Models (Mi), subject to Assumptions (A1–5), (Bi) and Remark 2 with regular initial conditions converge to the Keller–Segel model (1) in their drift-diffusion

In this paper, first we give a theorem which generalizes the Banach contraction principle and fixed point theorems given by many authors, and then a fixed point theorem for

A lemma of considerable generality is proved from which one can obtain inequali- ties of Popoviciu’s type involving norms in a Banach space and Gram determinants.. Key words

In section 2 we present the model in its original form and establish an equivalent formulation using boundary integrals. This is then used to devise a semi-implicit algorithm

[Mag3] , Painlev´ e-type differential equations for the recurrence coefficients of semi- classical orthogonal polynomials, J. Zaslavsky , Asymptotic expansions of ratios of