高次精度化&特性変数変換 柴山拓也 (名大、STEL) 2015年8月7日.

Slides:



Advertisements
Similar presentations
1 高速フーリエ変換 (fast Fourier transform). 2 高速フーリエ変換とは? – 簡単に言うとフーリエ変換を効率よく計算 する方法 – アルゴリズムの設計技法は分割統治法に基 づいている 今回の目的は? – 多項式の積を求める問題を取り上げ、高速 フーリエ変換のアルゴリズムを用いた解法.
Advertisements

Division of Process Control & Process Systems Engineering Department of Chemical Engineering, Kyoto University
1. 補間多項式 n 次の多項式 とは、単項式 の 線形結合の事である。 Definitions: ( 区間, 連続関数, abscissas (データ点、格子点、差分点), 多項 式 ) Theorem. (補間多項式の存在と一意性) 各 i = 0, …, n について、 をみたす、次数が高々 n.
1章 モーションプログラムの仕組みを理解する.
データ構造と アルゴリズム 第十二回 知能情報学部 知能情報学科 新田直也.
4章 制御の流れ-3.
FB6280E SSP パッケージ A3 BookEdge Scanner (NXScan SSP edition 添付版)
三重対角化アルゴリズムの性能評価 早戸拓也・廣田悠輔.
動機 目的 音声特徴比較システムの開発 結果を考察 サンプルデータ の収集及び統計 班員全員が音楽好き 歌の上手さを科学的に理解
24 両端単純支持梁に対する外乱抑制制御系の製作
プロセス制御工学 3.伝達関数と過渡応答 京都大学  加納 学.
スペクトル法による数値計算の原理 -一次元線形・非線形移流問題の場合-
数楽(微分方程式を使おう!) ~第5章 ラプラス変換と総仕上げ~
符号化のための重み付きジョイントバイラテラルフィルタを用いた 奥行き画像超解像
「データ学習アルゴリズム」 第3章 複雑な学習モデル 3.1 関数近似モデル ….. … 3層パーセプトロン
情報のディジタル化 情報量の単位(bit) 文字 数値 アナログ情報.
HLLD法に基づく 磁気流体方程式の差分解法
デジタル信号処理④
放射線計測エレクトロニクスの信号処理の為の アナログ電子回路の基礎 第十三回
(ラプラス変換の復習) 教科書には相当する章はない
電気回路学Ⅱ エネルギーインテリジェンスコース 5セメ 山田 博仁.
シミュレーション演習 G. 総合演習 (Mathematica演習) システム創成情報工学科
京都大学大学院医学研究科 画像応用治療学・放射線腫瘍学 石原 佳知
地理情報システム論 第14回 ラスタ形式による空間的演算 GISの応用(2)~土地利用の予測
システムモデルと伝達関数 1. インパルス応答と伝達関数 キーワード : 伝達関数、インパルス応答、 ステップ応答、ランプ応答
宇宙磁気流体・プラズマシミュレーション サマーセミナー ~三次元MHDコードの作成〜
音信号表現 音声波形のデジタル化(PCM) サンプリング、標本化定理、量子化 ソースフィルタモデル
・コンピュータのアナログデータの 扱いについて ・制御
電気回路学Ⅱ エネルギーインテリジェンスコース 5セメ 山田 博仁.
電気回路学Ⅱ 通信工学コース 5セメ 山田 博仁.
行列による画像処理 デジタル表現論 担当者:劉 雪峰 2017年6月1日.
「移流輸送の 高次精度かつ高解像度スキームの開発」
プログラミング序論 2. n人のインディアン.
アルゴリズムとデータ構造勉強会 B+tree.
6. ラプラス変換.
理化学研究所 重イオン核物理研究室 馬場 秀忠
磁気回転不安定性によるブラックホール降着流の角運動量輸送機構の 解明
ルンゲクッタ法 となる微分方程式の解を数値的に解く方法.
地理情報システム論(総)/ 国民経済計算論(商)
電気回路学Ⅱ コミュニケーションネットワークコース 5セメ 山田 博仁.
ソフトウェアを最適化 高い画像品質を維持しながら、 より高速なスキャンが可能
FEM勉強会 (第3回).
Mini-RT装置における 強磁場側からの異常波入射による 電子バーンシュタイン波の励起実験
Fourier 変換 Mellin変換 演習課題
1~15までの数字の中から、 1個の数字を選び、覚えて下さい。
UNIX演習 情報ネットワーク特論資料.
宇宙磁気流体・プラズマシミュレーションサマーセミナー
地域情報学 C言語プログラミング 第2回 変数・配列、型変換、入力 2017年10月20日
地理情報システム論(総)/ 国民経済計算論(商)
電気回路学Ⅱ エネルギーインテリジェンスコース 5セメ 山田 博仁.
第6回レポート解説 条件1 条件2 条件3 月の入力 月、日、曜日の表示 日の入力 曜日の入力
オブジェクト指向言語論 第三回 知能情報学部 新田直也.
電気回路学Ⅱ コミュニケーションネットワークコース 5セメ 山田 博仁.
ディジタル信号処理 Digital Signal Processing
電気回路学Ⅱ エネルギーインテリジェンスコース 5セメ 山田 博仁.
電気回路学Ⅱ コミュニケーションネットワークコース 5セメ 山田 博仁.
電気回路学Ⅱ 通信工学コース 5セメ 山田 博仁.
数値解析 第6章.
1~15までの数字の中から、 1個の数字を選び、覚えて下さい。
情報処理Ⅱ 第7回 2004年11月16日(火).
電気回路学Ⅱ 通信工学コース 5セメ 山田 博仁.
情報処理Ⅱ 2005年11月25日(金).
Q q 情報セキュリティ 第7回:2005年5月27日(金) q q.
電気回路学Ⅱ コミュニケーションネットワークコース 5セメ 山田 博仁.
信号データの変数代入と変数参照 フィードバック制御系の定常特性 フィードバック制御系の感度特性
復習 いろいろな変数型(2) char 1バイト → 英数字1文字を入れるのにぴったり アスキーコード → 付録 int
C言語講座 四則演算  if ,  switch 制御文.
アルゴリズム ~すべてのプログラムの基礎~.
プログラミング 2 静的変数.
情報処理3 第3回目講義         担当 鶴貝 達政 12/17/2019.
Presentation transcript:

