いわて駐在研究日誌2。

NEVER STAND BEHIND ME

includeAngleとresolveFeatureAngleについて

いまさらながらsHMねたです。 STLで出力した図形のresolveFeatureAngleをうまく設定するには、paraviewのFeature edges filterを使うと、角度毎に表示させて確認ができる。この値をsHMDictで指定すれば良い。 surfaceFeatureExtractDictに指定するincludeAng…

あらかじめ境界メッシュの動きが分かっている場合のdynamicMesh

境界パッチの動きが時間や空間の関数として分かっている場合にどうするかを調べたい。 境界パッチ形状の移動に伴う流体メッシュ変形については、例によってPenguintsさんのサイトにまとめられていた。ただし、パッチ形状が変形するわけではなく、あくまで剛…

mesh tutorial (OF4.x)

ということで、これまでmeshの剛体回転ばかりやってきたので、少し目先を変えてmesh morphing関係の調査。まずはチュートリアルから。 ① tutorials/incompressible/pimpleDyMFoam/movingCone 円筒領域(くさび型で2D化)の中をconeが平行移動する。メッシ…

OF-devのビルドエラー

学生に教える都合もあって、基本的に4.xの系列しか入れてないのだけど、岩手大の学生さんからの相談対応で、meshMorphing機能とかIB法関連の実装を調べるため、devとfoam-extend4.0を入れてみる。 extend4.0は問題なかったが、devがビルドエラーになり、/src…

新年度です。

結局、動的ライブラリが放置状態なので今年度こそはと思う。 今年度の目標は、 校務 授業は確実に準備をすすめること(前倒しでやること) 卒研・特研は学生との作業時間を増やすこと 地域連携を進め、共同研究3件を目指す(1件決定済み) 研究 マイクロO…

動的なsolidBodyMotionFvMeshのライブラリ作成(その2)

