著者

堀江 正信 森田 直樹 井原 遊 三目 直登

出版者

一般社団法人 日本計算工学会

雑誌

日本計算工学会論文集 (ISSN:13478826)

巻号頁・発行日

vol.2020, no.1, pp.20201005, 2020-11-13 (Released:2020-11-13)

参考文献数

26

---

メッシュは有限要素法や有限体積法で用いられる重要なデータ構造である.メッシュデータ構造はグラフと呼ばれるデータ構造の一種であるとみなせるため,メッシュを学習するために graph neural network (GNN) が広く用いられてきた.本研究では,GNN が有限要素解析の学習に有用なモデルであることを示す.提案手法では,メッシュの回転や並進に対して不変となるような形状の特徴量を入力として用いることによって学習をより効率的にしている.さらに,提案手法は学習データセットに含まれていない 100 万節点もの大きなメッシュに対しても精度を有意に落とすことなく推論ができることがわかった.

著者

松本 久也 井元 佑介 浅井 光輝 三目 直登

出版者

一般社団法人 日本計算工学会

雑誌

日本計算工学会論文集 (ISSN:13478826)

巻号頁・発行日

vol.2021, pp.20210017, 2021-10-29 (Released:2021-10-29)

参考文献数

33

---

粘性流れ問題に対し海底面形状に適合した座標変換を適用した,底面境界適合型 moving particle semi-implicit/simulation (MPS) 法を提案する.また,この座標変換の際に出現する混合偏微分の計算モデルを提案する.この混合偏微分モデルは,MPS 法の既存の微分作用素モデルを包括する形で導出する.提案手法の精度検証として,混合偏微分モデルおよび曲線座標系におけるラプラシアンモデルの収束性の評価を行う.加えて,底面が曲線的な容器内の静水圧問題を解析し,底面境界適合型 MPS 法の精度を検証するとともに,三角形状障害物を有するダムブレイク問題を解析し,動的な問題への適用性を議論する.

著者

松永 拓也 柴田 和也 室谷 浩平 越塚 誠一

出版者

日本計算工学会

雑誌

日本計算工学会論文集 (ISSN:13478826)

巻号頁・発行日

vol.2016, pp.20160012-20160012, 2016-05-11 (Released:2016-05-11)

参考文献数

20

---

本研究では, 粒子法を用いた非圧縮性流れの数値シミュレーションにおける圧力ポアソン方程式の解法として適するAMG法を用いた線形ソルバーについて検討を行った. 粒子法では, 格子法の場合と比較して, 係数行列が多くの非ゼロ要素をもち, なおかつ, 時間の進行と共に大きく変化することが特徴として挙げられるが, これに即する手法としてPA-AMG法ならびにK-cycleマルチグリッド法をBiCGSTAB法に対する前処理として適用することを考え, 数値実験に基づく性能検証を実施した. 基礎検証として2次元ポアソン方程式にMPS粒子間相互作用モデルを用いて離散化した線形方程式を対象に解析を行った. なお, 擬似乱数を用いることにより粒子分布のばらつきを模擬した. 検証の結果, 対象とする問題サイズに依らずおよそ一定のcoarsening ratio (3.9から4.0) を有する安定なコースグリッド生成が実現された. また, 求解にかかる計算時間について前処理なしBiCGSTAB法と比較したところ, 問題サイズnに対して, 前処理なしの場合ではnの1.5乗に, PA-AMG前処理付きの場合ではnの1乗に比例する傾向を示した. また, n > 10^4では前処理なしBiCGSTAB法より高速であり, n > 10^6では10倍以上の速度が得られた. 実問題に対する有用性を検証するため, MPS法を用いた2次元ダムブレイク問題の非圧縮流れの解析を行った. MPS法の計算において現れる圧力ポアソン方程式の解法として前処理なしBiCGSTAB法とPA-AMG前処理付きBiCGSTAB法を適用し, シミュレーションに要する全体の計算時間の比較を行った. その結果, 検証を行った粒子数範囲 (7000~130万) では総じてPA-AMG前処理付きBiCGSTAB法を用いた場合の方が高速であり, 最大で約16倍の高速化が得られた. 更に, 全体の計算時間に対する圧力ポアソン方程式解法が占める割合は粒子数が大きな問題に対してもおよそ一定に保たれる結果となった. 以上の結果より, PA-AMG法を前処理としたBiCGSTAB法は, 流体シミュレーションにおいて, 粒子分布の不均一化や複雑な境界形状変化に対する堅牢性を有し, 安定して高い計算速度を実現できることから, 粒子法の圧力ポアソン方程式に適する一解法であると結論付けることができた.

