createPatchによる2方向周期境界作成

[Tktk]

こんにちは. 大学でOpenFOAMによる流体解析をしているものです. お世話になります. 

当方, メッシュの作成にサードパーティツールを使用しているため, 周期境界をcreatePatchを用いて作成しようとしています. 

解析ケースはチャネル流れ, X（流れ）方向, Z（スパン）方向に周期境界を置きたいと考えています. 

おのおのでcyclicとなるようcreatePatchDictを作成し, 実行したところエラーとなり境界の作成ができません. 

以下エラー内容です. 

--> FOAM Warning :
From function syncPoints(const polyMesh&, pointField&, const CombineOp&, const point&)
in file createPatch.C at line 473
There are decomposed cyclics in this mesh with transformations.
This is not supported. The result will be incorrect
Points changed by average:0.3793571 max:1.414214



--> FOAM FATAL ERROR:
face 0 area does not match neighbour by 76.3932% -- possible face ordering problem.
patch:side0_0 my area:0.25 neighbour area:0.559017 matching tolerance:0.0001
Mesh face:20 fc:(0.25 0.75 0)
Neighbour fc:(0.25 0.75 0.5)
If you are certain your matching is correct you can increase the 'matchTolerance' setting in the patch dictionary in the boundary file.
Rerun with cyclic debug flag set for more information.

From function cyclicPolyPatch::calcTransforms()
in file meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C at line 221.

FOAM exiting

質問として, 以下の2点があります. 

・cfd-onlineを調べてみると, createPatchは複数方向の周期境界作成に対応していないというような内容の投稿が見つかりました. これは本当なのでしょうか. 

・メッシュ形状が複雑でメッシャに頼らざるを得ない状況です. 他に周期境界を置く方法はないでしょうか. 

以上何卒アドバイスいただけると幸いです. よろしくお願いいたします. 

[ohbuchi]

こんにちは。
メッシュから抽出したパッチ面が正しく周期性を満たしていないためのエラーの様です。
まずは、形状の正確さを検証した方が良いでしょう。

複数パッチを生成したいなら二度に分ければ良いのではないでしょうか？

cyclicではなくcyclicAMIを使った方が簡単な気がします。



2015年9月15日火曜日 15時33分46秒 UTC+9 Tktk:

[Tktk]

こんにちは. 返信, アドバイス有難うございます. 

ご指摘のところ確認してみたところ, 周期境界を置けたのでメモさせてください. 

流れ方向はcyclicとし, スパン方向にcyclicAMIを用いることで, 1度のcreatePatchで現象を解決できました. 

また, 流れ, スパンともcyclicとして2回に分けてcreatePatchを実行しても境界が置けないことを確認しました. 

このとき, 片一方がcyclicAMIならば境界を置けるか, ということは確認していません. 

投稿いたしましたエラー内容については, createPatchが2つ以上の, "cyclicを置く2つのパッチの法線ベクトル"を同時に扱うことができないことが原因, と推察しています. 

その他は未確認です. 

貴重なご意見有難うございました. 

2015年9月28日月曜日 14時22分09秒 UTC+9 ohbuchi: