OKWAVEのAI「あい」が美容・健康の悩みに最適な回答をご提案!
-PR-
解決
済み

偏微分方程式

  • すぐに回答を!
  • 質問No.104706
  • 閲覧数88
  • ありがとう数4
  • 気になる数0
  • 回答数3
  • コメント数0

お礼率 53% (212/400)

ある大学院の過去問なんですが、小問が3つあるんですが、まず1つ目について教えてください。

『次のような u(x, t) に関する偏微分方程式を考える。
    ∂u/∂t + u ∂u/∂x = ν ∂^2 u/∂x^2    (*)
ここで、νは正の定数である。この時以下の問いに答えよ。

(1) (*)式に ψ(x, t) と φ(x, t) を用いた2段階の変換
    i. u = ∂ψ/∂x
    ii. ψ = αlogφ
を行ない、定数αを適当に選ぶと φ(x, t) に関する線型方程式
    ∂φ/∂t = ν ∂^2 φ/∂x^2    (**)
が得られることを示せ。また、そのときのαを求めよ。』

と言う問題です。

変換i. ii. を行なってφの偏微分方程式には出来るのですが(合ってるかどうかは別問題)
αをどう取っても(**)になりそうにないんです。(ってことは合ってないって事?)

やったのは
    u = ∂ψ/∂x = (∂/∂x)(αlogφ) = (α/φ)(∂φ/∂x)
    ∂u/∂t = (∂/∂t){(α/φ)(∂φ/∂x)} = -(α/φ^2)(∂φ/∂t)(∂φ/∂x) + (α/φ)(∂^2 φ/∂t∂x)
    ∂u/∂x = (∂/∂x){(α/φ)(∂φ/∂x)} = -(α/φ^2)(∂φ/∂x)^2 + (α/φ)(∂^2 φ/∂x^2)
    u(∂u/∂x) = (α/φ)(∂φ/∂x){-(α/φ^2)(∂φ/∂x)^2 + (α/φ)(∂^2 φ/∂x^2)}
        = -(α^2/φ^3)(∂φ/∂x)^3 + (α/φ)^2(∂φ/∂x)(∂^2φ/∂x^2)
    ∂^2u/∂x^2 = (∂/∂x){(α/φ)(∂φ/∂x)} = -(2α/φ^2)(∂φ/∂x)^2 + (α/φ)(∂^2φ/∂x^2)
までなんですが、これを(*)に代入するとすごい事になってとても(**)にたどり着けそうにないんです。

ここまでで既に間違ってるんでしょうか?それともこの状態でαを適当に選べば(**)が導けるんでしょうか?

よろしくお願いします。
通報する
  • 回答数3
  • 気になる
    質問をブックマークします。
    マイページでまとめて確認できます。

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

  • 回答No.1

とりあえずできるところまでで、計算してみました。
偏微分の記号を(1)のように表現させていただきます。
(∂/∂t)=(Dt)
(∂/∂x)=(Dx)  (1)

問題の偏微分方程式は、

[Dt+u*Dx-ν*(Dx*Dx)]u=0 (2)

(2)式に、問題のヒントとして与えられている、変換式(3)を代入します。

u=(Dx)ψ (3)

すると(2)式は
(Dt)*(Dx)ψ+[(Dx)ψ][(Dx*Dx)ψ]-ν*(Dx)*(Dx)*(Dx)ψ=0

上式を整理すると

Dx*((Dt)ψ+0.5*[(Dx)ψ]^2-ν*(Dx*Dx)ψ)=0 (4)

(4)式を解くと
(Dt)ψ+0.5*[(Dx)ψ]^2-ν*(Dx*Dx)ψ=g(t) (5)

ここで、g(t)はtの任意の関数(当然0という定数でもよい)。
(5)式に変換式(6)を代入します。

ψ=α*log(φ) (6)

ここから先の計算ははしょりますと

α=2*νとすると、

(Dt)φ=ν*(Dx)*(Dx)φ+g(t) (7)

という式が導出できます。
問題では、g(t)=0の特別な場合になっているようです。

誤記、誤解があったらゴメンなさい。
-PR-
-PR-

その他の回答 (全2件)

  • 回答No.2

