数値計算により電流分布を算出
再分割アルゴリズムを用いた力学系の定常分布の数値計算
... 力学系の定常状態において状態空間の各点に軌道が訪れる確率分布を表し，系の複雑な挙動を特徴づける．定常 分布の計算方法として，状態空間を有限個の小領域に分割し，系の時間発展を有限マルコフ連鎖で近似し，定常 ...
8
斜め磁界中の超伝導平板内電流分布に関する考察
... Brandtモデルで（1.1）式を満たす斜め磁界中の平板 内電流前線パターンの右側の電流領域境界を与える前 線は，５個の区 曲線を連ねてできている 。この連続 区 曲線上の座標は，Brandtモデルの式により， に 対する 座標 として与えられている。そこで，右 側前線曲線上の座標点列を数値的に記録し，それらを ...
4
圧縮性LESを用いたエアリード楽器の発音機構の数値解析 (数値解析と数値計算アルゴリズムの最近の展開)
... たす． (2) 式の左辺は密度 $\rho$ に対する波動方程式であるので，右辺の非同次項は音源と解釈する事 が出来る．具体的に，音源を計算するには， Lighthill のテンソル $T_{ij}$ を計算し，その偏微分を取 ればよい．したがって，形式的にではあるが， Lighthill の方程式を使って流体の運動がどの様に ...
11
圧縮性LESを用いたエアリード楽器の発音機構の数値解析 (数値解析と数値計算アルゴリズムの最近の展開)
... 図 2: 小型エアリード楽器のモデル ( 長さの単位は mm, 角度の単位は度 ) この小論では，小型エアリード楽器の数値解析の結果の概要を報告するだけでなく，OpenFOAM にある他の圧縮流体ソルバー (LES, レイノルズ平均モデル (RANS)) との比較を行った結果を簡 単に紹介し，なぜ LES ...
11
イオン電流計算式に含まれる誤差関数の検討
... Ａｂｓｔｒａｃｔ Ｔｈｅ ｅ ｌｅｃ ｌ 口ｏｎ ｅｎｅｒｇｙ ｄｉ ｉｏｎ ｏｆ Ａ ｉｖｅ ＣｏｌｕｍＬｎ ｏｆｔｈｅ ｇｌｏｗ ｄｉ ｔ ｓｂｄｂｕｔ ｒ 血 ｔｈｅ Ｐｏｓｉ ｓｃｈａｒｇｅ ｃａｎ ｂｅ ｍｅａｓ‐ Ｌ ’ ｉｓ ｏｆ Ｄた 則ｒｅｄ ｏｎ 出ｅ ｂａｓ ｌｅｃｔ ｒ ｒｅｎｔｓ ａｒｅ ｔｈｅ ｓｕ ｕｍ ｏｆ ｅ ｙｖｅｓｔｅ～ｒｎ ｓ Ｐ[r] ...
10
台形則により得られるscaling係数, wavelet係数の誤差解析 (数値計算における前処理の研究)
... と変形出来る. これらの式に, 積分剰余項を持つ Taylor 展開式を適用し , scaling 関数, wavelet 関数の moment に関する性質を使うと誤差評価式を得る事が出来る . 3. scahing 係数, wavelet 係数等の別表現 . 本章では, 積分剰余項を持つ Taylor 展開 (文献 [2],[3]) を使って, ...
15
H-J法を用いた3次元渦電流計算 (計算力学の新解法と領域分割法)
... $\mathrm{c}\mathrm{u}\mathrm{r}\mathrm{r}$ en $\mathrm{t}$ $\mathrm{a}\mathrm{n}\mathrm{a}\mathrm{l}\mathrm{y}\mathrm{s}\mathrm{i}\mathrm{s}$ $\mathrm{u}\mathrm{s}\mathrm{i}$ ng $\mathrm{t}$ he Nede 1ec ...
13
密度不安定性により発生する塩水振動子の数値解析(計算流体力学に関わる数理的諸問題)
... の高さの時間変化の数値計算の結果を示した。内水面の高さは徐々に上昇してゆき、 3. 6 秒で内外水面の高さが等しくなっているが、 その高さは最終的な平衡状態に落ち着いた場 合の理論値の高さと – 致する。 この時点では、 密度差の分布はまだ不均–である。 にもか かわらず、 2 次元系であるため空間的広がりがないことから、 リミットサイクル的な周期 ...
12
渦運動に対する高速数値計算法 (数値計算における前処理の研究)
... クラスタと点の「遠近」 を判断する条件を定義しなけれ ばならない。 これは近似の精度や計算量に大きく関係するので、注意深く選ぶ必要がある。 これを以後「遠近条件」 と呼ぶ。具体的な条件については ...
23
Jordan 標準形の数値計算について(科学技術における数値計算の理論と応用II)
... 分の理論」 を用いることにより、 非対称で 2 次以上の Jordan ブロックをもつ行列につい て、 その重複固有値の近似値の情報から、 より精度の高い固有値を求め、 固有ベクトル 一般固有ベクトルを求める理論とアルゴリズムを得ることが出来たのでそれについて報告 する。 ここでの特徴としては、 固有値の近似値からの出発により 2 ...
10
Bessel関数を含む振動積分に対する数値積分公式(数値計算アルゴリズムの現状と展望II)
... これから反対称積分に対する Lagrange-Bessel 積分則を導出するが、 Bessel 関数の次数 $n$ が – 般の非負整数の場合計算が多少煩雑になるので、 $n=0$ の場合についてのみ具体的に 計算を示すことにする。 Lagrange-Bessel 補間 (2.2) の両辺に $\mathrm{s}\mathrm{g}\mathrm{n}x$ ...
11
DE積分変換を利用した水面波および孤立波の数値計算について (微分方程式の離散化手法と数値計算アルゴリズム)
... 2DE 変換公式を用いた数値解法 我々はまず, この方程式の数値計算を行うにあたり , DE 変換公式と $\arctan$ 関数によっ て生或される連続的な分点を用いた . これにより $\mu$ が非常に大きな値を取った時にも安 定して , かつ, Chandler and Graham ...
11
電子顕微鏡用電子銃の電位分布及び電子軌道の算出
... 電子銃の電子 光学的諸性質を論ずる時は，各電極の構造や静電界は軸 (光軸) i こ対して廻転対称と考え，主 l こ近軸光線を問題 にする.故l こ斯る近軸光線を形成する電子軌道を論ずる 際l とは先づ光軸上及びその附近の電位分布を知る事が必 要となる固 E 電位分布の算出 1.近似値計算法 上述の構造を有つ銃内の電位分布を詳細に正確に解析 的手段で求むる事は，出来ない[r] ...
10
概要 常温常圧のテトラヒドロフラン (THF) の液体について 系の作成と平衡化計算と本計算を実行し エネルギーとトラジェクトリの確認 比熱 圧縮率 動径分布関数 自己拡散係数の算出を行います 1-I. 平衡化計算エネルギー極小化 1-Il. 平衡化計算温度一定 MD 1-Ill. 平衡化計算温度
... 注意点： • 本チュートリアルでは、実施時間を短縮するため平衡化計算のステップ数を 短めに設定しています。 • 同様の理由で計算精度は落とし、ソルバ間で完全に計算条件を一致させる ことは難しいため、他のソルバで計算した結果と異なることがあります。 • 分子の種類、初期密度に応じて平衡化に必要なステップ数は変化します。 • ...
25
3次元渦電流問題の数値解法 (偏微分方程式の数値解法とその周辺II)
... 法の残差の履歴である . 縦軸は残差の対数スケールで横軸は $\mathrm{B}\mathrm{i}\mathrm{C}\mathrm{G}$ 法の反復回 数である . TEAM モデルと同様に補正を加えることによって $\mathrm{B}\mathrm{i}\mathrm{C}\mathrm{G}$ 法が収束するようになって いることがわかる . Table 2 は ...
8
数値計算法
... を最小にするようパラメータを決定する。 パラメータの数を として は自由度 = の 分布に従うことが期待される。 パラメータの誤差の推定: を最小にするパラメータ値 に対して、 を１だけ増加させる ...
28
数式処理と数値計算の融合
... 日本電信電話（株） NTT コミュニケーション科学基礎研究所 関川 浩 数式処理と数値計算双方の長所を取り入れた，信頼性が 高く効率のよい計算方法を実現する技術である数値数式 融合計算を紹介する．数式処理の長所は，数式の係数な どの入力値が誤差を含まなければ，得られる結果が常に ...
7
二次元空間内の砂山底面での圧力分布に関する数値計算 (複雑流体の数理II)
... Contact forces are calculated at each contact point, and a resulting force acting on each grain is the sum of the contact force vectors for that grain.. This enables us to solve the equa[r] ...
12
合流型超幾何関数の数値計算とその応用(科学技術における数値計算の理論と応用II)
... 特殊関数の計算は数理物理などでよく遭遇する。 Bessel 関数、 Gamma 関数 などよく出てくる関数を始めとして、 それぞれの計算法は独立に発展してきた。 ここで扱う合流型超幾何関数は次の Kummer の方程式の 1 っの解として定義さ れる。 ...
9
インピーダンス同定を伴う心電図逆問題の数値計算(数値計算アルゴリズムの研究)
... 方法である . 例題から, このアルゴリズムを用いて安定な数値解を得ることができると考えられる. しか し同定結果は不十分であった . 問題 (IEP) は–意可解ではない. このため汎関数の最小化問題に変換する にあたって汎関数の最小値が – 意であるよう設定をする必要がある . 本稿における汎関数最小化問題では , この点における考慮が不十分であるため最小化関数の – ...
8