質問編集履歴
2
書式の改善
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -66,11 +66,7 @@
|
|
|
66
66
|
|
|
67
67
|
t=dt*n
|
|
68
68
|
|
|
69
|
-
Q=
|
|
69
|
+
Q=40000000000*t*pow(2.718281828,-a*x)
|
|
70
|
-
|
|
71
|
-
|
|
72
|
-
|
|
73
|
-
|
|
74
70
|
|
|
75
71
|
|
|
76
72
|
|
|
@@ -126,13 +122,13 @@
|
|
|
126
122
|
|
|
127
123
|
|
|
128
124
|
|
|
129
|
-
#② FTCS差分式(右辺第一項に吸収熱を設定/0～50nsの間だけ照射、それ以降はQ＝0))
|
|
125
|
+
#② FTCS差分式(右辺第一項に吸収熱Q[i]を設定/0～50nsの間だけ照射、それ以降はQ＝0))
|
|
130
126
|
|
|
131
127
|
for i in range(1, nx-1):
|
|
132
128
|
|
|
133
129
|
if time < t_intv :
|
|
134
130
|
|
|
135
|
-
temp_new[i] = Q+temp[i] + dt * alpha * (temp[i+1] - 2.0 * temp[i] + temp[i-1]) / (dx * dx)
|
|
131
|
+
temp_new[i] = Q[i]+temp[i] + dt * alpha * (temp[i+1] - 2.0 * temp[i] + temp[i-1]) / (dx * dx)
|
|
136
132
|
|
|
137
133
|
else :
|
|
138
134
|
|
1
書式の改善
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -1,4 +1,4 @@
|
|
|
1
|
-
プログラミング未経験の者です。以下リンク、科学技術計算講座3-熱伝導方程式シミュレーションリンク内容を参考にして、厚さ1mmのSi材料にレーザー光を一定時間照射したときの試料の温度分布を求めたいと思っています。
|
|
1
|
+
プログラミング未経験の者です。以下リンク、科学技術計算講座3-熱伝導方程式シミュレーション[リンク内容](https://cattech-lab.com/science-tools/simulation-lecture-3-4/)を参考にして、厚さ1mmのSi材料にレーザー光を一定時間照射したときの試料の温度分布を求めたいと思っています。
|
|
2
2
|
|
|
3
3
|
試料は一次元（深さ）方向のみ考えており位置i、時間nとそれぞれを離散化しています。（照射）時間に比例、深さ方向に対しては指数関数的に減少するレーザー光の吸収熱Qを設定し、一次元非定常熱伝導方程式をもとにFTCS差分式の右辺第1項に足しました。
|
|
4
4
|
|