Twin Screw Simulator(Ver.6.0.0)
改良成果資料（発表用ダイジェスト版）

2018/11/15
株式会社HASLCopyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

① Template機能の改良② Job control form の新規実装③ 各種現象論的輸送方程式の陰解法化

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

1① Template機能の改良
・ Top, Flank分割数の一括更新

Division Number パネル内のUpdateボタンを押す
ことでTop, Flank分割数を一括更新
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved2

・ Kneading disk 軸方向分割数のチェック機能

Kneading disk 領域の軸方向分割数は，disk 枚数の整数倍に要請．この要請が満たされない入力情報に対し，Warningmessageを出力して分割数の訂正を促す．

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved3

・ スクリュ構成の一括置換，挿入，削除

B構成：Kneading disk ブロック
メッシュ生成コントロールファイル:B.tsmodel

A構成：標準スクリュ
メッシュ生成コントロールファイル:A.tsmodel

C構成：逆ねじスクリュ&標準スクリュ
メッシュ生成コントロールファイル:C.tsmodel

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

4新規実装された
Assemble import
ボタンをチェック状態とすると複数のtsmodelファ
イルの複合読込可能．

Importされた
tsmodelファイル
はリストボックスに登録され，置換，削除可能．

A+B+C 構成Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved5C+B+A 構成

B+C+A 構成

A+C+B 構成Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved6

② Job Control form の新規実装
新規実装されたMain menu/Tool/Job control form
プルダウン選択でポップアップ表示されるフォーム
内で複数のJobを登録可能

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

7

Job list

Job list 登録Jobを順次実行．
モデル構成，材料物性，計算コントロールデータ等を設定し，Job list に登録
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

8

Job control form でSave/Import するJobコントロールデータフ
ァイルは，Excelによる編集が可能なCSVファイル形式

スクリュ構成比較検討
材料物性比較検討

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

9

作業フォルダー名
Jobコメント
Job名=解析結果名
モデル構成
材料物性
計算コントロールフ
ァイル名③ 各種現象論的輸送方程式の陰解法化
輸送方程式

D uv
w
 f ( )
Dt t
x
y
z
表．TSSが解析対象とする現象論的輸送方程式一覧
解析モデル

ソース項 f

滞留時間

0

ひずみ

要素内一定平均ひずみ速度

繊維破断

非線形関数

液滴分裂

非線形関数

化学反応

非線形関数

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

10表．現象論的輸送方程式に対する解析法の比較
時間差分
処理

解析法

対流項処理

計算
負荷

時間刻制限

計算
回数

既往

陽解法

FVM

風上差分法

少CFL条件

大新規

陰解法

FEM

SUPG法

大無条件安定少既往：陽解法
uv
w
 f ( n )txyzn 1nnnn

新規：陰解法 n 1   ntCourant-Friedrichs-Lewy condition :

x
t u

流速（回転速） u
メッシュ幅（充満率） Δx

 n 1
 n 1
 n 1uv
w
 f ( n 1 )x
yz

時間刻Δtに制限無し

t  

当条件を外すとオーバーフロー（NAN：Not A Number 出力）定常方程式（非線形代数方程式）へ漸近
時間刻Δt

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

 n1 n1 n1uv
w f ( n1 )
x
yz

11

SUPG(Streamline Upwind/Petrov-Galerkin）法に基づくFEM離散化方程式

 n 1  nn 1
n 1n 1 n 1ˆˆuv
wd

f()d     txyz   : 要素内節点αの形状（補間）関数
標準的重み付き残差法：ˆ   

数値振動解SUPG法：

人工粘性項を付加した風上差分法と同様，解の数値振動を抑制   v  w  yz  xˆ      SUPG  u1  2 2  2 V 2  2 SUPG        ,  t   he   1 nen     V  u  v  w , he  2 uv
w
 V xyz1
222

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

121各種輸送方程式計算制御パネル内に陽解法/陰解法選択
ラジオボタンを実装

Explicit：陽解法
Implicit：陰解法

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved13陰解法を利用したひずみ解析結果

Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved

14表．解析法の比較
計算サイクル数/計算種別

計算時間

新規陰解法

定常解析（計算回数1)

1 sec未満

既往陽解法

非定常解析
(Calculation number:2000,Maximum calculation number:4000)

48 secStrain

解析法

Down stream distance (mm)
ひずみ解析結果の比較
Copyright© 2010 Hyper Advanced Simulation Laboratory Co., Ltd. All Rights Reserved15