趣味で計算流砂水理

Computational Sediment Hydraulics for Fun

R

Upwind Flux Scheme:不連続面の修正2

R

前回上げたものが若干間違っていたので、修正しました。 せっかくシンプルなスキームなのにどうしても複雑になってしまいます。 河床高と水位の関係より不連続面の判定をしていますが、dxが大きくなった場合、この判定が上手く使えるか微妙なところです。 概…

環境構築:bash on ubuntu on windowsのgcc,g++,gfortranのアップデート

R

bash on ubuntu on windows(Windows Subsystem for Linux:以下WSL)のgcc等は、5.4で、c++11以降あるいはf2003以降が使えません。なので、最新を入れました。 http://ftp.tsukuba.wide.ad.jp/software/gcc/releases/gcc-7.2.0/ から最新版の7.2をダウンロー…

QGISでSAGAが動かないときの対処法:windows10-64bit環境

R

最近は仕事モードになってきて、以前ぶっ飛んだPCの環境を整えてます。 めずらしく備忘録的にあげておきます。 私はなぜかSAGA(System for Automated Geoscientific Analyses)が好きで、QGISで処理をするときは、ほとんどSAGAプラグインを使います。 先日P…

seabornで日本語表示

R

何かとバタバタしていて全然書けないです。 今月は書きます。 seabornで日本語表示 windowsが変わるごとにフォントが変わるのでもう大変です。 多分この方法が一番ラクかなと思います。 ここを参考にしました。 ここからIPAexフォントをダウンロードします。…

Upwind Flux Schemeの問題点と修正

R

2d Upwind Flux Schemeの計算で厳しい条件の場合、全然回らなくて苦労していたのですが、1dでいろいろテストして、ようやく見つけました。 Upwind Flux Schemeに限ったことではないですが、水面の不連続面が出てきた場合、厄介なことになります。コロケート…

続:改めて不等間隔格子の河川流一次元計算

R

冷静に考えたら、座標変換ですね。 基礎式です。 実河川の問題では必須のような気がします。 今度テスト計算をしてみます。

改めて不等間隔格子の河川流一次元計算

R

改めて河川流一次元計算を考えていると、計算格子の定義点が測量点で決まるという一般的な流体解析ではあり得ない状況です。 顕著な不等間隔の場合、コロケート、スタッガードに関わらず、コントロールボリュームが以下のようになってしまいます。 なので、…

いろいろ

R

地理院地図とdeep learnning http://ccpn.gsi.go.jp/meeting_partners/data/20170608/4.pdf この前話そうと思っており、忘れてました。 deep learnningの部分を抜きにすると、公開データを一括処理でやっちゃうあたりは、我々がやりたいことと少し近いような…

勉強会での話

R

並列化 大きなプログラムを書きたいため、並列化の勉強をしているですが、気になる内容がありましたので、列挙しておきます。 ハイブリッド並列化 最近の主流っぽいですね。MPI+openmpが基本で。MPIで分割したものに、openmpで並列化をかけるようです。情報…

σ座標系に手を出そうかと

R

連投になってすいません。 VOFを書いたところで、自由表面流れにおけるデカルト座標系の限界を感じてしまいました。河川流には向いてないような気がします。周期境界条件も作りにくいし。 そこで、今まで敬遠していたσ座標系に手を出そうかと思ってます。 Wu…

ポインタとメモリリークについて

R

先日のアップしたプログラムを書いた際に初めてメモリリークという壁にぶち当たりました。 fortranでメモリリークって無いだろと思ってましたが、結構あるみたいです。 例えば、 http://jjoo.sakura.ne.jp/tips/f90/memory_leak.html ただし、今回陥ったもの…

似非VOF:like interfoam

R

以前Openfoamのinterfoamの話がでたので、気になって自分でコーディングしてみました。 延べ2週間以上かかりました。苦労しました。 いろんなことがわかったので書いておきます。 参考資料 このあたりです。 https://www.jstage.jst.go.jp/article/jscejoe/6…

windows環境

R

いろいろありがとうございます。 私の環境は、VS 2017 communityです。問題なく動いてます。intel compilerも入れてますが、VS 2017上で動きます。全ての日本語化されてます。 特に工夫もせずにクリックだけで環境が構築出来たので、よくわからないです。 VS…

基本的な計算機の作り方について

R

linux環境も揃えてるとは… コンパイラは何を使ってますか。gccですか。 基本的なことを教えていただきたいのですが、新規でPCを購入して計算機を作ろうと考えています。 環境はこんな感じをイメージしております。 言語:FortranとC++ CPU:Xeon phi 目的:…

渦度方程式

R

水工学における計算機利用の講習会の山坂先生の「渦度方程式を用いたワンド内流れの二次元計算」でした。 急に思い出しました。

