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

嫁さんがご飯を食べられないんだが、どうしたものか。 スイカだけどか、カブトムシのような食生活である。うーむ。

病院で確認したところ、7週2cmくらいだそうです。頑張れ。

最終的にやりたいことは「OpenFOAM+OpenMDAO」で流体計算による形状最適化なんだが、pythonでSTLのReader/Writerを作って（性格にはネットからの「かき集め」）みた。 もう1回流れを整理すると、以下の通り。 ・OpenMDAOの導入(done) ・FreeCADもしくはPytho…

最終的には、Linux環境でやるつもりだけれど、とりあえずWindows版でいろいろ確認。 環境 Windows 7 x86_64 + FreeCAD 0.13 32bit (Python 2.6.2 同梱) 準備 FreeCADを起動し、pythonコンソールを表示させる（表示されない場合は、表示 → ビュー → Pythonコ…

前のバージョンアップで、タムロンの望遠が動かなくなったのでガッカリしていたが、今回のL1.20で、再度動くようになったとの情報が！さっそくファームアップした。 今度の引越し先近くの川には、かわせみがいるそうなので、日曜日の朝はカメラもって出撃す…

式は挙げずに、入籍だけ。親戚同士の顔合わせはやりますが....。 こんな自分でも結婚してくれる人がいたことに感謝です。 これからは健康にもっと気をつけます。

注意： salome_mecaを既にインストール・設定ができていることが前提の記事です。 SL6.5 x86_64です。 OpenMPI, BLAS, libxml2などのシステムパッケージはインストール済みとします。（基本的にsalome_mecaに同梱のライブラリを使用しますが、libxml2はシス…

REFPROP 9.1を購入したので、ORCのシステム設計を開始。来年度前半くらいまでにシステムが組めると良いな。 まず、基本ランキンでタービン効率・発電機効率、ポンプ効率を考えたサイクルシミュレーションを実施してみた。媒体は水蒸気とR245FAで、文献値（エ…

計算をリスタートする時の処理がうまく行ってなかったんだが、shスクリプトとsed,awkをあれこれして何とかクリア。次は、並列計算の対応。 今回の場合、外部連携プログラム（Fortran）を動かすのに、OFに含まれるsystemCallを用いて起動shスクリプトを叩いて…

非定常入口流速境界条件の場合です。 慣性モーメントが小さいとレスポンス良く回ってくれます。当たり前の結果かもしれませんが。ただし、ローターが重いと安定して（風が落ちても回転数を維持して）回ってくれるという結果にもなっています。 非定常な境界…

昨日の動画をyoutubeにアップしてみました。ちょっと出力が荒すぎた感じがする...。 次のステップとして、入口の流速境界条件をswak4foamで非定常にしてみよう。 ここまで出来たら、とりあえず何処かに公開してみようかな。 2D Crossflow - YouTube

時間ステップを可変にすると10秒後ぐらいに発散するので（理由は不明）、時間ステップ一定で計算したところ計算はうまく進行した。 静止状態から無負荷でトルク0になるまでの回転数の変化が得られた（図1、図2）。 最大回転数は200rpmくらいで、従来のdynami…

~/OpenFOAM-2.2.x/src/dynamicFvMesh/solidBodyMotionFvMesh/solidBodyMotionFunctions/rotatingMotion/rotatingMotion.C(.H)を参考にライブラリを作ってみた。 オリジナルは、時刻tとdynamicMeshDictから読み込んだ角速度、回転中心を用いて、 // Rotation …

OpenFOAMでpimpleDyMFoamなんかの場合、回転数は固定になるため、例えば風車や水車のような動的回転数変化を考えることはできない。 回転メッシュのライブラリを見るとeulerAngleを計算する場合、時間×回転速度になっているようで、外部のプログラムと連携し…

スタンフォード大で出しているオープンソースCFDソフトウェアのSU２というものを使ってみる。 （１）ダウンロード http://su2.stanford.edu/download.html 現在のバージョンは3.0で、フルバージョンであるSU2 SUITEと、教育用に一部簡素化したSU２-Edu SUITE…

なにかしら１文書くことにします。

梅雨末期的様子ですが、セミが本格的に鳴き出し、昼はミンミンゼミ、夕方から夜にかけてはヒグラシが煩くなってきました。

はてなTシャツ2013ができました！ プレゼントキャンペーンを実施します - はてな広報ブログ あたるかな～。

OF1.6-extのgit版が動いているとして、 cd ~/OpenFOAM/OpenFOAM-1.6-extgit checkout -b feature/solidMechanics remotes/origin/feature/solidMechanics "application/solver/solidMechanics" にソルバーが追加されているので、./Allwmakeする。 IcoFSINonL…

大体こんなもの、やはり夜が問題ですな。 朝： ・おにぎり1個 ～200kcal ・サンドイッチ1個～300kcal程度のもの 昼： ・トマトジュース１パック（200ml)～40kcal ・おにぎり1個～200kcal ・サラダとか惣菜1品～300kcal 夜： ・All free シトラステイスト～0k…

こそこそやっていたDAKOTA+OFのサンプルができたので、とりあえず１段落 次のステップとしては、 ・GridgenやSWとの連携によるメッシュ生成の高度化（blockMeshでは限界） ・最適化理論をちゃんと勉強する といったところか。

SL6.3/VMware PlayerにDAKOTAをインストール http://dakota.sandia.gov/distributions/dakota/5.3/download.html - DAKOTA5.3 - JAGUAR2.1 ①準備 # yum install lesstif blass ②インストール > wget http://dakota.sandia.gov/distributions/dakota/5.3/dako…

新年度なので、家マシンでやりたいことを挙げておく。 OF＋DAKOTAで形状最適化 OFソルバー開発（PimpleDyMFoamの改良、interFSIFoam） ENSIS復活

試験とかの対応で日があいてしまった。研究生活のリハビリで更新されたSALOME-MECAをSL6.3へインストールしてみる。 SALOME-MECA＝Salome(モデラ)＋Code-Aster(構造解析ソルバ)なんですかね。 １．とりあえずダウンロード http://www.code-aster.org/V2/spip…

（昨日の続き） pgmon.outが作成されたので、pgprofで確認してみた。(pgprof -exe a.out) 計算時間がかかっているものから並べると main 92% is_raycross 2% geometres 2% rot_vec 1% check_face 1% : ということで、予想通りmainルーチンの中の計算時間が圧…

課題と目標 以前作った形態係数計算の計算プログラム Total number of nodes : 60734 Total number of boundary faces: 25704 Total number of elemnts : 326559のデータに対して、計算時間内訳は、現状で以下の通り [0:geometry] viewfactor preparation co…

@waku2011: そうか！安い所に外注すればいいんだ！！年収数千万円の会社員、個人の全業務を中国に外注していたのがバレて捕まる : ギズモード・ジャパン URL @gizmodojapanさんから2013-01-25 17:20:58 via Tweet Button @waku2011: 参考になる記事です。岩…

@waku2011: CUDA5では、CUDA_SAFE_CALL関数はいずこへ...。2013-01-24 12:13:18 via web