質問編集履歴
4
書式の改善
title
CHANGED
|
File without changes
|
body
CHANGED
|
@@ -1,25 +1,4 @@
|
|
|
1
1
|
```ここに言語を入力
|
|
2
|
-
コード
|
|
3
|
-
```### 前提・実現したいこと
|
|
4
|
-
プログラミング未経験の者です。以下リンク、科学技術計算講座3-熱伝導方程式シミュレーション-FTCS法のプログラミングのなかで
|
|
5
|
-
[リンク内容](http://cattech-lab.com/science-tools/simulation-lecture-3-4/)
|
|
6
|
-
指定したある時間幅だけに100℃の熱を加えるモデルを考えたとき、Pythonのプログラムはどのように設定するべきか教えていただきたいです。
|
|
7
|
-
具体的に、以下に示したソースコードでは、初期条件で表面温度（ここでは一次元的に考えてx=0の温度）temp_bc = 100.0で固定、計算期間が20000秒で0秒から0.1秒刻みで計算されています。計算期間をもっと短くしてある指定した時間にだけ表面に100℃の熱を与えたときの温度分布を出したいです。
|
|
8
|
-
|
|
9
|
-
サイトに問い合わせたところ特別に概要のみ以下のように教えて頂けました。
|
|
10
|
-
|
|
11
|
-
> プログラムの時間ループ（# time loop）のなかで、
|
|
12
|
-
> # Boundary condition
|
|
13
|
-
> temp[0] = temp_bc # Dirichlet @ x=0 の所の境界条件をif文を使い時間指定で切り替えればよいと思います。
|
|
14
|
-
|
|
15
|
-
|
|
16
|
-
簡単にできるとのことでしたが、プログラミング経験が全くないためif文について調べてもいまいち時間指定での切り替えまでは理解できませんでした。
|
|
17
|
-
|
|
18
|
-
|
|
19
|
-
### 発生している問題・エラーメッセージ
|
|
20
|
-
|
|
21
|
-
|
|
22
|
-
### 該当のソースコード
|
|
23
2
|
import matplotlib.pyplot as plt
|
|
24
3
|
import matplotlib.animation as animation
|
|
25
4
|
import numpy as np
|
|
@@ -89,8 +68,31 @@
|
|
|
89
68
|
anm = animation.ArtistAnimation(fig, ims, interval=50)
|
|
90
69
|
anm.save('animation.gif', writer='pillow')
|
|
91
70
|
plt.show()
|
|
71
|
+
``````ここに言語を入力
|
|
72
|
+
コード
|
|
73
|
+
```### 前提・実現したいこと
|
|
74
|
+
プログラミング未経験の者です。以下リンク、科学技術計算講座3-熱伝導方程式シミュレーション-FTCS法のプログラミングのなかで
|
|
75
|
+
[リンク内容](http://cattech-lab.com/science-tools/simulation-lecture-3-4/)
|
|
76
|
+
指定したある時間幅だけに100℃の熱を加えるモデルを考えたとき、Pythonのプログラムはどのように設定するべきか教えていただきたいです。
|
|
77
|
+
具体的に、以下に示したソースコードでは、初期条件で表面温度（ここでは一次元的に考えてx=0の温度）temp_bc = 100.0で固定、計算期間が20000秒で0秒から0.1秒刻みで計算されています。計算期間をもっと短くしてある指定した時間にだけ表面に100℃の熱を与えたときの温度分布を出したいです。
|
|
92
78
|
|
|
79
|
+
サイトに問い合わせたところ特別に概要のみ以下のように教えて頂けました。
|
|
93
80
|
|
|
81
|
+
> プログラムの時間ループ（# time loop）のなかで、
|
|
82
|
+
> # Boundary condition
|
|
83
|
+
> temp[0] = temp_bc # Dirichlet @ x=0 の所の境界条件をif文を使い時間指定で切り替えればよいと思います。
|
|
84
|
+
|
|
85
|
+
|
|
86
|
+
簡単にできるとのことでしたが、プログラミング経験が全くないためif文について調べてもいまいち時間指定での切り替えまでは理解できませんでした。
|
|
87
|
+
|
|
88
|
+
|
|
89
|
+
### 発生している問題・エラーメッセージ
|
|
90
|
+
|
|
91
|
+
|
|
92
|
+
### 該当のソースコード
|
|
93
|
+
|
|
94
|
+
|
|
95
|
+
|
|
94
96
|
### 試したこと
|
|
95
97
|
|
|
96
98
|
ここに問題に対して試したことを記載してください。
|
3
書式の改善
title
CHANGED
|
File without changes
|
body
CHANGED
|
@@ -18,9 +18,9 @@
|
|
|
18
18
|
|
|
19
19
|
### 発生している問題・エラーメッセージ
|
|
20
20
|
|
|
21
|
-
|
|
21
|
+
|
|
22
|
-
|
|
22
|
+
### 該当のソースコード
|
|
23
|
-
|
|
23
|
+
import matplotlib.pyplot as plt
|
|
24
24
|
import matplotlib.animation as animation
|
|
25
25
|
import numpy as np
|
|
26
26
|
|
|
@@ -90,12 +90,7 @@
|
|
|
90
90
|
anm.save('animation.gif', writer='pillow')
|
|
91
91
|
plt.show()
|
|
92
92
|
|
|
93
|
-
コード
|
|
94
|
-
```
|
|
95
93
|
|
|
96
|
-
### 該当のソースコード
|
|
97
|
-
|
|
98
|
-
|
|
99
94
|
### 試したこと
|
|
100
95
|
|
|
101
96
|
ここに問題に対して試したことを記載してください。
|
2
書式の改善
title
CHANGED
|
File without changes
|
body
CHANGED
|
@@ -20,10 +20,7 @@
|
|
|
20
20
|
|
|
21
21
|
```
|
|
22
22
|
エラーメッセージ
|
|
23
|
-
```
|
|
24
|
-
|
|
25
|
-
### 該当のソースコード
|
|
26
|
-
```
|
|
23
|
+
``````import matplotlib.pyplot as plt
|
|
27
24
|
import matplotlib.animation as animation
|
|
28
25
|
import numpy as np
|
|
29
26
|
|
|
@@ -93,7 +90,12 @@
|
|
|
93
90
|
anm.save('animation.gif', writer='pillow')
|
|
94
91
|
plt.show()
|
|
95
92
|
|
|
93
|
+
コード
|
|
94
|
+
```
|
|
96
95
|
|
|
96
|
+
### 該当のソースコード
|
|
97
|
+
|
|
98
|
+
|
|
97
99
|
### 試したこと
|
|
98
100
|
|
|
99
101
|
ここに問題に対して試したことを記載してください。
|
1
書式の改善
title
CHANGED
|
File without changes
|
body
CHANGED
|
@@ -1,4 +1,6 @@
|
|
|
1
|
+
```ここに言語を入力
|
|
2
|
+
コード
|
|
1
|
-
### 前提・実現したいこと
|
|
3
|
+
```### 前提・実現したいこと
|
|
2
4
|
プログラミング未経験の者です。以下リンク、科学技術計算講座3-熱伝導方程式シミュレーション-FTCS法のプログラミングのなかで
|
|
3
5
|
[リンク内容](http://cattech-lab.com/science-tools/simulation-lecture-3-4/)
|
|
4
6
|
指定したある時間幅だけに100℃の熱を加えるモデルを考えたとき、Pythonのプログラムはどのように設定するべきか教えていただきたいです。
|
|
@@ -21,7 +23,7 @@
|
|
|
21
23
|
```
|
|
22
24
|
|
|
23
25
|
### 該当のソースコード
|
|
24
|
-
import matplotlib.pyplot as plt
|
|
26
|
+
```こ```import matplotlib.pyplot as plt
|
|
25
27
|
import matplotlib.animation as animation
|
|
26
28
|
import numpy as np
|
|
27
29
|
|
|
@@ -91,6 +93,7 @@
|
|
|
91
93
|
anm.save('animation.gif', writer='pillow')
|
|
92
94
|
plt.show()
|
|
93
95
|
|
|
96
|
+
|
|
94
97
|
### 試したこと
|
|
95
98
|
|
|
96
99
|
ここに問題に対して試したことを記載してください。
|