著者

加藤 準治 谷地 大舜 西澤 峻祐 高瀬 慎介 寺田 賢二郎 京谷 孝史

出版者

一般社団法人 日本計算工学会

雑誌

日本計算工学会論文集 (ISSN:13478826)

巻号頁・発行日

vol.2016, pp.20160001, 2016-01-18 (Released:2016-01-18)

参考文献数

41

---

本研究は, 超惰性体二相複合材料を対象に, 分離型マルチスケール解析法を活用してミクロとマクロ構造の両方のトポロジーを同時に最適化する「マルチスケールトポロジー最適化手法」の提案を行うものである. ここでは, 使用材料体積量一定のもとマクロ構造のエンドコンプライアンス最小化を目的関数として定義した. マルチスケールトポロジー最適化は, 材料挙動が複雑な先端材料の設計に有用であるとして, 現在様々な分野で着目されている手法である. しかし, このようなマルチスケール解析を活用したトポロジー最適化は, 実際に生じる非線形の力学挙動を考慮すると, 計算コストが大きくなり, さらに目的関数の設計変数に対する感度の導出が複雑なることから, ほとんどの研究では線形弾性体を仮定した単純な問題のみを扱ってきた. そこで, 本研究ではこの現状を打破すべく, 二相複合材料の超弾性体を対象としたマルチスケールトポロジー最適化手法の確立に挑戦するものである. この問題を解くために, 当該研究では二変数境界値問題の局所化解析を活用した新しい感度解析法を構築した. これにより, ミクロ‐マクロ構造ともに最適な「理想的な構造物」の設計が可能となる. 最後にいくつかの数値シミュレーションによって本手法の妥当性を検証した.

著者

西野 崇行 加藤 準治 京谷 孝史

出版者

一般社団法人 日本計算工学会

雑誌

日本計算工学会論文集 (ISSN:13478826)

巻号頁・発行日

vol.2019, pp.20190004, 2019-04-10 (Released:2019-04-10)

参考文献数

28

---

本研究は,幾何学的非線形挙動を有する構造に対して,外荷重の方向が不確かである場合を想定し,そのロバスト性を向上させるためのトポロジー最適設計法を提案するものである.ここでは,エンドコンプライアンスの期待値と標準偏差の和を目的関数とするが,低い計算コストで精度よく期待値と標準偏差の推定を行うため,エンドコンプライアンスの2次近似式を用いる手法を提案する.そして,近似式による推定値の精度検証およびいくつかの最適化計算例を示し,本手法の妥当性を検証する.

著者

小野 謙二 古賀 壱成

出版者

一般社団法人 日本計算工学会

雑誌

日本計算工学会論文集 (ISSN:13478826)

巻号頁・発行日

vol.2020, no.1, pp.20201004, 2020-11-09 (Released:2020-11-09)

参考文献数

32

---

コンピュータ、観測、シミュレーション技術などの進歩により大量のデータが生成される時代になり、データに含まれる意味や知識をえることがますます重要になっている。この論文では、与えられたデータをシンボリック回帰問題として記述し,支配方程式を見つけるプロセスを構築した。提案手法では、遺伝的プログラミングに「偏微分関数」を導入して偏微分方程式を自動生成し、生成された方程式とデータを比較評価して誤差の少ない方程式を自動抽出する。数値実験を行い、流体シミュレーションデータから支配方程式を推定し、提案手法の有効性を評価した。その結果、元の方程式が高い確率で得られ、提案された方法がデータを表すための有用なモデルを見つけるための有効なツールになることがわかった。

著者

杉原 健太 青木 尊之

出版者

一般社団法人 日本計算工学会

雑誌

日本計算工学会論文集 (ISSN:13478826)

巻号頁・発行日

vol.2010, pp.20100018, 2010-12-03 (Released:2010-12-03)