高次精度化&特性変数変換 柴山拓也 (名大、STEL) 2015年8月7日

高次精度化の弊害 X 密度 基本変数 or 保存変数に補間をすると次数を上げるほど不連続の解像度は良くなるが数値振動が発生してしまう これじゃ精度上げた意味がない!!! Brio & Wu 1988

特性変数変換 基本変数: 開始 初期設定 CFL条件の評価 SSP-RK法 基本変数→特性変数: 特性変数変換 特性変数→基本変数: 高次精度補間 数値流束の評価 特性変数逆変換 時間積分 No t = tout Yes データ出力 Yes t < tend No 終了

振動を克服!! 不連続の解像度はWENOZが一番高いが振動が取りきれてない。 美しい!! もっとわかりやすい例 (Ryu & Jones) 密度 X 不連続の解像度はWENOZが一番高いが振動が取りきれてない。 美しい!!

Koren Limiter 密度 Koren Limiterは解像度は高いが特性変数補間すると変なアンダーシュートが。 X

まとめ HLL(1st), HLLD(1st), MUSCL(MM, MC, VL, KR), WENO(3, +3, 5, Z)とその特性変数変換のコードができた。 Koren Limiterは優秀だが特性変数補間をすると変なアンダーシュート Alfven波の減衰率を比較してみたかった。 オタクへの道を踏み出してしまった? 演算時間(Core i5-2500S 2.70GHz 1core, gfortran, 最適化オプションなし, Brio&Wu, Nx=800) 速くしようと意図してコーディングしてないので参考まで。  演算時間(s) hll hlld mm mc kr weno3 weno+3 weno5 wenoz 特性変数なし 0.47 0.51 2.04 2.46 2.52 2.36 2.30 3.46 3.45 特性変数あり 3.07 3.22 3.88 4.01 3.95 5.44 5.51