以下の前記事への意見です。 computational-sediment-hyd.hatenablog.jp 素晴らしい記事ありがとうございます。勉強になります。 の影響がかなり効きそうですね。 非保存形（一般的な非保存形と混合しそう）の方が流量が保存されるというとカオス感が。。。 …

久々に。。。 モチベーション ひょんなことからTVD-MacCormack法を実装しようと思い、山地河川における河床変動の数値計算法（砂防学会の「急流河川の１次元河床変動（その２）」）を見て実装しようと思ったのですが、川幅に変化がある場合の計算が上手くい…

Rouse分布の近似式であるLane-Kalinske分布って結構荒っぽくないかなと昔から気になっていたので検証してみました。 ついでに式の導出も手書きしか無かったのでTeXにまとめておきました。 土砂濃度の鉛直分布の導出 基本的な考え方 Rouse分布 参考：式展開 L…

jupyterでTeXを書くことが多いのですが、はてなブログ（markdown記法）に貼ると上手く行かないので変換用のスクリプトを書いておきました。 多分awkとかperlでもできるけど苦手なのでpythonで。 はてブロ用のTeXに対応するため、以下の処理を行う。 ※数式はa…

循環シフトとは 配列の最初と最後を循環してシフトするものです。 pythonのnumpy.rollで書くとこんな感じになります。 import numpy as np x = np.arange(5) x # array([0, 1, 2, 3, 4]) np.roll(x,2) # array([4, 0, 1, 2, 3]) np.roll(x,-1) # array([1, 2…

年に1回くらいこの記事が見たくなる。 karapaia.com 計算してみたいな。平面二次元で行けるのかな。三次元じゃないと無理かな。

前回の記事で椿先生の本1を参照してときに、ペラペラのめくっていると安定河道の話があったのでそれを元に考えてみました。 河道の安定条件のモデル化 簡単のため、一次元の不等流場、断面形状は矩形断面を想定する。また、河床材料は単一粒径で掃流運動のみ…

河道内貯留のついて勉強のためにも一度計算してみました。 河道内貯留とは 河道内貯留とは、下図1のように下流への伝播に伴いピーク流量（水深）が減衰する現象を示す。 このような伝播を計算するためには洪水波の移流拡散を考慮する必要があるため、次式のD…

頭の整理に書いてみました。ご指摘があれば Manningの公式の概要 Manningの公式への水理学的な意味付け 沖積河川の河床抵抗 河床波の変化を考慮した河床抵抗の評価方法 岸・黒木の方法による河床抵抗の評価 岸・黒木の方法の課題 河床波の遷移の履歴 河床波…

有名なEngelundによる有効掃流力の上図（元図はEngelund、上図は河村のp.226、実験値はUSGSのレポートらしい）は反砂堆領域の曲線の式が設定されてないです。 何かと不便なのでトレースしてみました。 こんな感じです。 トレースしたデータはこちら この図は…

MinGWのmakeはmingw32-make MinGWのmakeコマンドについて - （仮） 最近fortranを書く機会が減りましたが、年に数回書く際にmakeコマンドを忘れてしまいます。 windowsはやっぱり特殊ですね。個人的にはWSLを標準にして欲しいのですがなかなか完全には浸透し…

なんとなく気分転換に。この手の話って最近の本だとさらっとしか書いてないですよね。 開水路において水深平均化された対数則の導出 manning-stricklerによる近似形 対数則について思うこと 参考：対数則の導出 参考：対数則の数値計算上のテクニック 参考文…