参考文献数

20

---

GPUは演算性能が高いばかりでなくメモリバンド幅も広いため, CUDAのリリースによって格子系の流体計算などさまざまな分野でGPUを用いた高速計算の研究(GPGPU)が行われるようになっている. 本研究では流体の移流現象を記述する移流方程式に対して1次~6次精度有限差分法や5次精度WENO法を適用し, GPUによる実行性能の検証および評価, 複数GPUでの新たな並列計算方法の開発, 大規模高次精度移流計算の複数GPUを用いた並列計算での高速化と強スケーラビリティの検証を行う. 本研究は東京工業大学 学術国際情報センターのスーパーコンピュータ TSUBAME Grid Clusterを利用して行われた. NVIDIAのCUDAを用いてGPU上に実装する上で, 以下の方法を導入する. Global memoryへのアクセス回数を低減するためにx, y方向のデータをSM内のShared memoryをソフトウェアマネージド・キャッシュとして利用し, 各Thread内で閉じているz方向のデータはThread内の変数(Register ファイル)に格納し, 計算の高速化を図る. 3次元的に計算空間分割しBlockの数を増やすことでSMの並列実行効率を高める. 3次元領域分割法によるGPU並列計算により並列数の増加に伴い, 境界bufferのデータサイズを減少させ通信時間を短縮させる. 計算領域を7つのKernelに分け非同期実行することによる「計算Kernel, Device-Host間通信, MPI通信」3つのオーバーラップ技法を提案する.本論文では演算密度(flop/byte)と実行性能との関係に注目した評価を導入し, 低次から高次精度の差分を用いた移流計算を例に評価の妥当性を検証した. 高次精度手法は低次よりも演算密度が高く実行性能を引き出し易いことが明らかになった. また, GPU並列計算では3次元領域分割を使用し, 7個のstreamによる非同期実行によるkernel関数, MPI, Device-Host通信のオーバーラップ技法の提案をした. これらにより, 高次精度移流計算において60GPUで7.8TFlopsという非常に高い実行性能を達成した. 本研究によって以下の新しい知見が得られた. GPUの実行性能は演算密度(flop/byte)で有効な評価ができる. Shared memoryを用いたデータの再利用により演算密度は高くなり実行性能が向上する. 1GPUの実行性能が1CPU coreに比べて数十~百倍高速であり, 複数 GPUの場合はMPI通信に加えてDevice-Host間通信も必要になるため計算時間の中で通信時間の占める割合が高く, 高速化におけるボトルネックになる. 複数GPU計算において, 3次元領域分割が有効である. 本論文で提案した技法は移流計算に限らず直交格子でのさまざまなステンシル計算(例えば拡散方程,Jacobi法等の反復解法)にも適用可能であり, 汎用性が高いといえる.

著者

太田 悠誠 野口 悠暉 松島 慶 山田 崇恭

出版者

一般社団法人 日本計算工学会

雑誌

日本計算工学会論文集 (ISSN:13478826)

巻号頁・発行日

vol.2022, pp.20220014, 2022-10-19 (Released:2022-10-19)

参考文献数

18

---

本研究では,自立性を満たす構造設計案を創出するためのトポロジー最適化法を提案する.まず,自立性を満たさない構造を検出するために,仮想的な物理モデルを導入する.次に,レベルセット法に基づくトポロジー最適化法を導入し,自立性を実現するための制約条件を仮想的な物理場によって表現する.自立性が求められる代表例として,音響クローキング設計問題を対象に最適化問題の定式化を行う.トポロジー導関数の考え方に基づく設計感度の表式を述べ,複数の数値解析例を示すことによって提案手法の妥当性及び有用性を示す.具体的には,提案する仮想的な物理場の特徴を説明し,自立性を満たさない構造が排除された音響クローキング構造の最適設計案を示す.

著者

下野 智史 児玉 斎 鈴木 克幸

出版者

一般社団法人 日本計算工学会

雑誌

日本計算工学会論文集 (ISSN:13478826)

巻号頁・発行日

vol.2017, pp.20170009, 2017-07-24 (Released:2017-07-24)

参考文献数

10

---

