HSMAC法（Highly Simplified Marker And Cell method）で3次精度風上差分を用いてNavier-Stokes方程式を数値的に解いてみました。今回は高次の風上差分を入れてみました。以下も参照してください。 計算条件ですが、正方形のサイズは1×1、上部の壁は右へ速さ…

フラクショナルステップ法で3次精度風上差分を用いてNavier-Stokes方程式を数値的に解いてみました。今回は高次の風上差分を入れてみました。以下も参照してください。離散化は以下を見てください。河村哲也著『流体解析の基礎』のpp.75-76も非常に参考にな…

SMAC法（Simplified Maker And Cell method）で3次精度風上差分を用いてNavier-Stokes方程式を数値的に解いてみました。今回は高次の風上差分を入れてみました。以下も参照してください。離散化は以下を見てください。河村哲也著『流体解析の基礎』のpp.71-7…

HSMAC法（Highly Simplified Marker And Cell method）で風上差分を用いてNavier-Stokes方程式を数値的に解いてみました。前回は中心差分で解いたので、今回は風上差分を入れてみました。以下を参照してください。 計算条件ですが、正方形のサイズは1×1、上…

HSMAC法（Highly Simplified Marker And Cell method）を用いてNavier-Stokes方程式を数値的に解いてみました。テスト問題としてよく取り上げられる2次元キャビティー流れです。水で満たされた正方形があり、上部の壁だけ右へ一定の速さで動いていきます。最…

フラクショナルステップ法で風上差分を用いてNavier-Stokes方程式を数値的に解いてみました。前回は中心差分で解いたので、今回は風上差分を入れてみました。以下を参照してください。 離散化は以下を見てください。河村哲也著『流体解析の基礎』のpp.75-76…

フラクショナルステップ法を用いてNavier-Stokes方程式を数値的に解いてみました。テスト問題としてよく取り上げられる2次元キャビティー流れです。水で満たされた正方形があり、上部の壁だけ右へ一定の速さで動いていきます。最初は静止していた水が、上部…

フラクショナルステップ法（部分段階法）によるNavier-Stokes方程式の離散化について説明していきます。この手法は流速と圧力を直接未知数として計算する手法です。SMAC法では中間流速（予測子）を計算する際に圧力を使いましたが、フラクショナルステップ法…

SMAC法（Simplified Maker And Cell method）で風上差分を用いてNavier-Stokes方程式を数値的に解いてみました。前回は中心差分で解いたので、今回は風上差分を入れてみました。以下を参照してください。離散化は以下を見てください。河村哲也著『流体解析の…

SMAC法（Simplified Maker And Cell method）を用いてNavier-Stokes方程式を数値的に解いてみました。テスト問題としてよく取り上げられる2次元キャビティー流れです。水で満たされた正方形があり、上部の壁だけ右へ一定の速さで動いていきます。最初は静止…

MAC法（Maker And Cell method）を改良した、SMAC法（Simplified MAC法）によるNavier-Stokes方程式の離散化について説明していきます。この手法は流速と圧力を直接未知数として計算する手法です。MAC法ではPoisson方程式の右辺を計算するのに負荷がかかって…

MAC法（Maker And Cell method）で風上差分を用いてNavier-Stokes方程式を数値的に解いてみました。前回は中心差分で解いたので、今回は風上差分を入れてみました。以下を参照してください。 離散化は以下を見てください。河村哲也著『流体解析の基礎』のpp.…

MAC法（Marker And Cell method）を用いてNavier-Stokes方程式を数値的に解いてみました。テスト問題としてよく取り上げられる2次元キャビティー流れです。水で満たされた正方形があり、上部の壁だけ右へ一定の速さで動いていきます。最初は静止していた水が…

Navier-Stokes方程式に対する数値解法である、SIMPLE法とMAC法の解説リンクを見つけました。細かい離散化まで載っていてためになります。https://www.iitg.ac.in/director/files/Chapter-6.pdfこちらはMAC法だけ。スタガード格子における補間法がわかりやす…

