本研究では、NURBS ベースの等幾何解析と2次オーダーのStrang オペレータ分割法を組み合わせた手法を提案している。
まず、移流項はセミラグランジュ法で扱う。これにより、メッシュが歪むことなく固定されたままで計算を行えるため、リメッシュの必要がなく効率的である。
次に、拡散項は陰的な2次BDF法で離散化し、反応項は明示的なルンゲ・クッタ法で解く。これにより、非線形性への対処が容易になる。
提案手法の精度を検証するため、解析解を持つ問題や、Schnakenberg-Turing モデル、Gray-Scott モデルなどの非線形拡散反応方程式系を用いた数値実験を行った。その結果、複雑な幾何形状上でも高精度な解が得られ、パターン形成を正確に再現できることが示された。
特に、Schnakenberg-Turing モデルでは、幾何形状がTuring パターンに与える影響を明らかにできた。
다른 언어로
소스 콘텐츠 기반
arxiv.org
더 깊은 질문