[論文レビュー] Shoreline and Bathymetry Approximation in Mesh Generation for Tidal Renewable Simulations

Alexandros Avdis, Christian T. Jacobs+3 more|arXiv (Cornell University)|Oct 6, 2015

Tropical and Extratropical Cyclones Research参考文献 13被引用数 4

ひとこと要約

本論文は、潮流発電シミュレーションの非構造格子生成のための海岸線および水深の簡素化に、主成分分析（PCA）に基づく新規手法を提案する。地理的データを支配的空間モードに分解することで、計算コストを低減しつつ主要な大規模特徴を保持する滑らかで再現可能な簡素化が可能となる。英国の沿岸域、特にオキニーア諸島およびセーヴェルン・エストゥアリー周辺の高解像度格子を用いた実証により、有効性が示された。

ABSTRACT

Due to the fractal nature of the domain geometry in geophysical flow simulations, a completely accurate description of the domain in terms of a computational mesh is frequently deemed infeasible. Shoreline and bathymetry simplification methods are used to remove small scale details in the geometry, particularly in areas away from the region of interest. To that end, a novel method for shoreline and bathymetry simplification is presented. Existing shoreline simplification methods typically remove points if the resultant geometry satisfies particular geometric criteria. Bathymetry is usually simplified using traditional filtering techniques, that remove unwanted Fourier modes. Principal Component Analysis (PCA) has been used in other fields to isolate small-scale structures from larger scale coherent features in a robust way, underpinned by a rigorous but simple mathematical framework. Here we present a method based on principal component analysis aimed towards simplification of shorelines and bathymetry. We present the algorithm in detail and show simplified shorelines and bathymetry in the wider region around the North Sea. Finally, the methods are used in the context of unstructured mesh generation aimed at tidal resource assessment simulations in the coastal regions around the UK.

研究の動機と目的

• 関心の対象とは離れた地域で過剰な格子解像度が原因で生じる潮流シミュレーションの高コスト問題に対処すること。
• 海岸線や水深の複雑でマルチスケールな幾何構造を、強固で再現可能な方法で簡素化する手法の開発。
• 幾何的基準に基づく簡素化手法で一般的に見られる非滑らかでぎざぎざしたアーチファクトを回避することで、格子品質の向上。
• 簡素化された海岸線と水深を一貫して結合し、潮流資源評価のための正確な格子生成を可能にすること。
• 潮流エネルギー立地の水理的シミュレーションに用いられる非構造格子の自動的かつスケーラブルな格子生成を支援すること。

提案手法

• 主成分分析（PCA）をベクトル化された海岸線データおよびラスタライズドされた水深データに適用し、支配的構造を表す直交空間モードに分解する。
• データの分散を最も多く捉える主要な主成分のみを保持することで、低ランク近似を実現し、幾何構造を簡素化する。
• 簡素化された幾何構造は、モードの低次元集合から再構築され、滑らかさと大規模特徴の保持が保証されつつ、微小スケールのノイズがフィルタリングされる。
• 本手法は、海岸線の輪郭（ベクトルデータ）と水深ラスタ（ラスタデータ）の両方に別々に適用され、スケールと解像度の取り扱いが一貫している。
• 格子生成はGmshを用い、関心領域付近では細かい（50 m）から遠方では粗い（10 km）に段階的に要素サイズを調整する。
• Fluidity、OpenTidalFarm、Telemac、MIKEといった主要な海洋モデリングフレームワークとの統合が可能で、シームレスなシミュレーションワークフローを実現する。

実験結果

リサーチクエスチョン

• RQ1 PCAは、海岸線や水深といった複雑でマルチスケールな海岸幾何構造を、数学的に堅牢かつ計算的に効率的に簡素化するために効果的に適用可能か？
• RQ2 PCAに基づく簡素化は、従来の幾何的またはフィルタリング手法と比較して、大規模特徴の保持と微小スケールノイズの低減の両面で優れているか？
• RQ3 PCAに基づく簡素化は、潮流発電シミュレーションの非構造格子生成において、格子品質の向上と計算コストの低減にどの程度寄与するか？
• RQ4 本手法は、一貫したフレームワーク内で、ベクトル（海岸線）およびラスタ（水深）データの両方のデータタイプに適用可能か？
• RQ5 保持する主成分の数の制御が、幾何的忠実度と格子生成における計算効率のトレードオフにどのように影響するか？

主な発見

• PCAに基づく簡素化手法は、少数の主成分のみで海岸線および水深の支配的で大規模な特徴を的確に捉え、幾何構造の有効な低減を実現した。
• PCAを用いて生成された簡素化された海岸線および水深は、滑らかで連続的な輪郭を示し、幾何的基準に基づく簡素化で一般的に見られるぎざぎざしたアーチファクトを回避した。
• 簡素化された領域境界から生成された格子は、オキニーア諸島やセーヴェルン・エストゥアリーのような主要地域で高解像度（50 mの辺長）を達成し、最も細かい格子では最大340万頂点に達した。
• 本手法により、粗い（4,307頂点）、中程度（75,368頂点）、細かい（340万頂点）の3段階の異なる格子グレーデーションが作成可能であり、スケール間で一貫した境界表現が保証された。
• 本手法は、地理的データの再現可能でアルゴリズム的な処理を可能とし、格子生成パイプラインにおける手動介入への依存を低減した。
• 実装は、高解像度水深データと大域的な低解像度データの統合といった、マルチスケールデータ統合への拡張性を有する。