The English version is available here はじめに L-BFGS法の導出 L-BFGS法の実装 計算例 scipy.optimize.fmin_l_bfgs_bとの比較 メモ 最急降下法との比較 考察とまとめ Github はじめに 非線形問題の最適化に最もよく使用されるL-BFGS法(記憶制限BFGS法、BF…

モチベーション 横断面平均一次元河川流計算で実河川を対象とする場合、地形条件は河川横断測量データを使用することが一般的である。 このデータより、水位に対する河積、川幅、潤辺等を計算するが計算負荷が大きく、特に長期計算を行う場合には計算速度に…

勉強会で話が出たので整理しました。 速度比較 numba.jit最適化環境での速度比較 まとめ 速度比較 複数の戻り値をタプルとディクショナリで戻す場合の速度を比較してみました。 一応私のPCのスペックです。 ノートPCのスペック : VAIO Pro PG - 趣味で計算流…

ユーザーがかなり限定されますが... ※ テストが十分でないのでエラー等指摘いただけると助かります。 lemファイルとは pythonスクリプト Usage Exmaple 1 Exmaple 2 おまけ Environment Anaconda Google Colaboratory github 参考サイト lemファイルとは lem…

numba.jitによる高速化のクラス編です。公式にもほとんど情報が無いので結構苦労しました。 ※2024/2/25追記：以下も分かりやすいです。 qiita.com numba.jitを知らない方はまず関連記事を読んだほうがスムーズだと思います。 関連記事 computational-sedimen…

湯前町の雨量がすごいということだったので、湯前横谷（気象庁の地上雨量観測）のデータをみてみました。 累加雨量が500mmと半端じゃない雨ですが、近年豪雨続きであまり驚かないです。 台風じゃないので予測しづらいのが厄介ですね。 参考までに東京にこの…

球磨川の浸水推定図が公開されたので洪水ハザードマップと比べてみました。 浸水推定図（2020年7月4日20時作成） 洪水ハザードマップ:浸水想定区域図（想定最大規模） 今回の洪水が想定最大規模ほどは大きくないと思うのでこんな感じなのかな。

今回の記事はブログレベルの話ではないですが、ずっと悩んでいるのでとりあえずまとめておこうかなというところです。 砂州の定義は以下の記事を参照下さい。 computational-sediment-hyd.hatenablog.jp 支配方程式 流水 流砂、河床変動 計算条件 計算結果と…

砂州についてしゃべることは多いのですがまとめたことが無かったので整理しておこうかと。 面白い数値計算結果が出たので、その準備だったりもしもます。 はじめに 砂州の分類 砂州の分類 砂州の特徴 交互砂州(単列砂州) 多列砂州 網状流路 固定砂州 砂州の…

今までは、 import sympy as sp sp.init_printing() で、使えたのに何かのアップデートの関係で import sympy as sp sp.init_printing(use_latex='mathjax') じゃないと表示されなくなってしまいました。 原因は調査中です。

モチベーション このブログでもしばしば登場するnumba.jitを3D計算でどれほど使えるのかを 前記事(SIMPLE法でもpythonならシンプルに書けるのか ⇒ かなりシンプルになるがやっぱり...(ソースコードあり) - 趣味で計算流砂水理)のSIMPLE法による3Dキャビテー…

最近科学技術計算を書いてなかったのでウォーミングアップにこんな感じの3次元キャビティー流れをSIMPLE法で書いてみました。 大きいものはこちら はじめに 計算モデル 計算条件 ソースコード 仮の速度の計算式 圧力補正値の計算式 全体のループ コード全体 …

久しぶりにParaviewを使うおうとインストールしたら不具合が出たので。 元ネタIntel GPU driver update 26.20.100.6708 on windows breaks Qt GUI --> FIXED by UPGRADING intel drivers (#19364) · Issues · ParaView / ParaView · GitLab が英語しかなかっ…

CPU : intel core i7-8550U@ 1.80GHz 1.99GHz RAM : 16 GB OS : Windows 10 Pro 64bit

HHKBに慣れてしまうとノートPCを使ったときにCaps Lockを間違えてしまうので。 公式 Ctrl2cap | Microsoft Docs 参考サイト ５分でcapslockキーをctrlキーに変更する(windows10) | ジョイタスネット アップデートされてないしレジストリいじるから危ないけど…

なぜHoloviews？ コーディング：matplotlibとの比較 Holoviews matplotlib どうしても静的なグラフを書きたい場合 まとめ 参考図書 なぜHoloviews？ このブログに登場するグラフは、ほぼ全てHoloviewsを使っています。(正確には、PyViz-HoloVizファミリーのh…

以前の計算(多列砂州河道の再現計算（失敗例） - 趣味で計算流砂水理)をやり直してみました。 ましになりましたが、もう少しですね。 5000hrで計算機が落ちてしまいました。ちなみに計算時間は1ヶ月くらいです。

河道内に繁茂する木とか草についての話題です。河道内植生とか呼ばれており、河川分野では一つの研究テーマだったりします。 まず、大きな出水のあった2019年の荒川の熊谷周辺の経年変化をplanetでレンダリングしてみます。 画像を右クリックで再生を選択 大…