基礎にもどって、dynamicMeshの調査。 OF4.xでのdynamicメッシュは大きく分けると、 dynamicFvMesh ・・・Automatic mesh motion topoChangerFvMesh・・・Topological mesh changes に大別される(Ref.: A.O.Gonzalez,Mesh motion alternatives in OpenFOAM,…

動的なsolidBodyMotionFvMeshのライブラリ作成

OF、あるいは、動的メッシュ機能をもつたいていのCFDソフトは、メッシュの移動・回転速度を指定して計算するのが一般的だろう。 風車の研究をしているので、ブレードトルク→回転の運動方程式→角速度→回転角のような動的なメッシュ移動ができないか試行…

extBlockMesh for OF3.0 and higher

extBlockMeshをOF3.0以降で利用する場合はこちらから。 GitHub - HIKassem/extBlockMesh: extBlockMesh fork for OpenFOAM-dev version

SL6にdevtoolを入れる

SL6やCentOS6系はコンパイラのバージョンが古いため、OF4.xのビルドなどのためには、いろいろ支障が出てきている。そこで、devtoolと呼ばれるコンパイラ関係の更新パッケージをいれてみる(SL6.8)。 linux.web.cern.ch devtoolは現在4まであるみたいだが、…

chtMulitRegionFoam + DynamicMesh(fluidRegion) + interFoam + Evapolation/Condensation on OF3.0.x (その1)

というのをやりたいので、以下の順番に進める。 OF4.0が出たばっかりだが、ちょっとヘッダファイルが変わったようなので、とりあえず3.0.xで進める。 ①chtMulitRegionFoam + DynamicMesh(fluidRegion) = chtMultiReginDyMFoam

Immersed boundary method on foam-extend 3.2

(注意) foam-extend 3..2はgcc 4.9.2以上が必要なようです。 OpenFOAMのforkしたバージョンであるfoam-extend 3.2が昨年9月にリリースされ、前から気になっていたImmersed boundary methodが正式にリリースされていたので、知らべてみた。現在は、nextRele…

Adaptive Mesh Refinement(AMR) on OpenFOAM 3.0.x

AMRはDyMソルバーで利用可能。3Dメッシュのみ+Parallel計算可能。 interDyMFoamのチュートリアルに含まれているconstant/dynamicMeshDictでは以下の通りに記述されている。 dynamicFvMesh dynamicRefineFvMesh;dynamicRefineFvMeshCoeffs{ // How often to …

助成金申請

なかなか科研費があたらないので、助成金のほうで頑張って稼ごうと足掻いています。 本日2件提出しました。来月も2件提出予定です。 現状、 ・企業さんの申請に協力する形で2件(@30弱) ・某M社、E社の助成金(@100, @65) を提出して結果待ちで、来月〆…

Dakotaの並列計算機能について(その5)

Dakotaの並列計算機能について整理してみる。 計算機環境で整理すると、 1)マルチコアワークステーション1ノード向け dakotaのinstanceは1つだけ起動(シリアル) 例えばasynchrnous evaluation_concurrency 4 (4つ同時にeval.)を指定する。 1個1個のeva…

平成28年度開始しました。

年度初めのドタバタを終了し、なんとか毎日のリズムが出てきた。 昨年度は練習問題を一生懸命やっていたなぁというのが正直な所で、今年度はちゃんとアウトプットしていきたい。 研究リハビリ第1段として、昨秋にやっていたDakotaのRosenbrock関数(2D)の最…

OpenCV3.1.0+VS2015+Win10Pro

前から気になっていたOpenCVをいじってみようということで、とりあえずWebカメラのついている家PC(Win10Pro x86_64)に環境構築をトライしてみる。OpenCV3.1.0が年末に公開され、VS2015にも対応したとのことなので、C++の勉強もかねてVS2015 C++での開発を考…

blockMeshのmultiGrading 機能(2.3.x~)

blockMeshDGなんていうツールもあるが、ちょっと調べたら2.3.xから本家のblockMeshコマンドにmultiGrading機能が追加されていたとのことで調べてみた。 従来のblockMeshでは例えば矩形ダクト流れを計算するときに壁面にメッシュを寄せたい場合は、ブロックを…

Dakotaの並列計算機能について(その4)

ちょっと時間が開いてしまった。 4ノードの「なんちゃってクラスタ」でdakota+OF並列計算したいときのやり方をまとめておく。 クラスタでの利用にはいくつか(計算したい内容に応じて)方法があるが、ここでは、 個々の目的関数の評価はOpenFOAMの並列計算…

Dakotaの並列計算機能について(その3)

MOGAがやっぱり時間がかかるので、近似応答局面を導入した計算(Surrogate based optimization)の効率化を図ってみる。 テスト問題はBranin function(2D)で、3つの最小点がある。 f(-pi, 12.275)=f(pi,2.275)=f(9.42478,2.475)=0.397887 moga部分は共通で、以…

Dakotaの並列計算機能について(その2)

以下のような簡単な問題(Beale's function)を最適化で解いてみる。設計変数がx1, x2で、目的関数(単一)がf、制約条件(単一)がc>=0である。最小値の答えはf(3, 0.5)=0。 See: Test functions for optimization - Wikipedia, the free encyclopedia #!/bi…

Dakotaの並列計算機能について

ようやく復帰。Dakotaの並列計算機能について調べたのでまとめてみる。 概観と簡単なサマリー Dakota自体から起動する"simulator"はシリアル計算でも並列計算でももちろん構わない(たとえばOFの並列計算など)ので、ここではDakota本体の並列計算について調べ…

最近の状況

子育て中心で、なかなかブログを更新できておりません。 定時を過ぎたらばたばたと迎えにいって、おむつ替えて、遊びの相手して、風呂入れて、ミルク飲ませて、寝かしつける、で8時過ぎ。 飯を食って風呂入ると、疲れて眠くなり、机で気絶してたりする。少…

パスフレーズなしのSSHログイン(OpenMPI)と、NFSマウントの設定

パスフレーズなしのSSHログイン設定メモ (1)クライアント(ログインユーザはwaku)でssh2のRSA暗号鍵を生成 $ ssh-keygen -t rsa ※ パスフレーズを聞いてくるが、なしでよければそのままリターンでOK (2)公開鍵(*.pub)のほうを接続したいサーバ…

OF+外部ライブラリ(物性値ライブラリ)その3

水・水蒸気のライブラリであるfreesteamとOFの併用テスト freesteam 2.1: http://freesteam.sourceforge.net/ IAPWS-IF97-OF:https://github.com/romansCode/IAPWS-IF97-OF (freesteam 2.0~, OF 2.2~) 1)まず、freesteamをビルドして(sconsが必要)、.s…

VAWT翼型最適化(2)

メッシャーの問題がひと段落したので、いよいよ計算という感じであるが、ここから泥臭い作業が続くのだろうなと覚悟しておく。 その前に、Dakotaで計算する場合、特定のCFD計算が発散したりしてエラーになる場合、それを飛ばして次の計算に移る仕組みを調…

cfMeshで機能追加を希望するところ

cfMeshをいろいろ使っていて、やっぱり非常に使いやすい。 ちょっと不具合があってtwitterでぼやいていたら、わざわざ中の人が声を掛けてくれて直してくれたし(正確には、自分のSTLデータの造り方が想定外だったようだけど)、非常に好感触である。 さらに…

VAWT翼型最適化(1)

子育てしてると時間の感覚がなくなってしまっていかんな。もう6月半ばになってしまった。7月上旬に岩手大で発表するのに間に合わせるために、ちゃんとやらねば。 VAWT用2次元翼型の形状最適化を行いたいのだが、メッシュの自動生成でblockMeshではなかなか…

cfMesh v1.0.1 my tutorial

いろいろメモ(cartesian2DMesh関係) fms形式のジオメトリデータをstlに変換 FMSToSurface geom.fms geom.stl stlデータをfms形式に変換(同じファイル名) surfaceToFMS geom.stl 境界条件用に、パッチ毎のstl solidエンティティとする(1つのstlファイル…

C++のダイナミックライブラリの作成・使用

$ more hello.cpp #include<iostream>#include<stdio.h>using namespace std;int hello(){ cout<<"Hello,World!"<</stdio.h></iostream>

OF+外部ライブラリ(物性値ライブラリ)その2

去年の9月に調べていたOF+熱物性値ライブラリ関連を放置していたけど、μORCシステムを真面目に設計するため、再調査してみる。 使えそうなライブラリ freesteam 2.1: 国際標準IAPWSを利用した水・水蒸気熱物性値ライブラリ(C/C++) http://freesteam.source…