ソフトの設計ミス

液状化の検討（震度法）を実施する必要がありました。

数か月前に購入したサービスに含まれていたソフトで実施。どのような使い勝手か試してみたかった、というところです。

が、計算結果が何かおかしい。
よくよく見ていくと有効上載圧の計算結果が誤っています。なぜそのような計算結果になるのか調べてみると、わかりました。湿潤単位体積重量から 10kN/m3 を引く仕様になっていました。飽和・水中を入れる箇所がなく、どうにも手出しできません。

サポートに連絡しましたが、担当者不在とのことで、後から連絡するとのことでした。

その日の夕方に、担当者のオジサマから電話がかかってきました。
が、話が通じません。どうも、湿潤・飽和の違いがわからない土木素人サンのようでした。怖いですね~。実はこの会社、以前もCAD関連でこのような素人サン担当者が出てきて呆れたことがありました。有名な会社なのですが、ソフト開発は外に丸投げのようです。

最近、このような質の低下が目につきますね。市販ソフトでも「まず確認」「疑ってかかれ」が身に付きました。残念な流れです。

結局、いつものソフトで計算を終えました。
使い勝手が良かっただけに、残念ですね~。

DRASTIC Model

共立出版「地下水流動」を読んでいますと、DORASTIC Model といった評価法が紹介されていました。

地理、地形の分野で作られたのでしょうか？GISを利用した帯水層の脆弱性を表すツールだそうです。知りませんでしたね。ArcHydroとは考え方が異なりますが、同じようなにおいがします。
「脆弱性」が具体的にどのような物理現象を指すのか？それとDRASTICがどのように関連する理論なのかは不明ですが、何らかの答えを出せる点は良いと思います。

ま、こういう手法がある、ということは覚えておきましょう。

やり残し事項 2018

研究職の実績公開がありました。

色々なことに熱中されていました。基礎的な事項から応用まで。レベルの高いものから低いモノまで。自信が満ち溢れていました。
個人的には得られるものがあったり、気分的にへこんだり。

CIMの実績公開もありました。
コチラはまだ底上げ主体で、得るものはなし。

いずれにしても、業務で得られた利益をこのような玉石入り混じった内容に費やせる余裕があるということでしょう。

研究職と話をして、いくつかの進展がありました。
再度、整理しておきましょう。

道具
帯磁率計
ガンマ線測定器（携帯型）

技術
動的解析（耐震、液状化）
斜面設計

コード・ソフト

• 優先度低：DtransuのCUDA化・・・先月失敗。100万で購入可。
• 優先度低：PSIｎSAR・・・700万で購入可
• 優先度中：地表流＋地下水・・・GSFLOW、HydroGeoSphere、GETFLOWS・・・論文を読みましょう。
• 優先度中：GMS・・・J-SHISの取り込みまではできました。いつか続きを。
• 優先度低：AWS, Azure
• 優先度中：DualSPHysics・・・Foam Simulation が残っています。
• 優先度中：HEC-RAS、iRIC、河床変動計算
• 優先度低：DeepWater・・・よくできています。H2O4GPUで我慢するか。
• 優先度中：OpenMI・・・PHREEQCとDtransuの連成等

地下水流動

共立出版「地下水流動」谷口他

以前、CFCs の測定を行っていた時に書店で見かけました。関連する記述が出ていたレアな図書だったのですが、知りたい部分の内容が浅く、かつより良い資料があったので購入しませんでした。

余裕ができたので購入し、読んでいるところです。
内容は広く浅くといったところでしょうか。ある分野で最新の内容に特化したような図書ではありません。
それでも、報告書に引用するには良い文章が載っています。備忘録として残しておきましょう。

p106

モンスーンアジアの湿潤地域では、地形（によって高められた地下水面）は地下水の最も重要な駆動力を提供する。地質は地下水流動系の構造を修飾するが、全体としては地形の高低が地下水の流れを支配し、局地地下水流動系から地域地下水流動系に至る地下水流動系の階層構造を形成する。

p111