本研究では, ゴルフショットにおける再現性の指標として身体負荷量を提案し, 飛距離と再現性を同時に満たすためのゴルフシャフトの最適設計を試みた. 飛距離の指標となるヘッドスピードの算出には有限要素法, 再現性の指標となる身体負荷量の算出にはマルチボディダイナミクスによる逆動力学解析を用いた. 身体負荷量は, 閾値を考慮した上半身が発揮する関節トルクの総量と定義した. 最適化計算時には, 計算負荷を考慮して応答曲面法を用いた. パレート解より, ヘッドスピードを重視した場合, 身体負荷量を重視した場合それぞれの最適ゴルフシャフト及びそのときの動作が得られた. ヘッドスピードを重視した場合の最適シャフトを作成し実証実験を行った. コントロールシャフトと比較して, 飛距離が平均8yds, 再現性の評価項目とした左右方向における弾道のバラつきが平均16yds改善された.

著者

原田 隆宏 政家 一誠 越塚 誠一 河口 洋一郎

出版者

一般社団法人 日本計算工学会

雑誌

日本計算工学会論文集 (ISSN:13478826)

巻号頁・発行日

vol.2008, pp.20080016, 2008 (Released:2008-06-06)

参考文献数

27

---

本研究ではGraphics Processing Unit (GPU)上での粒子法シミュレーションをさらに高速化する手法を提案する.スライスグリッドは粒子法の計算のボトルネックである近傍粒子探索を効率化するだけでなく,計算効率も向上させるデータ構造であるが,既存研究はグラフィックスAPIを用いた実装を行なっており,一般的なストリームプロセッサ上でも有用であるかは不明確であった.そこで本研究ではまずスライスグリッドをCompute Unified Device Architecture (CUDA)を用いた実装方法を開発し,より一般的なストリームプロセッサ上での実装を示す.また本論文ではデータの時間軸上でのコヒレンシを利用したGPU上でのブロックトランジションソートを提案し,これを用いて粒子法シミュレーションのデータの空間局所性を高めて更なる高速化をはかる.そしてDistinct Element Method (DEM)を本手法を用いて実装し,近傍粒子探索を約3倍高速化し,計算全体では約1.5倍高速化した.

著者

尹 煕元 棚橋 隆彦

出版者

一般社団法人 日本計算工学会

雑誌

日本計算工学会論文集 (ISSN:13478826)

巻号頁・発行日

vol.2001, pp.20010048, 2001 (Released:2002-06-28)

参考文献数

8

---

We suggest expansion methodology for Black-Scholes equation in financial engineering field with hydrodynamics interpretation. Black-Scholes equation is famous for derivative pricing model, introducing from Ito’s process in theory of probability, and it is applying in financial area widely. Because of Black-Scholes equation based on ideal assumption, financial engineers have been devising parameters or have been creating parameter structure in the equation in order to explain financial market dynamics. They mainly have focused on diffusion term in stochastic process that is recognized as volatility on their studies. We apply the concept of advection term for financial market, and expand Black-Scholes equation. Obtaining equation can express market dynamics simply and explain derivative market distortion, volatility smile.

著者

車谷 麻緒 佐々木 浩武

出版者

一般社団法人 日本計算工学会

雑誌

日本計算工学会論文集 (ISSN:13478826)

巻号頁・発行日

vol.2018, pp.20180008-20180008, 2018-05-17 (Released:2018-05-17)

参考文献数

22

---

鉄筋コンクリートには,大小様々なひび割れが発生する.大きいひび割れは構造物の力学性能を低下させ,たとえ目に見えない小さなひび割れであっても,水やイオンの移動経路となり,鉄筋を腐食させるため,構造物の耐久性を低下させる.鉄筋コンクリートに生じるひび割れは3次元的に発生・進展を繰り返すため,実験的にも,解析的にも,精度よく捉えるのが困難という問題がある.鉄筋コンクリートに発生・進展するひび割れの3次元幾何形状を精度よく再現し,その結果をわかりやすく可視化することができれば,損傷状態の直観的な理解や,複雑な力学機構の解明に役立つと考えられる. 本論文では,物理シミュレーションを用いて,鉄筋コンクリート内部に形成するひび割れ進展挙動を3次元で詳細に再現し,その3次元幾何形状を3Dプリンタを用いて立体的に造形化する方法を提案する.物理シミュレーションには,破壊力学に基づく損傷モデルを用いて,著者らが既往の研究で開発した手法を用いる.本研究の新規性および可視化の特徴は,空気の部分に当たるひび割れを立体的に具現化し,通常は見えない部材内部のひび割れ進展挙動まで3次元で詳細に造形化する点にある. 第2節では,本研究で用いる物理シミュレーションの定式化について示す.第3節では,鉄筋コンクリートの引張と曲げに関する問題を対象に,部材内部および鉄筋まわりに進展するひび割れを再現し,3Dプリンタを用いてひび割れの3次元幾何形状を造形化した例を示す.最後に,本研究の総括を行い,課題と今後の予定について述べる.

