روش‌های تفاضل متناهی مرتبه‌بالا بر روی شبکه مکعب‌کره: کاربست به فرارفت کمیت نرده‌ای غیرفعال

نوع مقاله : مقاله پژوهشی

نویسندگان

1 گروه فیزیک فضا، موسسه ژئوفیزیک دانشگاه تهران

2 هیئت علمی موسسه ژئوفیزیک

3 مؤسسه ژئوفیزیک دانشگاه تهران

چکیده

در کار حاضر، پیاده‌سازی شبکه مکعب‌کره مرکزی و حل تفاضل متناهی مرتبه‌بالای معادله فرارفت بر روی آن ارائه شده‌اند. اگرچه این پیاده‌سازی به گونه‌ای است که به‌کارگیری عملگرهای تفاضل متناهی از نوع و مرتبه‌های مختلف را ممکن می‌کند، در این کار، صرفا به مقایسه روش‌های تفاضل متناهی صریح مرکزی مرتبه ۲، ۴ و ۶ بسنده شده است. انتگرال‌گیری زمانی به دو روش رونگه-کوتای کلاسیک مرتبه ۴ و لیپفراگ با پالایه زمانی وایازش چندجمله‌ای درجه یک با مرتبه دقت ۲ انجام شده است. این روش‌ها در معادله فرارفت با استفاده از آزمون گردش جسم صلب زنگوله‌های کسینوسی با رده‌های مشتق‌پذیری C1 و C3 مورد مقایسه قرار گرفته‌اند. روش‌های مکانی، زمانی و آزمون‌های موردی به‌گونه‌ای انتخاب شده‌اند که امکان بررسی عملکرد اجزای مختلف روش عددی و عوامل مؤثر بر دقت آنها تسهیل شود. نتایج آزمون‌ها نشان می‌دهد که اگرچه رده مشتق‌پذیری پایین یا مرتبه پایین روش انتگرال‌گیری زمانی می‌توانند مرتبه دقت دست‌یافتنی را محدود کنند، بااین‌حال به‌کارگیری عملگرهای تفاضل متناهی مرتبه‌بالا می‌تواند خواص پایستاری مدل را به‌طور قابل‌توجهی بهبود دهد. تمامی روش‌های به‌کاررفته در اینجا برای پایداری نیاز به مقداری میرایی مصنوعی مقیاس‌گزین دارند که در این کار به وسیله پالایه‌های فضایی تأمین می‌شود. نتایج آزمون بلندمدت نشان می‌دهد که روش‌های مرتبه‌بالا‌تر و تفکیک‌های بالاتر به میرایی مصنوعی کم‌تری برای پایداری محاسباتی نیاز دارند.

کلیدواژه‌ها

موضوعات


عنوان مقاله [English]

High-order finite difference methods on the cubed sphere: Application to passive scalar advection

نویسندگان [English]

  • Eliad Bagherzadegan 1
  • Alireza Mohebalhojeh 2
  • Sarmad Ghader 3
1 Department of Space Physics, Institute of Geophysics, University of Tehran
2 Space Physics Department, Institute of Geophysics, University of Tehran
3 Institute of Geophysics University of Tehran
چکیده [English]

The present work is part of our ongoing effort to make use of the cubed sphere for global atmospheric modeling applications. To this end, high-order numerical solutions of advection equation on the gnomonic cubed sphere grid are presented. Arbitrary-order one-dimensional interpolation based on spherical distance of the grid points have been used at the halo regions to bind the faces of the cubed sphere together. Although, the implementation of the model in this work is such that it makes using finite difference operators of various types and orders straightforward, we only discuss the results of 2nd-, 4th-, and 6th-order centered explicit operators. For time integration, the classical 4th-order Runge-Kutta method and the leapfrog method with a first-degree polynomial regression time filter with 2nd-order accuracy have been used. These methods have been compared by solving the advection equation using solid body rotations of C^1 and C^3 cosine bells. The long-range variant of this test is used to analyze the stability of the schemes. The methods and test cases have been chosen in such a way that facilitates the assessment of different components of the numerical scheme and factors affecting their accuracy. Results of the advection test cases, demonstrate that although low differentiability class of the advected field or low order of accuracy of the time integration scheme, could limit the actual order of accuracy of the solution, conservation properties of the model could, nevertheless, be improved significantly by the use of high-order finite difference operators. All of the methods used in this work need some scale-selective artificial dissipation for stability which has been supplied in the form of high-order explicit filters. Single-parameter fractional filters (i.e. filters that remove only a fraction of the shortest waves) are used in the long-range integrations to determine the required amount of artificial dissipation for each method. Results demonstrate that lower-order methods and lower resolutions require more artificial dissipation for stability. Since there is no systematic cascade of variance to smaller scales in the solid body rotation test-case, application of filter does not reduce the order of accuracy of the normalized variance. Also, the minimum required dissipation which has been calculated here might not be sufficient in cases where there is such a cascade. Time traces of 2-norm of error display a linear increase with time with no jumps or transient increases. Although such transient increases are present in time traces of normalized mean and variance, which might signify grid imprinting, this does not pose a problem since these non-conservation errors converge with or faster than the expected order. Overall, the results demonstrate the superiority of the high-order methods in terms of accuracy and performance. In long-range integrations, high-order methods exhibit convergence of errors at low resolutions even after 100 revolutions of the cosine bells.

کلیدواژه‌ها [English]

  • gnomonic cubed sphere
  • advection equation
  • finite difference method
  • filter
  • fractional filter