- 締切済み
Excel VBAでポテンシャル流れの計算プログラムを作成したい
- ExcelのVBAを使って流体のポテンシャル流れを差分法で計算するプログラムを作成しています。完成しているのですが、正しく動作しているか不安です。
- ポテンシャル流れのプログラムを作成したいので、ExcelのVBAを使用し、差分法を利用して計算しています。完成しているが、正確性に不安があります。
- ExcelのVBAを使って流体のポテンシャル流れを差分法で計算するプログラムを作成中です。完成していますが、動作の確認をお願いできますか?
- みんなの回答 (2)
- 専門家の回答
みんなの回答
関連するQ&A
- 循環とポテンシャル流れ
循環関連でやや混乱してしまったため質問させていただきます. 完全流体理論において,ダランベールの背理が示すように,循環Γを与えないと揚力が発生しません. それに関連して以下のような問題が出ました. 【問題】 ある二次元非圧縮一様流のポテンシャル流れW=Φ + iψ(Φおよびψは具体的に与えられている)に原点中心に半径aの円柱が置かれている. 今,流れ場全体に原点中心に時計回りの循環Γをさらに加えたとき,円柱に発生する揚力を求めよ. 【僕の考え】 円柱表面の流速がわかればベルヌーイ式より,円柱表面のpがわかるので それを円柱表面に沿って積分してやれば揚力が計算できる. よって,流速を知りたいから,重ね合わせの原理を用いて,循環Γを与える複素ポテンシャルを求め,(Φ’+iψ'とする) ∇(Φ+Φ’)よりvを求めて計算してやればいいと考えました. 【疑問点】 循環の定義は-Γ=∫VdR = ∫ω dS (時計回りなので) V R ω S はベクトル Rは線積分 Sは面積分 Vは流速 ωは渦度 であるはずだ. 疑問1.-Γ=∫VdRだが,dRは円周方向にそう積分なので ∫VdR=∫VθdR = 2pirVθ :VθはV(r,θ)のθ成分 であるはずだが,このrはaにすべきなのか,変数rなのかあるいは別なのかがよくわかりません. (つまり,流れ場全体にΓをあたえたということは,どの閉曲面で積分しているのかということがわからないということです) 疑問2.-Γ=∫ωdSだが,そもそもポテンシャル流れなのでωは今0ではないのか? そうすると右辺=0となるため,循環Γを与えたという事実に反するはず. しかし,定義が間違っているということもないはずです. これはいったいどういうことなのでしょうか? 混乱してしまって伝わりにくいかもしれませんが どうぞよろしくお願いいたします.
- 締切済み
- 物理学
- 平板を過ぎる流れについて
理想流体中の円柱を過ぎる流れを ジューコフスキー変換して平板を過ぎる 流れを求めると、 W=U(z*cosα-i√(z^2-4*a^2)sinα) と計算できました。 しかし、α=±π/2を代入すると実部が、 すなわち速度ポテンシャルが恒等的に0、 α=0,πを代入すると、虚部が、 つまり流関数が恒等的に0となって しまいます。 ジューコフスキー変換において、 α=±nπ/2を禁止する制約はなかったと 思いますが、なぜこのような矛盾が 生じてしまうのでしょうか?
- ベストアンサー
- 物理学
- 2次元ラプラス方程式を差分法で解くというプログラムなのですが、
2次元ラプラス方程式を差分法で解くというプログラムなのですが、 ・プログラムの流れ ・具体的にどのような計算をしているのか を教えていただけないでしょうか? 面倒だとは思いますがお願いします。 以下プログラムURL http://www.geocities.jp/laprog321/
- ベストアンサー
- C・C++・C#
- http://www.geocities.jp/laprog321/
http://www.geocities.jp/laprog321/ この上記URLのプログラムについて質問です。 このプログラムはラプラス方程式を差分法で解くプログラムと書いていたので、 自分はまずプログラムを理解しようとするために、 インターネットなどで、ラプラス方程式の差分法による解き方についていろいろ調べて、式を考えたりしました。 そして、このプログラムがいう差分法とは中央差分方程式を用いて解いているのかなと思ったのですが、 何卒知識が少ないもので、実際に合ってるかどうかが分かりません。 また、ラプラス方程式をコンピュータに解かせる場合、解き方の反復法として、 いろいろな種類(ガウスザイデル法、ヤコビ法、SOR法など)の方法があることがわかりました。 ここで質問させていただきます。 プログラムでは、実際にどの方法を利用してどのように計算しているのか教えてほしいです。 例えば、ここの部分(行数など)で「~式」、「~法」を用いてどのように計算してるかなど教えてほしいです。 プログラムの流れについては、違う質問でご教授いただいたので、上記した質問のように プログラム本文を使ってではなく、実際にこのプログラムでは、ラプラス方程式をどのような式を用いて、 また、どのような方法で解いているかを、矛盾しているようですが、プログラムと対比して教えていただきたいです。 同じような質問を繰り返し失礼いたします。 文章力のない質問で申し訳ありませんがよろしくお願いします。
- ベストアンサー
- C・C++・C#
- 流体解析のプログラムの作成で困っています…
現在「数値流体工学」著:荒川忠一 の本を参考に,差分法を用いた流体解析のプログラムを作成をしました. そして,解析をしたところ計算の途中で差分値が0になり,止まってしまいます. 原因を探してみたのですが,予想としましては, 1.モデルの格子幅の大きさに依存(幅が大きいと早く止まる) 2.収束率を速めるために用いている緩和係数の大きさに依存 と考えましたが,実際に何をどうすれば解決するのか という所までは考えが思い浮かばず,見当が付きません. もし,このような経験をした方,もしくはこの本をご存知の方がいらっしゃいましたら,なにかアドバイスをいただきたいのです. うまく問題事項が説明できませんので,ご質問とかありましたらお願いします. 大変困っております. ぜひ,小さな情報でも構いません.よろしくお願いします.
- ベストアンサー
- その他(学問・教育)
- 差分法の精度(電位の計算)
よろしくお願いします。 現在、差分法により電位を求めています。 ポアソン方程式を差分化してプログラムにより電位を求め、それと電位の式の計算結果をと比較しているのですが、差分法の計算結果があまり近い値となりません。 差分法の精度とはどのくらいなのでしょうか? 電荷がない場合では完璧に一致するのですが、電荷があると一致しません。ポアソン方程式の一般的な解はまだ導き出されておらず、電荷も常に一定ではないので十分な精度はない、みたいなことをちらっと聞いたこともあるので、もしやと思ったのですが、そもそも差分法の電位の計算とは一致しないものなのでしょうか? よろしくお願いします。
- ベストアンサー
- 物理学
- 井戸型ポテンシャルの数値解シミュレーション
井戸型ポテンシャルの数値解を「オイラー・クロマー法」を用いてシミュレーションで表す場合、どのようなプログラムにすれば良いのでしょうか?井戸の深さや幅をこちらで任意な数値を代入して、基底状態などにおける固有値や波動関数をグラフに表したいのです。 プログラムに関しては全くの初心者なのでどなたか具体的に教えてもらえれば幸いです。言語はC、グラフィクスの処理はPGPLOTのようなものは使用しないで、単純に結果をエクセルへ出力して求めたいです。
- ベストアンサー
- 物理学
- 1次元の井戸型ポテンシャル
以下のような1次元の井戸型ポテンシャル V(x)=0 (-L<x<L) =V (x<-L, x>L) ただし 0<E<V 中の質量mの粒子について・・・・ この問題でグラフの交点を求めることによって、固有値が求められますが単純な計算では出せません。よって、以下のような課題を出されました。 上の問題で適当なVとLについて、固有値をニュートン法などの簡単なプログラムを組んで、数値的に計算して求めよ。また、そのときの固有関数を求めてプロットせよ。 以上の問題なのですが、簡単なヒントなのでよろしいので分かる方がいたら教えてください。お忙しい中ありがとうございました。
- 締切済み
- 物理学
- 2階微分の差分法
こんにちわ。 現在、差分法にてプログラムを作っています。 やりたいことは、場所によって空間分解能を変え、重要な部分は細く計算し、そうでないところは分解能を下げて粗く計算したいと思っています。 その際、2階微分の空間分解能が変わる境界での差分はどのように書けばいいのでしょうか? 何か確立された方法をご存知であれば教えて下さい。 よろしくお願いします。
- ベストアンサー
- その他(プログラミング・開発)