No1でアドバイスさせていただきましたblue_monkeyです。 taropoo氏の質問を読み直すと、質問に対する直接的なアドバイス になっていないようです。 blue_monkeyも、始めにtaropoo氏のされたような方法で試みようとしたのですが、見通しが悪く、すごいことになりそうだと思いまして、少しやり方を変えてしまいました。このため、taropoo氏の質問からはずれてしまったような気がしま ...続きを読む
No1でアドバイスさせていただきましたblue_monkeyです。
taropoo氏の質問を読み直すと、質問に対する直接的なアドバイス
になっていないようです。
blue_monkeyも、始めにtaropoo氏のされたような方法で試みようとしたのですが、見通しが悪く、すごいことになりそうだと思いまして、少しやり方を変えてしまいました。このため、taropoo氏の質問からはずれてしまったような気がします。
とりあえず、問題を解く手がかりとしていただければ幸いです。
補足コメント
taropoo

お礼率 53% (212/400)

いえいえとんでもないです。
(4)への変形は目が覚める思いでした。これってさらっと積分してるんですよね?
同じく(5)も任意関数が出てくる所も納得です。
問題は(5)から(7)への変形なんですがどうもうまく行きません。

(6)より
(Dt)Ψ = (α/Φ)(Dt)Φ
(Dx)Ψ = (α/Φ)(Dx)Φ
(Dx)*(Dx)Ψ = -(α/Φ^2)(Dx)Φ + (α/Φ)(Dx)*(Dx)Φ
これらを(5)に代入すると
(α/Φ)(Dt)Φ + (1/2)(α/Φ)^2{(Dx)Φ}^2 - ν{-(α/Φ^2)(Dx)Φ + (α/Φ)(Dx)*(Dx)*Φ) = g(t)
この段階でg(t)=0としてしまい、全体を(α/Φ)で割ると
(Dt)Φ + (1/2)(α/Φ){(Dx)Φ}^2 + ν(1/Φ)(Dx)Φ - ν(Dx)*(Dx)*Φ) = 0
となるのですが、α=2νとしても
(Dt)Φ + (ν/Φ)[{(Dx)Φ}^2 + (Dx)Φ] - ν(Dx)*(Dx)*Φ) = 0
となってしまい、第2項が消えてくれません。
どこでミスってるんでしょうか?
投稿日時 - 2001-07-17 02:51:55
  • 回答No.3

(6)より (Dt)Ψ = (α/Φ)(Dt)Φ               (イ (Dx)Ψ = (α/Φ)(Dx)Φ                (ロ (Dx)*(Dx)Ψ = -(α/Φ^2)(Dx)Φ + (α/Φ)(Dx)*(Dx)Φ  (ハ と計算をされていますが、(ハの計算が異なるように思われます。 (Dx)Ψ = (α/Φ)(Dx)Φ   をさらにxで偏微分すると、 (D ...続きを読む
(6)より
(Dt)Ψ = (α/Φ)(Dt)Φ               (イ
(Dx)Ψ = (α/Φ)(Dx)Φ                (ロ
(Dx)*(Dx)Ψ = -(α/Φ^2)(Dx)Φ + (α/Φ)(Dx)*(Dx)Φ  (ハ
と計算をされていますが、(ハの計算が異なるように思われます。
(Dx)Ψ = (α/Φ)(Dx)Φ  
をさらにxで偏微分すると、
(Dx)*(Dx)Ψ=[(Dx)(α/Φ)](Dx)Φ+(α/Φ)(Dx)*(Dx)Φ
となり、
-(α/Φ^2)[(Dx)Φ][(Dx)Φ] + (α/Φ)(Dx)*(Dx)Φ  
となります。たぶんこれでうまくいくのでは?
ゴメンナサイ。
あとα=2*νではなくα=-2*νの誤記です。
            
お礼コメント
taropoo

お礼率 53% (212/400)

(Dx)*(Dx)Ψ = -(α/Φ^2)(Dx)Φ + (α/Φ)(Dx)*(Dx)Φ  (ハ

(Dx)*(Dx)Ψ = -(α/Φ^2){(Dx)Φ}^2 + (α/Φ)(Dx)*(Dx)Φ  (ハ '
の間違い、つまり
(fg)' = f'g + fg'
とすべき所を
(fg)' = f' + fg'
としてたという事ですね。

とても丁寧な解説、ありがとうございました。
今後ともよろしくお願いします。
投稿日時 - 2001-07-17 13:13:06
このQ&Aで解決しましたか?
関連するQ&A
-PR-
-PR-
このQ&Aにこう思った!同じようなことあった!感想や体験を書こう
このQ&Aにはまだコメントがありません。
あなたの思ったこと、知っていることをここにコメントしてみましょう。

その他の関連するQ&A、テーマをキーワードで探す

キーワードでQ&A、テーマを検索する
-PR-
-PR-
-PR-

特集


いま みんなが気になるQ&A

関連するQ&A

-PR-

ピックアップ

-PR-
ページ先頭へ