地盤災害システム論 前田研究室 maedalab Report 01

課題１： _Report １

表 ：検証 用いた関数 Function used for verification

関数 種類

Kind of Function

元 方程式

Source Equation

概観図 ‘Diffsion1D_000000000.thd’ ^出力 Graph(Output of ‘Diffsion1D_000000000.thd’)

Default Function

f _ sin _{ } 0 _{    }

0 0.2 0.4 0.6 0.8 1

0 0.2 0.4 0.6 0.8 1

x

u

Smoothing Function

(3D B-Spline Function)

^{ }

 





 























S

S

S

S

S

S

f

2

2

1

1

0

0

6 2

1

2

1

3

2

3 3 2

0 0.2 0.4 0.6 0.8 1

0 0.2 0.4 0.6 0.8 1

x

u

Triangle

f _ ax _ b

0 0.2 0.4 0.6 0.8 1

0 0.2 0.4 0.6 0.8 1

x

u

Step Function

f _ 0 . 8

0 0.2 0.4 0.6 0.8 1

0 0.2 0.4 0.6 0.8 1

x

u

表 ：各 テップ け グ フ形状 比較 Comparison of shape in graph in each step

STEP 0 STEP 250 STEP 500 STEP 750 STEP 1000

Default Function

Smoothing Function

Triangle

Step Function

※各関数 振幅を_0.8 計算 あ ま Each function has calculated the amplitude as 0.8.

計算 上手くい なくなる条件

CFL^条件 Courant-Friedrichs-Levy condition

CFL condition^：

 ¹

 ^{ x}

t

v

● 時間差分間隔

_{ t}

_ x

関係式 あ こ 条件 満たさ い 計算 不安定 こ 知 い

今回 課題 初期パ メ タ設定(Initial Condition of Parameter)

WavePropagation1D.idt ^{あ い 課題プ ント} How to use… Page2 ^を参考

 ¹ 2 ^{乗 値} 

0

.

10

0

.

100 _ v _ cwp

 

 ^{tinc 値} 

t _ 0 . 0001

 ^x ₁₀₀ ^{0 . 01}  ^xl ndivx ^値 

0

.

1 





い

1

1

.

01 0

.

0

0001

.

0

0

.

10 ^  

 ^{ x} 

t

v

(2)^式 (1) ^{条件を満た い 安定 た計算結果 得} (Good!)

そ _CFL条件を満たさ た う ？(When CFL condition is Broken, what happen?)

 1

 ^{ x}

t

v

_{ 1}

 ^{ x}

t

v

Default Function^{を用い パ メ タを以下 う 設定 た} (Setting Parameters)

変更 を

_{100 . 0}

₁₀₀₀₀





^{100 . 0}

^{10000 . 1}

計算結果_(Result) 以下 う ※_Fortran 計算時 パ メ タ設定 い 警告メッセ 出

STEP 25 STEP 50 STEP 75 STEP 100

う 波 伝播 こ 分 伝播 様子 明 く い こ 分

一方 波 振動 振幅 異常 増加 た _STEP100 時点 最大値_8.6×₁₀³⁵