RE3:pygridgen-c

R

pygridgen-cのwindows化 何とかコンパイルまでは行ったのですが、pythonから呼んだときにエラーが出ます。試しに32bit環境で動かすと何故か上手く行きます。 いろいろ調べているとctypesに32bitと64bitで差がありそうです。元のプログラムがpython側でしっか…

RE2:pygridgen-c

R

windows環境でインストールを試みましたが、どうも駄目です。というか、私がmsysを使いこなせてないです。 msys環境下で順番にconfigure - make installを行っていますが、 nearest neighbors、CSAをインストール後、gridutilsのconfigure時にそれらにリンク…

gridgen-c

R

久しぶりにgridgen-cを使ってみました. 同じものですが,githubに上がってました. その中にpythonのフロントエンドとしてpygridgenというものが追加されておりました. さらにその中で中段くらいに記載されているadd_focusという機能が追加されています.…

CQ2DF model + vegetation

R

面白いアイデアですね。 問題があるとしたら、渦動粘性係数の空間分布でしょうか。乱流モデルの影響がどこまで効くかですね。 時間をみてやってみます。 たまたまですが、先日ある方と清水先生の以下の論文を話していました。 https://www.jstage.jst.go.jp/…

デカルト座標系での平面二次元計算

R

昨日もお見せしてみましたが、とりあえずチャレンジしてみました。 計算条件 上流端境界に合わせて格子を設定しています。 下流端はあえて、横断方向と直交しないように設定しています。 一部堤防を切っています。これもあえてです。 初期条件は上手に設定で…

久しぶりになりますが。。。

R

年度末だったり、河川シンポだったりで疲れました。 いろいろな計算をやってるのですが、まとまっていないので文章だけ書きます。 変分法 お見せしているとおり結構良い結果がでたので、今後の展開を考えようかなと思ってます。 今のところは、支川を含む長…

2D Upwind Flux Scheme:Dry-Bed Dam-Break

R

※前々回の記事がまとまっていなかったので削除しました。 2D Upwind Flux Scheme:Dry-Bed Dam-Break ほぼ完成しました。一応ちょっとごまかしましたが、乱流項も入れました。 乱流項ありなしのgifを貼っておきます。 粗度・乱流なし 粗度・乱流あり ソースコ…

Cross section Quasi 2D Flow Model

R

githubに公開しました。 https://github.com/computational-sediment-hyd/CQ2DF windows+proxy環境下でのSourceTreeによるgithubへのアップロード やり方がわかれば一瞬ですが、結構大変でした。 忘れないためにもメモ的に置いておきます。 proxyの設定 prox…

Cross-sectionQuasi2DFlow

R

ギッハブ用のモデルとかを整理しました。DBにあるので見て頂けると助かります。windows環境があれば回してもらえるともっと助かります。 READMEとか結構頑張りました。 問題がなければ上げます。 英語にしたときに違和感があり、名前をCross-sectionQuasi2DF…

Section Quasi 2D Model:dll版が良いのか?

R

ライセンスはMITで行こうと思います。 リリース版として悩んでいるのが、現状のdll版がそうじゃないものかです。 計算の流れは、格子生成(python)⇒流速場の計算(fortran dll)⇒図化(python)ですが、個別の実行ファイルを作成し、テキストファイル経由でやり取…

EXCELのVBAが。。。

R

ふと気付いたのですが、Excel 2016に変更してからExcel 2007で作っていたマクロの中でユーザーフォームを使ってグラフをコントロールするものがすべて動かないです。 多分描画速度の低下(データ更新時処理のアルゴリズムが変わった?)の影響かと思いますが…

ライセンスの話

R

例のSection Quasi 2Dモデルについて、Githubで公開しようかなと考えているのですが、そろそろライセンスのことも気にしたほうが良いかもしれないと思ってます。 どうしましょうか。MITかGPLかと思いますが、GPLにするほど大したものじゃないですが悩ましい…

流砂の話とか

R

fortranとpythonですが、calSectionProfileとinitializeを呼んでます。いろいろ試していると、fortranの共有部分でallocateしたものが、ずっと残ってしまうので、計算が終わるごとにdeallocateしたほうが良さそうです。また、まとめて報告します。 XeonとGPU…

基礎水理シンポジウムとか

R

昨日所用で北海道に行ってました。寒すぎて生きている心地がしなかったです。 基礎水理シンポジウム 先日行われた基礎水理シンポジウムで気になった内容をまとめておこうと思っていましたが、すっかり忘れておりました。 テーマは、浮遊砂・ウォッシュロード…

Section Quasi2D model in natural river

R

できました。 とりあえず、結果だけ。いい感じです。 ptyhon+fortranです。ちょっと力作なのでコードを見てほしいです。 DBに置いておきます。