著者

佐藤 佳代 吉川 暢宏 楊 元稙 中埜 良昭

出版者

日本計算工学会

雑誌

日本計算工学会論文集 (ISSN:13478826)

巻号頁・発行日

vol.1999, pp.19990025-19990025, 1999 (Released:2004-11-11)

参考文献数

10

---

A new learning algorithm named whole learning algorithm is proposed for the feedforward neural network. Strictly speaking, the learning of the feedforward neural network is a kind of multi-objective optimization problem to minimize the errors of outputs for all the learning data sets with respect to the amount of weight modification. All the learning data sets are simultaneously taken into account to constitute the governing equation of the weight modification, which is formulated as linear simultaneous equations with rectangular matrix of coefficients in the proposed algorithm. The solution of the equation is determined by means of the Moore-Penrose generalized inverse to deal with the rectangular matrix. The efficiency of the proposed algorithm is demonstrated through the problem to learn the nolinear behavior described by the Ramberg-Osgood model. The applicability of the proposed algorithm is investigated in problem to learn the earthquake response of RC members.

著者

広瀬 孝三郎 江戸 孝昭 玉那覇 圭一 松原 仁 原 久夫

出版者

日本計算工学会

雑誌

日本計算工学会論文集 (ISSN:13478826)

巻号頁・発行日

vol.2016, pp.20160007-20160007, 2016

---

The fracture morphologies of mud pastes show significantly complicated patterns in nature. The mud pastes have initially fluid-like properties, but gradually change from "fluid" to porous "solid" in their drying process. However, although it has been considered that this phenomena may be induced by the differences between shrinkage ratios of soil skeletons, the mechanical details remain unknown. In this work, we attempt to develop a novel numerical model based on three dimensional finite element method for crack propagation phenomena in bentonite pastes. In order to validate the proposed method, Hausdorff's fractal dimension of the numerical cracking patterns was compared with those of the experiments. As a result, the numerical results showed reasonable cracking patterns, and the fractal dimension of the patterns by numerical simulation was almost consistent with the experimental results.

著者

古口 睦士 矢地 謙太郎 山田 崇恭 泉井 一浩 西脇 眞二

出版者

日本計算工学会

雑誌

日本計算工学会論文集 (ISSN:13478826)

巻号頁・発行日

vol.2015, pp.20150002-20150002, 2015-01-30 (Released:2015-01-30)

参考文献数

24

---

