Mathematical.jp

よいこの低学年向けすうがくひろば



偏微分方程式


偏微分方程式とは

線形偏微分方程式

2つ以上の独立変数とその偏導関数を含む微分方程式のことを偏微分方程式といいます。階とはその編微分方程式において含まれている最も高い値(最大の作用素)の階数に合わせて階偏微分方程式という言い方をし、このとき変数が2変数の場合についての係数をそれぞれ、そして未知関数を、既知の関数としてを定義すれば、2階の偏微分方程式は、

また、上記の式は常微分方程式のときと同様に偏微分方程式内における未知関数とその導関数、

に関して一次式であれば線形であるということを意味します。

通念的に解の偏微分方程式における一般解とは個の任意関数を含む解のことをいいます。これが常微分方程式の場合と違うところになります。

また上記偏微分式の既知の関数のときを同次、それ以外のときを非同次といいます。

2階線形偏微分方程式の一例

2階線形偏微分方程式の重要なものとしては以下のようなものがあります。

それぞれの偏微分方程式の特徴

  • 双曲型偏微分方程式
    • 波動方程式
    • 変数として時間変数tを含み、音や光などの波動を扱い次の方程式によって与えられる
    • 移流方程式
    • 速度場の下で保存スカラー量を記述する
  • 放物型偏微分方程式
    • 拡散方程式(熱伝導方程式)
    • 変数として時間変数を用いて与えられている領域内において時間とともに変化する物体の場を記述し以下の式で与えられる
  • 楕円型偏微分方程式
    • ラプラス方程式
    • 重力場や静電場ポテンシャルを記述しその解理論はポテンシャル論といい、の解では調和関数を与える
    • ポアソン方程式
    • ラプラス方程式の一般化による
    • ヘルムホルツ方程式
    • 電磁波放射など、次の方程式によって与えられる

境界値問題とは

一般的に無数の解を持つ微分方程式においてその定義されている領域の境界上で解、またはそれの導関数の値に対し何らかの条件を課すことによって解を指定します。この条件を境界条件と呼び、この境界条件を満たす解を求める問題を境界値問題といいます。

またさらにほかの場合においては、例えば時間が一つの変数として与えられている場合、あらかじめ定められている条件 ━ ある一定の初期条件によって決定される解を求める場合のものもあります。

解の重ね合わせ

同次型線形偏微分方程式において、ある領域内において関数が同時線形偏微分方程式の解ならば、

ラプラス方程式

それぞれの座標とした3次元座標空間において2階の偏微分作用素をとし、この作用素を次のようにおきます。

ここでナブラと呼ばれる次のもの、

これを使って表せば上記式の2階偏微分作用素式の右辺は次のようなものになります。

時間に依存しないある座標系におけるを変数とする関数をとおくとラプラス方程式は以下のようになります。

境界条件を課したある領域内においてなる条件を満たす解をラプラス方程式によって求めることができ、さまざまな解析を導くことが可能になります。

これにより、物体の温度分布などの状態を記述するのに重要な役割をします。

一階偏微分方程式

未知関数をとした場合これに対する一般的な1階偏微分方程式は次の形に書かれます。

この解を考える上で、1階偏微分方程式を、

としたときその中に含まれる独立な任意定数と置くと、

となるような解を完全解といいます。

例題

変数分離を行ってそれぞれ定数をと置きます。

まずより、

次にに関して、

これにより解は、

以上に挙げたものが偏微分方程式に関する簡単な一通りのものになります。このカテゴリでは境界値に関する問題、ラプラス方程式、その他1階偏微分方程式、定数係数偏微分方程式など、順次取り上げていく予定になっています。

境界値に関する問題

定義している領域の境界上で解を指定することができます。この条件を境界条件と呼び境界条件を満たす解を求める問題を境界値問題といいます。

ラプラス方程式

特殊な境界条件を与えた長方形プレートの平衡温度分布に関しての振る舞いを、2次元のラプラス方程式によって考察していきます。

Title Text

Lorem ipsum dolor sit amet, consectetur adipiscing elit, sed do eiusmod tempor incididunt ut labore et dolore magna aliqua.

Title Text

Lorem ipsum dolor sit amet, consectetur adipiscing elit, sed do eiusmod tempor incididunt ut labore et dolore magna aliqua.

コリオリの力とは

フランスの科学者で軍属でもあったガスパール・ギュスターブ・コリオリ

初歩的な力学の分野において慣性系に関する話の中にコリオリの力というものがあります。
この“コリオリ”とは人の名前であり地球が回転することによっておこる見かけの運動力を、回転座標上で移動したときの移動方向と垂直な方向に受ける慣性力の一種を数式で表現したものになります。

一般的にこのコリオリという人物は科学者という記述が多いのですが実は軍人でもあったことはあまり知られていないようです。

北半球において大砲を撃ったとき、その弾道が標的よりもわずかに右にずれることを数式によって説明した人物ということで自衛隊でも多少知られた存在らしく、話の小ネタとして話題になることもありました(おそらく重迫撃砲(ヘヴィモータ)などではコリオリの力による弾道補正が必要になるためだと思います)。

回転座標で運動する物体

円運動の角度を、周囲に沿った距離を、ボールの速度を、慣性力をとおくと、

上の図は、矢印の通りに回転させた円盤を真上からみたものであり、中心から外側へ投げた球があたかも右のほうへずれたように見えるその様子を表したものです。中心部分にいる人物は図のの方向へ投げたつもりがのほうへまるで曲がって投げたように見えてしまいます。

はは動径方向に垂直な方向に働く慣性力とし、一定加速度での移動距離はの形で表せるので、

この力をコリオリの力と呼び、回転座標系で運動する物体に加わる慣性力のことを言います。

この力を実際に数式を使って具体的に表現してみましょう。

左の座標系が3次元での回転、右の座標系の図が2次元での回転を表したものになります。

これより回転座標系において時間変化した角度をとすれば以下のように表せます。

または、

これを時間で微分します。

上記式を再度時間微分します。

の間の関係は次式で表されます。

さらに系(慣性系)では次のような運動方程式、

が成り立つので上の式の結果を用いてを表せば、

これらを代入し系の座標について整理すれば、

は回転系座標から見た加速度運動で、運動の原因となる力としてのほかに2つの力が加わった形であり、右辺第2項、第3項は見かけの力(慣性力)を示しています。

右辺のそれぞれの意味は、

右辺第2項(コリオリの成分)
右辺第3項(遠心力の成分)

を示しています。

回転座標系をまとめてみると次のようになります。

慣性系に対して運動する座標系

1.等速度運動する座標系は慣性系となる(ガリレイ変換)
2.加速度運動する座標系は非慣性系(見かけの力=慣性力)を考える
  • 加速度αで運動している系見かけの力見かけの力

3.回転運動をする座標系は2種類の見かけの力を考える
  • 遠心力:
  • コリオリの力:

これらを踏まえて今度は長距離射撃におけるコリオリ弾道軌道計算を考えていきましょう。
次回に続きます。


PAGE TOP