heatFluxの計算

勉強会で質問があって、ヌセルト数の計算について調べています。その前に、そもそもheatFluxの計算はと思い、ソースを眺めていたところ、OF23x付属のwallHeatfluxユーティリティは、どうも圧縮性を想定しているようで、thermoPhysicalPropertiesファイルを…

2014-11-27

CFDEM（OF+LIGGGHTS）その３

TUTORIAL(CFDEM) 以下のチュートリアルが付属している（OFのケースファイル）。 cfdemTestTUTコマンドで連続実行ができる...が、うまく動いている感じがしない。途中でevinceが起動したりするが、ファイルが見つかりませんというエラーになっている。 (11/28…

2014-11-26

CFDEM（OF+LIGGGHTS）その２

研究 OpenCAE

LIGGGHTSのビルド Installing LIGGGHTS on Ubuntu 12.x · CFDEMproject/LIGGGHTS-PUBLIC Wiki · GitHubを参考に、 serial版 cd LIGGGHTS_CFDEM/LIGGGHTS-PUBLIC/src cd STUBS make cd .. make serial ./lmp_serialLIGGGHTS (Version LIGGGHTS-PUBLIC 3.0.6, …

2014-11-25

CFDEM（OF+LIGGGHTS）

研究 OpenCAE 教育

ダウンロードここからインストラクションをダウンロードする。要登録。 http://www.cfdem.com/download インストール＆テストインストラクションに従って、gitレポジトリから、LIGGGHTSとCFDcouplingをgit cloneしてくる。環境変数は下記の通り、~/.bashr…

2014-11-25

OpenFOAM+DEM

教育研究 OpenCAE

しばらくぶりに調べてみる。１）OF本体に含まれるDEMソルバー DEM（個別要素法）のソルバー自体はOF2.0から含まれていて、lagrangianソルバーとして、幾つかのチュートリアルもすでにある。ソルバーとしては、以下の通り。 [waku@ensis10 lagrangian]$ pwd…

2014-11-16

OpenFOAMと外部プログラムの連携方法

OpenCAE 研究

OpenFOAM＋外部プログラムの連携方法について調査中。連携の形には、いくつかの形式がある。外部ソルバーから境界条件などを設定し、OFのケースを計算する。 OFケースの一部として例えばタイムステップ毎に外部プログラムを呼び出す。前者の例としては、D…

2014-11-01

OpenFOAMメモ

研究 OpenCAE

昨晩から調べてて、いまさら気がついたことをメモ (1) blockMeshDictで直接cellZone名を指定できる。 hex ( 0 1 2 3 4 5 6 7) rotor (20 40 1) simpleGrading ( 1 1 1) (2) mergeMeshesは、一度に２つしかマージできない。たとえば、stator, rotor1, rotor2…

2014-10-28

DAKOTA 6.0(その４：最適化手法のためのキーワード）

形状最適化研究

キーワード設計変数（ベクトル）：入力となる変数、寸法、座標、回転数、流量など制約条件：設計変数空間での制約条件（しばしば物理的な制約や、実際の製作時の制約条件）目的関数：設計変数に対する性能（出力）関数。複数の場合もあり。応答曲面：設計変…

2014-10-27

必ず1時間で終わらせる会議術

教育生活

必ず1時間で終わらせる会議術を自分の出ている会議にあてはめてみる | ICHINOHE Blog なるほどと思ったのでメモ。ちょっとだけ当方の現状に合わせてアレンジアジェンダは前日（できれば24時間以上前）のうちに関係者にメールする議題ごとに時間を割り振る…

2014-10-24

dakota 6.0（その３）

Joel Guerrero's Flubio Home Page ここにあるNACA4の形状最適化問題（設計変数は最大キャンバーと最大キャンバー位置（翼厚比一定）、目的関数：L/D最大化）は、dakota 5.4＋OF2.3＋openscadで計算されている。今後は、6.0系を使いたので、dakota_of.inを…

2014-10-24

久しぶりに何もない１日

研究形状最適化

プロコンやらイベントが続いて１０月中旬は何も出来ず...。・科研費申請書・某大学の共同研究計画書（こちらへの配分は０円だけど）を提出したので、さしあたりの締切ありの書類仕事もない。現状、タスクリストに並んでいるのは、・共同研究先から依頼…

2014-10-11

dakota 6.0（その２：DAKOTAとOFの連携TIPS）

OpenCAE 研究形状最適化

調べたところ、Dakota6.0/5.4/5.3のライブラリlibsampling.soと、OFのライブラリの名前が重複しているので、Dakotqのほうのライブラリ名をちょっと変えたものをソースからインストールする。１）準備(CentOS6.5) cmake 2.8.4以上(5.4)/2.8.8以上(6.0)→標準r…

2014-10-08

dakota 6.0（その１）

まずはユーザーマニュアルに従ってRosenblock関数の最小値の計算問題でテスト。 Rosenblock関数は、最適化問題のテストとしてよく出てくる。目的関数：f(x1,x2)=100*(x2-x1^2)^2+(1-x1)^2 設計変数：x1,x2 制約条件：|x1|≦2,|x2|≦2 ※最大値:f(-2,-2)=3609, …

2014-10-08

dakotaのインストールと6.0/5.4の違い

OpenCAE 研究形状最適化

dakotaのインストール SNLのサイトからバイナリパッケージをDL 必要なOS側パッケージをインストール yum install boost boost-system boost-signals boost-regex boost-filesystem blas lapack openmpi openmotif libX11 libXau libXext libXmu libXp libXpm…

2014-10-02

sedで特定の複数行を複数行で置換

研究

※ もっと賢いやり方があるはずです（確信）。こいつ↓を、 AMIr { type patch; nFaces 200; startFace 131854; } こうしたい AMIr { type cyclicAMI; inGroups 1(cyclicAMI); matchTolerance 0.0001; transform noOrdering; neighbourPatch AMIs; nFaces 200…

2014-10-01

python でコマンドライン引数をとり数値化する

python 研究

python でコマンドライン引数をとり、数値化(float)する。 argv[0]はスクリプト名（****.py）が入るので、argcは引数の数より１つ多い。引数は基本的にstringのため変換が必要。 import sys # for get argv argvs = sys.argvargc = len(argvs)print if (arg…

2014-09-30

指導方法について

教育研究

夏からいろいろ考えていたが、「研究する研究室」になれるよう、学生をどんどん”挑発する”ことにした。学生をあげつらったり、からかったりするのではなくて、どんどん質問するということです。「この寸法はどうやって計算したのか？」「この文献の意義は？…

2014-09-23

topoSetコマンド

OpenCAE 研究

OFでcellやfaceの集合を作るときは、setSetコマンドもしくはtopoSetコマンドを使用する。後者はsystem/topoSetDictで指定できるので、わかりやすいかも。 /*--------------------------------*- C++ -*----------------------------------*\| ========= | ||…

2014-09-16

筋肉痛の週明け

生活

今年は人手がたりないということで、嫁実家の稲刈りの手伝いを1日半ほど。コンバインではなく、歩きの稲刈り機で出てきた稲束をハサガケしていました。田んぼのぬかるみで踏ん張るため、大腿部内側が超絶筋肉痛です。よい運動になりました。今年は2枚休耕し…

2014-09-13

pythonでblockMeshDict

OpenCAE 研究

すごい力技だが、pythonの勉強もかねてblockMeshDictを出力するスクリプトを、ここ1週間ほどで書いてみた。 2次元のVAWTで、翼型はNACA4シリーズを出力できるようにしている。翼枚数は3枚前提だったり、はっきり行って汎用的ではないが、これをベースに改良…

2014-09-03

OF+外部ライブラリ

OpenCAE 研究

OFの計算で冷媒ガスの熱物性値を参照したいと考え調査してみた。事例としては、・OF+FreeSteam (水蒸気：IAPWSIF97) があっただけで、NIST REFPROP（最新版9.1）や、同様の機能を持つOSS版のCoolProp（最新版4.2.5?5系が開発中？？）とOFとの連携事例は、…

2014-08-31

1日体験入学

教育研究 OpenCAE

昨日、今日と学校行事の中学生向け体験入学を実施。午前中は説明役で、13：00～15:20までの自由見学は学生さんにトークをおまかせする。こういう世代の違う子供に話す力を養うのも教育の一環と思う。午前中、6回同じ内容をしゃべったせいか結構疲労感があ…

2014-08-25

8月の個人的イベント諸々が終了

生活

昨日の妹の結婚式をもって、8月に予定していた個人的イベントが終了したので、やれやれといったところ。二人ともお幸せに。 9月は職場関係の業務が多いのでメモメモ・今週末の体験入学・学会１＠仙台・報告会＠八王子・部活合宿・学会２＠仙台・エコラ…

2014-08-05

ダクト内間欠流の計算3(groovyBC+DirectEulerAngleRotatingMotion)

OpenCAE

入口・出口の圧力差と拙作DirectEulerAngleRotatingMotion(トルク＞回転数＞回転角を計算)の組み合わせてテストしてみた。うまく行っているかな。 ※ 画像は、はてなのアップロード容量を超えてしまったので、youtubeから。月たったの30MBというのはいかがな…

2014-08-01

ダクト内間欠流の計算2

OpenCAE

昨日の続きです。 inlet { type groovyBC; variables "p0=0.5;"; // =(0.5*U0*U0)/4 valueExpression "p0*(pow(sin(pi*time()/50),7))"; value uniform 0; } outlet { type groovyBC; variables "p0=0.5;"; // =(0.5*U0*U0)/4 valueExpression "p0*(-pow(sin…

2014-07-31

ダクト内間欠流の計算１

OpenCAE

とある用途のため、OpenFOAM(+swak4foam)で、ダクト内流れが間欠的に流れの向きが変わる計算をしているが、なかなかうまくいかない。現状の、境界条件は以下の通り。 0/Ｕ inlet { type pressureInletVelocity; value uniform (1 0 0); } outlet { type pres…