編集履歴

質問編集履歴

コードを修正

2020/10/19 00:47

投稿

退会済みユーザー

スコア0

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -218,8 +218,6 @@
 def reactor(N, γ, ΔT):
-    boolean=False
     P_log=[[]]
     P_list=[]

コードを修正しました。

2020/10/19 00:47

投稿

退会済みユーザー

スコア0

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -32,7 +32,7 @@
       return z
-print(division(2,0))
+division(2,0)
 ```

コードを修正しました。

2020/10/18 16:03

投稿

退会済みユーザー

スコア0

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -32,7 +32,7 @@
       return z
-print(2,0)
+print(division(2,0))
 ```

解決策を追記しました。

2020/10/18 16:02

投稿

退会済みユーザー

スコア0

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -154,6 +154,8 @@
 P_list: 反応器毎の攪拌消費電力を記録する(record the power consumption of each reactor)
+P_total: 最終的な計算結果を保存する(save the final result of calculation)
 """

解決策を追記しました。

2020/10/18 15:58

投稿

退会済みユーザー

スコア0

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -1,3 +1,5 @@
+returnで値が正しく返ってきません。
 ```Python
 def division(x,y):

解決策を追記しました。

2020/10/18 15:53

投稿

退会済みユーザー

スコア0

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -33,3 +33,279 @@
 print(2,0)
 ```
+多くの方にご回答いただき誠にありがとうございました。おかげさまで書きたかったプログラムが書けました。下のコードのように書くと、引数を更新して正しく計算されるまで関数を呼び直すという動作ができます。
+```Python
+Z=[]
+def division(x,y):
+  boolean=False
+  x=x
+  y=y
+  try:
+    z=x/y
+    boolean=True
+    Z.append(z)
+  except ZeroDivisionError:
+    y+=1
+    division(x,y)
+division(2,0)
+print(Z[0])
+```
+また、これを応用して、反応器段数、原料濃度、溶液と冷媒の温度差を変数にとり連続混合槽型反応器のシミュレーションを行うコードを作成しました。（大学で学んでいる化学工学の問題です。かなりマニアックな領域なので化学工学に馴染みのない方が見るのは難しいかもしれません。）
+```Python
+"""
+与えられた定数(prepared constants)
+Gp: 生産量(product flow rate) [kg/s]
+ζ: 反応転化率(conversion)
+γ: 原料中のブタジエン重量割合(weight fraction of monomer in the feed)
+τ: 滞留時間(residence time) [s]
+ρ: 溶液密度(density of mixture) [kg/m^3]
+α: 反応器の高さ/直径比(ratio of height to diameter of the reactor)
+ΔH: ブタジエンの反応熱(heat of polymerization of butadiene) [J/mol]
+ΔT: 反応溶液と冷媒の温度差(temperature difference between reactant and coolant) [℃]
+λ: 溶液の熱伝導率(thermal conductivity of solution) [W/(m・K)]
+μ: 溶液の粘度(viscosity of solution) [Pa・s]
+Cp: 溶液の比熱(specific heat of solution) [J/(kg・K)]
+M: ブタジエンの分子量(molecular weight of butadiene) [g/mol]
+"""
+"""
+計算により求められる値(values obtained by calculation)
+v: 溶液の体積流量(volumetric flow rate of solution)　[m^3/s]
+V: 反応器体積(volume of the reactor) [m^3]
+D: 反応器直径(diameter of the reactor) [m]
+Qp: 反応熱(heat of reaction) [J]
+P: 攪拌の消費電力(power consumption of stirring) [kW]
+h: 冷却装置の熱伝達率(heat transfer coefficiency of the coolant) [W/m^2・K]
+Nu: ヌッセルト数(Nusselt number)
+Pr: プラントル数(Prantdl number)
+Re: レイノルズ数(Reynols number)
+n: 1分あたり攪拌回数(number of stirring per minute) [/s]
+Np: 動力数(power number)
+"""
+"""
+値の記録に使用するリスト(list to record the value)
+P_log: i段目の反応器の消費電力を計算する際の挙動をi番目の要素に記録する
+       (record the behavior when calculating the power consumption of the i-th reactor in the i-th element)
+P_list: 反応器毎の攪拌消費電力を記録する(record the power consumption of each reactor)
+"""
+"""
+例外処理(exception handling)
+排熱を処理しきれない場合は冷媒の温度を1度ずつ下げる
+(if the exhaust heat cannot be processed, lower the temperature of the refrigerant by 1 degree.)
+"""
+import math
+Gp = 0.729
+ζ = 0.60
+γ_1 = 0.05
+τ_1 = 18000
+ρ = 850
+α = 1.3
+ΔH = 72800
+λ = 0.128
+Cp = 1680
+M = 54.0
+V_1=514.6 #V_1=v*τ_1
+v_1=0.02859 #v_1=Gp/(ρ*γ_1*ζ)
+#反応速度定数を求める(find the reaction rate constant)
+ζ_total=ζ
+k=v_1*ζ_total/(V_1*(1-ζ_total))
+P_log=[[]]
+P_list=[]
+P_total=[]
+def reactor(N, γ, ΔT):
+    boolean=False
+    P_log=[[]]
+    P_list=[]
+    N=N
+    γ=γ
+    ΔT=ΔT
+    #反応器のサイズを求める(find the size of each reactor)
+    τ=(1/k)*(1/math.pow(1-ζ_total,1/N)-1)
+    v=Gp/(ρ*γ*ζ_total)
+    V=v*τ
+    D=math.pow(4*V/(1.3*math.pi),1/3)
+    #各反応器の熱収支を計算する(calculate heat balance of each reactor)
+    ζ_each=1-math.pow(1-ζ_total,1/N)
+    try:
+        for i in range(N):
+            P_log.append([])
+            P_log[i+1].append(0)
+            Qp=(ΔH/(M*0.001))*k*ρ*V*(1-ζ_each)**(i+1)*γ
+            h=Qp/(1.3*math.pi*(D**2)*ΔT)
+            Nu=h*D/λ
+            μ=0.0010*math.pow(50,1.7)*math.pow(1-(1-ζ_each)**(i+1),2.5)*math.exp(21*γ)
+            Pr=μ*Cp/λ
+            Re=math.pow((1/0.5)*Nu*math.pow(Pr,-1/3),3/2)
+            n=μ*Re/(ρ*(D/2)**2)
+            Np=14.6*math.pow(Re,-0.28)
+            P=0.0010*Np*ρ*(n**3)*((D/2)**5)
+            P_log[i+1].append(P)
+            while abs((P-P_log[i+1][len(P_log[i+1])-2])/P)>1.00*math.pow(10, -3):
+                h=(Qp+P*1000)/(1.3*math.pi*(D**2)*ΔT)
+                Nu=h*D/λ
+                Pr=μ*Cp/λ
+                Re=math.pow((1/0.5)*Nu*math.pow(Pr,-1/3),3/2)
+                n=μ*Re/(ρ*(D/2)**2)
+                Np=14.6*math.pow(Re,-0.28)
+                P=0.0010*Np*ρ*(n**3)*((D/2)**5)
+                P_log[i+1].append(P)
+            P_list.append(P_log[i+1][-1])
+            P_total.append(sum(P_list))
+            P_total.append(ΔT)
+    except OverflowError:
+        if ΔT<=50+273:
+            ΔT+=1
+            reactor(N, γ, ΔT)
+        else:
+            print("Warning: The feed is too thick. Reduce γ.")
+reactor(1, 0.07, 65)
+print(P_total)
+```

コードを少し修正しました。

2020/10/18 15:52

投稿

退会済みユーザー

スコア0

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -20,12 +20,6 @@
     division(x,y)
-  else:
-    Z=z
-    pass
   if boolean==True:
       #print(z)は正しく値を返す
@@ -34,6 +28,8 @@
       #return zは値を返さない
-      return Z
+      return z
+print(2,0)
 ```

正しくマークダウン記法を使って書き直しました。また、プログラムを短いものに差し替えました。

2020/10/18 15:00

投稿

退会済みユーザー

スコア0

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -28,11 +28,11 @@
   if boolean==True:
+      #print(z)は正しく値を返す
       print(z)
-      #print(z)は正しく値を返す
-      return zは値を返さない
+      #return zは値を返さない
       return Z

正しくマークダウン記法で書き直しました。また、プログラムの重要な箇所だけ抽出して書き直しました。

2020/10/18 14:52

投稿

退会済みユーザー

スコア0

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -1,189 +1,39 @@
-### 前提・実現したいこと
+```Python
+def division(x,y):
+  boolean=False
-化学工学の課題で書いたプログラムです。
+  x=x
-反応器の攪拌に必要な消費電力Powerを計算したいです。
+  y=y
+  try:
+    z=x/y
-### 発生している問題・エラーメッセージ
+    boolean=True
+  except ZeroDivisionError:
+    y+=1
-printではきちんと値を表示できるのに、returnを使うと何も返ってきません。
+    division(x,y)
+  else:
+    Z=z
+    pass
+  if boolean==True:
-### 該当のソースコード
+      print(z)
+      #print(z)は正しく値を返す
+      return zは値を返さない
-import math
+      return Z
-Gp = 0.729
-ζ = 0.60
-γ_1 = 0.05
-τ_1 = 18000
-ρ = 850
-α = 1.3
-ΔH = 72800
-λ = 0.128
-Cp = 1680
-M = 54.0
-V_1=514.6 #V_1=v*τ_1
-v_1=0.02859 #v_1=Gp/(ρ*γ_1*ζ)
-#反応速度定数を求める(find the reaction rate constant)
-ζ_total=ζ
-k=v_1*ζ_total/(V_1*(1-ζ_total))
-P_log=[[]]
-P_list=[]
-def reactor(N, γ, ΔT):
-    boolean=False
-    P_log=[[]]
-    P_list=[]
-    N=N
+```
-    γ=γ
-    ΔT=ΔT
-    #反応器のサイズを求める(find the size of each reactor)
-    τ=(1/k)*(1/math.pow(1-ζ_total,1/N)-1)
-    v=Gp/(ρ*γ*ζ_total)
-    V=v*τ
-    D=math.pow(4*V/(1.3*math.pi),1/3)
-    #各反応器の熱収支を計算する(calculate heat balance of each reactor)
-    ζ_each=1-math.pow(1-ζ_total,1/N)
-    try:
-        for i in range(N):
-            P_log.append([])
-            P_log[i+1].append(0)
-            Qp=(ΔH/(M*0.001))*k*ρ*V*(1-ζ_each)**(i+1)*γ
-            h=Qp/(1.3*math.pi*(D**2)*ΔT)
-            Nu=h*D/λ
-            μ=0.0010*math.pow(50,1.7)*math.pow(1-(1-ζ_each)**(i+1),2.5)*math.exp(21*γ)
-            Pr=μ*Cp/λ
-            Re=math.pow((1/0.5)*Nu*math.pow(Pr,-1/3),3/2)
-            n=μ*Re/(ρ*(D/2)**2)
-            Np=14.6*math.pow(Re,-0.28)
-            P=0.0010*Np*ρ*(n**3)*((D/2)**5)
-            P_log[i+1].append(P)
-            while abs((P-P_log[i+1][len(P_log[i+1])-2])/P)>1.00*math.pow(10, -3):
-                h=(Qp+P*1000)/(1.3*math.pi*(D**2)*ΔT)
-                Nu=h*D/λ
-                Pr=μ*Cp/λ
-                Re=math.pow((1/0.5)*Nu*math.pow(Pr,-1/3),3/2)
-                n=μ*Re/(ρ*(D/2)**2)
-                Np=14.6*math.pow(Re,-0.28)
-                P=0.0010*Np*ρ*(n**3)*((D/2)**5)
-                P_log[i+1].append(P)
-            P_list.append(P_log[i+1][-1])
-            if len(P_list)==N:
-                boolean=True
-    except OverflowError:
-        if ΔT<=50+273:
-            ΔT+=1
-            reactor(N, γ, ΔT)
-        else:
-            print("Warning: The feed is too thick. Reduce γ.")
-    if boolean==True:
-        Power=P_list[-1]
-        print(Power)
-        return Power
-reactor(1, 0.07, 67)
-### 試したこと
-最後のreactor(1, 0.07, 67)をprint(reactor(1, 0.07, 67))とするとPowerの値ではなく"None"が返ってきました。
-### 補足情報（FW/ツールのバージョンなど）
-jupyter notebookで実行しています。
-printで表示すると確かに値を持っているのに、returnを使うとNoneになってしまう…。このようなことは初めてでとても困惑しています。