授業も2週目が終わり、新年度切り替え業務も概ね収まるところに収まってきたので、ちょっと仕事を整理してみよう。 教育活動 某数学科目の自学自習問題の更新 某工学のファシリテーション方法の検討（受講者4名なので柔軟に） 某実験の3Dプリンタ関連の資料 …

形状最適化にも様々なものが存在する（位相最適化、寸法最適化）が、ちょっと整理しておく。 形状最適化に当たって、以下の情報を求めることになる。 設計変数(Design parameters)：寸法、位置、境界条件.... 制約条件(Constraints)：設計変数の範囲、目的関…

// read patches for torque calculation from dynamicMeshDict IOdictionary dynamicMeshDict ( IOobject ( "dynamicMeshDict", runTime.constant(), mesh, IOobject::MUST_READ_IF_MODIFIED, IOobject::NO_WRITE ) ); // access subdict of dynamicMeshDic…

インストールがうまく行ったことを確認した上で、付属のTutorialの実行をしてみた。 cd tutorials/cartesian2DMesh/hatOctree FMSToSurface geom.fms geom.stl paraview geom.stl ./Allrun cd ../../cartesianMesh/bunnyOctree paraview bunnyWrapped.stl ./…

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, …

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

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

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

昨晩から調べてて、いまさら気がついたことをメモ (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…

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

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

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

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…

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

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

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

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

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

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

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

午前中4コマ分の実験（実質3コマぐらいだけど）の後、2章を読み始める。合間に事務雑用。今年も職員の懇親会組織の集金が3万円弱。どうすっかねー。 健康診断の結果が戻ってきた。今年から中年メニューです。安定の体重過多（180cm/100kg）で引っかかったの…

先週末に届いたWind Turbine Technology - Principles and Design -を読み始める。 主に水平軸風車の設計などに書かれた本らしく、インターネットで比較的好評だったのでアマゾンで購入してみた。今日読んだ1章はさらっと前置きが書かれているだけ。2章にレ…