しかし、いずれのメカニズムでも、地下環境が変化しなければ地下水中に高濃度でヒ素が溶出することはないので、モンスーンアジア地域における地下水中のヒ素濃度上昇の原因は、揚水など何らかの人間活動による地下環境の変化が原因と考えるのは妥当である。

地下水はサイエンスとしては面白い題材でしょう。
まだまだわからないことが多い分野だと思います。

昇格試験 その2

executive より管理職に関するヒアリングがありました。

「社内の昇格試験自体に興味はない。」「総監取得を昇格試験代わりにすれば良い。」等と伝えていますと、思いもよらぬ答えが。

executive 曰く「昇格試験は管理技術を見る試験ではなく、会社に忠誠を誓うかどうかの試験」「管理技術を問うなら、毎年同じような問題は出さない」

合点がいきました。

その前に、社内の組織運営上の問題を話していたのですが、それを招いたのは現在の管理職。要は監理能力の有無で管理職を選ばず、忠誠心の有無で管理職を認定していたため、問題がひき起ったのです。当然、その責任をとるという発想は出てきません。反省し監理能力を発揮する素質のある方を選抜する試験に変更する意思もありません。せいぜい、尻拭いを他者に押し付ける程度です。

個人的には、「まず問題点を認識する」「責任を取る」「管理技術を昇格試験で問う」等から始めるのが理想かと思います。が、無理でしょうね。ヒトの世ですから。

忠誠を誓うだけなら試験でなくて一筆ですませばお互い楽でしょうに。
どうしてヒトって形式や非効率なことを好むのでしょうか？

DualSPHysics 4.2

DualSPHysics が 5月の頭に Ver.UP しました。

あきらめていた foam simulation を確認しておきたかったのですが、これまで時間がありませんでした。ようやく時間ができたのでDL。

サンプル内に bat ファイルがあります。手早く確認するため、今回はプロンプトからそれを指定！

が、動きません。
CUDA ライブラリの Ver.が違うそうです。

VS の プロジュクトが同梱されていましたので、まず確認。
が、開けず。
CUDA 9.1 のファイルが見つからない、とのこと。

そういえば、PGI Community Edition 17.10 が CUDA 9.0 までしかサポートしていなかったので、9.0 しか入っていません（こういった時に、Docker の優秀さが再認されます）。
VS 上で 9.1 のファイル指定を 9.0 に書き換えてコンパイル。exe がどこにできたのか迷いましたが、bin フォルダの中の exe が更新されていました。配布元の配慮が伺えます。

今度は問題なく走りました。
01_DamBreak + GTX1060 で計算時間7分。短くなったような気がします。

ひとまず、手元の CUDA の Ver. にあわせてコンパイルすることで、環境は整いました。
次は foam simulation ですね。

河床変動計算の設定

砂防学会編「山地河川における河床変動の計算方法」の備忘録。

なぜ不等流？といった疑問に対する答えです。p15

河床変動の伝播の速さは水面変動のそれに比べてはるかに遅いため、流れの時間的な変化を無視しても良い場合が多くある。つまり、河床変動に着目して観察を行うと、流れの時間変動がはるかに速いため、ほとんど瞬間に終了してしまい、時間的に変化のない定常流として取り扱うことが可能だということである。こうして、河床変動に対しては時間的な変化を考えて式・・・・を用いるが、流れに対しては空間的な変化はあっても時間的に変化がないとして・・・・時間微分項を除いた、不等流の方程式を用いて議論することが多くある。

p51-54

設定
Case1：d=dmとした一様砂礫モデル
Case2：交換層厚さa=初期砂礫層の最大粒径
Case3：a=最大粒径の2倍
結果
Case1では河床勾配が全区間にわたって急速に一様化し、平均勾配に達する。Case2,3では遅い。洗堀区間での粗粒化と堆積区間の細粒化はaが小さいほど顕著。これが河床形状の変化傾向の相違の原因。
交換層の厚さの適正な設定法は確立されていない。

一様砂礫モデルでの急速な一様化は感じていたところです。
a=最大粒径として、粒度も複数階に見直しましょう。