構造最適化は, 数値解析による性能評価と数学的な最適化手法により, 最大限の性能を有する構造を求める手法で, 寸法最適化, 形状最適化, トポロジー最適化に大別される. このうちトポロジー最適化は, 構造の形状だけなく形態の変更も可能な最も自由度の高い手法で, 大幅な性能向上が期待できる. 構造最適化は, 当初構造問題への適用に限られていたが, 近年では様々な物理問題に適用されてきている. 流体問題への適用は, 電磁波問題などと比較するとやや古く, ストークス流を対象に流体に作用する抗力の最小化を目的とし, その形状感度を導出することにより最適形状を得る手法や, 変動拘束された境界を含む領域において粘性流体を対象とした散逸エネルギー最小化の構造最適化の手法が提案されている. しかしながら, これらの手法は, 対象とする設計領域における物体と流体の境界を変動させることにより最適構造を得る, いわゆる形状最適化の手法であるため大幅な性能向上が期待できないうえ, 最適構造が初期構造に強く依存する欠点を持つ. これに対して, 流体問題へのトポロジー最適化の展開も報告されている. その報告の例としては, 設計領域全体を多孔質体と仮定し, 物体領域と流体領域の境界で流速が零になるように仮想的な外力を与えることで物体と流体を区別し, 散逸エネルギーを最小化する最適構造を得る手法が挙げられる. トポロジー最適化の基本的な考え方は, 固定設計領域を導入し, 物体の有無を示す特性関数により構造最適化の問題を材料分布問題に置き換えることである. このため, 構造の形状だけでなく形態の変化も可能となるが, 特性関数が設計領域内のいたるところで無限小の範囲で不連続になる可能性を持つ不適切な問題(ill-posed problem)となる. この問題を解決するためには, 大域的な意味において不連続な関数を連続な関数に置き換える設計変数の緩和を行う. この緩和方法には, 均質化法や密度法などが提案されているが, どちらの手法においても, 最適構造の境界が明確にならない, いわゆるグレースケールを許容する欠点を持つ. 他方, 新たな構造最適化の手法として, レベルセット法に基づく形状最適化が提案されている. この手法では, レベルセット関数と呼ばれるスカラー関数を用いて, 設計領域中の物体の有無をその符号で示し, 零等位面を境界として表現するため, グレースケールを含まない明瞭な境界を有した最適構造が得られる利点がある. しかしながら, この方法は基本的には形状最適化の考え方に基づいており, 新たに境界が生成されるような構造の形態の変化を許容しない. この問題を本質的に解決する手法として, レベルセット法による境界表現を行いながら, トポロジー導関数の考え方に基づき設計変数を更新することにより, 形態変更を可能とした新しいトポロジー最適化の手法が提案されている. そこで本研究ではこの形態変更を可能とした新しい最適化の手法に基づき, ナビエ・ストークス方程式を支配方程式とする内部流れにおいて, エネルギー損失の最小化を目的としたトポロジー最適化の手法を構築する. トポロジー最適化は, 状態場および感度解析のための随伴場の計算に, 何らかの数値解析法を必要とする. その代表的な数値解析法として有限体積法が挙げられる. 有限体積法は, 離散化した各要素内において保存則を満たすように定式化した数値解析手法で, 状態方程式に存在する対流項や拡散項をガウスの積分定理により領域積分から境界積分に変換でき, 数値積分においても中点公式を用いる. また, 非構造格子を扱える汎用性の高さという利点もあり, 流れ場の数値解析法として広く利用されている. そこで本研究では, 状態場と感度解析には有限体積法を用いた新しい方法を開発することにより, 明確な形状表現を可能にしながら大規模問題への展開可能な方法論を構築する. 随伴方程式は連続系に基づき導出し, 有限体積法により離散化して数値安定性に優れている半陰解法のSIMPLE法(Semi-Implicit Method for Pressure Linked Equation) を用いて解析している. また, 最適化の過程におけるレベルセット関数を更新には, 時間発展方程式の数値解析法として有限体積法を適用した方法を構築する. これにより, 時間発展方程式の離散化の過程で弱形式による複雑な展開は必要としないうえ, 時間項については数値安定性のよいオイラー陰解法を適用することにより, 時間発展方程式の計算効率は状態場および随伴場と同様に向上させることができる.

著者

澤田 有弘 長濱 俊 佐々木 信也 手塚 明

出版者

日本計算工学会

雑誌

日本計算工学会論文集 (ISSN:13478826)

巻号頁・発行日

vol.2011, pp.20110003-20110003, 2011-03-18 (Released:2011-03-18)

参考文献数

48

---

本稿では,CAE(Computer-aided engineering)やSBD(Simulation-based design)に適した計算手法の開発という視点から,CAEに従事する人々がより自由に使え,詳細解析志向よりも,設計改善インスピレーションや設計の大まかな方向性・指針を得ることに注力したCAEソフトウェアの基盤技術に関する研究を報告する.そして,その構成要素の一つとして,X-FEMに基づく流体解析手法と流体構造連成解析手法の形状モデリングにおける幾何学的利点を活用した設計応用型の流体解析フレームワークを提案する. 本解析フレームワークは現段階では未完成の部分を多く含むが,物体の三次元形状測定やCADとの連結を図った解析例を示し,将来展望を含めた複数の活用例を示す.