今回は1次元移流方程式を有限体積法（流束制限関数）で解いていきます。C++コード付きです。1次元移流方程式とは のような一階の偏微分方程式のことをいいます。ここで、 は溶質の濃度、 は風などによる移流速度です。さて、slope limiterというのは、補間し…

今回は1次元移流方程式を有限体積法（流束制限関数）で解いていきます。C++コード付きです。1次元移流方程式とは のような一階の偏微分方程式のことをいいます。ここで、 は溶質の濃度、 は風などによる移流速度です。さて、slope limiterというのは、補間し…

今後キャビティー流れの計算をやりたいので、MAC法（Maker And Cell method）によるNavier-Stokes方程式の離散化について説明していきます。流速と圧力を直接未知数として計算する手法で陽解法です。以下の内容は『流れ解析のための有限要素法』のpp.154-155…

今回は1次元移流方程式を有限体積法（風上差分）で解いていきます。C++コード付きです。1次元移流方程式とは のような一階の偏微分方程式のことをいいます。ここで、 は溶質の濃度、 は風などによる移流速度です。移流方程式はその名前の通り、初期の波形（…

今回は1次元定常移流拡散方程式を有限体積法（ハイブリッド法）で解いていきます。C++コード付きです。1次元定常移流拡散方程式とは のような二階の常微分方程式のことをいいます。ここで、 は溶質の濃度、 は風などによる移流速度、 は拡散の効果を表す拡散…

今回は1次元定常移流拡散方程式を有限体積法（風上差分）で解いていきます。C++コード付きです。1次元定常移流拡散方程式とは のような二階の常微分方程式のことをいいます。ここで、 は溶質の濃度、 は風などによる移流速度、 は拡散の効果を表す拡散係数と…

今回は1次元定常移流拡散方程式を有限体積法（中心差分法）で解いていきます。C++コード付きです。ただし今回は風上化は入れていません。1次元定常移流拡散方程式とは のような二階の常微分方程式のことをいいます。ここで、 は溶質の濃度、 は風などによる…

有限体積法のほうもぼちぼち1次元の場合から解説していきます！今回は有限体積法の勉強で最初に解くであろう、1次元Poisson方程式を有限要素法（Galerkin法）で解くC++コードを公開します。普通の有限体積法で離散化します。二次元配列は一次元配列に収納し…

1次元Burgers方程式の解析解を求めます。『偏微分方程式の数値シミュレーション』のpp.219-220と以下のリンクを参照しています。偏微分方程式の数値シミュレーション作者: 登坂宣好,大西和栄出版社/メーカー: 東京大学出版会発売日: 1991/04メディア: 単行本…

今回は1次元Burgers方程式（バーガース方程式）を有限要素法（Galerkin法）で解いていきます。C++コード付きです。1次元Burgers方程式は、前に書いた記事「有限要素法のプログラミングを勉強するときにどの順序で学ぶのがよいか？」で紹介した、有限要素法を…

今回は1次元非定常移流拡散方程式を有限要素法（Galerkin法）で解いていきます。C++コード付きです。1次元非定常移流拡散方程式は、前に書いた記事「有限要素法のプログラミングを勉強するときにどの順序で学ぶのがよいか？」で紹介した、有限要素法を学ぶ場…

有限体積法を学びたいと考えている方はたくさんいると思いますが、如何せん有限体積法を解説している和書は少ないです（洋書まで含めればそれなりにあります）。それに加えて、ネット上にも有限体積法の情報は少ないように見受けられます。私が有限体積法を…

A. Siviglia著、"A short introduction to Finite volumes numerical methods for the solution of the Saint-Venant equations"のリンクです。浅水流方程式の数値計算手法とRiemann問題についての解説です。https://www.ethz.ch/content/dam/ethz/special-i…

夏です！Toro著"Shock-Capturing Methods for Free-Surface Shallow Flows"を読んで内容をまとめていこうと思います！夏といったら浅水流方程式でしょう！Shock-Capturing Methods for Free-Surface Shallow Flows作者: Eleuterio F. Toro出版社/メーカー: W…