• ベストアンサー

渦度方程式の導出 (特に非線形項)に困っています.

ナビエストークス方程式から,渦度方程式の導出を するときに非線形項の(u・∇)uに∇×(回転)を したいのですが,ベクトルの計算をどうしていいか わかりません. 本やネットで公式などで調べたのですがわからないのです.どうか教えてください.

noname#69940
noname#69940

質問者が選んだベストアンサー

  • ベストアンサー
回答No.1

(u・∇)u=1/2∇|u|^2-u×(∇×u) から ∇×((u・∇)u)=∇×(1/2∇|u|^2)-∇×(u×(∇×u)) =1/2∇×(∇|u|^2)-∇×(u×ω) |u|^2はスカラーなので上の第1項は消えて(同様の理由で圧力も消えます) ∇×((u・∇)u)=-∇×(u×ω) ここから先はお任せします。

noname#69940
質問者

お礼

すみません,わかりました. 公式 ∇(A・B)=(A・∇)B+A×(∇×B)+(B・∇)A+B×(∇×A) にA=B=uを代入すればでできました.これでなんとか導出できそうです.本当にありがとうございます.

noname#69940
質問者

補足

早速の回答ありがとうございます. あつかましいのですが,再び質問があります. (u・∇)u=1/2∇|u|^2-u×(∇×u) 上式がなぜこうなるのかわかりません. よろしければ教えて下さい.

関連するQ&A

  • 流れの中に渦が発生した場合の計算

    渦がない場合、ナビエストークス方程式を使って計算するとします。 もし渦が発生した場合、ナビエストークス方程式の中に、渦の項を追加して計算することもあるのでしょうか?

  • マクスウェル方程式から波動方程式の導出

    電磁気を知人からもらった本で勉強しているのですが,その本に、 「マクスウェルの方程式から、電磁波に関する波動方程式を導け。また、磁場ベクトルを発生させるのはどの方程式のどの項かを、理由と共に述べよ。」 という問題がありました。いろいろ調べて、なんとか波動方程式の導出はできたのですが、磁場ベクトルを発生させる項がわかりません。 わかる人がいれば、できるだけ詳しく教えてください。よろしくお願いします。

  • 円筒座標におけるナビエストークス方程式の導出を教えてください.

    デカルト座標でのナビエストークス方程式 を円筒座標系(r,φ,z)での式に変換するために 円筒座標系のナブラとラプラスを求めて, 速度を代入てみたのですが 計算しても出ない項があって困っています. 特にわからないのは r方向の式の左辺の中で4つ目の項です. -vφ^2/r  (vφは周方向の速度です.) 何かヒントやアドバイスがあれば教えてください.

  • ナビエストークス方程式

    最近ナビエストークス方程式というのを本で読んだのですが、この方程式が何を表しているのかいまいち分かりません。 よろしければ、分かりやすく教えてください。

  • 実質微分の定義

    ナビエ・ストークスの式の導出に当たり …これら3方向の運動量の式は一括してベクトルの式として書くことができる。 ∂(ρu)/∂t+∇・ρuu=-∇p+∇・τ+ρg …(*)(注:τ,g,および3つのuは太字,つまりベクトル) 左辺を展開して連続の式を用いると次式が得られる。(注:連続の式 ∂ρ/∂t+∇・(ρu)=0) ρDu/Dt=-∇p+∇・τ+ρg …(**) とテキストにあります。右辺が変わっていないので左辺がいずれも等しいはずなのですが 実質微分の定義を使ってρDu/Dtを普通の微分で表しても∂(ρu)/∂t+∇・ρuuになりません。(*)式は検査体積の運動量保存を考えることで出てきます。つまり,検査体積内の運動量変化は流入する流体の持つ運動量と圧力・応力による表面力,体積力の力積の合計として導かれます。導出過程を見ても恐らく(*)はあっているような気がします。ということは(**)が間違っている? ナビエ・ストークスの式の導出がわかる方は(*)か(**)のどちらが間違っているのでしょうか?(**)のτにニュートン流体の応力テンソルを代入するとナビエ・ストークスの式が出てくるはずなのですが,(**)が間違っていれば正しいナビエ・ストークスの式は導かれませんし,あっているなら反対に(*)が間違っているということになると思います。 それとも(*)と(**)は同値なのでしょうか。何度計算してもうまく変形できないのですが。よろしくお願いします。

  • ナビエストークス方程式中の表面張力項について

    初めて質問させてもらいまっす。 さて、非圧縮性流れの運動量保存式であるナビエストークス方程式中の表面張力項についての質問です。 この表面張力項は、(表面張力係数)*(曲率)*(界面法線方向ベクトル) / (密度) として定義されておりますよね。。 速度の時間微分項の次元はm/s^2に対して、表面張力項の次元がm^2/s^2になってしまいます。ここで示した表面張力項の記述が間違っているんでしょうか?それとも私の勘違いでしょうか?表面張力係数の次元はkg/s^2、曲率は1/m、法線ベクトルは無次元、密度はkg/m^3として計算しました。おそらく私の勘違いだと思いますが、どなたか、 お返事待ってます。

  • ナビエ・ストークスの方程式について

    ナビエ・ストークスの方程式はそれぞれ慣性項、体積力項、圧力項、粘性項によって成り立っていますが、それぞれはどうやって求めるのですか?慣性項の出し方、体積力項の出し方、圧力項の出し方、粘性項の出し方をそれぞれ教えていただきたいです。一生懸命勉強していますがよくわかりません。お願いします。

  • 円柱座標系によるナビエストークス運動方程式

    円柱座標系においてナビエ・ストークスの方程式 を出そうと思うのですが、 r成分では、たとえば τθθ/r の項(粘性項)が出ません。 扇形の収支式から考えているのですが、応力によるτ がうまく本と一致しません。 いろいろな文献を見てみたのですが、 ほとんどが直交座標系によるもので書いてありません。  説明するのがやや難しいかもしれませんが、お願いします。

  • 微分方程式 線形 非線形 その3

    y’+(x/y)=1が非線形であることは理解できました。 また、y’+y^2+xy=1も非線形ですね。 前回の質問で、 >(u-1)(v^2+v+1)w が、{ u,v } について 3 次であることも解りますか? は3次であることは理解できます。 { u,v,w}については4次ですね。 yについての微分方程式は、yについての次数を見れば、線形なのか 非線形なのかわかると思います。 例えば、 y’+yx^2=1やy''+2x^2=yなどは、線形微分方程式ですね。 ここまでで間違いはありますでしょうか?何度もすいませんm(__)m 間違いがある場合はご指摘よろしくお願い致します。

  • ナビエストークス方程式について

    こんばんは! 宿題で 『ナビエストークス方程式を3次元で各成分に展開せよ Du/Dt=-1/ρgradp+ν∇^2u』 という問題がでました。 手書きですが一応式を画像で貼っておきます。 同じような質問があって回答を見たのですが、私にはわかりませんでした。 できればわかりやすく回答よろしくお願いします!