著者

近藤 雅裕 越塚 誠一 滝本 正人

出版者

一般社団法人 日本計算工学会

雑誌

日本計算工学会論文集 (ISSN:13478826)

巻号頁・発行日

vol.2007, pp.20070021, 2007 (Released:2007-08-06)

参考文献数

10

---

粒子法は、界面を明示的に追跡することなく自由表面流れの計算ができる.自由表面を含んだマイクロスケールの流れでは、表面張力計算が重要となる.粒子法における表面張力の計算法は2種類計算されている.1つは、連続体力学に基づく支配方程式を離散化するもので、もう1つは、分子間力に似た粒子間ポテンシャル力による方法である.粒子間ポテンシャル力を用いたほうが計算効率はよい.しかし、これまで表面張力係数および接触角との関連が示されていなかった.本研究では、粒子法の1つであるMPS法において粒子間ポテンシャル力による表面張力の計算法を提案する.また、粒子間ポテンシャル力と表面張力係数および接触角との関連を示す.

著者

大村 浩之 三目 直登 浅井 光輝 磯部 大吾郎

出版者

一般社団法人 日本計算工学会

雑誌

日本計算工学会論文集 (ISSN:13478826)

巻号頁・発行日

vol.2021, pp.20210011, 2021-07-15 (Released:2021-07-15)

参考文献数

26

---

本稿では,壁領域の角を考慮し,Incompressible Smoothed Particle Hydrodynamics (ISPH) 法のような半陰解法の粒子法に適したポリゴン壁境界モデルを開発した.基本的な理論はExplicitly Represented Polygon (ERP) 壁境界モデルを踏襲しつつ,オリジナルのERPモデルで定義されている平面ポリゴンに対する壁面寄与計算を角に対しても拡張した.そのうえで,角に対して表現された壁領域の体積を角の平面角度および立体角度に基づく近似式で補正することで,計算効率を維持しながら精度の向上を図った.開発したモデルを改良ERP壁境界モデルとし,ISPH法へ導入した.開発したモデルの妥当性を検証するため,3次元の静水圧問題およびダムブレイク問題を解いた.計算結果から,改良ERPモデルを用いた場合,従来のERPモデルと比べて理論値もしくは実験値に近く,かつ滑らかな圧力分布が得られることが確認された.

著者

中山 茂

出版者

一般社団法人 日本計算工学会

雑誌

日本計算工学会論文集

巻号頁・発行日

vol.2013, pp.20130016, 2013

---

Adiabatic quantum computation has been proposed as quantum parallel processing with adiabatic evolution by using a superposition state to solve combinatorial optimization problem, then it has been applied to many problems like satisfiability problem. Among them, Deutsch and Deutsch-Jozsa problems have been tried to be solved by using adiabatic quantum computation. In our previous paper, it has been shown that the adiabatic quantum computation in Deutsch problem is modified by using a cubic step function instead of a linear step parameter. In this paper, it is proposed to solve Bernstein-Vazirani problem more efficiently by the same cubic method to obtain a solution with higher observation probability of 99.6%.

著者

今村 敬 古賀 浩二 荒牧 軍治

出版者

一般社団法人 日本計算工学会

雑誌

日本計算工学会論文集

巻号頁・発行日

vol.2001, pp.20010044, 2001

---

Object Oriented Programming (OOP) has been proposed for developing sophisticated software that needs future expansion. Successful applications of OOP in different fields such as graphic data processing and basic OS etc. is intensively reported in the literature. To apply OOP technique in softwares for structural analysis, the target structure is represented by the basic objects (beam, truss, etc.) using key concepts such as encapsulation, polymorphism and inheritance. These features facilitate easy development and management of the program. Yagawa and Kanto⁽¹⁾ proposed class structure in structural analysis program to apply OOP technique. They showed advantages of OOP for a two-dimensional linear problem. Referring their work, we have developed seismic design program for bridge, in which some new classes, class members, class member functions were developed. It has been shown that with newly developed classes, class members and class member functions and OOP techniques, the modification and expansion of program can be more easily achieved comparing with other researchers proposal, and of course traditional procedural programming technique.