Iform = 6
ブロックフォーマットキーワード この境界材料は、多相材料則/MAT/LAW51(Iform = 0、1、10または11)についてサイレント流出境界のシミュレーションを可能にします。本材料は、無限領域のモデル化に有益です。この材料則を割り当てられた要素はサイレント要素境界のように機能します。
これは、Iform = 3と使用される定式化を改善したもので、波形の反射と逆流の問題を最小にするのに役立ちます。計算領域の外側限界上の波形反射を回避するための音響処理および体積比率用の“連続”モデルの両方が含まれます。
フォーマット
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
---|---|---|---|---|---|---|---|---|---|
/MAT/LAW51/mat_ID/unit_ID | |||||||||
mat_title | |||||||||
空白行 | |||||||||
Iform | |||||||||
Pext | Tcp | Tca |
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
---|---|---|---|---|---|---|---|---|---|
空白行 |
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
---|---|---|---|---|---|---|---|---|---|
空白行 |
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
---|---|---|---|---|---|---|---|---|---|
空白行 |
定義
フィールド | 内容 | SI 単位の例 |
---|---|---|
mat_ID | 材料識別子 (整数、最大10桁) |
|
unit_ID | 単位識別子. (整数、最大10桁) |
|
mat_title | 材料のタイトル (文字、最大100文字) |
|
Iform |
(整数) |
|
Pext | 外圧 デフォルトはStarterによって計算されます(実数) |
|
Tcp | 圧力緩和のための特性時間 デフォルトはStarterによって計算されます(実数) |
|
Tca | 体積比率緩和のための特性時間 デフォルトはStarterによって計算されます(実数) |
|
副材料#j=1...3の初期体積比率 デフォルトはStarterによって計算されます(実数) |
||
副材料#j=1...3の初期密度 デフォルトはStarterによって計算されます(実数) |
||
副材料#j=1...3の初期体積エネルギー デフォルトはStarterによって計算されます(実数) |
||
副材料#j=1...3の圧力カットオフ デフォルトはStarterによって計算されます(実数) |
||
副材料#j=1...3の初期圧力 デフォルトはStarterによって計算されます(実数) |
||
副材料#j=1...3の初期音速 デフォルトはStarterによって計算されます(実数) |
例
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/LAW51/1
Boundary Element, AUTO-OUTLET
# IFORM
6
#---Global parameters------------------------------------------------------------------------------#
# PEXT TCP TC_ALPHA
0 0 0
#---Material#1:------------------------------------------------------------------------------------#
# ALPHA_1 RHO_0_1 E_0_1 P_MIN_1 P_0_1
0 0 0 0 0
# SSP_1
0
#---Material#2:------------------------------------------------------------------------------------#
# ALPHA_2 RHO_0_2 E_0_2 P_MIN_2 P_0_2
0 0 0 0 0
# SSP_2
0
#---Material#3:------------------------------------------------------------------------------------#
# ALPHA_3 RHO_0_3 E_0_3 P_MIN_3 P_0_3
0 0 0 0 0
# SSP_3
0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
コメント
- この圧力定式化は、Bayliss and Turkelによる方程式から導出されています。 8
(1) ここで、- 音響インピーダンス
- 流束が領域の内側からである、すなわち法線速度 <0(停滞圧力が代わりにかけられる)でない限り入力として与えられる遠距離場圧力。
- 材料の移流は次のように定義されます:
- 流出流束:
局所変数は単に更新されます。- 流入流束:
これは領域内で起こっていることが遠距離場でも起こると仮定しているため、局所変数は上述のケースのようには更新されません。しかしながら、流れが逆流する以前にそうであった状態と材料が局所的に類似していると考えるのは理にかなっています。緩和係数が追加され、流入材料が入力で与えられている設定された値へと次第に戻っていくことが可能となります。
- 定義されていない場合、全体パラメータPextと初期状態は、隣接要素から自動的に計算されます。
- 爆発副材料については、初期状態およびデータは常に、隣接要素からのデータを用いて自動的に計算されます。
- 各特性時間は、1桁の、圧力波(Tcp)または材料(Tca)が領域の最も近い限界に達するために必要な最短時間として計算されます。非圧縮性の場合、TcpとTcaは等しいとすることが可能です。爆轟または衝撃波のシミュレーションの場合、Tcpは1桁の波の伝播速度を用いて計算されます。
- 各流出要素は、ALE領域からの1つの要素と1つの面を共有することが推奨されます。
- 流出境界材料は、境界に達する波は局所的に平面波または球状波に近似されるべきことを意味する音響単極でできています。したがって、遠距離場の減少を正しくモデル化するには、流出要素は、爆発または衝突領域から、爆発または衝突領域の1倍から10倍の大きさ分のところに位置している必要があります。
- Alvin Bayliss and Eli Turkel, Outflow boundary condition for fluid dynamics, Report 80-21, Institute for Computer Applications in Science and Engineering, NASA Langley Research Center, Hampton, Virginia