ノイマン境界条件のフーリエ変換


8

2つの連立偏微分方程式系を数値的に解く必要があります。

x1t=c12x1+f1(x1,x2)x2t=c22x2+Kx1t

システムのドメインは正方形の領域です。

境界条件:

x=constantx1x=x2x=0y=constantx1y=x2y=0

このシステムをフーリエ変換で解いてみました。数回反復すると解が不安定になります。私は以前にこのシステムを有限差分法で解決しましたが、うまく機能したので、システムの定数は完全に問題ありません。

  • 私の質問は、フーリエ変換を使用してこれらの方程式を解くことができるかどうかです。
  • ノイマン境界条件のためにフーリエ変換を適用できないとどこかで読んだことがあります。これは正しいです?
  • はいの場合、代替手段は何ですか?

どのように定義しますか?ある周期関数?f 1f1(x1,x2)f1
ucsky 09/10/12

f1(x1,x2)=P(x1)arctan(x2)、ここでP(x)は多項式です。
チャトゥール

1
おそらく、数値の不安定性は関数ます。ドメインの制限でがゼロに等しくない場合、周期信号を再構築するときに境界条件でジャンプする可能性があります。 Pf1P
ucsky 2012

@異常、コメントありがとうございますが、それが何を意味するのか理解できません。たぶんFFTをもっと徹底的に勉強するべきだろう。しかし、すべての変更をシステムで解決している場合、境界の値が影響を受けないように、変更は中央部分のみに(大まかに)制限されている場合、ここでFFTを使用できますか?
チャトゥール

1
まあ私が言ったことは間違っているかもしれません、それはいくつかの滑らかな関数うまくいくかもしれませんが、とにかくそれは境界付近のプロセスが適切に解決されないのでおそらくそれはおそらく間違った解決策を与えます。おそらく、コサインモードまたはチェビシェフモードを使用する必要があります。本当にフーリエ変換を使いたいなら、ある種のペナルティ法を使うことができます。f1
ucsky 2012

回答:


4

FFTは、周期的な境界条件に使用できます。フォンノイマン境界条件は事実上「ミラー」境界条件であるため、FFTを適用する前に「ミラー継続」を行う必要があります。このアプローチの1つの欠点は、データ量を4倍に増やすことです(少し実験するだけの場合は重要ではありません)。コサイン変換を使用すると、暗黙的に「ミラーリングされた継続」が行われ、ファクター4のオーバーヘッドが回避されます。

境界の近くのグリッドポイントが配置されている場所に応じて、「離散ミラー継続」を実行するには2つの異なる方法があることに注意してください。したがって、FFTWのようなライブラリは、コサイン変換のさまざまなバリアントを提供していることがわかります(これらのさまざまな「離散ミラー継続」に対応しています)。

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.