前の一連の記事のプルーム上昇式について書いていた1983年８月の覚書きノートには続きがあります。それは水滴の成長と蒸発という項目です。
これも1000ページ近い"英文のANLモデル(アルゴンヌ国立研究所模型;後のSACTI)のドキュメントを全部翻訳して内容を完全に把握する"という当時のジョブの一環として勉強したものです。
冷却塔(cooling tower)からの温排気排出の影響には,水蒸気による白煙発生の他に,それが周辺大気で冷却され凝結した水滴の飛散による濃霧,また冷却水が海水のように塩分などを含む場合には飛散水滴による農作物等への塩害etc.も考えられるからです。
そこで,この項目についても覚書きとして記述しておきます。
(当時は,仕事の都合上,タイムリミットもあって,とにかく当座の計算アルゴリズム等に必要な式に到達するのが目的だったので,途中のプラス,マイナスなど細かいところはアバウトだったりテキトーだったりなところもあるので,あとで修正が必要なのです。)
以下は本文です。
§1.拡散流束による水滴質量の増減
水蒸気の密度をρv,その密度流束(局所流束密度)をｊvとします。
この局所流束密度ｊvは,移流風速(水蒸気を含む空気の平均流速)をｕとすると,そのゆらぎ(fluctuation)に相当する拡散流束を加えて,
ｊv＝ρvｕ－Ｄv∇ρv と表現されます。
ただし,Ｄvは拡散係数です。
この流束を用いた水蒸気の質量保存を示す連続の方程式は,
∂ρv/∂ｔ＋∇ｊv＝0 or ∂ρv/∂ｔ＝－∇ｊvです。
そして,ｊv＝ρvｕ－Ｄv∇ρvにより∇ｊv＝∇(ρvｕ－Ｄv∇ρv)です。
しかし,この水蒸気を含む空気は非圧縮(または,高々Boussinesq近似)を仮定しているため,∇ｕ＝0 を満たしますから,
これは,∇ｊv＝ｕ∇ρv－∇(Ｄv∇ρv) と書けます。
さらに,拡散係数Ｄvは空間位置によらない定数であると仮定すると,
∇ｊv＝ｕ∇ρv－Ｄv∇2ρvです。
結局,上記連続の方程式:∂ρv/∂ｔ＝－∇ｊvは,
∂ρv/∂ｔ＋ｕ∇ρv＝Ｄv∇2ρv (1) を意味します。
これは,Fick型の(移流)拡散方程式(diffusion equation)と呼ばれる方程式です。
一方,ρwを液体の水の密度として半径ａの球水滴の質量をｍとするとｍ＝∫ρwｄＶ＝4πａ3ρw/3です。
水滴の中心を原点:Ｏとするような極座標で考えると,ｒ＞ａの領域では水蒸気に対して,上述の拡散方程式:∂ρv/∂ｔ＋ｕ∇ρv＝Ｄv∇2ρv (1)が成立しますが,ｒ＜ａの球内は(1)が成立しない特異領域と考えられます。
そして,空気は水滴表面から内部に浸入することはないので,表面ｒ＝ａでは速度の表面への垂直成分:ｕr＝0 であるべきですが,現実の水滴表面ではより強い粘着条件:ｕ＝0 が満たされていると思われます。
つまり,表面ｒ＝ａではｕ＝0 なる境界条件が満たされると考えます。
水滴の質量保存から,明らかにｄｍ/ｄｔ＝－∫r=aｊvｎrｄＳですが,
上記境界条件を要求すると,ｒ＝ａでは
ｊv＝ρvｕ－Ｄv∇ρv＝－Ｄv∇ρvなので,
ｄｍ/ｄｔ＝Ｄv∫r=a(∂ρv/∂ｒ)r=aｄＳ (2)です。
特に,表面だけでなく全空間でｕ＝0 を満たす,水滴に対して静止した空気中でのｄｍ/ｄｔを(ｄｍ/ｄｔ)0と書けば,場は中心対称ですから,
(ｄｍ/ｄｔ)0＝4πａ2Ｄv(ｄρv/ｄｒ)r=a (3)と書けます。
他方,ｕ＝0 なら∂ρv/∂ｔ＝Ｄv∇2ρvです。さらにｔ→ ∞ の極限では拡散場は定常(∂ρv/∂ｔ＝0)になり,しかも球対称ですから,
∇2ρv＝ｒ-1ｄ2(ｒρv(ｒ)/ｄｒ2＝0 です。
これを,ｒ＝ａでρv＝ρv(ａ);ｒ＝∞でρv＝ρ∞の境界条件下で解けば
ρv(ｒ)＝ρ∞＋ａ{ρv(ａ)－ρ∞}/ｒ (4)を得ます。
これにより,(ｄρv/ｄｒ)r=a＝{ρ∞－ρv(ａ)}/ａ (5)ですから,
ｕ＝0の静止大気中での水滴の質量保存の式は,
(ｄｍ/ｄｔ)0＝4πａ2Ｄv(ｄρv/ｄｒ)r=a＝4πａＤv{ρ∞－ρv(ａ)}
と書けます。
水蒸気が理想気体なら,それが従う状態方程式:ｅ＝ρvＲＴ/Ｍwから,
ρ∞＝ｅ∞Ｍw/(ＲＴ∞),ρv(ａ)＝ｅaＭw/(ＲＴa)を代入して,
(ｄｍ/ｄｔ)0＝(4πａＤvＭw/Ｒ)(ｅ∞/Ｔ∞－ｅa/Ｔa) (6)
です。
または,ｍ＝4πａ3ρw/3なる式により,
方程式:ａ(ｄａ/ｄｔ)0＝{ＤvＭw/(ρwＲ)}(ｅ∞/Ｔ∞－ｅa/Ｔa) (7)
を得ます。
ただし,Ｔaは水滴表面(ｒ＝ａ)の温度,ｅは水蒸気圧,Ｒはモル気体定数(Ｒ～ 8.31J/(molＫ)),Ｍwは水の分子量(Ｍw～ 18ｇ/mol)です。
しかし,水滴半径ａが分子半径と比較してそれほど大きくない場合には,空気や水蒸気を連続体とみなす巨視的流体力学の適用範囲の外に出ると考えられるため,必要な小水滴に対する補正をします。
以下,Fuchs(1959)によって与えられた方法に従います。
拡散方程式とその解は水滴表面からの水蒸気の反跳長さ(recoiling length):Δvを超えないなら正しくないと考えます。
すなわち,ｒ≧ａ＋Δvではρv(ｒ)＝ρ∞＋Ａ/ｒ(8)が成立しますが,
ａ≦ｒ≦ａ＋Δｖの境界領域では水蒸気流束を構成する分子はMaxwell-Boltzmannの速度分布に従う粒子として気体運動論的に扱うことにします。
水分子の質量をｍwとし,分子の平均速度をｖ＝(ｖx,ｖy,ｖz),その大きさをｖ≡|ｖ|,Boltzmann定数をｋB≡Ｒ/ＮA(ＮAはAvogadro数)とすると,温度がＴのときの分子の平均速さは,
＜ｖ＞＝(8ｋBＴ)1/2/(πｍw)1/2＝(8ＲＴ)1/2/(πＭw)1/2 (9)
で与えられます。
※(注):＜ｖ＞＝∫ｖexp{－ｍwｖ2/(2ｋBＴ)}ｄ3ｖ/∫exp{－ｍwｖ2/(2ｋBＴ)}ｄ3ｖ＝∫0∞ｖ3exp{－ｍwｖ2/(2ｋBＴ)}ｄｖ/∫0∞ｖ2exp{－ｍwｖ2/(2ｋBＴ)}ｄｖ です。
ところが,∫0∞exp(－αｖ2)ｄｖ＝(1/2)(π/α)1/2ですから,両辺をαで微分することで,∫0∞ｖ2exp(－αｖ2)ｄｖ＝(1/4)π1/2/α3/2を得ます。
さらに,これをαで微分すると,∫0∞ｖ4exp(－αｖ2)ｄｖ＝(3/8)π1/2/α5/2です。
一方,∫0∞ｖ3exp(－αｖ2)ｄｖ＝(1/2)∫0∞ｕexp(－αｕ)ｄｕ＝(1/2)(1/α2)より,∫0∞ｖ3exp(－αｖ2)ｄｖ/∫0∞ｖ2exp(－αｖ2)ｄｖ＝2/(πα)1/2です。
＜ｖ＞＝∫0∞ｖ3exp{－ｍwｖ2/(2ｋBＴ)}ｄｖ/∫0∞ｖ2exp{－ｍwｖ2/(2ｋBＴ)}ｄｖの右辺は,上記のα＝ｍw/(2ｋBＴ)の場合に当たりますから,
＜ｖ＞＝(8ｋBＴ)1/2/(πｍw)1/2＝(8ＲＴ)1/2/(πＭw)1/2です。
一方,二乗平均速度を評価するのなら,
＜ｖ2＞1/2＝[∫ｖ2exp{－ｍwｖ2/(2ｋBＴ)}ｄ3ｖ]1/2＝[∫0∞ｖ4exp{－ｍwｖ2/(2ｋBＴ)}ｄｖ/∫0∞ｖ2exp{－ｍwｖ2/(2ｋBＴ)}ｄｖ]1/2
＝{3/(2α)}1/2＝(3ｋBＴ/ｍw)1/2＝(3ＲＴ/Ｍw)1/2です。
後者はエネルギーの等分配則＜ｍｖ2/2＞＝(3/2)ｋBＴから予想される結果に一致しています。
(8/π)1/2と31/2 には大差がないので,＜ｖ＞と＜ｖ2＞1/2は同じオーダーの量です。(注終わり)※
さて,ｎを分子濃度,つまり単位体積当たりの分子数とします。
分子の従うMaxwell-Boltzman分布を,
"ｖ＝(ｖx,ｖy,ｖz)とｖ＋ｄｖ＝(ｖx＋ｄｖx,ｖy＋ｄｖy,ｖz＋ｄｖz)を対角線とする速度空間の微小直方体内に,単位体積当たりＦ(ｖ2)ｄｖxｄｖyｄｖz個の分子があること"で表現すれば,
∫Ｆ(ｖ2)ｄｖxｄｖyｄｖz＝ｎ,∫ｖＦ(ｖ2)ｄｖxｄｖyｄｖz/ｎ＝＜ｖ＞
です。
すると,速度空間では速度ｖがｖ～ ｖ＋ｄｖの間にあり,かつ座標空間の任意の向きの平面の微小立体角:ｄΩ＝cosθｄθｄφに衝突する単位時間,単位面積当たりの分子数は,
{ｖsinθｄΩ/(4π)}Ｆ(ｖ2)ｄｖxｄｖyｄｖz
＝ｖＦ(ｖ2)ｄｖxｄｖyｄｖzsinθcosθｄθｄφ/(4π)
で与えられます。
ただし,この平面がｖとなす角がθであるように極座標軸を取ってます。
平面の一方の側の単位面積に衝突する単位時間当たりの分子の総数は,これを角度θについて 0 からπ/2まで,φについて 0 から2πまで,さらにｖについて 0 から∞まで積分した値です。
∫0π/2sinθcosθｄθ∫02πｄφ＝πですから,この値:[∫ｖＦ(ｖ2)ｄｖxｄｖyｄｖz∫0π/2sinθcosθｄθ∫0π/2ｄφ]/(4π)は,(1/4)∫ｖＦ(ｖ2)ｄｖxｄｖyｄｖz＝(1/4)ｎ＜ｖ＞に等しいことがわかります。
したがって,(分子数密度流束)＝(任意の面に衝突する単位時間,単位面積当たりの分子数)は,
(1/4)ｎ＜ｖ＞ (10) と書けます。
右辺の因子:＜ｖ＞は先に書いたように,
＜ｖ＞＝(8ｋBＴ)1/2/(πｍw)1/2＝(8ＲＴ)1/2/(πＭw)1/2 (9)
で与えられます。
また,水蒸気の温度を改めてＴの代わりにＴgと書くと,水蒸気はｅ＝ρwｋBＴg/ｍw＝ｎｋBＴgなる状態方程式に従います。
故に,水滴の周囲の水蒸気の分子数密度流束は
ｗ↓＝(1/4)ｎ＜ｖ＞＝(1/4){ｅ/(ｋBＴg)}(8ｋBＴg)1/2/(πｍw)1/2
＝ｅ/(2πwｋBＴg)1/2です。
他方,水滴表面に付着する水蒸気の分子数密度流束をｗ↑と書き,平衡状態にあるときのｗ↑の水滴の周囲の分子数密度流束ｗ↓に対する比を凝結定数と呼んでαcとします。
すると,ｗ↑＝αcｗ↓です。
そうすると,水滴表面の分子数密度流束は,
ｗ↑＝αcｅ/(2πｍwｋBＴa)1/2 (13)
で与えられることになります。
実験によれば,αcの値は,αc≒0.01～ 0.07です。
ただし,水蒸気分子の平均速さは,＜ｖ＞＝(8ＲＴg)1/2/(πＭw)1/2でも
分子速度ｖの分布は等方的で,そのベクトル平均はゼロと考えられるので
分子数密度流束ｗ↓の向きはあらゆる方向について対称です。
水滴表面の水蒸気圧:ｅが飽和水蒸気圧:ｅsatに等しいときに"定常状態＝蒸発平衡"に達するとすれば,
平衡時の水滴表面での飽和分子数流束:ｗ↑satは,
ｗ↑sat＝αcｅsat/(2πｍwｋBＴa)1/2 (14) です。
ただし,Ｔaは水滴表面の絶対温度です。
それ故,流入する分子数流束と流出する分子数流束は互いに独立であるとしてＴg～ Ｔa ＝Ｔと置けば,水滴に流入する総流束密度の大きさは
ｗsat＝ｗ↑－ｗ↑sat＝αc(ｅ－ｅsat)/(2πｍwｋBＴ)1/2(15)
と書けます。
この(15)式は,ｅ＜ｅsatなら水滴が蒸発し,ｅ≧ｅsatなら蒸発しないことを意味しています。
さて,水滴中心からの距離ｒにおける水蒸気の分子数密度をＮv(ｒ),密度をρv(ｒ)とすればρv(ｒ)＝ｍwＮv(ｒ)です。
したがって,ｒ＝ａ＋Δｖにおける水蒸気が,流出入する水蒸気を代表するとすれば,その水滴表面へと向かう凝結流束密度は,
αcｍwＮv(ａ＋Δｖ)＜ｖv＞/4＝αc＜ｖv＞ρv(ａ＋Δｖ)/4 です。
他方,水滴表面ｒ＝ａで蒸発に向かう蒸発流束密度は,
αc＜ｖv＞ρv(ａ)/4です。
そこで,表面での総蒸発流束Ｊv(ａ)は,
Ｊv(ａ)＝4πａ2ｊv(ａ)＝πａ2αc＜ｖv＞{ρv(ａ)－ρv(ａ＋Δc)} (16)
で与えられることがわかります。
一方,冒頭で書いたように,巨視的流体力学の見方では,
ｒ＝ａでの流速密度はｊv＝－Ｄv∇ρvなので,総流束は
Ｊv(ｒ)＝4πａｒ2ｊv(ｒ)＝－4πｒ2(ｄρv/ｄｒ)Ｄv
です。
ところで,今の気体運動論的考察では,
ｒ≧ａ＋Δｖでρv(ｒ)＝ρ∞＋Ａ/ｒ (8)が成立する,
と仮定しているので,
ｒ＝ａ＋Δｖにおいてはρv(ａ＋Δv)＝ρ∞＋Ａ/(ａ＋Δv)です。
故に,(ｄρv/ｄｒ)r=a+Δv＝－Ａ/(ａ＋Δc)2ですから,
Ｊv(ａ＋Δv)＝－4π(ａ＋Δv)2(ｄρv/ｄｒ)r=a+ΔvＤv＝4πＤvＡ (17)
を得ます。
"定常状態＝蒸発平衡"では,Ｊv(ａ＋Δv)＝Ｊv(ａ) (18)であるべきで,
ρv(ａ＋Δv)＝ρ∞＋Ａ/(ａ＋Δv)ですから(16),(17)より,
4πＤvＡ＝πａ2αc＜ｖv＞{ρv(ａ)－ρv(ａ＋Δc)}
＝πａ2αc＜ｖv＞{ρv(ａ)－ρ∞－Ａ/(ａ＋Δv)}
が成立します。
したがって,
[ａ2αc＜ｖv＞/(ａ＋Δv)＋4Ｄv]Ａ＝ａ2αc＜ｖv＞{ρv(ａ)－ρ∞},
および＜ｖv＞＝(8ｋaＴa)1/2/(πｍw)1/2＝(8ＲＴa)1/2/(πＭw)1/2 (9)
が成立します。
以上から,
Ａ＝ａ{ρv(ａ)－ρ∞}/[ａ/(ａ＋Δv)＋4Ｄv/(ａαc＜ｖv＞)]
＝ａ{ρv(ａ)－ρ∞}/[ａ/(ａ＋Δv)＋Ｄv{2πＭw/(ＲＴa)}1/2/(ａαc)] (19) を得ます。
これは,冒頭で述べた巨視的流体力学の考察:
(4)ρv(ｒ)＝ρ∞＋ａ{ρv(ａ)－ρ∞}/ｒ,
および(ｄｍ/ｄｔ)0＝4πａＤv{ρ∞－ρv(ａ)}から導かれた式,
(6)(ｄｍ/ｄｔ)0＝(4πａＤvＭw /Ｒ)(ｅ∞/Ｔ∞－ｅa/Ｔa),
または(7)ａ(ｄａ/ｄｔ)0＝{ＤvＭw /(ρwＲ)}(ｅ∞/Ｔ∞－ｅa/Ｔa)
を次のように補正することに相当します。
すなわち,(4)はａ{ρv(ａ)－ρ∞}→ Ａ＝ａ{ρv(ａ)－ρ∞}/[ａ/(ａ＋Δv)＋Ｄv{2πＭw/(ＲＴa)}1/2/(ａαc)]として,
ρv(ｒ)＝ρ∞＋ａ{ρv(ａ)－ρ∞}/ｒ→ (4)'ρv(ｒ)＝ρ∞＋Ａ/ｒ
と補正します。
そこで,(6),および(7)は拡散係数をＤv→ Ｄv'＝Ｄv/[ａ/(ａ＋Δv)＋Ｄv{2πＭw/(ＲＴa)}1/2/(ａαc)](20)として,
それぞれ(ｄｍ/ｄｔ)0＝(4πａＤv'Ｍw/Ｒ)(ｅ∞/Ｔ∞－ｅa/Ｔa)(6)',
およびａ(ｄａ/ｄｔ)0＝{Ｄv'Ｍw/(ρwＲ)}(ｅ∞/Ｔ∞－ｅa/Ｔa)(7)'
に補正します。
これで小水滴に対して補正された質量保存式を得ました。
次に,水蒸気を含む空気に対して運動している水滴に対する補正を考察します。
これは,水滴が運動する代わりに静止水滴に対して空気が運動する,つまりｕ≠0 の等価な問題として設定できます。
この場合も,ｒ＞ａにおける拡散は定常(∂ρv/∂ｔ＝0)に達しているとすれば,質量保存の拡散方程式はｕ∇ρv＝Ｄv∇2ρvです。
水滴表面での流束密度がｊv＝－Ｄv∇ρvなので,
やはりｄｍ/ｄｔ＝－∫r=aｊvｎrｄＳ＝∫r=aＤv(∂ρv/∂ｒ)ｄＳです。
ここで,平均通風係数と呼ばれる量ｆvを,
ｆv≡ｄｍ/ｄｔ/(ｄｍ/ｄｔ)0(21)で,
平均輸送係数ｋvをｋv≡－∫r=aＤv(∂ρv/∂ｒ)ｄＳ/[4πａ2{ρv(ａ)－ρ∞}] (22) で定義します。
すると,(2)'と(21),(22)からｋv＝Ｄvｆv/ａです。
さらに,無次元のSherwood数Ｓhを次式で定義します。
すなわち,Ｓh≡2ｋvａ/Ｄv＝2ｆv (23)です。
こう定義すると,ｄｍ/ｄｔ＝ＳhＤv2πａ(ρ∞－ρsat) (24)です。
ｕ＝0 のときにはｆv＝1 or Ｓh＝2です。
ただし,これらは単に定義した,または言葉を置き換えただけで何ら新しいことを述べていません。文献を読む際の便宜だけでしょう。
§2.水滴中の熱の移動
ｋを媒質中の熱伝導度,Ｔを絶対温度とすると,Fourierによれば熱流束ベクトル:ｊhはｊh＝－ｋ∇Ｔ (25)で与えられます。
そして,一般的な流体中のエネルギー保存の式は,粘性による散逸が無視できる場合には流体単位質量あたりのエントロピーをｓとしてρＴ(∂ｓ/∂ｔ＋ｕ∇ｓ)＝－∇ｊhと書けます。ρは流体の密度です。
今の場合のように,媒質流体(水,空気)が非圧縮と見なせるときには
∇ｕ＝0,かつＴｄｓ＝ｃpｄＴ(ｃpは単位質量当たりの定圧比熱)なので,
これはρｃp(∂Ｔ/∂ｔ＋ｕ∇Ｔ)＝κ∇2Ｔ
or ∂Ｔ/∂ｔ＋ｕ∇Ｔ＝κ∇2Ｔ(26) と変形されます。
ただし,κ≡ｋ/(ρｃp)でκは熱拡散係数と呼ばれる量です。
ｊh＝－ｋ∇Ｔ(25)より,水滴の熱量をｑとし,乾燥空気の熱伝導度をｋaとすると,
ｄｑ/ｄｔ＝－∫r=aｊhｎrｄＳ＝∫r=aｋa(∂Ｔ/∂ｒ)ｄＳ (27)
です。
特に,ｕ＝0(水滴に対して水蒸気静止)のときのｄｑ/ｄｔを(ｄｑ/ｄｔ)0と書きます。
ｕ＝0 の熱伝導方程式:∂Ｔ/∂ｔ＝κ∇2Ｔのｔ→∞での定常状態(∂Ｔ/∂ｔ＝0)の∇2Ｔ＝ｒ-1ｄ2(ｒＴ)/ｄｒ2＝0 の解:Ｔ＝Ｔ∞＋ａ(Ｔ∞－Ｔa)/ｒから,(ｄｑ/ｄｔ)0＝4πａｋa(Ｔ∞－Ｔa)(28)です。
ａが小さい小水滴のときには,先に拡散係数Ｄhの代わりにＤh'を用いるよう補正したのと同じく,熱伝導度ｋaの代わりにｋa'≡ｋa/[ａ/(ａ＋ΔT)＋ｋa{2πＭw/(ＲＴa)}1/2/(ａαT＋ρｃp)](29)を用います。
なお,ΔT,αTは拡散での反跳長さΔv,凝結定数αcに相当する熱伝導(熱拡散)での量です。
また,水滴に対して水蒸気が運動するときには,Nusselt数と呼ばれる数ＮuをＮu≡2ｄｑ/ｄｔ/(ｄｑ/ｄｔ)0で定義します。
(ｄｑ/ｄｔ)0＝4πａｋa(Ｔ∞－Ｔa) (28)からｄｑ/ｄｔ＝Ｎuｋa 2πａ(Ｔ∞－Ｔa)(30)と書けます。ｕ＝0 のときにはＮu＝2です。
今日はここまでですが,まだ続きます。
なお,もちろんオリジナルではなく当時,参考文献を見たはずなのですが,それをノートに書いてないので今は不明です。
確か1983年当時に勤めていた会社の蔵書だったと記憶してます。
PS:新しい試みとしてワードのルビ(ふりがな)を記号としてココログにコピーしてみたのですが取り合えずのところは失敗でした。おかげで文章を校正する時間を浪費しました。
貧乏ひまなしで師走は日曜日でも時間がありません。バスに遅れるので続きはあとで。
PS:民主党だけドタバタしてるように見えるのは与党だから仕方ないとはいえ,こんなとき野党は何のためにあるのか？責任追及だけやるしかないのか？そもそも代議士とは何をするのが仕事なのか？
政界再編はいずれ必要なんだろうけど,今やられてもそれが終わる頃には私自身の寿命などは終わってるだろうなあ。。。
イヤ,またまたマスコミや評論家の多くと同じく,自分のこと抜きで他人を揶揄しちまったか。
なんとか法とかの法律に賛成反対とか,国や公共団体の迷惑行為,大きなお世話はやめて欲しいけど,基本的に政治のお世話には期待してないアナーキーを自認している自分。。
自分の世話さえままならず,社会的にも無力で有象無象の自分がいまさら何の能書きをたれてるのか。お笑いだ。。。
PS2：機械文明や都会のコンビニエンスに依存し切っている身はいまさらどうしようもないという気分もあるけれど,もしもここ10年くらいの間,PCやネットがなかったら,既に廃人になってたかも。。
今のような状態でも,自分が生き生きと生きてるように感じることが多々あるのは悲しいことかも知れないけど,血の通ってないPCやネット依存症に身を委ねて合理化しているおかげでしょうネ。
。
最近のコメント